Abstract We investigate general relativistic magnetohydrodynamic simulations to determine the physical origin of the twisty patterns of linear polarization seen in spatially resolved black hole images and explain their morphological dependence on black hole spin. By characterizing the observed emission with a simple analytic ring model, we find that the twisty morphology is determined by the magnetic field structure in the emitting region. Moreover, the dependence of this twisty pattern on spin can be attributed to changes in the magnetic field geometry that occur due to the frame dragging. By studying an analytic ring model, we find that the roles of Doppler boosting and lensing are subdominant. Faraday rotation may cause a systematic shift in the linear polarization pattern, but we find that its impact is subdominant for models with strong magnetic fields and modest ion-to-electron temperature ratios. Models with weaker magnetic fields are much more strongly affected by Faraday rotation and have more complicated emission geometries than can be captured by a ring model. However, these models are currently disfavoured by the recent EHT observations of M87*. Our results suggest that linear polarization maps can provide a probe of the underlying magnetic field structure around a black hole, which may then be usable to indirectly infer black hole spins. The generality of these results should be tested with alternative codes, initial conditions, and plasma physics prescriptions.
more »
« less
The EB Correlation in Resolved Polarized Images: Connections to the Astrophysics of Black Holes
Abstract We present an in-depth analysis of the newly proposed correlation function in visibility space, between the E and B modes of linear polarization, hereafter the EB correlation, for a set of time-averaged general relativistic magnetohydrodynamical simulations compared with the phase map from different semianalytic models and the Event Horizon Telescope (EHT) 2017 data for M87*. We demonstrate that the phase map of time-averaged EB correlation contains novel information that might be linked to black hole (BH) spin, accretion state, and electron temperature. A detailed comparison with a semianalytic approach with different azimuthal expansion modes shows that to recover the morphology of real/imaginary part of the correlation function and its phase, we require higher orders of azimuthal modes. To extract the phase features, we use Zernike polynomial reconstruction developing an empirical metric to break degeneracies between models with different BH spins that are qualitatively similar. We use a set of geometrical ring models with various magnetic and velocity field morphologies, showing that both the image space and visibility-based EB -correlation morphologies in magnetically arrested disk simulations can be explained with simple fluid and magnetic field geometries as used in ring models. Standard and normal evolutions by contrast are harder to model, demonstrating that the simple fluid and magnetic field geometries of ring models are not sufficient to describe them owing to higher Faraday rotation depths. A qualitative comparison with the EHT data demonstrates that some of the features in the phase of EB correlation might be well explained by the current models for BH spins and electron temperatures, while others require larger theoretical surveys.
more »
« less
- PAR ID:
- 10465892
- Author(s) / Creator(s):
- ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more »
- Publisher / Repository:
- American Astronomical Society
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 955
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 6
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The spin of a newly formed black hole (BH) at the center of a massive star evolves from its natal value due to two competing processes: accretion of gas angular momentum that increases the spin and extraction of BH angular momentum by outflows that decreases the spin. Ultimately, the final, equilibrium spin is set by a balance between both processes. In order for the BH to launch relativistic jets and power aγ-ray burst (GRB), the BH magnetic field needs to be dynamically important. Thus, we consider the case of a magnetically arrested disk (MAD) driving the spin evolution of the BH. By applying the semianalytic MAD BH spin evolution model of Lowell et al. to collapsars, we show that if the BH accretes ∼20% of its initial mass, its dimensionless spin inevitably reaches small values,a≲ 0.2. For such spins, and for mass accretion rates inferred from collapsar simulations, we show that our semianalytic model reproduces the energetics of typical GRB jets,Ljet∼ 1050erg s−1. We show that our semianalytic model reproduces the nearly constant power of typical GRB jets. If the MAD onset is delayed, this allows powerful jets at the high end of the GRB luminosity distribution,Ljet∼ 1052erg s−1, but the final spin remains low,a≲ 0.3. These results are consistent with the low spins inferred from gravitational wave detections of binary BH mergers. In a companion paper by Gottlieb et al., we use GRB observations to constrain the natal BH spin to bea≃ 0.2.more » « less
-
ABSTRACT Event Horizon Telescope (EHT) observations of the core of the galaxy M87 suggest an observational appearance dominated by a ring of approximately 40 $$\mu$$as in diameter. The thickness of the ring is less certain: imaging efforts constrained it to be less than half the diameter (consistent with an imaging resolution of 20 $$\mu$$as), while visibility-domain modelling suggested a variety of fractional widths, including as low as $$10{{\ \rm per\ cent}}$$ on some days. The fractional width is very interesting as it has the potential to discriminate between different astrophysical scenarios for the source; in fact, the 10–$$20{{\ \rm per\ cent}}$$ range is so narrow as to be in tension with theoretical expectations. In the first of a series of papers on the width of the observed ring, we reproduce a subset of EHT visibility-domain modelling results and we explore whether alternative data analysis methods might favour thicker rings. We point out that the closure phase (and closure amplitude) likelihood function is not independent of residual station gain amplitudes, even at high signal-to-noise, and explore two approximations of practical interest: one standard in the field (and employed by the EHT collaboration), and a new one that we propose. Analysing the public data, we find that the new likelihood approximation prefers somewhat thicker rings, more in line with theoretical expectations. Further analysis is needed, however, to determine which approximation is better for the EHT data.more » « less
-
Context. Accreting black holes (BHs) may be surrounded by a highly magnetized plasma threaded by an organized poloidal magnetic field. Nonthermal flares and power-law spectral components at high energy could originate from a hot, collisionless, and nearly force-free corona. The jets we often observe from these systems are believed to be rotation-powered and magnetically driven. Aims. We study axisymmetric BH magnetospheres, where a fraction of the magnetic field lines anchored in a surrounding disk are connected to the event horizon of a rotating BH. For different BH spins, we identify the conditions and sites of magnetic reconnection within 30 gravitational radii. Methods. With the fully general relativistic particle-in-cell code GRZeltron , we solve the time-dependent dynamics of the electron–positron pair plasma and of the electromagnetic fields around the BH. The aligned disk is represented by a steady and perfectly conducting plasma in Keplerian rotation, threaded by a dipolar magnetic field. Results. For prograde disks around Kerr BHs, the topology of the magnetosphere is hybrid. Twisted open magnetic field lines crossing the horizon power a Blandford-Znajek jet, while open field lines with their footpoint beyond a critical distance on the disk could launch a magneto-centrifugal wind. In the innermost regions, coupling magnetic field lines ensure the transfer of significant amounts of angular momentum and energy between the BH and the disk. From the Y point at the intersection of these three regions, a current sheet forms where vivid particle acceleration via magnetic reconnection takes place. We compute the synchrotron images of the current sheet emission. Conclusions. Our estimates for jet power and BH–disk exchanges match those derived from purely force-free models. Particles are accelerated at the Y point, which acts as a heat source for the so-called corona. It provides a physically motivated ring-shaped source of hard X-rays above the disk for reflection models. Episodic plasmoid ejection might explain millisecond flares observed in Cygnus X-1 in the high-soft state, but are too fast to account for daily nonthermal flares from Sgr A * . Particles flowing from the Y point down to the disk could produce a hot spot at the footpoint of the outermost closed magnetic field line.more » « less
-
Abstract Filamentary structures in neutral hydrogen (Hi) emission are well aligned with the interstellar magnetic field, so Hiemission morphology can be used to construct templates that strongly correlate with measurements of polarized thermal dust emission. We explore how the quantification of filament morphology affects this correlation. We introduce a new implementation of the Rolling Hough Transform (RHT) using spherical harmonic convolutions, which enables efficient quantification of filamentary structure on the sphere. We use this Spherical RHT algorithm along with a Hessian-based method to construct Hi-based polarization templates. We discuss improvements to each algorithm relative to similar implementations in the literature and compare their outputs. By exploring the parameter space of filament morphologies with the Spherical RHT, we find that the most informative Histructures for modeling the magnetic field structure are the thinnest resolved filaments. For this reason, we find a ∼10% enhancement in theB-mode correlation with polarized dust emission with higher-resolution Hiobservations. We demonstrate that certain interstellar morphologies can produce parity-violating signatures, i.e., nonzeroTBandEB, even under the assumption that filaments are locally aligned with the magnetic field. Finally, we demonstrate thatBmodes from interstellar dust filaments are mostly affected by the topology of the filaments with respect to one another and their relative polarized intensities, whereasEmodes are mostly sensitive to the shapes of individual filaments.more » « less