Theoretical Study of Rb2 in HeN: Potential Energy Surface and Monte Carlo Simulations

Nonadiabatic Laser-Induced Alignment Dynamics of Molecules on a Surface

We demonstrate that a sodium dimer, Na_{2}(1^{3}Σ_{u}^{+}), residing on the surface of a helium nanodroplet, can be set into rotation by a nonresonant 1.0 ps infrared laser pulse. The time-dependent degree of alignment measured, exhibits a periodic, gradually decreasing structure that deviates qualitatively from that expected for gas-phase dimers. Comparison to alignment dynamics calculated from the time-dependent rotational Schrödinger equation shows that the deviation is due to the alignment dependent interaction between the dimer and the droplet surface. This interaction confines the dimer to the tangential plane of the droplet surface at the point where it resides and is the reason that the observed alignment dynamics is also well described by a 2D quantum rotor model.

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.

Structure, energetics, and spectroscopy of the K2+(X2Σ+g) interacting with the noble gas atoms Ar, Kr and Xe

The structure, energetic, and spectroscopy properties of the ionic system K₂⁺(X²Σ⁺_g) interacting with the noble gas atoms Argon, Krypton and Xenon are studied. The computations are done by an accurate ab initio approach based on the pseudo-potential technique, Gaussian basis sets, parameterized l-dependent polarization potentials and an analytic potential form for the K⁺Ar, K⁺Kr and K⁺Xe interactions. These interactions are added via the CCSD(T) potential taken from literature and fitted applying the analytical expression of Tang and Toennies. The application of the pseudo-potential approach reduces the number of active electrons of each complex to only one electron. The potential energy surfaces are analyzed on a large range of the Jacobi coordinates, R and θ. By the general interpolation outline based on the RKHS (Reproducing Kernel Hilbert Space) procedure, we have reproduced for each complex from our ab initio results the two-dimensional contour plots of an analytical potential. To evaluate the stability of each complex, we have determined from the potential energy surfaces the equilibrium distance (R_e), the well depth (D_e), the quantum energy (D₀), the zero-point-energy (ZPE) and the ZPE%. The results showed that the linear configurations, where the noble gas atom connected to the K₂⁺(X²Σ⁺_g) system, are the more stable.

Modelling interactions of cationic dimers in He droplets: microsolvation trends in HenK2+ clusters

We report the results of a detailed theoretical investigation of small K₂⁺-doped He clusters. The structural characteristics and stabilities of such cations are determined from ab initio electronic structure calculations at the MRCI+Q level of theory. The underlying interactions show a multireference character and such effects are analyzed. The interaction potentials are constructed employing an interpolation technique within the inverse problem theory method, while the nuclear quantum effects are computed for the trimers, their spatial arrangements are discussed, and information was extracted on the orientational anisotropy of the forces. We found that energetically the most stable conformer corresponds to linear arrangements that are taking place under large amplitude vibrations, with high zero-point energy. We have further looked into the behavior of higher-order species with various He atoms surrounding the cationic dopant. By using a sum of potentials approach and an evolutionary programming method, we analyzed the structural stability of clusters with up to six He atoms in comparison with interactions energies obtained from MRCI+Q quantum chemistry computations. Structures containing He_n motifs that characterize pure rare gas clusters, appear for the larger K₂⁺-doped He clusters, showing selective growth during the microsolvation process of the alkali-dimer cation surrounded by He atoms. Such results indicate the existence of local solvation microstructures in these aggregates, where the cationic impurity could get trapped for a short time, contributing to the slow ionic mobility observed experimentally in ultra-cold He-droplets.

Quantum effects on the stability of the He5 I2 van der Waals conformers

We present 15-dimensional quantum multiconfiguration time-dependent Hartree calculations of the vibrational levels of the He₅ I₂ van der Waals (vdW) complex employing an ab initio-based potential energy surface (PES). The energies and spatial features of such bound structures are analyzed, providing predictions on the structures and relative stabilities of its three lowest isomers. We found that the most stable isomer corresponds to all five He atoms encircling the I₂ molecule, indicating that in this case the anharmonic quantum effects do not stabilize the isomers involving a He atom in a linear configuration as reported previously for the smaller He_N I₂ systems. Such finding provides information on the overall structuring of the finite-size-solvent systems, highlighting the intriguing interplay between weak intermolecular interactions and quantum effects. © 2019 Wiley Periodicals, Inc.

Theoretical Study of Cationic Alkali Dimers Interacting with He: Li2+-He and Na2+-He van der Waals Complexes

We present a theoretical study on the potential energy surface and bound states of He-A₂⁺ complexes, where A is one of the alkali Li or Na atoms. The intermolecular interactions were systematically investigated by high-level ab initio electronic structure computations, and the corresponding raw data were then employed to reproduce accurate analytical expressions of the potential surfaces. In turn, we used these potentials to evaluate bound configurations of the trimers from nuclear quantum calculations and to extract information on the effect of orientational anisotropy of the forces and the interplay between repulsive and attractive interaction within the potential surfaces. The spatial features of the bound states are analyzed and discussed in detail. We found that both systems are going under large amplitude stretching and bending motions with high zero-point energies. Despite the large differences in the potential well-depths, the correct treatment of nuclear quantum effects provides insights on the effect of different strength of the ionic interaction on the spectral structure of such cationic alkali van der Waals complexes, related to the mobility of ions and the formation of cold-molecules in He-controlled environments.

An ab initio study of the ground and excited electronic states of the methyl radical

The ground and some excited electronic states of the methyl radical have been characterized by means of highly correlated ab intio techniques. The specific excited states investigated are those involved in the dissociation of the radical, namely the 3s and 3p_z Rydberg states, and the A₁ and B₁ valence states crossing them, respectively. The C-H dissociative coordinate and the HCH bending angle were considered in order to generate the first two-dimensional ab initio representation of the potential surfaces of the above electronic states of CH₃, along with the nonadiabatic couplings between them. Spectroscopic constants and frequencies calculated for the ground and bound excited states agree well with most of the available experimental data. Implications of the shape of the excited potential surfaces and couplings for the dissociation pathways of CH₃ are discussed in the light of recent experimental results for dissociation from low-lying vibrational states of CH₃. Based on the ab initio data some predictions are made regarding methyl photodissociation from higher initial vibrational states.

Vibrational Calculations of Higher-Order Weakly Bound Complexes: The He3,4I2 Cases

The structure and relative stability of higher-order He3,4I2 clusters are investigated by carrying out full-dimensional quantum calculations within the multiconfiguration time-dependent Hartree framework. The full interaction between the I2 molecule and the He atoms is based on analytical three-body ab initio He-I2 potentials obtained from high level ab initio calculations plus the He-He interaction. The low-lying minima on the potential surfaces are found to be very close in energy with the He atoms in a ring encircling the dopant for the global minimum structure, while for the local minima one or two of the He atoms prefer the linear arrangements along the I2-axis. Such classical description on the basis of the potential energy landscape is corrected by including anharmonic quantum effects, present in highly floppy systems, by carrying out full dimensional quantum calculations. The potential energy operator was constructed by natural potential fits, while a mode combination scheme was employed to optimize the computational cost of the improved relaxation calculations. The obtained results predict the relative stability of the He3,4I2 isomers at zero temperature and provide benchmark data on binding energies and structural properties of these van der Waals systems. The (2,1) and (2,2), involving two He atoms in the T-shape and one or two He atoms in the linear configurations, respectively, are found to be the most stable isomers, although extremely close in energy with the (3,0) and (4,0) ones as predicted by classical optimizations. Comparison with experimental data on similar systems at low temperatures is also discussed. This analysis indicates once more the importance of quantum delocalization and the need of accurate quantum-mechanical treatments to characterize such doped helium nanosystems.

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.

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.

Energetics and solvation structure of a dihalogen dopant (I2) in (4)He clusters

The energetics and structure of small HeNI2 clusters are analyzed as the size of the system changes, with N up to 38. The full interaction between the I2 molecule and the He atoms is based on analytical ab initio He-I2 potentials plus the He-He interaction, obtained from first-principle calculations. The most stable structures, as a function of the number of solvent He atoms, are obtained by employing an evolutionary algorithm and compared with CCSD(T) and MP2 ab initio computations. Further, the classical description is completed by explicitly including thermal corrections and quantum features, such as zero-point-energy values and spatial delocalization. From quantum PIMC calculations, the binding energies and radial/angular probability density distributions of the thermal equilibrium state for selected-size clusters are computed at a low temperature. The sequential formation of regular shell structures is analyzed and discussed for both classical and quantum treatments.

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

A new full-dimension potential energy surface of the three-body He-Rb₂(³Σ(u)(+)) complex and a quantum study of small (⁴He)(N)-Rb₂(³Σ(u)(+)) clusters, 1 ≤ N ≤ 4, are presented. We have accurately fitted the ab initio points of the interaction to an analytical form and addressed the dopant's vibration, which is found to be negligible. A Variational approach and a Diffusion Monte Carlo technique have been applied to yield energy and geometric properties of the selected species. Our quantum structure calculations show a transition in the arrangements of the helium atoms from N = 2, where they tend to be separated across the diatomic bond, to N = 4, in which a closer packing of the rare gas particles is reached, guided by the dominance of the He-He potential over the weaker interaction of the latter adatoms with the doping dimer. The deepest well of the He-Rb₂ interaction is placed at the T-shape configuration, a feature which causes the dopant to be located as parallel to the helium "minidroplet". Our results are shown to agree with previous findings on this and on similar systems.

Weakly bound finite systems: (⁴He)N-Rb₂(³Σu), clustering structures from a quantum Monte Carlo approach

We report here ((4)He)(N)-Rb(2)((3)Σ(u)) complexes, 2 ≤ N ≤ 20, analysed through a quantum diffusion Monte Carlo stochastic approach. The calculations show that the spin stretched dimer molecule is bound outside the pure He sub-complex, due to the stronger He-He potential as compared with the He-Rb(2) interaction, while the rare gas atom moiety presents, in turn, a shell-like structure with ten He adatoms completing the first shell. Our results agree with previous findings on this and similarly weakly interacting systems.