{"id":"5e896e41-6d6d-472e-a5df-c730c3dd3134","arxiv_id":"1908.09229","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Molecular gas discs in six low-excitation radio galaxies rotate in an ordered manner but show ubiquitous low-level distortions, with NGC 3557 yielding a black hole mass of 7.1e8 solar masses and NGC 3100 showing evidence of inflowing gas.","lead":"Six nearby radio galaxies observed with ALMA show cold molecular gas around their central black holes mostly rotating in orderly discs, but with small kinematic disturbances in every case. The results suggest the gas is not fully settled and may be flowing inward to feed the active nucleus in at least one galaxy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Residuals against the simple axisymmetric arctangent model in marginally resolved discs may not uniquely imply unsettled gas; the NGC 3557 SMBH-mass error is formal-only.","rationale":"The paper does many things well: it uses public KinMS/MCMC forward modelling, explicitly masks absorption, combines a higher-resolution cube for NGC 3557, and is unusually candid about limitations. The detection of (sub-)kpc molecular discs in LERGs is solid, and ordered rotation is securely established. The load-bearing step is the attribution of residual velocities to physical distortions in objects where the discs are only marginally resolved. The authors themselves say the data prevent strong conclusions, yet the discussion (Sections 5.1 and 5.2) draws exactly those conclusions, so this is not merely a stylistic caution. The reader's weakest-assumption analysis identifies the same issue; my independent reading agrees. Because the paper is already CONDITIONAL and the caveats are largely stated, I would not change the verdict.","tokens_in":29503,"tokens_out":7064,"duration_ms":73622,"concrete_test":"Generate mock ALMA cubes of intrinsically axisymmetric, relaxed discs with the best-fit parameters, resolution, channel width and noise of IC 1531, NGC 7075, IC 4296 and NGC 3557, and run them through the same KinMS/MCMC pipeline used in the paper; measure the residual moment-1 amplitudes. If the mock residuals are comparable to the observed ~30–60 km/s, those residuals are not diagnostic of unsettled gas or of a Keplerian upturn; if they are much smaller, the observed distortions are genuine.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's kinematic analysis (Sections 3–4) starts from a deliberately simple null model—thin, axisymmetric disc, single PA and inclination, arctangent rotation curve (Eq. 1), constant sigma_gas, Gaussian surface brightness—and interprets any residual above the channel width as a physical warp or non-circular motion. For IC 1531, NGC 7075 and IC 4296 the discs are only marginally resolved (Table 1 beams 0.6–0.9 arcsec against disc extents of a few arcsec), so the data are unlikely to distinguish a genuine kinematic distortion from a slightly more flexible model: a radial PA/inclination gradient, a non-arctangent rotation curve, or beam smearing of a thicker or ring-like disc can produce residual patterns of the same amplitude (~30–60 km/s). The paper explicitly concedes this in Section 5.1 ('differentiating between the two through disc modelling alone is not always straightforward') and in the abstract ('preventing us from drawing strong conclusions'). The 'not fully settled' conclusion in Section 5.1 and the NGC 3100 inflow inference in Section 5.2, however, both rely on the residuals being real non-circular motions rather than model artefacts. The NGC 3557 SMBH mass (Section 4.5) sharpens the issue: the Keplerian upturn is one-sided (the model predicts a symmetric upturn; only the positive-velocity side is seen), residuals of about 40 km/s remain in the centre, and the quoted +/-0.02 dex is only the random MCMC error, excluding MGE dust-masking, stellar M/L, inclination, and non-circular-motion systematics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the 3D kinematic modelling of the six molecular gas discs detected in CO(2-1) with ALMA in a volume- and flux-limited sample of eleven nearby low-excitation radio galaxies (Paper I). Using the KinMS forward-modelling code and an MCMC sampler, the authors fit a thin axisymmetric disc model with a Gaussian surface brightness distribution and an arctangent rotation curve (Eq. 1), allowing for warps and an inner hole where needed. They find that the bulk of the gas is in ordered rotation in all six systems, but with low-level residuals that they attribute to warps or non-circular motions. Two sources receive special treatment: NGC 3557, where a central velocity rise is modelled as a Keplerian upturn yielding M_SMBH = (7.10 +/- 0.02) x 10^8 M_sun, and NGC 3100, where a two-armed spiral plus PA and inclination warps are fitted and kinemetric harmonic analysis indicates non-circular motions, interpreted as possible inflows feeding the AGN. The paper also discusses the external origin of the cold gas and estimates Toomre Q parameters.","tokens_in":29936,"tokens_out":8895,"duration_ms":83818,"significance":"If the conclusions hold, the paper provides one of the few resolved molecular-gas kinematic studies of LERGs, supporting a picture in which cold gas is common in their centres and often not fully relaxed, with potential consequences for AGN fuelling. The NGC 3557 SMBH measurement adds a CO-based datum to the M_sigma relation, and the NGC 3100 analysis is a genuinely interesting case study. The authors are careful in several places, especially in acknowledging that the discriminating power of the data is limited; the use of public modelling tools and the explicit description of priors and fitting procedures is also a strength. However, the central \"unsettled gas\" claim and the two headline results rely on interpretations of residuals and of a one-sided kinematic feature that are not fully protected against model degeneracy and unquantified systematics.","major_comments":[{"comment":"The conclusion that the molecular gas discs are \"not fully settled\" in the stellar potential rests on residual moment-one maps obtained with a deliberately simple axisymmetric arctangent model (Eq. 1). For the marginally resolved targets (IC 1531, IC 4296, NGC 7075, and NGC 3557; beams 0.4–0.9 arcsec in Table 1), residual amplitudes of ~30–60 km/s are comparable to what beam smearing, a mild radial PA or inclination gradient, or a non-arctangent rotation curve would plausibly generate. Section 5.1 itself states that \"differentiating between the two through disc modelling alone is not always straightforward.\" The paper should either demonstrate with a synthetic-data or model-comparison test that these residuals cannot be absorbed by a slightly more flexible simple model, or restrict the unsettled-disc claim to the well-resolved cases (NGC 612, NGC 3100) and present the other cases as detections of ordered rotation with unclassified asymmetries.","section":"§5.1; Figures 1–4, Table 1"},{"comment":"The quoted SMBH mass of NGC 3557, log(M_SMBH) = 8.85 +/- 0.02, is not a supported measure of the total uncertainty. The text notes that the Keplerian rise is seen only on the positive-velocity side, that central residuals of ~40 km/s remain, and that \"we cannot definitively claim that a massive dark object is the cause.\" The error quoted in Table 4 is the statistical MCMC uncertainty only; the dependence on the dust-masked MGE model, the stellar M/L, the inclination, and the choice of the stellar potential is not quantified. Because the abstract and Section 6 present this value as a headline result and compare it with the M_sigma relation, the paper should supply a systematic error budget or explicitly downgrade the measurement to an indicative estimate with a conservatively enlarged uncertainty.","section":"§4.5; Table 4; Section 6"},{"comment":"The claim that radial motions in NGC 3100 are likely inflows is one interpretive step beyond what the current modelling measures. The Kinemetry analysis establishes a non-zero k5/k1 ratio, i.e., the presence of non-circular motions, but the sign and radial nature of the velocity component are inferred from the residual pattern and from the dust-absorption near-side determination rather than from a fit that includes an explicit radial velocity term. A model with a free radial velocity component, or a clear statement that the inflow direction is not directly constrained, is needed before the fuelling conclusion can be given the prominence it receives in the abstract and in Section 5.2.","section":"§4.6.1; §5.2"}],"minor_comments":[{"comment":"The note in Table 1 says the NGC 3557 observations are \"combined Cycle 3 and Cycle 6 data,\" while Section 2.1 states the archival data were taken during Cycle 4; please correct the inconsistency.","section":"Table 1; Section 2.1"},{"comment":"The text says that \"in two cases (IC 1531 and NGC 612) Q is below unity,\" but Table 6 lists IC 1531 only as a lower limit (>0.2); a lower limit does not establish Q < 1. Likewise, IC 4296 is listed as an upper limit (<29), so the statement that \"Q is well above unity\" should be qualified.","section":"§5.4; Table 6"},{"comment":"The captions describe panels as \"mock\" mean velocity maps; this should read \"model\" for consistency with Figures 1–4.","section":"Figures 5 and 7"},{"comment":"There are several typographical errors (e.g., \"thiss hypothesis\" in Section 5.1 and missing periods after \"e.g\" in Sections 5.1 and 5.2); a careful proofread would improve the presentation.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"This is a solid but not definitive paper; the referee report reflects the view that the authors' own caveats should be moved into the abstract and headline results, and that the systematic error budget for the SMBH mass should be quantified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a solid, workmanlike forward-modelling study of six CO(2-1) discs in low-excitation radio galaxies. It gives the first kinematic models for these specific discs, a CO-based SMBH mass for NGC 3557, and a reasonable, clearly argued case for cold-gas inflow in NGC 3100. It does not invent new methods; it applies the public KinMS tool to a new, well-defined sample. The authors are also unusually candid: the abstract itself says most discs are only marginally resolved and strong conclusions are not possible. That honesty carries through the text.\n\nWhat the paper does well: the per-source modelling is standard, transparent, and reproducible in principle—the data are public, the code is public, and the residual maps and parameter tables are clear. The NGC 3557 SMBH estimate uses an MGE stellar mass model, treats the black hole as a point mass, and checks against the M-sigma relation. The NGC 3100 analysis combines a warp model with kinemetry and makes a specific, testable claim that the radial motions are inflows connected to spiral streaming. The Toomre stability section is properly caveated about beam smearing and unresolved dispersion.\n\nSoft spots, in proportion: the main one is the interpretation of residuals from a deliberately simple axisymmetric arctangent model. In four of the six galaxies the discs are only marginally resolved. A radial PA/inclination gradient, a slightly different rotation curve, or beam smearing of a thicker or ring-like disc could plausibly create similar residual patterns. The authors acknowledge this in Section 5.1, but the summary statement that the gas is “not fully settled” leans on those residuals. For NGC 3557, the SMBH mass is driven by a one-sided central velocity upturn with ~40 km/s residuals right at the centre; the quoted ±0.02 dex is purely the MCMC random error and does not include systematics from dust masking, M/L, inclination, or non-circular motion. So I would treat that mass as provisional. The inflow direction in NGC 3100 depends on the dust-lane near-side assumption; it is plausible, not airtight.\n\nThese are caveats, not fatal flaws. The paper’s central claims are exactly as hedged as they should be, and the observational material is new and useful. This is a legitimate incremental contribution to a sparse literature.\n\nRecommendation: send it to a competent referee. It deserves peer review, and with the stated caveats it should be publishable. I would cite it if I were working on molecular gas in early-type galaxies or CO-based black hole masses. It is a reasonable reading-group choice for discussing model degeneracy in marginally resolved kinematic data.","headline":"A careful, honest kinematics paper that adds six LERGs to a small sample; the SMBH mass and inflow claims are plausible but rest on residuals from a simple model in marginally resolved discs.","tokens_in":30557,"tokens_out":2108,"would_cite":true,"duration_ms":25314,"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":"This paper claims that molecular gas discs in low-excitation radio galaxies rotate but are not settled, and in at least one object the gas is inflowing and can feed the black hole.","keywords":["molecular gas kinematics","low-excitation radio galaxies","AGN fuelling","CO(2-1) ALMA observations","supermassive black hole mass","galactic warps","non-circular motions","Toomre stability"],"falsifier":"Observe the same six galaxies at sub-100 pc resolution, or with a higher-density tracer such as CO(3-2), and ask whether the residuals between the data and the best simple rotating-disc model persist; if they vanish or shrink to the noise level, the conclusion that the gas is unsettled loses its support, and for NGC 3557 a map resolving the sphere of influence on both sides of the nucleus would show whether the velocity upturn is a Keplerian rise or a one-sided non-circular flow.","tokens_in":29373,"feed_emoji":"🕳️","tokens_out":5592,"duration_ms":52941,"temperature":0.7,"pith_summary":"This paper tries to establish what the cold molecular gas around nearby low-excitation radio galaxies is doing: is it feeding the central black hole, being pushed around by jets, or sitting in stable orbits? Using ALMA observations of six galaxies detected in carbon monoxide, the authors model each CO disc in three dimensions. They find that the bulk of the gas rotates in an ordered way, but every disc shows low-level kinematic distortions. They argue these distortions mean the gas is not fully settled into the host galaxy's gravitational potential, and in at least one case, NGC 3100, the radial motions are likely inflows that contribute to AGN fuelling.","feed_headline":"Unsettled gas discs could be feeding radio galaxy black holes","feed_subtitle":"ALMA maps of six low-excitation radio galaxies show rotating but distorted CO discs, including one galaxy with inflowing gas.","key_machinery":"The central tool is the KinMS forward-modelling code coupled to an MCMC fitter: it generates mock ALMA data cubes from a parametric gas distribution and velocity law, so observed cubes can be compared with models that include beam smearing, disc thickness, and velocity dispersion. The default velocity law is an arctangent rotation curve $v(R)=(2v_{\\rm flat}/\\pi)\\arctan(R/r_{\\rm turn})$. Deviations between data and model are then interpreted through added ingredients: position-angle and inclination warps, a central surface-brightness hole or ring, and, for NGC 3100, a two-armed logarithmic spiral brightness pattern. A harmonic decomposition of the velocity field (kinemetry) quantifies non-circular motions, and for NGC 3557 a multi-Gaussian expansion of the stellar light supplies the stellar mass model against which the black hole mass is fitted.","core_discovery":"The paper's central claim is that sub-kiloparsec molecular gas discs in low-excitation radio galaxies are common, mostly rotation-dominated, but not relaxed: deviations from a simple axisymmetric rotating disc appear in all six objects studied. For most the deviations are small and the discs are only marginally resolved, so the authors are cautious. The special cases carry the argument: NGC 3557 shows a central velocity upturn fitted as a Keplerian rise, giving a black hole mass of $(7.10\\pm0.02)\\times10^{8}\\,M_\\odot$; NGC 3100 requires a warped, two-armed spiral disc whose residual velocity field has non-negligible radial motions, most plausibly inflows toward the nucleus. On this basis the paper concludes that cold gas can contribute to fuelling the AGN in these objects, and that such disturbed discs distinguish radio-loud from radio-quiet early-type galaxies.","pith_inferences":["A testable extension the paper leaves implicit: compare CO residual maps with jet orientation and stellar kinematics in a larger LERG sample to see whether the asymmetries systematically align with jet axes, which would favour jet-gas interactions over unrelated warps.","If the central velocity rise in NGC 3557 were instead produced by a bar or a nuclear non-circular flow, the derived black hole mass would be an upper limit; a resolved map of the sphere of influence on both sides of the nucleus would settle this.","The relaxation timescales estimated for NGC 612 and NGC 3100 imply that externally acquired gas can survive several $10^8$ yr before settling, so young unsettled discs should be most common in galaxies with recent interaction signatures."],"forward_implications":["If the paper is right, low-level kinematic disturbances in CO discs should be regarded as the normal state of low-excitation radio galaxies, not as rare anomalies.","The NGC 3557 black hole mass joins the set of directly measured supermassive black hole masses in early-type galaxies and can be used to test the $M_{\\rm SMBH}$--$\\sigma_*$ relation at the high-mass end.","The inferred inflow in NGC 3100 supports externally accreted or chaotically cooling gas as a viable fuelling path for kinetic-mode AGN, not only hot halo accretion.","Disturbed, unrelaxed discs should be more common among radio-loud than radio-quiet early-type galaxies, a difference that can be tested with larger samples.","Most of the modelled discs are gravitationally stable (Toomre $Q>1$), so fuelling likely proceeds through slow, non-collapse channels rather than rapid fragmentation."],"supporting_citations":[{"why":"Paper I of this series; supplies the ALMA CO(2-1) detections, moment maps, and gas masses that this kinematic study models.","marker":"Ruffa et al. 2019"},{"why":"Introduces the KinMS forward-modelling tool used to fit the CO data cubes.","marker":"Davis et al. 2013"},{"why":"Establishes the method of deriving supermassive black hole masses from Keplerian upturns in molecular gas kinematics.","marker":"Davis et al. 2017"},{"why":"Provides the same Keplerian-modelling approach applied to CO discs in early-type galaxies.","marker":"Onishi et al. 2017"},{"why":"Supplies the comparison sample of mostly settled CO discs in radio-quiet early-type galaxies.","marker":"Alatalo et al. 2013"},{"why":"High-resolution ALMA study of radio-quiet early-type galaxies with settled CO discs; the contrast that makes LERGs appear disturbed.","marker":"Boizelle et al. 2017"},{"why":"Gives the kinemetry harmonic expansion used to quantify non-circular motions in NGC 3100.","marker":"Krajnović et al. 2006"},{"why":"Hydrodynamic simulations of externally accreted gas that motivate the settling-sequence interpretation and relaxation timescale estimates.","marker":"van de Voort et al. 2015"}],"fun_headline_variants":["Distorted gas discs hint at feeding paths for radio galaxies","Black hole mass weighed from gas motions in a radio galaxy","Inflowing gas in a warped disc feeds a nearby radio AGN","Distorted CO discs reveal gas feeding low-excitation radio galaxies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The cold-gas discs in most of the sample are only marginally resolved, and the central claim that the residual velocity features are genuine warps, non-circular motions, or inflows assumes those features are not artifacts of fitting a simple thin axisymmetric rotating-disc model to barely resolved data, including effects like beam smearing.","fun_headline_variants_meta":{"raw":{"variants":["Distorted gas discs hint at feeding paths for radio galaxies","Black hole mass weighed from gas motions in a radio galaxy","Inflowing gas in a warped disc feeds a nearby radio AGN","Distorted CO discs reveal gas feeding low-excitation radio galaxies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000615,"raw_usage":{"total_tokens":2904,"prompt_tokens":1041,"completion_tokens":1863,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":657,"completion_tokens_details":{"reasoning_tokens":1791}},"tokens_in":657,"tokens_out":1863,"duration_ms":14435,"temperature":1.0,"reasoning_tokens":1791,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:17:24.434031+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the same six galaxies at sub-100 pc resolution, or with a higher-density tracer such as CO(3-2), and ask whether the residuals between the data and the best simple rotating-disc model persist; if they vanish or shrink to the noise level, the conclusion that the gas is unsettled loses its support, and for NGC 3557 a map resolving the sphere of influence on both sides of the nucleus would show whether the velocity upturn is a Keplerian rise or a one-sided non-circular flow.","supporting_citations":[{"cited_title":"D., Barth A","cited_arxiv_id":null,"evidence_quote":"High-resolution ALMA study of radio-quiet early-type galaxies with settled CO discs; the contrast that makes LERGs appear disturbed."}],"review_version":1}