REVIEW 3 major objections 5 minor 2 cited by
Revision of two-temperature magnetically arrested flows onto a black hole
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Two-temperature magnetically arrested disk simulations around a black hole produce millimeter emission that is converged with grid resolution and practically indistinguishable from simple R(beta) electron-temperature models.
desk verdict A careful, honest two-temperature MAD study whose headline convergence claim is somewhat stronger than the evidence supports; still worth a serious referee. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is a two-temperature magnetically arrested disk simulation in which the electron temperature is evolved from an entropy equation with a subgrid heating split; the split, taken from a turbulent-cascade prescription, sends most dissipation to protons in high-$\beta$ plasma and to electrons in low-$\beta$ plasma. The electron-temperature model is compared against the parametric $R(\beta)$ prescription, in which the proton-to-electron temperature ratio is a function of plasma $\beta$ with limiting values $R_{\rm low}$ and $R_{\rm high}$. The argument is carried by post-processing the flows with polarized synchrotron radiative transfer to produce light curves and image metrics, notably the total-intensity modulation index, rotation measure, and fractional linear and circular polarization. The close similarity of these observable distributions between the turbulent-heating models and the R10 parametric models is the evidence for practical indistinguishability.
What would settle it
Run one of the zero-spin MAD models with roughly an order of magnitude more grid points in each direction and compare the distributions of the modulation index, rotation measure, and polarization fractions over the same time interval; if they move outside the run-to-run scatter among the four current resolutions, the claimed convergence plateau is false.
Extended reading notes
Core claim
The central claim is that the radiative output of two-temperature magnetically arrested disk (MAD) simulations around a black hole, with electrons heated by a prescribed turbulent cascade, is well converged at the grid resolutions used in current event-horizon-scale model libraries: total-intensity variability, rotation measure, and linear and circular polarization do not shift systematically as the grid is refined. The paper further claims that these self-consistent two-temperature models are in practice indistinguishable from the parametric $R(\beta)$ models with $R_{\rm high}=10$ in most observables, even though the spatial map of $T_p/T_e$ is not the same as R10. Radiative cooling and nonthermal electron distribution functions have only weak effects on the millimeter emission. When the models are scaled to Sgr A*, none reproduces all observed properties: the models are too variable and too optically thin, an external Faraday screen is needed to match the rotation measure, and among prograde models spin $a_*=0.5$ best recovers circular polarization at both observed frequencies.
Load-bearing premise
The whole comparison treats the numerical dissipation that appears at the grid scale in the ideal simulation as a stand-in for the physical turbulent cascade that heats electrons, and the paper does not independently validate that equivalence.
Editorial extensions
If this is right
- If the convergence result holds, existing model libraries built at these grid resolutions do not need to be regenerated at higher resolution to draw conclusions about millimeter emission statistics.
- If the turbulent-heating models are indeed practically indistinguishable from the R10 parametric models, parameter inference based on $R(\beta)$ electron temperatures remains a serviceable shortcut for interpreting current Sgr A* observations.
- Radiative cooling can be neglected when modeling Sgr A* millimeter emission from MADs with turbulent heating at these accretion rates.
- Circular polarization is the most discriminating observable: it distinguishes thermal from nonthermal electron distributions and favors an intermediate black hole spin ($a_*\approx 0.5$) for Sgr A*.
- None of the modeled MADs matches the observed variability or optical depth of Sgr A*, which points to physical ingredients beyond the thermal turbulent-heating picture.
Reading between the lines
- A practical consequence left implicit in the paper: if the turbulent-heating and R10 models are indistinguishable for current observables, then the remaining mismatch with Sgr A* data is probably not in the electron-temperature law but in other missing physics such as reconnection heating, resistivity, or anisotropic pressure.
- A natural next numerical test is to run the same convergence and comparison pipeline with a reconnection-based heating partition; if it drives electron temperatures lower while still mapping onto an $R(\beta)$-like model, it could resolve the too-variable and too-optically-thin discrepancies without abandoning the parametric approach.
- The weak effect of nonthermal electrons found at millimeter wavelengths may not persist at higher frequencies or for brighter sources; applying the kappa-distribution models to such cases could expose larger circular-polarization differences.
- The appendix's observation that raising $R_{\rm low}$ makes $R(\beta)$ images resemble those of a much larger $R_{\rm high}$ suggests that polarimetric scoring with $R_{\rm low}=1$ fixed may have biased inferred temperature ratios upward; verifying this could remove part of the apparent tension between turbulent-heating models and the best-fit $R(\beta)$ model.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports new 3D GRMHD simulations of two-temperature magnetically arrested disks (MADs) around black holes, using the ebhlight code and the Kawazura et al. (2019) turbulent electron heating prescription. The radiative output is post-processed with ipole to produce multifrequency light curves and polarimetric observables, which are then compared with the standard parametric R(β) electron temperature model, with models including radiative cooling, with a nonthermal f_kappa electron distribution, and finally with ALMA observations of Sgr A*. The principal claims are that the integrated radiative properties (M3, RM, LP, CP) are well converged with respect to grid resolution at a*=0, that the turbulent-heating K models are similar to, and in practice may be indistinguishable from, R(β) models (most closely R10), that radiative cooling has negligible effect at Sgr A* accretion rates, that nonthermal electrons affect mainly circular polarization, and that among the prograde models a*=0.5 at low viewing angles is closest to the ALMA constraints, although no model matches all observables.
Significance. If the convergence claim holds, it materially strengthens the reliability of grid-resolution choices used in EHT-scale model libraries and supports the use of integrated variability and polarization diagnostics for two-temperature MAD models. The claimed closeness of the turbulent-heating K models to the parametric R10 model is also practically significant because it suggests that simpler electron temperature prescriptions may suffice for current Sgr A* image and light-curve modeling. The paper benefits from long-duration simulations (up to 30,000M), a broad spin range, multifrequency polarized radiative transfer, and a direct comparison with ALMA data. The main weaknesses are that the convergence evidence is narrow and that the K-versus-R(β) similarity is presented through qualitative correlation and distribution comparisons rather than quantitative model-comparison statistics.
major comments (3)
- [Section 2.4, Figure 2, Table 1] The abstract's headline claim that radiative properties are 'well converged with respect to the numerical grid resolution' is not fully supported by the evidence shown. The four convergence runs are all at a*=0, and the refinement sequence is not nested: LR (120x120x128) to MR (240x120x128) doubles only N_r; MR to HR (240x240x128) doubles only N_theta; MR to HHR (360x120x192) increases N_r and N_phi but leaves N_theta unchanged. This is not a systematic convergence sequence. In addition, the reported comparison in Figure 2 is restricted to distributions of integrated quantities at 229 GHz, the 86/690 GHz agreement is asserted only in the text, and no image-space, spectral-index, or emissivity-weighted electron-temperature diagnostic is shown. Since the convergence claim is one of the paper's two central results, it should either be narrowed to 'integrated observables at 229 GHz for a*=0' or supported by additional diagnostics and additional spins.
- [Section 2.3] The physical interpretation of the K-models depends on an assumption that is stated but not validated in this work: the numerical truncation-error heating in the ideal GRMHD scheme is identified with the turbulent cascade heating that is partitioned between electrons and protons. The manuscript acknowledges that 'the total viscous heating is produced by truncation errors at the numerical grid level' and defers to Ressler et al. (2015) for the justification. This is not necessarily an error, but it is load-bearing for the claim that the K models represent physically motivated turbulent electron heating. The paper should either present a validation or sensitivity test for this identification, or explicitly frame the convergence and K-versus-R(β) results as properties of the numerical dissipation model rather than of the physical turbulent heating scenario.
- [Section 3.1, Figures 4-5] The second central claim, that the K models 'in practice, may be indistinguishable from the R(β) models,' is supported only by qualitative visual comparison of correlations and distributions. Figure 4 shows that K is strongly correlated with R10 for spectral index and CP but only weakly correlated for RM and LP, which the text itself identifies as the observables most sensitive to electron temperature. No quantitative distribution-overlap metric, likelihood, or model-selection statistic is provided. Given that the claim is stated in the abstract, a formal comparison (for example, a two-sample test on the observable distributions or a scoring against the same ALMA data used in Section 3.3) should be reported, or the wording should be weakened to reflect the qualitative agreement.
minor comments (5)
- [Abstract] The word 'self-consisitent' should be corrected to 'self-consistent'.
- [Section 2.4] The sentence 'Although not shown, this is also true for two other neighboring frequencies of 86 and 690 GHz' is an unsupported parenthetical in a section whose purpose is to establish convergence; the statement should either be removed or the supporting figure should be included.
- [Section 2.1] The text refers to 'mixed modified Kerr-Schield logarithmic coordinates'; the correct spelling is 'Kerr-Schild'.
- [Table 1] The 'Cooling' column contains both 'No' and 'no' entries; capitalization should be made uniform.
- [Section 2.2] The mass-scaling procedure is described as setting M to reproduce an average 229 GHz flux of 2-3 Jy, but the iterative nature of the procedure and the behavior of the 86 GHz flux under this scaling are not described; a sentence clarifying the procedure would improve reproducibility.
Circularity Check
No significant circularity: the K-versus-R(beta) comparison and the grid-resolution convergence claim are empirical results, not reductions to inputs.
full rationale
The paper's load-bearing comparisons are not constructed from one another. The turbulent-heating model K evolves electron temperatures from electron entropy using the Kawazura et al. (2019) heating partition, while the R(beta) models compute electron temperatures from the same GRMHD internal energy via Eqs. (3)-(4); the two prescriptions share only the underlying fluid simulation. The finding that K resembles R10 is presented as an emergent result of the simulations, with Figs. 4-7 showing that the Tp/Te maps of K and R10 actually differ, so the resemblance is not imposed by definition. The resolution-convergence claim is likewise an empirical grid study: Fig. 2 compares four distinct resolutions and several observables (M3, RM, LP, CP), and the paper itself cautions that 'we cannot exclude that our models are located at the resolution plateau.' That limitation is a robustness concern about single-spin, anisotropic refinement and integrated statistics, not circularity. The mass scaling M is fitted to the observed 229 GHz total flux, but the observables used for model comparison are normalized or dimensionless quantities (modulation index, fractional polarization, spectral index, RM), so the models are not 'predicting' the fitted flux. The author's self-citations (the ipole code, the R(beta) prescription from Moscibrodzka et al. 2016, and Yfantis et al. 2024) are code/tool and benchmark references; none is invoked as a uniqueness theorem or as the sole justification for a conclusion. Overall, the derivation chain is self-contained and no circular step can be exhibited.
Assumptions & free parameters
free parameters (3)
- Mass scaling M (per model, equivalent to accretion rate Mdot) =
1.0e-8 to 1.4e-8 Msun/yr for a*=0 convergence runs; 2.0e-9 to 1.4e-8 Msun/yr across fiducial models (Table 1)
- Rlow and Rhigh parameters of the R(beta) electron temperature model =
Rlow = 1 fixed; Rhigh scanned as 1, 10, 40, 160
- kappa shape parameter of the nonthermal f_kappa distribution =
kappa = 4.25
assumptions (5)
- domain assumption The grid-scale numerical dissipation in the ideal GRMHD code is equivalent to the physical turbulent cascade heating in a collisionless plasma.
- domain assumption The Kawazura et al. (2019) prescription for proton-to-electron heating ratio (Eq. 2) correctly describes turbulent dissipation in MADs.
- domain assumption The plasma is collisionless, two-temperature, and free of Coulomb coupling; radiative cooling is negligible for Sgr A* accretion rates.
- standard math Synchrotron emission, absorption, and internal Faraday rotation computed with ipole correctly model the millimeter radiation.
- domain assumption Ideal GRMHD with a constant adiabatic index gamma_ad = 13/9 adequately describes the MAD accretion flow.
Cite this review
Pith. "Pith review of Revision of two-temperature magnetically arrested flows onto a black hole." pith.science (2026). https://pith.science/paper/K4F6RVAE
@misc{pith2026241206492,
author = {Pith},
title = {Pith review of: Revision of two-temperature magnetically arrested flows onto a black hole},
year = {2026},
howpublished = {\url{https://pith.science/paper/K4F6RVAE}},
note = {Machine review of arXiv:2412.06492}
}
abstract
We revisit the radiative properties of 3D general relativistic magnetohydrodynamics (GRMHD) two-temperature magnetically arrested disk (MAD) models in which electrons are heated by a magnetic turbulent cascade. We focus on studying the model emission, whose characteristics include variability in both total intensity and linear/circular polarizations as well as rotation measures at energies around the synchrotron emission peak in millimeter waves. We find that radiative properties of MAD models with turbulent electron heating are well converged with respect to the numerical grid resolution, which has not been demonstrated before. We compare radiation from two-temperature simulations with turbulent heating to single-temperature models with electron temperatures calculated based on commonly used $R~(\beta)$ prescription. We find that the self-consisitent two-temperature models with turbulent heating do not significantly outperform the $R~(\beta)$ models and, in practice, may be indistinguishable from the $R~(\beta)$ models. Accounting for physical effects such as radiative cooling and nonthermal electron distribution function makes a weak impact on properties of millimeter emission. Models are scaled to Sgr A*, an accreting black hole in the center of our galaxy, and compared to the most complete observational datasets. We point out the consistencies and inconsistencies between the MAD models and observations of this source and discuss future prospects for GRMHD simulations.
Figures
Figures from the paper (14 more)
Forward citations
Cited by 2 Pith papers
-
Survey of Radiative, Two-Temperature Magnetically Arrested Simulations of the Black Hole M87* I: Turbulent Electron Heating
Radiative two-temperature MAD simulations of M87* with Kawazura et al. (2019) turbulent electron heating match EHT total intensity but over-predict beam-scale linear polarization (about 30% versus the observed under 1...
-
Two-temperature treatments in magnetically arrested disk GRMHD simulations more accurately predict light curves of Sagittarius A*
Adding separate electron thermodynamics and radiative cooling to magnetically arrested disk simulations of Sgr A* lowers predicted 230 GHz variability by nearly 50%, but still leaves it above observed levels.
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence := #2 '...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
- [3]
-
[4]
Chael , A., Narayan , R., & Johnson , M. D. 2019, , 486, 2873
2019
-
[5]
Chael , A., Rowan , M., Narayan , R., Johnson , M., & Sironi , L. 2018, , 478, 5209
work page 2018
-
[6]
A., Narayan , R., & Sadowski , A
Chael , A. A., Narayan , R., & Sadowski , A. 2017, , 470, 2367
work page 2017
-
[7]
Chan , C.-K., Psaltis , D., \"O zel , F., Narayan , R., & Sadowski , A. 2015, , 799, 1
work page 2015
-
[8]
Chan , H.-S., Chan , C.-k., Prather , B. S., Wong , G. N., & Gammie , C. 2024, , 964, 17
work page 2024
Show all 66 references
-
[9]
M., Mizuno , Y., et al
Cruz-Osorio , A., Fromm , C. M., Mizuno , Y., et al. 2022, Nature Astronomy, 6, 103
2022
-
[10]
2018, , 612, A34
Davelaar , J., Mo \'s cibrodzka , M., Bronzwaer , T., & Falcke , H. 2018, , 612, A34
2018
-
[11]
2019, , 632, A2
Davelaar , J., Olivares , H., Porth , O., et al. 2019, , 632, A2
2019
-
[12]
M., et al
Dexter , J., Jim \'e nez-Rosales , A., Ressler , S. M., et al. 2020 a , , 494, 4168
2020
-
[13]
2020 b , , 497, 4999
Dexter , J., Tchekhovskoy , A., Jim \'e nez-Rosales , A., et al. 2020 b , , 497, 4999
2020
-
[14]
K., Mizuno , Y., Fromm , C
Dihingia , I. K., Mizuno , Y., Fromm , C. M., & Rezzolla , L. 2023, , 518, 405
2023
-
[15]
C., Gammie , C
Dolence , J. C., Gammie , C. F., Mo \'s cibrodzka , M., & Leung , P. K. 2009, , 184, 387
2009
-
[16]
A., & Loeb , A
Emami , R., Anantua , R., Chael , A. A., & Loeb , A. 2021, , 923, 272
2021
-
[17]
2024, , 964, L26
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2024, , 964, L26
2024
-
[18]
2022, , 930, L16
Event Horizon Telescope Collaboration , Akiyama , K., Alberdi , A., et al. 2022, , 930, L16
2022
-
[19]
C., et al
Event Horizon Telescope Collaboration , Akiyama , K., Algaba , J. C., et al. 2021, , 910, L13
2021
-
[20]
Fishbone , L. G. & Moncrief , V. 1976, , 207, 962
1976
-
[21]
F., Quataert , E., & Tchekhovskoy , A
Foucart , F., Chandra , M., Gammie , C. F., Quataert , E., & Tchekhovskoy , A. 2017, , 470, 2240
2017
-
[22]
M., Cruz-Osorio , A., Mizuno , Y., et al
Fromm , C. M., Cruz-Osorio , A., Mizuno , Y., et al. 2022, , 660, A107
2022
-
[23]
2023, , 957, 103
Galishnikova , A., Philippov , A., & Quataert , E. 2023, , 957, 103
2023
-
[24]
F., McKinney , J
Gammie , C. F., McKinney , J. C., & T \'o th , G. 2003, , 589, 444
2003
-
[25]
C., Johnson , M
Gold , R., McKinney , J. C., Johnson , M. D., & Doeleman , S. S. 2017, , 837, 180
2017
-
[26]
Howes , G. G. 2010, , 409, L104
2010
-
[27]
Hunter , J. D. 2007, Computing in Science and Engineering, 9, 90
2007
-
[28]
V., Narayan , R., & Abramowicz , M
Igumenshchev , I. V., Narayan , R., & Abramowicz , M. A. 2003, , 592, 1042
2003
-
[29]
M., & Nathanail , A
Jiang , H.-X., Mizuno , Y., Fromm , C. M., & Nathanail , A. 2023, , 522, 2307
2023
-
[30]
Kawazura , Y., Barnes , M., & Schekochihin , A. A. 2019, Proceedings of the National Academy of Science, 116, 771
2019
-
[31]
Y., Krichbaum , T
Kim , J. Y., Krichbaum , T. P., Lu , R. S., et al. 2018, , 616, A188
2018
-
[32]
I., & Wilson , L
Livadiotis , G., Desai , M. I., & Wilson , L. B., I. 2018, , 853, 142
2018
-
[33]
C., Tchekhovskoy , A., & Blandford , R
McKinney , J. C., Tchekhovskoy , A., & Blandford , R. D. 2012, , 423, 3083
2012
-
[34]
M., Younsi , Z., et al
Mizuno , Y., Fromm , C. M., Younsi , Z., et al. 2021, , 506, 741
2021
-
[35]
& Falcke , H
Mo \'s cibrodzka , M. & Falcke , H. 2013, , 559, L3
2013
-
[36]
2016, , 586, A38
Mo \'s cibrodzka , M., Falcke , H., & Shiokawa , H. 2016, , 586, A38
2016
-
[37]
Mo \'s cibrodzka , M., Falcke , H., Shiokawa , H., & Gammie , C. F. 2014, , 570, A7
2014
-
[38]
& Gammie , C
Mo \'s cibrodzka , M. & Gammie , C. F. 2018, , 475, 43
2018
-
[39]
& Gammie , C
Mo \'s cibrodzka , M. & Gammie , C. F. 2024, , 968, 6
2024
-
[40]
F., Dolence , J
Mo \'s cibrodzka , M., Gammie , C. F., Dolence , J. C., & Shiokawa , H. 2011, , 735, 9
2011
-
[41]
2022, , 511, 3795
Narayan , R., Chael , A., Chatterjee , K., Ricarte , A., & Curd , B. 2022, , 511, 3795
2022
-
[42]
V., & Abramowicz , M
Narayan , R., Igumenshchev , I. V., & Abramowicz , M. A. 2003, , 55, L69
2003
-
[43]
Oliphant , T. E. 2007, Computing in Science and Engineering, 9, 10
2007
-
[44]
Palumbo , D. C. M., Wong , G. N., & Prather , B. S. 2020, , 894, 156
2020
-
[45]
& Gammie , C
Petersen , E. & Gammie , C. 2020, , 494, 5923
2020
-
[46]
Porth , O., Mizuno , Y., Younsi , Z., & Fromm , C. M. 2021, , 502, 2023
2021
-
[47]
& Begelman , M
Proga , D. & Begelman , M. C. 2003, , 592, 767
2003
-
[48]
M., Tchekhovskoy , A., Quataert , E., Chandra , M., & Gammie , C
Ressler , S. M., Tchekhovskoy , A., Quataert , E., Chandra , M., & Gammie , C. F. 2015, , 454, 1848
2015
-
[49]
M., Tchekhovskoy , A., Quataert , E., & Gammie , C
Ressler , S. M., Tchekhovskoy , A., Quataert , E., & Gammie , C. F. 2017, , 467, 3604
2017
-
[50]
Ricarte , A., Gammie , C., Narayan , R., & Prather , B. S. 2023, , 519, 4203
2023
-
[51]
E., Sironi , L., & Narayan , R
Rowan , M. E., Sironi , L., & Narayan , R. 2017, , 850, 29
2017
-
[52]
R., Dolence , J
Ryan , B. R., Dolence , J. C., & Gammie , C. F. 2015, , 807, 31
2015
-
[53]
R., Ressler , S
Ryan , B. R., Ressler , S. M., Dolence , J. C., Gammie , C., & Quataert , E. 2018, , 864, 126
2018
-
[54]
R., Ressler , S
Ryan , B. R., Ressler , S. M., Dolence , J. C., et al. 2017, , 844, L24
2017
-
[55]
2017, , 466, 705
Sadowski , A., Wielgus , M., Narayan , R., et al. 2017, , 466, 705
2017
-
[56]
Scepi , N., Dexter , J., & Begelman , M. C. 2022, , 511, 3536
2022
-
[57]
Tchekhovskoy , A., Narayan , R., & McKinney , J. C. 2011, , 418, L79
2011
-
[58]
T., Olivares , H., Cerutti , B., & Mo \'s cibrodzka , M
Vos , J. T., Olivares , H., Cerutti , B., & Mo \'s cibrodzka , M. 2024, , 531, 1554
2024
-
[59]
C., Hardee , P
Walker , R. C., Hardee , P. E., Davies , F. B., Ly , C., & Junor , W. 2018, , 855, 128
2018
-
[60]
2024, , 682, A97
Wielgus , M., Issaoun , S., Mart \' -Vidal , I., et al. 2024, , 682, A97
2024
-
[61]
2022 a , , 930, L19
Wielgus , M., Marchili , N., Mart \' -Vidal , I., et al. 2022 a , , 930, L19
2022
-
[62]
2022 b , , 665, L6
Wielgus , M., Moscibrodzka , M., Vos , J., et al. 2022 b , , 665, L6
2022
-
[63]
Wong , G. N. & Gammie , C. F. 2022, , 937, 60
2022
-
[64]
2006, Plasma Physics and Controlled Fusion, 48, 203
Xiao , F. 2006, Plasma Physics and Controlled Fusion, 48, 203
2006
-
[65]
I., Mo \'s cibrodzka , M
Yfantis , A. I., Mo \'s cibrodzka , M. A., Wielgus , M., Vos , J. T., & Jimenez-Rosales , A. 2024, , 685, A142
2024
-
[66]
2023, , 519, 340
Zhao , S.-S., Huang , L., Lu , R.-S., & Shen , Z. 2023, , 519, 340
2023
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.