⚛️ general relativity

Semi-analytic studies of accretion disk and magnetic field geometry in M87*

Using a semi-analytic radiatively inefficient accretion flow model combined with general relativistic ray tracing, this study demonstrates that M87* is most consistent with a poloidal magnetic field-dominated flow featuring radial inflow, a configuration that can be reliably distinguished from toroidal-dominated alternatives by Event Horizon Telescope observables.

Original authors: Saurabh, Maciek Wielgus, Arman Tursunov, Andrei P. Lobanov, Razieh Emami

Published 2026-02-04

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Saurabh, Maciek Wielgus, Arman Tursunov, Andrei P. Lobanov, Razieh Emami

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the center of the galaxy M87 as a cosmic whirlpool, a supermassive black hole swallowing everything around it. For a long time, we could only guess what the "plumbing" of this whirlpool looked like—specifically, how the invisible magnetic fields and swirling gas (plasma) were arranged.

This paper is like a team of cosmic detectives building a series of virtual simulations to figure out which "plumbing setup" matches the actual photos taken by the Event Horizon Telescope (EHT). They didn't just look at the black hole's shadow; they looked at the glowing ring of light around it and the invisible magnetic "threads" woven through that light.

Here is a breakdown of their investigation using simple analogies:

1. The Setup: Building Virtual Black Holes

The researchers built a "toy model" of the black hole's accretion disk (the swirling disk of hot gas). Instead of running a super-complex, slow-motion video game simulation for every single possibility, they used a semi-analytic model. Think of this as using a precise mathematical recipe to instantly generate thousands of different versions of the black hole, rather than waiting for a computer to simulate every drop of gas moving.

They changed four main ingredients in their recipe:

The Spin: How fast the black hole is spinning (like a top).
The Flow: How the gas moves (is it swirling in neat circles like a carousel, or falling straight in like a waterfall?).
The Thickness: Is the disk a thin, flat pancake or a puffy, thick donut?
The Magnetic Fields: This was the big variable. They tested different shapes for the magnetic fields, such as:
- Toroidal: Like rubber bands wrapped around a ball.
- Poloidal: Like the lines on a tennis ball or a globe, running from pole to pole.
- Dipole/Quadrupole: More complex patterns, like a bar magnet or a four-pole magnet.

2. The Experiment: Taking "Virtual Photos"

Once they built these virtual black holes, they used a technique called "ray-tracing." Imagine shooting millions of laser beams from the camera's eye, through the virtual black hole, and seeing how the light bends and changes color due to gravity and magnetism.

They then took these virtual photos and compared them to the real photos of M87* taken by the EHT. They looked for specific clues:

The Ring Size: How big is the glowing circle?
The Brightness: Is one side of the ring brighter than the other? (This happens because the gas moving toward us looks brighter, like a car headlight).
The Polarization: This is the "direction" of the light waves. It acts like a fingerprint for the magnetic fields. If the magnetic field is a rubber band, the light waves line up one way; if it's a globe, they line up another.

3. The Findings: What Fits and What Doesn't

The Magnetic Field Mystery
The most important discovery was about the magnetic fields.

The "Rubber Band" (Toroidal) vs. "Globe" (Poloidal): The team found they could clearly tell the difference between a magnetic field that wraps around the disk (toroidal) and one that runs through it (poloidal).
The Winner: The real photos of M87* look most like a model where the magnetic field is poloidal (running through the disk like a globe), mixed with some swirling motion. A purely "rubber band" style field didn't match the photos.

The Spin and the Flow

The Spin: The black hole is likely spinning in a "positive" direction (prograde), meaning the gas is swirling in the same direction the black hole spins. While they couldn't pinpoint the exact speed, they ruled out a slow or backward spin.
The Flow: The gas isn't just spinning in perfect circles. It's also falling inward. The models where the gas falls straight in (radial inflow) matched the real photos better than models where the gas just orbits perfectly.

The Thickness

Pancake vs. Donut: They tested if the disk was a thin pancake or a thick donut. Surprisingly, it didn't matter much. Whether the disk was thin or thick, the resulting "photo" looked very similar. This suggests that for the purpose of understanding the light we see, we can treat the disk as if it were flat without losing much accuracy.

The "Faraday" Problem
There was one snag. The real photos show that the light is somewhat "scrambled" (depolarized) as it passes through the gas, like looking through fog. The researchers' simple models were too "clear" (they didn't have enough "fog"). This suggests that the real black hole might have a more chaotic, turbulent structure or a jet of gas in front of the disk that is scrambling the light, which their simple models didn't fully capture.

4. The Conclusion

By comparing their virtual models to the real universe, the team concluded:

M87* is best described as a black hole with a poloidal magnetic field (like a globe) and gas that is falling inward while spinning.
The black hole is likely spinning moderately to fast in the forward direction.
The "thickness" of the gas disk is less important than we thought for these specific observations.

In short, the paper used a clever mix of math and computer simulations to narrow down the "personality" of the M87* black hole, telling us that its magnetic fields act more like the lines on a globe than rubber bands, and that the gas around it is in a hurry to fall in.

Technical Summary: Semi-analytic studies of the accretion disk and magnetic field geometry in M 87*

Problem and Context
The behavior of magnetic fields near black holes (BHs) and their role in accretion flows and jet formation remain central questions in BH astrophysics. While Event Horizon Telescope (EHT) observations of M 87* and Sagittarius A* have provided unprecedented insights into horizon-scale structures, interpreting these observations requires a theoretical framework linking polarized emission to system properties. General Relativistic Magnetohydrodynamics (GRMHD) simulations are the standard tool for this but are computationally expensive and rely on uncertain initial conditions and ideal MHD assumptions that may not fully capture collisionless plasma physics. Conversely, simple toy models lack the physical consistency to make quantitative predictions. This paper addresses the gap by employing a stationary, semi-analytic Radiatively Inefficient Accretion Flow (RIAF) model to systematically investigate how magnetic field geometry, black hole spin, plasma dynamics, and disk thickness affect observable synchrotron emission in M 87*.

Methodology
The authors constructed a stationary, axisymmetric RIAF model in Kerr spacetime using the open-source ray-tracing code ipole. The model assumes emission originates near the equatorial plane, neglecting the jet base contribution for the primary ring morphology.

Accretion Flow: The electron number density ( $n_e$ ) and temperature ( $T_e$ ) follow power laws of the spherical radius. The four-velocity of the flow is interpolated between Keplerian orbital motion and geodesic free-fall, parameterized by $\kappa_K$ (Keplerian) and $\kappa_{ff}$ (radial infall).
Magnetic Field Geometry: Seven distinct magnetic field configurations were tested: Radial (split monopole), Vertical, Dipole, Quadrupole, Parabolic, Combined (Wald solution + radial), and Toroidal. These were defined via vector potentials to ensure a divergence-free field.
Parameters: The study varied the dimensionless BH spin ( $a^*$ ) from -0.94 to +0.94, the geometric disk thickness ( $H$ ), and the flow dynamics. The fiducial model used $H=0.5$ and a mixed velocity profile ( $\kappa_{ff}=0.5, \kappa_K=0.5$ ).
Synthetic Observables: The authors generated synthetic images at 230 GHz and extracted EHT-constrained metrics: ring angular diameter ( $\alpha$ ), brightness asymmetry ( $I_r$ ), net fractional polarization ( $|m|_{net}$ ), resolved fractional polarization ( $\langle|m|\rangle$ ), and the second mode of the polarized image Fourier decomposition ( $\beta_2$ ).
Comparison: These semi-analytic results were compared against EHT observational constraints and a library of time-averaged GRMHD simulations (MAD and SANE states).

Key Results

Magnetic Field Geometry: The study confirms that magnetic field geometry strongly dictates the polarized image structure. Toroidal fields produce radial Electric Vector Position Angle (EVPA) patterns, while poloidal fields produce azimuthal patterns. The argument of the $\beta_2$ parameter ( $\arg \beta_2$ ) serves as a robust discriminator: it is near $0^\circ$ for azimuthal (toroidal) dominance and near $180^\circ$ for poloidal dominance. The EHT data, with $\arg \beta_2$ in the intermediate range ( $-163^\circ$ to $-127^\circ$ ), suggests a mixed or poloidal-dominated geometry rather than a purely toroidal one.
Plasma Dynamics and Asymmetry: The flow dynamics significantly impact brightness asymmetry ( $I_r$ ). Keplerian models (pure azimuthal motion) produce high asymmetry driven by Doppler boosting, often exceeding EHT constraints. In contrast, models with significant radial inflow (zero angular momentum) are too symmetric. The observed asymmetry in M 87* is best reproduced by models with sub-Keplerian velocities and moderate to high positive BH spins.
Black Hole Spin: While spin has a secondary effect on ring diameter, it is a primary driver of brightness asymmetry. Positive spins ( $a^* > 0$ ) align Doppler boosting and relativistic effects to enhance asymmetry, whereas negative spins often result in images that are too symmetric to match observations. The data favors moderate to high positive spin values.
Disk Thickness: The geometric thickness of the disk ( $H$ ) has a limited effect on most observables, validating the use of equatorial emission assumptions for global disk-dominated models, except in specific cases like dipole fields where off-equatorial regions contribute significantly to polarimetric signatures.
Polarization Discrepancy: A notable finding is that the semi-analytic RIAF models remain "Faraday-thin," resulting in high resolved polarization ( $\langle|m|\rangle$ ) and high $\beta_2$ amplitudes that are inconsistent with EHT observations. In contrast, GRMHD simulations (particularly with high ion-to-electron temperature ratios, $R_h$ ) exhibit significant Faraday depolarization. The authors suggest that the observed depolarization in M 87* likely arises from non-axisymmetric turbulent structures or off-equatorial jet base plasma, effects not captured in the stationary semi-analytic framework.

Significance and Conclusions
The paper demonstrates that horizon-scale observables depend systematically and physically on magnetic field geometry, plasma velocity, and BH spin. The authors conclude that the M 87* system is most consistent with a flow dominated by poloidal magnetic fields with partially radial inflow. While the BH spin remains elusive, the data disfavors zero or mildly negative spins, preferring moderate to high positive values.

The study highlights the limitations of purely semi-analytic models in reproducing the full complexity of EHT data, specifically regarding Faraday depolarization, which suggests that turbulent, non-axisymmetric features or jet contributions are critical. The work serves as a bridge between intuitive toy models and computationally exhaustive GRMHD, offering a parametric framework to interpret polarimetric constraints. The authors suggest that future work should explore mixed toroidal-poloidal field configurations and non-Kerr spacetimes within the RIAF framework to further refine these constraints.

1. The Setup: Building Virtual Black Holes

2. The Experiment: Taking "Virtual Photos"

3. The Findings: What Fits and What Doesn't

4. The Conclusion

More like this