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

The Journal of Chemical Physics

The interaction potential of a symmetric helium trimer

Aniket Bhattacharya, James B. Anderson

Published Jun 15, 1994• 19 citations
Publisher

Abstract

The interaction potential of three helium atoms in symmetric linear configurations is calculated under the Born–Oppenheimer approximation using an ‘‘exact’’ Green’s function quantum Monte Carlo method. The results are compared with those of ‘‘exact’’ quantum Monte Carlo calculations for the helium dimer and the Axilrod–Teller triple-dipole energy expression. For nearest neighbor separations of 5.6 bohrs, the equilibrium distance for a helium dimer, the corresponding trimer energy is −21.5±1.9 K, compared to −22.3±0.2 K, calculated for pairwise additive behavior. For all internuclear separations in the neighborhood of the van der Waals well (∼5.6 bohrs), the nonadditive contribution is found to be very small.

DOI
10.1063/1.467258

Related papers

Full-configuration-interaction calculation of three-body nonadditive contribution to helium interaction potential

Wojciech Cencek, Konrad Patkowski, Krzysztof Szalewicz

The three-body nonadditive interaction energy between helium atoms was calculated at 253 trimer configurations using the full-configuration-interaction (FCI) method. The analytic potential fitted to these energies is the best current representation of the three-body nonadditive interactions between helium atoms. At the equilateral triangle configuration with R=5.6 bohr, near the minimum of the total potential, the nonadditive three-body energy calculated at the FCI level amounts to -88.5 mK, compared to -98.5 mK at the coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] level. The uncertainty of the former result resulting from basis set incompleteness is estimated to be 1.5 mK. The relative uncertainty of our present complete three-body fit, including the uncertainties resulting from the fitting procedure, is estimated at 2%, a fivefold improvement over the previous best potential. Overall, the FCI contribution beyond CCSD(T) is rather important, being of the same order of magnitude as the uncertainty of the sum of two-body interactions. The inclusion of this contribution makes uncertainties of the total trimer interaction energies dominated by the uncertainties of the two-body component.

The Journal of Chemical Physics•Aug 11, 2009•60 citations

Structure and energetics of small helium clusters: Quantum simulations using a recent perturbational pair potential

Marius Lewerenz

Abstract Small helium clusters Hen (n = 2, 3, ... , 10) in molecular beams produced by supersonic expansion of gaseous helium have been observed in several different experiments. Among these clusters the dimer, which has a single very slightly bound state, and the trimer, which could have many bound states, have especially unusual properties which QMC calculations can explain in detail. In this paper, diffusion QMC calculations were used to determine ground state properties for clusters (n = 2, ... , 10) using two accurate, but slightly different pairwise-additive potential energy expressions, a perturbation-based potential by Tang et al.a and an empirical potential (HFD-B) proposed by Aziz et al.b Trial functions for importance sampling were based on those determined earlier by Rick et alc. Expectation values and probability distributions were determined by the descendent-weighting technique to give results independent of the trial functions. As expected, a single bound state of the dimer was found, and its average internuclear distance was predicted to be 51.5 A. Energies for clusters with n = 2, 7 with the HFD-B potential were in good agreement with those of earlier calculations. Comparisons of the results for the two different potential energy expressions showed energies differing by a few percent. Tests of possible nonadditive interactions showed these to be far less important in their effects on predicted energies.

•Jun 18, 2007•3 citations

Elusive structure of helium trimers

Petar Stipanović, Leandra Vranješ Markić, Jordi Boronat

Abstract Over the years many He–He interaction potentials have been developed, some very sophisticated, including various corrections beyond the Born–Oppenheimer approximation. Most of them were used to predict properties of helium dimers and trimers, examples of exotic quantum states, whose experimental study proved to be very challenging. Recently, detailed structural properties of helium trimers were measured for the first time, allowing a comparison with theoretical predictions and possibly enabling the evaluation of different interaction potentials. The comparisons already made included adjusting the maxima of both theoretical and experimental correlation functions to one, so the overall agreement between theory and experiment appeared satisfactory. However, no attempt was made to evaluate the quality of the interaction potentials used in the calculations. In this work, we calculate the experimentally measured correlation functions using both new and old potentials, compare them with experimental data and rank the potentials. We use diffusion Monte Carlo simulations at T = 0, which give within statistical noise exact results of the ground state properties. All models predict both trimers 4 He 3 and 4 He He to be in a quantum halo state.

arXiv:1607.05872•Journal of Physics B Atomic Molecular and Optical Physics•Aug 23, 2016•19 citations

Bounds for the scattering length of spin-polarized helium from high-accuracy electronic structure calculations

Michał Przybytek, Bogumił Jeziorski

We developed a series of correlation-consistent, polarized multiple zeta basis sets optimized specifically for the energy of the 2 3S state of helium atom. These basis sets were subsequently augmented with diffuse functions optimized for the van der Waals constants C6 through C14 which determine the asymptotic behavior of the second-order dispersion interaction between 2 3S helium atoms at large interatomic separation R. The resulting bases were applied to compute the Born-Oppenheimer (BO) potential for the lowest 5Sigmag+ state of the helium dimer. The coupled cluster and the full configuration-interaction techniques were employed to account for the electron correlation effects. The cardinal number extrapolation technique was used to obtain the complete-basis-set limit V(R) for the interaction potential and to find its lower VL(R) and upper VU(R) bounds. The resulting potentials were fitted to an analytical function containing accurate van der Waals constants C6 through C12 (including C11). We found that the complete-basis-set BO potential has a well depth De=1048.24+/-0.36 cm-1. The highest rotationless vibrational level is bound by D14=90.2+/-4.7 MHz, much stronger than the previous most accurate estimation of 15.2 MHz. The error bounds for De and D14 were obtained using the VL(R) and VU(R) potentials. The S-wave scattering length computed using the VL(R), V(R), and VU(R) potentials (assuming atomic masses) is aL=7.41 nm, a=7.54 nm, and aU=7.69 nm, respectively. We also computed the adiabatic, relativistic, and quantum electrodynamics (QED) corrections to the BO potential. When these corrections are taken into account the values of D14 and of a (both computed assuming nuclear masses) are 87.4+/-6.7 MHz and 7.64+/-0.20 nm; the error bounds reflect now also the uncertainty of the included adiabatic, relativistic, and QED corrections. The value of the scattering length resulting from our investigation lies outside the error bounds of all experimental determinations based on the properties of Bose-Einstein condensate of spin-polarized helium atoms.

The Journal of Chemical Physics•Sep 30, 2005•53 citations

Green’s function Monte Carlo calculation for the ground state of helium trimers

F. Cabral, Malvin H. Kalos

The ground state energy of weakly bound boson trimers interacting via Lennard-Jones (12,6) pair potentials is calculated using a Monte Carlo Green’s Function Method. Threshold coupling constants for self binding are obtained by extrapolation to zero binding.

The Journal of Chemical Physics•Feb 1, 1981•0 citations

Quantum Monte Carlo calculations of the potential energy curve of the helium dimer

Ryan Springall, Manolo C. Per, Salvy P. Russo +1 more

We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.

The Journal of Chemical Physics•Mar 19, 2008•17 citations
The interaction potential of three helium atoms in symmetric linear configurations is calculated under the Born–Oppenheimer approximation using an ‘‘exact’’ Green’s function quantum Monte Carlo method. The results are compared with those of ‘‘exact’’ quantum Monte Carlo calculations for the helium dimer and the Axilrod–Teller triple-dipole energy expression. For nearest neighbor separations of 5.6 bohrs, the equilibrium distance for a helium dimer, the corresponding trimer energy is −21.5±1.9 K, compared to −22.3±0.2 K, calculated for pairwise additive behavior. For all internuclear separations in the neighborhood of the van der Waals well (∼5.6 bohrs), the nonadditive contribution is found to be very small.
The three-body nonadditive interaction energy between helium atoms was calculated at 253 trimer configurations using the full-configuration-interaction (FCI) method. The analytic potential fitted to these energies is the best current representation of the three-body nonadditive interactions between helium atoms. At the equilateral triangle configuration with R=5.6 bohr, near the minimum of the total potential, the nonadditive three-body energy calculated at the FCI level amounts to -88.5 mK, compared to -98.5 mK at the coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] level. The uncertainty of the former result resulting from basis set incompleteness is estimated to be 1.5 mK. The relative uncertainty of our present complete three-body fit, including the uncertainties resulting from the fitting procedure, is estimated at 2%, a fivefold improvement over the previous best potential. Overall, the FCI contribution beyond CCSD(T) is rather important, being of the same order of magnitude as the uncertainty of the sum of two-body interactions. The inclusion of this contribution makes uncertainties of the total trimer interaction energies dominated by the uncertainties of the two-body component.
Abstract Small helium clusters Hen (n = 2, 3, ... , 10) in molecular beams produced by supersonic expansion of gaseous helium have been observed in several different experiments. Among these clusters the dimer, which has a single very slightly bound state, and the trimer, which could have many bound states, have especially unusual properties which QMC calculations can explain in detail. In this paper, diffusion QMC calculations were used to determine ground state properties for clusters (n = 2, ... , 10) using two accurate, but slightly different pairwise-additive potential energy expressions, a perturbation-based potential by Tang et al.a and an empirical potential (HFD-B) proposed by Aziz et al.b Trial functions for importance sampling were based on those determined earlier by Rick et alc. Expectation values and probability distributions were determined by the descendent-weighting technique to give results independent of the trial functions. As expected, a single bound state of the dimer was found, and its average internuclear distance was predicted to be 51.5 A. Energies for clusters with n = 2, 7 with the HFD-B potential were in good agreement with those of earlier calculations. Comparisons of the results for the two different potential energy expressions showed energies differing by a few percent. Tests of possible nonadditive interactions showed these to be far less important in their effects on predicted energies.
Abstract Over the years many He–He interaction potentials have been developed, some very sophisticated, including various corrections beyond the Born–Oppenheimer approximation. Most of them were used to predict properties of helium dimers and trimers, examples of exotic quantum states, whose experimental study proved to be very challenging. Recently, detailed structural properties of helium trimers were measured for the first time, allowing a comparison with theoretical predictions and possibly enabling the evaluation of different interaction potentials. The comparisons already made included adjusting the maxima of both theoretical and experimental correlation functions to one, so the overall agreement between theory and experiment appeared satisfactory. However, no attempt was made to evaluate the quality of the interaction potentials used in the calculations. In this work, we calculate the experimentally measured correlation functions using both new and old potentials, compare them with experimental data and rank the potentials. We use diffusion Monte Carlo simulations at T = 0, which give within statistical noise exact results of the ground state properties. All models predict both trimers 4 He 3 and 4 He He to be in a quantum halo state.
We developed a series of correlation-consistent, polarized multiple zeta basis sets optimized specifically for the energy of the 2 3S state of helium atom. These basis sets were subsequently augmented with diffuse functions optimized for the van der Waals constants C6 through C14 which determine the asymptotic behavior of the second-order dispersion interaction between 2 3S helium atoms at large interatomic separation R. The resulting bases were applied to compute the Born-Oppenheimer (BO) potential for the lowest 5Sigmag+ state of the helium dimer. The coupled cluster and the full configuration-interaction techniques were employed to account for the electron correlation effects. The cardinal number extrapolation technique was used to obtain the complete-basis-set limit V(R) for the interaction potential and to find its lower VL(R) and upper VU(R) bounds. The resulting potentials were fitted to an analytical function containing accurate van der Waals constants C6 through C12 (including C11). We found that the complete-basis-set BO potential has a well depth De=1048.24+/-0.36 cm-1. The highest rotationless vibrational level is bound by D14=90.2+/-4.7 MHz, much stronger than the previous most accurate estimation of 15.2 MHz. The error bounds for De and D14 were obtained using the VL(R) and VU(R) potentials. The S-wave scattering length computed using the VL(R), V(R), and VU(R) potentials (assuming atomic masses) is aL=7.41 nm, a=7.54 nm, and aU=7.69 nm, respectively. We also computed the adiabatic, relativistic, and quantum electrodynamics (QED) corrections to the BO potential. When these corrections are taken into account the values of D14 and of a (both computed assuming nuclear masses) are 87.4+/-6.7 MHz and 7.64+/-0.20 nm; the error bounds reflect now also the uncertainty of the included adiabatic, relativistic, and QED corrections. The value of the scattering length resulting from our investigation lies outside the error bounds of all experimental determinations based on the properties of Bose-Einstein condensate of spin-polarized helium atoms.
We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.