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
This content will become publicly available on July 1, 2027
Multifrequency Synthesis via CHIBI: Colorful Hierarchical Interferometric Bayesian Imaging
Abstract From magnetized plasma of relativistic jets to dust grains within protoplanetary disks, we study the emission mechanisms of radio sources via their rich spectral structure. Multifrequency synthesis (MFS) is a technique in which interferometric data at multiple frequencies are imaged simultaneously, resulting in a denser sampling of spatial scales, higher imaging fidelity, and tighter constraints on the source’s spectral structure and evolution. We describe a new method of MFS imaging reconstruction in a hierarchical interferometric Bayesian inference framework, CHIBI. The model parameterization is based on the spectral behavior of synchrotron radiation, the emission mechanism dominating the radio emission observed from galactic nuclei. We show results of this method on observations of jet sources from the MOJAVE catalog with the Very Long Baseline Array, and showcase the prospects for MFS imaging of M87* with simulated data from the Event Horizon Telescope (EHT) and future expansions such as the next-generation EHT and the black hole Explorer. These demonstrations highlight the benefit of MFS to reconstruct higher-fidelity images and spectral index maps, producing scientifically richer results in a statistically grounded framework, implemented inComrade.jl.
more »
« less
- PAR ID:
- 10701100
- Publisher / Repository:
- IOP
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 1005
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 132
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Abstract Radio observations provide powerful diagnostics of energy release, particle acceleration, and transport processes in solar flares. However, despite recent progress in radio interferometric imaging spectroscopy, current instruments still face limitations in image fidelity and resolution, restricting detailed spectroscopic studies of flaring regions. Here we present high-fidelity imaging spectroscopy of an M1.3 GOES class flare with MeerKAT, a precursor to the future-generation array SKA-Mid. Radio emissions at the observed frequencies typically originate in the low corona, offering valuable insights into magnetic reconnection and primary energy-release sites. The obtained images achieve an unprecedented dynamic range exceeding 103, enabling simultaneous analysis of bright coherent bursts and faint incoherent emission from the active region. Multiple spatially distinct coherent sources are identified, implying contributions from different populations of accelerated electrons. The incoherent emission extends beyond Atmospheric Imaging Assembly structures, highlighting MeerKAT’s ability to detect dilute but hot plasma invisible to extreme-ultraviolet instruments. Combined with cotemporal hard X-ray images and magnetic field extrapolations, the radio sources are located within distinct magnetic structures, further revealing their association with different populations of accelerated electrons. These results demonstrate MeerKAT imaging spectroscopy as a powerful diagnostic of solar flares and pave the way for future solar flare studies with SKA-Mid.more » « less
-
Abstract Very long baseline interferometry (VLBI) achieves the highest angular resolution in astronomy. VLBI measures corrupted Fourier components, known as visibilities. Reconstructing on-sky images from these visibilities is a challenging inverse problem, particularly for sparse arrays such as the Event Horizon Telescope (EHT) and the Very Long Baseline Array, where incomplete sampling and severe calibration errors introduce significant uncertainty in the image. To help guide convergence and control the uncertainty in image reconstructions, regularization on the space of images is utilized, such as enforcing smoothness or similarity to a fiducial image. Coupled with this regularization is the introduction of a new set of parameters that modulate its strength. We present a hierarchical Bayesian imaging approach (hierarchical interferometric Bayesian Imaging, HIBI) that enables the quantification of uncertainty for all parameters. Incorporating instrumental effects within HIBI is straightforward, allowing for simultaneous imaging and calibration of data. To showcase HIBI’s effectiveness and flexibility, we build a simple imaging model based on Markov random fields and demonstrate how different physical components can be included, e.g., black hole shadow size, and their uncertainties can be inferred. For example, while the original EHT publications were unable to constrain the ring width of M87*, HIBI measures a width of 9.3 ± 1.3μas. We apply HIBI to image and calibrate EHT synthetic data, real EHT observations of M87*, and multifrequency observations of OJ 287. Across these tests, HIBI accurately recovers a wide variety of image structures and quantifies their uncertainties. HIBI is publicly available in theComradeVLBI software repository.more » « less
-
ABSTRACT High-fidelity radio interferometric data calibration that minimizes spurious spectral structure in the calibrated data is essential in astrophysical applications, such as 21 cm cosmology, which rely on knowledge of the relative spectral smoothness of distinct astrophysical emission components to extract the signal of interest. Existing approaches to radio interferometric calibration have been shown to impart spurious spectral structure to the calibrated data if the sky model used to calibrate the data is incomplete. In this paper, we introduce BayesCal: a novel solution to the sky-model incompleteness problem in interferometric calibration, designed to enable high-fidelity data calibration. The BayesCal data model supplements the a priori known component of the forward model of the sky with a statistical model for the missing and uncertain flux contribution to the data, constrained by a prior on the power in the model. We demonstrate how the parameters of this model can be marginalized out analytically, reducing the dimensionality of the parameter space to be sampled from and allowing one to sample directly from the posterior probability distribution of the calibration parameters. Additionally, we show how physically motivated priors derived from theoretical and measurement-based constraints on the spectral smoothness of the instrumental gains can be used to constrain the calibration solutions. In a companion paper, we apply this algorithm to simulated observations with a HERA-like array and demonstrate that it enables up to four orders of magnitude suppression of power in spurious spectral fluctuations relative to standard calibration approaches.more » « less
An official website of the United States government
