Pith. sign in

REVIEW 4 major objections 3 minor 23 references

The boundary control approach to the Titchmarsh-Weyl $m-$function

T0 review · 4 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read One integral equation computes the Titchmarsh-Weyl m-function

desk verdict A genuinely new bound for the A-amplitude under weak decay, buried under a misstated integral equation and an unproved identification step; worth refereeing after revision. read the letter →

arxiv 2505.23332 v1 pith:JPG73OW6 submitted 2025-05-29 math.AP math.SP

S.A. Avdonin , V.S. Mikhaylov , A.V. Rybkin This is my paper
classification math.APmath.SP MSC 34B2034E0534L2534E4047B2081Q10
keywords SchrödingeroperatorTitchmarsh-Weylm-functionboundarycontrolA-amplituderesponseDirichlet-to-NeumannmapVolterraintegralequationinversespectraltheory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper shows that the Titchmarsh-Weyl m-function of a Schrödinger operator on the half-line, which encodes the spectral measure and is difficult to compute at finite spectral values, can be evaluated through the controlled wave equation associated with the operator. The m-function is the Laplace transform of the kernel of the Dirichlet-to-Neumann map (the response operator) of that wave equation, and that kernel is, up to scaling, the A-amplitude of the Gesztesy–Simon representation. The paper proves that the A-amplitude is the diagonal value of the solution of a linear Volterra integral equation, yielding an algorithm: solve once for the amplitude, then take Laplace transforms for each desired point. It also proves an explicit exponential bound for the A-amplitude under the condition that the L1-norm of the potential on unit intervals is uniformly bounded, answering a conjecture of Gesztesy and Simon.

What carries the argument

The central object is the response operator $R$ (the dynamic Dirichlet-to-Neumann map) of the wave equation associated with the Schrödinger operator, together with its kernel $r(s)=w_x(0,s)$, defined through the solution $w$ of the Goursat problem $w_{ss}-w_{xx}+q(x)w=0$, $w(x,0)=0$, $w(x,x)=-\frac12\int_0^x q$. The A-amplitude is the rescaled kernel $A(\alpha)=-2r(2\alpha)$. The argument's engine is the linear Volterra integral equation (3.2) for $A(x,y)$, obtained by a chain of variable changes in the Goursat problem; its diagonal value $A(\alpha,\alpha)$ is the A-amplitude. Iteration of the Volterra operator $K$, estimated with Lemma 1 and Stirling's bound, produces the explicit exponential estimate (4.5).

What would settle it

For a real potential satisfying (4.4) but with $\int|q|=\infty$ (for instance $q(x)=\sin(x^2)$), solve (3.2) by iteration to obtain $A(\alpha)=A(\alpha,\alpha)$, insert it into (2.15), and compare $m(-k^2)$ with a direct high-precision solution of (1.10)–(1.4) at a fixed $k$ with $\operatorname{Re}k>2\max\{\sqrt{2\|q\|},e\|q\|\}$; a disagreement beyond the bound (4.5) would show the integral equation does not yield the true A-amplitude.

Watch

Extended reading notes

Core claim

The paper establishes that the Dirichlet Titchmarsh–Weyl m-function associated with $H=-\partial_x^2+q(x)$ on $L^2(0,\infty)$ is the Laplace transform of the response function of the wave equation $u_{tt}-u_{xx}+q(x)u=0$ with zero initial data and boundary control $u(0,t)=f(t)$. Concretely, $m(-k^2)=-k+\int_0^\infty e^{-k\alpha}r(\alpha)\,d\alpha$, and equivalently $m(-k^2)=-k-\int_0^\infty A(\alpha)e^{-2\alpha k}\,d\alpha$, where $A(\alpha)=-2r(2\alpha)$ is the A-amplitude. The main new results are: first, $A(\alpha)$ equals the diagonal value $A(\alpha,\alpha)$ of the solution of the linear Volterra integral equation $A(x,y)=q(x)-\int_0^y\big(\int_x^v A(u,v)\,du\big)q(x-v)\,dv$, giving a direct algorithm to evaluate $m$ by solving once and then computing a Laplace transform; and second, under the condition $\|q\|:=\sup_x\int_x^{x+1}|q(s)|\,ds<\infty$, the error bound $|A(\alpha)-q(\alpha)|\le \frac12\big(\int_0^\alpha|q|\big)^2\big[e^{2\sqrt{2\|q\|\alpha}}+\tfrac1{\sqrt{2\pi}}e^{2e\|q\|\alpha}\big]$ holds, yielding absolute convergence of the Laplace integral for $\operatorname{Re} k>2\max\{\sqrt{2\|q\|},e\|q\|\}$ and settling the conjecture that the A-amplitude has an exponential bound for such potentials.

Load-bearing premise

The whole construction assumes that for every $q\in L^1_{\rm loc}$ the wave-equation solution has the representation (2.2) with $w$ solving the Goursat problem (2.3) and that $w_x(0,\cdot)$ is the response kernel; the paper does not prove existence or uniqueness of such $w$ for general $L^1_{\rm loc}$ potentials, and the convergence of the Volterra iteration is established only under the stronger uniform one-step $L^1$ bound (4.4).

Editorial extensions

If this is right

  • The m-function can be evaluated by solving one linear Volterra equation (3.2) and taking Laplace transforms; once $A(\alpha)$ is known, every new spectral point costs one transform rather than a separate ODE solve.
  • The A-amplitude acquires a concrete physical meaning: it is (minus twice) the impulse response kernel of the Dirichlet-to-Neumann map for the wave equation.
  • The representation (2.15) and its absolute convergence are extended to all potentials with uniformly bounded local $L^1$-norm, including potentials with no decay at infinity.
  • For nonnegative potentials the terms of the alternating series (5.1) are nonnegative, so the algorithm's series converges faster and truncation is easier to control.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Inverting the integral equation (3.2) would recover the potential $q$ from the boundary response kernel, suggesting a new inverse spectral algorithm for the half-line that bypasses Gelfand–Levitan–Marchenko machinery.
  • The Dirichlet-to-Neumann identification of the m-function is naturally multidimensional, so the same response-operator route could define operator-valued m-functions for $-\Delta+q$ in higher dimensions; the paper notes this possibility in Remark 2 but does not develop it.
  • The exponential estimate (4.5) doubles as an error bound for truncated Neumann series, giving a stopping criterion for numerical evaluation of the m-function when $\|q\|$ is moderate.
  • If the iteration of (3.2) converges for potentials outside the class (4.4), the Laplace representation may hold in a larger region, so the convergence condition rather than the Volterra equation itself is the current bottleneck.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 3 minor

Summary. The paper proposes a boundary-control approach to the Titchmarsh-Weyl m-function for the half-line Schrödinger operator with locally integrable potential. It introduces the response operator R of the wave equation (2.1), identifies its convolution kernel r with w_x(0,·), connects r to Simon's A-amplitude via (2.16), derives the linear Volterra equation (3.2) for A(x,y), and proves an exponential bound (4.5) for the diagonal A(α)=A(α,α) under the condition (4.4). It then presents Algorithm 1, which evaluates m by solving (3.2) and taking the Laplace transform (2.15).

Significance. If the identification between the formal solution of (3.2) and the response kernel was made rigorous, the paper would provide a new interpretation of the A-amplitude and a potentially efficient numerical procedure for the m-function. The explicit bound (4.5) with concrete constants is a substantive contribution and would answer a conjecture of Gesztesy-Simon for potentials with bounded L1 norm on unit intervals. The combination of the Goursat problem and boundary-control ideas is elegant and likely to be useful. However, the load-bearing analytic steps are currently asserted rather than proved, so the central claims are not yet established at the stated level of generality.

major comments (4)
  1. [§2, Eq. (2.2)–(2.3), (2.6)] The representation (2.2) of the weak solution to (2.1) in terms of a solution w of the Goursat problem (2.3), and the identification r(·)=w_x(0,·), are asserted without proof. For q merely in L1_loc, w is not shown to be differentiable up to the boundary x=0, nor is w_x(0,·) shown to be a locally integrable function whose Laplace transform can be used in (2.14). This is not a technicality: Theorem 1 and representation (2.15) both rely on this identification.
  2. [Theorem 1, Eq. (3.2)] Equation (3.2) is stated for x,y>0, but the term q(x-v) is evaluated at negative arguments whenever v>x, so the equation is not defined as written. The derivation from (3.5) is valid only on the triangular domain 0≤y≤x. The theorem should state that domain and prove existence and uniqueness of the solution A(x,y) there, for example by the Volterra iteration. Without this, the use of A(α,α) in Theorem 2 and Algorithm 1 is not justified.
  3. [Theorem 2, proof of (4.5) and of (2.15)] The proof estimates the formal Neumann series for (3.2) and concludes the bound for the A-amplitude. This conclusion requires that the actual A-amplitude defined by (2.16) satisfies (3.2) and equals the sum of the Neumann series. Under the stated assumption (4.4), neither the existence of a solution to (3.2) nor its coincidence with -2r(2α) is proved. The bound therefore applies to the formal series, not yet to the spectral object appearing in (2.15).
  4. [§5, Algorithm 1 and Remark 4] Algorithm 1 instructs to evaluate m(z) by (2.15) after computing A(α). The paper proves absolute convergence of the integral in (2.15) only for Re k > 2 max{√(2||q||), e||q||}, and it does not prove the identity m(-k^2) = -k - ∫_0∞ A(α)e^{-2αk}dα for all such k; the derivation via (2.11)–(2.14) inherits the gaps described above. Moreover, for a given finite z∈C+, the k with -k^2=z may have Re k below that threshold, so the algorithm's range of applicability to finite z is not demonstrated.
minor comments (3)
  1. [§2, Eq. (2.10)] Equation (2.10) appears to contain a typographical error: the intended statement is likely \(\widehat{Rf}(k)=\widehat{u_x}(0,k)\).
  2. [Lemma 1, Eq. (4.3)] The integration-by-parts proof uses the symbol b^n; for n=0 and b=0 the expression 0^0 should be avoided by treating n=0 separately or by taking n≥1.
  3. [References] In reference [12], the page range "491–436" appears to be printed in reverse order.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the derivation computes the A-amplitude from the potential via a Volterra integral equation and obtains m by a Laplace transform; no output quantity is used as an input.

full rationale

The derivation chain is self-contained and non-circular. The paper defines the Titchmarsh-Weyl m-function in the standard spectral way (1.4), introduces the response operator R of the associated wave equation, identifies its kernel r with w_x(0,·), and establishes the Laplace-domain relation (2.12)-(2.13); the A-amplitude is then defined by (2.16) as -2r(2α), and Theorem 1 derives the Volterra equation (3.2) for A from the Goursat problem (2.3) using only the potential q. Theorem 2 bounds the resulting object, and Algorithm 1 evaluates m by (2.15), i.e., by a Laplace transform of an object computed from q. No fitted parameter enters, no 'prediction' is a renaming of an input, and no part of the proof assumes the target m-function or the A-amplitude representation it claims to establish. The self-citations ([20], [21]) appear in introductory/comparative remarks and are not load-bearing for Theorem 1, Theorem 2, or Algorithm 1. The skeptical concern about existence and uniqueness of the Goursat problem and identification of the iterated solution of (3.2) with the response kernel for L1_loc potentials is a well-posedness and rigor gap, not circularity.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper relies on standard Titchmarsh-Weyl and boundary control facts, introduces no free parameters fitted to data, and invents no new entities. The main unproved structural assumption is the Goursat/wave representation for L1_loc potentials.

assumptions (4)
  • domain assumption The Schrödinger operator is in the limit point case at infinity, ensuring a unique Weyl solution u_+ for z∈C+.
    Invoked at the start, Eq (1.2)-(1.4), to define the Titchmarsh-Weyl m-function.
  • domain assumption The wave equation (2.1) with locally integrable q has a weak solution of the form (2.2), with w solving the Goursat problem (2.3).
    This representation is the basis for the response operator and the kernel r; the paper states it can be verified by direct computation but does not prove existence for the full L1_loc class.
  • standard math The relation bu_x(0,k) = m(-k^2) bf(k) between the Laplace-domain Dirichlet-to-Neumann map and the m-function holds.
    Used to derive (2.12) and the representation (2.14); standard in Titchmarsh-Weyl theory but not proved in the paper.
  • standard math Volterra iteration and the norm estimate in Lemma 1 are valid for the integral operator K defined in Section 4.
    These are used to prove Theorem 2 and the convergence of the series for A under condition (4.4).

how reviews work

0 comments
Cite this review

Pith. "Pith review of The boundary control approach to the Titchmarsh-Weyl $m-$function." pith.science (2026). https://pith.science/paper/JPG73OW6

@misc{pith2026250523332,
  author       = {Pith},
  title        = {Pith review of: The boundary control approach to the Titchmarsh-Weyl $m-$function},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JPG73OW6}},
  note         = {Machine review of arXiv:2505.23332}
}
abstract

We link the Boundary Control Theory and the Titchmarsh-Weyl Theory. This provides a natural interpretation of the $A-$amplitude due to Simon and yields a new efficient method to evaluate the Titchmarsh-Weyl $m-$function associated with the Schr\"{o}dinger operator $H=-\partial _{x}^{2}+q\left( x\right) $ on $L_{2}\left( 0,\infty \right) $ with Dirichlet boundary condition at $x=0.$

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    Amrein and D.B

    W.O. Amrein and D.B. Pearson,Moperators: a generalization of Weyl-Titchmarsh theory J. Comput. Appl. Math. 171 (2004), no. 1-2, 1–26

    work page 2004

  2. [2]

    Avdonin and M.I

    S.A. Avdonin and M.I. Belishev,Boundary control and dynamic inverse problem for non- selfadjoint Sturm–Liouville operator, Control and Cybernetics 25 (1996), 429–440

    work page 1996

  3. [3]

    Avdonin, M.I

    S.A. Avdonin, M.I. Belishev, and S. A. Ivanov,Matrix inverse problem for the equation utt −u xx +Q(x)u= 0, Math. USSR Sbornik 7 (1992), 287–310

    work page 1992

  4. [4]

    Borg,Uniqueness theorems in the spectral theory ofy ′′ + (λ−q(x))y= 0.In: Proc

    G. Borg,Uniqueness theorems in the spectral theory ofy ′′ + (λ−q(x))y= 0.In: Proc. 11th Scandinavian Congress of Mathematicians. Oslo: Johan Grundt Tanums Forlag, 1952, pp. 276–287

    work page 1952

  5. [5]

    Buslaev and V.B

    V.S. Buslaev and V.B. Matveev,Wave operators for the Schr¨ odinger equation with slowly decreasing potential, Theoret. and Math. Phys. 2 (1970), no. 3, 266–274

    work page 1970

  6. [6]

    Chadan and P.C

    K. Chadan and P.C. Sabatier,Inverse problems in quantum scattering theory. Second edition. Texts and Monographs in Physics. Springer-Verlag, New York, 1989. xxxii+499 pp

    work page 1989

  7. [7]

    Clark and F

    S. Clark and F. Gesztesy,Weyl-TitchmarshM-function asymptotics, local uniqueness re- sults, trace formulas, and Borg-type theorems for Dirac operators, Trans. Amer. Math. Soc. 354 (2002), no. 9, 3475–3534

    work page 2002

  8. [8]

    Freiling, V

    G. Freiling, V. Yurko,Inverse Sturm-Liouville problems and their applications.Nova Science Publishers, Inc., Huntington, NY, 2001. x+356 pp

    work page 2001

Show all 23 references
  1. [9]

    Gesztesy and B

    F. Gesztesy and B. Simon,A new approach to inverse spectral theory. II. General real po- tentials and the connection to the spectral measure,Ann. of Math. (2) 152 (2000), no. 2, 593–643. 13

    2000

  2. [10]

    B. J. Harris,An exact method for the calculation of certain Titchmarsh-Weylm-functions. Proc. Roy. Soc. Edinburgh Sect. A 106 (1987), no. 1-2, 137–142

    1987

  3. [11]

    Hinton, M

    D.B. Hinton, M. Klaus, and J.K. Shaw,Series representation and asymptotics for Titchmarsh-Weylm-functions.,Differential Integral Equations 2 (1989), no. 4, 419–429

    1989

  4. [12]

    Kachalov, Y

    A. Kachalov, Y. Kurylev, M. Lassas and N. Mandache,Equivalence of time-domain inverse problems and boundary spectral problems, Inverse Problems, 20 (2004), 491–436

    2004

  5. [13]

    Kaper and M.K

    H.G. Kaper and M.K. Kwong,Asymptotics of the Titchmarsh-Weylm-coefficient for in- tegrable potentials. II.Differential equations and mathematical physics (Birmingham, Ala., 1986), 222–229, Lecture Notes in Math., 1285, Springer, Berlin, 1987

    1986

  6. [14]

    Levitan, I.S

    B.M. Levitan, I.S. Sargsjan,Introduction to spectral theory: selfadjoint ordinary differen- tial operators.Translations of Mathematical Monographs, Vol. 39. American Mathematical Society, Providence, R.I., 1975. xi+525 pp

    1975

  7. [15]

    Marchenko,Certain problems in the theory of second-order differential operators.Dok- lady Akad

    V.A. Marchenko,Certain problems in the theory of second-order differential operators.Dok- lady Akad. Nauk SSSR 72, 457–460 (1950) (Russian)

    1950

  8. [16]

    Matveev and M.M

    V.B. Matveev and M.M. Skriganov,Scattering problem for radial Schr¨ odinger equation with slowly decreasing potential,Teor. Mat. Fiz. 10 (1972), no. 2, 238–248

    1972

  9. [17]

    Pavlov,S-Matrix and Dirichlet-to-Neumann Operators, In: Encyclopedia of Scattering, ed

    B. Pavlov,S-Matrix and Dirichlet-to-Neumann Operators, In: Encyclopedia of Scattering, ed. R. Pike, P. Sabatier, Academic Press, Harcourt Science and Tech. Company (2001) pp 1678-1688

    2001

  10. [18]

    Ramm,Recovery of the potential fromI−functionC

    A.G. Ramm,Recovery of the potential fromI−functionC. R. Math. Rep. Acad. Sci. Canada 9 (1987), no. 4, 177–182

    1987

  11. [19]

    Ramm and B

    A. Ramm and B. Simon,A new approach to inverse spectral theory. III. Short-range poten- tialsJ. Anal. Math. 80 (2000), 319–334

    2000

  12. [20]

    A. Rybkin,On a transformation of the Sturm-Liouville equation with slowly decaying poten- tials and the Titchmarsh-Weylm−functionSpectral methods for operators of mathematical physics, 185–201, Oper. Theory Adv. Appl., 154, Birkh¨ auser, Basel, 2004

    2004

  13. [21]

    Rybkin,Some new and old asymptotic representations of the Jost solution and the Weyl m-function for Schr¨ odinger operators on the line

    A. Rybkin,Some new and old asymptotic representations of the Jost solution and the Weyl m-function for Schr¨ odinger operators on the line. Bull. London Math. Soc. 34 (2002), no. 1, 61–72

    2002

  14. [22]

    Simon,A new approach to inverse spectral theory

    B. Simon,A new approach to inverse spectral theory. I. Fundamental formalism,Ann. of Math. 150 (1999), no. 2, 1029-1057

    1999

  15. [23]

    E. C. Titchmarsh,Eigenfunction expansions associated with second-order differential equa- tions. Part I,Second Edition, Clarendon Press, Oxford 1962, 203 pp. Department of Mathematics and Statistics, University of Alaska F airbanks, PO Box 756660, F airbanks, AK 99775 Email ad...

    1962

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.