The Event Horizon Telescope (EHT) has produced images of the plasma flow around the supermassive black holes in Sgr A* and M87* with a resolution comparable to the projected size of their event horizons. Observations with the next-generation Event Horizon Telescope (ngEHT) will have significantly improved Fourier plane coverage and will be conducted at multiple frequency bands (86, 230, and 345 GHz), each with a wide bandwidth. At these frequencies, both Sgr A* and M87* transition from optically thin to optically thick. Resolved spectral index maps in the near-horizon and jet-launching regions of these supermassive black hole sources can constrain properties of the emitting plasma that are degenerate in single-frequency images. In addition, combining information from data obtained at multiple frequencies is a powerful tool for interferometric image reconstruction, since gaps in spatial scales in single-frequency observations can be filled in with information from other frequencies. Here we present a new method of simultaneously reconstructing interferometric images at multiple frequencies along with their spectral index maps. The method is based on existing regularized maximum likelihood (RML) methods commonly used for EHT imaging and is implemented in the
VLBI imaging of black holes via second moment regularization
The imaging fidelity of the Event Horizon Telescope (EHT) is currently determined by its sparse baseline coverage. In particular, EHT coverage is dominated by long baselines, and is highly sensitive to atmospheric conditions and loss of sites between experiments. The limited short/mid-range baselines especially affect the imaging process, hindering the recovery of more extended features in the image. We present an algorithmic contingency for the absence of well-constrained short baselines in the imaging of compact sources, such as the supermassive black holes observed with the EHT. This technique enforces a specific second moment on the reconstructed image in the form of a size constraint, which corresponds to the curvature of the measured visibility function at zero baseline. The method enables the recovery of information lost in gaps of the baseline coverage on short baselines and enables corrections of any systematic amplitude offsets for the stations giving short-baseline measurements present in the observation. The regularization can use historical source size measurements to constrain the second moment of the reconstructed image to match the observed size. We additionally show that a characteristic size can be derived from available short-baseline measurements, extrapolated from other wavelengths, or estimated without complementary size constraints with parameter searches. We demonstrate the capabilities of this method for both static and movie reconstructions of variable sources.
more »
« less
- Award ID(s):
- 1716536
- NSF-PAR ID:
- 10156034
- Date Published:
- Journal Name:
- Astronomy & Astrophysics
- Volume:
- 629
- ISSN:
- 0004-6361
- Page Range / eLocation ID:
- A32
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract eht-imaging Python software library. We show results of this method on simulated ngEHT data sets as well as on real data from the Very Long Baseline Array and Atacama Large Millimeter/submillimeter Array. These examples demonstrate that simultaneous RML multifrequency image reconstruction produces higher-quality and more scientifically useful results than is possible from combining independent image reconstructions at each frequency. -
Abstract The extraordinary physical resolution afforded by the Event Horizon Telescope has opened a window onto the astrophysical phenomena unfolding on horizon scales in two known black holes, M87 * and Sgr A*. However, with this leap in resolution has come a new set of practical complications. Sgr A* exhibits intraday variability that violates the assumptions underlying Earth aperture synthesis, limiting traditional image reconstruction methods to short timescales and data sets with very sparse ( u , v ) coverage. We present a new set of tools to detect and mitigate this variability. We develop a data-driven, model-agnostic procedure to detect and characterize the spatial structure of intraday variability. This method is calibrated against a large set of mock data sets, producing an empirical estimator of the spatial power spectrum of the brightness fluctuations. We present a novel Bayesian noise modeling algorithm that simultaneously reconstructs an average image and statistical measure of the fluctuations about it using a parameterized form for the excess variance in the complex visibilities not otherwise explained by the statistical errors. These methods are validated using a variety of simulated data, including general relativistic magnetohydrodynamic simulations appropriate for Sgr A* and M87 * . We find that the reconstructed source structure and variability are robust to changes in the underlying image model. We apply these methods to the 2017 EHT observations of M87 * , finding evidence for variability across the EHT observing campaign. The variability mitigation strategies presented are widely applicable to very long baseline interferometry observations of variable sources generally, for which they provide a data-informed averaging procedure and natural characterization of inter-epoch image consistency.more » « less
-
more » « less
In April 2019, the Event Horizon Telescope (EHT) Collaboration reported the first-ever event-horizon-scale images of a black hole, resolving the central compact radio source in the giant elliptical galaxy M 87. These images reveal a ring with a southerly brightness distribution and a diameter of ∼42 μas, consistent with the predicted size and shape of a shadow produced by the gravitationally lensed emission around a supermassive black hole. These results were obtained as part of the April 2017 EHT observation campaign, using a global very long baseline interferometric radio array operating at a wavelength of 1.3 mm. Here, we present results based on the second EHT observing campaign, taking place in April 2018 with an improved array, wider frequency coverage, and increased bandwidth. In particular, the additional baselines provided by the Greenland telescope improved the coverage of the array. Multiyear EHT observations provide independent snapshots of the horizon-scale emission, allowing us to confirm the persistence, size, and shape of the black hole shadow, and constrain the intrinsic structural variability of the accretion flow. We have confirmed the presence of an asymmetric ring structure, brighter in the southwest, with a median diameter of 43.3−3.1+1.5 μas. The diameter of the 2018 ring is remarkably consistent with the diameter obtained from the previous 2017 observations. On the other hand, the position angle of the brightness asymmetry in 2018 is shifted by about 30° relative to 2017. The perennial persistence of the ring and its diameter robustly support the interpretation that the ring is formed by lensed emission surrounding a Kerr black hole with a mass ∼6.5 × 109
M ⊙. The significant change in the ring brightness asymmetry implies a spin axis that is more consistent with the position angle of the large-scale jet. -
General relativity predicts that images of optically thin accretion flows around black holes should generically have a “photon ring”, composed of a series of increasingly sharp subrings that correspond to increasingly strongly lensed emission near the black hole. Because the effects of lensing are determined by the spacetime curvature, the photon ring provides a pathway to precise measurements of the black hole properties and tests of the Kerr metric. We explore the prospects for detecting and measuring the photon ring using very long baseline interferometry (VLBI) with the Event Horizon Telescope (EHT) and the next-generation EHT (ngEHT). We present a series of tests using idealized self-fits to simple geometrical models and show that the EHT observations in 2017 and 2022 lack the angular resolution and sensitivity to detect the photon ring, while the improved coverage and angular resolution of ngEHT at 230 GHz and 345 GHz is sufficient for these models. We then analyze detection prospects using more realistic images from general relativistic magnetohydrodynamic simulations by applying “hybrid imaging”, which simultaneously models two components: a flexible raster image (to capture the direct emission) and a ring component. Using the Bayesian VLBI modeling package Comrade.jl, we show that the results of hybrid imaging must be interpreted with extreme caution for both photon ring detection and measurement—hybrid imaging readily produces false positives for a photon ring, and its ring measurements do not directly correspond to the properties of the photon ring.more » « less
-
Abstract In this paper we quantify the temporal variability and image morphology of the horizon-scale emission from Sgr A*, as observed by the EHT in 2017 April at a wavelength of 1.3 mm. We find that the Sgr A* data exhibit variability that exceeds what can be explained by the uncertainties in the data or by the effects of interstellar scattering. The magnitude of this variability can be a substantial fraction of the correlated flux density, reaching ∼100% on some baselines. Through an exploration of simple geometric source models, we demonstrate that ring-like morphologies provide better fits to the Sgr A* data than do other morphologies with comparable complexity. We develop two strategies for fitting static geometric ring models to the time-variable Sgr A* data; one strategy fits models to short segments of data over which the source is static and averages these independent fits, while the other fits models to the full data set using a parametric model for the structural variability power spectrum around the average source structure. Both geometric modeling and image-domain feature extraction techniques determine the ring diameter to be 51.8 ± 2.3 μ as (68% credible intervals), with the ring thickness constrained to have an FWHM between ∼30% and 50% of the ring diameter. To bring the diameter measurements to a common physical scale, we calibrate them using synthetic data generated from GRMHD simulations. This calibration constrains the angular size of the gravitational radius to be 4.8 − 0.7 + 1.4 μ as, which we combine with an independent distance measurement from maser parallaxes to determine the mass of Sgr A* to be 4.0 − 0.6 + 1.1 × 10 6 M ⊙ .more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1051/0004-6361/201936156
