Coulomb-Corrected Harmonic Oscillator Trap

• Coulomb-corrected harmonic oscillator traps are confined quantum systems where harmonic trapping is modified by long-range Coulomb or screened interactions to produce intricate correlation dynamics.
• They encompass models from two-electron harmonium to many-body scattering formalisms, utilizing Coulomb-dressed Green’s functions and reinterpreted Hamiltonians to address correlation effects.
• Recent studies reveal these systems yield tunable energy spectra, Wigner localization, and collective mode shifts, with validations from variational methods and numerical simulations.

Searching arXiv for the cited papers to ground the article in recent and relevant literature. A Coulomb-corrected harmonic oscillator trap is a confined system in which quadratic trapping is supplemented, deformed, or effectively reinterpreted by Coulomb, Coulomb-like, or screened-Coulomb interactions. Across the literature, this includes two-electron harmonium with explicit electron–electron repulsion, trapped few- and many-body systems with pairwise $1/r$ or $1/|x|$ interactions, Lorentz-symmetry-breaking constructions that induce an effective $1/\rho$ term for a neutral spin-$\tfrac12$ particle, and harmonic-trap quantization schemes used to extract Coulomb-modified scattering information from discrete spectra (Kroeze et al., 2015, Abraham et al., 2012, Bakke et al., 2012, Guo, 2021). The expression therefore does not identify a single canonical Hamiltonian; rather, it names a family of problems whose unifying feature is the competition between harmonic confinement and long-ranged electrostatic or electrostatic-like structure.

1. Scope and principal realizations

In bound-state problems, the simplest realization is the relative motion of two charged particles in an isotropic harmonic trap. For two electrons, the effective radial potential is the sum of a harmonic term, a Coulomb repulsion, and the centrifugal barrier, which is the standard harmonium setting (Kroeze et al., 2015). In few- and many-body formulations, the trap remains harmonic at the one-body level, while Coulomb effects enter through pairwise interactions such as $a/|\mathbf r_i-\mathbf r_j|$ or $g/|x_i-x_j|$, so the “correction” is not a perturbative modification of a one-particle oscillator spectrum but a correlation problem in relative coordinates and configuration space (Abraham et al., 2012, Cremon, 2012).

A second realization is formal rather than purely Hamiltonian. In trapped two-body scattering, the harmonic oscillator potential is introduced to discretize the continuum, and the Coulomb correction enters through Coulomb-dressed Green’s functions and Coulomb-modified quantization conditions rather than through a simple additive potential term (Guo, 2021, Bagnarol et al., 2024). A third realization is induced rather than fundamental: in a Lorentz-violating background, a neutral spin-$\tfrac12$ particle in a two-dimensional harmonic oscillator can acquire an effective Coulomb-like $1/\rho$ interaction, generated by the nonminimal coupling to external fields (Bakke et al., 2012).

Screened-Coulomb generalizations extend the same logic. A strongly coupled Yukawa plasma confined by a one-dimensional harmonic trap in the $z$-direction forms quasi-2D layers, and harmonically confined classical or quantum charges develop radial shells and intra-shell order at strong coupling (Pan et al., 2021, Wrighton et al., 2022). This suggests that “Coulomb-corrected harmonic oscillator trap” is best understood as an umbrella category spanning exact few-body models, many-body correlation problems, and trap-assisted scattering formalisms.

2. Canonical Hamiltonians and coordinate reduction

The two-body bound-state prototype is the relative-coordinate Schrödinger equation

$$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$

Here the center-of-mass sector is an ordinary 3D harmonic oscillator, so the nontrivial physics is entirely in the relative coordinate $1/|x|$0 (Kroeze et al., 2015). The same separation into center-of-mass and relative motion underlies trap problems more generally, including harmonically confined two-particle systems with other interparticle interactions (Chia et al., 2023).

For interacting trapped many-body systems, the Coulomb-corrected oscillator appears as a sum of one-body harmonic terms and pair interactions. In scaled units, the time-dependent Schrödinger equation studied for $1/|x|$1 identical particles in a one-dimensional harmonic trap with Coulomb interaction is

$1/|x|$2

with $1/|x|$3 interpreted as the ratio of interaction energy to confinement energy (Abraham et al., 2012). In the quasi-1D few-body case, the interaction was chosen explicitly as

$1/|x|$4

with quasi-one-dimensional regularization implemented by taking the transverse harmonic oscillator length to be $1/|x|$5 of the axial oscillator length (Cremon, 2012).

A related but anisotropic realization occurs in quasi-2D Yukawa matter. There the motion is free in the $1/|x|$6–$1/|x|$7 plane, harmonic only along $1/|x|$8,

$1/|x|$9

and the pair force is screened Coulomb,

$1/\rho$0

The control parameters are $1/\rho$1, where $1/\rho$2 measures confinement strength (Pan et al., 2021).

3. Correlation structure and ground-state descriptions

For a trapped few-body system in one dimension, a particularly compact description is provided by a Jastrow-type variational ansatz,

$1/\rho$3

with pair factor

$1/\rho$4

where $1/\rho$5 is the relative ground-state wavefunction of the two-particle problem and $1/\rho$6 is a fitted scaling parameter (Cremon, 2012). In this construction, even $1/\rho$7 yields bosonic symmetry and odd $1/\rho$8 yields fermionic antisymmetry; the special cases $1/\rho$9, $\tfrac12$0, and $\tfrac12$1 recover noninteracting bosons, a Slater determinant of the lowest $\tfrac12$2 orbitals for fermions, and Tonks–Girardeau bosons, respectively (Cremon, 2012).

The ansatz was benchmarked for $\tfrac12$3 against configuration interaction / exact diagonalization for bosons and spin-polarized fermions at $\tfrac12$4 and $\tfrac12$5. The energy as a function of $\tfrac12$6 has a clear minimum near $\tfrac12$7, and the reported residual error is of order $\tfrac12$8 of the gap to the first excited state (Cremon, 2012). At $\tfrac12$9, the repulsion is strong enough to produce Wigner localization, and bosons and fermions then have very similar energies (Cremon, 2012). The paper is explicit that this performance is established only for the very small systems tested, and that no general conclusion about $a/|\mathbf r_i-\mathbf r_j|$0-dependence of $a/|\mathbf r_i-\mathbf r_j|$1 can be drawn.

At larger $a/|\mathbf r_i-\mathbf r_j|$2 and in the classical or semiclassical regime, the dominant structural phenomenon is correlation-driven self-organization. For charges in an isotropic harmonic trap,

$a/|\mathbf r_i-\mathbf r_j|$3

the equilibrium density is obtained from density functional theory with excess free energy encoded by the direct pair correlation function $a/|\mathbf r_i-\mathbf r_j|$4. In the strongly coupled fluid phase, the adjusted hypernetted chain approximation applied to one-component-plasma correlations reproduces shell emergence, shell positions, shell occupancies, shell widths, and their evolution with $a/|\mathbf r_i-\mathbf r_j|$5 and $a/|\mathbf r_i-\mathbf r_j|$6, in good quantitative agreement with Monte Carlo and molecular dynamics simulations (Wrighton et al., 2022). At even stronger coupling, particles localize angularly within shells in a Thomson-like pattern, and the angular pair correlation function develops narrow peaks, including doublets and triplets (Wrighton et al., 2022).

The same mechanism appears in screened form for Yukawa interactions under one-dimensional harmonic confinement. Without interactions the density along $a/|\mathbf r_i-\mathbf r_j|$7 is Gaussian,

$a/|\mathbf r_i-\mathbf r_j|$8

but strong coupling reshapes this profile into multiple density maxima. Mean field remains monotonic and cannot produce layer splitting, whereas hypernetted chain theory reproduces the number and position of peaks seen in molecular dynamics (Pan et al., 2021). In that setting, $a/|\mathbf r_i-\mathbf r_j|$9 mainly determines how many layers form, while $g/|x_i-x_j|$0 mostly determines how pronounced they are (Pan et al., 2021).

4. Spectral crossover, collective modes, and strong-coupling regimes

The two-electron harmonium spectrum interpolates between two analytically distinct limits. For large trap frequency, the harmonic term dominates and the energies approach the free oscillator result

$g/|x_i-x_j|$1

For small trap frequency, Coulomb repulsion dominates and the energy scales as

$g/|x_i-x_j|$2

The problem is quasi-exactly solvable only at the Taut points, a countable set of trap frequencies for which the polynomial ansatz closes; outside those points the spectrum is obtained numerically or by harmonic approximations to the effective potential (Kroeze et al., 2015). The low-frequency regime is interpreted as strongly correlated and Wigner-crystal-like, and the dependence on angular momentum disappears asymptotically as $g/|x_i-x_j|$3 (Kroeze et al., 2015).

The corresponding many-particle picture is dynamical as well as structural. In a one-dimensional harmonic trap with Coulomb interaction, the breathing signal contains two frequencies: a universal center-of-mass mode

$g/|x_i-x_j|$4

and an interaction-dependent relative mode $g/|x_i-x_j|$5 (Abraham et al., 2012). The exact limiting values are $g/|x_i-x_j|$6 for $g/|x_i-x_j|$7 and $g/|x_i-x_j|$8 for $g/|x_i-x_j|$9, so

$\tfrac12$0

Configuration interaction benchmarks were carried out up to $\tfrac12$1, while time-dependent Hartree–Fock was used up to $\tfrac12$2; for $\tfrac12$3, $\tfrac12$4 decreases with $\tfrac12$5, reaches a minimum around $\tfrac12$6 in TDHF and $\tfrac12$7 in CI, and then increases again (Abraham et al., 2012). The TDHF spectra also contain unphysical artifacts, so peak assignment near $\tfrac12$8 requires care (Abraham et al., 2012).

A common simplification is to treat the Coulomb term as a small perturbation of the oscillator. The cited work shows that this is only regime-dependent. At high trap frequency the approximation is reasonable for harmonium, but at low frequency the Coulomb contribution changes the scaling law itself; in few-body quasi-1D systems it produces Wigner localization; and in classical or semiclassical many-body traps it generates shells, layers, and modified collective dynamics rather than a mere level shift (Kroeze et al., 2015, Cremon, 2012, Wrighton et al., 2022).

5. Quantization conditions and trap-based scattering formalisms

In trapped scattering theory, the Coulomb-corrected harmonic oscillator trap is used to convert continuum information into a discrete-spectrum problem. The central object is the Coulomb Green’s function in the trap,

$\tfrac12$9

which obeys the Dyson equation

$1/\rho$0

For separable zero-range short-range interactions, the trapped spectrum satisfies a Coulomb-modified BERW/Lüscher-type determinant condition,

$1/\rho$1

where $1/\rho$2 is built from the short-distance limit of the difference between the trapped and infinite-volume Coulomb Green’s functions (Guo, 2021). For the harmonic oscillator trap,

$1/\rho$3

the pure HO Green’s function is known explicitly in terms of Whittaker functions, but the paper does not derive a simple exact closed-form HO quantization condition with all-orders Coulomb corrections; instead it gives the HO Green’s function, the Dyson equation for Coulomb resummation, and a perturbative expansion strategy (Guo, 2021).

A subsequent numerical treatment of the generalized Coulomb-corrected BERW formula solved the required Green’s-function equations by two methods: successive approximation with stabilization and direct matrix inversion (Bagnarol et al., 2024). Applied to $1/\rho$4-$1/\rho$5 scattering with the phenomenological short-range interaction of Ali et al., these methods reproduce Numerov phase shifts for all $1/\rho$6 and $1/\rho$7 channels. For sufficiently weak traps,

$1/\rho$8

or equivalently

$1/\rho$9

the extracted phase shifts agree with Numerov results to within about $z$0 (Bagnarol et al., 2024). Numerical difficulties occur close to the poles of the HO Green’s function,

$z$1

where both iterative and inversion-based solvers become unstable (Bagnarol et al., 2024).

Within oscillator-basis scattering theory, the HORSE formalism has also been extended to the charged case by introducing an effective asymptotic diagonal Coulomb term in the three-term recurrence relation,

$z$2

This term is not the exact diagonal Coulomb matrix element; rather, it represents the influence of the full Coulomb interaction in the asymptotic oscillator region (Yanikov et al., 6 Apr 2025). The resulting method gives good convergence to direct Numerov solutions, is comparable in accuracy to the Coulomb-truncation approach of Bang et al., and is more accurate than the older summation-based treatment from the Kiev group in the example studied (Yanikov et al., 6 Apr 2025).

6. Induced Coulomb-like oscillators and conceptual limits

Not all Coulomb-corrected oscillator traps arise from literal charged particles. In the Lorentz-violating construction for a neutral spin-$z$3 particle, the nonminimal coupling

$z$4

is specialized to a fixed space-like vector field $z$5 and a uniform magnetic field $z$6, so that $z$7 generates an effective Coulomb-like interaction (Bakke et al., 2012). After separation of variables, the radial equation for the harmonic oscillator sector becomes

$z$8

which is a Coulomb-corrected harmonic oscillator with a $z$9 term generated by Lorentz symmetry breaking rather than by direct electrostatics (Bakke et al., 2012). After the substitution $$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$0, the problem reduces to a biconfluent Heun equation, and polynomial truncation quantizes the allowed frequencies $$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$1. The spectrum is

$$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$2

with the background contributing the additive shift $$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$3 (Bakke et al., 2012).

Several conceptual limits follow from the broader literature. First, a Coulomb-corrected harmonic oscillator trap is not generally an exactly solvable model: the HO scattering problem with Coulomb is handled through Green’s functions and Dyson equations rather than a universal closed-form spectrum (Guo, 2021), the two-electron bound-state problem is quasi-exact only at special Taut points (Kroeze et al., 2015), and the Jastrow description of the quasi-1D few-body case is explicitly validated only for $$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$4 (Cremon, 2012). Second, “Coulomb-corrected” need not mean bare $$-\frac{\hbar^2}{2\mu}\frac{d^2u}{dr^2} +\left( \frac{1}{2}\mu\omega^2r^2 +\frac{e^2}{4\pi\varepsilon_0 r} +\frac{\hbar^2l(l+1)}{2\mu r^2} \right)u=Eu.$$5: screened Yukawa repulsion produces analogous confinement-induced layering (Pan et al., 2021). Third, the correction is often nonperturbative in physical effect even when the formalism is perturbative or variational. This suggests that the most stable definition is operational: a Coulomb-corrected harmonic oscillator trap is any harmonic confinement problem whose spectral, structural, or scattering content is materially controlled by Coulombic correlations, Coulomb-dressed propagators, or an induced Coulomb-like interaction.