Physicists Solve Schrödinger Equation for an Exotic Langevin-Type Potential

A team of theoretical physicists has worked out, for the first time, how a quantum particle behaves inside a modified Langevin-type potential — deriving both the discrete energy levels of bound states and the continuous solutions that describe scattering. The study, published as open access research in Discover Physics on 6 January 2026, combines several competing analytical techniques with a numerical benchmark, and it uncovers striking scattering phenomena including resonances and the celebrated Ramsauer-Townsend effect.

The Langevin function, a mathematical form long valued for describing paramagnetic materials and the dielectric properties of insulators, takes the basic shape of a hyperbolic cotangent minus an inverse-distance term. The authors, E. Omugbe, T. A. Nahool, E. S. Eyube, E. M. Khokha, U. S. Okorie and C. A. Onate, modified this function into an attractive potential with two depth parameters, U0 and U1, plus a range parameter alpha, and appended a Coulomb-like term. Because the modified form contains the Coulomb function used to model hydrogen-like atoms, it inherits a singular behaviour at the origin while settling to a finite value at large distances — a combination that makes it physically interesting but mathematically stubborn.

That stubbornness is the crux of the work. The modified Langevin potential does not admit exact analytical solutions of the energy and wave functions, even without the centrifugal barrier that arises for particles with nonzero orbital angular momentum. To make progress, the researchers applied the Greene-Aldrich approximation scheme, which replaces the awkward inverse-square centrifugal term with an exponential-based surrogate valid for small values of the range parameter. With that substitution, the radial Schrödinger equation transforms into a hypergeometric differential equation whose solutions can be written in closed form.

From the asymptotic behaviour of the hypergeometric wave function, the team extracted the scattering phase shift — the angle-like quantity that encodes how an incoming particle’s wave is distorted by the potential — as an expression involving gamma functions of the asymptotic wave number and the potential parameters. They also derived the normalization factor of the continuous-state wave function and the total scattering cross section, the effective area a particle presents as a target. The bound-state energies then emerged elegantly from the poles of the scattering amplitude, where the gamma functions blow up at negative integers, yielding a closed-form energy equation for any orbital angular momentum quantum number.

To check these functional-analysis results, the authors brought in two independent approximation methods. The Wentzel-Kramers-Brillouin (WKB) approach, a semiclassical technique suited to slowly varying potentials, produced its own quantization condition and energy formula through a two-turning-point integral. Perturbation theory offered a third route, treating part of the potential as a small disturbance to an exactly solvable parent. Intriguingly, the choice of parent matters: when the Coulomb term dominates (U1 greater than U0), the Coulomb-plus-kinetic Hamiltonian serves as the unperturbed system; when U0 exceeds U1, that arrangement breaks down and the roles must be reversed, with the resulting correction integral evaluated numerically using Mathematica.

The decisive test came from the Matrix-Numerov method, a numerical scheme for solving differential equations that yields highly accurate bound-state eigenvalues. Comparing all four approaches across tabulated energy spectra in atomic units, the team found that the perturbation energies sit closest to the numerical values when the Coulomb term dominates and the orbital quantum number is zero, while the WKB energies converge best as the angular momentum increases. The eigenvalues themselves rise with the range parameter alpha before saturating, and they increase with the angular momentum quantum number — trends consistent across every method tested.

The scattering results carry the paper’s most visually dramatic findings. The phase shift grows with the angular momentum quantum number, but for low values it takes on a saw-tooth shape, with negative phase shifts at zero angular momentum — a direct signature of the strongly attractive potential in the absence of a centrifugal barrier. As the angular momentum climbs, the saw-tooth gives way to a pulsating waveform, which the authors attribute to quantum fluctuations of the centrifugal barrier that make the potential effectively less attractive and push scattering to higher energies. The phase shift likewise increases with the range parameter, the wave number and the Coulomb depth parameter.

The total cross section proved even richer. Plotted against the range parameter, it shows densely packed peaks at small alpha that spread into broader oscillations as alpha grows — peaks corresponding to resonant scattering, where the incident particles interact strongly with the potential. The zeros of the cross section, meanwhile, signal the Ramsauer-Townsend effect, the famous near-transparent scattering regime arising from minimal interaction between projectile and target. Against the asymptotic wave number, the cross section simply decays toward zero. But plotted against the potential depth parameters U0 and U1, the behaviour splits into three regimes: dense resonances interleaved with Ramsauer-Townsend minima at small values, more spread-out oscillations at intermediate values, and finally, beyond a critical depth, a saturated pulsating pattern in which both resonance and transparency effects are suppressed — a regime the authors link to the strongly attractive potential driving incident particles to ultralow energies.

The authors note that, to the best of their knowledge, the energy spectra and continuous-state solutions of the Langevin-type potential had not been studied before this work. Their results demonstrate that the centrifugal barrier and the attractive character of the potential jointly govern both the phase shift and the cross section, and they show that complementary approximation schemes — functional analysis, WKB and perturbation theory — can be cross-validated against numerical benchmarks even for potentials that resist exact solution. Because the modified Langevin potential may be used to describe binding energies and scattering states of atomic particles, the analytical formulas derived here offer a compact toolkit for future studies of quantum systems whose interactions blend Coulomb-like attraction with a saturating short-range form.

Subject of Research: Theoretical study of bound and scattering state solutions of the Schrödinger equation under a modified Langevin potential

Article Title: Energy spectra and scattering state solutions of a particle confined in a Langevin-type function

Article References: Omugbe, E., Nahool, T. A., Eyube, E. S., Khokha, E. M., Okorie, U. S., & Onate, C. A. (2026). Energy spectra and scattering state solutions of a particle confined in a Langevin-type function. Discover Physics, 2(1), Article 1. https://doi.org/10.1007/s44418-025-00003-7

Image Credits: AI Generated

DOI: 10.1007/s44418-025-00003-7

Keywords: quantum mechanics, Schrödinger equation, scattering theory, phase shift, Ramsauer-Townsend effect, WKB approximation, perturbation theory, Langevin potential, bound states, cross section, Matrix-Numerov method, theoretical physics