{"id":"2e33e613-e9f2-4b8a-9a31-ea4f1cda9fdf","arxiv_id":"2505.14781","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Multi-component emission line fitting of NGC 613 reveals bar-driven gas inflow along the dust lanes and an AGN-driven biconical ionised outflow aligned with the radio jet.","lead":"This paper analyzes MUSE spectra of the nearby barred galaxy NGC 613, separating multiple gas components and revealing gas flowing inward along the bar and a biconical outflow aligned with a radio jet. It also lays out a detailed pipeline for emission line fitting that could be reused on other integral-field surveys.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Residual-velocity zero point: the bar's non-circular streaming could masquerade as inflow, so the inflow claim needs an independent non-circular flow model.","rationale":"The reader's weakest assumption is precisely the one I identify: the single-disc circular model may not fully remove the gravitational signal in a strongly barred galaxy, so the residual velocity map is not automatically a measure of inflow. I agree with the reader's CONDITIONAL verdict because the inflow/outflow narrative is plausible and supported by prior work, but the inflow channel needs a more direct, model-independent test. The outflow detection and mass-rate estimates carry their own systematic uncertainties (density, C factor, metallicity calibration), but those are explicitly disclosed by the authors and do not affect the geometric alignment as severely. Thus the most load-bearing concern is the residual-velocity zero point. The proposed DiskFit test would settle whether the dust-lane residuals are genuine radial inflow or a projection of non-circular orbits. If the test confirms radial inflow, the paper's central claim is materially strengthened; if not, the claim should be downgraded to 'non-circular motions consistent with, but not uniquely requiring, inflow.' Since the reader already made the verdict CONDITIONAL, my read leaves the verdict unchanged.","tokens_in":37553,"tokens_out":3394,"duration_ms":36067,"concrete_test":"Fit the MUSE V_L1_gas map with an explicit non-circular kinematic model, e.g., DiskFit, allowing radial inflow and tangential streaming simultaneously with the bar potential's higher-order harmonics. If the best-fit radial velocity component along the bar dust lanes is inward with an amplitude comparable to the residual pattern in Fig. 14 (≳20–30 km/s at matched positions), the inflow claim is independently supported; if the pattern is absorbed by azimuthal or harmonic terms, the residual-based inflow claim is not yet established. This directly tests whether the single-disc circular model is the right zero point.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the translation of the residual map ΔV = V_L1_gas − V_circ (Fig. 14, middle) into evidence of gas inflow. V_circ is a Kinemetry model of a single inclined rotating disc, with inclination and PA adopted from Buta et al. (2015) and fitted to the stellar velocity field. NGC 613 hosts a strong bar with box/peanut, a nuclear disc/ring, and a possibly double-barred structure; in such a potential the stellar and gas orbits are not circular. Kinemetry on the stellar field will absorb some of the bar's non-circularity into a best-fit 'circular' model, and the residual map can contain line-of-sight signatures of x1/elliptical orbits and projection effects that resemble, but are not, radial inflow. The paper does not quantify how much of the residual is due to the bar distortion: the inflow interpretation therefore rests on the strong assumption that a single-disc circular model is the correct zero point. This is a correctness risk, not a contradiction: the dust-lane residual pattern is qualitatively what inflow produces, and the ALMA/torque analysis of Audibert et al. (2019) supports inward migration, so the claim is plausible. But the central narrative depends on this subtraction, and the same residual pattern could persist in a purely orbital non-circular model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a detailed analysis pipeline for MUSE TIMER data of the barred galaxy NGC 613, combining full-spectrum stellar fitting with multi-component Gaussian fitting of emission lines, AIC-based model selection, and spatially resolved diagnostics of extinction, ionisation, metallicity, kinematics, and outflows. The two central astrophysical claims are (i) that residual gas velocities relative to a single-disc circular model trace bar-driven inflows along the dust lanes, and (ii) that a biconical [O III] outflow aligned with the VLA radio jet is powered by the AGN, with an ionised outflow mass of log10 M_ion = 4.97 ± 0.01 M_sun and a mass outflow rate of about 0.04 ± 0.01 M_sun/yr. The methodology also includes a regularisation-setting test for SFH recovery and a multi-stage approach to avoid local minima in emission-line fitting. The paper is largely a case study whose central narrative is a complete gas cycle from bar inflow to nuclear feedback.","tokens_in":37692,"tokens_out":2966,"duration_ms":31413,"significance":"If the two claims hold, the paper provides a spatially resolved, multi-phase picture of gas inflow and AGN-driven outflow in a single nearby barred galaxy, and it strengthens the observational case for bar-driven secular feeding. The methodological contribution is also useful: the combination of global differential-evolution optimisation, outlier detection, RBF-based initialisation for spaxel-by-spaxel fitting, and AICc selection is a sensible and moderately novel recipe for multi-component emission-line work, and the scalar-field normalisation test for regularisation is a practical insight. The paper is honest about several caveats, stating that the outflow mass uncertainties are only statistical lower limits and explicitly questioning whether the low metallicity of the outflow is an artifact. However, the inflow claim depends on the residual-velocity zero point being a correct circular-disc model, and the outflow energetics rest on adopted n_e and C values; both are load-bearing assumptions that need more quantitative support before the central claims can be accepted at face value.","major_comments":[{"comment":"The inflow interpretation is built on the residual map V_L1_gas - V_circ, where V_circ is a Kinemetry model of a single inclined rotating disc with inclination and PA adopted from Buta et al. (2015) and fitted to the stellar velocity field. NGC 613 is strongly barred, has a box/peanut structure and a nuclear disc/ring, and the paper itself notes the nuclear disc and bar-related kinematic features. In such a potential, stellar and gas orbits are not circular, and Kinemetry on the stellar field can absorb part of the bar's non-circularity into a best-fit 'circular' model. The residual pattern along the dust lanes is qualitatively what inflow models produce, but the same pattern could in principle arise from x1/elliptical orbits or projection effects in a non-circular potential. The paper does not quantify how much of the residual amplitude is attributable to the bar distortion or to uncertainties in inclination and PA. Because this residual map is the primary evidence for the paper's signature inflow claim, I ask for a concrete robustness test: for example, a harmonic decomposition of the gas velocity field, a parametric non-circular flow model fitted to the gas kinematics, or at least a variation of i and PA over their published uncertainties and a statement of how the dust-lane residual pattern changes.","section":"Section 7.4, Fig. 14"},{"comment":"The outflow mass and mass outflow rate rest on three adopted quantities: a single average electron density n_e ≈ 112 cm^-3 derived from [S II] in the central kiloparsec, a condensation factor C = 1, and the O3N2-based oxygen abundance. The quoted uncertainty of ±0.01 dex in log10 M_ion is only the statistical propagation and is explicitly acknowledged in the text as a lower limit. However, the paper then presents M_ion ≈ 1e5 M_sun and Mdot_ion ≈ 0.04 M_sun/yr as central values without a systematic error budget. Since [S II] ratios in the field can range from saturation-limited regimes (n_e between roughly 50 and 2000 cm^-3) and C is typically between 0.1 and 1 in ionised outflows, the systematic uncertainty in M_ion is at least an order of magnitude. I request that the authors provide a propagated systematic range for M_ion, Mdot_ion, and E_dot_out, and explicitly state which values of n_e and C would be needed to change the conclusion that the AGN can power the outflow.","section":"Section 7.5, Eqs. (11)-(13)"},{"comment":"The metallicity used in Eq. (11) comes from O3N2, and the paper itself notes in Section 7.3 that the apparent metal-poor outflow seen with O3N2 may be an artifact of ionisation-parameter sensitivity; the alternative calibrations shown in Fig. 12 give different abundances in the outflow region, and for some indexes the outflow is outside the calibration range. Since the outflow mass scales as 10^-[O/H], this systematic uncertainty propagates directly into the mass and rate quoted in the abstract. I ask the authors to either recompute the outflow properties using the metallicity range allowed by the various calibrators, or to state clearly which calibrator was used and why, and how the outflow mass would change under the alternative estimates.","section":"Section 7.3 and 7.5"}],"minor_comments":[{"comment":"The word 'dissussed' should be 'discussed' in the sentence describing scalar vs field normalisation.","section":"Section 4.2"},{"comment":"The caption and text use 'scalar field normalisation' for what is elsewhere called 'field normalisation'; please make the terminology consistent (e.g., 'field-mode normalisation').","section":"Section 4.4.4 and Fig. 2"},{"comment":"The text contains the typo 'soucers' for 'sources' in the discussion of ionisation sources.","section":"Section 5.2"},{"comment":"The caption lists 'emission line-less retired and line-less retired (ELR and ELR, respectively)'; the second should be 'LLR'.","section":"Figure 11 caption"},{"comment":"The word 'redsifhted' should be 'redshifted' in the description of the counter-outflow.","section":"Section 7.5"},{"comment":"The claim that the multi-component approach is 'tested to be robust against local minima' is not supported by an explicit quantitative test in the text; the outlier-detection procedure is described, but no comparison of recovered chi-squared or parameters against known injected solutions is shown. Please either add such a test or soften the claim.","section":"Section 5.2/Summary item (ii)"},{"comment":"The data availability statement points to the ESO archive but no link or version is given for the custom fitting pipeline; for a methodology paper, providing a code repository or a clear statement of availability would improve reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is a workmanlike and generally honest case study, and the outflow morphology is genuinely interesting. My main reservation is that the inflow claim is built on a residual map whose zero point is a single circular-disc model, and the outflow mass/rate carry unquantified systematics from n_e, C, and metallicity. These are fixable within the scope of the paper by adding robustness tests and error budgets; they are not fatal. The paper could also be strengthened by releasing the pipeline code, since part of its stated contribution is methodological."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a solid, honest case study, and the pipeline description is genuinely useful. The new bits are the spatially resolved emission-line component maps for NGC 613, the residual-velocity evidence for inflow along the bar dust lanes, and the [O III] biconical outflow aligned with the VLA jet. The method assembly – global optimisation on binned spectra, LOF outlier rejection, RBF interpolation to seed spaxel fits, AICc for component counting – is sensible and clearly described; the field-normalisation trick for setting the regularisation parameter across the FoV is a practical detail worth stealing.\n\nThe outflow claim is on solid ground: the biconical [O III] morphology, the two kinematic components, and the alignment with the radio jet agree with earlier work (Davies et al. 2017; Audibert et al. 2019). The quantitative outflow mass and rate are much softer. They rest on an adopted n_e ≈ 112 cm^-3, C = 1, and metallicity calibrations that the authors themselves flag as possibly artifactual in the outflow; the ±0.01 dex error bar is a statistical lower limit. The authors say this plainly, so it is not a hidden flaw, but any use of the number should carry that caveat.\n\nThe inflow interpretation is the softest spot. The residual map is built against a single inclined-disc Kinemetry model fit to the stellar field. NGC 613 has a strong bar with box/peanut structure, a nuclear disc, and possibly an inner bar. In that potential, the stellar field itself is non-circular, and a circular best-fit model will absorb some of the bar's signature, leaving residuals that can look like radial streaming. The dust-lane coincidence and the sign of the residuals match inflow, and the ALMA torque analysis supports inward gas transport, so the claim is plausible. But the paper does not quantify how much of the residual could be produced by bar orbits alone. That is a moderate concern, not a fatal one.\n\nThe main practical gap is that no pipeline code or machine-readable maps are released; the description is detailed enough for a specialist group to re-implement, but it will cost someone weeks.\n\nOverall: worth a serious referee, and I would send it out. It is a careful, reproducible-in-principle case study that will be useful to people working on IFU emission-line fitting and bar-driven gas flows. The inflow conclusion needs a more explicit test against non-circular orbital models before it becomes a strong statement, but that is a revision-level issue, not a desk-reject issue.","headline":"A careful, honest case study with a reusable IFU fitting pipeline: the outflow detection is solid, the inflow claim is plausible but needs a stronger test against non-circular bar orbits.","tokens_in":38444,"tokens_out":2404,"would_cite":true,"duration_ms":22337,"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":"A study of the barred galaxy NGC 613 finds gas streaming inward along the bar's dust lanes and a biconical ionised outflow aligned with the radio jet, with a mass outflow rate near 0.04 solar masses per year.","keywords":["NGC 613","barred galaxies","emission-line fitting","integral field spectroscopy","gas inflows","AGN outflows","star formation history","galaxy evolution"],"falsifier":"Take the observed stellar mass distribution of NGC 613, including the bar, box/peanut, and nuclear disc, and run a hydrodynamical simulation with no net inflow; if it reproduces the same residual velocity pattern along the dust lanes, the inflow claim is falsified. A cheaper test is to recompute the residual map with the inclination and position angle varied within their published uncertainties and check whether the dust-lane inflow signature survives.","tokens_in":37236,"feed_emoji":"🌌","tokens_out":12252,"duration_ms":84997,"temperature":0.7,"pith_summary":"The paper develops a staged analysis pipeline for integral-field spectra—one spectrum per pixel—and uses it on the nearby barred galaxy NGC 613 to separate the many overlapping gas components in its centre. The authors claim that gas is flowing inward along the bar's dust lanes, reaching the nucleus, while a biconical ionised outflow, aligned with the radio jet, carries gas out of the galaxy at about $0.04\\,M_\\odot\\,\\mathrm{yr}^{-1}$. If correct, NGC 613 displays both ends of a gas cycle in one system: bar-driven inflow feeding the active nucleus and AGN feedback returning gas to the host. The methods—field-mode normalisation for regularised star-formation-history fits, global optimisation to avoid local minima, and Akaike-information-criterion model selection for the number of Gaussian components—are presented as reusable solutions to common problems in emission-line analysis.","feed_headline":"Gas streams down NGC 613's bar lanes into a central outflow","feed_subtitle":"MUSE spectra show bar-lane inflow feeding a biconical outflow, at 0.04 solar masses per year.","key_machinery":"The load-bearing object is a multi-component Gaussian emission-line model: lines are split into low- and high-ionisation groups, each fitted with up to three Gaussian components whose number is chosen by the Akaike Information Criterion, a model-selection rule that penalises extra free parameters. The fitting uses differential-evolution global optimisation on spatially binned spectra to escape local minima, then refines the solution pixel by pixel with an outlier-filtered radial-basis-function interpolation providing initial guesses. To turn kinematics into astrophysics, the authors subtract a Kinemetry circular-disc model from the gas velocity and overlay the residuals on HST dust-lane contours and VLA radio contours; those residual maps are the evidence for inflow and outflow.","core_discovery":"The central claim, stated on the paper's own terms, is that the residual velocity field of the ionised gas—after subtracting a best-fitting circular-rotation model built from the stellar kinematics with Kinemetry—shows non-circular motions along both HST-traced bar dust lanes that are directed toward the galaxy centre. The same analysis, applied to a second, broader low-ionisation component and to the high-ionisation [O III] lines, reveals a biconical (two-sided cone) outflow whose blueshifted and redshifted sides bracket the nucleus and align with the VLA radio jet. From the [O III] luminosity, a mean electron density of about $112\\,\\mathrm{cm}^{-3}$, and the O3N2 oxygen-abundance indicator, the authors derive an ionised outflow mass of $\\log_{10} M_{\\rm ion}=4.97\\,M_\\odot$, a bulk timescale of about $2.1\\,\\mathrm{Myr}$, a mass outflow rate of about $0.04\\,M_\\odot\\,\\mathrm{yr}^{-1}$, and a kinetic energy rate of $2.4\\times 10^{40}\\,\\mathrm{erg\\,s^{-1}}$, corresponding to about $1.5\\%$ of the AGN bolometric luminosity.","pith_inferences":["Beyond the paper: applying the same residual-velocity method to a sample of barred versus unbarred galaxies would test whether inflow strength scales with bar strength, which a single object cannot establish.","Beyond the paper: the apparent low metallicity of the outflow gas could be checked with temperature-sensitive auroral lines; if it survives, it would mean the outflow ejects comparatively metal-poor gas rather than processed nuclear gas.","Beyond the paper: the quoted $0.04\\,M_\\odot\\,\\mathrm{yr}^{-1}$ counts only ionised gas, so including molecular and neutral phases would probably raise the true mass-loss rate and could change the feedback-efficiency estimate.","Beyond the paper: a second-epoch observation of the same field could test whether the biconical outflow geometry is steady or flickering on megayear timescales."],"forward_implications":["Bar-driven inflow along dust lanes can deliver gas to the central kiloparsec, providing a direct observational link between bars and nuclear fuel supply.","The biconical outflow is consistent with being AGN-powered, since its kinetic power is about $1.5\\%$ of the AGN bolometric luminosity.","Diagnostics applied to individual Gaussian components reveal ionisation that a single integrated classification would misattribute, so future integral-field studies should separate components before classifying.","The field-normalisation and model-selection choices should make star-formation history and gas kinematics recovery more reliable in other nearby galaxies with complex centres.","The measured outflow is small but operates on a roughly $2\\,\\mathrm{Myr}$ timescale, matching a short feedback episode that can affect the central gas reservoir."],"supporting_citations":[{"why":"Establishes the theory that bars create dust-lane shocks and drive gas inflows, against which the residual velocities are interpreted.","marker":"Athanassoula (1992)"},{"why":"Supplies NGC 613's morphological classification and the disc inclination and orientation used for the circular-velocity model.","marker":"Buta et al. (2015)"},{"why":"Provides the bar's projected semi-major axis, position angle, ellipticity, and boxiness used for the bar contour overlays.","marker":"Kim et al. (2014)"},{"why":"Gives the 'galactic rivers' picture connecting bar dust lanes to inward gas transport.","marker":"Sormani et al. (2023)"},{"why":"Supplies the Kinemetry-based procedure for computing the circular-velocity model and residual maps.","marker":"Kolcu et al. (2023)"},{"why":"Detected the VLA radio jet and elongated [O III] emission that define the outflow axis.","marker":"Hummel et al. (1987)"},{"why":"Provides the [O III]-luminosity-based formula used to estimate the ionised outflow mass.","marker":"Carniani et al. (2015)"},{"why":"Provides the conventions and equations for outflow bulk velocity, timescale, and mass outflow rate.","marker":"Smethurst et al. (2021)"},{"why":"Presents the MUSE observations and data reduction for NGC 613 used throughout the analysis.","marker":"Gadotti et al. (2019)"}],"fun_headline_variants":["Bar dust lanes feed gas into NGC 613's biconical outflow","Gas inflow along bar lanes fuels biconical outflow in NGC 613","NGC 613: infalling gas meets outflow aligned with radio jet","Tracing gas cycle: bar lanes to central outflow in NGC 613","MUSE reveals gas inflow and outflow in barred galaxy NGC 613"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The inflow interpretation rests on the assumption that a single rotating stellar disc tilted at $39^\\circ$ with one adopted position angle fully accounts for the galaxy's gravitational circular motion, so the leftover gas velocities must be real inward streaming rather than artefacts of the bar's more complex potential or a mis-specified disc orientation.","fun_headline_variants_meta":{"raw":{"variants":["Bar dust lanes feed gas into NGC 613's biconical outflow","Gas inflow along bar lanes fuels biconical outflow in NGC 613","NGC 613: infalling gas meets outflow aligned with radio jet","Tracing gas cycle: bar lanes to central outflow in NGC 613","MUSE reveals gas inflow and outflow in barred galaxy NGC 613"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1576,"prompt_tokens":1100,"completion_tokens":476,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":382}},"tokens_in":716,"tokens_out":476,"duration_ms":5134,"temperature":1.0,"reasoning_tokens":382,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:29:24.294656+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the observed stellar mass distribution of NGC 613, including the bar, box/peanut, and nuclear disc, and run a hydrodynamical simulation with no net inflow; if it reproduces the same residual velocity pattern along the dust lanes, the inflow claim is falsified. A cheaper test is to recompute the residual map with the inclination and position angle varied within their published uncertainties and check whether the dust-lane inflow signature survives.","supporting_citations":[],"review_version":1}