Sturm-Liouville Hierarchy in Integrable Systems

• Sturm-Liouville hierarchy is a framework derived from isospectral deformations, unifying nonlinear evolution equations like KdV and Camassa-Holm.
• It employs zero-curvature conditions, Weyl function asymptotics, and Darboux transformations to construct varied spectral and fractional operator classes.
• Its generalizations extend to q-fractional, multi-parameter, and higher-order ladder problems, yielding novel exceptional orthogonal polynomials and integrable models.

The collected literature suggests that the expression “Sturm-Liouville hierarchy” denotes several mathematically distinct but related constructions built from Sturm-Liouville spectral problems. In one sense, it is a hierarchy of isospectral nonlinear evolution equations generated from the generalized spectral problem

$$-(p\varphi')' + q\varphi = \lambda y\varphi,\qquad y>0,$$

or, in the reduced form,

$-\varphi'' + q\varphi = \lambda y\varphi,$

with the Korteweg-de Vries and Camassa-Holm hierarchies appearing as special cases. In other senses, the same phrase refers to q-fractional and fractional generalizations, Darboux-generated hierarchies of rational Sturm-Liouville problems, higher-order shape-invariant families tied to Painlevé IV, and multi-parameter spectral constructions. Across these settings, the unifying theme is the persistence—often in modified form—of the core Sturm-Liouville structures of self-adjointness, orthogonality, variational characterization, spectral discreteness, and transform-based generation of solvable families (Johnson et al., 2014, Rubbioni et al., 24 Sep 2025, Turemuratova et al., 28 Apr 2025, Mansour, 2016, Natanson, 2013).

1. Zero-curvature and spectral formulations

A central usage of the term arises in the theory of integrable evolution equations. The hierarchy is built around the generalized Sturm-Liouville operator

$L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$

acting on the weighted Hilbert space $L^2(\mathbb R, y(x)\,dx)$, with $p(x), q(x), y(x)$ uniformly continuous and bounded, and $p,y>0$. The corresponding evolution equations are formulated through zero-curvature relations

$A_t - B_{r,x} + [A,B_r]=0,$

with $A$ constructed from the spectral problem and $B_r$ built from functions such as $U_g(x,\lambda)$, $-\varphi'' + q\varphi = \lambda y\varphi,$0, and $-\varphi'' + q\varphi = \lambda y\varphi,$1. In this framework, the hierarchy includes the Korteweg-de Vries hierarchy when $-\varphi'' + q\varphi = \lambda y\varphi,$2, and the Camassa-Holm hierarchy under different reductions and via Liouville transformations; the Hunter-Saxton equation also appears among the incorporated degenerate equations (Johnson et al., 2014).

A later formulation defines the hierarchy again through zero-curvature, but simultaneously through the asymptotic properties of the Weyl $-\varphi'' + q\varphi = \lambda y\varphi,$3-functions for

$-\varphi'' + q\varphi = \lambda y\varphi,$4

in $-\varphi'' + q\varphi = \lambda y\varphi,$5. For the half-line problems, the Weyl functions

$-\varphi'' + q\varphi = \lambda y\varphi,$6

satisfy the Riccati equation

$-\varphi'' + q\varphi = \lambda y\varphi,$7

and the diagonal Green’s function may be expressed as

$-\varphi'' + q\varphi = \lambda y\varphi,$8

The zero-curvature condition yields isospectral flows, and the evolution of the Weyl functions encodes the dynamics of the hierarchy. In this version, setting $-\varphi'' + q\varphi = \lambda y\varphi,$9, $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$0, $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$1 yields the classical KdV equation, while setting $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$2, $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$3, $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$4 gives the Camassa-Holm equation (Rubbioni et al., 24 Sep 2025).

The two formulations are compatible at the level of principle: both present the hierarchy as an isospectral deformation theory of Sturm-Liouville operators. This suggests that the hierarchy is best understood not as a single PDE sequence in isolation, but as a spectral mechanism in which zero-curvature, Riccati equations, and Weyl data are interchangeable descriptions of the same integrable structure.

2. Algebro-geometric and reflectionless sectors

For finite-gap spectra, explicit solutions are constructed through hyperelliptic Riemann surfaces. The motion of poles $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$5 associated to eigenvalues in the spectral gaps is governed by ODEs derived from the Riccati equation and the spectral data, and the potentials are reconstructed by trace formulas such as

$L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$6

and

$L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$7

The resulting solutions are algebro-geometric, almost periodic or quasi-periodic in $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$8 and $L_a = \frac{1}{y}\{-DpD+q\},\qquad D=\frac{d}{dx},$9, and preserve the spectrum under the hierarchy flows (Johnson et al., 2014).

The finite-gap construction extends to infinitely many gaps by taking limits under summability conditions on the gap endpoints. In that limit, one obtains reflectionless Sturm-Liouville potentials whose spectra may be a countable union of intervals, with accumulation at a finite point or at infinity. A defining property stated for these reflectionless potentials is

$L^2(\mathbb R, y(x)\,dx)$0

The limit procedure uses Weierstrass-Runge convergence factors for infinite products, and the trace formulas and pole-motion equations remain valid in the infinite-gap setting. A specific consequence emphasized in the literature is the construction of solutions to the Camassa-Holm hierarchy with reflectionless initial data corresponding to finite accumulation points in the spectrum (Johnson et al., 2014).

In the Weyl-function formulation, the same isospectral sector is described by the evolution of

$L^2(\mathbb R, y(x)\,dx)$1

which satisfies

$L^2(\mathbb R, y(x)\,dx)$2

The asymptotic Laurent expansion of $L^2(\mathbb R, y(x)\,dx)$3 at infinity provides recursive coefficients determined by $L^2(\mathbb R, y(x)\,dx)$4 and $L^2(\mathbb R, y(x)\,dx)$5, and the compatibility conditions extracted from this expansion generate the higher equations of the hierarchy. The literature states that full Laurent expansions correspond to algebro-geometric potentials and their limits, while truncated expansions correspond to other classes, including decaying and scattering potentials (Rubbioni et al., 24 Sep 2025).

3. Fractional, q-fractional, and graph-based extensions

A different but structurally related use of hierarchy language appears in fractional and q-fractional generalizations of Sturm-Liouville theory. For the regular q-fractional Sturm-Liouville problem, the operator

$L^2(\mathbb R, y(x)\,dx)$6

acts on a q-linear grid, with $L^2(\mathbb R, y(x)\,dx)$7, $L^2(\mathbb R, y(x)\,dx)$8, and $L^2(\mathbb R, y(x)\,dx)$9. The corresponding eigenvalue problem is

$p(x), q(x), y(x)$0

subject to boundary conditions involving $p(x), q(x), y(x)$1, $p(x), q(x), y(x)$2, and the right-sided q-fractional integral of $p(x), q(x), y(x)$3. The operator is stated to be self-adjoint in a q-weighted $p(x), q(x), y(x)$4-type Hilbert space; all eigenvalues are real; eigenfunctions corresponding to distinct eigenvalues are orthogonal; the q-fractional Wronskian is constant on the interval; and the geometric multiplicity of each eigenvalue is $p(x), q(x), y(x)$5. An explicit example has eigenfunctions $p(x), q(x), y(x)$6, where $p(x), q(x), y(x)$7 is the little q-Jacobi polynomial, with eigenvalues

$p(x), q(x), y(x)$8

The variational development for the q-fractional problem shows, for $p(x), q(x), y(x)$9 and Dirichlet conditions $p,y>0$0, the existence of a countable set of real eigenvalues and associated orthogonal eigenfunctions, together with a Rayleigh quotient whose minimum is the first eigenvalue (Mansour, 2016, Mansour, 2016).

The literature explicitly places this q-fractional problem in a hierarchy: it encompasses the classical regular Sturm-Liouville problem when $p,y>0$1 and $p,y>0$2, the regular q-Sturm-Liouville problem when $p,y>0$3, and the fractional Sturm-Liouville problem when $p,y>0$4 and $p,y>0$5. This is a precise sense in which “hierarchy” means a chain of increasingly general operator classes rather than a hierarchy of commuting flows (Mansour, 2016).

The fractional problem on metric graphs extends the same theme from intervals or q-grids to network domains. For a finite, connected metric graph $p,y>0$6 with edges $p,y>0$7, the operator is defined edgewise by

$p,y>0$8

using a left Riemann-Liouville fractional derivative and a right Caputo fractional derivative. Boundary vertices satisfy

$p,y>0$9

internal vertices satisfy continuity of $A_t - B_{r,x} + [A,B_r]=0,$0, and the transmission condition is

$A_t - B_{r,x} + [A,B_r]=0,$1

The main spectral theorem gives a discrete, real, non-negative spectrum

$A_t - B_{r,x} + [A,B_r]=0,$2

and a complete orthonormal system of eigenfunctions in $A_t - B_{r,x} + [A,B_r]=0,$3. The eigenvalues admit a min-max characterization, the quadratic form is non-negative, and

$A_t - B_{r,x} + [A,B_r]=0,$4

The paper also states the estimate

$A_t - B_{r,x} + [A,B_r]=0,$5

and identifies the classical Sturm-Liouville theory as the special case $A_t - B_{r,x} + [A,B_r]=0,$6 on a single interval (Turemuratova et al., 28 Apr 2025).

4. Darboux-generated rational hierarchies and exceptional orthogonal polynomials

In another established sense, a Sturm-Liouville hierarchy is produced by repeated Darboux transformations of rational Sturm-Liouville equations. The canonical starting point is a second-order Sturm-Liouville equation written in Schrödinger form,

$A_t - B_{r,x} + [A,B_r]=0,$7

where $A_t - B_{r,x} + [A,B_r]=0,$8 is the Bose invariant. Given a factorization function $A_t - B_{r,x} + [A,B_r]=0,$9 at $A$0, the first-order intertwiner

$A$1

generates a transformed problem whose Bose invariant differs by the logarithmic second derivative of $A$2. Recursive application of this canonical Darboux transformation produces hierarchies of rational extensions (Natanson, 2013).

The decisive structural condition in this setting is the presence of energy-independent characteristic exponents at the singular endpoints. The Darboux-Pöschl-Teller potential and the isotonic oscillator are presented as paradigmatic examples. When the exponents are energy-independent, the polynomial solution space is preserved under the transformations, and the transformed hierarchy remains within the domain of exceptional orthogonal polynomials. For the Darboux-Pöschl-Teller case, the partner Sturm-Liouville problem is solved by exceptional Jacobi polynomials; for the isotonic oscillator, the corresponding exceptional family is the exceptional Laguerre system. The literature states that, in this case, the exceptional sequences are infinite, whereas energy-dependent exponents permit only finite sets of exceptional polynomials (Natanson, 2013).

This framework also organizes a classification by seed type. The factorization functions are almost-everywhere holomorphic solutions, typically nodeless and lying below the ground state energy. The transformed potentials remain rational, but may acquire exceptional singularities. The paper further states that the trigonometric and hyperbolic Pöschl-Teller systems are unified under linear-fractional transformations and that their rational extensions generate the same $A$3 sets of GS Heine polynomials. It also identifies agreement with rational extensions previously constructed by Quesne, Odake & Sasaki, and with the theory of exceptional orthogonal polynomials due to Gómez-Ullate, Kamran, and Milson (Natanson, 2013).

A common misconception is that exceptional orthogonal polynomial families are merely isolated deformations of classical ones. In the Darboux framework summarized here, the relevant papers instead present them as systematically organized outputs of a rational Sturm-Liouville hierarchy generated by single-step and multi-step transformations.

5. Higher-order ladder hierarchies, Painlevé IV, and three-band decompositions

A further generalization constructs Sturm-Liouville problems from Hamiltonians satisfying a third-order shape-invariance condition. The Hamiltonian

$A$4

obeys

$A$5

with $A$6 and $A$7 third-order differential ladder operators. The potential is a rationally extended oscillator expressed through generalized Okamoto polynomials and rational solutions of the fourth Painlevé equation in the “$A$8” hierarchy:

$A$9

This produces a Sturm-Liouville spectral problem whose algebraic structure is controlled by Painlevé IV data (Hussin et al., 2021).

Because the annihilation operator is third order, the construction yields exactly three zero-modes,

$B_r$0

and the Hilbert space decomposes as

$B_r$1

The three spectral branches are

$B_r$2

and no eigenvalues from different branches overlap. Higher modes are generated by $B_r$3, while each eigenfunction has the form

$B_r$4

with

$B_r$5

The paper emphasizes that each branch satisfies its own three-term recurrence relation, with one initial condition per sequence (Hussin et al., 2021).

The significance of this result is not merely spectral multiplicity. The paper states that it provides, for the first time, a three-term recurrence relation for a special family of exceptional Hermite polynomials associated with “double partitions” $B_r$6, in contrast with much higher-order recurrences previously encountered in the literature. This suggests a new tier of Sturm-Liouville hierarchy in which higher-order ladder operators replace the classical single-chain orthogonal-polynomial picture by a finite direct sum of disjoint recurrence sectors (Hussin et al., 2021).

6. Multi-parameter, operator-pencil, and vessel-based generalizations

The hierarchy theme also appears when the spectral side itself is enlarged to several parameters. For

$B_r$7

explicit formulas are available for two linearly independent solutions as power series in $B_r$8. If $B_r$9 is a nonvanishing solution of the homogeneous equation, the solutions are

$U_g(x,\lambda)$0

and

$U_g(x,\lambda)$1

where the coefficients are recursively generated by indefinite integrals. These formulas reduce to the known single-parameter SPPS representation when $U_g(x,\lambda)$2, and they extend further to equations of the form

$U_g(x,\lambda)$3

The same paper notes that Sturm-Liouville pencils, or polynomial-type eigenvalue problems, fit into this framework directly (Porter, 2015).

In the language of explicit solution generation for nonlinear equations, the vessel formalism produces another hierarchy associated with the Sturm-Liouville equation

$U_g(x,\lambda)$4

The KdV vessel

$U_g(x,\lambda)$5

encodes both the spectral problem and the time evolution. The hierarchy is expressed for

$U_g(x,\lambda)$6

through

$U_g(x,\lambda)$7

with

$U_g(x,\lambda)$8

The tau function

$U_g(x,\lambda)$9

recovers the potential via

$-\varphi'' + q\varphi = \lambda y\varphi,$00

and explicit multisoliton solutions are obtained from exponential finite-dimensional realizations of $-\varphi'' + q\varphi = \lambda y\varphi,$01, $-\varphi'' + q\varphi = \lambda y\varphi,$02, and $-\varphi'' + q\varphi = \lambda y\varphi,$03 (Melnikov, 2012).

Taken together, these constructions show that “hierarchy” in the Sturm-Liouville setting can refer to at least three rigorous organizations: a sequence of commuting nonlinear flows, a nested sequence of generalized operator classes, and a recursively generated family of transformed or multi-parameter spectral problems. The literature does not collapse these meanings into one definition. A plausible implication is that the term functions less as a single taxonomy than as a recurrent structural motif: once the Sturm-Liouville problem is placed in a sufficiently rich algebraic, fractional, or spectral-analytic framework, hierarchies emerge naturally.