Event Horizon Scales Research Articles

Imaging the event horizon of M87* from space on different timescales

Imaging the event horizon of M87* from space on different timescales

Autoencoding Labeled Interpolator, Inferring Parameters from Image and Image from Parameters

The Event Horizon Telescope (EHT) provides an avenue to study black hole accretion flows on event-horizon scales. Fitting a semianalytical model to EHT observations requires the construction of synthetic images, which is computationally expensive. This study presents an image generation tool in the form of a generative machine-learning model, which extends the capabilities of a variational autoencoder. This tool can rapidly and continuously interpolate between a training set of images and can retrieve the defining parameters of those images. Trained on a set of synthetic black hole images, our tool showcases success in interpolating both black hole images and their associated physical parameters. By reducing the computational cost of generating an image, this tool facilitates parameter estimation and model validation for observations of black hole systems.

First Sagittarius A* Event Horizon Telescope Results. VIII. Physical Interpretation of the Polarized Ring

In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow.

First M87 Event Horizon Telescope Results. IX. Detection of Near-horizon Circular Polarization

Event Horizon Telescope (EHT) observations have revealed a bright ring of emission around the supermassive black hole at the center of the M87 galaxy. EHT images in linear polarization have further identified a coherent spiral pattern around the black hole, produced from ordered magnetic fields threading the emitting plasma. Here we present the first analysis of circular polarization using EHT data, acquired in 2017, which can potentially provide additional insights into the magnetic fields and plasma composition near the black hole. Interferometric closure quantities provide convincing evidence for the presence of circularly polarized emission on event-horizon scales. We produce images of the circular polarization using both traditional and newly developed methods. All methods find a moderate level of resolved circular polarization across the image (〈∣v∣〉 < 3.7%), consistent with the low image-integrated circular polarization fraction measured by the Atacama Large Millimeter/submillimeter Array (∣v int∣ < 1%). Despite this broad agreement, the methods show substantial variation in the morphology of the circularly polarized emission, indicating that our conclusions are strongly dependent on the imaging assumptions because of the limited baseline coverage, uncertain telescope gain calibration, and weakly polarized signal. We include this upper limit in an updated comparison to general relativistic magnetohydrodynamic simulation models. This analysis reinforces the previously reported preference for magnetically arrested accretion flow models. We find that most simulations naturally produce a low level of circular polarization consistent with our upper limit and that Faraday conversion is likely the dominant production mechanism for circular polarization at 230 GHz in M87*.

Polarimetry and astrometry of NIR flares as event horizon scale, dynamical probes for the mass of Sgr A*

We present new astrometric and polarimetric observations of flares from Sgr A* obtained with GRAVITY, the near-infrared interferometer at ESO’s Very Large Telescope Interferometer (VLTI), bringing the total sample of well-covered astrometric flares to four and polarimetric flares to six. Of all flares, two are well covered in both domains. All astrometric flares show clockwise motion in the plane of the sky with a period of around an hour, and the polarization vector rotates by one full loop in the same time. Given the apparent similarities of the flares, we present a common fit, taking into account the absence of strong Doppler boosting peaks in the light curves and the EHT-measured geometry. Our results are consistent with and significantly strengthen our model from 2018. First, we find that the combination of polarization period and measured flare radius of around nine gravitational radii (9Rg ≈ 1.5RISCO, innermost stable circular orbit) is consistent with Keplerian orbital motion of hot spots in the innermost accretion zone. The mass inside the flares’ radius is consistent with the 4.297 × 106 M⊙ measured from stellar orbits at several thousand Rg. This finding and the diameter of the millimeter shadow of Sgr A* thus support a single black hole model. Second, the magnetic field configuration is predominantly poloidal (vertical), and the flares’ orbital plane has a moderate inclination with respect to the plane of the sky, as shown by the non-detection of Doppler-boosting and the fact that we observe one polarization loop per astrometric loop. Finally, both the position angle on the sky and the required magnetic field strength suggest that the accretion flow is fueled and controlled by the winds of the massive young stars of the clockwise stellar disk 1–5″ from Sgr A*, in agreement with recent simulations.

Testing the Linear Relationship between Black Hole Mass and Variability Timescale in Low-luminosity AGNs at Submillimeter Wavelengths

The variability of submillimeter emission provides a useful tool to probe the accretion physics in low-luminosity active galactic nuclei. We accumulate four years of observations using the Submillimeter Array for Centaurus A, NGC 4374, NGC 4278, and NGC 5077, and one year of observations for NGC 4552 and NGC 4579. All sources are variable. We measure the characteristic timescale at which the variability is saturated by modeling these sources’ light curves as a damped random walk. We detect a timescale for all the sources except NGC 4552. The detected timescales are comparable to the orbital timescale at the event horizon scale for most sources. Combined with previous studies, we show a correlation between the timescale and the black hole mass over 3 orders of magnitude. This discovery suggests the submillimeter emission is optically thin with the emission originating from the event horizon. The mass scaling relationship further suggests that a group of radio sources with a broadband spectrum that peaks at submillimeter wavelengths have similar inner accretion physics. Sources that follow this relationship may be good targets for high-resolution imaging with the Event Horizon Telescope.

Dynamical analysis of the Tsallis holographic dark energy models with event horizon as cut-off and interaction with matter

The model of generalized Tsallis holographic dark energy (which is known to be particular representative of Nojiri–Odintsov HDE) with event horizon as cut-off is investigated using methods of dynamical analysis. We take into consideration possible interaction with dark energy and matter in various forms. Critical points are determined. Cosmological evolution of the universe depends from interaction parameters. If we use event horizon scale as cut-off quasi-de Sitter expansion is possible only for interaction of type [Formula: see text] (where [Formula: see text] is the Hubble parameter). For interactions [Formula: see text] and [Formula: see text] universe eventually stops ([Formula: see text]) or ends its existence in final singularity ([Formula: see text]). In the first case, fraction of dark energy tends to [Formula: see text] or constant value lesser than 1 because dynamical equilibrium between matter and dark energy is established on late times.

Dynamic and Polarimetric VLBI imaging with a multiscalar approach

Context. Due to the limited number of antennas and the limited observation time, an array of antennas in very long baseline interfer-ometry (VLBI) often samples the Fourier domain only very sparsely. Powerful deconvolution algorithms are needed to compute a final image. Multiscale imaging approaches such as DoG-HiT have recently been developed to solve the VLBI imaging problem and show promising performance: they are fast, accurate, unbiased, and automatic. Aims. We extend the multiscalar imaging approach to polarimetric imaging, to reconstructions of dynamically evolving sources, and finally to dynamic polarimetric reconstructions. Methods. These extensions (mr-support imaging) utilize a multiscalar approach. The time-averaged Stokes I image was decomposed by a wavelet transform into single subbands. We used the set of statistically significant wavelet coefficients, the multiresolution support (mr-support), computed by DoG-HiT as a prior in a constrained minimization manner; we fitted the single-frame (polarimetric) observables by only varying the coefficients in the multiresolution support. Results. The Event Horizon Telescope (EHT) is a VLBI array imaging supermassive black holes. We demonstrate on synthetic data that mr-support imaging offers ample regularization and is able to recover simple geometric dynamics at the horizon scale in a typical EHT setup. The approach is relatively lightweight, fast, and largely automatic and data driven. The ngEHT is a planned extension of the EHT designed to recover movies at the event horizon scales of a supermassive black hole. We benchmark the performance of mr-support imaging for the denser ngEHT configuration demonstrating the major improvements the additional ngEHT antennas will bring to dynamic polarimetric reconstructions. Conclusions. Current and upcoming instruments offer the observational possibility to do polarimetric imaging of dynamically evolving structural patterns with the highest spatial and temporal resolution. State-of-the-art dynamic reconstruction methods can capture this motion with a range of temporal regularizers and priors. With this work, we add an additional simpler regularizer to the list: constraining the reconstruction to the multiresolution support.

Collisionless Accretion onto Black Holes: Dynamics and Flares.

We study the accretion of collisionless plasma onto a rotating black hole from first principles using axisymmetric general-relativistic particle-in-cell simulations. We carry out a side-by-side comparison of these results to analogous general-relativistic magnetohydrodynamic simulations. Although there are many similarities in the overall flow dynamics, three key differences between the kinetic and fluid simulations are identified. Magnetic reconnection is more efficient, and rapidly accelerates a nonthermal particle population, in our kinetic approach. In addition, the plasma in the kinetic simulations develops significant departures from thermal equilibrium, including pressure anisotropy that excites kinetic-scale instabilities, and a large field-aligned heat flux near the horizon that approaches the free-streaming value. We discuss the implications of our results for modeling event-horizon scale observations of Sgr A* and M87 by GRAVITY and the Event Horizon Telescope.

Toward Horizon-scale Accretion onto Supermassive Black Holes in Elliptical Galaxies

We present high-resolution, three-dimensional hydrodynamic simulations of the fueling of supermassive black holes in elliptical galaxies from a turbulent medium on galactic scales, taking M87* as a typical case. The simulations use a new GPU-accelerated version of the Athena++ AMR code, and they span more than six orders of magnitude in radius, reaching scales similar to that of the black hole horizon. The key physical ingredients are radiative cooling and a phenomenological heating model. We find that the accretion flow takes the form of multiphase gas at radii less than about a kpc. The cold gas accretion includes two dynamically distinct stages: the typical disk stage in which the cold gas resides in a rotationally supported disk, and relatively rare chaotic stages (≲10% of the time) in which the cold gas inflows via chaotic streams. Though cold gas accretion dominates the time-averaged accretion rate at intermediate radii, accretion at the smallest radii is dominated by hot virialized gas at most times. The accretion rate scales with radius as when hot gas dominates, and we obtain near the event horizon, similar to what is inferred from EHT observations. The orientation of the cold gas disk can differ significantly on different spatial scales. We propose a subgrid model for accretion in lower-resolution simulations in which the hot gas accretion rate is suppressed relative to the Bondi rate by . Our results can also provide more realistic initial conditions for simulations of black hole accretion at the event horizon scale.

Probing plasma physics with spectral index maps of accreting black holes on event horizon scales

ABSTRACT The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of the supermassive black holes at the centre of our galaxy and at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. In this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code ipole to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically show that the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disc (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

The ngEHT’s Role in Measuring Supermassive Black Hole Spins

While supermassive black-hole masses have been cataloged across cosmic time, only a few dozen of them have robust spin measurements. By extending and improving the existing Event Horizon Telescope (EHT) array, the next-generation Event Horizon Telescope (ngEHT) will enable multifrequency, polarimetric movies on event-horizon scales, which will place new constraints on the space-time and accretion flow. By combining this information, it is anticipated that the ngEHT may be able to measure tens of supermassive black-hole masses and spins. In this white paper, we discuss existing spin measurements and many proposed techniques with which the ngEHT could potentially measure spins of target supermassive black holes. Spins measured by the ngEHT would represent a completely new sample of sources that, unlike pre-existing samples, would not be biased towards objects with high accretion rates. Such a sample would provide new insights into the accretion, feedback, and cosmic assembly of supermassive black holes.

Observational Signatures of Frame Dragging in Strong Gravity

Objects orbiting in the presence of a rotating massive body experience a gravitomagnetic frame-dragging effect, known as the Lense-Thirring effect, that has been experimentally confirmed in the weak-field limit. In the strong-field limit, near the horizon of a rotating black hole, frame dragging becomes so extreme that all objects must co-rotate with the black hole’s angular momentum. In this work, we perform general relativistic numerical simulations to identify observable signatures of frame dragging in the strong-field limit that appear when infalling gas is forced to flip its direction of rotation as it is being accreted. In total intensity images, infalling streams exhibit “S”-shaped features due to the switch in the tangential velocity. In linear polarization, a flip in the handedness of spatially resolved polarization ticks as a function of radius encodes a transition in the magnetic field geometry that occurs due to magnetic flux freezing in the dragged plasma. Using a network of telescopes around the world, the Event Horizon Telescope collaboration has demonstrated that it is now possible to directly image black holes on event horizon scales. We show that the phenomena described in this work would be accessible to the next-generation Event Horizon Telescope and extensions of the array into space, which would produce spatially resolved images on event horizon scales with higher spatial resolution and dynamic range.

Synthetic Images of Magnetospheric Reconnection-Powered Radiation around Supermassive Black Holes.

Accreting supermassive black holes can now be observed at the event-horizon scale at millimeter wavelengths. Current predictions for the image rely on hypotheses (fluid modeling, thermal electrons) which might not always hold in the vicinity of the black hole, so that a full kinetic treatment is in order. In this Letter, we describe the first 3D global general-relativistic particle-in-cell simulation of a black-hole magnetosphere. The system displays a persistent equatorial current sheet. Synthetic radio images are computed by ray-tracing synchrotron emission from nonthermal particles accelerated in this current sheet by magnetic reconnection. We identify several time-dependent features of the image at moderate viewing angles: a variable radius of the ring, and hot spots moving along it. In this regime, our model predicts that most of the flux of the image lies inside the critical curve. These results could help promote understanding of future observations of black-hole magnetospheres at improved temporal and spatial resolution.

Implications from the Velocity Profile of the M87 Jet: A Possibility of a Slowly Rotating Black Hole Magnetosphere

Motivated by the measured velocity profile of the M87 jet using the Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) by Park et al. indicating that the starting position of the jet acceleration is farther from the central engine of the jet than predicted in general relativistic magnetohydrodynamic simulations, we explore how to mitigate the apparent discrepancy between the simulations and the KaVA observation. We use a semi-analytic jet model proposed by Tomimatsu & Takahashi consistently solving the transmagnetic field structure but neglecting any dissipation effects. By comparing the jet model with the observed M87 jet velocity profile, we find that the model can reproduce the logarithmic feature of the velocity profile and can fit the observed data when choosing c/(100r g ) ≲ Ω F ≲ c/(70r g ), where r g is the gravitational radius. While a total specific energy () of the jet changes the terminal bulk Lorentz factor of the jet, a slower angular velocity of the black hole magnetosphere (Ω F ) makes a light-cylinder radius (r lc) larger, and it consequently pushes out a location of a starting point of the jet acceleration. Using the estimated Ω F we further estimate the magnetic field strength on the event horizon scale in M87 by assuming Blandford–Znajek process is in action. The corresponding magnetic flux threading the event horizon of M87 is in good agreement with a magnetically arrested disk regime.

How narrow is the M87* ring – II. A new geometric model

ABSTRACTThe 2017 Event Horizon Telescope (EHT) observations of M87* detected a ring-shaped feature ∼40 μas in diameter, consistent with the event horizon scale of a black hole of the expected mass. The thickness of this ring, however, proved difficult to measure, despite being an important parameter for constraining the observational appearance. In the first paper of this series, we asked whether the width of the ring was sensitive to the choice of likelihood function used to compare observed closure phases and closure amplitudes to model predictions. In this paper, we investigate whether the ring width is robust to changes in the model itself. We construct a more realistic geometric model with two new features: an adjustable radial falloff in brightness, and a secondary ‘photon ring’ component in addition to the primary annulus. This thin, secondary ring is predicted by gravitational lensing for any black hole with an optically thin accretion flow. Analysing the data using the new model, we find that the primary annulus remains narrow (fractional width ≤ 0.25) even with the added model freedom. This provides further evidence in favour of a narrow ring for the true sky appearance of M87*, a surprising feature that, if confirmed, would demand theoretical explanation. Comparing the Bayesian evidence for models with and without a secondary ring, we find no evidence for the presence of a lensed photon ring in the 2017 observations. However, the techniques we introduce may prove useful for future observations with a larger and more sensitive array.

Millimeter Light Curves of Sagittarius A* Observed during the 2017 Event Horizon Telescope Campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5–11 in the 1.3 mm wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very long baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare and exhibit strongly enhanced variability. All of the light curves are consistent with a red-noise process, with a power spectral density (PSD) slope measured to be between −2 and −3 on timescales between 1 minute and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 hr. The spectral energy distribution is flat at 220 GHz, and there are no time lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*’s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

Black hole magnetic fields and their imprint on circular polarization images

ABSTRACT The circular polarization of black hole accretion flows can encode properties of the underlying magnetic field structure. Using general relativistic magnetohydrodynamic (GRMHD) simulations, we study the imprint of magnetic field geometry on circular polarization images potentially observable by the Event Horizon Telescope (EHT). We decompose images into the different mechanisms that generate circular polarization in these models that are sensitive to both the line-of-sight direction and twist of the magnetic field. In these models, a stable sign of the circular polarization over time, as observed for several sources, can be attributed to a stability of these properties. We illustrate how different aspects of a generic helical magnetic field geometry become imprinted on a circular polarization image. We also identify novel effects of light bending that affect the circular polarization image on event horizon scales. One consequence is the sign flipping of successive photon rings in face-on systems, which if observable and uncorrupted by Faraday rotation, can directly encode the handedness of the approaching magnetic field.

An Event Horizon Imager (EHI) Mission Concept Utilizing Medium Earth Orbit Sub-mm Interferometry&amp;lt;sup&amp;gt;*&amp;lt;/sup&amp;gt;

Submillimeter interferometry has the potential to image supermassive black holes on event horizon scales, providing tests of the theory of general relativity and increasing our understanding of black hole accretion processes. The Event Horizon Telescope (EHT) performs these observations from the ground, and its main imaging targets are Sagittarius A* in the Galactic Center and the black hole at the center of the M87 galaxy. However, the EHT is fundamentally limited in its performance by atmospheric effects and sparse terrestrial (<italic>u,v</italic>)-coverage (Fourier sampling of the image). The scientific interest in quantitative studies of the horizon size and shape of these black holes has motivated studies into using space interferometry which is free of these limitations. Angular resolution considerations and interstellar scattering effects push the desired observing frequency to bands above 500 GHz. This paper presents the requirements for meeting these science goals, describes the concept of interferometry from Polar or Equatorial Medium Earth Orbits (PECMEO) which we dub the Event Horizon Imager (EHI), and utilizes suitable space technology heritage. In this concept, two or three satellites orbit at slightly different orbital radii, resulting in a dense and uniform spiral-shaped (<italic>u,v</italic>)-coverage over time. The local oscillator signals are shared via an inter-satellite link, and the data streams are correlated on-board before final processing on the ground. Inter-satellite metrology and satellite positioning are extensively employed to facilitate the knowledge of the instrument position vector, and its time derivative. The European space heritage usable for both the front ends and the antenna technology of such an instrument is investigated. Current and future sensors for the required inter-satellite metrology are listed. Intended performance estimates and simulation results are given.

First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon

Abstract Event Horizon Telescope (EHT) observations at 230 GHz have now imaged polarized emission around the supermassive black hole in M87 on event-horizon scales. This polarized synchrotron radiation probes the structure of magnetic fields and the plasma properties near the black hole. Here we compare the resolved polarization structure observed by the EHT, along with simultaneous unresolved observations with the Atacama Large Millimeter/submillimeter Array, to expectations from theoretical models. The low fractional linear polarization in the resolved image suggests that the polarization is scrambled on scales smaller than the EHT beam, which we attribute to Faraday rotation internal to the emission region. We estimate the average density n e ∼ 104–7 cm−3, magnetic field strength B ∼ 1–30 G, and electron temperature T e ∼ (1–12) × 1010 K of the radiating plasma in a simple one-zone emission model. We show that the net azimuthal linear polarization pattern may result from organized, poloidal magnetic fields in the emission region. In a quantitative comparison with a large library of simulated polarimetric images from general relativistic magnetohydrodynamic (GRMHD) simulations, we identify a subset of physical models that can explain critical features of the polarimetric EHT observations while producing a relativistic jet of sufficient power. The consistent GRMHD models are all of magnetically arrested accretion disks, where near-horizon magnetic fields are dynamically important. We use the models to infer a mass accretion rate onto the black hole in M87 of (3–20) × 10−4 M ⊙ yr−1.