The adiabatic, relativistic, and quantum electrodynamics (QED) contributions to the pair potential of helium were computed, fitted separately, and applied, together with the nonrelativistic Born-Oppenheimer (BO) potential, in calculations of thermophysical properties of helium and of the properties of the helium dimer. An analysis of the convergence patterns of the calculations with increasing basis set sizes allowed us to estimate the uncertainties of the total interaction energy to be below 50 ppm for interatomic separations R smaller than 4 bohrs and for the distance R = 5.6 bohrs. For other separations, the relative uncertainties are up to an order of magnitude larger (and obviously still larger near R = 4.8 bohrs where the potential crosses zero) and are dominated by the uncertainties of the nonrelativistic BO component. These estimates also include the contributions from the neglected relativistic and QED terms proportional to the fourth and higher powers of the fine-structure constant α. To obtain such high accuracy, it was necessary to employ explicitly correlated Gaussian expansions containing up to 2400 terms for smaller R (all R in the case of a QED component) and optimized orbital bases up to the cardinal number X = 7 for larger R. Near-exact asymptotic constants were used to describe the large-R behavior of all components. The fitted potential, exhibiting the minimum of -10.996 ± 0.004 K at R = 5.608 0 ± 0.000 1 bohr, was used to determine properties of the very weakly bound (4)He(2) dimer and thermophysical properties of gaseous helium. It is shown that the Casimir-Polder retardation effect, increasing the dimer size by about 2 Å relative to the nonrelativistic BO value, is almost completely accounted for by the inclusion of the Breit-interaction and the Araki-Sucher contributions to the potential, of the order α(2) and α(3), respectively. The remaining retardation effect, of the order of α(4) and higher, is practically negligible for the bound state, but is important for the thermophysical properties of helium. Such properties computed from our potential have uncertainties that are generally significantly smaller (sometimes by nearly two orders of magnitude) than those of the most accurate measurements and can be used to establish new metrology standards based on properties of low-density helium.

The variational method complemented with the use of explicitly correlated Gaussian basis functions is one of the most powerful approaches currently used for calculating the properties of few-body systems. Despite its conceptual simplicity, the method offers great flexibility, high accuracy, and can be used to study diverse quantum systems, ranging from small atoms and molecules to light nuclei, hadrons, quantum dots, and Efimov systems. The basic theoretical foundations are discussed, recent advances in the applications of explicitly correlated Gaussians in physics and chemistry are reviewed, and the strengths and weaknesses of the explicitly correlated Gaussians approach are compared with other few-body techniques.

The dissociation energies from all rovibrational levels of H2 and D2 in the ground electronic state are calculated with high accuracy by including relativistic and quantum electrodynamics (QED) effects in the nonadiabatic treatment of the nuclear motion. For D2, the obtained energies have theoretical uncertainties of 0.001 cm(-1). For H2, similar uncertainties are for the lowest levels, while for the higher ones the uncertainty increases to 0.005 cm(-1). Very good agreement with recent high-resolution measurements of the rotational v = 0 levels of H2, including states with large angular momentum J, is achieved. This agreement would not have been possible without accurate evaluation of the relativistic and QED contributions and may be viewed as the first observation of the QED effects, mainly the electron self-energy, in a molecular spectrum. For several electric quadrupole transitions, we still observe certain disagreement with experimental results, which remains to be explained.

The dissociation energy of molecular hydrogen is determined theoretically with a careful estimation of error bars by including nonadiabatic, relativistic, and quantum electrodynamics (QED) corrections. The relativistic and QED corrections were obtained at the adiabatic level of theory by including all contributions of the order α 2 and α 3 as well as the major (one-loop) α 4 term, where α is the fine structure constant. The computed α 0 , α 2 , α 3 , and α 4 components of the dissociation energy of the H 2 isotopomer are 36118.7978 (2) 2