First Principle Calculations of the Potential Energy Curves for Electronic States of the Lithium Dimer

An Intruder-Free Fock Space Coupled-Cluster Study of the Potential Energy Curves of LiMg+ within the (2,0) Sector

The potential energy curves (PECs) and spectroscopic constants of the ground and excited states of a LiMg+ molecular cation were investigated. We obtained accurate results for the fifteen lowest-lying states of the LiMg+ cation using the Intermediate Hamiltonian Fock Space Multireference Coupled Cluster (IH-FS-CC) method applied to the (2,0) sector. Relativistic corrections were accounted for using the third-order Douglas-Kroll method. In each instance, smooth PECs were successfully computed across the entire range of interatomic distances from equilibrium to the dissociation limit. The results are in good accordance with previous studies of this molecular cation. Notably, this study marks the first application of IH-FS-CC in investigating a mixed alkali and alkaline earth molecular cation, proving its usability in determining accurate PECs of such diatomics and their spectroscopic constants.

Benchmark Study of the Electronic States of the LiRb Molecule: Ab Initio Calculations with the Fock Space Coupled Cluster Approach

Accurate potential energy curves (PECs) are determined for the twenty-two electronic states of LiRb. In contrast to previous studies, the applied approach relies on the first principle calculations involving correlation among all electrons. The current methodology is founded on the multireference coupled cluster (CC) scheme constructed within the Fock space (FS) formalism, specifically for the (2,0) sector. The FS methodology is established within the framework of the intermediate Hamiltonian formalism and offers an intruder-free, efficient computational scheme. This method has a distinctive feature that, when applied to the doubly ionized system, provides the characteristics of the neutral case. This proves especially beneficial when investigating PECs in situations where a closed-shell molecule dissociates into open-shell fragments, yet its double positive ion forms closed-shell species. In every instance, we successfully computed continuous PECs spanning the entire range of interatomic distances, from the equilibrium to the dissociation limit. Moreover, the spectroscopic characteristic of various electronic states is presented, including relativistic effects. Relativistic corrections included at the third-order Douglas-Kroll level have a non-negligible effect on the accuracy of the determined spectroscopic constants.

On the nature of the two-positron bond: evidence for a novel bond type

The nature of the newly proposed two-positron bond in (PsH)2, which is composed of two protons, four electrons and two positrons, is considered in this contribution. The study is done at the multi-component-Hartree-Fock (MC-HF) and the Diffusion Monte Carlo (DMC) levels of theory by comparing ab initio data, analyzing the spatial structure of the DMC wavefunction, and applying the multi-component quantum theory of atoms in molecules and the two-component interacting quantum atoms energy partitioning schemes to the MC-HF wavefunction. The analysis demonstrates that (PsH)2 to a good approximation may be conceived of as two slightly perturbed PsH atoms, bonded through a two-positron bond. In contrast to the usual two-electron bonds, the positron exchange phenomenon is quite marginal in the considered two-positron bond. The dominant stabilizing mechanism of bonding is a novel type of classical electrostatic interaction between the positrons, which are mainly localized between nuclei, and the surrounding electrons. To emphasize its uniqueness, this mechanism of bonding is proposed to be called gluonic which has also been previously identified as the main driving mechanism behind formation of the one-positron bond in [H-,e+,H-]. We conclude that the studied two-positron bond should not be classified as a covalent bond and it must be seen as a brand-new type of bond, foreign to the electronic bonding modes discovered so far in the purely electronic systems.

An accurate potential model for the a3Σu+ state of the lithium dimer

An accurate Tang-Toennies (TT) model potential is introduced to describe the interatomic potential of the lithium dimer in the a³Σ_u⁺ state. With only one well-known parameter, the ionization energy, the new model potential compares favorably with the experimentally fitted Morse/Long-range (MLR) potential of Dattani and Le Roy [J. Mol. Spectrosc., 2011, 268, 199] and is in excellent agreement with the state-of-the-art ab initio potential of Lesiuk et al. [Phys. Rev. A, 2020, 102, 062806]. With the known dispersion coefficients and the ionization energy, the new potential requires only two experimental parameters, namely the depth of the potential well D_e and its location R_e. The new potential can be extended to the region of zero separation by the united atom limit.

Resonant Ionic, Covalent Bond, and Steric Characteristics Present in 1Σu+ States of Li2

The molecular bonding in the excited states of the alkali dimers involves the resonant ionic, covalent bond and steric interactions. We show here the case of the ¹Σ_u⁺ states of Li₂ by ab initio calculation. These interactions as functions of the internuclear distance lead to complex potential energy curves, providing an important application for high resolution laser spectroscopy. The spectroscopic constants for the 4 and 5 ¹Σ_u⁺ states are obtained for the first time.

Intermediate Hamiltonian Fock-space coupled-cluster theory for excitation energies, double ionization potentials, and double electron attachments with spin–orbit coupling

The intermediate Hamiltonian Fock-space coupled-cluster methods at the singles and doubles level (IHFSCCSD) for excitation energies in the (1p, 1h) sector, double ionization potentials in the (0p, 2h) sector, and double electron attachments in the (2p, 0h) sector of the Fock space are implemented based on the CCSD method with spin–orbit coupling (SOC) included in the post-Hartree–Fock treatment using a closed-shell reference in this work. The active space is chosen to contain those orbitals that have the largest contribution to principal ionized or electron-attached states obtained from the equation-of-motion coupled-cluster calculations. Both time-reversal symmetry and spatial symmetry are exploited in the implementation. Our results show that the accuracy of IHFSCCSD results is closely related to the active space, and the sufficiency of the active space can be assessed from the percentage of transitions within the active space. In addition, unreasonable results may be encountered when the ionized or electron-attached states with a somewhat larger contribution from double excitations are included to determine the active space and cluster operators in the (0p, 1h) or (1p, 0h) sector of the Fock space. A larger active space may be required to describe SO splitting reliably than that in the scalar-relativistic calculations in some cases. The IHFSCCSD method with SOC developed in this work can provide reliable results for heavy-element systems when a sufficient active space built upon the principal ionization potential/electron affinity states is adopted.

Signature of the neighbor's quantum nuclear dynamics in the electron transfer mediated decay spectra

We computed fully quantum nuclear dynamics, which accompanies electron transfer mediated decay (ETMD) in weakly bound polyatomic clusters. We considered two HeLi₂ clusters - with Li₂ being either in the singlet electronic ground state or in the triplet first excited state - in which ETMD takes place after ionization of He. The electron transfer from Li₂ to He⁺ leads to the emission of another electron from Li₂ into the continuum. Due to the weak binding of He to Li₂ in the initial states of both clusters, the involved nuclear wavepackets are very extended. This makes both the calculation of their evolution and the interpretation of the results difficult. We showed that despite the highly delocalized nature of the wavepackets the nuclear dynamics in the decaying state is imprinted on the ETMD electron spectra. The analysis of the latter helps understanding the effect which electronic structure and binding strength in the cluster produce on the quantum motion of the nuclei in the decaying state. The results produce a detailed picture of this important charge transfer process in polyatomic systems.

Electronic transitions in Rb2+ dimers solvated in helium

We have measured depletion spectra of the heteronuclear (85Rb87Rb+) dimer cation complexed with up to 10 He atoms. Two absorption bands are observed between 920 and 250 nm. The transition into the repulsive 12Σu+state of HeRb2+gives rise to a broad feature at 790 nm (12,650 cm−1); it exhibits a blueshift of 98 cm−1per added He atom. The transition into the bound 12Πustate of HeRb2+reveals vibrational structure with a band head at ≤ 15,522 cm−1, a harmonic constant of 26 cm−1, and a spin–orbit splitting of ≤ 183 cm−1. The band experiences an average redshift of − 38 cm−1per added He atom. Ab initio calculations rationalize the shape of the spectra and spectral shifts with respect to the number of helium atoms attached. For a higher number of solvating helium atoms, symmetric solvation on both ends of the Rb2+ion is predicted.

Electronic states of NaLi molecule: Benchmark results with Fock space coupled cluster approach

Accurate potential energy curves (PECs) are obtained for 20 lowest lying electronic states of the NaLi molecule. The computational scheme used here is based on the multireference coupled cluster theory formulated in the (2,0) sector of the Fock space. The latter sector provides the description of states obtained by attachment of two electrons to the reference system. This makes it possible to adopt the doubly ionized NaLi⁺² molecule as a Fermi vacuum. The latter has a very concrete advantage in calculations of the PECs since it dissociates into closed shell fragments (NaLi⁺² → Na⁺ + Li⁺); hence, the restricted Hartree-Fock method can be used within the whole range of interatomic distances. Computed PECs and spectroscopic constants stay very close to the experimental values (if the latter are available) with the accuracy exceeding the other theoretical approaches including those based on the effective core polarization potentials. Relativistic corrections included at the infinite-order two-component level have a non-negligible effect on the accuracy of computed excitation and dissociation energies with contributions up to 50 cm^-1.

Extension of the Fock-space coupled-cluster method with singles and doubles to the three-valence sector

The single-reference coupled-cluster method has proven very effective in the ab initio description of atomic and molecular systems, but its successful application is limited to states dominated by a single Slater determinant, which is used as the reference. In cases where several determinants are important in the wave function expansion, i.e., we have to deal with nondynamic correlation effects, a multi-reference version of the coupled-cluster method is required. The multi-reference coupled-cluster approaches are based on the effective Hamiltonian formulation providing a two-step procedure, in which dynamic correlation effects can be efficiently evaluated by the wave operator, while nondynamic correlation contributions are given by diagonalization of the effective Hamiltonian in the final step. There are two classical multi-reference coupled-cluster formulations. In this paper, the focus is on the so-called Fock-space coupled-cluster method in its basic version with one- and two-particle operators in the exponent. Computational schemes using this truncation of the cluster operator have been successfully applied in calculations in one- and two-valence sectors of the Fock space. In this paper, we show that the approach can be easily extended and effectively employed in the three-valence sector calculations.

A first principles study of the spin–orbit coupling effect in LiM (M = Na, K, Rb, Cs) molecules

Both fully relativistic and scalar-state based perturbation models provided the spin–orbit functions of the LiM (M = Na, K, Rb, Cs) molecules at almost experimental level of confidence.

Competition between charge migration and charge transfer induced by nuclear motion following core ionization: Model systems and application to Li2. J Chem Phys 2019;151:124108. [PMID: 31575213 DOI: 10.1063/1.5117246]

Attosecond and femtosecond spectroscopies present opportunities for the control of chemical reaction dynamics and products, as well as for quantum information processing; we address the somewhat unique situation of core-ionization spectroscopy which, for dimeric chromophores, leads to strong valence charge localization and hence tightly paired potential-energy surfaces of very similar shape. Application is made to the quantum dynamics of core-ionized Li₂ ⁺. This system is chosen as Li₂ is the simplest stable molecule facilitating both core ionization and valence ionization. First, the quantum dynamics of some model surfaces are considered, with the surprising result that subtle differences in shape between core-ionization paired surfaces can lead to dramatic differences in the interplay between electronic charge migration and charge transfer induced by nuclear motion. Then, equation-of-motion coupled-cluster calculations are applied to determine potential-energy surfaces for 8 core-excited state pairs, calculations believed to be the first of their type for other than the lowest-energy core-ionized molecular pair. While known results for the lowest-energy pair suggest that Li₂ ⁺ is unsuitable for studying charge migration, higher-energy pairs are predicted to yield results showing competition between charge migration and charge transfer. Central is a focus on the application of Hush's 1975 theory for core-ionized X-ray photoelectron spectroscopy to understand the shapes of the potential-energy surfaces and hence predict key features of charge migration.

Benchmark of correlation matrix renormalization method in molecule calculations

We report benchmark calculations of the correlation matrix renormalization (CMR) approach for 23 molecules in the well-established G2 molecule set. This subset represents molecules with spin-singlet ground state in a variety of chemical bonding and coordination environments. The QUAsi-atomic minimal basis-set orbitals (QUAMBOs) are used as local orbitals in both CMR and full configuration interaction (FCI) calculations for comparison. The results obtained from the calculations are also compared with available experimental data. It is shown that the CMR method produces binding and dissociation energy curves in good agreement with the QUAMBO-FCI calculations as well as experimental results. The CMR benchmark calculations yield a standard deviation of 0.09 Å for the equilibrium bond length and 0.018 Hartree/atom for the formation energy, with a gain of great computational efficiency which scales like Hartree-Fock method.

Electron transfer mediated decay in HeLi2 cluster: Potential energy surfaces and decay widths

Electron transfer mediated decay (ETMD) is a process responsible for double ionization of dopants in He droplets. It is initiated by producing He⁺ in the droplet, which is neutralized by ETMD, and has been shown to strongly enhance the dopant's double ionization cross section. The efficiency of ETMD, the spectra of emitted secondary electrons, and the character of the ionic products depend on the nuclear dynamics during the decay. To date, there has been no theoretical investigation of multimode dynamics which accompanies ETMD, which could help to understand such dynamics in a He droplet. In this article, we consider the He-Li₂ cluster where an ab initio examination of multimode dynamics during the electronic decay is feasible. Moreover, this cluster can serve as a minimal model for Li₂ adsorbed on the droplet's surface-a system where ETMD can be observed experimentally. In He droplets, Li₂ can be formed in both the ground X¹Σ_g ⁺ and the first excited a³Σ_u ⁺ states. In this article, we present ab initio potential energy surfaces of the electronic states of the He-Li₂ cluster involved in ETMD, as well as the respective decay widths. We show that the structure of these surfaces and expected nuclear dynamics strongly depend on the electronic state of Li₂. Thus, the overall decay rate and the appearance of the observable electron spectra will be dictated by the electronic structure of the dopant.

An accurate potential model for the a3Σu+ state of the alkali dimers Na2, K2, Rb2, and Cs2

A modified semi-empirical Tang-Toennies potential model is used to describe the a³Σ_u⁺ potentials of the alkali dimers. These potentials are currently of interest in connection with the laser manipulation of the ultracold alkali gases. The fully analytical model is based on three experimental parameters, the well depth D_e, well location R_e, and the harmonic vibrational frequency ω_e of which the latter is only slightly optimized within the range of the literature values. Comparison with the latest spectroscopic data shows good agreement for Na₂, K₂, Rb₂, and Cs₂, comparable to that found with published potential models with up to 55 parameters. The differences between the reduced potential of Li₂ and the conformal reduced potentials of the heavier dimers are analyzed together with why the model describes Li₂ less accurately. The new model potential provides a test of the principle of corresponding states and an excellent first order approximation for further optimization to improve the fits to the spectroscopic data and describe the scattering lengths and Feshbach resonances at ultra-low temperatures.

Note: Contamination-free loading of lithium metal into a nozzle source

This note describes a system for transferring a load of high purity lithium metal into a molecular or cluster beam source. A hot loading vessel is thoroughly baked out while empty and overpressured with argon. A clean Li rod is then dropped in through a long narrow tube. The thoroughly degassed interior of the vessel and the rapid melting of the inserted rod facilitate contamination-free transfer of the highly reactive liquid metal into the source oven.

Potential Energy Curves in the CASSCF/CASPT2 and FS-MR-CC Methods: The Role of Relativistic Effects

Ab initio CASSCF/CASPT2 calculations for the electronic ground and for a wide range of excited states of Li2 and Na2 dimers are presented. The computed spectroscopic parameters agree very well with the experimental data. This indicates that the old CASSCF/CASPT2 method can be as successfully applied to study excited states of molecules as recently developed the multireference Fock-space coulped-cluster method. The role of relativistic effects in the correct description of the potential energy curves has been investigated using as an example the SiAu molecule. The accuracy of the new infinite-order two-component relativistic method has been studied and its advantage over the Douglas-Kroll-Hess method demonstrated.

Spin Signature of Nonlocal Correlation Binding in Metal-Organic Frameworks

We develop a proper nonempirical spin-density formalism for the van der Waals density functional (vdW-DF) method. We show that this generalization, termed svdW-DF, is firmly rooted in the single-particle nature of exchange and we test it on a range of spin systems. We investigate in detail the role of spin in the nonlocal correlation driven adsorption of H_{2} and CO_{2} in the linear magnets Mn-MOF74, Fe-MOF74, Co-MOF74, and Ni-MOF74. In all cases, we find that spin plays a significant role during the adsorption process despite the general weakness of the molecular-magnetic responses. The case of CO_{2} adsorption in Ni-MOF74 is particularly interesting, as the inclusion of spin effects results in an increased attraction, opposite to what the diamagnetic nature of CO_{2} would suggest. We explain this counterintuitive result, tracking the behavior to a coincidental hybridization of the O p states with the Ni d states in the down-spin channel. More generally, by providing insight on nonlocal correlation in concert with spin effects, our nonempirical svdW-DF method opens the door for a deeper understanding of weak nonlocal magnetic interactions.