{"id":"e4469e3c-50bf-4ef5-bd84-ccc0ed833e09","arxiv_id":"2607.21966","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"A MAD+jet model reproduces the observed ν^{-1} radio core-size scaling of LLAGNs only if jet electrons are mostly non-thermal, and predicts strong suppression of power-law electrons in MADs above L/L_Edd ≈ 3–8×10^{-6}.","lead":"This paper models faint galaxy nuclei as a magnetically arrested disk plus a jet, and uses the model to explain why the radio 'core' of such nuclei shrinks with frequency. The result gives observers a way to read the electron content of jets and disks from radio images.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Synchrotron-boiler suppression threshold (L/L_Edd ≈ 3–8×10^-6) relies on an α-viscosity-based γ_cooling (Eq. 6) that is inconsistent with the MAD stress model; a self-consistent H/R could shift the threshold by an order of magnitude.","rationale":"Agree with the reader that Eq. (6) is the weakest link. This is the single most load-bearing concern because the abstract's final quantitative claim (the suppression threshold) and the fiducial MAD.th assumption for M104 both depend on it, and it is internally inconsistent with the MAD stress model: a MAD's inflow is not α-viscosity-driven, yet Eq. (5) uses V_R from α-viscosity. The authors flag IC neglect but not this inconsistency. The cubic H/R dependence makes the threshold fragile; a test that recomputes γ_cooling from the model's own H/R, β_tot, and mass-conservation V_R would adjudicate. If the threshold moves by an order of magnitude, the central conclusion about where PL electrons are suppressed changes, but the paper's conditional caveats and transparent 'crude estimation' warnings mean a CONDITIONAL verdict is appropriate; no further verdict change. The jet-size ξ_pl,jet diagnostic is real but explicitly conditional, so it is not the primary weak point.","tokens_in":22051,"tokens_out":11620,"duration_ms":108695,"concrete_test":"Recompute γ_cooling without Eq. (5) using the MAD model of Sec. 2.1 itself: for each radius, take H/R and β_tot from the self-consistent MAD solution (e.g., F.-G. Xie & Zdziarski 2019) and set V_R from mass conservation with the outflow profile Eq. (1) rather than the α-viscosity formula. Evaluate the L_bol/L_Edd at which γ_cooling falls below the thermal peak over R < 200 R_g. If the resulting threshold differs from (3–8)×10^-6 by more than a factor of ~3, the headline suppression claim needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing concern is the quantitative suppression threshold for PL electrons in MADs, featured as the abstract's final claim. Equation (6), derived from Eqs. (2)–(5), uses the α-viscosity radial velocity V_R ≈ −α_vis(H/R)^2 R Ω_K (Eq. 5), but Sec. 2.1 states that the MAD stress budget is dominated by ordered magnetic field stress B_z B_φ, with α_vis merely mimicking turbulent stress. The actual inflow in a MAD is governed by magnetic interchange instability, not by α-viscosity, so the accretion timescale in Eq. (2) may not represent the physical timescale. Moreover, H/R = 0.5 and β_tot = 6 are fixed without derivation; the MAD model of Sec. 2.1 yields H/R and β_tot that vary with radius and accretion rate. Since γ_cooling ∝ (H/R)^3 (Eq. 6), a plausible thinner flow with H/R = 0.3 shifts γ_cooling by a factor ~4.6, moving the threshold L/L_Edd upward by nearly an order of magnitude. This threshold underpins the fiducial MAD.th choice for M104 and the general conclusion that MAD.th is favored above the threshold; if it moves, the predicted high-frequency size slope and ν_shift in Figs. 5/6/9 could change. The authors acknowledge the neglect of inverse-Compton cooling but not the α-viscosity/MAD inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a coupled magnetically arrested disk (MAD) plus internal-shock jet model and uses it to predict the broadband spectrum and the frequency-dependent radio core size of LLAGNs. The model is first applied to M104, where the accretion rate is adjusted to reproduce the radio luminosity and the resulting SED and size are compared with observations. The authors then generalize to a 10^9 Msun black hole, normalizing the jet power through the empirical Merloni et al. (2003) Fundamental Plane. They report that a size–frequency slope d ∝ nu^{-0.94} is obtained when more than ~50% of the jet electrons follow a power-law distribution, that a jet-to-MAD transition occurs at a few tens of GHz, and that power-law electrons in the MAD are radiatively suppressed for L_bol/L_Edd ≳ (3–8)×10^-6, favoring a purely thermal MAD.","tokens_in":22379,"tokens_out":10273,"duration_ms":103784,"significance":"The paper is valuable because it makes concrete, falsifiable predictions that link the growing high-frequency VLBI/ALMA observations of LLAGNs to the physics of MADs and jets. If the central claims hold, the radio size–frequency slope becomes a diagnostic of the non-thermal electron fraction in jets, and the luminosity threshold identifies the regime in which thermal MAD models suffice. The M104 comparison is a genuine two-observable check (SED plus size), and the jet size–frequency slope is an actual model output rather than a fitted relation. The model is fully specified and the assumptions are largely transparent. The main weakness is that the quantitative suppression threshold, one of the headline conclusions, rests on a simplified cooling estimate that is not self-consistent with the MAD stress model used elsewhere in the paper.","major_comments":[{"comment":"The central suppression threshold quoted in the abstract and §5 is based on gamma_cooling computed with the viscous radial velocity V_R ≈ −alpha_vis(H/R)^2 R Omega_K (Eq. 5). This is inconsistent with the MAD stress model of §2.1, where accretion is driven by the ordered B_z B_phi stress and interchange instability, with alpha_vis only mimicking unresolved turbulent stress. The inflow time in a MAD need not be the alpha-viscosity time, so Eq. (6) is not a consequence of the model. Moreover, H/R = 0.5 and beta_tot = 6 are fixed in §4.2 without derivation, although §2.1 describes a MAD whose structure varies with radius and accretion rate. Because gamma_cooling ∝ alpha_vis^2 (1+beta_tot) (H/R)^3 R/Mdot, the threshold L_bol/L_Edd ≈ (3–8)×10^-6 can shift by a factor of several under plausible changes (e.g., lowering H/R from 0.5 to 0.3 changes gamma_cooling by a factor of 4.6). The acknowled","section":"§4.2, Eq. (6), Fig. 7"},{"comment":"The general predictions in Figs. 3–6, including the shift frequency nu_shift and the claim that MAD emission dominates at high radio frequencies with a flatter size slope, depend on normalizing the jet power through the empirical M03 Fundamental Plane. The jet power is manually adjusted to satisfy this relation, and for M104 the accretion rate is adjusted to fit the radio luminosity, with xi_pl,jet = 80% adopted because it is difficult to constrain. The size–frequency slope in Fig. 6 is a genuine model output and is not compromised by this normalization, but the relative jet/MAD fluxes, and hence nu_shift and the high-frequency size behavior, are. The paper should include a sensitivity test showing how nu_shift and the high-frequency size slope change when mdot_jet is varied within the FP scatter or deviates from the FP; without this, the claimed generality of the jet-to-MAD transition i","section":"§2.3, Fig. 2 caption, §4.1"},{"comment":"The abstract states that the model 'successfully reproduce[s] a size ∝ nu^{-1} scaling'. The model’s jet-dominated size for xi_pl,jet = 80% is d_jet ∝ nu^{-0.94} (Fig. 6, top panel), whereas the observed M104 size below ~30–40 GHz is d_core ∝ nu^{-1.13±0.04} (Sec. 3). The difference of ~0.2 in the power-law index is not discussed and the model slope lies outside the reported observational error. The claim of consistency is therefore quantitatively loose. Please provide a more careful comparison, e.g., fitting the model and the data over the same frequency range, including the MAD contribution, and stating the model uncertainty, before claiming successful reproduction of the size–frequency relation.","section":"§3 and §4.1, Fig. 6"}],"minor_comments":[{"comment":"Table 1 lists p_pl,mad = 2.2, but §2.1 states that the intrinsic MAD power-law index is fixed to p_mad = 2.5. Please reconcile this inconsistency.","section":"Table 1 vs §2.1"},{"comment":"Typo: 'an testbed' should be 'a testbed'.","section":"§3"},{"comment":"Typo: 'which is an good approximation' should be 'which is a good approximation'.","section":"§2.1"},{"comment":"The text says cooling is important 'as long as mdot_0 > 1×10^-7', while Fig. 7 and the abstract quote a threshold of L_bol/L_Edd ≈ (3–8)×10^-6. Clarify the conversion between mdot_0 and L_bol/L_Edd used to obtain the stated luminosity threshold.","section":"§4.2"},{"comment":"Mathematical notation such as 'νL ν ∝ν ∼+0.23' and 'νL ν ∝ν ∼−0.246' should be typeset as νL_ν ∝ ν^{+0.23} and νL_ν ∝ ν^{-0.246}.","section":"Fig. 8 caption"},{"comment":"The text says 'the jet size show no significant differences if xi_pl,jet ≳ 50%', but only the 10% and 80% cases are shown. It would be helpful to show a 50% case or state the range over which the slope is insensitive.","section":"§4.1, Fig. 6"},{"comment":"The size definition differs between MAD (diameter at half-maximum of cumulative flux) and jet (FWHM of dL/dz), and observed core sizes are usually obtained from Gaussian fits. A brief discussion of how these definitions compare to observational size measurements would improve the quantitative comparison.","section":"§2.4"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The genuinely new piece is the size/location calculation in a coupled MAD-jet model and the conclusion that the core-size slope below the jet/MAD transition is a practical readout of the jet's power-law electron fraction. That's a concrete, falsifiable diagnostic: ξ_pl,jet > 50% gives size ~ ν^{-1}, consistent with M104, while a thermal-dominated jet gives ~ν^{-0.2}. The M104 SED + size reproduction is a real two-observable comparison, even if the accretion rate is tuned to the radio luminosity.\n\nThe soft spots are real but addressable. The biggest one is the 'synchrotron boiler' suppression threshold (L/L_Edd ~ 3-8×10^{-6}) in Section 4.2. The estimate uses α-viscosity radial velocity V_R ~ -α_vis(H/R)^2 R Ω_K inside a MAD whose stress budget is explicitly magnetic (B_z B_phi). That's an internal tension; the authors don't flag it. Because γ_cooling ∝ (H/R)^3, the fixed H/R = 0.5 and β_tot = 6 choices are not innocuous — a thinner flow with H/R = 0.3 shifts the threshold by nearly an order of magnitude, which would undermine the abstract's final claim. This deserves a parameter scan or a self-consistent H/R from the MAD model itself.\n\nThe secondary issue is that the headline ν^{-1} reproduction is conditioned on ξ_pl,jet = 80%, chosen 'with our own preference.' The model slope -0.94 is quoted without uncertainty against the observed -1.13±0.04; a simple bootstrap or chi-square over the parameter range would sharpen the claim. The jet power normalization via the M03 FP is an empirical choice — fine, but it makes the general predictions in Figures 3-6 less than first-principles.\n\nNone of this is fatal. The paper is honest, the calculations are transparent, and the diagnostic is useful. I'd send it to peer review, with the expectation of a moderate revision: address the α-viscosity/MAD inconsistency, provide an uncertainty estimate on the slope, and possibly release the code.\n\nWho it's for: anyone working on LLAGN radio cores, jet composition, or MAD disk models. I'd cite it for the jet-composition diagnostic, though I'd wait until the threshold issue is cleaned up.\n\nRecommendation: engage — assign a serious referee.","headline":"A useful MAD-jet framework for radio core sizes, but the synchrotron-boiler threshold claim has a load-bearing internal inconsistency that should be fixed before publication.","tokens_in":23034,"tokens_out":2536,"would_cite":true,"duration_ms":26371,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A magnetically arrested disk plus jet model reproduces the observed radio core size scaling of low-luminosity active galactic nuclei and predicts that nonthermal electrons in the disk are radiatively suppressed above a narrow luminosity thr","keywords":["low-luminosity active galactic nuclei","magnetically arrested disk","radio core size","core shift","jet composition","power-law electrons","synchrotron boiler","M104"],"falsifier":"Measure the radio core size-frequency slope in a jet-dominated LLAGN whose spectrum independently shows nonthermal jet emission; if the slope is distinctly shallower than about −1 (for example near −0.2) for such a source, the claim that more than 50% power-law electrons produce size∝ν^(−1) is wrong. Conversely, for a LLAGN near L_bol/L_Edd ≈ 10^(−5)–10^(−6), a flat MAD-like size slope or a hard X-ray tail from power-law electrons would indicate that the synchrotron-boiler suppression threshold is too aggressive.","tokens_in":21805,"feed_emoji":"📡","tokens_out":5588,"duration_ms":59704,"temperature":0.7,"pith_summary":"The paper argues that in low-luminosity active galactic nuclei, the radio core size measured as a function of frequency can reveal where the radio emission originates and whether the jet's electrons are nonthermal. Using a coupled magnetically arrested disk and relativistic jet model, the authors reproduce the observed size∝ν^(−1) scaling for M104 between 1 and tens of GHz, provided more than half of the jet electrons follow a power-law energy distribution. They also show that at higher frequencies, disk emission overtakes jet emission and the size-frequency slope flattens. Their central prediction is that for systems with L_bol/L_Edd above about (3–8)×10^(−6), power-law electrons in the disk are rapidly cooled into a thermal distribution by the 'synchrotron boiler' effect, so a purely thermal disk model is sufficient. If correct, radio core-size measurements become a practical diagnostic of jet composition and of whether nonthermal electrons survive in accretion disks.","feed_headline":"Radio core size reveals whether AGN jet electrons are nonthermal","feed_subtitle":"A measured size-frequency slope distinguishes power-law from thermal jet electrons, separating jet from disk radio emission.","key_machinery":"The operative mechanism is the cooling-break Lorentz factor γ_cooling, computed by equating the synchrotron cooling time to the inflow time of the disk (or the dynamical time of the jet). It sets how much of any power-law electron population survives to radiate. Around this, the model assembles an analytic magnetically arrested disk (with thermal and optional power-law electrons) plus a conical internal-shock jet, and computes core size as a VLBI-like cumulative intensity radius for the disk and a FWHM of the jet emissivity along its axis. The relative fluxes and sizes of the two components then produce the frequency-dependent core size and core shift.","core_discovery":"On the paper's own terms, the central claim is that the MAD-jet model explains the size-frequency relation of LLAGN radio cores as a two-component story: below a shift frequency (~30–40 GHz for M104) the jet dominates, and its size scales as ν^(−0.9) to ν^(−1) only if at least half of the jet electrons are nonthermal; a thermal jet gives a much shallower slope (~ν^(−0.23)). Above the shift frequency the magnetically arrested disk dominates, with a flatter size slope. The same model reproduces M104's broadband spectrum and core size. The second claim is that in the disk, strong radiative cooling ('synchrotron boiler') thermalizes power-law electrons whenever L_bol/L_Edd ≳ (3–8)×10^(−6), so mo","pith_inferences":["I infer that the size-slope diagnostic could be tested immediately with existing VLBI data on a few nearby LLAGNs below their shift frequency, since the predicted difference between a −0.2 and a −1 slope is large enough to distinguish.","Because γ_cooling is independent of black hole mass, the same thermalization threshold should apply to stellar-mass black holes in the hard state; radio core-size and X-ray observations of X-ray binaries could extend the test across mass scales.","The quantitative threshold is sensitive to disk thickness; if a thinner MAD (H/R ≈ 0.3) is closer to reality, the luminosity above which power-law electrons are killed shifts up by almost an order of magnitude, leaving a wider population of sources in which nonthermal disk electrons might be visible.","The general LLAGN predictions inherit the empirical radio–X-ray–mass correlation for jet power, so sources that deviate from that correlation should have different shift frequencies and are natural targets for checking the assumed jet–disk coupling."],"forward_implications":["The measured slope of radio core size versus frequency below the shift frequency becomes a composition diagnostic: a slope near −1 means the jet is dominated by nonthermal electrons, while a slope near −0.2 suggests a mostly thermal jet.","Above the shift frequency, the MAD dominates and the size slope flattens, giving observers a direct way to identify the frequency at which the accretion disk overtakes the jet in radio.","For systems above L_bol/L_Edd ≈ (3–8)×10^(−6), the synchrotron-boiler effect predicts power-law electrons in the disk are strongly suppressed, so thermal-only disk models should suffice for most LLAGNs and hard-state X-ray binaries.","The shift frequency ν_shift depends mainly on black hole mass and the jet-to-disk mass-loss ratio, varying only weakly with X-ray luminosity, so high-frequency core-size observations can probe the accretion–ejection coupling.","M104's observed core size and spectrum, including the break near 30–40 GHz, are reproduced by the coupled model, supporting the MAD picture over a conventional weak-field hot flow for this source."],"fun_headline_variants":["Radio core size slope exposes jet electron type","MAD-jet model tells thermal from nonthermal jet electrons","Size-frequency slope reveals when disk outshines jet","Synchrotron boiler thermalizes AGN jet electrons","Jet core size scales as ν^-1 only with power-law electrons"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The claimed luminosity threshold for killing power-law electrons rests on an assumed disk thickness H/R = 0.5 and a viscous radial-velocity formula; since the threshold scales as (H/R)^3, a thinner magnetically arrested disk would move the threshold up by roughly an order of magnitude, weakening the headline suppression claim.","fun_headline_variants_meta":{"raw":{"variants":["Radio core size slope exposes jet electron type","MAD-jet model tells thermal from nonthermal jet electrons","Size-frequency slope reveals when disk outshines jet","Synchrotron boiler thermalizes AGN jet electrons","Jet core size scales as ν^-1 only with power-law electrons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000192,"raw_usage":{"total_tokens":1210,"prompt_tokens":800,"completion_tokens":410,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":544,"completion_tokens_details":{"reasoning_tokens":330}},"tokens_in":544,"tokens_out":410,"duration_ms":4960,"temperature":1.0,"reasoning_tokens":330,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T06:13:12.226291+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radio core size-frequency slope in a jet-dominated LLAGN whose spectrum independently shows nonthermal jet emission; if the slope is distinctly shallower than about −1 (for example near −0.2) for such a source, the claim that more than 50% power-law electrons produce size∝ν^(−1) is wrong. Conversely, for a LLAGN near L_bol/L_Edd ≈ 10^(−5)–10^(−6), a flat MAD-like size slope or a hard X-ray tail from power-law electrons would indicate that the synchrotron-boiler suppression threshold is too aggressive.","supporting_citations":[],"review_version":1}