Supersymmetric Expansion Algorithm

• Supersymmetric Expansion Algorithm is a systematic method in quantum mechanics that uses factorization and recursive expansions to derive spectra, eigenfunctions, and observables.
• It employs algebraic recurrences, logarithmic expansions, and quadrature formulations to bypass intermediate-state sums and ensure automatic normalization.
• Applied to potentials such as the Hulthén and anharmonic oscillator, SEA unifies analytic and numerical approaches for both exactly and non-exactly solvable models.

The Supersymmetric Expansion Algorithm (@@@@2@@@@) encompasses a distinct class of systematic expansions and algorithmic procedures arising in supersymmetric quantum mechanics, supersymmetric quantum field theory, and related domains of mathematical physics. Originating from attempts to generate analytical solutions and expansion schemes for quantum systems with non-trivial potentials, non-exactly solvable models, and complex algebraic structures, SEA exploits factorization properties, supersymmetric partner hierarchies, and recursive or algebraic expansions to produce explicit formulas for spectra, eigenfunctions, and observables, often circumventing intractable intermediate-state sums or inverting intractable operators. Over the past decade, several manifestations of SEA have been developed for both analytic and numerical applications in quantum mechanics, QFT, and mathematical computation, providing powerful tools for systems where traditional methods are cumbersome or inapplicable (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026, Zhou, 2023).

1. Supersymmetric Factorization and Edge-State Approach

A central theme of SEA is the representation of a Hamiltonian $H_0 = -\frac{d^2}{dx^2} + V_0(x, \lambda)$ via a factorization involving superpotentials and supersymmetric partner chains. For one-dimensional or radially reduced problems, SEA constructs a hierarchy of partner Hamiltonians: $$H_0 \longrightarrow H_1 \longrightarrow \dots \longrightarrow H_r \longrightarrow \cdots,$$ where each $H_r$ is related by intertwining differential operators

$$H_r = a_r^\dagger a_r + \varepsilon_r, \quad a_r = \frac{d}{dx} + W_r(x), \quad a_r^\dagger = -\frac{d}{dx} + W_r(x),$$

with $W_r(x)$ the superpotential. The "edge state" $u_r(x)$ of $H_r$ is the unique nodeless (usually ground) state obeying $a_r u_r(x) = 0$, i.e., $W_r(x) = -\frac{u_r'(x)}{u_r(x)}$. The Riccati equation for $W_r(x)$,

$W_r^2(x) - W_r'(x) = V_r(x) - \varepsilon_r,$

serves as the backbone of the expansion, with $W_{r+1}(x)$ and $V_{r+1}(x) = V_r(x) + 2 W_r'(x)$ generating the partner chain. Each excited state of $H_0$ is retrieved by applying the ladder of creation operators to the corresponding edge state of a higher partner: $$\phi_r(x) = a_0^\dagger a_1^\dagger \ldots a_{r-1}^\dagger\, u_r(x).$$ This structure ensures orthogonality and normalization throughout the spectrum without requiring an explicit Gram–Schmidt procedure (Napsuciale et al., 2024).

2. Logarithmic Expansion and Algebraic Recursion

When the potential $V_0(x, \lambda)$ possesses an analytic (often perturbative) dependence on a control parameter $\lambda$ (e.g., coupling or screening length), SEA systematically expands both the superpotential and energy: $$V_0(x, \lambda) = \sum_{k=0}^{\infty} v_{0k}(x) \lambda^k; \quad W_0(x, \lambda) = \sum_{k=0}^{\infty} w_{0k}(x) \lambda^k; \quad \varepsilon_0(\lambda) = \sum_{k=0}^{\infty} \varepsilon_{0k} \lambda^k.$$ Substituting these into the Riccati equation generates an infinite, hierarchical system of first-order ODEs for the $w_{0k}(x)$: $$\begin{aligned} O(\lambda^0):\quad & w_{00}^2 - w_{00}' = v_{00}(x) - \varepsilon_{00}, \ O(\lambda^k):\quad & 2 w_{00} w_{0k} - w_{0k}' = v_{0k}(x) - \varepsilon_{0k} - \sum_{m+n=k;\, m, n > 0} w_{0m} w_{0n}. \end{aligned}$$ For potentials where $v_{0k}(x)$ are (Laurent-)polynomials, one expands $w_{0k}(x)$ in a finite monomial basis, reducing the ODE to a finite, triangular algebraic system for each $k$: $$w_{0k}(x) = \sum_{\alpha} w_{0k,\alpha}\, x^{\alpha}.$$ The coefficients $w_{0k,\alpha}$ and $\varepsilon_{0k}$ are then solved algebraically order-by-order, enabling closed-form expressions for all expansion coefficients and providing analytic controllability to arbitrary order (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026).

3. Quadrature Formulation and Elimination of Intermediate-State Sums

A distinctive, practical advantage of SEA is the replacement of the infinite sums over intermediate eigenstates—ubiquitous in Rayleigh–Schrödinger perturbation theory—by univariate quadrature formulas. For each order $k$, the energy corrections and superpotential corrections for the edge state are given by: $$\varepsilon_{0k} = \int dx\, u_0(x,0)^2\; (v_{0k}(x) - B_{0k}(x)), \qquad w_{0k}(x) = u_0(x,0)^{-2} \int^{x} dt\, u_0(t,0)^2\; [\varepsilon_{0k} - v_{0k}(t) + B_{0k}(t)].$$ Here, $$B_{0k}(x) = \sum_{m+n=k;\, m, n > 0} w_{0m}(x)w_{0n}(x)$$ collects lower-order products. For excited states, the expansion is repeated for each partner $H_\ell$ and edge state $u_\ell(x,0)$, with corresponding formulas. For example, for the quartic oscillator, this yields the full Bender–Wu expansion for all levels, but without explicit reference to excited-state wavefunctions beyond the nodeless seeds (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026). This eliminates the computational bottleneck of crossings and sum convergence.

4. Algorithmic Steps and Implementation for General Potentials

The SEA, as a general method, can be outlined as follows:

1. Write $H_0$ in a dimensionless or scaled form and identify $V_0(x, \lambda)$, extracting $v_{0k}(x)$.
2. Obtain the $k=0$ (unperturbed) superpotential $w_{00}(x)$ and energy $\varepsilon_{00}$ from the Riccati equation.
3. For each order $k\geq 1$:
   • Compute the source term $B_{0k}(x)$ from lower-order $w_{0m}(x)$,
   • Convert the ODE for $w_{0k}(x)$ into an algebraic recurrence (via a monomial ansatz),
   • Solve for all $w_{0k,\alpha}$ and $\varepsilon_{0k}$, ensuring edge-state normalization.
4. Factorize $H_0$ and build $H_1$; repeat the expansion for $H_1$, $H_2$, ….
5. Construct all excited states of $H_0$ as $a_0^\dagger \cdots a_{r-1}^\dagger\, u_r(x)$, reconstructing the spectrum and eigenfunctions.

This approach is fully algebraic and recursive, amenable to algorithmic or symbolic computation, and applies universally to analytic, screened, and polynomial potentials (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026, Zhou, 2023).

5. Applications: Specific Potentials and Beyond

5.1 Hulthén and Screened Potentials

For the three-dimensional Hulthén potential, relevant in atomic, nuclear, and plasma problems: $$v_0(x, \lambda) = \frac{\ell(\ell + 1)}{x^2} - 2 \sum_{k=0}^\infty \frac{B_k^-}{k!}\, \lambda^k x^{k-1}$$ with $B_k^-$ the Bernoulli numbers, the SEA enables a closed-form, order-by-order analytic expansion for energy levels and radial wavefunctions for arbitrary quantum numbers, including nonzero angular momentum, previously inaccessible to perturbative methods. For $\ell=0$, the power series truncates and gives the known closed solution for $s$-states, while for $\ell>0$ it produces genuinely new analytic spectra and eigenfunctions (Napsuciale et al., 2024).

5.2 Anharmonic Oscillator

For polynomial cases such as the quartic oscillator,

$H_0 = -\frac{d^2}{dx^2} + x^2 + \lambda x^4,$

the SEA reconstructs the Bender–Wu expansion at arbitrary order for all levels. The edge-state algorithm sidesteps the need to handle nodal structures directly, with all nodality introduced solely via the intertwining operators. The resulting expansions can be resummed (e.g., Padé or Borel–Padé), yielding precision numerics at moderate coupling (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026).

5.3 Generalized Hypergeometric-like Operators

SEA has further been extended to completely general hypergeometric-like differential operators, including those with position-dependent kinetic terms and non-standard momentum maps. Active supersymmetrizations and shape-invariant recursions generate two distinct hierarchies of eigenfunctions and spectra: the principal and associated SUSYQM towers. These are determined recursively via algebraic relations on the coefficients of an expansion in the most suitable coordinate or factorization representation (Zhou, 2023).

6. Comparative Features and Algorithmic Advantages

Key technical and practical advantages of SEA include:

• Elimination of intermediate-state sums: All energy and wavefunction corrections are accessible via single or double quadratures against known edge-state densities, rather than needing infinite, slowly convergent sums over complete sets.
• Automatic orthogonality and normalization: All states are generated as ground states of properly intertwined SUSY partners, with the normalization and mutual orthogonality exact at each order.
• Applicability to both exactly and non-exactly solvable models: The method preserves full generality, treating analytic screening potentials and polynomially confined systems in a unified fashion, and admits extensions for potentials with shape-invariance or hidden symmetries.
• Closed-form recursive structure: All recursions are algebraic and suitable for symbolic or numerical automation to arbitrary perturbative order, or, in some potentials, to complete analytic solution (Napsuciale et al., 2024, Zhou, 2023).
• Flexibility for perturbation theory: SEA serves as the foundation for a supersymmetric reformulation of Rayleigh–Schrödinger perturbation theory, yielding compact quadrature-based formulas for corrections at all levels (Napsuciale et al., 14 Jan 2026).

7. Summary Table: Key Steps of the Supersymmetric Expansion Algorithm

Step	Description	Output
1. Factorization	$H_{r} = a_{r}^\dagger a_{r} + \varepsilon_r$ via $W_r(x)$	Riccati equation for $W_r(x)$
2. Logarithmic Expansion	Expand $V_0, W_0, \varepsilon_0$ in $\lambda$	Cascading ODEs for $w_{0k}(x)$
3. Algebraic Recurrences	Expand $w_{0k} = \sum_\alpha w_{0k,\alpha} x^\alpha$	Linear system for $w_{0k,\alpha}$, $\varepsilon_{0k}$
4. Construct Edge States	Integrate $W_{r}(x, \lambda)$ to get $u_r(x, \lambda)$	Nodeless eigenfunction
5. Build Spectrum	Apply $a_{0}^\dagger \cdots a_{r-1}^\dagger$ to $u_r(x)$	All states of $H_0$

SEA thus equips the supersymmetric, quantum, and mathematical physics communities with a robust analytic and algorithmic toolset for spectral problems that had remained recalcitrant to standard methods, enabling both exhaustive analytic control and practical numerical outputs (Napsuciale et al., 2024, Napsuciale et al., 14 Jan 2026, Zhou, 2023).