REVIEW 2 major objections 5 minor 72 references
The ALMA-ATOMS survey: Vibrationally excited HC$_3$N lines in hot cores
T0 review · 2 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Vibrationally excited cyanoacetylene in hot cores is lit up by mid-infrared dust radiation, not by collisions with hydrogen molecules.
desk verdict Useful HC3N* survey catalog, but the 'collisions ineffective' claim rests on an 8x underestimate of n(H2). read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument pivots on the critical density ncrit for collisional excitation of the lowest HC3N vibrational levels, taken as 4e8 $cm^{-3}$ for v7=1 at 300 K, with higher vibrational states requiring much larger densities. The survey derives H2 number densities from 3 mm dust continuum using a dust temperature adopted from complex organic molecule tracers, a gas-to-dust ratio, and a fixed dust opacity, then classifies each core by how many vibrational states are detected. Because nearly every fitted core has n(H2) below ncrit, the authors conclude that collisions cannot sustain the observed vibrational population and that mid-infrared pumping must be the dominant excitation route.
What would settle it
Observe a hot core at sub-arcsecond resolution, measure the local H2 density at the HC3N* line-forming region, and characterize the local mid-infrared field: a bright v7=1 line arising from gas shown to be denser than 4e8 $cm^{-3}$ while embedded in a weak mid-infrared environment would rule out infrared pumping as the dominant excitation mechanism.
Extended reading notes
Core claim
The central claim is that in hot cores the vibrational population of HC3N is pumped by radiation: molecules absorb mid-infrared photons from warm dust and climb from the ground vibrational state into states such as v7=1, v7=2, v6=1, v5=1, and v4=1. The observed support is that 52 of 60 hot cores show at least one HC3N* line, that hotter cores show more vibrational states, and that the derived H2 densities are below the critical densities for collisional excitation of the low-lying vibrational modes in 28 of the 29 cores that could be fitted. The paper also reports two column-density thresholds, about 1.9e16 and 7.6e16 $cm^{-2}$, that separate cores with few vibrational states from those with many, and it finds no significant difference in rotation temperature between cores with and without ultra-compact Hii regions.
Load-bearing premise
The comparison assumes that the H2 number density averaged over a 1.2-1.9 arcsec beam (about 2000-20000 au) is representative of the gas that actually emits the HC3N* lines, so that comparing that average with the 300 K critical density is a meaningful test of whether collisions matter.
Editorial extensions
If this is right
- HC3N* lines can serve as a diagnostic of the mid-infrared radiation field in hot cores rather than simply a probe of gas density.
- Counting how many HC3N* states appear in a core offers a quick observational classification of the core's heating environment.
- The two column-density thresholds, roughly 1.9e16 and 7.6e16 cm^-2, mark the stages at which higher-energy vibrational states begin to appear.
- In the single core whose H2 density exceeds the v7=1 critical density, collisional excitation may become relevant, making it a natural test case for the boundary between pumping regimes.
Reading between the lines
- If mid-IR pumping is the dominant route, HC3N* column densities could in principle be turned into estimates of the local dust radiation field, but doing so would require radiative-transfer modeling that accounts for line optical depth and beam dilution.
- Higher-resolution observations that resolve the line-forming region might uncover clumps denser than the beam-averaged value, which would raise the estimated collisional contribution.
- The same pumping logic should apply to other small linear molecules with low-lying bending vibrations, so vibrationally excited lines of related cyanopolyynes may also trace infrared fields rather than density alone.
- A clean test of the mechanism would be to compare HC3N* line intensities with the mid-infrared luminosity of the embedded protostars: a tight correlation would confirm radiative pumping, while the absence of one would point back toward collisions.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a survey of vibrationally excited HC3N (J=11-10) lines in 60 hot cores using ALMA-ATOMS Band 3 data. Detections are reported in 52 cores; for 29 cores with at least three available lines, XCLASS LTE fits yield rotation temperatures between 160 and 335 K and HC3N* column densities of about 7e15 to 1.7e18 cm^-2. The authors derive H2 masses, column densities, and number densities from the 3 mm dust continuum, find correlations between the number of detected HC3N* states and the derived column density and gas temperature, and conclude that HC3N* excitation is mainly driven by mid-IR pumping from hot dust, with collisional excitation ineffective because most cores have n(H2) below the v7=1 critical density of 4e8 cm^-3.
Significance. If the central conclusion were robust, this paper would establish vibrationally excited HC3N as a tracer of the mid-IR radiation field in hot cores rather than a purely density-sensitive tracer, with consequences for interpreting HC3N* line surveys. The observational catalog—52 detections, 29 fitted cores, line frequencies, XCLASS parameters, and continuum-derived core properties—is a useful community resource, and the authors provide an independent cross-check of their rotation temperatures via rotational temperature diagrams. The main weakness is that the load-bearing comparison between n(H2) and n_crit is affected by a systematic error in the radius definition and by beam averaging, so the unqualified conclusion is currently stronger than the data support.
major comments (2)
- [Sec. 3.2.1–3.2.2, Eq. (6), Tables C.1–C.4, Fig. 9] The definition of R_core in Sec. 3.2.1 as sqrt(theta_maj * theta_min) is a geometric mean of FWHM diameters, but Eq. (6) uses it directly as the radius of a uniform sphere. This underestimates the sphere volume by a factor of 8 (since (FWHM/2)^3 vs. FWHM^3) and therefore underestimates n(H2) by a factor of 8. With this correction, several entries in Tables C.1–C.4 approach or exceed the v7=1 critical density of 4e8 cm^-3: I18507+0110 is already above it (5.0e8 cm^-3), and I17016-4124c1 (2.5e8 -> 2.0e9), I19095+0930 (2.0e8 -> 1.6e9), and I15254-5621 (1.5e8 -> 1.2e9) would all cross it. The statement that 28/29 fitted cores are below n_crit is therefore not supported by the written formula. The authors should recompute n(H2) using a consistent radius (e.g., R_core/2) and propagate the effect into all density-dependent inferences and into the Abstract and Sect. 5 conclusions.
- [Sec. 3.3, Sec. 4.3, Fig. 3] Even with a corrected radius, Eq. (6) gives a source-averaged uniform-sphere density, whereas HC3N* emission arises from the central, densest part of the core. The paper itself notes in Sec. 3.3 that the 3 mm data are unable to resolve the emission distribution. For a plausible power-law density profile with p = 1.5–2, the density at the inner radius where v7=1 is excited is one to two orders of magnitude above the beam-averaged value, so the simple comparison of the averaged n(H2) with n_crit is not a secure test of whether collisions are effective in the emitting gas. I request a quantitative test: either use the observed continuum radial profiles to estimate central densities, or explicitly state that the conclusion applies only to beam-averaged densities and acknowledge that collisional excitation may be significant in the inner regions.
minor comments (5)
- [Sec. 4.3] The sentence 'All hot cores in our sample except for I18507+0110, for which n(H2) is lower than ncrit' is internally inconsistent with Table C.3, which lists n(H2) = 5.0e8 cm^-3 for I18507+0110, above 4e8 cm^-3; the text should read 'higher'. After the radius correction in the major comment, additional cores will also have n(H2) > n_crit, so the wording should be revised to reflect the actual exceptions.
- [Table C.1] The dust temperature for I18316-0602 is listed as 108 ± 204 K; the uncertainty is larger than the value and the lower bound is unphysical, suggesting a typographical error that should be corrected.
- [Abstract and Sect. 5, item (5)] The phrase 'collisional excitation is ineffective' is stated without qualification even though I18507+0110 already has n(H2) above the v7=1 critical density in the current tables; consider qualifying the conclusion to 'collisional excitation is ineffective for most cores' or to specific density regimes.
- [Sec. 3.1] The abbreviation HC3N* is used in the Abstract before it is formally defined in the Introduction; please define it at first use.
- [Fig. 7] The linear fits in the two panels are labeled only by the fitted equation; adding the correlation coefficient or reduced chi-squared would let the reader judge the claimed tightness of the relation for different excitation types.
Circularity Check
No significant circularity: the central mid-IR pumping conclusion is anchored to an external critical-density benchmark and to independently fitted dust-continuum densities, not to a self-referential fit.
full rationale
The paper's central inference (Abstract; Sect. 5 item 5) is that HC3N* excitation is mainly driven by mid-IR pumping because the measured H2 number densities lie below the v7=1 critical density of 4e8 cm^-3 (Sect. 4.3, Fig. 9). The two inputs to this comparison are independent of each other and of the conclusion: n(H2) is obtained from 3 mm dust continuum via Eqs. (3)-(6) (Hildebrand 1983; Ossenkopf & Henning 1994), while n_crit is taken from external quantum-chemical and experimental work (Wyrowski et al. 1999; Leach et al. 2014; CDMS/JPL line data). The LTE rotation temperatures and HC3N* column densities are fitted with XCLASS using externally catalogued line parameters and are cross-checked against an independent rotational temperature diagram (Sect. 3.2.2, Appendix E). Citations to Liu et al. (2020) and Qin et al. (2022) supply the ATOMS sample, distances, and adopted dust temperatures; these are project data papers with overlapping authors, but the key comparison does not reduce to them. The correlation between the number of detected states and fitted T_rot or column density is partly a selection/fit coupling, but it is not used as the derivation of the main conclusion. The skeptically noted volume convention for R_core (using the geometric-mean FWHM as a sphere radius in Eq. 6) is a physical or calibration concern that could bias n(H2), but it is not a case of a fitted input being renamed as a prediction, so it does not constitute circularity.
Assumptions & free parameters
free parameters (5)
- Rotation temperature T_rot per core =
160-335 K across 29 cores (Table B.1)
- HC3N* column density N per core =
6.9e15 to 1.7e18 cm^-2 (Table B.1)
- Dust temperature T_d per core =
100-230 K (Tables C.1-C.4)
- Distinguishing HC3N* column-density thresholds =
1.9e16 and 7.6e16 cm^-2
- Peak-separation threshold =
4500 au
assumptions (5)
- domain assumption LTE holds for the vibrationally excited HC3N lines, so a single rotation temperature and column density describe all observed states.
- domain assumption The 3 mm dust continuum is optically thin and can be converted to gas mass and H2 column density with a single dust opacity (0.24 cm^2 g^-1), a gas-to-dust ratio from Giannetti et al. (2017), and mean molecular weight 2.8.
- domain assumption Free-free emission from UC H ii regions is optically thin at 99 GHz with T_e=6000 K and N(He+)/N(H+)=0.08, so S_ff can be estimated from the H40alpha integrated intensity and subtracted.
- domain assumption If n(H2) is below n_crit for a given vibrational level, collisional excitation is ineffective, so the observed population must come from mid-IR pumping; the relevant critical densities are those for 300 K from Wyrowski et al. (1999).
- domain assumption The hot cores in the sample, selected via COM emission in Qin et al. (2022), are representative of hot cores and the continuum emission traces the same gas that emits HC3N*.
Cite this review
Pith. "Pith review of The ALMA-ATOMS survey: Vibrationally excited HC$_3$N lines in hot cores." pith.science (2026). https://pith.science/paper/TIYE76DL
@misc{pith2026241212546,
author = {Pith},
title = {Pith review of: The ALMA-ATOMS survey: Vibrationally excited HC$_3$N lines in hot cores},
year = {2026},
howpublished = {\url{https://pith.science/paper/TIYE76DL}},
note = {Machine review of arXiv:2412.12546}
}
abstract
Interstellar molecules are excellent tools for studying the physical and chemical environments of massive star-forming regions. In particular, vibrationally excited HC$_3$N (HC$_3$N*) lines are the key tracers for probing hot cores environments. We present the Atacama Large Millimeter/submillimeter Array (ALMA) 3 mm observations of HC$_3$N* lines in 60 hot cores, aiming to investigate how physical conditions affect the excitation of HC$_3$N* transitions. We have used the XCLASS for line identification. Under the assumption of local thermodynamic equilibrium (LTE), we derived the rotation temperature and column density of HC$_3$N* transitions in hot cores. Additionally, we calculated the H$_2$ column density and number density, along with the abundance of HC$_3$N* relative to H$_2$, to enable a comparison of the physical properties of hot cores with different numbers of HC$_3$N* states. We have detected HC$_3$N* lines in 52 hot cores, in which 29 cores showing more than one vibrationally excited state. Hot cores with higher gas temperatures have more detections of these vibrationally excited lines. The excitation of HC$_3$N* requires dense environments, with its spatial distribution influenced by the presence of UC Hii regions. The observed column density of HC$_3$N* contributes to the number of HC$_3$N* states in hot core environments. After analyzing the various factors influencing HC$_3$N* excitation in hot cores, we conclude that the excitation of HC$_3$N* is mainly driven by mid-IR pumping, while collisional excitation is ineffective.
Figures
Figures from the paper (6 more)
Reference graph
Works this paper leans on
- [1]
-
[2]
Armentrout , W. P. 2018, PhD thesis, West Virginia University
work page 2018
-
[3]
Belloche , A., M \"u ller , H. S. P., Menten , K. M., Schilke , P., & Comito , C. 2013, , 559, A47
work page 2013
- [4]
- [5]
-
[6]
Bonfand , M., Csengeri , T., Bontemps , S., et al. 2024, , 687, A163
work page 2024
-
[7]
2022, , 134, 114501
CASA Team , Bean , B., Bhatnagar , S., et al. 2022, , 134, 114501
2022
-
[8]
Chung , H. S., Osamu , K., & Masaki , M. 1991, Journal of Korean Astronomical Society, 24, 217
work page 1991
Show all 72 references
-
[9]
2002, , 40, 27
Churchwell , E. 2002, , 40, 27
2002
-
[10]
Condon , J. J. & Ransom , S. M. 2016, Essential Radio Astronomy
2016
-
[11]
& Aalto , S
Costagliola , F. & Aalto , S. 2010, , 515, A71
2010
-
[12]
2000, , 361, 1058
de Vicente , P., Mart \' n-Pintado , J., Neri , R., & Colom , P. 2000, , 361, 1058
2000
-
[13]
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
2016
-
[14]
B., Cernicharo , J., Viti , S., et al
Esplugues , G. B., Cernicharo , J., Viti , S., et al. 2013, , 559, A51
2013
-
[15]
M., Beltr \'a n , M
Frau , P., Girart , J. M., Beltr \'a n , M. T., et al. 2010, , 723, 1665
2010
-
[16]
A., & Goicoechea , J
Gerin , M., Neufeld , D. A., & Goicoechea , J. R. 2016, , 54, 181
2016
-
[17]
2017, , 606, L12
Giannetti , A., Leurini , S., K \"o nig , C., et al. 2017, , 606, L12
2017
-
[18]
F., Krotkov , R., & Snell , R
Goldsmith , P. F., Krotkov , R., & Snell , R. L. 1985, , 299, 405
1985
-
[19]
Goldsmith , P. F. & Langer , W. D. 1999, , 517, 209
1999
-
[20]
F., Snell , R
Goldsmith , P. F., Snell , R. L., Deguchi , S., Krotkov , R., & Linke , R. A. 1982, , 260, 147
1982
-
[21]
2021, , 253, 2
He , Y.-X., Henkel , C., Zhou , J.-J., et al. 2021, , 253, 2
2021
-
[22]
& van Dishoeck , E
Herbst , E. & van Dishoeck , E. F. 2009, , 47, 427
2009
-
[23]
Hildebrand , R. H. 1983, , 24, 267
1983
-
[24]
K., Belloche , A., & Garrod , R
J rgensen , J. K., Belloche , A., & Garrod , R. T. 2020, , 58, 727
2020
-
[25]
L., Evans , N
Kauffmann , J., Bertoldi , F., Bourke , T. L., Evans , N. J., I., & Lee , C. W. 2008, , 487, 993
2008
-
[26]
D., Lal , D
Khan , S., Pandian , J. D., Lal , D. V., et al. 2022, , 664, A140
2022
-
[27]
G., Aikin , A
Kunde , V. G., Aikin , A. C., Hanel , R. A., et al. 1981, , 292, 686
1981
-
[28]
A., Mahjoub , A., et al
Leach , S., Garcia , G. A., Mahjoub , A., et al. 2014, , 140, 174305
2014
-
[29]
& Goldsmith , P
Li , D. & Goldsmith , P. F. 2012, , 756, 12
2012
-
[30]
2017, , 466, 248
Li , S., Wang , J., Zhang , Z.-Y., et al. 2017, , 466, 248
2017
-
[31]
2021, , 505, 2801
Liu , H.-L., Liu , T., Evans , Neal J., I., et al. 2021, , 505, 2801
2021
-
[32]
J., Kim , K.-T., et al
Liu , T., Evans , N. J., Kim , K.-T., et al. 2020, , 496, 2790
2020
-
[33]
2024, , 271, 3
Liu , X., Liu , T., Shen , Z., et al. 2024, , 271, 3
2024
-
[34]
2020, , 894, L14
Lu , X., Cheng , Y., Ginsburg , A., et al. 2020, , 894, L14
2020
-
[35]
2024, Research in Astronomy and Astrophysics, 24, 065003
Luo , A.-X., Liu , H.-L., Li , G.-X., Pan , S., & Yang , D.-T. 2024, Research in Astronomy and Astrophysics, 24, 065003
2024
-
[36]
Mallinson , P. D. & Fayt , A. 1976, Molecular Physics, 32, 473
1976
-
[37]
G., Harada , N., et al
Mart \' n , S., Mangum , J. G., Harada , N., et al. 2021, , 656, A46
2021
-
[38]
Mauersberger , R., Henkel , C., & Sage , L. J. 1990, , 236, 63
1990
-
[39]
P., Waters , B., Schiebel , D., Young , W., & Golap , K
McMullin , J. P., Waters , B., Schiebel , D., Young , W., & Golap , K. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw , F. Hill , & D. J. Bell , 127
2007
-
[40]
2017, , 598, A7
M \"o ller , T., Endres , C., & Schilke , P. 2017, , 598, A7
2017
-
[41]
2018, , 56, 41
Motte , F., Bontemps , S., & Louvet , F. 2018, , 56, 41
2018
-
[42]
u ller , H. S. P., Schl \
M \"u ller , H. S. P., Schl \"o der , F., Stutzki , J., & Winnewisser , G. 2005, Journal of Molecular Structure, 742, 215
2005
-
[43]
M \"u ller , H. S. P., Thorwirth , S., Roth , D. A., & Winnewisser , G. 2001, , 370, L49
2001
-
[44]
& Henning , T
Ossenkopf , V. & Henning , T. 1994, , 291, 943
1994
-
[45]
F., et al
Pagani , L., Favre , C., Goldsmith , P. F., et al. 2017, , 604, A32
2017
-
[46]
2017, , 837, 49
Peng , Y., Qin , S.-L., Schilke , P., et al. 2017, , 837, 49
2017
-
[47]
M., Poynter , R
Pickett , H. M., Poynter , R. L., Cohen , E. A., et al. 1998, , 60, 883
1998
-
[48]
2022, , 511, 3463
Qin , S.-L., Liu , T., Liu , X., et al. 2022, , 511, 3463
2022
-
[49]
2018, , 98, 052708
Rankovi \'c , M., Nag , P., Zawadzki , M., et al. 2018, , 98, 052708
2018
-
[50]
Rico-Villas , F., Mart \' n-Pintado , J., Gonz \'a lez-Alfonso , E., Mart \' n , S., & Rivilla , V. M. 2020, , 491, 4573
2020
-
[51]
2021, , 502, 3021
Rico-Villas , F., Mart \' n-Pintado , J., Gonz \'a lez-Alfonso , E., et al. 2021, , 502, 3021
2021
-
[52]
Rosen , A. L. 2022, , 941, 202
2022
-
[53]
2022, , 516, 1983
Saha , A., Tej , A., Liu , H.-L., et al. 2022, , 516, 1983
2022
-
[54]
Schilke , P., Guesten , R., Schulz , A., Serabyn , E., & Walmsley , C. M. 1992, , 261, L5
1992
-
[55]
Shaver , P. A. 1970, , 5, 167
1970
-
[56]
2018, , 856, 141
Tang , M., Liu , T., Qin , S.-L., et al. 2018, , 856, 141
2018
-
[57]
K., & Minamidani , T
Taniguchi , K., Saito , M., Sridharan , T. K., & Minamidani , T. 2018, , 854, 133
2018
-
[58]
A., et al
Taniguchi , K., Sanhueza , P., Olguin , F. A., et al. 2023, , 950, 57
2023
-
[59]
Taniguchi , K., Tanaka , K. E. I., Zhang , Y., et al. 2022, , 931, 99
2022
-
[60]
Thorwirth , S., M \"u ller , H. S. P., & Winnewisser , G. 2000, Journal of Molecular Spectroscopy, 204, 133
2000
-
[61]
Turner , B. E. 1971, , 163, L35
1971
-
[62]
van Dishoeck , E. F. 2018, in IAU Symposium, Vol. 332, Astrochemistry VII: Through the Cosmos from Galaxies to Planets, ed. M. Cunningham , T. Millar , & Y. Aikawa , 3--22
2018
-
[63]
2015, , 575, A84
Velilla Prieto , L., S \'a nchez Contreras , C., Cernicharo , J., et al. 2015, , 575, A84
2015
-
[64]
F., et al
Wang , J., Li , D., Goldsmith , P. F., et al. 2020, , 889, 129
2020
-
[65]
Wyrowski , F., Schilke , P., & Walmsley , C. M. 1999, , 341, 882
1999
-
[66]
Yamada , K. M. T. & Creswell , R. A. 1986, Journal of Molecular Spectroscopy, 116, 384
1986
-
[67]
2019, , 489, 4497
Yu , N., Wang , J.-J., & Xu , J.-L. 2019, , 489, 4497
2019
-
[68]
2021, Research in Astronomy and Astrophysics, 21, 014
Yue , Y.-H., Qin , S.-L., Liu , T., et al. 2021, Research in Astronomy and Astrophysics, 21, 014
2021
-
[69]
& Li , G.-X
Zhang , C.-P. & Li , G.-X. 2017, , 469, 2286
2017
-
[70]
I., Dewangan , L
Zinchenko , I. I., Dewangan , L. K., Baug , T., Ojha , D. K., & Bhadari , N. K. 2021, , 506, L45
2021
-
[71]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sent...
-
[72]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.