{"id":"56b78911-7778-44a9-ba20-824b7a2e787b","arxiv_id":"2602.13114","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":20,"one_line_summary":"All three post-starburst galaxies show an earlier, weaker, slowly quenching starburst in their outskirts followed about 1 Gyr later by a stronger, faster-quenching central starburst.","lead":"This paper maps the star-formation histories of three recently quenched galaxies using integral-field spectra and a hierarchical Bayesian model. It finds that each galaxy's outer regions formed stars first and quenched slowly, while the centre starbursted about one billion years later and quenched fast—a pattern the authors interpret as the signature of a gas-rich merger.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed two-phase outside-in sequence rests on a one-burst-per-bin SFH; outer-bin 'slow quenching' could be an artifact of merging two bursts into one decline.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the per-bin SFH model contains exactly one recent starburst, and the radial gradient in burst age is interpreted as a two-episode galaxy-wide sequence. I agree that this is the most exposed part of the central claim. The paper's own Section 6.2 admits that the outer-region trends are also consistent with two distinct bursts per bin, so the 'slower quenching' and the clean first-pericentre/coalescence timeline are not uniquely determined by the current analysis. Other concerns (small sample, PSF convolution, posterior model selection) are real but less decisive: they affect precision and framing, whereas the one-burst-per-bin assumption could change the qualitative interpretation of the resolved SFHs. The stage-1 non-hierarchical fits provide some support for two spatially separated burst populations, but they share the same single-burst-per-bin restriction, so they cannot resolve this degeneracy. A mock-data stress test or a two-burst re-fit of the outer bins would settle whether the observed broad, slow decline is genuine or an artifact. Since the reader already assigns a CONDITIONAL verdict with moderate confidence, my read does not require changing that verdict; it reinforces the condition under which the central claim should be accepted.","tokens_in":41883,"tokens_out":6978,"duration_ms":70169,"concrete_test":"Construct mock MaNGA-like datacubes for a galaxy whose true resolved SFH has two recent bursts in every outer bin (a weak burst at t≈1.4 Gyr and a stronger one at t≈0.4 Gyr), with known radial mass/light profiles and PSF; run the full three-stage pipeline (Sections 4.2–4.4) on the mock. If the single-burst hierarchical model recovers a single outer burst with slow τ1/2 and an inner-outer age gap ~1 Gyr, the paper's interpretation is not identifiable from these data. Alternatively/additionally, re-fit the real outer Voronoi bins with a two-burst parametric SFH (e.g., French et al. 2018) and compare Bayesian evidence against the one-burst model; if the two-burst model is preferred in outer bins, the 'slower quenching' conclusion fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that each Voronoi bin's SFH contains at most one recent starburst. The model in Section 4.1.1/Table 3 uses an old exponential plus a single recent double power-law (β fixed at 250) for every bin. The recovered radial gradient in t_burst (Section 5.2, Fig. 8) is then read as two galaxy-wide episodes: an outer, slower-quenching burst followed by a central, faster one ~1 Gyr later. However, if an outer bin actually underwent two bursts—e.g., a weak first-pericentre burst and a later, coalescence-related burst—the single-burst model will tend to absorb both into one broad component with a slow apparent decline. The paper explicitly acknowledges this in Section 6.2: '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.' Thus the 'slower quenching' of the outer burst—one of the three qualitative features in the abstract—is exactly the observable most vulnerable to this degeneracy, and the ~1 Gyr age gap could also be biased if the fitted outer burst age is a luminosity-weighted blend of two episodes. The stage-1 non-hierarchical fits in Fig. 3 provide some independent support for two t_burst populations in 7965-1902, but those fits use the same one-burst-per-bin SFH, so they do not break the degeneracy. This is not a fatal flaw, but it is load-bearing: the merger-timeline interpretation (first pericentre vs coalescence) and the slow-vs-fast quenching contrast both depend on the outer bins' SFH shape being a single genuine burst.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":42505,"tokens_out":4059,"duration_ms":42368,"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":[{"comment":"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":"Section 4.1.1 / Table 3 / Section 6.2"},{"comment":"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":"Section 4.1.2 / Section 4.3"},{"comment":"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.","section":"Section 4.5 / Section 5"}],"minor_comments":[{"comment":"Galaxy ID typo: '12514-3792' should be '12514-3702'.","section":"Section 5.2"},{"comment":"Galaxy ID typo: '12607-3701' should be '12067-3701'.","section":"Section 6.4"},{"comment":"'in contract to conclusions' should read 'in contrast to conclusions'. Also '7976-1902' appears to be a typo for '7965-1902'.","section":"Section 6.5"},{"comment":"Minor wording: 'risingslope' should be 'rising slope'.","section":"Section 4.1.1"},{"comment":"The data availability statement uses placeholder 'url' for maps and scripts; actual repository links should be provided before publication.","section":"Data Availability"},{"comment":"Typo 'Voroni bin' should be 'Voronoi bin'.","section":"Figure 9 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is methodologically ambitious and transparent, but the central astrophysical claim is currently more model-dependent than the abstract suggests. The authors' own Section 6.2 caveat is important and needs to be addressed quantitatively before the two-phase merger interpretation can be accepted. The empirical-Bayes circularity in prior and radial-model construction is also worth a sensitivity test. I see no basis for rejection, but the requested additions (two-burst or non-parametric SFH tests, flexible radial model checks, and PSF bias quantification) are necessary to make the central claim robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a careful, honest methods paper on three MaNGA post-starburst galaxies, and the headline claim — an outer, weaker, slower-quenching starburst followed about 1 Gyr later by a central, stronger, faster-quenching one — is plausible but conditional. The paper basically says so itself in Section 6.2: the radial trends are also consistent with outer regions having experienced two distinct bursts that show up as one slow decline. That degeneracy sits right on the 'slower outer quenching' part of the claim, so I would not treat the two-phase sequence as settled until the SFH model allows more than one burst per bin.\n\nWhat is actually new and good. The hierarchical Bayesian machinery is properly presented — three-stage sampling via importance sampling and rejection sampling, with the collapsed-Gibbs justification spelled out. The stage-1 non-hierarchical fits in Fig. 3 are the most convincing piece: they show two distinct burst-age populations before any population model is applied, and that the global fit misses the outer, older burst entirely. That missed-burst result is the cleanest takeaway — integrated SED fits can lose a whole starburst episode. The cross-checks against pipe3D, MaNGA DAP, and Balmer-decrement dust are the right ones, and the paper is unusually open about its weaknesses (PSF-convolved gradients, conditional independence between bins, half-mass radii extrapolated beyond the observed range).\n\nSoft spots, in proportion. The one-burst-per-bin assumption is the load-bearing one, and the paper's own Section 6.2 is the right summary of it. The empirical-Bayes circularity — stage-0 global fits setting per-bin priors, and the radial functional forms chosen after inspecting the same galaxies' fits — is real but moderate, because the stage-1 bimodality shows the two-phase signal is not a construction of the hierarchical model. The more concrete problem is the Data Availability section: placeholders saying data and scripts are 'available at url'. For a methods-heavy paper without the code, the reproducibility burden is unmet.\n\nWho this is for: people working on PSBs, merger-driven quenching, or resolved spectral fitting. A referee should ask for the code and derived maps, and for a robustness test allowing two bursts in the outer bins. I would send it to review.","headline":"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.","tokens_in":43031,"tokens_out":4580,"would_cite":true,"duration_ms":40803,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["post-starburst galaxies","resolved star-formation histories","galaxy mergers","galaxy quenching","Bayesian hierarchical models","integral-field spectroscopy","stellar metallicity","starbursts"],"falsifier":"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.","tokens_in":41751,"feed_emoji":"💥","tokens_out":4555,"duration_ms":43288,"temperature":0.7,"pith_summary":"The paper tries to establish that the spatially resolved star-formation histories of three local post-starburst galaxies all follow the same two-phase sequence: an earlier, weaker, slowly quenching starburst in the outer regions, and a later, stronger, faster-quenching starburst in the centre that peaks about one billion years after the first. The authors argue this pattern matches what gas-rich galaxy mergers do — the first close passage lights up the outskirts, and the final coalescence drives a central burst. They further claim that the rapid quenching in these galaxies is best explained by gas consumption plus the stabilising effect of a growing spheroid, not necessarily by feedback from a supermassive black hole. A sympathetic reader would care because the result turns post-starburst galaxies into readable records of the merger timeline, and it challenges the default assumption that AGN feedback is the main quencher.","feed_headline":"Two starbursts a billion years apart drive these post-starburst galaxies","feed_subtitle":"Resolved star-formation histories tie the outer burst to first flyby and the central burst to final coalescence.","key_machinery":"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.","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"fun_headline_variants":["Two starbursts, one billion years apart, trace a galactic merger","First outer starburst, then central: a merger's double signature","Merger's first flyby sparks outer burst, coalescence ignites core","Galaxy's billion-year starburst gap reveals a wet merger","Outer starburst leads, central bursts a gigayear later: merger evidence"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Two starbursts, one billion years apart, trace a galactic merger","First outer starburst, then central: a merger's double signature","Merger's first flyby sparks outer burst, coalescence ignites core","Galaxy's billion-year starburst gap reveals a wet merger","Outer starburst leads, central bursts a gigayear later: merger evidence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000231,"raw_usage":{"total_tokens":1368,"prompt_tokens":833,"completion_tokens":535,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":577,"completion_tokens_details":{"reasoning_tokens":435}},"tokens_in":577,"tokens_out":535,"duration_ms":5539,"temperature":1.0,"reasoning_tokens":435,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T23:34:56.680445+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}