🔢 mathematics

An Accurate Zeroth-Order Perturbation Theory for Solving Power-Law Potentials within the Framework of the Asymptotic Iteration Method

This paper introduces an accurate zeroth-order perturbation theory using the Asymptotic Iteration Method with exponential and hyperbolic ansatz substitutions to solve the Schrödinger equation for spherically symmetric power-law potentials, effectively avoiding oscillatory behaviors and demonstrating high agreement with other methods for quartic and sextic singularities.

Original authors: A. J. Sous, Ibsal Assi, Nasser Saad

Published 2026-10-20

📖 4 min read🧠 Deep dive

CC0 1.0

Original authors: A. J. Sous, Ibsal Assi, Nasser Saad

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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

In the microscopic world of atoms and molecules, particles do not move along predictable tracks like planets orbiting a sun. Instead, they exist as clouds of probability, their behavior governed by a fundamental rule known as the Schrödinger equation. This mathematical law acts as the instruction manual for the universe at its smallest scale, dictating how energy is stored and how particles interact. For simple systems, scientists can solve this equation exactly to find precise answers. However, nature is often messy. When particles encounter complex forces, such as those found in molecules with multiple atoms or in exotic states of matter, the equations become incredibly difficult to crack. In these cases, the forces can become infinitely strong at certain points, creating mathematical "singularities" that break standard calculation methods. To understand these systems, physicists must rely on numerical techniques, using powerful computers to approximate the answers. Yet, even with modern computing, some of these problems remain stubbornly resistant, causing standard calculation tools to stumble and produce unreliable, oscillating results that fail to settle on a single answer.

A team of researchers has developed a new strategy to bypass these computational roadblocks, specifically targeting a class of difficult problems involving power-law potentials. These are forces that change strength based on distance, often appearing in models of molecular interactions. The researchers, working within a framework called the asymptotic iteration method, introduced a clever workaround that avoids the instability plaguing previous attempts. Instead of trying to solve the difficult equation directly, they proposed a two-step process. First, they replaced the original, unruly force with a slightly different, smoother version that is mathematically easier to handle. This new version is not a perfect match, but the difference between the two is treated as a tiny, manageable disturbance. By solving the easier version first, they obtain a solid foundation. They then apply a technique known as zeroth-order perturbation theory, which essentially means they accept this first, approximate solution as highly accurate without needing to perform complex, multi-step corrections. The key to their success lies in choosing a specific mathematical shape for this easier version, one that includes a small, adjustable knob. By turning this knob to a very small setting, the researchers ensure the approximation is close enough to the real thing to be useful, while remaining stable enough for the computer to solve without glitching.

The team tested this approach on three distinct types of physical problems, ranging from the well-understood behavior of charged particles to more complex scenarios involving sharp spikes in energy. In the first test, they applied their method to the Coulomb potential, which describes the attraction between an electron and an atomic nucleus. This is a problem with a known exact answer, making it a perfect benchmark. Using their new method, they calculated the energy levels of the system and found results that matched the known answers with remarkable precision, differing only in the eighth decimal place. They also compared two different mathematical shapes for their approximation—one based on exponential growth and another on hyperbolic curves. Both shapes worked well, producing nearly identical results and confirming that the method is robust regardless of the specific shape chosen.

Moving beyond simple cases, the researchers tackled a more challenging scenario known as the spiked isotropic oscillator. This system involves a particle moving in a potential that rises sharply, like a spike, due to a specific type of singularity. Standard methods often struggle here, but the team's approach yielded energy values that aligned closely with results from other advanced simulation techniques. They repeated this success with an even more complex system featuring a sextic singularity, a force that grows extremely rapidly at short distances. In every instance, the method provided stable, accurate energy spectra without the erratic fluctuations that typically plague these calculations. The researchers did not claim to have solved every possible problem in physics, but they demonstrated that their technique offers a reliable and efficient path through a landscape that was previously difficult to navigate. By treating the difference between the real world and their simplified model as a negligible perturbation, they turned a chaotic calculation into a steady, predictable process, offering a new tool for physicists studying the intricate dance of atoms and molecules.

Technical Summary: An Accurate Zeroth-Order Perturbation Theory for Solving Power-Law Potentials within the Framework of the Asymptotic Iteration Method

Problem Statement
The paper addresses the challenge of solving the radial Schrödinger equation for spherically symmetric singular power-law potentials, which are prevalent in molecular and atomic physics (e.g., Lennard-Jones potentials). While the Asymptotic Iteration Method (AIM) is a recognized numerical technique for eigenvalue problems, its direct application to these specific singular potentials is often hindered by oscillatory behavior in numerical computations. This oscillation leads to convergence difficulties, particularly when determining the energy spectrum. Previous attempts to apply perturbation theory within AIM were noted to be computationally demanding, sensitive to initial seed points, and reliant on higher-order corrections that were difficult to manage.

Methodology
The authors propose a zeroth-order perturbation theory formulated specifically to integrate with the AIM framework, aiming to bypass the oscillatory convergence issues. The core methodology involves the following steps:

Potential Decomposition: The original effective potential $V(r)$ is decomposed by adding and subtracting an Ansatz potential $V_0(\lambda, r)$ . This creates a new equation where the difference $U_\lambda(r) = V(r) - V_0(\lambda, r)$ is treated as a small perturbation.
Ansatz Substitutions: Two specific classes of Ansatz potentials are introduced to approximate the singular terms of $V(r)$ $V (r)$ :
- Exponential Ansatz (EA): Utilizes terms of the form $(e^{\lambda r} - 1)^{-k}$ .
- Hyperbolic Ansatz (HA): Utilizes terms of the form $(\sinh \lambda r)^{-k}$ .
  In both cases, $\lambda$ is a free parameter chosen to be sufficiently small. The coefficients of the Ansatz are selected such that $U_\lambda(r) \to 0$ as $\lambda \to 0$ .
Zeroth-Order Solution: The Schrödinger equation is solved for the unperturbed part involving only $V_0(\lambda, r)$ . This is achieved by applying coordinate transformations ( $y = y(r)$ ) and wavefunction factorizations ( $\psi = g(y)f(y)$ ) to map the equation into the standard form required for AIM.
AIM Implementation: The AIM recursive sequences ( $k_n, s_n$ ) are generated to find the roots of the termination condition $\Delta_n = 0$ , yielding the zeroth-order energy spectrum $E^{(0)}_\lambda$ and wavefunctions.
Parameter Control: The accuracy of the zeroth-order solution is controlled by the parameter $\lambda$ . The authors assert that by choosing $\lambda$ small enough, the zeroth-order solution provides high accuracy without requiring complex higher-order perturbation calculations.

Key Contributions and Results
The paper demonstrates the efficacy of this approach through three specific applications, comparing the results against exact solutions or other established numerical methods (such as the TRA method):

Coulomb Potential: The method was applied to the exactly solvable Coulomb problem. Using both EA and HA with $\lambda = 0.001$ , the authors obtained energy eigenvalues that showed excellent agreement with the exact analytical formula. The results indicated that both Ansatz types yield consistent and accurate spectra.
Spiked Isotropic Oscillator: The technique was applied to a potential with a quartic singularity ( $V(r) \propto r^2 + 1/r^4$ ). The calculated energy spectra for the first ten energy levels were compared with results from the TRA method. The AIM results using the Hyperbolic Ansatz showed good agreement with the TRA benchmarks.
Sextic Singularity Potential: To test the method's robustness on higher-order singularities, a potential with a $1/r^6$ term was solved. The authors successfully derived the discrete spectra using the proposed Ansatz substitutions.

In all cases, the authors utilized specific seed points ( $y_0$ ) and iteration counts (e.g., 400 iterations) to ensure convergence. The paper includes plots comparing the exact potentials with the Ansatz potentials, visually confirming the approximation quality for various parameter values.

Significance and Claims
The authors claim that this work introduces a "less demanding approach" compared to previous perturbative applications of AIM. The primary significance lies in the method's ability to:

Avoid Oscillatory Behavior: By reformulating the problem with a specific Ansatz, the method circumvents the numerical oscillations that typically prevent convergence in standard AIM applications for singular potentials.
Achieve Accuracy via Zeroth-Order: The approach demonstrates that high accuracy can be achieved using only zeroth-order perturbation theory, provided the free parameter $\lambda$ is carefully controlled. This reduces the computational burden associated with calculating higher-order corrections.
Broad Applicability: The method is shown to be effective for potentials ranging from the standard Coulomb interaction to complex singularities involving quartic and sextic terms.

The paper concludes that the proposed framework offers a reliable and computationally manageable alternative for solving the Schrödinger equation for a wide class of singular power-law potentials, yielding results that align well with exact solutions and other numerical techniques.

More like this