Energetic and Spectroscopic Properties of Astrophysically Relevant MgC4H Radicals Using High-Level Ab Initio Calculations.

Considering the importance of magnesium-bearing hydrocarbon molecules (MgCnH; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC4H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2Σ+), is the global minimum geometry on the MgC4H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ∼ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ3-magnesacycloprop-2-ene-2,3-diyl (2, 2A1, µ = 2.54 D), is 146 kJ mol-1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

Why an integrated approach between search algorithms and chemical intuition is necessary?

Though search algorithms are appropriate tools for identifying low-energy isomers, fixing several constraints seems to be a fundamental prerequisite to successfully running any structural search program. This causes some potential setbacks as far as identifying all possible isomers, close to the lowest-energy isomer, for any elemental composition. The number of explored candidates, the choice of method, basis set, and availability of CPU time needed to analyze the various initial test structures become necessary restrictions in resolving the issues of structural isomerism reasonably. While one could arrive at new structures through chemical intuition, reproducing or achieving those exact same structures requires increasing the number of variables in any given program, which causes further constraints in exploring the potential energy surface in a reasonable amount of time. Thus, it is emphasized here that an integrated approach between search algorithms and chemical intuition is necessary by taking the C12O2Mg2 system as an example. Our initial search through the AUTOMATON program yielded 1450 different geometries. However, through chemical intuition, we found eighteen new geometries within 40.0 kcal mol-1 at the PBE0-D3/def2-TZVP level. These results indirectly emphasize that an integrated approach between search algorithms and chemical intuition is necessary to further our knowledge in chemical space for any given elemental composition.

Why Are MgC3H Isomers Missing in the Interstellar Medium?

Considering the recent findings of linear doublet (2Σ+) MgCnH isomers (n = 2, 4, and 6) in the evolved carbon star IRC+10216, various structural isomers of MgC3H and MgC3H+ are theoretically investigated here. For MgC3H, 11 doublet and 8 quartet stationary points ranging from 0.0 to 71.8 and 0.0 to 110.1 kcal mol-1, respectively, have been identified initially at the UωB97XD/6-311++G(2d,2p) level. To get accurate relative energies, further energy evaluations are carried out for all isomers with coupled cluster methods and thermochemical modules such as G3//B3LYP, G4MP2, and CBS-QB3 methods. Unlike the even series, where the global minima are linear molecules with a Mg atom at one end, in the case of MgC3H, the global minimum geometry turns out to be a cyclic isomer, 2-magnesabicyclo[1.1.0]but-1,3,4-triyl (1, C2v, 2A1). In addition, five low-lying isomers, magnesium-substituted cyclopropenylidene (2, Cs, 2A'), 1-magnesabut-2,3-dien-1-yl-4-ylidene (3, Cs, 2A″), 1-magnesabut-2-yn-1-yl-4-ylidene (4, Cs, 2A″), 2λ3-magnesabicyclo[1.1.0]but-1,3-diyl-4-ylidene (5, C2v;, 2A1), and 1-magnesabut-2,3-dien-2-yl-4-ylidene (6, C∞v, 2Σ+), were also identified. The doublet linear isomer of MgC3H, 1-magnesabutatrienyl (10, C∞v, 2Σ+) turns out to be a minimum but lies 54.1 kcal mol-1 above 1 at the ROCCSD(T)/cc-pVTZ level. The quartet (4Σ+) electronic state of 10 was also found to be a minimum, but it lies 8.0 kcal mol-1 above 1 at the same level. Among quartets, isomer 10 is the most stable molecule. The next quartet electronic state (of isomer 11) is 34.4 kcal mol-1 above 10, and all other quartet electronic states of other isomers are not energetically close to low-lying doublet isomers 2 to 6. Overall, the chemical space of MgC3H contains more cyclic isomers (1, 2, and 3) on the low-energy side unlike their even-numbered MgCnH counterparts (n = 2, 4, and 6). Though the quartet electronic state of 10 is linear, it is not the global minimum geometry on the MgC3H potential energy surface. Isomerization pathways among the low-lying isomers (doublets of 1-4 and a quartet of 10) reveal that these molecules are kinetically stable. For the cation, MgC3H+, the cyclic isomers (1+, 2+, and 3+) are on the low-energy side. The singlet linear isomer, 10+, is a fourth-order saddle point. The low-lying cations are quite polar, with dipole moment values of >7.00 D. The current theoretical data would be helpful to both laboratory astrophysicists and radioastronomers for further studies on the MgC3H0/+ isomers.