REVIEW 3 major objections 6 minor 1 cited by
Spatially resolved star-formation histories of local post-starburst galaxies: Starburst and quenching spatial patterns consistent with recent mergers
T0 review · 3 major / 6 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read Three post-starburst galaxies each show an outer starburst followed, one billion years later, by a stronger central starburst — a sequence the paper attributes to a recent gas-rich merger.
desk verdict Resolved two-phase outside-in starburst sequence in three PSBs is a plausible new result, but the one-burst-per-bin SFH — whose degeneracy the paper itself concedes — makes the 'slower outer quenching' claim conditional. 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 central machinery is a hierarchical Bayesian spectral-fitting model applied to spatially binned integral-field spectra. Each bin is fitted with a two-component star-formation history — an old exponential component plus a recent double power-law starburst — while five properties (formed stellar mass, burst age, pre-burst metallicity, dust attenuation, and velocity dispersion) are assumed to follow smooth radial profiles, with burst age described by a logistic function of radius. Joint fitting lets the data reveal a radial gradient in burst timing while borrowing strength across bins to constrain low-surface-brightness outer regions.
What would settle it
Fit the outer-region spectra of one of the three galaxies with a star-formation history model that permits two distinct recent bursts (or with a non-parametric SFH) at sufficiently high signal-to-noise; if the outskirts require two separate burst peaks rather than a single slow decline, the claimed outer-then-central sequence is an artefact of the one-burst-per-bin assumption.
Extended reading notes
Core claim
The paper claims that all three galaxies first experienced an outer, weaker and slower-quenching starburst, followed by a central, stronger and faster-quenching starburst that peaked roughly 1 Gyr after the first. This spatial and temporal sequence matches binary merger simulations in which the first pericentre passage triggers star formation in the outer regions and the later coalescence triggers a stronger centralised starburst. The authors also find that the central starburst produced a significantly larger rise in stellar metallicity than the outer one, and that the rapid quenching is consistent with gas consumption and morphological stabilisation by a growing spheroid, without requiring
Load-bearing premise
The fitting model allows each spatial bin to have only one recent starburst, and the recovered radial gradient in burst age is then interpreted as two galaxy-wide bursts in sequence — an assumption the paper itself notes could be violated if outer regions actually had two bursts that blend into one slow decline.
Editorial extensions
If this is right
- Spatially integrated spectra of a post-starburst galaxy can completely miss an earlier, outer starburst, so the central burst alone may misrepresent the full recent star-formation history.
- The roughly 1 Gyr gap between outer and central burst peaks is a potential clock for the interval between first pericentre passage and final coalescence in a wet merger.
- Local post-starburst galaxies with extended PSB regions can be merger remnants even when they show no obvious morphological merger signatures.
- Rapid quenching in such galaxies does not require AGN feedback; gas consumption plus spheroid-driven morphological stabilisation can accomplish it.
- The common inside-out/outside-in quenching labels are too crude to describe the complex spatial and temporal patterns seen in resolved star-formation histories.
Reading between the lines
- If this two-burst pattern is common among post-starburst galaxies, the outer-burst age could be used as an observational estimate of the time since first pericentre, turning these galaxies into merger-timeline probes.
- The non-axisymmetric dust structures and burst-mass-fraction arcs would be promising targets for high-resolution cold-gas follow-up; molecular-gas kinematics could directly test their tidal origin.
- Because the model assumes one recent starburst per bin, fitting the same outer-region spectra with a two-burst or non-parametric star-formation history would test whether the slow outer decline is actually two unresolved bursts — a check that could either strengthen or revise the sequential picture.
- The same outside-in starburst sequence might be visible in high-redshift post-starbursts, linking local merger remnants to the bursty galaxies seen in the early universe.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a hierarchical Bayesian model for fitting spatially resolved MaNGA spectra of three local post-starburst galaxies, simultaneously constraining per-Voronoi-bin stellar population parameters and galaxy-wide radial trends in stellar mass, burst age, pre-burst metallicity, dust attenuation, and velocity dispersion. The analysis uses a three-stage importance/rejection-sampling scheme built on nested sampling. The authors report that all three galaxies show an outer, older, weaker and more slowly quenching starburst followed by a central, younger, stronger and more rapidly quenching starburst, with the two episodes separated by roughly 1 Gyr. They interpret this as evidence for a recent gas-rich merger, associating the outer burst with first pericentre passage and the central burst with coalescence, and argue that the quenching is more consistent with gas consumption plus morphological stabilization than with AGN feedback. The paper also maps non-axisymmetric features in burst mass fraction and dust attenuation, and validates selected results against pipe3D, MaNGA DAP, and Balmer-decrement dust estimates.
Significance. If the inferred two-phase outside-in starburst sequence is correct, the paper provides a rare, spatially resolved fossil record connecting the merger timeline (first pericentre vs. coalescence) to the resolved SFHs of local post-starburst galaxies, and it strengthens the case for gas-exhaustion and morphological quenching over AGN feedback. The methodological contribution is substantial: the hierarchical model is carefully derived, the nested-sampling implementation is modern, and the authors provide extensive validation and uncertainty maps. The analysis is also commendably transparent, with several important limitations acknowledged explicitly in the text. However, the central interpretation rests on a single-burst-per-bin parametric SFH, and the radial functional form for burst age is chosen after inspecting fits to the same galaxies. These choices make the main claim more model-dependent than the narrative suggests. With only three galaxies, the empirical basis is also narrow, although the paper is framed as an initial proof-of-concept study.
major comments (3)
- [Section 4.1.1 / Table 3 / Section 6.2] The per-Voronoi-bin SFH contains exactly one recent starburst (old exponential plus single double power-law). The recovered radial gradient in t_burst is then interpreted as two distinct galaxy-wide episodes, but the paper itself states in Section 6.2 that the radial trends are 'also consistent with the outer regions experiencing two distinct starbursts, manifesting as a burst with a slow apparent decline rate.' This is load-bearing because the 'slower quenching' of the outer burst and the ~1 Gyr age gap are central claims. The stage-1 non-hierarchical fits in Fig. 3 show two t_burst populations, but those fits use the same one-burst-per-bin model, so they do not break the degeneracy. I would like to see a quantitative test: either fit with a two-burst or non-parametric SFH in the outer bins, or use simulated SFHs to demonstrate that the single-burst model does not turn a two-burst input
- [Section 4.1.2 / Section 4.3] The hierarchical model uses the data twice in a way that can shape the main result. Stage-0 global fits are used to set informative priors for individual Voronoi bins (Section 4.1.2), and the radial functional forms—especially the logistic t_burst profile in Section 4.3.2—are chosen after inspecting non-hierarchical fits of the same galaxies ('patterns we observed from the fitted parameters', Section 4.3). The paper defends the prior update as admissible empirical Bayes, but the choice of the logistic radial model is not given the same scrutiny. A monotonically increasing logistic profile will tend to produce an 'outer-older, central-younger' sequence by construction. I recommend a sensitivity analysis with a more flexible radial model (e.g., spline or free per-annulus ages) or a hold-out-bin validation to show that the two-phase sequence is not imposed by the adopted functional form.
- [Section 4.5 / Section 5] The paper assumes conditional independence between Voronoi bins and applies the PSF correction only to the stellar mass surface density model, not to the burst-age or other radial gradients. As acknowledged in Section 4.5, the reported radial profiles other than stellar mass are PSF-convolved gradients. This is directly relevant to the quantitative claim of a ~1 Gyr offset between outer and central bursts: PSF scattering of the bright, young central population into outer bins could bias the inferred outer burst age, even if it is unlikely to create a spurious sign reversal. The authors caution that the sign of the gradient is probably robust, but the timing of the sequence, which is central to the merger-timeline interpretation, could be affected. A simple PSF-forward-modelling test for the burst-age map would help quantify this bias.
minor comments (6)
- [Section 5.2] Galaxy ID typo: '12514-3792' should be '12514-3702'.
- [Section 6.4] Galaxy ID typo: '12607-3701' should be '12067-3701'.
- [Section 6.5] 'in contract to conclusions' should read 'in contrast to conclusions'. Also '7976-1902' appears to be a typo for '7965-1902'.
- [Section 4.1.1] Minor wording: 'risingslope' should be 'rising slope'.
- [Data Availability] The data availability statement uses placeholder 'url' for maps and scripts; actual repository links should be provided before publication.
- [Figure 9 caption] Typo 'Voroni bin' should be 'Voronoi bin'.
Circularity Check
No significant circularity: the two-phase sequence is read directly from fitted radial gradients, with an acknowledged one-burst degeneracy and non-load-bearing self-citations.
full rationale
The claimed outside-in starburst sequence is a summary of the fitted per-bin SFH parameters (t_burst, tau_1/2) and their radial trends, not a quantity predicted from an independent input. Stage 1 non-hierarchical fits (Fig. 3, Section 4.2) already show two t_burst populations and a stacked two-episode SFH before the hierarchical radial model is introduced; the stage 2 logistic t_burst(R) is chosen after inspecting these same fits (Section 4.3), so the hierarchical result is a shrinkage/regularization of the same data rather than a construction. The priors on t0 and t∞ are identical, so the recovered t∞>t0 (outer older than center) is not forced by the functional form. Section 6.2 explicitly acknowledges the one-burst-per-bin parametric SFH cannot distinguish a slow outer decline from two distinct outer bursts ('the radial trends shown in Fig. 8 are therefore also consistent with the outer regions of the galaxies experiencing two distinct starbursts, manifesting as a burst with a slow apparent decline rate'); this is a genuine model degeneracy and a limitation, but it is not a circular reduction, and the authors propose future two-burst fitting. Self-citations (Leung et al. 2024; Wild et al. 2020) supply the SFH ansatz and fitting machinery, but the central interpretation is compared against external simulations (Petersson et al. 2023; Zheng et al. 2020) and independent observations (Cortijo-Ferrero et al. 2017; French et al. 2018), so no load-bearing self-citation chain forces the result. The stage-0 'more informative priors' are mildly data-informed but do not fix the burst-age gradient; the paper openly discusses this practice. Overall, no step reduces by definition to its inputs.
Assumptions & free parameters
free parameters (20)
- Re (half-mass radius) =
3.65, 8.92, 3.15 kpc for 7965-1902, 12067-3701, 12514-3702
- log10 Sigma_e (formed mass surface density at Re) =
8.23, 7.40, 8.58
- n_sersic =
3.91, 5.15, 3.63
- sigma_logM (scatter in log stellar mass) =
0.0840, 0.0996, 0.0952
- t0 (burst age at R=0) =
0.29, 0.98, 0.71 Gyr
- t_inf (burst age at large radius) =
1.25, 1.89, 1.58 Gyr
- Rc (inflection radius of logistic burst-age profile) =
1.62, 2.37, 1.98 kpc
- Rw (scale width of logistic burst-age profile) =
0.7, 0.6, 0.2 kpc
- sigma_t (scatter in burst age) =
0.127, 0.187, 0.187 Gyr
- m_logZ (pre-burst metallicity gradient) =
-0.069, -0.098, -0.023 dex/kpc
- Z_old0 (pre-burst metallicity at R=0) =
0.96, 1.31, 1.19 Zsun
- sigma_logZ (scatter in log metallicity) =
0.040, 0.049, 0.047 dex
- A_V0 (ISM dust attenuation at R=0) =
1.13, 0.17, 1.30 mag
- A_Vmax (ISM dust attenuation at Rmax) =
0.27, 0.27, 0.35 mag
- sigma_AV (scatter in A_V) =
0.217, 0.256, 0.286 mag
- m_log_sigma (velocity dispersion gradient) =
-0.097, -0.097, -0.031 dex/kpc
- sigma_disp0 (velocity dispersion at R=0) =
86.0, 125.8, 151.1 km/s
- sigma_log_sigma (scatter in log velocity dispersion) =
0.011, 0.050, 0.009 dex
- eta (birth-cloud dust factor) =
MAP from global fit per galaxy (not tabulated)
- Local per-Voronoi-bin spectral parameters (M*, tburst, Zold, Zburst, fburst, A_V, sigma_disp, SFH shape, redshift) =
Posterior per bin; maps in Figs. 5-7
assumptions (6)
- domain assumption Each Voronoi bin's spectrum is described by the two-component parametric SFH: an old exponential component plus a single recent double-power-law starburst, with fixed rising slope beta=250.
- domain assumption The Bruzual & Charlot (2003) stellar population models (MILES templates), Kroupa IMF, two-step metallicity evolution, and two-component dust law are adequate to separate age, metallicity and dust in these spectra.
- ad hoc to paper The radial population models are correct: Sersic for formed stellar mass, logistic for burst age, linear gradients in log Z, A_V and log sigma, with Gaussian azimuthal scatter.
- domain assumption Voronoi bins are conditionally independent given the population model; PSF and spatial covariance are neglected for all parameters except stellar mass.
- standard math The three-stage importance-sampling/rejection-sampling scheme produces samples equivalent to the joint posterior of the hierarchical model.
- domain assumption Absence of broad-line AGN features and non-Seyfert BPT/WHAN classifications rule out a significant AGN contribution to quenching.
Cite this review
Pith. "Pith review of Spatially resolved star-formation histories of local post-starburst galaxies: Starburst and quenching spatial patterns consistent with recent mergers." pith.science (2026). https://pith.science/paper/XCL4CUZ3
@misc{pith2026260213114,
author = {Pith},
title = {Pith review of: Spatially resolved star-formation histories of local post-starburst galaxies: Starburst and quenching spatial patterns consistent with recent mergers},
year = {2026},
howpublished = {\url{https://pith.science/paper/XCL4CUZ3}},
note = {Machine review of arXiv:2602.13114}
}
abstract
Post-starburst (PSB) galaxies, having recently experienced a starburst followed by rapid quenching, are excellent laboratories to probe physical mechanisms that drive starbursts and shutting down of star formation. Integral-field spectroscopy reveals the galaxies' spatially-resolved properties, where observed directional patterns can be linked to the galaxies' past evolution. We measure the resolved star-formation histories (SFHs), stellar metallicity evolution and dust properties of three local PSBs from the MaNGA survey, down to $0.5$" resolution ($\sim0.3\,$kpc) using a hierarchical Bayesian model. Local parameters were constrained simultaneously with parameters describing spatial trends. We found that all three galaxies first experienced an outer, weaker and slower quenching starburst, followed by a central, stronger and faster quenching starburst that peaked $\sim 1\,$Gyr after the first. The central starbursts induced a significantly stronger rise in stellar metallicity compared to the outer starbursts. These results are consistent with the effects of a recent gas-rich (wet) merger, where the first pericentre passage triggered starbursts in the outer regions, while the later coalescence triggers a stronger centralised starburst. We find non-axisymmetric features in the maps of burst mass fraction and dust attenuation in all galaxies, which could be caused by tidal effects during the recent merger. Comparisons with literature binary merger simulations suggests that the galaxies' rapid quenching was driven by gas consumption and the stabilisation against gas gravitational collapse by a growing spheroid, while AGN feedback was not necessarily a primary cause.
Figures
Figures from the paper (6 more)
Forward citations
Cited by 1 Pith paper
-
Towards Unbreaking the Universe: MINERVA Measurements of Color Gradients in Massive Quiescent Galaxies Can Help Ease Too-Early Star Formation Tensions
Resolved photometry of four high-redshift quiescent galaxies reveals negative color gradients that lower estimated stellar masses by 0.1 dex relative to slit measurements, reducing model tensions under an age-driven i...
Reference graph
Works this paper leans on
-
[1]
H., 2017, MNRAS, 466, 3272 Andreon S., Weaver B., 2015, Bayesian Methods for the Physical Sciences
AdamsM.R.,CornishN.J.,LittenbergT.B.,2012,Phys.Rev.D,86,124032 AlarconA.,SánchezC.,BernsteinG.M.,GaztañagaE.,2020,MNRAS,498, 2614 Alatalo K., et al., 2016, ApJS, 224, 38 Albán M., Wylezalek D., 2023, A&A, 674, A85 Alsing J., Heavens A., Jaffe A. H., 2017, MNRAS, 466, 3272 Andreon S., Weaver B., 2015, Bayesian Methods for the Physical Sciences. Learning fr...
arXiv 2012
-
[2]
(2004), and more recently Andreon & Weaver (2015) and Thrane & Talbot (2019)
For moredetailedderivations,thereaderisreferredtoBernardo&Smith (2000) and Gelman et al. (2004), and more recently Andreon & Weaver (2015) and Thrane & Talbot (2019). AswritteninSection3,thejointposteriordistributionof®𝛼and𝚯 (all®𝜃𝑖) is given by 𝑃(®𝛼,𝚯|D)= 𝑃(D|𝚯,®𝛼)𝑃(𝚯,®𝛼) 𝑃(D) .(A1) Using the definition of conditional probability𝑃(𝐴|𝐵)= 𝑃(𝐴,𝐵)/𝑃(𝐵), it f...
2000
-
[3]
Alsing et al
For examples applications of BHM in astronomy, the reader can refer to the growing literature in many fields, including in cosmol- ogy (e.g. Alsing et al. 2017; Hinton et al. 2019), gravitational waves (e.g. Adams et al. 2012; Thrane & Talbot 2019), stellar observa- tion and evolution (e.g. Leistedt & Hogg 2017; Hall et al. 2019; Lyttle et al. 2021), and ...
2017
-
[7]
This paper has been typeset from a TEX/LATEX file prepared by the author. MNRAS000, 1–24 (2026) Resolved star-formation histories of PSB galaxies25 log10(Σ ∗/M⊙ kpc□2) Burst Age / Gyr Burst mass fraction Half time τ1/2 / Gyr log10(Σ burst/M⊙ kpc□2) -5 0 5 ξ / arcsec -5 0 5 η / arcsec Zold/Z⊙ Zburst/Z⊙ AV / mag σdisp / km s□1 vLOS / km s□1 0.025 0.050 0.07...
2026
-
[8]
MNRAS000, 1–24 (2026)
2026
-
[11]
APPENDIX C: UNCERTAINTY MAPS For completeness, in Fig
Given that the sampler fromstage 2is converged and can be treated as a collapsed Gibbs sampler (van Dyk & Park 2008), it is reasonable to assume that the samples drawn through this rejection sampling method using the acceptanceprobabilityabovearethesameassamplesdrawndirectly from𝑃( ®𝜃𝑖|®𝐷𝑖,®𝛼𝑘). APPENDIX C: UNCERTAINTY MAPS For completeness, in Fig. C1 we...
2008
-
[2022]
Many studies have commented on the advantages of BHM in degeneracy-breaking (Juvela et al
and to study the stellar-halo mass relation with galaxy-galaxy lensing data (Sonnenfeld & Leauthaud 2018). Many studies have commented on the advantages of BHM in degeneracy-breaking (Juvela et al. 2013; Galliano 2018), uncertainty propagation (Curtis-Lake et al. 2021; Sandlesetal.2022),andreductionofbiasintheresults(Sonnenfeld & Leauthaud 2018). Leveragi...
2018
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.