A simple model for extracting astrophysics from black hole images

The Event Horizon Telescope (EHT) is providing unprecedented high-resolution images of supermassive black holes. These images are fundamentally related to properties of the luminous accretion disks, since black holes themselves produce no light. We develop a simple prescription to relate observational features of black hole images to a toy model for the intensity profile of the associated accretion disk. We apply our model to the original EHT image of M87*, as well as to the reanalyzed image from the PRIMO algorithm, providing generic, simultaneous constraints on the mass of the black hole and properties of the accretion disk emission. While current images lack the resolution to confidently detect the photon ring, the consideration of multiple image parameters are found to contain enough information to provide constraints on the inner edge of the accretion disk along with the black hole mass. Using observed features of the original EHT image, we constrain the mass of M87* to be $6.6^{+1.2}{-1.0}\times 10^9 M\odot$ to 68$%$ confidence, and find that emission may extend all the way to the black hole horizon. When instead using constraints from the PRIMO algorithm’s image along with constraints on the brightness asymmetry provided by the original EHT analysis, we find M87*’s mass to be $ 6.4^{+0.7}{-0.7}\times 10^9 M\odot$ to 68$%$ confidence, with the inner edge of the accretion disk between $3M$ and $5.3M$. Both analyses rule out an inner edge of the accretion disk coinciding with the innermost stable circular orbit for a Schwarzschild black hole. Furthermore, the narrow ring width reported in the PRIMO image also confidently rules out emission increasing all the way down to the black hole horizon. Further assumptions on the mass of M87* and connections between the accretion disk cutoff and physical radii allow for rudimentary black hole spin estimates.

💡 Research Summary

The paper presents a computationally inexpensive two‑parameter toy model that links observable features of Event Horizon Telescope (EHT) images of the supermassive black hole M87* to the underlying intensity profile of its accretion disk. The authors adopt the “just‑add‑one’’ approximation for photon trajectories in a Schwarzschild spacetime, which relates the impact parameter on the observer’s image plane to the emission radius in the disk with an error of ≤10 % for non‑spinning black holes and modestly larger errors for modest spins. Because the current EHT data lack the resolution to isolate the photon ring, the authors restrict themselves to a non‑spinning (Schwarzschild) geometry and demonstrate in an appendix that the ring diameter and width are largely insensitive to spin, justifying this simplification.

The accretion disk is modeled as geometrically thin, optically thin, and inclined by a fixed angle (θ₀≈17°). Its specific intensity depends only on radius and follows a broken power‑law with a sharp inner cutoff: I(r) ∝ r^β for r ≥ r_cut, and zero for r < r_cut, where β < 0. The cutoff radius r_cut = 2 M corresponds to emission reaching the event horizon; larger values imply the disk truncates before the horizon. This two‑parameter description (β, r_cut) is deliberately simple, yet it captures a wide range of plausible brightness distributions.

Ray tracing is performed by backward‑integrating null geodesics from the image plane (impact parameter b, azimuth α) to the disk, recording up to three intersections (direct image, lensing ring, photon ring). The observed specific intensity at each image pixel is the sum over intersections of the redshifted emissivity, I′(b,α)=∑_{i=1}^{3} g(b,α,r_i)^3 I(r_i), where g includes both gravitational redshift and, if a Keplerian velocity profile is assumed, the Doppler factor. The authors adopt a Keplerian velocity outside the ISCO (r=6 M) and a free‑fall continuation inside, noting that alternative velocity prescriptions would shift the quantitative results but not the qualitative conclusions.

To mimic the finite resolution of the 2017 EHT array, the generated images are convolved with a Gaussian kernel matching the nominal beam. The authors then extract three observable image parameters that were reported for the original EHT reconstruction and for the newer PRIMO reconstruction: (i) the ring diameter R, (ii) the ring width ΔR, and (iii) the brightness asymmetry A (ratio of peak intensities on opposite sides). These quantities are treated as Gaussian‑distributed measurements with known covariances.

A Bayesian inference framework is built with uniform priors on the black‑hole mass M, the power‑law index β, and the cutoff radius r_cut. The likelihood combines the three measured image parameters, and posterior distributions are sampled using Markov Chain Monte Carlo. Two separate analyses are performed: (a) using only the original EHT image constraints, and (b) using the PRIMO image constraints together with the original brightness‑asymmetry measurement.

Results from the original EHT image alone yield M = 6.6^{+1.2}_{‑1.0} × 10⁹ M⊙ (68 % confidence), with a broad allowed range for β (≈ ‑2 to ‑4) and r_cut that includes the horizon (r_cut≈2 M). Thus, with the current resolution, the data cannot exclude a disk that extends all the way to the event horizon.

When the PRIMO constraints are added, the posterior tightens: M = 6.4^{+0.7}_{‑0.7} × 10⁹ M⊙, and the inner edge is confined to 3 M ≤ r_cut ≤ 5.3 M. This range lies well inside the Schwarzschild ISCO (6 M), indicating that the emitting region does not reach the ISCO, contrary to many GRMHD‑based analyses that often assume emission down to the ISCO. Moreover, the narrow ring width reported by PRIMO rules out models with very negative β (i.e., strong central concentration of emission), implying that the intensity does not increase all the way to the horizon.

The authors discuss how, if an independent mass estimate is adopted, the mapping between r_cut and a physical radius (e.g., ISCO, photon sphere) can be used to infer a crude spin parameter, though they caution that spin remains poorly constrained because the brightness asymmetry is highly sensitive to the assumed velocity profile and to spin‑induced lensing effects.

Strengths of the work include: (1) a transparent analytic pipeline that avoids the computational expense of large GRMHD libraries; (2) a clear demonstration that multiple image observables together can break degeneracies between mass and disk structure; (3) an explicit quantification of how different image reconstructions (EHT vs PRIMO) affect astrophysical inference.

Limitations are also acknowledged: the Schwarzschild assumption neglects potentially important spin effects on the photon ring and on asymmetry; the sharp cutoff and single‑power‑law emissivity are idealizations that may not capture realistic temperature or magnetic field gradients; and the current data lack sufficient resolution to directly detect the photon ring, which would provide an independent geometric mass measurement.

The paper concludes that future higher‑resolution, multi‑frequency VLBI observations (e.g., next‑generation EHT) will enable detection of the photon ring and finer ring structure, allowing the toy model to be refined (e.g., multi‑segment power laws, smooth cutoffs) and extended to Kerr spacetimes. Such advances will permit simultaneous, model‑independent constraints on black‑hole mass, spin, inclination, and detailed accretion‑disk physics, moving beyond the present reliance on computationally intensive GRMHD simulations.