arxiv: 2605.08790 · v1 · submitted 2026-05-09 · 🌌 astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

The evolution of C4H and c-C3H2 in molecular cores

Yijia Liu , Junzhi Wang , Ningyu Tang , Yajiang Lu , Donghui Quan , Juan Li , Kai Yang , Shu Liu

show 3 more authors

Yuqiang Li Siqi Zheng Chao Ou

Authors on Pith no claims yet

Pith reviewed 2026-05-12 00:59 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords C4Hc-C3H2molecular coreschemical evolutionabundance ratiosstar-forming regionsinterstellar hydrocarbonsH13CO+ tracer

0 comments

The pith

Abundance ratios of C4H and c-C3H2 relative to H13CO+ decrease as molecular cores evolve.

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

The paper maps C4H 9-8 and c-C3H2 2-1 lines across 31 regions in 22 massive star-forming sources and combines the results with prior cold-core measurements. It reports that both C4H/H13CO+ and c-C3H2/H13CO+ ratios fall steadily from cold cores through more evolved stages. The decline is interpreted as evidence that these small unsaturated hydrocarbons are consumed and fed into pathways that build larger organic molecules. Emission is found concentrated at the edges of HII regions traced by H42, consistent with precursor chemistry occurring at the interfaces between molecular gas and ionized zones. A chemical model that tracks gas, grain-surface, and ice-mantle reactions is used to compare against the observed abundances.

Core claim

By mapping C4H 9-8, c-C3H2 2-1, H13CO+ 1-0, and H42 in massive star-forming regions and merging the data with cold-core observations, the C4H/H13CO+ and c-C3H2/H13CO+ ratios show strong decreasing trends as molecular cores evolve. This indicates that small unsaturated hydrocarbons are consumed and converted into other organic molecules during core evolution. The emission of C4H and c-C3H2 is concentrated at the edges of HII regions, supporting their role as precursors in interstellar chemical pathways.

What carries the argument

The evolutionary trend in the relative abundance ratios C4H/H13CO+ and c-C3H2/H13CO+, measured by on-the-fly mapping and interpreted with a gas-dust-ice chemical model.

If this is right

Small unsaturated hydrocarbons such as C4H and c-C3H2 are converted into more complex organic species as cores evolve.
The spatial concentration at HII-region edges marks the sites where precursor chemistry feeds complex-molecule formation.
A multi-phase chemical model is required to reproduce the observed drop in relative abundances.
The ratios can serve as evolutionary diagnostics for molecular cores in other sources.

Where Pith is reading between the lines

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

If the ratio decline is monotonic, it could be calibrated as a rough chemical clock for core age.
Repeating the measurements in low-mass star-forming regions would test whether the same consumption pathway operates across different mass regimes.
Incorporating these trends into larger chemical networks may help predict the timing of specific complex organics such as methanol or larger chains.

Load-bearing premise

The combined sample of massive star-forming regions and cold cores traces a single evolutionary sequence without major selection biases or unaccounted environmental differences that could alter the abundance ratios.

What would settle it

Finding that C4H/H13CO+ and c-C3H2/H13CO+ ratios do not decrease (or increase) when additional cold cores are compared to evolved regions in an unbiased sample would falsify the claimed consumption trend.

Figures

Figures reproduced from arXiv: 2605.08790 by Chao Ou, Donghui Quan, Juan Li, Junzhi Wang, Kai Yang, Ningyu Tang, Shu Liu, Siqi Zheng, Yajiang Lu, Yijia Liu, Yuqiang Li.

**Figure 1.** Figure 1: Velocity-integrated intensity maps and spatial averaged spectra of C4H 9–8, c-C3H2 2–1, H13CO+ 1–0, and H42α. The source names are presented in the maps and spectra. The grey scale colour at the right is in units of K km s−1 . (a) and (b) Velocity-integrated intensity maps of G015.03−00.67, where panel (a) shows C4H 9–8 (red contours) overlaid on c-C3H2 2–1 (blue contours and grey scale) and panel (b) show… view at source ↗

**Figure 2.** Figure 2: Relation between c-C3H2/H 13CO+ and C4H/H 13CO+ abundance ratios. Data are from 31 regions in 22 hot cores (red points) and 22 regions in 19 cold cores (blue points). Additionally, it has been established that C4H forms efficiently via gas-phase ion–molecule reactions, a process particArticle number, page 5 of 16 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Figure 3: Temporal evolution of c-C3H2, C4H, and HCO+ abundances predicted by models in cold core and hot core (warm-up). According to the model results, in cold cores, c-C3H2 is mainly formed via the dissociative recombination of C3H + 3 with e − , while C4H is primarily produced through the dissociative recombination of C4H + 2 with e − . During the collapse stage of the hot-core model, the formation pathways of … view at source ↗

**Figure 4.** Figure 4: Net percentage contributions of the formation and destruction pathways of c-C3H2 in cold and hot cores. (a) Net percentage contributions of the main formation pathways of c-C3H2 in cold molecular cloud cores. (b) Net percentage contributions of the main formation pathways of c-C3H2 during the hot-core collapse stage. (c) Net percentage contributions of the main formation pathways of c-C3H2 during the hot-c… view at source ↗

**Figure 5.** Figure 5: Net percentage contributions of the formation and destruction pathways of C4H in cold and hot cores. (a) Net percentage contributions of the main formation pathways of C4H in cold molecular cloud cores. (b) Net percentage contributions of the main formation pathways of C4H during the hot-core collapse stage. (c) Net percentage contributions of the main formation pathways of C4H during the hot-core warm-up … view at source ↗

read the original abstract

Linear C4H and cyclic c-C3H2, as small unsaturated hydrocarbons, are the key precursors to complex organic molecules and are critical components of the interstellar medium. We present on-the-fly mapping observations of C4H 9-8 lines, c-C3H2 2-1, H13CO+ 1-0, and H42 toward a sample of 22 massive star-forming regions using the IRAM 30m telescope. Our aim is to further explore the evolution of these carbon-chain molecules by combining observational results obtained in cold cores. We employed H13CO+ 1-0 and H42 as tracers to probe the positions of molecular cloud cores and ionised hydrogen regions (HII regions), respectively. One chemical model in particular, which includes gas, dust grain surface, and icy mantle phases for C4H and c-C3H2 molecules, was used to make comparisons with observed abundances. From mapping observations targeting 31 regions across 22 sources, C4H 9-8 (J = 19/2-17/2) and C4H 9-8 (J = 17/2-15/2) were detected in only 17 regions, while H13CO+ 1-0 and c-C3H2 2-1 were successfully detected in all 31 regions. We find that the emission of C4H 9-8 and c-C3H2 2-1 is concentrated at the edges of H42 emission regions. The C4H/H13CO+ and c-C3H2/H13CO+ relative abundance ratios range from 0.17 to 1.77 and 1.42 to 6.69, respectively, with a median C4H/c-C3H2 ratio of 0.13. By combining the observational results of cold cores, we find that C4H/H13CO+ and c-C3H2/H13CO+ ratios show a strong decreasing trend as molecular cores evolve. The decreasing trends in C4H/H13CO+ and c-C3H2/H13CO+ ratios imply that small unsaturated hydrocarbons can be consumed and converted into other organic molecules during the evolution of molecular cores. The spatial concentration of C4H and c-C3H2 emission at the edges of H42 regions further supports their role as precursors in the chemical pathways that lead to complex organic molecules in the interstellar medium.

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.

Desk Editor's Note private letter to a colleague

New IRAM maps give fresh C4H and c-C3H2 data in massive cores plus a combined declining-ratio trend, but the single evolutionary sequence assumption is the main loose thread.

read the letter

The main takeaway is that this work adds new on-the-fly IRAM 30m maps of the C4H 9-8 lines and c-C3H2 2-1 in 22 massive star-forming regions, and when these are put together with prior cold-core observations, the C4H/H13CO+ and c-C3H2/H13CO+ ratios show a decreasing trend with evolutionary stage. The observations are the strongest part. They targeted 31 regions, picked up C4H in 17 of them and c-C3H2 in every one, and used H13CO+ 1-0 and H42 to mark the cores and HII regions. The emission for both molecules sits at the edges of the ionized areas, which fits with their precursor role. The reported ratios span 0.17 to 1.77 for C4H over H13CO+ and 1.42 to 6.69 for c-C3H2, with a median C4H to c-C3H2 of 0.13. They bring in one chemical model covering gas, dust surfaces, and icy mantles to compare against the measured abundances. The combined dataset with cold cores then supports the idea that these small unsaturated hydrocarbons get used up and turned into larger organics as the cores develop. That gives a useful set of numbers and positions for people modeling carbon chemistry in star formation. The spatial concentration at HII edges is a nice observational result that adds to the picture. The part that feels less secure is the assumption that the massive regions and cold cores form a single, clean evolutionary sequence. The paper does not detail the specific metric used to order the massive cores or test whether differences in density, UV exposure, or initial conditions between the two samples could produce the same ratio drop without any chemical evolution. H13CO+ is a reasonable tracer, but if its abundance or excitation changes systematically, it could affect the normalized ratios. The model comparison is there, but without more on how the parameters were chosen or if separate runs were done for each core type, it is hard to judge how well it backs the trend. This is targeted at astrochemists who study small hydrocarbons in molecular clouds and star-forming regions. A reader looking for updated abundance data and spatial maps in massive cores will find value here, even if the broader evolutionary claim needs more support. I would recommend sending it to peer review. The new observations are worth the time, provided the authors address the staging and bias questions in revision.

Referee Report

3 major / 2 minor

Summary. The manuscript reports on-the-fly mapping observations with the IRAM 30m telescope of the C4H 9-8 lines, c-C3H2 2-1, H13CO+ 1-0 and H42 toward 22 massive star-forming regions (31 positions total). C4H is detected in 17 positions, the other lines in all positions; emission of C4H and c-C3H2 is concentrated at the edges of H42 regions. Column-density ratios C4H/H13CO+ (0.17–1.77) and c-C3H2/H13CO+ (1.42–6.69) are derived, together with a median C4H/c-C3H2 ratio of 0.13. Combining these data with an external cold-core sample yields a decreasing trend in both ratios with evolutionary stage; a single gas-grain chemical model is used for qualitative comparison. The authors interpret the trend as net consumption of the two hydrocarbons into more complex organics.

Significance. If the combined sample can be shown to lie on a single, bias-free evolutionary sequence and the trend is statistically robust, the work would supply direct observational support for the chemical processing of small unsaturated hydrocarbons during core evolution and their role as precursors to complex organic molecules. The mapped spatial anti-correlation with HII-region tracers adds useful morphological context.

major comments (3)

[Abstract and §3] Abstract and §3 (Results): the central claim of a 'strong decreasing trend' with evolutionary stage is based on merging the 22 massive sources with an external cold-core sample, yet no evolutionary staging metric (e.g., infrared luminosity, HII-region size, or chemical age indicator) is defined for the massive cores, nor is any statistical test (Spearman rank, linear regression with uncertainties) reported for the trend.
[Abstract and §4] Abstract and §4 (Chemical modeling): the model is invoked only for 'comparison' but the text supplies neither the initial abundances, density/temperature profiles, nor the quantitative metric (e.g., reduced χ² or abundance ratio residuals) used to judge agreement with the observed C4H/H13CO+ and c-C3H2/H13CO+ values across the two populations.
[Abstract] Abstract: the reported ratio ranges and median are given without uncertainties or discussion of how H13CO+ column density was derived (optical depth, excitation temperature assumptions) and whether these assumptions remain valid when the massive-star-forming and cold-core samples are combined.

minor comments (2)

[§2 and §3] The detection statistics (17/31 for C4H, 31/31 for the others) are stated clearly but would benefit from a table listing individual source positions, integrated intensities, and rms noise values.
[§3] Notation for the two C4H hyperfine components (J = 19/2–17/2 and 17/2–15/2) is correct but should be repeated once in the results text for readers unfamiliar with the 9–8 transition.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed comments. These have highlighted areas where the manuscript can be clarified and strengthened. We address each major comment below and outline the revisions we will make.

read point-by-point responses

Referee: [Abstract and §3] Abstract and §3 (Results): the central claim of a 'strong decreasing trend' with evolutionary stage is based on merging the 22 massive sources with an external cold-core sample, yet no evolutionary staging metric (e.g., infrared luminosity, HII-region size, or chemical age indicator) is defined for the massive cores, nor is any statistical test (Spearman rank, linear regression with uncertainties) reported for the trend.

Authors: We agree that the presentation would benefit from an explicit evolutionary staging metric for the massive cores and from a statistical quantification of the trend. In the revised manuscript we will define stages using the presence and spatial extent of HII regions (traced by H42) as the primary indicator, with the cold-core sample representing the earliest phase. We will also add a Spearman rank correlation analysis (including uncertainties) on the combined dataset to support the reported decreasing trend. revision: yes
Referee: [Abstract and §4] Abstract and §4 (Chemical modeling): the model is invoked only for 'comparison' but the text supplies neither the initial abundances, density/temperature profiles, nor the quantitative metric (e.g., reduced χ² or abundance ratio residuals) used to judge agreement with the observed C4H/H13CO+ and c-C3H2/H13CO+ values across the two populations.

Authors: The single gas-grain model is described qualitatively in §4, but we acknowledge that the manuscript lacks sufficient detail on the model setup and a quantitative measure of agreement. In the revision we will add the adopted initial abundances, the density and temperature profiles, and a simple quantitative comparison (e.g., residuals between observed and modeled ratios) for both the cold-core and massive-star-forming populations. revision: yes
Referee: [Abstract] Abstract: the reported ratio ranges and median are given without uncertainties or discussion of how H13CO+ column density was derived (optical depth, excitation temperature assumptions) and whether these assumptions remain valid when the massive-star-forming and cold-core samples are combined.

Authors: We will revise the abstract and the relevant results section to report uncertainties on the C4H/H13CO+ and c-C3H2/H13CO+ ratios and on the median C4H/c-C3H2 value. We will also include a brief discussion of the H13CO+ column-density derivation (LTE assumptions, optical-depth corrections, and excitation temperature) and explicitly address the applicability of these assumptions when combining the two samples. revision: yes

Circularity Check

0 steps flagged

No circularity; reported trends are direct column-density ratios from observations

full rationale

The paper's central results are abundance ratios (C4H/H13CO+ and c-C3H2/H13CO+) computed from detected line intensities in 31 regions across 22 sources, combined with external cold-core data. These ratios are presented as measured quantities that decrease with an assumed evolutionary ordering; no equations, fits, or models generate the trends as outputs. The chemical model is used solely for post-hoc comparison. No self-citations, uniqueness theorems, or ansatzes are invoked to justify the reported decline, and the quantities are not redefined in terms of themselves. The interpretive step of labeling the ordering 'evolutionary' is an assumption about sample merging, not a reduction of the data to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claim rests on standard radio-astronomy assumptions for line detection and abundance derivation plus comparison to one existing chemical model; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (1)

domain assumption Standard assumptions for converting integrated line intensities to column densities and abundances in molecular clouds
Implicit in all reported abundance ratios and comparisons to the chemical model.

pith-pipeline@v0.9.0 · 5803 in / 1277 out tokens · 39311 ms · 2026-05-12T00:59:49.217669+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

[1]

& Viti, S

Awad, Z. & Viti, S. 2022, MNRAS, 511, 3,

work page 2022
[2]

2001, ApJ, 552, 2,

doi:10.1093/mnras/stac297 Aikawa, Y ., Ohashi, N., Inutsuka, S.-ichiro, et al. 2001, ApJ, 552, 2,

work page doi:10.1093/mnras/stac297 2001
[3]

D., et al

doi:10.1086/320551 Canosa, A., Páramo, A., Le Picard, S. D., et al. 2007, Icarus, 187, 2,

work page doi:10.1086/320551 2007
[4]

M., Roueff, E., et al

doi:10.1016/j.icarus.2006.10.009 Fontani, F., Rivilla, V . M., Roueff, E., et al. 2025, A&A, 700, A245. doi:10.1051/0004-6361/202555901 Fossé, D., Cernicharo, J., Gerin, M., et al. 2001, ApJ, 552, 1,

work page doi:10.1016/j.icarus.2006.10.009 2006
[5]

A., Grassi, T., & Wester, R

doi:10.1086/320471 Gianturco, F. A., Grassi, T., & Wester, R. 2016, Journal of Physics B Atomic Molecular Physics, 49, 204003. doi:10.1088/0953-4075/49/20/204003 Ginsburg, A., Bressert, E., Bally, J., et al. 2012, ApJ, 758, 2, L29. doi:10.1088/2041-8205/758/2/L29 Harada, N. & Herbst, E. 2008, ApJ, 685, 1,

work page doi:10.1086/320471 2016
[6]

& van Dishoeck, E

doi:10.1086/590468 Herbst, E. & van Dishoeck, E. F. 2009, ARA&A, 47, 1,

work page doi:10.1086/590468 2009
[7]

& Osamura, Y

doi:10.1146/annurev-astro-082708-101654 Herbst, E. & Osamura, Y . 2008, ApJ, 679, 2,

work page doi:10.1146/annurev-astro-082708-101654 2008
[8]

M., Hoglund, B., Friberg, P., et al

doi:10.1086/587803 Irvine, W. M., Hoglund, B., Friberg, P., et al. 1981, ApJ, 248, L113. doi:10.1086/183637 Article number, page 7 of 16 A&A proofs:manuscript no. aa 101 102 103 104 105 106 Time (years) 0 20 40 60 80 100 120Percentage (%) C4H+ 2 + e -> C4H + H C5H+ + e -> C4H + C C + C3H2 -> C4H + H H + C4 -> C4H + e- O + C5H -> C4H + CO 10 12 10 11 10 10...

work page doi:10.1086/587803 1981
[9]

GlobularClusterIntrinsicIronAbundanceSpreads.I.Catalog,

doi:10.3847/1538- 4357/ad47a3 Loison, J.-C., Wakelam, V ., Hickson, K. M., et al. 2014, MNRAS, 437, 1,

work page doi:10.3847/1538- 2014
[10]

2017, MNRAS, 470, 4,

doi:10.1093/mnras/stt1956 Loison, J.-C., Agúndez, M., Wakelam, V ., et al. 2017, MNRAS, 470, 4,

work page doi:10.1093/mnras/stt1956 2017
[11]

& Liszt, H

doi:10.1093/mnras/stx1265 Lucas, R. & Liszt, H. S. 2000, A&A, 358,

work page doi:10.1093/mnras/stx1265 2000
[12]

C., Irvine, W

Madden, S. C., Irvine, W. M., Matthews, H. E., et al. 1989, AJ, 97,

work page 1989
[13]

doi:10.1086/115081 Müller H. S. P., Schlöder F., Stutzki J., Winnewisser G., 2005, JMoSt, 742,

work page doi:10.1086/115081 2005
[14]

u ller , H. S. P., Schl \

doi:10.1016/j.molstruc.2005.01.027 Olano, C. A., Walmsley, C. M., & Wilson, T. L. 1988, A&A, 196,

work page doi:10.1016/j.molstruc.2005.01.027 2005
[15]

J., Menten, K

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783,

work page 2014
[16]

J., Menten, K

doi:10.1088/0004-637X/783/2/130 Remijan, A., Scolati, H. N., Burkhardt, A. M., et al. 2023, ApJ, 944, L45. doi:10.3847/2041-8213/acb648 Rodgers, S. D. & Charnley, S. B. 2003, ApJ, 585, 1,

work page doi:10.1088/0004-637x/783/2/130 2023
[17]

D., & Charnley, S

doi:10.1086/345497 Ruaud, M., Wakelam, V ., & Hersant, F. 2016, MNRAS, 459, 4,

work page doi:10.1086/345497 2016
[18]

2008, ApJ, 672, 1,

doi:10.1093/mnras/stw88 Sakai, N., Sakai, T., Hirota, T., et al. 2008, ApJ, 672, 1,

work page doi:10.1093/mnras/stw88 2008
[19]

& Yamamoto, S

doi:10.1086/523635 Sakai, N. & Yamamoto, S. 2013, Chemical Reviews, 113, 12,

work page doi:10.1086/523635 2013
[20]

C., et al

doi:10.1021/cr4001308 Taniguchi, K., Gorai, P., Tan, J. C., et al. 2024, A&A, 692, A65. doi:10.1051/0004-6361/202451499 Takakuwa, S., Ohashi, N., & Aikawa, Y . 2011, ApJ, 728, 2,

work page doi:10.1021/cr4001308 2024
[21]

M., & Gottlieb, C

doi:10.1088/0004-637X/728/2/101 Thaddeus, P., Vrtilek, J. M., & Gottlieb, C. A. 1985, ApJ, 299, L63. doi:10.1086/184581 Wakelam, V ., Gratier, P., Loison, J.-C., et al. 2024, A&A, 689, A63. doi:10.1051/0004-6361/202450606 Article number, page 8 of 16 Liu et al.: The evolution of C4H andc-C 3H2 in molecular cores Appendix A: The spatial distribution maps a...

work page doi:10.1088/0004-637x/728/2/101 1985