An accurate potential energy curve for helium based on ab initio calculations

Reduced Radial Curves of Diatomic Molecules

The prospect of using the concept of a universal reduced potential energy curve (RPC) for a broader class of radial molecular functions is explored by performing appropriate model calculations for the electric dipole moment functions of the hydrogen halides HF, HCl, and HBr. The reduced radial functions of the model systems, constructed from their best available theoretical approximants, coincide so closely that they can be used as few-parameter universal representations of functions available in the literature. Given the mathematical nature of the problem addressed here, the results are not limited to the functions studied but can be applied equally well to all radial molecular functions that have similar shapes, such as electric quadrupole moment and dipole polarizability functions.

The Solid Phase of 4He: A Monte Carlo Simulation Study

The thermodynamics of solid (hcp) 4He is studied theoretically by means of unbiased Monte Carlo simulations at finite temperature, in a wide range of density. This study complements and extends previous theoretical work, mainly by obtaining results at significantly lower temperatures (down to 60 mK) and for systems of greater size, by including in full the effect of quantum statistics, and by comparing estimates yielded by different pair potentials. All the main thermodynamic properties of the crystal, e.g., the kinetic energy per atom, are predicted to be essentially independent of temperature below ∼ 1 K. Quantum-mechanical exchanges are virtually non-existent in this system, even at the lowest temperature considered. However, effects of quantum statistics are detectable in the momentum distribution. Comparison with available measurements shows general agreement within the experimental uncertainties.

A Glimpse into Quantum Triplet Structures in Supercritical 3He

A methodological study of triplet structures in quantum matter is presented. The focus is on helium-3 under supercritical conditions (4 < T/K < 9; 0.022 < ρN/Å-3 < 0.028), for which strong quantum diffraction effects dominate the behavior. Computational results for the triplet instantaneous structures are reported. Path integral Monte Carlo (PIMC) and several closures are utilized to obtain structure information in the real and the Fourier spaces. PIMC involves the fourth-order propagator and the SAPT2 pair interaction potential. The main triplet closures are: AV3, built as the average of the Kirkwood superposition and the Jackson-Feenberg convolution, and the Barrat-Hansen-Pastore variational approach. The results illustrate the main characteristics of the procedures employed by concentrating on the salient equilateral and isosceles features of the computed structures. Finally, the valuable interpretive role of closures in the triplet context is highlighted.

Dimerization dynamics of carboxylic acids in helium nanodroplets

The dimerization of molecules in helium nanodroplets is known to preferentially yield structures of higher energy than the global energy minimum structure for a number of quite different monomers. Here, we explore dimerization in this environment using an atomistic model within statistically converged molecular dynamics (MD) trajectories, treating the solvent implicitly through the use of a thermostat, or more explicitly by embedding one monomer in a He₁₀₀ cluster. The focus is on the two simplest carboxylic acids, formic and acetic, both of which have been studied experimentally. While the global minimum structure, which comprises two CO⋯HO hydrogen bonds, is predicted to be the most abundant dimer in the absence of the helium solvent, this is no longer the case once helium atoms are included. The simulations confirm the importance of kinetic trapping effects and also shed light on the occurrence of specific dynamical effects, leading to the occasional formation of high-energy structures away from minima, such as saddle configurations. Theoretically predicted infrared spectra, based on the MD statistics, are in good agreement with the experimental spectra.

Probing the Helium Dimer by Relativistic Highly Charged Projectiles

We study the fragmentation of He_{2} dimers into He^{+} ions by relativistic highly charged projectiles. We demonstrate that the interaction between an ultrafast projectile with an extremely extended object-the helium dimer-possesses interesting features that are absent in collisions with "normal" molecules. We also show that such projectiles, due to their enormous interaction range, can accurately probe the ground state of the dimer and even be used for a determination of its binding energy.

Fragmentation dynamics of Ar4He1000 upon electron impact ionization: Competition between ion ejection and trapping

The fragmentation upon electron impact ionization of Ar₄He₁₀₀₀ is investigated by means of mixed quantum-classical dynamics simulations. The Ar₄ ⁺ dopant dynamics is described by a surface hopping method coupled with a diatomics-in-molecules model to properly take into account the multiple Ar₄ ⁺ electronic surfaces and possible transitions between them. Helium atoms are treated individually using zero-point averaged dynamics, a method based on the building of an effective He-He potential. Fast electronic relaxation is observed from less than 2 ps to ∼30 ps, depending on initial conditions. The main fragments observed are Ar₂ ⁺He_q and Ar₃ ⁺He_q (q ≤ 1000), with a strong contribution of the bare Ar₂ ⁺ ion, and neither Ar⁺ nor Ar⁺He_q fragments are found. The smaller fragments (q ≤ 50) are found to mostly come from ion ejection, whereas larger fragments (q > 500) originate from long-term ion trapping. Although the structure of the trapped Ar₂ ⁺ ions is the same as in the gas phase, trapped Ar₃ ⁺ and Ar₄ ⁺ are rather slightly bound Ar₂ ⁺⋯Ar and Ar₂ ⁺⋯Ar⋯Ar structures (i.e., an Ar₂ ⁺ core with one or two argon atoms roaming within the droplet). These loose structures can undergo geminate recombination and release Ar₃ ⁺He_q or Ar₄ ⁺He_q (q ≤ 50) in the gas phase and/or induce strong helium droplet evaporation. Finally, the translational energy of the fragment center of mass was found to be suitable to provide a clear signature of the broad variety of processes at play in our simulations.

Inelastic scattering of photoelectrons from He nanodroplets

We present a detailed study of inelastic energy-loss collisions of photoelectrons emitted from He nanodroplets by tunable extreme ultraviolet (XUV) radiation. Using coincidence imaging detection of electrons and ions, we probe the lowest He droplet excited states up to the electron impact ionization threshold. We find significant signal contributions from photoelectrons emitted from free He atoms accompanying the He nanodroplet beam. Furthermore, signal contributions from photoionization and electron impact excitation/ionization occurring in pairs of nearest-neighbor atoms in the He droplets are detected. This work highlights the importance of inelastic electron scattering in the interaction of nanoparticles with XUV radiation.

Isomeric Broadening of C60+ Electronic Excitation in Helium Droplets: Experiments Meet Theory

Helium is considered an almost ideal tagging atom for cold messenger spectroscopy experiments. Although helium is bound very weakly to the ionic molecule of interest, helium tags can lead to shifts and broadenings that we recorded near 963.5 nm in the electronic excitation spectrum of C60+ solvated with up to 100 helium atoms. Dedicated quantum calculations indicate that the inhomogeneous broadening is due to different binding energies of helium to the pentagonal and hexagonal faces of C60+, their dependence on the electronic state, and the numerous isomeric structures that become available for intermediate coverage. Similar isomeric effects can be expected for optical spectra of most larger molecules surrounded by nonabsorbing weakly bound solvent molecules, a situation encountered in many messenger-tagging spectroscopy experiments.

Pair Potential with Submillikelvin Uncertainties and Nonadiabatic Treatment of the Halo State of the Helium Dimer

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].

Dirac bubble potential for He-He and inadequacies in the continuum: Comparing an analytic model with elastic collision experiments

We focus on the long-pending issue of the inadequacy of the Dirac bubble potential model in the description of He-He interactions in the continuum [L. L. Lohr and S. M. Blinder, Int. J. Quantum Chem. 53, 413 (1995)]. We attribute this failure to the lack of a potential wall to mimic the onset of the repulsive interaction at close range separations. This observation offers the explanation to why this excessively simple model proves incapable of quantitatively reproducing previous experimental findings of glory scattering in He-He, although being notorious for its capability of reproducing several distinctive features of the atomic and isotopic helium dimers and trimers [L. L. Lohr and S. M. Blinder, Int. J. Quantum Chem. 90, 419 (2002)]. Here, we show that an infinitely high, energy-dependent potential wall of properly calculated thickness r_c(E) taken as a supplement to the Dirac bubble potential suffices for agreement with variable-energy elastic collision cross section experiments for ⁴He-⁴He, ³He-⁴He, and ³He-³He [R. Feltgen et al., J. Chem. Phys. 76, 2360 (1982)]. In the very low energy regime, consistency is found between the Dirac bubble potential (to which our extended model is shown to reduce) and cold collision experiments [J. C. Mester et al., Phys. Rev. Lett. 71, 1343 (1993)]; this consistency, which in this regime lends credence to the Dirac bubble potential, was never noticed by its authors. The revised model being still analytic is of high didactical value while expected to increase in predictive power relative to other appraisals.

Imaging the He2 quantum halo state using a free electron laser

Quantum tunneling is a ubiquitous phenomenon in nature and crucial for many technological applications. It allows quantum particles to reach regions in space which are energetically not accessible according to classical mechanics. In this "tunneling region," the particle density is known to decay exponentially. This behavior is universal across all energy scales from nuclear physics to chemistry and solid state systems. Although typically only a small fraction of a particle wavefunction extends into the tunneling region, we present here an extreme quantum system: a gigantic molecule consisting of two helium atoms, with an 80% probability that its two nuclei will be found in this classical forbidden region. This circumstance allows us to directly image the exponentially decaying density of a tunneling particle, which we achieved for over two orders of magnitude. Imaging a tunneling particle shows one of the few features of our world that is truly universal: the probability to find one of the constituents of bound matter far away is never zero but decreases exponentially. The results were obtained by Coulomb explosion imaging using a free electron laser and furthermore yielded He2's binding energy of [Formula: see text] neV, which is in agreement with most recent calculations.

Transport coefficients of helium-argon mixture based on ab initio potential

The viscosity, thermal conductivity, diffusion coefficient, and thermal diffusion factor of helium-argon mixtures are calculated for a wide range of temperature and for various mole fractions up to the 12th order of the Sonine polynomial expansion with an ab initio intermolecular potential. The calculated values for these transport coefficients are compared with other data available in the open literature. The comparison shows that the obtained transport coefficients of helium-argon mixture have the best accuracy for the moment.

Onset of Casimir-Polder retardation in a long-range molecular quantum state

Two 4He atoms form a diatomic molecule with a significant vibrational wave function amplitude at interatomic separations R>100  Å, where the retardation switches the London R(-6) decay of the potential to the Casimir-Polder R(-7) form. It has been assumed that this effect of retardation on the long-range part of the potential is responsible for the 2 Å (4%) increase of the bond length <R> of 4He2. We show that <R> is, unexpectedly, insensitive to the potential at R>20  Å and its increase is due to quantum electrodynamics effects computed by us from expressions valid at short R--beyond the validity range of Casimir-Polder theory--that seamlessly extend this theory to distances relevant for properties of long molecules.

Approaching universality in weakly bound three-body systems

Atom-dimer scattering below the three-body breakup threshold is studied for a system of three identical bosons. The atom-dimer scattering length and the energy of the most weakly bound three-body state are shown to be strongly correlated. An appropriate rescaling of the observables reveals the subtlety of the correlation and serves to identify universal trends in the unitary limit of divergent two-body scattering length. The correlation provides a new quantitative measure of the degree of universality in three-body systems with short-ranged interactions, as well as a consistency check of effective field theories and other theoretical models.

Counterpoise correction is not useful for short and Van der Waals distances but may be useful at long range

This article investigates the errors in supermolecule calculations for the helium dimer. In a full CI calculation, there are two errors. One is the basis set superposition error (BSSE), the other is the basis set convergence error (BSCE). Both of the errors arise from the incompleteness of the basis set. These two errors make opposite contributions to the interaction energies. The BSCE is by far the largest error in the short range and larger than (but much closer to) BSSE around the Van der Waals minimum. Only at the long range, the BSSE becomes the larger error. The BSCE and BSSE largely cancel each other over the Van der Waals well. Accordingly, it may be recommended to not include the BSSE for the calculation of the potential energy curve from short distance till well beyond the Van der Waals minimum, but it may be recommended to include the BSSE correction if an accurate tail behavior is required. Only if the calculation has used a very large basis set, one can refrain from including the counterpoise correction in the full potential range. These results are based on full CI calculations with the aug-cc-pVXZ (X = D, T, Q, 5) basis sets.

Injection of atoms and molecules in a superfluid helium fountain: Cu and Cu2He(n) (n = 1, ..., ∞)

We introduce an experimental platform designed around a thermomechanical helium fountain, which is aimed at investigating spectroscopy and dynamics of atoms and molecules in the superfluid and at its vapor interface. Laser ablation of copper, efficient cooling and transport of Cu and Cu(2) through helium vapor (1.5 K < T < 20 K), formation of linear and T-shaped Cu(2)-He complexes, and their continuous evolution into large Cu(2)-He(n) clusters and droplets are among the processes that are illustrated. Reflection is the dominant quantum scattering channel of translationally cold copper atoms (T = 1.7 K) at the fountain interface. Cu(2) dimers mainly travel through the fountain unimpeded. However, the conditions of fountain flow and transport of molecules can be controlled to demonstrate injection and, in particular, injection into a nondivergent columnar fountain with a plug velocity of about 1 m/s. The experimental observables are interpreted with the aid of bosonic density functional theory calculations and ab initio interaction potentials.

Structures and energetics of helium cluster cations: equilibrium geometries revisited through the genetic algorithm approach

Equilibrium geometries and dissociation energies of He(N)(+) clusters have been calculated for N=3-35 using an extended genetic algorithm approach and a semiempirical model of intracluster interactions [P. J. Knowles, J. N. Murrell, and E. J. Hodge, Mol. Phys. 85, 243 (1995)]. A general aufbau principle is formulated for both ionic cores and neutral solvation shells, and the results are thoroughly compared with other theoretical data available for helium cluster cations in literature.

Computational study of the structures of gaseous helium-3 at low temperature

A computational study of gaseous helium-3 at T=5.23 K, for number densities rho N<0.0021 A(-3), analyzing the different pair and triplet structures in both the r and the k spaces, is presented. Structures in r space (i.e., instantaneous, total continuous linear response, and centroids) are determined via path-integral Monte Carlo simulations in the canonical ensemble by utilizing the Aziz-Slaman and the SAPT2 interatomic potentials. Additional results obtained with the application of two closures for triplets in r space, the Kirkwood superposition approximation and the Jackson-Feenberg convolution, are also reported. Besides, an analysis of the nonsuitability of quantum hard spheres for describing this system is included. The pair structures in k space are fixed via Ornstein-Zernike schemes appropriate for dealing with quantum diffraction effects in fluids. The effect on the quality of the computed isothermal compressibilities brought about by increasing the sample size in the simulations and by the subsequent application of a grand ensemble correction to the asymptotic behavior of the canonical pair radial correlation functions is also investigated. Furthermore, it is demonstrated analytically that the methods of classical statistical mechanics for dealing with the higher-order direct correlation functions remain fully valid for studying the higher-order correlations of path-integral centroids. By taking advantage of this result, the triplet structure factors for the centroid (also for the instantaneous) correlations are computed by following a number of distinct closures and methods that involve triplet direct correlation functions. The latter computations are intended to explore an alternative scheme to the very expensive fixing of triplet structure factors through direct path integral simulations, an alternative which is expected to yield the main features of these triplet quantities for this gas. Comparison with experiment is made wherever possible, and the results presented allow one to explain the substantial structural features existing in gaseous helium-3.

Fragmentation of ionized doped helium nanodroplets: theoretical evidence for a dopant ejection mechanism

We report a theoretical study of the effect induced by a helium nanodroplet environment on the fragmentation dynamics of a dopant. The dopant is an ionized neon cluster Ne(n) (+) (n=4-6) surrounded by a helium nanodroplet composed of 100 atoms. A newly designed mixed quantum/classical approach is used to take into account both the large helium cluster zero-point energy due to the light mass of the helium atoms and all the nonadiabatic couplings between the Ne(n) (+) potential-energy surfaces. The results reveal that the intermediate ionic dopant can be ejected from the droplet, possibly with some helium atoms still attached, thereby reducing the cooling power of the droplet. Energy relaxation by helium atom evaporation and dissociation, the other mechanism which has been used in most interpretations of doped helium cluster dynamics, also exhibits new features. The kinetic energy distribution of the neutral monomer fragments can be fitted to the sum of two Boltzmann distributions, one with a low kinetic energy and the other with a higher kinetic energy. This indicates that cooling by helium atom evaporation is more efficient than was believed so far, as suggested by recent experiments. The results also reveal the predominance of Ne(2) (+) and He(q)Ne(2) (+) fragments and the absence of bare Ne(+) fragments, in agreement with available experimental data (obtained for larger helium nanodroplets). Moreover, the abundance in fragments with a trimeric neon core is found to increase with the increase in dopant size. Most of the fragmentation is achieved within 10 ps and the only subsequent dynamical process is the relaxation of hot intermediate He(q)Ne(2) (+) species to Ne(2) (+) by helium atom evaporation. The dependence of the ionic fragment distribution on the parent ion electronic state reached by ionization is also investigated. It reveals that He(q)Ne(+) fragments are produced only from the highest electronic state, whereas He(q)Ne(2) (+) fragments originate from all the electronic states. Surprisingly, the highest electronic states also lead to fragments that still contain the original ionic dopant species. A mechanism is conjectured to explain this fragmentation inhibition.

Pair potential for helium from symmetry-adapted perturbation theory calculations and from supermolecular data

Symmetry-adapted perturbation theory (SAPT) was applied to the helium dimer for interatomic separations R from 3 to 12 bohrs. The first-order interaction energy and the bulk of the second-order contribution were obtained using Gaussian geminal basis sets and are converged to about 0.1 mK near the minimum and for larger R. The remaining second-order contributions available in the SAPT suite of codes were computed using very large orbital basis sets, up to septuple-zeta quality, augmented by diffuse and midbond functions. The accuracy reached at this level was better than 1 mK in the same region. All the remaining components of the interaction energy were computed using the full configuration interaction method in bases up to sextuple-zeta quality. The latter components, although contributing only 1% near the minimum, have the largest uncertainty of about 10 mK in this region. The total interaction energy at R=5.6 bohrs is -11.000+/-0.011 K. For R< or =6.5 bohrs, the supermolecular (SM) interaction energies computed by us recently turned out to be slightly more accurate. Therefore, we have combined the SM results for R< or =6.5 bohrs with the SAPT results from 7.0 to 12 bohrs to fit analytic functions for the potential and for its error bars. The potential fit uses the best available van der Waals constants C(6) through C(16), including C(11), C(13), and C(15), and is believed to be the best current representation of the Born-Oppenheimer (BO) potential for helium. Using these fits, we found that the BO potential for the helium dimer exhibits the well depth D(e)=11.006+/-0.004 K, the equilibrium distance R(e)=5.608+/-0.012 bohrs, and supports one bound state for (4)He(2) with the dissociation energy D(0)=1.73+/-0.04 mK, and the average interatomic separation R=45.6+/-0.5 A.

Electronic Spectroscopy of Nonalternant Hydrocarbons Inside Helium Nanodroplets

We have recorded the electronic spectra of three polycyclic aromatic hydrocarbons (acenaphtylene, fluoranthene, and benzo(k)fluoranthene) containing a five-member ring and their van der Waals complexes with argon and oxygen with a molecular beam superfluid helium nanodroplet spectrometer. Although the molecules, which differ by addition of one or two fused benzene rings to acenaphtylene, have the same point group symmetry, the spectral lineshapes show distinct differences in the number of zero phonon lines and shapes of the phonon wings. Whereas the smallest molecule (acenaphtylene) has the most complicated line shape, the largest molecule (benzo(k)fluoranthene) shows different lineshapes for different vibronic transitions. The van der Waals complexes of fluoranthene exhibit more peaks than the theoretically allowed number of isomeric complexes with argon/oxygen. The current models of molecular solvation in liquid helium do not adequately explain these discrepancies.

Electronic Spectroscopy of Biphenylene Inside Helium Nanodroplets

We have recorded the S1 <-- S0 electronic spectra of Biphenylene and its Ar and O2 van der Waals complexes inside helium nanodroplets using beam depletion detection. In general, the spectrum is similar to the previously reported high-resolution REMPI spectrum. The zero phonon lines, however, are split similar to the previously reported tetracene case. The calculated potential energy surface predicts that helium atoms can simultaneously occupy all equivalent global minima positions. Therefore, it appears that the splitting cannot be explained either by different isomers or by tunneling. Furthermore, surprisingly the splitting is retained for the Ar van der Waals complexes (and possibly for the O2 complex as well). This case suggests that the current models of the origin of zero phonon line splitting and the helium solvation are incomplete.

Photoionization Dynamics in Pure Helium Droplets

The photoionization and photoelectron spectroscopy of pure He droplets were investigated at photon energies between 24.6 eV (the ionization energy of He) and 28.0 eV. Time-of-flight mass spectra and photoelectron images were obtained at a series of molecular beam source temperatures and pressures to assess the effect of droplet size on the photoionization dynamics. At source temperatures below 16 K, where there is significant production of clusters with more than 10(4) atoms, the photoelectron images are dominated by fast electrons produced via direct ionization, with a small contribution from very slow electrons with kinetic energies below 1 meV arising from an indirect mechanism. The fast photoelectrons from the droplets have as much as 0.5 eV more kinetic energy than those from atomic He at the same photon energy. This result is interpreted and simulated within the context of a "dimer model", in which one assumes vertical ionization from two nearest-neighbor He atoms to the attractive region of the He2+ potential energy curve. Possible mechanisms for the slow electrons, which were also seen at energies below IE(He), are discussed, including vibrational autoionizaton of Rydberg states comprising an electron weakly bound to the surface of a large HeN+ core.

Accurate Pair Interaction Energies for Helium from Supermolecular Gaussian Geminal Calculations

Nonrelativistic clamped-nuclei pair interaction energy for ground-state helium atoms has been computed for 12 interatomic separations ranging from 3.0 to 9.0 bohr. The calculations applied the supermolecular approach. The major part of the interaction energy was obtained using the Gaussian geminal implementation of the coupled-cluster theory with double excitations (CCD). Relatively small contributions from single, triple, and quadruple excitations were subsequently included employing the conventional orbital coupled-cluster method with single, double, and noniterative triple excitations [CCSD(T)] and the full configuration interaction (FCI) method. For three distances, the single-excitation contribution was taken from literature Gaussian-geminal calculations at the CCSD level. The orbital CCSD(T) and FCI calculations used very large basis sets, up to doubly augmented septuple- and sextuple-zeta size, respectively, and were followed by extrapolations to the complete basis set limits. The accuracy of the total interaction energies has been estimated to be about 3 mK or 0.03% at the minimum of the potential well. For the attractive part of the well, the relative errors remain consistently smaller than 0.03%. In the repulsive part, the accuracy is even better, except, of course, for the region where the potential goes through zero. For interatomic separations smaller than 4.0 bohr, the relative errors do not exceed 0.01%. Such uncertainties are significantly smaller than the expected values of the relativistic and diagonal Born-Oppenheimer contributions to the potential.

Radiative correction to the helium dimer interaction energy

The leading-order radiative correction to the helium-helium interaction energy at the equilibrium internuclear distance has been calculated for the first time. The result is -1.27(2) mK. The calculations were performed using a new technique of evaluating expectation values of singular operators in connection with the most accurate wave functions of He(2) available today-the exponentially correlated Gaussian functions.