The Spin of M87*

TL;DR

This paper constrains the spin and magnetic flux of M87*'s black hole using relativistic jet measurements and models, indicating a moderate spin and a magnetically arrested disk state.

Contribution

It provides the first constraints on M87*'s spin and magnetic flux using jet power and accretion rate models, favoring a magnetically arrested disk.

Findings

Lower limit on spin a* > 0.4 (prograde) and > 0.5 (retrograde)

Lower limit on magnetic flux phi > 6 (prograde) and > 10 (retrograde)

Disfavors low phi SANE models

Abstract

Now that the mass of the central black hole in the galaxy M87 has been measured with great precision using different methods, the remaining parameter of the Kerr metric that needs to be estimated is the spin a*. We have modeled measurements of the average power of the relativistic jet and an upper limit to the mass accretion rate onto the black hole with general relativistic magnetohydrodynamic models of jet formation. This allows us to derive constraints on a* and the black hole magnetic flux phi. We find a lower limit on M87*'s spin and magnetic flux of a* > 0.4 and phi > 6 in the prograde case, and a* > 0.5 and phi > 10 in the retrograde case, otherwise the black hole is not able to provide enough energy to power the observed jet. These results indicate that M87* has a moderate spin at minimum and disfavor a variety of models typified by low values of phi known as ``SANE'', indicating that M87* prefers the magnetically arrested disk state. We discuss how different estimates of the jet power and accretion rate can impact a* and phi.

Equations

$$\eta ={\left(\frac{\varphi }{15}\right)}^{2}f(a_{\ast })$$

$${\eta }_{\mathrm{observed}}=\eta (a_{\ast },\varphi ).$$

Full text

accretion, accretion disks – black hole physics – galaxies: active – galaxies: individual (M87) – galaxies: jets

††journal: ApJL††facilities: Chandra X-ray Observatory, Submillimeter Array, Event Horizon Telescope††software: Jupyter (Pérez & Granger, 2007), mcerp, nmmn

1 Introduction

The shadow cast by the event horizon of a black hole (BH) has been imaged for the first time with the Event Horizon Telescope (EHT) for the supermassive black hole (SMBH) at the center of the galaxy M87, known as M87* (Event Horizon Telescope Collaboration et al. 2019a, hereafter EHTC1). The very long baseline interferometry (VLBI) observation at a wavelength of 1.3 mm of the asymmetric bright emission ring gives an angular diameter of $d=42\pm 3\ \mu$as, which allows an unprecedented constraint on the SMBH mass of $M=(6.5\pm 0.7)\times 10^{9}M_{\odot}$–the first fundamental parameter of the Kerr (1963) metric. The inferred mass is in agreement with–and hence strongly favors–the mass measurement based on stellar dynamics (Gebhardt et al., 2011).

However, the second parameter of the Kerr metric–the dimensionless spin $a_{*}\equiv Jc/GM^{2}$ where $J$ is the angular momentum of the black hole–is much harder to constrain using only shadow observations. One of the reasons is that the ring diameter has a very weak dependence on $a_{*}$ and the disk inclination, varying by only 4% in the range $a_{*}=0$ to $\approx 1$ (Takahashi, 2004; Johannsen & Psaltis, 2010).

Besides mass and spin, black hole accretion is also described by two other parameters: the mass accretion rate onto the BH $\dot{M}$ and the magnetic flux $\Phi$ crossing one hemisphere of the event horizon. Based on the properties of the asymmetric ring observed in 2017 with the EHT as well as other constraints, Event Horizon Telescope Collaboration et al. (2019b) (hereafter EHTC5) obtained some bounds on $\Phi$ and $a_{*}$. Models with high, retrograde spins and high $\Phi$ are rejected, as well as $a_{*}=0$ models which fail to provide enough power to the jet.

In this Letter, we set further bounds on the values of $a_{*}$ and $\Phi$ of M87* by modeling the energetics of its relativistic jet as being powered by the Blandford-Znajek process–the extraction of BH spin energy through electromagnetic torques. The observational input to our estimates are the measurements of $M$, $\dot{M}$ and the power carried by the jet. We assume a distance to M87* of 16.8 Mpc (e.g. Blakeslee et al. 2009; Event Horizon Telescope Collaboration et al. 2019c). At this distance, cosmological effects are negligible.

2 Observations

In this work, the fundamental quantity needed in order to constrain the BH spin in M87* is the efficiency of jet production $\eta\equiv P/\dot{M}c^{2}$, where $P$ is the jet power and $\dot{M}$ is the mass accretion rate onto the SMBH. Therefore, $P$ and $\dot{M}$ are the M87* observables that we need in order to apply our models of jet production to constrain the spin parameter.

There are different ways of measuring the jet power of M87, with different methods giving powers in the range $P\sim 10^{42}-10^{45}\ {\rm erg\ s}^{-1}$ (e.g., Reynolds et al. 1996; Allen et al. 2006; Abdo et al. 2009; de Gasperin et al. 2012; Nemmen et al. 2014; EHTC5 and references therein). Here, we use the jet-inflated X-ray cavities observed in the central regions of M87 with the Chandra X-ray Telescope as calorimeters to estimate the jet power (Russell et al., 2013). The jet power was estimated as $P=E_{\rm cav}/t_{\rm age}$, where $E_{\rm cav}$ is the energy required to create the observed cavities and $t_{\rm age}$ is the age of the cavity. The usual assumption in deriving $E_{\rm cav}$ is that the cavities are inflated slowly such that $E_{\rm cav}=4PV$ where $P$ is the thermal pressure of the surrounding X-ray emitting gas, $V$ is the volume of the cavity and the cavity is assumed to be filled up with relativistic plasma. This method of measuring $P$ is well-established and robust (e.g. Bîrzan et al. 2004; Dunn & Fabian 2004; McNamara & Nulsen 2012; Hlavacek-Larrondo et al. 2013), giving the jet power averaged over the timescale during which the central engine produces one continuous pair of jets ($\sim 10^{6}$ years for M87; Allen et al. 2006). We believe this is the most direct way of measuring the jet power and therefore accept the X-ray cavity power at face value as M87’s jet power, $\log P=42.9^{+0.27}_{-0.2}\ {\rm erg\ s}^{-1}$.

For the mass accretion rate onto the BH, we use the constraint obtained by Kuo et al. (2014) based on the Faraday rotation measure (RM) observed with the Submillimeter Array. Kuo et al. measured the M87*’s RM at four frequencies around 230 GHz to be in the range $−7.5\times 10^{5}-3.4\times 10^{5}\ {\rm rad\ m}^{-2}$. By making reasonable assumptions, Kuo et al. estimated the upper limit $\dot{M}(42r_{g})\leq 9.2\times 10^{-4}\ M_{\odot}\ {\rm yr}^{-1}$ where $\dot{M}(42r_{g})$ is the mass accretion rate at a distance of 42 gravitational radii from the SMBH ($r_{g}\equiv GM/c^{2}$). The corresponding accretion power is $\log(\dot{M}(42r_{g})c^{2})=43.7\ {\rm erg\ s}^{-1}$. The question is of course how to connect this measurement with the accretion rate $\dot{M}$ near the event horizon. We will come back to this question in section 3.

3 Models

The leading idea to understand how accreting Kerr BHs produce relativistic jets is the Blandford-Znajek process (Blandford & Znajek, 1977; Blandford et al., 2019). According to this mechanism, the rotating event horizon is threaded by large-scale magnetic field lines, which were brought in by accreted gas. The BH exerts a torque on the field lines and progressively transfers its rotational energy to the relativistic jet (e.g. McKinney & Gammie 2004; Semenov et al. 2004). According to the Blandford-Znajek model, the jet power depends on $a_{*}$ and the magnetic flux $\Phi$ threading the horizon as $P\propto(a_{*}\Phi/M)^{2}$ to first order (Blandford & Znajek, 1977). In reality, for rapidly spinning BHs the jet power depends in a more complicated nonlinear fashion on the spin as higher-order spin corrections become important (e.g. Tchekhovskoy et al. 2010).

We use the results of global, general relativistic magnetohydrodynamic (GRMHD) simulations of radiatively inefficient accretion flows (RIAFs; Yuan & Narayan 2014) around Kerr BHs to model the dependence of the jet power on $a_{*}$ and $\Phi$. Concretely, we model the jet power as $P=\eta(a_{*},\Phi)\dot{M}c^{2}$ where the jet production efficiency

[TABLE]

is a function of both the spin and dimensionless magnetic flux $\phi\equiv\Phi/(\dot{M}r_{g}^{2})^{1/2}$. We use the GRMHD results of Tchekhovskoy et al. (2012) which are based on the HARM code (Gammie et al., 2003) to set the spin-dependence of $\eta$. Tchekhovskoy et al. (2012) carried out RIAF simulations in the “magnetically arrested disk” (MAD) limit, for which the magnetic flux saturates at $\phi\sim 15$ (Narayan et al., 2003; Tchekhovskoy et al., 2011), with $h/r\approx 0.3$ where $h$ is the disk thickness. In order to consider the full range of astrophysically relevant magnetic fluxes, we also take into account the case of “standard and normal disk evolution” (SANE) with $\phi\sim 1$ (e.g. Narayan et al. 2012). Figure 1 illustrates the spin-dependence of $\eta$. Notice that this model encompasses both the prograde and retrograde cases in which the disk and BH are rotating in the same and opposite senses, respectively. Retrograde BHs produce less powerful jets (Tchekhovskoy & McKinney, 2012).

In order to connect the accretion rate at $r=42r_{g}$ constrained by the Faraday RMs of Kuo et al. (2014) with the BH rate, we adopt the simple radial scaling $\dot{M}(r)=\dot{M_{0}}(r/r_{0})^{s}$ originally proposed as an ansatz by Blandford & Begelman (1999) in their “ADIOS” model–and later supported by many global simulations (e.g. Yuan et al. 2015). This $\dot{M}(r)$ scaling corresponds to a density radial profile $\rho(r)\propto r^{-\beta}$ where $\beta=3/2-s$. We define the BH accretion rate as $\dot{M}(6r_{g})$ so we fix $r=6r_{g}$, $r_{0}=42r_{g}$ and $\dot{M}_{0}$ as the accretion rate constrained by Kuo et al. We want to be agnostic regarding the variety of possible density profiles in M87*, therefore we allow $\beta$ ($s$) to vary in the range $1.5-0.5$ ($0-1$), i.e. allowing for different levels of mass-loss in the RIAF. We should note that Russell et al. (2018) measured $\beta\approx 1.5$ at $r=(0.1-1)$ kpc in M87, which is outside but very close to the Bondi radius ($r_{B}=0.03$ kpc).

4 Results

Our first result is a model-independent estimate of the jet production efficiency from the SMBH in M87* from the observed jet power and mass accretion rate. Figure 2 shows this result allowing a variety of density profiles, with $\eta$ varying from $\approx 10\%$ ($1\sigma$ lower limit, $\beta=1.5$) up to about $200\%$ ($1\sigma$ upper limit, $\beta=0.5$) if the density profile flattens towards the BH.

With the considerations in the previous section, we have a model that provides a full mapping of the observed jet efficiency derived above to the spin and magnetic flux

[TABLE]

We proceed by solving this nonlinear equation, using the values of $\eta_{\rm observed}$ displayed in Figure 2 to constrain the physical parameters of the event horizon beyond its mass–assuming that the Kerr metric is the correct description of the spacetime. Figure 3 shows the inferred spin of M87* on the assumption that the SMBH is in the MAD state–i.e. with the maximum value of $\phi$.

The lessons behind Fig. 3 are the following: (i) If M87* is in the MAD state then the only allowed spins are $a_{*}\lesssim-0.5$ or $a_{*}\gtrsim 0.4$ (within the $1\sigma$ uncertainty bands). These are effectively lower limits on $a_{*}$. (ii) If the density profile of the accretion flow follows $\rho\propto r^{-1.5}$ as in RIAF models without mass-loss (e.g. Narayan & Yi 1994), then $a_{*}=0.45^{+0.12}_{-0.08}$ (prograde) or $a_{*}=-0.62^{+0.14}_{-0.30}$ (retrograde). (iii) If the BH is retrograde then values of $\beta\geq 0.8$ are favored otherwise M87* would not be able to power its jet through the Blandford-Znajek process; all values of $\beta$ are allowed if the BH is prograde. Notice that we limit the upper value of $|a_{*}|$ to one because the BH solution is not valid anymore above this limit.

Figure 4 shows the solutions of the equation 2 for $a_{*}$ and $\phi$ that are consistent with the mean values of $\eta_{\rm observed}$ (i.e. the values along the solid line in Fig. 2). As such, Fig. 4 gives us the observational constraints on M87*’s spin and magnetic flux. To begin with, the hatched area in the plot indicates the region of the parameter space which is forbidden for M87* on the assumption of the Kerr metric, because it would imply $|a_{*}|>0.998$ which is the astrophysical limiting value of the spin (Thorne, 1974). In other words, in the hatched region the SMBH does not provide enough energy to power the observed jet. We now describe separately the prograde and retrograde cases. The MAD state corresponds to the top part of the plot ($\phi\approx 15$) while the SANE mode with $\phi$ close to one–as considered among the models in EHTC5–corresponds to the bottom region.

Prograde. In this case, not all magnetic fluxes are accessible to the SMBH. For instance, only accretion flows with $\phi\gtrsim 5$ are permitted. In the extreme situation that $\beta=0.5$, the RIAF must be in the MAD state. All values of $\beta$ are allowed.

Retrograde. Retrograde BHs produce less powerful jets, therefore if the SMBH is retrograde this implies tighter constraints on M87*’s parameters compared to the prograde case. Only RIAFs with $\phi\gtrsim 10$ and $\beta\gtrsim 1.1$ are possible.

5 Discussion

Our results are fundamentally based on the modeling of $\eta_{\rm observed}$–the observationally constrained jet production efficiency of M87*. We compute $\eta_{\rm observed}$ from the jet power measured from the X-ray cavities and BH accretion rate upper limit measured from radio polarization. From our modeling of $\eta_{\rm observed}$, we rule out a considerable region of the BH accretion parameter space for M87*. We have found that most SANE models with $\phi\lesssim 5$–i.e. the family of SANE models considered in EHTC1, EHTC5–are inconsistent with the jet energetics since their magnetic flux gives very small jet efficiencies. Most of the MAD models considered in EHTC5 are consistent with our $\phi$-constraints. Reassuringly, the MAD models with $a_{*}=-0.94$ and $\phi=8.04$ rejected by EHTC5 because they fail to produce stable images are also rejected in our work. Further constraints on $a_{*}$ and $\phi$ can be obtained by using the image scoring in EHTC5 as a prior to the modeling performed here. This deserves further investigation in the future.

The spin bounds from this work encompass–but are less restrictive–than previous estimates in the literature based on semi-analytic spectral fits (Feng & Wu, 2017), TeV pair production (Li et al., 2009) and jet wobbling (Sob’yanin, 2018).

It is interesting to discuss one notable difference between our results and those of EHTC5. The SANE models with $\phi\approx 1$ and $a_{*}=-0.94$ simulated by EHTC5 are characterized by low jet efficiencies $\eta=5\times 10^{-3}$, therefore they do not agree with $\eta_{\rm observed}>0.08$ of this work. However, they produce jet powers in the range $P\sim 10^{42}-10^{43}\ {\rm erg\ s}^{-1}$ and pass the “jet power consistency test” of EHTC5. What is the origin of this apparent contradiction between the simulated large $P$ and low $\eta$ for these models? The answer is that $\dot{M}$ in the EHTC5 models is adjusted to produce the observed compact mm flux and is not bound by the upper limit provided by the RM of Kuo et al. (2014)–which we respect in our models. The exact value of $\dot{M}$ in the general relativistic ray-tracing simulations employed in EHTC5 depends on the electron thermodynamics and spans a wide range.

M87*’s mass accretion rate upper limit is obtained using a model for the density and magnetic field strength and geometry in the RIAF in which it is roughly spherical with $\rho(r)\propto r^{-\beta}$ and the magnetic field is well ordered, radial and of equipartition strength (cf. Marrone et al. 2006; Kuo et al. 2014). There are reasons to believe that the mapping between RM and $\dot{M}$ may not be as straightforward as the prediction by the simple model above (e.g. Mościbrodzka et al. 2017). For instance, since M87’s jet has a low viewing angle (e.g. Walker et al. 2018) it is not impossible that the line of sight does not pass through the accretion flow and the RM measured by Kuo et al. comes from the jet sheath. If this is the case, the RM would be consistent with accretion rates potentially much larger than considered by Kuo et al. (Mościbrodzka et al., 2017). In the extreme case that $\dot{M}$ is 100 times larger than the Kuo et al. value and assuming our fiducial $P$, then we find $|a_{*}|>0.1$, $\phi>0.5$ in the prograde case and $|a_{*}|>0.15$, $\phi>1$ in the retrograde case.

What would happen if the current value of $P$ is different than the measurement we have adopted due to e.g. variability? If $P$ is on the lower end of the estimates at $10^{42}\ {\rm erg\ s}^{-1}$ (e.g. Prieto et al. 2016) then the lower limits on $a_{*}$ and $\phi$ are somewhat alleviated: $|a_{*}|>0.2$, $\phi>2$ in the prograde case and $|a_{*}|>0.3$, $\phi>4$ in the retrograde one. On the other hand, if the actual $P$ lies on the high end of the estimates at $10^{45}\ {\rm erg\ s}^{-1}$ (e.g. de Gasperin et al. 2012) then the BH is unable to provide enough power to the jet via the Blandford-Znajek process given the $\dot{M}$ measured by Kuo et al. (2014). The maximum jet power that can be extracted from M87* fed at this $\dot{M}$ is $6\times 10^{43}\ {\rm erg\ s}^{-1}$, similarly to the most powerful jets in the EHTC5 models. If $P$ is established to be as high as $10^{45}\ {\rm erg\ s}^{-1}$, then either the $\dot{M}$ measurement of Kuo et al. is innacurate or there is a different mechanism (not Blandford-Znajek) powering the jet.

One assumption of our work is that the dependency of $\eta$ on the BH parameters is adequately described by the time-average of global ideal GRMHD simulations of relativistic jet formation. These models have achieved an impressive level of sophistication in the last few years and the results of different GRMHD codes broadly agree with each other (e.g. Porth et al. 2019). However, simulations which incorporate kinetic effects are suggesting that jet formation can be more complex than previously thought, with a significant population of particles with negative energies measured by a far observer that could lead to strong extraction of the BH spin energy via the Penrose process (Parfrey et al., 2019). This could alter expectations for the radiation spectrum from BH accretion. Furthermore, we have only considered the case of a BH spin vector parallel to the angular momentum vector of the accretion flow. The general case of arbitrary relative orientations of these vectors can lead to more complicated jet behaviors (Liska et al., 2018). These issues deserve further investigations.

Finally, as we already mentioned the BH shadow is weakly affected by the different values of $a_{*}$ and it is not clear whether further EHT observations will eventually be able to constrain the ring properties with enough sharpness to meaningfully constrain the spin. We are looking forward to advancements in VLBI techniques that will allow better measurements of the ring properties and potentially advance our understanding of M87*’s spin (e.g. going to shorter wavelengths, adding more telescopes, going to space-based interferometry; e.g. Roelofs et al. 2019) and also to upcoming polarimetric analysis of the EHT observations which will further constrain the magnetic flux and accretion rates.

6 Summary

We have compared measurements of the power carried by particles in the relativistic jet and the black hole mass accretion rate with the predictions of GRMHD models of jet formation, which allows us to derive constraints on the spin and magnetic flux of M87*’s supermassive black hole. The jet power comes from X-ray observations and corresponds to the time-average over $\sim 10^{6}$ years, while the accretion rate is estimated assuming the Faraday rotation measures come from external Faraday rotation due to the RIAF. Our main conclusions can be summarized as follows:

(i) The black hole in M87* is converting at least $\eta=10\%$ of the accreted rest mass energy to jet power, and up to $200\%$ depending on how shallow the density profile of the accretion flow is.

(ii) We derive a lower limit on M87*’s spin: $|a_{*}|\geq 0.4$ in the prograde case or $|a_{*}|\geq 0.5$ in the retrograde case. We are not able to distinguish between the prograde or retrograde scenarios based only on the data we have used.

(iii) We obtained lower limits on the BH magnetic flux, potentially ruling out a variety of models with low values of $\phi$ known as “SANE”. We find that $\phi\gtrsim 5$ in the prograde case and $\phi\gtrsim 10$ in the retrograde case. Therefore, the magnetically arrested disk state seems to be preferred by M87*.

A possibility that cannot currently be excluded is that the RM could be jet-dominated, given the low inclination angle to the relativistic jet. In this case, the accretion rate could potentially be much larger than we considered. This would be consistent with low spins ($|a_{*}|\geq 0.1$) or SANE accretion ($\phi\gtrsim 1$).

We hope that these constraints on the M87*’s black hole spin and magnetic flux will be useful in further extracting physical parameters from M87*’s BH shadow.

We would like to thank Sasha Tchekhovskoy for providing us with the electronic data for his GRMHD simulations in a convenient format. We thank Fabio Cafardo, Gustavo Soares and Sasha Tchekhovskoy for useful discussions, and the anonymous referee for a constructive report that led to improvements in the paper. This work was supported by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) under grant 2017/01461-2. The Black Hole Group at USP acknowledges the generous donation of a GPU from the NVIDIA Corporation.