Abstract A new image-reconstruction algorithm, Principal-component Interferometric Modeling (PRIMO), applied to the interferometric data of the M87 black hole collected with the Event Horizon Telescope (EHT), resulted in an image that reached the native resolution of the telescope array.PRIMOis based on learning a compact set of image building blocks obtained from a large library of high-fidelity, physics-based simulations of black hole images. It uses these building blocks to fill the sparse Fourier coverage of the data that results from the small number of telescopes in the array. In this paper, we show that this approach is readily justified. Since the angular extent of the image of the black hole and of its inner accretion flow is finite, the Fourier space domain is heavily smoothed, with a correlation scale that is at most comparable to the sizes of the data gaps in the coverage of Fourier space with the EHT. Consequently,PRIMOor other machine learning algorithms can faithfully reconstruct the images without the need to generate information that is unconstrained by the data within the resolution of the array. We also address the completeness of the eigenimages and the compactness of the resulting representation. We show thatPRIMOprovides a compact set of eigenimages that have sufficient complexity to recreate a broad set of images well beyond those in the training set.
more »
« less
PATOKA: Simulating Electromagnetic Observables of Black Hole Accretion
Abstract The Event Horizon Telescope (EHT) has released analyses of reconstructed images of horizon-scale millimeter emission near the supermassive black hole at the center of the M87 galaxy. Parts of the analyses made use of a large library of synthetic black hole images and spectra, which were produced using numerical general relativistic magnetohydrodynamics fluid simulations and polarized ray tracing. In this article, we describe thePATOKApipeline, which was used to generate the Illinois contribution to the EHT simulation library. We begin by describing the relevant accretion systems and radiative processes. We then describe the details of the three numerical codes we use,iharm,ipole, andigrmonty, paying particular attention to differences between the current generation of the codes and the originally published versions. Finally, we provide a brief overview of simulated data as produced byPATOKAand conclude with a discussion of limitations and future directions.
more »
« less
- PAR ID:
- 10486420
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Supplement Series
- Volume:
- 259
- Issue:
- 2
- ISSN:
- 0067-0049
- Format(s):
- Medium: X Size: Article No. 64
- Size(s):
- Article No. 64
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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 theeht-imagingPython 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.more » « less
-
Abstract The Event Horizon Telescope (EHT) has produced images of M87* and Sagittarius A*, and will soon produce time sequences of images, or movies. In anticipation of this, we describe a technique to measure the rotation rate, or pattern speed Ωp, from movies using an autocorrelation technique. We validate the technique on Gaussian random field models with a known rotation rate and apply it to a library of synthetic images of Sgr A* based on general relativistic magnetohydrodynamics simulations. We predict that EHT movies will have Ωp≈ 1° perGMc−3, which is of order 15% of the Keplerian orbital frequency in the emitting region. We can plausibly attribute the slow rotation seen in our models to the pattern speed of inward-propagating spiral shocks. We also find that Ωpdepends strongly on inclination. Application of this technique will enable us to compare future EHT movies with the clockwise rotation of Sgr A* seen in near-infrared flares by GRAVITY. Pattern speed analysis of future EHT observations of M87* and Sgr A* may also provide novel constraints on black hole inclination and spin, as well as an independent measurement of black hole mass.more » « less
-
Abstract The Event Horizon Telescope (EHT) images of the supermassive black hole at the center of the galaxy M87 provided the first image of the accretion environment on horizon scales. General relativity (GR) predicts that the image of the shadow should be nearly circular given the inclination angle of the black hole M87*. A robust detection of ellipticity in image reconstructions of M87* could signal new gravitational physics on horizon scales. Here we analyze whether the imaging parameters used in EHT analyses are sensitive to ring ellipticity, and measure the constraints on the ellipticity of M87*. We find that the top set is unable to recover ellipticity. Even for simple geometric models, the true ellipticity is biased low, preferring circular rings. Therefore, to place a constraint on the ellipticity of M87*, we measure the ellipticity of 550 synthetic data sets produced from GRMHD simulations. We find that images with intrinsic axis ratios of 2:1 are consistent with the ellipticity seen from EHT image reconstructions.more » « less
-
Abstract The Event Horizon Telescope (EHT) recently released the first linearly polarized images of the accretion flow around the supermassive black hole Messier 87*, hereafter M87*. The spiraling polarization pattern found in the EHT images favored magnetically arrested disks as the explanation for the EHT image. With next-generation improvements to very long baseline interferometry on the horizon, understanding similar polarized features in the highly lensed structure known as the “photon ring,” where photons make multiple half orbits about the black hole before reaching the observer, will be critical to the analysis of future images. Recent work has indicated that this image region may be depolarized relative to more direct emission. We expand this observation by decomposing photon half orbits in the EHT library of simulated images of the M 87* accretion system and find that images of magnetically arrested disk simulations show a relative depolarization of the photon ring attributable to destructive interference of oppositely spiraling electric field vectors; this antisymmetry, which arises purely from strong gravitational lensing, can produce up to ∼50% depolarization in the photon ring region with respect to the direct image. In systems that are not magnetically arrested and with the exception of systems with high spin and ions and electrons of equal temperature, we find that highly lensed indirect subimages are almost completely depolarized, causing a modest depolarization of the photon ring region in the complete image. We predict that next-generation EHT observations of M 87* polarization should jointly constrain the black hole spin and the underlying emission and magnetic field geometry.more » « less