Search for: All records

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. 

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.PRIMOis 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.PRIMOfollows Monte Carlo Markov Chains to fit a linear combination of principal components derived from an ensemble of simulated images to interferometric data. We show thatPRIMOcan 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. 

Abstract We demonstrate the use of an eigenbasis that is derived from principal component analysis (PCA) applied on an ensemble of random-noise images that have a “red” power spectrum; i.e., a spectrum that decreases smoothly from large to small spatial scales. The pattern of the resulting eigenbasis allows for the reconstruction of images with a broad range of image morphologies. In particular, we show that this general eigenbasis can be used to efficiently reconstruct images that resemble possible astronomical sources for interferometric observations, even though the images in the original ensemble used to generate the PCA basis are significantly different from the astronomical images. We further show that the efficiency and fidelity of the image reconstructions depends only weakly on the particular parameters of the red-noise power spectrum used to generate the ensemble of images. 

Abstract We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 s on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications. 

Context.The 2017 observing campaign of the Event Horizon Telescope (EHT) delivered the first very long baseline interferometry (VLBI) images at the observing frequency of 230 GHz, leading to a number of unique studies on black holes and relativistic jets from active galactic nuclei (AGN). In total, eighteen sources were observed, including the main science targets, Sgr A* and M 87, and various calibrators. Sixteen sources were AGN. Aims.We investigated the morphology of the sixteen AGN in the EHT 2017 data set, focusing on the properties of the VLBI cores: size, flux density, and brightness temperature. We studied their dependence on the observing frequency in order to compare it with the Blandford-Königl (BK) jet model. In particular, we aimed to study the signatures of jet acceleration and magnetic energy conversion. Methods.We modeled the source structure of seven AGN in the EHT 2017 data set using linearly polarized circular Gaussian components (1749+096, 1055+018, BL Lac, J0132–1654, J0006–0623, CTA 102, and 3C 454.3) and collected results for the other nine AGN from dedicated EHT publications, complemented by lower frequency data in the 2–86 GHz range. Combining these data into a multifrequency EHT+ data set, we studied the dependences of the VLBI core component flux density, size, and brightness temperature on the frequency measured in the AGN host frame (and hence on the distance from the central black hole), characterizing them with power law fits. We compared the observations with the BK jet model and estimated the magnetic field strength dependence on the distance from the central black hole. Results.Our observations spanning event horizon to parsec scales indicate a deviation from the standard BK model, particularly in the decrease of the brightness temperature with the observing frequency. Only some of the discrepancies may be alleviated by tweaking the model parameters or the jet collimation profile. Either bulk acceleration of the jet material, energy transfer from the magnetic field to the particles, or both are required to explain the observations. For our sample, we estimate a general radial dependence of the Doppler factorδ ∝ r^≤0.5. This interpretation is consistent with a magnetically accelerated sub-parsec jet. We also estimate a steep decrease of the magnetic field strength with radiusB ∝ r⁻³, hinting at jet acceleration or efficient magnetic energy dissipation. 

Context.Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. Aims.To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. Methods.We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. Our study is focused on the size and continuum spectrum of the innermost region containing the central engine and the footpoints of both jets. We employed a synchrotron-self absorption model to fit the continuum radio spectrum and we combined the size measurements from close to the central engine out to ∼1 pc to study the jet collimation. Results.For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 43 μas perpendicular to the jet axes, corresponding to just around 250 R_S(Schwarzschild radii). This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 3 × 10³ R_Sand 1 × 10⁴ R_Sdistances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 μas from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at ∼130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 2.6 × 10⁴ Gauss, which is consistent with previous measurements. Conclusions.The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength. 

Abstract In a companion paper, we present the first spatially resolved polarized image of Sagittarius A* on event horizon scales, captured using the Event Horizon Telescope, a global very long baseline interferometric array operating at a wavelength of 1.3 mm. Here we interpret this image using both simple analytic models and numerical general relativistic magnetohydrodynamic (GRMHD) simulations. The large spatially resolved linear polarization fraction (24%–28%, peaking at ∼40%) is the most stringent constraint on parameter space, disfavoring models that are too Faraday depolarized. Similar to our studies of M87*, polarimetric constraints reinforce a preference for GRMHD models with dynamically important magnetic fields. Although the spiral morphology of the polarization pattern is known to constrain the spin and inclination angle, the time-variable rotation measure (RM) of Sgr A* (equivalent to ≈46° ± 12° rotation at 228 GHz) limits its present utility as a constraint. If we attribute the RM to internal Faraday rotation, then the motion of accreting material is inferred to be counterclockwise, contrary to inferences based on historical polarized flares, and no model satisfies all polarimetric and total intensity constraints. On the other hand, if we attribute the mean RM to an external Faraday screen, then the motion of accreting material is inferred to be clockwise, and one model passes all applied total intensity and polarimetric constraints: a model with strong magnetic fields, a spin parameter of 0.94, and an inclination of 150°. We discuss how future 345 GHz and dynamical imaging will mitigate our present uncertainties and provide additional constraints on the black hole and its accretion flow. 

Abstract The Event Horizon Telescope observed the horizon-scale synchrotron emission region around the Galactic center supermassive black hole, Sagittarius A* (Sgr A*), in 2017. These observations revealed a bright, thick ring morphology with a diameter of 51.8 ± 2.3μas and modest azimuthal brightness asymmetry, consistent with the expected appearance of a black hole with massM≈ 4 × 10⁶M_⊙. From these observations, we present the first resolved linear and circular polarimetric images of Sgr A*. The linear polarization images demonstrate that the emission ring is highly polarized, exhibiting a prominent spiral electric vector polarization angle pattern with a peak fractional polarization of ∼40% in the western portion of the ring. The circular polarization images feature a modestly (∼5%–10%) polarized dipole structure along the emission ring, with negative circular polarization in the western region and positive circular polarization in the eastern region, although our methods exhibit stronger disagreement than for linear polarization. We analyze the data using multiple independent imaging and modeling methods, each of which is validated using a standardized suite of synthetic data sets. While the detailed spatial distribution of the linear polarization along the ring remains uncertain owing to the intrinsic variability of the source, the spiraling polarization structure is robust to methodological choices. The degree and orientation of the linear polarization provide stringent constraints for the black hole and its surrounding magnetic fields, which we discuss in an accompanying publication.