⚛️ high-energy theory

Stationary perturbation theory without sums over intermediate states: Supersymmetric Expansion Algorithm

This paper demonstrates that the supersymmetric expansion algorithm can efficiently derive Rayleigh-Schrödinger perturbation theory results for energy and eigenstate corrections without summing over intermediate states, instead expressing them as integrals weighted by the probability densities of edge states.

Original authors: M. Napsuciale, S. Rodríguez

Published 2026-01-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Napsuciale, S. Rodríguez

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a musical instrument will sound if you slightly change the tension of its strings. In the world of quantum physics, this is like calculating how the energy of an atom changes when you add a tiny bit of extra force (a "perturbation").

For a hundred years, physicists have used a standard method called Rayleigh-Schrödinger perturbation theory to do this. Think of this old method like trying to calculate the total weight of a backpack by adding up the weight of every single grain of sand inside it, one by one. It works, but it's messy. To get the answer, you have to sum up an infinite number of "intermediate states" (all the possible ways the system could wiggle in between). As you try to get more precise, the list of things you have to add gets longer and longer, making the math incredibly awkward and difficult to handle.

The New Approach: The "Supersymmetric Expansion Algorithm" (SEA)

The authors of this paper, M. Napsuciale and S. Rodríguez, propose a new way to solve this problem. They call it the Supersymmetric Expansion Algorithm. Instead of adding up an infinite list of grains of sand, they show you how to measure the backpack directly using a single, smooth calculation.

Here is how their method works, broken down into simple concepts:

1. The "Edge" of the Problem

In quantum mechanics, some states (like the ground state of an atom) are smooth and have no "kinks" or "nodes" (places where the wave goes flat). The authors realized that if you can solve the problem for these smooth, "nodeless" states first, you can use a special mathematical trick called supersymmetry to build the solutions for all the other, more complex states (the ones with kinks) from them.

Think of it like building a house. Instead of trying to build every room at once, you first build a perfect, solid foundation (the "edge state"). Once that foundation is solid, you can easily construct the rest of the house on top of it.

2. Turning a Sum into a Smooth Slide

The biggest breakthrough in this paper is how they handle the math.

The Old Way: To find the correction to the energy, you had to perform a "sum over intermediate states." Imagine trying to walk up a staircase where every step is a different, unknown height. You have to calculate the height of every single step before you can reach the top.
The New Way (SEA): The authors show that you can turn this staircase into a smooth slide. Instead of counting steps, you just calculate the area under a curve (an integral). In math terms, they reduce the problem to "quadrature forms."

This means the answer comes out as a single, clean integral (a type of area calculation) weighted by the probability of where the particle is likely to be. It's like measuring the total volume of water in a pool by looking at the shape of the pool, rather than counting every single drop of water.

3. The "Logarithmic" Shortcut

To get this smooth slide, the authors use a clever trick involving the "logarithmic form" of the Schrödinger equation.

Imagine the wave function (the description of the particle) is a complex, tangled rope.
The authors take the "logarithm" of this rope, which untangles it into a simpler shape called a superpotential.
They then expand this superpotential in a series, solving a cascade of simple, linear equations one after another. It's like peeling an onion layer by layer, where each layer is easy to handle once the previous one is removed.

4. What They Actually Did

The paper claims to do three main things:

Generalize a previous method: They took a method that worked well for simple, one-dimensional systems (where only the ground state is smooth) and expanded it to work for all states, including excited states (which have kinks/nodes), in any dimension.
Avoid the "Sum": They proved that you never need to sum over intermediate states again. You only need to perform integrals (area calculations) based on the probability densities of the "edge states."
Handle Complex Potentials: They showed this works even when the math isn't simple (polynomials). They tested their method on a "particle in a box" with a harmonic oscillator added to it. They successfully calculated the energy corrections up to the third order, getting the same results as the old, messy method but with a much cleaner process.

The Bottom Line

The authors aren't claiming to discover new particles or change how we build computers. They are offering a better calculator.

If the old method of quantum perturbation theory is like trying to solve a puzzle by gluing together thousands of tiny, jagged pieces, the Supersymmetric Expansion Algorithm is like having a template that lets you trace the whole picture in one smooth, continuous motion. It makes calculating energy corrections for quantum systems faster, cleaner, and avoids the "awkward sums" that have plagued physicists for a century.

Technical Summary: Stationary Perturbation Theory without Sums over Intermediate States: Supersymmetric Expansion Algorithm

Problem Statement
The paper addresses a fundamental limitation in conventional Rayleigh-Schrödinger perturbation theory (RSPT). While RSPT provides a systematic framework for calculating corrections to energy eigenvalues and eigenstates, it relies on sums over intermediate states. The first-order correction to eigenstates already requires an "awkward sum," and higher-order corrections necessitate increasingly complex sums over intermediate states at every order. This formalism becomes computationally cumbersome and often impractical beyond the lowest orders. Previous alternatives, such as sum rules or the logarithmic formulation of the Schrödinger equation proposed by Aharonov and Au, have successfully reduced ground-state calculations to quadrature forms (integrals) for one-dimensional systems but remain less convenient for excited states due to the presence of nodes.

Methodology: The Supersymmetric Expansion Algorithm (SEA)
The authors propose the Supersymmetric Expansion Algorithm (SEA) to obtain perturbative corrections to arbitrary order for all eigenstates (ground and excited) without summing over intermediate states. The method exploits the hidden supersymmetry inherent in non-relativistic quantum mechanical problems.

The core procedure involves the following steps:

Logarithmic Formulation: The Schrödinger equation is cast into a non-linear logarithmic form involving a superpotential $W(x, \lambda) = -\frac{d}{dx} \ln u(x, \lambda)$ , where $u$ is the wavefunction.
Perturbative Expansion: The superpotential, the potential, and the energy are expanded in powers of a small parameter $\lambda$ . This transforms the non-linear differential equation into a hierarchical cascade of coupled first-order linear differential equations for the expansion coefficients.
Nodeless "Edge" States: The algorithm first solves for the nodeless eigenstate (the "edge state") of the Hamiltonian at a given supersymmetric level. For one-dimensional systems, this is the ground state; for higher-dimensional systems reducible to effective one-dimensional potentials (like central potentials), these are the states with the highest angular momentum $l$ for a given principal quantum number $n$ .
Recursive Solution via Supersymmetry:
- The nodeless solution for the $r$ -th supersymmetric partner Hamiltonian $H_r$ is solved using the logarithmic expansion.
- The leading order term for the nodeless solution of $H_r$ is derived from the known unperturbed solutions of the original Hamiltonian $H_0$ by recursively applying the inverse of the supersymmetric raising operators ( $a^\dagger$ ).
- Once the nodeless solution for $H_r$ is established, the perturbative corrections for its energy and superpotential are calculated using integrals weighted by the unperturbed probability density of that specific nodeless state.
- The eigenstates of the original Hamiltonian $H_0$ with $r$ nodes are then reconstructed by applying the chain of raising operators ( $a^\dagger_0, a^\dagger_1, \dots$ ) to the nodeless solution of $H_r$ .

Key Contributions and Results

Elimination of Intermediate State Sums: The primary contribution is the reduction of perturbative corrections to "quadrature forms." Energy corrections and wavefunction corrections are expressed as integrals weighted by the probability densities of the relevant edge states, completely avoiding the summation over intermediate states required in RSPT.
Arbitrary Order and General Applicability: The formalism is applicable to arbitrary orders of perturbation theory and extends to all eigenstates, including excited states with nodes. It works for both polynomial and non-polynomial potentials.
Analytical Solutions for Specific Potentials: The authors demonstrate the method's efficacy by recovering known analytical solutions for the Yukawa, Hulthén, anharmonic, and Cornell potentials, where the superpotential coefficients are polynomials.
Non-Polynomial Example: The paper applies the method to a particle in a box with a harmonic oscillator perturbation. This case involves non-polynomial superpotentials. The authors successfully derive the energy corrections up to order $\lambda^3$ for both the ground state and the first excited state, providing explicit analytical expressions for the energy eigenvalues.

Significance and Claims
The paper claims that the SEA provides a "complete analytical solution" to perturbation problems in the form of a power series in the expansion parameter. By leveraging the hidden supersymmetry of quantum systems, the method offers a transparent and systematic alternative to RSPT. The authors emphasize that while previous methods (like Aharonov and Au's) were limited in handling excited states, the SEA naturally connects the solution of excited states to the supersymmetric hierarchy, making the calculation of corrections for all states straightforward. The work asserts that this approach is particularly powerful for systems where the coefficients in the potential expansion are not polynomials, as it does not rely on reducing the problem to algebraic recurrence relations but rather solves the differential hierarchy directly.