Pith. sign in

REVIEW 2 major objections 4 references

Spatially resolved data show multiple quenching pathways in ultra-massive galaxies at z approximately 3.5

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

2026-06-29 15:21 UTC pith:URODCNCM

load-bearing objection First spatially resolved [α/Fe] at z~3.5 for three galaxies is the actual advance, but the quenching mechanism mapping is preliminary. the 2 major comments →

arxiv 2605.27555 v1 pith:URODCNCM submitted 2026-05-26 astro-ph.GA

MAGAZ3NE: Spatially Resolved Ages and Chemical Abundances of Ultra-Massive Quiescent Galaxies at z sim 3.5 using JWST/NIRSpec IFU

Adit H. Edward , Jacqueline Antwi-Danso , Adam Muzzin , Ben Forrest , Ian McConachie , Aliza Beverage , Wenjun Chang , M. C. Cooper
show 9 more authors
Percy Gomez Massissilia L. Hamadouche Aur\'elien Henry Han Lei Danilo Marchesini Allison Noble Stephanie M. Urbano Stawinski Gillian Wilson M. E. Wisz
This is my paper
classification astro-ph.GA
keywords quiescent galaxieshigh-redshift galaxiesstellar ageschemical abundancesquenchingJWST observationsalpha-enhancement
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper measures stellar ages, iron abundances, and alpha-to-iron ratios at different locations within three ultra-massive quiescent galaxies at redshift around 3.5. The measurements come from integral field spectroscopy with the James Webb Space Telescope. All three galaxies have young centers with high alpha enhancement, pointing to quick star formation just before they stopped forming stars. Two of them show age increasing and alpha decreasing with radius, which matches expectations for quenching triggered by mergers, while the third has uniform ages suggesting a different process. These findings indicate that different ways of shutting down star formation were already in use early in the universe's history.

Core claim

The authors find that the central regions of the three galaxies are uniformly young at 0.6 to 0.7 billion years with elevated alpha to iron ratios of 0.2 to 0.5. Two galaxies show positive age gradients and negative alpha gradients outward, consistent with rapid merger-driven quenching, while the third has a flat age profile. Metallicities differ by 0.2 to 0.4 dex depending on whether alpha-enhanced models are used, underscoring the need for such models when studying these early systems.

What carries the argument

Spatially resolved stellar population parameters including age, [Fe/H], and [α/Fe] derived from JWST/NIRSpec IFU spectra, which preserve chemical signatures before mergers can erase them.

Load-bearing premise

The gradients in age and alpha-enhancement are caused by the quenching mechanisms themselves and not by choices in spectral modeling, dust, or post-quenching changes.

What would settle it

Repeating the analysis with solar-scaled templates or different dust corrections that removes the positive age gradients and negative alpha gradients would falsify the link to merger-driven quenching.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Quenching pathways were diverse already at z~3.5
  • Rapid merger-driven quenching operated in some massive galaxies within the first two billion years
  • Alpha-enhancement must be treated explicitly to interpret star-formation histories correctly
  • The data provide the first spatially-resolved alpha-enhancement constraints at this epoch

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If confirmed in larger samples, this would mean merger activity played a role in galaxy quenching very early on
  • Incorrect model choices could systematically bias age and metallicity estimates for high-redshift quiescent galaxies
  • Similar observations at even higher redshifts could trace the onset of these processes

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 0 minor

Summary. The paper presents the first spatially-resolved stellar population measurements (age, [Fe/H], [α/Fe]) for three ultra-massive (log M*/M⊙ > 11), compact (Re ≲ 2 kpc) quiescent galaxies at z ∼ 3.5 using JWST/NIRSpec IFU spectroscopy. All three show young (≈0.6–0.7 Gyr), α-enhanced (≈0.2–0.5) cores. Two exhibit positive age gradients and negative [α/Fe] gradients, interpreted as rapid merger-driven quenching; the third has a flat age profile, interpreted as uniform quenching. [Fe/H] gradients are noted as consistent with these trends, but metallicities differ by 0.2–0.4 dex between α-enhanced and solar-scaled models. The work concludes that quenching pathways are already diverse by z ∼ 3.5 and that explicit α-enhancement treatment is required for early quenched systems.

Significance. If the gradient-to-mechanism mapping and model comparisons are robust, the results would provide the first direct spatially-resolved α-enhancement constraints at this epoch, supporting early diversity in quenching and the necessity of α-enhanced templates. The N=3 sample and lack of quantitative robustness tests limit the strength of the diversity claim, but the observational advance on α at z∼3.5 would still be notable if the central interpretations hold after addressing systematics.

major comments (2)
  1. [Abstract] Abstract: The central claim that positive age + negative [α/Fe] gradients indicate rapid merger-driven quenching (while flat profiles indicate uniform quenching) is load-bearing for the diversity conclusion, yet the text provides no quantitative mapping, comparison to hydrodynamical simulations, or tests against alternatives (spectral-fitting degeneracies, dust, radial IMF variations, or post-quenching migration). The 0.2–0.4 dex [Fe/H] offset between model sets is noted but not propagated into gradient uncertainties or sign stability.
  2. [Abstract] Abstract: No error budgets, data-reduction details, or robustness tests (e.g., alternative fitting setups, different SSP libraries) are described for the reported gradients or the N=3 diversity statement; with such a small sample these omissions directly affect whether the quenching-pathway conclusions can be drawn from the observations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful and constructive report. The comments correctly identify areas where the manuscript's interpretive claims and robustness documentation can be strengthened. We address each point below and indicate the revisions that will be made.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that positive age + negative [α/Fe] gradients indicate rapid merger-driven quenching (while flat profiles indicate uniform quenching) is load-bearing for the diversity conclusion, yet the text provides no quantitative mapping, comparison to hydrodynamical simulations, or tests against alternatives (spectral-fitting degeneracies, dust, radial IMF variations, or post-quenching migration). The 0.2–0.4 dex [Fe/H] offset between model sets is noted but not propagated into gradient uncertainties or sign stability.

    Authors: We agree that the gradient-to-mechanism mapping is currently qualitative and would benefit from more explicit support. In the revised manuscript we will (i) add a dedicated paragraph in the discussion section that compares the observed gradient signs and amplitudes to published expectations from hydrodynamical simulations (e.g., those that track α-enhancement), (ii) explicitly list the main systematics (dust, IMF gradients, post-quenching migration) and state why they are unlikely to reverse the reported gradient signs given the data quality, and (iii) propagate the 0.2–0.4 dex [Fe/H] model offset into the quoted gradient uncertainties and re-evaluate sign stability. A full quantitative mapping or exhaustive alternative-model grid is beyond the scope of this first IFU study but will be flagged as future work. revision: partial

  2. Referee: [Abstract] Abstract: No error budgets, data-reduction details, or robustness tests (e.g., alternative fitting setups, different SSP libraries) are described for the reported gradients or the N=3 diversity statement; with such a small sample these omissions directly affect whether the quenching-pathway conclusions can be drawn from the observations.

    Authors: The full manuscript already contains data-reduction procedures (Section 2) and per-spaxel error budgets derived from the spectral fits. However, we acknowledge that additional robustness checks were not presented. We will add a new subsection (or appendix) that reports results from (a) two independent SSP libraries, (b) alternative regularization choices in the fitting code, and (c) a Monte-Carlo resampling of the IFU cubes. The diversity statement will be rephrased to emphasize that it is based on the variety seen in this pilot sample of three galaxies and that larger samples are required to quantify the relative frequency of quenching pathways. revision: yes

Circularity Check

0 steps flagged

No circularity: purely observational fitting with no self-referential derivations

full rationale

The paper reports direct measurements of age, [Fe/H], and [α/Fe] gradients from JWST/NIRSpec IFU spectra of three galaxies using standard stellar-population synthesis fitting. No equations, predictions, or uniqueness theorems are presented that reduce fitted quantities back to input parameters by construction. The interpretive claims about quenching mechanisms (merger-driven vs. uniform) are post-hoc attributions based on observed gradient signs, not derived from any self-citation chain or ansatz that loops to the data. This is a standard observational analysis with no load-bearing self-referential steps.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Measurements depend on the accuracy of stellar population synthesis models at high redshift and on the assumption that IFU spectra can be cleanly decomposed into age, metallicity, and alpha components without major systematics.

free parameters (1)
  • stellar population model parameters (age, [Fe/H], [α/Fe])
    These quantities are fitted to each spatial bin; the reported values and gradients are outputs of the fit.
axioms (1)
  • domain assumption Stellar population synthesis models with alpha-enhanced isochrones accurately reproduce the spectra of z~3.5 quiescent galaxies
    Invoked when comparing alpha-enhanced versus solar-scaled template results and when deriving ages and abundances.

pith-pipeline@v0.9.1-grok · 5939 in / 1401 out tokens · 50970 ms · 2026-06-29T15:21:06.484480+00:00 · methodology

0 comments
read the original abstract

We present spatially-resolved measurements of stellar age, [Fe/H], and [$\alpha$/Fe] in three ultra-massive ($\rm{log(M_{\ast}/M_{\odot})>11}$), compact ($\rm{R_e} \lesssim 2$ kpc) quiescent galaxies at $z\sim3.5$ using JWST/NIRSpec IFU spectroscopy. These observations provide the first spatially-resolved constraints on $\alpha$-enhancement at this epoch, enabling a direct test of quenching mechanisms before late-time assembly processes such as mergers can erase chemical signatures. The central regions of all three galaxies show both uniformly young ages ($\approx0.6-0.7$ Gyr) and elevated [$\alpha$/Fe] ($\approx0.2-0.5$), indicating rapid, enhanced star formation shortly before recent quenching. Beyond the cores, two galaxies display positive age gradients and negative [$\alpha$/Fe] gradients, consistent with rapid merger-driven quenching, while the third shows a flat age profile indicative of uniform quenching. The [Fe/H] gradients are also consistent with these trends, though we note that the metallicities reported by codes using $\alpha$-enhanced models differ significantly ($\approx0.2-0.4$ dex) from those reported using solar-scaled templates. These data demonstrate that quenching pathways are diverse by $z\sim3.5$, with rapid, merger-driven quenching already operating in a subset of massive quiescent galaxies in the first two billion years of cosmic time. Furthermore, these results establish that explicit treatment of $\alpha$-enhancement is essential for interpreting the star-formation histories of the earliest quenched systems.

Figures

Figures reproduced from arXiv: 2605.27555 by Adam Muzzin, Adit H. Edward, Aliza Beverage, Allison Noble, Aur\'elien Henry, Ben Forrest, Danilo Marchesini, Gillian Wilson, Han Lei, Ian McConachie, Jacqueline Antwi-Danso, Massissilia L. Hamadouche, M. C. Cooper, M. E. Wisz, Percy Gomez, Stephanie M. Urbano Stawinski, Wenjun Chang.

Figure 1
Figure 1. Figure 1: Left: False-color images of the MAGAZ3NE massive quiescent galaxies. We collapse the IFU data cubes along the wavelength axis and take the mean value for each pixel in the wavelength ranges of (6250˚A, 7500˚A), (4000˚A, 5600˚A) and (3735˚A, 4000˚A) for the R, G and B channels, respectively. We outline the effective radius (Re) of each galaxy in red and list the circularized Re in kpc from B. Forrest et al.… view at source ↗
Figure 2
Figure 2. Figure 2: Size vs. mass comparison of massive quiescent galaxies both from this work and in literature, colored by redshift. Galaxies used in C. M. Cheng et al. (2025) from JWST-SUSPENSE (M. Slob et al. 2024) are shown as dia￾monds, and the galaxy studied in F. D’Eugenio et al. (2024) is shown as a square, with reported error bars. Galaxies from L. Kawinwanichakij et al. (2026) are shown as hexagons, us￾ing the infe… view at source ↗
Figure 3
Figure 3. Figure 3: Left: Bagpipes star-formation histories. Solid lines show the non-parametric Gaussian Process SFH (K. G. Iyer et al. 2019), and dashed lines show SFHs adopting the continuity model (J. Leja et al. 2019). Right: Mass assembly histories based on the inferred SFHs. We compare these with the probability density function of the most massive galaxy expected in the XMM-VIDEO survey volume (assuming an area of 4.6… view at source ↗
Figure 4
Figure 4. Figure 4: Age (Top), [Fe/H] (Middle), and [α/Fe] (Bottom) gradients for all three galaxies in this study. We outline measure￾ments within R < 1Re and plot them at 1 2Re, slightly shifted for visual purposes. The remaining points are measured in elliptical annuli centered around the galaxy, in bins with a semi-major axis width of 1 pixel. We illustrate gradients using a least-squares fit to the radial measurements. A… view at source ↗
Figure 5
Figure 5. Figure 5: Same as [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Left: Radial SFHs of each massive quiescent galaxy, modelled by Bagpipes assuming a Gaussian Process SFH (K. G. Iyer et al. 2019) (solid lines) and continuity model (J. Leja et al. 2019) (dashed lines). Right: Mass assembly history as a function of cosmic time. Both are shown in units of surface density. Shaded regions correspond to 1σ confidence intervals. 5. DISCUSSION 5.1. Formation Processes for MQGs a… view at source ↗
Figure 7
Figure 7. Figure 7: Schematic illustrating different quenching pathways possible for massive quiescent galaxies, along with the expected gradients and timescales involved. Top row, (a): A rapid, merger driven quenching event. The expected gradients match those observed for XMM-VID1-2075 and XMM-VID3-1120. Upper middle row, (b): An inside-out quenching event, where star formation in the central regions terminates first. Lower … view at source ↗
Figure 8
Figure 8. Figure 8: Extracted spectra for XMM-VID1-2075, within R < Re (top, black) and elliptical radial bins, normalized and scaled arbitrarily for visual purposes. We show the results of this testing in [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Same as [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Same as [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Results from fitting a model galaxy made in Bagpipes with alfα, using the MILES+IRTF models from C. Conroy et al. (2018) (orange) and sMILES (A. T. Knowles et al. 2023) (blue) for age (top left), [Fe/H] (top right), metallicity (bottom left), and [α/Fe] (bottom right). Posterior fitting distributions are presented as probability density functions. The model galaxy is made with an age of 1.1 Gyr, Z∗ = Z⊙, … view at source ↗
Figure 12
Figure 12. Figure 12: Top: Red: Bagpipes model spectrum of a 1.1 Gyr galaxy with Z∗ = Z⊙, [α/Fe] = 0, and σ∗ = 200 km/s. MILES+IRTF (orange) and sMILES (blue) fits from alfα. Bottom: The residuals between the MILES+IRTF and sMILES fits and the model galaxy. Significant differences between the models are seen, particularly for Balmer absorption line depths. The MILES+IRTF spectrum has deeper Balmer lines than either sMILES or B… view at source ↗
Figure 13
Figure 13. Figure 13: A comparison of different stellar population parameters inferred from alfα (y-axis) and Bagpipes (x-axis) for our integrated spectra (R < Re). The one-to-one line is indicated as a black dashed line. Left: SSP-equivalent ages from alfα compared to the mass-weighted ages from Bagpipes. Center-left: [Z/H] inferred from alfα versus log(Z∗/Z⊙) from Bagpipes. Center-right: [Fe/H] versus log(Z∗/Z⊙) for the thre… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

4 extracted references · 3 canonical work pages · 2 internal anchors

  1. [1]

    E., Brammer, G., et al

    Akhshik, M., Whitaker, K. E., Brammer, G., et al. 2020, ApJ, 900, 184, doi: 10.3847/1538-4357/abac62 Akhshik, M., Whitaker, K. E., Leja, J., et al. 2023, ApJ, 943, 179, doi: 10.3847/1538-4357/aca677 Antwi-Danso, J., Papovich, C., Esdaile, J., et al. 2025, ApJ, 978, 90, doi: 10.3847/1538-4357/ad8b30 Araya-Araya, P., Cochrane, R. K., Hayward, C. C., et al. ...

    work page doi:10.3847/1538-4357/abac62 2020

  2. [2]

    (2019) Non-parametric SFH SFR (M⊙ yr−1) (1×10 −8, 2×10

    Iyer et al. (2019) Non-parametric SFH SFR (M⊙ yr−1) (1×10 −8, 2×10

    2019

  3. [3]

    M., Lim , P

    log 10 αλ<5000˚A (−2,4) Gaussian (µ=−1.5, σ= 0.5) αλ>5000˚A (−2,2) Gaussian (µ= 0.5, σ= 0.5) aOnly for XMM-VID3-2457. T able 3.Priors used inBagpipesspectral fitting. Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 336, Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, ed...

    work page internal anchor Pith review doi:10.3847/1538-4357/ac7c74 2005

  4. [4]

    Fast and scalable Gaussian process modeling with applications to astronomical time series

    (blue) for age (top left), [Fe/H] (top right), metallicity (bottom left), and [α/Fe] (bottom right). Posterior fitting distributions are presented as probability density functions. The model galaxy is made with an age of 1.1 Gyr, Z ∗ = Z⊙, [α/Fe] = 0, andσ ∗ = 200 km/s. These values are indicated by a black dashed line. Foreman-Mackey, D., Agol, E., Ambik...

    work page internal anchor Pith review doi:10.3847/1538-3881/aa9332 2017