BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Plamen Stefanov 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
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.
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.
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.
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 $α$.
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.
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.
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.
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.
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.
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.
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 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 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.
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 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.
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.
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.
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.
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.
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 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 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.
We study the problem of recovery both the attenuation $a$ and the source $f$ in the attenuated X-ray transform in the plane. We study the linearization as well. It turns out that there are natural Hamiltonian flow that determines which singularities we can recover. If the perturbations $δa$, $δf$ are supported in a compact set that is non-trapping for that flow, then the problem is well posed. Otherwise, it may not be, and least in the case of radial $a$, $f$, it is not. We present uniqueness and non-uniqueness results for both the linearized and the non-linear problem; as well as a Hölder stability estimate.
In the inverse stationary transport problem through anisotropic attenuating, scattering, and refractive media, the albedo operator stably determines the gauge equivalent class of the attenuation and scattering coefficients.
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.
For anisotropic attenuating media, the albedo operator determines the scattering and the attenuation coefficients up to a gauge transformation. We show that such a determination is stable.
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.