🔢 mathematics

Field-Modified Quantum Potentials from Tridiagonal Representations: Analytical Spectra and Galerkin Simulations

This paper develops an effective radial framework for charged spinless particles in collinear electric and magnetic fields under a weak-field approximation, utilizing Tridiagonal representation approaches in Laguerre and Jacobi bases to derive new analytical energy spectra and numerical benchmarks for field-modified potentials.

Original authors: Tunde Joseph Osunmusanmi, Berihu Teklu

Published 2026-09-22

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Tunde Joseph Osunmusanmi, Berihu Teklu

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 a tiny, charged particle, like an electron stripped of its spin, floating in a vacuum. In the quiet of the quantum world, this particle usually follows a predictable path, circling a center of force much like a planet orbits a star. But what happens when we introduce the invisible hands of nature: electric and magnetic fields? These fields are not just background noise; they are powerful forces that can stretch, squeeze, and twist the very shape of the particle's possible existence. When scientists try to predict how a particle behaves under the combined influence of both an electric field and a magnetic field, the mathematics becomes incredibly difficult. The equations that describe this dance are so complex that they often cannot be solved exactly, forcing researchers to rely on approximations or massive computer simulations. Understanding this behavior is crucial for everything from designing better atomic clocks to controlling ions in advanced quantum computers, yet the full picture has remained elusive for many specific setups.

A team of researchers at Khalifa University has now carved out a new, manageable path through this mathematical wilderness. They did not attempt to solve the entire, messy three-dimensional problem at once. Instead, they devised a clever way to simplify the scene. They imagined the particle moving in a very specific, aligned configuration where the electric and magnetic fields run parallel to each other. By focusing on a narrow slice of the particle's possible orientations and making a careful assumption that the magnetic field is weak enough not to cause certain secondary effects, they were able to strip away the most confusing parts of the equation. What remained was a streamlined, one-dimensional model that captures the essential physics of the situation without the overwhelming complexity. This new model acts as a precise map, allowing them to see exactly how the particle's energy levels and its wave-like shape change when the fields are turned on.

Using a sophisticated mathematical toolkit known as the Tridiagonal Representation Approach, the team translated this simplified model into two different languages of mathematics, each revealing a different kind of potential energy landscape. In the first approach, they discovered that the particle behaves as if it were trapped in a valley formed by two competing forces: one that pulls it inward like gravity and another that pushes it outward like a stretched spring. This creates a stable pocket where the particle can settle into distinct, quantized energy states. The researchers were able to write down exact formulas for these energy levels and the shapes of the waves that describe the particle's position. They found that the energy of these states shifts in a predictable way depending on the strength of the magnetic field, a result that can be calculated directly without needing a computer.

In the second approach, the researchers mapped the particle's movement onto a finite range, effectively trapping it between two walls. Here, the potential energy landscape looks different, rising sharply at the boundaries to keep the particle confined. This setup led to a more complex set of rules for the particle's energy, which could not be solved with a single neat formula. Instead, the team used a numerical method, essentially building a giant grid of numbers to simulate the system and find the answers. Their calculations showed that even in this more complicated scenario, the energy levels stabilize quickly as the simulation becomes more detailed, giving them high confidence in the results. They found that the particle's behavior in this confined space is governed by a unique set of mathematical patterns that had not been fully explored before.

The beauty of this work lies in its clarity and its utility. The researchers did not just find new numbers; they created a set of benchmark models that other scientists can use to test their own theories and simulations. Because they kept their assumptions explicit and their methods transparent, anyone can check their work or use their formulas to study similar systems. The paper explicitly states that these results are not the exact solution to the most general, chaotic version of the problem, but rather a highly accurate and useful approximation for a specific, controlled environment. This distinction is vital: it means the findings are reliable within the boundaries they set, offering a solid foundation for future exploration.

Looking ahead, the authors suggest that these models could be tested in real-world experiments, particularly with single ions trapped in electromagnetic fields, a technology that is already advanced enough to measure these subtle shifts. The work also opens doors for studying similar effects in condensed matter physics, such as in tiny structures made of graphene or other materials where electrons are confined to small spaces. By providing a clear, analytical way to understand how electric and magnetic fields reshape the quantum world, this research offers a new lens through which scientists can view the fundamental building blocks of matter, turning a previously intractable problem into a tractable and insightful story of confinement and energy.

Technical Summary: Field-Modified Quantum Potentials from Tridiagonal Representations

Problem Statement
The paper addresses the challenge of analyzing the response of atomic and few-body quantum systems to combined external electric and magnetic fields. While exact treatments of the three-dimensional Stark–Zeeman problem are generally nonseparable and computationally demanding, the authors propose an effective mathematical formulation. The core problem is to derive a reduced, analytically tractable radial Schrödinger equation that captures the dominant field-induced modifications to bound-state spectra without solving the full nonseparable 3D problem. The study focuses on a charged spinless particle in uniform, parallel electric ( $E$ ) and magnetic ( $B$ ) fields, operating under a weak-magnetic-field approximation where the diamagnetic term is neglected.

Methodology
The authors employ the Tridiagonal Representation Approach (TRA), a method closely related to the J-matrix program, which utilizes square-integrable orthogonal polynomial bases to convert the Schrödinger equation into three-term recurrence relations.

Effective Model Construction:
- The authors start with the standard Hamiltonian for a particle in uniform $E$ and $B$ fields (parallel to the $z$ -axis).
- They introduce a specific approximation: replacing the angular factor $\cos \theta$ in the electric dipole term with a constant parameter $C = \cos \theta$ , where $0 < \theta < 90^\circ$ . This projects the anisotropic coupling onto a restricted class of angularly localized states.
- This results in an effective radial Schrödinger equation where the potential includes a linear electric term ( $-qECr$ ) and a linear magnetic energy shift ( $qB m_\phi \hbar / 2m$ ).
Basis Set Analysis:
The study applies the TRA in two distinct bases to identify solvable potential forms:
- Laguerre Basis: Used for the radial half-line $r \in [0, \infty)$ . The authors utilize the "indirect mode" (or search mode) of TRA, imposing a tridiagonal constraint on the Hamiltonian matrix elements. This algebraic constraint determines the functional form of the unknown effective potential $V(r)$ .
- Jacobi Basis: Used to map the radial half-line to a finite interval $x \in [-1, 1]$ via the transformation $x = (r-\lambda)/(r+\lambda)$ . This approach treats the resulting singular confining potential through a truncated generalized eigenvalue problem (Galerkin simulation).

Key Contributions and Results

Identification of Solvable Potentials (Laguerre Basis):
By enforcing the tridiagonal constraint, the authors derive two specific forms for the effective central potential $V(r)$ :
1. Linear-plus-Inverse Potential: $V(r) = ar + b/r$. This potential combines a linear confinement term (analogous to a constant external field) and an inverse-distance term (Coulomb-like). The authors show that for positive coefficients, this potential supports bound states with a minimum at a finite radius.
2. Linear-plus-Quadratic Potential: $V(r) = \kappa r + \eta r^2$ . This represents a superposition of linear and quadratic terms, describing a confined system under a constant external force. Stability requires the quadratic coefficient $\eta$ to be positive.
Analytical Spectra and Wavefunctions:
- For both Laguerre-based potentials, the recurrence relations for the expansion coefficients are identified as belonging to the Meixner–Pollaczek and Meixner polynomial families.
- Through analytic continuation of the Meixner–Pollaczek parameters, the authors derive closed-form expressions for the discrete energy spectra ( $\xi_k$ ) and the corresponding wavefunctions.
- The total energy spectrum is explicitly coupled to the magnetic field via a linear shift: $\xi_{total} = \xi_k + \frac{qB m_\phi \hbar}{2m}$ .
Numerical Spectra (Jacobi Basis):
- The Jacobi basis construction yields a singular, boundary-confined effective potential.
- The resulting recurrence relation does not map to a standard classical orthogonal polynomial family. Consequently, the energy spectrum is not obtained via a closed-form analytic formula but is computed numerically.
- The authors formulate a generalized eigenvalue problem ( $H \mathbf{f} \approx \xi \Omega \mathbf{f}$ ) and demonstrate convergence of the low-lying eigenvalues as the basis size $N$ increases, with errors scaling as $O(N^{-2})$ .
Visualizations:
The paper provides plots of the effective potentials $V(r)$ and discrete energy levels for various orientation angles $\theta$ , illustrating how the potential shape and spectrum evolve with the parameter $C = \cos \theta$ . Wavefunctions for localized states ( $l=0$ ) and angularly excited states ( $l=1$ ) are also plotted.

Significance and Claims
The paper claims to provide new analytical and numerical benchmark models for field-modified potential functions. Its significance lies in:

Analytical Tractability: It offers exact solutions for field-modified potentials that are typically intractable in full 3D, preserving the self-adjoint radial structure while remaining solvable.
Explicit Assumptions: The authors emphasize that these are spectra of projected radial models, not exact hydrogenic Stark–Zeeman spectra. The validity is explicitly limited to weak magnetic fields and angularly localized states.
Benchmarking: The derived models serve as testbeds for investigating how external fields modify central potentials, bound-state spectra, and confinement properties.
Experimental Relevance: The authors suggest that the linear-plus-quadratic branch of their model is particularly suitable for experimental validation using single trapped ions in Penning traps or gate-defined quantum dots, where equilibrium position shifts and level displacements can be measured.

The work concludes that the Tridiagonal Representation Approach provides a systematic algebraic route to constructing solvable quantum systems in external fields, offering a flexible framework for exploring field-modified dynamics in various physical settings.

More like this