The pair potential for helium is computed with accuracy improved by an order of magnitude relative to the best previous determination. For the well region, its uncertainties are now below 1 millikelvin. The main improvement is due to the use of explicitly correlated wave functions at the nonrelativistic Born-Oppenheimer (BO) level of theory. The diagonal BO and the relativistic corrections are obtained from large full configuration interaction calculations. Nonadiabatic perturbation theory is used to predict the properties of the halo state of the helium dimer. Its binding energy and the average value of the interatomic distance are found to be 138.9(5) neV and 47.13(8) Å. The binding energy agrees with its first experimental determination of 151.9(13.3) neV [Zeller et al., Proc. Natl. Acad. Sci. U.S.A. 113, 14651 (2016)PNASA60027-842410.1073/pnas.1610688113].
The pair potential for helium is computed with accuracy improved by an order of magnitude relative to the best previous determination. For the well region, its uncertainties are now below 1 millikelvin. The main improvement is due to the use of explicitly correlated wave functions at the nonrelativistic Born-Oppenheimer (BO) level of theory. The diagonal BO and the relativistic corrections are obtained from large full configuration interaction calculations. Nonadiabatic perturbation theory is used to predict the properties of the halo state of the helium dimer. Its binding energy and the average value of the interatomic distance are found to be 138.9(5) neV and 47.13(8) Å. The binding energy agrees with its first experimental determination of 151.9(13.3) neV [Zeller et al., Proc. Natl. Acad. Sci. U.S.A. 113, 14651 (2016)PNASA60027-842410.1073/pnas.1610688113].
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.
The precision of the photoassociation spectroscopy of the Yb dimer in degenerate gases is enough to improve the constraints on short-range gravitylike forces if the theoretical knowledge of the Born-Oppenheimer interatomic potential and non-Born-Oppenheimer interactions is refined [Borkowski et al., Sci. Rep. 9, 14807 (2019)]. The ground-state interaction potential of the ytterbium dimer is investigated at the exact two-component core-correlated coupled-cluster method with singles, doubles, and noniterative triples [CCSD(T)] level of ab initio theory in the complete-basis-set limit with extensive augmentation by diffuse functions. For the small basis set the comparison is made with the four-component relativistic finite-nuclei CCSD(T) calculations to identify the contraction of the dimer bond length as the main unrecoverable consequence of the scalar-relativistic approximation. The empirical constraint on the number of bound vibrational energy levels of the 174Yb2 dimer is accounted for by representing the global ab initio--based Born-Oppenheimer potential with the model semianalytical function containing the scale and shift parameters. The results support the previous evaluation of the Yb dimer potentials from the photoassociation spectroscopy data and provide an accurate and flexible reference for future refinement of the constraints on short-range gravitylike forces by ultracold atomic spectroscopy.
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.
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 To describe long-range behaviour of one particle removed from a few- or a many-body system, a hyperspherical cluster model has been developed. It has been applied to the ground and first excited states of helium drops with five, six, eight and ten atoms interacting via a two-body soft gaussian potential. Convergence of the hyperspherical cluster harmonics expansion is studied for binding energies, root-mean-squared radii and overlaps of the wave functions of two helium drops differing by one atom. It was shown that with increasing model space the functional form of such overlaps at large distances converges to the correct asymptotic behaviour. The asymptotic normalization coefficients that quantify the overlaps’ amplitudes in this region are calculated. It was also shown that in the first excited state one helium atom stays far apart from the rest forming a two-body molecule, or a halo. The probability of finding the halo atom in the classically-forbidden region of space depends on the definition of the latter and on the valence atom binding energy. The total norm of the overlap integrals, the spectroscopic factor, represents the number of partitions of a many-body state into a chosen state of the system with one particle removed. The spectroscopic factors have been calculated and their sum rules are discussed giving a further insight into the structure of helium drops.
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.