Plamen Stefanov

Inverse Boundary Problem for the Two Photon Absorption Transport Equation

This work studies the inverse boundary problem for the two photon absorption radiative transport equation. We show that the absorption coefficients and scattering coefficients can be uniquely determined from the \emph{albedo} operator. If scattering is absent, we do not require smallness of the incoming source and the reconstructions of the absorption coefficients are explicit.

Recovery of a cubic non-linearity in the wave equation in the weakly non-linear regime

We study the inverse problem of recovery a compactly supported non-linearity in the semilinear wave equation $u_{tt}-Δu+ α(x) |u|^2u=0$ , in two and three dimensions. We probe the medium with complex-valued harmonic waves of wavelength $h$ and amplitude $h^{-1/2}$ , then they propagate in the weakly non-linear regime; and measure the transmitted wave when it exits the support of $α$ . We show that one can extract the Radon transform of $α$ from the phase shift of such waves, and then one can recover $α$ . We also show that one can probe the medium with real-valued harmonic waves and obtain uniqueness for the linearized problem.

The Radon transform with finitely many angles

We study the Radon transform in the plane in parallel geometry possibly undersampled in the angular variables. We study resolution, aliasing artifacts, and edge recovery.

Weakly nonlinear geometric optics for the Westervelt equation and recovery of the nonlinearity

We study the non-diffusive Westervelt equation in the weakly nonlinear regime. We show that the leading profile equation is of Burgers' type. We show that a compactly supported nonlinearity $α$ can be reconstructed from the tilt of the transmitted high frequency wave packets sent from different directions since those tilts are proportional to the X-ray transform of $α$ .

preprint2020arXiv

Inverse problems for real principal type operators

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.

preprint2019arXiv

The Light Ray transform on Lorentzian manifolds

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.

preprint2016arXiv

Multiwave tomography with reflectors: Landweber's iteration

We use the Landweber method for numerical simulations for the multiwave tomography problem with a reflecting boundary and compare it with the averaged time reversal method. We also analyze the rate of convergence and the dependence on the step size for the Landweber iterations on a Hilbert space.

preprint2016arXiv

Thermo and photoacoustic Tomography with variable speed and planar detectors

We analyze the mathematical model of multiwave tomography with a variable speed with integrating measurements on planes tangent to a sphere surrounding the source. We prove sharp uniqueness and stability estimates with full and partial data and propose a time reversal algorithm which recovers the visible singularities.

Boundary rigidity with partial data

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.

Multiwave tomography in a closed domain: averaged sharp time reversal

We study the mathematical model of multiwave tomography including thermo and photoacoustic tomography with a variable speed for a fixed time interval $[0,T]$ . We assume that the waves reflect from the boundary of the domain. We propose an averaged sharp time reversal algorithm. In case of measurements on the whole boundary, we give an explicit solution in terms of a Neumann series expansion. When the measurements are taken on a part of the boundary, we show that the same algorithm produces a parametrix.

On the inverse problem of finding cosmic strings and other topological defects

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.

On the stable recovery of a metric from the hyperbolic DN map with incomplete data

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.

Support theorems for the Light Ray transform on analytic Lorentzian manifolds

We study the weighted ray transform of integrating functions on a Lorentzian manifold over lightlike geodesics. We prove support theorems if the manifold and the weight are analytic.

preprint2014arXiv

Inverting the local geodesic X-ray transform on tensors

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.

preprint2014arXiv

Modulated Luminescent Tomography

We propose and analyze a mathematical model of Modulated Luminescent Tomography. We show that when single X-rays or focused X-rays are used as an excitation, the problem is similar to the inversion of weighted X-ray transforms. In particular, we give an explicit inversion in the case of Dual Cone X-ray excitation.

preprint2014arXiv

The geodesic ray transform on Riemannian surfaces with conjugate points

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.

preprint2013arXiv

Stability of Coupled-Physics Inverse Problems with internal measurements

In this paper, we develop a general approach to prove stability for the non linear second step of hybrid inverse problems. We work with general functionals of the form $σ|\nabla u|^p$ , $0 < p \leq 1$ , where $u$ is the solution of the elliptic partial differential equation $\nabla\cdot σ\nabla u =0$ on a bounded domain $Ω$ with boundary conditions $u|_{\partial Ω} = f$ . We prove stability of the linearization and Hölder conditional stability for the non-linear problem of recovering $σ$ from the internal measurement.

preprint2013arXiv

Weyl asymptotics of the transmission eigenvalues for a constant index of refraction

We prove Weyl type of asymptotic formulas for the real and the complex internal transmission eigenvalues when the domain is a ball and the index of refraction is constant.

preprint2012arXiv

Instability of the linearized problem in multiwave tomography of recovery both the source and the speed

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.

preprint2012arXiv

Is a curved flight path in SAR better than a straight one?

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.

A New Numerical Algorithm for Thermoacoustic and Photoacoustic Tomography with Variable Sound Speed

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.

Recovery of a source term or a speed with one measurement and applications

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.

Stability for the multi-dimensional Borg-Levinson theorem with partial spectral data

We prove a stability estimate related to the multi-dimensional Borg-Levinson theorem of determining a potential from spectral data: the Dirichlet eigenvalues and the normal derivatives of the eigenfunctions on the boundary of a bounded domain. The estimate is of Hölder type, and we allow finitely many eigenvalues and normal derivatives to be unknown. We also show that if the spectral data is known asymptotically only, up to $O(k^{-a})$ with $a>>1$ , then we still have Hölder stability.