The pair potential for helium has been 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 were obtained from large full configuration interaction calculations. The nonadiabatic perturbation theory was used to predict the properties of the halo state of helium dimer. Its binding energy and the average value of 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., PNAS 113, 14651 (2016)].Helium is expected to become an important medium in determining thermodynamic metrology standards and the future system of SI units [1,2]. Several elements of such standards will be established by ab initio quantum mechanical calculations [3][4][5][6][7]. An important theory input is the helium pair potential. Its knowledge is required to account for the imperfection of helium gas and the necessary extrapolations to zero pressure [2]. The more accurate this potential is, the smaller will be the uncertainties of the resulting standards.There are other reasons of interest in the helium pair potential. The dimer composed of 4 He atoms, 4 He 2 , has a single very weakly bound vibrational state-an example of a quantum halo state-where atoms move mainly in the classically forbidden tunneling region of the configuration space [8]. This state was the subject of several experimental investigations [9][10][11][12][13][14]. We present here the development of a new potential with uncertainties reduced by an order of magnitude compared to the previous most accurate determination [15]. This potential and the nonadiabatic perturbation theory [16], accounting for the coupling of the electronic and nuclear motion, are used to obtain an accurate theoretical prediction of the properties of the halo state.The potential of Ref.[15] contained the BornOppenheimer (BO) component from Ref. [17]. Its uncertainty, amounting to several millikelvin (mK) in the well region, was due to the slow convergence of a part of the wave function expanded in terms of orbital products. Since it is impossible to converge the orbital expansion sufficiently well [18], we now follow Refs. [19,20] and expand the BO wave function using the four-electron explicitly correlated Gaussian (ECG) basis. Several improvements to the approach of Refs. [19,20] In Ref.[15], the BO potential of Ref.[17] was combined with the adiabatic (diagonal BO), relativistic, and quantum electrodynamics (QED) contributions, as well as with an appropriate retardation correction [24]. Its uncertainties were almost entirely determined by the uncertainties of the BO component. With the much improved BO potential computed in the present work, the accuracy of the adiabatic and relativistic co...