1. Introduction
The Event Horizon Telescope (EHT) Collaboration has recently published polarization images of the event-horizon-scale emission around the supermassive black hole at the center of the M87 galaxy (M87
; [
1,
2]). The images reveal that a significant fraction of the emission is linearly polarized, as expected for synchrotron emission from relativistic electrons gyrating around well organized magnetic field lines in the M87 jet and/or the surrounding accretion disk. The sign of the Faraday Rotation Measure (RM) derived from these images determines the direction of the magnetic field that threads the disk: by convention, negative values correspond to an average magnetic field pointing
away from the observer, while positive values
towards the observer. The M87 jet points towards us at
to our line-of-sight, and the disk around M87
rotates in the clockwise direction in the plane of the sky [
3]. This offers a unique opportunity to test the Cosmic Battery model for the origin of magnetic fields around astrophysical black holes which surmises that the large scale axial magnetic field in the immediate vicinity of the M87
black hole points away from us (along the direction of
of the disk), while further away it points towards us
1.
In the next Section, we will present a rough prediction for the distribution of the axial magnetic field around the central black hole according to the Cosmic Battery. We will then provide the general features of an RM map of the average axial field whenever it will become reliably available. In
Section 3 we discuss that, in order to confirm the prediction of the Cosmic Battery, we are only interested to know whether the average line-of-sight magnetic field changes direction at some particular distance from the center. Episodic enhancement of the radiation in such regions of opposite line-of-sight magnetic field direction could give rise to the observed rapid flips in the sign of the average RM [
5]. However, these RM sign flips may also have to do with the turbulent events observed in the General Relativistic Magnetohydrodynamic (GRMHD) numerical simulations of [
6]. Further monitoring of the polarization emission of M87
as well as polarization angle observations over a wider frequency range will help disentangle between these possible scenarios. We conclude in
Section 4 with a discussion of the prospects for testing the Cosmic Battery.
2. RM Map According to the Cosmic Battery
The Cosmic Battery model for the generation of poloidal magnetic fields around astrophysical black holes was first proposed by Contopoulos Kazanas [
7]. According to that mechanism, the aberrated radiation pressure on the electrons of the electron-proton plasma in orbit around the black hole decelerates the electrons and thus induces an azimuthal electric field opposite to the direction of rotation. The rotation (curl) of this electric field generates poloidal magnetic field flux which, near the black hole, is along the direction of disk-hole rotation
(‘magnetism along spin’, [
8]), and in the surrounding disk is opposite to
(e.g., [
9,
10]). In a disk that is observed almost face-on with
pointing away from us (as is the case with the clockwise rotating disk of M87
), one expects that the central parts of the image will show generally negative RM values, whereas the outer parts of the image will consistently show generally positive RM values.
In order to generate theoretical RM maps, we will use a simple model for the distribution of electron density
and axial magnetic field
with distance
R in the equatorial plasma disk of M87
. We will thus assume that
drops exponentially inside the distance of field reversal
, and inversely with distance outside. Similarly, we will assume that
stays roughly constant inside
, and also drops inversely with distance outside, according to the observational conclusions of Fukumura et al. [
11] on self-similar magnetic winds from astrophysical accretion disks. The distance
is currently a big unknown that will be determined with more numerical simulations analogous to those of [
12]. Finally, we will assume that most of the contribution to the RM comes from a thin disk with scale-height/thickness that increases proportionally with distance
R from the center. We will thus consider the ad hoc expressions
If we were to observe such a magnetic field/plasma distribution face-on from below, this would yield
Here, the normalization value
is taken to be equal to
. For simplicity, these approximate expressions ignore the internal Faraday depolarization although several models considered in EHT VIII [
2] (see [
6]) indicate that Faraday depolarization in the midplane is strong. We also ignore any contribution to the line-of-sight magnetic field and to the RM from the turbulent magnetic field in the disk. However, we note that the GRMHD numerical simulations of [
6] showed that there are periods during the evolution of the system that the turbulent magnetic field may dominate the RM, thus our expressions may be valid only on average, and during a period of quiescent (non-turbulent) disk activity.
The M87
disk is not observed face-on, thus photons emitted from the disk that reach the EHT observer cross the disk at some nonzero angle
from the vertical to the disk (see
Figure 1). For simplicity, we will only consider the line-of-sight contribution of the axial field
in the RM map, and ignore the line-of-sight contribution of the field in the equatorial plane of the disk. As we just acknowledged, this is not always true. Under this simplifying approximation, the local contribution to the RM from a path length
through the thin disk is equal to
, thus
In order to obtain RM maps on the sky, we will follow the approximate analytical formulation provided by Beloborodov [
13] and most recently by Narayan et al. [
14] to account for the geodesics of photons emitted from the M87
disk that reach the EHT observer on earth. According to Narayan et al. [
14], each point
in the disk is mapped onto
on the plane of the sky, where
are expressed in units of
(
M is the mass of the black hole,
G is Newton’s constant, and
c is the speed of light).
x is measured along the line of nodes of the plane of the disk along
, and
y is measured perpendicularly to it along
.
is the inclination of the M87
disk with respect to our line of sight. Note that the line of nodes is tilted with respect to the East-West horizontal on the sky by
clockwise. We have tacitly assumed here that the direction of the large scale jet in the galaxy M87 coincides with the z-direction perpendicular to the accretion disk around the central black hole in M87
. Obviously, the two may differ.
As we acknowledged above, the distance
is currently a big unknown. If we choose the distance of field reversal to concide with the position of the light emitting ring in the EHT image, namely
[
14], we obtain the theoretical RM map shown in
Figure 2.
3. Rapid RM Reversals
In order to generate RM maps, we need simultaneous knowledge of the Electric Vector Polarization Angles (EVPAs) in the high and low bands of the EHT polarization maps. The RM value at each point of the polarization map is determined as
Here, EVPAs are measured in degrees counter-clockwise from North (up), and “low”/“high” refer to the so-called low/high band observations of EHT at 227.1/229.1 GHz respectively. Although the raw EVPA data is not yet available from the EHT collaboration, a direct visual inspection of the fiducial M87
daily average maps shown in Figure 28 of EHT VII [
1] reveals EVPA differences between the low and high bands of up to about
, corresponding to RM values on the order of
, i.e., much higher than the expected contribution of foreground Faraday screens and the contribution of our own Galaxy which are all expected to be on the order of only
. Therefore, at the resolution of EHT the sign (positive/negative) of the RM values calculated directly through Equation (
8) represent the direction of the line-of-sight magnetic field (towards us/away from us respectively). The above visual inspection revealed various regions around M87
with either positive or negative differences between the EVPA of the two bands (
), corresponding to positive and negative RM values, respectively, as described in Equation (
8). However, the uncertainty of the EVPA measurements is probably of the order of
(e.g., see Figure 8 of EHT VII [
1]), which is similar to the EVPA differences betwen the low and high bands, and hence the RM sign cannot be reliably constrained by the visual inspection of the EHT polarization images of M87
described above. Nevertheless, we plan to perform detailed RM investigations as soon as the raw EVPA data become available.
One way to improve this unfortunate situation is to look at the average RM values. The published average RM values for the four days 5, 6, 10, and 11 April 2017 of EHT observations, as constrained with ALMA observations at 1.3 mm and 3 mm are
,
,
,
rad/m
2 respectively Goddi et al. [
5]. Since we are primarily interested in the sign of the line-of-sight magnetic fields, i.e., the sign in the RM maps, the above average RM values suggest that the RM sign of the polarization image alternates within a few days. This is a suggestion that the accretion disk in the immediate vicinity of the central black hole is threaded by both positive and negative vertical magnetic field
over extended parts of the inner disk. Similar transient structures have also been observed in the canonical MAD GRMHD simulations performed by the EHT Collaboration to account for the overall features of the EHT image of M87
[
6]. The initial conditions in these simulations assume a large scale axial magnetic field in the accretion disk that is unidirectional, i.e.,
without field reversals. Other researchers favor initial configurations
with field reversals (e.g., [
15]). Nevertheless, the fact that the turbulence in the disk can induce transient axial field reversals even in numerical simulations without initial field reversals, leads us to accept that the observation of RM sign flips in the EHT polarization images does not yet consist confirmation of the Cosmic Battery.
4. Conclusions
Confirmation of the Cosmic Battery would consist of observing a large-scale pattern of alternating field polarities similar to the one shown in
Figure 2 over several EHT observing periods with no transient RM sign flips. One possibility is that the interface between the two
polarities does not have the well-defined compact shape shown in
Figure 2, with “tongues” of alternating sign RM contributing to the observed RM values. In other words, it is not only the instantaneous RM map morphology that matters, but also its day-to-day variability. The determination of a large scale geometry such as that of
Figure 2 may require either a longer term monitoring, in order to determine an average value of the RM as a function of distance from the black hole, or even the synergy with lower resolution RM observations that can capture the spatial average of such a topology at both small and large scales.
The published EHT observations in the low and high bands cover a very narrow region of the RM space and hence they are insensitive to lower RM values which can be manifested by observations over a wider range of frequencies. This is more relevant for regions of low density plasma and/or low magnetic field strengths which give rise to the lower RM values. Such regions may belong to larger spatial scales where the field pattern predicted by the Cosmic Battery may be easier to detect. If indeed we observe an average pattern similar to the one shown in
Figure 2, this would be a hint that the standard scenario according to which the magnetic field that threads the disk is brought in from large distances and saturates in the immediate vicinity of the central black hole with one and the same polarity
may not work. The alternative is that dipolar magnetic flux of both polarities is generated locally by some physical mechanism like the one we propose in this paper. While it appears currently difficult to establish the average large-scale geometry of the near black hole axial field, we believe these observations present an interesting first step in this direction, with the hope future refinements will allow us to reach this goal.
It is interesting that the reversal of the magnetic field polarity at some distance from the central black hole will correspond to a
large scale current sheet. This configuration is expected to lead to several important features such as dynamic intermittent behavior with continuous reconnection and plasmoid formation along it. This may be a fundamental element of astrophysical jets according to several theoretical models of magnetically driven jets and winds (e.g., [
16]). The current sheet may naturally account for several characteristic features of AGN jets such as particle acceleration and enhanced radiation emission along its direction (as is the case with the pulsar curent sheet [
17]),
polarization angle swings as it sweeps the observer’s line-of-sight (e.g., [
18]), etc. Another important feature is that since there is a limited amount of magnetic flux threading the innermost accretion disk around the central black hole, the magnetic flux of the opposite (return) polarity threading the large scale disk is also limited. Therefore, for the return flux to fill the whole accretion disk, the magnetic field immediately outside the distance of field reversal will drop. We thus predict a region of weaker magnetic field (and thus also weaker synchrotron emission) surrounding the jet [
19]. This may be one of the reasons AGN jets have this characteristic sheath structure.
Time and more observations will tell if we will ever be able to observe the Cosmic Battery magnetic field configuration and the footpoint of the purported current sheet along the surface of the core jet in M87*.