Free Access
Volume 652, August 2021
Article Number L9
Number of page(s) 7
Section Letters to the Editor
Published online 18 August 2021

© ESO 2021

1. Introduction

The recent discovery with the QUIJOTE1 line survey of abundant hydrocarbons such as the propargyl radical (CH2CCH), vinylacetylene (CH2CHCCH), and ethynylallene (H2CCCHCCH) and cyclic hydrocarbons such as cyclopentadiene (c-C5H6) and indene (c-C9H8) in the cold dark cloud TMC-1 (Agúndez et al. 2021a; Cernicharo et al. 2021a,b,c; Burkhardt et al. 2021), as well as the detection of their cyano derivatives (McGuire et al. 2018, 2021) in the same source, points to a new and interesting chemistry that had been completely unforeseen by the most sophisticated state-of-the-art chemical models. The detection of indene, the first non-functionalised polycyclic aromatic hydrocarbons (PAHs) detected in space, in the unexpected environment of a cold starless core is surprising and introduces new challenges to our understanding of the place were these PAHs are formed in space.

Many efforts have been devoted over the last 40 years to understanding the chemical processes that lead to the formation of these molecular species (see e.g. Jones et al. 2011; Parker et al. 2012; Joblin & Cernicharo 2018). Circumstellar envelopes around carbon-rich asymptotic giant branch (AGB) stars have been suggested as the factories of PAHs (Cherchneff et al. 1992). However, PAHs are only detected in well-UV-illuminated regions of the interstellar medium. The detection of benzene in the carbon-rich protoplanetary nebula CRL 618 (Cernicharo et al. 2001) suggests a bottom-up approach in which the small hydrocarbons that formed during the AGB phase, such as C2H2 and C2H4, interact with the UV radiation produced by the star in its evolution to the white dwarf phase (Woods et al. 2002; Cernicharo 2004). Other hypotheses involve the processing of dust grains around evolved stars, either through UV photons (Pilleri et al. 2015) or by chemical processes (Martínez et al. 2020). Hence, it has been surprising to see that pure PAHs and their cyanide derivatives have been found in TMC-1, which is well protected against UV radiation (Cernicharo et al. 2021c; Burkhardt et al. 2021; McGuire et al. 2018, 2021). It is unlikely that these molecules arise from a reservoir existing since the early stages of the cloud because these relatively small PAHs would not have survived the diffuse cloud stage. An in situ formation mechanism for benzene, cyclopentadiene, indene, and naphthalene must involve abundant hydrocarbons containing from two to five carbon atoms (Kaiser & Hansen 2021). Moreover, some of these species should allow an easy cyclisation in bimolecular reactions to efficiently form the first aromatic ring – benzene (c-C6H6), benzyne (c-C6H4), the phenyl radical (c-C6H5), or any other species – from which larger PAHs can grow (Jones et al. 2011; Parker et al. 2012; McCabe et al. 2020; Doddipatla et al. 2021). The progargyl radical and the closed-shell hydrocarbons methyl-, vinyl-, and ethynylallene have been detected in TMC-1 with high abundances (Agúndez et al. 2021a; Cernicharo et al. 2021b,c). Detecting all possible products of the reactions between these hydrocarbons and radicals is a mandatory step towards understand the chemistry of PAHs in cold dark clouds.

The name benzyne refers to three isomeric species with molecular formula c-C6H4, which are six-membered rings with four hydrogen atoms. They can be considered as derivatives of benzene molecules due to the loss of two hydrogen atoms. The most stable isomer is 1,2-didehydrobenzene (ortho-benzyne), followed by 1,3-didehydrobenzene (meta-benzyne) and the 1,4-didehydrobenzene (para-benzyne), placed at 147 and 220 kJ mol−1 from the ortho-benzyne, respectively (Olsen 1971). The molecular structure of ortho-benzyne (hereinafter referred to as benzyne or as o-C6H4) is interesting because it contains a formal triple bond (‘-yne’) within a six-membered ring (see Fig. 1). This peculiar structure confers a high reactivity, and thus it is a reactive intermediate in organic chemistry that plays an important role in aromatic nucleophilic substitution reactions. In addition, it is a potential key species in the formation of PAHs and soot in pyrolysis experiments of organic molecules (Hirsch et al. 2018).

thumbnail Fig. 1.

Scheme of the ortho, para, and meta isomers of benzyne.

In this Letter we report the discovery of the ortho isomer of benzyne, o-C6H4, and we discuss the relative abundances of pure hydrocarbons, chains and cycles, in the context of chemical models that include new radical-neutral reactions involving the abundant hydrocarbons recently found in TMC-1 (Agúndez et al. 2021a; Cernicharo et al. 2021a,b,c). We also discuss the evolution of these molecules when the volume density increases due to the gravitational collapse of the core to form protostars and their dusty planetary disks.

2. Observations

New receivers, built within the Nanocosmos project2 and installed at the Yebes 40 m radio telescope, were used for the observations of TMC-1 (αJ2000 = 4h41m41.9s and δJ2000 = +25° 41′27.0″). A detailed description of the system is given by Tercero et al. (2021). The receiver consists of two cold high electron mobility transistor amplifiers that cover the 31.0−50.3 GHz band with horizontal and vertical polarisations. Receiver temperatures in the runs achieved in 2020 vary from 22 K at 32 GHz to 42 K at 50 GHz. Some power adaptation in the down-conversion chains have reduced the receiver temperatures in 2021 to 16 K at 32 GHz and 25 K at 50 GHz. The back ends are 2 × 8 × 2.5 GHz fast Fourier transform spectrometers with a spectral resolution of 38.15 kHz, providing the whole coverage of the Q band in both polarisations. All observations were performed in the frequency switching mode. The main beam efficiency varies from 0.6 at 32 GHz to 0.43 at 50 GHz. The intensity scale used in this work, antenna temperature (), was calibrated using two absorbers at different temperatures and the atmospheric transmission model ATM (Cernicharo 1985; Pardo et al. 2001). Calibration uncertainties were adopted to be 10%. All data were analysed using the GILDAS package3.

The QUIJOTE line survey was performed in several sessions between 2019 and 2021. Previous QUIJOTE results on the detection of C3N and C5N (Cernicharo et al. 2020a), HC5NH+ (Marcelino et al. 2020), HC4NC (Cernicharo et al. 2020b), and HC3O+ (Cernicharo et al. 2020c) were based on two observing runs performed in November 2019 and February 2020. Two different frequency coverages were achieved during these runs, 31.08−49.52 GHz and 31.98−50.42 GHz, in order to check that no spurious spectral ghosts are produced in the down-conversion chain. In these observations the frequency throw was 10 MHz.

Additional data were taken in October 2020, December 2020, and January−May 2021. In these observations the selected frequency coverage was 31.08−49.52 GHz, and the frequency throw for the frequency switching observations was 8 MHz. The two different frequency throws allow unambiguously negative features produced in the folding of the frequency switching data to be identified. These data allowed the detection of the acetyl cation, CH3CO+ (Cernicharo et al. 2021d), HC3S+ (Cernicharo et al. 2021e), the isomers of C4H3N (Marcelino et al. 2021), l-H2C5 (Cabezas et al. 2021a), vinyl and allenyl acetylene (Cernicharo et al. 2021a,b), ethynyl cyclopropenylidene, cyclopentadiene, and indene (Cernicharo et al. 2021c). Several sulphur-bearing species (NCS, HCCS, H2CCS, H2CCCS, C4S, HCSCN, and HCSCCH) have also been detected with these data (Cernicharo et al. 2021f,g). In addition, some molecules typical of hot cores and corinos (C2H3CHO, C2H3OH, HCOOCH3, CH3OCH3, CH3CH2CN) and species typical of circumstellar envelopes (C5S, HCCN, HC4N) were also detected (Agúndez et al. 2021b; Cernicharo et al. 2021a). Last but not least, a significant number of singly and doubly deuterated species, such as HDCCN (Cabezas et al. 2021b), CH2DC3N (Cabezas et al. 2021c), and the D and 13C doubly substituted isotopologues of HC3N (Tercero et al., in prep.), have also been observed for the first time in space thanks to QUIJOTE.

The total observing time on the source by May 2021 was 240 h. QUIJOTE is a living line survey, and we expect a total observing time on the source of 450 hours by the end of 2021. The final goal of QUIJOTE, a sensitivity of 0.1 mK (lines of 0.5 mK detected at the 5σ level), will require around 1000 hours of observing time and will be reached by the winter of 2023−2024.

The sensitivity of our TMC-1 data is better than that of previously published line surveys of this source at the same frequencies (Kaifu et al. 2004) by a factor 20−30. Below 40 GHz the sensitivity is 0.2 mK, increasing slowly up to 0.5 mK at 49.5 GHz. Thanks to this achievement, it has been possible to detect many individual lines from molecules (Cernicharo et al. 2021a) that have been reported only by stacking techniques (Burkhardt et al. 2021). The level of sensitivity we have reached transforms TMC-1 into a source in which, although still far from the spectral confusion limit, the density of lines below 1 mK in intensity is particularly high. TMC-1 can no longer be considered as a line-poor source. Hence, the identification of new species requires a systematic exploration of the rotational transitions of single and double isotopologues of abundant species (see e.g. Cabezas et al. 2021b,c; Tercero et al., in prep.) before the detection of a new species can be claimed. In this context, line identification in this work was done using the catalogues MADEX (Cernicharo 2012), CDMS (Müller et al. 2005), JPL (Pickett et al. 1998), and Splatalogue4. By May 2021, the MADEX code contained 6369 spectral entries corresponding to the ground and vibrationally excited states, together with the corresponding isotopologues, of 1691 molecules.

3. Results

Our recent discovery of pure hydrocarbon cycles such as cyclopentadiene and indene (Cernicharo et al. 2021c) prompted us to search for other chemically related hydrocarbon cycles, benzyne (o-C6H4) being an obvious search candidate. The benzyne molecule is an asymmetric top with a b dipole directed along the C2v symmetry axis, estimated from ab initio calculations to be 1.38 D (Kraka & Cremer 1993). Laboratory data for this species have been taken from Brown et al. (1986), Kukolich et al. (2003), and Robertson et al. 2003. The rotational levels of benzyne are subjected to spin statistics due to the presence of two pairs of equivalent hydrogen nuclei. As shown in Bunker & Jensen (1998), the nuclear statistical weight for the rotational levels with Ka + Kc even and odd is five and three, respectively (ortho and para species; these names refer to the symmetry states and are unrelated to the same names used to distinguish the different isomers). The first para level, 10, 1, is 0.42 K above the ground ortho state, 000. Ortho-benzyne was previously searched for towards CRL618 without success (Widicus Wwaver et al. 2007).

An examination of the QUIJOTE data quickly indicates that all ortho lines are detected with antenna temperatures between 0.5 and 1.8 mK. A total of seven ortho lines and four para lines are detected. Figure 2 shows the nine lines detected above 3σ. Two of the lines of para benzyne are detected only at the 3σ level and are not shown in this figure (see Table A.1). Nevertheless, their intensities are compatible with our model predictions. Table A.1 provides the line parameters of the observed lines. One of the para lines, the 33, 0 − 22, 1, is fully blended with a line of 2 mK arising from D13CCCN (Tercero et al., in prep.). Three additional lines of the para species are below the 3σ detection limit at the corresponding frequencies (see Table A.1). Some interfering features produce partial blends with three lines of benzyne, as shown in Fig. 2. One of these features arises from a rotational transition of the A species of propylene and overlaps in frequency with the 70, 7 − 61, 6 transition of benzyne. From the derived column density of propylene by Marcelino et al. (2007), the expected contribution from the transition of benzyne is 0.5 mK. A couple of unidentified lines also have some partial blend with the 70, 7 − 61, 6 and 61, 5 − 52, 4 transitions of benzyne. However, the line parameters for these transitions can be derived reasonably well. All remaining lines in the Q band of benzyne are well below our present detection limit. At the present level of sensitivity of QUIJOTE, the number of unidentified lines becomes a concern when any stacking technique is used to detect molecules that produce weak emission. Our detection of benzyne based on the observation of several ortho and para lines is robust, reliable, and provides an accurate estimation of its column density.

thumbnail Fig. 2.

Subset of the observed lines of o-C6H4 in the 31−50 GHz frequency range towards TMC-1 (see Table A.1 for the complete list of detected benzyne lines). The abscissa corresponds to the rest frequency assuming a local standard of rest velocity of 5.83 km s−1. The ordinate is the antenna temperature corrected for atmospheric and telescope losses, in mK. The red line shows the synthetic spectrum obtained from a fit to the observed line profiles, which provides Tr = 7 ± 1 K and N(o-C6H4) = (5.0 ± 1.0) × 1011 cm−2. The rotational quantum numbers are indicated in each panel. Blanked channels correspond to negative features produced in the folding of the frequency switching data. Cyan numbers indicate the multiplicative factor applied to the best model, when needed, to match the observations.

An analysis of the observed intensities using a line profile fitting method (Cernicharo et al. 2021e) provides a rotational temperature of 7 ± 1 K and a total (ortho plus para) column density for benzyne of (5.0 ± 1.0) × 1011 cm−2. In the fitting procedure we considered the ortho and para levels as those of a single molecule with the corresponding statistical weights (5/3). A separate fit considering the two symmetry species separately and fixing the rotational temperature to 7 K provides similar results and an ortho/para ratio of 1.4 ± 0.6, which is consistent with the expected value of 5/3.

Adopting an H2 column density of 1022 cm−2 for TMC-1 (Cernicharo & Guélin 1987), the abundance of o-C6H4 relative to H2 is (5.0 ± 1) × 10−11. This is significantly below those of other cycles, such as c-C3H2, c-C5H6, and c-C9H8 (Cernicharo et al. 2021a).

4. Chemistry

The detection of o-C6H4 in TMC-1 raises interesting questions about its formation. The possible reactions leading to this molecule are discussed in Appendix B. To investigate the overall contribution of reactions (1) and (7) to the formation o-C6H4 in TMC-1, astrochemical models were exploited (see Fig. 3 and Appendix B). We built a pseudo-time-dependent gas-phase chemical model. We adopted typical parameters of cold dense clouds, namely a gas kinetic temperature of 10 K, a volume density of H2 of 2 × 104 cm−3, a cosmic-ray ionisation rate of H2 of 1.3 × 10−17 s−1, a visual extinction of 30 mag, and the set of so-called ‘low-metal’ elemental abundances (see e.g. Agúndez & Wakelam 2013). We note that the volume density of H2 in TMC-1 could be somewhat higher than that in our adopted value. The main effect of a slight increase in the density is that the abundance curves shift to shorter times, while peak abundances may increase or decrease slightly. The gas-phase chemical network is based on the RATE12 network from the UMIST database (McElroy et al. 2013), with updates relevant for the chemistry of CH2CCH, CH2CHCCH, H2CCCHCCH, and c-C5H6 (Agúndez et al. 2021a; Cernicharo et al. 2021a,b,c). Benzyne was included and assumed to be formed by reactions (1) and (7) and destroyed by reactions with C atoms and C+ and H+ ions. The peak fractional abundance of o-C6H4 is 3 × 10−13 relative to H2 (see model A in Fig. 3), which is well below the observed value, (5 ± 1) × 10−11. This suggests that alternative reactions besides (1) and (7) may contribute to o-C6H4.

thumbnail Fig. 3.

Calculated fractional abundance of o-C6H4 as a function of time in models A and B (see text). The horizontal grey band corresponds to the observed abundance in TMC-1.

Among these reactions, (2), (4), (5), (8), (9), (11), and (12) are unlikely to form o-C6H4 (Appendix B). Reactions (3) and (6) are difficult to implement since the radicals C5H3 and C4H5 are not included in the UMIST database. Adding ill-defined radical reactions would add more uncertainty to the output of the models. However, the remaining reactions, (13) and (14) and potentially (10), are plausible routes to o-C6H4. Moreover, these reactions involve reactants that are abundant in TMC-1 (Fossé et al. 2001; Cabezas et al. 2021b; Agúndez et al. 2021a). Including these three reactions results in a peak abundance for o-C6H4 of 3 × 10−11 relative to H2 (model B in Fig. 3), which is in good agreement with the observed value, (5 ± 1) × 10−11. Therefore, reactions (13), (14), and (10) deserve to be studied experimentally and theoretically in the future to understand how this highly reactive cyclic and aromatic molecule, o-C6H4, can be formed in TMC-1. It should be noted that benzyne could also be formed through ionic routes. For example, the dissociative recombination of the c-C6H ion with electrons could produce o-C6H4. The precursor ion, c-C6H, could in turn be formed through different ion-neutral reactions, such as C4H + C2H2 (Anicich 1993; see also Woods et al. 2002) or C3H + CH2CCH.

It would be interesting to investigate whether o-C6H4 could be formed upon the dissociative recombination of C6H or C6 ions with more than five H atoms.

5. Discussion and conclusions

The high abundance of hydrocarbon chains and cycles in TMC-1 suggests that one of the main sites of production of the seeds for PAH formation could be cold dark clouds and that this process could be a transitory one. Figure 3 shows that the peak of formation of o-C6H4 and other hydrocarbons is produced in a few 105 years. When the density increases due to the collapse of the cloud, most of these large hydrocarbons will stick onto dust grains. The chemistry of ices of cyclopentadiene, benzene, indene, and other species is poorly known as most experiments have been done with large PAHs. The dust grains enriched in these molecules, together with their pristine carbon composition produced during their early formation, will evolve with time and will become the constituents of the solid blocks that produce asteroids and planetesimals during the first stages of star and planet formation. Eventually, some of these grains will be exposed to UV radiation from external and/or internal sources producing photodissociation regions. A significant fraction of the PAH emission observed in these regions could be due to PAHs formed from hydrocarbon-rich ices (Abplanalp et al. 2018; Kaiser & Roessler 1997), and ices of cyclopentadiene, benzene, and indene deposited in a transitory phase of between 105 and 106 yr during the evolution of the cloud (we note that chemical times depend on the volume density and adopted initial physical and chemical conditions). None of these species have been found in the deep Q-band line survey of IRC +10216 performed with the Yebes 40 m radio telescope (Pardo et al. 2021). While in the external layers of this carbon-rich circumstellar envelope the abundance of radicals is very high (Agúndez et al. 2017), the chemical time corresponding to these regions is relatively short, ∼104 yr, implying that the first cycles may not have had enough time to be formed. A different context is that of carbon-rich protoplanetary nebulae such as CRL 618, where photodissociation and ionisation is enhanced with respect to the AGB phase and where benzene has been found from infrared observations with the Infrared Space Observatory (Cernicharo et al. 2001). Further, distinct temperature ranges from 10 K to a few thousand K facilitate diverse low- and high-temperature molecular mass growth processes (Cernicharo 2004; Kaiser & Hansen 2021).

To summarise, we report in this work the detection of the ortho isomer of benzyne, o-C6H4, in TMC-1. The observed abundance, (5 ± 1) × 10−11 relative to H2, can be nicely reproduced with fractional abundances of 3 × 10−11 considering that o-C6H4 is formed by the neutral-neutral reactions C2H + CH2CHCCH and C + C5H5 and possibly through the reactions C4H + C2H4, C3H + CH2CCH2, and C3H2 + C3H3. With a sensitive line survey such as QUIJOTE, the discovery of new molecules can be performed in the classical line by line detection method. Hence, QUIJOTE provides unambiguous detection of new molecular species and reliable column density determinations.


Q-band Ultrasensitive Inspection Journey to the Obscure TMC-1 Environment.


We thank ERC for funding through grant ERC-2013-Syg-610256-NANOCOSMOS. M.A. thanks MICIU for grant RyC-2014-16277. We also thank Ministerio de Ciencia e Innovación of Spain (MICIU) for funding support through projects AYA2016-75066-C2-1-P, PID2019-106110GB-I00, PID2019-107115GB-C21/AEI/10.13039/501100011033, and PID2019-106235GB-I00.


  1. Abplanalp, M. J., Jones, B. M., & Kaiser, R. I. 2018, Phys. Chem. Chem. Phys., 20, 5435 [CrossRef] [Google Scholar]
  2. Agúndez, M., & Wakelam, V. 2013, Chem. Rev., 113, 8710 [Google Scholar]
  3. Agúndez, M., Cernicharo, J., Quintana-Lacaci, G., et al. 2017, A&A, 601, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Agúndez, M., Cabezas, C., Tercero, B., et al. 2021a, A&A, 647, L10 [EDP Sciences] [Google Scholar]
  5. Agúndez, M., Marcelino, N., Tercero, B., et al. 2021b, A&A, 649, L4 [CrossRef] [EDP Sciences] [Google Scholar]
  6. Anicich, V. G. 1993, J. Phys. Chem. Ref. Data, 22, 1469 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bunker, P. R., & Jensen, P. 1998, Molecular Symmetry and Spectroscopy, 2nd edn. (Ottawa: NRC Research Press) [Google Scholar]
  8. Burkhardt, A. M., Lee, L. K., Bryan Changala, P., et al. 2021, ApJ, 913, L18 [Google Scholar]
  9. Brown, R. D., Godfrey, P. D., & Rodler, M. 1986, J. Am. Chem. Soc., 108, 1296 [CrossRef] [Google Scholar]
  10. Cabezas, C., Tercero, B., Agúndez, M., et al. 2021a, A&A, 650, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Cabezas, C., Endo, Y., Roueff, E., et al. 2021b, A&A, 646, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cabezas, C., Roueff, E., Tercero, B., et al. 2021c, A&A, 650, L15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Canosa, A., Sims, I. R., Travers, D., et al. 1997, A&A, 323, 644 [Google Scholar]
  14. Cernicharo, J. 1985, Internal IRAM Report (Granada: IRAM) [Google Scholar]
  15. Cernicharo, J. 2004, ApJ, 608, L41 [NASA ADS] [CrossRef] [Google Scholar]
  16. Cernicharo, J. 2012, EAS Publ. Ser., 251 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cernicharo, J., & Guélin, M. 1987, A&A, 176, 299 [Google Scholar]
  18. Cernicharo, J., Heras, A. M., Tielens, A. G. G. M., et al. 2001, ApJ, 546, L123 [NASA ADS] [CrossRef] [Google Scholar]
  19. Cernicharo, J., Marcelino, N., Pardo, J. R., et al. 2020a, A&A, 641, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Cernicharo, J., Marcelino, N., Agúndez, M., et al. 2020b, A&A, 642, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Cernicharo, J., Marcelino, N., Agúndez, M., et al. 2020c, A&A, 642, L17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Cernicharo, J., Agúndez, M., Cabezas, C., et al. 2021a, A&A, 647, L2 [EDP Sciences] [Google Scholar]
  23. Cernicharo, J., Cabezas, C., Agúndez, M., et al. 2021b, A&A, 647, L3 [EDP Sciences] [Google Scholar]
  24. Cernicharo, J., Agúndez, M., Cabezas, C., et al. 2021c, A&A, 649, L15 [EDP Sciences] [Google Scholar]
  25. Cernicharo, J., Cabezas, C., Bailleux, S., et al. 2021d, A&A, 646, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Cernicharo, J., Cabezas, C., Endo, Y., et al. 2021e, A&A, 646, L3 [EDP Sciences] [Google Scholar]
  27. Cernicharo, J., Cabezas, C., Agúndez, M., et al. 2021f, A&A, 648, L3 [EDP Sciences] [Google Scholar]
  28. Cernicharo, J., Cabezas, C., Agúndez, M., et al. 2021g, A&A, 650, L14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Chastaing, D., James, P. L., Sims, I. R., & Smith, I. W. M. 1998, Faraday Discuss., 109, 165 [NASA ADS] [CrossRef] [Google Scholar]
  30. Cherchneff, I., Barker, J. R., & Tielens, A. G. G. M. 1992, ApJ, 401, 269 [NASA ADS] [CrossRef] [Google Scholar]
  31. Doddipatla, S., Galimova, G. R., Wei, H., et al. 2021, Sci. Adv., 7, eabd4044 [Google Scholar]
  32. Fossé, D., Cernicharo, J., Gerin, M., & Cox, P. 2001, ApJ, 552, 168 [Google Scholar]
  33. Hahndorf, I., Lee, Y. T., & Mebel, A. M. 2000, J. Chem. Phys., 113, 9622 [CrossRef] [Google Scholar]
  34. Hahndorf, I., Lee, Y. T., Kaiser, R. I., et al. 2002, J. Chem. Phys., 116, 3248 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hirsch, F., Reusch, E., Constantinidis, P., et al. 2018, J. Phys. Chem. A, 122, 9563 [CrossRef] [Google Scholar]
  36. Joblin, C., & Cernicharo, J. 2018, Science, 359, 156 [CrossRef] [Google Scholar]
  37. Jones, B. M., Zhang, F., Kaiser, R. I., et al. 2011, Proc. Natl. Acad. Sci. USA, 108, 452 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kaifu, N., Ohishi, M., Kawaguchi, K., et al. 2004, PASJ, 56, 69 [Google Scholar]
  39. Kaiser, R. I. 2002, Chem. Rev., 102, 1309 [Google Scholar]
  40. Kaiser, R. I., & Hansen, N. 2021, J. Phys. Chem. A, 125, 3826 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kaiser, R. I., & Roessler, K. 1997, ApJ, 475, 144 [NASA ADS] [CrossRef] [Google Scholar]
  42. Kraka, E., & Cremer, D. 1993, Chem. Phys. Lett., 216, 333 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kukolich, S. G., Tanjaroon, C., McCarthy, M. C., & Thaddeus, P. 2003, J. Chem. Phys., 119, 4353 [NASA ADS] [CrossRef] [Google Scholar]
  44. Marcelino, N., Cernicharo, J., & Agúndez, M. 2007, ApJ, 665, L127 [Google Scholar]
  45. Marcelino, N., Agúndez, M., Tercero, B., et al. 2020, A&A, 643, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Marcelino, N., Tercero, B., Agúndez, M., & Cernicharo, J. 2021, A&A, 646, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Martínez, L., Santoro, G., Merino, P., et al. 2020, Nat. Astron., 4, 97 [NASA ADS] [CrossRef] [Google Scholar]
  48. McCabe, M. N., Hemberger, P., Reusch, E., et al. 2020, J. Phys. Chem. Lett., 11, 2859 [CrossRef] [Google Scholar]
  49. McElroy, D., Walsh, C., Markwick, A. J., et al. 2013, A&A, 550, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. McGuire, B. A., Burkhardt, A. M., Kalenskii, S., et al. 2018, Science, 359, 202 [Google Scholar]
  51. McGuire, B. A., Loomis, R. A., Burkhardt, A. M., et al. 2021, Science, 371, 1265 [Google Scholar]
  52. Mebel, A. M., & Kaiser, R. I. 2015, Int. Rev. Phys. Chem., 17, 21564 [Google Scholar]
  53. Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215 [Google Scholar]
  54. Olsen, J. F. 1971, J. Mol. Struct., 8, 307 [NASA ADS] [CrossRef] [Google Scholar]
  55. Pardo, J. R., Cernicharo, J., & Serabyn, E. 2001, IEEE Trans. Antennas Propag., 49, 12 [Google Scholar]
  56. Pardo, J. R., Cernicharo, J., Tercero, B., et al. 2021, A&A, submitted [Google Scholar]
  57. Parker, D. S. N., Zhang, F., Kim, Y. S., et al. 2012, Proc. Natl. Acad. Sci. USA, 109, 53 [NASA ADS] [CrossRef] [Google Scholar]
  58. Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spectrosc. Radiat. Transf., 60, 883 [Google Scholar]
  59. Pilleri, P., Joblin, C., Boulanger, F., & Onaka, T. 2015, A&A, 577, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Robertson, E. G., Godfrey, P. D., & McNaughton, D. 2003, J. Mol. Spectrosc., 217, 123 [NASA ADS] [CrossRef] [Google Scholar]
  61. Tercero, F., López-Pérez, J. A., Gallego, J. D., et al. 2021, A&A, 645, A37 [EDP Sciences] [Google Scholar]
  62. Widicus Wwaver, S., Remijan, A. J., McMahon, R. J., & McCall, B. J. 2007, ApJ, 671, L153 [NASA ADS] [CrossRef] [Google Scholar]
  63. Woods, P. M., Millar, T. J., & Zijlstra, A. A. 2002, ApJ, 574, L167 [NASA ADS] [CrossRef] [Google Scholar]
  64. Zhang, F., Parker, D., Kim, Y. S., et al. 2011, ApJ, 728, 141 [Google Scholar]

Appendix A: Benzyne line parameters

Line parameters were derived from a Gaussian fit to the observed lines using the GILDAS package. A velocity coverage of ±15 km s−1 was selected for each line. Observed frequencies were derived assuming a local standard of rest velocity of 5.83 km s−1 (Cernicharo et al. 2020b). Observed velocities for the detected lines were determined assuming that rest frequencies are those predicted from the rotational and distortion constants provided by the fit to the laboratory data (Brown et al. 1986; Kukolich et al. 2003; Robertson et al. 2003). The predicted frequencies in MADEX (Cernicharo 2012) agree within 1-5 kHz with those of the CDMS (Müller et al. 2005). Differences between predicted and observed frequencies never exceed 1.5σΔν, where Δν is the estimated frequency uncertainty of the observed lines. The averaged VLSR from the 11 observed lines of benzyne is 5.80±0.05 km s−1, a value similar to that derived from all lines of HC5N and its isotopologues in the Q band (Cernicharo et al. 2020b).

Table A.1.

Observed line parameters of benzyne.

Appendix B: The chemical network for the formation of benzyne

The low number densities of the reactants along with the typical 10 K translational temperature of a cold starless core such as TMC-1 only support barrierless and exoergic bimolecular reactions. Focusing on neutral-neutral reactions, the scheme is A + B →[AB]* → C + D, where A and B represent the reactants, [AB]* is the reaction intermediate(s), and C and D are the products (Kaiser 2002). In TMC-1, the newly detected o-C6H4 defines the product C; atomic hydrogen represents one of the most prominent counter fragments, D, that carry away a significant amount of the total angular momentum. This concludes that [AB]* has to fulfil the stoichiometry of [C6H5]*. Since the intermediate holds six carbon atoms, this can be achieved in bimolecular reactions in C1-C5, C2-C4, and C3-C3 systems with reactants at various degrees of hydrogenation. In principle, the following reactions, (1) to (14), could yield rovibrationally excited [C6H5]* reaction intermediates, which undergo unimolecular decomposition via atomic hydrogen loss to form the newly detected o-C6H4 molecule.

(1) C + C5H5o-C6H4 + H

(2) CH + C5H4o-C6H4 + H

(3) CH2 + C5H3o-C6H4 + H

(4) CH3 + C5H2o-C6H4 + H

(5) CH4 + C5H → o-C6H4 + H

(6) C2 + C4H5o-C6H4 + H

(7) C2H + C4H4o-C6H4 + H

(8) C2H2 + C4H3o-C6H4 + H

(9) C2H3 + C4H2o-C6H4 + H

(10) C2H4 + C4H → o-C6H4 + H

(11) C2H5 + C4o-C6H4 + H

(12) C3 + C3H5o-C6H4 + H

(13) C3H + C3H4o-C6H4 + H

(14) C3H2 + C3H3o-C6H4 + H.

Among these reactions, only the reaction between the ethynyl radical (C2H) and vinylacetylene (CH2CHCCH), reaction (7), has been studied via crossed molecular beam experiments and electronic structure calculations (Zhang et al. 2011). This reaction is barrierless and exoergic, and thus it is expected to be fast at low temperatures. This finding is in line with the known high reactivity of C2H radicals with unsaturated hydrocarbons at low temperatures (Chastaing et al. 1998), revealing rates of up to 4×10−10 cm3s−1. Moreover, Zhang et al. (2011) predicted that o-C6H4 is formed with a branching ratio of up to 20% at 10 K, thus providing an effective rate of 8×10−11 cm3s−1. Reaction (7) is thus a plausible route to benzyne in TMC-1. Further, in analogy to the C(3P)–benzene system (Hahndorf et al. 2002), the reaction of atomic carbon (C(3P)) with cyclopropadienyl (c-C5H5) – although studied neither experimentally nor computationally – is expected to form o-C6H4 plus atomic hydrogen (reaction (1)). Hahndorf et al. (2000) predicted computationally that the propargyl radical (C3H3) can react barrierlessly and exoergically with the vinyl radical (C2H3) on the singlet surface to form c-C5H5 plus atomic hydrogen, thus providing the c-C5H5 required for reaction (1).

B.1. Discounted reactions to benzyne

Reaction (2) has never been studied in the laboratory with any isomer such as methyldiacetylene (HCCCCCH3), penta-1,4-diyne (CH2(C2H)2), or 1,2,3,4-pentatetraene (H2CCCCCH2). Rate constant measurements of methylidyne (CH) with unsaturated hydrocarbons by Canosa et al. (1997) reveal barrierless pathways with rates of a few 10−10 cm3s−1. However, crossed molecular beam reactions indicate that the CH radical adds to the double or triple bond of, for example, acetylene, ethylene, diacetylene, allene, and methylacetylene, followed by isomerisation through the ring opening and/or hydrogen shifts, followed by hydrogen atom losses. Six-membered ring formation only occurs through the ring expansion of an existing five-membered ring (i.e. a hypothetical cyclic c-C5H4 isomer). However, no confirmed formation routes to cyclic c-C5H4 are known so far; therefore, we do not include reaction (2) as a contributor to benzyne.

Reaction (4), which has never been studied either computationally nor experimentally, involves the addition of a methyl radical (CH3) to a triplet carbene (HCCCCCH). The addition of a methyl radical to triple bonds of the ethynyl moiety involves barriers of up to 50 kJ mol−1, which cannot be overcome in TMC-1. Methyl radicals could add barrierlessly to the central carbon atom, which essentially carries the unpaired electrons. However, this requires small, nearly zero-impact parameters and an extensive reorganisation of the carbon skeleton; therefore, although exoergic by 207 kJ mol−1 overall, the formation of benzyne is unlikely.

Reaction (5) represents a direct reaction that leads, via hydrogen abstraction, to a methyl radical plus C5H2 molecules; although never studied experimentally, this hydrogen abstraction is expected to have a barrier. Further, this reaction is endoergic by 27 kJ mol−1 and hence cannot take place at 10 K. Reaction (5) cannot form o-C6H4 + H at 10 K since hydrogen abstraction represents the only feasible pathway.

Reaction (8) has an entrance barrier of addition of the C4H3 with its radical centre to the carbon-carbon triple bond of acetylene 25–37 kJ mol−1. Hence, its reaction is closed at 10 K as well.

Reaction (9) involves the reaction of a vinyl radical with diacetylene; reactions of vinyl radicals with closed-shell olefines and alkynes have entrance barriers of 10 to 35 kJ mol−1, which cannot be overcome in TMC-1. This has been confirmed computationally (Zhang et al. 2011). Finally, reaction (11) and, in particular, reaction (12) likely have entrance barriers (Mebel & Kaiser 2015).

All Tables

Table A.1.

Observed line parameters of benzyne.

All Figures

thumbnail Fig. 1.

Scheme of the ortho, para, and meta isomers of benzyne.

In the text
thumbnail Fig. 2.

Subset of the observed lines of o-C6H4 in the 31−50 GHz frequency range towards TMC-1 (see Table A.1 for the complete list of detected benzyne lines). The abscissa corresponds to the rest frequency assuming a local standard of rest velocity of 5.83 km s−1. The ordinate is the antenna temperature corrected for atmospheric and telescope losses, in mK. The red line shows the synthetic spectrum obtained from a fit to the observed line profiles, which provides Tr = 7 ± 1 K and N(o-C6H4) = (5.0 ± 1.0) × 1011 cm−2. The rotational quantum numbers are indicated in each panel. Blanked channels correspond to negative features produced in the folding of the frequency switching data. Cyan numbers indicate the multiplicative factor applied to the best model, when needed, to match the observations.

In the text
thumbnail Fig. 3.

Calculated fractional abundance of o-C6H4 as a function of time in models A and B (see text). The horizontal grey band corresponds to the observed abundance in TMC-1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.