{"id":"6cddc8ae-175c-4e42-99ec-e18165866e10","arxiv_id":"2412.11440","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Synchrotron cooling becomes dynamically important in magnetically arrested disks above roughly 10^-5.5 Eddington accretion rates, altering the disk magnetic flux and jet efficiency by about a factor of two.","lead":"This paper finds a critical mass accretion rate around 10^-5.5 times the Eddington rate, above which synchrotron radiation efficiently cools magnetically arrested disks around black holes and changes their dynamics and jet efficiency. The result gives a simple analytical threshold that can help interpret simulations and observations of low-luminosity black hole systems.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mdot_crit normalization rests on Te=Tp/3 and the no-cooling virial Te; in a MAD with strong magnetic energy extraction, this could shift the quoted 10^-5.5 threshold by about an order of magnitude.","rationale":"The reader's weakest assumption is exactly the electron-temperature normalization: the no-cooling virial estimate and Te=Tp/3 set the scale of Mdot_crit. This is load-bearing because Eq. 6 depends quadratically on tau=Te/Tp, and the paper's own discussion notes tau may reach 8-10 at higher Mdot. The numerical confirmation uses the same Te=Tp/3, so it cannot settle the absolute normalization. The concern does not invalidate the qualitative claim that a cooling transition exists around low Eddington rates; it only questions the precise value of the threshold. The reader's CONDITIONAL verdict already captures this, so no verdict change is needed. The proposed two-temperature simulation or a self-consistent analytic Te balance would directly test whether the threshold shifts materially. I found no fatal internal inconsistency, and the analytical mass-independence is plausible once the proper dimensionless phi_B is used; the main caveat remains the electron temperature physics.","tokens_in":14526,"tokens_out":40961,"duration_ms":346733,"concrete_test":"Run a two-temperature GRMHD simulation suite over the same Mdot range (1e-7 to 1e-4 Mdot_Edd) for a=0.94 and a=0, evolving Te separately with Coulomb and turbulent heating prescriptions, then identify where phi_B and eta begin to change. If the transition shifts by more than 0.5 dex relative to the Te=Tp/3 run, the quoted Mdot_crit is not robust. Alternatively, re-derive Mdot_crit by solving H = C(T_e) self-consistently with Te set by the two-temperature energy balance, and compare the resulting gamma_e to Eq. 5.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytic threshold in Eq. 6 is fixed by Eq. 5, which assumes (i) the full gravitational binding energy is converted to gas internal energy, ug = n mp GM/r, and (ii) Te = Tp/3. For a MAD, however, a large fraction of accretion power is carried off magnetically: Fig. 4 shows jet efficiencies near 100% or higher for a=0.94, so the gas internal energy is not simply n mp GM/r. Even in the two-temperature picture, Te can be substantially below Tp if electron heating is inefficient (e.g., Kawazura et al. 2019; Rowan et al. 2019). Since Eq. 6 scales as (Te/Tp)^-2, an uncertainty of a factor of 3 in Te/Tp shifts Mdot_crit by nearly an order of magnitude. The numerical confirmation inherits the same Te=Tp/3 prescription, so it cannot independently validate the absolute value of Mdot_crit. If the true electron temperature in the cold phase is lower, the transition moves to higher Mdot; if Te/Tp is closer to unity at the threshold, it moves lower. Thus the central quantitative claim, 10^-5.5 Eddington, is not yet robust to the assumed electron thermodynamics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper argues that in magnetically arrested disks (MADs), synchrotron cooling becomes dynamically important above a critical Eddington-scaled accretion rate. The authors equate the synchrotron cooling rate (Eq. 2) with the gravitational energy gain rate (Eq. 3), use MAD flux saturation to express the magnetic field through the accretion rate, and derive Mdot_crit ≈ 2.8×10^-6 (A α)(φ_B/20)^-2 (τ/0.333)^-2 Mdot_Edd (Eq. 6), independent of black-hole mass. They then present GRMHD simulations with radiative cooling for five spins and target accretion rates 10^-7 to 10^-4 Mdot_Edd, reporting that the MAD parameter φ_B and jet efficiency η change by about a factor of 2 above this threshold (Fig. 4), and that the radial force balance shifts from thermal-pressure support toward magnetic and inertial support as cooling increases (Figs. 5–6).","tokens_in":14787,"tokens_out":8572,"duration_ms":79904,"significance":"The paper offers a compact, transparent heuristic for a potentially important transition in MAD dynamics. Its strengths are the explicit analytic parameterization, the multi-spin simulation suite, and the cross-checks against published radiative efficiencies at low accretion rate. If the threshold were robust to the assumed electron thermodynamics, it would provide a useful organizing result for low-luminosity AGN and X-ray binaries. As it stands, the absolute normalization of Mdot_crit is conditional on the assumed Te/Tp and on a no-cooling virial temperature estimate; the simulations adopt the same electron-temperature prescription and therefore do not independently validate the absolute threshold. The existence of a transition, and its rough Eddington scaling, are nevertheless credible and worth publishing once the main caveats are quantified.","major_comments":[{"comment":"The central quantitative claim, Mdot_crit ≈ 10^-5.5 Mdot_Edd, is fixed by two assumptions that the simulations do not test: (i) the full gravitational binding energy is converted to gas internal energy, u_g = n m_p GM/r, even in the cooling-dominated regime, and (ii) Te/Tp = 1/3. Equation (6) scales as τ^-2, so if Te/Tp at the threshold is closer to unity the threshold shifts down by nearly an order of magnitude, while if Te lies well below Tp it shifts up. The authors themselves call the no-cooling estimate \"crude\" immediately after Eq. (5), and Section 4 acknowledges that τ can increase toward 8–10 at high Mdot, but this is not folded into the error budget. Since Section 3.1 sets Te = Tp/3 in the cooling prescription, the numerical confirmation cannot break this degeneracy. I request either a two-temperature treatment or an explicit conditional statement of the threshold with a propagated uncertainty in τ.","section":null},{"comment":"The empirical position of the transition is measured against the targeted accretion rate, but Fig. 1 (bottom right) shows that the effective, measured accretion rate deviates from the target by up to a factor of about 4 in the cooled runs; this alone shifts the apparent threshold by roughly 0.6 dex. Fig. 4 reports time-averaged φ_B and η without error bars, and each parameter point is a single realization, so the claimed factor-of-2 change and the non-monotonic behavior for negative spins (φ_B increases before decreasing) are not quantitatively supported. Please report variability intervals (e.g., 2σ or interquartile ranges) and, at minimum, consider plotting against the effective accretion rate rather than the targeted one.","section":null},{"comment":"There is an internal tension in the quoted location of the threshold: Eq. (6) and the abstract give Mdot_crit ≈ 10^-5.5 Mdot_Edd, while the Fig. 4 caption states that φ_B changes around 10^-6 Mdot_Edd, a difference of more than half a decade. Because the numerical confirmation is visual against sparse points, the manuscript should either reconcile these numbers or present the transition as a range with an uncertainty derived from the scatter in Fig. 4.","section":null}],"minor_comments":[{"comment":"The definition of τ is inconsistent: Section 4 writes \"τ = Tp/Te = 3,\" while Eqs. (5)–(6) and Section 3.1 use τ = Te/Tp; please correct the notation.","section":null},{"comment":"The abstract states that the surveyed accretion rates range from 10^-7 to 10^-4 Mdot_Edd, but Section 3.1 says the simulations cover 10^-7 to 10^-3.5 Mdot_Edd; these numbers should be harmonized.","section":null},{"comment":"The no-cooling simulations are plotted at a nominal accretion rate of 10^-8 Mdot_Edd; this should be explicitly labeled as a reference value rather than a physically evolved accretion rate.","section":null},{"comment":"The text says that at low accretion rate φ_B and η are \"roughly constants,\" but the right panel of Fig. 4 shows sizable variation even at the lowest rates for the highest spins; please specify the tolerance within which they are considered constant.","section":null},{"comment":"The force-balance analysis is shown only for a = 0.94, while the text asserts that the results are consistent across spins; since the spin dependence of the transition is one of the paper's predictions, a compact cross-spin comparison (or an explicit statement that it is deferred to the companion paper) would strengthen the claim.","section":null}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its assumptions and the heuristic is useful, but the central numerical value of the threshold is not yet robust to the electron-temperature model and to the missing variability estimates in the simulations. I would support publication after a major revision that either reframes Mdot_crit as a conditional scaling with an explicit uncertainty or adds a two-temperature test of the normalization. There is also a risk of overclaiming universality based on a single cooling prescription and one realization per parameter point."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Know this: the paper gives a clean analytical estimate for when synchrotron cooling starts to matter in MADs, and backs it with a spin-resolved GRMHD survey. The claimed transition at ~1e-5.5 Mdot_Edd is physically plausible and consistent with earlier hints from Liska et al. (2024). Read the number as an order-of-magnitude estimate, though; the normalization is sensitive to assumptions about electron temperature that the current simulations cannot independently calibrate.\n\nWhat's genuinely new: the closed-form Mdot_crit expression in terms of phi_B, tau, and Gamma (Eq. 6), and a systematic look at how phi_B and jet efficiency vary across five spins and accretion rates from 1e-7 to 1e-4. The force-balance decomposition (Figs. 5-6) is a useful addition, showing thermal pressure gradient yielding to magnetic and inertial terms as cooling increases. The authors are also upfront about the crudeness of their electron temperature prescription, which I appreciate.\n\nSoft spots, in order of concern. First, the analytic threshold is built on Eq. 5, which assumes the full gravitational binding energy goes into gas internal energy and fixes Te=Tp/3. Since Eq. 6 scales as tau^-2, an uncertainty of a factor of 3 in Te/Tp moves Mdot_crit by nearly an order of magnitude. The simulations use the same tau, so they can show that a transition exists in that model family but cannot pin down the absolute Eddington rate. Second, the numerical evidence is thinner than the abstract suggests: one realization per parameter point, no error bars on phi_B or eta in Fig. 4, and a single BH mass (Sgr A*). The effective accretion rate also deviates from the target by up to a factor of 4 at the high end, which blurs the placement of the threshold. Third, the full numerics are deferred to a companion paper, making independent verification harder.\n\nNone of this sinks the central claim. The qualitative transition is robust across spins, and the factor-of-2 changes are much larger than the typical fluctuation shown for Mdot in Fig. 1. The paper is honest about its assumptions and frames the result as an estimate rather than a precision measurement.\n\nWho is this for? Anyone modeling low-luminosity AGN, Sgr A*, or X-ray binaries with MADs, and people interpreting spin-dependent jet efficiencies in GRMHD. It deserves a serious referee; I'd recommend acceptance after revisions that add error bars, test at least one additional BH mass, and expand the discussion of how the threshold depends on the electron heating prescription.","headline":"A plausible Eddington-scaled cooling threshold for MADs, with an analytic estimate and a spin-resolved GRMHD survey; the qualitative transition holds, but the absolute normalization is hostage to Te/Tp and deserves revision.","tokens_in":15372,"tokens_out":3586,"would_cite":true,"duration_ms":34259,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A critical accretion rate near 10^-5.5 times Eddington separates two dynamical regimes of magnetically arrested disks, with the MAD parameter and jet efficiency changing by about a factor of two across the transition.","keywords":["magnetically arrested disk","radiative cooling","synchrotron emission","GRMHD simulation","jet efficiency","MAD parameter","black hole accretion","Eddington accretion rate"],"falsifier":"Run the same GRMHD setup at $\\dot M = 10^{-5}\\dot M_{\\rm Edd}$ for spin $a=0.94$ while evolving the electron temperature separately instead of fixing $T_e = T_p/3$; if $\\phi_B$ and $\\eta$ stay at their no-cooling values rather than changing by about a factor of two, the threshold is an artifact of the temperature prescription.","tokens_in":14305,"feed_emoji":"🕳️","tokens_out":9991,"duration_ms":79216,"temperature":0.7,"pith_summary":"This paper argues that magnetically arrested disks—accretion flows in which the magnetic field near the black hole is so strong that it regulates the inflow—have a sharp transition set by how fast matter falls in. Below about $\\dot M_{\\rm crit} \\approx 10^{-5.5}\\dot M_{\\rm Edd}$, synchrotron cooling is too weak to matter and the disk behaves as in the standard no-cooling picture. Above that rate, cooling radiates the thermal energy as fast as accretion supplies it, so the disk reconfigures: the MAD parameter and jet efficiency move by about a factor of two, and magnetic forces take over from gas pressure in holding the disk in balance. The authors derive the critical rate analytically and confirm it with GRMHD simulations for five black hole spins and accretion rates from $10^{-7}$ to $10^{-4}\\dot M_{\\rm Edd}$. If correct, the result is a mass-independent, Eddington-scaled dividing line for a whole class of accretion flows.","feed_headline":"Magnetic black-hole disks switch regime at ~10^-5.5 Eddington","feed_subtitle":"Synchrotron cooling kicks in above that rate, changing jet power and flux by about a factor of two.","key_machinery":"The load-bearing relation is the MAD saturation condition $\\phi_B = \\Phi_B/\\sqrt{\\dot M} \\approx 4\\pi B r/\\sqrt{4\\pi \\rho u^r}$, which ties the magnetic field strength near the horizon to the accretion rate; because $\\phi_B$ saturates near a value of roughly 15–30 in a MAD, the magnetic field $B$ can be eliminated from the cooling rate and written in terms of $\\dot M$. The electron Lorentz factor $\\gamma_e$ is estimated from the no-cooling internal energy $u_g = n m_p (GM/r)$, giving $\\gamma_e \\approx (3/2)(\\tau/(m_e c^2))(\\Gamma-1)(GM/r)m_p$, with $\\tau = T_e/T_p = 1/3$. Equating synchrotron emission with gravitational energy release then produces the closed-form critical rate of Eq. (6). This cancellation of the magnetic field strength through the MAD parameter is what makes the threshold independent of the black hole mass.","core_discovery":"The paper's central claim is that synchrotron cooling sets a threshold in magnetically arrested disks (MADs): once the mass accretion rate exceeds $\\dot M_{\\rm crit} \\approx 10^{-5.5}\\dot M_{\\rm Edd}$, the synchrotron loss rate equals the gravitational energy gain rate of the accreting gas, so cooling stops being a small correction and starts shaping the disk. The threshold is derived by equating the synchrotron emissivity $Q_s = (4/3)c\\sigma_T\\gamma_e^2\\beta_e^2 U_B n_e$ to the accretion heating rate $(GM/r^2)\\rho u^r$, using the MAD saturation relation $\\phi_B = \\Phi_B/\\sqrt{\\dot M}$ to eliminate the magnetic field. The result depends on the saturated flux and on the electron-to-proton temperature ratio but not on the black hole mass. The simulations show $\\phi_B$ and the jet efficiency $\\eta$ changing by about a factor of two as $\\dot M$ crosses $\\dot M_{\\rm crit}$, with the direction of the change depending on spin. The force balance shifts correspondingly: the thermal pressure gradient weakens and magnetic stresses take over.","pith_inferences":["A population-level prediction follows: observing jet power or radiative efficiency as a function of Eddington ratio across low-luminosity active nuclei and X-ray binaries should show a break near $\\sim 10^{-5.5}$, because the critical rate does not scale with mass.","If the electron-to-proton temperature ratio rises to 8–10 near $10^{-4}\\dot M_{\\rm Edd}$, as the paper notes it may, the same derivation predicts a shifted threshold there; a two-temperature version of the calculation would quantify that shift.","The same flux-saturation logic could be applied to inverse Compton cooling at higher accretion rates, predicting a second dynamical transition where Compton losses become competitive.","The spin dependence of the jet-efficiency change suggests the sign of the black hole spin could be inferred, in principle, from how a source's jet responds to luminosity changes across the threshold."],"forward_implications":["At accretion rates above $\\sim 10^{-5.5}\\dot M_{\\rm Edd}$, the MAD parameter and jet efficiency deviate by roughly a factor of two from their no-cooling values, so radiative cooling must be included in models of such disks.","The transition is independent of black-hole mass, so the same Eddington-scaled switch should appear for stellar-mass black holes and for supermassive black holes.","Above the threshold, the thermal pressure gradient no longer dominates the radial force balance; magnetic stresses, and at the highest rates inertia, become the main support against gravity.","Below the threshold, cooling is dynamically negligible, so the standard non-radiative MAD results remain valid.","The radiative efficiency grows with accretion rate, from about 4% at $10^{-7}\\dot M_{\\rm Edd}$ to higher values above the critical rate."],"supporting_citations":[{"why":"Establishes that the MAD parameter saturates near 15–30, the empirical input that ties magnetic field strength to $\\dot M$ in the derivation.","marker":"Tchekhovskoy et al. 2011"},{"why":"Provides the baseline radial force balance of MADs (pressure-gradient dominated) that the cooling runs are compared against.","marker":"Chatterjee & Narayan 2022"},{"why":"Supplies the synchrotron emissivity formula used to compute the cooling rate in Eq. (2).","marker":"Rybicki & Lightman 1986"},{"why":"Gives the synchrotron and Bremsstrahlung cooling prescriptions adopted in the simulations.","marker":"Esin et al. 1996"},{"why":"Presents the GRMHD code with which the simulations were run.","marker":"Bégué et al. 2023"},{"why":"Provides the spin dependence of $\\phi_B$ used to predict that $\\dot M_{\\rm crit}$ should be smaller for positive spins.","marker":"Narayan et al. 2022"}],"fun_headline_variants":["Cooling flips black-hole disk dynamics at critical accretion rate","Synchrotron cooling sets MAD disk switch at ~10^-5.5 Eddington","Above ~10^-5.5 Eddington, cooling reshapes black-hole disks","Critical accretion rate turns on cooling effects in black-hole disks","MADs: radiative cooling matters above ~10^-5.5 Eddington"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The threshold formula is built on the no-cooling electron temperature estimate with electrons fixed at one third of the proton temperature, so if the real electrons are cooler the predicted switch rate shifts.","fun_headline_variants_meta":{"raw":{"variants":["Cooling flips black-hole disk dynamics at critical accretion rate","Synchrotron cooling sets MAD disk switch at ~10^-5.5 Eddington","Above ~10^-5.5 Eddington, cooling reshapes black-hole disks","Critical accretion rate turns on cooling effects in black-hole disks","MADs: radiative cooling matters above ~10^-5.5 Eddington"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000307,"raw_usage":{"total_tokens":1814,"prompt_tokens":1057,"completion_tokens":757,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":655}},"tokens_in":673,"tokens_out":757,"duration_ms":6748,"temperature":1.0,"reasoning_tokens":655,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:57:18.237559+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same GRMHD setup at $\\dot M = 10^{-5}\\dot M_{\\rm Edd}$ for spin $a=0.94$ while evolving the electron temperature separately instead of fixing $T_e = T_p/3$; if $\\phi_B$ and $\\eta$ stay at their no-cooling values rather than changing by about a factor of two, the threshold is an artifact of the temperature prescription.","supporting_citations":[],"review_version":1}