Abstract We present a set of eleven two-temperature, radiative, general relativistic magnetohydrodynamic (2TGRRMHD) simulations of the black hole M87* in the magnetically arrested (MAD) state, surveying different values of the black hole spin a*. Our 3D simulations self-consistently evolve the temperatures of separate electron and ion populations under the effects of adiabatic compression/expansion, viscous heating, Coulomb coupling, and synchrotron, bremsstrahlung, and inverse Compton radiation. We adopt a sub-grid heating prescription from gyrokinetic simulations of plasma turbulence. Our simulations have accretion rates $\dot{M}=(0.5-1.5)\times 10^{-6}\dot{M}_{\rm Edd}$ and radiative efficiencies εrad = 3 − 35 %. We compare our simulations to a fiducial set of otherwise identical single-fluid GRMHD simulations and find no significant changes in the outflow efficiency or black hole spindown parameter. Our simulations produce an effective adiabatic index for the two-temperature plasma of Γgas ≈ 1.55, larger than the Γgas = 13/9 value often adopted in single-fluid GRMHD simulations. We find moderate ion-to-electron temperature ratios in the 230 GHz emitting region of R = Ti/Te ≈ 5. While total intensity 230 GHz images from our simulations are consistent with Event Horizon Telescope (EHT) results, our images have significantly more beam-scale linear polarization (〈|m|〉 ≈ 30 %) than is observed in EHT images of M87* (〈|m|〉 < 10 %). We find a trend of the average linear polarization pitch angle ∠β2 with black hole spin consistent with what is seen in single-fluid GRMHD simulations, and we provide a simple fitting function for ∠β2(a*) motivated by the wind-up of magnetic field lines by black hole spin in the Blandford-Znajek mechanism.
Read full abstract