REVIEW 3 major objections 4 minor 1 cited by
The contribution of stars, dust, neutral gas and supermassive black holes in galaxies to the cosmic baryon inventory
T0 review · 3 major / 4 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read Using spectral energy distribution fits to about 800,000 galaxies, this paper argues that stars, dust, neutral gas, and supermassive black holes inside galaxies account for only about 5% of the universe's baryons.
desk verdict A genuinely useful self-consistent baryon census from GAMA+DEVILS, but the dust-mass scale calibration is too fragile to trust the quantitative CDMH/CGMH without caveats. 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 ProSpect, an SED-fitting code that enforces energy balance between stellar attenuation (Charlot & Fall) and dust re-emission (Dale et al.), yielding a dust mass for every galaxy. Because ProSpect's default constant dust-to-hydrogen ratio produces dust masses about 2.5 times higher than the MAGPHYS code, the authors apply a wavelength-dependent dust-to-hydrogen weighting with a PAH+very-small-grain mass fraction of q=0.14, calibrated on 218 bright far-infrared-detected galaxies, to bring ProSpect into agreement with MAGPHYS. A metallicity-dependent dust-to-gas ratio (from Rémy-Ruyer et al.) then converts dust to neutral gas, and double Schechter function fits to the
What would settle it
Measure dust masses for a sample of faint, high-redshift galaxies using far-infrared or submillimeter data that directly constrain the cold-dust component (e.g., ALMA), without imposing the q=0.14 scaling. If these directly constrained dust masses are systematically about 2.5 times higher than the scaled ProSpect values, the global correction is over-applied and the 5% figure is too low; if they agree, the absolute scale is supported.
Extended reading notes
Core claim
The central claim is that galaxies as bounded systems contain only about 5% of the universe's baryons. The authors derive this by fitting FUV-to-FIR spectral energy distributions with the ProSpect code to obtain stellar and dust masses for ~800,000 galaxies in GAMA and DEVILS at 0<z<3, converting dust masses to neutral gas masses using a metallicity-dependent dust-to-gas ratio, and folding in a supermassive black hole mass history inferred from the same data. The resulting cosmic dust mass history peaks at z≈1–2 and falls more gently than the cosmic star formation history, indicating that dust growth and destruction timescales lengthen as star formation declines. The neutral gas density infe
Load-bearing premise
The entire absolute scale of the dust, neutral gas, and the 5% census rests on one global correction: dust masses from ProSpect are divided by about 2.5 so that they match MAGPHYS on a sample of 218 bright, far-infrared-detected galaxies, using an assumed mass fraction of q=0.14 for PAHs and very small grains; if that correction is wrong or does not apply to the faint, high-redshift galaxies that dominate the sample, every derived mass density shifts by the same factor.
Editorial extensions
If this is right
- The cosmic dust mass history peaks at z≈1–2 and declines roughly half as steeply as star formation, implying that dust growth and destruction processes respond to the declining star formation rate over several gigayears.
- The dust-traced neutral gas density is ~0.6 dex lower than 21-cm HI measurements, confirming that HI extends well beyond the optical radii within which the dust emission is measured.
- The baryon fraction locked in galaxies grows from ~1% at z≈3 to ~5% today, meaning the overwhelming majority of baryons have never been incorporated into galaxies at any time in the last 12 billion years.
- With roughly 95% of baryons omitted, the census directly motivates searches for warm-hot ionized gas in haloes and the cosmic web, such as via fast radio burst dispersion measures or X-ray absorption.
- Comparisons with galaxy formation simulations show broad agreement at z≲0.4 but reveal ~0.1 dex higher observed dust density at higher redshift, providing a constraint on feedback and dust destruction implementations.
Reading between the lines
- The q=0.14 dust-mass correction is calibrated on 218 bright, mostly low-redshift galaxies but is applied uniformly to all ~800,000 galaxies; spatially resolved dust maps from ALMA or JWST across a wider redshift range could test whether the correction varies with environment or redshift.
- If the global dust correction is over-applied at high redshift, the neutral gas history and the 5% figure would also shift, since the neutral gas masses are dust masses divided by a metallicity-dependent dust-to-gas ratio.
- The paper's self-consistent design (all four mass densities from the same data and SED fits) means future observational or theoretical comparisons can isolate astrophysical evolution without worrying about cross-survey systematic offsets.
- One could re-scale the neutral gas history using CO-based molecular gas observations without redoing the SED fits, since the underlying dust mass distributions are already provided.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper uses ProSpect SED fits of approximately 800,000 galaxies from GAMA and DEVILS at 0<z<3 to construct stellar, dust, and neutral gas mass distribution functions, integrate them into cosmic mass densities, and combine them with a previously derived SMBH mass density to present a galaxy baryon census. The authors report that the cosmic dust mass history (CDMH) peaks at z≈1–2 and declines more slowly than the cosmic star formation history; that dust-traced neutral gas is on average about 0.6 dex below 21cm HI measurements; and that stars, neutral gas, SMBHs, and dust within optical radii make up approximately 5 per cent of the cosmic baryon budget. The central caveats are that the absolute dust scale is set by calibrating ProSpect dust masses to MAGPHYS using 218 low-redshift FIR-bright galaxies, and that high-redshift DEVILS dust masses rely on a Gaussian prior on the Dale et al. (2014) alpha parameter.
Significance. If the dust-scale uncertainty is addressed, the paper makes a valuable contribution: it provides a homogeneous census of galaxy-bound baryons at 0<z<3, with tabulated Schechter fits, public code and data, and a falsifiable headline statement (the 5 per cent galactic baryon fraction) that is robust because stars dominate the budget and dust plus neutral gas contribute less than about 1 per cent. The comparisons of the CDMH with simulations and of the CGMH with 21cm measurements are also useful. However, the CDMH amplitude and shape, and the reported 0.6 dex CGMH deficit relative to HI surveys, are not as secure as the quoted ~0.05 dex statistical errors imply, because they inherit the global dust-mass calibration and the alpha prior. These issues do not undermine the 5 per cent census, but they do affect two of the paper's other headline results.
major comments (3)
- [§3.2, Eq. (4); Table A5] The dust-mass scale is set by a global q=0.14 variable-DTH correction calibrated on 218 low-redshift, FIR-bright GAMA galaxies and then applied uniformly to all ~800,000 galaxies, including DEVILS galaxies at z>2 where almost none have FIR detections. The paper itself notes that q=0.1 gives ~3.1 times lower dust masses, and §3.3 shows that plausible changes in the opacity coefficient can change dust masses by factors of a few. These are systematic uncertainties, but Table A5 quotes CDMH and CGMH uncertainties of only ~0.05 dex and appears not to include them. Because the CGMH is just the dust mass divided by a metallicity-dependent DTG, the reported 0.6 dex offset from 21cm measurements could absorb a rescaling of this order. I ask for a per-object or at least redshift-dependent correction, and for a propagated systematic error envelope in Table A5 and Fig. 9 covering the plausible range
- [§3.1, Fig. 2; Fig. 9, Table A5] For z>2, Fig. 1 shows that only about 1 per cent of DEVILS sources have MIR/FIR photometry. Their dust masses are therefore controlled by the Gaussian prior on alpha with mean 2 and sigma 1, rather than by direct FIR constraints. Fig. 2 demonstrates that including FIR data shifts the dust-mass posterior by about 0.7 dex because the fit prefers alpha≈3. The high-redshift CDMH points in Fig. 9 and Table A5 (e.g. z=2.311 and 3.090) and the location of the peak at z≈1–2 are thus prior-dominated to a significant degree. Please quantify this sensitivity by rerunning the high-redshift CDMH with an alpha prior centred at 3, or with no prior, and reporting the resulting shape; if this is not possible, the claims about the peak and the high-z decline should be softened.
- [§6, Fig. 12; references to D'Silva et al. (2023, 2025)] The term 'self-consistent' overstates the case for the dust and SMBH components. The dust scale is calibrated to make ProSpect agree with MAGPHYS on 218 galaxies, so the subsequent agreement of the CDMH with MAGPHYS-based literature values is partly by construction. The SMBH mass density is not derived in this paper but is imported from D'Silva et al. (2023, 2025). The 5 per cent baryon census remains a useful aggregation, but the text should describe it as homogeneous in data and fitting infrastructure, while acknowledging that the dust absolute scale is inherited from a MAGPHYS-based calibration and that the SMBH curve is taken from the authors' earlier analysis.
minor comments (4)
- [Abstract vs. §6] The abstract states the neutral gas deficit is on average ≈0.7 dex lower than 21cm measurements, while the full text and Section 6 say ≈0.6 dex. Please make these consistent.
- [Caption of Fig. 12] The caption says 'baron density' in the bottom panel; this should be 'baryon density'.
- [§7, item (iii)] The conclusions say ProSpect gives '≈3 times higher dust masses' than MAGPHYS, while §3.2 says the correction is '≈2.5 times'. Please reconcile these numbers or clarify that 3 refers to the uncorrected code comparison and 2.5 to the adopted correction.
- [§3.3, opacity discussion] The opacity coefficient discussion would benefit from a quantitative statement of how much of the factor-of-three offset remains after applying the variable DTH correction, since the paper currently leaves the reader with two partially degenerate explanations for the ProSpect-MAGPHYS offset.
Circularity Check
No significant circularity: the dust-mass rescaling is an openly described calibration, the SMBH density is an independent prior result, and the 5% census is an aggregation of measured components.
full rationale
The central derivation chain is linear and not self-referential: stellar masses come from ProSpect SED fits, which give the SMFs and CSMH; dust masses come from the same fits, giving the DMFs and CDMH; neutral gas masses are dust masses divided by a metallicity-dependent DTG, giving the CGMH; and the SMBH density is imported from prior published work. The only step that could look circular is the Section 3.2 dust-mass rescaling, where q = 0.14 is chosen so that ProSpect dust masses agree with MAGPHYS on 218 FIR-bright galaxies, and the resulting factor of about 2.5 is then applied globally. But the paper presents this explicitly as a calibration, not as a prediction: it states that the scaling is sensitive to q and that q = 0.1 would give about 3.1 times lower dust masses. The redshift-dependent shape of the CDMH and CGMH is produced by the SED fits themselves, not by this multiplicative constant, so the comparison to MAGPHYS-based literature is a consistency check rather than a derivation forced by construction. The SMBH density is taken from D'Silva et al. 2023, 2025, which is a self-citation, but those are previously published, externally reviewed measurements using the same data and methods; they are not derived from equations in this paper, and the SMBH term is a minor contributor to the final 5% figure. No self-definition, no fitted parameter renamed as a prediction, no ansatz smuggled in via citation, no imported uniqueness theorem, and no renaming of a known result was found. The paper's own caveat about the q-sensitivity of the absolute dust/gas scale is a robustness limitation, not a circularity.
Assumptions & free parameters
free parameters (4)
- PAH+VSG mass fraction q =
0.14
- Global dust-mass correction factor =
≈2.5
- Gaussian prior on Dale et al. α parameters =
mean=2, σ=1
- SMBH radiative efficiency =
0.1
assumptions (6)
- domain assumption ProSpect energy-balance SED fitting maps absorbed starlight to reradiated dust emission, with dust mass from template mass-to-light ratios.
- domain assumption Dale et al. (2014) dust emission templates and Charlot & Fall (2000) attenuation describe the dust SED; the α prior constrains the radiation field.
- domain assumption A solar-metallicity DTH=0.0073, modified by a q-weighted wavelength-dependent DTH, sets the dust mass scale.
- domain assumption The Rémy-Ruyer et al. (2014) broken power-law dust-to-gas ratio applies as a function of ProSpect Zfinal to convert dust to neutral gas for all galaxies and redshifts.
- domain assumption The SMBH mass history from D'Silva et al. (2023, 2025), based on a 10% radiative efficiency and integration of the AGN luminosity history, is adopted unchanged.
- ad hoc to paper Quadratic polynomials in redshift for the double-Schechter parameters provide the smooth underlying evolution whose residual versus direct fits gives the LSS correction.
Cite this review
Pith. "Pith review of The contribution of stars, dust, neutral gas and supermassive black holes in galaxies to the cosmic baryon inventory." pith.science (2026). https://pith.science/paper/JMLBACYE
@misc{pith2026260108112,
author = {Pith},
title = {Pith review of: The contribution of stars, dust, neutral gas and supermassive black holes in galaxies to the cosmic baryon inventory},
year = {2026},
howpublished = {\url{https://pith.science/paper/JMLBACYE}},
note = {Machine review of arXiv:2601.08112}
}
abstract
We compute the cosmic stellar, dust and neutral gas mass history at $0<z\lesssim3$ using ProSpect spectral energy distribution modelling of $\approx 800 \, 000$ galaxies in the Galaxy and Mass Assembly (GAMA) survey and the Deep Extragalactic VIsible Legacy Survey (DEVILS). The cosmic dust mass history broadly follows the shape of the cosmic star formation history; though, the decline is slower, suggestive of a slowing rate of dust growth and destruction as the star formation declines past its peak at $z\approx 2$. Neutral gas masses were estimated by scaling the dust masses by the metallicity-dependent dust-to-gas ratio. The neutral gas mass density as traced by the dust is an average of $\approx 0.7$ dex lower than that measured from $21$cm experiments, most likely due to differences in the spatial scales inhabited by dust and HI. Folding in measurements of the supermassive black hole mass density obtained previously with similar data and methods, we present a self-consistent census of the baryons confined to galaxies. Stars, neutral gas, SMBHs and dust contained within the optical radii of galaxies account for $\approx 5$ per cent of the baryons. Most of the remaining $\approx 95$ per cent of baryons must be ionised and dispersed throughout the interstellar, circumgalactic and intergalactic media within, around and between galaxies.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 1 Pith paper
-
The evolution of galaxy dust scaling relations in the COLIBRE simulations
The COLIBRE simulations reproduce the broad evolution of galaxy dust scaling relations from z=15 to 0, but overproduce the z<1 cosmic dust mass density and miss the most extreme sub-millimeter galaxies.
Reference graph
Works this paper leans on
-
[1]
Aghanim N., et al., 2020, Astronomy & Astrophysics, 641, A6 Aniano G., et al., 2020, The Astrophysical Journal, 889, 150 BaldryI.K.,etal.,2012,MonthlyNoticesoftheRoyalAstronomicalSociety, 421, 621 Beeston R. A., Gomez H. L., Dunne L., Maddox S., Eales S. A., Smith M. W. L., 2024, Monthly Notices of the Royal Astronomical Society, 535, 3162 BellagambaF.,Ma...
arXiv 2020
-
[12]
MNRAS000, 1–14 (2026)
Each of these columns and the uncertainties factor in the LSScorrectionsandtheirassociateduncertainties,whicharethemselvespresentedintherightmostcolumn.TheCDMHandCGMHincludecorrectionsforthe metallicity-dependent DTH and have been lowered by≈2.5as per the discussion in Section 3.2. MNRAS000, 1–14 (2026)
2026
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.