{"id":"d2af9c00-767d-4e05-bbef-bde9b9fb3e06","arxiv_id":"2607.11439","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Future VLBI with dynamic range 10^4 from 86–345 GHz can distinguish turbulent vs reconnection electron heating and thermal vs kappa distributions in M87 via spectral indices and jet structure.","lead":"Computer simulations of the black hole in M87 show that future radio telescope arrays can tell apart different ways electrons get heated in the jet and disk. This gives observers concrete targets for distinguishing microphysics with multi-frequency images.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"The claimed distinguishability of heating models is not quantified against realistic VLBI imaging systematics, so the abstract/§6 claim may overstate what the idealised maps actually demonstrate.","rationale":"The Reader correctly flags the unvalidated sub-grid heating fractions as a medium-risk assumption and lands on CONDITIONAL. That concern is real but secondary: even if the PIC-calibrated recipes were perfect, the paper never shows that the spectral-index contrasts survive realistic VLBI sampling and imaging. The strongest claim is therefore only half-supported—idealised maps differ, but the observational pathway that would actually “tell them apart” is asserted rather than demonstrated. A single synthetic-observation recovery test would settle whether the claim holds or needs to be softened to “idealised maps show potential differences.” Because the methods are otherwise standard and no derivation error appears, the verdict remains CONDITIONAL; the concrete test simply tightens the condition that must be met before the claim can be regarded as observationally robust.","tokens_in":21517,"tokens_out":564,"duration_ms":5994,"concrete_test":"Take one pair of snapshots (turbulent vs reconnection, both with hybrid kappa eDF), generate synthetic 86/230 GHz visibilities for a representative ngEHT array at the paper’s DR=10^4, image with standard regularisers, and recompute the spectral-index map. If the jet-axis α difference of ~0.3–0.5 (Fig. 9) falls below ~0.1 or becomes spatially unresolved, the distinguishability claim is not supported by the presented evidence.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (abstract, §6) is that future VLBI with DR ~ 10^4 and 86–345 GHz can tell turbulent vs reconnection heating and thermal vs kappa eDFs apart. Support comes from idealised, time-averaged GRRT maps (Figs. 6, 8–11) that show systematic differences in spectral index and optically-thin slope along the jet axis and in the sheath. Those maps are never passed through a realistic array response, thermal noise, sparse uv-coverage, or imaging regularisation. The paper itself notes that total-intensity images are largely degenerate (§5) and that the strongest differences appear only in spectral-index structure; without a controlled recovery test it remains unproven that those differences survive the actual ngEHT/ngVLA pipeline at the quoted dynamic range. The sub-grid heating recipes (Eqs. 2–3, 8–9) are a secondary, already-flagged limitation; the more immediate load-bearing gap is the missing end-to-end observability demonstration.","agreement_with_reader":"partial"},"referee_report":{"model":"grok-4.5","summary":"This paper presents long-term 3D GRMHD simulations of MAD accretion onto a high-spin Kerr black hole (a*=0.9375) that evolve electron entropy under two sub-grid heating prescriptions (turbulent from Kawazura et al. 2019 and reconnection from Rowan et al. 2017). The snapshots are post-processed with GRRT (BHOSS) using thermal Maxwell–Jüttner and hybrid thermal–kappa electron distribution functions (with kappa and width parametrized from PIC results of Ball et al. 2018 and Meringolo et al. 2023). For M87 parameters the authors compute broadband spectra, 230 GHz images with disk/jet decompositions, 86–230 GHz spectral-index maps, optical-depth profiles, and turnover-frequency/flux maps. They conclude that, while total-intensity images and integrated spectra are largely degenerate, spectral-index structure and optically-thin slopes differ systematically between heating models and eDFs, so that future VLBI arrays with dynamical range ~10^4 over 86–345 GHz have the potential to distinguish them.","tokens_in":21837,"tokens_out":1235,"duration_ms":30561,"significance":"If the reported spectral-index and turnover distinctions survive realistic observing conditions, the work supplies concrete, falsifiable predictions that link sub-grid electron microphysics to multi-frequency VLBI observables of M87 on scales from the photon ring to ~1000 µas. This is timely for ngEHT and ngVLA planning. Strengths include the consistent two-temperature GRMHD implementation, PIC-calibrated heating and kappa recipes, careful component decompositions, optical-depth analysis, radial profiles, and the parameter explorations already present in the appendices. The calculations are performed with publicly documented codes (BHAC, BHOSS) and therefore in principle reproducible.","major_comments":[{"comment":"The abstract and §6 claim that a dynamical range of 1×10^4 over 86–345 GHz is sufficient for future VLBI to distinguish turbulent versus reconnection heating and thermal versus kappa eDFs. This rests entirely on idealised, time-averaged GRRT images and maps (Figs. 6, 8–11) that are never passed through a realistic array response, thermal noise, sparse uv-coverage or imaging regularisation. The paper itself emphasises that total-intensity images are largely degenerate (§5) and that the discriminatory power resides in spectral-index structure; without a controlled recovery test it remains unproven that those differences remain detectable at the quoted dynamic range. Either synthetic ngEHT/ngVLA observations or an explicit qualification of the claim is required.","section":"Abstract and Section 6"},{"comment":"The free parameters ε = 0.5 and r_inj = 10 M are fixed to values previously chosen to fit the broadband SED and jet width, while ṁ is iterated to match 1 Jy of compact 230 GHz flux. Appendix B demonstrates that the spatial extent of the steep-to-flat spectral-index transition depends on both parameters. The main-text results and the abstract claim should quantify how the distinguishability between heating models persists (or degrades) across a plausible range of ε and r_inj rather than presenting primarily the single best-fit case.","section":"Section 3 and Appendix B"},{"comment":"The electron-heating fractions (Eqs. 2–3) and kappa parametrizations (Eqs. 8–9) are taken from idealised PIC simulations and applied as sub-grid models. The authors correctly flag the limitations (omission of compressive fluctuations, full reconnection geometries and non-local effects). Because the quantitative spectral-index differences that underpin the central claim are generated by these prescriptions, a short robustness test—e.g., modest variations of the functional forms within published uncertainties, or a direct comparison with a simple R–β model—would strengthen that the reported distinctions are not artefacts of the specific sub-grid choices.","section":"Section 2, Eqs. (2)–(3) and (8)–(9)"}],"minor_comments":[{"comment":"Several typographical errors remain (e.g., “Appdenix B”, “moti-vated”, inconsistent spacing in “eDF : thermal”). A careful proof-reading pass is needed.","section":"Throughout"},{"comment":"Figure 2 and Figure 4 colour bars use non-standard symbols (æ, Øp, £e); these should be replaced by conventional σ, β_p, Θ_e for readability.","section":"Figures 2 and 4"},{"comment":"The viewing angle is stated as ϑ = 160° in the text but the jet-axis orientation in the images is not explicitly related to the observer’s line of sight; a short clarifying sentence would help.","section":"Section 3"},{"comment":"The floor model and σ_cut = 3 are standard, yet the residual North–South asymmetries after time/azimuthal averaging (Fig. 2) are attributed partly to floors. A quantitative statement of how much of the jet-spine emission is discarded would be useful.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid, carefully executed contribution that expands earlier work by the same group. The missing end-to-end synthetic-observation test is the only load-bearing gap; once addressed (or the claim appropriately qualified) the paper will be suitable for A&A. No concerns about novelty disclosure or citation practices."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a careful extension of the two-temperature MAD pipeline (Fromm, Zhang, Chael, Mościbrodzka) that actually delivers new, usable diagnostics: 86–230 GHz spectral-index maps, turnover-frequency maps, and component-decomposed spectra for turbulent vs reconnection heating plus thermal vs κ eDFs at a*=0.9375. The differences are real in the idealised images—reconnection flattens α closer in, turbulence does better in the outer sheath, κ extends the jet—and the optical-depth and radial-profile checks are clean. That is the part worth citing if you are planning ngEHT/ngVLA strategy papers.\n\nWhat it does well: the GRMHD setup (BHAC, MAD torus, electron entropy evolution) and BHOSS post-processing are standard and well-documented. They separate disk/forward/counter-jet properly, show that total intensity is largely degenerate while spectral index is not, and include a short appendix on r_inj and ε. The free parameters (ṁ normalised to 1 Jy at 230 GHz, ε=0.5 from prior SED fits, σ_cut=3) are conventional and do not force the spectral-index contrast that is the main result.\n\nSoft spots, in proportion. The abstract and §6 claim that DR~10^4 and 86–345 GHz “have the potential to tell the models apart.” The maps support that potential, but they never go through a realistic array response, thermal noise, sparse uv-coverage or imaging regularisation. The paper itself notes total-intensity degeneracy; without a recovery test the claim is a bit ahead of the evidence. Secondary and already flagged: the sub-grid heating fractions and κ(σ,β) recipes are idealised PIC extrapolations; floors force the σ-cut. Neither overturns the maps, but both limit how strongly you can read them as unique fingerprints.\n\nThis is for people who already live in the EHT/ngEHT modelling ecosystem and need concrete multi-frequency targets. It is not a breakthrough, but it is honest computational work with clear figures. I would send it to peer review; a referee can ask for a simple mock-observation test or a clearer statement of the idealised nature of the claim. Worth engaging if you care about heating diagnostics; skip if you only want total-intensity horizon images.","headline":"Solid multi-frequency GRMHD+GRRT diagnostics for M87 heating models; the abstract claim is a bit ahead of the idealised maps, but the work is still useful and referee-ready.","tokens_in":22461,"tokens_out":600,"would_cite":true,"duration_ms":6751,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Future mm-VLBI of M87 can separate turbulent from reconnection electron heating and thermal from non-thermal particle distributions.","keywords":["M87","black-hole accretion","GRMHD","electron heating","kappa distribution","VLBI","spectral index","jets"],"falsifier":"A multi-frequency VLBI campaign of M87 that reaches dynamic range ~10^4 between 86 GHz and 345 GHz and measures spectral-index maps on 50–150 µas scales; if those maps show no systematic difference between the jet sheath and the highly magnetised interior, the claimed separability fails.","tokens_in":22432,"feed_emoji":"🌌","tokens_out":908,"duration_ms":8658,"temperature":0.7,"pith_summary":"The paper asks whether next-generation millimetre VLBI arrays can distinguish how electrons are heated near M87’s black hole and whether those electrons follow a thermal or a hybrid thermal-plus-power-law distribution. It runs three-dimensional general-relativistic magnetohydrodynamic simulations of a magnetically arrested disk around a rapidly spinning Kerr black hole, applying two sub-grid heating prescriptions (turbulent cascade and magnetic reconnection) and both Maxwell–Jüttner and kappa electron distributions. Synthetic multi-frequency images and spectra from 86 GHz to 345 GHz show that total-intensity snapshots alone look similar, yet spectral-index maps, component decompositions, and optically-thin slopes differ systematically between the models. With a dynamic range of order 10^4 these differences become observationally accessible, giving a concrete route to diagnose microphysics that single-frequency images cannot resolve.","feed_headline":"Future VLBI can tell how M87 heats its electrons","feed_subtitle":"Spectral maps from 86–345 GHz separate turbulent vs reconnection heating and thermal vs power-law particles","key_machinery":"Two-temperature MAD GRMHD runs with explicit electron-entropy evolution under turbulent (Kawazura-type) and reconnection (Rowan-type) heating fractions, post-processed with GRRT that includes both Maxwell–Jüttner and kappa distributions whose parameters are set by local magnetisation and plasma beta.","core_discovery":"With a dynamical range of 10^4 and simultaneous coverage from 86 GHz to 345 GHz, future VLBI arrays can observationally separate turbulent versus magnetic-reconnection electron heating and thermal versus hybrid kappa electron distributions in M87, even though the integrated spectra and single-frequency total-intensity images remain largely degenerate.","pith_inferences":["If the same spectral-index contrast appears in other low-luminosity AGN, the method becomes a general diagnostic rather than an M87-specific tool.","Polarimetric extensions of the same multi-frequency campaign would further break remaining degeneracies between heating models once Faraday rotation and ordered-field geometry are included.","The requirement for dynamic range 10^4 sets a concrete performance target for ngEHT and ngVLA array design and calibration strategies."],"forward_implications":["Spectral-index maps between 86 GHz and 230 GHz become the primary observable for distinguishing heating mechanisms, not single-frequency total-intensity images.","Hybrid thermal-plus-kappa models are required to match both the radio-to-NIR spectrum and the extended jet emission at millimetre wavelengths.","The radial transition from steep (disk-dominated) to flat (jet-dominated) spectral index can locate the non-thermal particle injection radius.","Turnover frequency and turnover flux density are largely insensitive to heating model and electron distribution, while the optically thin spectral index remains diagnostic."],"fun_headline_variants":["Future VLBI can separate turbulent vs reconnection heating in M87","86–345 GHz maps with 10^4 range distinguish M87 electron models","ngEHT/ngVLA imaging tells thermal vs kappa electrons apart in M87","Multi-frequency VLBI probes radiation microphysics of M87 jets","Dense Fourier sampling reveals M87 electron heating signatures"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The electron-heating fractions and kappa slopes taken from idealised particle-in-cell simulations remain valid when applied as sub-grid prescriptions to the global, time-dependent magnetically arrested flow.","fun_headline_variants_meta":{"raw":{"variants":["Future VLBI can separate turbulent vs reconnection heating in M87","86–345 GHz maps with 10^4 range distinguish M87 electron models","ngEHT/ngVLA imaging tells thermal vs kappa electrons apart in M87","Multi-frequency VLBI probes radiation microphysics of M87 jets","Dense Fourier sampling reveals M87 electron heating signatures"]},"model":"grok-4.5","effort":"low","cost_usd":0.006816,"raw_usage":{"total_tokens":1712,"prompt_tokens":816,"num_sources_used":0,"completion_tokens":97,"cost_in_usd_ticks":68160000,"prompt_tokens_details":{"text_tokens":816,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":799,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":816,"tokens_out":97,"duration_ms":6792,"temperature":1.0,"reasoning_tokens":799,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T05:35:09.771430+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A multi-frequency VLBI campaign of M87 that reaches dynamic range ~10^4 between 86 GHz and 345 GHz and measures spectral-index maps on 50–150 µas scales; if those maps show no systematic difference between the jet sheath and the highly magnetised interior, the claimed separability fails.","supporting_citations":[],"review_version":1}