Helmholtz-Schrödinger Equation Overview

• Helmholtz–Schrödinger equation is a family of stationary wave equations that interpolates between the Helmholtz operator and the time-independent Schrödinger operator.
• Key methodologies include analytical transformations, similarity reductions, and canonical form derivations to construct exactly solvable potentials.
• Its applications span inverse boundary value problems, quantum simulations, and advanced numerical techniques for high-frequency wave propagation.

The Helmholtz–Schrödinger equation denotes a family of stationary wave equations that interpolate between the Helmholtz operator and the time-independent Schrödinger operator, often at fixed frequency or fixed energy. In the literature covered here, the term is not tied to a single canonical PDE: it includes the semilinear fixed-frequency boundary value problem

$-\Delta u-k^2u+q(x)u^2=0,$

the one-dimensional stationary Schrödinger equation

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$

and variable-coefficient formulations of

$-\Delta \psi-k^2(x)\psi=f$

with $k^2(x)=E-V(x)$ or $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$ (Lu, 2022, Dong et al., 2017, Bonazzoli et al., 2017, Yu et al., 1 May 2025). A related paraxial usage arises when the Helmholtz equation is reduced, by a multiple-scaling ansatz, to a Schrödinger equation for a slowly varying envelope (Klumpp et al., 2021).

1. Terminological scope and canonical forms

A central feature of the terminology is that “Helmholtz” and “Schrödinger” describe the same stationary operator from different viewpoints. In one dimension, the time-independent Schrödinger equation

$\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$

is explicitly of Helmholtz form, with local effective wave number $\kappa^2(x)=k^2-V(x)$ (Dong et al., 2017). In higher dimensions, the general Helmholtz equation is

$\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$

so the stationary Schrödinger equation is recovered by identifying the potential with an energy shift in the coefficient $\kappa^2$ (Dong et al., 2017).

In PDE and inverse problems, the term is used for fixed-frequency elliptic equations whose linear part is the Helmholtz operator $-\Delta-k^2$ and whose interpretation is Schrödinger-like after rewriting $\psi''(x)+[k^2-V(x)]\psi(x)=0,$0 as $\psi''(x)+[k^2-V(x)]\psi(x)=0,$1 with $\psi''(x)+[k^2-V(x)]\psi(x)=0,$2. The nonlinear model studied in the partial-data Calderón setting is

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$3

with $\psi''(x)+[k^2-V(x)]\psi(x)=0,$4, $\psi''(x)+[k^2-V(x)]\psi(x)=0,$5, bounded with smooth boundary, $\psi''(x)+[k^2-V(x)]\psi(x)=0,$6 fixed outside a discrete resonance set, and $\psi''(x)+[k^2-V(x)]\psi(x)=0,$7 or more generally $\psi''(x)+[k^2-V(x)]\psi(x)=0,$8 (Lu, 2022).

In computational work, the same identification is used to transfer Helmholtz solvers to stationary quantum problems. The stationary Schrödinger equation

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$9

is written as

$-\Delta \psi-k^2(x)\psi=f$0

or, in nonrelativistic normalization,

$-\Delta \psi-k^2(x)\psi=f$1

making it a variable-coefficient Helmholtz problem (Bonazzoli et al., 2017, Yu et al., 1 May 2025).

A common source of ambiguity is therefore terminological rather than mathematical. Inverse problems, spectral theory, numerical linear algebra, and scattering theory all use closely related stationary operators, but emphasize different structures: fixed frequency, fixed energy, semilinearity, variable refractive index, or potential recovery. This suggests that “Helmholtz–Schrödinger equation” is best understood as a class label for stationary elliptic wave equations rather than a unique normal form.

2. Transformations, invariants, and operator reductions

A major analytical theme is the reduction of Helmholtz–Schrödinger equations to alternative normal forms. For the one-dimensional equation

$-\Delta \psi-k^2(x)\psi=f$2

a similarity transformation $-\Delta \psi-k^2(x)\psi=f$3 with $-\Delta \psi-k^2(x)\psi=f$4 yields

$-\Delta \psi-k^2(x)\psi=f$5

Its invariant is

$-\Delta \psi-k^2(x)\psi=f$6

and the Schwarzian derivative of the coordinate change,

$-\Delta \psi-k^2(x)\psi=f$7

enters the identity

$-\Delta \psi-k^2(x)\psi=f$8

which the paper calls the Schrödinger invariant (Dong et al., 2017). This furnishes a constructive route to exactly solvable potentials: specify a target ODE such as the Heun equation, choose $-\Delta \psi-k^2(x)\psi=f$9, and recover $k^2(x)=E-V(x)$0 from $k^2(x)=E-V(x)$1.

The same paper uses the general Heun equation as the target canonical form,

$k^2(x)=E-V(x)$2

and derives families of solvable potentials from the quadratic relation

$k^2(x)=E-V(x)$3

Special choices produce exponential, hyperbolic, trigonometric, quadratic, and linear changes of variables, with resulting potentials described as Heun–Eckart type, Heun–Pöschl–Teller I type, Heun–Pöschl–Teller II type, rational, and Coulomb-plus-centrifugal type (Dong et al., 2017).

A different reduction appears in paraxial wave theory. Starting from the Helmholtz equation

$k^2(x)=E-V(x)$4

and treating $k^2(x)=E-V(x)$5 as the evolution variable, a multiple-scaling ansatz

$k^2(x)=E-V(x)$6

with $k^2(x)=E-V(x)$7, leads at leading order to the Schrödinger equation

$k^2(x)=E-V(x)$8

The approximation is justified in Sobolev norms with an $k^2(x)=E-V(x)$9 error on propagation distances $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$0, despite the fact that the Helmholtz equation is ill-posed as an initial-value problem in $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$1 without a low-transverse-frequency restriction (Klumpp et al., 2021).

In axially symmetric settings, the stationary equation

$k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$2

admits a nonlocal Darboux transformation after the substitution $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$3, with

$k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$4

The resulting first-order system involves an auxiliary nonlocal variable $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$5, and in the special case $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$6 the transformation reduces to a generalized Moutard transformation adapted to cylindrical coordinates (Kudryavtsev, 2021).

3. Inverse boundary value theory

The most developed inverse-problem formulation in the supplied corpus is the partial-data Calderón problem for the nonlinear equation

$k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$7

in a bounded smooth domain $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$8 (Lu, 2022). For $k^2(x)=\frac{2m}{\hbar^2}(E-V(x))$9, a discrete exceptional set $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$0, and sufficiently small boundary data in $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$1, the forward problem

$\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$2

has a unique solution $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$3 with estimate

$\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$4

The nonlinear Dirichlet-to-Neumann map is

$\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$5

and the partial map restricts Dirichlet data to $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$6 and Neumann observations to $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$7. The principal uniqueness theorem states that if

$\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$8

on $\frac{d^2}{dx^2}\psi(x)+[k^2-V(x)]\psi(x)=0$9 for all sufficiently small $\kappa^2(x)=k^2-V(x)$0, $\kappa^2(x)=k^2-V(x)$1, $\kappa^2(x)=k^2-V(x)$2, then $\kappa^2(x)=k^2-V(x)$3 in $\kappa^2(x)=k^2-V(x)$4.

The method is higher-order linearization. With boundary data $\kappa^2(x)=k^2-V(x)$5, the first derivatives

$\kappa^2(x)=k^2-V(x)$6

solve the linear Helmholtz equation

$\kappa^2(x)=k^2-V(x)$7

The second derivatives

$\kappa^2(x)=k^2-V(x)$8

satisfy

$\kappa^2(x)=k^2-V(x)$9

After integration by parts with a third linear solution $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$0, one obtains the key identity

$\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$1

The density argument needed to force $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$2 has two levels. First, a local density theorem, proved by a Dos Santos Ferreira–Kenig–Sjöstrand–Uhlmann type microlocal construction, shows that if

$\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$3

for all solutions of

$\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$4

then $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$5 near a boundary point of $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$6. Second, Runge approximation and unique continuation propagate this vanishing along curves to the whole domain. The same strategy also yields simultaneous identification of unknown cavities and coefficients, and recovery of unknown outer boundaries together with the nonlinear coefficient (Lu, 2022).

4. Exact solutions, separability, and solvable geometries

Closed-form solutions are available in several highly structured settings. For the inhomogeneous Helmholtz equation in $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$7,

$\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$8

with the Sommerfeld radiation condition, the outgoing solution can be written explicitly in terms of $\nabla^2\psi+\kappa^2(\mathbf r)\psi=0,$9: $\kappa^2$0 The same paper derives a closed formula for the free Schrödinger evolution with initial data $\kappa^2$1, expressed by complex error functions, and extends these exact ball-data solutions to approximation formulas for general radial $\kappa^2$2-data by shell decomposition (Kirkeby, 2021).

Separation of variables can also be generalized beyond classical orthogonal systems. Polyelliptic coordinates are built by gluing several local elliptic coordinate systems associated with the sides of a convex polygon, producing a global orthogonal coordinate system outside the polygonal “protected area.” The metric assumes Stäckel form,

$\kappa^2$3

which yields separated ODEs for the 2D Helmholtz equation and, by extension, for the time-independent Schrödinger equation with piecewise constant potential (Kovalev, 2014). This framework is proposed for polygonal scattering and finite-wall polygonal quantum wells, including the square well.

Exactly solvable potential construction via Heun functions fits naturally into the same landscape. Choosing the Heun equation as canonical target and using the invariant relation

$\kappa^2$4

produces explicit potential families with Heun-function wavefunctions

$\kappa^2$5

(Dong et al., 2017).

5. Numerical solution methods and computational formulations

High-frequency Helmholtz–Schrödinger problems are dominated by indefiniteness, oscillation, and dispersion error. Domain decomposition methods address this through absorbed local problems and coarse corrections. For the interior Helmholtz problem

$\kappa^2$6

two-level ORAS/ImpRAS preconditioners combine local solves with a coarse space $\kappa^2$7, either from an absorbed coarse mesh or from local Dirichlet-to-Neumann eigenproblems (Bonazzoli et al., 2017). In the DtN construction, interface modes satisfy

$\kappa^2$8

and modes with $\kappa^2$9 are selected. The experiments use $-\Delta-k^2$0, subdomain diameter $-\Delta-k^2$1, and absorbed preconditioners with $-\Delta-k^2$2, where $-\Delta-k^2$3 performs markedly better than $-\Delta-k^2$4 (Bonazzoli et al., 2017).

A more recent nonoverlapping spectral additive Schwarz framework attacks the ill-posedness of local Dirichlet problems by splitting local interior modes through generalized eigenproblems

$-\Delta-k^2$5

and promoting small eigenmodes to a global coarse space (Yu et al., 1 May 2025). The resulting preconditioners $-\Delta-k^2$6 and $-\Delta-k^2$7 are purely algebraic and applicable to heterogeneous Helmholtz coefficients, with theoretical spectral clustering when thresholds are chosen close to zero.

Other solver architectures make the Helmholtz–Schrödinger relation explicit. A preconditioned iterative solver for scattering solutions of the Schrödinger equation uses partial-wave expansion, converting a high-dimensional driven Schrödinger equation into coupled 2D Helmholtz blocks and proposing a quadrant-definite preconditioner

$-\Delta-k^2$8

for the Helmholtz subproblems (Zubair et al., 2010). A time-domain preconditioner instead constructs a matrix recurrence whose steady state equals the discrete Helmholtz solution, then uses the resulting map as a GMRES preconditioner; this allows compact-stencil Helmholtz discretizations for which no native time-domain solver exists (Stolk, 2020). Far-field computation can also be reformulated on a complex contour, converting the far-field integral to an integral over a rotated domain where the scattered wave solves a damped problem equivalent to a Complex Shifted Laplacian or Complex Stretched Grid discretization, making multigrid effective in 2D and 3D Helmholtz and Schrödinger scattering (Cools et al., 2012).

Data-driven solvers follow the same structure. A plane-wave-activation neural network represents the solution as

$-\Delta-k^2$9

thereby generalizing the plane wave partition of unity method

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$00

by learning both amplitudes and directions (Wang et al., 2020). The loss enforces

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$01

at collocation points, and the method is positioned as directly relevant to constant-potential Schrödinger problems and as a starting point for variable-potential cases.

6. Quantum and operator-theoretic extensions

Quantum algorithms turn the discrete Helmholtz system

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$02

into a Schrödinger-type Hamiltonian evolution by “Schrödingerization.” The paper first forms the damped dynamical system

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$03

whose steady state is the desired solution. After rewriting it as a first-order system

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$04

and splitting the homogeneous lift $\psi''(x)+[k^2-V(x)]\psi(x)=0,$05 into $\psi''(x)+[k^2-V(x)]\psi(x)=0,$06, a warped phase transformation produces a higher-dimensional Schrödinger-type equation with Hamiltonian

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$07

The stated query complexity is $\psi''(x)+[k^2-V(x)]\psi(x)=0,$08, and for Helmholtz a simple preconditioner reduces the complexity to $\psi''(x)+[k^2-V(x)]\psi(x)=0,$09 at the level claimed in the paper (Gu et al., 31 Jul 2025).

A related, but conceptually different, operator-theoretic extension appears in compressible Navier–Stokes flow written in Schrödinger-type variables. There the velocity is Helmholtz-decomposed,

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$10

in two dimensions, and logarithmic transforms of the compressive and vortical potentials yield imaginary-time Schrödinger-type equations

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$11

together with a mixed density amplitude

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$12

satisfying a vector-potential-coupled equation

$\psi''(x)+[k^2-V(x)]\psi(x)=0,$13

(Beattie et al., 29 Apr 2026). This is not a Helmholtz–Schrödinger equation in the stationary scattering sense, but it shows that Helmholtz decomposition and Schrödinger-type evolution can be combined into an exact change of variables for a nonlinear dissipative system.

Across these usages, a stable conclusion is that the Helmholtz–Schrödinger equation is less a single model than a structural nexus. It joins fixed-frequency wave propagation, stationary quantum mechanics, inverse boundary value theory, special-function solvability, multigrid and domain decomposition, neural trial spaces built from plane waves, and quantum simulation frameworks that convert indefinite elliptic systems into Schrödinger-type dynamics (Lu, 2022, Dong et al., 2017, Gu et al., 31 Jul 2025).