Structure and photoabsorption properties of cationic alkali dimers solvated in neon clusters

Potential Energy Surface and Bound States of Ne-Li2+(X2Σg+) van der Waals Complex Based on Ab Initio Calculations

Theoretical studies of the potential energy surface and vibrational bound states calculations were performed for the ground state of the Ne-Li2+(X2Σg+) van der Waals (vdW) complex. The intermolecular interactions were investigated by using an accurate monoconfigurational RCCSD(T) method and large basis sets (aug-cc-pVnZ, n = T, Q, 5), extrapolated to the complete basis set (CBS) limit. In turn, the obtained raw data from RCCSD(T)/CBS(Q5) calculations were numerically interpolated using the Morse + vdW model and the Reproducing Kernel Hilbert Space (RKHS) polynomial method to generate analytic expressions for the 2D-PES. The RKHS interpolated PES was then used to assess the bound states of the Ne-Li2+(X2Σg+) system through nuclear quantum calculations. By studying the aspect of the potential energy surface, the analysis sheds light on the behavior of the Ne-Li2+(X2Σg+) complex and its interactions between repulsive and attractive forces with other particles. By examining the vibrational states and wave functions of the system, the researchers were able to gain a better understanding of the behavior of the Ne-Li2+(X2Σg+) complex. The calculated radial and angular distributions for all even and odd symmetries are discussed in detail. We observe that the radial distributions exhibit a more complicated nodal structure, representing stretching vibrational behavior in the neon atom along its radial coordinate. For the highest bound states, the situation is very different, and the energies surpass the angular barrier.

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 the matrix shift on the vibrational frequency of ThO by DFT-D3 calculations

Benchmark calculations with a goal to find dispersion-corrected DFT-D3 methods suitable for a reliable estimation of matrix shifts on the vibrational frequency were carried out on the ThO molecule in three rare gas (Rg = Ne, Ar, and Kr) matrices. The matrices were modelled by the explicit approach, in which a single and a double shell of Rg atoms around ThO was considered. The selection of exchange-correlation functionals was based on test calculations on triatomic ThO⋯Rg models. The B3LYP, PBE0, CAM-B3LYP, and LC-ωPBE functionals were found to be the best suited for the estimation of matrix shifts. The single shell of Rg's around ThO accounted for a major part of the shifts; the addition of a second Rg shell resulted only in a minor improvement. Continuum solvation models considerably overestimated the effect of Rg matrices both when the whole matrix was treated by the model and when the first shell was treated explicitly and the rest with a continuum solvation model.

Nonadiabatic molecular dynamics of photoexcited Li2(+) Ne(n) clusters

We investigate the relaxation of photoexcited Li(2)(+) chromophores solvated in Ne(n) clusters (n = 2-22) by means of molecular dynamics with surface hopping. The simplicity of the electronic structure of these ideal systems is exploited to design an accurate and computationally efficient model. These systems present two series of conical intersections between the states correlated with the Li+Li(2s) and Li+Li(2p) dissociation limits of the Li(2)(+) molecule. Frank-Condon transition from the ground state to one of the three lowest excited states, hereafter indexed by ascending energy from 1 to 3, quickly drives the system toward the first series of conical intersections, which have a tremendous influence on the issue of the dynamics. The states 1 and 2, which originate in the Frank-Condon area from the degenerated nondissociative 1(2)Π(u) states of the bare Li(2)(+) molecule, relax mainly to Li+Li(2s) with a complete atomization of the clusters in the whole range of size n investigated here. The third state, which originates in the Frank-Condon area from the dissociative 1(2)Σ(u)(+) state of the bare Li(2)(+) molecule, exhibits a richer relaxation dynamics. Contrary to intuition, excitation into state 3 leads to less molecular dissociation, though the amount of energy deposited in the cluster by the excitation process is larger than for excitation into state 1 and 2. This extra amount of energy allows the system to reach the second series of conical intersections so that approximately 20% of the clusters are stabilized in the 2(2)Σ(g)(+) state potential well for cluster sizes n larger than 6.