Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2505.10452 · The Journal of Chemical Physics

Exotic Harmonium model: Exploring correlation effects of attractive Coulomb interaction

Nahid Sadat Riyahi, Mohammad Goli, Shant Shahbazian

Published Aug 15, 2025• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.1063/5.0284258

Related papers

Development of Exotic Harmonium Model to Investigate Electron-Positively Charged Particle Correlation in Two-Component Quantum Systems

Nahid Sadat Riyahi

The problem of calculating the electron-positively charged particle correlation energy poses a challenge in the field of quantum chemistry beyond the adiabatic approximation. In this study, a toy model called Exotic Harmonium is developed to enhance our understanding of this type of correlation. Since the analytical methods to solve the eigenvalue equation for Exotic Harmonium lead to the excited rather than the ground state of the system, we employ the variation method in this study to obtain an approximate ground state wave function. By considering the asymptotic behavior of the wave function, we derive a compact but highly accurate variational wave function. Using this wave function, we are able to determine various properties of the system, including energy and its components, as well as single-particle densities. Additionally, we extend key concepts and quantities specifically tailored to study electron correlation in the Harmonium model to the Exotic Harmonium model. For instance, the "correlation hill" concept effectively highlights the distinct nature of electron-positively charged particle correlation compared to electron-electron correlation. Furthermore, we utilize Exotic Harmonium to assess the accuracy of five electron-positively charged particle correlation functionals developed within the context of the two-component density functional theory. Through a regression process, we obtained a compact yet precise analytical form for wave function which explicitly depends on the two crucial variables: oscillator field frequency and positively charged particle mass. This analytical form provides valuable insights into the behavior of the system. Also, the Exotic Harmonium model enables the investigation of electron-positively charged particle correlation in a vast range of particle masses.

arXiv:2504.18118•Apr 25, 2025•0 citations

Contactium: A strongly correlated model system

Jerzy Ciosłowski, Berthold‐Georg Englert, Martin-Isbjörn Trappe +1 more

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.

arXiv:2303.14982•The Journal of Chemical Physics•May 12, 2023•3 citations

Potential harmonics expansion method for trapped interacting bosons: Inclusion of two-body correlation

Tapan Kumar Das, Barnali Chakrabarti

We study a system of $A$ identical interacting bosons trapped by an external field by solving ab initio the many-body Schr\"odinger equation. A complete solution by using, for example, the traditional hyperspherical harmonics (HH) basis develops serious practical problems due to the large degeneracy of HH basis. Symmetrization of the wave function, calculation of the matrix elements, etc., become an immensely formidable task as $A$ increases. Instead of the HH basis, here we use a new basis, called ``potential harmonics'' (PH) basis, which is a subset of HH basis. We assume that the contribution to the orbital and grand orbital [in $3(A\ensuremath{-}1)$-dimensional space of the reduced motion] quantum numbers comes only from the interacting pair. This implies inclusion of two-body correlations only and disregard of all higher-body correlations. Such an assumption is ideally suited for the Bose-Einstein condensate (BEC), which is required, for experimental realization of BEC, to be extremely dilute. Hence three and higher-body collisions are almost totally absent. Unlike the $(3A\ensuremath{-}4)$ hyperspherical variables in HH basis, the PH basis involves only three active variables, corresponding to three quantum numbers---the orbital $l$, azimuthal $m$, and the grand orbital $2K+l$ quantum numbers for any arbitrary $A$. It drastically reduces the number of coupled equations and calculation of the potential matrix becomes tremendously simplified, as it involves integrals over only three variables for any $A$. One can easily incorporate realistic atom-atom interactions in a straightforward manner. We study the ground and excited state properties of the condensate for both attractive and repulsive interactions for various particle number. The ground state properties are compared with those calculated from the Gross-Pitaevskii equation. We notice that our many-body results converge towards the mean field results as the particle number increases.

arXiv:quant-ph/0408061•Physical Review A•Dec 1, 2004•54 citations

Relativistic effects on entangled single-electron traps

Marko Toroš, Patrick Andriolo, Martine Schut +2 more

The manipulation of individual charged particles has been deeply explored in physics’s theoretical and experimental domains during the past decades. It is the pillar of several existing devices used for metrology and sensing and is a promising platform for realizing future technologies, such as quantum computers. It is also known that in the relativistic regime, interactions between charged particles become affected by post-Coulombian corrections, with the dominant couplings encoded in the . The Darwin term has been extensively studied in atomic physics, where the interaction range is confined to the sub-angstrom scale. Still, there is a lack of understanding about whether (and when) Darwin’s contributions are relevant at larger scales. In this paper, we explore the effects of these corrections in a system of two harmonically trapped electrons, where we look into the behavior of quantum entanglement present in the static and dynamical regimes. We explore the parameter space of the developed model and seek frequencies, distances, and squeezing parameters for which relativistic effects become relevant for the generation of entanglement. Published by the American Physical Society 2024

arXiv:2406.17848•Physical review. D/Physical review. D.•Sep 16, 2024•7 citations

Trapped electrons in the quantum degenerate regime

Ronen M. Kroeze, Jom Luiten, S. J. J. M. F. Kokkelmans

A full strength Coulomb interaction between trapped electrons can be felt only in absence of a neutralizing background. In order to study quantum degenerate electrons without such a background, an external trap is needed to compensate for the strong electronic repulsion. As a basic model for such a system, we study a trapped electron pair in a harmonic trap with an explicit inclusion of its Coulomb interaction. We find the eigenenergy of the ground state, confirming earlier work in the context of harmonium. We extend this to a complete set of properly scaled energies for any value of the trapping strength, including the excited states. The problem is solved either numerically or by making harmonic approximations to the potential. As function of the trapping strength a crossover can be made from the strongly to the weakly-coupled regime, and we show that in both regimes perturbative methods based on a pair-wise electron description would be effective for a many-particle trapped electron system, which resembles a Wigner crystal in the ground state of the strongly coupled limit.

arXiv:1508.00365•Aug 3, 2015•0 citations

Many-body physics in small systems: Observing the onset and saturation of correlation in linear atomic chains

Emily Townsend, Tomáš Neuman, Alex Debrecht +2 more

The exact study of small systems can guide us toward relevant measures for extracting information about many-body physics as we move to larger and more complex systems capable of quantum information processing or quantum analog simulation. We use exact diagonalization to study many electrons in short 1-D atom chains represented by long-range extended Hubbard-like models. We introduce a novel measure, the Single-Particle Excitation Content (SPEC) of an eigenstate and show that the dependence of SPEC on eigenstate number reveals the nature of the ground state (ordered phases), and the onset and saturation of correlation between the electrons as Coulomb interaction strength increases. We use this SPEC behavior to identify five regimes as interaction is increased: a non-interacting single-particle regime, a regime of perturbative Coulomb interaction in which the SPEC is a nearly universal function of eigenstate number, the onset and saturation of correlation, a regime of fully correlated states in which hopping is a perturbation and SPEC is a different universal function of state number, and the regime of no hopping. In particular, the behavior of the SPEC shows that when electron-electron correlation plays a minor role, all of the lowest energy eigenstates are made up primarily of single-particle excitations of the ground state, and as the Coulomb interaction increases, the lowest energy eigenstates increasingly contain many-particle excitations. In addition, the SPEC highlights a fundamental, distinct difference between a non-interacting system and one with minute, very weak interactions. While SPEC is a quantity that can be calculated for small exactly diagonalizable systems, it guides our intuition for larger systems, suggesting the nature of excitations and their distribution in the spectrum. Thus, this function, like correlation functions or order parameters, provides us with a window of intuition about the behavior of a physical system.

arXiv:2011.07113•Physical review. B./Physical review. B•May 20, 2021•5 citations
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 problem of calculating the electron-positively charged particle correlation energy poses a challenge in the field of quantum chemistry beyond the adiabatic approximation. In this study, a toy model called Exotic Harmonium is developed to enhance our understanding of this type of correlation. Since the analytical methods to solve the eigenvalue equation for Exotic Harmonium lead to the excited rather than the ground state of the system, we employ the variation method in this study to obtain an approximate ground state wave function. By considering the asymptotic behavior of the wave function, we derive a compact but highly accurate variational wave function. Using this wave function, we are able to determine various properties of the system, including energy and its components, as well as single-particle densities. Additionally, we extend key concepts and quantities specifically tailored to study electron correlation in the Harmonium model to the Exotic Harmonium model. For instance, the "correlation hill" concept effectively highlights the distinct nature of electron-positively charged particle correlation compared to electron-electron correlation. Furthermore, we utilize Exotic Harmonium to assess the accuracy of five electron-positively charged particle correlation functionals developed within the context of the two-component density functional theory. Through a regression process, we obtained a compact yet precise analytical form for wave function which explicitly depends on the two crucial variables: oscillator field frequency and positively charged particle mass. This analytical form provides valuable insights into the behavior of the system. Also, the Exotic Harmonium model enables the investigation of electron-positively charged particle correlation in a vast range of particle masses.
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.
We study a system of AAA identical interacting bosons trapped by an external field by solving ab initio the many-body Schr\"odinger equation. A complete solution by using, for example, the traditional hyperspherical harmonics (HH) basis develops serious practical problems due to the large degeneracy of HH basis. Symmetrization of the wave function, calculation of the matrix elements, etc., become an immensely formidable task as AAA increases. Instead of the HH basis, here we use a new basis, called ``potential harmonics'' (PH) basis, which is a subset of HH basis. We assume that the contribution to the orbital and grand orbital [in $3(A\ensuremath{-}1)$-dimensional space of the reduced motion] quantum numbers comes only from the interacting pair. This implies inclusion of two-body correlations only and disregard of all higher-body correlations. Such an assumption is ideally suited for the Bose-Einstein condensate (BEC), which is required, for experimental realization of BEC, to be extremely dilute. Hence three and higher-body collisions are almost totally absent. Unlike the $(3A\ensuremath{-}4)$ hyperspherical variables in HH basis, the PH basis involves only three active variables, corresponding to three quantum numbers---the orbital lll, azimuthal mmm, and the grand orbital 2K+l2K+l2K+l quantum numbers for any arbitrary AAA. It drastically reduces the number of coupled equations and calculation of the potential matrix becomes tremendously simplified, as it involves integrals over only three variables for any AAA. One can easily incorporate realistic atom-atom interactions in a straightforward manner. We study the ground and excited state properties of the condensate for both attractive and repulsive interactions for various particle number. The ground state properties are compared with those calculated from the Gross-Pitaevskii equation. We notice that our many-body results converge towards the mean field results as the particle number increases.
The manipulation of individual charged particles has been deeply explored in physics’s theoretical and experimental domains during the past decades. It is the pillar of several existing devices used for metrology and sensing and is a promising platform for realizing future technologies, such as quantum computers. It is also known that in the relativistic regime, interactions between charged particles become affected by post-Coulombian corrections, with the dominant couplings encoded in the . The Darwin term has been extensively studied in atomic physics, where the interaction range is confined to the sub-angstrom scale. Still, there is a lack of understanding about whether (and when) Darwin’s contributions are relevant at larger scales. In this paper, we explore the effects of these corrections in a system of two harmonically trapped electrons, where we look into the behavior of quantum entanglement present in the static and dynamical regimes. We explore the parameter space of the developed model and seek frequencies, distances, and squeezing parameters for which relativistic effects become relevant for the generation of entanglement. Published by the American Physical Society 2024
A full strength Coulomb interaction between trapped electrons can be felt only in absence of a neutralizing background. In order to study quantum degenerate electrons without such a background, an external trap is needed to compensate for the strong electronic repulsion. As a basic model for such a system, we study a trapped electron pair in a harmonic trap with an explicit inclusion of its Coulomb interaction. We find the eigenenergy of the ground state, confirming earlier work in the context of harmonium. We extend this to a complete set of properly scaled energies for any value of the trapping strength, including the excited states. The problem is solved either numerically or by making harmonic approximations to the potential. As function of the trapping strength a crossover can be made from the strongly to the weakly-coupled regime, and we show that in both regimes perturbative methods based on a pair-wise electron description would be effective for a many-particle trapped electron system, which resembles a Wigner crystal in the ground state of the strongly coupled limit.
The exact study of small systems can guide us toward relevant measures for extracting information about many-body physics as we move to larger and more complex systems capable of quantum information processing or quantum analog simulation. We use exact diagonalization to study many electrons in short 1-D atom chains represented by long-range extended Hubbard-like models. We introduce a novel measure, the Single-Particle Excitation Content (SPEC) of an eigenstate and show that the dependence of SPEC on eigenstate number reveals the nature of the ground state (ordered phases), and the onset and saturation of correlation between the electrons as Coulomb interaction strength increases. We use this SPEC behavior to identify five regimes as interaction is increased: a non-interacting single-particle regime, a regime of perturbative Coulomb interaction in which the SPEC is a nearly universal function of eigenstate number, the onset and saturation of correlation, a regime of fully correlated states in which hopping is a perturbation and SPEC is a different universal function of state number, and the regime of no hopping. In particular, the behavior of the SPEC shows that when electron-electron correlation plays a minor role, all of the lowest energy eigenstates are made up primarily of single-particle excitations of the ground state, and as the Coulomb interaction increases, the lowest energy eigenstates increasingly contain many-particle excitations. In addition, the SPEC highlights a fundamental, distinct difference between a non-interacting system and one with minute, very weak interactions. While SPEC is a quantity that can be calculated for small exactly diagonalizable systems, it guides our intuition for larger systems, suggesting the nature of excitations and their distribution in the spectrum. Thus, this function, like correlation functions or order parameters, provides us with a window of intuition about the behavior of a physical system.