Probing the Black Hole Metric. I. Black Hole Shadows and Binary Black-Hole Inspirals

In General Relativity, the spacetimes of black holes have three fundamental properties: (i) they are the same, to lowest order in spin, as the metrics of stellar objects; (ii) they are independent of mass, when expressed in geometric units; and (iii) they are described by the Kerr metric. In this paper, we quantify the upper bounds on potential black-hole metric deviations imposed by observations of black-hole shadows and of binary black-hole inspirals in order to explore the current experimental limits on possible violations of the last two predictions. We find that both types of experiments provide correlated constraints on deviation parameters that are primarily in the tt-components of the spacetimes, when expressed in areal coordinates. We conclude that, currently, there is no evidence for a deviations from the Kerr metric across the 8 orders of magnitudes in masses and 16 orders in curvatures spanned by the two types of black holes. Moreover, because of the particular masses of black holes in the current sample of gravitational-wave sources, the correlations imposed by the two experiments are aligned and of similar magnitudes when expressed in terms of the far field, post-Newtonian predictions of the metrics. If a future coalescing black-hole binary with two low-mass (e.g., ~3 Msun) components is discovered, the degeneracy between the deviation parameters can be broken by combining the inspiral constraints with those from the black-hole shadow measurements.

Gravitational Test Beyond the First Post-Newtonian Order with the Shadow of the M87 Black Hole

The 2017 Event Horizon Telescope (EHT) observations of the central source in M87 have led to the first measurement of the size of a black-hole shadow. This observation offers a new and clean gravitational test of the black-hole metric in the strong-field regime. We show analytically that spacetimes that deviate from the Kerr metric but satisfy weak-field tests can lead to large deviations in the predicted black-hole shadows that are inconsistent with even the current EHT measurements. We use numerical calculations of regular, parametric, non-Kerr metrics to identify the common characteristic among these different parametrizations that control the predicted shadow size. We show that the shadow-size measurements place significant constraints on deviation parameters that control the second post-Newtonian and higher orders of each metric and are, therefore, inaccessible to weak-field tests. The new constraints are complementary to those imposed by observations of gravitational waves from stellar-mass sources.

Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

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

Discretization and Filtering Effects on Black Hole Images Obtained with the Event Horizon Telescope

Interferometers, such as the Event Horizon Telescope (EHT), do not directly observe the images of sources but rather measure their Fourier components at discrete spatial frequencies up to a maximum value set by the longest baseline in the array. Construction of images from the Fourier components or analysis of them with high-resolution models requires careful treatment of fine source structure nominally beyond the array resolution. The primary EHT targets, Sgr A* and M87, are expected to have black-hole shadows with sharp edges and strongly filamentary emission from the surrounding plasma on scales much smaller than those probed by the currently largest baselines. We show that for aliasing not to affect images reconstructed with regularized maximum likelihood methods and model images that are directly compared to the data, the sampling of these images (i.e., their pixel spacing) needs to be significantly finer than the scale probed by the largest baseline in the array. Using GRMHD simulations of black-hole images, we estimate the maximum allowable pixel spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT targets, this corresponds to an angular pixel size of <0.5 microarcseconds. With aliasing under control, we then advocate use of the second-order Butterworth filter with a cut-off scale equal to the maximum array baseline as optimal for visualizing the reconstructed images. In contrast to the traditional Gaussian filters, this Butterworth filter retains most of the power at the scales probed by the array while suppressing the fine image details for which no data exist.

The Scattering and Intrinsic Structure of Sagittarius A* at Radio Wavelengths

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 pm 0.013) lambda_{rm cm}^2,{rm mas}$ along the major axis, $(0.703 pm 0.013) lambda_{rm cm}^2,{rm mas}$ along the minor axis, and a position angle of $81.9^circ pm 0.2^circ$. We estimate that the turbulent dissipation scale is at least $600,{rm km}$, with tentative support for $r_{rm in} = 800 pm 200,{rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($beta=11/3$), and we estimate $beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{rm in} = 800,{rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

Testing General Relativity with the Event Horizon Telescope

The Event Horizon Telescope is a millimeter VLBI array that aims to take the first pictures of the black holes in the center of the Milky Way and of the M87 galaxy, with horizon scale resolution. Measurements of the shape and size of the shadows cast by the black holes on the surrounding emission can test the cosmic censorship conjecture and the no-hair theorem and may find evidence for classical effects of the quantum structure of black holes. Observations of coherent structures in the accretion flows may lead to accurate measurements of the spins of the black holes and of other properties of their spacetimes. For Sgr A*, the black hole in the center of the Milky Way, measurements of the precession of stellar orbits and timing monitoring of orbiting pulsars offer complementary avenues to the gravitational tests with the Event Horizon Telescope.

A Model for Anisotropic Interstellar Scattering and its Application to Sgr A*

Scattering in the ionized interstellar medium is commonly observed to be anisotropic, with theories of magnetohydrodynamic (MHD) turbulence explaining the anisotropy through a preferred magnetic field direction throughout the scattering regions. In particular, the line of sight to the Galactic Center supermassive black hole, Sgr A*, exhibits strong and anisotropic scattering, which dominates its observed size at wavelengths of a few millimeters and longer. Therefore, inferences of the intrinsic structure of sgra at these wavelengths are sensitive to the assumed scattering model. In addition, extrapolations of the scattering model from long wavelengths, at which its parameters are usually estimated, to 1.3 mm, where the Event Horizon Telescope (EHT) seeks to image Sgr A* on Schwarzschild-radius scales, are also sensitive to the assumed scattering model. Past studies of Sgr A* have relied on simple Gaussian models for the scattering kernel that effectively presume an inner scale of turbulence far greater than the diffractive scale; this assumption is likely violated for Sgr A* at 1.3 mm. We develop a physically motivated model for anisotropic scattering, using a simplified model for MHD turbulence with a finite inner scale and a wandering transverse magnetic field direction. We explore several explicit analytic models for this wandering and derive the expected observational properties --- scatter broadening and refractive scintillation --- for each. For expected values of the inner scale, the scattering kernel for all models is markedly non-Gaussian at 1.3 mm but is straightforward to calculate and depends only weakly on the assumed model for the wandering of the magnetic field direction. On the other hand, in all models, the refractive substructure depends strongly on the wandering model and may be an important consideration in imaging Sgr A* with the EHT.

Event Horizon Telescope Observations as Probes for Quantum Structure of Astrophysical Black Holes

The need for a consistent quantum evolution for black holes has led to proposals that their semiclassical description is modified not just near the singularity, but at horizon or larger scales. If such modifications extend beyond the horizon, they influence regions accessible to distant observeration. Natural candidates for these modifications behave like metric fluctuations, with characteristic length and time scales set by the horizon radius. We investigate the possibility of using the Event Horizon Telescope to observe these effects, if they have a strength sufficient to make quantum evolution consistent with unitarity. We find that such quantum fluctuations can introduce a strong time dependence for the shape and size of the shadow that a black hole casts on its surrounding emission. For the black hole in the center of the Milky Way, detecting the rapid time variability of its shadow will require non-imaging timing techniques. However, for the much larger black hole in the center of the M87 galaxy, a variable black-hole shadow, if present with these parameters, would be readily observable in the individual snapshots that will be obtained by the Event Horizon Telescope.

Model-Independent Inference of Neutron Star Radii from Moment of Inertia Measurements

A precise moment of inertia measurement for PSR J0737-3039A in the double pulsar system is expected within the next five years. We present here a new method of mapping the anticipated measurement of the moment of inertia directly into the neutron star structure. We determine the maximum and minimum values possible for the moment of inertia of a neutron star of a given radius based on physical stability arguments, assuming knowledge of the equation of state only at densities below the nuclear saturation density. If the equation of state is trusted up to the nuclear saturation density, we find that a measurement of the moment of inertia will place absolute bounds on the radius of PSR J0737-3039A to within $pm$1 km. The resulting combination of moment of inertia, mass, and radius measurements for a single source will allow for new, stringent constraints on the dense-matter equation of state.

A Quantitative Test of the No-Hair Theorem with Sgr A* using stars, pulsars, and the Event Horizon Telescope

The black hole in the center of the Milky Way, Sgr A*, has the largest mass-to-distance ratio among all known black holes in the Universe. This property makes Sgr A* the optimal target for testing the gravitational no-hair theorem. In the near future, major developments in instrumentation will provide the tools for high-precision studies of its spacetime via observations of relativistic effects in stellar orbits, in the timing of pulsars, and in horizon-scale images of its accretion flow. We explore here the prospect of measuring the properties of the black-hole spacetime using all these three types of observations. We show that the correlated uncertainties in the measurements of the black-hole spin and quadrupole moment using the orbits of stars and pulsars are nearly orthogonal to those obtained from measuring the shape and size of the shadow the black hole casts on the surrounding emission. Combining these three types of observations will, therefore, allow us to assess and quantify systematic biases and uncertainties in each measurement and lead to a highly accurate, quantitative test of the gravitational no-hair theorem.