The iron line and the continuum-fitting methods are currently the two leading techniques for measuring black hole spins with electromagnetic radiation. They can be naturally extended for probing the spacetime geometry around black holes and testing general relativity in the strong field regime. In the past couple of years, there has been significant work to use the iron line method to test the nature of black holes. In this Letter, we use the continuum-fitting method and we show its capability of constraining the spacetime geometry around black holes by analyzing 17 RXTE data of the X-ray binary LMC X-1.

A fraction of the electromagnetic radiation emitted from the surface of a geometrically thin and optically thick accretion disk of a black hole returns to the disk because of the strong light bending in the vicinity of the compact object (returning radiation). While such radiation clearly affects the observed spectrum of the source, it is often neglected in theoretical models. In the present paper, we study the impact of the returning radiation on relativistic reflection spectra. Assuming neutral material in the disk, we estimate the systematic uncertainties on the measurement of the properties of the system when we fit the data with a theoretical model that neglects the returning radiation. Our simulations show that the inclination angle of the disk and the black hole spin parameter tend to be overestimated for low viewing angles, while no clear bias is observed for high viewing angles. The iron abundance of the disk is never overestimated. In the most extreme cases (in particular, for maximally rotating black holes) the returning radiation flattens the radial emissivity beyond a few gravitational radii. In such cases, it also produces residuals that cannot be compensated by adjusting the parameters of models neglecting the returning radiation. This may be an important issue for interpretation of data from future X-ray missions.

The analysis of the thermal spectrum of geometrically thin and optically thick accretion disks of black holes, the so-called continuum-fitting method, is one of the leading techniques for measuring black hole spins. Current models normally approximate the disk as infinitesimally thin, while in reality the disk thickness is finite and increases as the black hole mass accretion rate increases. Here we present an XSPEC model to calculate the multi-temperature blackbody spectrum of a thin accretion disk of finite thickness around a Kerr black hole. We test our new model with an RXTE observation of the black hole binary GRS 1915+105. We find that the spin value inferred with the new model is slightly higher than the spin value obtained with a model with an infinitesimally thin disk, but the difference is small and the effect is currently subdominant with respect to other sources of uncertainties in the final spin measurement.

X-ray reflection spectroscopy (or iron line method) is a powerful tool to probe the strong gravity region of black holes, and currently is the only technique for measuring the spin of the supermassive ones. While all the available relativistic reflection models assume thin accretion disks, we know that several sources accrete near or above the Eddington limit and therefore must have thick accretion disks. In this work, we employ the Polish donut model for the description of thick disks. We thus estimate the systematic error on the spin measurement when a source with a thick accretion disk is fitted with a thin disk model. Our results clearly show that spin measurements can be significantly affected by the morphology of the accretion disk. Current spin measurements of sources with high mass accretion rate are therefore not reliable.

Asymptotically safe quantum gravity is a promising candidate scenario to provide a UV extension for the effective quantum field theory of Einstein's gravity. The theory has its foundations on the very successful framework of quantum field theory, which has been extensively tested for electromagnetic and nuclear interactions. However, observational tests of asymptotically safe quantum gravity are more challenging. Recently, a rotating black hole metric inspired by asymptotically safe quantum gravity has been proposed, and this opens the possibility of astrophysical tests of the theory. In the present paper, we show the capabilities of X-ray reflection spectroscopy to constrain the inverse dimensionless fixed-point value $\gamma$ from the analysis of a Suzaku observation of the X-ray binary GRS 1915+105. We compare these constraints with those obtained from black hole imaging.

We construct a relativistic reflection model in a non-Kerr spacetime in which, depending on the value of the deformation parameter of the metric, there are black hole solutions with spin parameter $|a_*| > 1$. We apply our model to fit Suzaku data of four Seyfert galaxies (Ton S180, Ark 120, 1H0419-577, and Swift J0501.9-3239). These galaxies host at the center supermassive black holes that were previously interpreted as near-extremal Kerr black holes. For Ton S180 and 1H0419-577, our measurements are still consistent with the Kerr hypothesis. For Ark 120 and Swift J0501.9-3239, the Kerr solution is not recovered at 3-$\sigma$. We discuss our results and possible systematic uncertainties in the model.

Signatures of X-ray reprocessing (reflection) out of an accretion disk are commonly observed in the high-energy spectrum of accreting black holes, and can be used to probe the strong gravity region around these objects. In this paper, we extend previous work in the literature and we employ a full emission model for relativistic reflection in non-Kerr spacetime to demonstrate an approach that tests the Kerr black hole hypothesis. We analyze a composite spectrum obtained with the Proportional Counter Array in the Rossi X-ray Timing Explorer (RXTE), of the stellar-mass black hole GX 339-4 in its brightest hard state. With a remarkable sensitivity of ~0.1% and 40 million counts in the 3-45 keV band to capture the faint features in the reflection spectrum, we demonstrate that it is possible with existing data and an adequate model to place constraints on the black hole spin $a_*$ and the deformation parameter that quantifies the departure from the Kerr metric. Our measurement obtained with the best fit model, which should be regarded as principally a proof of concept, is $a_*=0.92^{+0.07}_{-0.12}$ and $\alpha_{13}=-0.76^{+0.78}_{-0.60}$ with a 90% confidence level and is consistent with the hypothesis that the compact object in GX 339-4 is a Kerr black hole. We also discuss how the physical model choice and the emissivity profile adopted could make an impact on the constraints of $\alpha_{13}$ and spin. To enable Kerr metric test using X-ray reflection spectroscopy, it is essential to improve our astrophysical understanding of accreting black holes, e.g., the natures of accretion flow and corona.

X-ray reflection spectroscopy is potentially a powerful tool to probe the spacetime geometry around astrophysical black holes and test general relativity in the strong field regime. However, precision tests of general relativity are only possible if we employ the correct astrophysical model and we can limit the systematic uncertainties. It is thus crucial to select the sources and the observations most suitable for these tests. In this work, we analyze simultaneous observations of XMM-Newton and NuSTAR of the supermassive black hole in Fairall 9. This source has a number of properties that make it a promising candidate for tests of general relativity using X-ray reflection spectroscopy. Nevertheless, we find that with the available data there is not a unique interpretation of the spectrum of Fairall 9, which prevents, for the moment, to use this source for robust tests of general relativity. This issue may be solved by future X-ray missions with a higher energy resolution near the iron line.

X-ray reflection spectroscopy is a powerful tool for probing the strong gravity region of black holes and can be used for testing general relativity in the strong field regime. Simplifications of the available relativistic reflection models limit the capability of performing accurate measurements of the properties of black holes. In this paper, we present an extension of the model RELXILL_NK in which the accretion disk has a finite thickness rather than being infinitesimally thin. We employ the accretion disk geometry proposed by Taylor & Reynolds (2018) and we construct relativistic reflection models for different values of the mass accretion rate of the black hole. We apply the new model to high quality Suzaku data of the X-ray binary GRS 1915+105 to explore the impact of the thickness of the disk on tests of the Kerr metric.

In a previous paper, one of us calculated iron K$\alpha$ line profiles emitted from possible accretion disks around traversable wormholes as a first step to use X-ray reflection spectroscopy to search for astrophysical wormholes in active galactic nuclei. In the present paper, we extend that work and construct an XSPEC model for the whole relativistic reflection spectrum. We apply our model to XMM-Newton and NuSTAR observations of the supermassive object in MCG-6-30-15 and we check whether these observations prefer the hypothesis that the central body is a Kerr black hole or a traversable wormhole. We find that the two models provide equally good fits, so with the available data we cannot distinguish the black hole and wormhole scenarios.

There are many simplifications in current relativistic reflection models and this introduces systematic uncertainties in the measurement of the properties of the source. In this paper, we study the impact of the radiation crossing the equatorial plane between the black hole and the inner edge of the accretion disk; that is, the radiation produced by the other side of the disk or circling the black hole one or more times (higher order disk images). For slow-rotating black holes with a larger plunging region, the effect is not very sensitive to the exact inclination angle of the disk. For fast-rotating black holes with a smaller plunging region, the effect is very weak for low inclination angles and becomes more important as the angle increases. We simulate some observations without and with higher order disk images and fit the data with the reflection model RELXILL_NK to check its capability of recovering the correct input parameters. Our results suggest that the effect of higher order disk images can be safely ignored for observations with present and near future X-ray missions, even for tests of the Kerr hypothesis.

Relativistic reflection features are commonly observed in the X-ray spectra of stellar-mass and supermassive black holes and originate from illumination of the inner part of the accretion disk by a hot corona. All the available relativistic reflection models assume that the disk is infinitesimally thin and the inner edge is at the innermost stable circular orbit or at a larger radius. However, we know that several sources, especially among supermassive black holes, have quite high mass accretion rates. In such a case, the accretion disk becomes geometrically thick and the inner edge of the disk is expected to be inside the innermost stable circular orbit. In this work, we employ the Polish donut model to describe geometrically thick disks and we study the iron line shapes from similar systems. We also simulate full reflection spectra and we analyze the simulated observations with a thin disk relativistic reflection model to determine the impact of the disk structure on the estimation of the model parameters, in particular in the case of tests of the Kerr hypothesis.

In a previous paper, we tried to test the Kerr nature of the stellar-mass black hole in GRS 1915+105 by analyzing NuSTAR data of 2012 with our reflection model RELXILL_NK. We found that the choice of the intensity profile of the reflection component is crucial and eventually we were not able to get any constraint on the spacetime metric around the black hole in GRS 1915+105. In the present paper, we study the same source with Suzaku data of 2007. We confirm that the intensity profile plays an important role, but now we find quite stringent constraints consistent with the Kerr hypothesis. The key differences with respect to our previous study are likely the lower disk temperature in the Suzaku observation and the higher energy resolution near the iron line of the Suzaku data. We also apply different RELXILL flavors (different descriptions of the coronal spectrum and variable disk electron density) obtaining essentially the same results. We thus conclude that this choice is not very important for our tests of the Kerr hypothesis while the intensity profile does play an important role, and that with high quality data it is possible to measure both the spacetime metric and the intensity profile.

In this paper, we briefly present RELXILL_NK, the first and currently only readily available model of the relativistic reflection spectrum of black hole accretion disks that includes non-Kerr solutions for the black hole spacetime, thus allowing for tests of the Kerr hypothesis and GR. RELXILL_NK makes use of a general relativistic ray-tracing code to calculate the relativistic effects of any well-behaved, stationary, axisymmetric, and asymptotically flat black hole spacetime, while the disk physics is handled through the non-relativistic X-ray reflection code XILLVER. A number of different flavors are available within RELXILL_NK; we summarize and compare these flavors using the Johannsen metric for the black hole spacetime.

Can the observation of the "shadow" allow us to distinguish a black hole from a more exotic compact object? We study the motion of photons in a class of vacuum static axially-symmetric space-times that is continuously linked to the Schwarzschild metric through the value of one parameter that can be interpreted as a measure of the deformation of the source. We investigate the lensing effect and shadow produced by the source with the aim of comparing the expected image with the shadow of a Schwarzschild black hole. In the context of astrophysical black holes we found that it may not be possible to distinguish an exotic source with small deformation parameter from a black hole. However, as the deformation increases noticeable effects arise. Therefore, the future more precise measurement of the shadow of astrophysical black hole candidates would in principle allow to put constraints on the deviation of the object from spherical symmetry.

Einstein's gravity has undergone extensive tests in the weak field gravitational limit, with results in agreement with theoretical predictions. There exist theories beyond general relativity (GR) which modify gravity in the strong field regime but agree with GR in the weak field. Astrophysical black holes are believed to be described by the Kerr metric and serve as suitable candidates to test strong gravity with electromagnetic radiation. We perform such a test by fitting one Suzaku dataset of the narrow-line Seyfert 1 (NLS1) galaxy Mrk 335 with X-ray reflection spectroscopy, using the Johannsen metric to model the black hole spacetime and test for deviations from Kerr. We find the data is best modeled with a hybrid model that includes both partial covering absorption and a reflection component. This is the first time such a model has been proposed for a high-flux (low reflection) Mrk 335 dataset. We constrain the Johannsen deformation parameter $\alpha_{13}$ to $-1.5<\alpha_{13}<0.6$ with spin parameter $a_{*}>0.8$, and the $\alpha_{22}$ parameter to $-0.4<\alpha_{22}<2.1$ with $a_{*}>0.7$, both at the 99% confidence level. Although additional solutions at large deviations from the Kerr metric show statistical similarity with the ones above, further analysis suggests these solutions may be manifestations of uncertainties beyond our control and do not represent the data. Hence, our results are in agreement with the idea that the supermassive compact object at the center of Mrk 335 is described by the Kerr metric.

We present the public release version of relxill_nk, an X-ray reflection model for testing the Kerr hypothesis and general relativity. This model extends the relxill model that assumes the black hole spacetime is described by the Kerr metric. We also present relxilllp_nk, the first non-Kerr X-ray reflection model with a lamppost corona configuration, as well as all other models available in the full relxill_nk package. In all models the relevant relativistic effects are calculated through a general relativistic ray-tracing code that can be applied to any well-behaved, stationary, axisymmetric, and asymptotically flat black hole spacetime. We show that the numerical error introduced by using a ray-tracing code is not significant as compared with the observational error present in current X-ray reflection spectrum observations. In addition, we present the reflection spectrum for the Johannsen metric as calculated by relxill_nk.

We continue exploring the constraining capabilities of X-ray reflection spectroscopy to test the Kerr-nature of astrophysical black holes and we present the results of our analysis of two NuSTAR observations of Cygnus X-1 in the soft state. We find that the final measurement can strongly depend on the assumption of the intensity profile. We conclude that Cygnus X-1 is not suitable for tests of general relativity using X-ray reflection spectroscopy and we discuss the desired properties of a source to be a good candidate for our studies.

In this paper, we review our current efforts to test General Relativity in the strong field regime by studying the reflection spectrum of supermassive black holes. So far we have analyzed 11 sources with observations of NuSTAR, Suzaku, Swift, and XMM-Newton. Our results are consistent with general relativity, according to which the spacetime metric around astrophysical black holes should be well approximated by the Kerr solution. We discuss the systematic uncertainties in our model and we present a preliminary study on the impact of some of them on the measurement of the spacetime metric.

We present NKBB, an XSPEC model for the thermal spectrum of thin accretion disks in parametric black hole spacetimes. We employ the Novikov-Thorne model for the description of the accretion disk and the formalism of the transfer function proposed by Cunningham for storing all the relativistic effects of the spacetime metric. The current version of the model assumes the Johannsen metric, but it can be easily extended to any stationary, axisymmetric, and asymptotically-flat spacetime without pathological properties. The model can be used within the XSPEC package to test the Kerr nature of astrophysical black holes with the continuum-fitting method.

Astrophysical black hole systems are the ideal laboratories for testing Einstein's theory of gravity in the strong field regime. We have recently developed a framework which uses the reflection spectrum of black hole systems to perform precision tests of general relativity by testing the Kerr black hole hypothesis. In this paper, we analyze XMM-Newton and NuSTAR observations of the supermassive black hole in the Seyfert 1 galaxy MCG-06-30-15 with our disk reflection model. We consider the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$, which quantify deviations from the Kerr metric. For $\alpha_{22} = 0$, we obtain the black hole spin $0.928 < a_* < 0.983$ and $-0.44 < \alpha_{13} < 0.15$. For $\alpha_{13} = 0$, we obtain $0.885 < a_* < 0.987$ and $-0.12 < \alpha_{22} < 1.05$. The Kerr solution is recovered for $\alpha_{13} = \alpha_{22} = 0$. Thus, our results include the Kerr solution within statistical uncertainties. Systematic uncertainties are difficult to account for, and we discuss some issues in this regard.

We employ the accretion disk reflection model RELXILL_NK to test the spacetime geometry around the stellar-mass black hole in GRS 1915+105. We adopt the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$, for which the Kerr solution is recovered when $\alpha_{13} = \alpha_{22} = 0$. We analyze a NuSTAR observation of 2012, obtaining vanishing and non-vanishing values of the deformation parameters depending on the astrophysical model adopted. Similar difficulties were not found in our previous tests with other sources. The results of this work can shed light on the choice of sources suitable for testing the Kerr metric using X-ray reflection spectroscopy and on the parts of our reflection models that more urgently require improvement.

We extend our previous work on tests of the Kerr black hole hypothesis using X-ray reflection spectroscopy and we infer observational constraints on the Johannsen deformation parameter $\epsilon_3$. We analyze a NuSTAR observation of the stellar-mass black hole in GS 1354-645, Suzaku data of the supermassive black holes in Ark 564, 1H0419-577, and PKS 0558-504, and some simultaneous XMM-Newton and NuSTAR observations of the supermassive black hole in MCG-6-30-15. All our measurements of $\epsilon_3$ are consistent with the Kerr metric and we find quite strong constraints from Ark 564 and MCG-6-30-15. We discuss the implications of our results and the next steps to improve our tests.

We present the study of a sample of seven "bare" active galactic nuclei (AGN) observed with Suzaku. We interpret the spectrum of these sources with a relativistic reflection component and we employ our model RELXILL_NK to test the Kerr nature of their supermassive black holes. We constrain the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$, in which the Kerr metric is recovered when $\alpha_{13} = \alpha_{22} = 0$. All our measurements are consistent with the hypothesis that the spacetime geometry around these supermassive objects is described by the Kerr solution. For some sources, we obtain quite strong constraints on $\alpha_{13}$ and $\alpha_{22}$ when compared to those found in our previous studies. We discuss the systematic uncertainties in our tests and the implications of our results.

X-ray reflection spectroscopy is a promising technique for testing general relativity in the strong field regime, as it can be used to test the Kerr black hole hypothesis. In this context, the parametrically deformed black hole metrics proposed by Konoplya, Rezzolla \& Zhidenko (Phys. Rev. D93, 064015, 2016) form an important class of non-Kerr black holes. We implement this class of black hole metrics in \textsc{relxill\_nk}, which is a framework we have developed for testing for non-Kerr black holes using X-ray reflection spectroscopy. We compare the results from this new model in the Kerr limit with existing codes, illustrate the effect of the leading deformation parameters on the X-ray reflection spectrum and present constraints on some KRZ parameters using data from the supermassive black hole in Ark 564.

It is thought that the spacetime metric around astrophysical black holes is well described by the Kerr solution of Einstein's gravity. However, a robust observational evidence of the Kerr nature of these objects is still lacking. Here we fit the X-ray spectrum of the stellar-mass black hole in GS 1354-645 with a disk reflection model beyond Einstein's gravity in order test the Kerr black hole hypothesis. We consider the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$. The Kerr metric is recovered for $\alpha_{13} = \alpha_{22} = 0$. For $\alpha_{22} = 0$, our measurements of the black hole spin and of the deformation parameter $\alpha_{13}$ are $a_* > 0.975$ and $-0.34 < \alpha_{13} < 0.16$, respectively. For $\alpha_{13} = 0$, we find $a_* > 0.975$ and $-0.09 < \alpha_{22} < 0.42$. All the reported uncertainties are at 99% of confidence level for two relevant parameters.

The observation of the shadows cast by the event horizon of black holes on the light emitted in its neighborhood is the target of current very-long-baseline-interferometric observations. When considering supermassive black holes, the light source is the black hole's accretion disk, and therefore, the observation of the shadow may reveal information about the black hole and the accretion flow. We here study the shadows cast by stellar-mass black holes that are illuminated not by an accretion disk but by a stellar companion in a wide binary orbit. We calculate the shadows produced in such a configuration for the first time and show snapshots of the time-dependent shadow "movie" that is generated. We also study the minimal criteria for detecting and resolving such shadows with very-long-baseline-interferometric observations. We find that one would need telescopes capable of resolving apparent magnitudes greater than $33$ with baselines larger than $10^{6}$ km. Therefore, although very-long-baseline-interferometric efforts such as the Event Horizon Telescope would not be able to detect these shadows, their detection is possible with future telescopes in the next few decades.

Recently, a rotating black hole solution was found in Vector-Tensor Galileons modified gravity that has some significant differences from the Kerr black hole solution. We study the iron line shape that is part of the reflection spectrum of accretion disks around black holes in this new black hole solution. We simulate and compare the iron lines of this solution with those of the Kerr solution to see if the technique of iron line spectroscopy can be used as a tool to test Vector-Tensor Galileons modified gravity. Our analysis shows that current X-ray facilities can, in principle, be used to test and place constraints on Vector-Tensor Galileons modified gravity in the strong coupling and ultraspinning regimes.

Einstein's theory of general relativity has been extensively tested in weak gravitational fields, mainly with experiments in the Solar System and observations of radio pulsars, and current data agree well with the theoretical predictions. Nevertheless, there are a number of scenarios beyond Einstein's gravity that have the same predictions for weak fields and present deviations only when gravity becomes strong. Here we try to test general relativity in the strong field regime. We fit the X-ray spectrum of the supermassive black hole in Ark 564 with a disk reflection model beyond Einstein's gravity, and we are able to constrain the black hole spin $a_*$ and the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$ separately. For $\alpha_{22} = 0$, we find $a_* > 0.96$ and $-1.0 < \alpha_{13} < 0.2$ with a 99% confidence level. For $\alpha_{13} = 0$, we get $a_* > 0.96$ and $-0.1 < \alpha_{22} < 0.9$ with a 99% confidence level. Our measurements are thus consistent with the hypothesis that the supermassive compact object in Ark 564 can be described by the Kerr metric.

Observations of the black hole shadow of supermassive black holes, such as Sagittarius A* at the center of our Milky Way galaxy, allow us to study the properties of black holes and the nature of strong-field gravity. According to the Kerr hypothesis, isolated, stationary, and axisymmetric astrophysical black holes are described by the Kerr metric. The Kerr hypothesis holds in General Relativity and in some modified gravity theories, but there are others in which it is violated. In principle, black hole shadow observations can be used to determine if the Kerr metric is the correct description for black holes, and in turn, they could be used to place constraints on modified gravity theories that do not admit the Kerr solution. We here investigate whether black hole shadow observations can constrain deviations from general relativity, focusing on two well-motivated modified quadratic gravity theories: Einstein-dilaton-Gauss-Bonnet gravity and dynamical Chern-Simons gravity. We show that current constraints on Einstein-dilaton-Gauss-Bonnet gravity are stronger than any that could be placed with black hole shadow observations of supermassive black holes. We also show that the same holds for dynamical Chern-Simons gravity through a systematic bias and a likelihood analysis when considering slowly-rotating supermassive black holes. However, observations of more rapidly-rotating black holes, with dimensionless spins $|\vec{J}|/M^{2}\simeq0.5$, could be used to better constrain dynamical Chern-Simons gravity.

Recently, two of us have found a family of singularity-free rotating black hole solutions in Einstein's conformal gravity. These spacetimes are characterized by three parameters: the black hole mass $M$, the black hole spin angular momentum $J$, and a parameter $L$ that is not specified by the theory but can be expected to be proportional to the black hole mass $M$. The Kerr black hole solution of Einstein's gravity is recovered for $L = 0$. In a previous paper, we showed that X-ray data of astrophysical black holes require $L/M < 1.2$. In the present paper, we report the results of a more sophisticated analysis. We apply the X-ray reflection model {\sc relxill\_nk} to \textsl{NuSTAR} and \textsl{Swift} data of the supermassive black hole in 1H0707-495. We find the constraint $L/M < 0.45$ (90% confidence level).

Einstein's theory of general relativity was proposed over 100 years ago and has successfully passed a large number of observational tests in the weak field regime. However, the strong field regime is largely unexplored, and there are many modified and alternative theories that have the same predictions as Einstein's gravity for weak fields and present deviations when gravity becomes strong. RELXILL_NK is the first relativistic reflection model for probing the spacetime metric in the vicinity of astrophysical black holes and testing Einstein's gravity in the strong field regime. Here we present our current constraints on possible deviations from Einstein's gravity obtained from the black holes in 1H0707-495, Ark 564, GX 339-4, and GS 1354-645.

Observations of the continuum spectrum emitted by accretion disks around black holes allows us to infer their properties, including possibly whether black holes are described by the Kerr metric. Some modified gravity theories do not admit the Kerr metric as a solution, and thus, continuum spectrum observations could be used to constrain these theories. We here investigate whether current and next generation X-Ray observations of the black hole continuum spectrum can constrain such deviations from Einstein's theory, focusing on two well-motivated modified quadratic gravity theories: dynamical Chern-Simons gravity and Einstein-dilaton-Gauss-Bonnet gravity. We do so by determining whether the non-Kerr deviations in the continuum spectrum introduced by these theories are larger than the observational error intrinsic to the observations. We find that dynamical Chern-Simons gravity cannot be constrained better than current bounds with current or next generation continuum spectrum observations. Einstein-dilaton-Gauss-Bonnet gravity, however, may be constrained better than current bounds with next generation telescopes, as long as the systematic error inherent in the accretion disk modeling is decreased below the predicted observational error.

Future electromagnetic observations of black holes may allow us to test General Relativity in the strong-field regime. Such tests, however, require knowledge of rotating black hole solutions in modified gravity theories, a class of which does not admit the Kerr metric as a solution. Several rotating black hole solutions in modified theories have only been found in the slow-rotation approximation (i.e. assuming the spin angular momentum is much smaller than the mass squared). We here investigate whether the systematic error due to the approximate nature of these black hole metrics is small enough relative to the observational error to allow their use in electromagnetic observations to constrain deviations from General Relativity. We address this by considering whether electromagnetic observables constructed from a slow-rotation approximation to the Kerr metric can fit observables constructed from the full Kerr metric with systematic errors smaller than current observational errors. We focus on black hole shadow and continuum spectrum observations, as these are the least influenced by accretion disk physics, with current observational errors of about 10%. We find that the fractional systematic error introduced by using a second-order, slowly-rotating Kerr metric is at most 2% for shadows created by black holes with dimensionless spins $\chi\leq0.6$. We also find that the systematic error introduced by using the slowly-rotating Kerr metric as an exact metric when constructing continuum spectrum observables is negligible for black holes with dimensionless spins of $\chi \lesssim 0.9$. Our results suggest that the modified gravity solutions found in the slow-rotation approximation may be used to constrain realistic deviations from General Relativity with continuum spectrum and black hole shadow observations.