{"id":"239604da-5c20-477f-9b60-55e948ee815e","arxiv_id":"2507.10662","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"First sphere-of-influence-resolving CO dynamical measurement gives the NGC 1574 black hole a mass of (6.2 +/- 1.2) x 10^7 solar masses.","lead":"Astronomers used ALMA's sharpest millimeter-wave vision to map gas swirling near the center of galaxy NGC 1574, catching the first clear gravitational signature of its supermassive black hole. They measure the black hole at about 62 million solar masses, lower than a previous estimate and a sign that lower-resolution measurements can be biased.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unmodeled non-circular motions (±55 km/s residuals) can bias M_BH by tens of percent via M_BH∝v², comparable to the quoted 20% uncertainty; the paper only argues by amplitude comparison, not by fitting a non-circular model.","rationale":"Agree with reader's weakest_assumption. The strongest claim—accurate M_BH from resolving the SoI—requires that the circular-orbit model is a good description of the kinematics within the SMBH-dominated region. The residual map shows that it is not perfect. Because M_BH is quadratic in the deprojected velocity, a 15% systematic velocity error translates to ~30% mass error, larger than the reported ~20% statistical uncertainty. The paper's argument based on amplitude comparison is insufficient; a proper test requires fitting non-circular terms. Thus I recommend CONDITIONAL acceptance: the measurement is plausible and well-executed, but the central accuracy claim should be contingent on demonstrating that M_BH is robust to non-circular motions. The concern is not an accusation of error; it is a request for a quantitative systematic check that is feasible with the existing code and data. The secondary SoI-definition issue is noted but does not affect the mass value; it only tempers the 'first to resolve' novelty statement.","tokens_in":21069,"tokens_out":6460,"duration_ms":77166,"concrete_test":"Re-run the GAStimator MCMC fit with an added parameter for a radial velocity component (e.g., v_R = f × v_circ, or a power-law radial flow) in the KinMS model, holding all other settings identical. If the marginalized posterior of M_BH shifts by more than 1σ (i.e., central value moves by >1.2×10^7 M_sun) when v_R is allowed, the circular-orbit assumption is not innocuous and the accuracy claim is overstated. Alternatively (or in addition), mask the central 0.15 arcsec region where residuals are largest and refit; a significant M_BH change would confirm the bias. The no-SMBH comparison should also be re-run with all other parameters refitted to verify the 'unambiguous' detection claim.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The dynamical model assumes purely circular orbits (Section 3.2). Section 4.1 reports residual spiral-arm-like structures of up to ±55 km/s (deprojected) in the first-moment map, which are not reproduced by the model and are attributed to bar-induced non-circular motions. The inferred SMBH mass obeys M_BH ∝ [v_circ(R)/sin i]², so a coherent non-circular velocity component of ~15% of the local circular speed can alter the mass estimate by ~30%, exceeding the listed 1σ uncertainty (6.2±1.2×10^7 M_sun). The paper's rebuttal compares the residual amplitude to the maximum circular velocity (~380 km/s) rather than to the circular velocity at the radii that dominate the M_BH likelihood, and does not test a model that includes radial or streaming motions. Because the central claim emphasizes improved accuracy from resolving the SoI, an unquantified systematic of this size is load-bearing. A related but secondary concern is that the 'first to spatially resolve the SoI' claim relies on the R_eq definition; the beam (0.074 arcsec) is larger than the traditional R_SoI=GM/σ²≈0.059 arcsec, though the data do trace the Keplerian rise at R_min≈0.043 arcsec.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents new ALMA 12CO(2-1) observations of the lenticular galaxy NGC 1574 with a synthesised beam of 0.078\" x 0.070\" (≈7.5 x 6.7 pc), roughly a factor of two higher resolution than the previous intermediate-resolution data. The observations reveal a prominent central velocity rise in the position-velocity diagram, which the authors interpret as the first spatial resolution of the SMBH's sphere of influence in this galaxy. Using the KinMS forward-modelling code and GAStimator MCMC sampler, they model the full data cube and infer a SMBH mass of (6.2 ± 1.2) x 10^7 M_sun, consistent with but slightly smaller than the earlier value of (1.0 ± 0.2) x 10^8 M_sun from Ruffa et al. (2023). They also characterise a position-angle warp, an M/L ratio, and place the galaxy on the M_BH-σ_e and M_BH-M_* scaling relations.","tokens_in":21402,"tokens_out":4354,"duration_ms":53739,"significance":"If the measurement is robust, the paper provides a valuable high-spatial-resolution molecular gas dynamical mass measurement that directly tests the accuracy of a previous, marginally resolved measurement and contributes to the WISDOM sample. The strengths of the work include the forward modelling of the complete data cube (not just moment maps), the explicit comparison with a no-SMBH model that fails to reproduce the central velocity rise, the use of publicly available and widely used tools (KinMS, GAStimator, JAM), and a detailed discussion of statistical and systematic uncertainties. The main weakness is the treatment of non-circular motions, which are detected as ±55 km/s residual spiral-arm features but are only dismissed by an amplitude comparison rather than by a quantitative bias estimate. A secondary weakness is the definition of the sphere of influence used to support the 'first to resolve' claim, which depends on R_eq and R_min rather than the traditional R_SoI that is slightly smaller than the beam.","major_comments":[{"comment":"The manuscript argues that the ±55 km/s deprojected residuals are negligible because they are smaller than the maximum circular velocity of ~380 km/s within the SoI. This comparison is not appropriate for assessing the mass bias: the SMBH mass is constrained primarily by the innermost Keplerian region, where the circular velocities are considerably lower than 380 km/s, and M_BH scales as the square of the deprojected velocity. A coherent non-circular component of ~15% of the local circular speed at the radii that dominate the likelihood would shift M_BH by ~30%, which exceeds the quoted 1σ statistical uncertainty of ±1.2 x 10^7 M_sun. To support the claim of improved accuracy, the authors should test the sensitivity of M_BH to non-circular motions, for example by (i) fitting a model that includes radial or streaming motions, (ii) masking the residual spiral-arm regions and refitting, or (iii) explicitly adding a systematic uncertainty term derived from the residual velocity field. As it stands, the potential bias is unquantified and comparable to the reported error budget.","section":"§4.2, Eq. (3)"},{"comment":"The claim that these are 'the first to spatially resolve the SMBH's sphere of influence' rests on the choice of R_eq (0.12\") and R_min (0.043\") rather than the traditional R_SoI = G M_BH / σ_e^2 ≈ 0.059\". Since the synthesised beam FWHM of 0.074\" is actually slightly larger than R_SoI, the statement is definition-dependent. The authors do motivate R_eq and R_min convincingly, and the no-SMBH model failure provides strong internal support, but the paper should explicitly acknowledge that under the classical definition the beam does not resolve the SoI, and explain why R_eq and R_min are the more meaningful metrics for this dataset. This is important because the 'first to resolve' wording appears in the abstract and is used to argue that the new measurement is more accurate than the previous one.","section":"§4.2, Eq. (3)"}],"minor_comments":[{"comment":"The cleaning-depth flux correction is described thoroughly, but the resulting corrected total flux (10.1 Jy km/s) is substantially lower than the uncorrected intermediate-resolution value (17.9 Jy km/s); it would help the reader if the Jorsater & van Moorsel correction were sketched in a sentence or a reference to the exact equation, rather than only cited.","section":"§2.3.3"},{"comment":"The statement that the inner PA differs from the Ruffa et al. (2023) value by 8σ would be more informative if the quoted uncertainty of the previous measurement were given explicitly in the text (the abstract of that paper already does so, but the value is not repeated here).","section":"§3.4"},{"comment":"There are minor typographical issues, including missing spaces ('the12CO(2-1)'), 'mean line-of-light velocities' instead of 'line-of-sight velocities' in §2.3.2, and inconsistent use of 'circumnuclear' vs 'circum-nuclear' in the introduction. These should be corrected in a final proof.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the analysis is broadly sound, but the non-circular motion systematic is a load-bearing issue that needs to be addressed before the accuracy-improvement claim can be accepted. The paper's own admission of ±55 km/s residuals, together with the M_BH ∝ v² scaling, makes this a real concern rather than a mere stylistic quibble. I believe the authors can address it within the scope of a revision, either by adding a non-circular model or by a careful masking test. The 'first to resolve the SoI' claim is also somewhat definition-dependent; given the prominence of this claim in the abstract, it should be softened or rigorously justified. The self-citation pattern is not problematic for a series paper, but the over-reliance on the same group's previous analysis of the same galaxy is acceptable here because the new data are genuinely independent at higher resolution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: this paper delivers the first spatially resolved SMBH mass for NGC 1574 and the Keplerian detection is convincing, but the unmodeled non-circular motions are a bigger systematic than the paper concedes and could shift the mass by tens of percent.\n\nThe new ALMA CO(2-1) observations (0.078\" x 0.070\") clearly resolve the central velocity rise. The forward modelling with KinMS is standard but well executed; the no-SMBH model fails, the new model reproduces the PVD better than the previous best fit, and the derived mass (6.2 +- 1.2 x 10^7 Msun) is consistent with but slightly lower than the earlier unresolved measurement. The larger central PA warp is a genuine improvement. The paper is honest about limitations and the uncertainty budget, and the data products are publicly available.\n\nThe main soft spot is the treatment of non-circular motions. Residuals show spiral-like features up to +-55 km/s after deprojection, while the model assumes pure circular orbits. Since M_BH scales as v^2, a coherent 14% error in the circular speed would shift the mass by ~30%, larger than the quoted 1-sigma error. The paper's rebuttal compares the residual amplitude to the maximum circular velocity (~380 km/s) in the centre, but that does not rule out a coherent systematic component at the radii that dominate the fit. A more convincing test would be to fit a model with radial or streaming motions, or at least to quantify the effect of the residuals on M_BH. This is not necessarily fatal: the Keplerian rise is prominent and the no-SMBH model is strongly rejected. But the claim of improved accuracy should be tempered until the streaming-motion bias is addressed.\n\nThe 'first to resolve the SoI' claim leans on the R_eq definition; by the classic R_SoI = GM/sigma^2 definition the beam is slightly larger, though the data do trace emission down to R_min < R_SoI. That is a semantic point and the paper discusses it openly.\n\nOverall this is a solid measurement in an established series. The reader's ACCEPT verdict is reasonable, but the non-circular motion systematics deserve a closer look. I would send this to peer review and ask the authors to explicitly test for streaming-motion bias. The paper will be of value to anyone working on molecular gas dynamical SMBH masses and the M_BH scaling relations.","headline":"First SoI-resolving SMBH mass for NGC 1574 is convincing, but unmodeled non-circular motions could bias the result by tens of percent—worth peer review with a requested sensitivity test.","tokens_in":21985,"tokens_out":3758,"would_cite":true,"duration_ms":40633,"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 reports the first NGC 1574 black hole mass measurement from CO gas kinematic data that spatially resolve the black hole's sphere of influence, yielding $M_{\\rm BH}=(6.2\\pm1.2)\\times10^7\\,M_\\odot$ and indicating that the earlier…","keywords":["supermassive black hole mass","molecular gas dynamics","ALMA CO(2-1) observations","sphere of influence","Keplerian rotation","NGC 1574","position-angle warp","galaxy scaling relations"],"falsifier":"Fit the same data cube with a dynamical model that adds radial or other non-circular motions matching the observed $\\pm55$ km/s spiral-arm residuals: if the best-fitting $M_{\\rm BH}$ shifts by more than $1\\sigma$ from $6.2\\pm1.2\\times10^7\\,M_\\odot$, the circular-orbit assumption is falsified. Alternatively, a deeper observation that resolves radii below $R_{\\rm min}\\approx4$ pc should show the projected velocity continuing to rise roughly as $v\\propto r^{-1/2}$ for the lower mass; if the rise instead follows the curve expected for $1.0\\times10^8\\,M_\\odot$, the resolved measurement is wrong.","tokens_in":20880,"feed_emoji":"🕳️","tokens_out":11718,"duration_ms":124473,"temperature":0.7,"pith_summary":"The paper sets out to show that a black hole mass measured from molecular gas kinematics becomes trustworthy only when the observations spatially resolve the black hole's sphere of influence. Using new high-resolution ALMA observations of carbon monoxide in the lenticular galaxy NGC 1574, it reports the first spatially resolved view of that sphere of influence and a clear Keplerian rise in the rotation curve toward the nucleus. Forward modelling of the full data cube gives a black hole mass of $(6.2\\pm1.2)\\times10^7\\,M_\\odot$, slightly smaller than but consistent with the earlier intermediate-resolution measurement. The authors argue that the value of the new observation is accuracy rather than precision: resolving the sphere of influence breaks the degeneracy between black hole mass and stellar mass-to-light ratio, and reveals a central position-angle warp that previous data could not trace. If correct, the result implies that molecular-gas black hole masses from data that only marginally resolve the sphere of influence can be systematically overestimated.","feed_headline":"Resolved black hole sphere sets NGC 1574 mass at 62 million Suns","feed_subtitle":"Resolving the black hole's reach gives a mass nearly 40 percent below the previous estimate.","key_machinery":"The central machinery is the forward modelling of the entire CO data cube rather than just its moment maps. A model galaxy is built from an exponential molecular gas disc and a stellar mass distribution derived from a multi-Gaussian expansion of the Hubble image, deprojected at the fitted inclination and scaled by a constant mass-to-light ratio, with a point mass at the centre representing the black hole. From this combined mass model a circular velocity curve is computed, a constant gas velocity dispersion is added, a position-angle warp that varies linearly with radius is applied, and the resulting velocity field is projected, convolved with the synthesised beam, and spectrally binned to produce a synthetic cube that is compared pixel-by-pixel with the data in a Markov chain Monte Carlo fit. The key quantity carrying the argument is the resolved Keplerian rise in the position-velocity diagram: the innermost detected gas radius is $R_{\\rm min}\\approx4.1$ pc, about a third of the equality radius $R_{\\rm eq}\\approx11$ pc where stellar and black hole masses are equal, so the central velocities are dominated by the black hole and are nearly independent of the stellar mass-to-light ratio.","core_discovery":"On its own terms, the discovery is that the previously unresolved central few parsecs of NGC 1574 contain a Keplerian velocity rise that unambiguously detects the supermassive black hole's gravity, and that modelling this rise gives $M_{\\rm BH}=(6.2\\pm1.2)\\times10^7\\,M_\\odot$. The observations reach a synthesised beam of $0\\farcs078\\times0\\farcs070$ (about $7.5\\times6.7$ pc), and while this beam is larger than the traditional $GM_{\\rm BH}/\\sigma_e^2$ sphere-of-influence radius, it resolves the physically meaningful equality radius $R_{\\rm eq}\\approx0\\farcs12$ (11 pc) at which the enclosed stellar mass equals the black hole mass; the innermost detected gas orbits at $R_{\\rm min}\\approx4.1$ pc, well inside the region where the black hole dominates. The best-fitting model also finds that the isovelocity twist is a position-angle warp that steepens toward the centre, with an inner position angle of $326^\\circ$ rather than the previously inferred $342^\\circ$, and the velocity residuals contain spiral-arm-like non-circular motions of up to $\\pm55$ km s$^{-1}$. The new mass is $1.7\\sigma$ smaller than the previous value of $(1.0\\pm0.2)\\times10^8\\,M_\\odot$, so the unresolved data appear to have been slightly high.","pith_inferences":["If marginally resolved molecular-gas masses are often biased high, part of the observed scatter in the $M_{\\rm BH}$--$\\sigma_e$ relation may be measurement scatter; a systematic comparison of resolved and unresolved measurements of the same galaxies would test this.","The spiral-arm-like residuals suggest bar-driven radial gas flows down to the circumnuclear disc; modelling these flows explicitly, for example with a barred potential or radial-velocity terms, could both correct the residuals and quantify gas inflow toward the black hole.","Because the galaxy is nearly face-on, inclination dominates the statistical error; an independent stellar-dynamical measurement of the same black hole would provide a cross-check of both $M_{\\rm BH}$ and the velocity dispersion.","A future observation with an even smaller beam, about 3 pc or less, could test whether the central Keplerian rise continues along the $6.2\\times10^7\\,M_\\odot$ curve or turns over, revealing additional central mass or a breakdown of the linear warp model."],"forward_implications":["NGC 1574's black hole mass is $(6.2\\pm1.2)\\times10^7\\,M_\\odot$, about $1.7\\sigma$ below the previous unresolved measurement, so the older data overestimated the mass slightly.","Resolving the sphere of influence breaks the usual degeneracy between black hole mass and stellar mass-to-light ratio, so the black hole mass is constrained by the inner Keplerian rise rather than by the stellar light model.","The position-angle warp in the central $0\\farcs2$ region is larger than previously inferred, with the inner position angle at $326^\\circ$ rather than $342^\\circ$; a two-component warp model does not improve the fit, so a single linear warp describes the disc.","NGC 1574 sits slightly below the $1\\sigma$ scatter of the $M_{\\rm BH}$--$\\sigma_e$ relation but within the $3\\sigma$ scatter, and within $1\\sigma$ of the $M_{\\rm BH}$--$M_*$ relation, suggesting the offset is more likely due to an overestimated velocity dispersion.","The measurement implies that molecular-gas dynamical masses from data that only marginally resolve the sphere of influence can be systematically inaccurate, reinforcing the need for high-resolution observations in black hole scaling-relation studies."],"supporting_citations":[{"why":"Provides the previous intermediate-resolution CO dynamical mass measurement and the multi-Gaussian expansion stellar light model that this paper improves upon.","marker":"Ruffa et al. 2023"},{"why":"Supplies the forward-modelling code used to construct simulated data cubes from an input gas distribution and circular velocity field.","marker":"Davis et al. 2013a"},{"why":"Provides the routine used to compute the circular velocity curve from the combined stellar and black hole mass distribution.","marker":"Cappellari 2008"},{"why":"Provides the multi-Gaussian expansion fitting method used to deproject the Hubble image into the stellar mass model.","marker":"Cappellari 2002"},{"why":"Supplies the stellar velocity dispersion of 216 km/s used to estimate the sphere of influence and to place the galaxy on the $M_{\\rm BH}$--$\\sigma_e$ relation.","marker":"Bernardi et al. 2002"},{"why":"Supplies the surface-brightness-fluctuation distance of 19.9 Mpc that sets the physical scale of the beam, radii, and inferred masses.","marker":"Tonry et al. 2001"},{"why":"Provides the $M_{\\rm BH}$--$\\sigma_e$ and $M_{\\rm BH}$--$M_*$ scaling relations against which the new black hole mass is compared.","marker":"van den Bosch 2016"},{"why":"Source of the chi-squared rescaling used to inflate uncertainties and produce the quoted $1\\sigma$ confidence intervals.","marker":"van den Bosch & van de Ven 2009"},{"why":"Companion high-resolution black hole measurement whose methodology for estimating the innermost detected kinematic tracer radius is adopted here.","marker":"Zhang et al. 2025"}],"fun_headline_variants":["ALMA resolves black hole sphere in NGC 1574: mass 62 million Suns","Keplerian rise resolves NGC 1574 black hole mass: 62 million Suns","NGC 1574 black hole mass cut to 62 million Suns by resolved sphere","ALMA sees Keplerian rise, sets NGC 1574 black hole at 62 million Suns","Resolving black hole sphere gives NGC 1574 mass 62 million Suns"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The gas particles are assumed to move on circular orbits when the mass model is converted into kinematics, and if the mild non-circular motions seen as $\\pm55$ km/s spiral-arm residuals are stronger near the centre than the models assume, the black hole mass could be biased.","fun_headline_variants_meta":{"raw":{"variants":["ALMA resolves black hole sphere in NGC 1574: mass 62 million Suns","Keplerian rise resolves NGC 1574 black hole mass: 62 million Suns","NGC 1574 black hole mass cut to 62 million Suns by resolved sphere","ALMA sees Keplerian rise, sets NGC 1574 black hole at 62 million Suns","Resolving black hole sphere gives NGC 1574 mass 62 million Suns"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00083,"raw_usage":{"total_tokens":3737,"prompt_tokens":1172,"completion_tokens":2565,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":788,"completion_tokens_details":{"reasoning_tokens":2460}},"tokens_in":788,"tokens_out":2565,"duration_ms":20213,"temperature":1.0,"reasoning_tokens":2460,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:27:54.596826+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same data cube with a dynamical model that adds radial or other non-circular motions matching the observed $\\pm55$ km/s spiral-arm residuals: if the best-fitting $M_{\\rm BH}$ shifts by more than $1\\sigma$ from $6.2\\pm1.2\\times10^7\\,M_\\odot$, the circular-orbit assumption is falsified. Alternatively, a deeper observation that resolves radii below $R_{\\rm min}\\approx4$ pc should show the projected velocity continuing to rise roughly as $v\\propto r^{-1/2}$ for the lower mass; if the rise instead follows the curve expected for $1.0\\times10^8\\,M_\\odot$, the resolved measurement is wrong.","supporting_citations":[{"cited_title":"V., da Costa L","cited_arxiv_id":null,"evidence_quote":"Supplies the stellar velocity dispersion of 216 km/s used to estimate the sphere of influence and to place the galaxy on the $M_{\\rm BH}$--$\\sigma_e$ relation."}],"review_version":1}