Publications

Light bending by the strong gravity around the black hole will form the so-called black hole shadow, the shape of which can shed light on the structure of the near-horizon geometry to possibly reveal novel physics of strong gravity and black hole. In this work, we adopt both analytical and ray-tracing methods to study the black hole shadow in the presence of the infrared structure of gravity theory, which manifests the asymptotic symmetries of spacetime as the supertranslation soft hairs of the black hole. Though the black hole metrics with and without the soft hair are related by large gauge transformations, the near horizon geometries relevant for the shape of the shadow are quite different. Moreover, the Hamiltonian for the geodesic seems intrinsically different, i.e., the loss of separability due to the breaking of spherical symmetry by soft hair. By applying ray-tracing computations, we find that the soft hair, although not affecting the shape of the shadow, may change the average size and position of the shadow. Images resulting from soft hair black holes with surrounding accretion flows are also discussed.
 

Quadratic gravity constitutes a prototypical example of a perturbatively renormalizable quantum theory of the gravitational interactions. In this work, we construct the associated phase space of static, spherically symmetric, and asymptotically flat spacetimes. It is found that the Schwarzschild geometry is embedded in a rich solution space comprising horizonless, naked singularities and wormhole solutions. Characteristically, the deformed solutions follow the Schwarzschild solution up outside of the photon sphere while they differ substantially close to the center of gravity. We then carry out an analytic analysis of observable signatures accessible to the Event Horizon Telescope, comprising the size of the black hole shadow as well as the radiation emitted by infalling matter. On this basis, we argue that it is the brightness within the shadow region which constrains the phase space of solutions. Our work constitutes the first step towards bounding the phase space of black hole type solutions with a clear quantum gravity interpretation based on observational data.

The Event Horizon Telescope (EHT) collaboration has produced the first resolved images of M87*, the supermassive black hole at the centre of the elliptical galaxy M87. As both technology and analysis pipelines improve, it will soon become possible to produce spectral index maps of black hole accretion flows on event horizon scales. In this work, we predict spectral index maps of both M87* and Sgr A* by applying the general relativistic radiative transfer (GRRT) code {\sc ipole} to a suite of general relativistic magnetohydrodynamic (GRMHD) simulations. We analytically explore how the spectral index increases with increasing magnetic field strength, electron temperature, and optical depth. Consequently, spectral index maps grow more negative with increasing radius in almost all models, since all of these quantities tend to be maximized near the event horizon. Additionally, photon ring geodesics exhibit more positive spectral indices, since they sample the innermost regions of the accretion flow with the most extreme plasma conditions. Spectral index maps are sensitive to highly uncertain plasma heating prescriptions (the electron temperature and distribution function). However, if our understanding of these aspects of plasma physics can be tightened, even the spatially unresolved spectral index around 230 GHz can be used to discriminate between models. In particular, Standard and Normal Evolution (SANE) flows tend to exhibit more negative spectral indices than Magnetically Arrested Disk (MAD) flows due to differences in the characteristic magnetic field strength and temperature of emitting plasma.

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.

Inference is crucial in modern astronomical research, where hidden astrophysical features and patterns are often estimated from indirect and noisy measurements. Inferring the posterior of hidden features, conditioned on the observed measurements, is essential for understanding the uncertainty of results and downstream scientific interpretations. Traditional approaches for posterior estimation include sampling-based methods and variational inference. However, sampling-based methods are typically slow for high-dimensional inverse problems, while variational inference often lacks estimation accuracy. In this paper, we propose alpha-DPI, a deep learning framework that first learns an approximate posterior using alpha-divergence variational inference paired with a generative neural network, and then produces more accurate posterior samples through importance re-weighting of the network samples. It inherits strengths from both sampling and variational inference methods: it is fast, accurate, and scalable to high-dimensional problems. We apply our approach to two high-impact astronomical inference problems using real data: exoplanet astrometry and black hole feature extraction.

The black hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5–61 days) is comparable to the 6 day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure-phase measurements on all six linearly independent nontrivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of ∼3°–5°. The only triangles that exhibit substantially higher variability (∼90°–180°) are the ones with baselines that cross the visibility amplitude minima on the u–v plane, as expected from theoretical modeling. We used two sets of general relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas.

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 the PATOKA pipeline, 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, and igrmonty, 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 by PATOKA and conclude with a discussion of limitations and future directions.

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.

The Event Horizon Telescope recently captured images of the supermassive black hole in the center of the M87 galaxy, which shows a ring-like emission structure with the south side only slightly brighter than the north side. This relatively weak asymmetry in the brightness profile along the ring has been interpreted as a consequence of the low inclination of the observer (around 17° for M87), which suppresses the Doppler beaming and boosting effects that might otherwise be expected due to the nearly relativistic velocities of the orbiting plasma. In this work, we use a large suite of general relativistic magnetohydrodynamic simulations to reassess the validity of this argument. By constructing explicit counterexamples, we show that low inclination is a sufficient but not necessary condition for images to have low brightness asymmetry. Accretion flow models with high accumulated magnetic flux close to the black hole horizon (the so-called magnetically arrested disks) and low black hole spins have angular velocities that are substantially smaller than the orbital velocities of test particles at the same location. As a result, such models can produce images with low brightness asymmetry even when viewed edge on.

Recent developments in very long baseline interferometry (VLBI) have made it possible for the Event Horizon Telescope (EHT) to resolve the innermost accretion flows of the largest supermassive black holes on the sky. The sparse nature of the EHT's (u, v)-coverage presents a challenge when attempting to resolve highly time-variable sources. We demonstrate that the changing (u, v)-coverage of the EHT can contain regions of time over the course of a single observation that facilitate dynamical imaging. These optimal time regions typically have projected baseline distributions that are approximately angularly isotropic and radially homogeneous. We derive a metric of coverage quality based on baseline isotropy and density that is capable of ranking array configurations by their ability to produce accurate dynamical reconstructions. We compare this metric to existing metrics in the literature and investigate their utility by performing dynamical reconstructions on synthetic data from simulated EHT observations of sources with simple orbital variability. We then use these results to make recommendations for imaging the 2017 EHT Sgr A* data set.