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 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 Dirac bubble potential, previously used to model the helium dimer, is applied to the trimer. It is shown that the quantum mechanical kinetic energy operator for a three-body system contains terms over and above the analogues of classical pairwise contributions. The additional terms, the ‘Borromean couplings’, are responsible for a dramatic increase in the binding energy of the trimer. For example, 4He3 is bounded by two orders of magnitude more strongly than 4He2, 279 mK versus 1.31 mK, respectively, according to our calculations. Moreover, the trimer is considerably more compact than the dimer, with (r12) decreasing to 9.01 Å from 51.9 Å.
The
4He3 system is studied using the adiabatic hyperspherical representation. We adopt the current state-of-the-art helium interaction potential including retardation and the nonadditive three-body term to calculate all low-energy properties of the triatomic
4He system. The bound state energies of the
4He trimer are computed as well as the
4He+4He2 elastic scattering cross sections, the three-body recombination, and collision induced dissociation rates at finite temperatures. We also treat the system that consists of two
4He and one
3He atoms, and compute the spectrum of the isotopic trimer
4He23He, the
3He+4He2 elastic scattering cross sections, the rates for three-body recombination, and the collision induced dissociation rate at finite temperatures. The effects of retardation and the nonadditive three-body term are investigated. Retardation is found to be significant in some cases, while the three-body term plays only a minor role for these systems.
We theoretically investigate long-range interactions between an excited
L-state He atom and two identical
S-state He atoms for the cases of the three atoms all in spin-singlet states or all in spin-triplet states, denoted by He(${n}_{0}^{\ensuremath{\lambda}}S$)-He(${n}_{0}^{\ensuremath{\lambda}}S$)-He(${n}_{0}^{\ensuremath{'}\ensuremath{\lambda}}L$), with
n0 and ${n}_{0}^{\ensuremath{'}}$ principal quantum numbers, $\ensuremath{\lambda}=1$ or 3 the spin multiplicity, and
L the orbital angular momentum of a He atom. Using degenerate perturbation theory for the energies up to second-order, we evaluate the coefficients
C3 of the first-order dipolar interactions and the coefficients
C6 and
C8 of the second-order additive and nonadditive interactions. Both the dipolar and dispersion interaction coefficients, for these three-body degenerate systems, show dependences on the geometrical configurations of the three atoms. The nonadditive interactions start to appear in second-order. To demonstrate the results and for applications, the obtained coefficients
Cn are evaluated with highly accurate variationally generated nonrelativistic wave functions in Hylleraas coordinates for $\text{He}(1{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(1{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{1}S)$, $\text{He}(1{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(1{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{1}P)$, $\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{1}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{1}P)$, and $\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{3}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{3}S)\text{\ensuremath{-}}\text{He}(2{\phantom{\rule{0.16em}{0ex}}}^{3}P)$. The calculations are given for three like nuclei for the cases of hypothetical infinite mass He nuclei, and of real finite mass
4He or
3He nuclei. The special cases of the three atoms in equilateral triangle configurations are explored in detail, and for the cases in which one of the atoms is in a
P state, we also present results for the atoms in an isosceles right triangle configuration or in an equally spaced collinear configuration. The results can be applied to construct potential energy surfaces for three helium atom systems.
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].
We present detailed calculations of bound and scattering states of dimers and trimers of He to produce highly accurate data and to test a non-relativistic three-body code currently in development for public distribution. For these systems, uncertainties and inaccuracies in the fundamental constants frequently used in published works can substantially exceed numerical errors. Our benchmark calculations include specific estimates of the numerical accuracy of the calculations and also explore sensitivity to fundamental constants and their uncertainties. The use of an inexact coupling constant in the previous calculations leads to 0.08% error for the ground state energy, 0.3% error for the excited state energy and up to 0.15% error for the atom–dimer scattering length in the system of three 4 He atoms. The corresponding errors for the asymmetric 4 He 2 3 He system are 0.3% for the bound state energy and 0.03% for the atom–dimer scattering length.