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.

The reflection spectroscopic model RELXILL is commonly implemented in studying relativistic X-ray reflection from accretion disks around black holes. We present a systematic study of the model's capability to constrain the dimensionless spin and ionization parameters from $\sim$6,000 NuSTAR simulations of a bright X-ray source employing the lamppost geometry. We employ high count spectra to show the limitations in the model without being confused with limitations in signal-to-noise. We find that both parameters are well-recovered at 90% confidence with improving constraints at higher reflection fraction, high spin, and low source height. We test spectra across a broad range - first at 10$^6-$10$^7$ and then $\sim$10$^5$ total source counts across the effective 3-79 keV band of NuSTAR, and discover a strong dependence of the results on how fits are performed around the starting parameters, owing to the complexity of the model itself. A blind fit chosen over an approach that carries some estimates of the actual parameter values can lead to significantly worse recovery of model parameters. We further stress on the importance to span the space of nonlinear-behaving parameters like $log~\xi$ carefully and thoroughly for the model to avoid misleading results. In light of selecting fitting procedures, we recall the necessity to pay attention to the choice of data binning and fit statistics used to test the goodness of fit by demonstrating the effect on the photon index $\Gamma$. We re-emphasize and implore the need to account for the detector resolution while binning X-ray data and using Poisson fit statistics instead while analyzing Poissonian data.