© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2303.14982 · The Journal of Chemical Physics

Contactium: A strongly correlated model system

Jerzy Ciosłowski, Berthold‐Georg Englert, Martin-Isbjörn Trappe, Jun Hao Hue

Published May 12, 2023• 3 citations

Abstract

At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.

DOI: 10.1063/5.0150523

Related papers

One-electron reduced density matrices of strongly correlated harmonium atoms

Jerzy Ciosłowski

Explicit asymptotic expressions are derived for the reduced one-electron density matrices (the 1-matrices) of strongly correlated two- and three-electron harmonium atoms in the ground and first excited states. These expressions, which are valid at the limit of small confinement strength ω, yield electron densities and kinetic energies in agreement with the published values. In addition, they reveal the ω(5/6) asymptotic scaling of the exchange components of the electron-electron repulsion energies that differs from the ω(2/3) scaling of their Coulomb and correlation counterparts. The natural orbitals of the totally symmetric ground state of the two-electron harmonium atom are found to possess collective occupancies that follow a mixed power/Gaussian dependence on the angular momentum in variance with the simple power-law prediction of Hill's asymptotics. Providing rigorous constraints on energies as functionals of 1-matrices, these results are expected to facilitate development of approximate implementations of the density matrix functional theory and ensure their proper description of strongly correlated systems.

The Journal of Chemical Physics•Mar 17, 2015•10 citations

Wigner molecules: Natural orbitals of strongly correlated two-electron harmonium

Jerzy Ciosłowski, Marcin Buchowiecki

Explicit asymptotic expressions for natural orbitals and their occupancies are derived for the harmonium atom at the strong-correlation limit at which the confinement strength omega tends to zero. Unlike in systems with moderate correlation effects, the occupancies at the omega-->0 limit (derived from occupation amplitudes with alternating sign patterns) are vanishingly small and asymptotically independent of the angular momentum, forming a geometric progression with the scale factor proportional to omega(1/3) and the common ratio of ca. 0.0186. The radial components of the natural orbitals are given by products of polynomials and Gaussian functions that, as expected, peak at approximately half of the equilibrium interelectron distance.

The Journal of Chemical Physics•Aug 10, 2006•29 citations

Exotic Harmonium model: Exploring correlation effects of attractive Coulomb interaction

Nahid Sadat Riyahi, Mohammad Goli, Shant Shahbazian

Simple few-body systems often serve as theoretical laboratories across various branches of theoretical physics. A prominent example is the two-electron Harmonium model, which has been widely studied over the past three decades to gain insights into the nature of the electron-electron correlations in many-electron quantum systems. Building on our previous work [Riyahi et al., Phys. Rev. B 108, 245155 (2023)], we introduce an analogous model consisting of an electron and a positively charged particle (PCP) with variable mass, interacting via Coulomb forces while confined by external harmonic potentials. Termed the exotic Harmonium model, this provides insights into the electron-PCP correlations, a cornerstone of the emerging field of the ab initio study of multi-component many-body quantum systems. Through a systematic exploration of the parameter space and numerical solutions of the corresponding Schrödinger equation, we identify two extreme regimes: the atom-like and the particle-in-trap-like behavior. The electron-PCP correlation dominates in the atom-like regime, significantly influencing physical observables, while its role diminishes in the particle-in-trap-like limit. Between these two extremes lies a complex intermediate regime that challenges qualitative interpretation. Overall, the exotic Harmonium model offers a powerful framework to unravel the electron-PCP correlations across diverse systems, spanning particles of varying masses and conditions, from ambient to high-pressure environments.

arXiv:2505.10452•The Journal of Chemical Physics•Aug 15, 2025•1 citation

Entanglement in N-harmonium: bosons and fermions

C L Benavides-Riveros, I V Toranzo, J S Dehesa

The ground-state entanglement of a single particle of the N-harmonium system (i.e., a completely-integrable model of $N$ particles where both the confinement and the two-particle interaction are harmonic) is shown to be analytically determined in terms of $N$ and the relative interaction strength. For bosons, we compute the von Neumann entropy of the one-body reduced density matrix by using the corresponding natural occupation numbers. There exists a critical number $N_c$ of particles so that below it, for positive values of the coupling constant, the entanglement grows when the number of particles is increasing; the opposite occurs for $N > N_c$. For fermions, we compute the one-body reduced density matrix for the closed-shell spinned case. In the strong coupling regime, the linear entropy of the system decreases when $N$ is growing. For fixed $N$, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Moreover, the spatial and spin contributions to the total entanglement are found to be of comparable size.

arXiv:1404.4447•Journal of Physics B Atomic Molecular and Optical Physics•Sep 24, 2014•35 citations

Note: The weak-correlation limit of the three-electron harmonium atom

Jerzy Ciosłowski, Eduard Matito

Asymptotic energy expressions for the weak-correlation limits of the two lowest energy states of the three-electron harmonium atom are obtained in closed forms. When combined with the known results for the strong-correlation limit, these expressions, which are correct through the second order of perturbation theory, yield robust Padé approximants that allow accurate estimation of energies in question for all magnitudes of the confinement strength.

The Journal of Chemical Physics•Mar 17, 2011•13 citations

Hooke’s law correlation in two-electron systems

Pierre‐François Loos

We study the properties of the Hooke's law correlation energy (${E}_{\mathrm{c}}$), defined as the correlation energy when two electrons interact via a harmonic potential in a $D$-dimensional space. More precisely, we investigate the ${}^{1}S$ ground-state properties of two model systems: the Moshinsky atom (in which the electrons move in a quadratic potential) and the spherium model (in which they move on the surface of a sphere). A comparison with their Coulombic counterparts is made that highlights the main differences of the ${E}_{\mathrm{c}}$ in both the weakly and strongly correlated limits. Moreover, we show that the Schr\"odinger equation of the spherium model is exactly solvable for two values of the dimension ($D=1 \text{and} 3$) and that the exact wave function is based on Mathieu functions.

arXiv:1002.3396•Physical Review A•Mar 10, 2010•34 citations

At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.

Explicit asymptotic expressions are derived for the reduced one-electron density matrices (the 1-matrices) of strongly correlated two- and three-electron harmonium atoms in the ground and first excited states. These expressions, which are valid at the limit of small confinement strength ω, yield electron densities and kinetic energies in agreement with the published values. In addition, they reveal the ω(5/6) asymptotic scaling of the exchange components of the electron-electron repulsion energies that differs from the ω(2/3) scaling of their Coulomb and correlation counterparts. The natural orbitals of the totally symmetric ground state of the two-electron harmonium atom are found to possess collective occupancies that follow a mixed power/Gaussian dependence on the angular momentum in variance with the simple power-law prediction of Hill's asymptotics. Providing rigorous constraints on energies as functionals of 1-matrices, these results are expected to facilitate development of approximate implementations of the density matrix functional theory and ensure their proper description of strongly correlated systems.

Explicit asymptotic expressions for natural orbitals and their occupancies are derived for the harmonium atom at the strong-correlation limit at which the confinement strength omega tends to zero. Unlike in systems with moderate correlation effects, the occupancies at the omega-->0 limit (derived from occupation amplitudes with alternating sign patterns) are vanishingly small and asymptotically independent of the angular momentum, forming a geometric progression with the scale factor proportional to omega(1/3) and the common ratio of ca. 0.0186. The radial components of the natural orbitals are given by products of polynomials and Gaussian functions that, as expected, peak at approximately half of the equilibrium interelectron distance.

Simple few-body systems often serve as theoretical laboratories across various branches of theoretical physics. A prominent example is the two-electron Harmonium model, which has been widely studied over the past three decades to gain insights into the nature of the electron-electron correlations in many-electron quantum systems. Building on our previous work [Riyahi et al., Phys. Rev. B 108, 245155 (2023)], we introduce an analogous model consisting of an electron and a positively charged particle (PCP) with variable mass, interacting via Coulomb forces while confined by external harmonic potentials. Termed the exotic Harmonium model, this provides insights into the electron-PCP correlations, a cornerstone of the emerging field of the ab initio study of multi-component many-body quantum systems. Through a systematic exploration of the parameter space and numerical solutions of the corresponding Schrödinger equation, we identify two extreme regimes: the atom-like and the particle-in-trap-like behavior. The electron-PCP correlation dominates in the atom-like regime, significantly influencing physical observables, while its role diminishes in the particle-in-trap-like limit. Between these two extremes lies a complex intermediate regime that challenges qualitative interpretation. Overall, the exotic Harmonium model offers a powerful framework to unravel the electron-PCP correlations across diverse systems, spanning particles of varying masses and conditions, from ambient to high-pressure environments.

The ground-state entanglement of a single particle of the N-harmonium system (i.e., a completely-integrable model of

N

particles where both the confinement and the two-particle interaction are harmonic) is shown to be analytically determined in terms of

N

and the relative interaction strength. For bosons, we compute the von Neumann entropy of the one-body reduced density matrix by using the corresponding natural occupation numbers. There exists a critical number

N_c

of particles so that below it, for positive values of the coupling constant, the entanglement grows when the number of particles is increasing; the opposite occurs for

N > N_c

. For fermions, we compute the one-body reduced density matrix for the closed-shell spinned case. In the strong coupling regime, the linear entropy of the system decreases when

N

is growing. For fixed

N

, the entanglement is found (a) to decrease (increase) for negatively (positively) increasing values of the coupling constant, and (b) to grow when the energy is increasing. Moreover, the spatial and spin contributions to the total entanglement are found to be of comparable size.

We study the properties of the Hooke's law correlation energy (

{E}_{\mathrm{c}}

), defined as the correlation energy when two electrons interact via a harmonic potential in a

D

-dimensional space. More precisely, we investigate the

{}^{1}S

ground-state properties of two model systems: the Moshinsky atom (in which the electrons move in a quadratic potential) and the spherium model (in which they move on the surface of a sphere). A comparison with their Coulombic counterparts is made that highlights the main differences of the

{E}_{\mathrm{c}}

in both the weakly and strongly correlated limits. Moreover, we show that the Schr\"odinger equation of the spherium model is exactly solvable for two values of the dimension (

D=1 \text{and} 3

) and that the exact wave function is based on Mathieu functions.