Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U, and V, respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

The Event Horizon Telescope observed the horizon-scale synchrotron emission region around the Galactic center supermassive black hole, Sagittarius A* (Sgr A*), in 2017. These observations revealed a bright, thick ring morphology with a diameter of 51.8 ± 2.3 μas and modest azimuthal brightness asymmetry, consistent with the expected appearance of a black hole with mass M ≈ 4 × 10⁶M_⊙. From these observations, we present the first resolved linear and circular polarimetric images of Sgr A*. The linear polarization images demonstrate that the emission ring is highly polarized, exhibiting a prominent spiral electric vector polarization angle pattern with a peak fractional polarization of ∼40% in the western portion of the ring. The circular polarization images feature a modestly (∼5%–10%) polarized dipole structure along the emission ring, with negative circular polarization in the western region and positive circular polarization in the eastern region, although our methods exhibit stronger disagreement than for linear polarization. We analyze the data using multiple independent imaging and modeling methods, each of which is validated using a standardized suite of synthetic data sets. While the detailed spatial distribution of the linear polarization along the ring remains uncertain owing to the intrinsic variability of the source, the spiraling polarization structure is robust to methodological choices. The degree and orientation of the linear polarization provide stringent constraints for the black hole and its surrounding magnetic fields, which we discuss in an accompanying publication.

Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here we present the first analysis of circular polarization using EHT data, acquired in 2017, which can potentially provide additional insights into the magnetic fields and plasma composition near the black hole. Interferometric closure quantities provide convincing evidence for the presence of circularly polarized emission on event-horizon scales. We produce images of the circular polarization using both traditional and newly developed methods. All methods find a moderate level of resolved circular polarization across the image (〈∣v∣〉 < 3.7%), consistent with the low image-integrated circular polarization fraction measured by the Atacama Large Millimeter/submillimeter Array (∣v_int∣ < 1%). Despite this broad agreement, the methods show substantial variation in the morphology of the circularly polarized emission, indicating that our conclusions are strongly dependent on the imaging assumptions because of the limited baseline coverage, uncertain telescope gain calibration, and weakly polarized signal. We include this upper limit in an updated comparison to general relativistic magnetohydrodynamic simulation models. This analysis reinforces the previously reported preference for magnetically arrested accretion flow models. We find that most simulations naturally produce a low level of circular polarization consistent with our upper limit and that Faraday conversion is likely the dominant production mechanism for circular polarization at 230 GHz in M87*.