Next Article in Journal
Pharmacoinformatics-Based Approach for Uncovering the Quorum-Quenching Activity of Phytocompounds against the Oral Pathogen, Streptococcus mutans
Previous Article in Journal
A Multiple Stimuli–Responsive Ag/P/S Complex Showing Solvochromic and Mechanochromic Photoluminescence
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Laboratory of Interfaces and Advanced Materials, Physics Department, Faculty of Sciences of Monastir, Avenue de l’Environnement, Monastir 5019, Tunisia
2
Department of Physics, College of Science and Humanities in Al-Kharj, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
3
Mathematics and Natural Sciences Department, School of Arts and Sciences, American University of Ras Al Khaimah, Ras Al-Khaimah 10021, United Arab Emirates
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(14), 5512; https://doi.org/10.3390/molecules28145512
Submission received: 23 May 2023 / Revised: 10 July 2023 / Accepted: 17 July 2023 / Published: 19 July 2023
(This article belongs to the Section Computational and Theoretical Chemistry)

Abstract

:
We report a computational study of the potential energy surface (PES) and vibrational bound states for the ground electronic state of L i 2 + K r . 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 L i 2 + K r 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 K r 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.

1. Introduction

Over the past few decades, rare gas media and matrices have attracted great attention as fascinating environments for investigating the spectroscopy, structure, and dynamics embedded inside or deposited on their surfaces. Rare gas atoms are renowned as non-reactive solvents [1,2,3] and provide a model system for both theoretical and experimental investigations involving a large number of atoms, representing a high degree of freedom. Furthermore, as their condensates exhibit transparency across a broad range of radiation energies, inert gas matrices offer the advantage of selectively exciting chromophores embedded within them. Considering the case of alkali atoms, these matrices have proven useful in studying the spectroscopy, photo-association, and optical absorption spectra of alkali atoms [2,4,5]. Matrix isolation spectroscopy (MIS) is a widely utilized technique in various research fields, particularly for investigating alkali atoms in solid matrices of noble gases at cryogenic temperatures [6], as well as studying the high-spin states of alkali molecules [7,8,9,10].
The exceptional characteristics of rare gas atoms, combined with the growing interest in studying alkali atoms following the observation of the photo-association of cold alkali atoms and their Bose–Einstein condensation [4], have made the interaction between alkali atoms and rare gas atoms a prominent research theme in numerous studies. Notably, both experimental and theoretical works [4,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29] have been dedicated to exploring the interaction between alkali neutral/ionic dimers and rare gas atoms. The primary focus of these studies is calculating the potential energy surfaces (PES) of dimers with a single rare gas atom. Subsequently, these surfaces are fitted with suitable analytical forms, which are then utilized for dynamic studies such as bound state calculations, as well as for structural investigations, such as the geometric minimization of large-sized clusters to examine the micro-solvation of the embedded dimers.
From an experimental perspective, recent work by Kristensen et al. [21] introduced the use of Coulomb explosion induced by an intense fs laser pulse as a means of studying alkali dimers ( R b 2 , L i 2 , N a 2 and K 2 ) on H e droplets. They demonstrated that Coulomb explosion enables the determination, within a single measurement, of whether alkali dimers, A k 2 , on the surface of H e nanodroplets, are formed in either the X 1 g + ground state or in the lowest-lying triplet state a 3 u + . Fuchs et al. [19] experimentally investigated the collisional cross-sections of lithium dimer L i 2   in selected vibrational states with H e and K r atoms. An der Lan et al. [13] presented high-resolution mass spectra of N a + H e n , K + H e n , N a 2 + H e n and K 2 + H e n , formed via the electronionization of doped helium droplets. They observed two distinct anomalies in N a 2 + H e n clusters at n = 2 and n = 6.
On the theoretical side, Guillon et al. [11] recently performed diffusion and path integral quantum Monte Carlo studies of the rubidium dimer in its ground triplet state in a helium environment [11]. The focus of their work was the influence of helium atoms on the rotational motion of the R b 2 dimer. Alharzali et al. [12,22], Bodo et al. [17,23,29], and Marinetti et al. [26] focused on the structure and dynamic properties of cationic dimers ( L i 2 + ,   N a 2 +   a n d   K 2 + ) interacting with helium atoms using RCCSDT and Post Hartree–Fock approaches. They found that for all three species, the helium atom preferred linear attachment to the cationic dimer rather than a T-shaped configuration. In the case of larger clusters, they observed the solvation of cationic dimers within the helium clusters. Douady et al. [20] investigated the solvation of ionic sodium dimers in argon clusters. They confirmed that the cationic dimer enters the clusters rather than resides at the surface. This behavior is attributed to the relatively strong interaction between the closed-shell neon atoms and the cationic dimer. Zanuttini et al. [15] employed a pseudopotential technique and molecular dynamics with a surface hopping approach to investigate the structure and optical absorption of L i 2 + , N a 2 + , and K 2 + alkali dimers in neon clusters. They concluded that the modification of the PES by surrounding neon atoms indicated the strong confinement of the lowest energy states of the three cationic dimers. Saidi et al. [18] explored the structure and stability of the lithium dimer with xenon atoms, while Ghanmi et al. [25] performed a structural, energetic, and spectroscopic study of the potassium cationic dimer interacting with rare gas atoms, such as A r ,   K r , and X e .
In this paper, we present a computational investigation of the structure and dynamics of the lithium cationic L i 2 + dimer interacting with a K r atom. Section 3 provides an overview of the computational method details employed in this study, including the basis sets and extrapolation schemes used to calculate the interaction energies of the L i 2 + K r complex. This potential energy surface is, then, employed in the calculation of bound state levels after being reproduced via the RKHS method [30]. The results are presented and discussed in Section 2. The used methodologies are detailed in Section 3.

2. Results and Discussion

The resulting interaction energies were employed to calculate the spectroscopic constants R e   and D e for each configuration corresponding to angles ranging from 0 °   t o   90 ° . The found equilibrium distances ( R e ) and the well depths ( D e ) for the L i 2 + K r using the RCCSD(T) method and different basis are summarized in Table 1.
From Table 1, it can be observed that the well depth of the L i 2 + K r complex at all geometries increases when moving from AVQZ to CBS[Q5]. Taken the example of their collinear arrangement ( θ   = 0 ° )   D e varies from 1730 cm−1 to 1737 cm−1 for the RCCSD(T)/AVQZ and RCCSD(T)/AV5Z calculations, respectively. However, upon computing the RCCSD(T)/CBS[Q5] using the RCCSD(T)/AVQZ/AV5Z results, the energies increase to 1745 cm−1. This indicates that the RCCSD(T)/CBS[Q5] energies are lower than those obtained from the RCCSD(T)/AVQZ and RCCSD(T)/AV5Z estimations.
A comparison of the RCCSD(T)/AVQZ, RCCSD(T)/AV5Z, and RCCSD(T)/CBS[Q5] extrapolations is depicted in Figure 1, considering all angular orientations. It can be observed that despite a slight shift between the three curves, the energies obtained from the RCCSD(T)/AVQZ/AV5Z calculations are lower than those from the RCCSD(T)/CBS[Q5] extrapolation. Therefore, the RCCSD(T)/CBS[Q5] energies will be considered as reference data in the subsequent sections.
It is worth noting that as the angle increases from θ   =   0 ° to θ   =   90 ° , the well depth decreases from D e   = 1745 cm−1 to D e   = 317 cm−1. This implies that the linear configuration is the most stable. This arrangement can be explained by considering that the L i 2 + dimer in its ground state can be approximated as two closed-shell Li+ cores with an electron cloud located between them. Consequently, the short-range repulsion between the electron and the krypton atom, combined with the attraction between the two Li+ cores and the closed-shell K r atom, favors the positioning of the krypton atom at the ends of the dimer.
It is important to mention that the interaction of the lithium dimer with noble gas atoms has been extensively studied recently due to its simplicity and significance. Studies conducted by Zanuttini et al. [15] and Alharzali et al. [22] have contributed to this research area. In Table 2, we have grouped the findings from these studies with the current results obtained for L i 2 + K r . It is evident that the repulsive interactions decrease in their impact compared to the attractive long-range interaction forces when moving from L i 2 + H e to L i 2 + K r . In fact, the well depths in their linear configuration increase from 380 cm−1 for L i 2 + H e to 700 cm−1 for L i 2 + N e and further to 1745 cm−1 for L i 2 + K r . Hence, the interactions of the three rare gas atoms ( H e ,   N e ,   K r ) with the lithium dimer become increasingly more attractive, exhibiting larger attractive cores as one moves from helium to krypton.
The RKHS potential curves together with the ab-initio RCCSD(T)/CBS[Q5] interaction energies along the R coordinate for each θ angle from 0 ° to 90 ° are illustrated in Figure 2. A good agreement is observed between the RKHS potential and the RCCSD(T)/CBS[Q5] calculations for all orientations.
To further verify the quality of the fit we compute the relative error E ( % ) between the original ab initio RCCSD(T)/CBS[Q5] and the values of the RKHS potential for all orientations. The results are presented in Table 3. We note that the relative error does not exceed 0.175%.
The angular minimum energy path for all configurations of the L i 2 + K r complex, obtained from the RKHS and the RCCSD(T)/CBS[Q5] potentials, are depicted in Figure 3. Several observations could be concluded from the behavior of the angular minimum energy path. Firstly, the RKHS fitting method demonstrates good performance in describing the ab-initio RCCSD(T)/CBS[Q5] PES for all orientations. Secondly, as the configuration changes from linear to T-shaped, the attractive effect decreases compared to the increasing potential values. Furthermore, there is a non-equally spaced energy between successive orientations, and the differences become more pronounced at angles between 30° and 60°. This phenomenon can be attributed to the emergence of repulsive effects, which become more important than the attractive ones as the K r atom approaches the center of mass of the ionic dimer L i 2 + .
Figure 4 displays a two-dimensional contour plot in the ( θ ,   R ) plane, representing the fitted RCCSD(T)/CBS[Q5] potential of L i 2 + K r . In this plot, we observe the presence of two symmetrical minima, which correspond to the linear orientations. These minima have a well depth of 1746 cm⁻1 located at an internuclear distance R = 4.17 Å.
Based on the parity of the quantum number j and the total wave function, with respect to the exchange of the two identical L i + atoms, the results will be divided into two types (even/odd). It is expected to obtain a significant number of bound states when we consider the energy of the linear arrangement, 1746 cm−1. In our calculations, we obtained 272/260 states for even and odd symmetry, respectively. In Table 4, we list only the energies of the lowest bound states and some selected states that lie below, around, and above the T-shaped 90° barrier well for both even and odd symmetry. The other values are reported in the Supplementary Materials.
From the full table, we observe that up to n = 78, the even and odd vibrational states are fully degenerated. However, the higher-lying even and odd parity states exhibit different energies, and the energy difference becomes more pronounced for levels above the potential barrier in the T-shaped configuration. Furthermore, Table 4 shows that the energy of the lowest vibrational level is −1666 cm−1, with a zero-point energy ( Z P E ) of about 88 cm−1, which corresponds to only 5% of the well depth.
Figure 5, Figure 6 and Figure 7 illustrate the radial and angular distributions for selected states lying below, around, and above the T-shaped 90° barrier for both even and odd symmetries, respectively.
For the lowest bound states wave functions, such as n = 0–5, there is a minimal angular population in the region of the 90° barrier well, indicating a strong localization of the   K r atom in the collinear geometric configuration. The ground state wave function exhibits nodal structures, with angular distributions symmetrically positioned at θ   =   0 ° and 180°, while the radial distribution peaks at R = 4.3 Å, which is close to the equilibrium intermolecular distance of the linear geometry. The low-lying excited states can be attributed to the stretching and bending motions of the K r atom. The bending motions correspond to the n = 1 and n = 2 levels, as evidenced by the presence of nodes in their angular distributions. The n = 3 state clearly displays nodes in its radial distribution, indicating stretching motions.
For the bound states that lie around the 90° barrier well (just below and above it), such as n = 114–119, Figure 6 illustrates more complex nodal patterns in both the radial and angular distributions. These distributions extend to larger R distances and configurations, and some population appears in the region of the T-shaped barrier. This behavior can be associated with the mixed stretch-bending vibrational “mode” of the   K r atom along the radial and angular coordinates. Additionally, noticeable differences between even and odd states become apparent in both radial and angular probability distributions for levels n = 114 and above.
For the highest bound states, as depicted in Figure 7, the shape of the distributions becomes even more complicated, with extended angular distributions and prominent differences between even and odd parity states.

3. Computational Methods

3.1. Potential Energy Surface

3.1.1. Ab Initio Calculations

To describe the intermolecular interactions between the krypton atom and the diatomic molecule L i 2 + , Jacobi coordinates ( r ,   R ,   θ ) were employed, as illustrated in Figure 8. Here, r represents the equilibrium distance of L i 2 + , R denotes the distance between the center of mass of L i 2 + and the krypton atom, and θ represents the Jacobi angle between the vectors r and R .
The calculations were carried out by keeping the diatomic L i 2 + frozen at its experimental equilibrium distance, r e = 3.11 Å, obtained via pulsed optical–optical double resonance spectroscopy [16]. The internuclear distance R was varied from 2 Å to 20 Å with an irregular step, while the Jacobi angle θ was varied from 0° to 90° with a step of 10°.
In this present work, the Molpro 2010 Package [31] was used to perform all ab initio calculations. The potential energy surface of L i 2 + K r was calculated using the restricted Hartree–Fock calculation followed by a single-reference restricted open-shell coupled cluster method with single, double, and no-iterative triple excitations RCCSD(T) method [32,33]. This method is known for its high level of accuracy in describing electronic correlation effects.
To assess the validity of the RCCSD(T) method employed, correlation factors   D 1 and T 1 were calculated for all investigated geometries. The T 1 and   D 1 diagnostics, as defined in reference [34], were used for this purpose. The T 1 diagnostic is given by T 1 = i o c c a v i r t i a 2 n , where n is the number of electrons, and t i a represents the single excitation amplitudes. The   D 1 diagnostic is defined as D 1 C C S = T 2 , where T is the Euclidean norm of the matrix T   calculated from the CCSD wave function. The results of the calculations are summarized in Table 5.
The obtained results indicate that both, T 1 and D 1 factors, are below the threshold values of 0.02 and 0.025, respectively, as suggested by Lee and Taylor [34]. This confirms the suitability of the mono-configurational approach RCCSD(T) chosen for this study.
For both L i and K r atoms, the aug-cc-pVn (n = 4, 5) Z basis sets [35] were utilized. The energies obtained were then extrapolated to the complete basis set (CBS) limit using a two-parameter expression [36] applied to the correlation energies as follows:   E n =   E C B S + A n 3 . Here, n = 4 and 5 correspond to AVQZ and AV5Z basis sets, respectively, E n is the computed total energy corresponding to that basis set, E C B S is the CBS extrapolated energy, and A is a fitting parameter. The energies from the AVQZ and AV5Z basis set calculations were employed for the extrapolation.
In the calculations of the internuclear interaction energies V ( r e ,   R ,   θ ) between K r and L i 2 + , the standard counterpoise method of Boys and Bernardi [37] was used to correct the basis set superposition error (BSSE) at all configurations. This correction is applied according to the equation: V r e ,   R ,   θ = E L i 2 + K r r e , R , θ [ E L i 2 + ( r e ,   R ,   θ ) + E K r ( r e ,   R ,   θ ) ] , where E L i 2 + K r ( r e ,   R ,   θ ) represents the total energy of the complex, and E L i 2 + ( r e ,   R ,   θ ) and E K r ( r e ,   R ,   θ ) are the energies of the L i 2 + and K r monomers, respectively, calculated in the full basis set of the system.

3.1.2. Analytical Representation of the Ab Initio Surface

In order to conduct dynamic investigations, an analytical representation of the ab-initio potential energy surface (PES) is necessary. In this study, the 2D-PES of L i 2 + K r was interpolated using the reproducing kernel Hilbert space (RKHS) procedure developed by Ho and Rabitz [30]. The two-dimensional potential function is given by the following:
V r e , R , θ = i = 1 N R j = 1 N θ v i j q 1 2,3 ( R i , R ) q 2 ( y j , y )
V R K H S r e , R , θ = i = 1 N R j = 1 N θ V i j q 1 n , m R i , R q 2 y j , y
In the above equations, y   =   c o s θ , and N R and N θ represent the number of calculated ab initio energies in the R and θ coordinates, respectively. The v i j coefficients are determined by solving Equation (1), where V ( r e , R i , θ j ) represents the ab initio RCCSD(T)/CBS[Q5] energy at each ( R i ,   θ j , r e ) grid point.
The one-dimensional distance-like q 1 n , m and angle-like q 2 reproducing kernels are expressed as follows:
q 1 n , m x , x = n 2 x > m + 1 B m + 1 , n 2 F 1 n + 1 , m + 1 ; n + m + 1 ; x < x >
q 2 y , y = l 2 l + 1 2 P l ( y ) P l ( y )
where x > and x < refer to the maximum and minimum values of x and x , respectively. B   a n d   2 F 1 represent the Beta and Gauss hypergeometric functions [38], respectively. P l denotes the Legendre polynomials with l = 0, 2, 4, 6, 8, 10, 12.
Based on the previous work by Alharzali et al. [12,24] for similar interactions ( L i 2 + H e , N a 2 + H e and K 2 + H e ), it is assumed that the 2D-PES V ( R ,   θ ) of L i 2 + K r is a smooth function with derivatives up to the second order, n = 2, in both R and θ . Additionally, to account for the dominant dispersion interaction between the K r atom and the ionic dimer L i 2 + , a weighting factor w x = x m with m = 3 is introduced.

3.2. Bound States Calculation

The calculations of bound states energies and wave functions are crucial steps for studying dynamics, such as vibrational predissociation and photoionization. In this study, variational quantum bound states calculations were performed using the fitted RKHS potential energies. The bound states energies were obtained by diagonalizing the Hamiltonian expressed in Jacobi coordinates, as follows:
H ^ = 2 2 μ 1 2 2 R + j ^ 2 2 μ 2 r e 2 + l ^ 2 2 μ 1 R 2 + V ( r e , R , θ )
In the above equation, 1 μ 1 = 1 m K r + 1 2 m L i 2 + and 1 μ 2 = 1 m L i + 1 m L i are the reduced masses of the L i 2 + K r complex and the L i 2 + dimer, respectively, where m L i and m K r   are the atomic masses of L i and K r atoms. j ^ and l ^ represent the angular momenta associated with the vectors r (rotational momentum of the dimer) and R (orbital angular momentum), respectively. The sum of these angular momenta gives the total angular momentum J ^ , which was taken as zero in these calculations. V ( r e , R , θ ) represents the fitted 2D-RKHS potential, where r is the fixed equilibrium bond length of the diatomic L i 2 + .
In the framework of zero total angular momentum J = 0 , a product of radial and angular basis functions is used to represent the Hamiltonian. For the angular coordinate, the employed basis function is an orthonormalized Legendre polynomial P j c o s ( θ ) , where j ranges up to 40 for even and odd symmetry, respectively. For the radial coordinate, a discrete variable representation (DVR) basis set, based on particle in a box eigenfunctions [39], was employed. A 300 points DVR was used over a range from R   = 2.5 to 20 Å, and a convergence criterion of 10 6 was established.

4. Conclusions

In this study, we conducted structural and dynamic investigations of the lithium cationic dimer, L i 2 + , interacting with the K r atom. The two-dimensional potential energy surface was computed using the RCCSD(T) method and the aug-cc-pVnZ (n = 4, 5) basis sets, and then extrapolated to its CBS[Q5] limit. Subsequently, the RKHS method was employed for numerical interpolation to generate the RCCSD(T)/CBS[Q5] analytical potentials. Remarkably, this interpolation method accurately matched the numerical curves for all orientations used to determine the two main spectroscopic constants, R e and D e , for each configuration. The analysis of the results revealed that the linear configuration, where the krypton atom is linearly attached to the lithium dimer L i 2 + is found to be more stable than the T-shaped configuration.
The fitted 2D-RKHS potential was employed to calculate the vibrational bound state energies, resulting in a significant number of bound states as expected. It is worth noting the limited contribution of the zero-point energy ( Z P E ), which accounts for only 5% of the well depth. Additionally, the system exhibits very large amplitude motions, both as stretching and bending modes in its lower states, but rapidly transformed into non-regular features of their wave function in the excited states.
The obtained bound states, along with the potential energy surfaces, have various applications. They can be utilized to investigate vibrational predissociation or photoionization processes of L i 2 K r complex. Furthermore, an analytical fitting of the obtained 2D-PES can be employed to explore the micro-solvation process of the lithium cationic dimer in a krypton matrix. These valuable results could assist in experimental settings and help in the interpretation of observations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28145512/s1, Table S1: Vibrational energies (in cm−1) of all bound states ( J = 0) of the L i 2 + K r complex.

Author Contributions

Conceptualization, S.S. and J.D.; methodology, S.S., J.D. and N.M.; investigation, H.B., S.S. and J.D.; data curation, N.M.; writing—original draft preparation, S.S. and J.D.; writing—review and editing, S.S., J.D. and H.B.; project administration, S.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia (Project No. IF-PSAU-2021/01/18968).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IF-PSAU-2021/01/18968).

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not applicable.

References

  1. Liu, Q.; Wang, J.K.; Zewail, A.H. Femtosecond dynamics of dissociation and recombination in solvent cages. Nature 1993, 364, 427–430. [Google Scholar] [CrossRef]
  2. Apkarian, V.A.; Schwentner, N. Molecular Photodynamics in Rare Gas Solids. Chem. Rev. 1999, 99, 1509–1514. [Google Scholar] [CrossRef]
  3. Greenblatt, B.J.; Zanni, M.T.; Neumark, D.M. Photodissociation of I 2 A r n Clusters Studied with Anion Femtosecond Photoelectron Spectroscopy. Science 1997, 276, 1675–1678. [Google Scholar] [CrossRef]
  4. Kokoouline, V.; Dulieu, O.; Kosloff, R.; Masnou-Seeuws, F. Mapped Fourier methods for long-range molecules: Application to perturbations in the Rb2( O u + ) photoassociation spectrum. J. Chem. Phys. 1999, 110, 9865–9876. [Google Scholar] [CrossRef] [Green Version]
  5. Weyhmann, W.; Pipkin, F.M. Optical Absorption Spectra of Alkali Atoms in Rare-Gas Matrices. Phys. Rev. 1965, 137, A490. [Google Scholar] [CrossRef]
  6. Braggio, C.; Calabrese, R.; Carugno, G.; Fiscelli, G.; Guarise, M.; Khanbekyan, A.; Noto, A.; Passante, R.; Rizzuto, L.; Ruoso, G.; et al. Spectroscopy of Alkali Atoms in Solid Matrices of Rare Gases: Experimental Results and Theoretical Analysis. Appl. Sci. 2022, 12, 6492. [Google Scholar] [CrossRef]
  7. Rodríguez-Cantano, R.; López-Durán, D.; González-Lezana, T.; Delgado-Barrio, G.; Villarreal, P.; Yurtsever, E.; Gianturco, F.A. Spin-polarized Rb2 interacting with bosonic He atoms: Potential energy surface and quantum structures of small clusters. J. Phys. Chem. A. 2012, 116, 2394–2404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Higgins, J.; Callegari, C.; Reho, J.; Stienkemeier, F.; Ernst, W.E.; Gutowski, M.; Scoles, G. Helium Cluster Isolation Spectroscopy of Alkali Dimers in the Triplet Manifold. J. Phys. Chem. A. 1998, 102, 4952–4965. [Google Scholar] [CrossRef]
  9. Schulz, C.P.; Claas, P.; Schumacher, D.; Stienkemeier, F. Formation and Stability of High-Spin Alkali Clusters. Phys. Rev. Lett. 2004, 92, 013401–013404. [Google Scholar] [CrossRef]
  10. Vongehr, S.; Kresin, V.V. Unusual pickup statistics of high-spin alkali agglomerates on helium nanodroplets. J. Chem. Phys. 2003, 119, 11124–11129. [Google Scholar] [CrossRef]
  11. Guillon, G.; Zanchet, A.; Leino, M.; Viel, A.; Zillich, R.E. Theoretical Study of Rb2 in HeN: Potential Energy Surface and Monte Carlo Simulations. J. Phys. Chem. A. 2011, 115, 6918–6926. [Google Scholar] [CrossRef] [PubMed]
  12. Alharzali, N.; Rodríguez-Segundo, R.; Prosmiti, R. Modelling interactions of cationic dimers in He droplets: Microsolvation trends in Hen K 2 + clusters. Phys. Chem. Chem. Phys. 2021, 23, 7849–7859. [Google Scholar] [CrossRef]
  13. An der Lan, L.; Bartl, P.; Leidlmair, C.; Jochum, R.; Denifl, S.; Echt, O.; Scheier, P. Solvation of Na+, K+ and their dimers in helium. Chem. Eur. J. 2012, 18, 4411–4418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Prosmiti, R.; Delgado-Barrio, G.; Villarreal, P.; Yurtsever, E.; Coccia, E.; Gianturco, F.A. Structuring a Quantum Solvent around a Weakly Bound Dopant: The HeCs2(3u) Complex. J. Phys. Chem. A. 2009, 113, 14718–14729. [Google Scholar] [CrossRef]
  15. Zanuttini, D.; Douady, J.; Jacquet, E.; Giglio, E.; Gervais, B. Structure and photoabsorption properties of cationic alkali dimers solvated in neon clusters. J. Chem. Phys. 2010, 133, 174503. [Google Scholar] [CrossRef]
  16. Bernheim, R.A.; Gold, L.P.; Tipton, T. Rydberg states of 7Li2 by pulsed optical–optical double resonance spectroscopy: Molecular constants of 7 L i 2 + . J. Chem. Phys. 1983, 78, 3635–3646. [Google Scholar] [CrossRef]
  17. Bodo, E.; Gianturco, F.A.; Yurtsever, E.; Yurtsever, M. Neutral and ionic dopants in helium clusters: Interaction forces for the Li2(a3 u + ) − He and L i 2 + (X2 g + ) − He complexes. Mol. Phys. 2005, 103, 3223–3231. [Google Scholar] [CrossRef] [Green Version]
  18. Saidi, S.; Ghanmi, C.; Hassen, F.; Berriche, H. Ab initio Study of the Potential Energy Surface and Stability of the L i 2 + (X2 g + ) Alkali Dimer in Interaction with a Xenon Atom. Prog. Theor. Chem. Phys. 2012, 26, 321–330. [Google Scholar] [CrossRef]
  19. Fuchs, M.; Toennies, J.P. Scattering of highly vibrationally excited Li2 from He and Kr. J. Chem. Phys. 1986, 85, 7062–7076. [Google Scholar] [CrossRef]
  20. Douady, J.; Jacquet, E.; Giglio, E.; Zanuttini, D.; Gervais, B. Solvation of N a 2 + in Arn clusters. I. Structures and spectroscopic properties. J. Chem. Phys. 2008, 129, 184303–184313. [Google Scholar] [CrossRef]
  21. Kristensen, H.H.; Kranabetter, L.; Schouder, C.A.; Stapper, C.; Arlt, J.; Mudrich, M.; Stapelfeldt, H. Quantum-State-Sensitive Detection of Alkali Dimers on Helium Nanodroplets by Laser-Induced Coulomb Explosion. Phys. Rev. Lett. 2022, 128, 93201. [Google Scholar] [CrossRef] [PubMed]
  22. Alharzali, N.; Berriche, H.; Villarreal, P.; Prosmiti, R. Theoretical Study of Cationic Alkali Dimers Interacting with He: L i 2 + He and N a 2 + He van der Waals Complexes. J. Phys. Chem. A 2019, 123, 7814–7821. [Google Scholar] [CrossRef] [PubMed]
  23. Bodo, E.; Sebastianelli, F.; Gianturco, F.A.; Yurtsever, E.; Yurtsever, M. Ab initio quantum dynamics with very weak van der Waals interactions: Structure and stability of small L i 2 + (X1 g + ) − (He)n clusters. J. Chem. Phys. 2004, 120, 9160–9166. [Google Scholar] [CrossRef] [PubMed]
  24. Grebenev, S.; Toennies, J.P.; Vilesov, A.F. Superfluidity within a small helium-4 cluster: The microscopic andronikashvili experiment. Science 1998, 279, 2083–2086. [Google Scholar] [CrossRef]
  25. Ghanmi, C.; Nakbi, H.; Al-Qarni, H.J.; Alharzali, N.; Berriche, H. Structure, energetics, and spectroscopy of the K 2 + (X2 g + ) interacting with the noble gas atoms Ar, Kr and Xe. J. Mol. Graph. Model. 2023, 120, 108413–108418. [Google Scholar] [CrossRef]
  26. Marinetti, F.; Uranga-Piña, L.; Coccia, E.; López-Durán, D.; Bodo, E.; Gianturco, F.A. Microsolvation of Cationic Dimers in 4He Droplets: Geometries of A 2 + (He)N (A = Li, Na, K) from Optimized Energies. J. Phys. Chem. A 2007, 111, 12289–12294. [Google Scholar] [CrossRef] [Green Version]
  27. Stienkemeier, F.; Ernst, W.E.; Higgins, J.; Scoles, G. On the use of liquid helium cluster beams for the preparation and spectroscopy of the triplet states of alkali dimers and other weakly bound complexes. J. Chem. Phys. 1995, 102, 615–617. [Google Scholar] [CrossRef]
  28. Albertini, S.; Martini, P.; Schiller, A.; Schöbel, H.; Ghavidel, E.; Ončák, M.; Echt, O.; Scheier, P. Electronic transitions in R b 2 + dimers solvated in helium. Theor. Chem. Acc. 2021, 140, 29–40. [Google Scholar] [CrossRef]
  29. Bodo, E.; Yurtsever, E.; Yurtsever, M.; Gianturco, F.A. Ionic dimers in He droplets: Interaction potentials for L i 2 + He, N a 2 + He, and K 2 + He and stability of the smaller clusters. J. Chem. Phys. 2006, 124, 074320–074332. [Google Scholar] [CrossRef]
  30. Ho, T.S.; Rabitz, H. A general method for constructing multidimensional molecular potential energy surfaces from ab initio calculations. J. Chem. Phys. 1996, 104, 2584–2597. [Google Scholar] [CrossRef]
  31. Werner, H.J.; Knowles, P.J.; Knizia, G.; Manby, F.R.; Schütz, M. Molpro: A general-purpose quantum chemistry program package. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 242–253. [Google Scholar] [CrossRef]
  32. Čížek, J. On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods. J. Chem. Phys. 1966, 45, 4256–4266. [Google Scholar] [CrossRef]
  33. Čížek, J. On the Use of the Cluster Expansion and the Technique of Diagrams in Calculations of Correlation Effects in Atoms and Molecules. In Advances in Chemical Physics; LeFebvre, R., Moser, C., Eds.; Wiley: Hoboken, NJ, USA, 1969; Volume 14, pp. 35–89. [Google Scholar] [CrossRef]
  34. Lee, T.J.; Taylor, P.R. A diagnostic for determining the quality of single-reference electron correlation methods. Int. J. Quantum Chem. 1989, 36, 199–207. [Google Scholar] [CrossRef] [Green Version]
  35. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef] [Green Version]
  36. Schwartz, C. Importance of Angular Correlations between Atomic Electrons. Phys. Rev. 1962, 126, 1015–1019. [Google Scholar] [CrossRef]
  37. Boys, S.F.; Bernardi, F. The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol. Phys. 1970, 19, 553–566. [Google Scholar] [CrossRef]
  38. Abramowitz, M.; Irene, S. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, 9th ed.; Dover: New York, NY, USA, 1964; p. 470. [Google Scholar]
  39. Prosmiti, R.; Cunha, C.; Villarreal, P.; Delgado-Barrio, G. The van der Waals potential energy surfaces and structures of HeICl and NeICl clusters. J. Chem. Phys. 2002, 117, 7017–7023. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Minimum energy values obtained from RCCSD(T) level of theory using the AVQZ and AV5Z basis sets and CBS extrapolation as a function of θ (deg).
Figure 1. Minimum energy values obtained from RCCSD(T) level of theory using the AVQZ and AV5Z basis sets and CBS extrapolation as a function of θ (deg).
Molecules 28 05512 g001
Figure 2. Comparison between the 2D-PES RKHS analytical fitting and the RCCSD(T)/CBS(Q5) ab-initio surfaces of L i 2 + K r .
Figure 2. Comparison between the 2D-PES RKHS analytical fitting and the RCCSD(T)/CBS(Q5) ab-initio surfaces of L i 2 + K r .
Molecules 28 05512 g002
Figure 3. Minimum energy values obtained from the RKHS analytical interpolation with the corresponding RCCSD(T)/CBS[Q5] interaction energies as a function of angular orientations θ.
Figure 3. Minimum energy values obtained from the RKHS analytical interpolation with the corresponding RCCSD(T)/CBS[Q5] interaction energies as a function of angular orientations θ.
Molecules 28 05512 g003
Figure 4. Two-dimensional contour plot of the analytical potential using the RKHS analytical interpolation in the ( θ ,   R ) plane.
Figure 4. Two-dimensional contour plot of the analytical potential using the RKHS analytical interpolation in the ( θ ,   R ) plane.
Molecules 28 05512 g004
Figure 5. Radial (left) and angular (right) distributions of both even (e) and odd (o) states of L i 2 + K r complex for n = 0–5.
Figure 5. Radial (left) and angular (right) distributions of both even (e) and odd (o) states of L i 2 + K r complex for n = 0–5.
Molecules 28 05512 g005
Figure 6. Computed radial (left) and angular (right) distributions for n = 114 to n = 119 bound states of L i 2 + K r for both even (right line) and odd (dashed line) states.
Figure 6. Computed radial (left) and angular (right) distributions for n = 114 to n = 119 bound states of L i 2 + K r for both even (right line) and odd (dashed line) states.
Molecules 28 05512 g006
Figure 7. Computed radial (left) and angular (right) distributions for n = 240, 250, 260, and 272 bound states of L i 2 + K r for both even (right line) and odd (dashed line) states.
Figure 7. Computed radial (left) and angular (right) distributions for n = 240, 250, 260, and 272 bound states of L i 2 + K r for both even (right line) and odd (dashed line) states.
Molecules 28 05512 g007
Figure 8. Jacobi coordinates ( r ,   R ,   θ ) of L i 2 + K r complex.
Figure 8. Jacobi coordinates ( r ,   R ,   θ ) of L i 2 + K r complex.
Molecules 28 05512 g008
Table 1. Equilibrium distance ( R e ) and the well depth ( D e ) for the L i 2 + K r using the RCCSD(T) method and different basis sets. Re are in Å and De in cm−1.
Table 1. Equilibrium distance ( R e ) and the well depth ( D e ) for the L i 2 + K r using the RCCSD(T) method and different basis sets. Re are in Å and De in cm−1.
Basis SetAVQZAV5ZCBS[Q5]
θ   ( ° ) R e D e R e D e R e D e
04.2217304.1717374.171745
104.1417154.1317224.131730
204.0116984.0117064.011715
303.8815643.8715723.871580
403.6913913.6814003.681408
503.5211223.5111303.511139
603.457753.457833.44791
703.614793.604853.59492
803.853403.833463.82353
903.933043.923103.90317
Table 2. Trends of equilibrium distances ( R e ) and the well depths ( D e ) for L i 2 + ( X 2 g + ) alkali dimer in interaction with H e ,   N e and K r rare gas atoms. R e are in Å and D e in cm−1.
Table 2. Trends of equilibrium distances ( R e ) and the well depths ( D e ) for L i 2 + ( X 2 g + ) alkali dimer in interaction with H e ,   N e and K r rare gas atoms. R e are in Å and D e in cm−1.
Complex R e D e Method/References
L i 2 + ( X 2 g + ) H e 3.53342RCCSD(T)/CBS[Q5] [24]
L i 2 + ( X 2 g + ) N e -700Pseudopotential [15]
L i 2 + ( X 2 g + ) K r 4.011745RCCSD(T)/CBS[Q5]
Table 3. The equilibrium distance ( R e ) and the well depth ( D e ) of L i 2 + K r complex obtained with the RCCSD(T)/CBS[Q5] and the RKHS method for θ ranging from 0° to 90°.
Table 3. The equilibrium distance ( R e ) and the well depth ( D e ) of L i 2 + K r complex obtained with the RCCSD(T)/CBS[Q5] and the RKHS method for θ ranging from 0° to 90°.
RCCSD(T)/CBS[Q5]RKHS E   ( % )
θ   ( ° ) R e   ( Å ) D e   ( cm 1 ) R e   ( Å ) D e   ( cm 1 )
4.1717454.1717460.057
10°4.1317304.1317300.000
20°4.0117154.0317150.000
30°3.8715803.8715810.063
40°3.6814083.6814070.071
50°3.5111393.5011410.175
60°3.447913.447900.126
70°3.594923.594920.000
80°3.823534.823530.000
90°3.903174.903170.000
Table 4. Vibrational energies (in cm−1). of the lowest and some selected bound states ( J = 0) that lie below, around, and above the T-shaped 90° barrier well of the L i 2 + K r complex.
Table 4. Vibrational energies (in cm−1). of the lowest and some selected bound states ( J = 0) that lie below, around, and above the T-shaped 90° barrier well of the L i 2 + K r complex.
n j = even/odd
0−1633.440613/−1633.440613
1−1562.620803/−1562.620803
2−1492.719188/−1492.719188
3−1481.617393/−1481.617393
4−1422.606448/−1422.606448
5−1405.444466/−1405.444466
6−1352.065571/−1352.065571
7−1339.178576/−1339.178576
8−1331.488421/−1331.488421
9−1281.201931/−1281.201931
10−1263.686846/−1263.686846
114−335.655266/−335.655351
115−335.190171/−335.190064
116−326.764948/−326.764891
117−322.031270/−322.031221
118−318.591568/−318.591555
119−316.216834/−316.216848
240−25.004181/−18.144948
250−14.747134/−8.049236
260−6.805258/−0.477206
272−0.298997/-
Table 5. T 1 and   D 1 diagnostics for L i 2 + K r complex around the equilibrium positions.
Table 5. T 1 and   D 1 diagnostics for L i 2 + K r complex around the equilibrium positions.
θ   ( ° ) T 1 D 1
0.003378670.00727516
10°0.003461170.00761505
20°0.009730380.0087485
30°0.004071420.00994863
40°0.004556820.01166247
50°0.005102740.01350919
60°0.005469150.01473819
70°0.005085290.01356188
80°0.004300780.01101894
90°0.003995470.00998581
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saidi, S.; Mabrouk, N.; Dhiflaoui, J.; Berriche, H. Structural, Spectroscopic, and Dynamic Properties of Li2+(X2g+) in Interaction with Krypton Atom. Molecules 2023, 28, 5512. https://doi.org/10.3390/molecules28145512

AMA Style

Saidi S, Mabrouk N, Dhiflaoui J, Berriche H. Structural, Spectroscopic, and Dynamic Properties of Li2+(X2g+) in Interaction with Krypton Atom. Molecules. 2023; 28(14):5512. https://doi.org/10.3390/molecules28145512

Chicago/Turabian Style

Saidi, Samah, Nesrine Mabrouk, Jamila Dhiflaoui, and Hamid Berriche. 2023. "Structural, Spectroscopic, and Dynamic Properties of Li2+(X2g+) in Interaction with Krypton Atom" Molecules 28, no. 14: 5512. https://doi.org/10.3390/molecules28145512

APA Style

Saidi, S., Mabrouk, N., Dhiflaoui, J., & Berriche, H. (2023). Structural, Spectroscopic, and Dynamic Properties of Li2+(X2g+) in Interaction with Krypton Atom. Molecules, 28(14), 5512. https://doi.org/10.3390/molecules28145512

Article Metrics

Back to TopTop