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.

We analyze eleven NuSTAR and Swift observations of the black hole X-ray binary GX 339-4 in the hard state, six of which were taken during the end of the 2015 outburst, five during a failed outburst in 2013. These observations cover luminosities from 0.5%-5% of the Eddington luminosity. Implementing the most recent version of the reflection model relxillCp, we perform simultaneous spectral fits on both datasets to track the evolution of the properties in the accretion disk including the inner edge radius, the ionization, and temperature of the thermal emission. We also constrain the photon index and electron temperature of the primary source (the "corona"). We find the disk becomes more truncated when the luminosity decreases, and observe a maximum truncation radius of $37R_g$. We also explore a self-consistent model under the framework of coronal Comptonization, and find consistent results regarding the disk truncation in the 2015 data, providing a more physical preferred fit for the 2013 observations.