Pith. Sign in

REVIEW

Resolving the physics of Quasar Ly$\alpha$ Nebulae (RePhyNe): II. The dense and clumpy CGM of Quasars at $z \sim$ 3.5

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Titouan Lazeyras, Sebastiano Cantalupo, Andrea Travascio, Nicolas Ledos, Joop Schaye, Gabriele Pezzulli, Federico Baraggioni, Alessandro Lupi
show 5 more authors
Andrea V. Macciò, Marta Galbiati, Antonio Pensabene, Giada Quadri, Weichen Wang

Astrophysics of Galaxies · Submitted 1 Oct 2026 · 230 references

Pith's verdict: Solid qualitative comparison of Ly-alpha SB distributions vs. simulations, but the headline quantitative claims (C >= 8x, sigma >= 2.6, densities >= 10 cm^-3) rest on unpublished companion papers and an arbitrary median-density fix that make the specific numbers fragile; the qualitative dichotomy survives.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 2, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 6 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.01400 (version 1)
Pith ID3MIKIYOZ, Pith's permanent record of this paper and its review
CitationAlessandro Lupi, Andrea Travascio, Andrea V. Macci\`o, Antonio Pensabene, Federico Baraggioni, Gabriele Pezzulli, Giada Quadri, Joop Schaye, Marta Galbiati, Nicolas Ledos, Sebastiano Cantalupo, Titouan Lazeyras, Weichen Wang. "Resolving the physics of Quasar Ly$\alpha$ Nebulae (RePhyNe): II. The dense and clumpy CGM of Quasars at $z \sim$ 3.5." arXiv:2610.01400 (2026). https://pith.science/paper/3MIKIYOZ
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

This paper argues that the cold gas around quasars at z~3.5 is far denser and clumpier than any current cosmological simulation produces.

The reading

The paper tries to establish that current cosmological simulations fail to build the dense, clumpy cold gas that surrounds quasars at z~3.5. By comparing MUSE observations of quasar Ly$\alpha$ nebulae with mock observations from three simulation suites under maximal fluorescence, it finds that observed nebulae contain far more bright pixels — a covering fraction reaching 100% inside 100–200 ckpc — than simulations, implying gas densities of at least 10 cm$^{-3}$ must be present in quasar CGM. It then shows every simulated cold-gas density distribution, a skewed lognormal with clumping factor 16–102, is too narrow: the observed HeII/Ly$\alpha$ ratio requires a clumping factor at least 8 times larger at fixed median density. This matters because the cold CGM is the reservoir for star formation, so getting its small-scale structure wrong distorts our picture of accretion and feedback at high redshift.

What carries the argument

The Ly$\alpha$ surface brightness function — the count of pixels per logarithmic SB bin in the 100–200 ckpc annulus — used for the first time as a full distribution rather than a mean, since the emissivity under maximal fluorescence is set by the cold-gas density squared, $\epsilon_{\mathrm{Ly}\alpha} \propto n^2 \alpha_{\mathrm{eff}}(T)$. The clumping factor $C = \langle n^2 \rangle / \langle n \rangle^2$ connects SB to the unresolved density distribution, with $\mathrm{SB} = B(T)\, f_V\, C_{\mathrm{cold}}\,\langle n \rangle_{\mathrm{cold}}^2$ for a cold, low-filling-factor phase. The observed HeII/Ly$\alpha$ ratio supplies the independent, maximal-fluorescence-free constraint on the PDF width, converted into a required $\sigma$ via CLOUDY photoionization models in the companion paper.

Extended reading notes

Core claim

The central claim is that the average density, clumpiness, or volume filling factor of the cold Ly$\alpha$-emitting gas in the CGM — or all three — is too low in current simulations at $z\sim 3.5$. Under maximal fluorescence the Ly$\alpha$ surface brightness scales as $n^2$, so the shape of the SB distribution directly encodes the cold-gas density distribution; the observed SB function, peaking more than an order of magnitude above the detection limit and reaching full covering fraction within 100–200 ckpc, demands densities of at least 10 cm$^{-3}$ that only the highest-resolution COLIBRE run approaches. Independently of the maximal-fluorescence assumption, the measured HeII/Ly$\alpha$ ratio, modeled via CLOUDY with a lognormal density PDF, requires a width $\sigma \ge 2.6$ at fixed median density $\log_{10} n_0 = -2$, i.e. a clumping factor near 860, against the 16–102 found in simulations — at least a factor of 8 broader. The paper also reports that all simulated cold CGM density PDFs are skewed lognormals, the first systematic statement of that shape for the CGM, and that averaged SB quantities hide the discrepancy while distribution-level comparisons expose it.

Load-bearing premise

The quantitative claims — densities of at least 10 cm$^{-3}$ and a clumping factor at least 8 times the simulated value — rest on the HeII/Ly$\alpha$ photoionization modeling of unpublished companion papers and on fixing the median cold-gas density to $\log_{10} n_0 = -2$, a choice the authors themselves call arbitrary.

Editorial extensions

If this is right

  • Simulations that impose an effective equation of state at a density threshold near 10 cm$^{-3}$ should raise that threshold or drop the equation of state entirely so the CGM itself can form dense clumps.
  • Averaged comparisons — mean SB or circularly averaged profiles, which agree within a factor of 3 — cannot reveal the real discrepancy; distribution-level quantities like the SB function and covering fraction are required.
  • Matching the observations will need higher resolution combined with additional physics, likely including self-shielding prescriptions, cooling functions, and small-scale turbulence rather than resolution alone.
  • All simulated cold CGM phases are under-pressurized relative to the hot phase by factors of 5–18, offering a plausible origin for the missing dense clumps if that imbalance is physical.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The '8 times' figure and the '≥10 cm$^{-3}$' requirement are not independent: the $\sigma \ge 2.6$ bound follows from fixing $\log_{10} n_0 = -2$, which the authors themselves call arbitrary, and they note that fixing the simulated width instead forces $\log_{10} n_0 > -1$ with mean density above 10 cm$^{-3}$ — so the tension is real but its magnitude hinges on a normalization the paper flags as u
  • The COLIBRE m5 exception, at the same resolution as DaLya HR which shows no improvement, suggests resolution-dependent subgrid physics — self-shielding and cooling permitting gas above ~0.1 cm$^{-3}$ — matters more than raw resolution, an implication the paper hints at but does not fully resolve.
  • A testable extension would be to re-derive the required $(n_0, \sigma)$ pair jointly across the full halo mass range $10^{12}{-}10^{12.5}\,M_\odot$ to see whether any (median density, width) combination simultaneously reproduces both the Ly$\alpha$ SB function and the HeII/Ly$\alpha$ ratio, bracketing the clumpiness the CGM actually needs.
Original abstract

The small scale density distribution of the cold circumgalactic medium (CGM) plays a central role in galaxy evolution, yet it is still poorly understood. We provide new constraints on the cold CGM density distribution by comparing the CGM emission from MUSE observations of quasar (QSO) nebulae at $z \sim 3.8$ with mock observations generated from the new DaLya simulation suite under the assumption of maximal fluorescence. These are further complemented by mock observations generated from the COLIBRE (m5 and m6 resolutions) and HELLO simulations. We first investigate the distribution of the observed and simulated Ly$\alpha$ surface brightness (SB) values, which directly depend on the cold gas density distribution under our assumption, finding stark differences between observations and upper limits provided by simulations with the exception of the COLIBRE higher-resolution run (m5). The observations generally show a much larger covering fraction of regions with relatively large Ly$\alpha$ SB values and suggest that densities reaching at least 10 cm$^{-3}$ must be present in the CGM of QSOs within a radial distance of 50 to 200 ckpc. To understand the origin of this discrepancy, we study the cold CGM density distribution finding in all cases simulated PDFs consistent with skewed lognormal ones with medians $\sim 10^{-2}$ cm$^{-3}$ and clumping factors between 16 and 102. We compare the broadness of these distributions to the one implied by the observed HeII/Ly$\alpha$ CGM ratio, which is independent on the assumption of maximum fluorescence. Assuming lognormal PDFs with median values similar to the simulations, the observed ratio requires a broader PDF with a clumping factor at least 8 times larger than predicted by simulations. These results suggest that higher resolution and/or additional physical mechanisms are needed to reproduce the observed CGM of QSOs at these redshifts.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-02: 1 request, 7,320 input tokens.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.87 · human participants 0.01 · controlled comparison 0.28 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.86 · mentions code or data availability 0.04.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

The small scale density distribution of the cold circumgalactic medium (CGM) plays a central role in galaxy evolution, yet it is still poorly understood. We provide new constraints on the cold CGM density distribution by comparing the CGM emission from MUSE observations of quasar (QSO) nebulae at $z \sim 3.8$ with mock observations generated from the new DaLya simulation suite under the assumption of maximal fluorescence. These are further complemented by mock observations generated from the COLIBRE (m5 and m6 resolutions) and HELLO simulations. We first investigate the distribution of the observed and simulated Ly$\alpha$ surface brightness (SB) values, which directly depend on the cold gas density distribution under our assumption, finding stark differences between observations and upper limits provided by simulations with the exception of the COLIBRE higher-resolution run (m5). The observations generally show a much larger covering fraction of regions with relatively large Ly$\alpha$ SB values and suggest that densities reaching at least 10 cm$^{-3}$ must be present in the CGM of QSOs within a radial distance of 50 to 200 ckpc. To understand the origin of this discrepancy, we study the cold CGM density distribution finding in all cases simulated PDFs consistent with skewed lognormal ones with medians $\sim 10^{-2}$ cm$^{-3}$ and clumping factors between 16 and 102. We compare the broadness of these distributions to the one implied by the observed HeII/Ly$\alpha$ CGM ratio, which is independent on the assumption of maximum fluorescence. Assuming lognormal PDFs with median values similar to the simulations, the observed ratio requires a broader PDF with a clumping factor at least 8 times larger than predicted by simulations. These results suggest that higher resolution and/or additional physical mechanisms are needed to reproduce the observed CGM of QSOs at these redshifts.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 The small scale density distribution of the cold circumgalactic medium (CGM) plays a central role in galaxy evolution, yet it is still poorly understood. background 0.64 0.04 0.02 asserted 1.00
s1 We provide new constraints on the cold CGM density distribution by comparing the CGM emission from MUSE observations of quasar (QSO) nebulae at $z \sim 3.8$ with mock observations generated from the new DaLya simulation suite under the assumption of maximal fluorescence. contribution 0.63 0.95 0.06 asserted 1.00
s2 These are further complemented by mock observations generated from the COLIBRE (m5 and m6 resolutions) and HELLO simulations. method 0.98 0.79 0.02 asserted 0.96
s3 We first investigate the distribution of the observed and simulated Ly$\alpha$ surface brightness (SB) values, which directly depend on the cold gas density distribution under our assumption, finding stark differences between observations and upper limits provided by simulations with the exception of the COLIBRE higher-resolution run (m5). result 0.94 0.95 0.11 asserted 1.00
s4 The observations generally show a much larger covering fraction of regions with relatively large Ly$\alpha$ SB values and suggest that densities reaching at least 10 cm$^{-3}$ must be present in the CGM of QSOs within a radial distance of 50 to 200 ckpc. result 0.99 0.90 0.38 hedged possibility 0.85
s5 To understand the origin of this discrepancy, we study the cold CGM density distribution finding in all cases simulated PDFs consistent with skewed lognormal ones with medians $\sim 10^{-2}$ cm$^{-3}$ and clumping factors between 16 and 102. result 0.94 0.95 0.10 asserted 0.96
s6 We compare the broadness of these distributions to the one implied by the observed HeII/Ly$\alpha$ CGM ratio, which is independent on the assumption of maximum fluorescence. method 0.91 0.80 0.05 asserted 0.90
s7 Assuming lognormal PDFs with median values similar to the simulations, the observed ratio requires a broader PDF with a clumping factor at least 8 times larger than predicted by simulations. result 0.94 0.86 0.29 asserted 0.70
s8 These results suggest that higher resolution and/or additional physical mechanisms are needed to reproduce the observed CGM of QSOs at these redshifts. interpretation 0.78 0.47 0.68 hedged possibility 1.00

What kind of paper the classifier reports this is: empirical study 0.57 · evidence from observation or experiment 0.44.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. A visual impression of the DaLya simulation (fiducial). The center panel shows the full base simulation with the two re-simulated regions highlighted. The inserts on the left and right present different quantities for two halos, one in each region, in a 2x2x2 cMpc volume. From top to bottom: projected dark matter density, gas column density, cold gas column density, and the Lyα nebula as extracted by CubEx (Cantalupo et al. 2019) from the mock cube generated using the P2C output (see section 3.1 for details). and mb = 9.0·104M⊙. We refer to these as DaLya HR. DaLya HR contains only…
Figure 2. Mosaic of selected mock pseudo-narrow band Lyα SB maps (without noise for clarity) at z = 3.5 in a 1×1 cMpc field of view for DaLya fiducial. Each panel shows a randomly selected LOS for halos spanning the mass range 12.0 ≤ log M ≤ 12.5. The white circles indi￾cate the radial annulus 100 – 200 ckpc. In addition to ionizing the CGM, the central quasar also heats the surrounding medium through photo-electric heating. Due to the gas cooling when it recombines and emits radiation, the tem￾perature of the gas will eventually reach an equilibrium between photo-heating and cooling. We…
Figure 3. Left panel: The Lyα SB function (i.e. the count of intrinsic SB values per logarithmic SB bins) at z = 3.5 in the 100-200 ckpc annulus and comparison with observations from Travascio et al. in prep.. Results for DaLya fiducial are shown in blue while results from observations are in black. Each time we show the mean number of pixels per LOS. The dashed vertical lines indicate the mean of each histogram, and the dotted vertical line indicates the average detection limit in observations. Central panel: The covering fraction of Lyα pixels, corresponding to the cumulative version of the…
Figure 4. Average PDFs of the density of the cold (103 ≤ T ≤ 105K) CGM in DaLya fiducial at z = 3.5 projected in two radial annuli. The solid histograms show the density distributions of simulations particles. For the 50-160ckpc annulus, we quote the median density n0 and the 1σ equivalent width of the distribution of lnnc (see appendix D). The solid gray line shows the inferred lognormal PDF from observations of the HeII/Lyα ratio from Travascio et al. in prep. in the 50-160ckpc radial bin when fixing log10 n0 = −2 and the corresponding minimum width σ = 2.6 compatible with the observed ratio
Figure 5. Comparison of the Lyα SB function with results from COL￾IBRE m6 (resolution similar to DaLya fiducial) and from other sim￾ulations at 8 times higher resolution: DaLya HR, COLIBRE m5, and HELLO. DaLya , COLIBRE and HELLO are presented in blue, red and green respectively. The linestyle (solid or dashed) indicates the resolu￾tion (∼ 106 and ∼ 105M⊙ respectively). All simulations with the excep￾tion of COLIBRE m5 present PDF shapes which are strikingly different than the observed one, i.e. they predict a larger number of faint pixels than observed. We stress again that the simulated SB…
Figure 6. Same as figure 4 but comparing to DaLya HR (left panel), COLIBRE m6 and m5 (central panel), and HELLO (right panel), focusing on the 50–160 ckpc annulus only. Again, n0 denotes the median density while σ is the 1σ equivalent width of ln nc . The corresponding mean density and clumping factor are listed in table 2. The solid gray line shows the inferred PDF from the observations of the HeII/Lyα ratio from Travascio et al. in prep. in the same radial bin when fixing log10 n0 = −2. The grey arrows illustrate that these are lower limit distributions. Sim. log10 n0 σ n¯ [cm−3 ] Ccold…
Figure 7. The temperature – density diagrams (mass weighted) at z = 3.5 for particles within the 50 –160 kpc annulus in each simulations. The color-coding indicates the number density of particles in each T–n bin. We focus on the temperature range 3 ≤ log T[K] ≤ 5. The black solid line indicates T=104K for reference. cooler and denser in m5 than m6. Additional differences in the cooling prescription bring the median temperature of dense gas (n ≳ 1cm−3 ) below 104K in COLIBRE. While a complete pic￾ture of the impact of resolution and different physical mecha￾nisms on the Lyα SB values would…

Figures reproduced from arXiv:2610.01400 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.