Open Access
Issue
A&A
Volume 693, January 2025
Article Number A277
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202452210
Published online 23 January 2025

© The Authors 2025

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1 Introduction

Interstellar polyynes are a large class of chemical compounds that can be described with the common formula R1–(C)2n–R2, where n > 1. Due to the presence of the conjugated unsaturated triple bonds, these chemical compounds are typically linear in shape and, along with the simplest alkynes HCCH (C2H2, acetylene) and HCCCH3 (C3H4 propyne), they are often referred as “carbon chains”. Polyynes and other carbon chain molecules comprise the large fraction of species identified to date in the inter- and circumstellar medium (ISM and CSM). Taniguchi et al. (2024) reported the detection of more than 100 carbon- chain species. This accounts for about a third of species currently identified in the ISM.

Polyynes and carbon chain molecules are ubiquitous in the ISM. They have been detected in extragalactic sources, photodissociation regions (PDRs), planetary nebulae, and the envelopes around carbon-rich AGB stars (Pardo et al. 2005, 2022; Cuadrado et al. 2015; Martín et al. 2021; Taniguchi et al. 2024). Observations of the Taurus Molecular Cloud-1 (TMC-1) using the 40-meter Yebes telescope and the 100-meter Green Bank telescope resulted in the identification of many new gas phase carbon chain constituents. These developments have raised questions about their associated formation and destruction routes (Loomis et al. 2021; Cernicharo et al. 2021a,b, 2022; Fuentetaja et al. 2022a,b; Taniguchi et al. 2016, 2017, 2019). A steady growth in the number of carbon chain molecules identified towards protostars and protoplanetary disks can be outlined (Sakai & Yamamoto 2013; Bergner et al. 2018; Ilee et al. 2021; Guzmán et al. 2021; van Dishoeck et al. 2023).

It is generally assumed that carbon chain molecules in starforming regions are primarily produced in the gas phase through various ion-molecule and neutral-neutral reaction sequences. Astrochemical numerical simulations have a clear focus on the introduction of new gas-phase formation routes explaining the abundance of the newly observed species. Sometimes, the introduction of new gas-phase consumption routes is required to avoid the overproduction of one of the observed intermediate species. Even though solid-state reactions on the surface of interstellar dust grains are often required to produce the initial precursors for these carbon chain molecules gas-phase formation routes, systematic studies of carbon chain reactivity on the grain surface are lacking.

The dust temperature in the densest regions of molecular dark clouds drops to as low as ~ 10 K, resulting in the continuous freeze-out of the species, including atoms and molecules previously formed in the gas phase. The most prominent example is the freeze out of CO molecules produced in a sequence of gas-phase reactions starting from C+ ions. A CO freeze-out is the initial point for the formation of various complex organic molecules (COMs) as well as carbon chain oxides (Gerakines & Moore 2001; Öberg et al. 2009; Modica & Palumbo 2010; Maity et al. 2015; Fedoseev et al. 2015, 2022; Chuang et al. 2017; Qasim et al. 2019a; Urso et al. 2019). Surface interactions of accreting CO molecules and H-atoms result in the formation of CH3OH ice and organic molecules as complex as glycerol, a key component of lipids and the membranes of living cells (Fedoseev et al. 2017). Unsaturated triple bonds comprising polyynes have high reactivities with accreting H atoms and other open-shell species, such as free radicals, even at 10 K temperature (Molpeceres & Rivilla 2022). Indeed, dedicated laboratory studies performed by Hiraoka et al. (2000) and Kobayashi et al. (2017) revealed high reactivity of the simplest alkyne C2H2 with H atoms, resulting in the formation of C2H4 (CH2CH2, ethylene) and C2H6 (CH3CH3, ethane). Subsequently, Qasim et al. (2019b) reported similar results for the surface hydrogenation of C3H4. The formation of C3H6 (CH2CHCH3, propene) followed by the ultimate formation of the fully saturated C3H8 (CH3CH2CH3, propane) at 10 K. Saturation of triple bonds is also observed in the presence of other radicals upon (non)energetic processing. Chuang et al. (2020, 2021) and Qasim et al. (2019b) investigated the solid-state reactivity of C2H2 and C3H4 with OH (hydroxyl) radicals, a well-known intermediate in forming H2O ice. It was found that the addition of OH radicals to the triple CC bonds can proceed at 10 K. The activation barrier for this reaction is expected to be lower than that for the similar addition of H atoms. Moreover, the addition of OH radicals to the double CC bonds proceeds with even lower activation barriers, see Qasim et al. (2019b) and Molpeceres & Rivilla (2022) for more details. The final products obtained after the simultaneous exposure of C2H2 and C3H4 to OH radicals and H atoms are various alcohols and aldehydes (ketones). These new reaction routes have been used to explain recent detections of i-C3H7OH (2-propanol) in Sgr B2(N2) in Belloche et al. (2022). They used a modified rate equation code accounting for the non-diffusive reactions, such as the interactions of species with immobile OH radicals. Similarly the routes provided by Qasim et al. (2019a) have been suggested to explain the first detection of n-C3H7OH (n-propanol; Jiménez-Serra et al. 2022). So far, the simplest alkynes HCCH and HCCCH3 are the only representatives of H–(C)2n–H and H–(C)2n–CH3 rows of alkynes and polyynes, whose surface reactivity with H atoms and OH radicals are systematically investigated.

This work aims to provide the first experimental confirmations for surface hydrogenation of H–(C)2nH (n > 1) polyynes under the conditions relevant to the early stages of star formation. Based on the outcome of the aforementioned experimental and theoretical studies, we expect that hydrogenation of H–(C)2n–H (n > 1) polyynes results in the intermediate formation of various semi-saturated and fully saturated linear alkanes described by the formula C2nH4n+2. This will result in the presence of linear alkanes on the surface of icy grains already at prestellar stages. The partial survival of the pristine ice through the later stages of star formation would enrich celestial bodies with alkanes. Linear alkanes, including C2H6, C4H10, and C6H14, have been identified in comet 67P/Churyumov–Gerasimenko with high abundances of ~5 × 10−1 (C2H6) and ~5 × 10−3 (C4H10) with respect to the H2O (Schuhmann et al. 2019; Altwegg et al. 2017). Moreover, the Hayabusa2 mission reported the presence of the saturated aliphatic chains identified in Ryugu’s organic matter (Yabuta et al. 2023).

The surface chemistry of H–(C)2n–H polyynes with n > 1 has not been investigated as systematically in the laboratory due to the difficulties associated with the acquisition and preservation of polyyne samples. Experimental samples synthesized in the laboratory are kept at liquid-N2 temperature to avoid spontaneous explosive decomposition (see Khlifi et al. 1995 and Kim & Kaiser 2009, for example). In the present study, a twostep approach has been utilized. For the first time, a mixture of polyynes was obtained in situ on the substrate and the reactivity of this mixed polyyne ice with H-atoms was investigated in the same experimental setup.

Section 2 describes the applied experimental methodology and utilised ultra-high vacuum (UHV) cryogenic setup. In Sect. 3, the experimental results are presented and discussed. Section 4 focuses on the astrochemical implications and offers our conclusions.

2 Experimental

All experiments are performed using an ultra-high vacuum (UHV) setup, SURFRESIDEUV, with a UV-photon source module. The original SURFRESIDE design is optimized to study various astrochemically relevant solid-state atom- and radicaladdition reactions at temperatures as low as 10 K, as described in Ioppolo et al. (2013). The implementation of a UV-photon source allows us to make a quantitative comparison between the outcome of these “non-energetic” reaction routes, i.e., the reaction routes between the species in the thermal equilibrium with the surface, and the outcome of UV photon-induced “energetic” reaction routes (Chuang et al. 2017). The description of the upgrade of the setup is presented in Fedoseev et al. (2016) and Chuang et al. (2018a). In this work, a UV-photon source is used to produce the thermodynamically unstable species directly in the ice samples. This allows for the study of the chemical reactivity between species that cannot be introduced into UHV setup by other means (see, e.g., Butscher et al. 2015, 2016, 2017). Following this latter approach, a mixture of unsaturated hydrocarbons, including unstable C4H2 and C6H2, can be synthesized in situ by UV irradiation of C2H2 ice. Then the reactivity of the newly produced hydrocarbons with H atoms is investigated following the exposure of the ice samples to an H-atom beam. It should be noted that this experimental approach does not aim to mimic the exact composition of ice mantles in the ISM or their complex interactions with impinging UV-photons and H-atoms. Instead, it is aimed at an experimental verification that the hydrogenation of polyynes occurs on a 10 K surface by accreting H atoms using astronomically relevant H-atom fluences.

2.1 Experimental setup

The SURFRESIDEUV setup is composed of a UHV main chamber with a typical base pressure of ~ 10−10 mbar. The background residual gas mainly comprises H2, while the accretion rate of background H2O is negligible (<6.3 × 1010 molecules cm−2 s−1, Chuang et al. 2018b) for the period of the experiment. A rotatable gold-plated copper substrate is mounted on the cold tip of a helium closed-cycle cryostat that is positioned in the center of the UHV main chamber. The temperature of the substrate is regulated in the range between 8 and 450 K by means of resistive heating and monitored using two silicone diodes with the absolute accuracy of ±0.5 K. The C2H2:He gas mixture or Ar gas is admitted to the main chamber through high-precision all-metal leak valves with the 22° angle to the substrate surface normal. Two independent dosing lines are used to keep these gases for the whole extent of the experiment. A mixture of C2H2 in He “bathing” gas is used to avoid the risk of spontaneous C2H2 decomposition. The substrate temperature is maintained far above the He condensation point. This prevents the formation of He ice in all experiments. The mixed beam of H atoms and undissociated H2 molecules is obtained by thermal cracking of H2 molecules, using a Hydrogen Atom Beam Source (HABS; Tschersich 2000). HABS is mounted in an independent UHV chamber separated from the main chamber by a UHV shutter. A U-bent quartz pipe is placed along the path of H-atom beam to efficiently quench the high translational or ro-vibrational excited states of H-atoms and H2 -molecules prior to their impact on the substrate under the incident angle of 45°. The application of two independent commercial H-atom beam sources allows for the crosscheck of the obtained results. UV-photons are generated by a microwave-discharged hydrogen flowing lamp (MDHL) attached to a MgF2 viewport (with a cutoff at 115 nm) on the main chamber. The MgF2 viewport is positioned in front of the substrate parallel to the substrate plane. This configuration allows incident photons to cover uniformly the whole area of the 2.5 × 2.5 cm2 substrate. The H-atom and UV-photon fluxes are calibrated in situ, and the details of this procedure are available from Ioppolo et al. (2013) as well as Fedoseev et al. (2016) and Ligterink et al. (2015), respectively. The H-atom flux utilised in the present study amounts to 2.4 × 1013 atoms cm−2 s−1. The lower limit on UV-photon flux at the substrate plane is estimated to be ~1.3 × 1013 cm−2 s−1. The example of UV emission spectra obtained for the used F-type MDHL is available in Chen et al. (2014) and Ligterink et al. (2015). The following commercially available reagents are used in the experiments: a 5% mixture of HCCH in He (Praxair), H2 (Linde, 5.0), and Ar (Linde, 5.0).

2.2 Experimental methods

The experiments are performed in the following procedure. First, a “pillow” of Ar ice (10 Langmuir1) is deposited on top of a bare gold substrate. Then a ~10 monolayer (ML) thick C2H2 ice layer is deposited on top of the predeposited Ar ice layer. The layered ice sample is exposed to the UV-photon beam generated by MDHL to achieve partial polymerization of C2H2 into its C2nH2 (n > 1) polyynes derivatives. This polymerization has been extensively studied by Cuylle et al. (2014) using a combination of mid-infrared and ultraviolet-visible (UV-VIS) spectroscopy (see also Compagnini et al. 2009; Wu et al. 2010; Puglisi et al. 2014; Abplanalp & Kaiser 2020). The thickness of the grown C2H2 ice is well below the penetration depth of UV-photons. This aims to secure the uniform distribution of the formed products in the photo-processed ice. On the other hand, the interaction of UV-photons passing through the ice with the gold substrate may result in the production of photoelectrons. Hopkins & Riviere (1964) reported the work function of >4.7 eV for gold. The deposited ten Langmuir layer of Ar ice aims to prevent the interaction of the HCCH ice with the photoelectrons due to the limited penetration depth of photoelectrons in condensed phases (Jo & White 1991).

The ice composition is monitored in situ by means of Fourier-transform reflection absorption infrared spectroscopy (RAIRS). Upon completion of the UV-exposure stage, a quadrupole mass spectrometer-temperature programmed desorption (QMS-TPD) measurement is performed. The substrate is rotated by 135° to face directly the ion source of the QMS. Then, during the TPD measurements, the ice is gently warmed up with a 2 K/min rate from 10 to 50 K to remove the “pillow” of Ar ice. Subsequently, the 5 K/min heating rate is used in the range of 50 to 250 K to increase the sensitivity of the applied TPD techniques. The sublimating species are continuously monitored by means of QMS. The formed species are then identified through the presence of their characteristic m/z (mass-to-charge) values at characteristic desorption temperatures. Thus, two independent parameters are used to identify the obtained products.

The described experiment is then repeated using the very same experimental conditions followed by irradiation of the obtained ice with an overabundance of H atoms. Upon completion of the H-atom exposure stage, QMS-TPD measurement is performed by applying the same heating routine. The difference in the outcome of the performed experiments is used to identify the presence of four- and six-carbon bearing products of C2nH2 (n > 1) polyynes hydrogenation, formed along with the expected products of well-investigated C2H2 hydrogenation (Hiraoka et al. 2000; Kobayashi et al. 2017; Chuang et al. 2020).

Only the top monolayer of C2H2 ice is assumed to be available for hydrogenation by incident H-atoms (Ioppolo et al. 2013; Chuang et al. 2018b). This assumption is in agreement with the results of the pure HCCH ice hydrogenation experiment reported by Kobayashi et al. (2017). The hydrogenation of 12 ML of pure HCCH ice at 10 K resulted in a relative HCCH consumption of 0.06 at an equilibrium state (see Fig. 7 of Kobayashi et al. 2017). This observation allows for the derivation of RAIRS setupspecific C2H2 band strength value. With this goal, pure C2H2 ice is hydrogenated upon reaching the maximum C2H2 consumption. Then, the difference in the total area for the selected C2H2 absorption features is obtained at the end of the experiment. This area corresponds to the absorption by ~1 ML of C2H2 ice (where 1 ML surface coverage is assumed to be 1 × 1015 cm−2). With the geometry utilized in our setup, the absorbance area equal to 0.028 cm−1 was obtained for the C-H stretching mode of C2H2 (v3, 3244 cm−1). This corresponds to 0.064 cm−1 in optical depth, resulting in the setup-specific RAIRS band strength values equal to 6.5 × 10−17 cm molecule−1.

3 Results and discussion

3.1 Application of QMS-TPD technique

The QMS-TPD spectra obtained after the UV-exposure of the deposited C2H2 ice are presented in the left panel of Fig. 1. It is important to note that the ordinate axis is presented as a logarithm of 10 to highlight the appearance of the small desorption signals following the deviations from the exponential decays. Only the most relevant m/z values are selected and shown for better clarity. Three distinct TPD peaks corresponding to the desorption of three individual H–(C)2n–H alkynes or polyynes can be clearly observed in the spectrum. These are C2H2 (m/z = 26), C4H2 (m/z = 50), and C6H2 (m/z = 74). Under our experimental conditions, the desorption of C2H2 starts at about 55 K and peaks at 75 K, while the desorption of C4H2 occurs at 95 K and peaks at 120 K. These temperature ranges are in agreement with the previous experimental results (see Zhou et al. 2009, 2010; Abplanalp & Kaiser 2017). To our knowledge, no literature value has been reported for pure C6H2 ice. Under our experimental conditions, the corresponding m/z = 74 peak appears at 125 K and reaches the maximum at about 150 K. Several additional m/z values are presented in the spectra for reference. The m/z = 30 and m/z = 58 correspond to the “parent masses” of C2H6 and C4H10, i.e., the masses of ions obtained by the direct non-dissociative ionisation of C2H6 and C4H10, respectively. While the m/z = 29 and m/z = 43 spectra correspond to the signals of C2H5+ and C3H7+, the two most likely ions obtained by dissociative ionization of various short linear aliphatic hydrocarbons in the ion head of the QMS. Appendix B in Schuhmann et al. (2019) offers typical examples of mass-spectra obtained for the short linear hydrocarbons. The QMS-TPD spectra presented in the left panel of Fig. 1 reveal the traces of C2H6 ice through the rise of m/z =30 and m/z = 29 signals in the range from 55 to 80 K. No signs of other saturated hydrocarbons can be found.

The right panel of Fig. 1 presents a general view of the QMS-TPD spectra obtained after the hydrogenation of the UV- exposed C2H2 ice under the same experimental conditions. Here, a minor drop in the intensity of C2H2 , C4H2, and C6H2 desorption features can be distinguished (see Fig. 2 for more details). On the other hand, the appearance of new desorption features corresponding to the saturated aliphatic hydrocarbons can be observed. The rise of m/z = 30 (C2H6+) and m/z = 29 (C2H5+) signals in the temperature range between 55 and 85 K can be attributed to the desorption of C2H6. The rise of m/z = 58 (C4H10+), m/z = 43 (C3H7+) and m/z = 29 (C2H5+) in the range from 95 to 115 K can be attributed to the desorption of C4H10. The third desorption feature in the range from 130 to 150 K can be tentatively assigned to semi-saturated or fully saturated hydrocarbons with six carbon atoms through the rise of m/z = 29 (C2H5+). More details are given in Fig. 3.

The detailed comparison between the QMS-TPD curves obtained for the H−(C)2n−H polyynes with and without the hydrogenation of the ice is presented in Fig. 2. The linear scale is used for all of the plots for the direct comparison of the corresponding areas. The m/z signals assigned to C2H2 (HCCH, m/z = 26) and C4H2 (HCCCCH, m/z =50) demonstrate about 0.1 drop in the intensity upon hydrogenation of the ice. In the case of C6H2 (HCCCCCCH, m/z=74) molecules, the observation of such a drop can only be considered tentative because of the low signal-to-noise ratio. These drops in the areas of TPD curves in the right panel of Fig. 2 are due to the consumption of C2H2, C4H2 and C6H2 molecules in the upper layer of the UV-exposed ice through the reactions with H atoms. These experimental results show that the reactivity of unsaturated triple bonds of C4H2 (and, likely, C6H2) with H-atoms is similar to the reactivity of the C2H2 triple bond (see also Hiraoka et al. 2000; Kobayashi et al. 2017). The shape and profile changes in the TPD curves obtained after the hydrogenation of the ice can be attributed to the more complex (non-uniform) composition of the ice. Indeed, the binding energies for polyynes are higher than those for their corresponding saturated counterparts (see, e.g., Abplanalp & Kaiser 2017). This can be explained by the interactions between π-electrons of their triple bonds. This difference in the binding energies can affect the shape and profile of polyyne desorption features obtained for the mixed ice.

In addition to the data presented in Fig. 2, we present the comparison between the QMS-TPD curves corresponding to the fully saturated aliphatic hydrocarbons obtained with and without H-atom exposure to the ice in Fig. 3. The dissociative ionisation of various linear saturated aliphatic hydrocarbons in the ion source of the QMS results in the efficient cleavage of both C-C and C-H bonds. The upper panels in Fig. 3 show QMS-TPD curves obtained for the “parent masses” of C2H6 (CH3CH3+, m/z = 30) and C4H10 (CH3CH2CH2CH3+, m/z = 58). At the same time, the m/z = 29 (CH3CH2+) and m/z = 43 (CH3CH2CH2+) correspond to the common fragments obtained after dissociative ionization of various linear aliphatic hydrocarbons. The peak in the range from 55 to 85 K can be confidently assigned to C2H6. The m/z = 29 to m/z = 30 signal ratio of 0.8 is consistent with the ratio reported for C2H6 in the NIST database2. The clear peaks for the characteristic m/z signals of linear C4H10 equal to 58, 43 and 29 are observed in the range from 90 to 105 K. These desorption temperatures are in close agreement with the values reported previously for C4H10 (Abplanalp & Kaiser 2017, the structural isomer is not specified). In addition to the two resolved desorption peaks, a weaker desorption feature can be observed in the range from 125 to 150 K for m/z = 29 (CH3CH2+) and m/z = 43 (CH3CH2CH2+). Observation of these dissociative ionization products is in agreement with the possible presence of C6H14 (CH3CH2CH2CH2CH2CH3); namely, a fully saturated derivative of C6H2 (HCCCCCCH). The observed desorption temperature range overlaps with that reported by Abplanalp & Kaiser (2017) for various hydrocarbons consisting of six carbon atoms. However, without the identification of m/z = 86 “parent mass”, this attribution is considered tentative. The comparison between the observed intensities of m/z signals and the standard mass spectra available from the NIST database is presented in Appendix B.

It should be noted that the m/z signals associated with the presence of semi-saturated hydrocarbons can also be identified in the mass spectra. However, the unambiguous identification of semi-saturated hydrocarbons becomes increasingly difficult due to the high number of possible structure isomers and the existing overlap between the m/z signals of fully saturated and semi-saturated hydrocarbons. The combination of C4H4, C4H6, and C4H8 alone has more than 20 stable structural isomers, excluding isotopes. Meanwhile, C6Hn semi-saturated hydrocarbons have more than 100 isomers. Some observations can be made based on the obtained laboratory data, see Fig. 4. It is seen that prior to the hydrogenation of the UV-exposed ice, the peaks of m/z = 52 (C4H4+) and m/z = 54 (C4H6+) signals are already present in the spectra. The desorption starts from 100 K and peaks in the range from 115 to 120 K. The m/z = 52 value corresponds to the “parent” and most intense signal of vinylacetylene (C4H4, CH2CHCCH). The desorption temperature registered for the m/z = 52 signal is higher than that observed for the m/z = 58 signal of C4H10 and is consistent with the literature value reported for the pure CH2CHCCH ice by Kim & Kaiser (2009) (see Fig. 4 of Kim & Kaiser 2009). Generally, a gradual shift towards lower desorption temperature is expected with the increase of hydrogenation degree, namely, Tdep(C4H2) > Tdep(C4H4–8) > Tdep(C4H10), as, for example, described in Abplanalp & Kaiser (2017). The formation of CH2CHCCH upon UV exposure of C2H2 ice is in agreement with the experimental results of Cuylle et al. (2014). The exact identification of produced semi-saturated hydrocarbons requires other experimental techniques and is beyond the scope of this work.

A weak desorption feature observed for the m/z = 54 (C4H6+) signal can be attributed to the traces of one of the C4H6 isomers in the irradiated ice. The most probable candidates are ethylacetylene (CH3CH2CCH) and divinyl (CH2CHCHCH2). Here, the preference can be given to CH3CH2CCH. The observed m/z = 54 signal is the most intense peak in the mass spectra of CH3CH2CCH and corresponds to its non-dissociative ionization in the ion source of the QMS. Previous laboratory studies showed that surface hydrogenation of HCCH and HCCCH3 triple bonds proceeds with a higher activation barrier than the hydrogenation of CH2CH2 and CH2CHCH3 double bonds, see Kobayashi et al. (2017) and Qasim et al. (2019b). Similar behavior can be expected for the hydrogenation of CH2CHCCH triple and double bonds. Preferable formation of CH3CH2CCH can be suggested in this case. Although CH2CHCHCH2 is thermodynamically more stable than CH3CH2CCH due to the presence of the conjugated double bonds, the preference should be given to the kinetic argument. Only traces of the C4H8+ m/z desorption feature can be observed in the QMS-TPD spectrum of UV-exposed ice prior to hydrogenation.

The m/z curves obtained after the hydrogenation of the ice reveal the overall consumption of C4H4 and the formation of C4H6 and C4H8. This observation is consistent with the increase of hydrocarbon hydrogenation degree following various H-atom addition reactions in the mixed ice. The simultaneous rise of m/z = 56 (C4H8+) and m/z = 27 (CH2CH+) signals in the range from 100 to 115 K is consistent with the desorption of 1-butene (CH2CHCH2CH3). The lower desorption temperature registered for m/z = 56 (C4H8+) in comparison to m/z = 52 (C4H4+) and m/z=54 (C4H6+) is in line with the higher hydrogenation degree of CH2CHCH2CH3. It should be noted that the rise of m/z = 56 and m/z = 27 signals is also observed in the range from 125 to 150 K. In this range, the desorption of 6-carbon bearing hydrocarbons should occur, see Abplanalp & Kaiser (2017). Both of these m/z values are present in the mass spectra of fully saturated linear C6H14.

thumbnail Fig. 1

Left: QMS-TPD spectra in the range from 50 to 225 K for the selected m/z values obtained after the UV-exposure of 10 ML of pure HCCH ice with 4.6 × 1016 photons cm−2 at 10 K. Right: QMS-TPD spectrum obtained after the hydrogenation of UV-exposed ice with 1.3 × 1017 cm−2 atoms at 10 K. The ordinate axis is presented as a logarithm of 10. An exponential decay of ion intensity appears as a straight line.

thumbnail Fig. 2

Direct comparison between the fragments of QMS-TPD spectra in the range from 50 to 225 K obtained after the UV-exposure of 10 ML thick pure HCCH ice at 10 K with >4.6 × 1016 cm−2 photons (black empty circles) and the subsequent exposure of this photoprocessed ice with 1.3 × 1017 cm−2 H atoms at 10 K (red filled circles). The results obtained upon subtraction of the initial mass spectrum from the mass spectrum obtained upon hydrogenation of the ice for each of the presented m/z values are presented in Appendix A.

thumbnail Fig. 3

Similar to Fig. 2. M/z signals corresponding to linear saturated aliphatic hydrocarbons are shown.

thumbnail Fig. 4

Similar to Figs. 2 and 3. m/z signals associated with the semisaturated aliphatic hydrocarbons are presented.

thumbnail Fig. 5

Left: RAIR spectra obtained after the deposition of 10 ML of HCCH ice at 10 K (a, black line), after the exposure of this obtained ice with 4.6 × 1016 cm−2 photons (b, red line) and by subtracting the initial spectrum (a) from spectrum (b) (green line). Right: RAIR spectra obtained after the UV-exposure of 10 ML thick pure HCCH ice at 10 K with 4.6 × 1016 cm−2 photons (a, red line), the consequent exposure of this photoprocessed ice with 1.3 × 1017 cm−2 H atoms at 10 K (b, blue line), and by subtracting the initial spectrum (a) from spectrum (b) (green line). The grey dashed and dotted lines are the reference RAIRS spectra obtained for the pure C4H2 and C4H4 ices by Zhou et al. (2009) and Kim & Kaiser (2009). All spectra are offset for clarity.

3.2 RAIRS and kinetic analysis

Complementarily to the QMS-TPD data, RAIRS results provide information about the ice composition at 10 K prior to the sample thermal processing and consequent destruction during the TPD experiment. Moreover, information about the phase changes and kinetic data can be acquired during the UV-processing and H- atom exposure of the ice. In the left panel of Fig. 5, the RAIR spectra obtained after the deposition of pure C2H2 ice at 10 K on top of the 10 Langmuir Ar “pillow” (spectrum a), and the consequent UV-exposure of C2H2 ice (spectrum b) are presented. The obtained RAIR difference spectrum; namely, (b)-(a), is added for comparison. A clear decrease in the total area of the absorption features associated with the CH stretching mode (ν3, 3265–3210 cm−1), CCH bending mode (ν5, 825–725 cm−1), and the combinational mode (ν4+ν5, 1430–1350 cm−1) is observed, indicating the consumption of HCCH. Moreover, the morphological changes in HCCH ice can be observed following the disappearance of C-H stretching and ν4+ν5 combinational mode absorption features associated with the contribution of “polycrystalline” HCCH ice (i.e., 3226 and 1390 cm−1) and the corresponding growth of the C–H stretching mode absorption feature (i.e., 3240 cm−1) associated with “amorphous” HCCH ice (see Fig. 8 of Hudson et al. 2014). The “amorphous” and “polycrystalline” terms are used in quotes, as this band profile modification can also be attributed to the different orientation of HCCH molecules located in the interphase with vacuum or with the Ar ice “pillow”, in comparison to the orientation of HCCH molecules in the bulk of HCCH ice. The effects in the interphase should impact only a few adjunct MLs of the ice. Thus, the contribution of the interphase is consistent with the low 10 ML thickness of the ices used in present study, in comparison to the thick µm ice used by Hudson et al. (2014). A clear growth of the absorption feature in the range from 3308 to 3265 cm−1 centered around 3285 cm−1 is observed. This feature can be successfully assigned to the C-H stretching modes of both HCC- CCH and CH2CHCCH (Tørneng et al. 1980; Zhou et al. 2009; Kim & Kaiser 2009), one of the most intensive modes for both of the species. The assignment of CH2CHCCH can be further supported by the observation of a broad merged absorption feature in the range from 1020 to 900 cm−1 with two peaks centered around 975 and 935 cm−1. These features can be tentatively assigned to the ν14 and ν15 modes of CH2CHCCH (Tørneng et al. 1980; Kim & Kaiser 2009), along with the ν7 mode of C2H4. These observations further support the assignments done by the QMS-TPD technique, confirming the formation of both species already at 10 K.

In the right panel of Fig. 5, the RAIRS spectra (a) and (b) obtained prior to and after the H-atom exposure of the ice are presented. In a similar way, the RAIR difference spectrum (b)–(a) is shown for comparison. The expected decrease in the intensity of absorption features associated with unsaturated hydrocarbons (i.e., HCCH, HCCCCH and CH2CHCCH) is observed. Simultaneously, multiple absorption features appear in the range from 2990 to 2870 cm−1 and can be assigned to the C-H stretching modes of newly formed saturated hydrocarbons. The zoom-in spectral region of interest is presented in Fig. 6. Three characteristic absorption features of CH3CH3 are observed at 2975, 2943, and 2882 cm−1 corresponding to the v10, v8+v11 and ν5 modes, respectively (Comeford & Gould 1961; Kim et al. 2010; Turner et al. 2018). However, in addition to the indicated absorption features of CH3 CH3, the traces of two distinct absorption features of CH3CH2CH2CH3 can be tentatively identified at 2959 (v12, v27) and 2925 cm−1 (v13) corresponding to the asymmetric CH stretching modes of CH3– and –CH2– respectively (Tørneng et al. 1980; Kim et al. 2010; Turner et al. 2018). To further outline the plausible presence of –CH2– asymmetric stretching mode of CH3CH2CH2CH3 the reference RAIRS spectrum of C2H6 obtained from Öberg et al. (2009) is plotted versus the data obtained in the hydrogenation experiment. The asymmetric stretching mode of −CH2− at 2925 cm−1 can be considered a characteristic absorption feature for all saturated linear aliphatic hydrocarbons (see, e.g., Yabuta et al. 2023). The RAIR difference spectrum obtained after hydrogenation of the ice does not reveal the growth of any absorption features that can be assigned to the semi-saturated hydrocarbons with C=C double bonds.

In first approximation the column density of the species is directly proportional to the total absorbance areas of their corresponding RAIRS absorption features. The “in situ” RAIRS technique allows for tracking the evolution of the total areas of the assigned absorbance features over the time of the H-exposure experiment. Figure 7 presents the obtained kinetic curves for the species identified in RAIR data. The total areas of absorption features assigned to the ≡CH stretching modes of HC≡CH (3265−3210 cm−1) or HC≡CC≡CH and CH2≡CHC≡CH (3308–3265 cm−1) were used to track the consumption of these species following hydrogenation of their C≡C triple bonds. The formation of C2H6 is tracked through the absorption of its CH asymmetric stretching mode (v10, 2995−2964 cm−1). The integration of the peak (2929−2909 cm−1) is used to track C4H10 following the tentative assignment of this absorption feature to the CH asymmetric stretch of –CH2–. The modulus of absorbance is taken for C2H2 and C4H2 consumed in the experiments. This allows for the direct comparison with produced species. The presented kinetic curves show that the total areas of all selected features follow a very similar profile and reach the plateau around the same hydrogenation time. This indicates that the rates (or activation barriers) of H-atom addition reactions to the triple bonds of C2H2, C4H2, and C4H4 have similar values. Otherwise, a delay in the appearance of a plateau could be observed for the kinetic curve of a molecule with a slower H-atom addition rate. The lack of shift between the C2H2 and C2H6 kinetic curves is observed. This is in agreement with the faster rate for hydrogenation of C=C double bonds of the intermediate C2H4 previously reported by Kobayashi et al. (2017). The slower hydrogenation rate of the intermediate would result in a delay in the appearance of the fully saturated hydrocarbons with respect to the consumption of the initial reagent. Similarly, there is no shift between the C4H2 and C4H10 kinetic curves. Thus, a faster rate for hydrogenation of the intermediate C4H4−8 double bonds in comparison to the hydrogenation of C4H2 triple bonds can be inferred. This observation is also consistent with the lack of absorption features that can be assigned to the semi-saturated hydrocarbons with double bonds. The plateau appears in the range of H-atom fluence from 3 × 1016 to 6 × 1016 cm−2. Assuming the H-atom flux in dark clouds is ~104 atoms cm−2 s−1, these values of H- atom fluences are obtained at ~1 × 105 and ~2 × 105 years, respectively. It should be noted that these periods fall within the characteristic timescales of the dark clouds stage (Boogert et al. 2015; Chevance et al. 2020).

thumbnail Fig. 6

Zoom-in into the CH stretching vibration modes region for the RAIR difference spectrum obtained after hydrogenation of UV-exposed C2H2 ice (green line). See the right panel of Fig. 5 for more details. The spectrum of pure C2H6 ice at 20 K is adapted from Öberg et al. (2009) (black dots). The spectrum of pure C4H10 at 10 K is obtained by Strazzulla & Baratta (1992) (blue line). All spectra are offset for clarity.

thumbnail Fig. 7

Kinetic curves obtained for C2H2, C6H2, C4H2, and C4H10. Note: the area of the (ν4) absorption feature of C4H2 cannot be deconvoluted from the contribution of the (ν1) absorption feature of C4H4. Thus, the combined area is used for the integration. The modulus of absorbance is taken for C2H2 and C4H2 consumed in the experiments. This allows for the direct comparison with produced species. The ordinate axis is presented as a logarithm of 10.

4 Astrochemical implications and conclusions

Our laboratory results show that fully saturated linear hydrocarbons (n-alkanes) such as n-C4H10 and plausibly n-C6H14 can be produced by the surface hydrogenation of the corresponding C2nH2 (H–(C)2n–H) polyynes at 10 K. This work further extends the experimental finding obtained by Hiraoka et al. (2000), Kobayashi et al. (2017) and Qasim et al. (2019b), which demonstrates the formation of C2H4, C2H6, C3H6, and C3H8 by surface hydrogenation of simplest alkynes C2H2 and C3H8, under the physical conditions resembling those founds in dark molecular clouds. Kobayashi et al. (2017) reported a four times higher rate for the hydrogenation of the HCCH triple bond than for the hydrogenation of CO molecules. This makes HCCH hydrogenation efficient at dark molecular cloud conditions. Comparable rates for the hydrogenation of HCCH, HCCCCH, and H2CCHCCH triple bonds can be proposed based on the analysis of the presented kinetic data. Even higher rates are expected for hydrogenation of C=C double bonds of intermediate semi-saturated hydrocarbons produced by hydrogenation of triple bonds (see also Kobayashi et al. 2017; Qasim et al. 2019b).

These experimental findings have a direct application to the surface chemistry at the early stages of star formation and, in particular, dark molecular clouds. This evolutionary stage is characterized by low pressures (104–106 cm−3) and temperatures (~10 K). It sets up the initial chemical composition for the future evolutionary stages of star formation, providing material for the formation of protostars and protoplanetary disks. The data obtained by various space exploration missions such as Rosetta and Hayabusa2 suggests that part of the chemical inventory produced during the dark clouds stage may survive the harsh later stages of star formation and provide pristine material for comets and asteroids (Rubin et al. 2015; Yabuta et al. 2023). The chemistry of dark clouds is driven by the interplay between gas and dust. Ubiquitous sub-µm dust grains provide the surface on which the simplest gas-phase chemical species (e.g., H, O, N, CO, etc.) accrete and react. Atom addition reactions among accreting species result in the formation of “H2 O-rich” and “COrich” ice mantles on the surface of the dust grains (Öberg et al. 2011; Boogert et al. 2015).

Interstellar carbon chains and polyynes are efficiently produced early in dark molecular clouds through cold gas-phase reaction routes, see Taniguchi et al. (2024) and references therein. The produced polyynes and carbon chains should accrete on the surface of dust grains and become available for reactions with other accreting species, such as H atoms. The efficient hydrogenation of polyynes demonstrated in the present study indicates that the formation of various linear aliphatic hydrocarbons can be expected at the beginning of the dark cloud stage. Moreover, the incorporation of other species (e.g., OH, O, N, etc.) into the polyyne hydrogenation chain will result in the formation of even more complex organic species (see Qasim et al. 2019a; Molpeceres & Rivilla 2022; Santos et al. 2024; Chuang et al. 2020, 2021, 2024).

The formation of n-alkanes and semi-saturated aliphatic hydrocarbons early in dark clouds is important for the preservation of these species during the later stages of star formation. The “H2 O-rich” and “CO-rich” ice mantles along with coagulation of the icy grains can provide limited protection from UV-photons originating both from the cosmic ray induced UV- field and the protostar. Preservation of produced linear alkanes and the intermediate semi-saturated hydrocarbons followed by their integration into the pristine material of comets is consistent with the detection of alkanes released by comet 67P/Churyumov- Gerasimenko (see Schuhmann et al. 2019; Altwegg et al. 2017). In addition, Schuhmann et al. (2019) reported that saturated alkanes account for 42% of the total amount of detected hydrocarbons in Comet 67P/Churyumov-Gerasimenko, while the contribution of unsaturated hydrocarbons is 56%. The cycloalkanes contribution is calculated to be at most 2%. A minimum contribution of cycloalkanes is consistent with the production of aliphatic hydrocarbons by the hydrogenation of linear polyynes. It should be noted that even though the reference spectra of linear alkanes were used in the study by Schuhmann et al. (2019), the contribution of branched alkanes could not be deconvoluted in that work.

Besides observations in comet 67P, hydrocarbons have been recently identified on the surface of Ceres (De Sanctis et al. 2017) and in the samples from the asteroid Ryugu delivered by the Hayabusa2 mission. Observation of asymmetric stretching modes of CH3- (2960 cm−1, 3.38 µm) and –CH2– (2930 cm−1, 3.41 µm) along with the other bands in the infrared range is used to constrain the presence of aliphatic hydrocarbons in Ryugu samples by Yabuta et al. (2023). Identification of characteristic –CH2– modes is consistent with the presence of linear alkanes comprising three or more carbon atoms or with the presence of linearly saturated hydrocarbon segments incorporated into the organic residue. Moreover, Yabuta et al. (2023) report a relatively high –CH2– to –CH3 peak intensity ratio equal to 1.9. Based on these observations and on the comparison with the ratios observed in Murchison and Ivuna meteorites, Yabuta et al. (2023) inferred that Ryugu’s organic matter may contain longer aliphatic chains than in the meteoritic insoluble organic matters (IOMs).

Gas-phase observations of saturated linear hydrocarbons are challenging because of the low dipole moment of these molecules. Alternative method for alkanes identification in the ISM can be suggested. If alkanes are present on the grain surface early during the dark cloud stage, these species can become a target for JWST observations in the mid-infrared range. Albeit the abundance of each individual alkane is expected to be low, their absorption features overlap (Tørneng et al. 1980; Turner et al. 2018). This will result in the cumulative add-up of individual absorption of each species into the combined optical depth. In this way the whole class of chemical compounds can be observed simultaneously. Besides the CH3– (2960 cm−1, 3.38 µm) and –CH2– (2930 cm−1, 3.41 µm) bands mentioned earlier, the CH3–/–CH2– deformation/scissoring modes at 1460 cm−1 (6.85 µm) and CH3– symmetric deformation mode at 1380 cm−1 (7.25 µm) are good potential candidates for such “family” spectral indicators.

Semi-saturated hydrocarbons can be observed in the gas phase. Such observations would require a mechanism responsible for the partial grain-to-gas transfer of produced intermediate species at low temperatures, such as reactive desorption (see, e.g., Minissale et al. 2016; Oba et al. 2018; Chuang et al. 2018b). Cernicharo et al. (2021b) reported the first detection of C4H4 (vinylacetylene) in TMC-1 through observation of its 6 transitions. C4H6 (ethylacetylene) was also tentatively identified in their study. So far, there have been no reports about successful detections of 1-butene, despite the fact that its structural isomer isobutene was observed in TMC-1; see Fatima et al. (2023). If solid-state formation routes contribute to the gas-phase enrichment of these species, then other semi-saturated hydrocarbons can be set as targets for future gas-phase observations.

The precise identification of intermediate four- and six-carbon bearing semi-saturated hydrocarbons in the present work was challenging. This is mainly because of the high number of possible isomers and an overlap between their characteristic m/z values in the mass spectra. The determination of the exact C4H2 and C6H2 hydrogenation sequences could serve as the subject of future experimental and theoretical studies. Nevertheless, we suggest that the preference should be given to the observations of semi-saturated hydrocarbons with triple bonds due to their higher stability.

Finally, the present work uses the two-step experimental technique to investigate the polyynes reactions with H atoms. The mixture of unstable polyynes was synthesized in situ on the substrate in the main UHV chamber of the setup by UV- irradiation of C2H2 ice. Then, the hydrogenation process of the obtained mixed ice using an H-atom beam was performed. Simultaneous hydrogenation of the mixed ice allowed us to compare the relative reactivity of involved species. This approach can be used to study the hydrogenation of other species, which might not otherwise be so easily obtained in the lab. This may include cyanopolyynes or oxygen- and sulphur-bearing carbon chains. The former, for example, can be synthesized by UV-exposure to HCN/C2H2 mixed ice. Then, their overall reactivity with H- atoms can be investigated. Analogously to the present study and the works of Kobayashi et al. (2017) and Qasim et al. (2019b), we can expect that C≡C triple bonds of cyanopolyynes can participate in hydrogenation reactions that result in the formation of aliphatic nitriles or even amines (see also Theule et al. 2011; Nguyen et al. 2019). Similarly, the hydrogenation of (H)CnO and CnS carbon chains can result in the formation of various aliphatic aldehydes, alcohols, and thiols. We give more details in Fig. 8. These substituted hydrocarbons have direct astrobiologi- cal implications as they are amphiphilic and may play a part in the formation of micelles in aqueous solutions, the most primitive prototypes of cells and organelles (Segré et al. 2001; Trevors & Psenner 2001; Stevenson et al. 2015).

thumbnail Fig. 8

Schematic representation of the overall suggested mechanism. Various polyynes observed in dark clouds (e.g., TMC-1) are presented in the left part of the reaction scheme. The accretion of polyynes on the surface of the grains is followed by hydrogenation of their unsaturated CC triple bonds, resulting in the formation of their semi-saturated and fully saturated counterparts. Besides CH2 CHCCCHO, all of the first- generation products obtained by such hydrogenation have also been observed in the gas phase in TMC-1 or other dark clouds (Brown 1981; Irvine & Hjalmarson 1984; Agúndez et al. 2021; Cernicharo et al. 2001, 2020, 2021a,b, 2024; Lee et al. 2021; McGuire et al. 2017). In the right part of the reaction scheme, we present the fully saturated aliphatic organic compounds. Most of these molecules have been observed in the comet 67P/Churyumov-Gerasimenko. The heavy molecules whose identifications are still lacking are marked with asterisks. It should be noted that UV- photolysis and cosmic ray processing of fully saturated aliphatic compounds result in their efficient dehydrogenation. The high saturation degree of the detected aliphatic compounds hints at their partial formation through surface hydrogenation mechanisms (Altwegg et al. 2017; Schuhmann et al. 2019; Hänni et al. 2021, 2022, 2023). The image of comet 67P/Churyumov-Gerasimenko is taken by ESA Rosetta on September 2, 20143.

Acknowledgements

This work is dedicated to the memory of Prof. Dr. Harold V.J. Linnartz. This research was funded through financial support from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie actions grant agreement no. 664931, and the Netherlands Research School for Astronomy (NOVA). GF thanks the Xinjiang Tianchi Talent Program (2024). GAB and MEP acknowledge Progetto Premiale 2012 ‘iALMA’ grant (CUP C52I13000140001) and INAF Grant “Ricerca Fondamentale 2022” (CUP C63C22000900005). K.J.C. is grateful for support from NWO via a VENI fellowship (VI.Veni.212.296). We would like to thanks Juan Tuo for her help in preparation of this manuscript.

Appendix A QMS-TPD spectra subtraction

Figures A.1, A.2, and A.3 show the results of subtraction between the mass spectra obtained prior and after hydrogenation of the ice for each of the individual m/z values presented through Figs. 2, 3, and 4.

thumbnail Fig. A.1

Results of subtraction for each of the m/z values presented in Fig. 2

thumbnail Fig. A.2

Results of subtraction for each of the m/z values presented in Fig. 3

thumbnail Fig. A.3

Results of subtraction for each of the m/z values presented in Fig. 4

Appendix B Fragmentation pattern comparison

In Fig. B.1, examples of comparison between the m/z signals registered at 105 and 138 K with the mass spectra available from the NIST database are presented. It is shown that the ratio between m/z≥43 signals demonstrates relatively good agreement with the NIST spectra of linear n-C4H10 and n-C6H14. However, co-desorption of entrapped C2H6 (or other species) along with the main constituents of the ice is required to account for the high relative intensity of m/z= 27, 29, 30 signals registered at these temperatures.

thumbnail Fig. B.1

Left: Comparison between the observed intensities of m/z signals at 105 K for 70 eV electron ionization energy of the QMS and the mass fragmentation patterns of C4H10 (CH3CH2CH2CH3), C4H8 (CH2CHCH2CH3), and C2H6 (CH3CH3) acquired from NIST database. Right: Similar comparison between intensities obtained at 138 K with the mass fragmentation patterns of C6H14 (CH3CH2CH2CH2CH2CH3) and C2H6 (CH3 CH3) acquired from the NIST database. Experimental values are baseline subtracted. Only the selected m/z values with resolved peaks are presented.

References

  1. Abplanalp, M. J., & Kaiser, R. I. 2017, ApJ, 836, 195 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abplanalp, M. J., & Kaiser, R. I. 2020, ApJ, 889, 3 [Google Scholar]
  3. Agúndez, M., Marcelino, N., Tercero, B., et al. 2021, A&A, 649, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Altwegg, K., Balsiger, H., Berthelier, J. J., et al. 2017, MNRAS, 469, S130 [Google Scholar]
  5. Belloche, A., Garrod, R. T., Zingsheim, O., Müller, H. S. P., & Menten, K. M. 2022, A&A, 662, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Bergner, J. B., Guzmán, V. G., Öberg, K. I., Loomis, R. A., & Pegues, J. 2018, ApJ, 857, 69 [Google Scholar]
  7. Boogert, A. C. A., Gerakines, P. A., & Whittet, D. C. B. 2015, ARA&A, 53, 541 [Google Scholar]
  8. Brown, R. L. 1981, ApJ, 248, L119 [NASA ADS] [CrossRef] [Google Scholar]
  9. Butscher, T., Duvernay, F., Theule, P., et al. 2015, MNRAS, 453, 1587 [NASA ADS] [CrossRef] [Google Scholar]
  10. Butscher, T., Duvernay, F., Danger, G., & Chiavassa, T. 2016, A&A, 593, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Butscher, T., Duvernay, F., Rimola, A., Segado-Centellas, M., & Chiavassa, T. 2017, Phys. Chem. Chem. Phys. (Incorp. Faraday Trans.), 19, 2857 [NASA ADS] [CrossRef] [Google Scholar]
  12. Cernicharo, J., Heras, A. M., Tielens, A. G. G. M., et al. 2001, ApJ, 546, L123 [Google Scholar]
  13. Cernicharo, J., Marcelino, N., Agúndez, M., et al. 2020, A&A, 642, L8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Cernicharo, J., Agúndez, M., Cabezas, C., et al. 2021a, A&A, 656, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Cernicharo, J., Cabezas, C., Agúndez, M., et al. 2021b, A&A, 648, L3 [EDP Sciences] [Google Scholar]
  16. Cernicharo, J., Fuentetaja, R., Cabezas, C., et al. 2022, A&A, 663, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cernicharo, J., Tercero, B., Agúndez, M., et al. 2024, A&A, 686, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Chen, Y. J., Chuang, K. J., Muñoz Caro, G. M., et al. 2014, ApJ, 781, 15 [Google Scholar]
  19. Chevance, M., Kruijssen, J. M. D., Vazquez-Semadeni, E., et al. 2020, Space Sci. Rev., 216, 50 [CrossRef] [Google Scholar]
  20. Chuang, K. J., Fedoseev, G., Qasim, D., et al. 2017, MNRAS, 467, 2552 [NASA ADS] [Google Scholar]
  21. Chuang, K. J., Fedoseev, G., Qasim, D., et al. 2018a, A&A, 617, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Chuang, K. J., Fedoseev, G., Qasim, D., et al. 2018b, ApJ, 853, 102 [NASA ADS] [CrossRef] [Google Scholar]
  23. Chuang, K. J., Fedoseev, G., Qasim, D., et al. 2020, A&A, 635, A199 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Chuang, K. J., Fedoseev, G., Scirè, C., et al. 2021, A&A, 650, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Chuang, K. J., Jäger, C., Santos, J. C., & Henning, T. 2024, A&A, 687, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Comeford, J., & Gould, J. H. 1961, J. Mol. Spectrosc., 5, 474 [Google Scholar]
  27. Compagnini, G., D’Urso, L., Puglisi, O., Baratta, G., & Strazzulla, G. 2009, Carbon, 47, 1605 [Google Scholar]
  28. Cuadrado, S., Goicoechea, J. R., Pilleri, P., et al. 2015, A&A, 575, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Cuylle, S. H., Zhao, D., Strazzulla, G., & Linnartz, H. 2014, A&A, 570, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. De Sanctis, M. C., Ammannito, E., McSween, H. Y., et al. 2017, Science, 355, 719 [NASA ADS] [CrossRef] [Google Scholar]
  31. Fatima, M., Müller, H. S. P., Zingsheim, O., et al. 2023, A&A, 680, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Fedoseev, G., Cuppen, H. M., Ioppolo, S., Lamberts, T., & Linnartz, H. 2015, MNRAS, 448, 1288 [NASA ADS] [CrossRef] [Google Scholar]
  33. Fedoseev, G., Chuang, K. J., van Dishoeck, E. F., Ioppolo, S., & Linnartz, H. 2016, MNRAS, 460, 4297 [Google Scholar]
  34. Fedoseev, G., Chuang, K. J., Ioppolo, S., et al. 2017, ApJ, 842, 52 [Google Scholar]
  35. Fedoseev, G., Qasim, D., Chuang, K.-J., et al. 2022, ApJ, 924, 110 [NASA ADS] [CrossRef] [Google Scholar]
  36. Fuentetaja, R., Agúndez, M., Cabezas, C., et al. 2022a, A&A, 667, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Fuentetaja, R., Cabezas, C., Agúndez, M., et al. 2022b, A&A, 663, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Gerakines, P. A., & Moore, M. H. 2001, Icarus, 154, 372 [NASA ADS] [CrossRef] [Google Scholar]
  39. Guzmán, V. V., Bergner, J. B., Law, C. J., et al. 2021, ApJS, 257, 6 [CrossRef] [Google Scholar]
  40. Hänni, N., Altwegg, K., Balsiger, H., et al. 2021, A&A, 647, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Hänni, N., Altwegg, K., Combi, M., et al. 2022, Nat. Commun., 13, 3639 [CrossRef] [Google Scholar]
  42. Hänni, N., Altwegg, K., Baklouti, D., et al. 2023, A&A, 678, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Hiraoka, K., Takayama, T., Euchi, A., Handa, H., & Sato, T. 2000, ApJ, 532, 1029 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hopkins, B. J., & Riviere, J. C. 1964, Br. J. Appl. Phys., 15, 941 [NASA ADS] [CrossRef] [Google Scholar]
  45. Hudson, R. L., Ferrante, R. F., & Moore, M. H. 2014, Icarus, 228, 276 [NASA ADS] [CrossRef] [Google Scholar]
  46. Ilee, J. D., Walsh, C., Booth, A. S., et al. 2021, ApJS, 257, 9 [NASA ADS] [CrossRef] [Google Scholar]
  47. Ioppolo, S., Fedoseev, G., Lamberts, T., Romanzin, C., & Linnartz, H. 2013, Rev. Sci. Instrum., 84, 073112 [NASA ADS] [CrossRef] [Google Scholar]
  48. Irvine, W. M., & Hjalmarson, A. 1984, Origins Life, 14, 15 [NASA ADS] [CrossRef] [Google Scholar]
  49. Jiménez-Serra, I., Rodríguez-Almeida, L. F., Martín-Pintado, J., et al. 2022, A&A, 663, A181 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Jo, S. K., & White, J. M. 1991, J. Chem. Phys., 94, 5761 [NASA ADS] [CrossRef] [Google Scholar]
  51. Khlifi, M., Paillous, P., Delpech, C., et al. 1995, J. Mol. Spectrosc., 174, 116 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kim, Y. S., & Kaiser, R. I. 2009, ApJS, 181, 543 [NASA ADS] [CrossRef] [Google Scholar]
  53. Kim, Y. S., Bennett, C. J., Chen, L.-H., O’Brien, K., & Kaiser, R. I. 2010, ApJ, 711, 744 [NASA ADS] [CrossRef] [Google Scholar]
  54. Kobayashi, H., Hidaka, H., Lamberts, T., et al. 2017, ApJ, 837, 155 [Google Scholar]
  55. Lee, K. L. K., Loomis, R. A., Burkhardt, A. M., et al. 2021, ApJ, 908, L11 [CrossRef] [Google Scholar]
  56. Ligterink, N. F. W., Paardekooper, D. M., Chuang, K. J., et al. 2015, A&A, 584, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Loomis, R. A., Burkhardt, A. M., Shingledecker, C. N., et al. 2021, Nat. Astron., 5, 188 [NASA ADS] [CrossRef] [Google Scholar]
  58. Maity, S., Kaiser, R. I., & Jones, B. M. 2015, Phys. Chem. Chem. Phys. (Incorp. Faraday Trans.), 17, 3081 [NASA ADS] [CrossRef] [Google Scholar]
  59. Martín, S., Mangum, J. G., Harada, N., et al. 2021, A&A, 656, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. McGuire, B. A., Burkhardt, A. M., Shingledecker, C. N., et al. 2017, ApJ, 843, L28 [NASA ADS] [CrossRef] [Google Scholar]
  61. Minissale, M., Moudens, A., Baouche, S., Chaabouni, H., & Dulieu, F. 2016, MNRAS, 458, 2953 [Google Scholar]
  62. Modica, P., & Palumbo, M. E. 2010, A&A, 519, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Molpeceres, G., & Rivilla, V. M. 2022, A&A, 665, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Nguyen, T., Fourré, I., Favre, C., et al. 2019, A&A, 628, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Oba, Y., Tomaru, T., Lamberts, T., Kouchi, A., & Watanabe, N. 2018, Nat. Astron., 2, 228 [NASA ADS] [CrossRef] [Google Scholar]
  66. Öberg, K. I., Garrod, R. T., van Dishoeck, E. F., & Linnartz, H. 2009, A&A, 504, 891 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Öberg, K. I., Boogert, A. C. A., Pontoppidan, K. M., et al. 2011, ApJ, 740, 109 [Google Scholar]
  68. Pardo, J. R., Cernicharo, J., & Goicoechea, J. R. 2005, ApJ, 628, 275 [NASA ADS] [CrossRef] [Google Scholar]
  69. Pardo, J. R., Cernicharo, J., Tercero, B., et al. 2022, A&A, 658, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Puglisi, O., Compagnini, G., D’Urso, L., et al. 2014, Nucl. Instrum. Methods Phys. Res. B, 326, 2 [Google Scholar]
  71. Qasim, D., Fedoseev, G., Chuang, K. J., et al. 2019a, A&A, 627, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Qasim, D., Fedoseev, G., Lamberts, T., et al. 2019b, ACS Earth Space Chem., 3, 986 [NASA ADS] [CrossRef] [Google Scholar]
  73. Rubin, M., Altwegg, K., Balsiger, H., et al. 2015, Science, 348, 232 [NASA ADS] [CrossRef] [Google Scholar]
  74. Sakai, N., & Yamamoto, S. 2013, Chem. Rev., 113, 8981 [Google Scholar]
  75. Santos, J. C., Enrique-Romero, J., Lamberts, T., Linnartz, H., & Chuang, K.-J. 2024, ACS Earth Space Chem., 8, 1646 [NASA ADS] [CrossRef] [Google Scholar]
  76. Schuhmann, M., Altwegg, K., Balsiger, H., et al. 2019, A&A, 630, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Segré, D., Ben-Eli, D., Deamer, D. W., & Lancet, D. 2001, Origins of Life and Evolution of the Biosphere, 31, 119 [CrossRef] [Google Scholar]
  78. Stevenson, J., Lunine, J., & Clancy, P. 2015, Sci. Adv., 1, 1400067 [NASA ADS] [CrossRef] [Google Scholar]
  79. Strazzulla, G., & Baratta, G. A. 1992, A&A, 266, 434 [NASA ADS] [Google Scholar]
  80. Taniguchi, K., Saito, M., & Ozeki, H. 2016, ApJ, 830, 106 [NASA ADS] [CrossRef] [Google Scholar]
  81. Taniguchi, K., Ozeki, H., & Saito, M. 2017, ApJ, 846, 46 [NASA ADS] [CrossRef] [Google Scholar]
  82. Taniguchi, K., Herbst, E., Ozeki, H., & Saito, M. 2019, ApJ, 884, 167 [Google Scholar]
  83. Taniguchi, K., Gorai, P., & Tan, J. C. 2024, Ap&SS, 369, 34 [NASA ADS] [CrossRef] [Google Scholar]
  84. Theule, P., Borget, F., Mispelaer, F., et al. 2011, A&A, 534, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Tørneng, E., Nielsen, C. J., Klaeboe, P., Hopf, H., & Priebe, H. 1980, Spectroch. Acta A: Mol. Spectrosc., 36, 975 [CrossRef] [Google Scholar]
  86. Trevors, J. T., & Psenner, R. 2001, FEMS Microbiol. Rev., 25, 573 [CrossRef] [Google Scholar]
  87. Tschersich, K. G. 2000, J. Appl. Phys., 87, 2565 [NASA ADS] [CrossRef] [Google Scholar]
  88. Turner, A. M., Abplanalp, M. J., Blair, T. J., Dayuha, R., & Kaiser, R. I. 2018, ApJS, 234, 6 [NASA ADS] [CrossRef] [Google Scholar]
  89. Urso, R. G., Palumbo, M. E., Ceccarelli, C., et al. 2019, A&A, 628, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. van Dishoeck, E. F., Grant, S., Tabone, B., et al. 2023, Faraday Discuss., 245, 52 [NASA ADS] [CrossRef] [Google Scholar]
  91. Wu, Y.-J., Lin, M.-Y., Chou, S.-L., et al. 2010, ApJ, 721, 856 [NASA ADS] [CrossRef] [Google Scholar]
  92. Yabuta, H., Cody, G. D., Engrand, C., et al. 2023, Science, 379, abn9057 [NASA ADS] [CrossRef] [Google Scholar]
  93. Zhou, L., Kaiser, R. I., & Tokunaga, A. T. 2009, Planet. Space Sci., 57, 830 [Google Scholar]
  94. Zhou, L., Zheng, W., Kaiser, R. I., et al. 2010, ApJ, 718, 1243 [NASA ADS] [CrossRef] [Google Scholar]

1

1 Langmuir (L) corresponds to exposure of substrate at 1 × 10−6 Torr pressure during one second at room temperature. This results in a surface coverage of about 1 monolayer. The usage of Langmuir estimation for the Ar ice is caused by the lack of absorption in the IR range.

2

NIST Mass Spec Data Center, S. E. Stein, director, “Mass Spectra” in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, eds. P. J. Linstrom and W. G. Mallard, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA.

3

Emily Lakdawalla, Comet jets! NavCam view of comet Churyumov- Gerasimenko on September 2, 2014, ESA, Rosetta, https://www.planetary.org/space-images/comet-cg-20140902-jets

All Figures

thumbnail Fig. 1

Left: QMS-TPD spectra in the range from 50 to 225 K for the selected m/z values obtained after the UV-exposure of 10 ML of pure HCCH ice with 4.6 × 1016 photons cm−2 at 10 K. Right: QMS-TPD spectrum obtained after the hydrogenation of UV-exposed ice with 1.3 × 1017 cm−2 atoms at 10 K. The ordinate axis is presented as a logarithm of 10. An exponential decay of ion intensity appears as a straight line.

In the text
thumbnail Fig. 2

Direct comparison between the fragments of QMS-TPD spectra in the range from 50 to 225 K obtained after the UV-exposure of 10 ML thick pure HCCH ice at 10 K with >4.6 × 1016 cm−2 photons (black empty circles) and the subsequent exposure of this photoprocessed ice with 1.3 × 1017 cm−2 H atoms at 10 K (red filled circles). The results obtained upon subtraction of the initial mass spectrum from the mass spectrum obtained upon hydrogenation of the ice for each of the presented m/z values are presented in Appendix A.

In the text
thumbnail Fig. 3

Similar to Fig. 2. M/z signals corresponding to linear saturated aliphatic hydrocarbons are shown.

In the text
thumbnail Fig. 4

Similar to Figs. 2 and 3. m/z signals associated with the semisaturated aliphatic hydrocarbons are presented.

In the text
thumbnail Fig. 5

Left: RAIR spectra obtained after the deposition of 10 ML of HCCH ice at 10 K (a, black line), after the exposure of this obtained ice with 4.6 × 1016 cm−2 photons (b, red line) and by subtracting the initial spectrum (a) from spectrum (b) (green line). Right: RAIR spectra obtained after the UV-exposure of 10 ML thick pure HCCH ice at 10 K with 4.6 × 1016 cm−2 photons (a, red line), the consequent exposure of this photoprocessed ice with 1.3 × 1017 cm−2 H atoms at 10 K (b, blue line), and by subtracting the initial spectrum (a) from spectrum (b) (green line). The grey dashed and dotted lines are the reference RAIRS spectra obtained for the pure C4H2 and C4H4 ices by Zhou et al. (2009) and Kim & Kaiser (2009). All spectra are offset for clarity.

In the text
thumbnail Fig. 6

Zoom-in into the CH stretching vibration modes region for the RAIR difference spectrum obtained after hydrogenation of UV-exposed C2H2 ice (green line). See the right panel of Fig. 5 for more details. The spectrum of pure C2H6 ice at 20 K is adapted from Öberg et al. (2009) (black dots). The spectrum of pure C4H10 at 10 K is obtained by Strazzulla & Baratta (1992) (blue line). All spectra are offset for clarity.

In the text
thumbnail Fig. 7

Kinetic curves obtained for C2H2, C6H2, C4H2, and C4H10. Note: the area of the (ν4) absorption feature of C4H2 cannot be deconvoluted from the contribution of the (ν1) absorption feature of C4H4. Thus, the combined area is used for the integration. The modulus of absorbance is taken for C2H2 and C4H2 consumed in the experiments. This allows for the direct comparison with produced species. The ordinate axis is presented as a logarithm of 10.

In the text
thumbnail Fig. 8

Schematic representation of the overall suggested mechanism. Various polyynes observed in dark clouds (e.g., TMC-1) are presented in the left part of the reaction scheme. The accretion of polyynes on the surface of the grains is followed by hydrogenation of their unsaturated CC triple bonds, resulting in the formation of their semi-saturated and fully saturated counterparts. Besides CH2 CHCCCHO, all of the first- generation products obtained by such hydrogenation have also been observed in the gas phase in TMC-1 or other dark clouds (Brown 1981; Irvine & Hjalmarson 1984; Agúndez et al. 2021; Cernicharo et al. 2001, 2020, 2021a,b, 2024; Lee et al. 2021; McGuire et al. 2017). In the right part of the reaction scheme, we present the fully saturated aliphatic organic compounds. Most of these molecules have been observed in the comet 67P/Churyumov-Gerasimenko. The heavy molecules whose identifications are still lacking are marked with asterisks. It should be noted that UV- photolysis and cosmic ray processing of fully saturated aliphatic compounds result in their efficient dehydrogenation. The high saturation degree of the detected aliphatic compounds hints at their partial formation through surface hydrogenation mechanisms (Altwegg et al. 2017; Schuhmann et al. 2019; Hänni et al. 2021, 2022, 2023). The image of comet 67P/Churyumov-Gerasimenko is taken by ESA Rosetta on September 2, 20143.

In the text
thumbnail Fig. A.1

Results of subtraction for each of the m/z values presented in Fig. 2

In the text
thumbnail Fig. A.2

Results of subtraction for each of the m/z values presented in Fig. 3

In the text
thumbnail Fig. A.3

Results of subtraction for each of the m/z values presented in Fig. 4

In the text
thumbnail Fig. B.1

Left: Comparison between the observed intensities of m/z signals at 105 K for 70 eV electron ionization energy of the QMS and the mass fragmentation patterns of C4H10 (CH3CH2CH2CH3), C4H8 (CH2CHCH2CH3), and C2H6 (CH3CH3) acquired from NIST database. Right: Similar comparison between intensities obtained at 138 K with the mass fragmentation patterns of C6H14 (CH3CH2CH2CH2CH2CH3) and C2H6 (CH3 CH3) acquired from the NIST database. Experimental values are baseline subtracted. Only the selected m/z values with resolved peaks are presented.

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.