BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Gunther Uhlmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider an inverse problem arising in nonlinear ultrasound imaging. The propagation of ultrasound waves is modeled by a quasilinear wave equation. We make measurements at the boundary of the medium encoded in the Dirichlet-to-Neumann map, and we show that these measurements determine the nonlinearity.
We consider the ultrasound imaging problem governed by a nonlinear wave equation of Westervelt type with variable wave speed. We show that the coefficient of nonlinearity can be recovered uniquely from knowledge of the Dirichlet-to-Neumann map. Our proof is based on a second order linearization and the use of Gaussian beam solutions to reduce the problem to the inversion of a weighted geodesic ray transform. We propose an inversion algorithm and report the results of a numerical implementation to solve the nonlinear ultrasound imaging problem in a transmission setting in the frequency domain.
We consider an isotropic elastic medium occupying a bounded domain D whose density and Lamé parameters are piecewise smooth. In the elastic wave initial value inverse problem, we are given the solution operator for the elastic wave equation, but only outside the domain D and only for initial data supported outside D, and we study the recovery of the density and Lamé parameters. For known density, results have recently been obtained using the scattering control method to recover wave speeds. Here, we extend this result to include the recovery of the density in addition to the Lamé parameters under certain geometric conditions using techniques from microlocal analysis and a connection to local tensor tomography.
We revisit the problem of recovering wave speeds and density across a curved interface from reflected wave amplitudes. Such amplitudes have been exploited for decades in (exploration) seismology in this context. However, the analysis in seismology has been based on linearization and mostly flat interfaces. Here, we present a nonlinear analysis allowing curved interfaces, establish uniqueness and provide a reconstruction, while making the notion of amplitude precise through a procedure rooted in microlocal analysis.
We study an inverse problem for the fractional Schrödinger equation (FSE) with a local perturbation by a linear partial differential operator (PDO) of order smaller than the order of the fractional Laplacian. We show that one can uniquely recover the coefficients of the PDO from the Dirichlet-to-Neumann (DN) map associated to the perturbed FSE. This is proved for two classes of coefficients: coefficients which belong to certain spaces of Sobolev multipliers and coefficients which belong to fractional Sobolev spaces with bounded derivatives. Our study generalizes recent results for the zeroth and first order perturbations to higher order perturbations.
We consider an inverse boundary value problem for a semilinear wave equation on a time-dependent Lorentzian manifold with time-like boundary. The time-dependent coefficients of the nonlinear terms can be recovered in the interior from the knowledge of the Neumann-to-Dirichlet map. Either distorted plane waves or Gaussian beams can be used to derive uniqueness.
We consider an inverse boundary value problem for a nonlinear elastic wave equation which was studied in [de Hoop, Uhlmann, Wang. Math. Ann. (2019) doi:10.1007/s00208-018-01796-y]. We show that all the parameters appearing in the equation can be uniquely determined from boundary measurements under certain geometric assumptions. The proof is based on second order linearization and Gaussian beams.
In this work we study the optimality of increasing stability of the inverse boundary value problem (IBVP) for Schrödinger equation. The rigorous justification of increasing stability for the IBVP for Schrödinger equation were established by Isakov \cite{Isa11} and by Isakov, Nagayasu, Uhlmann, Wang of the paper \cite{INUW14}. In \cite{Isa11}, \cite{INUW14}, the authors showed that the stability of this IBVP increases as the frequency increases in the sense that the stability estimate changes from a logarithmic type to a Hölder type. In this work, we prove that the instability changes from an exponential type to a Hölder type when the frequency increases. This result verifies that results in \cite{Isa11}, \cite{INUW14} are optimal.
We consider the mixed ray transform of tensor fields on a three-dimensional compact simple Riemannian manifold with boundary. We prove the injectivity of the transform, up to natural obstructions, and establish stability estimates for the normal operator on generic three dimensional simple manifold in the case of 1+1 and 2+2 tensors fields. We show how the anisotropic perturbations of averaged isotopic travel-times of qS-polarized elastic waves provide partial information about the mixed ray transform of 2+2 tensors fields. If in addition we include the measurement of the shear wave amplitude, the complete mixed ray transform can be recovered. We also show how one can obtain the mixed ray transform from an anisotropic perturbation of the Dirichlet-to-Neumann map of an isotropic elastic wave equation on a smooth and bounded domain in three dimensional Euclidean space.
We consider inverse boundary value problems for general real principal type differential operators. The first results state that the Cauchy data set uniquely determines the scattering relation of the operator and bicharacteristic ray transforms of lower order coefficients. We also give two different boundary determination methods for general operators, and prove global uniqueness results for determining coefficients in nonlinear real principal type equations. The article presents a unified approach for treating inverse boundary problems for transport and wave equations, and highlights the role of propagation of singularities in the solution of related inverse problems.
On a time-oriented Lorentzian manifold $(M,g)$ with non-empty boundary satisfying a convexity assumption, we show that the topological, differentiable, and conformal structure of suitable subsets $S\subset M$ of sources is uniquely determined by measurements of the intersection of future light cones from points in $S$ with a fixed open subset of the boundary of $M$; here, light rays are reflected at $\partial M$ according to Snell's law. Our proof is constructive, and allows for interior conjugate points as well as multiply reflected and self-intersecting light cones.
We consider an inverse problem for the Boltzmann equation with nonlinear collision operator in dimensions $n\geq 2$. We show that the kinetic collision kernel can be uniquely determined from the incoming-to-outgoing mappings on the boundary of the domain provided that the kernel satisfies a monotonicity condition. Furthermore, a reconstruction formula is also derived. The key methodology is based on the higher-order linearization scheme to reduce a nonlinear equation into simpler linear equations by introducing multiple small parameters into the original equation.
In this paper, we consider the boundary rigidity problem on a cylindrical domain in $\mathbb R^{1+n}$, $n\geq 2$, equipped with a stationary (time-invariant) Lorentzian metric. We show that the time separation function between pairs of points on the boundary of the cylindrical domain determines the stationary spacetime, up to some time-invariant diffeomorphism, assuming that the metric satisfies some a-priori conditions.
We show global uniqueness in the fractional Calderón problem with a single measurement and with data on arbitrary, possibly disjoint subsets of the exterior. The previous work \cite{GhoshSaloUhlmann} considered the case of infinitely many measurements. The method is again based on the strong uniqueness properties for the fractional equation, this time combined with a unique continuation principle from sets of measure zero. We also give a constructive procedure for determining an unknown potential from a single exterior measurement, based on constructive versions of the unique continuation result that involve different regularization schemes.
We study the weighted light ray transform $L$ of integrating functions on a Lorentzian manifold over lightlike geodesics. We analyze $L$ as a Fourier Integral Operator and show that if there are no conjugate points, one can recover the spacelike singularities of a function $f$ from its the weighted light ray transform $Lf$ by a suitable filtered back-projection.
Consider a Riemannian manifold in dimension $n\geq 3$ with strictly convex boundary. We prove the local invertibility, up to potential fields, of the geodesic ray transform on tensor fields of rank four near a boundary point. This problem is closely related with elastic \textit{qP}-wave tomography. Under the condition that the manifold can be foliated with a continuous family of strictly convex hypersurfaces, the local invertibility implies a global result. One can straightforwardedly adapt the proof to show similar results for tensor fields of arbitrary rank.
We show global uniqueness in an inverse problem for the fractional Schrödinger equation: an unknown potential in a bounded domain is uniquely determined by exterior measurements of solutions. We also show global uniqueness in the partial data problem where the measurements are taken in arbitrary open, possibly disjoint, subsets of the exterior. The results apply in any dimension $\geq 2$ and are based on a strong approximation property of the fractional equation that extends earlier work. This special feature of the nonlocal equation renders the analysis of related inverse problems radically different from the traditional Calderón problem.
We consider inverse problems in space-time $(M, g)$, a $4$-dimensional Lorentzian manifold. For semilinear wave equations $\square_g u + H(x, u) = f$, where $\square_g$ denotes the usual Laplace-Beltrami operator, we prove that the source-to-solution map $L: f \rightarrow u|_V$, where $V$ is a neighborhood of a time-like geodesic $μ$, determines the topological, differentiable structure and the conformal class of the metric of the space-time in the maximal set where waves can propagate from $μ$ and return back. Moreover, on a given space-time $(M, g)$, the source-to-solution map determines some coefficients of the Taylor expansion of $H$ in $u$.
We present two range characterizations for the attenuated geodesic X-ray transform defined on pairs of functions and one-forms on simple surfaces. Such characterizations are based on first isolating the range over sums of functions and one-forms, then separating each sub-range in two ways, first by implicit conditions, second by deriving new inversion formulas for sums of functions and one-forms.
We survey some results on travel time tomography. The question is whether we can determine the anisotropic index of refraction of a medium by measuring the travel times of waves going through the medium. This can be recast as geometry problems, the boundary rigidity problem and the lens rigidity problem. The boundary rigidity problem is whether we can determine a Riemannian metric of a compact Riemannian manifold with boundary by measuring the distance function between boundary points. The lens rigidity problem problem is to determine a Riemannian metric of a Riemannian manifold with boundary by measuring for every point and direction of entrance of a geodesic the point of exit and direction of exit and its length. The linearization of these two problems is tensor tomography. The question is whether one can determine a symmetric two-tensor from its integrals along geodesics. We emphasize recent results on boundary and lens rigidity and in tensor tomography in the partial data case.
It is well-known that reconstruction algorithms in quantitative susceptibility mapping often contain streaking artifacts. In [1], the cause of the artifacts is identified as propagation of singularities. In this work, we analyze such singularities carefully and propose some strategies to reduce the artifacts.
Motivated by non-destructive testing of optical fiber, we consider the problem of determining the index of refraction of a two-dimensional medium from magnitude of the total field resulting from known incident plane waves at a fixed frequency. The finiteness of the fiber cross section allows us to assume that the index of refraction is variable over a compact set. The measurement set is a closed curve, e.g., a circle of fixed radius, containing the fiber cross section. We use complex geometric optics (CGO) solutions to show that the inverse problem for the unknown index of refraction is unique in $C^{k}$ for large enough $k$. The key tool in this work is the multi-term expansion of the CGO solutions which we rely on to address the uniqueness question.
We show that the spectrum of a Schrödinger operator on $\mathbb{R}^n$, $n\ge 3$, with a periodic smooth Riemannian metric, whose conformal multiple has a product structure with one Euclidean direction, and with a periodic electric potential in $L^{n/2}_{\text{loc}}(\mathbb{R}^n)$, is purely absolutely continuous. Previously known results in the case of a general metric are obtained in [12], see also [8], under the assumption that the metric, as well as the potential, are reflection symmetric.
We study the boundary rigidity problem with partial data consisting of determining locally the Riemannian metric of a Riemannian manifold with boundary from the distance function measured at pairs of points near a fixed point on the boundary. We show that one can recover uniquely and in a stable way a conformal factor near a strictly convex point where we have the information. In particular, this implies that we can determine locally the isotropic sound speed of a medium by measuring the travel times of waves joining points close to a convex point on the boundary. The local results lead to a global lens rigidity uniqueness and stability result assuming that the manifold is foliated by strictly convex hypersurfaces.
This paper concerns thermoacoustic tomography and photoacoustic tomography, two couple-physics imaging modalities that attempt to combine the high resolution of ultrasound and the high contrast capabilities of electromagnetic waves. We give sufficient conditions to recover both the sound speed of the medium being probed and the source.
In this paper, we consider two regularized transformation-optics cloaking schemes for electromagnetic (EM) waves. Both schemes are based on the blowup construction with the generating sets being, respectively, a generic curve and a planar subset. We derive sharp asymptotic estimates in assessing the cloaking performances of the two constructions in terms of the regularization parameters and the geometries of the cloaking devices. The first construction yields an approximate full-cloak, whereas the second construction yields an approximate partial-cloak. Moreover, by incorporating properly chosen conducting layers, both cloaking constructions are capable of nearly cloaking arbitrary EM contents. This work complements the existing results in [5--7] on approximate EM cloaks with the generating set being a singular point, and it also extends [9,25] on regularized full and partial cloaks for acoustic waves governed by the Helmholtz system to the more challenging EM case governed by the full Maxwell system.
In the recent articles \cite{PSU1,PSU3}, a number of tensor tomography results were proved on two-dimensional manifolds. The purpose of this paper is to extend some of these methods to manifolds of any dimension. A central concept is the surjectivity of the adjoint of the geodesic ray transform, or equivalently the existence of certain distributions that are invariant under geodesic flow. We prove that on any Anosov manifold, one can find invariant distributions with controlled first Fourier coefficients. The proof is based on subelliptic type estimates and a Pestov identity. We present an alternative construction valid on manifolds with nonpositive curvature, based on the fact that a natural Beurling transform on such manifolds turns out to be essentially a contraction. Finally, we obtain uniqueness results in tensor tomography both on simple and Anosov manifolds that improve earlier results by assuming a condition on the terminator value for a modified Jacobi equation.
We show that the knowledge of the Dirichlet-to-Neumann map on the boundary of a bounded open set in $R^n$ for the perturbed polyharmonic operator $(-Δ)^m +q$ with $q\in L^{n/2m}$, $n>2m$, determines the potential $q$ in the set uniquely. In the course of the proof, we construct a special Green function for the polyharmonic operator and establish its mapping properties in suitable weighted $L^2$ and $L^p$ spaces. The $L^p$ estimates for the special Green function are derived from $L^p$ Carleman estimates with linear weights for the polyharmonic operator.
We consider the linearization of the Dirichlet-to-Neumann (DN) map as a function of the potential. We show that it is injective at a real analytic potential for measurements made at an open subset of analyticity of the boundary. More generally, we relate the analyticity up to the boundary of the variations of the potential to the analyticity of the symbols of the corresponding variations of the DN-map.
The aim of this work is to derive sharp quantitative estimates of the qualitative convergence results developed in [28] for regularized full- and partial-cloaks via the transformation-optics approach. Let $Γ_0$ be a compact set in $\mathbb{R}^3$ and $Γ_δ$ be a $δ$-neighborhood of $Γ_0$ for $δ\in\mathbb{R}_+$. $Γ_δ$ represents the virtual domain used for the blow-up construction. By incorporating suitably designed lossy layers, it is shown that if the generating set $Γ_0$ is a generic curve, then one would have an approximate full-cloak within $δ^2$ to the perfect full-cloak; whereas if $Γ_0$ is the closure of an open subset on a flat surface, then one would have an approximate partial-cloak within $δ$ to its perfect counterpart. The estimates derived are independent of the contents being cloaked; that is, the cloaking devices are capable of nearly cloaking an arbitrary content. Furthermore, as a significant byproduct, our argument allows the relaxation of the convexity requirement on $Γ_0$ in [28], which is critical for the Mosco convergence argument therein.
We consider how microlocal methods developed for tomographic problems can be used to detect singularities of the Lorentzian metric of the Universe using measurements of the Cosmic Microwave Background radiation. The physical model we study is mathematically rigorous but highly idealized.
We study the microlocal properties of the geodesic X-ray transform $\mathcal{X}$ on a manifold with boundary allowing the presence of conjugate points. Assuming that there are no self-intersecting geodesics and all conjugate pairs are nonsingular we show that the normal operator $\mathcal{N} = \mathcal{X}^t \circ \mathcal{X}$ can be decomposed as the sum of a pseudodifferential operator of order $-1$ and a sum of Fourier integral operators. We also apply this decomposition to prove inversion of $\mathcal{X}$ is only mildly ill-posed in dimension three or higher.
We show that given two hyperbolic Dirichlet to Neumann maps associated to two Riemannian metrics of a Riemannian manifold with boundary which coincide near the boundary are close then the lens data of the two metrics is the same. As a consequence, we prove uniqueness of recovery a conformal factor (sound speed) locally under some conditions on the latter.
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.
We extend a global uniqueness result for the Calderón problem with partial data, due to Kenig-Sjöstrand-Uhlmann, to the case of less regular conductivities. Specifically, we show that in dimensions $n\ge 3$, the knowledge of the Diricihlet-to-Neumann map, measured on possibly very small subsets of the boundary, determines uniquely a conductivity having essentially $3/2$ derivatives in an $L^2$ sense.
We consider integral geometry inverse problems for unitary connections and skew-Hermitian Higgs fields on manifolds with negative sectional curvature. The results apply to manifolds in any dimension, with or without boundary, and also in the presence of trapped geodesics. In the boundary case, we show injectivity of the attenuated ray transform on tensor fields with values in a Hermitian bundle (i.e. vector valued case). We also show that a connection and Higgs field on a Hermitian bundle are determined up to gauge by the knowledge of the parallel transport between boundary points along all possible geodesics. The main tools are an energy identity, the Pestov identity with a unitary connection, which is presented in a general form, and a precise analysis of the singularities of solutions of transport equations when there are trapped geodesics. In the case of closed manifolds, we obtain similar results modulo the obstruction given by twisted conformal Killing tensors, and we also study this obstruction.
Given a bounded domain $M$ in $\mathbb{R}^n$ with a conformally Euclidean metric $g=ρ\,dx^2$, in this paper we consider the inverse problem of recovering a semigeodesic neighborhood of a domain $Γ\subset \partial M$ and the conformal factor $ρ$ in the neighborhood from the travel time data (defined below) and the Cartesian coordinates of $Γ$. We develop an explicit reconstruction procedure for this problem. The key ingredient is the relation between the reconstruction and a Cauchy problem of the conformal Killing equation.
In this paper, we prove in two dimensions global identifiability of the viscosity in an incompressible fluid by making boundary measurements. The main contribution of this work is to use more natural boundary measurements, the Cauchy forces, than the Dirichlet-to-Neumann map previously considered in \cite{IY} to prove the uniqueness of the viscosity for the Stokes equations and for the Navier-Stokes equations.
We consider inverse problems for the coupled Einstein equations and the matter field equations on a 4-dimensional globally hyperbolic Lorentzian manifold $(M,g)$. We give a positive answer to the question: Do the active measurements, done in a neighborhood $U\subset M$ of a freely falling observed $μ=μ([s_-,s_+])$, determine the conformal structure of the spacetime in the minimal causal diamond-type set $V_g=J_g^+(μ(s_-))\cap J_g^-(μ(s_+))\subset M$ containing $μ$? More precisely, we consider the Einstein equations coupled with the scalar field equations and study the system $Ein(g)=T$, $T=T(g,ϕ)+F_1$, and $\square_gϕ-\mathcal V^\prime(ϕ)=F_2$, where the sources $F=(F_1,F_2)$ correspond to perturbations of the physical fields which we control. The sources $F$ need to be such that the fields $(g,ϕ,F)$ are solutions of this system and satisfy the conservation law $\nabla_jT^{jk}=0$. Let $(\hat g,\hat ϕ)$ be the background fields corresponding to the vanishing source $F$. We prove that the observation of the solutions $(g,ϕ)$ in the set $U$ corresponding to sufficiently small sources $F$ supported in $U$ determine $V_{\hat g}$ as a differentiable manifold and the conformal structure of the metric $\hat g$ in the domain $V_{\hat g}$. The methods developed here have potential to be applied to a large class of inverse problems for non-linear hyperbolic equations encountered e.g. in various practical imaging problems.
We prove the local invertibility, up to potential fields, and stability of the geodesic X-ray transform on tensor fields of order 1 and 2 near a strictly convex boundary point, on manifolds with boundary of dimension n>=3. We also present an inversion formula. Under the condition that the manifold can be foliated with a continuous family of strictly convex surfaces, we prove a global result which also implies a lens rigidity result near such a metric. The class of manifolds satisfying the foliation condition includes manifolds with no focal points, and does not exclude existence of conjugate points.
We study the Einstein equations coupled with the scalar field equations, $\hbox{Ein}(g)=T$, $T=T(g,ϕ)+F^1$, and $\square_gϕ^\ell-m^2ϕ^\ell= F^2$, where the sources $F=(F^1, F^2)$ correspond to perturbations of the physical fields which we control. Here $ϕ=(ϕ^\ell)_{\ell=1}^L$ and $(M,g)$ is a 4-dimensional globally hyperbolic Lorentzian manifold. The sources $F$ need to be such that the fields $(g,ϕ,F)$ satisfy the conservation law $\hbox{div}_g(T)=0$. If $(g_ε,ϕ_ε)$ solves the above equations, $\dot g=\partial_εg_ε|_{ε=0}$, $\dotϕ=ϕ_ε|_{ε=0}$, and $f=(f^1,f^2)= \partial_εF_ε|_{ε=0}$ solve the linearized Einstein equations and the linearized conservation law $$ \frac 12 \hat g^{pk}\hat \nabla_p f^1_{kj}+ \sum_{\ell=1}^L f^2_\ell \, \partial_j\hatϕ_\ell=0, $$ where $\hat g= g_ε|_{ε=0}$ and $\hat ϕ= ϕ_ε|_{ε=0}$. Then $(\hat g,\hat ϕ)$ and $f$ have the linearization stability property. Here ask the converse: If $\dot g$, $\dot ϕ$, and $f$ solve the linearized Einstein equations and the linearized conservation law, are there $F_ε=(F^1_ε,F^2_ε)$ and $(g_ε,ϕ_ε)$ depending on $ε\in [0,ε_0)$, $ε_0>0$, such that $(g_ε,ϕ_ε)$ solves the Einstein-scalar field equations and the conservation law. When $\hat g$ and $\hat ϕ$ vary enough and $L\geq 5$, we prove a microlocal version of this: When $Y\subset M$ is a 2-surface and $(y,η)\in N^*Y$, there is $f$ that is a conormal distibutions wrt. the surface $Y$ with a given principal symbol at $(y,η)$ such that $(\hat g,\hat ϕ)$ and $f$ have the linearization stability property.
We consider transformation-optics based cloaking in acoustic and electromagnetic scattering. The blueprints for an ideal cloak use singular acoustic and electromagnetic materials, posing server difficulties to both theoretical analysis and practical fabrication. In order to avoid the singular structures, various regularized approximate cloaking schemes have been developed. We survey these developments in this paper. We also propose some challenging issues for further investigation.
This article considers inverse problems on closed Riemannian surfaces whose geodesic flow is Anosov. We prove spectral rigidity for any Anosov surface and injectivity of the geodesic ray transform on solenoidal 2-tensors. We also establish surjectivity results for the adjoint of the geodesic ray transform on solenoidal tensors. The surjectivity results are of independent interest and imply the existence of many geometric invariant distributions on the unit sphere bundle. In particular, we show that on any Anosov surface $(M,g)$, given a smooth function $f$ on $M$ there is a distribution in the Sobolev space $H^{-1}(SM)$ that is invariant under the geodesic flow and whose projection to $M$ is the given function $f$.
We study the geodesic X-ray transform $X$ on compact Riemannian surfaces with conjugate points. Regardless of the type of the conjugate points, we show that we cannot recover the singularities and therefore, this transform is always unstable (ill-posed). We describe the microlocal kernel of $X$ and relate it to the conjugate locus. We present numerical examples illustrating the cancellation of singularities. We also show that the attenuated X-ray transform is well posed if the attenuation is positive and there are no more than two conjugate points along each geodesic; but still ill-posed, if there are three or more conjugate points. Those results follow from our analysis of the weighted X-ray transform.
We consider the inverse problem of recovering a potential by measuring the response at a point to a source located at the same point and then varying the point on the surface of a sphere. This is a similar to the inverse back-scattering problem. We show that if the angular derivatives of the difference of two potentials having the same data is controlled by the L^2 norm of the difference of the potentials they must be equal. In particular this shows injectivity of the inverse problem for radial potentials.
We consider the problem of recovering a smooth, compactly supported potential on R^3 from its backscattering data. We show that if two such potentials have the same backscattering data and the difference of the two potentials has controlled angular derivatives then the two potentials are identical. In particular, if two potentials differ by a finite linear combination of spherical harmonics with radial coefficinets and have the same backscattering data then the two potentials are identical.
We consider inverse problems for the Einstein equation with a time-depending metric on a 4-dimensional globally hyperbolic Lorentzian manifold $(M,g)$. We formulate the concept of active measurements for relativistic models. We do this by coupling the Einstein equation with equations for scalar fields and study the system $Ein(g)=T$, $T=T(g,ϕ)+F_1$, and $\square_g ϕ=F_2+S(g,ϕ,F_1,F_2)$. Here $F=(F_1,F_2)$ correspond to the perturbations of the physical fields which we control and $S$ is a secondary source corresponding to the adaptation of the system to the perturbation so that the conservation law $div_g(T)=0$ will be satisfied. The inverse problem we study is the question, do the observation of the solutions $(g,ϕ)$ in an open subset $U\subset M$ of the space-time corresponding to sources $F$ supported in $U$ determine the properties of the metric in a larger domain $W\subset M$ containing $U$. To study this problem we define the concept of light observation sets and show that these sets determine the conformal class of the metric. This corresponds to passive observations from a distant area of the space which is filled by light sources (e.g. we see light from stars varying in time). One can apply the obtained result to solve inverse problems encountered in general relativity and in various practical imaging problems.
We develop an algorithm for the computation of general Fourier integral operators associated with canonical graphs. The algorithm is based on dyadic parabolic decomposition using wave packets and enables the discrete approximate evaluation of the action of such operators on data in the presence of caustics. The procedure consists in the construction of a universal operator representation through the introduction of locally singularity-resolving diffeomorphisms, enabling the application of wave packet driven computation, and in the construction of the associated pseudo-differential joint-partition of unity on the canonical graphs. We apply the method to a parametrix of the wave equation in the vicinity of a cusp singularity.
We prove uniform $L^p$ estimates for resolvents of higher order elliptic self-adjoint differential operators on compact manifolds without boundary, generalizing a corresponding resul of [3] in the case of Laplace-- Beltrami operators on Riemannian manifolds. In doing so, we follow the methods, developed in [1] very closely. We also show that spectral regions in our $L^p$ resolvent estimates are optimal.
We describe the range of the attenuated ray transform of a unitary connection on a simple surface acting on functions and 1-forms. We use this to determine the range of the ray transform acting on symmetric tensor fields.
We consider the mathematical model of photoacoustic and thermoacoustic tomography in media with a variable sound speed. When the sound speed is known, the explicit reconstruction formula by P. Stefanov and G. Uhlmann (Inverse Problems, 25(7):075011, 16, 2009) can be used. We study how a modelling error in the sound speed affects the reconstruction formula and quantify the effect in terms of a stability estimate.
We develop a very general theory on the regularized approximate invisibility cloaking for the wave scattering governed by the Helmholtz equation in any space dimensions via the approach of transformation optics.
We survey recent progress in the problem of recovering a tensor field from its integrals along geodesics. We also propose several open problems.
We show that the knowledge of the set of the Cauchy data on the boundary of a $C^1$ bounded open set in $\R^n$, $n\ge 3$, for the Schrödinger operator with continuous magnetic and bounded electric potentials determines the magnetic field and electric potential inside the set uniquely. The proof is based on a Carleman estimate for the magnetic Schrödinger operator with a gain of two derivatives.
In this paper we show that for metrics with conormal singularities that correspond to class C^{1,α} with α>0, the reflected wave is more regular than the incident wave in a Sobolev sense. This is helpful in the analysis of the multiple scattering series since higher order terms can be effectively `peeled off'.
In this paper we consider the linearized problem of recovering both the sound speed and the thermal absorption arising in thermoacoustic and photoacoustic tomography. We show that the problem is unstable in any scale of Sobolev spaces.
In the plane, we study the transform $R_γf$ of integrating a unknown function $f$ over circles centered at a given curve $γ$. This is a simplified model of SAR, when the radar is not directed but has other applications, like thermoacoustic tomography, for example. We study the problem of recovering the wave front set $\WF(f)$. If the visible singularities of $f$ hit $γ$ once, we show that the "artifacts" cannot be resolved. If $γ$ is a closed curve, we show that this is still true. On the other hand, if $f$ is known a priori to have singularities in a compact set, then we show that one can recover $\WF(f)$, and moreover, this can be done in a simple explicit way, using backpropagation for the wave equation.
For the isotropic Lamé system, we prove in dimensions three or larger that both Lamé coefficients are uniquely recovered from partial Cauchy data on an arbitrary open subset of the boundary provided that the coefficient $μ$ is a priori close to a constant.
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.
We show that on simple surfaces the geodesic ray transform acting on solenoidal symmetric tensor fields of arbitrary order is injective. This solves a long standing inverse problem in the two-dimensional case.
We show that for a simple surface with boundary the attenuated ray transform in the presence of a unitary connection and a skew-Hermitian Higgs field is injective modulo the natural obstruction for functions and vector fields. We also show that the connection and the Higgs field are uniquely determined by the scattering relation modulo a gauge transformation. The proofs involve a Pestov type energy identity for connections together with holomorphic gauge transformations which arrange the curvature of the connection to have definite sign.
Under a convexity assumption on the boundary we solve a local inverse problem, namely we show that the geodesic X-ray transform can be inverted locally in a stable manner; one even has a reconstruction formula. We also show that under an assumption on the existence of a global foliation by strictly convex hypersurfaces the geodesic X-ray transform is globally injective. In addition we prove stability estimates and propose a layer stripping type algorithm for reconstruction.
We show that the knowledge of the set of the Cauchy data on the boundary of a bounded open set in $\R^n$, $n\ge 3$, for the magnetic Schrödinger operator with $L^\infty$ magnetic and electric potentials determines the magnetic field and electric potential inside the set uniquely. The proof is based on a Carleman estimate for the magnetic Schrödinger operator with a gain of two derivatives.
We present a new algorithm for reconstructing an unknown source in Thermoacoustic and Photoacoustic Tomography based on the recent advances in understanding the theoretical nature of the problem. We work with variable sound speeds that might be also discontinuous across some surface. The latter problem arises in brain imaging. The new algorithm is based on an explicit formula in the form of a Neumann series. We present numerical examples with non-trapping, trapping and piecewise smooth speeds, as well as examples with data on a part of the boundary. These numerical examples demonstrate the robust performance of the new algorithm.
We consider an operator $Δ^2 + A(x)\cdot D+q(x)$ with the Navier boundary conditions on a bounded domain in $R^n$, $n\ge 3$. We show that a first order perturbation $A(x)\cdot D+q$ can be determined uniquely by measuring the Dirichlet--to--Neumann map on possibly very small subsets of the boundary of the domain. Notice that the corresponding result does not hold in general for a first order perturbation of the Laplacian.
In this work we study the inverse boundary value problem of determining the refractive index in the acoustic equation. It is known that this inverse problem is ill-posed. Nonetheless, we show that the ill-posedness decreases when we increase the frequency and the stability estimate changes from logarithmic type for low frequencies to a Lipschitz estimate for large frequencies.
We prove identification of coefficients up to gauge by Cauchy data at the boundary for elliptic systems on oriented compact surfaces with boundary or domains of $\mathbb{C}$. In the geometric setting, we fix a Riemann surface with boundary, and consider both a Dirac-type operator plus potential acting on sections of a Clifford bundle and a connection Laplacian plus potential (i.e. Schrödinger Laplacian with external Yang-Mills field) acting on sections of a Hermitian bundle. In either case we show that the Cauchy data determines both the connection and the potential up to a natural gauge transformation: conjugation by an endomorphism of the bundle which is the identity at the boundary. For domains of $\mathbb{C}$, we recover zeroth order terms up to gauge from Cauchy data at the boundary in first order elliptic systems.
We show that a first order perturbation $A(x)\cdot D+q(x)$ of the polyharmonic operator $(-Δ)^m$, $m\ge 2$, can be determined uniquely from the set of the Cauchy data for the perturbed polyharmonic operator on a bounded domain in $R^n$, $n\ge 3$. Notice that the corresponding result does not hold in general when $m=1$.
In this paper we study inverse boundary value problems with partial data for the magnetic Schrödinger operator. In the case of an infinite slab in $R^n$, $n\ge 3$, we establish that the magnetic field and the electric potential can be determined uniquely, when the Dirichlet and Neumann data are given either on the different boundary hyperplanes of the slab or on the same hyperplane. This is a generalization of the results of [41], obtained for the Schrödinger operator without magnetic potentials. In the case of a bounded domain in $R^n$, $n\ge 3$, extending the results of [2], we show the unique determination of the magnetic field and electric potential from the Dirichlet and Neumann data, given on two arbitrary open subsets of the boundary, provided that the magnetic and electric potentials are known in a neighborhood of the boundary. Generalizing the results of [31], we also obtain uniqueness results for the magnetic Schrödinger operator, when the Dirichlet and Neumann data are known on the same part of the boundary, assuming that the inaccessible part of the boundary is a part of a hyperplane.
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.
We study the problem of recovery the source $a(t,x)F(x)$ in the wave equation in anisotropic medium with $a$ known so that $a(0,x)\not=0$ with a single measurement. We use Carleman estimates combined with geometric arguments and give sharp conditions for uniqueness. We also study the non-linear problem of recovery the sound speed in the equation $u_{tt} -c^2(x)Δu =0$ with one measurement. We give sharp conditions for stability, as well. An application to thermoacoustic tomography is also presented.
The advent of transformation optics and metamaterials has made possible devices producing extreme effects on wave propagation. Here we give theoretical designs for devices, Schrödinger hats, acting as invisible concentrators of waves. These exist for any wave phenomenon modeled by either the Helmholtz or Schrödinger equations, e.g., polarized waves in EM, pressure waves in acoustics and matter waves in QM, and occupy one part of a parameter space continuum of wave-manipulating structures which also contains standard transformation optics based cloaks, resonant cloaks and cloaked sensors. For EM and acoustic Schrödinger hats, the resulting centralized wave is a localized excitation. In QM, the result is a new charged quasiparticle, a \emph{quasmon}, which causes conditional probabilistic illusions. We discuss possible solid state implementations.
We study the asymptotic behavior of an incompressible fluid around a bounded obstacle. The problem is modeled by the stationary Navier-Stokes equations in an exterior domain in $\R^n$ with $n\ge 2$. We will show that, under some assumptions, any nontrivial velocity field obeys a minimal decaying rate $\exp(-Ct^2\log t)$ at infinity. Our proof is based on appropriate Carleman estimates.
It is generally believed that transformation optics based cloaking, besides rendering the cloaked region invisible to detection by scattering of incident waves, also shields the region from those same waves. We demonstrate a coupling between the cloaked and uncloaked regions, exposing a difference between cloaking for rays and waves. Interior resonances allow this coupling to be amplified, and careful choice of parameters leads to effective cloaks with degraded shielding. As one application, we describe how to use transformation optics to hide sensors in the cloaked region and yet enable the sensors to efficiently measure waves incident on the exterior of the cloak, an effect similar to the plasmon based approach of Alu' and Engheta.
Consider a real-analytic orientable connected complete Riemannian manifold $M$ with boundary of dimension $n\ge 2$ and let $k$ be an integer $1\le k\le n$. In the case when $M$ is compact of dimension $n\ge 3$, we show that the manifold and the metric on it can be reconstructed, up to an isometry, from the set of the Cauchy data for harmonic $k$-forms, given on an open subset of the boundary. This extends a result of [13] when $k=0$. In the two-dimensional case, the same conclusion is obtained when considering the set of the Cauchy data for harmonic $1$-forms. Under additional assumptions on the curvature of the manifold, we carry out the same program when $M$ is complete non-compact. In the case $n\ge 3$, this generalizes the results of [12] when $k=0$. In the two-dimensional case, we are able to reconstruct the manifold from the set of the Cauchy data for harmonic $1$-forms.
We prove that a potential $q$ can be reconstructed from the Dirichlet-to-Neumann map for the Schrodinger operator $-Δ_g + q$ in a fixed admissible 3-dimensional Riemannian manifold $(M,g)$. We also show that an admissible metric $g$ in a fixed conformal class can be constructed from the Dirichlet-to-Neumann map for $Δ_g$. This is a constructive version of earlier uniqueness results by Dos Santos Ferreira et al. on admissible manifolds, and extends the reconstruction procedure of Nachman in Euclidean space. The main points are the derivation of a boundary integral equation characterizing the boundary values of complex geometrical optics solutions, and the development of associated layer potentials adapted to a cylindrical geometry.
We show that the attenuated geodesic ray transform on two dimensional simple surfaces is injective. Moreover we give a stability estimate and develop a reconstruction procedure.
We establish a decoupling result for the $P$ and $S$ waves of linear, isotropic elasticity, in the setting of twice-differentiable Lamé parameters. Precisely, we show that the $P\leftrightarrow S$ components of the wave propagation operator are regularizing of order one on $L^2$ data, by establishing the diagonalization of the elastic system modulo a $L^2$-bounded operator. Effecting the diagonalization in the setting of twice-differentiable coefficients depends upon the symbol of the conjugation operator having a particular structure.
We give a detailed microlocal study of X-ray transforms over geodesics-like families of curves with conjugate points of fold type. We show that the normal operator is the sum of a pseudodifferential operator and a Fourier integral operator. We compute the principal symbol of both operators and the canonical relation associated to the Fourier integral operator. In two dimensions, for the geodesic transform, we show that there is always a cancellation of singularities to some order, and we give an example where that order is infinite; therefore the normal operator is not microlocally invertible in that case. In the case of three dimensions or higher if the canonical relation is a local canonical graph we show microlocal invertibility of the normal operator. Several examples are also studied.
We study the mathematical model of thermoacoustic and photoacoustic tomography when the sound speed has a jump across a smooth surface. This models the change of the sound speed in the skull when trying to image the human brain. We derive an explicit inversion formula in the form of a convergent Neumann series under the assumptions that all singularities from the support of the source reach the boundary.
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.
We study the mathematical model of thermoacoustic tomography in media with a variable speed for a fixed time interval, greater than the diameter of the domain. In case of measurements on the whole boundary, we give an explicit solution in terms of a Neumann series expansion. We give necessary and sufficient conditions for uniqueness and stability when the measurements are taken on a part of the boundary.