Ground‐state properties of alkali dimers XY (X, Y=Li to Cs)

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.

Theoretical study of the Coriolis effect in LiNa, LiK, and LiRb molecules

The non-adiabatic electronic matrix elements, LΠΣ(R), that arise from the spin-conserving electron-rotational interactions between all mΣ+ and mΠ states, where multiplicity m = 1, 3, converging to the lowest three dissociation limits of Li-containing alkali diatomics, LiM (M = Na, K, Rb), were calculated ab initio up to large internuclear distances, R. The required electronic wavefunctions were obtained within the framework of the multi-reference configuration interaction treatment of the two-valence-electron problem constructed using small-core scalar-relativistic effective core potentials and l-independent core-polarization potentials. A least squares analysis of the ab initio functions at large internuclear distances in conjunction with long-range perturbation theory (LRPT) revealed three different asymptotic behaviors of the LΠΣ(R → +∞)-functions: const. + β[n]/Rn, characterized by n = -1, 3 and 6. The asymptotic coefficients β[n], extracted from the point-wise ab initio data, were found to be in agreement with their LRPT counterparts, which were evaluated analytically using the relevant atomic parameters. The mass dependence of the LΠΣ matrix elements was investigated analytically and numerically. To confirm the reliability of the LΠΣ(R)-functions and interatomic potentials at small and intermediate distances, the empirical q-factors available for the D1Π-states of all LiM molecules studied were compared with their theoretical counterparts derived from the present ab initio data.

Measurement of the permanent electric dipole moment of ultracold ground state 85Rb133Cs molecules by microwave coherent spectroscopy

We demonstrate measurement of the permanent electric dipole moment (EDM) of ⁸⁵Rb¹³³Cs molecules in the absolute vibrational ground state by microwave (MW) coherent spectroscopy. The rotational states of the considered molecules, which are formed from short-range photoassociation of mixed cold atoms, are nondegenerated under external electric field. To measure the EDM based on electric-field-induced shifts of the sublevels of X¹Σ⁺(v = 0, J = 1) rotational state, we utilized a MW coherent spectroscopy, which has a higher resolution than depletion spectroscopy and one-photon MW spectroscopy and can also eliminate the influence from Stark shift of the excited state existing in both spectroscopies above. In order to acquire accurate electric intensity, electromagnetic induced transparency spectroscopy of ⁸⁵Rb Rydberg atoms is used to implement the calibration. The permanent EDM of ⁸⁵Rb¹³³Cs molecules is finally determined to be 1.266(15) D, which agrees with the theoretical calculations and is comparable with the value of its isotopic molecule.

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.

Ab initio studies on the spin-forbidden cooling transitions of the LiRb molecule

The spin-forbidden cooling of the LiRb molecule is investigated based on ab initio quantum chemistry calculations. The multireference configuration interaction method is used to generate the potential energy curves (PECs) of the ground state X(1)Σ(+) and the low-lying excited states a(3)Σ(+), B(1)Π, and b(3)Π. The spin-orbit coupling effects for the PECs and the transition dipole moments (TDMs) between the X(1)Σ(+), b(3)Π and a(3)Σ(+) states are also calculated. The analytical functions for the PECs are deduced. The rovibrational energy levels, the spectroscopic parameters and the Franck-Condon factors (FCF) are determined by solving the Schrödinger equation of nuclear movement with the obtained analytical functions. The b(3)Π0 ↔ X(1)Σ(+) and b(3)Π1 ↔ X(1)Σ(+) transitions have highly diagonal distributed FCFs and non-zero TDMs, demonstrating that the LiRb molecule could be a very promising candidate for laser cooling. Therefore, a three-cycle laser cooling scheme for the molecule has been proposed based on these two spin-forbidden transitions. Using the radiative lifetime and linewidth calculated from the obtained TDM functions, we present further analysis of the cooling of LiRb and the corresponding KRb molecule. The transition b(3)Π0 ↔ X(1)Σ(+) is found to be a practical transition to cool the LiRb molecule, and a sub-microkelvin cool temperature could be reached for the KRb molecule using a similar laser cooling scheme.

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.

Creation of an Ultracold Gas of Ground-State Dipolar ^{23}Na^{87}Rb Molecules

We report the successful production of an ultracold sample of absolute ground-state ^{23}Na^{87}Rb molecules. Starting from weakly bound Feshbach molecules formed via magnetoassociation, the lowest rovibrational and hyperfine level of the electronic ground state is populated following a high-efficiency and high-resolution two-photon Raman process. The high-purity absolute ground-state samples have up to 8000 molecules and densities of over 10^{11}  cm^{-3}. By measuring the Stark shifts induced by external electric fields, we determined the permanent electric dipole moment of the absolute ground-state ^{23}Na^{87}Rb and demonstrated the capability of inducing an effective dipole moment over 1 D. Bimolecular reaction between ground-state ^{23}Na^{87}Rb molecules is endothermic, but we still observed a rather fast decay of the molecular sample. Our results pave the way toward investigation of ultracold molecular collisions in a fully controlled manner and possibly to quantum gases of ultracold bosonic molecules with strong dipolar interactions.

Analytic functions for potential energy curves, dipole moments, and transition dipole moments of LiRb molecule

The analytic potential energy functions (APEFs) of the X(1)Σ(+), 2(1)Σ(+), a(3)Σ(+), and 2(3)Σ(+) states of the LiRb molecule are obtained using Morse long-range potential energy function with damping function and nonlinear least-squares method. These calculations were based on the potential energy curves (PECs) calculated using the multi-reference configuration interaction (MRCI) method. The reliability of the APEFs is confirmed using the curves of their first and second derivatives. By using the obtained APEFs, the rotational and vibrational energy levels of the states are determined by solving the Schrödinger equation of nuclear movement. The spectroscopic parameters, which are deduced using Dunham expansion, and the obtained rotational and vibrational levels are compared with the reported theoretical and experimental values. The correlation effect of the electrons of the inner shell remarkably improves the results compared with the experimental spectroscopic parameters. For the first time, the APEFs for the dipole moments and transition dipole moments of the states have been determined based on the curves obtained from the MRCI calculations.

Relativistic state-specific multireference perturbation theory incorporating improved virtual orbitals: Application to the ground state single-bond dissociation

Using four-component (4c) relativistic spinors, we present a computationally economical relativistic ab initio method for molecular systems employing our recently proposed second-order state-specific multireference perturbation theory (SSMRPT) incorporating the improved virtual orbital-complete active space configuration interaction (IVO-CASCI) reference wavefunction. The resulting method, 4c-IVO-SSMRPT [calculate one state at a time] is tested in pilot calculations on the homonuclear dimers including Li(2), Na(2), K(2), Rb(2), F(2), Cl(2), and Br(2) through the computations of the ground state potential energy curves (PECs). As SSMRPT curbs intruder effects, 4c-IVO-SSMRPT is numerically stable. To our knowledge, the SSMRPT in the 4c relativistic framework has not been explored in the past. Selective spectroscopic constants that are closely related to the correct shape and accuracy of the energy surfaces have been extracted from the computed PECs. For the halogen molecules, a relativistic destabilization of the bond has been found. Relativistic and electron correlation effects need to be incorporated to get reliable estimates. Our results are in good accordance with reference theoretical and experimental data which manifests the computational accuracy and efficiency of the new 4c-IVO-SSMRPT method. The method opens for an improved description of MR systems containing heavy elements. The inexpensiveness of IVO-CASCI makes 4c-IVO-SSMRPT method promising for studies on large systems of heavy elements.

Theoretical study of the CsNa molecule: adiabatic and diabatic potential energy and dipole moment

The adiabatic and diabatic potential energy curves of the low-lying electronic states of the NaCs molecule dissociating into Na (3s, 3p) + Cs (6s, 6p, 5d, 7s, 7p, 6d, 8s, 4f) have been investigated. The molecular calculations are performed using an ab initio approach based on nonempirical pseudopotential, parametrized l-dependent polarization potentials and full configuration interaction calculations through the CIPCI quantum chemistry package. The derived spectroscopic constants (Re, De, Te, ωe, ωexe, and Be) of the ground state and lower excited states are compared with the available theoretical and experimental works. Moreover, accurate permanent and transition dipole moment have been determined as a function of the internuclear distance. The adiabatic permanent dipole moment for the first nine (1)Σ(+) electronic states have shown both ionic characters associated with electron transfer related to Cs(+)Na(-) and Cs(-)Na(+) arrangements. By a simple rotation, the diabatic permanent dipole moment is determined and has revealed a linear behavior, particularly at intermediate and large distances. Many peaks around the avoided crossing locations have been observed for the transition dipole moment between neighbor electronic states.

Comment on "Going beyond the frozen core approximation: development of coordinate-dependent pseudopotentials and application to Na2(+)" [J. Chem. Phys. 138, 054110 (2013)]

The new coordinate-dependent pseudopotential for Na2(+) by Kahros and Schwartz [J. Chem. Phys. 138, 054110 (2013)] is assessed and compared to the pseudopotential approach by Fuentealba et al. [Chem. Phys. Lett. 89, 418 (1982)] which incorporates the coordinate-dependent core-polarization potential by Müller and Meyer [J. Chem. Phys. 80, 3311 (1984)]. In contrast to the latter approach, the one by Kahros and Schwartz does not reproduce the accurately known experimental data and∕or high level theoretical results for Na2(+). The treatment of core polarization by Kahros and Schwartz neglects the dynamic polarization of atomic cores which is much more important for Na2(+) than the static one. On the other hand, the Kahros and Schwartz method heavily overestimates frozen-core corrections at the Hartree-Fock level by compounding them with artifacts of a superposition of non-norm-conserving pseudopotentials.

Relations between harmonic frequencies of diatomic molecules

The relations between the harmonic frequencies of different molecules are revealed with the aid of the spring constants of atoms in molecules. Using the atomic spring constants in the related molecules, the force constants for a new molecule can be estimated. The simplest scheme to obtain the force constant of a given molecule is similar to a simple chemical reaction formula, such as A(2) + B(2) → AB, and the corresponding relation between the molecular force constants is k(AB)(-1) = (2k(A(2)))(-1) + (2k(B(2)))(-1). For a given molecule, one can design numerous schemes to obtain its force constant from the atomic spring constants in other molecules. A high degree of periodical regularity appears in the application of different kinds of schemes to the ground states. The reliable schemes for the ground electronic states can be adopted for the excited states. Over two hundred molecules with experimental data available for comparison have been tested. The discrepancies between the calculated and the experimental harmonic frequencies can reach 1% and better; the results show that the present approach is simple in theory and handy to use. The harmonic frequencies for dozens of hetero-nuclear molecules of the transition-metal elements are also predicted.

Collisional control of ground state polar molecules and universal dipolar scattering

We explore the impact of the short-range interaction on the scattering of ground state polar molecules and study the transition from a weak to strong dipolar scattering over an experimentally reasonable range of energies and electric field values. In the strong dipolar limit, the scattering scales with respect to a dimensionless quantity defined by mass, induced dipole moment, and collision energy. The scaling has implications for all quantum mechanical dipolar scattering. Furthermore the universal scattering regime will readily be achieved with polar molecules at ultracold temperatures.

Theoretical investigation of RbCs via two-component spin-orbit pseudopotentials: Spectroscopic constants and permanent dipole moment functions

Potential energy curves for the 28 lowest LambdaSigma states and 49 Omega states of RbCs are obtained from large-scale multireference configuration interaction calculations using both spin-averaged and two-component spin-orbit energy-consistent effective core potentials. Spectroscopic properties of all states are compared across available data in literature to date. Variations of the permanent dipole moments on the internuclear separation (R) for the (1)Sigma(+), (3)Sigma(+), (1)Pi, and (3)Pi states are evaluated over a wide range of R. The most important effects of the spin-orbit interaction on the dipole moment distribution are discussed.

Theoretical calculation of the excited states of the KCs molecule including the spin-orbit interaction

For the molecule KCs the potential energy has been calculated for the 72 lowest molecular states Omega. Using an ab initio method the calculation is based on nonempirical pseudopotentials within the range of 5.0a0-34.0a0 of the internuclear distance R. Gaussian basis sets have been used for both atoms and spin-orbit effects have been taken into account through a semiempirical spin-orbit pseudopotential added to the electrostatic Hamiltonian. The spectroscopic constants for 60 states have been calculated by fitting the calculated energy values to a polynomial in terms of the internuclear distance R. The components of the spin-orbit splitting for (1,2,5,6) 3Pi and (1) 3delta have been identified. The comparison of the present results with those available in the literature shows a very good agreement, while the other results, to the best of our knowledge, are given here for the first time.

Ab initio relativistic calculation of the RbCs molecule

We apply the relativistic configuration-interaction valence-bond method to calculate various characteristics of the alkali-metal RbCs dimer. These include the electronic potentials and transition dipole moments between the ground and first excited states and permanent dipole moments of the X 1sigma+ and a 3sigma+ states of the ground configuration. In addition, we estimate the lifetime of the rovibrational levels of the X state due to blackbody radiation. These data can help experimentalists to optimize photoassociative formation of ultracold RbCs molecules and their longevity in a trap or in an optical lattice. Extended basis sets, constructed from Dirac-Fock and Sturm's orbitals, have been used to ensure convergence of our calculations. We compare our data with other theoretical and experimental results when they were available.

Optical production of ultracold polar molecules

We demonstrate the production of ultracold polar RbCs molecules in their vibronic ground state, via photoassociation of laser-cooled atoms followed by a laser-stimulated state transfer process. The resulting sample of X1Sigma+ (nu = 0) molecules has a translational temperature of approximately 100 microK and a narrow distribution of rotational states. With the method described here it should be possible to produce samples even colder in all degrees of freedom, as well as other bialkali species.

Calculation of accurate permanent dipole moments of the lowest Σ+1,3 states of heteronuclear alkali dimers using extended basis sets

Obtaining ultracold samples of dipolar molecules is a current challenge which requires an accurate knowledge of their electronic properties to guide the ongoing experiments. In this paper, we systematically investigate the ground state and the lowest triplet state of mixed alkali dimers (involving Li, Na, K, Rb, Cs) using a standard quantum chemistry approach based on pseudopotentials for atomic core representation, Gaussian basis sets, and effective terms for core polarization effects. We emphasize on the convergence of the results for permanent dipole moments regarding the size of the Gaussian basis set, and we discuss their predicted accuracy by comparing to other theoretical calculations or available experimental values. We also revisit the difficulty to compare computed potential curves among published papers, due to the differences in the modelization of core-core interaction.

Ground-state properties and static dipole polarizabilities of the alkali dimers from K2n to Fr2n(n=0,+1) from scalar relativistic pseudopotential coupled cluster and density functional studies

The newly adjusted energy-consistent nine-valence-electron pseudopotentials for K to Fr are used to calculate spectroscopic properties for the neutral and positively charged alkali dimers using coupled cluster and density functional theory. For the neutral dimers the static dipole polarizability was calculated. The coupled cluster results are all in excellent agreement with experimental values. The density functionals used can give quite different spectroscopic properties especially for the dipole polarizability, with the Perdew-Wang PW91 functional performing best.

Quantum computation with trapped polar molecules

We propose a novel physical realization of a quantum computer. The qubits are electric dipole moments of ultracold diatomic molecules, oriented along or against an external electric field. Individual molecules are held in a 1D trap array, with an electric field gradient allowing spectroscopic addressing of each site. Bits are coupled via the electric dipole-dipole interaction. Using technologies similar to those already demonstrated, this design can plausibly lead to a quantum computer with greater, approximately > or = 10(4) qubits, which can perform approximately 10(5) CNOT gates in the anticipated decoherence time of approximately 5 s.

Theoretical Study of the Electronic States of the Rb(2) Molecule

We have calculated the electronic states of Rb(2) by multireference configuration interactions using the averaged relativistic effective small-core potential and the core-polarization potential. The potential energy curves for a large number of states dissociating into from 5s+5s up to 7s+5s asymptotic limits are calculated and the spectroscopic constants are reported. The spin-orbit effects for the states dissociating into 5p+5s and 4d+5s are calculated using the effective spin-orbit potential. The results are compared with available experimental data and other theoretical works. Copyright 2001 Academic Press.

The Ground State X1Sigma+ of NaK Revisited

Laser-induced fluorescence (LIF) optical spectra of the D1Pi-X1Sigma+ system of the NaK molecule have been reinvestigated. In this new analysis, quantum numbers have been assigned to 3000 selected fluorescence lines for a wide range of vibrational and rotational levels (J </= 149, v' </= 23, v" </= 47). A new set of constants of the X1Sigma state (Dunham-type coefficients) and term values of the D1Pi state have been determined by a weighted least-squares fit of the LIF optical lines of the D-X system, some lines of the C-X system, and microwave spectral data reported recently. Constants for X1Sigma (in cm-1) are as follows: omegae = 123.993, omegaexe = 0.3045, Be = 0.095229, and alphae = 4.48 x 10(-4). The equilibrium internuclear distance is re = 3.49903 Å, very little changed from the microwave value. Copyright 1998 Academic Press.