The diversity of the circumgalactic medium around z = 0 Milky Way-mass galaxies from the Auriga simulations
Pith reviewed 2026-05-25 00:00 UTC · model grok-4.3
The pith
Simulations of Milky Way-mass galaxies show the circumgalactic medium varies by 3-4 orders of magnitude in column density even at similar masses.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L* galaxies is extremely diverse: column densities of commonly observed species span ~3-4 dex and their covering fractions range from ~5 to 90 per cent. Despite this diversity, we identify the following correlations: the covering fractions of hydrogen and metals positively correlate with stellar mass, the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample also表明
What carries the argument
Analysis of 28 isolated Milky Way-mass galaxy haloes from the Auriga cosmological magneto-hydrodynamical zoom-in simulation suite, used to compute column densities and covering fractions of gas species such as H I, C IV, and Si II in the circumgalactic medium.
If this is right
- Covering fractions of hydrogen and metals positively correlate with stellar mass.
- Covering fractions of H I, C IV, and Si II anticorrelate with AGN luminosity due to ionization effects.
- Covering fractions of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows.
- Long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.
Where Pith is reading between the lines
- Observational studies that rely on a limited number of single lines-of-sight will face systematic challenges when trying to reconstruct typical CGM properties.
- The results suggest that internal galaxy properties like disc fraction and AGN activity can drive CGM variation more than assembly history.
- Future work could test whether adding environmental factors beyond halo mass improves agreement between simulated and observed CGM diversity.
Load-bearing premise
The physical processes and feedback mechanisms implemented in the Auriga magneto-hydrodynamical simulations accurately capture the real-world behavior of the circumgalactic medium.
What would settle it
Multi-sightline absorption observations of several real Milky Way-mass galaxies at z=0 that find column density spreads much narrower than 3-4 dex across the sample.
read the original abstract
Galaxies are surrounded by massive gas reservoirs (i.e. the circumgalactic medium; CGM) which play a key role in their evolution. The properties of the CGM, which are dependent on a variety of internal and environmental factors, are often inferred from absorption line surveys which rely on a limited number of single lines-of-sight. In this work we present an analysis of 28 galaxy haloes selected from the Auriga project, a cosmological magneto-hydrodynamical zoom-in simulation suite of isolated Milky Way-mass galaxies, to understand the impact of CGM diversity on observational studies. Although the Auriga haloes are selected to populate a narrow range in halo mass, our work demonstrates that the CGM of L* galaxies is extremely diverse: column densities of commonly observed species span ~3-4 dex and their covering fractions range from ~5 to 90 per cent. Despite this diversity, we identify the following correlations: 1) the covering fractions (CF) of hydrogen and metals of the Auriga haloes positively correlate with stellar mass, 2) the CF of H I, C IV, and Si II anticorrelate with active galactic nucleus luminosity due to ionization effects, and 3) the CF of H I, C IV, and Si II positively correlate with galaxy disc fraction due to outflows populating the CGM with cool and dense gas. The Auriga sample demonstrates striking diversity within the CGM of L* galaxies, which poses a challenge for observations reconstructing CGM characteristics from limited samples, and also indicates that long-term merger assembly history and recent star formation are not the dominant sculptors of the CGM.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes the circumgalactic medium (CGM) in 28 Milky Way-mass haloes from the Auriga magneto-hydrodynamical zoom-in simulations. Despite the narrow halo-mass selection, it reports that column densities of commonly observed ions span ~3-4 dex and covering fractions range from ~5% to 90%. It identifies positive correlations of covering fractions with stellar mass and disc fraction, and anticorrelations with AGN luminosity, concluding that long-term merger history and recent star formation are not the dominant drivers of CGM properties.
Significance. If the Auriga subgrid model produces CGM statistics representative of real L* galaxies, the result is significant because it quantifies the intrinsic scatter that limited observational sightlines must contend with and supplies concrete correlations that can be tested against absorption surveys. The use of a controlled sample of 28 haloes at fixed mass is a clear strength for isolating diversity from mass dependence.
major comments (2)
- [Abstract, §1] Abstract and opening of §1: the central claim that 'the CGM of L* galaxies is extremely diverse' is presented as applying to observed systems, yet the entire result rests on the Auriga implementation (cooling, star-formation threshold, stellar/AGN feedback, metal diffusion). No quantitative comparison of the simulated column-density distributions or covering-fraction statistics to observational samples (e.g., COS-Halos or other quasar-absorption surveys) is provided; this comparison is load-bearing for the extrapolation beyond the model.
- [Results section on covering fractions] Results on covering fractions (figures and associated text): the reported ranges (5–90 per cent) are given without accompanying discussion of the number of sightlines per halo, the impact of projection effects, or convergence tests with respect to resolution or subgrid parameter variations. If these ranges are sensitive to the specific numerical choices, the claimed diversity could be partly numerical rather than physical.
minor comments (2)
- [Abstract] The abstract states the sample size (28) only after the diversity claim; moving the sample description earlier would improve clarity.
- Notation for column densities and covering fractions should be defined consistently on first use in the text and figures.
Simulated Author's Rebuttal
We thank the referee for their constructive comments on our manuscript. We respond point-by-point to the major comments below.
read point-by-point responses
-
Referee: [Abstract, §1] Abstract and opening of §1: the central claim that 'the CGM of L* galaxies is extremely diverse' is presented as applying to observed systems, yet the entire result rests on the Auriga implementation (cooling, star-formation threshold, stellar/AGN feedback, metal diffusion). No quantitative comparison of the simulated column-density distributions or covering-fraction statistics to observational samples (e.g., COS-Halos or other quasar-absorption surveys) is provided; this comparison is load-bearing for the extrapolation beyond the model.
Authors: The manuscript's primary goal is to quantify the diversity of CGM properties within a controlled sample of 28 Auriga haloes at fixed halo mass. The abstract and §1 language is intended to highlight the implications of this simulated diversity for observational studies of L* galaxies. We agree the wording could be tightened to avoid implying a direct claim about observed systems. In revision we will modify the abstract and opening of §1 to state explicitly that the Auriga simulations exhibit this diversity. A full quantitative comparison to COS-Halos or similar surveys is outside the scope of the present work, which focuses on internal scatter at fixed mass rather than model calibration; we therefore do not plan to add such a comparison. revision: partial
-
Referee: [Results section on covering fractions] Results on covering fractions (figures and associated text): the reported ranges (5–90 per cent) are given without accompanying discussion of the number of sightlines per halo, the impact of projection effects, or convergence tests with respect to resolution or subgrid parameter variations. If these ranges are sensitive to the specific numerical choices, the claimed diversity could be partly numerical rather than physical.
Authors: We will expand the relevant results section to specify the number of sightlines per halo used to compute covering fractions and to discuss projection effects. The Auriga project has already demonstrated resolution convergence for halo gas properties in earlier papers; we will add an explicit reference to those tests. Full convergence tests that vary subgrid parameters would require an entirely new simulation suite and are therefore beyond the scope of this study. revision: yes
Circularity Check
No circularity: direct analysis of simulation outputs
full rationale
The paper reports column-density spans and covering-fraction ranges measured directly from the 28 Auriga zoom-in runs. These quantities are computed from the simulation gas cells and sightlines; no parameter is fitted to a subset and then re-predicted, no quantity is defined in terms of itself, and no uniqueness theorem or ansatz is imported via self-citation to close the argument. The central claim is therefore an empirical statement about the chosen simulation suite rather than a derivation that reduces to its own inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Auriga simulations correctly model CGM physics at z=0 for MW-mass galaxies
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.