{"id":"29668be8-5b28-4002-8f3b-3be28b9ceef9","arxiv_id":"2412.04541","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"First resolved mass-metallicity relation for 24 local extremely metal-poor galaxies reveals a metal-poor star-forming branch, attributed to episodic star formation fueled by dilute gas infall.","lead":"Using Subaru IFU observations of 24 extremely metal-poor galaxies, the authors map gas metallicity spaxel by spaxel and find that patches with more intense star formation are more metal-poor at fixed stellar density; they interpret this as pristine gas falling in and triggering episodic starbursts. The work extends resolved galaxy scaling relations into a regime that matters for understanding the earliest galaxies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"R3 metallicity calibration's EW(Hβ) proxy shares a parameter with Σ_SFR; residual bias could create the reported anti-correlation.","rationale":"The central claim is that spatially resolved metallicity decreases with increasing Σ_SFR at fixed Σ_* (Sect. 3.3, Fig. 4, Eq. 2) and that the metal-poor horizontal branch reflects metal-poor gas infall (Sect. 3.5). The most load-bearing step is the metallicity calibration in Sect. 3.1: R3 is corrected for ionization using EW(Hβ), while Σ_SFR is derived from dust-corrected Hα in Sect. 3.3. Because both Balmer lines trace recent star formation, EW(Hβ) and Σ_SFR are strongly correlated, so any residual calibration dependence on ionization or SFR will produce a spurious anti-correlation. The paper's Figure 1 validation uses integrated central spectra and cannot rule out a per-spaxel bias. The added 0.15 dex systematic uncertainty in the binned fit does not address a correlated bias. The proposed test—checking the metallicity–Σ_SFR relation within narrow EW(Hβ) bins—would cleanly separate a calibration artifact from a physical trend. Given the plausibility of the infall scenario, the partial Te-based validation, and the kinematic support from Isobe et al. (2023) and Xu et al. (2024), the concern is not fatal but is severe enough to warrant conditional acceptance pending this direct test. The reader's weakest assumption identifies exactly this issue, and our analysis agrees.","tokens_in":35882,"tokens_out":3736,"duration_ms":85261,"concrete_test":"Restrict the spaxel sample to narrow bins of EW(Hβ) (e.g., Δlog EW = 0.2 dex) and, within each bin at fixed log Σ_*, test whether 12+log(O/H) still decreases with log Σ_SFR. If the anti-correlation disappears or weakens substantially within fixed EW(Hβ) bins, the global trend is an artifact of the EW(Hβ)–Σ_SFR coupling rather than evidence for metal-poor gas infall.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.1 calibrates R3 = [OIII]λ5007/Hβ using EW(Hβ) as an ionization proxy (Nakajima et al. 2022), while Section 3.3 derives Σ_SFR from dust-corrected Hα following Kennicutt (1998). Hα and Hβ both trace ongoing star formation, so EW(Hβ) correlates strongly with Σ_SFR. If the R3 calibration retains any residual dependence on ionization or SFR that is not fully removed by the EW(Hβ) correction, the inferred metallicity will be systematically lower in high-Σ_SFR spaxels, directly generating the claimed anti-correlation. The validation in Figure 1 compares only integrated central spectra and does not test the per-spaxel residual trend. The 'remarkably tight' offset–Σ_SFR correlation (Fig. 4c) and the derived α = 0.66 in Equation (2) could be inflated by this systematic. In the extreme, a compressed low-metallicity calibration end could produce flat metal-poor branches without true metal-poor gas infall. This is an internal correctness risk, not a mere disagreement with external consensus.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents Subaru/FOCAS IFU observations of 24 extremely metal-poor galaxies (EMPGs), mapping gas-phase metallicity (via the R3 index with an EW(Hβ) ionization correction), stellar mass surface density (from 6500–7000 Å continuum), and SFR surface density (from dust-corrected Hα). The authors report that the resolved mass–metallicity relation (rMZR) shows large scatter at low Σ* and that this scatter is strongly correlated with Σ_SFR, with metallicity decreasing as Σ_SFR increases at fixed Σ*. They further identify a 'metal-poor horizontal branch' in 11 of 22 spatially resolved EMPGs, located at the peaks of Σ* and Σ_SFR and surrounded by more metal-enriched gas, which they interpret as evidence for episodic star formation fueled by metal-poor gas infall. Four of the most metal-poor galaxies show only a metal-poor clump, proposed as candidates for a first phase of chemical evolution. The paper extends the rMZR to metallicities 12+log(O/H) ≈ 6.9–7.9 and includes a comparison of the classification with kinematic properties from Isobe et al. (2023).","tokens_in":36179,"tokens_out":6053,"duration_ms":60023,"significance":"If the central results hold, this is the first systematic spatially resolved study of the mass–metallicity relation in the extremely metal-poor regime, a parameter space barely explored by MaNGA/CALIFA and high-redshift IFU surveys. The reported anti-correlation between metallicity and Σ_SFR at fixed Σ*, and the existence of metal-poor horizontal branches, would provide direct observational evidence for metal-poor gas infall driving episodic star formation in low-mass galaxies. The four Category B objects (Z ≲ 0.03 Z☉) are potentially valuable local analogs of early galaxy formation. The paper benefits from careful data reduction, a cross-check of R3 metallicities against direct Te abundances (Figure 1), BPT-based exclusion of evolved AGN, and the use of existing high-resolution kinematics to support the inflow scenario. However, the astrophysical conclusions rest heavily on the assumption that spaxel-level metallicity differences are not dominated by systematic trends in the strong-line calibration with SFR or ionization.","major_comments":[{"comment":"The R3 metallicity calibration uses EW(Hβ) as the ionization proxy, while Σ_SFR is independently derived from dust-corrected Hα. Because Hα and Hβ both trace recent star formation, EW(Hβ) and Σ_SFR are strongly correlated. If the R3 calibration retains any residual dependence on ionization or SFR beyond the EW(Hβ) correction, the inferred metallicities will be systematically lower in high-Σ_SFR spaxels, artificially producing or inflating the anti-correlation in Figure 4(c) and the derived α = 0.66 in Equation (2). The validation in Figure 1 uses only integrated central spectra and does not test the per-spaxel calibration trend. The authors should quantify this residual dependence, for example by comparing R3 metallicities to Te-based metallicities in spaxels where [OIII]λ4363 is detected, or by demonstrating that the offset–Σ_SFR correlation persists when the sample is restricted to a narrow EW(Hβ) range or when an alternative ionization correction (e.g., from [OIII]/[OII] where available) is adopted.","section":"§3.1, Figure 4(c), Eq. (2)"},{"comment":"The statistical significance of the offset–Σ_SFR correlation and the uncertainties on the fitted parameters are computed by treating the 9,177 spaxels as independent samples. However, spaxels within a single galaxy are spatially correlated, so the effective number of independent measurements is far smaller than the spaxel count. This likely overstates the significance of the trend and underestimates the uncertainties on α in Equation (2). Please re-analyze the correlation and fit using a bootstrap resampled at the galaxy level, or a hierarchical model, and report the significance at the number of independent galaxies (24, or 22 for the resolved subsample).","section":"§3.3, Figure 4"},{"comment":"The classification of the 22 resolved EMPGs into Categories A–D, and particularly the identification of the 'metal-poor horizontal branch' in 11 objects, appears to be based on visual inspection of rMZR plots. Since this classification is the basis for the central claim of episodic star formation fueled by metal-poor gas infall, the criteria should be made quantitative and reproducible. Please provide a clear operational definition of the branch (e.g., a minimum number of spaxels forming a flat sequence at low metallicity, a threshold on the metallicity drop near the Σ_SFR peak, or a quantitative measure of the contrast with surrounding spaxels) and demonstrate that the 11/22 classification is robust to reasonable variations of those criteria.","section":"§3.5, Figures 7–12"}],"minor_comments":[{"comment":"The sentence beginning 'Note that we do not apply slit-loss corrections...' would be clearer if the authors stated exactly which lines are used for the R3 metallicity and the continuum, since R3 = [OIII]λ5007/Hβ is a line ratio but EW(Hβ) also requires continuum.","section":"§3.1"},{"comment":"The phrase 'The FOCAS IFU data was reduced' should read 'The FOCAS IFU data were reduced'.","section":"§2.2"},{"comment":"The color bars in several figures (e.g., Figures 7–11) appear to be labeled 'log SFR' rather than 'log Σ_SFR'; please make the notation consistent and ensure the units (M☉ yr⁻¹ kpc⁻²) are always shown.","section":"Figure 4 and related panels"},{"comment":"Please state explicitly how the uncertainties on the best-fit parameters a, b, c, and α were computed, and whether the binning of spaxels or the added 0.15 dex systematic uncertainty was included in those uncertainties.","section":"§3.3, Eq. (2)"},{"comment":"The dependent variable y is not explicitly defined in the text; please state that y = 12 + log(O/H) in both Equations (1) and (2).","section":"§3.3, Eq. (1)"},{"comment":"The classification of metallicity gradients into negative, positive, and flat appears to be qualitative; please specify the threshold (e.g., slope uncertainty) used to assign each category.","section":"§3.4"},{"comment":"The histogram in Figure 13 uses literature central metallicities from Table 1, not the resolved metallicity maps; this is not immediately clear from the caption and should be stated more prominently.","section":"§4.1, Figure 13"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a very timely and important question with a unique dataset, and the overall interpretation is plausible. The main risk is the strong dependence of the central claim on the R3 calibration's EW(Hβ) ionization correction, which is correlated with the key physical quantity Σ_SFR. The authors should be encouraged to perform the additional tests requested in Major Comment 1, as these are feasible within the scope of the current dataset (e.g., using spaxels with [OIII]4363 detections, or split-sample tests in EW(Hβ)). I would also suggest that the authors make the spaxel catalogs (metallicity, Σ*, Σ_SFR, EW(Hβ)) publicly available, as this would greatly facilitate independent verification of the reported trends."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline result is genuinely new: this is the first systematic, spatially resolved mass–metallicity relation for extremely metal-poor galaxies, built from 24 EMPGs and 9,177 spaxels in a parameter space (Z ~ 6.9–7.9, Σ* ~ 1e5–1e7 Msun/kpc2) that large IFU surveys have not reached. The metal-poor horizontal branch seen in 11 of 22 resolved systems, and the four-category evolutionary sequence (A through D), are interesting and likely to be cited. The data reduction is careful and the paper is honest about several of its own limitations, including the two-valued R3 calibration near Z = 8.0, the 0.3 dex offset for the most metal-poor galaxy, and the possibility that some Category B classifications are just detection limits.\n\nThe soft spot is the calibration–SFR coupling. The R3 metallicity uses EW(Hβ) as an ionization proxy, while Σ_SFR is derived from Hα; Hα and Hβ both trace current star formation, so any residual dependence of the calibration on ionization or SFR could inflate or even create the reported anti-correlation between metallicity and Σ_SFR. The Te-based validation in Fig. 1 is only for integrated central spectra, not per spaxel, so it does not directly test this. That is a real internal-correctness risk, and the paper should be asked to address it with spaxel-level Te measurements or photoionization-model tests before the result is taken as established. The stress-test note is on point here, but I would not call the paper hopeless: the central trend is shown directly in Fig. 4c before any model fitting, the data are new, and the α = 0.66 from the fitted relation matches independent global results, which is reassuring.\n\nMinor issues: the definition of the metal-poor horizontal branch and the category assignments are visual and not yet quantitative; no data or code are released, which makes independent checks harder. These are fixable in revision.\n\nWho is this for? Anyone working on low-mass galaxy chemical evolution, local analogs of high-redshift systems, or resolved scaling relations. It deserves a serious referee: the observation is substantial, the interpretation is physically reasonable, and the main concern is testable rather than terminal. I would recommend sending it to review, with the explicit request that the authors validate the per-spaxel metallicity calibration against ionization- and SFR-related systematics.","headline":"First resolved mass-metallicity relation for extremely metal-poor galaxies, with a plausible but not yet bulletproof anti-correlation with star-formation surface density.","tokens_in":36986,"tokens_out":1539,"would_cite":true,"duration_ms":22856,"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":"Half of extremely metal-poor galaxies show fresh metal-poor cores inside metal-enriched surroundings, evidence that infalling gas drives their star formation in bursts.","keywords":["extremely metal-poor galaxies","spatially resolved metallicity","mass-metallicity relation","star formation surface density","metal-poor gas infall","episodic star formation","integral field spectroscopy","local dwarf galaxies"],"falsifier":"Measure spaxel-level metallicities in the same galaxies with an independent, direct method, using electron temperatures from $[\\mathrm{O\\,III}]\\.\\lambda4363$ where it is detected, and compare the metallicity versus $\\Sigma_{\\rm SFR}$ trend at fixed $\\Sigma_\\star$; if the trend and the metal-poor branch disappear under direct metallicities, the central claim is a calibration artifact.","tokens_in":35683,"feed_emoji":"🌌","tokens_out":7403,"duration_ms":71968,"temperature":0.7,"pith_summary":"The paper maps gas-phase oxygen abundance pixel by pixel across 24 of the most metal-poor galaxies known, then asks whether those maps match the usual picture of chemically young galaxies. It finds instead that metallicity falls where star-formation surface density rises at fixed stellar-mass surface density, and that in 11 of 22 resolved galaxies the metal-poorest gas sits at the star-forming core while the surrounding gas is more enriched. The authors read this \"metal-poor horizontal branch\" pattern as evidence that infalling metal-poor gas dilutes galactic centers and triggers episodic star formation. If correct, the resolved mass-metallicity relation reaches $12+\\log(\\mathrm{O/H})\\approx 6.9$–$7.9$, and local extremely metal-poor galaxies split into at least two chemical histories: recent-infall bursts and possibly first-light systems.","feed_headline":"Half of extremely metal-poor galaxies show infall-fed starburst cores","feed_subtitle":"Spatially resolved maps show metal deficits at star-formation peaks, a signature of episodic growth fed by pristine gas.","key_machinery":"The analysis rests on the R3 metallicity indicator, $R_3=[\\mathrm{O\\,III}]\\,\\lambda5007/\\mathrm{H}\\beta$, applied per spaxel and corrected for ionization using the H$\\beta$ equivalent width; this is the tool that turns emission-line maps into metallicity maps in a regime where $[\\mathrm{O\\,II}]$ and $[\\mathrm{O\\,III}]\\.\\lambda4363$ are often undetected. The other central object is the \"metal-poor horizontal branch\" on the resolved mass-metallicity diagram, a flat or slightly rising sequence of low-metallicity spaxels that reaches the $\\Sigma_\\star$ and $\\Sigma_{\\rm SFR}$ peaks and is surrounded by more metal-rich gas; this spatial pattern is the diagnostic through which infall and episodic star formation are inferred.","core_discovery":"Using the R3 strong-line index, $R_3=[\\mathrm{O\\,III}]\\,\\lambda5007/\\mathrm{H}\\beta$, calibrated on extremely metal-poor galaxies with the H$\\beta$ equivalent width as an ionization correction, the paper maps metallicities across 9,177 spatial pixels (spaxels) in 24 EMPGs. At fixed local stellar-mass surface density, spaxel metallicity decreases as local star-formation surface density increases, and the scatter of the resolved mass-metallicity relation is largely absorbed by a relation of the form $y=a+b(\\mu_\\alpha-c)e^{-(\\mu_\\alpha-c)}$ with $\\mu_\\alpha=\\log\\Sigma_\\star-\\alpha\\log\\Sigma_{\\rm SFR}$ and $\\alpha=0.66\\pm0.04$. Half of the resolved galaxies show a distinct, nearly horizontal metal-poor branch at the peaks of $\\Sigma_\\star$ and $\\Sigma_{\\rm SFR}$, surrounded by gas enriched to about 0.1–0.2 solar metallicity; four galaxies with $Z\\lesssim0.03\\,Z_\\odot$ show only the metal-poor branch. The paper interprets the branch-plus-envelope pattern as recent metal-poor gas infall that dilutes the center and fuels a starburst, consistent with turbulence-dominated kinematics and gas fractions near unity, and identifies the isolated-branch galaxies as likely first-phase chemical evolution systems.","pith_inferences":["A direct test would compare spaxel-level R3 metallicities against auroral-line metallicities in galaxies where $[\\mathrm{O\\,III}]\\.\\lambda4363$ is detected: the reported metallicity–$\\Sigma_{\\rm SFR}$ anti-correlation should persist if it is physical rather than a calibration artifact.","The Category A/B split predicts correlated abundance signatures, such as low N/O and high gas fraction concentrated in the metal-poor branch regions; the paper notes that forthcoming N/O work will test this.","The same rMZR–$\\Sigma_{\\rm SFR}$ dependence could be sought in JWST/NIRSpec resolved observations of $z>3$ galaxies, where rapid gas accretion should produce similar cold, low-metallicity cores surrounded by more enriched envelopes.","Mapping velocity and metallicity simultaneously in these systems could reveal whether the metal-poor branch gas is actually moving inward, giving a kinematic test of the infall interpretation."],"forward_implications":["The resolved mass-metallicity relation now reaches $12+\\log(\\mathrm{O/H})\\approx6.9$–$7.9$ and $\\Sigma_\\star\\approx10^5$–$10^7\\,M_\\odot\\,\\mathrm{kpc}^{-2}$, and its low-mass scatter is nearly closed by $\\Sigma_{\\rm SFR}$.","If infalling metal-poor gas dilutes cores and fuels bursts, many EMPGs are not chemically young; selection on strong Balmer lines preferentially catches post-infall starbursts.","The four Category B galaxies with isolated metal-poor clumps at $Z\\lesssim0.03\\,Z_\\odot$ are the best local analogues of early-universe galaxies.","Kinematic evidence ($v_{\\rm rot}/\\sigma<1$, gas fraction $\\approx0.9$–$1.0$) supports the infall-triggered starburst picture and implies detectable kinematic discontinuities at branch locations."],"supporting_citations":[{"why":"Calibrates the R3 strong-line metallicity indicator used for every spaxel, including the Hβ-equivalent-width ionization correction.","marker":"Nakajima et al. (2022)"},{"why":"Supplies the high-metallicity resolved mass-metallicity reference relation and its functional form used for extrapolation and offset measurements.","marker":"Barrera-Ballesteros et al. (2016)"},{"why":"Provides kinematics, inclinations, gas fractions, and rotation-to-dispersion ratios used to support the infall-triggered starburst interpretation.","marker":"Isobe et al. (2023)"},{"why":"Companion kinematic study reporting low vrot/σ values and gas-rich properties, including category-dependent trends.","marker":"Xu et al. (2024)"},{"why":"Provides EMPG stellar masses from SED fitting and the mass-to-light calibration used to derive Σ⋆ per spaxel.","marker":"Kojima et al. (2020)"},{"why":"Compiles the high-metallicity resolved relation data used for comparison and joint fitting.","marker":"Maiolino & Mannucci (2019)"},{"why":"Documents metallicity drops at starburst heads in tadpole galaxies, the local analogue used to support gas-inflow dilution.","marker":"Sánchez Almeida et al. (2014b)"}],"fun_headline_variants":["Pristine gas infall fuels starbursts in dwarf galaxy cores","Metal-poor gas infall drives episodic star formation in EMPGs","Half of EMPGs show metal-poor starburst cores from infall","Spatial maps link metal-poor gas to episodic starbursts","Ultra-metal-poor galaxies may mirror early galaxy evolution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that the metallicity calibration returns unbiased oxygen abundances for every pixel, with no residual hidden dependence on star-formation intensity; if star-forming pixels are systematically read as more metal-poor, the anti-correlation and the flat branch could be artifacts.","fun_headline_variants_meta":{"raw":{"variants":["Pristine gas infall fuels starbursts in dwarf galaxy cores","Metal-poor gas infall drives episodic star formation in EMPGs","Half of EMPGs show metal-poor starburst cores from infall","Spatial maps link metal-poor gas to episodic starbursts","Ultra-metal-poor galaxies may mirror early galaxy evolution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000919,"raw_usage":{"total_tokens":4085,"prompt_tokens":1227,"completion_tokens":2858,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":843,"completion_tokens_details":{"reasoning_tokens":2763}},"tokens_in":843,"tokens_out":2858,"duration_ms":19742,"temperature":1.0,"reasoning_tokens":2763,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:25:04.412995+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure spaxel-level metallicities in the same galaxies with an independent, direct method, using electron temperatures from $[\\mathrm{O\\,III}]\\.\\lambda4363$ where it is detected, and compare the metallicity versus $\\Sigma_{\\rm SFR}$ trend at fixed $\\Sigma_\\star$; if the trend and the metal-poor branch disappear under direct metallicities, the central claim is a calibration artifact.","supporting_citations":[],"review_version":1}