pith. sign in

arxiv: 2605.06767 · v1 · submitted 2026-05-07 · 🌌 astro-ph.GA

Complex organic molecules and cosmic ray ionisation rate towards the massive protostar Cepheus A HW2

Pith reviewed 2026-05-11 01:05 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords cosmic ray ionizationcomplex organic moleculeshigh-mass protostarCepheus A HW2molecular spectroscopygas-grain chemistryspectral survey
0
0 comments X

The pith

Observations show locally enhanced cosmic ray ionization rates in the complex organic molecule emitting gas towards Cepheus A HW2.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper reports a wide-band spectral survey detecting ions and complex organic molecules such as methanol, methyl cyanide, and methyl formate towards the high-mass protostar Cepheus A HW2. Column densities, abundances, deuterium fractions, and molecular hydrogen density are derived for two kinematic components using rotational diagrams and radiative transfer modeling. Analytic chemistry applied to ion abundances yields a cosmic ray ionization rate of 6.8×10^{-17} s^{-1} in the -11 km/s component that dominates COM emission, compared to an upper limit of 9.2×10^{-19} s^{-1} in the -5 km/s component. The Nautilus gas-grain code reproduces the observed abundances when run at the measured density, temperature, and ionization rate. This finding matters because cosmic rays are known to drive molecular chemistry and release organics from grain surfaces, yet their role had mainly been tested in cold low-mass cores rather than high-mass environments.

Core claim

The central claim is that the cosmic ray ionization rate of the kinematic component associated with most of the COMs' emission in the region is locally enhanced. This is shown by deriving CRIR values from observed ion abundances via analytic chemistry expressions, with the -11 km/s component at 6.8×10^{-17} s^{-1} while the -5 km/s component is at most 9.2×10^{-19} s^{-1}. The Nautilus chemical model matches the measured abundances of CH3OH, CH3CN, HCO+, and N2H+ at the observed density of 2.6×10^5 cm^{-3}, temperature, and the higher CRIR.

What carries the argument

Analytic chemistry expressions that convert observed abundances of ions such as HCO+ and N2H+ into cosmic ray ionization rates, validated by comparison to abundances predicted by the Nautilus gas-grain chemical code.

If this is right

  • The Nautilus model reproduces the observed abundances of CH3OH, CH3CN, HCO+, and N2H+ at the derived density, temperature, and CRIR within model uncertainty.
  • Deuterium fractions in the range 0.002-0.3 are measured in the -11 km/s component along with a molecular hydrogen density of 2.6×10^5 cm^{-3}.
  • Abundances of additional COMs including t-HCOOH, H2CCO, CH3CHO, and CH3OCHO are reported for the main emitting component.
  • The contrast in CRIR between components points to spatially varying ionization linked to the presence of COM emission.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The enhancement may result from the ionized jet accelerating cosmic rays locally, an effect that uniform-CRIR models of high-mass regions would miss.
  • Higher-resolution mapping could test whether the CRIR increase correlates directly with jet position or with COM column density peaks.
  • If the pattern holds in other high-mass protostars, chemical models should incorporate locally boosted ionization when predicting organic inventories near jets.

Load-bearing premise

The analytic chemistry expressions that convert observed ion abundances into a cosmic-ray ionization rate accurately reflect the dominant ionization processes and are not significantly affected by the nearby ionized jet or other unmodeled sources.

What would settle it

An independent measurement of the cosmic ray ionization rate, for instance via gamma-ray observations or additional ion ratios unaffected by the jet, that yields values inconsistent with 6.8×10^{-17} s^{-1} in the -11 km/s component and ≤9.2×10^{-19} s^{-1} in the -5 km/s component.

Figures

Figures reproduced from arXiv: 2605.06767 by Anna Punanova, A. O. Henrik Olofsson, Brandt Gaches, Emma W. Nielsen, Eva Wirstr\"om, Jonathan C. Tan, Paola Caselli, Prasanta Gorai.

Figure 1
Figure 1. Figure 1: Gaussian fits to the integrated intensities observed at different positions along RA and Dec axes. The (0,0) position corresponds to Cep A HW2. age pyspeckit (Ginsburg & Mirocha 2011; Ginsburg et al. 2022). Assuming the excitation temperature is the same for all hyper￾fine components, pyspeckit calculates the HFS from the relative velocities and intensities of the hyperfine components. The HFS fits constra… view at source ↗
Figure 2
Figure 2. Figure 2: Rotational diagrams at VLSR = −11 km s−1 (left) and −5 km s−1 (right). Top: Methanol: CH3OH-E (blue dots) and CH3OH-A (red crosses). Bottom: Methyl cyanide: p-CH3CN (blue dots) and o￾CH3CN (red crosses). densities of methanol A may be overestimated2 . Because of the non-thermal nature of methanol excitation and to probe volume density and kinetic temperature, we also modelled N(CH3OH) with RADEX and used t… view at source ↗
Figure 3
Figure 3. Figure 3: Rotational diagrams at VLSR = −11 km s−1 . Top: Acetaldehyde: CH3CHO-A (blue dots) and CH3CHO-E (red crosses). Individual fits (blue dashed and red dotted lines for CH3CHO-A and CH3CHO-E, re￾spectively) and a combined fit (black line) are shown. Bottom: Ketene (H2CCO) (left) and trans-formic acid (t-HCOOH) (right). 3.2.2. Optically thin/thick transitions When only one transition was detected for a given mo… view at source ↗
Figure 4
Figure 4. Figure 4: RADEX grid search using low-energy (Eu < 50 K) methanol E lines at −11 km s−1 . The coloured lines (and 1σ shaded regions) show what combinations of parameters are consistent with each of the measured integrated intensities. The grey regions show the 1, 2 and 3σ regions of the fitted parameters in the plane of the best fit Tkin, n(H2) and N(CH3OH) (from left to right). The excitation temperatures for the l… view at source ↗
Figure 5
Figure 5. Figure 5: RADEX grid search using low-energy (Eu < 50 K) methanol E/A lines at −5 km s−1 . Coloured lines and shaded areas are as in [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: RADEX grid search using high-energy (Eu > 50 K) methanol E/A lines at −5 km s−1 (excluding the line at 95.17 GHz). Coloured lines and shaded areas are as in [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Column densities of deuterated against hydrogenated species. whereas Luo et al. (2024) consider the formation of N2H + and H3O + . Caselli et al. (1998) trace this chemistry from measure￾ments of HCO+ , DCO+ and CO, and Luo et al. (2024) use N2H + and HCO+ as tracers. We inferred N(HCO+ ) from the measured HC18O + (1–0) transition and N(CO) from the measured C18O (1–0) transition (see Table B.2). The Casel… view at source ↗
Figure 8
Figure 8. Figure 8: CRIRs (left) and electron fractions (right) for Cep A HW2 (at VLSR = −11 km s−1 ) compared to LMCs, IRDCs, HMCs, HMPOs and UCHIIs. Medians, minimum and maximum values are shown with black lines; 16% and 84% percentiles are shown with dotted black lines. The CRIR measurements for HMCs, HMPOs and UCHIIs are from Luo et al. (2024). The measurements for Cep A HW2 were calculated using the method from Luo et al… view at source ↗
Figure 9
Figure 9. Figure 9: Abundance of COMs for Cep A HW2 (at VLSR = −11 km s−1 ) compared to LMCs, IRDCs, HMCs, HMPOs and UCHIIs. Medians, minimum and maximum values are shown with black lines. If fewer than 5 measurements are reported, the individual measurements are shown instead of a violin plot. The measurements for IRDCs, HMCs, HMPOs and UCHIIs are from Nummelin et al. (1998); Ikeda et al. (2001); Gerner et al. (2014); Vasyun… view at source ↗
Figure 10
Figure 10. Figure 10: Top and middle: 0D Chemical model abundances for HCO+ , N2H + , CH3OH, and CH3CN as a function of the H2 CRIR. The hori￾zontal dashed lines indicate the observed inferred abundances. The solid lines with bands show the model results with a 1-dex spread. Bottom: Abundance ratio of CH3OH and CH3CN as a function of H2 CRIR. The horizontal dashed cyan line shows the observed abundance ratio, and the black lin… view at source ↗
Figure 11
Figure 11. Figure 11: Deuterium fractions for Cep A HW2 (at VLSR = −11 km s−1 ) compared to LMCs, IRDCs, HMCs, HMPOs and UCHIIs. Medians, minimum and maximum values are shown with black lines. The measurements for IRDCs, HMCs, HMPOs and UCHIIs are from Pillai et al. (2007); Miettinen et al. (2011); Fontani et al. (2011, 2014, 2015); Gerner et al. (2015). The measurements for LMCs are from Hatchell (2003); Roueff & Gerin (2003)… view at source ↗
read the original abstract

Cosmic rays (CRs) are important drivers for molecular chemistry in star-forming regions, and laboratory experiments have shown that CRs can stimulate the release of complex organic molecules (COMs) such as methanol. Observationally, this has primarily been tested in cold, low-mass cores, so studying how CRs affect COM formation in a high-mass star-forming environment is of great interest. We performed a high-sensitivity wide-band spectral line survey with the Onsala 20 m telescope towards the high-mass protostar Cepheus A HW2, which is known to host an ionised jet. Consistent with previous studies, two primary velocity components ($-11$ km s$^{-1}$ and $-5$ km s$^{-1}$) were identified. Column densities and relative abundances of the detected ions and COMs were estimated from rotational diagrams, single transitions and RADEX grid searches (CH$_3$OH: $1.6\times10^{-9}$, CH$_3$CN: $5.9\times10^{-11}$, t-HCOOH: $7.9\times10^{-11}$, H$_2$CCO: $1.7\times10^{-11}$, CH$_3$CHO: $1.9\times10^{-11}$, CH$_3$OCHO: $7.6\times10^{-10}$ at $-11$ km s$^{-1}$). Deuterium fractions were also estimated (in range $0.002-0.3$ at $-11$ km s$^{-1}$), and the volume density of molecular hydrogen ($2.6\times10^5$ cm$^{-3}$ at $-11$ km s$^{-1}$) was constrained from the RADEX grid searches. Electron fractions and CR ionisation rates (CRIR, $6.8\times10^{-17}$ s$^{-1}$ at $-11$ km s$^{-1}$, $\leq9.2\times10^{-19}$ s$^{-1}$ at $-5$ km s$^{-1}$) were estimated through analytic chemistry using different ions as probes. The gas-grain chemical code Nautilus reproduced the observed abundances of CH$_3$OH, CH$_3$CN, HCO$^+$, N$_2$H$^+$ at the observed density, temperature and CRIR within the uncertainty of the model. The results indicate that the CR ionisation rate of the kinematic component associated with most of the COMs' emission in the region is locally enhanced.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript reports results from a wide-band spectral line survey of the high-mass protostar Cepheus A HW2 with the Onsala 20 m telescope. Two kinematic components are identified at approximately -11 km s^{-1} and -5 km s^{-1}. Column densities and abundances of detected COMs (e.g., CH3OH at 1.6e-9, CH3CN at 5.9e-11) and ions are derived using rotational diagrams, single-line estimates, and RADEX non-LTE grids; H2 volume density is constrained at 2.6e5 cm^{-3} for the -11 km s^{-1} component. Electron fractions and cosmic-ray ionization rates (CRIR) are obtained via analytic ion chemistry (CRIR = 6.8e-17 s^{-1} at -11 km s^{-1} versus an upper limit of 9.2e-19 s^{-1} at -5 km s^{-1}). The Nautilus gas-grain code is then run at the observed density, temperature, and derived CRIR to reproduce the abundances of CH3OH, CH3CN, HCO+, and N2H+. The central claim is that the CRIR is locally enhanced in the velocity component associated with the majority of the COM emission.

Significance. If the analytic CRIR derivation is shown to be robust, the work would extend studies of CR-driven COM chemistry from cold low-mass cores to a high-mass protostellar environment and provide a quantitative link between elevated CRIR and COM abundances. The multi-method approach (rotational diagrams + RADEX + analytic chemistry + Nautilus validation) and the successful reproduction of key abundances within model uncertainties are positive features that support the internal consistency of the analysis.

major comments (1)
  1. [CRIR estimation via analytic ion chemistry (methods and results sections)] The analytic chemistry expressions used to convert observed HCO+, N2H+, and electron abundances into CRIR values (reported as 6.8×10^{-17} s^{-1} at -11 km s^{-1}) assume cosmic-ray ionization is the dominant process balancing recombination and ion-molecule reactions. The source is known to host an ionized jet; any UV or X-ray leakage from the jet could produce the same ion enhancements without an increase in true CR flux. The Nautilus runs fix the CRIR at the analytically derived value and do not test whether an additional non-CR ionization term at a canonical CRIR can reproduce the observed abundances, which directly affects the load-bearing claim of local CRIR enhancement in the -11 km s^{-1} component.
minor comments (2)
  1. Explicit 1σ uncertainties or error bars should be reported on all derived quantities (column densities, abundances, CRIR values, and deuterium fractions) so that the statistical significance of the difference between the two velocity components can be evaluated.
  2. A supplementary table listing all detected transitions, integrated intensities, line widths, and the specific method used for each column-density estimate would improve reproducibility and allow independent verification of the rotational-diagram and RADEX results.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We appreciate the referee's detailed and constructive comments on our manuscript. The feedback highlights an important assumption in our analysis of the cosmic ray ionization rate. Below we provide a point-by-point response to the major comment and indicate the revisions we will make to the manuscript.

read point-by-point responses
  1. Referee: The analytic chemistry expressions used to convert observed HCO+, N2H+, and electron abundances into CRIR values (reported as 6.8×10^{-17} s^{-1} at -11 km s^{-1}) assume cosmic-ray ionization is the dominant process balancing recombination and ion-molecule reactions. The source is known to host an ionized jet; any UV or X-ray leakage from the jet could produce the same ion enhancements without an increase in true CR flux. The Nautilus runs fix the CRIR at the analytically derived value and do not test whether an additional non-CR ionization term at a canonical CRIR can reproduce the observed abundances, which directly affects the load-bearing claim of local CRIR enhancement in the -11 km s^{-1} component.

    Authors: We thank the referee for this insightful comment. The analytic expressions do indeed assume that cosmic-ray ionization is the primary driver. We have revised the manuscript to state this assumption more clearly and to discuss the ionized jet as a potential source of UV or X-ray ionization. We note that the localization of the high CRIR and COM emission to one kinematic component supports a local process rather than widespread jet leakage. Regarding the chemical models, the Nautilus runs were intended to check consistency with the derived CRIR. We have added text acknowledging that models with canonical CRIR plus extra ionization were not explored in the current work and suggest this as an avenue for future investigation. This does not change our main conclusion but highlights a limitation in distinguishing ionization sources. revision: partial

Circularity Check

0 steps flagged

No circularity: CRIR from standard analytic ion chemistry; Nautilus is post-hoc consistency check

full rationale

The paper derives CRIR values (6.8e-17 s^-1 and upper limit 9.2e-19 s^-1) directly from observed ion abundances via standard analytic chemistry expressions that balance cosmic-ray ionization against recombination and ion-molecule reactions. These expressions are applied separately to the two velocity components using HCO+, N2H+ and electron fractions. Nautilus is then run forward at the fixed observed density, temperature and this independently-derived CRIR to verify that CH3OH, CH3CN, HCO+ and N2H+ abundances fall within model uncertainty; it is not used to fit or back-solve for CRIR. No self-citation chain, ansatz smuggling, or renaming of known results is present in the derivation steps. The central claim therefore rests on external analytic formulae applied to new observations rather than reducing to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard but unverified assumptions in observational astrochemistry techniques and chemical modeling; no new physical entities are postulated.

axioms (2)
  • domain assumption Steady-state assumptions in analytic chemistry relating ion abundances to cosmic-ray ionization rate
    Invoked to derive CRIR from ions such as HCO+ and N2H+.
  • domain assumption Validity of RADEX non-LTE modeling for constraining H2 density and molecular abundances from observed lines
    Used to obtain the volume density 2.6e5 cm^-3 and COM abundances.

pith-pipeline@v0.9.0 · 5798 in / 1474 out tokens · 54624 ms · 2026-05-11T01:05:05.917423+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

103 extracted references · 103 canonical work pages

  1. [1]

    2021, , 649, L4

    Ag \'u ndez , M., Marcelino , N., Tercero , B., et al. 2021, , 649, L4

  2. [2]

    & Faure , A

    Bacmann , A. & Faure , A. 2016, , 587, A130

  3. [3]

    2012, , 541, L12

    Bacmann , A., Taquet , V., Faure , A., Kahane , C., & Ceccarelli , C. 2012, , 541, L12

  4. [4]

    Bettens , R. P. A., Hansen , T. A., & Collins , M. A. 1999, , 111, 6322

  5. [5]

    Bobylev, V. V. & Bajkova, A. T. 2010, Monthly Notices of the Royal Astronomical Society, 408, 1788–1795

  6. [6]

    L., Chandler , C

    Brogan , C. L., Chandler , C. J., Hunter , T. R., Shirley , Y. L., & Sarma , A. P. 2007, , 660, L133

  7. [7]

    M., et al

    Cabedo , V., Maury , A., Girart , J. M., et al. 2023, , 669, A90

  8. [8]

    M., Tercero , B., Berm \'u dez , C., & Cernicharo , J

    Cabezas , C., Neeman , E. M., Tercero , B., Berm \'u dez , C., & Cernicharo , J. 2021, , 645, A75

  9. [9]

    2025, , 698, L24

    Cabezas , C., V \'a vra , K., Molpeceres , G., et al. 2025, , 698, L24

  10. [10]

    2008, , 492, 703

    Caselli , P., Vastel , C., Ceccarelli , C., et al. 2008, , 492, 703

  11. [11]

    M., Terzieva, R., & Herbst, E

    Caselli, P., Walmsley, C. M., Terzieva, R., & Herbst, E. 1998, The Astrophysical Journal, 499, 234

  12. [12]

    M., Zucconi, A., et al

    Caselli, P., Walmsley, C. M., Zucconi, A., et al. 2002, The Astrophysical Journal, 565, 331–343

  13. [13]

    2025, pythonradex

    Cataldi, G. 2025, pythonradex

  14. [14]

    2019, , 622, A141

    Chac \'o n-Tanarro , A., Caselli , P., Bizzocchi , L., et al. 2019, , 622, A141

  15. [15]

    M., et al

    Crapsi , A., Caselli , P., Walmsley , C. M., et al. 2005, , 619, 379

  16. [16]

    C., & Tafalla , M

    Crapsi , A., Caselli , P., Walmsley , M. C., & Tafalla , M. 2007, , 470, 221

  17. [17]

    Dagdigian , P. J. 2024, , 527, 2209

  18. [18]

    2020, , 634, A103

    Dartois , E., Chabot , M., Bacmann , A., et al. 2020, , 634, A103

  19. [19]

    M., Liu , M., Tan , J

    De Buizer , J. M., Liu , M., Tan , J. C., et al. 2017, , 843, 33

  20. [20]

    Draine , B. T. 1978, , 36, 595

  21. [21]

    P., Schlemmer , S., Schilke , P., Stutzki , J., & M \"u ller , H

    Endres , C. P., Schlemmer , S., Schilke , P., Stutzki , J., & M \"u ller , H. S. P. 2016, Journal of Molecular Spectroscopy, 327, 95

  22. [22]

    C., Cosentino , G., et al

    Entekhabi , N., Tan , J. C., Cosentino , G., et al. 2022, , 662, A39

  23. [23]

    D., Cloutier , P., Sanche , L., & Huels , M

    Esmaili , S., Bass , A. D., Cloutier , P., Sanche , L., & Huels , M. A. 2018, , 148, 164702

  24. [24]

    C., Telkamp , Z., et al

    Fedriani , R., Tan , J. C., Telkamp , Z., et al. 2023, , 942, 7

  25. [25]

    2015, , 575, A87

    Fontani , F., Busquet , G., Palau , A., et al. 2015, , 575, A87

  26. [26]

    2011, , 529, L7

    Fontani , F., Palau , A., Caselli , P., et al. 2011, , 529, L7

  27. [27]

    2014, , 440, 448

    Fontani , F., Sakai , T., Furuya , K., et al. 2014, , 440, 448

  28. [28]

    2015, , 584, A124

    Furuya , K., Aikawa , Y., Hincelin , U., et al. 2015, , 584, A124

  29. [29]

    Gaches, B. A. L. & Viti, S. 2026, ACS Earth and Space Chemistry

  30. [30]

    T., Wakelam , V., & Herbst , E

    Garrod , R. T., Wakelam , V., & Herbst , E. 2007, , 467, 1103

  31. [31]

    2014, , 563, A97

    Gerner , T., Beuther , H., Semenov , D., et al. 2014, , 563, A97

  32. [32]

    L., Beuther , H., et al

    Gerner , T., Shirley , Y. L., Beuther , H., et al. 2015, , 579, A80

  33. [33]

    2023, , 676, A78

    Giers , K., Spezzano , S., Caselli , P., et al. 2023, , 676, A78

  34. [34]

    & Mirocha , J

    Ginsburg , A. & Mirocha , J. 2011, PySpecKit: Python Spectroscopic Toolkit , Astrophysics Source Code Library, record ascl:1109.001

  35. [35]

    2022, The Astronomical Journal, 163, 291

    Ginsburg, A., Sokolov, V., de Val-Borro, M., et al. 2022, The Astronomical Journal, 163, 291

  36. [36]

    A., Pipher , J

    Goetz , J. A., Pipher , J. L., Forrest , W. J., et al. 1998, , 504, 359

  37. [37]

    Goldsmith, P. F. & Langer, W. D. 1999, The Astrophysical Journal, 517, 209

  38. [38]

    2021, , 506, 3734

    Guti \'e rrez-Quintanilla , A., Layssac , Y., Butscher , T., et al. 2021, , 506, 3734

  39. [39]

    D., & van Dishoeck , E

    Hacar , A., Bosman , A. D., & van Dishoeck , E. F. 2020, , 635, A4

  40. [40]

    2017, , 600, A61

    Harju, J., Daniel, F., Sipilä, O., et al. 2017, , 600, A61

  41. [41]

    2024, , 688, A117

    Harju , J., Vastel , C., Sipil \"a , O., et al. 2024, , 688, A117

  42. [42]

    2003, , 403, L25

    Hatchell , J. 2003, , 403, L25

  43. [43]

    & Klemperer , W

    Herbst , E. & Klemperer , W. 1973, , 185, 505

  44. [44]

    V., Pety , J., et al

    Hern \'a ndez-Vera , C., Guzm \'a n , V. V., Pety , J., et al. 2026, arXiv e-prints, arXiv:2601.20146

  45. [45]

    Hughes , V. A. & Wouterloot , J. G. A. 1984, , 276, 204

  46. [46]

    L., Henkel , C., & Mauersberger , R

    Huttemeister , S., Wilson , T. L., Henkel , C., & Mauersberger , R. 1993, , 276, 445

  47. [47]

    2001, , 560, 792

    Ikeda , M., Ohishi , M., Nummelin , A., et al. 2001, , 560, 792

  48. [48]

    V., Pellegrin , T., et al

    Indriolo , N., Ivlev , A. V., Pellegrin , T., et al. 2026, , 997, 123

  49. [49]

    I., Caselli , P., et al

    Jim \'e nez-Serra , I., Vasyunin , A. I., Caselli , P., et al. 2016, , 830, L6

  50. [50]

    I., Spezzano , S., et al

    Jim \'e nez-Serra , I., Vasyunin , A. I., Spezzano , S., et al. 2021, , 917, 44

  51. [51]

    C., Wakelam , V., & Sipil \"a , O

    Kong , S., Caselli , P., Tan , J. C., Wakelam , V., & Sipil \"a , O. 2015, , 804, 98

  52. [52]

    H., Knacke , R., Geballe , T

    Lacy , J. H., Knacke , R., Geballe , T. R., & Tokunaga , A. T. 1994, , 428, L69

  53. [53]

    I., et al

    Lattanzi , V., Bizzocchi , L., Vasyunin , A. I., et al. 2020, , 633, A118

  54. [54]

    2017, , 466, 248

    Li , S., Wang , J., Zhang , Z.-Y., et al. 2017, , 466, 248

  55. [55]

    L., Draine , B

    Linsky , J. L., Draine , B. T., Moos , H. W., et al. 2006, , 647, 1106

  56. [56]

    G., Padovani, M., & Gaches, B

    Luo, G., Bisbas, T. G., Padovani, M., & Gaches, B. A. L. 2024, , 690, A293

  57. [57]

    2023, , 519, 1601

    Meg \' as , A., Jim \'e nez-Serra , I., Mart \' n-Pintado , J., et al. 2023, , 519, 1601

  58. [58]

    2011, , 534, A134

    Miettinen , O., Hennemann , M., & Linz , H. 2011, , 534, A134

  59. [59]

    J., Menten , K

    Moscadelli , L., Reid , M. J., Menten , K. M., et al. 2009, , 693, 406

  60. [60]

    M \"u ller , H. S. P., Thorwirth , S., Roth , D. A., & Winnewisser , G. 2001, , 370, L49

  61. [61]

    Müller, H. S. P., Menten, K. M., & Mäder, H. 2004, , 428, 1019–1026

  62. [62]

    2019, , 630, A136

    Nagy , Z., Spezzano , S., Caselli , P., et al. 2019, , 630, A136

  63. [63]

    E., Bergman , P., et al

    Nummelin , A., Dickens , J. E., Bergman , P., et al. 1998, , 337, 275

  64. [64]

    Oberg , K. I. 2016, arXiv e-prints, arXiv:1609.03112

  65. [65]

    2016, , 590, A8

    Padovani , M., Marcowith , A., Hennebelle , P., & Ferri \`e re , K. 2016, , 590, A8

  66. [66]

    G., Zinchenko , I

    Pazukhin , A. G., Zinchenko , I. I., Trofimova , E. A., & Henkel , C. 2022, Astronomy Reports, 66, 1302

  67. [67]

    G., Zinchenko , I

    Pazukhin , A. G., Zinchenko , I. I., Trofimova , E. A., Henkel , C., & Semenov , D. A. 2023, , 526, 3673

  68. [68]

    2026, , 997, 296

    Petrashkevich , I., Punanova , A., Caselli , P., et al. 2026, , 997, 296

  69. [69]

    V., Punanova , A

    Petrashkevich , I. V., Punanova , A. F., Caselli , P., et al. 2024, , 528, 1327

  70. [70]

    M., Poynter , R

    Pickett , H. M., Poynter , R. L., Cohen , E. A., et al. 1998, , 60, 883

  71. [71]

    G., & Thompson , M

    Pillai , T., Wyrowski , F., Hatchell , J., Gibb , A. G., & Thompson , M. A. 2007, , 467, 207

  72. [72]

    E., Sipilä, O., Segura-Cox, D

    Pineda, J. E., Sipilä, O., Segura-Cox, D. M., et al. 2024, , 686, A162

  73. [73]

    2016, , 587, A118

    Punanova , A., Caselli , P., Pon , A., Belloche , A., & Andr \'e , P. 2016, , 587, A118

  74. [74]

    2014, , 440, 3557

    Reboussin , L., Wakelam , V., Guilloteau , S., & Hersant , F. 2014, , 440, 3557

  75. [75]

    2019, , 629, A15

    Redaelli , E., Bizzocchi , L., Caselli , P., et al. 2019, , 629, A15

  76. [76]

    2024, , 685, A67

    Redaelli , E., Bovino , S., Lupi , A., et al. 2024, , 685, A67

  77. [77]

    2025, , 702, A210

    Redaelli , E., Bovino , S., Sabatini , G., et al. 2025, , 702, A210

  78. [78]

    & Gerin , M

    Roueff , E. & Gerin , M. 2003, , 106, 61

  79. [79]

    C., van der Tak , F

    Roueff , E., Lis , D. C., van der Tak , F. F. S., Gerin , M., & Goldsmith , P. F. 2005, , 438, 585

  80. [80]

    2016, , 459, 3756

    Ruaud , M., Wakelam , V., & Hersant , F. 2016, , 459, 3756

Showing first 80 references.