Source author record

Guillaume Bal

Guillaume Bal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

60works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

60 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.

preprint2016arXiv

A dynamical polynomial chaos approach for long-time evolution of SPDEs

We propose a Dynamical generalized Polynomial Chaos (DgPC) method to solve time-dependent stochastic partial differential equations (SPDEs) with white noise forcing. The long-time simulation of SPDE solutions by Polynomial Chaos (PC) methods is notoriously difficult as the dimension of the stochastic variables increases linearly with time. Exploiting the markovian property of white noise, DgPC [1] implements a restart procedure that allows us to expand solutions at future times in terms of orthogonal polynomials of the measure describing the solution at a given time and the future white noise. The dimension of the representation is kept minimal by application of a Karhunen--Loeve (KL) expansion. Using frequent restarts and low degree polynomials on sparse multi-index sets, the method allows us to perform long time simulations, including the calculation of invariant measures for systems which possess one. We apply the method to the numerical simulation of stochastic Burgers and Navier--Stokes equations with white noise forcing. Our method also allows us to incorporate time-independent random coefficients such as a random viscosity. We propose several numerical simulations and show that the algorithm compares favorably with standard Monte Carlo methods.

preprint2016arXiv

Dynamical polynomial chaos expansions and long time evolution of differential equations with random forcing

Polynomial chaos expansions (PCE) allow us to propagate uncertainties in the coefficients of differential equations to the statistics of their solutions. Their main advantage is that they replace stochastic equations by systems of deterministic equations. Their main challenge is that the computational cost becomes prohibitive when the dimension of the parameters modeling the stochasticity is even moderately large. We propose a generalization of the PCE framework that allows us to keep this dimension as small as possible in favorable situations. For instance, in the setting of stochastic differential equations (SDEs) with Markov random forcing, we expect the future evolution to depend on the present solution and the future stochastic variables. We present a restart procedure that precisely allows PCE to depend only on that information. The computational difficulty then becomes the construction of orthogonal polynomials for dynamically evolving measures. We present theoretical results of convergence for our Dynamical generalized Polynomial Chaos (DgPC) method. Numerical simulations for linear and nonlinear SDEs show that it adequately captures the long-time behavior of their solutions as well as their invariant measures when the latter exist.

preprint2016arXiv

Photo-acoustic tomography in a rotating measurement setting

Photo-acoustic tomography (PAT) aims to leverage the photo-acoustic coupling between optical absorption of light sources and ultrasound (US) emission to obtain high contrast reconstructions of optical parameters with the high resolution of sonic waves. Quantitative PAT often involves a two-step procedure: first the map of sonic emission is reconstructed from US boundary measurements; and second optical properties of biological tissues are evaluated. We consider here a practical measurement setting in which such a separation does not apply. We assume that the optical source and an array of ultrasonic transducers are mounted on a rotating frame (in two or three dimensions) so that the light source rotates at the same time as the US measurements are acquired. As a consequence, we no longer have the option to reconstruct a map of sonic emission corresponding to a given optical illumination. We propose here a framework where the two steps are combined into one and an absorption map is directly reconstructed from the available ultrasound measurements.

preprint2015arXiv

Imaging of complex-valued tensors for two-dimensional Maxwell's equations

This paper concerns the imaging of a complex-valued anisotropic tensor γ = σ+ιωε from knowledge of several inter magnetic fields H where H satisfies the anisotropic Maxwell system on a bounded domain with prescribed boundary conditions. We show that γ can be uniquely reconstructed with a loss of two derivatives from errors in the acquisition H. A minimum number of five well-chosen functionals guaranties a local reconstruction of γ in dimension two. The explicit inversion procedure is presented in several numerical simulations, which demonstrate the influence of the choice boundary conditions on the stability of the reconstruction. This problem finds applications in the medical imaging modalities Current Density Imaging and Magnetic Resonance Electrical Impedance Tomography.

preprint2015arXiv

Reconstruction of a fully anisotropic elasticity tensor from knowledge of displacement fields

We present explicit reconstruction algorithms for fully anisotropic unknown elasticity tensors from knowledge of a finite number of internal displacement fields, with applications to transient elastography. Under certain rank-maximality assumptions satified by the strain fields, explicit algebraic reconstruction formulas are provided. A discussion ensues on how to fulfill these assumptions, describing the range of validity of the approach. We also show how the general method can be applied to more specific cases such as the transversely isotropic one.

preprint2015arXiv

Reconstruction of Cloud Geometry from High-Resolution Multi-Angle Images

We consider the reconstruction of the interface of compact, connected "clouds" from satellite or airborne light intensity measurements. In a two dimensional setting, the cloud is modeled by an interface, locally represented as a graph, and an outgoing radiation intensity that is consistent with a diffusion model for light propagation in the cloud. Light scattering inside the cloud and the optical internal parameters of the cloud are not modeled. The main objective is to understand what can or cannot be reconstructed in such a setting from intensity measurements in a finite (on the order of 10) number of directions along the path of a satellite. Numerical simulations illustrate the theoretical predictions.

preprint2015arXiv

Ultrasound modulated bioluminescence tomography and controllability of the radiative transport equation

We propose a method to reconstruct the density of an optical source in a highly scattering medium from ultrasound-modulated optical measurements. Our approach is based on the solution to a hybrid inverse source problem for the radiative transport equation (RTE). A controllability result for the RTE plays an essential role in the analysis.

preprint2014arXiv

Fluctuations of Parabolic Equations with Large Random Potentials

In this paper, we present a fluctuation analysis of a type of parabolic equations with large, highly oscillatory, random potentials around the homogenization limit. With a Feynman-Kac representation, the Kipnis-Varadhan's method, and a quantitative martingale central limit theorem, we derive the asymptotic distribution of the rescaled error between heterogeneous and homogenized solutions under different assumptions in dimension $d\geq 3$. The results depend highly on whether a stationary corrector exits.

preprint2014arXiv

Homogenization of Parabolic Equations with Large Time-dependent Random Potential

This paper concerns the homogenization problem of a parabolic equation with large, time-dependent, random potentials in high dimensions $d\geq 3$. Depending on the competition between temporal and spatial mixing of the randomness, the homogenization procedure turns to be different. We characterize the difference by proving the corresponding weak convergence of Brownian motion in random scenery. When the potential depends on the spatial variable macroscopically, we prove a convergence to SPDE.

preprint2014arXiv

Imaging of anisotropic conductivities from current densities in two dimensions

We consider the imaging of anisotropic conductivity tensors $γ=(γ_{ij})_{1\leq i,j\leq 2}$ from knowledge of several internal current densities $\mathcal{J}=γ\nabla u$ where $u$ satisfies a second order elliptic equation $\nabla\cdot(γ\nabla u)=0$ on a bounded domain $X\subset\mathbb{R}^2$ with prescribed boundary conditions on $\partial X$. We show that $γ$ can be uniquely reconstructed from four {\em well-chosen} functionals $\mathcal{J}$ and that noise in the data is differentiated once during the reconstruction. The inversion procedure is local in the sense that (most of) the tensor $γ(x)$ can be reconstructed from knowledge of the functionals $\mathcal{J}$ in the vicinity of $x$. We obtain the existence of an open set of boundary conditions on $\partial X$ that guaranty stable reconstructions by using the technique of complex geometric optics (CGO) solutions. The explicit inversion procedure is presented in several numerical simulations, which demonstrate the influence of the choice boundary conditions on the stability of the reconstruction. This problem finds applications in the medical imaging modality called Current Density Imaging or Magnetic Resonance Electrical Impedance Tomography.

preprint2014arXiv

Reconstruction of constitutive parameters in isotropic linear elasticity from noisy full-field measurements

Within the framework of linear elasticity we assume the availability of internal full-field measurements of the continuum deformations of a non-homogeneous isotropic solid. The aim is the quantitative reconstruction of the associated moduli. A simple gradient system for the sought constitutive parameters is derived algebraically from the momentum equation, whose coefficients are expressed in terms of the measured displacement fields and their spatial derivatives. Direct integration of this system is discussed to finally demonstrate the inexpediency of such an approach when dealing with noisy data. Upon using polluted measurements, an alternative variational formulation is deployed to invert for the physical parameters. Analysis of this latter inversion procedure provides existence and uniqueness results while the reconstruction stability with respect to the measurements is investigated. As the inversion procedure requires differentiating the measurements twice, a numerical differentiation scheme based on an ad hoc regularization then allows an optimally stable reconstruction of the sought moduli. Numerical results are included to illustrate and assess the performance of the overall approach.

preprint2014arXiv

Weak Convergence Approach for Parabolic Equations with Large, Highly Oscillatory, Random Potential

This paper concerns the macroscopic behavior of solutions to parabolic equations with large, highly oscillatory, random potential. When the correlation function of the random potential satisfies a specific integrability condition, we show that the random solution converges, as the correlation length of the medium tends to zero, to the deterministic solution of a homogenized equation in dimension $d\geq3$. Our derivation is based on a Feynman-Kac probabilistic representation and the Kipnis-Varadhan method applied to weak convergence of Brownian motions in random sceneries. For sufficiently mixing coefficients, we also provide an optimal rate of convergence to the homogenized limit using a quantitative martingale central limit theorem. As soon as the above integrability condition fails, the solution is expected to remain stochastic in the limit of a vanishing correlation length. For a large class of potentials given as functionals of Gaussian fields, we show the convergence of solutions to stochastic partial differential equations (SPDE) with multiplicative noise. The Feynman-Kac representation and the corresponding weak convergence of Brownian motions in random sceneries allows us to explain the transition from deterministic to stochastic limits as a function of the correlation function of the random potential.

preprint2013arXiv

Homogenization of elliptic equations with large random potential

We consider an elliptic equation with purely imaginary, highly heterogeneous, and large random potential with a sufficiently rapidly decaying correlation function. We show that its solution is well approximated by the solution to a homogeneous equation with a real-valued homogenized potential as the correlation length of the random medium $\varepsilon\rightarrow 0$ and estimate the size of the random fluctuations in the setting $d\geq3$.

preprint2013arXiv

Hybrid inverse problems and redundant systems of partial differential equations

Hybrid inverse problems are mathematical descriptions of coupled-physics (also called multi-waves) imaging modalities that aim to combine high resolution with high contrast. The solution of a high-resolution inverse problem, a first step that is not considered in this paper, provides internal information combining unknown parameters and solutions of differential equations. In several settings, the internal information and the differential equations may be described as a redundant system of nonlinear partial differential equations. We propose a framework to analyze the uniqueness and stability properties of such systems. We consider the case when the linearization of the redundant system is elliptic and with boundary conditions satisfying the Lopatinskii conditions. General theories of elliptic systems then allow us to construct a parametrix for such systems and derive optimal stability estimates. The injectivity of the nonlinear problem or its linearization is not guaranteed by the ellipticity condition. We revisit unique continuation principles, such as the Holmgren theorem and the uniqueness theorem of Calderón, in the context of redundant elliptic systems of equations. The theory is applied to the case of power density measurements, which are internal functionals of the form $γ|\nabla u|^2$ where $γ$ is an unknown parameter and $u$ is the solution to the elliptic equation $\nabla\cdotγ\nabla u=0$ on a bounded domain with appropriate boundary conditions.

preprint2013arXiv

Hybrid Inverse Problems for a System of Maxwell's Equations

This paper concerns the quantitative step of the medical imaging modality Thermo-acoustic Tomography (TAT). We model the radiation propagation by a system of Maxwell's equations. We show that the index of refraction of light and the absorption coefficient (conductivity) can be uniquely and stably reconstructed from a sufficiently large number of TAT measurements. Our method is based on verifying that the linearization of the inverse problem forms a redundant elliptic system of equations. We also observe that the reconstructions are qualitatively quite different from the setting where radiation is modeled by a scalar Helmholtz equation as in [10].

preprint2013arXiv

Inverse anisotropic conductivity from internal current densities

This paper concerns the reconstruction of an anisotropic conductivity tensor $γ$ from internal current densities of the form $J = γ\nabla u$, where $u$ solves a second-order elliptic equation $\nabla\cdot(γ\nabla u) = 0$ on a bounded domain $X$ with prescribed boundary conditions. A minimum number of such functionals equal to $n + 2$, where $n$ is the spatial dimension, is sufficient to guarantee a local reconstruction. We show that $γ$ can be uniquely reconstructed with a loss of one derivative compared to errors in the measurement of $J$. In the special case where $γ$ is scalar, it can be reconstructed with no loss of derivatives. We provide a precise statement of what components may be reconstructed with a loss of zero or one derivatives.

preprint2013arXiv

Linearized internal functionals for anisotropic conductivities

This paper concerns the reconstruction of an anisotropic conductivity tensor in an elliptic second-order equation from knowledge of the so-called power density functionals. This problem finds applications in several coupled-physics medical imaging modalities such as ultrasound modulated electrical impedance tomography and impedance-acoustic tomography. We consider the linearization of the nonlinear hybrid inverse problem. We find sufficient conditions for the linearized problem, a system of partial differential equations, to be elliptic and for the system to be injective. Such conditions are found to hold for a lesser number of measurements than those required in recently established explicit reconstruction procedures for the nonlinear problem.

preprint2013arXiv

Local Inversions in Ultrasound Modulated Optical Tomography

Ultrasound modulated optical tomography, also called acousto-optics tomography, is a hybrid imaging modality that aims to combine the high contrast of optical waves with the high resolution of ultrasound. We follow the model of the influence of ultrasound modulation on the light intensity measurements developed in [Bal Schotland PRL 2010]. We present sufficient conditions ensuring that the absorption and diffusion coefficients modeling light propagation can locally be uniquely and stably reconstructed from the corresponding available information. We present an iterative procedure to solve such a problem based on the analysis of linear elliptic systems of redundant partial differential equations.

preprint2013arXiv

Propagation of seismic waves through a spatio-temporally fluctuating medium: Homogenization

Measurements of seismic wave travel times at the photosphere of the Sun have enabled inferences of its interior structure and dynamics. In interpreting these measurements, the simplifying assumption that waves propagate through a temporally stationary medium is almost universally invoked. However, the Sun is in a constant state of evolution, on a broad range of spatio-temporal scales. At the zero wavelength limit, i.e., when the wavelength is much shorter than the scale over which the medium varies, the WKBJ (ray) approximation may be applied. Here, we address the other asymptotic end of the spectrum, the infinite wavelength limit, using the technique of homogenization. We apply homogenization to scenarios where waves are propagating through rapidly varying media (spatially and temporally), and derive effective models for the media. One consequence is that a scalar sound speed becomes a tensorial wavespeed in the effective model and anisotropies can be induced depending on the nature of the perturbation. The second term in this asymptotic two-scale expansion, the so-called corrector, contains contributions due to higher-order scattering, leading to the decoherence of the wavefield. This decoherence may be causally linked to the observed wave attenuation in the Sun. Although the examples we consider here consist of periodic arrays of perturbations to the background, homogenization may be extended to ergodic and stationary random media. This method may have broad implications for the manner in which we interpret seismic measurements in the Sun and for modeling the effects of granulation on the scattering of waves and distortion of normal-mode eigenfunctions.

preprint2013arXiv

Reconstruction of complex-valued tensors in the Maxwell system from knowledge of internal magnetic fields

This paper concerns the reconstruction of a complex-valued anisotropic tensor $γ=σ+ıω\varepsilon$ from knowledge of several internal magnetic fields $H$, where $H$ satisfies the anisotropic Maxwell system on a bounded domain with prescribed boundary conditions. We show that $γ$ can be uniquely reconstructed with a loss of two derivatives from errors in the acquisition of $H$. A minimum number of 6 such functionals is sufficient to obtain a local reconstruction of $γ$. In the special case where $γ$ is close to a scalar tensor, boundary conditions are chosen by means of complex geometric optics (CGO) solutions. For arbitrary symmetric tensors $γ$, a Runge approximation property is used to obtain partial results. This problem finds applications in the medical imaging modalities Current Density Imaging and Magnetic Resonance Electrical Impedance Tomography.

preprint2012arXiv

A corrector theory for diffusion-homogenization limits of linear transport equations

This paper concerns the diffusion-homogenization of transport equations when both the adimensionalized scale of the heterogeneities $α$ and the adimensionalized mean-free path $\eps$ converge to 0. When $α=\eps$, it is well known that the heterogeneous transport solution converges to a homogenized diffusion solution. We are interested here in the situation where $0<\eps\llα\ll1$ and in the respective rates of convergences to the homogenized limit and to the diffusive limit. Our main result is an approximation to the transport solution with an error term that is negligible compared to the maximum of $α$ and $\frac\epsα$. After establishing the diffusion-homogenization limit to the transport solution, we show that the corrector is dominated by an error to homogenization when $α^2\ll\eps$ and by an an error to diffusion when $\eps\llα^2$. Our regime of interest involves singular perturbations in the small parameter $η=\frac\epsα$. Disconnected local equilibria at $η=0$ need to be reconnected to provide a global equilibrium on the cell of periodicity when $η>0$. This reconnection between local and global equilibria is shown to hold when sufficient {\em no-drift} conditions are satisfied. The Hilbert expansion methodology followed in this paper builds on corrector theories for the result developed in \cite{NBAPuVo}.

preprint2012arXiv

Boundary control of elliptic solutions to enforce local constraints

We present a constructive method to devise boundary conditions for solutions of second-order elliptic equations so that these solutions satisfy specific qualitative properties such as: (i) the norm of the gradient of one solution is bounded from below by a positive constant in the vicinity of a finite number of prescribed points; and (ii) the determinant of gradients of $n$ solutions is bounded from below in the vicinity of a finite number of prescribed points. Such constructions find applications in recent hybrid medical imaging modalities. The methodology is based on starting from a controlled setting in which the constraints are satisfied and continuously modifying the coefficients in the second-order elliptic equation. The boundary condition is evolved by solving an ordinary differential equation (ODE) defined so that appropriate optimality conditions are satisfied. Unique continuations and standard regularity results for elliptic equations are used to show that the ODE admits a solution for sufficiently long times.

preprint2012arXiv

Cauchy problem for Ultrasound Modulated EIT

Ultrasound modulation of electrical or optical properties of materials offers the possibility to devise hybrid imaging techniques that combine the high electrical or optical contrast observed in many settings of interest with the high resolution of ultrasound. Mathematically, these modalities require that we reconstruct a diffusion coefficient $σ(x)$ for $x\in X$, a bounded domain in $\Rm^n$, from knowledge of $σ(x)|\nabla u|^2(x)$ for $x\in X$, where $u$ is the solution to the elliptic equation $-\nabla\cdotσ\nabla u=0$ in $X$ with $u=f$ on $\partial X$. This inverse problem may be recast as a nonlinear equation, which formally takes the form of a 0-Laplacian. Whereas $p-$Laplacians with $p>1$ are well-studied variational elliptic non-linear equations, $p=1$ is a limiting case with a convex but not strictly convex functional, and the case $p<1$ admits a variational formulation with a functional that is not convex. In this paper, we augment the equation for the 0-Laplacian with full Cauchy data at the domain's boundary, which results in a, formally overdetermined, nonlinear hyperbolic equation. The paper presents existence, uniqueness, and stability results for the Cauchy problem of the 0-Laplacian. In general, the diffusion coefficient $σ(x)$ can be stably reconstructed only on a subset of $X$ described as the domain of influence of the space-like part of the boundary $\partial X$ for an appropriate Lorentzian metric. Global reconstructions for specific geometries or based on the construction of appropriate complex geometric optics solutions are also analyzed.

preprint2012arXiv

Convergence to SPDE of the Schrodinger equation with large, random potential

We study the asymptotic behavior of solutions to the Schr{ö}dinger equation with large-amplitude, highly oscillatory, random potential. In dimension $d<\mathfrak{m}$, where $\mathfrak{m}$ is the order of the leading operator in the Schrödinger equation, we construct the heterogeneous solution by using a Duhamel expansion and prove that it converges in distribution, as the correlation length $\varepsilon$ goes to 0, to the solution of a stochastic differential equation, whose solution is represented as a sum of iterated Stratonovich integral, over the space $C([0,+\infty),\mathcal{S}')$. The uniqueness of the limiting solution in a dense space of $L^2(Ω\times\mathbb{R}^d)$ is shown by verifying the property of conservation of mass for the Schrödinger equation. In dimension $d>\mathfrak{m}$, the solution to the Schr{ö}dinger equation is shown to converge in $L^2(Ω\times\mathbb{R}^d)$ to a deterministic Schr{ö}dinger solution in \cite{ZB-12}.

preprint2012arXiv

Explicit Reconstructions in QPAT, QTAT, TE, and MRE

Photo-acoustic Tomography (PAT) and Thermo-acoustic Tomography (TAT) are medical imaging modalities that combine the high contrast of radiative properties of tissues with the high resolution of ultrasound. In both modalities, a first step concerns the reconstruction of the radiation-induced source of ultrasound. Transient Elastography (TE) and Magnetic Resonance Elastography (MRE) combine the high elastic contrast of tissues with the high resolution of ultrasound and magnetic resonance, respectively. In both modalities, a first step concerns the reconstruction of the elastic displacement. The result of this first step, which is not considered in this paper, is the availability of internal functionals of the unknown tissue properties. All imaging modalities are recast as the reconstruction of parameters in elliptic equations from knowledge of solutions to such equations. This paper provides a characterization of the parameters that may or may not be reconstructed from such internal functionals. We provide explicit reconstruction procedures and indicate how stable they are with respect to errors in the available measurements. The modalities PAT, TAT, TE, and MRE allow us to reconstruct high-contrast optical and elastic properties of tissues with the high resolution of ultrasound or magnetic resonance imaging. They provide a means to reconstruct second-order tensors modeling tissue anisotropy as well as complex-valued coefficients modeling absorbing and dissipative effects.

preprint2012arXiv

Homogenization of the Schroedinger equation with large, random potential

We study the behavior of solutions to a Schr{ö}dinger equation with large, rapidly oscillating, mean zero, random potential with Gaussian distribution. We show that in high dimension $d>\mathfrak{m}$, where $\mathfrak{m}$ is the order of the spatial pseudo-differential operator in the Schr{ö}dinger equation (with $\mathfrak{m}=2$ for the standard Laplace operator), the solution converges in the $L^2$ sense uniformly in time over finite intervals to the solution of a deterministic Schr{ö}dinger equation as the correlation length $\varepsilon$ tends to 0. This generalizes to long times the convergence results obtained for short times and for the heat equation. The result is based on a careful decomposition of multiple scattering contributions. In dimension $d<\mathfrak{m}$, the random solution converges to the solution of a stochastic partial differential equation.

preprint2012arXiv

Inverse anisotropic conductivity from power densities in dimension $n\ge 3$

We investigate the problem of reconstructing a fully anisotropic conductivity tensor $γ$ from internal functionals of the form $\nabla u\cdotγ\nabla u$ where $u$ solves $\nabla\cdot(γ\nabla u) = 0$ over a given bounded domain $X$ with prescribed Dirichlet boundary condition. This work motivated by hybrid medical imaging methods covers the case $n\ge 3$, following the previously published case $n=2$ \cite{Monard2011}. Under knowledge of enough such functionals, and writing $γ= τ\tilde γ$ ($\det \tildeγ= 1$) with $τ$ a positive scalar function, we show that all of $γ$ can be explicitely and locally reconstructed, with no loss of scales for $τ$ and loss of one derivative for the anisotropic structure $\tildeγ$. The reconstruction algorithms presented require rank maximality conditions that must be satisfied by the functionals or their corresponding solutions, and we discuss different possible ways of ensuring these conditions for $\C^{1,α}$-smooth tensors ($0<α<1$).

preprint2012arXiv

Inverse anisotropic diffusion from power density measurements in two dimensions

This paper concerns the reconstruction of an anisotropic diffusion tensor $γ=(γ_{ij})_{1\leq i,j\leq 2}$ from knowledge of internal functionals of the form $γ\nabla u_i\cdot\nabla u_j$ with $u_i$ for $1\leq i\leq I$ solutions of the elliptic equation $\nabla \cdot γ\nabla u_i=0$ on a two dimensional bounded domain with appropriate boundary conditions. We show that for I=4 and appropriately chosen boundary conditions, $γ$ may uniquely and stably be reconstructed from such internal functionals, which appear in coupled-physics inverse problems involving the ultrasound modulation of electrical or optical coefficients. Explicit reconstruction procedures for the diffusion tensor are presented and implemented numerically.

preprint2012arXiv

Inverse diffusion from knowledge of power densities

This paper concerns the reconstruction of a diffusion coefficient in an elliptic equation from knowledge of several power densities. The power density is the product of the diffusion coefficient with the square of the modulus of the gradient of the elliptic solution. The derivation of such internal functionals comes from perturbing the medium of interest by acoustic (plane) waves, which results in small changes in the diffusion coefficient. After appropriate asymptotic expansions and (Fourier) transformation, this allow us to construct the power density of the equation point-wise inside the domain. Such a setting finds applications in ultrasound modulated electrical impedance tomography and ultrasound modulated optical tomography. We show that the diffusion coefficient can be uniquely and stably reconstructed from knowledge of a sufficient large number of power densities. Explicit expressions for the reconstruction of the diffusion coefficient are also provided. Such results hold for a large class of boundary conditions for the elliptic equation in the two-dimensional setting. In three dimensions, the results are proved for a more restrictive class of boundary conditions constructed by means of complex geometrical optics solutions.

preprint2012arXiv

Inverse diffusion problems with redundant internal information

This paper concerns the reconstruction of a scalar diffusion coefficient $σ(x)$ from redundant functionals of the form $H_i(x)=σ^{2α}(x)|\nabla u_i|^2(x)$ where $α\in\Rm$ and $u_i$ is a solution of the elliptic problem $\nabla\cdot σ\nabla u_i=0$ for $1\leq i\leq I$. The case $α=\frac12$ is used to model measurements obtained from modulating a domain of interest by ultrasound and finds applications in ultrasound modulated electrical impedance tomography (UMEIT) as well as ultrasound modulated optical tomography (UMOT). The case $α=1$ finds applications in Magnetic Resonance Electrical Impedance Tomography (MREIT). We present two explicit reconstruction procedures of $σ$ for appropriate choices of $I$ and of traces of $u_i$ at the boundary of a domain of interest. The first procedure involves the solution of an over-determined system of ordinary differential equations and generalizes to the multi-dimensional case and to (almost) arbitrary values of $α$ the results obtained in two and three dimensions in \cite{CFGK} and \cite{BBMT}, respectively, in the case $α=\frac12$. The second procedure consists of solving a system of linear elliptic equations, which we can prove admits a unique solution in specific situations.

preprint2012arXiv

Radiative Transport Limit of Dirac Equations with Random Electromagnetic Field

This paper concerns the kinetic limit of the Dirac equation with random electromagnetic field. We give a detailed mathematical analysis of the radiative transport limit for the phase space energy density of solutions to the Dirac equation. Our derivation is based on a martingale method and a perturbed test function expansion. This requires the electromagnetic field to be a space-time random field. The main mathematical tool in the derivation of the kinetic limit is the matrix-valued Wigner transform of the vector-valued Dirac solution. The major novelty compared to the scalar (Schrödinger) case is the proof of convergence of cross-modes to 0 weakly in space and almost surely in probability. The propagating modes are shown to converge in an appropriate strong sense to their deterministic limit.

preprint2012arXiv

Reconstructions for some coupled-physics inverse problems

This letter announces and summarizes results obtained in arXiv:1111.5051 and considers several natural extensions. The aforementioned paper proposes a procedure to reconstruct coefficients in a second-order, scalar, elliptic equation from knowledge of a sufficiently large number of its solutions. We present this derivation and extend it to show which parameters may or may not be reconstructed for several hybrid (also called coupled physics) imaging modalities including photo-acoustic tomography, thermo-acoustic tomography, transient elastography, and magnetic resonance elastography. Stability estimates are also proposed.

preprint2012arXiv

The Levenberg-Marquardt Iteration for Numerical Inversion of the Power Density Operator

In this paper we develop a convergence analysis in an infinite dimensional setting of the Levenberg-Marquardt iteration for the solution of a hybrid conductivity imaging problem. The problem consists in determining the spatially varying conductivity $σ$ from a series of measurements of power densities for various voltage inductions. Although this problem has been very well studied in the literature, convergence and regularizing properties of iterative algorithms in an infinite dimensional setting are still rudimentary. We provide a partial result under the assumptions that the derivative of the operator, mapping conductivities to power densities, is injective and the data is noise-free. Moreover, we implemented the Levenberg-Marquardt algorithm and tested it on simulated data.

preprint2011arXiv

A Hybrid (Monte-Carlo/Deterministic) Approach for Multi-Dimensional Radiation Transport

A novel hybrid Monte Carlo transport scheme is demonstrated in a scene with solar illumination, scattering and absorbing 2D atmosphere, a textured reflecting mountain, and a small detector located in the sky (mounted on a satellite or a airplane). It uses a deterministic approximation of an adjoint transport solution to reduce variance, computed quickly by ignoring atmospheric interactions. This allows significant variance and computational cost reductions when the atmospheric scattering and absorption coefficient are small. When combined with an atmospheric photon-redirection scheme, significant variance reduction (equivalently acceleration) is achieved in the presence of atmospheric interactions.

preprint2011arXiv

Corrector theory for MsFEM and HMM in random media

We analyze the random fluctuations of several multi-scale algorithms such as the multi-scale finite element method (MsFEM) and the finite element heterogeneous multiscale method (HMM), that have been developed to solve partial differential equations with highly heterogeneous coefficients. Such multi-scale algorithms are often shown to correctly capture the homogenization limit when the highly oscillatory random medium is stationary and ergodic. This paper is concerned with the random fluctuations of the solution about the deterministic homogenization limit. We consider the simplified setting of the one dimensional elliptic equation, where the theory of random fluctuations is well understood. We develop a fluctuation theory for the multi-scale algorithms in the presence of random environments with short-range and long-range correlations. What we find is that the computationally more expensive method MsFEM captures the random fluctuations both for short-range and long-range oscillations in the medium. The less expensive method HMM correctly captures the fluctuations for long-range oscillations and strongly amplifies their size in media with short-range oscillations. We present a modified scheme with an intermediate computational cost that captures the random fluctuations in all cases.

preprint2011arXiv

Hybrid inverse problems and internal functionals

This paper reviews recent results on hybrid inverse problems, which are also called coupled-physics inverse problems of multi-wave inverse problems. Inverse problems tend to be most useful in, e.g., medical and geophysical imaging, when they combine high contrast with high resolution. In some settings, a single modality displays either high contrast or high resolution but not both. In favorable situations, physical effects couple one modality with high contrast with another modality with high resolution. The mathematical analysis of such couplings forms the class of hybrid inverse problems. Hybrid inverse problems typically involve two steps. In a first step, a well-posed problem involving the high-resolution low-contrast modality is solved from knowledge of boundary measurements. In a second step, a quantitative reconstruction of the parameters of interest is performed from knowledge of the point-wise, internal, functionals of the parameters reconstructed during the first step. This paper reviews mathematical techniques that have been developed in recent years to address the second step. Mathematically, many hybrid inverse problems find interpretations in terms of linear and nonlinear (systems of) equations. In the analysis of such equations, one often needs to verify that qualitative properties of solutions to elliptic linear equations are satisfied, for instance the absence of any critical points. This paper reviews several methods to prove that such qualitative properties hold, including the method based on the construction of complex geometric optics solutions.

preprint2011arXiv

Importance Sampling and Adjoint Hybrid Methods in Monte Carlo Transport with Reflecting Boundaries

Adjoint methods form a class of importance sampling methods that are used to accelerate Monte Carlo (MC) simulations of transport equations. Ideally, adjoint methods allow for zero-variance MC estimators provided that the solution to an adjoint transport equation is known. Hybrid methods aim at (i) approximately solving the adjoint transport equation with a deterministic method; and (ii) use the solution to construct an unbiased MC sampling algorithm with low variance. The problem with this approach is that both steps can be prohibitively expensive. In this paper, we simplify steps (i) and (ii) by calculating only parts of the adjoint solution. More specifically, in a geometry with limited volume scattering and complicated reflection at the boundary, we consider the situation where the adjoint solution "neglects" volume scattering, whereby significantly reducing the degrees of freedom in steps (i) and (ii). A main application for such a geometry is in remote sensing of the environment using physics-based signal models. Volume scattering is then incorporated using an analog sampling algorithm (or more precisely a simple modification of analog sampling called a heuristic sampling algorithm) in order to obtain unbiased estimators. In geometries with weak volume scattering (with a domain of interest of size comparable to the transport mean free path), we demonstrate numerically significant variance reductions and speed-ups (figures of merit).

preprint2011arXiv

Inverse transport with isotropic time-harmonic sources

This paper concerns the reconstruction of the scattering coefficient in a two-dimensional transport equation from angularly averaged measurements when the probing source is isotropic and time-harmonic. This is a practical setting in the medical imaging modality called Optical Tomography. As the modulation frequency of the source increases, we show that the reconstruction of the scattering coefficient improves. More precisely, as the frequency $ω$ increases, we show that all frequencies of the scattering coefficient lower than $b$ are reconstructed stably with an accuracy that improves as $ω$ increases and $b$ decreases. The proofs are based on an analysis of the single scattering singularities of the transport equation and on careful analyses of oscillatory integrals by stationary phase arguments.

preprint2011arXiv

Reconstruction of coefficients in scalar second-order elliptic equations from knowledge of their solutions

This paper concerns the reconstruction of possibly complex-valued coefficients in a second-order scalar elliptic equation posed on a bounded domain from knowledge of several solutions of that equation. We show that for a sufficiently large number of solutions and for an open set of corresponding boundary conditions, all coefficients can be uniquely and stably reconstructed up to a well characterized gauge transformation. We also show that in some specific situations, a minimum number of such available solutions equal to $I_n=\frac12n(n+3)$ is sufficient to uniquely and globally reconstruct the unknown coefficients. This theory finds applications in several coupled-physics medical imaging modalities including photo-acoustic tomography, transient elastography, and magnetic resonance elastography.

preprint2010arXiv

Large Deviation Theory for a Homogenized and "Corrected" Elliptic ODE

We study a one-dimensional elliptic problem with highly oscillatory random diffusion coefficient. We derive a homogenized solution and a so-called Gaussian corrector. We also prove a "pointwise" large deviation principle (LDP) for the full solution and approximate this LDP with a more tractable form. These results allow one to access the limits of Gaussian correctors. In general, the corrector does not capture the large deviation behavior. Applications to uncertainty quantification are considered.

preprint2010arXiv

Ray Transforms on a Conformal Class of Curves

We introduce a technique for recovering a sufficiently smooth function from its ray transform over a wide class of curves in a general region of Euclidean space. The method is based on a complexification of the underlying vector fields defining the initial transport and recasting the problem in terms of complex-analytic function theory. Explicit inversion formulae are then given in a unified form. The method is then used to give inversion formulae for the attenuated ray transform.

preprint2010arXiv

The Inversion of Ray Transforms on a Conformal Class of Curves

We introduce a technique for recovering a sufficiently smooth function from its ray transform over a wide class of curves in a general region of Euclidean space. The method is based on a complexification of the underlying vector fields defining the initial transport and recasting the problem in terms of complex-analytic function theory. Explicit inversion formulae are then given in a unified form. The method is then used to give inversion formulae for the attenuated ray transform.

preprint2009arXiv

Inverse Diffusion Theory of Photoacoustics

This paper analyzes the reconstruction of diffusion and absorption parameters in an elliptic equation from knowledge of internal data. In the application of photo-acoustics, the internal data are the amount of thermal energy deposited by high frequency radiation propagating inside a domain of interest. These data are obtained by solving an inverse wave equation, which is well-studied in the literature. We show that knowledge of two internal data based on well-chosen boundary conditions uniquely determines two constitutive parameters in diffusion and Schroedinger equations. Stability of the reconstruction is guaranteed under additional geometric constraints of strict convexity. No geometric constraints are necessary when $2n$ internal data for well-chosen boundary conditions are available, where $n$ is spatial dimension. The set of well-chosen boundary conditions is characterized in terms of appropriate complex geometrical optics (CGO) solutions.

preprint2009arXiv

Inverse Transport Theory of Photoacoustics

We consider the reconstruction of optical parameters in a domain of interest from photoacoustic data. Photoacoustic tomography (PAT) radiates high frequency electromagnetic waves into the domain and measures acoustic signals emitted by the resulting thermal expansion. Acoustic signals are then used to construct the deposited thermal energy map. The latter depends on the constitutive optical parameters in a nontrivial manner. In this paper, we develop and use an inverse transport theory with internal measurements to extract information on the optical coefficients from knowledge of the deposited thermal energy map. We consider the multi-measurement setting in which many electromagnetic radiation patterns are used to probe the domain of interest. By developing an expansion of the measurement operator into singular components, we show that the spatial variations of the intrinsic attenuation and the scattering coefficients may be reconstructed. We also reconstruct coefficients describing anisotropic scattering of photons, such as the anisotropy coefficient $g(x)$ in a Henyey-Greenstein phase function model. Finally, we derive stability estimates for the reconstructions.

preprint2009arXiv

Time-dependent angularly averaged inverse transport

This paper concerns the reconstruction of the absorption and scattering parameters in a time-dependent linear transport equation from knowledge of angularly averaged measurements performed at the boundary of a domain of interest. We show that the absorption coefficient and the spatial component of the scattering coefficient are uniquely determined by such measurements. We obtain stability results on the reconstruction of the absorption and scattering parameters with respect to the measured albedo operator. The stability results are obtained by a precise decomposition of the measurements into components with different singular behavior in the time domain.

preprint2008arXiv

Convergence to SPDEs in Stratonovich form

We consider the perturbation of parabolic operators of the form $\partial_t+P(x,D)$ by large-amplitude highly oscillatory spatially dependent potentials modeled as Gaussian random fields. The amplitude of the potential is chosen so that the solution to the random equation is affected by the randomness at the leading order. We show that, when the dimension is smaller than the order of the elliptic pseudo-differential operator $P(x,D)$, the perturbed parabolic equation admits a solution given by a Duhamel expansion. Moreover, as the correlation length of the potential vanishes, we show that the latter solution converges in distribution to the solution of a stochastic parabolic equation with a multiplicative term that should be interpreted in the Stratonovich sense. The theory of mild solutions for such stochastic partial differential equations is developed. The behavior described above should be contrasted to the case of dimensions that are larger than or equal to the order of the elliptic pseudo-differential operator $P(x,D)$. In the latter case, the solution to the random equation converges strongly to the solution of a homogenized (deterministic) parabolic equation as is shown in the companion paper [2]. The stochastic model is therefore valid only for sufficiently small space dimensions in this class of parabolic problems.