{"id":"48dbdceb-a004-4072-bff2-884cb77e3367","arxiv_id":"2501.09827","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In three local elliptical galaxies, the inner component formed early and rapidly and the outer component assembled later, consistent with the two-phase formation scenario.","lead":"Astronomers split three nearby elliptical galaxies into two light components and found that their inner parts are old and metal-rich, while the outer parts are younger and more metal-poor. The pattern supports the idea that the inner cores formed early and quickly, while the outer regions were added later by merging with smaller galaxies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 50% inner-component stellar mass fraction in the Abstract is asserted without derivation and is not reproduced by Table 3's r-band light fractions (34–45%); the central co-dominant claim needs a stated mass-to-light conversion.","rationale":"The reader's weakest assumption correctly identifies the BUDDI decomposition as model-dependent and lacking model-selection statistics, and I agree that the unreliable J205050 outer component is a genuine caveat. However, I find a more specific and more decisive gap: even if the decomposition is accepted, the headline 50% stellar-mass fraction does not follow from the published numbers. Table 3 is explicitly a flux ratio at 6166 Å, Table 5 is normalized within each component, and no equation for converting component fluxes to component stellar masses appears anywhere in the text. A constant-M/L assumption would itself conflict with the paper's own age and metallicity differences between the components, while a non-constant conversion is not documented. This makes the central quantitative claim unauditable, which is a stronger reason for caution than decomposition uncertainty alone. The qualitative two-phase and inside-out narrative is plausible and broadly consistent with the literature, so I would not reject the paper outright; the appropriate outcome is the same CONDITIONAL verdict the reader reached, with the added requirement that the mass-fraction calculation be explicitly derived and the abstract's 50% be either substantiated or revised. My recommendation therefore remains UNCHANGED relative to the reader's verdict.","tokens_in":35551,"tokens_out":5675,"duration_ms":56779,"concrete_test":"Reconstruct the stellar mass of each component from the BUDDI extracted spectra: integrate each component SED weighted by the pPXF mass weights, or equivalently combine the best-fit M/L at the 6166 Å BUDDI flux with Table 3, to compute M_inner/(M_inner+M_outer) for each galaxy, propagating M/L uncertainties from the pPXF fits. State the equation used and compare to 50%: if the derived fractions differ by more than 5–10 percentage points for any galaxy, the abstract and Sect. 4.4.2 must be revised to the correct co-dominant fraction; if the conversion cannot be unambiguously reconstructed from the current text, the claim is unreviewable as written. A useful secondary check is to rerun pPXF on the same BUDDI spectra with the regularization parameter changed by one order of magnitude and recompute τ50 and τ90 to test whether the inner-first assembly ordering survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim that the inner component contributes ~50% of the stellar mass is asserted, not derived, and it is not supported by the only tabulated quantity that bears on it. Table 3 gives r-band (6166 Å) flux ratios F_inner/F_total = 0.34 (J020536), 0.45 (J205050), and 0.34 (J225546); these are light fractions, not mass fractions. Section 4.4.2 states that the relative contribution of the inner component to the total stellar mass was found to be 50% for all three galaxies, but no equation, mass-to-light conversion, or uncertainty is provided. Table 5 normalizes the pPXF weights within each component, so it cannot supply the inter-component mass ratio. The equality would hold only under identical mass-to-light ratios in both components, but Table 4 shows the components differ substantially in age and metallicity (e.g., J205050 mass-weighted inner age 13.80 Gyr at [M/H]=0.245 versus outer age 8.83 Gyr at [M/H]=-0.426), so a constant-M/L assumption is not stated or justified. If the 50% value was obtained from a more elaborate conversion of the BUDDI component spectra, that step is absent from the text and cannot be checked. Because the abstract and conclusions are framed around a co-dominant ~50% outer component, this missing derivation is the most load-bearing gap. The paper's own caveats (weaker J205050 outer fit, old-age degeneracy) compound but do not cause the problem.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a pilot study of three isolated local elliptical galaxies (J020536, J205050, J225546; z<0.06, M* ~ 10^11 M_sun) observed with MUSE. The authors decompose each datacube into two Sérsic components (inner and outer) using BUDDI, a spectro-photometric decomposition code, and fit the extracted component spectra with pPXF to derive mass- and light-weighted ages, metallicities, and star formation histories, including cumulative SFHs and tau50/tau90 assembly timescales. The central claims are that (i) the inner components are old and metal-rich, assembled rapidly at early epochs via dissipative collapse or gas-rich major mergers; (ii) the outer components are co-dominant, contributing approximately 50% of the total stellar mass, and assembled their mass a few gigayears later via dry mergers and possible gas accretion; and (iii) these results support the two-phase formation scenario with inside-out growth. The paper also compares the BUDDI-based results with Voronoi-binned 2D stellar population maps and Lick indices as consistency checks.","tokens_in":35851,"tokens_out":10831,"duration_ms":100036,"significance":"If supported, the paper's qualitative results--metal-rich old inner components, more extended and less enriched outer components, inside-out growth--would add an IFS-based, spectro-photometric decomposition perspective to the two-phase formation scenario at intermediate mass scales. The methodology is clearly described, and the paper is honest about its caveats: lunar contamination for J205050's outer component (Sect. 3.2), the old-age degeneracy (Sect. 4.2), and the statement in Sect. 5.3 that the fossil-record method does not directly identify in-situ versus ex-situ populations. Strengths include the tabulated structural parameters, the decomposed spectra in the appendix, and the qualitative agreement across all three galaxies and with the Voronoi-binned maps. However, the quantitative backbone of the paper--the approximately 50% inner stellar mass fraction and the 1-5 Gyr assembly delays--rests on an undocumented mass calculation and on the least reliable component in the sample, so the significance is currently conditional on a revised, fully derived mass analysis.","major_comments":[{"comment":"The central claim of the Abstract and Section 6, that the inner component 'contributes to ~50% of the galaxy stellar mass' and that the outer component is therefore 'co-dominant', is asserted in Sect. 4.4.2 ('The relative contribution of the inner component to the total stellar mass was found to be 50% for all three galaxies') without any derivation. The only tabulated quantity bearing on this is the r-band (6166 Å) flux ratio F_inner/F_total in Table 3, which is 0.34, 0.45, and 0.34 for J020536, J205050, and J225546; these are light fractions, not mass fractions. Converting them to a mass fraction requires a mass-to-light ratio for each component, and Table 4 shows the components differ by up to ~5 Gyr in age and ~0.7 dex in metallicity (e.g., J205050: inner age 13.80 Gyr at [M/H]=0.245 versus outer 8.83 Gyr at [M/H]=-0.426), under which a constant M/L cannot be assumed without explicit justification. Table 5 cannot supply the inter-component ratio because the pPXF weights are normalized within each component. The authors should state the M/L conversion, present the resulting inner mass fractions with uncertainties, and either support or temper the 'co-dominant' wording in the Abstract.","section":"Sect. 4.4.2; Table 3; Abstract"},{"comment":"The two-component Sérsic model is the foundation of the entire analysis, since all stellar population and SFH results derive from the BUDDI-decomposed component spectra, yet its selection over single- and three-component alternatives is justified only by 'visual inspection of the fit residuals and the estimated structural parameters, along with the extracted spectra' (Sect. 4.1), and the Chebyshev polynomial orders in Table 2 are likewise chosen by visual inspection. BUDDI/GalfitM computes goodness-of-fit statistics, so the authors should report a quantitative model comparison (e.g., delta-chi^2 or BIC between one-, two-, and three-component fits) or otherwise demonstrate that the decomposition is not driven by the polynomial-order choices, sky model, and mask geometry. Without this, the physical separation into 'inner' and 'outer' components, and hence the in-situ/ex-situ interpretation, rests on an unquantified modelling assumption.","section":"Sect. 3.1, Sect. 4.1"},{"comment":"The quantitative assembly-delay result depends most heavily on the object the authors themselves caution against. Sect. 3.2 warns that the outer component of J205050 has the strongest lunar contamination and that 'we caution the direct interpretation of the fits and the resulting stellar populations for this component'. Yet this is the component that produces the largest inner-outer tau50 delay in Table 4 (4.74 Gyr versus 2.09 Gyr and 0.67 Gyr for the other two galaxies) and the largest intermediate-age mass fraction in Table 5 (35.8%), and it drives the 'delay of several gigayears' narrative in Sect. 4.4.2. With a sample of only three galaxies, the paper should either demonstrate through a sky-subtraction sensitivity test that the lunar contamination does not bias the derived ages and metallicities, or down-weight or exclude this component when stating the quantitative delays.","section":"Sect. 3.2; Table 4; Table 5"}],"minor_comments":[{"comment":"The final sentence of Sect. 4.1 is grammatically incomplete and contains an unmatched parenthesis: 'the 1D profiles rarely capture mismatches with the surface brightness models (depending on the choice of sampling rate and that are clearly visible spatially in the 2D profiles.'; it should be reworded.","section":"Sect. 4.1"},{"comment":"The phrase 'the inner component had already formed half of their mass by ~13 Gyr ago' has a singular/plural agreement error, and 'a delay of several gigayears' is broader than the tabulated tau50 delays, which range from 0.67 Gyr (J225546) to 4.74 Gyr (J205050); please rephrase to match Table 4.","section":"Sect. 4.4.2"},{"comment":"Several tau50 and tau90 values in Table 4 exceed the age of the Universe under the adopted Planck 2016 cosmology (e.g., J205050 inner tau50 = 14.00 Gyr and single-component tau50 = 13.96 Gyr versus ~13.8 Gyr), consistent with the old-age degeneracy noted in Sect. 4.2; this caveat should be restated in Sect. 4.4.2 where the assembly timescales are used quantitatively.","section":"Sect. 4.2; Table 4"},{"comment":"The 1D surface brightness profiles are extracted along the position angle of the single-Sérsic model, but Table 3 gives appreciably different PAs for the two components (e.g., J020536: inner 6.05 deg versus outer -8.47 deg); the authors should state whether the conclusions drawn from the 1D profiles (e.g., the ~3.5-4 kpc inflection discussed in Sect. 5.1) are robust to this choice.","section":"Sect. 4.1; Fig. 1; Sect. 5.1"},{"comment":"The colour scales of the 2D SFH grids differ between sub-panels and between galaxies (e.g., the inner-component mass-fraction scale reaches 0.4 for J205050 but only 0.1 for J020536), making visual cross-comparison of absolute mass fractions difficult; consider using a common colour scale or stating explicitly that each panel is individually normalized.","section":"Sect. 4.4; Figs. 7-9"},{"comment":"The identification of the inner component with in-situ star formation and the outer component with ex-situ accretion is an interpretive step; the authors themselves note in Sect. 5.3 that the fossil-record method does not directly characterise a stellar population as formed in situ or ex situ, so this caveat should be restated in the Abstract or Conclusions to avoid overclaiming the two-phase framing.","section":"Sect. 5.3; Abstract"}],"recommendation":"major_revision","confidential_remarks":"The paper is a reasonable fit for A&A as a pilot study, but the abstract overstates what the body supports: the '~50%' co-dominance claim is not derived anywhere, and the quantitative delay signal is carried mostly by the one component the authors flag as less reliable. Both issues are fixable in revision. I would encourage the editors to ask for either a release of the BUDDI component spectra or a quantitative model-selection table, since the decomposition code is developed partly by the authorship team and the two-component choice is otherwise hard to check independently. The citation practice appears appropriate and no novelty-disclosure concerns arose."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper gives you genuinely new component-resolved stellar populations for three isolated ellipticals from MUSE, and the qualitative picture—old, metal-rich inner component; older but more spread outer component—is believable. But the headline claim that the inner component contributes ~50% of the stellar mass is asserted, not shown. Table 3 lists r-band light fractions of 34–45% for the inner component; Table 5 normalizes the pPXF weights within each component. There is no stated mass-to-light conversion, no equation, no uncertainty. That gap matters because the abstract and conclusions are built around the “co-dominant” 50/50 picture. The components differ substantially in age and metallicity, so you cannot quietly assume equal M/L.\n\nCredit where it is due: the BUDDI decomposition plus pPXF is a clean methodological demonstration. They also do useful due diligence—Lick indices as a sanity check, Voronoi maps, explicit masking of a nearby dwarf and a background galaxy, and an honest discussion of the bright-moon conditions. The J205050 outer component is explicitly flagged as less robust, which is good practice. The reconstructed SFHs, taken qualitatively, are a reasonable read.\n\nSoft spots: beyond the missing mass derivation, the choice of two Sérsic components rests on visual inspection rather than any model comparison. That is acceptable for a pilot, but it is a free parameter. The 50% number appears in the abstract as if measured; it needs to be derived or demoted. Also, the outer component of J205050 was observed under strong lunar contamination and the paper itself cautions against direct interpretation—yet the sample-wide conclusions include it. That does not sink the paper, but the authors should show the results with and without that galaxy.\n\nWho this is for: people using spectro-photometric decomposition, and observers working on elliptical assembly. It is a pilot study, not a definitive result. I would send it to review, but require the mass fraction derivation and a clearer separation of supported versus inferred quantities.","headline":"The component-resolved data are new and the qualitative picture is plausible, but the paper's central ~50% inner-mass claim is asserted rather than derived, and the tabulated light fractions do not support it.","tokens_in":36450,"tokens_out":1961,"would_cite":false,"duration_ms":20657,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Three local ellipticals show a two-phase, inside-out assembly history in their stellar populations.","keywords":["elliptical galaxies","stellar populations","star formation histories","integral field spectroscopy","galaxy decomposition","two-phase formation","Sersic profiles","MUSE"],"falsifier":"Re-observing the galaxy with the strongest lunar contamination and repeating the decomposition under dark-sky conditions would settle whether its outer component's young, metal-poor population is real, since the paper itself flags this component as less robust. More generally, if a single-Sérsic fit or a three-component fit with a formal model comparison outperforms the two-component fit on the same datacubes, the inside-out two-phase interpretation would lose its structural basis.","tokens_in":35327,"feed_emoji":"🌌","tokens_out":7298,"duration_ms":67576,"temperature":0.7,"pith_summary":"This paper tries to show that the formation of elliptical galaxies can be read as two distinct phases, visible in the stars themselves. Using integral-field spectroscopy of three isolated, roughly $10^{11} M_\\odot$ ellipticals at $z<0.06$, the authors decompose each galaxy into an inner and an outer component and extract each component's spectrum separately. They find that the inner component formed early and rapidly, is old and metal-rich, and contains roughly half of the galaxy's stellar mass. The outer component assembled most of its mass shortly afterward, through dry mergers and possible gas accretion, and hosts a mix of old and intermediate-age stars with lower metallicities. The result supports the two-phase scenario of elliptical assembly and demonstrates a method for separating in-situ from ex-situ stellar populations.","feed_headline":"Elliptical galaxies assembled inside-out in two phases","feed_subtitle":"Inner cores formed early and fast; outer envelopes later accreted half the stellar mass.","key_machinery":"The load-bearing tool is a spectro-photometric decomposition code, BUDDI, that fits two Sérsic profiles to every wavelength slice of the MUSE datacube, thereby separating the integrated light into a clean one-dimensional spectrum for each component. These component spectra are then fitted with a full spectral fitting routine (pPXF) that combines simple stellar population templates to return mass-weighted and light-weighted ages, metallicities, and the full distribution of stellar populations on an age–metallicity grid, from which cumulative star formation histories and $\\tau_{50}$/$\\tau_{90}$ formation and quenching timescales are derived. Voronoi-binned 2D maps provide a consistency check of the radial trends. The key identity is the decomposition itself: the two Sérsic components are interpreted as the in-situ and ex-situ phases, so the physical conclusions depend on the decomposition being faithful.","core_discovery":"On the paper's own terms, the central discovery is that in all three galaxies the stellar populations separate cleanly into two structurally distinct components with different assembly histories. The inner component (Sérsic half-light radius $\\sim 1.4\\text{--}2.8$ kpc) is composed entirely of stars older than about 8 Gyr, is metal-rich, and built up half of its mass by a look-back time of roughly 13 Gyr, implying a rapid, early dissipative collapse or gas-rich major merger. The outer component ($R_e \\sim 6.7\\text{--}19$ kpc) contributes about 50% of the total stellar mass, but assembled its mass more slowly, reaching $\\tau_{90}$ several gigayears later, with a substantial intermediate-age population and sub-solar metallicities, consistent with accretion of smaller systems via dry mergers. The reconstructed star formation histories show both components quiescent over the last several gigayears, with the inner core quenched first. The authors interpret this as inside-out mass growth under the two-phase scenario, where the inner component is in-situ and the outer component is largely ex-situ.","pith_inferences":["A testable extension would be to apply the same two-component decomposition to a larger, mass- and environment-matched sample; if the roughly 50% inner-mass fraction is universal, it would set a strong constraint on merger histories.","The interpretation of the outer component as ex-situ rests on matching stellar populations to merger simulations; a direct check would be to search for tidal streams or shells around these galaxies, which should be present if dry mergers built the envelope.","Because the fossil record cannot distinguish stars formed in situ from those accreted, the same data could be reinterpreted with a single, radially continuous star formation history; the two-phase claim would be strengthened by showing a real discontinuity in mass-weighted age and metallicity at the component boundary."],"forward_implications":["If the two-phase picture holds, the inner component of a local elliptical is the surviving high-redshift core, so its stellar population should match the compact red nuggets seen at $z \\gtrsim 1.5$.","The outer component carrying roughly half of the stellar mass means that dry merging after quenching is not a minor addition but a co-dominant channel of mass growth.","The short $\\tau_{50}$–$\\tau_{90}$ interval in the inner component implies that star formation was intense and quenched early, so later gas accretion did not reignite star formation.","The negative metallicity gradient between components (metal-rich inner, metal-poor outer) traces the in-situ/ex-situ transition rather than a continuous disk-like gradient."],"supporting_citations":[{"why":"Supplies the two-phase simulation framework that the interpretation invokes.","marker":"Oser et al. 2010"},{"why":"Extends the two-phase simulations to predict stellar population properties of accreted versus in-situ components.","marker":"Johansson et al. 2012"},{"why":"Established multi-component decomposition of ellipticals and linked inner components to in-situ formation.","marker":"Huang et al. 2013a"},{"why":"Presents the BUDDI code used for spectro-photometric decomposition.","marker":"Johnston et al. 2017"},{"why":"Presents the pPXF full spectral fitting method used to derive stellar populations.","marker":"Cappellari & Emsellem 2004"},{"why":"Provides the MILES simple stellar population template library used for age and metallicity estimation.","marker":"Vazdekis et al. 2010"},{"why":"Connects structural components to in-situ and accreted origins through surface brightness and stellar population analysis.","marker":"Spavone et al. 2017"}],"fun_headline_variants":["Ellipticals assembled inside-out in two phases","Two-phase assembly: inner core early, outer halo late","MUSE reveals dual assembly histories in ellipticals","Inside-out growth: cores first, halos accreted later","Fast early core, slow accreted halo in ellipticals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim collapses if the two-Sérsic decomposition does not yield physically clean, uncontaminated spectra for each component, since all stellar population and star formation history results are derived from those spectra.","fun_headline_variants_meta":{"raw":{"variants":["Ellipticals assembled inside-out in two phases","Two-phase assembly: inner core early, outer halo late","MUSE reveals dual assembly histories in ellipticals","Inside-out growth: cores first, halos accreted later","Fast early core, slow accreted halo in ellipticals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000337,"raw_usage":{"total_tokens":1944,"prompt_tokens":1104,"completion_tokens":840,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":762}},"tokens_in":720,"tokens_out":840,"duration_ms":8964,"temperature":1.0,"reasoning_tokens":762,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T19:37:54.820625+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-observing the galaxy with the strongest lunar contamination and repeating the decomposition under dark-sky conditions would settle whether its outer component's young, metal-poor population is real, since the paper itself flags this component as less robust. More generally, if a single-Sérsic fit or a three-component fit with a formal model comparison outperforms the two-component fit on the same datacubes, the inside-out two-phase interpretation would lose its structural basis.","supporting_citations":[{"cited_title":"H., Naab, T., & Ostriker, J","cited_arxiv_id":null,"evidence_quote":"Extends the two-phase simulations to predict stellar population properties of accreted versus in-situ components."},{"cited_title":"R., et al","cited_arxiv_id":null,"evidence_quote":"Connects structural components to in-situ and accreted origins through surface brightness and stellar population analysis."}],"review_version":1}