REVIEW 3 major objections 2 minor 1 cited by
Quasars with the most massive molecular gas reservoirs show the dimmest central Lyα nebulae, while bright central Lyα sources lack CO detections.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.3
2026-07-01 01:43 UTC pith:XBI7K5XT
load-bearing objection New ALMA CO data on 37 MUSE quasars at z~3 shows an anti-correlation with central Lyα brightness, but the non-detection interpretation carries the main uncertainty. the 3 major comments →
ALMA visits the QSO MUSEUM: connecting molecular gas and the cool circumgalactic medium around 37 z~3 quasars
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Among the 37 quasars, those with the highest molecular gas masses are found with the centrally dimmest Lyα nebulae, whereas quasars hosting the centrally brightest Lyα nebulae are generally undetected in CO(4-3). This pattern indicates that the molecular gas and dust in the quasar hosts regulate Lyα escape and thereby shape the observed emission from the circumgalactic medium.
What carries the argument
ALMA CO(4-3) line observations that measure molecular gas masses and compare them directly to the central surface brightness of the Lyα nebulae mapped by MUSE.
Load-bearing premise
That the absence of CO detection means the quasar truly has little molecular gas rather than the line being missed because of sensitivity, excitation, or dust effects that differ between the bright and dim Lyα groups.
What would settle it
Finding strong CO(4-3) emission in several of the quasars that currently show the brightest central Lyα nebulae would remove the reported anti-correlation.
If this is right
- Quasars with Eddington ratios below about 0.9 retain larger molecular gas reservoirs.
- Quasars accreting near or above Eddington ratio 0.9 appear to have depleted their gas through outflows.
- Even the CO-detected low-Eddington quasars show gas fractions around 0.10, lower than typical for inactive star-forming galaxies at similar redshifts.
- The six marginally resolved CO sources extend up to 8 kpc, and 14 companion galaxies are detected, pointing to overdense fields around the quasars.
Where Pith is reading between the lines
- If the anti-correlation holds, models of AGN feedback must include how host gas content modulates the visibility of halo Lyα emission.
- Follow-up observations at higher sensitivity or different CO transitions could test whether the non-detections in bright Lyα systems are truly gas-poor.
- The reported cross-correlation length of nearly 10 comoving megaparsecs suggests these quasars sit in the same large-scale structures as their CO companions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents ALMA CO(4-3) observations of 37 z~3 quasars from the QSO MUSEUM survey that were previously mapped in Lyα with MUSE. 21 sources are detected in CO with derived molecular gas masses M_gas ≈ (3-40)×10^9 M_⊙. The central claim is a qualitative anti-correlation: quasars with the most massive molecular gas reservoirs show the centrally dimmest Lyα nebulae, while those with the brightest central Lyα are generally CO non-detections; this is interpreted as evidence that host gas and dust regulate Lyα escape into the halo. Secondary results include a trend with Eddington ratio, low gas fractions (median M_gas/M_* ~0.10) even for CO-detected low-Eddington systems, six marginally resolved CO sources, and 14 high-fidelity companions implying overdense fields with a measured cross-correlation length.
Significance. If the anti-correlation survives quantitative scrutiny and bias checks, the result would furnish a direct observational connection between the cold molecular ISM of quasar hosts and the cool CGM traced by extended Lyα at cosmic noon, with implications for AGN feedback and Lyα radiative transfer. The sample size (37 objects), joint ALMA+MUSE coverage, and companion detections are concrete strengths; the reported cross-correlation length is a falsifiable, quantitative prediction.
major comments (3)
- [Abstract and results section on CO detections/non-detections] The central claim (abstract and §4) that the 16 CO(4-3) non-detections correspond to intrinsically low molecular gas reservoirs (M_gas ≪ 3×10^9 M_⊙) permitting bright central Lyα escape is load-bearing, yet the manuscript provides no quantitative demonstration that the ALMA upper limits are uniform across the sample or that excitation conditions (e.g., r_{4-3/1-0}) do not differ systematically between the bright-Lyα and dim-Lyα subsets. If sensitivity, beam dilution, or excitation biases correlate with Lyα central brightness, the inferred physical link is undermined.
- [Abstract and §4 (correlation analysis)] No error bars are reported on the derived gas masses, no statistical measure of the anti-correlation significance (Spearman rank, p-value, or survival analysis accounting for upper limits) is given, and no completeness or selection-bias corrections are discussed. These omissions make it impossible to assess whether the claimed trend is robust or could arise from the QSO MUSEUM selection function.
- [§4 (Eddington ratio trends)] The statement that strongly accreting quasars (λ_Edd ≳ 0.9) deplete their gas via outflows is presented without supporting kinematic evidence or comparison to control samples; the low gas fractions for CO-detected low-Eddington systems are noted but not placed in context with matched inactive galaxies at the same stellar mass and redshift.
minor comments (2)
- [Throughout] Notation for molecular gas mass (M_gas) and gas fraction should be defined once with explicit conversion factors and assumptions (e.g., α_CO, r_{4-3/1-0}) rather than repeated in the abstract and multiple sections.
- [Results on resolved sources] The effective radii of the six marginally resolved CO sources (~8 kpc) would benefit from a table listing individual sizes, beam sizes, and S/N to allow assessment of resolution claims.
Simulated Author's Rebuttal
We thank the referee for the constructive report. We agree that several aspects of the central claim require quantitative support and will revise the manuscript to address the points raised. Below we respond point-by-point.
read point-by-point responses
-
Referee: [Abstract and results section on CO detections/non-detections] The central claim (abstract and §4) that the 16 CO(4-3) non-detections correspond to intrinsically low molecular gas reservoirs (M_gas ≪ 3×10^9 M_⊙) permitting bright central Lyα escape is load-bearing, yet the manuscript provides no quantitative demonstration that the ALMA upper limits are uniform across the sample or that excitation conditions (e.g., r_{4-3/1-0}) do not differ systematically between the bright-Lyα and dim-Lyα subsets. If sensitivity, beam dilution, or excitation biases correlate with Lyα central brightness, the inferred physical link is undermined.
Authors: We agree that the current text does not explicitly verify uniformity of the upper limits or test for systematic excitation differences. In the revision we will add a supplementary table listing the 3σ upper limits for all non-detections (converted to M_gas assuming the same r_{4-3/1-0}=0.5 and α_CO=0.8 used for detections) together with the rms noise and beam size for each target. We will also include a short discussion noting that the ALMA integration times were chosen to reach a uniform sensitivity goal and that no trend between rms and Lyα central surface brightness is present in the data; we will flag the assumption on excitation as a caveat. revision: yes
-
Referee: [Abstract and §4 (correlation analysis)] No error bars are reported on the derived gas masses, no statistical measure of the anti-correlation significance (Spearman rank, p-value, or survival analysis accounting for upper limits) is given, and no completeness or selection-bias corrections are discussed. These omissions make it impossible to assess whether the claimed trend is robust or could arise from the QSO MUSEUM selection function.
Authors: We acknowledge these omissions. The revised manuscript will report 1σ uncertainties on all M_gas values (propagating the CO flux errors and the 0.3 dex systematic uncertainty on α_CO). We will add a quantitative assessment of the anti-correlation using a Spearman rank test on the detected sources and a survival-analysis treatment (e.g., Kaplan–Meier or censored Spearman) that incorporates the upper limits. A brief section on possible selection biases will be included, noting that the parent QSO MUSEUM sample was selected on Lyα properties and that the ALMA follow-up was not biased by expected gas mass. revision: yes
-
Referee: [§4 (Eddington ratio trends)] The statement that strongly accreting quasars (λ_Edd ≳ 0.9) deplete their gas via outflows is presented without supporting kinematic evidence or comparison to control samples; the low gas fractions for CO-detected low-Eddington systems are noted but not placed in context with matched inactive galaxies at the same stellar mass and redshift.
Authors: The claim is interpretive and rests only on the observed anti-correlation between λ_Edd and M_gas within the sample. We will rephrase the text to “suggest” rather than “likely deplete” and will add a short comparison to literature gas fractions for star-forming galaxies at z≈3 and similar M_* (citing Tacconi et al. 2018 and similar works). No kinematic evidence for outflows is present in the current ALMA data cubes, so we cannot strengthen that part of the argument; the statement will be presented as a possible interpretation rather than a firm conclusion. revision: partial
Circularity Check
No circularity: central claim is direct empirical comparison of independent ALMA and MUSE datasets
full rationale
The paper reports observational associations between CO(4-3) detections/non-detections (yielding M_gas estimates) and Lyα nebula properties measured from separate VLT/MUSE data. No equations, fitted parameters, or self-citations are invoked to derive the main trend; the statement that massive reservoirs correlate with dimmest central Lyα (and bright Lyα with non-detections) is presented as a direct data comparison without reduction to inputs by construction. Assumptions about non-detections are interpretive and do not create definitional or fitted-input circularity. The derivation chain is self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption CO(4-3) luminosity can be converted to molecular hydrogen mass using a standard conversion factor appropriate for high-redshift quasar hosts
- domain assumption Non-detections in CO(4-3) correspond to genuinely lower molecular gas content rather than observational or excitation biases
read the original abstract
Extended Ly$\alpha$ emission is ubiquitous around quasars and traces cool circumgalactic gas, providing insight into halo gas dynamics and active galactic nuclei (AGN) feedback. However, its connection to the cold molecular gas of the host galaxies remains largely unexplored. We aim to characterize the molecular gas reservoirs in quasars at cosmic noon and investigate how they are linked to extended Ly$\alpha$ emission. To this end, we present ALMA CO(4-3) observations of 37 quasars at $z\sim3$ from the QSO MUSEUM survey, previously mapped in Ly$\alpha$ with VLT/MUSE. We derive molecular gas masses and gas fractions, explore correlations with Ly$\alpha$ nebula and quasar properties, and search for CO-emitting companions. Of 37 quasars, 21 are detected in CO(4-3), with gas masses $M_\mathrm{gas}\approx(3-40) \times10^9\,\mathrm{M_\odot}$. Quasars with the most massive molecular gas reservoirs are associated with the centrally dimmest Ly$\alpha$ nebulae, while those hosting the centrally brightest Ly$\alpha$ nebulae are generally not detected in CO. This suggests that gas and dust in the hosts regulate Ly$\alpha$ escape and consequently affect the emission from halo gas. We find evidence that quasars with lower Eddington ratios harbor more massive gas reservoirs, whereas strongly accreting quasars ($\lambda_\mathrm{Edd} \gtrapprox 0.9$) likely deplete their gas through quasar-driven outflows. Despite their higher molecular gas masses within the sample, CO-detected low-Eddington quasars exhibit low gas fractions, with a median $M_\mathrm{gas}/M_* \sim 0.10$, below those typical of inactive star-forming galaxies. Six quasars are marginally resolved in CO, with effective radii up to $\sim 8\,\mathrm{kpc}$. In addition, we detect 14 high-fidelity companion galaxies, indicating overdense quasar fields with a quasar-galaxy cross-correlation length of $9.81^{+2.22}_{-2.05}\,h^{-1}\mathrm{cMpc}$.
Figures
Forward citations
Cited by 1 Pith paper
-
Even central galaxies feel their environment: cosmic siphoning of cool gas accretion
Isolated central galaxies near massive halos are systematically more gas-poor than group centrals of the same black hole mass and star-formation offset, evidence for environment-driven gas removal.
Reference graph
Works this paper leans on
-
[1]
Akritas, M. G. & Siebert, J. 1996, MNRAS, 278, 919 Arrigoni Battaia, F., Chen, C.-C., Liu, H.-Y . B., et al. 2022, ApJ, 930, 72 Arrigoni Battaia, F., Hennawi, J. F., Prochaska, J. X., & Cantalupo, S. 2015, ApJ, 809, 163 Article number, page 16 Jelena Ritter et al.: Molecular gas reservoirs ofz∼3 quasars Arrigoni Battaia, F., Hennawi, J. F., Prochaska, J. ...
work page 1996
-
[2]
2010, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol
Bacon, R., Accardo, M., Adjali, L., et al. 2010, in Society of Photo-Optical In- strumentation Engineers (SPIE) Conference Series, V ol. 7735, Ground-based and Airborne Instrumentation for Astronomy III, ed. I. S. McLean, S. K. Ram- say, & H. Takami, 773508
work page 2010
-
[3]
Banerji, M., Carilli, C. L., Jones, G., et al. 2017, MNRAS, 465, 4390
work page 2017
-
[4]
Banerji, M., Jones, G. C., Wagg, J., et al. 2018, MNRAS, 479, 1154
work page 2018
-
[5]
Banerji, M., McMahon, R. G., Hewett, P. C., et al. 2012, MNRAS, 427, 2275
work page 2012
- [6]
-
[7]
Becker, R. H., White, R. L., & Helfand, D. J. 1994, in Astronomical Society of the Pacific Conference Series, V ol. 61, Astronomical Data Analysis Software and Systems III, ed. D. R. Crabtree, R. J. Hanisch, & J. Barnes, 165
work page 1994
-
[8]
Bertola, E., Circosta, C., Ginolfi, M., et al. 2024, A&A, 691, A178
work page 2024
-
[9]
Bischetti, M., Feruglio, C., Piconcelli, E., et al. 2021, A&A, 645, A33
work page 2021
-
[10]
Bolatto, A. D., Wolfire, M., & Leroy, A. K. 2013, ARA&A, 51, 207
work page 2013
- [11]
- [12]
- [13]
- [14]
-
[15]
2017, ApJ, 837, 71 Calistro Rivera, G., Alexander, D
Cai, Z., Fan, X., Yang, Y ., et al. 2017, ApJ, 837, 71 Calistro Rivera, G., Alexander, D. M., Rosario, D. J., et al. 2021, A&A, 649, A102
work page 2017
-
[16]
Cantalupo, S., Arrigoni-Battaia, F., Prochaska, J. X., Hennawi, J. F., & Madau, P. 2014, Nature, 506, 63
work page 2014
-
[17]
Cantalupo, S., Porciani, C., Lilly, S. J., & Miniati, F. 2005, ApJ, 628, 61
work page 2005
-
[18]
Carilli, C. L. & Walter, F. 2013, ARA&A, 51, 105
work page 2013
-
[19]
Chartas, G., Davidson, E., Brusa, M., et al. 2020, MNRAS, 496, 598
work page 2020
-
[20]
Chen, C.-C., Arrigoni Battaia, F., Emonts, B. H. C., Lehnert, M. D., & Prochaska, J. X. 2021, ApJ, 923, 200
work page 2021
-
[21]
Cicone, C., Mainieri, V ., Circosta, C., et al. 2021, A&A, 654, L8
work page 2021
-
[22]
Circosta, C., Mainieri, V ., Lamperti, I., et al. 2021, A&A, 646, A96
work page 2021
-
[23]
Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693
work page 1998
- [24]
- [25]
-
[26]
Costa, T., Rosdahl, J., Sijacki, D., & Haehnelt, M. G. 2018b, MNRAS, 479, 2079 da Cunha, E., Walter, F., Smail, I. R., et al. 2015, ApJ, 806, 110 D’Amato, Q., Gilli, R., Vignali, C., et al. 2020, A&A, 636, A37
work page 2079
-
[27]
Y ., Seebeck, J., Veilleux, S., et al
Dan, K. Y ., Seebeck, J., Veilleux, S., et al. 2025, ApJ, 979, 68
work page 2025
-
[28]
I., Baes, M., Bianchi, S., et al
Davies, J. I., Baes, M., Bianchi, S., et al. 2017, PASP, 129, 044102
work page 2017
-
[29]
2024, A&A, 689, A263 de Beer, S., Cantalupo, S., Travascio, A., et al
Davies, R., Shimizu, T., Pereira-Santaella, M., et al. 2024, A&A, 689, A263 de Beer, S., Cantalupo, S., Travascio, A., et al. 2023, MNRAS, 526, 1850 De Breuck, C., Downes, D., Neri, R., et al. 2005, A&A, 430, L1 De Breuck, C., Neri, R., & Omont, A. 2003, New A Rev., 47, 285
work page 2024
- [30]
-
[31]
Decarli, R., Falomo, R., Treves, A., et al. 2010, MNRAS, 402, 2453
work page 2010
-
[32]
Decarli, R., Pensabene, A., Venemans, B., et al. 2022, A&A, 662, A60
work page 2022
-
[33]
Decarli, R., Walter, F., Gónzalez-López, J., et al. 2019, ApJ, 882, 138
work page 2019
- [34]
- [35]
-
[36]
Dekel, A. & Birnboim, Y . 2006, MNRAS, 368, 2 Di Matteo, T., Springel, V ., & Hernquist, L. 2005, Nature, 433, 604 Díaz-Santos, T., Assef, R. J., Blain, A. W., et al. 2018, Science, 362, 1034
work page 2006
- [37]
-
[38]
Saas-Fee Lecture Notes: Physics of Lyman Alpha Radiative Transfer
Dijkstra, M. 2017, arXiv e-prints, arXiv:1704.03416
work page internal anchor Pith review Pith/arXiv arXiv 2017
- [39]
- [40]
- [41]
-
[42]
Drake, A. B., Neeleman, M., Venemans, B. P., et al. 2022, ApJ, 929, 86
work page 2022
-
[43]
Emonts, B. H. C., Cai, Z., Prochaska, J. X., Li, Q., & Lehnert, M. D. 2019, ApJ, 887, 86
work page 2019
-
[44]
Emonts, B. H. C., Lehnert, M. D., Lebowitz, S., et al. 2023, ApJ, 952, 148
work page 2023
-
[45]
Fabian, A. C. 2012, ARA&A, 50, 455
work page 2012
- [46]
-
[47]
P., Arrigoni-Battaia, F., Costa, T., et al
Farina, E. P., Arrigoni-Battaia, F., Costa, T., et al. 2019, ApJ, 887, 196 Faucher-Giguère, C.-A., Kereš, D., Dijkstra, M., Hernquist, L., & Zaldarriaga, M. 2010, ApJ, 725, 633
work page 2019
-
[48]
Feruglio, C., Fiore, F., Carniani, S., et al. 2018, A&A, 619, A39
work page 2018
-
[49]
Feruglio, C., Fiore, F., Carniani, S., et al. 2015, A&A, 583, A99
work page 2015
-
[50]
Feruglio, C., Maiolino, R., Piconcelli, E., et al. 2010, A&A, 518, L155
work page 2010
-
[51]
Fiore, F., Feruglio, C., Shankar, F., et al. 2017, A&A, 601, A143
work page 2017
-
[52]
Fluetsch, A., Maiolino, R., Carniani, S., et al. 2021, MNRAS, 505, 5753
work page 2021
-
[53]
Flury, S. R., Jaskot, A. E., Ferguson, H. C., et al. 2022, ApJ, 930, 126
work page 2022
-
[54]
Fogasy, J., Knudsen, K. K., Drouart, G., Lagos, C. D. P., & Fan, L. 2020, MN- RAS, 493, 3744
work page 2020
- [55]
-
[56]
2012, MNRAS, 424, 2276 García-Vergara, C., Hennawi, J
Fumagalli, M., Dessauges-Zavadsky, M., Furniss, A., et al. 2012, MNRAS, 424, 2276 García-Vergara, C., Hennawi, J. F., Barrientos, L. F., & Arrigoni Battaia, F. 2019, ApJ, 886, 79 García-Vergara, C., Hennawi, J. F., Barrientos, L. F., & Rix, H.-W. 2017, ApJ, 848, 7 García-Vergara, C., Rybak, M., Hodge, J., et al. 2022, ApJ, 927, 65
work page 2012
- [57]
- [58]
-
[59]
Genzel, R., Price, S. H., Übler, H., et al. 2020, ApJ, 902, 98 González Lobos, J., Arrigoni Battaia, F., Chang, S.-J., et al. 2023, A&A, 679, A41 González Lobos, J., Arrigoni Battaia, Fabrizio, Obreja, Aura, et al. 2026, A&A, 707, A380 González-López, J., Bauer, F. E., Aravena, M., et al. 2017, A&A, 608, A138 González-López, J., Decarli, R., Pavesi, R., e...
work page 2020
- [60]
-
[61]
Gronke, M. & Oh, S. P. 2018, MNRAS, 480, L111
work page 2018
-
[62]
Gronke, M. & Oh, S. P. 2020, MNRAS, 492, 1970
work page 2020
- [63]
- [64]
-
[65]
Heckman, T. M., Lehnert, M. D., Miley, G. K., & van Breugel, W. 1991, ApJ, 381, 373
work page 1991
-
[66]
Hennawi, J. F., Prochaska, J. X., Cantalupo, S., & Arrigoni-Battaia, F. 2015, Science, 348, 779
work page 2015
-
[67]
Herrera-Camus, R., Tacconi, L., Genzel, R., et al. 2019, ApJ, 871, 37
work page 2019
-
[68]
Herwig, E., Arrigoni Battaia, F., González Lobos, J., et al. 2024, A&A, 691, A210
work page 2024
-
[69]
Herwig, E., Battaia, F. A., Chen, C.-C., et al. 2025, A&A, 703, A17
work page 2025
-
[70]
Hill, R., Chapman, S. C., Scott, D., et al. 2019, MNRAS, 485, 753
work page 2019
-
[71]
Hopkins, P. F., Hernquist, L., Cox, T. J., & Kereš, D. 2008, ApJS, 175, 356
work page 2008
-
[72]
Hunter, T. R., Indebetouw, R., Brogan, C. L., et al. 2023, PASP, 135, 074501
work page 2023
-
[73]
J., Smail, I., Amblard, A., et al
Ivison, R. J., Smail, I., Amblard, A., et al. 2012, MNRAS, 425, 1320
work page 2012
-
[74]
Jarvis, M. E., Harrison, C. M., Mainieri, V ., et al. 2020, MNRAS, 498, 1560
work page 2020
-
[75]
C., Maiolino, R., Circosta, C., et al
Jones, G. C., Maiolino, R., Circosta, C., et al. 2023, MNRAS, 518, 691
work page 2023
- [76]
- [77]
-
[78]
Kinney, A. L., Schmitt, H. R., Clarke, C. J., et al. 2000, ApJ, 537, 152
work page 2000
-
[79]
Kormendy, J. & Ho, L. C. 2013, ARA&A, 51, 511
work page 2013
-
[80]
J., Strittmatter, B., Lamperti, I., et al
Koss, M. J., Strittmatter, B., Lamperti, I., et al. 2021, ApJS, 252, 29
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.