REVIEW 4 major objections 6 minor 64 references
Discovery of ionized circumstellar gas emission around the long-period Cepheid $\ell$ Carinae with ALMA
T0 review · 4 major / 6 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read ALMA observations of the long-period Cepheid ℓ Carinae reveal ionized circumstellar gas: a millimeter continuum 2.5 times the photospheric flux and a hydrogen H29α recombination line at the stellar rest velocity.
desk verdict Likely real detection of ionized gas around ℓ Car, but the H29α line rests on a single stacked dataset and the abstract oversells the profile symmetry. 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 is carried by two complementary tracers of ionized gas seen at millimeter wavelengths. The first is a free-free continuum excess over the Rayleigh-Jeans photospheric prediction, quantified by the spectral index $S_ u \propto \nu^{+1.26}$, which distinguishes an optically thin wind, an optically thick blackbody-like source, and a partially optically thick chromosphere with a steep density gradient. The second is the hydrogen radio recombination line H29$\alpha$ (the $n=29\to 28$ transition at 256.3 GHz), an unambiguous tracer of ionized hydrogen whose line-to-continuum ratio yields an electron temperature estimate. The photospheric baseline itself is set by Eq. 1, a Rayleigh-Jeans blackbody using interferometric uniform-disk angular diameters and literature effective temperatures.
What would settle it
Measure the angular size of the 1 mm source with long-baseline interferometry: if the ~3.5 mJy continuum is confined to the ~2.9 mas photospheric disk rather than extended, the ionized circum-stellar interpretation fails. A complementary check is to detect additional hydrogen recombination lines, such as H30α or H40α, and verify that their centroids and widths match the stellar rest velocity seen in H29α.
Extended reading notes
Core claim
The paper establishes that ℓ Carinae has a millimeter emission component that cannot be photospheric. Combining ALMA continuum measurements at 221.6 and 264.9 GHz with published angular diameters and effective temperatures, the authors predict photospheric fluxes of 1.23 and 1.77 mJy; the observed fluxes of 3.31 and 3.99 mJy exceed these predictions by factors of 2.7 and 2.2. The spectral index of the combined flux density is +1.26±0.44, characteristic of partially optically thick ionized gas rather than an optically thin wind. The detection of the H29α radio recombination line at 256.3 GHz, with a line flux about 35% of the continuum and a width of 55.3±7.5 km/s, identifies the emitting material as ionized hydrogen. From the line-to-continuum ratio the authors estimate an electron temperature on the order of $10^4$ K, and they interpret the emission as a chromosphere of ionized gas similar to what is observed around red supergiants.
Load-bearing premise
The predicted photospheric millimeter flux assumes the continuum forms at the same radius as the optical and infrared photosphere, using literature values of effective temperature and angular diameter; if the true millimeter-emitting radius or temperature is larger, the reported excess shrinks, although the H29α line independently supports the presence of ionized gas.
Editorial extensions
If this is right
- ℓ Carinae becomes the first Cepheid with a direct recombination-line tracer of ionized circumstellar gas, making chromospheric envelopes an observed rather than hypothetical feature of Cepheids.
- At millimeter wavelengths the star is about 2.5 times brighter than its photosphere, so Cepheid spectral energy distributions that omit ionized gas will underpredict radio and millimeter fluxes.
- The spectral index of +1.26 favors a partially optically thick chromosphere with a steep density gradient over an optically thin ionized wind, which means the emission cannot yet be converted into a mass-loss rate.
- The lack of significant variability across the sampled pulsation phases suggests the ionized gas is not a pulsation-driven wind, in contrast to the behavior reported for δ Cephei.
- If other long-period Cepheids resemble ℓ Carinae, extragalactic Cepheids observed by JWST could host similar ionized envelopes, and whether that affects the period-luminosity relation remains an open question.
Reading between the lines
- A survey of Cepheids across pulsation periods could test whether the millimeter excess scales with luminosity or convection activity; shorter-period Cepheids would be predicted to show weaker chromospheric emission if convection is the heating mechanism.
- The electron temperature estimate of roughly 23,500 K rests on an optically thin assumption and a single line; observing two recombination lines simultaneously would check this by requiring a consistent temperature from both.
- Calibration systematics between separate spectral setups contribute about ±0.4 to the spectral index; observing the same source in one ALMA band with multiple widely spaced spectral windows could sharpen the density-gradient diagnosis.
- Continuum subtraction for future Cepheid line studies at millimeter wavelengths will need to account for this excess, since the ionized component is a large fraction of the total flux.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports ALMA Band 6 observations of the long-period Cepheid ℓ Carinae in two spectral setups, centered near 221 and 261 GHz. It claims (1) a millimeter continuum flux of about 3.5 mJy, roughly 2.5 times the predicted photospheric contribution, (2) a spectral index α = +1.26 ± 0.44, and (3) a detection of the H29α recombination line at the stellar rest velocity with an FWHM of 55.3 ± 7.5 km/s. The authors interpret these as evidence of free-free emission from an ionized gas envelope, likely a chromosphere, around the Cepheid, and argue this is the first direct evidence of ionized gas around a Cepheid. The letter includes a journal of observations, an SED assembled from optical to millimeter data, and a comparison to the earlier ionized-gas envelope model of Hocdé et al. (2021).
Significance. If correct, the H29α line detection would be the first direct spectroscopic evidence of ionized gas physically associated with a Cepheid, with implications for Cepheid mass loss, atmospheric structure, and the period-luminosity relation. The paper has several strengths: the photospheric baseline is derived from external interferometric angular diameters and literature effective temperatures rather than from the target result; the comparison to the earlier Hocdé et al. (2021) model is a genuine test; and the authors explicitly provide an alternative spectral-index fit in Appendix C. The main weaknesses are that the decisive H29α line is presented from a single stacked dataset without per-epoch verification, and that the continuum/spectral-index claims are sensitive to a single high-flux epoch whose calibration is argued to be acceptable but not quantitatively demonstrated.
major comments (4)
- [§3.2 and Appendix D] The H29α line detection is the pivotal evidence for the claim of ionized gas, yet it is derived only from a single stacked Setup 2 dataset. Appendix D lists the 2024-10-02 epoch with Sν = 4.653 ± 0.032 mJy and a beam of 0.87″ × 0.77″, while the two other Setup 2 epochs have 3.265 and 3.821 mJy with beams of about 0.2″; the outlier is 22–40% above the others. Because the line and continuum are measured from the same visibilities, this epoch may dominate the stacked line profile, and the upper limit of about 0.2″ on the line extent is not supported by any of the quoted beam sizes (the Setup 2 beam in Table 1 is 0.37″ × 0.30″). I request per-epoch H29α spectra, or at least a re-stacking that excludes the 2024-10-02 epoch, with the resulting line flux, centroid, and FWHM reported. Without this, the abstract's 'direct evidence' wording is premature.
- [§3.1 and Table D.1] The spectral index α = +1.26 ± 0.06 is not robust to the exclusion of the 2024-10-02 continuum epoch. Averaging only the two high-resolution Setup 2 epochs gives a Setup 2 flux of about 3.54 mJy at 264.86 GHz, versus 3.306 mJy at 221.61 GHz for Setup 1; the two-point spectral index then becomes α ≈ 0.4, consistent with optically thin free-free emission rather than the claimed +1.26. The manuscript attributes the 40% excess in the 2024-10-02 epoch to 'more realistic statistical errors' from Francis et al. (2020), but this is an assertion rather than a demonstrated error budget. Please provide α computed with and without that epoch and reconcile the discrepancy.
- [§3.1 and Appendix C] The two spectral-index fitting methods disagree at a level larger than the quoted systematic error. The CASA mt-mfs fit gives α = +1.26 ± 0.06 (Table 1), while the direct log-log fit in Appendix C yields α = +1.78 ± 0.03; the 0.5 difference exceeds the 0.4 systematic uncertainty assigned to cross-calibration. The manuscript prefers the CASA result because of calibration sensitivity, but the reader is not shown why the CASA method is less affected. Since the physical interpretation (density slope, comparison with RSG chromospheres) relies on the value of α, this choice should be justified quantitatively or the quoted uncertainty should encompass both estimates.
- [§3, Eq. (1)] Equation (1) computes the photospheric contribution using the H/K-band uniform-disk angular diameter, assuming the 1 mm continuum forms at the same radius as the optical/IR photosphere. The authors acknowledge this approximation and note that a radio-photosphere can be larger in evolved stars; the maximum-θUD test (θ = 3.2 mas) only covers the pulsation-phase variation of the adopted angular diameter, not a frequency-dependent radius. A 30% increase in the radius would raise the predicted photospheric flux by about 70% and shrink the claimed excess from roughly 2.5× to about 1.5×. The H29α line independently supports the presence of ionized gas, so I do not view this as fatal, but the '2.5 times' figure in the abstract and conclusions should be qualified as model-dependent.
minor comments (6)
- [Abstract vs §3.2] The abstract describes the H29α line as having a 'symmetric profile', while §3.2 explicitly states that the data do not allow a conclusive determination of the line symmetry; the abstract should be revised to match this caveat.
- [Abstract vs §2/Table A.1] The abstract lists the observing frequencies as 'near 212 and 253 GHz', but the text and Table A.1 use 221.61 and 264.86 GHz for the continuum reference frequencies; please make the numbers consistent.
- [§3.1] The notation 'α = +1.26 ± 0.44 (3σ)' is ambiguous: if ±0.44 is a 1σ uncertainty, then the significance is about 2.9σ and should be stated as such; if it is a 3σ uncertainty, the relationship to the 0.4 systematic error should be explained. The abstract's '∼3σ error' makes this ambiguity worse.
- [Appendix D] The two Setup 1 observations on 2023-12-16 are only about one hour apart and share the same baseline configuration; treating them as independent data points in the mean and in the variability check may understate the uncertainty. Please quantify the correlation or merge them.
- [§3.2] The statement that the line is 'smaller in spatial extent than about 0.2″' is not supported by the quoted beams; please specify how this upper limit was derived (e.g., deconvolved Gaussian fit, uniform-weighting image) and report the effective beam used for the line map.
- [§3.1 and Abstract] Please check whether Matthews et al. (2023) already reported a spectral index for δ Cep; if so, the 'For the first time' claim in §3.1 and the abstract needs adjustment.
Circularity Check
No significant circularity: the ALMA measurements are new, the photospheric baseline uses external angular diameters and temperatures, and the H29α line is an independent detection.
full rationale
The paper's central claims rest on new ALMA data, not on a fitted parameter that is then renamed as a prediction. The predicted stellar continuum (Eq. 1) uses the uniform-disk angular diameter from PIONIER/VINCI (Anderson et al. 2016; Kervella et al. 2004) and effective temperatures from Luck (2018) and SPIPS interpolation; none of these inputs are fitted to the ALMA flux densities. The reported excess is measured flux minus this externally anchored baseline, so it is not circular by construction. Even if the photospheric baseline were revised upward, the H29α recombination line is an independent spectroscopic detection whose existence does not depend on the continuum prediction. The spectral index α = +1.26 ± 0.44 is measured directly from the two ALMA bands, with a stated systematic uncertainty from calibration; the comparison to the earlier Hocdé et al. (2021) model is explicitly falsifying, since the paper concludes that the model does not explain the millimeter data. Self-citations appear as motivation and as a model to be tested, but they are not load-bearing for the detection claim. The concerns raised about the Appendix D outlier epoch and the overstatement of line symmetry are robustness/correctness issues, not circularity. No equation reduces to its own input, no fitted parameter is relabeled a prediction, and no uniqueness theorem is imported from the authors' prior work.
Assumptions & free parameters
assumptions (6)
- domain assumption The stellar millimeter continuum follows the Rayleigh-Jeans blackbody law (Eq. 1) with literature Teff and UD angular diameter.
- domain assumption The millimeter continuum is formed at a similar photospheric radius as the optical/IR photosphere.
- standard math Free-free emission spectral index theory for ionized envelopes spans -0.1 to +2 depending on density gradient (Wright & Barlow 1975; Panagia & Felli 1975).
- standard math RRL line formation and the electron temperature formula (Gordon & Sorochenko 2002, Eq. 2.124) apply.
- domain assumption The H29α line identification at 256.302035 GHz is correct and the feature is not an artifact.
- domain assumption ALMA absolute flux calibration accuracy is ~10% (Francis et al. 2020).
Cite this review
Pith. "Pith review of Discovery of ionized circumstellar gas emission around the long-period Cepheid $\ell$ Carinae with ALMA." pith.science (2026). https://pith.science/paper/AEX6RO4G
@misc{pith2026250119011,
author = {Pith},
title = {Pith review of: Discovery of ionized circumstellar gas emission around the long-period Cepheid $\ell$ Carinae with ALMA},
year = {2026},
howpublished = {\url{https://pith.science/paper/AEX6RO4G}},
note = {Machine review of arXiv:2501.19011}
}
abstract
Cepheid circumstellar emissions have previously been detected via both infrared excess and infrared interferometric observations at a few stellar radii. Those studies have shown that these circumstellar emission can be produced by ionized gas, however there is no direct observational evidence to confirm this hypothesis. In this letter we explore the continuum emission and a spectrum of the bright and long-period Cepheid, $\ell$~Car ($P=35.56\,$day) at millimeter-wavelengths in order to detect possible effects of ionized gas emission. We presented ALMA observations of $\ell$ Car in two spectral setups in Band~6 (near 212 and 253\,GHz, respectively) and we compared the measured flux density to one expected for the stellar continuum. We also derived the spectral index and probed the presence of Radio Recombination Lines (RRL). We report statistically significant emission of about 3.5$\,$mJy in the two spectral ranges, which is about 2.5 times the stellar continuum emission. For the first time, we are also able to derive the spectral index of the flux density ($S_\nu \propto \nu^\alpha$), $\alpha=+1.26\pm$0.44 ($\sim$3$\sigma$ error), which is characteristic of partially optically thick ionized gas emission. Additionally, we discovered an emission line from a RRL of hydrogen H29$\alpha$ centered on the stellar rest velocity, smaller in spatial extent than about 0\farcs2 ($\lesssim 100\,$AU), with a symmetric profile with a width at half power of 55.3$\pm$7.5\,\kms (1$\sigma$ error). It confirms the presence of ionized gas emission near $\ell$~Car. The millimeter emission detected from $\ell$ Car can be attributed to ionized gas emission from the Cepheid's chromosphere. Further radio interferometric observations are necessary to confirm the occurrence of these ionized gas envelopes around Cepheids of different pulsation periods.
Figures
Reference graph
Works this paper leans on
-
[1]
, " * 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.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
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 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
-
[3]
Altenhoff , W. J., Thum , C., & Wendker , H. J. 1994, , 281, 161
work page 1994
-
[4]
I., M \'e rand , A., Kervella , P., et al
Anderson , R. I., M \'e rand , A., Kervella , P., et al. 2016, , 455, 4231
work page 2016
-
[5]
Anderson , R. I., Viviani , G., Shetye , S. S., et al. 2024, , 686, A177
work page 2024
-
[6]
Andrievsky , S. M., Kovtyukh , V. V., & Korotin , S. A. 2023, , 671, A27
work page 2023
-
[7]
M., Sip o cz , B
Astropy Collaboration , Price-Whelan , A. M., Sip o cz , B. M., et al. 2018, , 156, 123
2018
- [8]
Show all 64 references
-
[9]
& Love , S
Bohm-Vitense , E. & Love , S. G. 1994, , 420, 401
1994
-
[10]
& Seaton , M
Brocklehurst , M. & Seaton , M. J. 1972, , 157, 179
1972
-
[11]
2022, , 134, 114501
CASA Team , Bean , B., Bhatnagar , S., et al. 2022, , 134, 114501
2022
-
[12]
& Kurucz , R
Castelli , F. & Kurucz , R. L. 2003, in IAU Symposium, Vol. 210, Modelling of Stellar Atmospheres, ed. N. Piskunov , W. W. Weiss , & D. F. Gray , A20
2003
-
[13]
2022, , 511, 2125
Cs \"o rnyei , G., Szabados , L., Moln \'a r , L., et al. 2022, , 511, 2125
2022
-
[14]
Deasy , H. P. 1988, , 231, 673
1988
-
[15]
Dent , W. R. F., Harper , G. M., Richards , A. M. S., Kervella , P., & Matthews , L. D. 2024, , 966, L13
2024
-
[16]
G., Guinan , E
Engle , S. G., Guinan , E. F., Harper , G. M., et al. 2017, , 838, 67
2017
-
[17]
G., Guinan , E
Engle , S. G., Guinan , E. F., Harper , G. M., Neilson , H. R., & Remage Evans , N. 2014, , 794, 80
2014
-
[18]
2007, , 476, 73
Fouqu \'e , P., Arriagada , P., Storm , J., et al. 2007, , 476, 73
2007
-
[19]
R., & Harsono , D
Francis , L., Johnstone , D., Herczeg , G., Hunter , T. R., & Harsono , D. 2020, , 160, 270
2020
-
[20]
Freedman , W. L. & Madore , B. F. 2024, in IAU Symposium, Vol. 376, IAU Symposium, ed. R. de Grijs , P. A. Whitelock , & M. Catelan , 1--14
2024
-
[21]
Gaia Collaboration , Vallenari , A., Brown , A. G. A., et al. 2023, , 674, A1
2023
-
[22]
2013, , 558, A140
Gallenne , A., M \'e rand , A., Kervella , P., et al. 2013, , 558, A140
2013
-
[23]
1968, , 17, 445
Goldwire , Henry C., J. 1968, , 17, 445
1968
-
[24]
Gordon , M. A. & Sorochenko , R. L. 2002, Radio Recombination Lines. Their Physics and Astronomical Applications , Vol. 282
2002
-
[25]
2024, VizieR Online Data Catalog: Herschel/SPIRE point source catalog (HSPSC) (Herschel team+, 2017) , VizieR On-line Data Catalog: VIII/112
Herschel Team , Schulz , B., Marton , G., et al. 2024, VizieR Online Data Catalog: Herschel/SPIRE point source catalog (HSPSC) (Herschel team+, 2017) , VizieR On-line Data Catalog: VIII/112. Originally published in: Herschel catalogs (2017)
2017
-
[26]
2024, arXiv e-prints, arXiv:2501.00373
Hocd \'e , V., Matter , A., Nardetto , N., et al. 2024, arXiv e-prints, arXiv:2501.00373
2024 arXiv
-
[27]
2020 a , , 641, A74
Hocd \'e , V., Nardetto , N., Borgniet , S., et al. 2020 a , , 641, A74
2020
-
[28]
2020 b , , 633, A47
Hocd \'e , V., Nardetto , N., Lagadec , E., et al. 2020 b , , 633, A47
2020
-
[29]
2021, , 651, A92
Hocd \'e , V., Nardetto , N., Matter , A., et al. 2021, , 651, A92
2021
-
[30]
R., et al
Kervella , P., Gallenne , A., Evans , N. R., et al. 2019, , 623, A116
2019
-
[31]
2009, , 498, 425
Kervella , P., M \'e rand , A., & Gallenne , A. 2009, , 498, 425
2009
-
[32]
2006, , 448, 623
Kervella , P., M \'e rand , A., Perrin , G., & Coud \'e du Foresto , V. 2006, , 448, 623
2006
-
[33]
2004, , 416, 941
Kervella , P., Nardetto , N., Bersier , D., Mourard , D., & Coud \'e du Foresto , V. 2004, , 416, 941
2004
-
[34]
Kraft , R. P. 1957, , 125, 336
1957
-
[35]
Luck , R. E. 2018, , 156, 171
2018
-
[36]
2024, , 688, A203
Marton , G., Gezer , I., Madar \'a sz , M., et al. 2024, , 688, A203
2024
-
[37]
D., Evans , N
Matthews , L. D., Evans , N. R., & Rupen , M. P. 2020, in American Astronomical Society Meeting Abstracts, Vol. 235, American Astronomical Society Meeting Abstracts \#235, 106.03
2020
-
[38]
D., Evans , N
Matthews , L. D., Evans , N. R., & Rupen , M. P. 2023, , 165, 92
2023
-
[39]
D., Marengo , M., & Evans , N
Matthews , L. D., Marengo , M., & Evans , N. R. 2016, , 152, 200
2016
-
[40]
D., Marengo , M., Evans , N
Matthews , L. D., Marengo , M., Evans , N. R., & Bono , G. 2012, , 744, 53
2012
-
[41]
P., Kervella , P., et al
M \'e rand , A., Aufdenberg , J. P., Kervella , P., et al. 2007, , 664, 1093
2007
-
[42]
2015, , 584, A80
M \'e rand , A., Kervella , P., Breitfelder , J., et al. 2015, , 584, A80
2015
-
[43]
2006, , 453, 155
M \'e rand , A., Kervella , P., Coud \'e du Foresto , V., et al. 2006, , 453, 155
2006
-
[44]
H., Kraus , S., Millour , F., & Gillet , D
Nardetto , N., Groh , J. H., Kraus , S., Millour , F., & Gillet , D. 2008, , 489, 1263
2008
-
[45]
R., Ngeow , C.-C., Kanbur , S
Neilson , H. R., Ngeow , C.-C., Kanbur , S. M., & Lester , J. B. 2009, , 692, 81
2009
-
[46]
R., Ngeow , C.-C., Kanbur , S
Neilson , H. R., Ngeow , C.-C., Kanbur , S. M., & Lester , J. B. 2010, , 716, 1136
2010
-
[47]
J., van Duinen , R., et al
Neugebauer , G., Habing , H. J., van Duinen , R., et al. 2019, IRAS Point Source Catalog v2.1 (PSC) , NASA IPAC DataSet, IRSA4
2019
-
[48]
M., Brown , A., et al
O'Gorman , E., Harper , G. M., Brown , A., et al. 2015, , 580, A101
2015
-
[49]
M., Ohnaka , K., et al
O'Gorman , E., Harper , G. M., Ohnaka , K., et al. 2020, , 638, A65
2020
-
[50]
Olnon , F. M. 1975, , 39, 217
1975
-
[51]
& Felli , M
Panagia , N. & Felli , M. 1975, , 39, 1
1975
-
[52]
N., & Dullemond , C
Peters , T., Longmore , S. N., & Dullemond , C. P. 2012, , 425, 2352
2012
-
[53]
Reid , M. J. & Menten , K. M. 1997, , 476, 327
1997
-
[54]
G., Anand , G
Riess , A. G., Anand , G. S., Yuan , W., et al. 2024, , 962, L17
2024
-
[55]
Rodgers , A. W. & Bell , R. A. 1968, , 138, 23
1968
-
[56]
Sasselov , D. D. & Lester , J. B. 1994 a , , 423, 795
1994
-
[57]
Sasselov , D. D. & Lester , J. B. 1994 b , , 423, 777
1994
-
[58]
Schmidt , E. G. & Parsons , S. B. 1982, , 48, 185
1982
-
[59]
L., et al
Scowcroft , V., Seibert , M., Freedman , W. L., et al. 2016, , 459, 1170
2016
-
[60]
J., Price , S
Smith , B. J., Price , S. D., & Baker , R. I. 2004, , 154, 673
2004
-
[61]
2016, , 200, 1
Wedemeyer , S., Bastian , T., Braj s a , R., et al. 2016, , 200, 1
2016
-
[62]
Welch , D. L. & Duric , N. 1988, , 95, 1794
1988
-
[63]
& Thompson , A
Wootten , A. & Thompson , A. R. 2009, IEEE Proceedings, 97, 1463
2009
-
[64]
Wright , A. E. & Barlow , M. J. 1975, , 170, 41
1975
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.