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:2504.18118

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

Nahid Sadat Riyahi

Published Apr 25, 2025• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2504.18118

Related papers

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

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

Electron correlations. II. Ground-state results at low and metallic densities

R. F. Bishop, K.H. Lührmann

In this second work in a series of papers the coupled-cluster or $\mathrm{exp}(S)$ formalism is applied to the problem of ground-state correlations in the one-component Fermi plasma. The main approximation employed is the so-called SUB2 approximation for the two-body subsystem correlation operator ${S}_{2}$ which provides a measure of the two-particle---two-hole component in the true ground-state wave function, and in which the coupling to larger subsystems is neglected. The SUB2 equations for Fermi systems are brought into tractable form by a "state-averaging" procedure which we develop by analogy with the comparable Bose equations, and which is shown by comparison with earlier exact results to be accurate in the metallic density regime at about the 1% level. Our final results for metallic densities include the completely integrated and self-consistent effects of the terms which by themselves generate (i) the random-phase approximation (RPA) and its associated exact long-range screening effects, (ii) the extra random-phase approximation exchange terms necessary to keep the RPA explicitly antisymmetric, (iii) the self-consistent particle-particle ladders (LAD) that describe two-particle scattering within the many-body medium and which describe the exact short-range behavior, (iv) a class of particle-hole ladder terms, and (v) the self-consistent hole-potential terms. Particular attention is paid to the important effects caused by the interference at intermediate separations of the long-range RPA and the short-range LAD effects. By comparison with recent and essentially exact stochastic simulations of the many-body Schr\"odinger equation, our results are seen to be accurate to about 1% for metallic densities, and hence to provide what is probably the currently best microscopic description available for this system. We show also that even in the low-density limit the SUB2 approximation provides a good "translationally-invariant-solid" description of both charged Fermi and Bose systems in this exact Wigner crystal regime.

Physical review. B, Condensed matter•Nov 15, 1982•142 citations

Development and applications of density matrix functional theory of generalized Hubbard Models

Matthieu Saubanère

In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.

•Jan 25, 2012•0 citations

Putting a Temperature on Electron Correlation

Paul G. Graf, Florian Matz, Markus Penz +1 more

It is common wisdom that a better understanding of how to quantify electron correlation should lead to more efficient treatments of the correlation problem. Here, we not only substantiate this view but turn it into a concrete, systematic framework for ground states. Since prominent correlation measures such as the nonfreeness depend only on the one-particle reduced density matrix, modeling the correlation energy as a function of such a measure naturally places the framework within reduced density matrix functional theory, while keeping the cost at mean-field level. Strikingly, the variable conjugate to the measure emerges as a ``correlation temperature'', thereby connecting correlation measures, functional theory, and effective temperatures. The framework rationalizes methods such as i-DMFT, and we illustrate its key concepts analytically for the Hubbard dimer and numerically for small molecular systems. This brings a fundamental question within reach: how much of the correlation problem can be solved at mean-field cost?

cond-mat.str-elCross-listedarXiv:2610.22134•Oct 16, 2026•0 citations

The three-electron harmonium atom: The lowest-energy doublet and quadruplet states

Jerzy Ciosłowski, Krzysztof Strasburger, Eduard Matito

Calculations of sub-μhartree accuracy employing explicitly correlated Gaussian lobe functions produce comprehensive data on the energy E(ω), its components, and the one-electron properties of the two lowest-energy states of the three-electron harmonium atom. The energy computations at 19 values of the confinement strength ω ranging from 0.001 to 1000.0, used in conjunction with a recently proposed robust interpolation scheme, yield explicit approximants capable of estimating E(ω) and the potential energy of the harmonic confinement within a few tenths of μhartree for any ω ≥ 0.001, the respective errors for the kinetic energy and the potential energy of the electron-electron repulsion not exceeding 2 μhartrees. Thanks to the correct ω → 0 asymptotics incorporated into the approximants, comparable accuracy is expected for values of ω smaller than 0.001. Occupation numbers of the dominant natural spinorbitals and two different measures of electron correlation are also computed.

The Journal of Chemical Physics•May 21, 2012•30 citations
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.
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.
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.
In this second work in a series of papers the coupled-cluster or exp(S)\mathrm{exp}(S)exp(S) formalism is applied to the problem of ground-state correlations in the one-component Fermi plasma. The main approximation employed is the so-called SUB2 approximation for the two-body subsystem correlation operator S2{S}_{2}S2​ which provides a measure of the two-particle---two-hole component in the true ground-state wave function, and in which the coupling to larger subsystems is neglected. The SUB2 equations for Fermi systems are brought into tractable form by a "state-averaging" procedure which we develop by analogy with the comparable Bose equations, and which is shown by comparison with earlier exact results to be accurate in the metallic density regime at about the 1% level. Our final results for metallic densities include the completely integrated and self-consistent effects of the terms which by themselves generate (i) the random-phase approximation (RPA) and its associated exact long-range screening effects, (ii) the extra random-phase approximation exchange terms necessary to keep the RPA explicitly antisymmetric, (iii) the self-consistent particle-particle ladders (LAD) that describe two-particle scattering within the many-body medium and which describe the exact short-range behavior, (iv) a class of particle-hole ladder terms, and (v) the self-consistent hole-potential terms. Particular attention is paid to the important effects caused by the interference at intermediate separations of the long-range RPA and the short-range LAD effects. By comparison with recent and essentially exact stochastic simulations of the many-body Schr\"odinger equation, our results are seen to be accurate to about 1% for metallic densities, and hence to provide what is probably the currently best microscopic description available for this system. We show also that even in the low-density limit the SUB2 approximation provides a good "translationally-invariant-solid" description of both charged Fermi and Bose systems in this exact Wigner crystal regime.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
It is common wisdom that a better understanding of how to quantify electron correlation should lead to more efficient treatments of the correlation problem. Here, we not only substantiate this view but turn it into a concrete, systematic framework for ground states. Since prominent correlation measures such as the nonfreeness depend only on the one-particle reduced density matrix, modeling the correlation energy as a function of such a measure naturally places the framework within reduced density matrix functional theory, while keeping the cost at mean-field level. Strikingly, the variable conjugate to the measure emerges as a ``correlation temperature'', thereby connecting correlation measures, functional theory, and effective temperatures. The framework rationalizes methods such as i-DMFT, and we illustrate its key concepts analytically for the Hubbard dimer and numerically for small molecular systems. This brings a fundamental question within reach: how much of the correlation problem can be solved at mean-field cost?
Calculations of sub-μhartree accuracy employing explicitly correlated Gaussian lobe functions produce comprehensive data on the energy E(ω), its components, and the one-electron properties of the two lowest-energy states of the three-electron harmonium atom. The energy computations at 19 values of the confinement strength ω ranging from 0.001 to 1000.0, used in conjunction with a recently proposed robust interpolation scheme, yield explicit approximants capable of estimating E(ω) and the potential energy of the harmonic confinement within a few tenths of μhartree for any ω ≥ 0.001, the respective errors for the kinetic energy and the potential energy of the electron-electron repulsion not exceeding 2 μhartrees. Thanks to the correct ω → 0 asymptotics incorporated into the approximants, comparable accuracy is expected for values of ω smaller than 0.001. Occupation numbers of the dominant natural spinorbitals and two different measures of electron correlation are also computed.