skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.
Attention:The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 7:00 AM ET to 7:30 AM ET on Friday, April 24 due to maintenance. We apologize for the inconvenience.


Title: Variational Image Feature Extraction for the Event Horizon Telescope
Abstract Imaging algorithms form powerful analysis tools for very long baseline interferometry (VLBI) data analysis. However, these tools cannot measure certain image features (e.g., ring diameter) by their nonparametric nature. This is unfortunate since these image features are often related to astrophysically relevant quantities such as black hole mass. This paper details a new general image feature-extraction technique that applies to a wide variety of VLBI image reconstructions called variational image domain analysis . Unlike previous tools, variational image domain analysis can be applied to any image reconstruction regardless of its structure. To demonstrate its flexibility, we analyze thousands of reconstructions from previous Event Horizon Telescope synthetic data sets and recover image features such as diameter, orientation, and ellipticity. By measuring these features, our technique can help extract astrophysically relevant quantities such as the mass and orientation of the central black hole in M87.  more » « less
Award ID(s):
1935980
PAR ID:
10328067
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
925
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
122
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for the EHT is how well key parameters, such as trapped magnetic flux and the associated disk models, can be extracted from present and future EHT VLBI data products. The process of modeling visibilities and analyzing them is complicated by the fact that the data are sparsely sampled in the Fourier domain while most of the theory/simulation is constructed in the image domain. Here we propose a data-driven approach to analyze complex visibilities and closure quantities for radio interferometric data with neural networks. Using mock interferometric data, we show that our neural networks are able to infer the accretion state as either high magnetic flux (MAD) or low magnetic flux (SANE), suggesting that it is possible to perform parameter extraction directly in the visibility domain without image reconstruction. We have applied VLBInet to real M87 EHT data taken on four different days in 2017 (April 5, 6, 10, 11), and our neural networks give a score prediction 0.52, 0.4, 0.43, 0.76 for each day, with an average score 0.53, which shows no significant indication for the data to lean toward either the MAD or SANE state. 
    more » « less
  2. The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for the EHT is how well key parameters, such as trapped magnetic flux and the associated disk models, can be extracted from present and future EHT VLBI data products. The process of modeling visibilities and analyzing them is complicated by the fact that the data are sparsely sampled in the Fourier domain while most of the theory/simulation is constructed in the image domain. Here we propose a data- driven approach to analyze complex visibilities and closure quantities for radio interferometric data with neural networks. Using mock interferometric data, we show that our neural networks are able to infer the accretion state as either high magnetic flux (MAD) or low magnetic flux (SANE), suggesting that it is possible to perform parameter extraction directly in the visibility domain without image reconstruction. We have applied VLBInet to real M87 EHT data taken on four di↵erent days in 2017 (April 5, 6, 10, 11), and our neural networks give a score prediction 0.52, 0.4, 0.43, 0.76 for each day, with an average score 0.53, which shows no significant indication for the data to lean toward either the MAD or SANE state. 
    more » « less
  3. A bstract We construct a family of non-supersymmetric extremal black holes and their horizonless microstate geometries in four dimensions. The black holes can have finite angular momentum and an arbitrary charge-to-mass ratio, unlike their supersymmetric cousins. These features make them and their microstate geometries astrophysically relevant. Thus, they provide interesting prototypes to study deviations from Kerr solutions caused by new horizon-scale physics. In this paper, we compute the gravitational multipole structure of these solutions and compare them to Kerr black holes. The multipoles of the black hole differ significantly from Kerr as they depend non-trivially on the charge-to-mass ratio. The horizonless microstate geometries (that are comparable in size to a black hole) have a similar multipole structure as their corresponding black hole, with deviations to the black hole multipole values set by the scale of their microstructure. 
    more » « less
  4. Abstract We report results from the first Earth-space VLBI observations of the Galactic Center supermassive black hole, Sgr A*. These observations used the space telescope Spektr-R of the RadioAstron project together with a global network of 20 ground telescopes, observing at a wavelength of 1.35 cm. Spektr-R provided baselines up to 3.9 times the diameter of the Earth, corresponding to an angular resolution of approximately 55 μ as and a spatial resolution of 5.5 R Sch at the source, where R Sch ≡ 2 GM / c 2 is the Schwarzschild radius of Sgr A*. Our short ground baseline measurements ( ≲ 80 M λ ) are consistent with an anisotropic Gaussian image, while our intermediate ground baseline measurements (100–250 M λ ) confirm the presence of persistent image substructure in Sgr A*. Both features are consistent with theoretical expectations for strong scattering in the ionized interstellar medium, which produces Gaussian scatter-broadening on short baselines and refractive substructure on long baselines. We do not detect interferometric fringes on any of the longer ground baselines or on any ground–space baselines. While space-VLBI offers a promising pathway to sharper angular resolution and the measurement of key gravitational signatures in black holes, such as their photon rings, our results demonstrate that space-VLBI studies of Sgr A* will require sensitive observations at submillimeter wavelengths. 
    more » « less
  5. Abstract The sparse interferometric coverage of the Event Horizon Telescope (EHT) poses a significant challenge for both reconstruction and model fitting of black hole images. PRIMO is a new principal components analysis-based algorithm for image reconstruction that uses the results of high-fidelity general relativistic, magnetohydrodynamic simulations of low-luminosity accretion flows as a training set. This allows the reconstruction of images that are consistent with the interferometric data and that live in the space of images that is spanned by the simulations. PRIMO follows Monte Carlo Markov Chains to fit a linear combination of principal components derived from an ensemble of simulated images to interferometric data. We show that PRIMO can efficiently and accurately reconstruct synthetic EHT data sets for several simulated images, even when the simulation parameters are significantly different from those of the image ensemble that was used to generate the principal components. The resulting reconstructions achieve resolution that is consistent with the performance of the array and do not introduce significant biases in image features such as the diameter of the ring of emission. 
    more » « less