{"id":"db45676d-8b09-4d8a-9927-478ac1f87d8c","arxiv_id":"2507.21241","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":10,"one_line_summary":"The radio jet of the low-luminosity AGN M84 changes from a slowly widening flow to a cone at only about 16,000 Schwarzschild radii, much closer to the black hole than the Bondi radius.","lead":"Astronomers measured the shape of the radio jet from the nearby galaxy M84 and found that it transitions from a modestly focused stream to a wide cone about 16,000 black hole radii from the center. Because M84 is one of the faintest actively accreting black holes known, this result tests how jets form around weakly feeding black holes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The broken power-law break at 1.67e4 rs is not robust against combining data with different frequencies, resolutions, and epochs; re-fitting per epoch would settle it.","rationale":"The reader's weakest assumption correctly identified the unmodeled combination of multi-frequency, multi-epoch data as the main vulnerability. My stress-test sharpens this into a specific mechanism: the break radius may coincide with a transition in instrumental resolution and observing frequency, such that the fitted index change (0.71 to 1.16) is driven by the switch from high-frequency VLBI to low-frequency/VLA measurements rather than by a genuine change in jet shape. The paper does not provide the per-point source-to-beam ratios or separate fits by epoch, so this possibility is untested. Because the same concern motivated the reader's CONDITIONAL verdict, my finding does not change the verdict, but it strengthens the need for the conditional analysis. If the proposed split-by-epoch fits show inconsistent break radii, the paper's headline conclusion about the Bondi-radius challenge would be substantially weakened; if they agree, the conclusion is robust. The paper's other elements—careful phase-referencing, core-shift analysis, and comparison with other LLAGNs—are independent and largely unaffected by this concern, so no outright rejection is warranted without the test.","tokens_in":91,"tokens_out":7004,"duration_ms":106615,"concrete_test":"Fit the broken power-law (Eq. 2) separately to (a) the 2014 VLBA I widths, (b) the 2020–2021 VLBA II + 1980 VLA widths, and (c) only points with Phi0/Phi_b > 1.5. If the best-fit r0, au, and ad from (a) and (b) differ by more than the 1σ errors of the combined fit, or if the break disappears in (c), then the reported transition is an artifact of the heterogeneous dataset rather than a property of the M84 jet.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—a parabolic-to-conical transition at r0 = 15.8 ± 3.0 mas—rests on the assumption that the deconvolved Gaussian width Wj = sqrt(Phi0^2 - Phi_b^2) measured at six VLBA frequencies (VLBA I, 2014), two additional epochs (VLBA II, 2020–2021), and a 1980 VLA 1.4 GHz image all trace the same steady-state jet cross-section (Section 3.1, Figure 3). This assumption is not tested. The inner (r < a few mas) points are observed mainly at high frequencies (15–43 GHz) with sub-mas beams, while the outer points are observed at low frequencies (1.4–8.4 GHz) and, at the largest scales, with a 3.86-arcsec VLA beam. Because synchrotron opacity makes the apparent jet photosphere frequency-dependent (the core shifts by up to ~1 mas between 1.4 and 43 GHz per Section 3.2), low-frequency widths can be systematically biased even after Gaussian deconvolution. The break at 15.8 mas is close to the radius where the data switch from high-frequency VLBI to low-frequency/VLA observations, so it may be an instrumental transition rather than a physical collimation break. The paper neither tabulates the individual Wj measurements nor reports the ratio Phi0/Phi_b per point, so there is no way to check whether the break is caused by one or two low-frequency points with large beams.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper measures the jet collimation profile of the low-luminosity AGN M84 by combining multi-frequency VLBA observations from 2014 (VLBA I), supplementary VLBA data from 2020-2021 (VLBA II), and a 1980 VLA 1.4 GHz archival image. Using Gaussian deconvolution of transverse slices, the authors derive jet widths from ~10^2 to ~10^7 r_s and fit single and broken power-law models. They report a transition from a semi-parabolic profile (W ∝ r^0.71) to a conical shape (W ∝ r^1.16) at r0 = 15.8 ± 3.0 mas ≈ 1.67 × 10^4 r_s, with a reduced chi-square improvement from 2.59 to 1.09. They also measure a core shift between 1.4 and 43 GHz of ≲1 mas, use it to place an upper limit on the black hole position, and estimate a magnetic field strength of ~6.5 mG at 1 pc. The paper then compares M84 with other LLAGNs, finding a tentative anticorrelation between the upstream power-law index and normalized accretion rate, and concludes that the jet break occurs almost two orders of magnitude inside the Bondi radius, challenging a universal Bondi-radius transition.","tokens_in":17261,"tokens_out":3373,"duration_ms":37821,"significance":"If the central measurement is robust, the paper extends jet collimation studies to the lowest-accretion-rate regime yet probed and provides a concrete counterexample to the idea that the parabolic-to-conical transition occurs near the Bondi radius. The paper's strengths include a detailed error budget (Table 3), the use of atmosphere-corrected phase-referencing, a careful core-shift analysis, and an explicit comparison with M87. The authors also report all fitting parameters and make their statistical improvement quantitatively clear (χ²/dof from 2.59 to 1.09). The finding of a break at ~1.67 × 10^4 r_s, well inside the Bondi radius, is astrophysically interesting and, if confirmed, would support models where the collimation break is governed by the accretion-flow pressure profile rather than the Bondi scale alone. The paper is appropriately cautious in presenting the accretion-rate correlations as tentative (p-values are quoted), and the comparison across LLAGNs is a useful compilation for the community.","major_comments":[{"comment":"The width measurements used for the broken power-law fit are not tabulated, and the fit combines VLBA I (six frequencies, 2014), VLBA II (5–88 GHz, 2020–2021), and a 1980 VLA 1.4 GHz image as a single steady-state jet profile. The break radius r0 = 15.8 mas is close to the spatial scale where high-frequency VLBI measurements (sub-mas beams) give way to low-frequency and VLA measurements with much larger beams (Table 1 shows beams from 0.34 mas to 3860 mas). Because the radio core shifts by up to ~1 mas between 1.4 and 43 GHz (Section 3.2), the apparent jet width at a given de-projected radius may depend on observing frequency through opacity effects. The paper does not report Φ0 and Φb for each slice, nor does it test whether the break survives when the data are split by epoch or by frequency. I request that the authors tabulate all width measurements with their fitted Gaussian widths and beam sizes, and that they perform subset fits (e.g., VLBA I only, VLBA II only, high-frequency-only, low-frequency-only) to demonstrate that the broken power-law and its break radius are not an artifact of combining heterogeneous data.","section":"Section 3.1, Figure 3"},{"comment":"The statement 'Despite many years of difference, the data set shows consistency and hence constrains the result' is an assertion that is not quantified. The reduced chi-square improvement from 2.59 to 1.09 is encouraging, but it does not by itself establish that multi-epoch and multi-frequency data trace the same intrinsic jet structure. To support this claim, the authors should compare width measurements at overlapping radii from different epochs (e.g., VLBA I 2014 versus VLBA II 2020/2021 at 5 and 22/24 GHz) and show that the residuals of the broken power-law fit do not correlate with observing frequency or epoch. Without such a test, the fitted indices and break radius may be biased by systematic differences between datasets.","section":"Section 3.1, Section 4.1"},{"comment":"The description of how the jet radius is referenced to the central engine is incomplete. The text says the jet radii are plotted 'with respect to the location of the central engine by assuming de-projected distance due to the inclination angle i = 74°' (Section 3.1), but it does not explicitly state whether the frequency-dependent core shift measured in Section 3.2 is applied to shift the origin for each frequency. The core shift is ≲1 mas, which is small compared with r0 = 15.8 mas but could be significant for the innermost points (r ≲ a few mas) that anchor the upstream power-law index a_u = 0.71. The authors should state clearly how the absolute position of the black hole was set for each dataset and, if no correction was applied, quantify the effect of a ~1 mas origin shift on a_u and r0.","section":"Section 3.1, Section 3.2"}],"minor_comments":[{"comment":"The title contains an apparent typesetting artifact: 'Jet F ormation' should read 'Jet Formation'.","section":"Title"},{"comment":"The broken power-law function is written with a sharpness parameter n, but the fit returns n = 75. The sharpness is not discussed in the text; a sentence explaining why such a large n is preferred (i.e., an almost sharp break) would improve readability.","section":"Section 3.1"},{"comment":"The figure and table do not include the number of data points or the degrees of freedom explicitly; reporting dof alongside χ²/dof would allow the reader to judge the fit quality more directly.","section":"Figure 3"},{"comment":"The Kendall tau p-values for the two correlations (0.23 and 0.48) are quoted, but the sample size (N ≈ 5–6) is not displayed in the figure; adding it would clarify the statistical weight of the tentative trends.","section":"Section 4.2"},{"comment":"The footnotes in Table Appendix II.1 are helpful, but the superscript '3' on n_e(r_B) and T(r_B) is not explained in the table itself; please ensure the note appears in the table caption or as a footnote.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"To the editor: The paper is technically sound in its calibration and error treatment, and the science case is interesting. However, the central measurement—the broken power-law transition—rests on a heterogeneous dataset whose internal consistency is not demonstrated. The requested robustness checks (tabulated widths, subset fits, per-frequency/per-epoch comparisons) are feasible with the existing data and would either strengthen or overturn the main conclusion. I recommend major revision rather than rejection because the issue is a missing analysis, not a fundamental flaw in the technique. The statistical significance of the accretion-rate correlations is appropriately downplayed by the authors, and I do not see a circularity problem that would invalidate the central measurement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper gives us the first good collimation profile for an LLAGN at the very low accretion end: M84, from roughly 100 to 10^7 r_s. That alone is worth having. The broken power-law fit is clean on its face (χ2/dof from 2.59 to 1.09), the atmosphere-corrected phase-referencing is careful, and the core-shift analysis is a real step forward. Placing M84 next to M87 and other LLAGNs makes the claim tangible: the break at about 1.7×10^4 r_s, well inside the Bondi radius, challenges the simple universal-break picture. That is an interesting data point even if not decisive.\n\nThe soft spot is exactly what the reader flagged. The collimation profile combines data at six frequencies over a decade plus a 1980 VLA 1.4 GHz image, and all points are fitted to one steady-state power law. The deconvolved Gaussian width is treated as the local jet radius, but if opacity pushes the low-frequency/VLA widths above the true cross-section at those radii, the outer points are biased upward and the break could be an instrumental feature where the data switch from VLBI to VLA. The paper says the data are consistent, but it does not show the individual Wj values or per-frequency fits, so I cannot verify. That is the main thing to fix.\n\nMinor points: the magnetic-field estimate comes without error bars and uses the fitted opening angle as input; it is suggestive, not a robust constraint. The trend with accretion rate is weak (Kendall tau p-values around 0.2–0.5) and the authors are appropriately cautious. The Bondi-radius comparison depends on assumed BH mass and Bondi parameters, but that is not their fault.\n\nOverall this is a solid, honest measurement paper. It deserves serious peer review. I would want a response to the frequency/epoch robustness concern plus a table of the individual widths before accepting it; if that check holds, the result stands.\n\nRecommendation: yes, send it to review.","headline":"A useful VLBI measurement of the M84 jet collimation profile at the low-accretion extreme, but the parabolic-to-conical break radius needs a robustness check against frequency/epoch systematics before I trust it.","tokens_in":17856,"tokens_out":2834,"would_cite":true,"duration_ms":32403,"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 M84 jet switches from a semi-parabolic to a conical width profile at about 17,000 black-hole radii, far inside its Bondi radius.","keywords":["active galactic nuclei","relativistic jets","jet collimation","low-luminosity AGN","M84","very long baseline interferometry","core shift","radio galaxies"],"falsifier":"Measure the M84 jet width at a single frequency (43 GHz) over the full range from $10^2$ to $10^7\\,r_s$ with one epoch and matched resolution; if the width profile shows no break near $1.67\\times10^4\\,r_s$ within the quoted errors, the claimed transition is an artifact of blending multi-frequency data.","tokens_in":16721,"feed_emoji":"🔭","tokens_out":14463,"duration_ms":130308,"temperature":0.7,"pith_summary":"This study sets out to measure how the jet of M84, one of the least luminous active galactic nuclei for which a jet collimation profile has been measured, widens with distance from its black hole. Combining phase-referenced VLBA images at six frequencies, supplementary archival VLBA data, and a 1980 VLA image, it finds that the jet width does not follow one power law: inside about $1.67\\times10^4$ Schwarzschild radii the width grows as $r^{0.71}$, and outside it grows as $r^{1.16}$. The break is significant because it sits roughly 35 times closer to the black hole than M84's Bondi radius, undercutting the idea that the collimation transition happens where gravity captures the surrounding gas. If correct, M84 becomes the lowest-accretion example of a jet that loses collimation close to the black hole, and the tightest constraint yet on how accretion rate controls jet shaping.","feed_headline":"M84 jet widens into a cone far inside its Bondi radius","feed_subtitle":"New radio maps place the jet's collimation break at 17,000 black-hole radii, about 35 times closer than the Bondi radius.","key_machinery":"The load-bearing object is the jet collimation profile $W_j(r)$, the deconvolved Gaussian width of transverse radio slices plotted against deprojected distance from the black hole, fitted with a broken power law $W_j(r)=W_0\\left[\\left(r/r_0\\right)^{n a_u}+\\left(r/r_0\\right)^{n a_d}\\right]^{1/n}$. The upstream index $a_u$ and break radius $r_0$ are what carry the argument: they distinguish a collimated parabolic jet from a freely expanding cone and locate where the ambient pressure stops shaping the flow. The profile is anchored to the black hole by a frequency-dependent core-shift measurement, which sets an upper limit of $\\lesssim14\\,r_s$ on the distance from the 43 GHz core to the central engine.","core_discovery":"The central claim is that the M84 jet has a two-part geometry: a semi-parabolic inner region $W(r)\\propto r^{0.71\\pm0.03}$ inside $r_0=15.8\\pm3.0$ mas ($\\approx1.67\\times10^4\\,r_s$) and a conical shape $W(r)\\propto r^{1.16\\pm0.01}$ beyond it, with the width at the break $W_0\\approx1.46\\times10^3\\,r_s$. A single power law fits the data poorly ($\\chi^2/{\\rm dof}=2.59$), while the broken power law fits well ($\\chi^2/{\\rm dof}=1.09$). Compared with M87, whose break occurs near $2.5\\times10^5\\,r_s$, M84's break is almost two orders of magnitude closer, and since M84's Bondi radius is $\\sim5.9\\times10^5\\,r_s$, the transition cannot be universally tied to the Bondi radius. The paper also derives a small core shift ($\\lesssim1$ mas between 1.5 and 43 GHz) that places the black hole within $\\lesssim14\\,r_s$ of the 43 GHz core, and a magnetic field strength of $\\sim6.5$ mG at 1 pc, which it argues is sufficient for a magnetically arrested accretion state.","pith_inferences":["If the multi-frequency data are biased by opacity, the true break could be frequency-dependent; a single-frequency 43 GHz monitoring campaign would reveal whether the $r^{0.71}$ to $r^{1.16}$ transition is a structural feature or a spectral artifact.","The paper's interpretation implies that the ambient pressure profile at $\\sim10^4\\,r_s$ must fall steeply in M84; X-ray surface-brightness deprojection of the Bondi sphere could test this independently.","The same analysis applied to Sgr A*, an even lower-accretion black hole with no persistent jet, would be a direct test of the low-accretion trend, predicting that any jet there should break within a few hundred Schwarzschild radii."],"forward_implications":["M84 becomes the least collimated LLAGN jet measured to date, extending the known relation between jet collimation and accretion rate down to an Eddington ratio of about $5\\times10^{-7}$.","The parabolic-to-conical break at $\\sim1.67\\times10^4\\,r_s$, well inside the Bondi radius, rules out a universal Bondi-radius transition and points instead to a steep drop in the confining gas density near the black hole.","If the tentative negative correlation between upstream power-law index and normalized accretion rate holds, lower-accretion jets should be systematically wider at a given radius, a prediction for the next generation of VLBI surveys.","The inferred magnetic flux approaching the magnetically arrested disk threshold implies that even very weakly accreting black holes can launch jets through a well-ordered magnetic field, giving low-power jets the same fundamental engine as powerful ones."],"supporting_citations":[{"why":"supplies the adopted black hole mass of $8.5\\times10^8\\,M_\\odot$, setting the Schwarzschild-radius scale used throughout","marker":"Walsh et al. 2010"},{"why":"gives the Bondi radius $\\sim5.9\\times10^5\\,r_s$ and Bondi accretion rate used to show the break is far inside","marker":"Bambic et al. 2023"},{"why":"provides the M87 jet profile whose break at $\\sim2.5\\times10^5\\,r_s$ is the main comparison","marker":"Asada & Nakamura 2012"},{"why":"adds the 2020-2021 VLBA data at 5, 22, 44 and 88 GHz and reports jet kinematics consistent with acceleration in the parabolic zone","marker":"Wang et al. 2022"},{"why":"contributes the 1980 VLA 1.4 GHz image that extends the profile to about $10^7\\,r_s$","marker":"Laing & Bridle 1987"},{"why":"supplies the $74^\\circ$ viewing angle used to deproject the measured widths","marker":"Meyer et al. 2018"},{"why":"establishes the phase-referencing core-shift method and the M87 reference position used to anchor M84's core","marker":"Hada et al. 2011"}],"fun_headline_variants":["M84 jet widens into a cone at 17,000 Schwarzschild radii","M84 jet's collimation break sits 35x inside its Bondi radius","M84 jet becomes conical 15 times closer than M87's","M84: low-luminosity jet turns cone at 17,000 r_s","M84 jet flares to cone far inside its Bondi radius"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The deconvolved Gaussian width of each radio slice is assumed to be the true jet diameter at that radius, although the data combine six frequencies and four epochs without a model for opacity or time variability.","fun_headline_variants_meta":{"raw":{"variants":["M84 jet widens into a cone at 17,000 Schwarzschild radii","M84 jet's collimation break sits 35x inside its Bondi radius","M84 jet becomes conical 15 times closer than M87's","M84: low-luminosity jet turns cone at 17,000 r_s","M84 jet flares to cone far inside its Bondi radius"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000422,"raw_usage":{"total_tokens":2252,"prompt_tokens":1112,"completion_tokens":1140,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":728,"completion_tokens_details":{"reasoning_tokens":1050}},"tokens_in":728,"tokens_out":1140,"duration_ms":12897,"temperature":1.0,"reasoning_tokens":1050,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:58:28.092740+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the M84 jet width at a single frequency (43 GHz) over the full range from $10^2$ to $10^7\\,r_s$ with one epoch and matched resolution; if the width profile shows no break near $1.67\\times10^4\\,r_s$ within the quoted errors, the claimed transition is an artifact of blending multi-frequency data.","supporting_citations":[{"cited_title":"L., Barth, A","cited_arxiv_id":null,"evidence_quote":"supplies the adopted black hole mass of $8.5\\times10^8\\,M_\\odot$, setting the Schwarzschild-radius scale used throughout"},{"cited_title":"2023, Monthly Notices of the Royal Astronomical Society, 522, 4374","cited_arxiv_id":null,"evidence_quote":"gives the Bondi radius $\\sim5.9\\times10^5\\,r_s$ and Bondi accretion rate used to show the break is far inside"},{"cited_title":"2012, The Astrophysical Journal Letters, 745, L28","cited_arxiv_id":null,"evidence_quote":"provides the M87 jet profile whose break at $\\sim2.5\\times10^5\\,r_s$ is the main comparison"},{"cited_title":"2022, The Astrophysical Journal, 941, 140","cited_arxiv_id":null,"evidence_quote":"adds the 2020-2021 VLBA data at 5, 22, 44 and 88 GHz and reports jet kinematics consistent with acceleration in the parabolic zone"},{"cited_title":"A., & Bridle, A","cited_arxiv_id":null,"evidence_quote":"contributes the 1980 VLA 1.4 GHz image that extends the profile to about $10^7\\,r_s$"},{"cited_title":"T., Petropoulou, M., Georganopoulos, M., et al","cited_arxiv_id":null,"evidence_quote":"supplies the $74^\\circ$ viewing angle used to deproject the measured widths"}],"review_version":1}