BZPEERReading, trust and collaboration for researchDiscover
Workspaces
Network
Opportunities
Account
Researcher profile
Gunther Uhlmann contributes to research discovery and scholarly infrastructure.
Trust snapshot
Actions
Identity and collaboration
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Claim this author entity first to unlock direct invitations.
Research graph
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph 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 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 survey recent progress in the problem of recovering a tensor field from its integrals along geodesics. We also propose several open problems.
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.
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 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$.
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.
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.
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 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.