Pith. sign in

REVIEW 5 minor 27 references

Sensitivities of Black Hole Images from GRMHD Simulations

T0 review · 0 major / 5 minor · reviewed 2026-07-12 · grok-4.5

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

Pith's one-line read Pixel-wise image derivatives from GRMHD simulations can guide black-hole parameter recovery even with blur and noise.

desk verdict First solid AD sensitivities for real GRMHD black-hole images; feasibility claim holds, limitations are stated honestly. read the letter →

arxiv 2604.11869 v2 pith:NJZK2IAY submitted 2026-04-13 astro-ph.HE

Pedro Naethe Motta , Mário Raia Neto , Cora Prather , Alejandro Cárdenas-Avendaño This is my paper · ORCID
classification astro-ph.HE
keywords blackholeimagingGRMHDradiativetransferautomaticdifferentiationimagesensitivitiesparameterrecoveryEventHorizonTelescopeelectronheating
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

High-fidelity black-hole images are expensive to produce from GRMHD simulations, so comparing models to data across many imaging parameters is costly. This paper shows that a differentiable radiative-transfer code can compute, for each pixel, how the intensity changes with key post-processing parameters such as observer inclination and electron-heating strength. Those sensitivities form a local map from parameter space to image space. The authors map the resulting error surface and find it is structured—with anisotropies and local minima—yet still navigable when the gradients are used. In mock recovery tests, including cases with realistic blur and noise, automatic-differentiation gradients successfully steer a simple optimizer back to the injected parameters. The result is a concrete basis for gradient-informed model–data comparison rather than pure library search.

What carries the argument

Image sensitivities (the Jacobian of the forward model): pixel-wise derivatives dI/dP obtained by integrating the differentiated geodesic and radiative-transfer equations with automatic differentiation, giving a local map from parameter space to image space.

What would settle it

A side-by-side recovery experiment on the same GRMHD snapshot in which automatic-differentiation gradients systematically fail to reduce image error while finite-difference gradients (or a dense library search) succeed under identical blur and noise.

Watch

Extended reading notes

Core claim

Automatic differentiation can produce stable, physically informative pixel-wise derivatives of GRMHD black-hole images with respect to post-processing parameters, and those derivatives remain useful for parameter recovery even after the images are blurred and contaminated with noise.

Load-bearing premise

The radiative-transfer and geodesic equations remain continuous and differentiable enough, after the authors’ chosen step-size rule and optional magnetization cutoff, that the automatic-differentiation derivatives stay both numerically stable and physically meaningful.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

0 major / 5 minor

Summary. The paper introduces and validates the first computation of pixel-wise image sensitivities (Jacobians dI/dP) for GRMHD-based black-hole images using the differentiable radiative-transfer code Jipole. Focusing on two post-processing parameters—observer inclination θ_o and electron-heating parameter R_high—the authors derive the sensitivity equations from the geodesic and covariant transfer equations, integrate them with automatic differentiation, and validate both the forward images (NMSE ~10^{-13} vs. ipole) and the AD derivatives (excellent agreement for R_high; acceptable for θ_o after removing step-size and magnetization-cutoff artifacts). They map the resulting NMSE error landscape, revealing a symmetry-induced local minimum near 180°-θ_o and anisotropic gradients in R_high, then demonstrate that a simple conjugate-gradient scheme guided by these sensitivities recovers the injected parameters under idealized, blurred, and SNR=15 noisy conditions. The work is framed as a feasibility study establishing that AD-computed gradients remain stable and informative for GRMHD imaging, thereby motivating their later incorporation into full Bayesian pipelines.

Significance. If the result holds, the paper supplies a concrete, publicly available technical foundation (Jipole on GitHub/Zenodo) for gradient-informed model–data comparison in EHT-style analyses. The careful diagnosis of numerical discontinuities (Sec. 3.3), the AD-vs-FD validation, and the explicit mapping of the structured error landscape are genuine advances over purely library-based or finite-difference approaches. The mock recoveries under blur and noise, while idealized, demonstrate that the sensitivities remain useful even when the image is degraded, which is a necessary first step before integration into production samplers such as Comrade.jl. The limitations (Stokes I only, two parameters, illustrative CG) are clearly stated and do not undermine the feasibility claim.

minor comments (5)
  1. Throughout the manuscript many figure captions and section headings appear as strings of black squares (e.g., “Figure 1 shows…” followed by garbled text). These are almost certainly encoding artifacts from the draft PDF; they should be cleaned before final production so that every caption is readable.
  2. Sec. 3.2: the NMSE for dI/d heta_o is quoted as ~0.3. While the text correctly attributes the discrepancy to the different geodesic sampling of AD versus FD, a short quantitative statement of the typical absolute residual (or a zoomed inset of the photon-ring region) would help readers judge whether residual differences remain negligible for optimization.
  3. Sec. 5.1: the basin-hopping-inspired stagnation detector is described only in prose. A brief algorithmic box or pseudocode listing the probing directions, step-size reduction factor (0.8), and maximum rounds (15) would improve reproducibility.
  4. Eq. (19) defines NMSE with a sum over pixels; it would be useful to state explicitly whether the sum is restricted to the field of view shown in the figures or includes the full 160 µas camera plane.
  5. The paper cites the authors’ prior Jipole work extensively; a single sentence in the introduction clarifying what is new relative to Naethe Motta et al. (2025) (namely the first GRMHD application and the error-landscape analysis) would help readers unfamiliar with that paper.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: AD sensitivities are derived from the geodesic and transfer equations, validated against an independent code and finite differences, and used only for mock recovery of injected ground-truth parameters.

full rationale

The paper’s central claim is a feasibility result: that automatic-differentiation image sensitivities of GRMHD post-processing remain stable and informative enough to guide local parameter recovery under blur and noise. The derivation chain is self-contained. Equations (1)–(3) and the differentiated forms (8)–(12) are the standard geodesic and covariant radiative-transfer equations; the partials are obtained by automatic differentiation (ForwardDiff) along a single geodesic, not by fitting. Validation is external: forward images match the independent code ipole to NMSE ~ 10^{-13} (Sec. 3.1), and AD derivatives match finite-difference estimates (NMSE ~ 10^{-14} for dI/dR_high and ~0.3 for dI/dθ_o, Sec. 3.2). Numerical artifacts from the original ipole step-size prescription and magnetization cutoffs are diagnosed and replaced by a continuous prescription (Sec. 3.3). The NMSE landscapes (Sec. 4) are computed by direct image comparison and reveal known geometric features (supplementary-angle local minimum) without circular construction. The conjugate-gradient experiments (Sec. 5) recover injected ground-truth values (θ_o = 163°, R_high = 20) under idealized, blurred, and noisy conditions; they are explicitly labeled “illustrative” and “not a final inference strategy.” The only self-citation is to the authors’ prior Jipole paper for the method itself; that citation is not load-bearing for the GRMHD results, which rest on independent validation. No fitted input is renamed a prediction, no uniqueness theorem is imported from the authors, and no ansatz is smuggled in via citation. Score 0 is therefore appropriate.

Assumptions & free parameters 6 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard GRMHD+GRRT physics plus a small set of numerical and modeling choices that the authors explicitly vary and validate. No new physical entities are postulated; free parameters are the usual post-processing knobs plus a few integration-control constants.

free parameters (6)
  • R_high
    Electron-to-ion temperature ratio in weakly magnetized regions; free post-processing parameter whose derivative is computed and fitted.
  • θ_o
    Observer inclination; free geometric parameter whose derivative is computed and fitted.
  • β_crit
    Transition threshold in the R-model, fixed to 1.0 for all runs.
  • σ_cut
    Magnetization cutoff that zeros emission in the funnel; varied to diagnose artifacts.
  • ε_res / step-size prescription
    Geodesic integration step-size control (Δλ_ipole vs. Δλ_new); chosen to restore differentiability.
  • mass unit M
    Density scaling set to produce ~0.5 Jy compact flux.
assumptions (4)
  • domain assumption The geodesic and covariant radiative-transfer equations are continuous and differentiable with respect to the post-processing parameters.
    Stated in Sec. 2.1 as the prerequisite for writing the sensitivity ODEs (Eqs. 8–12).
  • domain assumption Thermal synchrotron emissivity and the Moscibrodzka et al. (2016) R-model correctly describe the electron thermodynamics of the GRMHD snapshot.
    Standard EHT library assumptions used throughout Secs. 2–5.
  • domain assumption Linear interpolation of fluid primitives on the GRMHD grid is adequate for both intensity and its derivatives.
    Implicit in the ray-tracing implementation (Sec. 2).
  • domain assumption Kerr metric with a* = 0.9375 and the chosen SANE snapshot at t = 4500 r_g/c are representative for the sensitivity study.
    Single-snapshot validation and all landscape/fitting experiments (Sec. 3).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Sensitivities of Black Hole Images from GRMHD Simulations." pith.science (2026). https://pith.science/paper/NJZK2IAY

@misc{pith2026260411869,
  author       = {Pith},
  title        = {Pith review of: Sensitivities of Black Hole Images from GRMHD Simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NJZK2IAY}},
  note         = {Machine review of arXiv:2604.11869}
}
abstract

The advent of high-fidelity imaging of supermassive black holes calls for efficient and robust data-analysis methods. In this work, we use $\texttt{Jipole}$, a differentiable, $\texttt{ipole}$-based radiative transfer code, to enable gradient-based analyses of images generated from state-of-the-art general relativistic magnetohydrodynamic (GRMHD) simulations. We compute image sensitivities, i.e., pixel-wise derivatives of the intensity with respect to model parameters, which form the Jacobian of the forward model and define a local map from parameter space to image space. Using these sensitivities in a mock data analysis, we find that GRMHD-based images generate a structured error landscape for parameter fitting, with anisotropies and local minima, making parameter exploration nontrivial but still tractable when guided by gradient information. We characterize this landscape through the Jacobian and assess the feasibility of gradient-based recovery under idealized, blurred, and noisy conditions. Our results show that automatic differentiation-computed image gradients can guide parameter exploration effectively even in the presence of noise. These findings establish a basis for efficient, high-precision model--data comparisons in black hole imaging and motivate the integration of these sensitivities into advanced inference frameworks.

Figures

Figures reproduced from arXiv: 2604.11869 by the authors.

Figure 1
Figure 1. Image comparison between Jipole and ipole using a GRMHD snapshot. Panel (a): Intensity image computed by Jipole for the parameters specified in Section 3. Panel (b): Intensity image computed by ipole for the parameters specified in Section 3. For each image, the total unpolarized compact flux is shown on top of the image. Panel (c): The pixel-wise absolute difference between Jipole and ipole, defined as |Iipole − IJ… view at source ↗
Figure 2
Figure 2. Comparison of the intensity derivative dI/dRhigh for each pixel computed with Jipole. Panel (a) shows the result obtained using automatic differentiation (AD), while panel (b) displays the result using finite-difference (FD). The values of this sensitivities are presented on a symmetric logarithmic scale in both cases. Panel (c) illustrates the logarithmic relative difference between the two methods, with the corres… view at source ↗
Figure 3
Figure 3. Comparison of the intensity derivative dI/dθo for each pixel as computed by Jipole. Panel (a) shows the result obtained using automatic differentiation (AD), while panel (b) displays the finite-difference (FD) estimate; both are presented on a symmetric logarithmic scale. Panel (c) illustrates the logarithmic relative difference between the two methods, with the corresponding NMSE (Equation (19)) reported in the top… view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: The sensitivity dI/dθ computed using different combinations of step size prescriptions and magnetization cutoffs. The top and middle rows correspond to the automatic differentiation (AD) and finite differences (FD) algorithms, respectively, while the bottom row shows t…
Figure 5
Figure 5. Figure 5: The normalized mean squared error (NMSE) between images at an angle θo and an image at 60◦ , 90◦ , 163◦ in panels (a), (b) and (c), respectively, all with Rhigh = 20. The presence of a local minimum is denoted as a black cross in every panel. For panels (a) and (c), we…
Figure 6
Figure 6. Figure 6: The normalized mean squared error (NMSE) between images at a given Rhigh and the ground truth image at Rhigh = 5, 20, and 60 in panels (a), (b), and (c), respectively, all with θo = 163◦ . The presence of a local minimum is denoted as a black cross in every panel. On e…
Figure 7
Figure 7. Figure 7: The behavior (landscape) of the NMSE across the parameter space of electron temperature ratio (Rhigh) and observer angle (θo). (a) A 2D logarithmic contour map illustrating the error distribution. The global minimum is in￾dicated by the red star at Rhigh = 21.0 and θo …
Figure 8
Figure 8. Figure 8: Panel (a): The evolution of the observer’s incli￾nation, θo, as a function of the iteration step during the CG optimization. The black dashed line corresponds the true reference value. The blue line corresponds to an initial guess of 175◦ , while the red line correspon…
Figure 11
Figure 11. Figure 11: The resulting image after the inclusion of blur￾ring and noise, as described in Section 5.3. This corresponds to [PITH_FULL_IMAGE:figures/full_fig_p013_11.png]
Figure 12
Figure 12. Figure 12: Joint fitting of θo and Rhigh under noise and blur using the CG optimization. The transparent lines refer to the fitting without blur and noise portrayed in [PITH_FULL_IMAGE:figures/full_fig_p014_12.png]

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 1 canonical work pages

  1. [1]

    2018, A&A, 613, A2, doi: 10.1051/0004-6361/201732149

    Bronzwaer, T., Davelaar, J., Younsi, Z., et al. 2018, A&A, 613, A2, doi: 10.1051/0004-6361/201732149

    doi:10.1051/0004-6361/201732149 2018

  2. [2]

    2013, The Astrophysical Journal, 777, 13, doi: 10.1088/0004-637X/777/1/13

    Chan, C.-k., Psaltis, D., & ¨Ozel, F. 2013, The Astrophysical Journal, 777, 13, doi: 10.1088/0004-637X/777/1/13

    doi:10.1088/0004-637x/777/1/13 2013

  3. [3]

    2015, The Astrophysical Journal, 799, 1, doi: 10.1088/0004-637x/799/1/1

    Chan, C.-K., Psaltis, D., ¨Ozel, F., Narayan, R., & Sa¸dowski, A. 2015, The Astrophysical Journal, 799, 1, doi: 10.1088/0004-637x/799/1/1

    doi:10.1088/0004-637x/799/1/1 2015

  4. [4]

    N., & Gammie, C

    Dhruv, V., Prather, B., Wong, G. N., & Gammie, C. F. 2025, The Astrophysical Journal Supplement Series, 277, 16, doi: 10.3847/1538-4365/adaea6

    doi:10.3847/1538-4365/adaea6 2025

  5. [5]

    C., Gammie, C

    Dolence, J. C., Gammie, C. F., Mo´ scibrodzka, M., & Leung, P. K. 2009, The Astrophysical Journal Supplement Series, 184, 387, doi: 10.1088/0067-0049/184/2/387 Event Horizon Telescope Collaboration, Akiyama, K.,

    doi:10.1088/0067-0049/184/2/387 2009

  6. [6]

    2019a, 875, L5, doi: 10.3847/2041-8213/ab0f43 Event Horizon Telescope Collaboration, Akiyama, K.,

    Alberdi, A., et al. 2019a, 875, L5, doi: 10.3847/2041-8213/ab0f43 Event Horizon Telescope Collaboration, Akiyama, K.,

    doi:10.3847/2041-8213/ab0f43

  7. [7]

    2019b, ApJL, 875, L6, doi: 10.3847/2041-8213/ab1141 Event Horizon Telescope Collaboration, Akiyama, K.,

    Alberdi, A., et al. 2019b, ApJL, 875, L6, doi: 10.3847/2041-8213/ab1141 Event Horizon Telescope Collaboration, Akiyama, K.,

    doi:10.3847/2041-8213/ab1141

  8. [8]

    C., et al

    Algaba, J. C., et al. 2021, The Astrophysical Journal Letters, 910, L13, doi: 10.3847/2041-8213/abe4de Event Horizon Telescope Collaboration, Akiyama, K.,

    doi:10.3847/2041-8213/abe4de 2021

Show all 27 references
  1. [9]

    2022, 930, L16, doi: 10.3847/2041-8213/ac6672

    Alberdi, A., et al. 2022, 930, L16, doi: 10.3847/2041-8213/ac6672

    2022 doi

  2. [10]

    W., Broderick, A

    Georgiev, B., Pesce, D. W., Broderick, A. E., et al. 2022, The Astrophysical Journal Letters, 930, L20, doi: 10.3847/2041-8213/ac65eb

    2022 doi

  3. [11]

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

    Johnson, M. D., Akiyama, K., Blackburn, L., et al. 2023, Galaxies, 11, 61, doi: 10.3390/galaxies11030061

    2023 doi

  4. [12]

    D., Akiyama, K., Baturin, R., et al

    Johnson, M. D., Akiyama, K., Baturin, R., et al. 2024, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 13092, Space Telescopes and Instrumentation 2024: Optical, Infrared, and Millimeter Wave, ed. L. E. Coyle, S. Matsuura, & M. D. Perrin, 1...

    2024 doi

  5. [13]

    K., Gammie, C

    Leung, P. K., Gammie, C. F., & Noble, S. C. 2011, ApJ, 737, 21, doi: 10.1088/0004-637X/737/1/21

    2011 doi

  6. [14]

    2022, Universe, 8, doi: 10.3390/universe8020085 Mo´ scibrodzka, M., & Gammie, C

    Mizuno, Y. 2022, Universe, 8, doi: 10.3390/universe8020085 Mo´ scibrodzka, M., & Gammie, C. F. 2018, MNRAS, 475, 43, doi: 10.1093/mnras/stx3162

    2022 doi

  7. [15]

    A., & Yfantis, A

    Moscibrodzka, M. A., & Yfantis, A. I. 2023, The Astrophysical Journal Supplement Series, 265, 22, doi: 10.3847/1538-4365/acb6f9 Mo´ scibrodzka, M., Falcke, H., & Shiokawa, H. 2016, Astronomy & Astrophysics, 586, A38, doi: 10.1051/0004-6361/201526630 Naethe Motta, P., mraia...

    2023 doi arXiv

  8. [16]

    2025, GRMHD File for EHT Polarized Code

    Prather, B. 2025, GRMHD File for EHT Polarized Code

    2025

  9. [17]

    Comparison, V1, Harvard Dataverse, doi: 10.7910/DVN/XZECPF

    doi

  10. [18]

    2021, The Journal of Open Source Software, 6, 3336, doi: 10.21105/joss.03336

    Prather, B., Wong, G., Dhruv, V., et al. 2021, The Journal of Open Source Software, 6, 3336, doi: 10.21105/joss.03336

    2021 doi

  11. [19]

    S., Dexter, J., Moscibrodzka, M., et al

    Prather, B. S., Dexter, J., Moscibrodzka, M., et al. 2023, The Astrophysical Journal, 950, 35, doi: 10.3847/1538-4357/acc586

    2023 doi

  12. [20]

    2016, ApJ, 820, 105, doi: 10.3847/0004-637X/820/2/105

    Pu, H.-Y., Yun, K., Younsi, Z., & Yoon, S.-J. 2016, ApJ, 820, 105, doi: 10.3847/0004-637X/820/2/105

    2016 doi

  13. [21]

    2016, arXiv:1607.07892 [cs.MS]

    Revels, J., Lubin, M., & Papamarkou, T. 2016, arXiv:1607.07892 [cs.MS]. https://arxiv.org/abs/1607.07892

    2016 arXiv

  14. [22]

    N., et al

    Sharma, A., Medeiros, L., Wong, G. N., et al. 2025, Astrophys. J., 985, 40, doi: 10.3847/1538-4357/adc104

    2025 doi

  15. [23]

    Tchekhovskoy, A., Narayan, R., & McKinney, J. C. 2011, Monthly Notices of the Royal Astronomical Society: Letters, 418, L79–L83, doi: 10.1111/j.1745-3933.2011.01147.x

    2011 doi

  16. [24]

    2022, Journal of Open Source Software, 7, 4457, doi: 10.21105/joss.04457

    Tiede, P. 2022, Journal of Open Source Software, 7, 4457, doi: 10.21105/joss.04457

    2022 doi

  17. [25]

    2018, Automatic differentiation of ODE integration

    Willkomm, J. 2018, Automatic differentiation of ODE integration. https://arxiv.org/abs/1802.02247

    2018 arXiv

  18. [26]

    N., et al

    Wong, G. N., et al. 2022, Astrophys. J. Supp., 259, 64, doi: 10.3847/1538-4365/ac582e

    2022 doi

  19. [27]

    B., et al

    Yoon, D., Chatterjee, K., Markoff, S. B., et al. 2020, MNRAS, 499, 3178, doi: 10.1093/mnras/staa3031

    2020 doi

Pith tools

Reviewed July 12, 2026 · model on record in the stance chip above.