{"id":"499e9d64-ae7d-48c1-8d8b-1cac038655f1","arxiv_id":"2506.13093","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"NGC 4261's twin jets are collimated and accelerated in the same sub-parsec region, with a maximum Lorentz factor of about 2.6, evidence for a compact acceleration and collimation zone.","lead":"Astronomers used archival Very Long Baseline Array observations to map the twin jets of the nearby galaxy NGC 4261 and found that the jets are narrowed and sped up in the same sub-parsec region. This provides a rare, close-up look at how a low-power black hole jet is collimated and accelerated.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The kinematic part of the ACZ claim hinges on interpreting the jet-to-counterjet brightness ratio as purely Doppler, but residual free-free absorption in the counterjet could inflate R_B and mimic the inferred acceleration; a quantitative multifrequency R_B test would settle this.","rationale":"The collimation side of the paper is robust: the parabolic-to-conical transition is confirmed by three independent methods and prior work (N18, Y23), and the break locations are consistent across approaches. The kinematic side, however, is the fragile leg. The velocity profile is derived solely from R_B and α via Eq. (2), and because α is nearly constant, βint(r) is essentially the R_B profile rescaled by 1/cos(68°). The R_B peak at ~4–5×10³ R_s sets Γmax≈2.6, and the subsequent dip/rise becomes the 'deceleration/reacceleration' that the paper aligns with the collimation break. This makes the ACZ claim directly dependent on R_B being purely Doppler. The paper itself documents FFA in the counterjet base and excludes only the region inside the minor-axis beam; it also discards the 2.3 GHz 2002 data because they are not optically thin until >10 mas. These exclusions are post-hoc and are not replaced with a quantitative opacity model. The assertion that FFA/SSA are negligible in the measured region rests only on the qualitative similarity of R_B trends at 5 and 8.4 GHz. A quantitative check of R_B frequency dependence at fixed radii would either validate or falsify the Doppler interpretation. If R_B proves frequency-dependent, the inferred acceleration, deceleration, and Γmax are biased, and the co-spatial ACZ conclusion weakens to a collimation-only result. The reader's CONDITIONAL verdict is appropriate; my concern sharpens the condition but does not move the verdict.","tokens_in":34966,"tokens_out":7060,"duration_ms":78823,"concrete_test":"Compute R_B from the quasi-simultaneous 1.4/2.3/5.0/8.4 GHz maps (BS094b/d) at fixed deprojected radii in the acceleration zone (r ≈ 3, 4, 5, 6, 8 ×10³ R_s), convolving to a common beam and aligning on the core-shift-corrected SMBH position. Under pure Doppler boosting, R_B should be independent of frequency for constant α; residual FFA would make R_B rise toward lower frequency (τ ∝ ν^-2.1). If R_B(1.4 GHz) − R_B(8.4 GHz) exceeds the propagated ~10% intensity uncertainties, FFA contaminates the ratio and Γmax is an overestimate; if R_B is flat within ~20%, the Doppler interpretation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim — a co-spatial ACZ with Γmax≈2.6 — rests on the velocity field derived from Eq. (2) in §3.2.2, which converts the jet-to-counterjet brightness ratio R_B into βint assuming the asymmetry is purely relativistic Doppler boosting, with a single viewing angle 68°±4° and negligible FFA/SSA. Because the spectral index is nearly constant (§3.2.2), the resulting βint(r) is essentially a rescaling of R_B(r). The R_B profile rises to 8–10 at ~4–5 mas (~4–5×10³ R_s), then dips and rises again before declining; this peak is what sets Γmax≈2.6. However, the counterjet is known to suffer free-free absorption near the base (§3.1, Fig. 1, Appendix D), and the authors exclude only data inside the minor-axis beam and the 2.3 GHz 2002 data because they are not optically thin until >10 mas. If foreground ionized gas absorbs the counterjet preferentially at the distances where R_B is measured, the inferred speeds are overestimated, and the peak/dip/rise pattern could be an opacity artifact rather than real acceleration/deceleration. The paper's claim that FFA/SSA are 'negligible' in the measured region is supported only by the qualitative frequency-independence of R_B in Fig. 5, not by a quantitative opacity estimate. This is the load-bearing weak point: the collimation break survives even if the velocity field is wrong, but the co-spatial ACZ conclusion does not.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes archival VLBA observations of the nearby LLAGN NGC 4261 at 1.4–43 GHz to investigate the jet collimation profile and velocity field on sub-parsec scales. The authors combine deconvolved jet widths, model-fitted component sizes, and multifrequency core-size measurements to infer a parabolic-to-conical structural transition in both the jet and counterjet, at (1.23±0.24) pc and (0.97±0.29) pc respectively, corresponding to ~(6–8)×10^3 Schwarzschild radii. Using the jet-to-counterjet brightness ratio and spectral indices under a Doppler-boosting assumption (Eq. 2), they derive an intrinsic velocity field showing overall acceleration to a maximum Lorentz factor Γmax≈2.6 near the collimation break, followed by gradual deceleration. They interpret the spatial coincidence of acceleration and collimation as evidence for a co-spatial sub-parsec acceleration and collimation zone (ACZ), and compare the source with M 87 and NGC 1052. The paper also reports proper motions from radial intensity profiles and discusses possible recollimation signatures.","tokens_in":35401,"tokens_out":4266,"duration_ms":48486,"significance":"If the kinematic interpretation holds, this would add NGC 4261 as the fifth AGN (third LLAGN) with observational evidence for a co-spatial ACZ, strengthening the comparative study of jet acceleration and collimation across the AGN population. The collimation analysis is a genuine strength: the parabolic-to-conical transition is derived from two independent width methods (deconvolved transverse profiles and model-fitted component sizes), is consistent with previous work by Nakahara et al. (2018) and Yan et al. (2023), and is supported by well-documented fits with reduced χ² near unity. The archival data are carefully tabulated, and the authors are explicit about caveats such as low signal-to-noise, beam effects, and limited physical coverage. However, the central claim of a co-spatial ACZ rests on the velocity field, which depends on the assumption that the measured brightness asymmetry is purely Doppler in origin. The paper provides only a qualitative frequency-independence argument against free-free absorption and synchrotron self-absorption, and does not quantitatively demonstrate that residual opacity in the counterjet cannot mimic the inferred acceleration pattern.","major_comments":[{"comment":"The reconstruction of the intrinsic speed β_int assumes that the jet-to-counterjet brightness ratio R_B is entirely due to Doppler boosting, with negligible free-free absorption (FFA) and synchrotron self-absorption (SSA) in the measured region. The only support offered is the qualitative frequency independence of R_B in Fig. 5, whereas the counterjet is known to suffer FFA near the base (§3.1, Fig. 1, Appendix D). Because the R_B peak at ~4–5 mas sets Γmax≈2.6 and drives the acceleration/deceleration pattern, a residual frequency-dependent FFA screen on the counterjet could mimic the inferred acceleration. I request a quantitative test — for example, fitting a free-free absorption model to the multifrequency R_B values, or showing that the spectral turnovers of the jet and counterjet are consistent with zero extra opacity — before the co-spatial ACZ claim can be accepted.","section":"§3.2.1/§3.2.2, Eq. (2)"},{"comment":"The velocity field for epochs without spectral maps is derived using the average spectral index <α> = −0.95 ± 0.47 obtained from the same dataset, and binned values with <α> > −0.5 are excluded. This introduces a mild self-reference and a potential selection bias. The authors should demonstrate that the inferred acceleration and Γmax are robust to (i) using different spectral index values within the quoted uncertainty and (ii) including the excluded bins. If the acceleration pattern changes under these tests, the co-spatial ACZ claim would need to be softened.","section":"§3.2.2"},{"comment":"The claim of an overall acceleration from ~10^3 to ~8×10^3 R_s rests largely on the brightness-ratio peak at ~4–5 mas and the subsequent dip and rise. The paper’s own caveats in §4.4.2 note low S/N and beam-shape effects, but no quantitative estimate is given for how opacity, beam-smearing, or local intensity enhancements would alter the R_B profile. Without such an estimate, the peak could reflect a local emission feature rather than a relativistic-speed maximum. The collimation break is robust, but the kinematic part of the ACZ conclusion is not.","section":"§3.2.2 / Figure 7"},{"comment":"The causality argument Γ φ_open / 2 ≤ sqrt(σ_m) assumes σ_m ≥ 1 and φ_open ≲ 5°, but the magnetization parameter is not independently constrained. This is an interpretive step rather than a measurement, and the conclusion of a causally connected ACZ should be presented with that dependence made explicit. A sensitivity check using a range of σ_m would clarify how much of the ACZ claim rests on this assumption.","section":"§4.2"}],"minor_comments":[{"comment":"The broken power-law fit fixes the sharpness parameter at s=10 following N18; a sentence reporting the sensitivity of the fitted break location and indices to s (e.g., s=5 or s=20) would strengthen confidence in the collimation results.","section":"§3.1.1, Eq. (1)"},{"comment":"The colored solid sticks on the horizontal axes of the right panel are described only in the caption; the text should explicitly define that the brightness-ratio measurements are restricted to regions to the right of these sticks, and the excluded inner region should be marked consistently in the left panel.","section":"§3.2.1 / Figure 5"},{"comment":"The sign convention in Eq. (C1) should be stated explicitly for the counterjet: the reader cannot tell from the formula alone whether the plus sign applies to the receding or approaching side, and the paper’s Table 6 reports positive β_int for both sides.","section":"Appendix C, Eq. (C1)"},{"comment":"The simple power-law fits to the model-fitted width data have reduced χ² of 2.6 and 3.2, which are not ideal; a sentence acknowledging this scatter and its possible origin (e.g., component identification or intrinsic width variations) would be helpful.","section":"§3.1.2"},{"comment":"Minor language issue: “will play an unique role” should be “will play a unique role”; similar small grammatical fixes are needed in a few other places throughout the text.","section":"§4.4.4"}],"recommendation":"major_revision","confidential_remarks":"The collimation analysis is solid and publishable, but the co-spatial ACZ conclusion depends on the Doppler-only interpretation of the brightness ratio. The authors should be asked to provide a quantitative opacity test or to clearly reframe the kinematic claim as tentative. The paper’s claim to be the 'first robust evidence' also deserves scrutiny given that Y23 already suggested a coincidence of acceleration and collimation; the novelty should be stated more precisely."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead of Yan et al. on NGC 4261. The collimation result is the strongest part: with archival VLBA data (some epochs not used before), they measure the jet and counterjet width profiles and get a parabolic-to-conical transition at 0.97-1.23 pc, consistent with N18 and Y23. Two independent width methods—deconvolved Gaussians and model-fit component sizes—give the same break, and the symmetry between the two sides is a good argument for a global confining environment. That section is solid.\n\nThe kinematic part is the soft spot. The velocity field is mostly derived from the jet-to-counterjet brightness ratio via Eq. (2), which assumes the asymmetry is purely Doppler, adopts a single viewing angle of 68± 4 degrees, and requires FFA/SSA to be negligible in the measured region. The stress-test worry is legitimate: if residual FFA eats into the counterjet, R_B is inflated and the acceleration/deceleration/re-acceleration pattern could be an artifact. The paper's defense—R_B is roughly frequency-independent in Fig. 5—is reasonable evidence against strong FFA, but it is not a quantitative opacity bound. They also use an average spectral index from the same dataset for epochs without spectral maps, and they exclude some binned alpha > -0.5 values after the fact. Those are minor choices, but they should be defended more explicitly.\n\nTo be fair, the independent proper-motion points (including counterjet features) are broadly consistent with the brightness-ratio speeds, so the picture is not built on one fragile method. And the collimation break survives even if the velocity field shifts. The co-spatial ACZ conclusion is plausible, but I'd want the authors to either produce a multifrequency opacity estimate or soften the claim to \"suggestive\" until then.\n\nFor the AGN jet community this is a useful addition: the fifth AGN and third LLAGN with a candidate ACZ, and the counterjet proper motions are a nice bonus. The writing is clear and the caveats are out in the open.\n\nI'd send it to a referee. With a quantitative FFA check and a slightly more cautious abstract, it's a solid publication.\n\nBest.","headline":"Solid collimation result, plausible but model-dependent kinematics; the co-spatial ACZ claim needs a quantitative FFA check before it carries the weight the abstract gives it.","tokens_in":35915,"tokens_out":4488,"would_cite":true,"duration_ms":51811,"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":"The twin jets of NGC 4261 accelerate to relativistic speeds and change their collimation shape in the same sub-parsec region, providing the first evidence for a co-spatial acceleration and collimation zone in this low-luminosity AGN.","keywords":["low-luminosity active galactic nuclei","relativistic jets","very long baseline interferometry","high angular resolution","NGC 4261","jet acceleration","jet collimation","Doppler boosting"],"falsifier":"Measure the jet-to-counterjet brightness ratio at 22 and 43 GHz with a high-sensitivity array such as the Very Long Baseline Array or the Global mm-VLBI Array: at these frequencies free-free absorption is negligible, so if the ratio implies a different acceleration pattern from the one derived at 1.4--8.4 GHz, the Doppler-boosting assumption breaks down. Alternatively, a direct measurement of the jet viewing angle that falls outside 64--72 degrees would shift the entire speed curve and could erase the coincidence with the collimation break.","tokens_in":34753,"feed_emoji":"🔭","tokens_out":7626,"duration_ms":75211,"temperature":0.7,"pith_summary":"This paper argues that the two-sided relativistic jet of the nearby galaxy NGC 4261 undergoes most of its acceleration in exactly the same sub-parsec region where its shape changes from a parabola to a cone. Using archival Very Long Baseline Array images between 1.4 and 43 GHz, the authors measure the jet width on both sides and the jet-to-counterjet brightness ratio, from which they reconstruct the jet's speed out to about 20,000 Schwarzschild radii. They find the jet accelerates from roughly 0.3 c near 1,000 Schwarzschild radii to a maximum Lorentz factor of about 2.6 near 8,000 Schwarzschild radii, then slowly decelerates. Because this acceleration region overlaps the parabolic-to-conical collimation break at about 1.2 pc on the approaching side and 0.97 pc on the receding side, they conclude that acceleration and collimation are co-spatial and causally connected, forming a compact acceleration and collimation zone (ACZ). If correct, NGC 4261 becomes the third low-luminosity AGN, after M87 and NGC 315, with direct evidence that both processes act together, which matters for testing magnetohydrodynamic jet-launching models.","feed_headline":"Twin jets of NGC 4261 accelerate and collimate in one sub-parsec zone","feed_subtitle":"The jet peaks at Lorentz factor 2.6 near 8,000 Schwarzschild radii, exactly where its shape turns from parabolic to conical.","key_machinery":"The central object is the jet-to-counterjet brightness ratio $R_B = I_{\\rm jet}/I_{\\rm cjet}$ combined with the spectral index $\\alpha$ and the assumed viewing angle $\\phi_{\\rm view}=68^\\circ\\pm4^\\circ$, converted into an intrinsic speed $\\beta_{\\rm int}$ through the Doppler-boosting formula (Eq. 2). The collimation side is carried by the deconvolved jet width $W(r)$ fitted with a broken power law whose break radius $r_b$ locates the parabolic-to-conical transition; a possible dip in the unbinned width profile at the transition is read as a hint of a recollimation shock. These two independent measurements, kinematic via brightness ratio and structural via width profile, are brought together by comparing $r_b$ with the radius where $\\Gamma\\beta_{\\rm int}$ peaks, establishing the co-spatial ACZ.","core_discovery":"NGC 4261's twin jets accelerate to relativistic speeds and change their collimation geometry in the same sub-parsec volume, providing the first robust evidence for a co-spatial acceleration and collimation zone in this source. The paper derives the collimation profile from deconvolved VLBA jet widths and model-fitted component sizes, finding a clear break from a parabolic width profile ($W \\propto r^{0.48\\text{--}0.51}$) to a conical one ($W \\propto r^{1.06\\text{--}1.15}$) at about $8.1\\times10^{3} R_{\\rm s}$ (1.23 pc) on the approaching side and $6.4\\times10^{3} R_{\\rm s}$ (0.97 pc) on the counterjet side. Using the jet-to-counterjet brightness ratio with the spectral index and a viewing angle of $68^\\circ \\pm 4^\\circ$, it reconstructs an intrinsic-speed field showing overall acceleration from about $0.3c$ at $\\sim10^{3} R_{\\rm s}$ to $\\Gamma_{\\max} \\approx 2.6$ near $8\\times10^{3} R_{\\rm s}$, followed by deceleration to sub-relativistic speeds. The spatial coincidence of the acceleration and collimation break, together with the fulfilment of the condition $\\Gamma\\,\\phi_{\\rm open}/2 \\lesssim \\sqrt{\\sigma_m}$, is interpreted as evidence that the jet is accelerated by magnetic-to-kinetic energy conversion while being confined by external pressure, in a zone confined to $r \\lesssim 1.5$ pc, well inside the Bondi radius of about 99 pc. The authors note that NGC 4261 may be a scaled-down analogue of M87's ACZ, with a lower maximum Lorentz factor suggesting a less magnetized jet base.","pith_inferences":["The brightness-ratio method used here offers a way to measure acceleration zones in other two-sided jets that lack long monitoring campaigns; if applied to a larger sample, it could test whether sub-parsec ACZs are generic in low-luminosity AGN rather than rare.","The frequency-independent behaviour of the brightness-ratio radial profile claimed by the paper suggests that a single well-calibrated epoch at several frequencies could anchor the velocity field, which would make ACZ surveys much cheaper; this is an extension the authors do not explicitly develop.","If the tentative recollimation-shock signature at the transition is real, the deceleration-reacceleration pattern in the speed field may be the kinematic fingerprint of a standing shock, a prediction that future millimeter-VLBI imaging of the same region could check.","A direct test worth pursuing is combining brightness-ratio speeds with actual multi-epoch proper motions of individual jet components on the same scales; this would tell whether the velocity field is steady or time-variable."],"forward_implications":["NGC 4261 becomes the fifth radio-loud AGN, and the third low-luminosity AGN after M87 and NGC 315, with evidence that jet acceleration and collimation occur simultaneously.","The ACZ lies at $r \\lesssim 1.5$ pc (about $10^{4} R_{\\rm s}$), far inside the Bondi radius of roughly 99 pc, implying that the jet is collimated by external pressure from an advection-dominated accretion flow or disk wind rather than by the Bondi-scale environment.","The modest maximum Lorentz factor $\\Gamma_{\\max} \\approx 2.6$ combined with a shallow acceleration profile $\\Gamma \\propto r^{0.21}$ in the inner region suggests the jet base is only weakly magnetized, similar to NGC 315 and unlike M87.","Downstream of the collimation break, the jet decelerates to 0.3--0.5 c at $(2\\text{--}3)\\times10^{4} R_{\\rm s}$, indicating in-situ dissipation or a recollimation shock, possibly seen as a width dip in the transition zone."],"supporting_citations":[{"why":"Supplies the core-shift measurements that locate the central black hole, align the twin jets, and identify the counterjet core used throughout the analysis.","marker":"T. Haga et al. (2015)"},{"why":"Established the earlier parabolic-to-conical transition in NGC 4261 and provided the kiloparsec-scale VLA width measurement that anchors the outer collimation profile.","marker":"S. Nakahara et al. (2018)"},{"why":"Provided four-epoch 15 GHz kinematics and previous collimation-profile measurements that this paper expands with additional epochs and frequencies.","marker":"X. Yan et al. (2023)"},{"why":"Supplied the deconvolved-width method and the brightness-ratio velocity method, as well as the NGC 315 comparison ACZ.","marker":"J. Park et al. (2021)"},{"why":"Provided the M87 acceleration profile and the framework for deriving jet kinematics from alternative methods.","marker":"J. Park et al. (2019)"},{"why":"Measured an apparent jet speed from radial intensity profiles, used as a cross-check for the brightness-ratio speeds.","marker":"B. G. Piner et al. (2001)"},{"why":"Applied the brightness-ratio technique to NGC 315, serving as the methodological template for deriving the velocity field.","marker":"L. Ricci et al. (2022)"},{"why":"Supplies the magnetization condition $\\Gamma\\phi_{\\rm open}/2 \\lesssim \\sqrt{\\sigma_m}$ used to argue that acceleration and collimation are causally connected.","marker":"S. S. Komissarov et al. (2009)"},{"why":"Provided the M87 parabolic-to-conical transition and ACZ comparison that frames the scaled-down analogue interpretation.","marker":"K. Asada & M. Nakamura (2012)"}],"fun_headline_variants":["Same sub-parsec zone accelerates and focuses NGC 4261's twin jets","Jet acceleration and collimation break coincide in NGC 4261","NGC 4261's jet shape change marks its speed peak","One sub-parsec zone controls NGC 4261 jet speed and collimation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derived velocity field assumes that the brighter jet is brighter only because of Doppler boosting: no significant free-free or synchrotron self-absorption contaminates the jet-to-counterjet brightness ratio in the regions used, and the jet keeps a single viewing angle of 68 degrees plus or minus 4 degrees.","fun_headline_variants_meta":{"raw":{"variants":["Same sub-parsec zone accelerates and focuses NGC 4261's twin jets","Jet acceleration and collimation break coincide in NGC 4261","NGC 4261's jet shape change marks its speed peak","One sub-parsec zone controls NGC 4261 jet speed and collimation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000722,"raw_usage":{"total_tokens":3430,"prompt_tokens":1328,"completion_tokens":2102,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":944,"completion_tokens_details":{"reasoning_tokens":2037}},"tokens_in":944,"tokens_out":2102,"duration_ms":15959,"temperature":1.0,"reasoning_tokens":2037,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:37:10.345064+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the jet-to-counterjet brightness ratio at 22 and 43 GHz with a high-sensitivity array such as the Very Long Baseline Array or the Global mm-VLBI Array: at these frequencies free-free absorption is negligible, so if the ratio implies a different acceleration pattern from the one derived at 1.4--8.4 GHz, the Doppler-boosting assumption breaks down. Alternatively, a direct measurement of the jet viewing angle that falls outside 64--72 degrees would shift the entire speed curve and could erase the coincidence with the collimation break.","supporting_citations":[{"cited_title":"G., Jones , D","cited_arxiv_id":null,"evidence_quote":"Measured an apparent jet speed from radial intensity profiles, used as a cross-check for the brightness-ratio speeds."}],"review_version":1}