Stellar Surface Density Modulates MgII Cool-gas Outflow Absorption in DESI Star-forming Galaxies
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 kernel 2026-07-01 06:30 UTCgrok-4.3pith:DR6JZMAArecord.jsonopen to challenge →
The pith
Stellar surface density increases MgII cool-gas outflow absorption in star-forming galaxies beyond effects of total mass or star-formation rate.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In AGN-clean samples of massive star-forming galaxies at 0.35 < z < 1.0 drawn from DESI Data Release 1, the MgII outflow equivalent width rises monotonically with stellar surface density Σ★ = M★ / (2π R_e²) in every redshift bin. The increase from lowest to highest Σ★ tertile is 0.37–0.61 Å when samples are matched on stellar mass alone and persists in a stricter sample also matched on a Balmer-line SFR proxy. The absolute outflow velocity changes only weakly with Σ★, while the structural trend is interpreted as evidence that galaxy compactness alters the absorbing velocity distribution or the effective covering fraction of cool outflowing gas.
What carries the argument
Stellar surface density Σ★ = M★ / (2π R_e²), used as an independent empirical coordinate that modulates the strength of down-the-barrel MgII cool-gas absorption after mass and SFR matching.
Load-bearing premise
That matching samples on stellar mass and Balmer-line SFR proxy plus AGN cleaning isolates stellar surface density as the independent variable without residual biases from other galaxy properties or measurement systematics.
What would settle it
A new stacked or individual-galaxy sample, matched identically on stellar mass and Balmer SFR proxy, in which MgII equivalent width shows no systematic rise across stellar surface density tertiles.
Figures
read the original abstract
Galaxy outflows are usually ordered by stellar mass and star-formation rate (SFR), but the same feedback budget may couple differently to gas in diffuse and compact galaxies. We use Dark Energy Spectroscopic Instrument (DESI) Data Release 1 stacked spectra of massive star-forming galaxies at $0.35<z<1.0$ to test whether stellar surface density, $\Sigma_\star=M_\star/(2\pi R_e^2)$, is an independent empirical coordinate of down-the-barrel singly ionized magnesium (MgII) cool-gas absorption. In AGN-clean samples matched in stellar mass, and in a stricter sample matched in both stellar mass and a Balmer-line SFR proxy, the MgII outflow equivalent width (EW) rises monotonically with $\Sigma_\star$ in every redshift bin. From the lowest to highest $\Sigma_\star$ tertile, EW increases by 0.37-0.61 Angstrom, while the absolute outflow velocity changes only weakly. DESI therefore shows that cool-gas outflow strength in massive star-forming galaxies is not set only by how much stellar mass or star formation a galaxy has, but also by how tightly the galaxy is built. The structural dependence points to changes in the absorbing velocity distribution and/or the effective covering fraction of cool outflowing gas.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes DESI DR1 stacked spectra of massive star-forming galaxies at 0.35<z<1.0 to show that MgII cool-gas outflow equivalent width increases monotonically with stellar surface density Σ*=M*/(2πRe²) after AGN cleaning and matching on stellar mass (and more strictly on both M* and a Balmer-line SFR proxy). EW rises 0.37-0.61 Å across Σ* tertiles in every redshift bin while outflow velocity changes only weakly, implying that galaxy compactness modulates cool-gas outflow properties independently of total stellar mass or SFR.
Significance. If the matching isolates Σ* without residual biases, the result supplies a clear empirical demonstration that structural compactness is an additional coordinate for cool-gas outflow strength in massive star-forming galaxies, based on a large spectroscopic dataset with redshift binning. This strengthens the case for incorporating surface density into feedback models beyond the usual M*–SFR scaling.
major comments (2)
- [Sample construction and matching] The sample-matching procedure (described in the methods section on AGN-clean samples) must demonstrate that residuals in dust attenuation, Sérsic index, or gas-phase metallicity are balanced across Σ* tertiles; otherwise the monotonic EW trend could arise from these correlated properties rather than Σ* itself, undermining the claim that Σ* is an independent driver.
- [Results on velocity and EW trends] Given the paper's note that absolute outflow velocity changes only weakly, the EW increase is plausibly driven by covering fraction or optical-depth variations rather than changes in outflow strength or velocity distribution; this distinction should be quantified (e.g., via line-profile modeling) to support the interpretation of modulated outflow strength.
minor comments (2)
- [Methods] Clarify the exact definition and measurement of the Balmer-line SFR proxy used for the stricter matching, including any aperture or dust corrections applied.
- [Results] Add error bars or bootstrap uncertainties to the reported EW differences (0.37-0.61 Å) and state the number of galaxies per tertile and redshift bin.
Simulated Author's Rebuttal
We thank the referee for their constructive comments, which help clarify the robustness of our results. We address each major comment below and indicate where revisions will be made.
read point-by-point responses
-
Referee: [Sample construction and matching] The sample-matching procedure (described in the methods section on AGN-clean samples) must demonstrate that residuals in dust attenuation, Sérsic index, or gas-phase metallicity are balanced across Σ* tertiles; otherwise the monotonic EW trend could arise from these correlated properties rather than Σ* itself, undermining the claim that Σ* is an independent driver.
Authors: We agree that explicit checks on these residuals are needed to support the independence claim. In the revised manuscript we will add supplementary material (tables and/or violin plots) showing the distributions and medians of dust attenuation (A_V from SED fits or Balmer decrement), Sérsic index (from DESI imaging or external catalogs), and gas-phase metallicity (where emission-line measurements are available) across the Σ* tertiles in both the M*-matched and M*+SFR-proxy-matched samples. Preliminary internal checks indicate these quantities are well balanced after matching, but we will quantify any small residuals and discuss their possible influence on the EW trend. revision: yes
-
Referee: [Results on velocity and EW trends] Given the paper's note that absolute outflow velocity changes only weakly, the EW increase is plausibly driven by covering fraction or optical-depth variations rather than changes in outflow strength or velocity distribution; this distinction should be quantified (e.g., via line-profile modeling) to support the interpretation of modulated outflow strength.
Authors: We already note in the manuscript that the weak velocity change implies the EW trend is driven by covering fraction or optical-depth variations rather than bulk velocity shifts. Full line-profile modeling of stacked spectra would require detailed assumptions about geometry, ionization, and multiple components that are not feasible within the scope of this empirical study. We will expand the discussion to more explicitly frame the result in terms of increased effective covering fraction or column density in compact galaxies, while retaining the current interpretation that Σ* modulates cool-gas outflow properties independently of M* and SFR. revision: partial
Circularity Check
No circularity: purely observational empirical result
full rationale
The paper reports an empirical trend from stacked DESI spectra: after matching AGN-clean samples on stellar mass (and stricter matching on both M* and Balmer-line SFR proxy), MgII outflow EW increases monotonically with Σ* across tertiles in each redshift bin. No derivations, model equations, parameter fits, or predictions are present that could reduce to inputs by construction. No self-citations, ansatzes, or uniqueness theorems are invoked. The central claim is a direct observational measurement of a correlation after sample matching, with no load-bearing steps that equate output to input by definition. This is the expected non-finding for a data-driven galaxy survey analysis.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
- [1]
- [2]
-
[3]
Bordoloi, R., Lilly, S. J., Hardmeier, E., et al. 2014, ApJ, 794, 130
work page 2014
- [4]
-
[5]
A., Leitherer, C., Chen, Y ., Wofford, A., & Lundgren, B
Chisholm, J., Tremonti, C. A., Leitherer, C., Chen, Y ., Wofford, A., & Lundgren, B. 2015, ApJ, 811, 149
work page 2015
-
[6]
A., Leitherer, C., Chen, Y ., Wofford, A
Chisholm, J., Tremonti, C. A., Leitherer, C., Chen, Y ., Wofford, A. 2016, MNRAS, 457, 3133
work page 2016
- [7]
- [8]
-
[9]
M., Moustakas, J., Tremonti, C
Diamond-Stanic, A. M., Moustakas, J., Tremonti, C. A., et al. 2012, ApJL, 755, L26 DESI Collaboration et al. 2022, AJ, 164, 207 DESI Collaboration et al. 2026, AJ, 171, 285
work page 2012
- [10]
-
[11]
Erb, D. K., Quider, A. M., Henry, A. L., & Martin, C. L. 2012, ApJ, 759, 26
work page 2012
-
[12]
Franx, M., van Dokkum, P. G., F¨orster Schreiber, N. M., Wuyts, S., Labb´e, I., Toft, S. 2008, ApJ, 688, 770
work page 2008
-
[13]
Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Armus, L. 2000, ApJS, 129, 493
work page 2000
- [14]
-
[15]
Kauffmann, G., Heckman, T. M., Tremonti, C., et al. 2003, MNRAS, 346, 1055
work page 2003
- [16]
- [17]
-
[18]
J., Groves, B., Kauffmann, G., & Heckman, T
Kewley, L. J., Groves, B., Kauffmann, G., & Heckman, T. 2006, MNRAS, 372, 961
work page 2006
-
[19]
Kornei, K. A., Shapley, A. E., Martin, C. L., Coil, A. L., Lotz, J. M., Schiminovich, D., Bundy, K., Noeske, K. G. 2012, ApJ, 758, 135
work page 2012
-
[20]
Martin, C. L. 2005, ApJ, 621, 227
work page 2005
-
[21]
Martin, C. L., Shapley, A. E., Coil, A. L., et al. 2012, ApJ, 760, 127
work page 2012
-
[22]
Martin, C. L., Shapley, A. E., Coil, A. L., Kornei, K. A., Murray, N., & Pancoast, A. 2013, ApJ, 770, 41
work page 2013
-
[23]
L., Kereˇs, D., Faucher-Gigu`ere, C.-A., et al
Muratov, A. L., Kereˇs, D., Faucher-Gigu`ere, C.-A., et al. 2015, MNRAS, 454, 2691
work page 2015
-
[24]
Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569
work page 2005
-
[25]
Naab, T., & Ostriker, J. P. 2017, ARA&A, 55, 59
work page 2017
-
[26]
F., Genzel, R., F¨orster-Schreiber, N
Newman, S. F., Genzel, R., F¨orster-Schreiber, N. M., et al. 2012, ApJ, 761, 43
work page 2012
- [27]
-
[28]
J., Zhang, H.-X., Cao, T., Puzia, T
Rong, Y ., Zhu, K., Johnston, E. J., Zhang, H.-X., Cao, T., Puzia, T. H., Galaz, G. 2020, ApJL, 899, L12
work page 2020
-
[29]
Rubin, K. H. R., Weiner, B. J., Koo, D. C., et al. 2010, ApJ, 719, 1503
work page 2010
-
[30]
Rubin, K. H. R., Prochaska, J. X., M´enard, B., et al. 2011, ApJ, 728, 55
work page 2011
- [31]
-
[32]
S., Veilleux, S., & Sanders, D
Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005, ApJS, 160, 115
work page 2005
-
[33]
Sell, P. H., Tremonti, C. A., Hickox, R. C., et al. 2014, MNRAS, 441, 3417
work page 2014
-
[34]
Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65
work page 2003
- [35]
- [36]
-
[37]
A., Moustakas, J., & Diamond-Stanic, A
Tremonti, C. A., Moustakas, J., & Diamond-Stanic, A. M. 2007, ApJL, 663, L77
work page 2007
-
[38]
Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARA&A, 43, 769
work page 2005
-
[39]
Veilleux, S., Maiolino, R., Bolatto, A. D., & Aalto, S. 2020, A&A Rev., 28, 2
work page 2020
- [40]
-
[41]
Wuyts, S., F¨orster Schreiber, N. M., van der Wel, A., et al. 2011, ApJ, 742, 96
work page 2011
-
[42]
Stellar feedback drives the baryon deficiency in low-mass galaxies
Yu, H., et al. 2025, arXiv:2512.05584
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [43]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.