🔢 mathematics

High-Precision Computation and PSLQ Identification of Stokes Multipliers for Anharmonic Oscillators

This paper presents a high-precision computational study using arbitrary-precision arithmetic and the PSLQ algorithm to derive exact closed-form expressions for Stokes multipliers in anharmonic oscillators, revealing a new number-theoretic pattern linking their existence to Euler's totient function while identifying the $x^8$ case as a likely new transcendental constant.

Original authors: Jian Zhou

Published 2026-03-31

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jian Zhou

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 the weather. You have a perfect mathematical model, but when you try to calculate the future using that model, the numbers get bigger and bigger until they explode into infinity. This is a common problem in quantum physics called a divergent series.

For decades, physicists have known that even though these calculations "explode," they hide a secret message. If you look at the pattern of how the numbers explode, you can extract a hidden constant called the Stokes Multiplier. Think of this multiplier as the "secret code" or the "DNA" of the system that tells you how the universe behaves when things get really weird (like particles tunneling through walls).

This paper is a massive computational detective story where the author, Jian Zhou, tries to crack the code for a family of quantum systems called Anharmonic Oscillators.

Here is the breakdown of the adventure, explained simply:

1. The Challenge: The Exploding Numbers

The author is studying different types of "quantum springs."

The Simple Spring (M=2): This is the standard spring everyone knows. We already know its secret code.
The Weird Springs (M=3 to M=11): These are springs that get stiffer or weirder as you stretch them. For these, nobody knew the exact secret code.

The author had to calculate thousands of terms of a mathematical series that explodes to infinity. To do this, he used a super-precise calculator (300 digits of precision! That's like measuring the distance to the moon with an error smaller than a grain of sand) to get 1,200 terms of the series.

2. The Tools: The Mathematical Swiss Army Knife

To find the hidden code in the exploding numbers, the author used three powerful tools:

The Ladder (Rayleigh-Schrödinger Recursion): A step-by-step method to climb up the mathematical ladder and generate the exploding numbers.
The Time Machine (Richardson Extrapolation): Since the numbers explode slowly at first, this tool acts like a time machine. It looks at the pattern of the explosion and "fast-forwards" the calculation to see what the number would be if we had infinite terms. It turns a blurry, slow-moving picture into a sharp, high-definition image.
The Pattern Matcher (PSLQ Algorithm): This is the star of the show. Imagine you have a number like 3.14159.... The PSLQ algorithm is like a super-smart detective that asks: "Is this number actually just a combination of famous constants like Pi ( $\pi$ ), the number 2, 3, or special Gamma functions?" It searches for a hidden mathematical equation that equals zero.

3. The Discoveries: Finding the Hidden Formulas

The author ran this pipeline for springs ranging from M=2 to M=11. Here is what happened:

The Confirmed Cases (M=2, 3, 5, 7): For these springs, the Pattern Matcher found beautiful, simple formulas.
- For the M=3 spring, the secret code turned out to be a clean mix of Pi and square roots.
- For the M=5 and M=7 springs, the codes involved complex "Gamma" numbers (which are like advanced factorials), but they still fit a perfect, neat mathematical pattern.
- Analogy: It's like finding that a complex lock opens with a key made of standard Lego bricks.
The Mystery Case (M=4): This is the most exciting part. The author looked at the M=4 spring (the "Octic" oscillator).
- He calculated the number to 30 digits of precision.
- He asked the Pattern Matcher to find a formula using Pi, 2, 3, and Gamma numbers.
- Result: Nothing. The algorithm searched billions of combinations and found no match.
- Analogy: It's like trying to build a house out of standard Lego bricks, but the blueprint requires a piece of the house that is made of a material that doesn't exist in your box. The author suspects the M=4 spring has a "secret ingredient" that is a completely new, unique type of mathematical number that we haven't named yet.

4. The Number Theory Connection: The "Prime" vs. "Composite" Puzzle

The author noticed a strange pattern.

The springs that had neat formulas were M=2, 3, 5, and 7. These are all Prime Numbers.
The spring that didn't have a neat formula was M=4. This is a Composite Number (2 x 2).
The Trap: The author thought, "Maybe only Prime numbers have neat formulas!"
The Twist: He checked M=11 (which is Prime). It didn't have a neat formula either!
The Real Rule: The author realized the rule isn't about primes. It's about a specific number theory concept called Euler's Totient Function.
- Think of this function as a "complexity meter."
- For M=2, 3, 5, 7, the complexity meter was low, allowing for neat formulas.
- For M=4, the meter was low, but something went wrong (it's a unique anomaly).
- For M=6, 8, 9, 10, 11, the meter was high, meaning the math is too messy for a simple formula.

5. Why Does This Matter?

This isn't just about solving math puzzles.

New Physics: Finding these exact formulas helps physicists understand how the universe works at a fundamental level, specifically how particles "tunnel" through energy barriers.
New Math: The fact that M=4 has no known formula suggests there might be a brand new type of number hiding in nature that we haven't discovered yet.
The "M=4" Mystery: The author leaves us with a cliffhanger. Why is the M=4 spring so stubborn? Is it a glitch in the math, or is it a sign of a deeper, stranger reality?

Summary

The author built a super-precise digital microscope to look at the "DNA" of quantum springs. He found that most of them have simple, elegant genetic codes made of known mathematical constants. However, the M=4 spring is an outlier—it seems to carry a genetic code made of a mysterious, unknown substance. This discovery bridges the gap between quantum physics and the ancient world of number theory, suggesting that the universe's deepest secrets might be written in a language we are only just beginning to learn.

1. Problem Statement

The paper addresses the challenge of determining exact closed-form expressions for Stokes multipliers ( $C_M$ ) in the context of one-dimensional anharmonic oscillators defined by the Hamiltonian:
$H = \frac{p^2}{2} + \frac{x^2}{2} + g x^{2M}$
where $M = 2, 3, \dots, 11$ .

While the large-order behavior of perturbation coefficients ( $a_k$ ) is known to be governed by instanton physics and the Stokes multiplier (connecting perturbative and non-perturbative sectors via resurgence theory), exact closed-form expressions for $C_M$ are largely unknown for $M > 2$ . The quartic oscillator ( $M=2$ ) is the only case with a known closed form ( $C_2 = -\sqrt{6}/\pi^{3/2}$ ). For higher $M$ , the complexity of instanton solutions and fluctuation determinants makes analytical derivation difficult. The goal is to use high-precision numerical computation to discover exact identities for $C_M$ and identify cases where no such simple form exists.

2. Methodology

The author employs a three-stage computational pipeline combining arbitrary-precision arithmetic, sequence acceleration, and integer relation detection:

Rayleigh–Schrödinger Recursion:
- Perturbation coefficients $a_k$ are computed up to order $k=1200$ using the harmonic oscillator basis.
- Calculations are performed using the mpmath library in Python with working precision set to 300 decimal digits (plus guard digits) to ensure high accuracy.
- The recursion handles matrix elements of the perturbation $V=x^{2M}$ , scaling as $O(Mk)$ per order.
Parameter Extraction via Richardson Extrapolation:
- The Stokes multiplier $C_M$ is extracted from the asymptotic behavior of the coefficients:
  $a_k \sim C_M \cdot (-1)^{k+1} \cdot A_M^{-k} \cdot \Gamma((M-1)k + b + 1)$
- The instanton action ( $A_M$ ) is derived analytically via a Beta-function parametrization and verified numerically.
- The Gamma-shift parameter ( $b$ ) is confirmed numerically to be $-1/2$ for all $M$ .
- Richardson extrapolation (orders 40–100) is applied to the sequence of estimates for $C_M$ to eliminate $O(1/k)$ subleading corrections, accelerating convergence to 17–30 reliable significant digits.
PSLQ Integer Relation Detection:
- The PSLQ algorithm is used to search for integer relations between $\ln|C_M|$ and a basis of logarithms of fundamental constants: $\ln \Gamma(a/(M-1))$ , $\ln \pi$ , $\ln 2$ , and $\ln 3$ .
- The search is exhaustive, testing coefficient bounds up to 2000 (and up to $10^6$ in detection limit terms) at 30-digit precision.

3. Key Contributions & Results

A. New Exact Closed Forms

The study discovers exact closed-form expressions for Stokes multipliers in cases where none were previously known:

$M=3$ (Sextic): $C_3^2 \pi^4 = 32 \implies C_3 = -4\sqrt{2}/\pi^2$ .
$M=5$ (Decic): $C_5^4 \Gamma(1/4)^4 \pi^5 = 2^{12} \cdot 3^2$ .
$M=7$ (Tetradecic): $C_7^6 \Gamma(1/3)^9 \pi^6 = 2^{20} \cdot 3^3$ .
$M=2$ (Quartic): Confirmed the known result $C_2^2 \pi^3 = 6$ .

B. The $M=4$ Anomaly (Octic Oscillator)

Negative Result: Exhaustive PSLQ searches for $M=4$ (using 1200 terms at 300-digit precision) found no closed form involving $\Gamma(1/3)$ , $\pi$ , and algebraic numbers within the searched bounds.
Significance: This suggests that $C_4$ introduces a genuinely new transcendental number, distinct from the standard Gamma-function values and $\pi$ , likely arising from the complexity of the fluctuation determinant for the $x^8$ potential.

C. Instanton Action and Gamma Shift

Universal Instanton Action: The paper provides a unified Beta-function formula for the instanton action $A_M$ :
$A_M = \left[ \frac{2^{-1/(M-1)} B(\frac{1}{M-1}, \frac{3}{2})}{M-1} \right]^{M-1}$
This is verified numerically to relative errors $< 10^{-30}$ .
Universal Shift: The parameter $b$ in the Gamma function argument is confirmed to be exactly $-1/2$ for all $M \in [2, 11]$ , consistent with one-loop fluctuation determinant arguments.

D. Number-Theoretic Patterns

The paper establishes a correlation between the existence of closed forms and Euler's totient function $\phi(n)$ :

The number of algebraically independent Gamma transcendentals at denominator $M-1$ is $\phi(M-1)/2$ .
Conjecture: Closed forms exist when $\phi(M-1)/2 \le 1$ (with the exception of $M=4$ ).
Refutation of Primality: The existence of a closed form is not determined by whether $M$ is prime. For instance, $M=11$ (prime) yielded no closed form, while $M=4$ (composite) was the specific anomaly among low-totient cases.

4. Significance and Implications

Resurgence Theory: The results provide precise benchmarks for the resurgent trans-series structure of anharmonic oscillators, validating the connection between large-order perturbation theory and instanton physics.
Computational Physics: The work demonstrates the power of combining high-precision arithmetic with PSLQ to discover exact mathematical identities in quantum mechanics, a field traditionally dominated by analytical methods.
Number Theory Connection: It bridges quantum mechanics with classical number theory, suggesting that the "complexity" of the Stokes multiplier is governed by the arithmetic properties of the denominator in the Gamma function arguments ( $\phi(M-1)$ ).
The $M=4$ Mystery: The failure to find a closed form for the octic oscillator ( $M=4$ ) highlights a unique mathematical obstruction in the fluctuation determinant of the $x^8$ potential, potentially pointing to new classes of transcendental numbers in quantum field theory.

5. Data Availability

The author provides full reproducibility, including all Python code (Rayleigh–Schrödinger recursion, Richardson extrapolation, PSLQ scripts) and pre-computed perturbation coefficients, hosted on GitHub under the MIT license.

Conclusion:
This paper successfully extends the knowledge of Stokes multipliers from the well-known quartic case to higher-order anharmonic oscillators. It uncovers new exact identities, identifies a specific anomaly at $M=4$ , and proposes a number-theoretic framework (based on Euler's totient function) to predict the existence of closed forms, fundamentally linking the solvability of quantum mechanical asymptotics to the arithmetic properties of the Gamma function.