Ground state potential energy curves for He–Kr, Ne–Kr, Ar–Kr, and Kr2: Coupled-cluster calculations and comparison with experiment

A Simple Model for the Pauli Repulsion with Possible Utility in QM, MM and Chemical Education

The Pauli repulsion is the intermolecular force responsible for the volume and low compressibility of condensed-phase matter at normal conditions. A simple model for this force is presented, wherein per-atom electron densities are represented as spherical charge distributions that are prevented from significantly overlapping. In the example of two noble gas atoms approaching one another beyond their van der Waals radii, the distance between the centers of the electronic charge distributions becomes larger than the distance between the nuclei, giving rise to an unfavorable electrostatic interaction. For the purpose of calculating this interaction, the model is further simplified by representing the per-atom electron density as a negative point charge, loosely inspired by the classical Drude oscillator. The dispersion interaction is simplified to an R-6 term, centered on the aforementioned point charges. Despite the gross simplicity of the resulting formalism, near-quantitative agreement with high-level QM interaction energies of noble gas dimers is achieved. Accordingly, the present model is thought to have utility in force fields, in post-HF and post-DFT dispersion corrections, and in chemical education.

Tunnelling of electrons via the neighboring atom

As compared to the intuitive process that the electron emits straight to the continuum from its parent ion, there is an alternative route that the electron may transfer to and be trapped by a neighboring ionic core before the eventual release. Here, we demonstrate that electron tunnelling via the neighboring atomic core is a pronounced process in light-induced tunnelling ionization of molecules by absorbing multiple near-infrared photons. We devised a site-resolved tunnelling experiment using an Ar-Kr+ ion as a prototype system to track the electron tunnelling dynamics from the Ar atom towards the neighboring Kr+ by monitoring its transverse momentum distribution, which is temporally captured into the resonant excited states of the Ar-Kr+ before its eventual releasing. The influence of the Coulomb potential of neighboring ionic cores promises new insights into the understanding and controlling of tunnelling dynamics in complex molecules or environment.

Trapping and thermal migration of the first- and second-row atoms in Ar, Kr and Xe crystals

Trapping and temperature-induced migration (TIM) of the first- and second-row atoms A from H to Ne in the face-centered cubic rare gas RG = Ar, Kr and Xe crystals are investigated within the classical crystal model parameterized by the empirically modified pairwise potentials. New ab initio coupled cluster A-RG potentials computed in a uniform way for all the atoms A are used to represent the atom-crystal interactions. Absolute and relative stabilities of the substitutional and interstitial trapping sites, their structures, interstitial migration pathways, related activation energies and rough estimates of the TIM rates are obtained. The isotropic model, which neglects non-zero atomic electronic orbital momentum, reveals that migration of interstitial atoms along the network of conjugated fcc octahedral voids is the generic case for atomic mobility. Anisotropic interactions with a crystal inherent to P-state atoms B, C, O and F are accounted for using the non-relativistic diatomics-in-molecule method. Depending on its sign, interaction anisotropy can alter the structures of interstitial trapping sites and transition states remarkably. This, in turn, can dramatically affect the TIM rates. Comparison with reliable experimental data available for oxygen and hydrogen indicates a systematic overestimation of the measured activation energies, by 30% at worst. A comprehensive literature review accomplished for other atoms reveals a lack of information on the TIM processes and rates, though makes it possible to verify a part of the present results on the trapping site energies and structures.

Universal Pairwise Interatomic van der Waals Potentials Based on Quantum Drude Oscillators

Repulsive short-range and attractive long-range van der Waals (vdW) forces play an appreciable role in the behavior of extended molecular systems. When using empirical force fields, the most popular computational methods applied to such systems, vdW forces are typically described by Lennard-Jones-like potentials, which unfortunately have a limited predictive power. Here, we present a universal parameterization of a quantum-mechanical vdW potential, which requires only two free-atom properties─the static dipole polarizability α1 and the dipole-dipole C6 dispersion coefficient. This is achieved by deriving the functional form of the potential from the quantum Drude oscillator (QDO) model, employing scaling laws for the equilibrium distance and the binding energy, and applying the microscopic law of corresponding states. The vdW-QDO potential is shown to be accurate for vdW binding energy curves, as demonstrated by comparing to the ab initio binding curves of 21 noble-gas dimers. The functional form of the vdW-QDO potential has the correct asymptotic behavior at both zero and infinite distances. In addition, it is shown that the damped vdW-QDO potential can accurately describe vdW interactions in dimers consisting of group II elements. Finally, we demonstrate the applicability of the atom-in-molecule vdW-QDO model for predicting accurate dispersion energies for molecular systems. The present work makes an important step toward constructing universal vdW potentials, which could benefit (bio)molecular computational studies.

Probing Resonant Photoionization Time Delay by Self-Referenced Molecular Attoclock

Attosecond time-resolved electron tunneling dynamics have been investigated by using attosecond angular streaking spectroscopy, where a clock reference to the laser field vector is required in atomic strong-field ionization and the situation becomes complicated in molecules. Here we reveal a resonant ionization process via a transient state by developing an electron-tunneling-site-resolved molecular attoclock in Ar-Kr^{+}. Two distinct deflection angles are observed in the photoelectron angular distribution in the molecular frame, corresponding to the direct and resonant ionization pathways. We find the electron is temporally trapped in the Coulomb potential wells of the Ar-Kr^{+} before finally releasing into the continuum when the electron tunnels through the internal barrier. By utilizing the direct tunneling ionization as a self-referenced arm of the attoclock, the time delay of the electron trapped in the resonant state is revealed to be 3.50±0.04  fs. Our results give an impetus to exploring the ultrafast electron dynamics in complex systems and also endow a semiclassical presentation of the electron trapping dynamics in a quantum resonant state.

Ultrafast temporal evolution of interatomic Coulombic decay in NeKr dimers

We investigate interatomic Coulombic decay in NeKr dimers after neon inner-valence photoionization [Ne⁺(2s^-1)] using a synchrotron light source. We measure with high energy resolution the two singly charged ions of the Coulomb-exploding dimer dication and the photoelectron in coincidence. By carefully tracing the post-collision interaction between the photoelectron and the emitted ICD electron we are able to probe the temporal evolution of the state as it decays. Although the ionizing light pulses are 80 picoseconds long, we determine the lifetime of the intermediate dimer cation state and visualize the contraction of the nuclear structure on the femtosecond time scale.

Conformal Analytical Potential for All the Rare Gas Dimers over the Full Range of Internuclear Distances

An analytical model for the potential between two rare gas atoms at distances between R=0 to R→∞ is assumed to be conformal with the previously published potential for He_{2} [J. Chem. Phys. 142, 131102 (2015)JCPSA60021-960610.1063/1.4916740]. The potential curves of the rare gas dimers all have the same shape and only depend on the well parameters D_{e} and R_{e}. The potentials and the vibrational levels for the 11 homonuclear and heteronuclear dimers for which recent ab initio calculations are available agree, within several percent, with the ab initio results. For the other rare gas dimers, the new potential provides the first realistic estimates for the potentials.

Molecular dynamics calculations of collision-induced absorption in a gas mixture of neon and krypton

We continue the development of the in-house molecular dynamics software package SpaCIAL and test it for the computation of the collision-induced absorption coefficients for a neon (Ne) and krypton (Kr) gas mixture. An apodization procedure for the dipole autocorrelation function is implemented and tested. We also carry out a statistical study of the convergence rate with respect to ensemble size. The resulting absorption coefficients show a good accordance with quantum mechanical results. Comparison with laboratory measurements shows agreement within 10%-20% at T = 295 K. At T = 480 K, a larger difference of 40%-80% is observed, which can presumably be explained by experimental uncertainties. For the study, an empirical (Barker, Fisher, and Watts) interaction-potential [Mol. Phys. 21, 657 (1971)] for Ne-Kr has been developed. Ab initio {coupled cluster with singles and doubles (triples) [CCSD(T)]} potentials for Ne-Ne, Kr-Kr, and Ne-Kr have been computed, as well as the CCSD(T) interaction-induced Ne-Kr dipole moment curve.

Determination of the Residual Entropy of Simple Mixtures by Monte Carlo Simulation

Extensive Monte Carlo simulations were performed for (neon + krypton) mixtures for temperatures between 200 K and 600 K and pressures up to 1 GPa, using Lennard-Jones potentials to describe the intermolecular interactions. The residual entropies were obtained via Widom's insertion method, as well as via an integration technique. At high pressures, the residual entropy is, to a very good approximation, a linear function of λ_a^-1, which is the reciprocal value of the average Monte Carlo displacement parameter that gives the acceptance ratio a for translational moves. The slope of this linear function varies linearly with the mole fraction and is related to the effective collision diameters of the molecules. As the displacement parameter is available during Monte Carlo simulations of fluids, its linearity with the residual entropy can be used to compute this properties with negligible computational effort at high densities, when particle insertion methods become unreliable.

Thermophysical properties of krypton-helium gas mixtures from ab initio pair potentials

A new potential energy curve for the krypton-helium atom pair was developed using supermolecular ab initio computations for 34 interatomic distances. Values for the interaction energies at the complete basis set limit were obtained from calculations with the coupled-cluster method with single, double, and perturbative triple excitations and correlation consistent basis sets up to sextuple-zeta quality augmented with mid-bond functions. Higher-order coupled-cluster excitations up to the full quadruple level were accounted for in a scheme of successive correction terms. Core-core and core-valence correlation effects were included. Relativistic corrections were considered not only at the scalar relativistic level but also using full four-component Dirac-Coulomb and Dirac-Coulomb-Gaunt calculations. The fitted analytical pair potential function is characterized by a well depth of 31.42 K with an estimated standard uncertainty of 0.08 K. Statistical thermodynamics was applied to compute the krypton-helium cross second virial coefficients. The results show a very good agreement with the best experimental data. Kinetic theory calculations based on classical and quantum-mechanical approaches for the underlying collision dynamics were utilized to compute the transport properties of krypton-helium mixtures in the dilute-gas limit for a large temperature range. The results were analyzed with respect to the orders of approximation of kinetic theory and compared with experimental data. Especially the data for the binary diffusion coefficient confirm the predictive quality of the new potential. Furthermore, inconsistencies between two empirical pair potential functions for the krypton-helium system from the literature could be resolved.

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.

Ab initio simulation of transport phenomena in rarefied gases

Ab initio potentials are implemented into the direct simulation Monte Carlo (DSMC) method. Such an implementation allows us to model transport phenomena in rarefied gases without any fitting parameter of intermolecular collisions usually extracted from experimental data. Applying the method proposed by Sharipov and Strapasson [Phys. Fluids 24, 011703 (2012)], the use of ab initio potentials in the DSMC requires the same computational efforts as the widely used potentials such as hard spheres, variable hard sphere, variable soft spheres, etc. At the same time, the ab initio potentials provide more reliable results than any other one. As an example, the transport coefficients of a binary mixture He-Ar, viz., viscosity, thermal conductivity, and thermal diffusion factor, have been calculated for several values of the mole fraction.

Consistent van der Waals radii for the whole main group

Atomic radii are not precisely defined but are nevertheless widely used parameters in modeling and understanding molecular structure and interactions. The van der Waals radii determined by Bondi from molecular crystals and data for gases are the most widely used values, but Bondi recommended radius values for only 28 of the 44 main-group elements in the periodic table. In the present Article, we present atomic radii for the other 16; these new radii were determined in a way designed to be compatible with Bondi's scale. The method chosen is a set of two-parameter correlations of Bondi's radii with repulsive-wall distances calculated by relativistic coupled-cluster electronic structure calculations. The newly determined radii (in A) are Be, 1.53; B, 1.92; Al, 1.84; Ca, 2.31; Ge, 2.11; Rb, 3.03; Sr, 2.49; Sb, 2.06; Cs, 3.43; Ba, 2.68; Bi, 2.07; Po, 1.97; At, 2.02; Rn, 2.20; Fr, 3.48; and Ra, 2.83.

London dispersion forces by range-separated hybrid density functional with second order perturbational corrections: The case of rare gas complexes

A satisfactory account of the van der Waals (vdW) (London dispersion) forces is, in general not possible by the Kohn-Sham method using standard local, semilocal generalized gradient approximation (GGA), or meta-GGA density functionals. The recently proposed range-separated hybrid (RSH) approach, supplemented by second order perturbational corrections (MP2) to include long-range dynamic correlation effects, offers a physically consistent, seamless description of dispersion forces. It is based on a rigorous generalization of the Kohn-Sham method, where long-range exchange and correlation effects are treated by wave function methods, while short-range electron exchange and correlation are handled by local or semilocal functionals. The method is tested on a series of rare gas dimers in comparison with standard wave function theory and density functional theory approaches. In contrast to the most successful exchange correlation functionals, which describe at best the vdW minimum, the RSH+MP2 approach is valid also in the asymptotic region and the potential curve displays the correct 1/R(6) behavior at large internuclear separations. In contrast to usual MP2 calculations, the basis set superposition error is considerably reduced, making RSH+MP2 an ideal tool for exploring the potential energy surface of weakly bound molecular complexes.

Rydberg states of the rare gas dimers

This work is dedicated to Gerhard Herzberg and his prodigious contributions to molecular spectroscopy. Of particular relevance here is Herzberg's seminal 1987 article (Annu. Rev. Phys. Chem. 38, 27 (1987)) in which he discussed the electronic structures of several groups of molecules he termed "Rydberg molecules". Among these are the rare gas dimers (Rg2), a group whose study has benefited significantly from recent advances in laser excitation and synchrotron-based spectroscopies, as well as in theory. Following the spirit of Herzberg's 1987 article, this paper reviews some of the more prominent features of Rydberg states in the Rg2 family as viewed from the current perspective.Key words: rare gas dimers, rare gas dimer excited states, Rg2 electronic structures, Rydberg states of He2, Ne2, Ar2, Kr2, and Xe2, electronic states of He2, Ne2, Ar2, Kr2, and Xe2.