Discovery of two new magnesium-bearing species in IRC+10216: MgC3N and MgC4H

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.

Spectroscopic Study of [Mg, H, N, C, O] Species: Implications for the Astrochemical Magnesium Chemistry

Magnesium is an abundant metal element in space, and magnesium chemistry has vital importance in the evolution of interstellar medium (ISM) and circumstellar regions, such as the asymptotic giant branch star IRC+10216 where a variety of Mg compounds bearing H, C, N, and O have been detected and proposed as the important components in the gas-phase molecular clouds and solid-state dust grains. Herein, we report the formation and infrared spectroscopic characterization of the Mg-bearing molecules HMg, [Mg, N, C], [Mg, H, N, C], [Mg, N, C, O], and [Mg, H, N, C, O] from the reactions of Mg/Mg+ and the prebiotic isocyanic acid (HNCO) in the solid neon matrix. Based on their thermal diffusion and photochemical behavior, a complex reactivity landscape involving association, decomposition, and isomerization reactions of these Mg-bearing molecules is developed, which can not only help understand the chemical processes of the magnesium (iso)cyanides in astrochemistry but also provide implications on the presence of magnesium (iso)cyanates in the ISM and the chemical model for the dust grain surface reactions. It also provides a new paradigm of the key intermediate nature of the cationic complexes in the formation of neutral interstellar species.

The role of small molecular cations in the chemical flow of the interstellar environments

Molecular ions have been ubiquitous in a variety of environments in the interstellar medium, from Circumstellar Envelopes to Dark Molecular Clouds and to Diffuse Clouds. Their role in the multitude of molecular processes which have been found to occur in those environments has been the subject of many studies over the years, so that we have acquired by now a complex body of data on their chemical structures, their possible function within chemical reactions and their most likely paths to formation. In the present work we review a broad range of such molecular ions, focusing exclusively on positive ions involving the smallest and simplest cations which have been either detected or conjectured as present in the interstellar medium (ISM). We therefore consider mainly molecular cations formed with components like H, H+, He and He+, atomic species which are by far the most abundant baryons in the ISM in general. Their likely structures and their roles in a variety of chemical energy flow paths are discussed and presented within the context of their interstellar environments.

Icy Grain Mantle Surface Astrochemistry of MgNC: The Emergence of Metal Ion Catalysis Studied via Model Ice Cluster Calculations

One of a small number of known magnesium-containing astromolecules, magnesium isocyanide (MgNC) was first detected in 1986. MgNC is an intriguing reactant to consider: it is an open-shell radical in which its metal atom forms a bond with CN that is a mixture of ionic and covalent character. While its gas phase astrochemistry has received prior attention, the grain surface chemistry of MgNC has never been studied. Because of its ionic character, MgNC is found to interact far more strongly with an ice surface than molecules with a greater degree of covalency. As a radical, it may react with closed-shell molecules deposited from the gas phase. In this work, cluster calculations treated with density functional theory and correlation consistent basis sets were used to model the deposition of MgNC on clusters containing 17 and 24 water molecules, which were then allowed to react with acetylene (HCCH) and hydrogen cyanide (HCN) as well as with H atoms. The addition of H to MgNC-nH₂O yields hydromagnesium isocyanide (HMgNC), a known astromolecule that may be ejected into the gas phase. HCCH and HCN bind to MgNC-nH₂O to form intermediate radical compounds that may then also react with H atoms. There is enough reaction energy from H addition to eject fragments of the intermediates into the gas phase: the vinyl radical (C₂H₃) for HCCH and the methaniminyl radical (H₂CN) for HCN. That leaves MgNC-nH₂O to perform further catalytic activity. Alternatively, various hydrogenated divalent Mg compounds may also be stabilized and frozen into the ice or potentially ejected into the gas phase. Benchmark coupled cluster theory calculations in limited systems were used to characterize the submerged reaction barriers present when HCCH or HCN add to MgNC in the gas phase.

Why Are MgC3H Isomers Missing in the Interstellar Medium?

Considering the recent findings of linear doublet (²Σ⁺) MgC_nH isomers (n = 2, 4, and 6) in the evolved carbon star IRC+10216, various structural isomers of MgC₃H and MgC₃H⁺ are theoretically investigated here. For MgC₃H, 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 MgC₃H, the global minimum geometry turns out to be a cyclic isomer, 2-magnesabicyclo[1.1.0]but-1,3,4-triyl (1, C_2v, ²A₁). In addition, five low-lying isomers, magnesium-substituted cyclopropenylidene (2, C_s, ²A'), 1-magnesabut-2,3-dien-1-yl-4-ylidene (3, C_s, ²A″), 1-magnesabut-2-yn-1-yl-4-ylidene (4, C_s, ²A″), 2λ³-magnesabicyclo[1.1.0]but-1,3-diyl-4-ylidene (5, C_2v;, ²A₁), and 1-magnesabut-2,3-dien-2-yl-4-ylidene (6, C_∞v, ²Σ⁺), were also identified. The doublet linear isomer of MgC₃H, 1-magnesabutatrienyl (10, C_∞v, ²Σ⁺) turns out to be a minimum but lies 54.1 kcal mol^-1 above 1 at the ROCCSD(T)/cc-pVTZ level. The quartet (⁴Σ⁺) 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 MgC₃H contains more cyclic isomers (1, 2, and 3) on the low-energy side unlike their even-numbered MgC_nH counterparts (n = 2, 4, and 6). Though the quartet electronic state of 10 is linear, it is not the global minimum geometry on the MgC₃H 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, MgC₃H⁺, 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 MgC₃H^0/+ isomers.

Theoretical Rovibrational Spectroscopy of Magnesium Tricarbide-Multireference Character Thwarts a Full Analysis of All Isomers

Magnesium tricarbide isomers are studied herein with coupled cluster theory and multireference configuration interaction to support their possible detection in astrochemical environments such as the circumstellar envelope surrounding the star IRC +10216 or in terrestrial laboratories. Magnesium-bearing species may abound in the interstellar medium (ISM), but only eight (MgNC, MgCN, HMgNC, MgC₂H, MgC₃N, MgC₄H, MgC₅N, and MgC₆H) have been directly identified thus far. Several possible isomers for the related MgC₃ system are explored in their singlet and triplet spin multiplicities. Overall, this work offers quantum chemical insight of rovibrational spectroscopic data for MgC₃ using quartic force fields (QFFs) based on the CCSD(T) and CCSD(T)-F12 levels of theory at the complete basis set (CBS) limit. Additional corrections with small basis set CCSDT(Q) and scalar relativistic effects are also included in the analysis. Salient multireference character is found in the singlet diamond electronic state, which makes a definitive assignment of the ground state challenging. Nevertheless, coupled cluster-based composite energies and multireference configuration interaction both predict that the ¹A₁ diamond isomer is 1.6-2.2 kcal mol^-1 lower in energy than the ³A₁ diamond isomer. Furthermore, highly accurate binding energies of various isomers MgC₃ are provided for comparison to photodetachment experiments. Dipole moments along with harmonic infrared intensities will guide efforts for astronomical and spectroscopic characterization.

Rotational spectroscopy of the large saturated dinitriles hexanedinitrile and heptanedinitrile

Dinitriles with a saturated hydrocarbon skeleton and a -C≡N group at each end can have large electric dipole moments. Their formation can be related to highly reactive radicals such as CH₂CN, C₂N, or CN. Thus, these saturated dinitriles are potential candidates to be observed in the interstellar medium. In this work, two members of this family, hexanedinitrile and heptanedinitrile, have been investigated through their rotational spectra. The jet-cooled broadband chirped-pulse Fourier transform microwave spectra of both molecules were measured in the 2-8 GHz frequency region. Three and six conformers of hexanedinitrile and heptanedinitrile, respectively, were detected and assigned based on the rotational and quadrupole coupling constants.

Theoretical investigation of the structures, stabilities, and vibrational and rotational spectroscopic parameters of linear HOMgNC and HMgNCO molecules by density functional theory and coupled-cluster method

Linear HOMgNC and HMgNCO molecules: two appropriate candidates for interstellar observation and experimental preparation.

TMC-1, the starless core sulfur factory: Discovery of NCS, HCCS, H2CCS, H2CCCS, and C4S and detection of C5S

We report the detection of the sulfur-bearing species NCS, HCCS, H2CCS, H2CCCS, and C4S for the first time in space. These molecules were found towards TMC-1 through the observation of several lines for each species. We also report the detection of C5S for the first time in a cold cloud through the observation of five lines in the 31-50 GHz range. The derived column densities are N(NCS) = (7.8±0.6)×1011 cm-2, N(HCCS) = (6.8±0.6)×1011 cm-2, N(H2CCS) = (7.8±0.8)×1011 cm-2, N(H2CCCS) = (3.7±0.4)×1011 cm-2, N(C4S) = (3.8±0.4)×1010 cm-2, and N(C5S) = (5.0±1.0)×1010 cm-2. The observed abundance ratio between C3S and C4S is 340, that is to say a factor of approximately one hundred larger than the corresponding value for CCS and C3S. The observational results are compared with a state-of-the-art chemical model, which is only partially successful in reproducing the observed abundances. These detections underline the need to improve chemical networks dealing with S-bearing species.

Discovery of allenyl acetylene, H2CCCHCCH, in TMC-1. A study of the isomers of C5H4

We present the discovery in TMC-1 of allenyl acetylene, H2CCCHCCH, through the observation of nineteen lines with a signal-to-noise ratio ~4-15. For this species, we derived a rotational temperature of 7±1K and a column density of 1.2±0.2×1013 cm-2. The other well known isomer of this molecule, methyl diacetylene (CH3C4H), has also been observed and we derived a similar rotational temperature, Tr=7.0±0.3 K, and a column density for its two states (A and E) of 6.5±0.3×1012 cm-2. Hence, allenyl acetylene and methyl diacetylene have a similar abundance. Remarkably, their abundances are close to that of vinyl acetylene (CH2CHCCH). We also searched for the other isomer of C5H4, HCCCH2CCH (1.4-Pentadiyne), but only a3σ upper limit of 2.5×1012 cm-2 to the column density can be established. These results have been compared to state-of-the-art chemical models for TMC-1, indicating the important role of these hydrocarbons in its chemistry. The rotational parameters of allenyl acetylene have been improved by fitting the existing laboratory data together with the frequencies of the transitions observed in TMC-1.

Space and laboratory discovery of HC3 . Astron Astrophys 2021;646:L3. [PMID: 33824540 PMCID: PMC7610522 DOI: 10.1051/0004-6361/202040013]

We report the detection in TMC-1 of the protonated form of C3S. The discovery of the cation HC3S+ was carried through the observation of four harmonically related lines in the Q band using the Yebes 40m radiotelescope, and is supported by accurate ab initio calculations and laboratory measurements of its rotational spectrum. We derive a column density N(HC3S+) = (2.0 ± 0.5) × 1011 cm-2, which translates to an abundance ratio C3S/HC3S+ of 65 ± 20. This ratio is comparable to the CS/HCS+ ratio (35 ± 8) and is a factor of about ten larger than the C3O/HC3O+ ratio previously found in the same source. However, the abundance ratio HC3O+/HC3S+ is 1.0 ± 0.5, while C3O/C3S is just ~ 0.11. We also searched for protonated C2S in TMC-1, based on ab initio calculations of its spectroscopic parameters, and derive a 3σ upper limit of N(HC2S+)≤ 9×1011 cm-2 and a C2S/HC2S+ ≥ 60. The observational results are compared with a state-of-the-art gas-phase chemical model and conclude that HC3S+ is mostly formed through several pathways: proton transfer to C3S, reaction of S+ with c-C3H2, and reaction between neutral atomic sulfur and the ion C3H+ 3.

Discovery of the acetyl cation, CH3CO+, in space and in the laboratory

Using the Yebes 40m and IRAM 30m radiotelescopes, we detected two series of harmonically related lines in space that can be fitted to a symmetric rotor. The lines have been seen towards the cold dense cores TMC-1, L483, L1527, and L1544. High level of theory ab initio calculations indicate that the best possible candidate is the acetyl cation, CH3CO+, which is the most stable product resulting from the protonation of ketene. We have produced this species in the laboratory and observed its rotational transitions Ju = 10 up to Ju = 27. Hence, we report the discovery of CH3CO+ in space based on our observations, theoretical calculations, and laboratory experiments. The derived rotational and distortion constants allow us to predict the spectrum of CH3CO+ with high accuracy up to 500 GHz. We derive an abundance ratio N(H2CCO)/N(CH3CO+)~44. The high abundance of the protonated form of H2CCO is due to the high proton affinity of the neutral species. The other isomer, H2CCOH+, is found to be 178.9 kJ mol-1 above CH3CO+. The observed intensity ratio between the K=0 and K=1 lines, ~2.2, strongly suggests that the A and E symmetry states have suffered interconversion processes due to collisions with H and/or H2, or during their formation through the reaction of H 3 + with H2CCO.

Yebes 40 m radio telescope and the broad band NANOCOSMOS receivers at 7 mm and 3 mm for line surveys

CONTEXT

Yebes 40m radio telescope is the main and largest observing instrument at Yebes Observatory and it is devoted to Very Long Baseline Interferometry (VLBI) and single dish observations since 2010. It has been covering frequency bands between 2 GHz and 90 GHz in discontinuous and narrow windows in most of the cases, to match the current needs of the European VLBI Network (EVN) and the Global Millimeter VLBI Array (GMVA).

AIMS

Nanocosmos project, a European Union funded synergy grant, opened the possibility to increase the instantaneous frequency coverage to observe many molecular transitions with single tunnings in single dish mode. This reduces the observing time and maximises the output from the telescope.

METHODS

We present the technical specifications of the recently installed 31.5 - 50GHz (Q band) and 72 - 90.5 GHz (W band) receivers along with the main characteristics of the telescope at these frequency ranges. We have observed IRC+10216, CRL 2688 and CRL 618, which harbour a rich molecular chemistry, to demonstrate the capabilities of the new instrumentation for spectral observations in single dish mode.

RESULTS

The results show the high sensitivity of the telescope in the Q band. The spectrum of IRC+10126 offers a signal to noise ratio never seen before for this source in this band. On the other hand, the spectrum normalised by the continuum flux towards CRL 618 in the W band demonstrates that the 40 m radio telescope produces comparable results to those from the IRAM 30 m radio telescope, although with a smaller sensitivity. The new receivers fulfil one of the main goals of Nanocosmos and open the possibility to study the spectrum of different astrophysical media with unprecedented sensitivity.

Discovery of HC4NC in TMC-1: A study of the isomers of HC3N, HC5N, and HC7N

We present a study of the isocyano isomers of the cyanopolyynes HC₃N, HC₅N, and HC₇N in TMC-1 and IRC+10216 carried out with the Yebes 40m radio telescope. This study has enabled us to report the detection, for the first time in space, of HCCCCNC in TMC-1 and to give upper limits for HC₆NC in the same source. In addition, the deuterated isotopologues of HCCNC and HNCCC were detected, along with all ¹³C substitutions of HCCNC, also for the first time in space. The abundance ratios of HC₃N and HC₅N, with their isomers, are very different in TMC-1 and IRC+10216, namely, N(HC₅N)/N(HC₄NC) is ~300 and ≥2100, respectively. We discuss the chemistry of the metastable isomers of cyanopolyynes in terms of the most likely formation pathways and by comparing observational abundance ratios between different sources.

Discovery of HC3O+ in space: The chemistry of O-bearing species in TMC-1

Using the Yebes 40m and IRAM 30m radio telescopes, we detected a series of harmonically related lines with a rotational constant B 0=4460.590±0.001 MHz and a distortion constant D 0=0.511 ±0.005 kHz towards the cold dense core TMC-1. High-level-of-theory ab initio calculations indicate that the best possible candidate is protonated tricarbon monoxide, HC3O+. We have succeeded in producing this species in the laboratory and observed its J u -J l = 2-1 and 3-2 rotational transitions. Hence, we report the discovery of HC3O+ in space based on our observations, theoretical calculations, and laboratory experiments. We derive an abundance ratio N(C3O)/N(HC3O+)~7. The high abundance of the protonated form of C3O is due to the high proton affinity of the neutral species. The chemistry of O-bearing species is modelled, and predictions are compared to the derived abundances from our data for the most prominent O-bearing species in TMC-1.

MgC6H2 Isomers: Potential Candidates for Laboratory and Radioastronomical Studies

Eighty three stationary points of MgC₆H₂ isomers spanning from 0 to 215 kcal mol^-1 have been theoretically identified using density functional theory at the B3LYP/6-311++G(2d,2p) level of theory. Among them, four low-lying isomers lying within 23.06 kcal mol^-1 (1 eV) have been further characterized in detail using high-level coupled-cluster (CC) methods. The thermodynamically most stable isomer turns out to be 1-magnesacyclohepta-4-en-2,6-diyne (1). The other three isomers, 3-magnesahepta-1,4,6-triyne (2), 1-magnesacyclohepta-2,3,4-trien-6-yne (3), and 1-magnesahepta-2,4,6-triyne (4) lie 8.24, 19.76, and 21.36 kcal mol^-1, respectively, above 1 at the ae-CCSD(T)/cc-pCVTZ level of theory. All the four isomers are polar with a permanent electric dipole moment (μ ≠ 0). Hence, they are potential candidates for rotational spectroscopic studies. Considering the recent identification of magnesium-bearing hydrocarbons such as, MgC₂H and MgC₄H in IRC+10216, it is believed that the current theoretical data may be of relevance to laboratory molecular spectroscopic and radioastronomical studies on MgC₆H₂ isomers. The energetic and spectroscopic information gathered in this study would aid the detection of low-lying MgC₆H₂ isomers in the laboratory, which are indispensable for radioastronomical studies. It is also noted here that neither the National Institute of Standards and Technology Chemistry WebBook nor the Kinetic Database for Astrochemistry lists any isomer of MgC₆H₂ at the moment. Therefore, these isomers are studied here theoretically for the very first time.

Interstellar nitrile anions: Detection of C3N- and C5N- in TMC-1. Astron Astrophys 2020;641:L9. [PMID: 33173234 PMCID: PMC7116340 DOI: 10.1051/0004-6361/202039231]

We report on the first detection of C3N- and C5N- towards the cold dark core TMC-1 in the Taurus region, using the Yebes 40 m telescope. The observed C3N/C3N- and C5N/C5N- abundance ratios are ~140 and ~2, respectively; that is similar to those found in the circumstellar envelope of the carbon-rich star IRC +10216. Although the formation mechanisms for the neutrals are different in interstellar (ion-neutral reactions) and circumstellar clouds (photodissociation and radical-neutral reactions), the similarity of the C3N/C3N- and C5N/C5N- abundance ratios strongly suggests a common chemical path for the formation of these anions in interstellar and circumstellar clouds. We discuss the role of radiative electronic attachment, reactions between N atoms and carbon chain anions C n -, and that of H- reactions with HC3N and HC5N as possible routes to form C n N-. The detection of C5N- in TMC-1 gives strong support for assigning to this anion the lines found in IRC +10216, as it excludes the possibility of a metal-bearing species, or a vibrationally excited state. New sets of rotational parameters have been derived from the observed frequencies in TMC-1 and IRC +10216 for C5N- and the neutral radical C5N.

Detection of vibrationally excited HC7N and HC9N in IRC+10216 ⋆

Observations of IRC +10216 with the Yebes 40m telescope between 31 and 50 GHz have revealed more than 150 unidentified lines. Some of them can be grouped into a new series of 26 doublets, harmonically related with integer quantum numbers ranging from J up=54 to 80. The separation of the doublets increases systematically with J, i.e., as expected for a linear species in one of its bending modes. The rotational parameters resulting from the fit to these data are B = 290.8844 ± 0.0004 MHz, D = 0.88 ± 0.04 Hz, q = 0.1463 ± 0.0001 MHz. The rotational constant is very close to that of the ground state of HC9N. Ab initio calculations show an excellent agreement between these parameters and those predicted for the lowest energy vibrationally excited state, ν 19=1, of HC9N. This is the first detection, and complete characterization in space, of vibrationally excited HC9N. An energy of 41.5 cm-1 is estimated for the ν 19 state. In addition, 17 doublets of HC7N in the ν 15=1 state, for which laboratory spectroscopy is available, have been detected for the first time in IRC+10216. Several doublets of HC5N in its ν 11=1 state have been also observed. The column density ratio between the ground and the lowest excited vibrational states are ≈127, 9.5, and 1.5 for HC5N, HC7N, and HC9N, respectively. We find that these lowest-lying vibrational states are most probably populated via infrared pumping to vibrationally excited states lying at ≈600 cm-1. The lowest vibrationally excited states thus need to be taken into account to precisely determine absolute abundances and abundanceratios for long carbon chains. The abundance ratios N(HC5N)/N(HC7N) and N(HC7N)/N(HC9N) are 2.4 and 7.7 respectively.