---
title: 'Sturm–Liouville: Impedance Form'
url: https://www.emergentmind.com/topics/sturm-liouville-equation-in-impedance-form
type: topic
---

# Sturm–Liouville: Impedance Form

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}}$ [1210.7628].

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 [2508.12369].

## 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$ [2508.12369]. 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 [1204.5458].

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 [1204.5458].

## 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 [1307.1924, 2508.12369]. 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 [1307.1924].

### 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 [1210.7628].

## 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 [2508.12369].

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 [2508.12369].

### 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 [2506.16604]. In weakened regularity, the bounded operator property persists, allowing Volterra-type representations and completeness in $L^2$ [2506.16604, 2511.12094].

### 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 [2601.04513]. 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}$ [2601.04513].

## 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 [1307.1924].
- 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 [1210.7628].

Local uniqueness holds via local Borg–Marchenko-type conditions on the Weyl $m$-functions, providing recovery of the impedance on subintervals from spectral asymptotics [1210.7628].

## 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$ [2506.16604, 2511.12094]. 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$ [1204.5458]. 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 [2601.04513]. The completeness of formal powers ensures dense approximation property in $L^2$ and Sobolev spaces, fundamental for efficient solution approximation [2506.16604].

---

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.

Source: https://www.emergentmind.com/topics/sturm-liouville-equation-in-impedance-form