Explicit inversion for variable-speed wave equations on bounded domains
Pith reviewed 2026-05-24 07:09 UTC · model grok-4.3
The pith
Explicit formulas recover spectral coefficients of initial pressure from boundary measurements for variable-speed waves.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within a unified framework, explicit formulas are presented that recover the spectral coefficients ⟨f,ϕ_k^B⟩ of f with respect to the eigenfunction bases of the operator −c(⋅)Δx for boundary types B∈{D,R} from time-resolved boundary measurements of p or its normal derivative under the respective boundary conditions.
What carries the argument
The explicit inversion formulas that extract the inner products ⟨f,ϕ_k^B⟩ directly from the time-resolved boundary traces by integrating against the eigenfunction expansion of the wave solution.
If this is right
- The initial pressure f can be reconstructed by summing the recovered coefficients against the eigenfunctions.
- The same set of formulas applies uniformly to both Dirichlet and Robin boundary conditions.
- Direct coefficient recovery becomes possible without iterative solution of the inverse problem.
- Variable sound speed c(x) is incorporated explicitly into the inversion without additional approximations.
Where Pith is reading between the lines
- Numerical implementation on simple domains such as intervals or disks could verify the formulas' accuracy for specific c(x).
- The direct recovery might lower computational demands compared to optimization-based methods in imaging applications.
- The framework could be examined for adaptation to other linear hyperbolic equations with similar spectral structures.
Load-bearing premise
The time-resolved boundary measurements contain the information needed to extract the spectral coefficients directly via the formulas, relying on the eigenfunctions of -c(x)Δ forming a basis under the given boundary conditions.
What would settle it
Simulate boundary data from the wave equation for a known f and c(x), apply the formulas to recover the coefficients, and check whether they match the true inner products ⟨f,ϕ_k^B⟩; mismatch would show the formulas do not hold.
read the original abstract
We study the reconstruction of the initial pressure $f(x)=p(x,0)$ for the wave model \[ \partial_t^2 p(x,t)=c(x)\Delta_{x}p(x,t)\qquad (x,t)\in\Omega\times[0,\infty), \] posed on a bounded domain $\Omega$ with variable sound speed $c(\cdot)$. From time-resolved boundary measurements, we consider two settings: (i) measurement of $p|_{\partial\Omega\times[0,\infty)}$ under a Robin boundary condition $p+\alpha\,\partial_\nu p=0$ on $\partial\Omega\times[0,\infty)$ with $\alpha\gneq 0$, and (ii) measurement of $\partial_\nu p|_{\partial\Omega\times[0,\infty)}$ under a Dirichlet boundary condition $p=0$ on $\partial\Omega\times[0,\infty)$. Within a unified framework, we present explicit formulas that recover the spectral coefficients $\langle f,\phi_k^B\rangle$ of $f$ with respect to the eigenfunction bases of the operator $-c(\cdot)\Delta_{x}$ for boundary types $B\in\{D,R\}$. The framework integrates variable sound speed with Dirichlet/Robin boundary conditions in a single setting, enabling direct coefficient-level recovery from boundary data.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives explicit formulas recovering the spectral coefficients ⟨f, ϕ_k^B⟩ of the initial pressure f with respect to the eigenfunctions of the operator −c(x)Δ (for B ∈ {D, R}) directly from time-resolved boundary traces (p on ∂Ω for Robin; ∂_ν p on ∂Ω for Dirichlet) of the solution to the variable-speed wave equation ∂_t²p = c(x)Δp on a bounded domain Ω.
Significance. If the derivations hold, the result supplies a direct, non-iterative recovery of individual eigen-coefficients from boundary data within a single framework that treats both Dirichlet and Robin conditions uniformly. This is a modest but concrete advance for coefficient-level inversion in the variable-coefficient wave equation, resting on standard self-adjoint spectral theory (completeness of the eigenbasis in L²(Ω) under the stated boundary conditions and positivity/smoothness assumptions on c).
minor comments (3)
- The abstract and introduction state that the formulas are 'explicit,' but the precise sense (closed-form expressions involving only the boundary traces, known c, and the eigenvalues, without auxiliary solves) should be clarified in §2 or §3 with an example computation for the constant-c case to illustrate the reduction.
- Notation for the two boundary conditions is introduced in the abstract but the precise Robin parameter α (constant or function) and the domain regularity (C^∞ or C^{2,α}) are not restated in the main theorems; add a short paragraph in §1.2 listing all standing assumptions on Ω and c.
- The time-harmonic expansion used to isolate the coefficients is standard, yet the passage from the boundary trace to the coefficient formula (likely via Laplace transform or Fourier sine transform in time) should be written with an explicit reference to the relevant identity or lemma number.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of the manuscript and the recommendation for minor revision. The report correctly identifies the core contribution as explicit, non-iterative recovery of individual eigen-coefficients from boundary traces for the variable-speed wave equation under both boundary conditions.
Circularity Check
No significant circularity identified
full rationale
The paper presents a mathematical derivation of explicit inversion formulas recovering spectral coefficients ⟨f, ϕ_k^B⟩ directly from boundary traces of the variable-speed wave equation. The supporting elements—self-adjointness of −c(x)Δ, completeness of the eigenbasis in L²(Ω), and time-harmonic expansion—are standard results in spectral theory for elliptic operators on bounded domains with the stated boundary conditions; they are not derived from or equivalent to the target formulas. No fitted parameters are renamed as predictions, no self-citations form a load-bearing chain, and no ansatz is smuggled via prior work. The logical chain from the PDE and boundary data to the coefficient formulas is therefore self-contained.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Eigenfunctions of the operator -c(x)Δ_x form a complete basis allowing spectral expansion and coefficient recovery from boundary data.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We provide explicit formulas that recover the spectral coefficients ⟨f,ϕ_k^B⟩ of f with respect to the eigenfunction bases of the operator −c(⋅)Δ_x for boundary types B∈{D,R}.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Proposition 2. ... T is a compact self-adjoint operator on L²(Ω,c(x)⁻¹dx).
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
M. Agranovsky and P. Kuchment. Uniqueness of reconstructio n and an inversion procedure for ther- moacoustic and photoacoustic tomography with variable sound spe ed. Inverse Problems , 23(5):2089, 2007
work page 2089
- [2]
- [3]
-
[4]
H Ammari, H Kang, and S Kim. Sharp estimates for the neumann fun ctions and applications to quanti- tative photo-acoustic imaging in inhomogeneous media. Journal of Differential Equations , 253(1):41–72, 2012
work page 2012
-
[5]
M. A Anastasio, J. Zhang, D. Modgil, and P. J La Rivi` ere. Applicat ion of inverse source concepts to photoacoustic tomography. Inverse Problems , 23(6):S21, 2007
work page 2007
-
[6]
A.G. Bell. On the production and reproduction of sound by light. American Journal of Science , 20:305–324, October 1880
-
[7]
H. Brezis. Functional Analysis, Sobolev Spaces and Partial Differenti al Equations. Universitext. Springer New York, 2010
work page 2010
-
[8]
F. Dreier and M. Haltmeier. Explicit inversion formulas for the two- dimensional wave equation from neumann traces. SIAM Journal on Imaging Sciences , 13(2):589–608, 2020
work page 2020
-
[9]
A. Erdelyi. Tables of Integral Transforms, Vols. I and II . Batemann Manuscript Project, McGrawHill, New York, 1954
work page 1954
-
[10]
L.C. Evans. Partial Differential Equations . Graduate studies in mathematics. American Mathematical Society, 2010
work page 2010
-
[11]
M. Haltmeier, O. Scherzer, P. Burgholzer, and G. Paltauf. The rmoacoustic computed tomography with large planar receivers. Inverse Problems , 20(5):1663–1673, 2004-10-01T00:00:00
work page 2004
-
[12]
Temporal backward projection of optoacoustic pressure transients using fourier transform methods
K P K¨ ostli, M Frenz, H Bebie, and H P Weber. Temporal backward projection of optoacoustic pressure transients using fourier transform methods. Physics in Medicine & Biology , 46(7):1863, 2001
work page 2001
- [13]
- [14]
- [15]
-
[16]
M. Moon, I. Hur, and S. Moon. Singular value decomposition of th e wave forward operator with radial variable coefficients. SIAM Journal on imaging science , To appear
-
[17]
P. Stefanov and G. Uhlmann. Thermoacoustic tomography with variable sound speed. Inverse Problems, 25(7):075011, 2009
work page 2009
-
[18]
M. Xu and L.V. Wang. Universal back-projection algorithm for p hotoacoustic computed tomography. Physical Review E , 71(1):016706, 2005. 7
work page 2005
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.