Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability

Pressure-anisotropy-induced nonlinearities in the kinetic magnetorotational instability

Squire, J.; Quataert, E.; Kunz, M. W.

In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are very small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low-collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the `parasitic instabilities’ that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in magnetohyrodynamics). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large-amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.

 

Stochastic Electron Acceleration by the Whistler Instability in a Growing Magnetic Field

Stochastic Electron Acceleration by the Whistler Instability in a Growing Magnetic Field

Riquelme, Mario; Osorio, Alvaro; Quataert, Eliot

We use 2D particle-in-cell simulations to study the effect of the saturated whistler instability on the viscous heating and nonthermal acceleration of electrons in a shearing, collisionless plasma with a growing magnetic field, {\boldsymbol{B}}. In this setup, an electron pressure anisotropy with {p}\perp ,e> {p}| | ,e naturally arises due to the adiabatic invariance of the electron magnetic moment ({p}| | ,e and {p}\perp ,e are the pressures parallel and perpendicular to {\boldsymbol{B}}). If the anisotropy is large enough, then the whistler instability arises, efficiently scattering the electrons and limiting {{∆ }}{p}e (\equiv {p}\perp ,e-{p}| | ,e). In this context, {{∆ }}{p}e taps into the plasma velocity shear, producing electron heating by the so-called anisotropic viscosity. In our simulations, we permanently drive the growth of | {\boldsymbol{B}}| by externally imposing a plasma shear, allowing us to self-consistently capture the long-term, saturated whistler instability evolution. We find that besides the viscous heating, the scattering by whistler modes can stochastically accelerate electrons to nonthermal energies. This acceleration is most prominent when initially {β }ẽ 1, gradually decreasing its efficiency for larger values of {β }e (\equiv 8π {p}e/| {\boldsymbol{B}}{| }2). If initially {β }ẽ 1, then the final electron energy distribution can be approximately described by a thermal component, plus a power-law tail with a spectral index of ̃3.7. In these cases, the nonthermal tail accounts for ̃ 5 % of the electrons and for ̃ 15 % of their kinetic energy. We discuss the implications of our results for electron heating and acceleration in low-collisionality astrophysical environments, such as low-luminosity accretion flows.

 

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

O. Porth et al.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole

First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole

Event Horizon Telescope collaboration et al.

We present measurements of the properties of the central radio source in M87 using Event Horizon Telescope data obtained during the 2017 campaign. We develop and fit geometric crescent models (asymmetric rings with interior brightness depressions) using two independent sampling algorithms that consider distinct representations of the visibility data. We show that the crescent family of models is statistically preferred over other comparably complex geometric models that we explore. We calibrate the geometric model parameters using general relativistic magnetohydrodynamic (GRMHD) models of the emission region and estimate physical properties of the source. We further fit images generated from GRMHD models directly to the data. We compare the derived emission region and black hole parameters from these analyses with those recovered from reconstructed images. There is a remarkable consistency among all methods and data sets. We find that >50% of the total flux at arcsecond scales comes from near the horizon, and that the emission is dramatically suppressed interior to this region by a factor >10, providing direct evidence of the predicted shadow of a black hole. Across all methods, we measure a crescent diameter of 42 ± 3 μas and constrain its fractional width to be <0.5. Associating the crescent feature with the emission surrounding the black hole shadow, we infer an angular gravitational radius of GM/Dc 2 = 3.8 ± 0.4 μas. Folding in a distance measurement of 16.8-0.7+0.8 Mpc gives a black hole mass of M=6.5+/- 0.2| stat+/- 0.7| sys× 109 M. This measurement from lensed emission near the event horizon is consistent with the presence of a central Kerr black hole, as predicted by the general theory of relativity.

First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring

First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring

Event Horizon Telescope collaboration et al.

The Event Horizon Telescope (EHT) has mapped the central compact radio source of the elliptical galaxy M87 at 1.3 mm with unprecedented angular resolution. Here we consider the physical implications of the asymmetric ring seen in the 2017 EHT data. To this end, we construct a large library of models based on general relativistic magnetohydrodynamic (GRMHD) simulations and synthetic images produced by general relativistic ray tracing. We compare the observed visibilities with this library and confirm that the asymmetric ring is consistent with earlier predictions of strong gravitational lensing of synchrotron emission from a hot plasma orbiting near the black hole event horizon. The ring radius and ring asymmetry depend on black hole mass and spin, respectively, and both are therefore expected to be stable when observed in future EHT campaigns. Overall, the observed image is consistent with expectations for the shadow of a spinning Kerr black hole as predicted by general relativity. If the black hole spin and M87’s large scale jet are aligned, then the black hole spin vector is pointed away from Earth. Models in our library of non-spinning black holes are inconsistent with the observations as they do not produce sufficiently powerful jets. At the same time, in those models that produce a sufficiently powerful jet, the latter is powered by extraction of black hole spin energy through mechanisms akin to the Blandford-Znajek process. We briefly consider alternatives to a black hole for the central compact object. Analysis of existing EHT polarization data and data taken simultaneously at other wavelengths will soon enable new tests of the GRMHD models, as will future EHT campaigns at 230 and 345 GHz.

First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole

First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole

Event Horizon Telescope collaboration et al.

We present the first Event Horizon Telescope (EHT) images of M87, using observations from April 2017 at 1.3 mm wavelength. These images show a prominent ring with a diameter of ~40 micro-as, consistent with the size and shape of the lensed photon orbit encircling the “shadow” of a supermassive black hole. The ring is persistent across four observing nights and shows enhanced brightness in the south. To assess the reliability of these results, we implemented a two-stage imaging procedure. In the first stage, four teams, each blind to the others’ work, produced images of M87 using both an established method (CLEAN) and a newer technique (regularized maximum likelihood). This stage allowed us to avoid shared human bias and to assess common features among independent reconstructions. In the second stage, we reconstructed synthetic data from a large survey of imaging parameters and then compared the results with the corresponding ground truth images. This stage allowed us to select parameters objectively to use when reconstructing images of M87. Across all tests in both stages, the ring diameter and asymmetry remained stable, insensitive to the choice of imaging technique. We describe the EHT imaging procedures, the primary image features in M87, and the dependence of these features on imaging assumptions.

First M87 Event Horizon Telescope Results. III. Data Processing and Calibration

First M87 Event Horizon Telescope Results. III. Data Processing and Calibration

Event Horizon Telescope collaboration et al.

We present the calibration and reduction of Event Horizon Telescope (EHT) 1.3mm radio wavelength observations of the supermassive black hole candidate at the center of the radio galaxy M87 and the quasar 3C 279, taken during the 2017 April 5-11 observing campaign. These global very long baseline interferometric observations include for the first time the highly sensitive Atacama Large Millimeter/submillimeter Array (ALMA); reaching an angular resolution of 25 micro-as, with characteristic sensitivity limits of ~1 mJy on baselines to ALMA and ~10 mJy on other baselines. The observations present challenges for existing data processing tools, arising from the rapid atmospheric phase fluctuations, wide recording bandwidth, and highly heterogeneous array. In response, we developed three independent pipelines for phase calibration and fringe detection, each tailored to the specific needs of the EHT. The final data products include calibrated total intensity amplitude and phase information. They are validated through a series of quality assurance tests that show consistency across pipelines and set limits on baseline systematic errors of 2% in amplitude and 1 degree in phase. The M87 data reveal the presence of two nulls in correlated flux density at ~3.4 and ~8.3 giga-lambda and temporal evolution in closure quantities, indicating intrinsic variability of compact structure on a timescale of days, or several light-crossing times for a few billion solar-mass black hole. These measurements provide the first opportunity to image horizon-scale structure in M87.

First M87 Event Horizon Telescope Results. II. Array and Instrumentation

First M87 Event Horizon Telescope Results. II. Array and Instrumentation

Event Horizon Telescope collaboration et al.

The Event Horizon Telescope (EHT) is a very long baseline interferometry (VLBI) array that comprises millimeter- and submillimeter-wavelength telescopes separated by distances comparable to the diameter of the Earth. At a nominal operating wavelength of ∼1.3 mm, EHT angular resolution (λ/D) is ∼25 μas, which is sufficient to resolve nearby supermassive black hole candidates on spatial and temporal scales that correspond to their event horizons. With this capability, the EHT scientific goals are to probe general relativistic effects in the strong-field regime and to study accretion and relativistic jet formation near the black hole boundary. In this Letter we describe the system design of the EHT, detail the technology and instrumentation that enable observations, and provide measures of its performance. Meeting the EHT science objectives has required several key developments that have facilitated the robust extension of the VLBI technique to EHT observing wavelengths and the production of instrumentation that can be deployed on a heterogeneous array of existing telescopes and facilities. To meet sensitivity requirements, high-bandwidth digital systems were developed that process data at rates of 64 gigabit s−1, exceeding those of currently operating cm-wavelength VLBI arrays by more than an order of magnitude. Associated improvements include the development of phasing systems at array facilities, new receiver installation at several sites, and the deployment of hydrogen maser frequency standards to ensure coherent data capture across the array. These efforts led to the coordination and execution of the first Global EHT observations in 2017 April, and to event-horizon-scale imaging of the supermassive black hole candidate in M87.

First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole

First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole

Event Horizon Telescope collaboration et al.

When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42 ± 3 μas, which is circular and encompasses a central depression in brightness with a flux ratio ≳10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M = (6.5 ± 0.7) × 109 M . Our radio-wave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

Numerical Evaluation of the Relativistic Magnetized Plasma Susceptibility Tensor and Faraday Rotation Coefficients

Numerical Evaluation of the Relativistic Magnetized Plasma Susceptibility Tensor and Faraday Rotation Coefficients

Pandya, Alex; Chandra, Mani; Joshi, Abhishek; Gammie, Charles F.

Polarized models of relativistically hot astrophysical plasmas require transport coefficients as input: synchrotron absorption and emission coefficients in each of the four Stokes parameters, as well as three Faraday rotation coefficients. Approximations are known for all coefficients for a small set of electron distribution functions, such as the Maxwell-Jüttner relativistic thermal distribution, and a general procedure has been obtained by Huang & Shcherbakov for an isotropic distribution function. Here we provide an alternative general procedure, with a full derivation, for calculating absorption and rotation coefficients for an arbitrary isotropic distribution function. Our method involves the computation of the full plasma susceptibility tensor, which in addition to absorption and rotation coefficients may be used to determine plasma modes and the dispersion relation. We implement the scheme in a publicly available library (https://github.com/afd-illinois/symphony) with a simple interface, thus allowing for easy incorporation into radiation transport codes. We also provide a comprehensive survey of the literature and comparison with earlier results.