⚛️ quantum physics

Novel method for evaluating the eigenvalues of the Heun differential equation with an application to the Breit equation

This paper presents a novel method for accurately calculating the eigenvalues of the Breit equation by deriving a second-order Heun differential equation and utilizing both continued fraction and Green function-based semi-infinite determinant approaches, which yield consistent results with 25-digit precision and extend previous findings in the literature.

Original authors: P. J. Rijken, Th. A. Rijken

Published 2026-01-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: P. J. Rijken, Th. A. Rijken

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 the universe as a giant, complex machine where tiny particles like electrons and protons dance around each other. To predict exactly how they move and what energy they have, scientists use a set of rules called equations. One of these famous rules is the Breit equation, which describes how two particles interact, taking into account both their electric charge and their "spin" (a quantum property like a tiny internal magnet).

For decades, solving this equation has been like trying to find a needle in a haystack. The math involved is so messy that it creates a "monster" equation with four tricky points where the math breaks down (called singularities), rather than the usual three. This monster equation is known as the Heun differential equation.

Here is what this paper does, broken down into simple concepts:

1. The Problem: A Math Maze

Think of the Breit equation as a maze. In the past, scientists tried to guess the answer (the energy levels of the particles) by approximating or using trial-and-error methods. They could get close, but they couldn't get the exact answer, especially for very high precision. They were stuck because they didn't have a direct map to the exit.

2. The New Map: Two Different Keys

The authors of this paper, P.J. and Th.A. Rijken, built two new "keys" to unlock the maze and find the exact energy levels.

Key #1: The Infinite Ladder (Continued Fractions)
Imagine a ladder that goes up forever. To find the right energy level, you have to climb this ladder. The paper shows how to calculate the rungs of this ladder using a method called "continued fractions." They found a special algorithm (like a super-fast elevator) that lets them climb this ladder with incredible speed and precision, reaching over 25 digits of accuracy. This is like measuring the height of a mountain to the width of a single atom.
Key #2: The Giant Grid (The Determinant)
The authors also tried a completely different approach. They turned the maze into a giant, infinite grid (a matrix). Imagine a spreadsheet with infinite rows and columns. The "answer" to the energy problem is hidden in the condition where this giant grid collapses or "flattens" in a specific way. They call this a "semi-infinite determinant."
- Why two keys? Having two completely different ways to solve the same puzzle is powerful. If both keys open the door to the exact same spot, you know for sure you haven't made a mistake. The paper proves that both methods agree with each other perfectly.

3. The "Spooky" Ghost State

While solving the maze, the authors noticed something strange. Their math predicted a "ghost" state—a solution that looks like a real particle but isn't. It appears near a specific distance from the center.

The Analogy: Imagine you are looking at a reflection in a mirror. Sometimes, the reflection looks so real you think it's a person, but it's just light. The paper identifies this "ghost" state and explains that it's a mathematical artifact of the specific equation they are using, not a real physical particle. This helps other scientists avoid getting confused by it.

4. Connecting the Dots

The paper also shows that if you tweak the settings of their complex "monster" equation, it turns into a simpler, well-known equation called the Dirac equation (which describes a single electron).

The Metaphor: It's like showing that a complex, high-tech robot can be simplified into a standard bicycle if you remove certain parts. This proves their new, complex math is correct because it matches the known, simpler math when the situation allows it.

5. What They Actually Found

The authors didn't just build the tools; they used them to calculate the energy levels for the first few "steps" (energy states) of the particles.

They calculated these levels for different types of particle spins (like when the particles are spinning in opposite directions or the same direction).
Their results match everything scientists have found in the past but with much higher precision (adding several more correct digits to the numbers).
They confirmed that their method works for the simplest case (where the two particles have the same mass) and laid the groundwork for more complex cases.

Summary

In short, this paper is about building better tools to solve a very old, very difficult math puzzle.

They created two distinct, highly accurate methods (an infinite ladder and a giant grid) to find the exact energy of interacting particles.
They proved these methods work by showing they give the same answer and match known results.
They identified and explained a "ghost" solution that appears in the math.
They showed how their complex math connects back to simpler, established physics.

The paper is a technical triumph in pure mathematics and theoretical physics, offering a clearer, more precise way to understand how these fundamental particles behave, without venturing into medical or engineering applications.

Technical Summary: Novel Method for Evaluating the Eigenvalues of the Heun Differential Equation with an Application to the Breit Equation

Problem Statement
The Breit equation, describing the interaction of two fermions (e.g., electron and proton) via a static Coulomb potential, leads to a second-order differential equation with four singular points. This classifies the equation as a Heun differential equation, which is more singular than the standard confluent hypergeometric equations typically encountered in mathematical physics. While previous studies (e.g., Yamaguchi; Tursunov et al.) have analyzed the Breit equation using variational methods or by matching series solutions, no determination of energy levels currently exists that is directly based on the second-order Heun equation derived from the Breit equation. Furthermore, no explicit expression exists for the expansion of energy levels in powers of the fine-structure constant squared ( $\alpha^2$ ) beyond the first order.

Methodology
The authors develop a theoretical framework to solve the Breit equation for spin-singlet ( $S=0$ ) and specific spin-triplet ( $S=1$ ) states. The methodology proceeds through two distinct, mutually validating approaches:

Continued Fraction (CF) Approach:
- The radial Breit equation is transformed into a standard form of the confluent Heun equation.
- The eigenfunctions are expanded in terms of associated Laguerre polynomials.
- This expansion yields a three-term recurrence relation for the expansion coefficients.
- The recurrence relation is converted into a continued fraction. The eigenvalues are determined by finding the zeros of this continued fraction.
- Numerical evaluation utilizes Miller's backward recurrence algorithm for wave function normalization and Peter Wynn's $\rho$ -algorithm (a non-linear sequence accelerator) to handle the slow convergence of the CF near eigenvalues.
Green Function and Semi-Infinite Determinant Approach:
- The differential equation is reformulated as a homogeneous Fredholm integral equation of the second kind using a Green function constructed from the eigenfunctions of the associated Laguerre operator.
- Expanding the solution to this integral equation in associated Laguerre polynomials leads to a condition on a semi-infinite determinant ( $\det \Delta_{m,n} = 0$ ).
- The authors introduce specific special functions ( $I_{m,n}$ , $S^m_n$ , $\hat{S}^m_n$ ) to simplify the determinant elements.
- To ensure numerical stability and convergence, a "converging factor" is introduced to modify the diagonal elements of the determinant.
- The determinant is solved numerically using three methods: direct calculation (using Dodgson's condensation method), forward recursion, and the CF method for cross-verification.

Key Contributions

Derivation of the Heun Equation: The paper explicitly derives the second-order Heun differential equation from the Breit equation for the spin-singlet case ( $m=M$ ) and generalizes the parameterization for various spin-triplet states.
Novel Determinant Method: The authors present a new method for calculating eigenvalues based on a semi-infinite determinant derived from a Green function integral equation. This provides an alternative to the continued fraction method.
High-Precision Numerical Results: The study calculates energy levels for the first five principal quantum numbers ( $N=1$ to $5$) for spin-singlet and spin-triplet states ( $J=0$ and $J=L$ ). The results achieve accuracy up to 25–30 digits.
Consistency Check: The two distinct mathematical approaches (CF and Determinant) are shown to be consistent within 25 digits of accuracy, validating the numerical methods.
Reproduction of Known Physics: The methods successfully reproduce known results from the literature (e.g., Scott et al., Yamaguchi et al.) and correctly identify the "spurious state" behavior near $\rho \sim 1/y$ previously discussed by Patterson.
Dirac Limit: The framework demonstrates that the Dirac equation solution can be recovered from the determinant condition via a suitable mapping of variables, providing insight into the structure of the determinant.

Results

Eigenvalues: The paper provides tabulated values for the eigenvalues ( $z_n$ ) and energy levels ( $E_n$ ) for the first five levels. For the spin-singlet case, the results match literature values up to 18–25 digits.
Wave Functions: The radial wave functions ( $F, K, G$ ) are efficiently calculated using Miller's backward recurrence algorithm. The paper visualizes these components, confirming the expected asymptotic behavior.
Convergence Analysis: The study analyzes the convergence of the series expansion, proving absolute convergence for $\rho \ge 0$ . It also addresses the numerical challenges of the determinant method for large physical values of the parameter $y$ (related to the fine-structure constant), showing that forward recursion is superior to direct calculation for high precision in this regime.
Analytic Structure: The analysis of the determinant suggests that the series expansion of eigenvalues may involve non-analytical terms, specifically logarithms of the fine-structure constant ( $\ln(1/\alpha)$ ), contrasting with some existing literature but aligning with recent high-order calculations for positronium.

Significance and Claims
The authors claim that their work extends previous calculations by several digits of additional accuracy and provides a robust, alternative theoretical framework for solving the Breit equation. The primary significance lies in:

Providing a direct method to calculate eigenvalues based on the Heun equation without relying on matching series solutions at boundaries.
Offering a semi-infinite determinant condition that serves as a rigorous consistency check for the continued fraction method.
Laying the theoretical groundwork for future work to obtain explicit expansions of energy levels in powers of $\alpha^2$ (beyond first order) and potentially including the full Breit interaction (photon exchange) to resolve state ordering issues (e.g., the reversal of $2^1P_1$ and $2^1S_0$ states).

The paper remains modest regarding the analytic expansion, noting that while the determinant structure suggests the presence of $\ln(1/\alpha)$ terms, an explicit expansion in $\alpha^2$ has not yet been successfully derived in this work but is a target for future research.