Researcher profile

Guillaume Bal

Guillaume Bal contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2023arXiv

Approximations of interface topological invariants

This paper concerns the asymmetric transport observed along interfaces separating two-dimensional bulk topological insulators modeled by (continuous) differential Hamiltonians and how such asymmetry persists after numerical discretization. We first demonstrate that a relevant edge current observable is quantized and robust to perturbations for a large class of elliptic Hamiltonians. We then establish a bulk edge correspondence stating that the observable equals an integer-valued bulk difference invariant depending solely on the bulk phases. We next show how to extend such results to periodized Hamiltonians amenable to standard numerical discretizations. A form of no-go theorem implies that the asymmetric transport of periodized Hamiltonians necessarily vanishes. We introduce a filtered version of the edge current observable and show that it is approximately stable against perturbations and converges to its quantized limit as the size of the computational domain increases. To illustrate the theoretical results, we finally present numerical simulations that approximate the infinite domain edge current with high accuracy, and show that it is approximately quantized even in the presence of perturbations.

preprint2022arXiv

Magnetic slowdown of topological edge states

We study the propagation of wavepackets along curved interfaces between topological, magnetic materials. Our Hamiltonian is a massive Dirac operator with a magnetic potential. We construct semiclassical wavepackets propagating along the curved interface as adiabatic modulations of straight edge states under constant magnetic fields. While in the magnetic-free case, the wavepackets propagate coherently at speed one, here they experience slowdown, dispersion, and Aharonov - Bohm effects. Several numerical simulations illustrate our results.

preprint2022arXiv

Topological charge conservation for continuous insulators

This paper proposes a classification of elliptic (pseudo-)differential Hamiltonians describing topological insulators and superconductors in Euclidean space by means of domain walls. Augmenting a given Hamiltonian by one or several domain walls results in confinement that naturally yields a Fredholm operator, whose index is taken as the topological charge of the system. A Fedosov-Hörmander formula implementing in Euclidean spaces an Atiyah-Singer index theorem allows for an explicit computation of the index in terms of the symbol of the Fredholm operator. For Hamiltonians admitting an appropriate decomposition in a Clifford algebra, the index is given by the easily computable degree of a naturally associated map. A practically important property of topological insulators is the asymmetric transport observed along one-dimensional lines generated by the domain walls. This asymmetry is captured by a line conductivity, a physical observable of the system. We prove that the line conductivity is quantized and given by the index of a second Fredholm operator of Toeplitz type. We also prove a topological charge conservation stating that the two aforementioned indices agree. This result generalizes to higher dimensions and higher-order topological insulators the bulk-edge correspondence of two-dimensional materials. We apply this procedure to evaluate the topological charge of several classical examples of (standard and higher-order) topological insulators and superconductors in one, two, and three spatial dimensions.

preprint2021arXiv

An inversion algorithm for P-functions with applications to Multi-energy CT

Multi-energy computed tomography (ME-CT) is an x-ray transmission imaging technique that uses the energy dependence of x-ray photon attenuation to determine the elemental composition of an object of interest. Mathematically, forward ME-CT measurements are modeled by a nonlinear integral transform. In this paper, local conditions for global invertibility of the ME-CT transform are studied, and explicit stability estimates quantifying the error propagation from measurements to reconstructions are provided. Motivated from the inverse problem of image reconstruction in ME-CT, an iterative inversion algorithm for the so-called P-functions is proposed. Numerical simulations for ME-CT, in two and three materials settings with an equal number of energy measurements, confirm the theoretical predictions.

preprint2021arXiv

Multiscale Invariants of Floquet Topological Insulators

This paper analyzes Floquet topological insulators resulting from the time-harmonic irradiation of electromagnetic waves on two dimensional materials such as graphene. We analyze the bulk and edge topologies of approximations to the evolution of the light-matter interaction. Topologically protected interface states are created by spatial modulations of the drive polarization across an interface. In the high-frequency modulation regime, we obtain a sequence of topologies that apply to different time scales. Bulk-difference invariants are computed in detail and a bulk-interface correspondence is shown to apply. We also analyze a high-frequency high-amplitude modulation resulting in a large-gap effective topology topologically that remains valid only for moderately long times.

preprint2021arXiv

Topological invariants for interface modes

We consider topologically non-trivial interface Hamiltonians, which find several applications in materials science and geophysical fluid flows. The non-trivial topology manifests itself in the existence of topologically protected, asymmetric edge states at the interface between two two-dimensional half spaces in different topological phases. It is characterized by a quantized interface conductivity. The objective of this paper is to compute such a conductivity and show its stability under perturbations. We present two methods. The first one computes the conductivity using the winding number of branches of absolutely continuous spectrum of the interface Hamiltonian. This calculation is independent of any bulk properties but requires a sufficient understanding of the spectral decomposition of the Hamiltonian. In the fluid flow setting, it also applies in cases where the so-called bulk-interface correspondence fails. The second method establishes a bulk-interface correspondence between the interface conductivity and a so-called bulk-difference invariant. We introduce the bulk-difference invariants characterizing pairs of half spaces. We then relate the interface conductivity to the bulk-difference invariant by means of a Fedosov-Hörmander formula, which computes the index of a related Fredholm operator. The two methods are used to compute invariants for representative $2\times2$ and $3\times3$ systems of equations that appear in the applications.

preprint2020arXiv

Homogenization of hydrodynamic transport in Dirac fluids

Large-scale electrical and thermal currents in ordinary metals are well approximated by effective medium theory: global transport properties are governed by the solution to homogenized coupled diffusion equations. In some metals, including the Dirac fluid of nearly charge neutral graphene, microscopic transport is not governed by diffusion, but by a more complicated set of linearized hydrodynamic equations, which form a system of degenerate elliptic equations coupled with the Stokes equation for fluid velocity. In sufficiently inhomogeneous media, these hydrodynamic equations reduce to homogenized diffusion equations. We re-cast the hydrodynamic transport equations as the infimum of a functional over conserved currents, and present a functional framework to model and compute the homogenized diffusion tensor relating electrical and thermal currents to charge and temperature gradients. We generalize to this system two well-known results in homogenization theory: Tartar's proof of local convergence to the homogenized theory in periodic and highly oscillatory media, and sub-additivity of the above functional in random media with highly oscillatory, stationary and ergodic coefficients.

preprint2020arXiv

Pencil-beam approximation of stationary Fokker-Planck

Solutions of stationary Fokker-Planck equations in the narrow beam regime are commonly approximated by either ballistic linear transport or by a Fermi pencil-beam equation. We present a rigorous approximation analysis of these three models in a half-space geometry. Error estimates are obtained in a 1-Wasserstein sense, which is an adapted metric to quantify beam spreading. The required well-posedness and regularity results for the stationary Fokker-Planck equation with singular internal and boundary sources are also presented in detail.

preprint2020arXiv

Uniqueness criteria in multi-energy CT

Multi-Energy Computed Tomography (ME-CT) is a medical imaging modality aiming to reconstruct the spatial density of materials from the attenuation properties of probing x-rays. For each line in two- or three-dimensional space, ME-CT measurements may be written as a nonlinear mapping from the integrals of the unknown densities of a finite number of materials along said line to an equal or larger number of energy-weighted integrals corresponding to different x-ray source energy spectra. ME-CT reconstructions may thus be decomposed as a two-step process: (i) reconstruct line integrals of the material densities from the available energy measurements; and (ii) reconstruct densities from their line integrals. Step (ii) is the standard linear x-ray CT problem whose invertibility is well-known, so this paper focuses on step (i). We show that ME-CT admits stable, global inversion provided that (a well-chosen linear transform of) the differential of the transform in step (i) satisfies appropriate orientation constraints that makes it a P-matrix. We introduce a notion of quantitative P-function that allows us to derive global stability results for ME-CT in the determined as well as over-determined (with more source energy spectra than the number of materials) cases. Numerical simulations based on standard material properties in imaging applications (of bone, water, contrast agents) and well accepted models of source energy spectra show that ME-CT is often (always in our simulations) either (i) non-globally injective because it is non-injective locally (differential not of full rank), or (ii) globally injective as soon as it is locally injective (differentials satisfy our proposed constraints).

preprint2019arXiv

Photo-acoustic tomography in a rotating setting

Photo-acoustic tomography is a coupled-physics (hybrid) medical imaging modality that aims to reconstruct optical parameters in biological tissues from ultrasound measurements. As propagating light gets partially absorbed, the resulting thermal expansion generates minute ultrasonic signals (the photo-acoustic effect) that are measured at the boundary of a domain of interest. Standard inversion procedures first reconstruct the source of radiation by an inverse ultrasound (boundary) problem and second describe the optical parameters from internal information obtained in the first step. This paper considers the rotating experimental setting. Light emission and ultrasound measurements are fixed on a rotating gantry, resulting in a rotation-dependent source of ultrasound. The two-step procedure we just mentioned does not apply. Instead, we propose an inversion that directly aims to reconstruct the optical parameters quantitatively. The mapping from the unknown (absorption and diffusion) coefficients to the ultrasound measurement via the unknown ultrasound source is modeled as a composition of a pseudo-differential operator and a Fourier integral operator. We show that for appropriate choices of optical illuminations, the above composition is an elliptic Fourier integral operator. Under the assumption that the coefficients are unknown on a sufficiently small domain, we derive from this a (global) injectivity result (measurements uniquely characterize our coefficients) combined with an optimal stability estimate. The latter is the same as that obtained in the standard (non-rotating experimental) setting.