Sturm–Liouville: Impedance Form

• The impedance-form Sturm–Liouville equation is a second-order differential operator defined via a quasi-derivative, extending both self-adjoint and non-self-adjoint problem frameworks.
• It employs the Liouville transform to convert the operator to a Schrödinger form, preserving spectral data and enabling explicit inverse spectral analysis.
• Advanced computational techniques such as the SPPS and NSBF methods yield high-accuracy eigenvalue estimates and efficient numerical approximations.

The Sturm–Liouville equation in impedance form is a second-order linear differential equation central to modern spectral theory, functional analysis, and mathematical physics. Distinct from the classical form by the presence of the impedance (or weight) function acting through a quasi-derivative, the impedance formulation captures a wider class of self-adjoint and non-self-adjoint problems, accommodates distributional potentials, facilitates explicit transforms to Schrödinger operators, and serves as a robust foundation for analytic, spectral, and computational advances.

1. Differential Operators and Quasi-derivative Formulation

The general Sturm–Liouville differential expression admits a quasi-derivative framework: $$\tau f(x) := \frac{1}{r(x)} \left( -\frac{d}{dx}\big( p(x)[f'(x) + s(x)f(x)] \big) + s(x)p(x)[f'(x) + s(x)f(x)] + q(x)f(x) \right)$$ where $p(x)$, $q(x)$, $r(x)$, and $s(x)$ are measurable, real-valued coefficients, $p(x) \neq 0$, $r(x) > 0$ a.e., $p^{-1}, q, r, s \in L^1_{\mathrm{loc}}$, and $f$ is locally absolutely continuous with $p(f' + s f) \in AC_{\mathrm{loc}}$ (Eckhardt et al., 2012).

The impedance form is derived by specializing:

• $q(x) = 0$,
• $s(x) = 0$,
• $r(x) = p(x) =: \kappa(x) > 0$,

yielding

$-(\kappa(x) u'(x))' = \lambda \kappa(x) u(x),$

which constitutes the hallmark impedance-form Sturm–Liouville equation and is frequently normalized so that $\kappa(0) = 1$. More generally, formulations with quasi-derivative $D_s[y](x) = y'(x) + s(x)y(x)$ and

$$-\frac{d}{dx}[p_1(x) D_s[y](x)] + p_0(x)y(x) = \lambda w(x) y(x),$$

cover a broader operator class, including non-self-adjoint cases (Bandyopadhyay et al., 17 Aug 2025).

2. Boundary Conditions and Self-Adjoint Realization

Sturm–Liouville equations in impedance form admit separated boundary conditions: $$\alpha_1 y(a) + \alpha_2 D_s[y](a) = 0, \qquad \beta_1 y(b) + \beta_2 D_s[y](b) = 0,$$ subject to $\alpha_1^2 + \alpha_2^2 > 0$, $\beta_1^2 + \beta_2^2 > 0$ (Bandyopadhyay et al., 17 Aug 2025). Classical Robin (impedance-type) conditions, such as $p(a) u'(a) + Z_a u(a) = 0$, arise both in time-independent and evolutionary (PDE) settings, modeling finite-mass or dissipative boundaries (Picard et al., 2012).

In the Hilbert space $L^2((a,b); r(x)\,dx)$, the operator is self-adjoint on the domain of functions in $H^2$ satisfying the boundary conditions. The spectrum is real and discrete on bounded intervals and incorporates the impedance boundary contributions as mass-like terms in the energy form: $$\langle Au \mid u \rangle_{L^2(r\,dx)} = \int_a^b \kappa(x) |u'(x)|^2\,dx + Z_b |u(b)|^2 + Z_a |u(a)|^2,$$ with regularizing effects on both spectral and evolutionary problems (Picard et al., 2012).

3. Spectral Theory and Transformations

Liouville Transformation and Equivalence to Schrödinger Form

The impedance operator is unitarily equivalent to a Schrödinger operator via the Liouville transformation. Setting $q(x) = \kappa'(x)/\kappa(x)$ and $f(x) = y(x)/\kappa(x)^{1/2}$, the change of variable

$v(x) = \exp\left(\int_a^x s(t)\,dt\right) y(x)$

converts the quasi-derivative problem to

$$-\bigl(p_1(x) v'(x)\bigr)' + Q(x) v(x) = \lambda w(x) v(x)$$

with $Q(x) = p_0(x) - (p_1(x)s(x))' + p_1(x) s^2(x)$. For the canonical impedance equation,

$-\kappa^{-2} (\kappa^2 u')' = \lambda u,$

the Liouville transform $y(x) = \kappa(x) u(x)$ yields the Schrödinger equation

$$-y''(x) + p(x) y(x) = \lambda y(x), \quad p(x) = q'(x) + q(x)^2,$$

with preservation of spectral data and boundary conditions (Isozaki et al., 2013, Bandyopadhyay et al., 17 Aug 2025). The transformation establishes a real-analytic isomorphism between the Banach manifolds of impedance functions (with $q \in H^1$) and Schrödinger potentials (with $p \in L^2$), ensuring bijective correspondence of inverse problems (Isozaki et al., 2013).

Fundamental Solutions, Weyl–Titchmarsh Theory, and m-Function

With separated boundary conditions, spectral theory proceeds via construction of a real-entire fundamental system $(\varphi_z, \theta_z)$ for $(\tau - z)u = 0$, with the Weyl–Titchmarsh solution and associated m-function

$$m(z) = -\frac{\varphi_z[1](b)\cos\beta - \varphi_z(b)\sin\beta}{\theta_z[1](b)\cos\beta - \theta_z(b)\sin\beta}$$

serving as the analytic parameterization of spectral data. The spectral measure is generated via the Herglotz representation for $m(z)$ and yields a unitary spectral transform for the impedance operator (Eckhardt et al., 2012).

4. Analytic and Algebraic Tools: Prüfer Transformation, Formal Powers, and Transmutation

Generalized Prüfer Substitution and Oscillation Theory

The generalized Prüfer transformation, suitable for impedance/quasi-derivative operators, is defined via

$$y(x) = \rho(x) \cos\theta(x), \qquad D_s[y](x) = \rho(x) \sin\theta(x),$$

leading to a nonlinear amplitude-phase system, e.g.,

$$\theta'(x) = -\frac{\cos^2\theta(x)}{p_1(x)} - s(x) \sin\theta(x)\cos\theta(x) + p_0(x) \sin^2\theta(x) - \frac{s(x)}{p_1(x)} \rho(x) \sin^2\theta(x)\cos\theta(x)$$

for the phase variable (Bandyopadhyay et al., 17 Aug 2025).

Major consequences:

• Comparison theorems: ordering of potential coefficients translates to phase monotonicity and classical Sturm oscillation structure.
• Monotonicity in the spectral parameter: the zeros of eigenfunctions shift strictly with $\lambda$.
• Variational bounds: sharp n-th eigenvalue estimates using minimax techniques.
• Phase-function eigenvalue characterization: root-finding for $\theta(b;\lambda) = n\pi$ replaces two-point boundary value solvers (Bandyopadhyay et al., 17 Aug 2025).

Formal Powers, Transmutation Operators, and SPPS Method

The spectral parameter power series (SPPS) method constructs “formal powers” $\varphi_a^{(k)}(x)$ recursively, forming a complete basis in $L^p(I; a^2)$ and even $W^{m,p}(I)$ provided $a, 1/a \in L^2$ and sufficient regularity. The method produces the “transmutation operator”

$$(T_a[u])(x) = u(x) - \int_{x_0}^x K_a(x,t) u'(t) dt,$$

K_a(x, t) continuous and satisfying the Goursat boundary data $K_a(x, x-x_0) = 1-1/a(x), K_a(x, -[x-x_0])=0$, such that $T_a[x^n] = \varphi_a^{(n)}(x)$ and $T_a$ intertwines the impedance and Laplacian operators (Vicente-Benítez, 19 Jun 2025). In weakened regularity, the bounded operator property persists, allowing Volterra-type representations and completeness in $L^2$ (Vicente-Benítez, 19 Jun 2025, Vicente-Benítez, 15 Nov 2025).

Neumann Series of Spherical Bessel Functions (NSBF)

For $\kappa \in W^{1,2}$ with no zeros, solutions admit uniform convergent NSBF representations: $$e_\kappa(\rho,x) = e^{i \rho x} - i \rho \sum_{n=0}^\infty i^n \alpha_n(x) j_n(\rho x),$$ with $\alpha_n(x)$ given via recursive single-integral quadratures and $j_n$ denoting the spherical Bessel functions (Márquez-Hernández et al., 8 Jan 2026). Series for Dirichlet or impedance boundary value solutions systematically translate spectral problems into root-finding on explicitly constructed analytic functions, with truncation errors decaying as $O(N^{-(p+1/2)})$ for $\kappa \in C^{p+1}$ (Márquez-Hernández et al., 8 Jan 2026).

5. Inverse Problems and Uniqueness Results

The analytic isomorphism induced by the Liouville transformation ensures that the inverse spectral problem for the impedance Sturm–Liouville operator is globally well-posed and equivalent to the established theory for Schrödinger operators:

• The mapping from the impedance function (parametrized by $q \in H^1$) and boundary parameters to spectral data (eigenvalues and norming constants) is real-analytic and bijective (Isozaki et al., 2013).
• The inverse procedure reconstructs $q$, and thus $\kappa(x)$, from spectral data by pullback through the analytic inverse of the Liouville map.
• Uniqueness is guaranteed up to shift and scaling by two- or three-spectral data: if two sets of impedance data $(p, r)$ produce the same spectra for Dirichlet/Neumann (or internal subinterval) conditions, they are equivalent under explicit rescaling and translation (Eckhardt et al., 2012).

Local uniqueness holds via local Borg–Marchenko-type conditions on the Weyl $m$-functions, providing recovery of the impedance on subintervals from spectral asymptotics (Eckhardt et al., 2012).

6. Darboux Transformation and Operator Intertwining

The operator-level Darboux transformation applies to impedance-form Sturm–Liouville equations, paralleling the Schrödinger case. The first-order operator intertwines original and Darboux-transformed equations: $$D_a \circ \mathbf{L}_a = \mathbf{L}_{1/a} \circ D_a,$$ where $D_a[u] = a^2 u'$ and $\mathbf{L}_{1/a}$ is the impedance operator with reciprocal function. The transmutation operator kernels for $a$ and $1/a$ are related by explicit integral identities involving both $a^2$ and derivatives in $x, t$ (Vicente-Benítez, 19 Jun 2025, Vicente-Benítez, 15 Nov 2025). These relations enable the inheritance of structural and spectral properties across Darboux partners.

7. Computational and Evolutionary Aspects

Impedance-form Sturm–Liouville equations appear naturally in evolutionary (PDE) contexts, such as hyperbolic equations with impedance-type boundary conditions: $$r(x)\partial_t^2 u + \eta(x) \partial_t u + q(x)u - \partial_x(p(x) \partial_x u) = f(t,x),$$ with boundary conditions $p(a) \partial_x u(t,a) + Z_a u(t,a) = 0$ (Picard et al., 2012). The structure ensures existence and regularity (via semigroup theory) and preserves self-adjointness or maximal dissipativity required for well-posedness.

Contemporary numerical spectral methods for impedance-form problems leverage the SPPS and NSBF expansions for high-accuracy eigenvalue computations, with uniform error control that does not deteriorate for high spectral indices, thus supporting the computational tractability of large-scale or high-frequency problems (Márquez-Hernández et al., 8 Jan 2026). The completeness of formal powers ensures dense approximation property in $L^2$ and Sobolev spaces, fundamental for efficient solution approximation (Vicente-Benítez, 19 Jun 2025).

---

The impedance form encapsulates both structural generality and computational accessibility, unifies direct, spectral, and inverse results (including under weak regularity or distributional data), and enables effective algebraic, analytic, and numerical methodologies foundational for modern operator theory and applications in mathematical physics.