Spin-polarized Rb2 interacting with bosonic He atoms: potential energy surface and quantum structures of small clusters

Computational Modeling: Up-to-Date Approaches and Cutting-Edge Applications from Clusters, Nanostructures to Bulk Systems

The contributions in this special theme collection, in honor to Prof. P. Villarreal, cover a broad variety of computational methodologies and experimental techniques, containing studies on gas phase, clusters and condensed phase systems.

Structural, Spectroscopic, and Dynamic Properties of Li2+(X2∑g+) in Interaction with Krypton Atom

We report a computational study of the potential energy surface (PES) and vibrational bound states for the ground electronic state of Li2+Kr. The PES was calculated in Jacobi coordinates at the Restricted Coupled Cluster method RCCSD(T) level of calculation and using aug-cc-pVnZ (n = 4 and 5) basis sets. Afterward, this PES is extrapolated to the complete basis set (CBS) limit for correction. The obtained interaction energies were, then, interpolated numerically using the reproducing kernel Hilbert space polynomial (RKHS) approach to produce analytic expressions for the 2D-PES. The analytical PES is used to solve the nuclear Schrodinger equation to determine the bound states' eigenvalues of Li2+Kr for a J = 0 total angular momentum configuration and to understand the effects of orientational anisotropy of the forces and the interplay between the repulsive and attractive interaction within the potential surface. In addition, the radial and angular distributions of some selected bound state levels, which lie below, around, and above the T-shaped 90° barrier well, are calculated and discussed. We note that the radial distributions clearly acquire a more complicated nodal structure and correspond to bending and stretching vibrational motions "mode" of the Kr atom along the radial coordinate, and the situation becomes very different at the highest bound states levels with energies higher than the T-shaped 90° barrier well. The shape of the distributions becomes even more complicated, with extended angular distributions and prominent differences between even and odd states.

A Path Integral Molecular Dynamics Simulation of a Harpoon-Type Redox Reaction in a Helium Nanodroplet

We present path integral molecular dynamics (PIMD) calculations of an electron transfer from a heliophobic Cs2 dimer in its (3Σu) state, located on the surface of a He droplet, to a heliophilic, fully immersed C60 molecule. Supported by electron ionization mass spectroscopy measurements (Renzler et al., J. Chem. Phys.2016, 145, 181101), this spatially quenched reaction was characterized as a harpoon-type or long-range electron transfer in a previous high-level ab initio study (de Lara-Castells et al., J. Phys. Chem. Lett.2017, 8, 4284). To go beyond the static approach, classical and quantum PIMD simulations are performed at 2 K, slightly below the critical temperature for helium superfluidity (2.172 K). Calculations are executed in the NVT ensemble as well as the NVE ensemble to provide insights into real-time dynamics. A droplet size of 2090 atoms is assumed to study the impact of spatial hindrance on reactivity. By changing the number of beads in the PIMD simulations, the impact of quantization can be studied in greater detail and without an implicit assumption of superfluidity. We find that the reaction probability increases with higher levels of quantization. Our findings confirm earlier, static predictions of a rotational motion of the Cs2 dimer upon reacting with the fullerene, involving a substantial displacement of helium. However, it also raises the new question of whether the interacting species are driven out-of-equilibrium after impurity uptake, since reactivity is strongly quenched if a full thermal equilibration is assumed. More generally, our work points towards a novel mechanism for long-range electron transfer through an interplay between nuclear quantum delocalization within the confining medium and delocalized electronic dispersion forces acting on the two reactants.

Many-Body Dispersion

A broad range of approaches to many-body dispersion are discussed, including empirical approaches with multiple fitted parameters, augmented density functional-based approaches, symmetry adapted perturbation theory, and a supermolecule approach based on coupled cluster theory. Differing definitions of "body" are considered, specifically atom-based vs molecule-based approaches.

Quantum Features of Anionic Species He*⁻ and He₂*⁻ in Small He(N) Clusters

We present variational calculations on systems containing a few boson helium atoms attached to electronically excited atomic and molecular helium anions He*⁻ and He₂*⁻ and characterize their structures and energetics. Previously reported high-level ab initio results [Huber, S. E.; Mauracher, A. Mol. Phys. 2014, 112, 794] to describe the interactions between excited (metastable) anions and a neutral He atom have been employed. For the case of the atomic species He*⁻, the corresponding interaction with He suggests large anharmonicity effects due to the presence of a deep well of ∼17,500 cm⁻¹ at short distances, together with a more external shallow secondary well of ∼4 cm⁻¹, both supporting bound levels. Moreover, when a sum of pairwise interactions is assumed to describe the full PES corresponding to the presence of several neutral He atoms, geometrical constraints already predict the complete solvation of the anionic impurity by six helium atoms, giving rise to a bipyramidal structure. In turn, for the anisotropic weak interaction He-He₂*⁻, where the anionic dimer is considered as a rigid rotor, the obtained structures show the tendency of the helium atoms to pack themselves together and largely far away from the dopant, thereby confirming the heliophobic character of He₂*⁻.

Reactive scattering calculations for (87)Rb+(87)RbHe→Rb2((3)Σ(u)(+),v)+He from ultralow to intermediate energies

We investigate atom-diatom reactive collisions, as a preliminary step,in order to assess the possibility of forming Rb(2) molecules in their lowest triplet electronic state by cold collisions of rubidium atoms on the surface of helium nanodroplets [corrected]. A simple model related to the well-known Rosen treatment of linear triatomic molecules [N. Rosen, J. Chem. Phys. 1, 319 (1933)] in relative coordinates is used, allowing to estimate reactive probabilities for different values of the total angular momentum. The best available full dimensional potential energy surface [Guillon et al., J. Chem. Phys. 136, 174307 (2012)] is employed through the calculations. Noticeable values of the probabilities in the ultracold regime, which numerically fulfill the Wigner threshold law, support the feasibility of the process. The rubidium dimer is mainly produced at high vibrational states, and the reactivity is more efficient for a bosonic helium partner than when the fermion species is considered.

Solvation and spectral line shifts of chromium atoms in helium droplets based on a density functional theory approach

The interaction of an electronically excited, single chromium (Cr) atom with superfluid helium nanodroplets of various size (10 to 2000 helium (He) atoms) is studied with helium density functional theory. Solvation energies and pseudo-diatomic potential energy surfaces are determined for Cr in its ground state as well as in the y⁷P, a⁵S, and y⁵P excited states. The necessary Cr–He pair potentials are calculated by standard methods of molecular orbital-based electronic structure theory. In its electronic ground state the Cr atom is found to be fully submerged in the droplet. A solvation shell structure is derived from fluctuations in the radial helium density. Electronic excitations of an embedded Cr atom are simulated by confronting the relaxed helium density (ρ_He), obtained for Cr in the ground state, with interaction pair potentials of excited states. The resulting energy shifts for the transitions z⁷P ← a⁷S, y⁷P ← a⁷S, z⁵P ← a⁵S, and y⁵P ← a⁵S are compared to recent fluorescence and photoionization experiments.

Low-temperature rate coefficients for vibrational relaxation of (3)Σ(u)(+)Rb2 molecules by (3)He and (4)He atoms

We present quantum-scattering calculations of (4)He and (3)He colliding with (87)Rb2. For both helium isotopes, the elastic and inelastic rate coefficients are strongly influenced by the J = 1 partial wave. For the lighter isotope, a strong resonance feature of the J = 1 partial wave is responsible for an extremely efficient vibrational relaxation process. We also perform bound-state calculations of the Rb2He complex for even Rb permutation symmetry and nonzero total angular momentum. The global Rb2He (3)Σu(+) potential-energy surface used supports four bound states for (4)He and a single one for (3)He. We propose an analysis of the (87)Rb2(4)He spectrum separating the contributions of Rb2 rotation and helium motion.

Binding energies of the ground triplet state a(3) Σ(u)(+) of Rb2 and Cs2 in terms of the generalized Le Roy-Bernstein near-dissociation expansion

Formulae of Le Roy-Bernstein near-dissociation theory are derived in a general isotope-invariant form, applicable to any term in the rotational expansion of a diatomic ro-vibrational term value. It is proposed to use the generalized Le Roy-Bernstein expansion to describe the binding energies (ro-vibrational term values) of the ground triplet state a(3) Σ(u)(+) of alkali metal dimers. The parameters of this description are determined for Rb2 and Cs2 molecules. This approach gives a recipe to calculate the whole variety of the binding energies with characteristic accuracies from ∼1 × 10(-3) to 1 × 10(-2) cm(-1) using a relatively simple algebraic equation.