General Relativistic Magnetohydrodynamic Simulations Research Articles

The impact of resistivity on the variability of black hole accretion flows

The accretion of magnetized plasma onto black holes is a complex and dynamic process in which the magnetic field plays a crucial role. The amount of magnetic flux that is accumulated near the event horizon significantly impacts the accretion flow behavior. Resistivity, which is a measure of how easily magnetic fields can dissipate, is thought to be a key factor influencing this process. This work explores the influence of resistivity on the accretion flow variability. We investigated simulations that reached the limit of the magnetically arrested disk (MAD) and simulations with an initial multi-loop magnetic field configuration. We employed 3D resistive general relativistic magnetohydrodynamic (MHD) simulations to model the accretion process under various regimes, where resistivity is globally constant (uniform resistivity). Our findings reveal distinct flow behaviors depending on resistivity. High-resistivity simulations never achieved the MAD state, which indicates a disturbed magnetic-flux accumulation process. Conversely, low-resistivity simulations converged toward the ideal MHD limit. The key results are that i) for the standard MAD model, resistivity plays a minimum role in flow variability, suggesting that flux eruption events dominate the dynamics. ii) High-resistivity simulations exhibit strong magnetic field diffusion into the disk that rearranges the efficient magnetic flux accumulation from the accretion flow. iii) In multi-loop simulations, resistivity significantly reduces the flow variability, which was not expected. However, magnetic flux accumulation becomes more variable as a result of frequent reconnection events at very low resistivity values. This study shows that resistivity affects how much the flow is distorted as a result of the magnetic field dissipation. Our findings provide new insights into the interplay between magnetic field accumulation, resistivity, variability, and the dynamics of black hole accretion.

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a “multizone” technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH (a * = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r −1 inside the Bondi radius R B , which is located at R B = 2 × 105 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ∼ unity, indicating an extended magnetically arrested state; the mass accretion rate Ṁ is ≈1% of the analytical Bondi accretion rate ṀB; and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate ≈0.02Ṁc2 . Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ∼ 102−107 r g, and find that Ṁ/ṀB≈(RB/6rg)−0.5 .

GR-RMHD Simulations of Super-Eddington Accretion Flows onto a Neutron Star with Dipole and Quadrupole Magnetic Fields

Although ultraluminous X-ray pulsars (ULXPs) are believed to be powered by super-Eddington accretion onto a magnetized neutron star (NS), the detailed structures of the inflow–outflow and magnetic fields are still not well understood. We perform general relativistic radiation magnetohydrodynamics (GR-RMHD) simulations of super-Eddington accretion flows onto a magnetized NS with dipole and/or quadrupole magnetic fields. Our results show that an accretion disk and optically thick outflows form outside the magnetospheric radius, while inflows aligned with magnetic field lines appear inside. When the dipole field is more prominent than the quadrupole field at the magnetospheric radius, accretion columns form near the magnetic poles, whereas a quadrupole magnetic field stronger than the dipole field results in the formation of a belt-like accretion flow near the equatorial plane. The NS spins up as the angular momentum of the accreting gas is converted into the angular momentum of the electromagnetic field, which then flows into the NS. Even if an accretion column forms near one of the magnetic poles, the observed luminosity is almost the same on both the side with the accretion column and the side without it, because the radiation energy is transported to both sides through scattering. Our model suggests that galactic ULXP Swift J0243.6+6124 has a quadrupole magnetic field of 2 × 1013 G and a dipole magnetic field of less than 4 × 1012 G.

Broadband Spectral Modeling of the M87 Nucleus

We study spectra produced by weakly accreting black hole (BH) systems using the semianalytic advection-dominated accretion flow (ADAF) model and the general-relativistic magnetohydrodynamic (GRMHD) simulation. We find significant differences between these two approaches related to a wider spread of the flow parameters as well as a much steeper radial distribution of the magnetic field in the latter. We apply these spectral models to the broadband spectral energy distribution (SED) of the nucleus of the M87 galaxy. The standard (in particular, 1D) formulation of the ADAF model does not allow us to explain it; previous claims that this model reproduces the observed SED suffer from an inaccurate treatment of the Compton process. The spectra based on GRMHD simulation are in a much better agreement with the observed data. In our GRMHD model, in which we assumed the BH spin a = 0.9, the bulk of radiation observed between the millimeter and the X-ray range is produced in the disk area within 4 gravitational radii from the BH. In this solution, the synchrotron component easily reproduces the spectral data between the millimeter and the UV range, and the Compton component does not violate the X-rays constraints, for Ṁ≲0.01M⊙ yr−1 and a relatively strong magnetic field, with the plasma β ∼ 1 in the region where radiation is produced. However, the Compton component cannot explain the observed X-ray spectrum. Instead, the X-ray spectrum can be reproduced by a high-energy tail of the synchrotron spectrum if electrons have a hybrid energy distribution with a ∼5% nonthermal component.

Testing Bayesian inference of GRMHD model parameters from VLBI data

ABSTRACT Recent observations by the Event Horizon Telescope (EHT) of supermassive black holes M87* and Sgr A* offer valuable insights into their space–time properties and astrophysical conditions. Utilizing a library of model images ($\sim 2$ million for Sgr A*) generated from general-relativistic magnetohydrodynamic (GRMHD) simulations, limited and coarse insights on key parameters such as black hole spin, magnetic flux, inclination angle, and electron temperature were gained. The image orientation and black hole mass estimates were obtained via a scoring and an approximate rescaling procedure. Lifting such approximations, probing the space of parameters continuously, and extending the parameter space of theoretical models is both desirable and computationally prohibitive with existing methods. To address this, we introduce a new Bayesian scheme that adaptively explores the parameter space of ray-traced, GRMHD models. The general relativistic radiative transfer code IPOLE is integrated with the EHT parameter estimation tool THEMIS. The pipeline produces a ray-traced model image from GRMHD data, computes predictions for very long baseline interferometric (VLBI) observables from the image for a specific VLBI array configuration and compares to data, thereby sampling the likelihood surface via a Markov chain Monte Carlo scheme. At this stage we focus on four parameters: accretion rate, electron thermodynamics, inclination, and source position angle. Our scheme faithfully recovers parameters from simulated VLBI data and accommodates time-variability via an inflated error budget. We highlight the impact of intrinsic variability on model fitting approaches. This work facilitates more informed inferences from GRMHD simulations and enables expansion of the model parameter space in a statistically robust and computationally efficient manner.

Mildly Super-Eddington Accretion onto Slowly Spinning Black Holes Explains the X-Ray Weakness of the Little Red Dots

JWST has revealed a population of low-luminosity active galactic nuclei at z > 4 in compact, red hosts (the “Little Red Dots,” or LRDs), which are largely undetected in X-rays. We investigate this phenomenon using General Relativistic Radiation Magnetohydrodynamics simulations of super-Eddington accretion onto a supermassive black hole (SMBH) with M • = 107 M ⊙ at z ∼ 6, representing the median population; the spectral energy distributions (SEDs) that we obtain are intrinsically X-ray weak. The highest levels of X-ray weakness occur in SMBHs accreting at mildly super-Eddington rates (1.4 < f Edd < 4) with zero spin, viewed at angles >30° from the pole. X-ray bolometric corrections in the observed 2–10 keV band reach ∼104 at z = 6, ∼5 times higher than the highest constraint from X-ray stacking. Most SEDs are extraordinarily steep and soft in the X-rays (median photon index Γ = 3.1, mode of Γ = 4.4). SEDs strong in the X-rays have harder spectra with a high-energy bump when viewed near the hot (>108 K) and highly relativistic jet, whereas X-ray weak SEDs lack this feature. Viewing an SMBH within 10° of its pole, where beaming enhances the X-ray emission, has a ∼1.5% probability, matching the LRD X-ray detection rate. Next-generation observatories like AXIS will detect X-ray-weak LRDs at z ∼ 6 from any viewing angle. Although many SMBHs in the LRDs are already estimated to accrete at super-Eddington rates, our model explains 50% of their population by requiring that their masses are overestimated by a mere factor of ∼3. In summary, we suggest that LRDs host slowly spinning SMBHs accreting at mildly super-Eddington rates, with large covering factors and broad emission lines enhanced by strong winds, providing a self-consistent explanation for their X-ray weakness and complementing other models.

Measuring Black Hole Light Echoes with Very Long Baseline Interferometry

Light passing near a black hole can follow multiple paths from an emission source to an observer due to strong gravitational lensing. Photons following different paths take different amounts of time to reach the observer, which produces an echo signature in the image. The characteristic echo delay is determined primarily by the mass of the black hole, but it is also influenced by the black hole spin and inclination to the observer. In the Kerr geometry, echo images are demagnified, rotated, and sheared copies of the direct image and lie within a restricted region of the image. Echo images have exponentially suppressed flux, and temporal correlations within the flow make it challenging to directly detect light echoes from the total light curve. In this Letter, we propose a novel method to search for light echoes by correlating the total light curve with the interferometric signal at high spatial frequencies, which is a proxy for indirect emission. We explore the viability of our method using numerical general relativistic magnetohydrodynamic simulations of a near-face-on accretion system scaled to M87-like parameters. We demonstrate that our method can be used to directly infer the echo delay period in simulated data. An echo detection would be clear evidence that we have captured photons that have circled the black hole, and a high-fidelity echo measurement would provide an independent measure of fundamental black hole parameters. Our results suggest that detecting echoes may be achievable through interferometric observations with a modest space-based very long baseline interferometry mission.

She’s Got Her Mother’s Hair: Unveiling the Origin of Black Hole Magnetic Fields through Stellar to Collapsar Simulations

Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a low net magnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so.

Chaotic magnetic disconnections trigger flux eruptions in accretion flows channeled onto magnetically saturated Kerr black holes

Magnetized accretion flow onto a black hole (BH) may lead to the accumulation of poloidal magnetic flux across its horizon, which for high BH spin can power far-reaching relativistic jets. The BH magnetic flux is subject to a saturation mechanism by means of magnetic flux eruptions involving relativistic magnetic reconnection. Such accretion flows have been described as magnetically arrested disks (MAD) or magnetically choked accretion flows (MCAF). The main goal of this work is to describe the onset of relativistic reconnection and initial development of magnetic flux eruption in accretion flow onto magnetically saturated BHs. We analyzed the results of 3D general relativistic ideal magnetohydrodynamic (GRMHD) numerical simulations in the Kerr metric, starting from weakly magnetized geometrically thick tori rotating either prograde or retrograde. We integrate d large samples of magnetic field lines in order to probe magnetic connectivity with the BH horizon. The boundary between magnetically connected and disconnected domains coincides roughly with enthalpy equipartition. The geometrically constricted innermost part of the disconnected domain develops a rigid structure of magnetic field lines -- rotating slowly and insensitive to the BH spin orientation. The typical shape of innermost disconnected lines is a double spiral converging to a sharp inner tip anchored at the single equatorial current layer. The foot points of magnetic flux eruptions are found to zip around the BH along with other azimuthal patterns. Magnetic flux eruptions from magnetically saturated accreting BHs can be triggered by minor density gaps in the disconnected domain, resulting from the chaotic disconnection of plasma-depleted magnetospheric lines. Accretion flow is effectively channeled along the disconnected lines toward the current layer, and further toward the BH by turbulent cross-field diffusion. Rotation of flux eruption foot points may contribute to the variability of BH crescent images.

Constraints on the accretion properties of quasi-periodic erupters from GRMHD simulations

Some apparently quiescent supermassive black holes (BHs) at centers of galaxies show quasi-periodic eruptions (QPEs) in the X-ray band, the nature of which is still unknown. A possible origin for the eruptions is an accretion disk. However, the properties of such disks are restricted by the timescales of recurrence and the duration of the flares. In this work, we test the possibility that the temporal properties of known QPEs can be explained by accretion from a compact accretion disk with an outer radius out g and we focus on a particular object, GSN 069. We ran several 3D general relativistic magnetohydrodynamic (GRMHD) simulations with the H-AMR code of thin and thick disks and studied how the initial disk parameters such as thickness, magnetic field configuration, magnetization, and Kerr parameter affect the observational properties of QPEs. We show that accretion onto a slowly rotating BH through a small, moderately thin accretion disk with an initially low plasma beta can explain the observed time between outbursts and the lack of evidence for a variable jet emission. In order to form such a disk, the accreting matter should have a low net angular momentum. A potential source for such low angular momentum matter with a quasi-periodic feeding mechanism might be a tight binary of wind-launching stars. Apart from their primary application, our results can also be useful for general studies of systems with small accretion disks, in which evolution occurs very rapidly so that the disks cannot be considered stationary. For such systems, it is important to understand how the initial conditions affect the results.

Energy-dependent and Energy-integrated Two-moment General-relativistic Neutrino Transport Simulations of a Hypermassive Neutron Star

We compare two-moment-based energy-dependent and three variants of energy-integrated neutrino transport general-relativistic magnetohydrodynamics simulations of a hypermassive neutron star. To study the impacts due to the choice of the neutrino transport schemes, we perform simulations with the same setups and input neutrino microphysics. We show that the main differences between energy-dependent and energy-integrated neutrino transport are found in the disk and ejecta properties, as well as in the neutrino signals. The properties of the disk surrounding the neutron star and the ejecta in energy-dependent transport are very different from the ones obtained using energy-integrated schemes. Specifically, in the energy-dependent case, the disk is more neutron-rich at early times and becomes geometrically thicker at later times. In addition, the ejecta is more massive and, on average, more neutron-rich in the energy-dependent simulations. Moreover, the average neutrino energies and luminosities are about 30% higher. Energy-dependent neutrino transport is necessary if one wants to better model the neutrino signals and matter outflows from neutron star merger remnants via numerical simulations.

Accreting Neutron Stars in 3D General-relativistic Magnetohydrodynamic Simulations: Jets, Magnetic Polarity, and the Interchange Slingshot

Abstract Accreting neutron stars differ from black holes by the presence of the star’s own magnetic field, whose interaction with the accretion flow is a central component in understanding these systems’ disk structure, outflows, jets, and spin evolution. It also introduces an additional degree of freedom, as the stellar dipole can have any orientation relative to the inner disk’s magnetic field. We present a suite of 3D general-relativistic magnetohydrodynamic simulations in which we investigate the two extreme polarities, with the dipole field being either parallel or antiparallel to the initial disk field, in both the accreting and propeller states. When the magnetosphere truncates the disk near or beyond the corotation radius, most of the system’s properties, including the relativistic jet power, are independent of the star–disk relative polarity. However, when the disk extends well inside the corotation radius, in the parallel orientation the jet power is suppressed and the inner disk is less dense and more strongly magnetized. We suggest a physical mechanism that may account for this behavior—the interchange slingshot—and discuss its astrophysical implications. When the star is in the rapidly accreting regime, which in most cases will be associated with strong spin-up, we expect large observational differences between the two magnetic orientations. This may be reflected in increased variability as the accretion flow drags in successive magnetic structures of varying polarity.

General relativistic magnetohydrodynamic simulations with bam: Implementation and code comparison

General relativistic magnetohydrodynamic simulations with bam: Implementation and code comparison

Bayesian Black Hole Photogrammetry

We propose an analytic dual-cone accretion model for horizon-scale images of the cores of low-luminosity active galactic nuclei, including those observed by the Event Horizon Telescope (EHT). Our model is of synchrotron emission from an axisymmetric, magnetized plasma, constrained to flow within two oppositely oriented cones that are aligned with the black hole’s spin axis. We show this model can accurately reproduce images of a variety of time-averaged general relativistic magnetohydrodynamic simulations and that it accurately recovers the black hole spin, orientation, emission scale height, peak emission radius, and fluid flow direction from these simulations within a Bayesian inference framework using radio interferometric data. We show that nontrivial topologies in the images of relativistic accretion flows around black holes can result in nontrivial multimodal solutions when applied to observations with a sparse array, such as the EHT 2017 observations of M87*. The presence of these degeneracies underscores the importance of employing Bayesian techniques to adequately sample the posterior space for the interpretation of EHT measurements. We fit our model to the EHT observations of M87* and find a 95% highest posterior density interval for the mass-to-distance ratio of θ g ∈ (2.84, 3.75) μas, and give an inclination of θ o ∈ (11°, 24°). These new measurements are consistent with mass measurements from the EHT and stellar dynamical estimates and with the spin axis inclination inferred from properties of the M87* jet.

Shock-induced Partial Alignment in Geometrically Thick Tilted Accretion Disks Around Black Holes

We carry out idealized three-dimensional general-relativistic magnetohydrodynamic simulations of prograde, weakly magnetized, and geometrically thick accretion flows where the gas distribution is misaligned from the black hole (BH) spin axis. We evolve the disk for three BH spins: a = 0.5, 0.75, and 0.9375, and we contrast them with a standard aligned disk simulation with a = 0.9375. The tilted disks achieve a warped and twisted steady-state structure, with the outer disk misaligning further away from the BH and surpassing the initial 24° misalignment. However, closer to the BH, there is evidence of partial alignment, as the inclination angle decreases with radius in this regime. Standing shocks also emerged in proximity to the BH, roughly at ∼6 gravitational radii. We show that these shocks act to partially align the inner disk with the BH spin. The rate of alignment increases with increasing BH spin magnitude, but in all cases is insufficient to fully align the gas before it accretes. Additionally, we present a toy model of orbit crowding that can predict the location of the shocks in moderate-to-fast rotating BHs, illustrating a potential physical origin for the behavior seen in simulations—with possible applications in determining the positions of shocks in real misaligned astrophysical systems.

General Relativistic Radiation Magnetohydrodynamics Simulations of Precessing Tilted Super-Eddington Disks

We perform a three-dimensional general relativistic radiation magnetohydrodynamics simulation of a tilted super-Eddington accretion disk around a spinning black hole (BH). The disk, which tilts and twists as it approaches the BH, precesses while maintaining its shape. The gas is mainly ejected around the rotation axis of the outer part of the disk rather than around the spin axis of the BH. The disk precession changes the ejection direction of the gas with time. The radiation energy is also released in approximately the same direction as the outflow, so the precession is expected to cause a quasiperiodic time variation of the observed luminosity. The timescale of the precession is about 10 s for a 10 M ⊙ BH and for the radial extent of a disk of several tens of gravitational radii, where M ⊙ is the solar mass. This timescale is consistent with the frequency of the low-frequency quasiperiodic oscillation (0.01–1 Hz) observed in some ultraluminous X-ray sources.

GRMHD Simulations of Accretion Structures with Different Angular Momentum Profiles

In this article, we explore the dynamics of accretion structures encircling spherically symmetric black holes, comparing three accretion disk models with distinct angular momentum profiles: (i) the geometrically thin Keplerian disk, (ii) the Fishbone–Moncrief torus; and (iii) the Polish Doughnut. Employing general relativistic magnetohydrodynamics simulations with the High Accuracy Relativistic Magnetohydrodynamics code, we investigate these three models, considering the magnetic field’s influence on the accretion disk angular momentum redistribution. We show that the magnetic field is a key factor in accretion disk structures, especially in regions with lower mass density. Our investigation verifies the well-established fact that the presence of a magnetic field significantly influences the accretion rate and its temporal variability.

In LIGO’s Sight? Vigorous Coherent Gravitational Waves from Cooled Collapsar Disks

We present the first numerical study of gravitational waves (GWs) from collapsar disks, using state-of-the-art 3D general relativistic magnetohydrodynamic simulations of collapsing stars. These simulations incorporate a fixed Kerr metric for the central black hole (BH) and employ simplified prescriptions for disk cooling. We find that cooled disks with an expected scale height ratio of H/R ≳ 0.1 at ∼10 gravitational radii induce Rossby instability in compact, high-density rings. The trapped Rossby vortices generate vigorous coherent emission regardless of disk magnetization and BH spin. For BH mass of ∼10 M ⊙, the GW spectrum peaks at ∼100 Hz, with some breadth due to various nonaxisymmetric modes. The spectrum shifts toward lower frequencies as the disk viscously spreads and the circularization radius of the infalling gas increases. Weaker-cooled disks with H/R ≳ 0.3 form a low-density extended structure of spiral arms, resulting in a broader, lower-amplitude spectrum. Assuming an optimistic detection threshold with a matched-filter signal-to-noise ratio of 20 and a rate similar to Type Ib/c supernovae, LIGO–Virgo–KAGRA (LVK) could detect ≲1 event annually, suggesting that GW events may already be hidden in observed data. Third-generation GW detectors could detect dozens to hundreds of collapsar disks annually, depending on the cooling strength and the disk formation rate. The GW amplitudes from collapsar disks are ≳100 times higher with a substantially greater event rate than those from core-collapse supernovae, making them potentially the most promising burst-type GW class for LVK and Cosmic Explorer. This highlights the importance of further exploration and modeling of disk-powered GWs, promising insights into collapsing star physics.

Dynamics and emission properties of flux ropes from two-temperature GRMHD simulations with multiple magnetic loops

Context. Flux ropes erupting in the vicinity of a black hole are thought to be a potential model for the flares observed in Sagittarius A*. Aims. In this study, we examine the radiative properties of flux ropes that emerged in the vicinity of a black hole. Methods. We performed three-dimensional two-temperature general relativistic magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows with alternating multiple magnetic loops and general relativistic radiation transfer (GRRT) calculations. In the GRMHD simulations, we implemented two different sizes of initial magnetic loops. Results. In the small loop case, magnetic dissipation leads to a weaker excitement of magneto-rotational instability inside the torus, which generates a lower accretion rate compared to the large loop case. However, the small loop case generates more flux ropes due to frequent reconnection by magnetic loops with different polarities. By calculating the thermal synchrotron emission, we found that the variability of light curves and the emitting region are tightly related. At 230 GHz and higher, the emission from the flux ropes is relatively stronger compared to the background, which is responsible for the filamentary structure in the images. At lower frequencies (e.g., 43 GHz), emission comes from more extended regions, which have a less filamentary structure in the image. Conclusions. Our study shows that self-consistent electron temperature models are essential for the calculation of thermal synchrotron radiation and the morphology of the GRRT images. Flux ropes contribute considerable emission at frequencies ≳230 GHz.

Bayesian Constraints on the Ring Ellipticity of M87* 2017 Using Themis

Measuring the properties of black hole images has the potential to constrain deviations from general relativity on horizon scales. Of particular interest is the ellipticity of the ring that is sensitive to the underlying spacetime. In 2019, the Event Horizon Telescope (EHT) produced the first-ever image of a black hole on horizon scales. Here, we reanalyze the M87* EHT 2017 data using Bayesian imaging (BI) techniques, constructing a posterior of the ring shape. We find that BI recovers the true on-sky ring shape more reliably than the original imaging methods used in 2019. As a result, we find that M87*'s ring ellipticity is 0.09−0.06+0.07 and is consistent with the measured ellipticity from general relativistic magnetohydrodynamic simulations.