Accurate computations of the rovibrational spectrum of the He–HF van der Waals complex

Quantum study of the rovibrational relaxation of HF by collision with 4He on a new potential energy surface

The HF molecule is considered the main reservoir of fluorine in the interstellar medium (ISM). Also, the interactions of this molecule with the most common atoms and molecules in the ISM have attracted great interest from the astrochemical community. Collisions between HF and helium have recently caused controversy following a study using a two-dimensional SAPT potential energy surface (PES) that exhibited large discrepancies with previous scattering calculations based on more recent ab initio potentials. To address this issue, our current work aims to develop the most precise three-dimensional PES for the HF+He system. We employ the size-consistent CCSD(T) method in conjunction with the aug-cc-pV6Z basis set. The main features of the new PES as well as the bound states of the He-HF complex are compared to the existing data. The new PES is then utilised to conduct close coupling calculations that demonstrate He-HF as a good instance of vibration-rotation near resonant energy transfer. The novel rate coefficients will be accessible via the BASECOL database, and the use of the new PES is advised when describing HF in helium droplets.

Complexes of helium with neutral molecules: Progress toward a quantitative scale of bonding character

The complexes of helium with nearly 30 neutral molecules (M) were investigated by various techniques of bonding analysis and symmetry-adapted perturbation theory (SAPT). The main investigated function was the local electron energy density H(r), analyzed, in particular, so to estimate the degree of polarization (DoP) of He in the various He(M). As we showed recently (Borocci et al., J. Comput. Chem., 2019, 40, 2318-2328), the DoP is a quantitative index that is generally informative about the role of polarization (induction plus charge transfer [CT]) and dispersion in noncovalent noble gas complexes. As further evidence in this regard, we presently ascertained quantitative correlations between the DoP(He) of the He(M) and indices based on the electron density ρ(r), including the molecular electrostatic potential at the HeM bond critical point, as well as the percentage contributions of induction and dispersion to the SAPT binding energies. Based also on the explicit evaluation of the CT, accomplished through the study of the charge-displacement function, we derived a quantitative scale that ranks the He(M) according to their dispersive, inductive, and CT bonding character. Our taken approach could be conceivably extended to other types of noncovalent complexes.

Rotational energy relaxation quantum dynamics of a diatomic molecule in a superfluid helium nanodroplet and study of the hydrogen isotopes case

The rotational energy relaxation (RER) of a molecule X2(j,mj) in a 4He superfluid nanodroplet [HeND or (4He)N; T = 0.37 K] has been investigated using a hybrid quantum dynamics approach recently proposed by us. As far as we know, this is the first theoretical study about rotational relaxation inside HeNDs, and here several (real and hypothetical) isotopes of H2 have been examined, in order to analyze the influence of the rotational constant Be of these fast rotors on the dynamics. The structure of the nanodroplet practically does not change during the RER process, which approximately takes place according to a cascade mechanism j → j - 2; j - 2 → j - 4; …; 2 → 0, and mj is conserved. The results are consistent with the very scarce estimated experimental data available. The lifetime of an excited rotational state (≈1.0-7.6 ns) increases when: (a) Be increases; (b) j increases; and (c) N decreases (above N = 100 there is a small influence of N on the lifetime). This also applies to the global relaxation time and transition time. The analysis of the influence of the coupling between the j and j - 2 rotational states (due to the X2-helium interaction) and the X2 angular velocity on the lifetime and related properties has been helpful to better understand the dynamics. In contrast to the RER results, for the vibrational energy relaxation (VER) in HeNDs, when the quantum number v increases a decrease is observed in the lifetime of the excited vibrational state. This difference can be interpreted taking into account that RER and VER are associated with very different types of motion. Besides, in VER the intermediate excited states show metastability, differing from the RER case. We hope that the present study will encourage more studies to be developed on the RER dynamics in HeNDs, a basic, interesting and difficult to study physical phenomenon about which we still know very little.

Noncovalent Complexes of the Noble-Gas Atoms: Analyzing the Transition from Physical to Chemical Interactions

The bonding character of the noncovalent complexes of the noble-gas (Ng) atoms ranges from nearly purely dispersive contacts to interactions featuring appreciable contributions of induction and charge transfer. In this study, we discuss a new quantitative index that seems peculiarly informative about these diverse bonding situations. This index was termed as the degree of polarization (DoP) of Ng, as it measures, in essence, the Ng polarization promoted by the binding partner. The definition of the DoP(Ng) relies on the analysis of the local electron energy density H(r), and its physical meaning was best appreciated by studying also the charge-displacement function and the molecular electrostatic potential of the investigated benchmark species, that include nearly 60 Ngs complexes of different bonding character. The DoP(Ng) appears of general applicability, and is also positively correlated with other bonding character indices. © 2019 Wiley Periodicals, Inc.

van der Waals potential energy surfaces from the exchange-hole dipole moment dispersion model

The potential energy surfaces (PESs) of 28 simple van der Waals complexes, each consisting of a rare-gas (Rg) atom interacting with a linear molecule, are calculated using the exchange-hole dipole moment (XDM) dispersion model in conjunction with three base density functionals (HFPBE, PW86PBE, and a commensurate hybrid functional). Results are compared with literature coupled-cluster reference data. The quality of the computed PESs is assessed based on the positions of the global minima and the corresponding binding energies. Only the hybrid functional is found to provide generally reliable PESs. Dispersion-corrected HFPBE strongly underestimates the equilibrium intermolecular separations and predicts different global minima than the reference PESs for Rg–HCl, Rg–HBr, and two of the Rg–HCN complexes. Analysis of the binding-energy errors reveals that the performance of HFPBE degrades as the size of the Rg atoms increase down the group, while the performance of PW86PBE is significantly worse for strongly-polar molecules. PW86PBE, and to a lesser extent the hybrid, strongly overbind Kr–HF due to charge-transfer error. Despite this, the XDM-corrected hybrid functional displays the best overall error statistics and provides binding energies to within ca. 10 cm–1 of the coupled-cluster reference data at a greatly reduced computational cost.

Quantum solvation dynamics of HCN in a helium-4 droplet

Ultracold nanodroplets of helium-4, containing several thousands of He atoms, offer considerable promise as microscopic cryogenic chambers. Potential applications include the creation of tailor-made chemical or biomolecular complexes and studies of superfluidity in nanoscale systems. Recent experiments have succeeded in interrogating droplets of quantum solvent which consist of as few as 1-20 helium-4 atoms and which contain a single solute molecule. This allows the transition from a floppy, but essentially molecular, complex to a dissolved molecule to be followed and, surprisingly, the transition is found to occur quite rapidly, in some cases for as few as N = 7-20 solvent atoms. For example, in experiments on helium-4 droplets seeded with CO molecules [Tang and McKellar, J. Chem. Phys. 119, 754 (2003)], two series of transitions are observed which correlate with the a-type (Delta K = 0) and b-type (Delta K = +/-1) lines of the binary complex, CO-He (K is the quantum number associated with the projection of the total angular momentum onto the vector connecting the atom and the molecular center of mass). The a-type series, which evolves from the end-over-end rotational motion of the CO-He binary complex, saturates to the nanodroplet limit for as few as 10-15 helium-4 atoms, i.e., the effective moment of inertia of the molecule converges to its asymptotic (solvated) value quite rapidly. In contrast, the b-type series, which evolves from the free-molecule rotational mode, disappears altogether for N approximately 7 atoms. Similar behavior is observed in recent computational studies of HCN(4He)N droplets [Paolini et al., J. Chem. Phys. 123, 114306 (2005)]. In this article the quantum solvation of HCN in small helium-4 droplets is studied using a new fixed-node diffusion Monte Carlo (DMC) procedure. In this approach a Born-Oppenheimer-type separation of radial and angular motions is introduced as a means of computing nodal surfaces of the many-body wave functions which are required in the fixed-node DMC method. Excited rotational energies are calculated for HCN(4He)N droplets with N = 1-20: the adiabatic node approach also allows concrete physical mechanisms to be proposed for the predicted disappearance of the b-type series as well as the rapid convergence of the a-type series to the nanodroplet limit with increasing N. The behavior of the a-type series is traced directly to the mechanics of angular momentum coupling-and decoupling-between identical bosons and the molecular rotor. For very small values of N there exists significant angular momentum coupling between the molecule and the helium atoms: at N approximately 10 solvation appears to be complete as evidenced by significant decoupling of the molecule and solvent angular momenta. The vanishing of the b-type series is predicted to be a result of increasing He-He repulsion as the number of solvent atoms increases.