Principal-Component Interferometric Modeling (PRIMO), an Algorithm for EHT Data I: Reconstructing Images from Simulated EHT Observations

TL;DR

PRIMO is a novel principal components analysis-based algorithm that reconstructs black-hole images from EHT data by leveraging simulation-based training, effectively handling sparse coverage and reducing biases.

Contribution

It introduces a new image reconstruction method that uses simulation-derived principal components and Monte Carlo Markov Chains for accurate EHT data analysis.

Findings

Successfully reconstructs synthetic EHT data with high fidelity.

Achieves resolution consistent with array performance.

Reduces biases in key image features.

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 both 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.