Pith. sign in

REVIEW 2 major objections 4 minor 1 cited by

Radiative cooling cools the dense black-hole disk far more than the jet sheath, dimming the ring, brightening extended jets, and cutting 230 GHz flux at fixed accretion rate.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-12 16:44 UTC pith:RXOUPILJ

load-bearing objection Solid, useful cooling-vs-noncooling comparison for MAD M87* with hybrid κ imaging; the midplane Te drop is robust, the relative jet brightening is the part that needs the nonthermal-cooling caveat the authors already flag. the 2 major comments →

arxiv 2605.15502 v2 pith:RXOUPILJ submitted 2026-05-15 astro-ph.HE

Impacts of radiative cooling on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

Mingyuan Zhang , Yosuke Mizuno , Indu K. Dihingia , Christian M. Fromm , Ziri Younsi , Hai Yang , Alejandro Cruz-Osorio This is my paper
classification astro-ph.HE
keywords black hole shadowradiative coolingtwo-temperature GRMHDelectron temperaturesynchrotron jetsM87*nonthermal κ distributionspectral energy distribution
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Event Horizon Telescope images of M87* are usually interpreted with hot accretion models that ignore radiative losses. The paper shows that when radiative cooling is included self-consistently in two-temperature GRMHD simulations, the electron temperature collapses in the dense inner disk near the midplane while the jet sheath cools only modestly. The resulting synthetic 230 GHz images therefore display a dimmer ring, more extended and relatively brighter jets, lower total flux, and a spectral energy distribution that peaks at lower frequency. Variability still comes mainly from the midplane and weakens at higher accretion rates. These morphological and spectral shifts are currently below EHT dynamic range but should become visible with next-generation arrays, so cooling cannot be treated as optional when modelling low-luminosity AGN shadows and jets.

Core claim

Radiative cooling sharply lowers electron temperature in the dense equatorial disk (r ≲ 10 rg) while only slightly cooling the jet sheath; consequently the 230 GHz image becomes a dimmer disk with more extended, brighter jets and reduced total flux, and the high-frequency SED is suppressed at fixed mass accretion rate.

What carries the argument

Two-temperature GRMHD simulations that evolve electron entropy with Coulomb coupling plus radiative cooling (bremsstrahlung, cyclo-synchrotron, Compton), followed by GRRT with a hybrid thermal-plus-variable-κ electron distribution, allow direct before/after comparison of density, temperature, images, SEDs and light curves.

Load-bearing premise

Cooling by nonthermal electrons is omitted from the fluid evolution even though nonthermal electrons are used later to compute the images; if those electrons cool efficiently the jet-brightening result could reverse.

What would settle it

A simultaneous 230 GHz image of the M87* ring and jet base with dynamic range high enough to measure the jet-to-disk flux ratio and turnover frequency; if the jet remains faint relative to the ring at the accretion rates required by the total flux, the cooling-induced jet brightening is ruled out.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 4 minor

Summary. The paper compares two-temperature MAD GRMHD simulations (BHAC) of a spinning black hole (a=0.9375) without and with radiative cooling plus Coulomb coupling, for mass accretion rates ṁ=(1–10)×10^{-6} Ṁ_Edd and two electron heating prescriptions (turbulent Kawazura and reconnection Rowan). Cooling is implemented as a local source term (thermal bremsstrahlung, cyclo-synchrotron, multiple Compton). GRRT (BHOSS) at 230 GHz and i=163° for M87* uses a hybrid thermal+variable-κ eDF (PIC-based κ matched to heating, ε=0.5, σ_cut=1). The central claims, supported by density/Θ_e maps, decomposed images, SEDs and light curves (§§3.1–3.6, Appendices A–C), are that cooling sharply lowers electron temperature in the dense midplane (r≲10 r_g) while only slightly cooling the jet sheath, producing a dimmer disk, relatively more extended/brighter jets, lower total 230 GHz flux, reduced high-frequency SED flux at fixed ṁ, and midplane-dominated variability that decreases with rising ṁ.

Significance. If the differential cooling and resulting morphology hold, the work is significant for EHT/ngEHT modeling of M87* and other MAD LLAGNs: it shows that the observed high radiative efficiency requires cooling in electron thermodynamics, that fixed-ṁ comparisons already alter ring/jet contrast and SEDs, and that horizon-scale jet emission (currently below EHT dynamic range) is a concrete ngEHT target. Strengths include the systematic multi-rate, multi-heating design, self-consistent two-temperature evolution, decomposed GRRT images that isolate midplane vs. jet contributions, matched PIC κ prescriptions, and explicit appendices on σ_cut and heating differences. The limitations (no nonthermal cooling, models not flux-matched) are stated, making the results usable as controlled experiments rather than direct fits.

major comments (2)
  1. [§2.1, §2.2 Eqs. (7)–(8), §3.4] §2.1 explicitly omits radiative cooling from nonthermal electrons (“ignored for simplicity”) while §2.2 and Eqs. (7)–(8) still inject a hybrid thermal+κ eDF with nonthermal efficiency η→0.5 (and PIC-based κ) that dominates jet-sheath emission. Midplane emission is thermal-dominated (η→0), so the sharp Θ_e drop at r≲10 r_g is robust; jet-sheath emission is not. Nonthermal electrons cool more efficiently via synchrotron and inverse-Compton, so the selective under-cooling of the sheath relative to the midplane is load-bearing for the paper’s main imaging claim (Abstract; §3.4; Summary points 1–2) of a dimmer disk plus relatively brighter/more extended jets. §4 correctly flags this as future work, but the present morphology and fractional-jet results require a stronger quantitative caveat or an order-of-magnitude estimate of the missing sink before the claim can be taken at face value.
  2. [§3.4–3.5, Figs. 5–7] All image and SED comparisons are performed at fixed ṁ rather than fixed 230 GHz flux. Resulting total fluxes (Figs. 5–6, 8–9) range from ~0.4 Jy (lowest cooling) to tens of Jy (non-cooling), far above the observed ~0.5 Jy for most models. Consequently the optical-depth regime, self-absorption turnover locations (§3.5), and the relative jet contribution that underpins the “brighter jets” morphology are not guaranteed to map onto M87*. The paper notes that models are “not precisely tuned,” yet the Abstract and Summary still present the morphology change as directly relevant to EHT/ngEHT. Either a flux-matched rescaling discussion or an explicit statement that the reported jet fractions apply only at the simulated (higher) optical depths is needed.
minor comments (4)
  1. [Fig. 1 caption] Figure 1 and several captions contain residual box characters (□) in place of minus signs or subscripts; these should be cleaned for production.
  2. [§2.2] The free parameters ε=0.5, σ_min=0.01, r_inj=10 r_g and the precise polar-angle cuts used for image decomposition (57.3°/122.7°) are stated but not varied; a short sensitivity sentence or pointer to Appendix A would help readers assess robustness.
  3. [§3.3, Fig. 4] In §3.3 the angular-thickness definition of the disk is given; it would be clearer to mark the corresponding θ boundaries consistently on all panels of Fig. 4 rather than only for the non-cooling case.
  4. [Table 1] Table 1 reports modulation indices; adding the absolute mean fluxes (already in Fig. 9) in the same table would make the relative contributions easier to read without cross-referencing.

Circularity Check

0 steps flagged

No significant circularity: cooling-on vs cooling-off morphology and SED differences are independent GRMHD+GRRT outputs, not forced by definition or fit.

full rationale

The paper's load-bearing claims (sharp Te drop in the dense midplane r≲10 rg, milder drop in the jet sheath, consequent dimmer disk + relatively brighter/more extended jets, reduced total 230 GHz flux and high-frequency SED at fixed Ṁ) are obtained by running otherwise identical two-temperature MAD GRMHD simulations with and without the radiative-cooling + Coulomb source terms, then performing GRRT with a hybrid thermal+variable-κ eDF. Density, Θe maps, decomposed images, SEDs and modulation indices are direct numerical outputs of those runs (Sects. 3.2–3.6, Figs. 2–9, Table 1); they are not algebraically identical to any input parameter. Accretion rates are chosen a priori from EHT estimates rather than fitted to the 0.5 Jy flux, and the authors explicitly note that none of the models is tuned to the data. Self-citations (Zhang et al. 2024 for the prior non-cooling setup, Dihingia et al. 2023 / Mizuno et al. 2021 for the two-temperature + cooling modules, Ball/Meringolo/Cruz-Osorio for the κ prescriptions) supply standard numerical methods and context; they do not constitute a uniqueness theorem or ansatz that forces the reported temperature contrast or jet-fraction increase. The deliberate omission of nonthermal cooling is an explicit modeling assumption, not a circular reduction of the claimed result to its inputs. Hence the derivation chain is self-contained.

Axiom & Free-Parameter Ledger

6 free parameters · 5 axioms · 0 invented entities

The central claim rests on standard GRMHD/GRRT machinery plus several hand-chosen microphysical and numerical cutoffs. No new physical entity is invented; the free parameters control how much nonthermal emission and how much highly magnetized plasma contribute to the images. The load-bearing modeling choice is cooling only thermal electrons while imaging with nonthermal electrons.

free parameters (6)
  • ε (κ width magnetic contribution)
    Set to 0.5 by hand to include magnetic energy in w for r > r_inj; directly affects nonthermal jet brightness.
  • σ_cut (magnetization ceiling for emission)
    Default σ_cut=1 zeros emission in highly magnetized regions; Appendix A shows total flux and nearside jet fraction change with σ_cut up to 25.
  • η nonthermal efficiency parameters (ϵ=0.5, σ_min=0.01)
    Control the thermal/κ mix; η→ϵ in the jet and η→0 in the disk by construction of Eq. (8).
  • r_inj injection radius
    Fixed at 10 r_g as a stagnation-surface proxy; enters w and thus the nonthermal distribution.
  • mass accretion rates ṁ
    Three discrete values (1e-6, 5e-6, 1e-5) chosen from EHT estimates rather than fit to 0.5 Jy; absolute fluxes are not observation-matched.
  • black hole spin a and inclination i
    a=0.9375 and i=163° fixed to M87*-like values; morphology conclusions are for this geometry only.
axioms (5)
  • domain assumption Electron thermodynamics evolve via a two-temperature entropy equation with grid-scale heating (turbulent Kawazura or reconnection Rowan) plus Spitzer Coulomb coupling when cooling is on.
    §2.1; standard in the Ressler/Mizuno/Dihingia line of work, not re-derived here.
  • ad hoc to paper Radiative cooling is applied as a local source term (bremsstrahlung, thermal cyclo-synchrotron, multiple Compton) without solving radiation fields or cooling nonthermal electrons.
    §2.1 explicitly ignores nonthermal cooling 'for simplicity'; this is the main modeling idealization.
  • domain assumption Hybrid Maxwell–Jüttner plus variable-κ eDF with PIC-based κ(σ,β) (PIC-TURB or PIC-CS) and emission mix Eq. (7)–(8).
    §2.2; taken from Ball, Meringolo, EHT, Cruz-Osorio, Fromm prior work.
  • domain assumption Single-loop MAD torus (Fishbone–Moncrief, rin=20 rg, rmax=40 rg, γ=4/3) reaches quasi-steady MAD by t≳10 000 tg.
    §2.1; standard MAD setup; results apply to MAD, not SANE.
  • domain assumption Synchrotron GRRT at 230 GHz with FoV 760 µas, M_BH=6.5e9 M_⊙, D=16.8 Mpc, and emission cut at σ>σ_cut.
    §2.2; standard EHT-style imaging assumptions for M87*.

pith-pipeline@v1.1.0-grok45 · 30421 in / 3669 out tokens · 39779 ms · 2026-07-12T16:44:54.090811+00:00 · methodology

0 comments
read the original abstract

The recent 230 GHz observations from the Event Horizon Telescope collaboration have successfully imaged the supermassive black hole shadow of the M87 galaxy. However, the relatively high radiative efficiency observed in the hot accretion flow suggests that radiative cooling is non-negligible and should be considered when calculating the electron temperature. In this study, we compare accretion models without and with radiative cooling across a range of mass accretion rates, $\dot{M}_{\mathrm{BH}} = (1.0 - 10) \times 10^{-6}\,\dot{M}_{\mathrm{Edd}}$, aiming to assess the impact of cooling on the disk structure, electron temperature distribution (eDF), black hole shadow morphology, broadband spectral energy distributions (SEDs), and flux variability. We performed general relativistic radiative transfer (GRRT) calculations on two-temperature, radiative, general relativistic magnetohydrodynamic (GRMHD) simulations, employing different electron heating prescriptions and nonthermal eDFs, analyzing the radiation transfer due to synchrotron emission at 230 GHz with inclination angle of $163^\circ$. These simulations are targeted toward M87$^{*}$. By comparing density profiles, eDFs, GRRT images, SEDs, and time variability between models, we find that the radiative cooling sharply decreases the electron temperature in the dense inner disk around the equatorial plane ($r\lesssim 10\,r_\mathrm{g}$), while slightly reducing jet sheath temperature. Cooling leads to a dimmer disk, more extended and brighter jets, and reduced total flux. For a given accretion rate, cooling reduces the high-frequency flux. Time variability originates primarily from the midplane in both non-cooling and cooling cases and decreases as accretion rates rise. Although currently below the dynamic range of EHT observations, the features identified in this study could be resolved by next-generation arrays such as the ngEHT.

Figures

Figures reproduced from arXiv: 2605.15502 by Alejandro Cruz-Osorio, Christian M. Fromm, Hai Yang, Indu K. Dihingia, Mingyuan Zhang, Yosuke Mizuno, Ziri Younsi.

Figure 1
Figure 1. Figure 1: Accretion rates measured at the event horizon (top) and normalized magnetic flux at the horizon (bottom). The curves in different colors correspond to the different electron heating prescriptions, radiative cooling, and time-averaged accretion rates: without cooling of ˙m = 1 × 10−6 (black), the turbulent heating model with cooling of ˙m = 1 × 10−6 (blue), 5 × 10−6 (green), and 1 × 10−5 (red), and reconnec… view at source ↗
Figure 2
Figure 2. Figure 2: Logarithmic density distribution averaged in time and azimuth over the interval t = 12 000 tg to 15 000 tg. From left to right: without cooling (a), turbulent heating with cooling at ˙m = 1×10−6 (b), ˙m = 5×10−6 (c), and ˙m = 1×10−5 (e), and reconnection heating with cooling at ˙m = 5×10−6 (d). The dashed white and solid black curves represent the magnetization σ = 0.1 and 1, respectively. cooling of the e… view at source ↗
Figure 3
Figure 3. Figure 3: Panels (a) − (f) show the logarithm of the dimensionless electron temperature averaged in time and azimuth over the interval t = 12 000 tg to 15 000 tg. Panels (g) − (j) highlight the differences in linear scale by subtracting the dimensionless electron temperature in the corresponding non-cooling case. The solid black curves represent σ = 1. The dashed skyblue thin to thick curves represent Θe = 10, 32, a… view at source ↗
Figure 4
Figure 4. Figure 4: Angular distribution of time- and azimuthally-averaged dimen￾sionless electron temperature at the given radii on a logarithmic scale. From top to bottom: The radius increases from 7 rg to 20 rg. The curves in different colors correspond to the non-cooling or radiative cooling under different normalized mass accretion rates: the turbulent heating model without cooling (black), with cooling of ˙m = 1 × 10−6 … view at source ↗
Figure 5
Figure 5. Figure 5: Time-averaged GRRT decomposed images from MAD simulations in the interval t = 12 000 tg to 15 000 tg, assuming a black hole spin of a = 0.9375, observed at 230 GHz with an inclination angle of 163◦ . From top to bottom: the accretion rates are M˙ BH/M˙ Edd = 1 × 10−6 , 5 × 10−6 , 5 × 10−6 , and 1 × 10−5 , respectively. The electron heating prescriptions are used for turbulent heating and reconnection heati… view at source ↗
Figure 6
Figure 6. Figure 6: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Light curves of flux at 230 GHz with a 163◦ inclination angle and spin a = 0.9375. Curves are plotted using accretion rates M˙ BH/M˙ Edd = 1×10−6 (top), 5×10−6 (middle), and 1×10−5 (bottom). The solid curves represent the cases without cooling, and the dashed curves correspond to those with radiative cooling. All curves are in turbulent heating and adopt the thermal/variable κ eDF. of time variability on t… view at source ↗
Figure 9
Figure 9. Figure 9: Total flux variation in turbulent heating without cooling (dots) and with radiative cooling (squares) at 230 GHz with a 163◦ inclination angle and spin a = 0.9375. The different colors correspond to the different accretion rates: M˙ BH/M˙ Edd = 1 × 10−6 (black), 5 × 10−6 (red), and 1 × 10−5 (blue). The labels on the x-axis denote the emissions that come from every region: the whole region is depicted first… view at source ↗

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Probing radiation micro-physics in M 87 I. Total intensity and broad-band spectra

    astro-ph.HE 2026-07 conditional novelty 5.0

    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.

Reference graph

Works this paper leans on

69 extracted references · 1 linked inside Pith · cited by 1 Pith paper

  1. [1]

    2025, Living Reviews in Relativity, 28, 4

    Ayzenberg, D., Blackburn, L., Brito, R., et al. 2025, Living Reviews in Relativity, 28, 4

    2025

  2. [2]

    2018, ApJ, 862, 80

    Ball, D., Sironi, L., & Özel, F. 2018, ApJ, 862, 80

    2018

  3. [3]

    2019, ARA&A, 57, 467

    Blandford, R., Meier, D., & Readhead, A. 2019, ARA&A, 57, 467

    2019

  4. [4]

    2025, MNRAS, 537, 2496

    Chael, A. 2025, MNRAS, 537, 2496

    2025

  5. [5]

    Chael, A., Narayan, R., & Johnson, M. D. 2019, MNRAS, 486, 2873

    2019

  6. [6]

    2018, MNRAS, 478, 5209

    Chael, A., Rowan, M., Narayan, R., Johnson, M., & Sironi, L. 2018, MNRAS, 478, 5209

    2018

  7. [7]

    1984, ApJ, 280, 319

    Colpi, M., Maraschi, L., & Treves, A. 1984, ApJ, 280, 319

    1984

  8. [8]

    M., Mizuno, Y ., et al

    Cruz-Osorio, A., Fromm, C. M., Mizuno, Y ., et al. 2022, Nature Astronomy, 6, 103

    2022

  9. [9]

    M., et al

    Cruz-Osorio, A., Meringolo, C., Fromm, C. M., et al. 2026, ApJ, 1001, 227

    2026

  10. [10]

    2018, A&A, 612, A34

    Davelaar, J., Mo´scibrodzka, M., Bronzwaer, T., & Falcke, H. 2018, A&A, 612, A34

    2018

  11. [11]

    2019, A&A, 632, A2

    Davelaar, J., Olivares, H., Porth, O., et al. 2019, A&A, 632, A2

    2019

  12. [12]

    2023, ApJ, 959, L3 De Villiers, J.-P

    Davelaar, J., Ripperda, B., Sironi, L., et al. 2023, ApJ, 959, L3 De Villiers, J.-P. & Hawley, J. F. 2003, ApJ, 589, 458

    2023

  13. [13]

    C., & McKinney, J

    Dexter, J., Agol, E., Fragile, P. C., & McKinney, J. C. 2010, ApJ, 717, 1092

    2010

  14. [14]

    M., et al

    Dexter, J., Jiménez-Rosales, A., Ressler, S. M., et al. 2020, MNRAS, 494, 4168

    2020

  15. [15]

    C., Markoff, S., & Dexter, J

    Dibi, S., Drappeau, S., Fragile, P. C., Markoff, S., & Dexter, J. 2012, MNRAS, 426, 1928

    2012

  16. [16]

    Dihingia, I. K. & Fendt, C. 2025, in New Frontiers in GRMHD Simulations, ed. C. Bambi, Y . Mizuno, S. Shashank, & F. Yuan, 327–360

    2025

  17. [17]

    K., Mizuno, Y ., Fromm, C

    Dihingia, I. K., Mizuno, Y ., Fromm, C. M., & Rezzolla, L. 2023, MNRAS, 518, 405

    2023

  18. [18]

    K., Mizuno, Y ., Fromm, C

    Dihingia, I. K., Mizuno, Y ., Fromm, C. M., & Younsi, Z. 2025, J. Cosmology Astropart. Phys., 2025, 152

    2025

  19. [19]

    K., Vaidya, B., & Fendt, C

    Dihingia, I. K., Vaidya, B., & Fendt, C. 2021, MNRAS, 505, 3596

    2021

  20. [20]

    2010, ApJ, 708, 1545

    Ding, J., Yuan, F., & Liang, E. 2010, ApJ, 708, 1545

    2010

  21. [21]

    A., Narayan, R., Ostriker, E., & Yi, I

    Esin, A. A., Narayan, R., Ostriker, E., & Yi, I. 1996, ApJ, 465, 312 Event Horizon Telescope Collaboration, Akiyama, K., Albentosa-Ruíz, E., et al. 2025, A&A, 693, A265 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al. 2024, A&A, 681, A79 Event Horizon Telescope Collaboration, Akiyama, K., Alberdi, A., et al. 2022, ApJ, 930, L16 Even...

    1996

  22. [22]

    Fishbone, L. G. & Moncrief, V . 1976, ApJ, 207, 962

    1976

  23. [23]

    M., Cruz-Osorio, A., Mizuno, Y ., et al

    Fromm, C. M., Cruz-Osorio, A., Mizuno, Y ., et al. 2022, A&A, 660, A107

    2022

  24. [24]

    Gammie, C. F. 2025, ApJ, 980, 193

    2025

  25. [25]

    F., McKinney, J

    Gammie, C. F., McKinney, J. C., & Tóth, G. 2003, ApJ, 589, 444

    2003

  26. [26]

    2013, ApJ, 773, 118

    Hoshino, M. 2013, ApJ, 773, 118

    2013

  27. [27]

    1977, ApJ, 214, 840

    Ichimaru, S. 1977, ApJ, 214, 840

    1977

  28. [28]

    2025, ApJ, 990, L33

    Imbrogno, M., Meringolo, C., Cruz-Osorio, A., et al. 2025, ApJ, 990, L33

    2025

  29. [29]

    2024, Astrophys

    Imbrogno, M., Meringolo, C., Servidio, S., et al. 2024, Astrophys. J. Lett., 972, L5

    2024

  30. [30]

    D., Akiyama, K., Blackburn, L., et al

    Johnson, M. D., Akiyama, K., Blackburn, L., et al. 2023, Galaxies, 11, 61

    2023

  31. [31]

    Kawazura, Y ., Barnes, M., & Schekochihin, A. A. 2019, Proceedings of the Na- tional Academy of Science, 116, 771

    2019

  32. [32]

    K., Gammie, C

    Leung, P. K., Gammie, C. F., & Noble, S. C. 2011, ApJ, 737, 21

    2011

  33. [33]

    C., Tchekhovskoy, A., & Blandford, R

    McKinney, J. C., Tchekhovskoy, A., & Blandford, R. D. 2012, MNRAS, 423, 3083

    2012

  34. [34]

    2023, Astrophys

    Meringolo, C., Cruz-Osorio, A., Rezzolla, L., & Servidio, S. 2023, Astrophys. J., 944, 122

    2023

  35. [35]

    M., Younsi, Z., et al

    Mizuno, Y ., Fromm, C. M., Younsi, Z., et al. 2021, MNRAS, 506, 741

    2021

  36. [36]

    M., et al

    Mizuno, Y ., Younsi, Z., Fromm, C. M., et al. 2018, Nature Astronomy, 2, 585 Mo´scibrodzka, M. 2025, ApJ, 981, 145 Mo´scibrodzka, M., Falcke, H., & Shiokawa, H. 2016, A&A, 586, A38 Mo´scibrodzka, M., Falcke, H., Shiokawa, H., & Gammie, C. F. 2014, A&A, 570, A7 Mo´scibrodzka, M., Gammie, C. F., Dolence, J. C., & Shiokawa, H. 2011, ApJ, 735, 9 Mo´scibrodzka...

    2018

  37. [37]

    2018, ApJ, 868, 146

    Nakamura, M., Asada, K., Hada, K., et al. 2018, ApJ, 868, 146

    2018

  38. [38]

    V ., & Abramowicz, M

    Narayan, R., Igumenshchev, I. V ., & Abramowicz, M. A. 2003, PASJ, 55, L69

    2003

  39. [39]

    F., & Kulkarni, A

    Narayan, R., SÄ dowski, A., Penna, R. F., & Kulkarni, A. K. 2012, MNRAS, 426, 3241

    2012

  40. [40]

    Narayan, R. & Yi, I. 1994, ApJ, 428, L13

    1994

  41. [41]

    C., Leung, P

    Noble, S. C., Leung, P. K., Gammie, C. F., & Book, L. G. 2007, Classical and Quantum Gravity, 24, S259

    2007

  42. [42]

    2019, A&A, 629, A61

    Olivares, H., Porth, O., Davelaar, J., et al. 2019, A&A, 629, A61

    2019

  43. [43]

    Pandya, A., Zhang, Z., Chandra, M., & Gammie, C. F. 2016, ApJ, 822, 34

    2016

  44. [44]

    2019, ApJS, 243, 26

    Porth, O., Chatterjee, K., Narayan, R., et al. 2019, ApJS, 243, 26

    2019

  45. [45]

    2017, Computational Astrophysics and Cosmology, 4, 1

    Porth, O., Olivares, H., Mizuno, Y ., et al. 2017, Computational Astrophysics and Cosmology, 4, 1

    2017

  46. [46]

    A., Fernández-Ontiveros, J

    Prieto, M. A., Fernández-Ontiveros, J. A., Markoff, S., Espada, D., & González- Martín, O. 2016, MNRAS, 457, 3801

    2016

  47. [47]

    2026, MNRAS, 546, stag148

    Raha, R., Mukhopadhyay, B., & Chatterjee, K. 2026, MNRAS, 546, stag148

    2026

  48. [48]

    M., Tchekhovskoy, A., Quataert, E., Chandra, M., & Gammie, C

    Ressler, S. M., Tchekhovskoy, A., Quataert, E., Chandra, M., & Gammie, C. F. 2015, MNRAS, 454, 1848

    2015

  49. [49]

    M., Tchekhovskoy, A., Quataert, E., & Gammie, C

    Ressler, S. M., Tchekhovskoy, A., Quataert, E., & Gammie, C. F. 2017, MNRAS, 467, 3604

    2017

  50. [50]

    & Zanotti, O

    Rezzolla, L. & Zanotti, O. 2013, Relativistic Hydrodynamics

    2013

  51. [51]

    D., Kovalev, Y

    Ricarte, A., Johnson, M. D., Kovalev, Y . Y ., Palumbo, D. C. M., & Emami, R. 2023, Galaxies, 11, 5

    2023

  52. [52]

    E., Sironi, L., & Narayan, R

    Rowan, M. E., Sironi, L., & Narayan, R. 2017, ApJ, 850, 29

    2017

  53. [53]

    R., Ressler, S

    Ryan, B. R., Ressler, S. M., Dolence, J. C., Gammie, C., & Quataert, E. 2018, ApJ, 864, 126

    2018

  54. [54]

    Salas, L. D. S., Liska, M. T. P., Markoff, S. B., et al. 2025, MNRAS, 538, 698

    2025

  55. [55]

    V ., Penna, R

    Shcherbakov, R. V ., Penna, R. F., & McKinney, J. C. 2012, ApJ, 755, 133

    2012

  56. [56]

    2025, ApJ, 981, L11

    Singh, A., Bégué, D., & Pe’er, A. 2025, ApJ, 981, L11

    2025

  57. [57]

    2026, ApJ, submitted, arXiv:2605.09326 S˛ adowski, A., Narayan, R., McKinney, J

    Singh, A., Begue, D., & Pe’er, A. 2026, ApJ, submitted, arXiv:2605.09326 S˛ adowski, A., Narayan, R., McKinney, J. C., & Tchekhovskoy, A. 2014, MN- RAS, 439, 503 S˛ adowski, A., Narayan, R., Penna, R., & Zhu, Y . 2013, MNRAS, 436, 3856 S˛ adowski, A., Wielgus, M., Narayan, R., et al. 2017, MNRAS, 466, 705

    Pith/arXiv arXiv 2026

  58. [58]

    1965, Physics of fully ionized gases

    Spitzer, L. 1965, Physics of fully ionized gases

    1965

  59. [59]

    & McKinney, J

    Tchekhovskoy, A. & McKinney, J. C. 2012, MNRAS, 423, L55

    2012

  60. [60]

    Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, MNRAS, 418, L79

    2011

  61. [61]

    W., Narayan, R., et al

    Tsunetoe, Y ., Pesce, D. W., Narayan, R., et al. 2025, ApJ, 984, 35 V ourellis, C., Fendt, C., Qian, Q., & Noble, S. C. 2019, ApJ, 882, 2

    2025

  62. [62]

    2006, Plasma Physics and Controlled Fusion, 48, 203

    Xiao, F. 2006, Plasma Physics and Controlled Fusion, 48, 203

    2006

  63. [63]

    B., et al

    Yoon, D., Chatterjee, K., Markoff, S. B., et al. 2020, MNRAS, 499, 3178

    2020

  64. [64]

    M., & Olivares, H

    Younsi, Z., Porth, O., Mizuno, Y ., Fromm, C. M., & Olivares, H. 2020, in Perseus in Sicily: From Black Hole to Cluster Outskirts, ed. K. Asada, E. de Gouveia Dal Pino, M. Giroletti, H. Nagai, & R. Nemmen, V ol. 342, 9–12

    2020

  65. [65]

    2023, ApJ, 942, 47

    Younsi, Z., Psaltis, D., & Özel, F. 2023, ApJ, 942, 47

    2023

  66. [66]

    Younsi, Z., Wu, K., & Fuerst, S. V . 2012, A&A, 545, A13

    2012

  67. [67]

    & Narayan, R

    Yuan, F. & Narayan, R. 2014, ARA&A, 52, 529

    2014

  68. [68]

    A., Poutanen, J., Mikolajewska, J., et al

    Zdziarski, A. A., Poutanen, J., Mikolajewska, J., et al. 1998, MNRAS, 301, 435

    1998

  69. [69]

    M., Younsi, Z., & Cruz-Osorio, A

    Zhang, M., Mizuno, Y ., Fromm, C. M., Younsi, Z., & Cruz-Osorio, A. 2024, A&A, 687, A88 Article number, page 12 of 15 Zhang et al.: Impacts of radiative cooling on the images of black hole shadow and jets Appendix A: Exclusion of magnetized region Due to the density, pressure, and internal energy in simulations may reach the floor value in highly magnetiz...

    2024