{"id":"12a464b8-08f7-45b7-9eb9-a2c5cf40b8d0","arxiv_id":"2502.02080","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Simulated galaxies in the CIELO suite frequently show broken (non-linear) metal abundance profiles, and the type of break traces specific gas accretion and feedback processes.","lead":"A new fitting algorithm called DB-A was used to detect kinks, or breaks, in the metal abundance profiles of 45 simulated galaxies, and it found that most of them have broken rather than smooth profiles. The paper links each break type to a physical process, such as gas inflows, star formation, feedback outflows, or cold gas accretion, using two galaxies as detailed case studies.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The physical process attributions in the abstract are inferred from only two case-study galaxies, with no quantitative linkage between break type and measured inflow/outflow/accretion properties; the generality of the diagnostic claim is therefore not established.","rationale":"The reader's weakest-assumption concern about the subjective manual revision in DB-A is valid and important, and it is one source of uncertainty in the break classification. However, the more load-bearing concern is that even a perfectly objective break classifier would not establish the paper's central claim: the physical process attributions come from only two galaxies and are based on qualitative inspection of temporal sequences, accretion fractions, and δr-δz diagrams. The z=0 sample of 45 galaxies provides the taxonomy of profile shapes, but no quantitative connection between those shapes and the proposed physical drivers. The paper explicitly defers a statistical study to future work, yet the abstract and conclusions present the attributions as general findings. This is a scope-of-evidence problem rather than an internal inconsistency. The proposed test would directly assess whether the case-study associations hold across a larger sample. If they do, the diagnostic claim is strengthened; if not, the conclusions should be limited to the two galaxies. Since the reader already recommended CONDITIONAL and my concern reinforces that recommendation, the verdict remains unchanged, with the condition being that the process attributions be rephrased as case-study-based or extended to a statistically meaningful sample.","tokens_in":26901,"tokens_out":2619,"duration_ms":30878,"concrete_test":"Apply the same time-resolved break classification and gas-source tracking used for LG1-4469 and LG1-4337 to a larger sample of at least 10 additional CIELO galaxies (ideally all 45), and compute, for each classified break type (inner rise, inner drop, outer rise, outer drop), the distribution of measured physical quantities: dominant accreted gas source fraction (ISM inflow, CGM, unbound/cold flow, satellite), SFR surface density in the relevant region, and outflow fraction (δr<0 gas). Then perform a statistical test (e.g., permutation test or logistic regression) of whether break type predicts the dominant source and SFR enhancement at, say, p<0.05. If the associations do not replicate or are not statistically significant, the abstract's diagnostic claims must be rephrased as case-study illustrations rather than general results.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim extends the case-study findings to general diagnostics: 'Most inner breaks show central oxygen enhancement, linked to gas accretion and star formation. Inner drops result from disrupted gas due to feedback-driven outflows. Outer breaks with high metallicities arise from re-accreted material...' These attributions are drawn from a detailed temporal analysis of exactly two galaxies, LG1-4469 and LG1-4337 (Sections 4-5). The z=0 statistics from 45 galaxies (Section 3) quantify the frequency of profile shapes, but they do not measure the physical processes responsible for those shapes. The process links are established qualitatively: for example, inner rises are associated with 'gas inflows' and 'efficient star formation' based on one or a few temporal sequences (Fig. 10), and outer drops are linked to 'metal-poor gas accretion from cold flows' using two intervals (Fig. 12 upper panels). No quantitative test is provided that the classified break type predicts the dominant accretion source, outflow strength, or SFR enhancement across the sample. Since the paper's stated purpose is to individualize the physical processes, a two-galaxy demonstration is suggestive but insufficient to support the abstract's population-level diagnostic statements. The authors themselves acknowledge this scope in Section 5 ('These results will be used to define a strategy to perform a statistical study in a forthcoming paper'), but the abstract and conclusions present the associations as established findings. This is a load-bearing concern because if the two case studies are not representative, the central claim that profile shape can be used as a diagnostic of gas accretion, feedback, and CGM mixing is unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the CIELO cosmological zoom simulations to characterize the gas-phase oxygen abundance profiles of 45 central galaxies at z=0 and to interpret profile breaks physically. The authors introduce a new fitting algorithm, DB-A, which fits up to two broken power laws and classifies profiles as linear, inner-broken, outer-broken, or doubly broken. At z=0 they report 38 per cent linear, 13 per cent inner, 40 per cent outer, and 9 per cent double breaks, with gradient values consistent with MUSE/CALIFA/MaNGA observations. They then follow two galaxies in a Local Group analog across about 87 snapshots to z~6, finding that inner rises are associated with gas inflows and subsequent star formation, inner drops with feedback-driven outflows, outer rises with re-accretion and CGM mixing, and outer drops with metal-poor cold-flow accretion. The paper concludes that metallicity-profile shapes can serve as diagnostics of these processes.","tokens_in":27166,"tokens_out":6145,"duration_ms":59909,"significance":"If the claims hold, the paper provides a simulation-based interpretation of broken metallicity profiles that ties break types to specific baryon-cycle processes, and it offers a new automatic fitting tool that is benchmarked against Theil-Sen single-slope fits and against observational break statistics. The authors are candid about the main limitations: no AGN feedback, restriction to sub-Milky-Way masses, and the case-study nature of the process attributions. The strength of the paper is its use of full cosmological simulations with metal-dependent cooling, star formation, and SN feedback, and its detailed temporal tracking of two galaxies, including a quantitative decomposition of accretion sources. The main risks are the subjectivity in the break-detection step and the extrapolation from two case studies to population-level statements.","major_comments":[{"comment":"The description of the DB-A algorithm leaves a load-bearing subjective step unspecified. The text states that \"conflictive cases (i.e. very small or very large break radius or inner and outer break radius too close to each other) were manually revised, decreasing the number of breaks used to fit the profiles,\" but no quantitative thresholds are given for \"very small,\" \"very large,\" or \"too close.\" Because every downstream result — the z=0 profile fractions (Section 3), the temporal break statistics (Section 4: 26 per cent inner, 43 per cent outer, 72 per cent inner rises, 68 per cent outer rises), and the case-study attributions (Section 5) — depends on the number and location of breaks, this manual intervention could bias all of them. The authors should specify the criteria in advance and should show that the reported fractions are stable when the revision step is replaced by an objective rule (for example, requiring the break radius to lie within a fixed range and requiring separated inner and outer breaks).","section":"2.2"},{"comment":"The fitted gradients and break radii are presented without any uncertainty estimates. Table A.1 lists single values for each galaxy, and the time series in Figure 8 are plotted without error bars, although the underlying profiles are noisy, as Figure 1 shows broad abundance distributions at fixed radius. Without uncertainties one cannot assess whether the 38/40/13/9 per cent distribution is distinguishable from, say, a 30/40/15/15 distribution, nor whether the 72 per cent inner-rise fraction or the redshift dependence of outer drops is statistically significant. The authors should provide bootstrap or resampling errors on the gradients and break radii and propagate them to the quoted fractions.","section":"2.2 and Table A.1"},{"comment":"The robustness of DB-A to its hand-set parameters is not demonstrated. The running-median window N=40, the KDE outlier threshold 0.05, the tolerance tau=0.99, and nmax=2 are chosen once and used everywhere; the only validation reported is a statement that BIC and AIC give similar results. Since the break-detection method is the paper's central tool, the authors should show how the z=0 fractions and the temporal classifications change under plausible variations of N, tau, and the outlier threshold (for example, N=20 and N=80), and they should describe the BIC/AIC comparison in enough detail to be checked.","section":"2.2 and Section 3"},{"comment":"The abstract and conclusions generalize the physical attributions to the population: \"Most inner breaks show central oxygen enhancement, linked to gas accretion and star formation. Inner drops result from disrupted gas due to feedback-driven outflows. Outer breaks with high metallicities arise from re-accreted material...\" However, the quantitative link between break type and measured inflow/outflow/accretion properties is established only for the two case-study galaxies LG1-4469 and LG1-4337 (Sections 4-5). The z=0 sample (Section 3) provides frequencies of profile shapes but not measurements of the physical processes responsible for those shapes. The authors themselves note in Section 5 that these results \"will be used to define a strategy to perform a statistical study in a forthcoming paper.\" As it stands, the paper should either soften the population-level wording or add a quantitative test, such as comparing the dominant accretion source, outflow strength, or SFR enhancement between break types across the full sample.","section":"Abstract and Sections 4-5"}],"minor_comments":[{"comment":"The abbreviation B17 is used in the phrase \"B17 and SM18\" but is never defined; it presumably refers to Belfiore et al. 2017, which should be cited explicitly at that point.","section":"2.2"},{"comment":"There are several typos that should be corrected in a revision: \"Form the upper panels\" (Sec. 5.4), \"ouflows\" (Sec. 4.2), \"statical study\" (Sec. 5), \"high redshit\" (Sec. 5.1), \"shirring sphere\" (Sec. 2.1), \"understating\" (Conclusions), and a duplicated \"showing showing\" in the caption of Fig. 12.","section":"Throughout"},{"comment":"The table header labels the columns as \"Inner, Middle, Outer\" but the caption does not explicitly say that the middle value is the DB-A middle-region gradient and that linear profiles use the full [epsilon, 2 r83] range; please state this in the caption.","section":"Table A.1"},{"comment":"The observational comparison would be more informative if the SM18 break categories (linear, inner, outer, double) were shown in the same panels; currently only their linear-gradient values are overplotted in the lower panels.","section":"Figure 3"},{"comment":"The binomial probability quoted for the difference in break fractions between CIELO and SM18 should be accompanied by the mass-range caveat already mentioned in the text, and a chi-square or Anderson-Darling test on the category counts would be more natural than a binomial test on mutually exclusive fractions of one sample.","section":"Section 3"},{"comment":"The sign convention delta_r < 0 for outward flow is opposite to the usual radial-coordinate convention and is easy to misread; please state it explicitly in the captions of Figs. 9 and 12 as well as in the text.","section":"4.1.3 and Figs. 9/12"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of A&A and the authors are transparent about the simulations' limitations. My main editorial concern is the calibration of the abstract and conclusions to the evidence: the population-level diagnostic claims should be softened or supported by a quantitative multigalaxy test. The manual revision step in DB-A and the absence of uncertainties on the fitted parameters are fixable with additional analysis and should be addressed before publication. I see no issue with citation practice or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a read if you work on metallicity gradients: the DB-A algorithm is a sensible, tested way to fit broken profiles, and the z=0 census of 45 CIELO galaxies (38% linear, 40% outer, 13% inner, 9% double) is new. The two-galaxy temporal analysis is well done, with careful separation of accreted gas sources (ISM inflow/outflow, CGM, unbound/cold flows, satellites) and clear figures showing how inflow plus starburst makes inner rises, feedback makes inner drops, recycled CGM gas makes outer rises, and metal-poor cold flows make outer drops. That is the real content.\n\nThe main soft spot is scope. The abstract states these associations as if they were population-level findings, but they are inferred from exactly two galaxies. The authors themselves note in Section 5 that a statistical study is forthcoming. The associations are plausible and consistent with the literature, but they are not yet a demonstrated diagnostic tool. That should be fixed by qualifying the abstract and conclusions.\n\nSecond, DB-A has a manual revision step with no stated thresholds. The authors say 'conflictive cases' were manually revised and the number of breaks reduced, but do not say what counts as conflictive. This matters because the z=0 break fractions could shift. I would want objective criteria or a sensitivity test. Third, the fitted slopes and break radii are reported without uncertainties; a bootstrap or jackknife would help, though for a first presentation this is minor.\n\nCredit where due: the comparison with observations (SM18, CALIFA, etc.) is careful, the robustness check against BIC/AIC is good, and the paper is transparent about missing AGN feedback and the sub-Milky Way mass range. The citation pattern is solid.\n\nFor simulators and observers interested in broken metallicity profiles, the algorithm and census are directly usable. Overall, a solid, honest paper with a useful new tool and suggestive physical results, but the abstract overreaches the evidence. A serious referee can fix that. I would accept it for peer review with the expectation of major revision.","headline":"Useful new fitting algorithm and honest case-study physics, but the abstract sells a two-galaxy result as a general diagnostic.","tokens_in":27814,"tokens_out":2245,"would_cite":true,"duration_ms":23745,"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":"Broken shapes of galaxy metallicity profiles can be read as a record of gas inflow, star formation, and feedback-driven outflows, not as fitting artifacts.","keywords":["galaxies: abundances","galaxies: evolution","galaxies: formation","galaxies: ISM","metallicity gradients","galactic feedback","cosmological simulations","oxygen abundance profiles"],"falsifier":"Re-fit the 45 simulated oxygen profiles with an independent, fully automated piecewise-linear fitter using a fixed model-selection rule and no manual case-by-case adjustments; if the resulting fractions of linear, inner-broken, outer-broken, and doubly-broken profiles differ materially from 38/13/40/9 per cent, the break taxonomy and the physical interpretations attached to it are not robust.","tokens_in":26675,"feed_emoji":"🌌","tokens_out":13816,"duration_ms":124506,"temperature":0.7,"pith_summary":"The paper aims to show that the radial profile of oxygen abundance in a galaxy's star-forming gas is a structured record of how gas enters, enriches, and leaves the galaxy, not just a smooth gradient. Analysing 45 simulated central galaxies with masses between $10^{8.5}$ and $10^{10.5}$ solar masses at $z=0$, the authors find that only about 38 per cent of profiles are single power laws; the rest contain inner breaks, outer breaks, or both, and the gradient values match observations. Tracking two Local Group analogue galaxies back to $z\\sim6$, they tie inner rises to gas inflow followed by star formation, inner drops to feedback-driven outflows of enriched gas, outer rises to re-accretion of enriched circumgalactic gas and extended star formation, and outer drops to metal-poor cold-flow accretion at high redshift. The payoff is that the shape of a metallicity profile can be used as a diagnostic of the gas cycle inside galaxies.","feed_headline":"Trace galaxy metallicity breaks to gas inflows and outflows","feed_subtitle":"Each kind of break is tied to gas inflow, star formation, or feedback-driven outflows.","key_machinery":"The load-bearing tool is the DB-A algorithm, a piecewise-linear fitter that removes outliers with a Gaussian kernel density estimate, allows up to two breaks in the oxygen profile, and selects the number of breaks by comparing relative fit errors against a tolerance parameter ($\\tau=0.99$); it defines the inner, middle, and outer slopes that the analysis maps onto physical processes. The second mechanism is the gas bookkeeping applied to the two case-study galaxies: gas particles entering each radial region are classified by origin (unbound cold flows, circumgalactic gas, satellite gas, inflowing or outflowing interstellar gas) and by radial and vertical displacement, with star formation and accretion rates computed per region. This combination turns a one-dimensional abundance profile into an accounting of where the gas comes from, where it travels, and what it does when it arrives.","core_discovery":"The central discovery is a break taxonomy with a physical cause attached to each entry. Using the DB-A algorithm, which fits piecewise power laws with up to two breaks to the running-median oxygen profile, the paper classifies the $z=0$ sample into 38 per cent linear, 40 per cent outer-broken, 13 per cent inner-broken, and 9 per cent doubly-broken profiles, with middle-region slopes that agree with single-gradient measurements in MUSE and CALIFA galaxies. In the two case studies traced across time, inner breaks occur in 26 per cent of snapshots (72 per cent inner rises, 28 per cent inner drops) and outer breaks in 43 per cent (68 per cent outer rises, 32 per cent outer drops). The authors attribute inner rises to cold gas falling into the centre followed by a starburst that raises central oxygen, inner drops to supernova-driven ejection of enriched gas, outer rises to the return of enriched material from the circumgalactic reservoir plus extended star formation and tidally driven mixing, and outer drops, found almost exclusively at $z>1$, to accretion of pristine cold-flow gas. The claim is that each break shape advertises a specific episode in the cycle of accretion, star formation, and outflow.","pith_inferences":["The same break taxonomy could be applied to stellar-population abundance maps rather than gas; if the gas-cycle interpretation is right, starlight-weighted profiles should show smoothed echoes of the same breaks, shifted by the ages of the stars.","Because these simulations lack AGN feedback and stop at $10^{10.5}$ solar masses, the near absence of inner drops may be a mass-range effect; running DB-A on simulations with AGN feedback or on more massive galaxies would test whether inner drops become as common as observed.","The tolerance parameter and the two-break maximum are modeling choices; varying them systematically would reveal whether the 38/13/40/9 per cent split is stable or an artifact of the selection rule, and would calibrate the method for noisier observational data.","Observational outer breaks are likely under-counted because of surface-brightness limits; if profile shape is a true gas-cycle diagnostic, deeper IFU observations at fixed mass should find a higher outer-break fraction than current surveys report."],"forward_implications":["Middle-region gradients from DB-A agree with single linear fits over $[0.5,1.5]\\,r_{50}$, so past observational comparisons based on one gradient remain meaningful even when the true profile is broken.","Inner breaks last about 1.5 Gyr in the massive case study versus about 0.3 Gyr in the lower-mass one, so break persistence carries information about how deep the galaxy's potential well is.","Inner drops and inverted central gradients are produced by the same feedback-driven ejection of enriched gas, so they should appear together and be accompanied by disturbed gas morphology.","Outer drops are a high-redshift feature in these simulations, tied to metal-poor cold flows; mapping them at $z>1$ with deep spectroscopy would test whether the simulated cold-flow channel operates in real galaxies.","Breaks migrate outward in units of $r_{50}$ as enrichment proceeds, so the location of a break can serve as a rough clock for recent accretion and star-formation episodes."],"supporting_citations":[{"why":"Supplies the observed MUSE-based classification of inner, outer, and double breaks in local galaxies that the simulated break fractions and radii are compared against.","marker":"Sánchez-Menguiano et al. 2018"},{"why":"Establishes the mass dependence of MaNGA metallicity gradients and the central/outer flattenings that motivate the re-accretion and feedback interpretations.","marker":"Belfiore et al. 2017"},{"why":"Simulation study showing gas inflows first dilute central abundances and then re-steepen the gradient via star formation; underpins the inner-rise scenario.","marker":"Perez et al. 2011"},{"why":"Shows how interaction-driven inflows, star formation, and outflows modify metallicity gradients and open spiral arms, used for outer mixing and outer rises.","marker":"Sillero et al. 2017"},{"why":"Provides the EAGLE-based expectation of near-zero median gradients with growing scatter at high redshift and gradient recovery timescales that frame the evolutionary analysis.","marker":"Tissera et al. 2022"},{"why":"Defines the cold-flow gas accretion channel invoked to explain high-redshift outer drops of metal-poor gas.","marker":"Ceverino et al. 2016"},{"why":"Supplies the baryon-cycle framework of inflows, outflows, and circumgalactic gas that organizes the interpretation of re-accretion and galactic fountains.","marker":"Péroux & Howk 2020"},{"why":"Reports observed broken profiles in barred galaxies whose inner and outer slopes are used to validate the simulated gradient ranges.","marker":"Chen et al. 2023"},{"why":"Provides CALIFA inner-break statistics and the inner-drop correlation with specific star formation rate used to compare the simulated inner features.","marker":"Cardoso et al. 2025"}],"fun_headline_variants":["Galaxy metallicity breaks reveal gas inflow and outflow history","Metallicity breaks trace gas accretion, star formation, and outflows","Break taxonomy ties galaxy metallicity profile shapes to gas flows","Galaxy chemical breaks flag gas inflows and outflows"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classification of every profile as linear, inner-broken, outer-broken, or doubly-broken rests on the DB-A algorithm's decision about the true number and position of breaks, with ambiguous 'conflictive' cases settled by manual revision rather than a fixed objective rule; if those decisions are biased, the reported break fractions and the process attributions built on them shift.","fun_headline_variants_meta":{"raw":{"variants":["Galaxy metallicity breaks reveal gas inflow and outflow history","Metallicity breaks trace gas accretion, star formation, and outflows","Break taxonomy ties galaxy metallicity profile shapes to gas flows","Galaxy chemical breaks flag gas inflows and outflows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001481,"raw_usage":{"total_tokens":6043,"prompt_tokens":1130,"completion_tokens":4913,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":746,"completion_tokens_details":{"reasoning_tokens":4844}},"tokens_in":746,"tokens_out":4913,"duration_ms":32577,"temperature":1.0,"reasoning_tokens":4844,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T13:24:54.490248+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the 45 simulated oxygen profiles with an independent, fully automated piecewise-linear fitter using a fixed model-selection rule and no manual case-by-case adjustments; if the resulting fractions of linear, inner-broken, outer-broken, and doubly-broken profiles differ materially from 38/13/40/9 per cent, the break taxonomy and the physical interpretations attached to it are not robust.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Simulation study showing gas inflows first dilute central abundances and then re-steepen the gradient via star formation; underpins the inner-rise scenario."},{"cited_title":"B., Lambas, D","cited_arxiv_id":null,"evidence_quote":"Shows how interaction-driven inflows, star formation, and outflows modify metallicity gradients and open spiral arms, used for outer mixing and outer rises."}],"review_version":1}