Research ArticlePHYSICS

Quantum computation solves a half-century-old enigma: Elusive vibrational states of magnesium dimer found

See allHide authors and affiliations

Science Advances  03 Apr 2020:
Vol. 6, no. 14, eaay4058
DOI: 10.1126/sciadv.aay4058

Abstract

The high-lying vibrational states of the magnesium dimer (Mg2), which has been recognized as an important system in studies of ultracold and collisional phenomena, have eluded experimental characterization for half a century. Until now, only the first 14 vibrational states of Mg2 have been experimentally resolved, although it has been suggested that the ground-state potential may support five additional levels. Here, we present highly accurate ab initio potential energy curves based on state-of-the-art coupled-cluster and full configuration interaction computations for the ground and excited electronic states involved in the experimental investigations of Mg2. Our ground-state potential unambiguously confirms the existence of 19 vibrational levels, with ~1 cm−1 root mean square deviation between the calculated rovibrational term values and the available experimental and experimentally derived data. Our computations reproduce the latest laser-induced fluorescence spectrum and provide guidance for the experimental detection of the previously unresolved vibrational levels.

INTRODUCTION

The weakly bound alkaline-earth dimers (AE2) have emerged as probes of fundamental physics relevant to ultracold collisions (1), doped helium nanodroplets (2), coherent control of binary reactions (3), and even fields rarely associated with molecular science, such as optical lattice clocks (4) and quantum gravity (5). The magnesium dimer is especially important, since it has several desirable characteristics that can be useful in the above applications, such as the absence of hyperfine structure in the most abundant 24Mg isotope that facilitates the analysis of binary collisions involving laser-cooled and trapped atoms, it helps us understand heavier AE2 diatomics, and, unlike its lighter Be2 analog, it is nontoxic (6). Unfortunately, the status of Mg2 as a prototype heavier AE2 species is complicated by the fact that its high-lying vibrational levels and, consequently, the long-range part of its ground-state potential energy curve (PEC) have eluded experimental characterization for half a century. In this regard, the magnesium dimer is even more challenging than its celebrated beryllium counterpart, whose elusive 12th vibrational level near the dissociation threshold (7, 8), which we also found in (9), was confirmed in 2014 (10) after reanalyzing the spectra obtained in stimulated emission pumping experiments (11).

Experimentally, probing vibrational manifold of the magnesium dimer in its ground, X1Σg+, electronic state has to involve excited electronic states, since Mg2, being a homonuclear diatomic, is infrared inactive. The first high-resolution photoabsorption spectra of Mg2, corresponding to a transition from the ground state to the electronically excited A1Σu+ state, were reported in 1970 by Balfour and Douglas (12). Their spectroscopic analysis resulted in 285 G(v″, J″) and 656 G(v′, J′) rovibrational term values of 24Mg2 involving 13 (v″ = 0 to 12) X1Σg+ and 24 (v′ = 1 to 24) A1Σu+ vibrational levels, respectively. Here, we are using the notation in which the vibrational, v, and rotational, J, quantum numbers in the ground electronic state are designated by a double prime, whereas those corresponding to the excited A1Σu+ state are marked with a prime. In their pioneering work, Balfour and Douglas constructed a Rydberg-Klein-Rees (RKR) (13–16) X1Σg+ PEC in the 3.25 to 7.16 Å range and located the last experimentally resolved v″ = 12 level about 25 cm−1 below the dissociation threshold, pointing to the existence of extra vibrational states with v″ > 12. It did not take long to detect one of such states. In 1973, Li and Stwalley (17) identified X1Σg+→A1Σu+ transitions involving the v″ = 13 level in the spectra reported in (12). They accomplished this by extending the original RKR PEC of Balfour and Douglas to the asymptotic region beyond 7.16 Å using theoretical values of C6 and C8 van der Waals coefficients (18, 19). The resulting PEC supported 19 vibrational levels, i.e., five levels more than what was observed experimentally (17). Four decades later, in an effort to characterize states with v″ > 13, Knöckel et al. (20, 21) examined the A1Σu+→X1Σg+ transition using laser-induced fluorescence (LIF), repeating and refining the earlier LIF experiment by Scheingraber and Vidal (22). They improved and expanded the original 24Mg2 dataset of Balfour and Douglas by reporting a total of 333 G(v″, J″) and 1,351 G(v′, J′) rovibrational term values involving v″ = 0 to 13 and v′ = 1 to 46, respectively, and constructed a few experimentally derived analytical forms of the X1Σg+ PEC, extrapolated to the asymptotic region using the theoretical C6 (23), C8 (24), and C10 (24) coefficients, which support the discrete spectral data in the 3.27 to 8.33 Å range (20). Although these refined PECs supported 19 24Mg2 vibrational levels, reinforcing the initial prediction of Li and Stwalley (17), Knöckel et al. (20) were unable to identify A1Σu+(v′,J′)→X1Σg+(v″,J″) transitions involving the elusive high-lying vibrational levels with v″ > 13 in their LIF spectra.

Typically, high-lying vibrational states near dissociation constitute a small fraction of the entire vibrational manifold, but this is not the case for the weakly bound magnesium dimer, which has a shallow minimum on the ground-state PEC at re = 3.89039 Å (20) and a tiny dissociation energy De of 430.472(500) cm−1 (20, 21). If the five extra levels, which have been speculated about, truly existed, they would represent more than a quarter of the entire vibrational manifold in the ground electronic state. Furthermore, without precise knowledge of the ground-state PEC of Mg2, especially its long-range part that determines the positions of the high-lying vibrational states near the dissociation threshold, one cannot accurately interpret the aforementioned ultracold and collisional phenomena involving interacting magnesium atoms. It is intriguing why a seemingly docile main group diatomic continues to challenge state-of-the-art spectroscopic techniques. The experimental difficulties in detecting the elusive v″ > 13 states of the magnesium dimer originate from several factors, including small energy gaps between high-lying vibrations that are comparable to rotational spacings (12, 25), resulting in overlapping spectral lines, and unfavorable signal-to-noise ratio in the existing LIF spectra (20). Rotational effects complicate the situation even more, since, in addition to affecting line intensities (20, 22, 25), they may render the high-lying vibrational states of Mg2 unbound. All of these and similar difficulties prompted Knöckel et al. (20, 21) to conclude that experimental work alone is insufficient and that accurate theoretical calculations are needed to guide further analysis of the ground-state PEC and rovibrational states of Mg2, especially the elusive v″ > 13 levels near the dissociation threshold.

Unfortunately, there have only been a handful of theoretical investigations attempting to determine the entire vibrational manifold of the magnesium dimer. This is, at least in part, related to the intrinsic complexity of the underlying electronic structure and difficulties with obtaining an accurate representation of the ground-state PEC using purely ab initio quantum-chemical means. At the Hartree-Fock theory level, which neglects electron correlation and dispersion interactions, Mg2 remains unbound. As demonstrated in this work, one needs to go to much higher theory levels, incorporate high-order many-electron correlation effects, including valence as well as inner-shell electrons, and use large, carefully calibrated, one-electron basis sets to accurately capture the relevant physics and obtain a reliable description of the X1Σg+ potential and of the corresponding rovibrational manifold [see (26) for a detailed discussion and historical account, including references to the earlier quantum chemistry computations for the magnesium dimer]. Ab initio quantum mechanical calculations for the A1Σu+ PEC, the rovibrational states supported by it, and the X1Σg+−A1Σu+ electronic transition dipole moment function, needed to interpret and aid the photoabsorption and LIF experiments using purely theoretical means, are similarly challenging, and the present study shows this too.

The initial theoretical estimates of the number of vibrational states supported by the X1Σg+ potential ranged from 18 to 20 (27), while the more recent ab initio quantum chemistry computations based on the various levels of coupled-cluster (CC) theory (28), reported in (26, 29), suggested that the highest vibrational level of 24Mg2 is v″ = 18. Among the previous theoretical studies, only Amaran et al. (29) considered the A1Σu+ state involved in the photoabsorption and LIF experiments and included rotational effects, but they have not provided any information about the calculated rovibrational term values other than the root mean square deviations (RMSDs) relative to the experimental data of Balfour and Douglas (12). Furthermore, as demonstrated in our recent benchmark study (26), where a large number of CC methods were tested using the X1Σg+ PEC of the magnesium dimer and the rotationless term values of 24Mg2 as examples, and consistent with the earlier calculations (30, 31), the popular CCSD(T) approximation (32) exploited in (29) could not possibly produce the small RMSD value reported in (29), of 1.3 cm−1, for the rovibrational manifold of Mg2 in its ground electronic state; the value on the order of a dozen cm−1 would be more appropriate (26) (see fig. S1). Similar remarks apply to the A1Σu+ state, which was treated in (29) using the low-level variant of the linear-response CC theory (33), resulting in noticeable deviations from the experimentally derived A1Σu+ potential shown in figure 4 of (21). To simulate and properly interpret the A1Σu+→X1Σg+ LIF spectra obtained in (20) using purely theoretical means, one needs much higher accuracy levels in the computations of line positions and robust information about line intensities, which has not been obtained in the previous quantum chemistry studies.

The call for a reliable ab initio computation of the ground-state PEC and rovibrational states of Mg2, including the v″ > 13 levels that have eluded experimentalists for decades, expressed by Knöckel et al. (20, 21), is answered in the present work. We report the highly accurate PECs for the ground, X1Σg+, and excited, A1Σu+, electronic states, involved in the previous experimental investigations of Mg2 (12, 20–22), obtained with state-of-the-art ab initio quantum chemistry, and use them to determine the corresponding rovibrational manifolds. Consistent with the conclusions of our recent benchmark study (26), to obtain a highly accurate representation of the ground-state PEC, we combine the numerically exact description of the valence electron correlation effects, provided by the full configuration interaction (CI) approach, with the nearly exact description of subvalence correlations involving all electrons but the 1s shells of Mg atoms offered by the CC theory with a full treatment of singly, doubly, and triply excited clusters, abbreviated as CCSDT (34, 35). Our computational protocol for the A1Σu+ excited state, which we did not consider in (26), is similar, except that to capture subvalence electron correlation effects in this state we adopt one of the carefully chosen approximations to the equation-of-motion (EOM) CC theory (36) with singles, doubles, and triples (EOMCCSDT) (37, 38) belonging to the completely renormalized CR-EOMCCSD(T) family (39), which is considerably more affordable than EOMCCSDT without substantial loss of accuracy. As in the case of the X1Σg+ state, the remaining electron correlations originating from the valence shells are captured using full CI. After thorough examination of basis set effects (see the Supplementary Materials for details), to obtain reasonably well converged rovibrational and LIF spectra, needed to correctly interpret the available experimental or experimentally derived data (12, 17, 20, 21), and to accurately describe the relevant many-electron correlation effects, the electronic structure calculations reported in this work rely on the carefully calibrated augmented polarized valence and weighted core-valence correlation-consistent bases of quadruple-ζ quality developed in (40), designated as aug-cc-pV(Q + d)Z and aug-cc-pwCVQZ, respectively. We used these basis sets in our earlier benchmark calculations for the ground state of the magnesium dimer (26) but only for the methods up to CCSDT, i.e., the important post-CCSDT electron correlation effects were treated in (26) with smaller, less saturated, bases.

To make our comparisons with the experiment more complete, for each of the two electronic potentials considered in this study, we examine both the most abundant 24Mg2 species and the 24Mg25Mg, 24Mg26Mg, 25Mg2, 25Mg26Mg, and 26Mg2 isotopologs (to our knowledge, rovibrational levels of the Mg2 species other than 24Mg2 have not been calculated using ab initio potentials before). We combine the above information with the X1Σg+−A1Σu+ electronic transition dipole moment function resulting from the same valence full CI computations as used in the PEC determination to accurately simulate the LIF spectra reported in (20, 21), including line positions and the corresponding line intensities, as defined via the Einstein coefficients, and provide the long-awaited theoretical guidance for the possible experimental detection of the A1Σu+(v', J') → X1Σg+(v″, J″) rovibronic transitions involving the v″ > 13 levels.

RESULTS

The most essential numerical information, generated in the present study using the computational protocol described in Materials and Methods, is summarized in Figs. 1 to 3 and Tables 1 to 3. All of the numerical data supporting the content and conclusions of this work are included in the main text and compiled in the Supplementary Materials document and data files S1 and S2 attached to it. In describing and discussing our results, we begin with the PECs and rovibrational term values characterizing the X1Σg+ and A1Σu+ states of the magnesium dimer, focusing on a comparison of our ab initio calculations with the available experimental and experimentally derived data reported in (12, 20, 21). Next, we compare the experimental LIF spectra reported in (20, 21) with those resulting from our computations and suggest potential avenues for detection of the elusive v″ > 13 levels of the magnesium dimer. Auxiliary information, which complements the discussion in this section, including further comments on the accuracy and convergence characteristics of the computational protocol used in the present study, the effect of isotopic substitution on the calculated rovibrational term values, the discussion of the validity of the Franck-Condon analysis adopted in (20) to examine the LIF spectra reported in (20, 21), and the lifetimes for predissociation by rotation characterizing quasi-bound rovibrational states supported by the X1Σg+ potential, is provided in the Supplementary Materials.

Fig. 1 The wave functions of the high-lying, purely vibrational, states of 24Mg2 and the underlying X1Σg+ potential.

The last experimentally observed v″ = 13 level is marked in blue, the predicted v″ = 14 to 18 levels are marked in green, and the ab initio X1Σg+ PEC obtained in this study is marked by a long-dashed black line. The inset is a Birge-Sponer plot comparing the rotationless G(v″ + 1) − G(v″) energy differences as functions of v″ + ½ obtained in this work (black circles) with their experimentally derived counterparts (red open squares) based on the data reported in (12) (v″ = 0 to 12) and (17) (v″ = 13; cf. also Table 1). The red solid line is a linear fit of the experimental points.

Fig. 2 Schematics of the pump, X1Σg+(v″ = 5, J″ = 10) → A1Σu+(v′ = 3, J′ = 11), and fluorescence, A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12), processes resulting in the LIF spectrum for 24Mg2 shown in figure 3 of (20).

The X1Σg+ and A1Σu+ PECs and the corresponding X1Σg+(v″ = 5, J″ = 10) and A1Σu+(v′ = 3, J′ = 11) rovibrational wave functions were calculated in this work. The A1Σu+ PEC was shifted to match the experimentally determined adiabatic electronic excitation energy Te of 26,068.9 cm−1 (21) (see Materials and Methods for the details).

Fig. 3 The A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) LIF spectrum of 24Mg2.

(A) Comparison of the experimental A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) fluorescence progression [black solid lines; adapted from figure 3 of (20) with the permission of AIP Publishing] with its ab initio counterpart obtained in this work (red dashed lines). The theoretical line intensities were normalized such that the tallest peaks in the calculated and experimental spectra corresponding to the v″ = 5 P12 line match. (B) Magnification of the low-energy region of the LIF spectrum shown in (A), with red solid lines representing the calculated transitions. The blue arrows originating from the v″ = 13 label indicate the location of the experimentally observed v″ = 13 P12/R10 doublet. The blue arrows originating from the v″ = 14 and 15 labels point to the most probable locations of the corresponding P12/R10 doublets. Spectral lines involving v″ = 16 and 17 are buried in the noise (see also Table 3).

Table 1 Comparison of the ab initio (Calc.) and experimentally derived (Expt.) rovibrational G(v″, J″) energies for selected values of J″ characterizing 24Mg2 in the ground electronic state (in reciprocal centimeter), along with the corresponding dissociation energies De (in reciprocal centimeter) and equilibrium bond lengths re (in angstrom).

The G(v″, J″) energies calculated using the ab initio X1Σg+ PEC defined by Eq. 1 are reported as errors relative to experiment, whereas De and re are the actual values of these quantities. If the experimental G(v″, J″) energies are not available, we provide their calculated values in square brackets. Quasi-bound rovibrational levels are given in italics. Horizontal bars indicate term values not supported by the X1Σg+ PEC.

View this table:
Table 2 Comparison of the rovibrational G(v″, J″) energies obtained using the ab initio X1Σg+ PEC defined by Eq. 1 (Calc.) and its X-representation counterpart constructed in (20) (X-rep.) for selected values of J″ characterizing 24Mg2 in the ground electronic state (in reciprocal centimeter), along with the corresponding dissociation energies De (in reciprocal centimeter) and equilibrium bond lengths re (in angstrom).

The G(v″, J″) energies calculated using the ab initio X1Σg+ PEC are reported as errors relative to the X-representation data, whereas De and re are the actual values of these quantities. If a given G(v″, J″) state corresponding to our ab initio X1Σg+ PEC is not supported by the X-representation potential of (20), we provide its energy in square brackets. Quasi-bound rovibrational levels are given in italics. Horizontal bars indicate term values not supported by the X1Σg+ PEC.

View this table:
Table 3 Comparison of the theoretical line positions of the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) fluorescence progression in the LIF spectrum of 24Mg2 calculated in this work with experiment.

All line positions are in reciprocal centimeter. The available experimental values are the actual line positions, whereas our calculated results are errors relative to experiment. If the experimentally determined line positions are not available, we provide their calculated values in square brackets. Horizontal bars indicate term values not supported by the X1Σg+ PEC.

View this table:

PECs and rovibrational states

As shown in Table 1, our ab initio X1Σg+ PEC reproduces the experimentally derived dissociation energy De and equilibrium bond length re of Mg2 (20, 21) to within 0.9 cm−1 (0.2%) and 0.003 Å (0.07%), respectively. These high accuracies in describing De and re are reflected in our calculated rovibrational term values of 24Mg2 and its isotopologs, which are in very good agreement with the available experimental information (12, 20, 21). As shown in the spreadsheets included in data file S1, the RMSDs characterizing our ab initio G(v″, J″) values for 24Mg2 relative to their experimentally determined counterparts, reported in (12) for v″ < 13 and (20, 21) for v″ < 14, are 1.1 cm−1, when the spectroscopic data from (12) are used, and 1.5 cm−1, when we rely on (20, 21) instead. At the same time, the maximum unsigned errors in our calculated G(v″, J″) values relative to the experiment do not exceed ~2 cm−1, even when the quasi-bound states above the potential asymptote arising from centrifugal barriers are considered. Although the experimental information about the G(v″, J″) values characterizing other Mg2 isotopologs is limited to 24Mg25Mg, 24Mg26Mg, and 26Mg2 and includes very few v″ values (20, 21), the RMSDs relative to the experiment resulting from our calculations are similarly small (1.0 cm−1 for 24Mg25Mg, 1.2 cm−1 for 24Mg26Mg, and 0.6 cm−1 for 26Mg2; cf. the Supplementary Materials).

Further insights into the quality of our ab initio calculations for the ground-state PEC can be obtained by comparing the resulting rovibrational term values with their counterparts determined using the most accurate, experimentally derived, analytical forms of the X1Σg+ potential to date constructed in (20). In the discussion below, we focus on the so-called X-representation of the ground-state PEC developed in (20), which the authors of (20) regard as a reference potential in their analyses (see Table 2). We recall that the X-representation of the ground-state PEC of the magnesium dimer was obtained by simultaneously fitting the X1Σg+ and A1Σu+ PECs to a large number of the experimentally determined A1Σu+(v′,J′)→X1Σg+(v″,J″) rovibronic transition frequencies and extrapolating the resulting X1Σg+ PEC to the asymptotic region using the theoretical C6 (23), C8 (24), and C10 (24) coefficients. As shown in Table 2, our ab initio G(v″, J″) energies characterizing the most abundant 24Mg2 isotopolog are in very good agreement with those generated using the X-representation of the ground-state PEC developed in (20). When all of the rovibrational bound states supported by both potentials are considered, the RMSD and the maximum unsigned error characterizing our ab initio G(v″, J″) values for 24Mg2 relative to their counterparts arising from the X-representation are 1.3 and 2.0 cm−1, respectively. What is especially important in the context of the present study is that our ab initio ground-state PEC and the state-of-the-art analytical fit to the experimental data defining the X-representation, constructed in (20), bind the v″ = 18 level if the rotational quantum number J″ is not too high (see the discussion below).

The high quality of our calculated G(v″, J″) values and spacings between them, which can also be seen in Tables 1 and 2 and Fig. 1, allows us to comment on the existence of the v″ > 13 levels that have escaped experimental detection for decades. As already alluded to above and as shown in Table 2 and Fig. 1, our ab initio X1Σg+ PEC supports the same number of rotationless vibrational levels as the latest experimentally derived PEC defining the X-representation (20), which for the most abundant 24Mg2 isotopolog is 19 (see the Supplementary Materials for the information about the remaining Mg2 species). Table 1, which compares the rovibrational term values of 24Mg2 resulting from our ab initio calculations for the representative rotational quantum numbers ranging from 0 to 80 with the available experimental data, shows that the elusive high-lying states with v″ > 13 quickly become unbound as J″ increases, so by the time J″ = 20, the v″ = 15 to 18 levels are no longer bound (see fig. S3 for a graphical representation of the J″ = 20, 40, 60, and 80 effective potentials including centrifugal barriers characterizing the rotating 24Mg2 molecule, along with the corresponding vibrational wave functions and information about the lifetimes for predissociation by rotation associated with tunneling through centrifugal barriers characterizing quasi-bound states). In fact, according to our ab initio data compiled in the Supplementary Materials, the maximum rotational quantum number that allows for at least one bound rovibrational state decreases with v″, from J″ = 68 for v″ = 0 to J″ = 4 for v″ = 18, with all states becoming quasi-bound or unbound when J″ ≥ 70, when the most abundant 24Mg2 isotopolog is considered. In general, as shown in fig. S3 and the lifetime data compiled in the Supplementary Materials, the mean lifetimes for predissociation by rotation characterizing quasi-bound states with a given J″ rapidly decrease as v″ becomes larger. They decrease equally fast when J″ increases and v″ is fixed. These observations imply that the spectroscopic detection of the high-lying vibrational states of Mg2 can only be achieved if the molecule does not rotate too fast (cf. Table 1 and fig. S3). We could not find any information regarding the timescales involved in the LIF experiments carried out by Knöckel et al. (20). However, a comparison of our ab initio–determined quasi-bound rovibrational states, including their energies and lifetimes compiled in the Supplementary Materials, with the observed rovibronic transitions reported in the supplementary material of (21) suggests that the mean lifetimes for predissociation by rotation characterizing quasi-bound states seen in the experimentally resolved LIF spectral lines are on the order of 0.1 ns or longer.

As shown in Fig. 1, where we plot the wave functions of the high-lying, purely vibrational, states of 24Mg2, starting with the last experimentally observed v″ = 13 level, along with the X1Σg+ PEC obtained in our ab initio calculations, the v″ = 18 state, located only 0.2 cm−1 below the potential asymptote, is barely bound (see also Table 1). This makes the existence of an additional, v″ = 19, level for the most abundant isotopolog of the magnesium dimer unlikely. Further insights into the number of purely vibrational bound states of 24Mg2 supported by the X1Σg+ PEC are provided by the inset in Fig. 1, where we plot the rotationless G(v″ + 1) − G(v″) energy differences, resulting from the ab initio calculations reported in this work and the experiment, as a function of v″ + ½ (the Birge-Sponer plot). Fitting the experimental data to a line, i.e., assuming a Morse potential, results in v″ = 16 being the last bound vibrational level of 24Mg2. Although the deviation from the Morse potential, as predicted by our ab initio calculations, is not as severe as in the case of Be2 (11), it is large enough to result in the v″ = 17 and 18 states becoming bound, emphasizing the importance of properly describing the long-range part of the PEC.

As shown in Table 1 and Fig. 1, the G(v″ + 1) − G(v″) vibrational spacings rapidly decrease with increasing v″, from 47.7 cm−1 or 68.6 K for v″ = 0 to 11.7 cm−1 or 16.8 K for v″ = 12, and to 0.8 cm−1 or 1.2 K for v″ = 17, when 24Mg2 is considered. This means that at regular temperatures all vibrational levels of the magnesium dimer, which is a very weakly bound system, are substantially populated, making selective probing of the closely spaced higher-energy states, including those with v″ > 13, virtually impossible, since practically every molecular collision (e.g., with another dimer) may result in a superposition of many rovibrational states, with some breaking the dimer apart. At room temperature, for example, the cumulative population of the v″ > 13 states of 24Mg2, determined using the normalized Boltzmann distribution involving all rotationless levels bound by the X1Σg+ potential, of about 12%, is comparable to the populations of the corresponding low-lying states (16% for v″ = 0, 13% for v″ = 1, and 10% for v″ = 2). The situation changes in the cold/ultracold regime, where the available thermal energies, which are on the order of millikelvin or even microkelvin, are much smaller than the vibrational spacings, even when the high-lying states with v″ > 13 near the dissociation threshold are considered, suppressing collisional effects and allowing one to probe the long-range part of the ground-state PEC, where the v″ > 13 states largely localize (cf. Fig. 1). This makes the accurate characterization of the v″ > 13 bound and quasi-bound states provided by the high-level ab initio calculations reported in this work relevant to the applications involving cold/ultracold Mg atoms separated by larger distances in magneto-optical traps [see, e.g., (6)].

The accuracy of our ab initio description of the more strongly bound A1Σu+ electronic state [De = 9414 cm−1 and re = 3.0825 Å (21); cf. Fig. 2 for the corresponding PEC], which we need to consider to simulate the LIF spectra, is consistent with that obtained for the weakly bound ground state. For example, the errors relative to the experiment (21) resulting from our calculations of the dissociation energy De and equilibrium bond length re are 0.91% (86 cm−1) and 0.2% (0.006 Å), respectively (see the Supplementary Materials). This high accuracy of our ab initio A1Σu+ PEC is reflected in the excellent agreement between the 24Mg2 G(v′, J′) values obtained in this work and their experimentally derived counterparts reported in (12, 21). In particular, the RMSDs characterizing our rovibrational term values in the A1Σu+ state relative to the data of Balfour and Douglas (12) and Knöckel et al. (21) are only 3.2 and 4.5 cm−1, respectively, which is a major improvement over the RMSD of 30 cm−1 reported in (29). They are similarly small for the rovibrational states supported by the A1Σu+ potential that characterize the remaining, experimentally observed, 24Mg25Mg, 24Mg26Mg, and 26Mg2 isotopologs examined in (21) (3.7, 4.1, and 4.1 cm−1, respectively; cf. the Supplementary Materials). According to our ab initio calculations using the computational protocol described in Materials and Methods, the total number of vibrational states supported by the A1Σu+ potential well for the most abundant 24Mg2 species is 169 (see data file S1).

LIF: Ab initio theory versus experiment

The most compelling evidence for the predictive power of our ab initio electronic structure and rovibrational calculations is the nearly perfect reproduction of the experimental A1Σu+→X1Σg+ LIF spectrum reported in (20, 21), shown in Fig. 3 and Table 3, with further information provided in the Supplementary Materials. Figure 2 uses our calculated X1Σg+ and A1Σu+ PECs and the corresponding rovibrational wave functions to illustrate the photoexcitation and fluorescence processes that resulted in the experimental LIF spectrum shown in figure 3 of (20), which is reproduced in Fig. 3A. This particular spectrum represents the fluorescence progression from the A1Σu+(v′ = 3, J′ = 11) state of 24Mg2, populated by laser excitation from the X1Σg+(v″ = 5, J″ = 10) state, to all accessible X1Σg+(v″, J″) rovibrational levels, resulting in the P12/R10 doublets that correspond to J″ = 12 for the P branch and J″ = 10 for the R branch. Figure 3 and Table 3 compare the experimentally observed A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) transitions with the corresponding line positions (Fig. 3 and Table 3) and intensities (Fig. 3) resulting from our ab initio calculations. The only adjustment that we made to produce the theoretical LIF spectrum shown in Fig. 3 and Table 3 was a uniform shift of the entire A1Σu+ PEC obtained in our ab initio computations to match the experimentally determined adiabatic electronic excitation energy Te of 26,068.9 cm−1 (21) (see Materials and Methods for the details). Other than that, the theoretical LIF spectrum in Fig. 3 and Table 3 relies on the raw ab initio electronic structure and rovibrational data.

The notable agreement between the theoretical and experimental LIF spectra shown in Fig. 3A and Table 3, with differences in line positions not exceeding 1 to 1.5 cm−1 and with virtually identical intensity patterns, suggests that our predicted transition frequencies involving the elusive v″ > 13 states are very accurate, allowing us to provide guidance for their potential experimental detection in the future. Before discussing our suggestions in this regard, we note that owing to our ab initio calculations, we can now locate the previously unidentified P12/R10 doublets involving the v″ > 13 states within the experimental LIF spectrum reported in figure 3 of (20). Indeed, as shown in Fig. 3 and Table 3, the LIF spectrum corresponding to the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) transitions contains the P12/R10 doublets involving the v″ = 0 to 16 states and the R10 line involving the v″ = 17 state. The A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 17, J″ = 12) and A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 18, J″ = 10,12) transitions are absent, since the v″ = 17, J″ = 12 and v″ = 18, J″ = 10 and 12 states are unbound (see the Supplementary Materials), but they could potentially be observed if one used different initial A1Σu+(v′, J′) states (see the discussion below).

As one can see by inspecting data file S2 and Fig. 3, and consistent with the remarks made by Knöckel et al. in (20), the experimental detection of the P12/R10 doublets involving v″ > 13, when transitioning from the A1Σu+(v′ = 3, J′ = 11) state, was hindered by the unfavorable signal-to-noise ratio (transitions to the v″ = 16 and 17 states exhibit low Einstein coefficients) and the presence of overlapping lines outside the P12/R10 progression, originating from collisional relaxation effects (20) and having similar (v″ = 15) or higher (v″ = 14) intensities. To fully appreciate this, in Fig. 3B, we magnified the region of the LIF spectrum recorded in (20) that contains the calculated A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 13 to 16, J″ = 10,12) and A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 17, J″ = 10) transitions. As shown in Fig. 3 and Table 3, the identification of the P12/R10 doublets corresponding to the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 0 to 13, J″ = 10,12) transitions is unambiguous. The observed and calculated line positions and intensities and line intensity ratios within every doublet match each other very closely. Figure 3B demonstrates that the identification of the remaining doublets in the P12/R10 progression is much harder. On the basis of our ab initio work and taking into account the fact that our calculated line positions may be off by about 1 cm−1 (cf. Table 3), the v″ = 14 P12/R10 doublet, marked in Fig. 3B by the blue arrows originating from the v″ = 14 label, is largely hidden behind the higher-intensity feature that does not belong to the P12/R10 progression and that most likely originates from collisional relaxation (20). Because of our calculations, we can also point to the most likely location of the v″ = 15 P12/R10 doublet in the LIF spectrum recorded in (20) (see the blue arrows originating from the v″ = 15 label in Fig. 3B). Doing this without backing from the theory is virtually impossible due to the presence of other lines near the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 15, J″ = 10,12) transitions having similar intensities. As shown in Fig. 3B, the situation with the remaining A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 16, J″ = 10,12) and A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 17, J″ = 10) transitions is even worse, since they have very low Einstein coefficients that hide them in the noise.

Knöckel et al. (20) also suggested that the difficulties with detecting the P12/R10 doublets involving the v″ = 14 and 15 states, which have higher Franck-Condon factors than those characterizing the experimentally observed A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 0, J″ = 10,12) transitions, might be related to the variation of the X1Σg+−A1Σu+ electronic transition dipole moment function μz(X−A)(r) with the internuclear separation r and limitations of the Franck-Condon principle, but our ab initio calculations do not confirm this. As shown in the Supplementary Materials (see fig. S2), in the region of r values where the respective rovibrational wave functions, χv",J"(X)(r), with v″ = 14,15 and J″ = 10,12 for the X1Σg+ state, and χv',J'(A)(r), with v′ = 3 and J′ = 11 for the A1Σu+ state, overlap, changes in μz(X−A)(r) do not exceed 3%, i.e., the Frank-Condon analysis is well justified. Furthermore, μz(X−A)(r) does not vary too much, even when the entire r = 2.2 to 100.0 Å region examined in this work is considered. In agreement with the analysis presented in (20), our calculated Franck-Condon factors for the v″ = 14 and 15 P12/R10 doublets are higher than those characterizing the analogous v″ = 0 lines, but, as demonstrated in data file S2, the same holds for the respective Einstein coefficients. Thus, it is the presence of densely spaced and overlapping lines outside the P12/R10 progression having similar or higher intensities than the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 14,15, J″ = 10,12) transitions that makes the experimental identification of the v″ = 14 and 15 P12/R10 doublets very hard.

Theory-inspired avenues for detection of elusive states

In general, our ab initio calculations carried out in this work indicate that under the constraints of the LIF experiments reported in (20, 21), where the authors populated the A1Σu+(v′, J′) states with v′ = 1 to 46, the X1Σg+(v″, J″) states with v″ = 14 to 18 cannot be realistically detected because of very small Franck-Condon factors and Einstein coefficients characterizing the corresponding A1Σu+(v′, J′) → X1Σg+(v″, J″) transitions (see data file S2). As shown in Fig. 1, the v″ = 14 to 18 states are predominantly localized in the long-range r = 8 to 16 Å region. At the same time, as illustrated in Fig. 2, the potential well characterizing the electronically excited A1Σu+ state is much deeper and shifted toward shorter internuclear separations compared to its X1Σg+ counterpart. Thus, the only way to access the X1Σg+(v″, J″) states with v″ = 14 to 18 via fluorescence from A1Σu+ is by populating the high-lying A1Σu+(v′, J′) levels with v′ ≫ 46.

In an effort to assist the experimental community in detecting the elusive v″ = 14 to 18 vibrational levels, we searched for the A1Σu+(v′, J′) → X1Σg+(v″ = 14 to 18, J″ = J′ ± 1) transitions in the most abundant isotopolog of the magnesium dimer, 24Mg2, that would result in spectral lines of maximum intensity based on the Einstein coefficients compiled in data file S2. To ensure the occurrence of allowed transitions involving the last, v″ = 18, level, which for 24Mg2 becomes unbound when J″ > 4, we focused on the J″ values not exceeding 4, i.e., the fluorescence from the A1Σu+(v′, J′) states with J′ = 1, 3, and 5. According to our calculations, the optimum v′ values for observing the v″ = 14 to 18, J″ ≤ 4 states via the LIF spectroscopy are in the neighborhood of v′ = 60, 66 to 69, and 74 to 84 for v″ = 14; 72 to 75 and 80 to 91 for v″ = 15; 79 to 82 and 88 to 100 for v″ = 16; 88, 89, and 97 to 111 for v″ = 17; and 109 to 129 for v″ = 18 [see data file S2 for the details of all allowed rovibronic transitions in 24Mg2 involving the X1Σg+ and A1Σu+ states, including, in particular, the relevant X1Σg+(v″, J″ ≤ 4) → A1Σu+(v′, J′) pump and A1Σu+(v′, J′ = 1,3,5) → X1Σg+(v″ = 14 to 18, J″ ≤ 4) fluorescence processes]. In determining these optimum v′ values, we chose the cutoff value of 1.0 × 107 Hz in the Einstein coefficients, which is similar to the Einstein coefficients calculated for the most intense v″ = 5 P12/R10 doublet in the experimental LIF spectrum shown in figure 3 of (20), reproduced in Fig. 3A. Our predicted A1Σu+(v′, J′ = 1,3,5) → X1Σg+(v″ = 14 to 18, J″ ≤ 4) fluorescence frequencies resulting from the aforementioned optimum v′ ranges, which might allow one to detect the v″ = 14 to 18 states of 24Mg2 via a suitably designed LIF experiment, are estimated at about 33,360, 33,740 to 33,910, and 34,150 to 34,530 cm−1 for v″ = 14; 34,050 to 34,190 and 34,390 to 34,710 cm−1 for v″ = 15; 34,350 to 34,460 and 34,640 to 34,880 cm−1 for v″ = 16; 34,630 to 34,660 and 34,830 to 35,000 cm−1 for v″ = 17; and 34,990 to 35,100 cm−1 for v″ = 18 [given the 86 cm−1 error in the calculated De characterizing the A1Σu+ state and the RMSD of ~3 to 5 cm−1 in our 24Mg2 G(v′, J′) values relative to the spectroscopic data of (12, 21), the above frequency ranges may have to be shifted by a dozen or so cm−1].

DISCUSSION

We used state-of-the-art ab initio quantum-mechanical methodologies to address a half-century-old enigma regarding the v″ = 14 to 18 vibrational states of the magnesium dimer. We provided the highly accurate ground-state PEC and rovibrational term values of 24Mg2 and its less abundant 24Mg25Mg, 24Mg26Mg, 25Mg2, 25Mg26Mg, and 26Mg2 isotopologs. We demonstrated that the X1Σg+ PEC supports rovibrational levels of 24Mg2 up to v″ = 18, although the elusive v″ > 13 states become unbound as the rotational quantum number J″ increases, which contributes to difficulties with their experimental detection. We also obtained an accurate representation of the A1Σu+ potential, which, according to our calculations, supports 169 vibrational states of 24Mg2, and, with the help of the ab initio electronic transition dipole moment function, determined in this study as well, accurately simulated the LIF spectra recorded in (20, 21), including line positions and intensities. Our work provides the long-awaited guidance for possible experimental identification of rovibronic transitions involving the v″ > 13 levels that have eluded scientists for five decades.

We hope that this study will fuel new spectroscopic investigations of the challenging Mg2 species and its heavier group 2 analogs, which are important in a variety of phenomena at the intersection of chemistry and atomic, molecular, and optical physics. A few years ago, ab initio calculations (8) combined with spectroscopic analyses (7, 10) led to the discovery of the elusive 12th vibrational level of the beryllium dimer. By dealing with five similarly challenging states in a system three times larger than Be2, we demonstrated that the predictive power of modern ab initio quantum chemistry is no longer limited to small few-electron species.

MATERIALS AND METHODS

Ab initio electronic structure calculations

The goal of the ab initio electronic structure calculations performed in this study was to obtain highly accurate X1Σg+ and A1Σu+ PECs of the magnesium dimer and the corresponding X1Σg+−A1Σu+ transition dipole moment function μz(X−A)(r) involved in the photoabsorption and LIF experiments reported in (12, 20–22). In the case of the ground-state PEC, we combined the numerically exact description of the valence electron correlation effects provided by full CI with the high-level description of subvalence correlations involving all electrons but the 1s shells of Mg atoms obtained using CCSDT (34, 35). Thus, the X1Σg+ PEC of Mg2 reported in this work was obtained by adopting the composite schemeEX1Σg+=EX1Σg+(CCSDT/AwCQZ)+(EX1Σg+(Full CI/A(Q+d)Z)−EX1Σg+(CCSDT/A(Q+d)Z))(1)

The first term on the right-hand side of Eq. 1 denotes the total electronic energy obtained in the full CCSDT calculations correlating all electrons other than the 1s shells of the Mg monomers and using the aug-cc-pwCVQZ basis set developed in (40), abbreviated in this section and in the Supplementary Materials as AwCQZ. The second and third terms on the right-hand side of Eq. 1, which represent the difference between the frozen-core full CI and CCSDT energies obtained using the aug-cc-pV(Q + d)Z basis of (40), abbreviated in this section and in the Supplementary Materials as A(Q + d)Z, correct the nearly all-electron CCSDT/AwCQZ energy for the valence correlation effects beyond CCSDT. The A(Q + d)Z and AwCQZ basis sets were taken from the Peterson group’s website (41). We used these bases rather than their standard aug-cc-pVnZ and aug-cc-pCVnZ counterparts, since it has been demonstrated that the aug-cc-pV(n + d)Z and aug-cc-pwCVnZ basis set families, including A(Q + d)Z and AwCQZ, accelerate the convergence of bond lengths, dissociation energies, and spectroscopic properties of magnesium compounds (26, 40). The aug-cc-pV(T + d)Z, aug-cc-pwCVTZ, and aug-cc-pwCV5Z bases (40), abbreviated in this section and in the Supplementary Materials as A(T + d)Z, AwCTZ, and AwC5Z, respectively, and used in the auxiliary calculations discussed in section S1 to demonstrate the convergence of our computational protocol with respect to the basis set size (see tables S1 and S2), were taken from the Peterson group’s website (41) as well.

As shown in section S1, the AwCQZ and A(Q + d)Z bases are large and rich enough to provide spectroscopic properties of the magnesium dimer that can be regarded as reasonably well converged with respect to the basis set size, to within ~0.1 to 2 cm−1 for the experimentally observed v″ ≤ 13 levels and ~3 to 5 cm−1 for the remaining high-lying vibrational states and De (see, e.g., table S2). Ideally, one would like to improve these results further by extrapolating, for example, the nearly all-electron CCSDT energetics in Eq. 1, which are responsible for the bulk of the many-electron correlation effects in Mg2, to the complete basis set (CBS) limit. Unfortunately, a widely used two-point CBS extrapolation (42) based on the subvalence CCSDT/AwCTZ and CCSDT/AwCQZ data, which are the only CCSDT data of this type available to us, to determine the CBS counterpart of the first term on the right-hand side of Eq. 1 would not be reliable enough. As demonstrated in (26) and as elaborated on in section S1 (see table S2), a CBS extrapolation using the AwCTZ and AwCQZ basis sets worsens, instead of improving, the De, re, and vibrational term values of the magnesium dimer compared to the unextrapolated results using the AwCQZ basis. As shown in table S2, the CBS extrapolation using the AwCQZ and AwC5Z basis sets would be accurate enough, but the CCSDT/AwC5Z calculations for the magnesium dimer correlating all electrons but the 1s shells of Mg atoms are prohibitively expensive. One could try to address this concern by replacing CCSDT in Eq. 1 by the more affordable CCSD(T) approach (32), resulting inE˜X1Σg+=EX1Σg+(CCSD(T)/AwCQZ)+(EX1Σg+(Full CI/A(Q+d)Z)−EX1Σg+(CCSD(T)/A(Q+d)Z))(2)but, as explained in section S2, the computational protocol defined by Eq. 2 is not sufficiently accurate for the spectroscopic considerations reported in this work due to the inadequate treatment of triples by the baseline CCSD(T) approximation (cf. fig. S1). For all these reasons, we have to rely on Eq. 1, in which we use CCSDT, not CCSD(T), and finite (albeit large and carefully optimized) AwCQZ and A(Q + d)Z basis sets rather than the poor-quality CBS extrapolation from the CCSDT/AwCTZ and CCSDT/AwCQZ information.

In principle, one could extend the above composite scheme, given by Eq. 1, to the electronically excited A1Σu+ state by replacing CCSDT in Eq. 1 with its EOMCCSDT counterpart (37, 38), but the nearly all-electron full EOMCCSDT calculations using the large AwCQZ basis set turned out to be prohibitively expensive for us. To address this problem, we resorted to one of the CR-EOMCCSD(T) approximations to EOMCCSDT, namely, CR-EOMCCSD(T),IA (39), which is capable of providing highly accurate excited-state PECs of near-EOMCCSDT quality at the small fraction of the cost. Thus, our composite scheme for the calculations of the A1Σu+ PEC was defined asEA1Σu+=EA1Σu+(CR-EOMCCSD(T),IA/AwCQZ)+(EA1Σu+(Full CI/A(Q+d)Z)−EA1Σu+(CR-EOMCCSD(T),IA/A(Q+d)Z))(3)where the first term on the right-hand side of Eq. 3 is the total electronic energy of the A1Σu+ state obtained in the CR-EOMCCSD(T), IA/AwCQZ calculations correlating all electrons other than the 1s shells of the Mg monomers and the next two terms correct the nearly all-electron CR-EOMCCSD(T),IA/AwCQZ calculations for the valence correlation effects beyond the CR-EOMCCSD(T),IA level using the difference of the full CI and CR-EOMCCSD(T),IA energies obtained with the A(Q + d)Z basis. Before deciding on the use of CR-EOMCCSD(T),IA, we tested other CR-EOMCC schemes (43) by comparing the resulting A1Σu+ potentials obtained using Eq. 3 and the corresponding rovibrational term G(v′, J′) values with the available experimentally derived data reported in (21, 44). Although all of these schemes worked well, the computational protocol defined by Eq. 3, with the CR-EOMCCSD(T),IA approach serving as a baseline method, turned out to produce the smallest maximum unsigned errors and RMSD values relative to experiment.

While the X1Σg+ and A1Σu+ PECs obtained in this study appear to be accurate enough for reproducing and interpreting the experimental A1Σu+→X1Σg+ LIF spectra reported in (20, 21), one might wonder whether the neglect of the post–Born-Oppenheimer and relativistic effects in our ab initio calculations could substantially affect our main conclusions. According to (20, 21), the non-adiabatic Born-Oppenheimer corrections (BOCs) for the X1Σg+ and A1Σu+ states and the mass-dependent adiabatic BOC for the X1Σg+ state are negligible. The adiabatic BOC for the A1Σu+ state, as defined in (21), may have to be accounted for, but, based on the numerical data reported in (21), its magnitude is well within the uncertainty of the ab initio calculations reported in this work. According to (30), special relativity reduces the dissociation energy De characterizing the X1Σg+ PEC by 4.3 cm−1, i.e., the relativistic effects change the De by about 1%. However, our preliminary analysis using the modified version of the ab initio protocol adopted in the present work, in which the valence full CI and CCSDT calculations using the A(Q + d)Z basis set and the nearly all-electron CCSDT/AwCQZ computations are replaced by their scalar-relativistic counterparts employing the third-order Douglas-Kroll (DK) Hamiltonian (45, 46) and the triple-ζ aug-cc-pV(T + d)Z-DK and aug-cc-pwCVTZ-DK bases (40), demonstrates that the number of bound vibrational states supported by the relativity-corrected X1Σg+ potential is exactly the same as in the case of the analogous nonrelativistic calculations using the A(T + d)Z and AwCTZ bases [the small negative differences between the relativity-corrected and nonrelativistic rotationless G(v″) values vary from <1 cm−1 or 0.8% for v″ = 0 to 2 to ~1% for the highest vibrational states near the corresponding dissociation thresholds]. Similar applies to the A1Σu+ PEC, where the effect of relativity on the De value, estimated using the triple-ζ DK analog of the quadruple-ζ nonrelativistic computational protocol adopted in this work, is 0.2%, but the total number of bound vibrational states supported by the nonrelativistic and relativity-corrected potentials remains the same. The ab initio vibrational spectra corresponding to the X1Σg+ and A1Σu+ electronic states obtained using the triple-ζ DK modification of the nonrelativistic protocol used in the present study also show that the effects of relativity on the rotationless G(v″ + 1) − G(v″) and G(v′ + 1) − G(v′) energy spacings do not exceed 0.4 cm−1 in the former case and 0.3 cm−1 in the case of the latter energy differences. Thus, while our preliminary findings regarding the small, but nonnegligible, effects of relativity need a thorough reexamination using both the larger basis sets, such as aug-cc-pV(Q + d)Z-DK and aug-cc-pwCVQZ-DK, and the various truncations in the DK Hamiltonian expansions, which may influence the calculated spectra too (46), and we will return to these issues in the future work, the X1Σg+ and A1Σu+ PECs obtained in the present study are sufficiently accurate to interpret and analyze the LIF spectra reported in (20, 21) and to comment on the corresponding rovibrational manifolds, especially for the ground electronic state.

All electronic structure calculations for Mg2 performed in this study, summarized in tables S3 to S5, were based on the tightly converged restricted Hartree-Fock (RHF) reference functions (the convergence criterion for the RHF density matrix was set up at 10−9). The valence full CI calculations for the X1Σg+ and A1Σu+ states were performed using the GAMESS package (47), whereas the valence and subvalence CCSDT computations for the X1Σg+ state were carried out with NWChem (48). The valence and subvalence CR-EOMCCSD(T),IA calculations for the A1Σu+ state were executed using the RHF-based CR-EOMCCSD(T) routines developed in (39), which take advantage of the underlying ground-state CC codes described in (49) and which are part of GAMESS as well. The GAMESS RHF-based CC routines (49) were also used to perform the CCSD(T) calculations needed to explore the basis set convergence and the viability (or the lack thereof) of the alternative to the CCSDT-based composite scheme given by Eq. 1, defined by Eq. 2 (see sections S1 and S2, especially table S2 and fig. S1). The convergence thresholds used in the post-RHF steps of the CC and EOMCC computations reported in this work were set up at 10−7 for the relevant excitation amplitudes and 10−7 hartree (0.02 cm−1) for the corresponding electronic energies. The default GAMESS input options that were used to define our full CI calculations guaranteed energy convergence to 10−10 hartree.

The grid of Mg-Mg separations r, at which the electronic energies of the X1Σg+ and A1Σu+ states reported in this study (cf. tables S3 and S4) were determined, was as follows: 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.1, 3.2, 3.3, 3.4, 3.6, 3.7, 3.8, 3.9, 4.0, 4.1, 4.2, 4.4, 4.6, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0, 6.4, 6.8, 7.2, 7.6, 8.0, 8.4, 8.8, 9.2, 9.6, 10.0, 11.0, 12.0, 13.0, 15.0, 20.0, 25.0, 30.0, and 100.0 Å. We adopted the same set of r values to determine the electronic transition dipole moment function μz(X−A)(r) between the X1Σg+ and A1Σu+ electronic states, needed to calculate LIF line intensities using the Einstein coefficients. The μz(X−A)(r) calculations reported in this work were performed using the valence full CI approach, as implemented in GAMESS, adopting the A(Q + d)Z basis set of (40) (see fig. S2 and table S5).

It is worth pointing out that our ab initio data points representing the X1Σg+ PEC calculated on the above grid of r values are consistent with the expected long-range physics. One can see this, for example, by comparing our ab initio electronic energies for the X1Σg+ state compiled in the last column of table S3 with their X-representation counterparts obtained using the potential parameters provided in table 1 of (20). Indeed, if we align the X-representation and our ab initio potentials such that the energies at r = 100.0 Å are identical (without this alignment, the X-representation and our ab initio energies at r = 100.0 Å calculated relative to the corresponding potential minima differ by 0.9 cm−1), the differences between the two PECs in the r > 8.5 Å region, where the X-representation potential has the form V(r)=De−Σm=03C6+2m/r6+2m, do not exceed 0.8 cm−1, rapidly approaching zero as r increases (r = 100.0 Å is large enough to define the asymptotic region; for example, the difference between the X-representation energies at r = 30.0 and 100.0 Å is only 0.004 cm−1; our ab initio calculations at the same two r values produce the numerically identical energy difference). In the r ≥ 20 Å region, where the X-representation energies are flat to within about 0.05 cm−1, the X-representation and shifted ab initio energies, as described above, differ by less than 0.01 cm−1.

Calculations of rovibrational term values and rovibronic transitions

The rovibrational term values, including bound and quasi-bound states supported by our ab initio X1Σg+ and A1Σu+ PECs defined by Eqs. 1 and 3, were computed by numerically integrating the radial Schrödinger equation from 2.2 to 100.0 Å using the Numerov-Cooley algorithm (50) available in the LEVEL16 code (51) [LEVEL16 uses the Airy function approach described in (52) to locate quasi-bound states]. The widths and the tunneling lifetimes for predissociation by rotation characterizing the quasi-bound rovibrational states supported by the X1Σg+ potential were calculated using LEVEL16 as well. In this case, we followed the semiclassical procedure described in (51) and implemented in LEVEL16, which requires numerical integrations between turning points in the classically allowed and classically forbidden regions of the relevant effective potentials including centrifugal barriers shown, for example, in fig. S3 [see (51) for further details]. To produce electronic energies V(r) on a dense grid of internuclear distances r with the step size of 0.001 Å, needed to perform the required numerical integrations and determine the corresponding equilibrium bond lengths re and dissociation energies De, we proceeded as follows. To obtain V(r) values every 0.001 Å in the r = 2.3 to 30.0 Å region, which excludes the innermost and outermost PEC parts defined by the 2.2 to 2.3 Å and 30.0 to 100.0 Å intervals, we used cubic splines available in LEVEL16, interpolating between pairs of nearest-neighbor r values used in the ab initio electronic structure calculations, starting from (2.3 Å, 2.4 Å) and ending up with (25.0 Å, 30.0 Å). To generate the equally densely spaced electronic energies in the innermost and outermost segments of each of the two PECs considered in this work, we resorted to analytical potential fits provided by the LEVEL16 code. In the case of the innermost parts of the X1Σg+ and A1Σu+ potentials, corresponding to the 2.2 to 2.3 Å interval, we used the formula V(r) = A + Be−Cr, where parameters A, B, and C were determined by fitting the respective electronic energies calculated at r = 2.2, 2.3, and 2.4 Å. For the outermost, 30.0 to 100.0 Å, PEC segments, we adopted the appropriate long-range forms of the X1Σg+ and A1Σu+ potentials, which areV(r)=De−Σm=02C6+2m/r6+2m(4)in the X1Σg+ case andV(r)=De−Σm=05C3+m/r3+m(5)in the case of the A1Σu+ PEC. The C6, C8, and C10 coefficients entering the former formula were obtained by fitting the X1Σg+ electronic energies calculated at r = 25.0, 30.0, and 100.0 Å to Eq. 4, in which De was defined as the relevant energy difference between r = 100.0 Å and re, with re representing the previously determined equilibrium internuclear separation in the ground electronic state. The six coefficients C3 through C8 entering the latter expression were obtained by fitting the A1Σu+ electronic energies calculated at r = 13.0, 15.0, 20.0, 25.0, 30.0, and 100.0 Å to Eq. 5, in which De was set as the energy difference between r = 100.0 Å and the corresponding re.

The quality of the potential fits generated by LEVEL16 is very high. We illustrate it here by summarizing the results of two of the several numerical tests that we carried out for the ground-state PEC. In one of the tests, we computed the electronic energy of the X1Σg+ state at the internuclear distance r = 3.893 Å, which is the equilibrium bond length determined by the potential fit V(r) produced by LEVEL16, using our ab initio quantum chemistry protocol defined by Eq. 1. The resulting energy, determined relative to the asymptotic value of the X1Σg+ potential corresponding to r = 100.0 Å, matched the value of V(r) at r = 3.893 Å obtained with LEVEL16 to within 0.0001 cm−1. In another test, aimed at examining the ability of the interpolation scheme used by LEVEL16 to reproduce the ab initio energetics obtained using Eq. 1, we removed the electronic energies calculated at r = 4.6, 5.2, 5.8, and 6.4 Å, which is the region of the X1Σg+ PEC where V(r) changes its curvature, and regenerated the potential fit using the remaining ab initio points. The new potential fit, based on fewer ab initio energies, reproduced the equilibrium bond length re = 3.893 Å resulting from the original fit, constructed using all values of r in our grid, to within 0.001 Å (which is the step for numerical integration in LEVEL16). The mean unsigned error characterizing the X1Σg+ energies at the removed points r = 4.6, 5.2, 5.8, and 6.4 Å resulting from the new fit relative to their ab initio values obtained using Eq. 1 was 0.05 cm−1. On the basis of these and similar analyses, including the X1Σg+ and A1Σu+ PECs, we can safely conclude that the potential fits generated by LEVEL16 faithfully represent our ab initio data.

We also used LEVEL16 to determine the rovibrational term values characterizing the experimentally derived analytical X-representation potential developed in (20), which we used to assess the accuracy of our ab initio–determined X1Σg+ PEC in Table 2. To be consistent with our LEVEL16 calculations for the ground-state PEC resulting from the ab initio protocol based on Eq. 1, we first determined the energies corresponding to the X-representation potential on the grid of 47 internuclear distances r adopted in our ab initio work. We then followed the same numerical procedure as described above for the X1Σg+ PEC resulting from the ab initio quantum chemistry calculations.

Last, but not least, we used LEVEL16 to compute the line positions of all allowed A1Σu+(v′, J′) → X1Σg+(v″, J″) rovibronic transitions and, with the help of our ab initio transition dipole moment function μz(X−A)(r), the corresponding line intensities, as defined by the Einstein coefficients. The only adjustment that we had to make to be able to compare our calculated line positions and intensities for the allowed A1Σu+(v′, J′) → X1Σg+(v″, J″) transitions with the LIF data reported in (20, 21) was a uniform downward shift of the entire A1Σu+ PEC resulting from our ab initio computations by 1,543.2 cm−1, needed to match the experimentally determined adiabatic electronic gap Te of 26,068.9 cm−1 (21). Other than that, all of the calculated spectroscopic properties, including the De, re, and rovibrational term values corresponding to the X1Σg+ and A1Σu+ states and the line positions and intensities characterizing the A1Σu+(v′, J′) → X1Σg+(v″, J″) transitions reported in this study, rely on the raw ab initio data compiled in tables S3 to S5 and data files S1 and S2, combined with the LEVEL16 processing, as described above. To produce Fig. 3, we superimposed our theoretical A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″, J″ = 10,12) LIF spectrum on top of the experimental one reported in figure 3 of (20). The theoretical line intensities shown in Fig. 3 were normalized such that the tallest peaks in the calculated and experimental LIF spectra corresponding to the v″ = 5 P12 line representing the A1Σu+(v′ = 3, J′ = 11) → X1Σg+(v″ = 5, J″ = 12) transition match.

SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/6/14/eaay4058/DC1

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.

REFERENCES AND NOTES

  1. ↵
    1. J.-C. Zhang,
    2. J.-F. Sun,
    3. Z.-L. Zhu,
    4. Y.-F. Liu
    , Elastic collision properties of ultracold calcium atoms. Phys. Scr. 87, 025302 (2013).
    OpenUrl
  2. ↵
    1. E. Krotscheck,
    2. R. E. Zillich
    , Solvation of Mg in helium-4: Are there meta-stable Mg dimers? J. Chem. Phys. 145, 244317 (2016).
    OpenUrl
  3. ↵
    1. L. Rybak,
    2. S. Amaran,
    3. L. Levin,
    4. M. Tomza,
    5. R. Moszynski,
    6. R. Kosloff,
    7. C. P. Koch,
    8. Z. Amitay
    , Generating molecular rovibrational coherence by two-photon femtosecond photoassociation of thermally hot atoms. Phys. Rev. Lett. 107, 273001 (2011).
    OpenUrlPubMed
  4. ↵
    1. P. Wolf,
    2. P. Lemonde,
    3. A. Lambrecht,
    4. S. Bize,
    5. A. Landragin,
    6. A. Clairon
    , From optical lattice clocks to the measurement of forces in the Casimir regime. Phys. Rev. A 75, 063608 (2007).
    OpenUrl
  5. ↵
    1. C. Simon,
    2. D. Jaksch
    , Possibility of observing energy decoherence due to quantum gravity. Phys. Rev. A 70, 052104 (2004).
    OpenUrl
  6. ↵
    1. M. Machholm,
    2. P. S. Julienne,
    3. K.-A. Suominen
    , Calculations of collisions between cold alkaline-earth-metal atoms in a weak laser field. Phys. Rev. A 64, 033425 (2001).
    OpenUrl
  7. ↵
    1. P. F. Bernath
    , Extracting potentials from spectra. Science 324, 1526–1527 (2009).
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. K. Patkowski,
    2. V. Špirko,
    3. K. Szalewicz
    , On the elusive twelfth vibrational state of beryllium dimer. Science 326, 1382–1384 (2009).
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. I. Magoulas,
    2. N. P. Bauman,
    3. J. Shen,
    4. P. Piecuch
    , Application of the CC(P;Q) hierarchy of coupled-cluster methods to the beryllium dimer. J. Phys. Chem. A 122, 1350–1368 (2018).
    OpenUrl
  10. ↵
    1. V. V. Meshkov,
    2. A. V. Stolyarov,
    3. M. C. Heaven,
    4. C. Haugen,
    5. R. J. Le Roy
    , Direct-potential-fit analyses yield improved empirical potentials for the ground X1Σg+ state of Be2. J. Chem. Phys. 140, 064315 (2014).
    OpenUrl
  11. ↵
    1. J. M. Merritt,
    2. V. E. Bondybey,
    3. M. C. Heaven
    , Beryllium dimer–caught in the act of bonding. Science 324, 1548–1551 (2009).
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. W. J. Balfour,
    2. A. E. Douglas
    , Absorption spectrum of the Mg2 molecule. Can. J. Phys. 48, 901–914 (1970).
    OpenUrl
  13. ↵
    1. R. Rydberg
    , Graphische Darstellung einiger bandenspektroskopischer Ergebnisse. Z. Phys. 73, 376–385 (1932).
    OpenUrlCrossRef
    1. R. Rydberg
    , Über einige Potentialkurven des Quecksilberhydrids. Z. Phys. 80, 514–524 (1933).
    OpenUrl
    1. O. Klein
    , Zur Berechnung von Potentialkurven für zweiatomige Moleküle mit Hilfe von Spektraltermen. Z. Phys. 76, 226–235 (1932).
    OpenUrl
  16. ↵
    1. A. L. G. Rees
    , The calculation of potential-energy curves from band-spectroscopic data. Proc. Phys. Soc. 59, 998–1008 (1947).
    OpenUrl
  17. ↵
    1. K. C. Li,
    2. W. C. Stwalley
    , Vibrational levels near dissociation in Mg2 and long-range forces. J. Chem. Phys. 59, 4423–4427 (1973).
    OpenUrl
  18. ↵
    1. W. C. Stwalley
    , Long-range analysis of the internuclear potential of Mg2. Chem. Phys. Lett. 7, 600–602 (1970).
    OpenUrl
  19. ↵
    1. W. C. Stwalley
    , Polarizability and long-range interactions of magnesium atoms. J. Chem. Phys. 54, 4517–4518 (1971).
    OpenUrl
  20. ↵
    1. H. Knöckel,
    2. S. Rühmann,
    3. E. Tiemann
    , The X1Σg+ ground state of Mg2 studied by Fourier-transform spectroscopy. J. Chem. Phys. 138, 094303 (2013); 138, 189901 (2013) [Erratum].
    OpenUrlPubMed
  21. ↵
    1. H. Knöckel,
    2. S. Rühmann,
    3. E. Tiemann
    , The A1Σu++(1)1Πu system of Mg2. Eur. Phys. J. D 68, 293 (2014).
    OpenUrl
  22. ↵
    1. H. Scheingraber,
    2. C. R. Vidal
    , Discrete and continuous Franck Condon factors of the Mg2 A1Σu+−X1Σg+ system and their J dependence. J. Chem. Phys. 66, 3694–3704 (1977).
    OpenUrl
  23. ↵
    1. S. G. Porsev,
    2. A. Derevianko
    , High-accuracy relativistic many-body calculations of van der Waals coefficients C6 for alkaline-earth-metal atoms. Phys. Rev. A 65, 020701(R) (2002).
    OpenUrl
  24. ↵
    1. S. G. Porsev,
    2. A. Derevianko
    , High-accuracy calculations of dipole, quadrupole, and octupole electric dynamic polarizabilities and van der Waals coefficients C6, C8, and C10 for alkaline-earth dimers. J. Exp. Theor. Phys. 102, 195–205 (2006).
    OpenUrl
  25. ↵
    1. W. J. Balfour,
    2. R. F. Whitlock
    , Rotational dependence of Franck–Condon factors in the A1Σu+←X1Σg+ system of 24Mg2. Can. J. Phys. 50, 1648–1651 (1972).
    OpenUrl
  26. ↵
    1. S. H. Yuwono,
    2. I. Magoulas,
    3. J. Shen,
    4. P. Piecuch
    , Application of the coupled-cluster CC(P;Q) approaches to the magnesium dimer. Mol. Phys. 117, 1486–1506 (2019).
    OpenUrl
  27. ↵
    1. E. Tiesinga,
    2. S. Kotochigova,
    3. P. S. Julienne
    , Scattering length of the ground-state Mg + Mg collision. Phys. Rev. A 65, 042722 (2002).
    OpenUrl
  28. ↵
    1. J. Čížek
    , 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. 45, 4256–4266 (1966).
    OpenUrlCrossRefWeb of Science
  29. ↵
    1. S. Amaran,
    2. R. Kosloff,
    3. M. Tomza,
    4. W. Skomorowski,
    5. F. Pawłowski,
    6. R. Moszynski,
    7. L. Rybak,
    8. L. Levin,
    9. Z. Amitay,
    10. J. M. Berglund,
    11. D. M. Reich,
    12. C. P. Koch
    , Femtosecond two-photon photoassociation of hot magnesium atoms: A quantum dynamical study using thermal random phase wavefunctions. J. Chem. Phys. 139, 164124 (2013).
    OpenUrl
  30. ↵
    1. K. Patkowski,
    2. R. Podeszwa,
    3. K. Szalewicz
    , Interactions in diatomic dimers involving closed-shell metals. J. Phys. Chem. A 111, 12822–12838 (2007).
    OpenUrl
  31. ↵
    1. H. Li,
    2. H. Feng,
    3. W. Sun,
    4. Y. Zhang,
    5. Q. Fan,
    6. K. A. Peterson,
    7. Y. Xie,
    8. H. F. Schaefer III.
    , The alkaline earth dimer cations (Be2+, Mg2+, Ca2+, Sr2+, and Ba2+). Coupled cluster and full configuration interaction studies. Mol. Phys. 111, 2292–2298 (2013).
    OpenUrl
  32. ↵
    1. K. Raghavachari,
    2. G. W. Trucks,
    3. J. A. Pople,
    4. M. Head-Gordon
    , A fifth-order perturbation comparison of electron correlation theories. Chem. Phys. Lett. 157, 479–483 (1989).
    OpenUrlCrossRefPubMedWeb of Science
  33. ↵
    1. H. J. Monkhorst
    , Calculation of properties with the coupled-cluster method. Int. J. Quantum Chem., Symp. 11, 421–432 (1977).
    OpenUrl
  34. ↵
    1. J. Noga,
    2. R. J. Bartlett
    , The full CCSDT model for molecular electronic structure. J. Chem. Phys. 86, 7041–7050 (1987); 89, 3401 (1988) [Erratum].
    OpenUrl
  35. ↵
    1. G. E. Scuseria,
    2. H. F. Schaefer III.
    , A new implementation of the full CCSDT model for molecular electronic structure. Chem. Phys. Lett. 152, 382–386 (1988).
    OpenUrl
  36. ↵
    1. J. F. Stanton,
    2. R. J. Bartlett
    , The equation of motion coupled-cluster method. A systematic biorthogonal approach to molecular excitation energies, transition probabilities, and excited state properties. J. Chem. Phys. 98, 7029–7039 (1993).
    OpenUrlCrossRefWeb of Science
  37. ↵
    1. K. Kowalski,
    2. P. Piecuch
    , The active-space equation-of-motion coupled-cluster methods for excited electronic states: Full EOMCCSDt. J. Chem. Phys. 115, 643–651 (2001).
    OpenUrl
  38. ↵
    1. S. A. Kucharski,
    2. M. Włoch,
    3. M. Musiał,
    4. R. J. Bartlett
    , Coupled-cluster theory for excited electronic states: The full equation-of-motion coupled-cluster single, double, and triple excitation method. J. Chem. Phys. 115, 8263–8266 (2001).
    OpenUrl
  39. ↵
    1. K. Kowalski,
    2. P. Piecuch
    , New coupled-cluster methods with singles, doubles, and noniterative triples for high accuracy calculations of excited electronic states. J. Chem. Phys. 120, 1715–1738 (2004).
    OpenUrlPubMed
  40. ↵
    1. B. P. Prascher,
    2. D. E. Woon,
    3. K. A. Peterson,
    4. T. H. Dunning Jr..,
    5. A. K. Wilson
    , Gaussian basis sets for use in correlated molecular calculations. VII. Valence, core-valence, and scalar relativistic basis sets for Li, Be, Na, and Mg. Theor. Chem. Acc. 128, 69–82 (2011).
    OpenUrlCrossRef
  41. ↵
    K. A. Peterson, Correlation consistent basis sets; http://tyr0.chem.wsu.edu/~kipeters/basis.html.
  42. ↵
    1. T. Helgaker,
    2. W. Klopper,
    3. H. Koch,
    4. J. Noga
    , Basis-set convergence of correlated calculations on water. J. Chem. Phys. 106, 9639–9646 (1997).
    OpenUrlCrossRefWeb of Science
  43. ↵
    1. P. Piecuch,
    2. J. A. Hansen,
    3. A. O. Ajala
    , Benchmarking the completely renormalised equation-of-motion coupled-cluster approaches for vertical excitation energies. Mol. Phys. 113, 3085–3127 (2015).
    OpenUrl
  44. ↵
    1. C. R. Vidal,
    2. H. Scheingraber
    , Determination of diatomic molecular constants using an inverted perturbation approach: Application to the A1Σu+−X1Σg+ system of Mg2. J. Mol. Spectrosc. 65, 46–64 (1977).
    OpenUrl
  45. ↵
    1. M. Douglas,
    2. N. M. Kroll
    , Quantum electrodynamical corrections to the fine structure of helium. Ann. Phys. 82, 89–155 (1974).
    OpenUrlCrossRefWeb of Science
  46. ↵
    1. T. Nakajima,
    2. K. Hirao
    , The Douglas–Kroll–Hess approach. Chem. Rev. 112, 385–402 (2012).
    OpenUrlPubMed
  47. ↵
    M. S. Gordon, M. W. Schmidt, Advances in electronic structure theory: GAMESS a decade later. Eds., C. E. Dykstra, G. Frenking, K. S. Kim, G. E. Scuseria, in Theory and Applications of Computational Chemistry: The First Forty Years (Elsevier, 2005), pp. 1167–1189.
  48. ↵
    1. M. Valiev,
    2. E. J. Bylaska,
    3. N. Govind,
    4. K. Kowalski,
    5. T. P. Straatsma,
    6. H. J. J. Van Dam,
    7. D. Wang,
    8. J. Nieplocha,
    9. E. Apra,
    10. T. L. Windus,
    11. W. A. de Jong
    , NWChem: A comprehensive and scalable open-source solution for large scale molecular simulations. Comput. Phys. Commun. 181, 1477–1489 (2010).
    OpenUrlCrossRefWeb of Science
  49. ↵
    1. P. Piecuch,
    2. S. A. Kucharski,
    3. K. Kowalski,
    4. M. Musiał
    , Efficient computer implementation of the renormalized coupled-cluster methods: The R-CCSD[T], R-CCSD(T), CR-CCSD[T], and CR-CCSD(T) approaches. Comp. Phys. Commun. 149, 71–96 (2002).
    OpenUrlCrossRef
  50. ↵
    1. J. W. Cooley
    , An improved eigenvalue corrector formula for solving the Schrödinger equation for central fields. Math. Comput. 15, 363–374 (1961).
    OpenUrl
  51. ↵
    1. R. J. Le Roy
    , LEVEL: A computer program for solving the radial Schrödinger equation for bound and quasibound levels. J. Quant. Spectrosc. Radiat. Transf. 186, 167–178 (2017).
    OpenUrl
  52. ↵
    1. R. J. Le Roy,
    2. R. B. Bernstein
    , Shape resonances and rotationally predissociating levels: The atomic collision time-delay functions and quasibound level properties of H2(X1Σg+). J. Chem. Phys. 54, 5114–5126 (1971).
    OpenUrlCrossRefWeb of Science
Acknowledgments: Funding: This work has been supported by the Chemical Sciences, Geosciences and Biosciences Division, Office of Basic Energy Sciences, Office of Science, U.S. Department of Energy (grant no. DE-FG02-01ER15228) and National Science Foundation (grant no. CHE-1763371). Author contributions: I.M. and P.P. conceived the project. S.H.Y. and I.M. performed the electronic structure calculations under P.P.’s supervision. S.H.Y. carried out the rovibrational analysis and computed the rovibronic spectrum, which was subsequently used to simulate the relevant LIF spectra. All authors contributed equally and extensively to the discussion and interpretation of the data. All authors worked closely together on preparing the manuscript. P.P. supervised the entire project. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions of this study are included in the paper and/or the Supplementary Materials. Any additional information related to this paper may be requested from the authors.
View this article with LENS

Article Tools

My saved folders

Stay Connected to Science Advances

Related Content

Similar Articles in:

Citing Articles in:

Navigate This Article