---
abstract: 'X-ray reflection models are used to constrain the properties of the accretion disk, such as the degree of ionization of the gas and the elemental abundances. In combination with general relativistic ray tracing codes, additional parameters like the spin of the black hole and the inclination to the system can be determined. However, current reflection models used for such studies only provide angle-averaged solutions for the flux reflected at the surface of the disk. Moreover, the emission angle of the photons changes over the disk due to relativistic light bending. To overcome this simplification, we have constructed an angle-dependent reflection model with the [xillver]{} code and self-consistently connected it with the relativistic blurring code [relline]{}. The new model, [relxill]{}, calculates the proper emission angle of the radiation at each point on￼ the accretion disk, and then takes the corresponding reflection spectrum into account. We show that the reflected spectra from illuminated disks follow a limb-brightening law highly dependent on the ionization of disk and yet different from the commonly assumed form $I \propto \ln(1+1/\mu)$. A detailed comparison with the angle-averaged model is carried out in order to determine the bias in the parameters obtained by fitting a typical relativistic reflection spectrum. These simulations reveal that although the spin and inclination are mildly affected, the Fe abundance can be over-estimated by up to a factor of two when derived from angle-averaged models. The fit of the new model to the [*Suzaku*]{} observation of the Seyfert galaxy Ark 120 clearly shows a significant improvement in the constrain of the physical parameters, in particular by enhancing the accuracy in the inclination angle and the spin determinations.'
author:
- 'J. García, T. Dauser, A. Lohfink, T. R. Kallman, J. Steiner, J. E. McClintock, L. Brenneman J. Wilms, W. Eikmann, C. S. Reynolds, F. Tombesi'
bibliography:
- 'my-references.bib'
title: 'Improved Reflection Models of Black-Hole Accretion Disks: Treating the Angular Distribution of X-rays'
---

Introduction
============

Despite the huge difference in mass, Galactic black holes (GBH) in binary systems and super-massive black holes in active galactic nuclei (AGN) show similar X-ray properties. In particular, in many sources an emission component is present, which is generated by high-energy coronal photons that are reprocessed in an optically-thick accretion disk. The most prominent feature in this component – commonly referred to as the “reflection” spectrum – is the Fe K emission line at 6–7 keV, which is produced by fluorescence. This reflection spectrum from the inner disk region carries information on the physical composition and condition of the matter in strong gravitational fields [e.g., @fab00; @fab03; @rey03; @dov04; @mil08]. Doppler effects, light bending and gravitational redshift skew the Fe line (and the fluorescence lines of other elements) and can, for rapidly spinning black holes, extend the red wing of the line to very low energies [@fab82; @lao91; @dab97; @dov04; @bre06; @dau10]. Hence, in order to fit observations, the pure reflection spectrum has to be convolved by relativistic blurring algorithms that account for these relativistic effects [e.g., @mil08; @ste11; @rey12; @fab12b; @Dauser2012a]. This relativistic smearing highly depends on the parameters of the system, like the black hole spin and the inclination of the system. However, most current reflection models used for such studies only provide an angle-averaged solution for the flux reflected at the surface of the disk, which can systematically affect the inferred disk parameters. Furthermore, an isotropic distribution or a particular limb-darkening law is often chosen [e.g. @chi12; @wal13], despite the fact that radiative transfer calculations predict an emission excess near grazing angles [known as limb brightening, e.g., @roz08; @svo09; @roz11]. In fact, the pioneer and widely used relativistic line convolution model [laor]{} [@lao91] intrinsically assumes a limb-darkening law.

Several authors have studied relativistic effects on the emission of X-rays from accretion disks, while considering the angular distribution of the radiation. [@mar00] performed calculations of the relativistic effects acting on both the reflection continuum and the iron line in a Kerr-metric, and they showed that for rapidly spinning black holes the line equivalent width is substantially enhanced. [@rey00] discussed (within the context the analysis of [*ASCA*]{} data for NGC 4258) how gravitational light bending can affect the Fe line emission in high-inclination sources. [@bec04] and [@dov04] presented similar calculations emphasizing how the assumed angular emissivity law (limb darkening or brightening) can significantly affect the line profile when relativistic effects are important. These studies were extended by [@nie08] in order to explain the Fe line profile and variability observed in MCG$-$6-30-15. More recently, [@svo09] performed relativistic radiative transfer calculations of X-ray irradiated disk atmospheres with the [noar]{} code to determine the impact of light bending on X-rays in the 2–10 keV band. They concluded that the uncertainty in the angular distribution of reflected radiation can translate into a $\sim 20\%$ uncertainty in the determination of $R_\mathrm{in}$, and thus also in the spin parameter $a_*$. The emphasis of their work is on the relativistic effects and the estimation of uncertainties, and they omit a detailed ionization-balance calculation of the reflection spectrum.

In the past, X-ray reflection from a cold, neutral slab for the K lines of heavy elements has been discussed extensively in the context of X-ray binaries [@bas78] and AGN [@geo91; @mat91]. Additional earlier work treats the angular-dependence of electron scattering in cold material [@ghi94], and in hot thermal plasmas [@haa93], while neglecting photoelectric absorption and line emission. [@mat96] carried out detailed Monte Carlo calculations of the Fe K$\alpha$ emission from X-ray photoionized accretion disks, showing that the line strength and shape depend on both the ionization stage of the material and on the disk inclination angle. [@vrt93] computed the disk corona structure in X-ray binaries under the assumptions of ionization, thermal, and hydrostatic balance of gas illuminated by the central continuum source. They explained the Fe K emission and the absorption edge in terms of the inclination angle and the shape of the incident source spectrum. More sophisticated ionized reflection calculations have been presented by [@roz08] in the context of accretion disks around supermassive black holes; their models depict the effects of a limb-brightening law on the reflected X-rays. The reflection model that has been most widely used by observers, for both general application and for measuring black hole spin via the Fe-line method, is [reflionx]{} [@ros05]. This industry-standard model solves the radiation transfer problem using a diffusion equation, which imposes the limitation that the model can only deliver an angle-averaged spectrum.

We overcome the simplifications adopted by previous models – i.e., limited ionization balance calculations and angle-averaged reflected flux – by exploiting the full capabilities of our reflection code [xillver]{} [@gar10; @gar13a], to treat the angular distribution of the reflected X-rays. The Feautrier method we employ [@fea64] enables [xillver]{} to calculate the specific intensity of the radiation field as a function of energy, position in the slab, and viewing angle. This allows us to construct a grid of reflection models in which the inclination angle is included as an explicit fitting parameter. Furthermore, using this approach we are able to address a further complication. Due to general relativistic light-bending, photons are emitted over a wide range of angles depending on their location on the accretion disk. This complication can be solved by directly linking the reflection code to a relativistic smearing kernel. Accordingly, we use the angle-dependent [xillver]{} reflection code to assign a proper reflected spectrum for each point on the disk, which then depends on the actual emission angle derived from general relativity. These single spectra are then relativistically smeared using [relline]{} [@dau10; @dau13] before being integrated to yield the total reflection spectrum. As each of these emission points has to be treated separately, a simple combination of [xillver]{} and [relline]{} will yield inconsistent results.

In this paper we present for the first time a complete description of the spectrum reflected by an ionized accretion disk around a black hole that is complete, self-consistent, and takes into account the angular distribution of the reflected X-rays. The new model, [relxill]{}, is provided in the standard format, enabling its use with commonly employed fitting packages such as [xspec]{} [@arn96] and [isis]{} [@hou00]; the model’s input parameters are the same as that of its parent model [relconv]{}. An additional version of the new model, [relxill\_lp]{}, which simulates the reflected spectra assuming a lamppost geometry, is also supplied.

This paper is organized as follows: In Section \[sec:xill\], we describe our angle-dependent solution for the reflected spectrum and compare it with the commonly used angle-averaged solution. Sections \[sec:acc\] and \[sec:ang\] first explore how relativistic effects modify the reflected spectrum, and they then discuss how we integrated the reflection code [xillver]{} with the relativistic blurring code [relline]{}. In Section \[sec:sim\], we assess, via a detailed error analysis, how the angle-averaged solution yields biased values of the various model parameters. An application of our model to an analysis of [*Suzaku*]{} data for Ark 120 is discussed in Section \[sec:obs\], and we offer our conclusions in Section \[sec:con\].

Angular Solution of the Reflected Spectrum {#sec:xill}
==========================================

The X-ray reflection calculations presented in this paper are based upon those from [@gar13a], with the difference being that we now exploit the full capabilities of our reflection code [xillver]{} by extracting the solution of the emergent spectra for individual viewing angles. These models and the details of the numerical code have been extensively described by [@gar10], [@gar11], and [@gar13a]. Here we will discuss the theoretical aspects concerning the angular dependence of the radiation fields, and their effects on the observed spectrum.

The basic layout of the reflection problem is shown in the diagram of Figure \[cartoon\]. The gas in the atmosphere of the accretion disk (i.e., within a few Thomson depths) is assumed to be at a density of $n_{\mathrm H}=10^{15}$ cm$^{-3}$, in a plane-parallel geometry, where the coordinate of importance is shown in the vertical direction perpendicular to the slab ($z$-axis). The surface of the disk is set at $\tau_{\mathrm T}=0$, and the base of the atmosphere at $\tau_\mathrm{max}=10$, where $d\tau_{\mathrm T}=1.2n_{\mathrm H} \sigma_{\mathrm T}dz$ is the Thomson optical depth, and $\sigma_{\mathrm T}=6.65\times 10^{-25}$ cm$^2$ is the Thomson classical cross section for electron scattering.

The disk is illuminated by a primary source of X-ray radiation from above. The radiation enters the slab at $\tau_{\mathrm T}=0$, and is reprocessed by the atmosphere. The intensity of the radiation field $I(z,\mu,{\varepsilon})$ at each position $z$ the atmosphere, per energy ${\varepsilon}$, unit area, and steradian is determined by solving the time-independent, one-dimensional plane-parallel radiation transfer equation $$\label{ertz}
\mu \frac{\partial I_{{\varepsilon}}(z,\mu)}{\partial z} = \eta_{{\varepsilon}}(z) - \chi_{{\varepsilon}}(z) I_{{\varepsilon}}(z,\mu)$$ where $\mu$ is the cosine of the angle $\theta$ with respect to the disk normal, and $\eta_{{\varepsilon}}(z)$ and $\chi_{{\varepsilon}}(z)$ are the total emissivity and opacity, respectively.

At each point in the slab one can distinguish between the incoming and the outgoing components of the radiation for $\pm \mu$ (see Figure \[cartoon\]), and thus write two separate transfer equations $$\label{ertp}
\mu\frac{\partial I^+_{{\varepsilon}}(\tau_{{\varepsilon}},\mu)}{\partial \tau_{{\varepsilon}}} = I^+_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) - S_{{\varepsilon}}(\tau_{{\varepsilon}})$$ $$\label{ertm}
- \mu\frac{\partial I^-_{{\varepsilon}}(\tau_{{\varepsilon}},\mu)}{\partial \tau_{{\varepsilon}}} = I^-_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) - S_{{\varepsilon}}(\tau_{{\varepsilon}})$$ where $I^+_{{\varepsilon}}$ and $I^-_{{\varepsilon}}$ refer to positive and negative values of $\mu$, respectively, and $\tau_{{\varepsilon}} = -\chi_{{\varepsilon}}(z)dz$ is the total optical depth. By defining the symmetric and anti-symmetric averages $$\label{eu}
u_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) = \frac{1}{2}\left[ I^+_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) + I^-_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) \right]$$ $$\label{ev}
v_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) = \frac{1}{2}\left[ I^+_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) - I^-_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) \right]$$ the transfer equation (\[ertz\]) can be rewritten as $$\label{ertu}
\mu^2 \frac{\partial^2 u_{{\varepsilon}}(\tau_{{\varepsilon}},\mu)}{\partial \tau_{{\varepsilon}}^2} = u_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) - S_{{\varepsilon}}(\tau_{{\varepsilon}}),$$ where $S_{{\varepsilon}}(\tau_{{\varepsilon}}) = \eta_{{\varepsilon}}/\chi_{{\varepsilon}}$ is the source function. The line plus continuum emissivities and opacities depend on the structure of the gas, which is determined at each point in the slab by the photoionization code [xstar]{} [@kal01]. Details on the calculation of these quantities and on the atomic data implemented in these calculations are fully described in [@gar10] and [@gar13a]. Expression (\[ertu\]) is a second-order differential equation subject to two boundary conditions. At the top ($\tau_{{\varepsilon}}=0$), the incoming radiation field $I^-_{{\varepsilon}}(0,\mu)$ is known. Since $u_{{\varepsilon}} - v_{{\varepsilon}} = I^-_{{\varepsilon}}$, and $v_{{\varepsilon}} = \mu \frac{\partial u_{{\varepsilon}}}{\partial \tau_{{\varepsilon}}}$, the boundary condition at the surface of the slab can be expressed as $$\label{ebc1}
u_{{\varepsilon}}(0,\mu) - \mu \left( \frac{\partial u_{{\varepsilon}}}{\partial \tau_{{\varepsilon}}}\right)_{0} = I_\mathrm{inc}.$$ Assuming no irradiation from below, we set $I^+_{{\varepsilon}}(\tau_\mathrm{max},\mu)=0$ at the bottom of the atmosphere ($\tau_{{\varepsilon}}=\tau_\mathrm{max}$). Thus, using $u_{{\varepsilon}} + v_{{\varepsilon}} = I^+_{{\varepsilon}}$, the boundary condition at the bottom of the slab becomes $$\label{ebc2}
u_{{\varepsilon}}(\tau_\mathrm{max},\mu) + \mu \left( \frac{\partial u_{{\varepsilon}}}{\partial \tau_{{\varepsilon}}}\right)_{\tau_\mathrm{max}} = 0$$ and $\tau_\mathrm{max} = 10 \tau_\mathrm{T}$ is chosen for all the calculations.

The transfer equation (\[ertu\]) and its boundary conditions (\[ebc1\]) and (\[ebc2\]) are converted into a set of difference equations via the discretization of all variables, in the form $$\label{erti}
\mu^2 u_i'' = u_i - S_i$$ where $i=1,\ldots,N$ denotes a particular position in a grid of $N$ points $\tau_i$. Notice that, for simplicity, we have intentionally suppressed the explicit dependence on ${\varepsilon}$ and $\mu$.

The derivatives of $u_i$ in a non-uniform mesh can be found using the Taylor’s expansions $$\label{efor}
u_{i+1} = u_i + u_i' \Delta\tau_i + u_i'' \frac{\Delta\tau_i^2}{2} + ...$$ $$\label{ebac}
u_{i-1} = u_i - u_i' \Delta\tau_{i-1} + u_i'' \frac{\Delta\tau_{i-1}^2}{2} + ...$$ which are usually referred to as the forward and backward finite differences. Solving for $u_i'$ in (\[efor\]) and substituting in (\[ebac\]), the second derivative can be expressed as $$\label{esed}
u_i'' = \frac{2}{\Delta\tau_i + \Delta\tau_{i-1}} \left[
u_{i-1} \left(\frac{1}{\Delta\tau_{i-1}} \right)
- u_i \left(\frac{\Delta\tau_i + \Delta\tau_{i-1}}{\Delta\tau_i\Delta\tau_{i-1}}  \right)
u_{i+1} \left(\frac{1}{\Delta\tau_i}\right) 
\right].$$

This can now be inserted into (\[erti\]) to write $$\begin{aligned}
u_{i-1}\left[\frac{-2\mu^2}{\Delta\tau_{i-1}(\Delta\tau_i + \Delta\tau_{i-1})} \right]
 &+& u_{i}\left[\frac{2\mu^2}{\Delta\tau_{i}\Delta\tau_{i-1}} \right] \\
 &+& u_{i+1}\left[\frac{-2\mu^2}{\Delta\tau_{i}(\Delta\tau_i + \Delta\tau_{i-1})} \right]
 = S_i \nonumber\end{aligned}$$ Equation (\[esed\]) can be inserted into (\[efor\]) to get a similar expression for the first derivative in terms of $u_{i-1}, u_i,$ and $u_{i+1}$. However, for the boundary condition at the top $$\label{ebc}
u_1 - \mu u_1' = I_\mathrm{inc}$$ a different expression is needed. An alternative is to rely upon the second derivative provided by equation (\[erti\]). Thus $$u_i' = u_{i+1}\left(\frac{1}{\Delta\tau_i}\right) 
- u_i \left(\frac{1}{\Delta\tau_i} + \frac{\Delta\tau_i}{2\mu^2} \right)
       + S_i \frac{\Delta\tau_i}{2\mu^2}$$ and (\[ebc\]) becomes $$u_1 \left( 1 + \frac{\mu}{\Delta\tau_1} + \frac{\Delta\tau_1}{2\mu} \right)
- u_2 \left(\frac{\mu}{\Delta\tau_1}\right)
= I_\mathrm{inc} + S_1 \frac{\Delta\tau_1}{2\mu}$$ A similar expression can be found for the inner boundary at $\tau=\tau_\mathrm{max}$.

The numerical solution in this approach proceeds by a forward elimination plus backward substitution scheme. We start using the boundary condition at the surface to express $u_1$ in terms of $u_2$. The result is then applied to variables $u_1$, $u_2$ and $u_3$ in order to express $u_2$ in terms of $u_3$. Iteratively, and in the same way, we can express each $u_i$ in terms of $u_{i+1}$. We finally reach the last point in the grid, and the boundary condition for $\tau_\mathrm{max}$ gives $u_N$. Back substitutions then produce $u_N \rightarrow u_{N-1} \rightarrow u_{N-2} \rightarrow \cdots \rightarrow u_2
\rightarrow u_1$, so all $u_i$ are known.

Once the solution for all $u_i$ is known, it is trivial to obtain the outgoing intensity $I^+(0)$ at the surface of the slab (the reflected component), by simply using (\[eu\]) $$I^+_1 = 2u_1 - I^-_1$$ where, again, $I^-_1 = I_\mathrm{inc}$. This is the solution of the reflected intensity at the illuminated surface of the slab for every energy $E$ and every angle $\mu$. For the calculations presented here, the grid in $\mu$ angles is defined as $$\mu_j = \frac{j-0.5}{M} \ \ \ \ \ \ \ \ \ (j = 1, 2,\ldots, M)$$ with $M=10$, and $\Delta\mu = 1/M$. Note that for the [xillver]{} models previously published [@gar10; @gar11; @gar13a], the reflected spectra were provided in terms of the [*angle-averaged*]{} quantity, $$\label{eave}
F_{{\varepsilon}}^+(0) = \frac{1}{2} \int_0^1 I^+_{{\varepsilon}}(0,\mu) \mu d\mu$$ defined as the first moment of the radiation field (although we are only taking the component of the field for positives $\mu$, the one that reaches the observer).

The left column in Figure \[xillver\] shows the reflection spectra across all viewing angles for different values of the ionization parameter, as indicated in each panel. These are displayed in different colors (see right column for an identification). At first glance the spectra look very similar to each other. However, the ratio between spectra for individual $\mu$ and the angle-averaged spectrum reveals differences that can reach up to a factor 2.5. Especially interesting is that the spectral features behave differently than the continuum. This makes sense because the viewing angle affects the intensity of the emergent radiation by changing the [*effective optical depth*]{}, given by the geometrical projection of the optical depth along the line of sight according to $$\label{eteff}
\tau_\mathrm{eff} = \tau/\mu \quad.$$ An observer looking at grazing angles (small $\mu$) will effectively see a larger optical depth. Or in other words, photons emitted at a particular depth in the slab will see a larger effective depth to escape out of the disk at grazing angles. Therefore, looking at the energy distribution of the emergent radiation, spectral features that are originally produced at large $\tau$ (such as those from Fe and Ni) will be more attenuated than those that are produced near the surface (such as C, N, and O lines). This is also the case for photons in the continuum, which is the reason for the Compton hump above $20$ keV being strongly affected by the viewing angle [@lig80; @lig81; @lig88; @mag95]. In fact, high energy photons are the most affected, since they experience a larger number of scatterings (thus, they will be very sensitive to any small change in the optical depth).

The angular dependence of the reflected spectra for energies above $\sim
20$ keV is the same for all ionization parameters, since Compton scattering is the main source of opacity, and these energies are insensitive to the gas temperature. Conversely, the impact of the viewing angle at lower energies is more pronounced where the photoelectric opacity becomes important. Because low-ionization models are dominated by photoelectric absorption, angular effects are accordingly stronger the lower the ionization parameter. An important effect is the change in the slope of the reflected spectrum, as the the angular dependence of the solution is different at different energies. In general, models with an inclination of $60{^{\circ}}$-$70{^{\circ}}$ seem to be the closest to the angle-averaged solution. Therefore, we expect the models presented here to be most useful for systems with either very low or very high inclination angles.

The energy dependence of the angular effects on the reflected spectrum becomes very relevant when combined with relativistic blurring models, such as [relline]{}. Commonly, in these models one can choose between a limb-brightening or a limb-darkening law which behaves uniformly for all energies [@svo09; @dau10; @dau13]. In the right column of Figure \[xillver\], we compare the intensity predicted by our angular dependent simulation with these empirical models. It is clear that our simulations tend towards limb brightening, i.e., the integrated flux increases when the inclination angle increases. For low ionization models, the predicted trend departs significantly from the usual limb-brightening law, showing a much flatter profile. As the ionization grows, the integrated flux predicted by [xillver]{} becomes steeper and approaches the limb-brightening law. However, for log $\xi \gtrsim
3$, the profile becomes flatter, resembling isotropic behavior. This is due to the fact that for high-ionization models, the ionization structure is dominated by electron scattering, rather than photoionization and recombination. In this regime, photons suffer a large number of scatterings in a hot atmosphere before escaping towards the observer. The gas remains close to the Compton temperature for a large range of optical depths [@gar13a], meaning that an observer will see a nearly isothermal atmosphere at any inclination. These results are of great relevance when contrasted with earlier relativistic convolution models (e.g., [kdblur, kerrconv]{}), which arbitrarily assume a limb-darkening law [@lao91]. Instead, all our models predict a limb-brightening behavior highly dependent on the particular ionization state of the gas.

We have also explored the possibility of including an angular dependence in the scattering term of the source function of Equation (\[erti\]). In the absence of absorption and intrinsic emission (i.e., pure scattering), the source function is simply reduced to $S_{{\varepsilon}}(\tau_{{\varepsilon}})=J_{{\varepsilon}}(\tau_{{\varepsilon}})$, where $J_{{\varepsilon}}(\tau_{{\varepsilon}})=\int u_{{\varepsilon}}(\mu,\tau_{{\varepsilon}})d\mu$ is the zeroth moment of the radiation field (i.e., the mean intensity), which is an angle-averaged quantity. One can, however, introduce a phase function to describe the angular properties of the electron scattering. A simple approximation in the non-relativistic limit is the Rayleigh scattering phase function, for which the source function can be expressed as $$S_{{\varepsilon}}(\tau_{{\varepsilon}},\mu) = \frac{3}{8} [(3-\mu^2) J_{{\varepsilon}}(\tau_{{\varepsilon}}) + (3\mu^2-1) K_{{\varepsilon}}(\tau_{{\varepsilon}})]$$ [see Eq. 103, Chap. 1 of @cha60], where $K_{{\varepsilon}}(\tau_{{\varepsilon}}) = \int
u_{{\varepsilon}}(\mu,\tau_{{\varepsilon}})\mu^2d\mu$ is the second moment of the radiation field. Using this form of the source function we have calculated the radiation transfer on an isothermal, pure scattering atmosphere. Figure \[scattering\] shows that the angular distribution of the integrated flux behaves very closely to the isotropic case. This result demonstrates that the treatment of the Compton scattering in [xillver]{} is adequate for these kind of calculations and it is upheld by 3-dimensional Monte Carlo simulations which have shown that the reflection intensity of a Compton scattering dominated atmosphere appears nearly isotropic for all viewing angles.

Accretion Disks in General Relativity {#sec:acc}
=====================================

In this section we discuss the effects of strong gravity on photons emitted from the surface of an accretion disk around compact objects, and the numerical methods employed to predict the emitted spectrum perceived by a distant observer. In the following we briefly summarize how we calculate the general relativistic radiative transfer around a rotating black hole. Accordingly, the equations of motion in the @Kerr1963 metric are used to describe the photon and particle trajectories [@Bardeen1972]. The actual transfer is solved by using “Cunningham Transfer Function” [@Cunningham1975], as implemented by @Speith1995. This approach is used by the [relline]{} code, which was employed to obtain the results presented in the following.[^1]

In the following, we use the dimensionless spin parameter $a_*$, which takes values from $a_*=1$ (maximally rotating) to $a_*=0$ (non-spinning) and $a_*=-1$ (negatively rotating, i.e., the black hole and the accretion disk are counter rotating). Moreover, all lengths are given in units of the gravitational radius, which is defined as ${\ensuremath{r_\mathrm{g}}}= GM/c^2$.

Using the approach outlined above, we are able to calculate the apparent image of the accretion disk on the sky a distant observer would see [e.g. @Cunningham1973]. Figure \[disk1\] shows how light-bending creates apparent asymmetries in the shape and warps the disk towards the observer. Moreover, the energy shift the photons experience from the disk to the observer is illustrated as it varies over the disk [see @dau10 for more details and explicit formulae]. The main effects which impact the energy are the Doppler shift (which leads to the blue and red division of the disk) and the gravitational redshift, which increases with proximity to the black hole.

Besides the energy shift, general relativistic effects also alter the direction of the photon, i.e., the photon is likely to be emitted at a different angle than it is observed. From the momentum of the emitted photon $\vec{p}_\mathrm{e}$ [see @Bardeen1972], the emission angle can be easily calculated via $$\label{eq:emis_angle}
  \cos({\ensuremath{\theta_\mathrm{e}}}) =  \frac{\vec{p}_{\mathrm{e}\perp}}{|\vec{p}_\mathrm{e}|} \quad.$$ Figure \[disk2\] shows a map of ${\ensuremath{\theta_\mathrm{e}}}$ on the surface of a disk seen at inclinations of $\theta =40{^{\circ}}$ and $\theta=80{^{\circ}}$. In both cases, the emission angle converges as expected towards the viewing angle $\theta$ at sufficiently large distances from the central black hole. This is expected, as this implies that the photon travels on a straight trajectory to the observer. In particular for $\theta=40{^{\circ}}$, slight differences at larger distances are visible between the left and the right sides, due to the rotation of the accretion disk. Looking at the zoomed images in Figure \[disk2\], it is obvious that essentially every emission angle from $0{^{\circ}}$ and $90{^{\circ}}$ provides some contribution to the observed image. Specifically, ${\ensuremath{\theta_\mathrm{e}}}$ is significantly lower than the viewing angle for photons emitted behind the black hole, as light bending effects are stronger in this case. Even for $\theta=80{^{\circ}}$, the emission angle approaches $0{^{\circ}}$ for a small part of the disk. Again, due to the rotation of the disk and also the rotation of the black hole, this effect is stronger for particles moving towards the observer and in the rotational direction of the black hole. Therefore, the region of low emission angles is in both cases shifted asymmetrically to the left. Additionally, there is also a region of relatively high emission angles at the innermost radii of the disk. This is best seen for $\theta = 40{^{\circ}}$, visible as a blue ring in Figure \[disk2\]. Note that although most of the emission angles in the figure are close to the viewing angle, the inner parts radiate with much stronger intensity and generally dominate the observation. In detail, for a standard @Shakura1973a disk the emitted intensity behaves like $I \propto
r^{-3}$, while for a disk irradiated by a source on the black hole’s rotational axis, generally, an even steeper profile is expected [e.g., @Fukumura2007a; @Wilkins2012a; @dau13].

Angle-Dependent Blurred Reflection {#sec:ang}
==================================

As discussed in Section \[sec:xill\], the reflection spectra change significantly depending on the inclination angle. Moreover, in the previous section we showed that the emission angles are not necessarily equal (or even close) to the inclination angle of the system. Therefore, it is imperative to take into account the complete relativistic transfer function for each region on the accretion disk, in order to properly predict the emitted spectrum.

General Properties (`relxill`)
------------------------------

For this purpose, we have extended the relativistic blurring code [relline]{} and incorporated angle-dependent reflection from the [xillver]{} code. Accordingly, the model behaves similarly to a convolution of the [xillver]{} tables with `relconv`. But instead of smearing the angle averaged reflection, this new model, called `relxill`, combines the spectra from all points in the image plane according to their observed emission angles, all while properly smearing each point relativistically. Although more physically rigorous, under this new approach, no additional parameters enter the model. In fact, because the angular solution of the reflected spectrum is self-consistently calculated, no assumptions have to be made for the limb-brightening/darkening laws, and thus the [limb]{} parameter is removed from the model. However, we have included a flag variable called [angleon]{} to switch between the new model with the full angular solution ([angleon=1]{}), and the old angle-average version ([angleon=0]{}). This model is provided in the appropriate format to be used in combination with other models included in the commonly-used X-ray fitting packages such as [xspec]{} [@arn96] and [isis]{} [@hou00].

It is important to clarify that a simple convolution of the angle-dependent [xillver]{} spectra with a relativistic blurring code will yield incorrect results. As described in Section \[sec:acc\], in presence of relativistic effects photons emitted at many different angles can contribute to the spectrum observed at one particular inclination, due to strong light-bending. This effect is different at different locations on the surface of the accretion disk. Therefore, the convolution of the reflected spectra needs to be done differently at each particular radius, taking the correct contribution of the multiple angular solutions to the given viewing angle. Figure \[fig:comparison\] illustrates the differences in the final spectrum obtained from an simple convolution of [relconv]{} on the angle-dependent [xillver]{} spectra with the inclination angles linked together, and the one calculated self-consistently with our new model [relxill]{}.

In Figure \[fig:relconv\_xill\] we show the results of the new model `relxill` with typical parameter settings for different inclinations of the system. As is well known and evident in the spectra in the left panel, different inclinations primarily result in a net energy shift for the spectrum. However, when compared to the usual, averaged relativistic convolution (middle panel) it is obvious that a deviation is present which is highly dependent on the inclination. In general, these deviations are as large as 20% and distinct for lines and for the continuum. The reason for such differences becomes evident in the right panel, which shows the distribution of the flux per emission angle. In this particular case, there is a flat distribution at large inclinations ($\theta=80{^{\circ}}$), i.e., photons from almost all angles are being observed with a similar contribution. Consequently, the resulting spectra are also very similar. This is equivalent to an angle-averaged solution, which is assumed by common reflection codes. In contrast, for small inclination angles the distribution is peaked at the actual inclination angle. This leads to an intrinsically different reflection spectrum and hence the relativistically smeared spectrum also differs significantly. However, these results depend upon various parameter settings, like the ionization of the disk ($\xi$), the abundances, the photon index $\Gamma$, the emissivity $\epsilon$, or the spin $a_*$.

Note that while for low spin the relativistic effects are generally less pronounced, this does not mean that the difference between the average and the proper angle treatment vanishes. As it can be seen in Figure \[xillver\], even in the non-relativistic case, the averaged spectrum coincides best with the spectrum emitted at $\theta=65{^{\circ}}$ and usually differs in shape for other smaller and larger inclinations.

Differences to the Angle-Average (`relxill_lp`)
-----------------------------------------------

Recent measurements [@wilms2001a; @fabian2002a; @bre06; @Dauser2012a; @fab12a; @gal11; @pon10; @bre11; @bre13; @dur11; @mil06; @rei08] have revealed a steeper emissivity at the inner parts of the accretion disk than the canonically assumed $r^{-3}$ dependence [@Shakura1973a]. This is well understood in the so-called “lampost” geometry [@mat91; @Martocchia1996a]. For this scenario, instead of coming from a corona around the inner regions of the disk [@haa93; @Dove1997a], the primary radiation illuminating the disk comes from a source placed on the rotational axis of the black hole [e.g., @dau13 for a detailed discussion]. Recent analyzes show that this model is capable of explaining the observed spectra in some systems [e.g., @Wilkins2011a; @dur11; @Dauser2012a].

Due to the success of this model, we investigate the effect of a rigorous and full treatment of the emission angle on the parameters describing the reflection and the relativistic blurring. Under lampost geometry, the emissivity profile is specified by the height $h$ of the primary source. While a large height usually implies a flat emissivity in the inner region, it steepens dramatically for decreasing height [e.g., @dau13]. This means that for a source located near the black hole (low $h$), the inner portion of the accretion disk more strongly dominates the reflection spectrum than for a source located far away (high $h$).

Figure \[fig:rel\_refl\_profiles\] shows the result of incorporating the angle-dependent reflection models from [xillver]{} into lampost geometry, synthesized with the relativistic blurring code into the end-product [relxill\_lp]{} [@dau13]. Each row represents a different inclination angle of the system (as indicated in the left panels). In all of the panels, each color represent a different height of the illumination source (as indicated in the middle panels). The left column shows the integrated spectra reflected from the disk over all considered heights of the source. Intuitively and evidently, the closer the proximity of the source to the black hole, the stronger the relativistic effects that distort the shape of the lines. The middle column contains the ratio of each spectrum to the solution predicted using the angle-averaged case. Significant differences can be seen across the entire energy spectrum. For both very low and very high inclinations, the discrepancies in the Fe K line region can be as large as 40%. As expected, the differences with respect to the angle-averaged solution are smaller for intermediate inclinations. The right column shows the distribution of emission angles for each configuration. In other words, these panels show the relative contribution of photons emitted from different angles to a spectrum corresponding to an observation at a particular inclination. As expected, the distribution peaks at the viewing angle, in particular when the source is located far away from the center. For cases where the height of the illumination source is low ($h\sim 2-10~{\ensuremath{r_\mathrm{g}}}$), the distribution can become very flat, and even more so when the disk is observed at large inclination angles. A very similar result can be seen for the models with $\log\xi=3$ (Figure \[fig:rel\_refl\_profiles\], middle column). Notice that when the distribution of the emission angle becomes very flat, the resultant spectrum is very similar to the angle-averaged solution (since in this case the observer is detecting photons coming from all angles with similar contributions, which has the effect of averaging the signal).

An important caveat to notice is that all the [xillver]{} models have been calculated implementing an incident irradiation at $45{^{\circ}}$ with respect to the disk normal. As discussed in [@gar10] and [@dau13], the incidence angle affects the specific intensity at the illuminated boundary of the slab, which influences to certain degree its ionization structure. The assumption is that the differences introduced by the incidence angle can be mimicked by a change in the ionization parameter. However, in the context of a lampost scenario these effects could have a larger impact on the simulated spectrum. A proper treatment of variable incidence angles requires a much larger calculation of reflection models which is beyond the scope of the present work, but it will be implemented in future versions of the model.

Estimating Bias from Angle-Averaged Reflection Results {#sec:sim}
======================================================

As evident from the previous section, significant differences exist when using the angle-averaged emissivity instead of the proper relativistically lensed treatment. However, prior to this work, the only available models for data fitting were angle averaged. Hence, we aim to determine to what extent this approximate treatment affects the fitted parameters as well as what can be gained by using the fully relativistic emissivity angle distribution.

In order to estimate the degree of bias from a typical observation, we produce a spectrum of an AGN showing relativistic reflection by simulating a 100ksec observation with *XMM-Newton* assuming an unblurred reflection component with flux comparable to [@bre06]. We employ the lampost version of our new model (`relxill_lp`) for the intrinsic, emitted reflection from which the observation is obtained. Subsequently, the simulated data is then fit by a commensurate (and common) angle-averaged model for X-ray reflection, by setting the flag [angleon=0]{}[^2].

Figure \[fig:sim\_refl\] shows the bias estimates obtained from the simulation. In general, the behavior of the fitting parameters can be divided into two classes, depending on the irradiation of the accretion disk. These classes are delineated by the emissivity profile as (a): showing strong, centrally-concentrated emissivity which predominantly fluoresces the innermost accretion disk (a high inner emissivity index), and (b): showing a flattened emissivity profile with a scaling below $r^{-3}$. Class (a) is readily understood as resulting from a low source height, whereas (b) represents a larger source height which provides a more uniform illumination of the inner disk. In the following we will interpret these results in detail. The values of our bias estimates, for both classes, can be found in Table \[tab:sim\].

In order to study class (a) in which the emissivity is steep (Figure \[fig:sim\_refl\], left panels), we use a source height of $h=3\,{\ensuremath{r_\mathrm{g}}}$. In general, it is evident that the iron abundance is most affected. In almost all the simulations the abundance is over-predicted and generally the effect is largest for low inclination angles. For a large spin ($a_*=0.99$) and an almost neutral disk ($\xi=1$), this effect can be more than a factor of two. For high inclinations this effect is reduced for all parameter configurations, while for $a_*=0.0$ and $\xi=1$ this turns into the opposite trend. Especially interesting is how the angle-averaged treatment influences spin determination. For the usual AGN observation (large spin, low ionization disk) there seem to be only very slight deviation from the input parameters. If the ionization is larger ($\xi=10^{3}$), then stronger biases arise: A large spin ($a_*=0.99$) will lead to an underestimate when fitting with angle-averaged models for low inclinations ($\theta<50^\circ$), while for larger inclinations ($\theta>60^\circ$), low spin values ($a_*=0.0$) can be overestimated (but with large uncertainties). The inclination itself is essentially unaffected by the angle-averaged treatment. Only in the case of very low inclinations ($\theta<40^\circ$) and high spin, the inferred inclination can be slightly overestimated (low ionization) or underestimated (high ionization), but the change is still quite small.

In the other instance (class b), Figure \[fig:sim\_refl\] (right panels) shows the results for a primary irradiating X-ray source placed at $h=50\,{\ensuremath{r_\mathrm{g}}}$. For this particular choice of parameters the spin cannot be very well constrained [see @dau13], and thus we analyze the effect of the averaged-angle treatment on the fitted height of the primary source, which is equivalent to determining the steepness of the emissivity profile. As in class (a), iron abundance is generally overestimated. However, the absolute bias is only around 20% which much lower for (b) compared to the bias obtained from the steep emissivity of (a). As before, this effect diminishes for higher inclinations. For the fitted inclination and the fitted height of the primary source, the degree of bias is very similar: A significant difference is only obtained for large inclinations ($\theta>60^\circ$), in the sense that at increasing values of (input) inclination, the determination of height and the inclination will be underestimated by angle-averaged models. This effect is again much stronger and more significant for a low ionization of the accretion disk.

Application to Observational Data {#sec:obs}
=================================

To illustrate how the discussed changes in the model influence the spectroscopic results of a real dataset, we study the [*Suzaku*]{} spectrum of the ‘bare’ Seyfert 1 galaxy Ark120. The term ‘bare’ refers to the weak amount of any intrinsic absorption, which simplifies the spectral analysis considerably [@pat11b]. The $\sim 100$ks exposure (ObsID:702014010) was taken in April 2007 and reduced following the instructions of the [*Suzaku*]{} ABC Guide using the newest available calibration. For the spectral analysis, the two front-illuminated CCDs were co-added and binned to a signal-to-noise of 10. The energies $0.7-10$ keV were noticed for fitting. The PIN data was used from $15-45$ keV and binned to a signal-to-noise of 5. The back-illuminated XIS chip (XIS1) was not considered for this analysis. @nar11 and @wal13 already showed that the spectrum of Ark 120 can be well described by a simple reflection model, including both a cold and blurred ionized reflector. Utilizing a similar model set-up to @wal13, as well as applying the same restrictions to the parameters, we first describe the spectrum with the angle-averaged version of our new model [relxill]{} (by setting the flag [angleon=0]{}) for the blurred ionized reflector, and the cold reflector with a simple [xillver]{} table (with $\xi$ fixed to 1). The key parameters of this fit are presented in Table 2. In particular, we derive a spin of $a_*=0.674_{-0.203}^{+0.118}$ and an inclination of $i=45.0_{-2.7}^{+5.3}$ deg. To see the influence of the angle-dependence we then set the flag [angleon=1]{} for the blurred reflector, and replace the cold reflector with the angle dependent version of the [xillver]{} table. The inclination angles of the two reflectors were linked for fitting. Similar to the angle-averaged fit the quality is very good and the spectrum is well described by this model (see Figure \[ark120\]). The angle-dependence improves the spin constraint to $a_*=0.655_{-0.126}^{+0.122}$ and also improves the constraint on the inclination slightly ($i=45.3_{-2.4}^{+4.8}$ deg). This is in agreement with the difference from the angle-averaged model predicted in Section 4.1, which is of order 10%, and modest compared to the measurement errors. Error bars were calculated to a 90% confidence level.

To better illustrate the improvement achieved by our new model in constraining the physical parameters, we have produced the contour plots for the spin parameter $a_*$ and the inclination angle. In Figure \[fig:contour\] solid lines show the 69%, 90%, and 99% confidence regions of these two parameters for the fit using the angle-resolved ([angleon=1]{}) option in the new model [relxill]{}. Dashed-lines show the same contours obtained while fitting with the angle-averaged solution ([angleon=0]{}). The crosses show the position of the minimum value of $\chi^2$ for each case. Although both fit recover similar best-fit values, the accuracy in the determination of both the spin and inclination is significantly higher with the angle-resolved model. One can see that both parameters are better determined at 90% level with the new model than a 69% level with the old model. With this particular test case we have shown the capabilities of the new model in the interpretation of X-ray spectra from accreting sources, as it constitutes a more consistent and physical description of the reflected spectra from accretion disks around black holes.

Evidently, this is just one particular example of the potential of the new model in fitting X-ray observations significantly improved: the errors in both spin and inclination are better

Summary and Conclusions {#sec:con}
=======================

The solution of the angular distribution of X-ray reflection from an accretion disk has been calculated with our reflection code [xillver]{}. The reflected spectrum at any given angle shows clear discrepancies of up to $\sim 25 \%$ when compared to the angle-averaged solution. This is observed to occur for all angles, and at any ionization parameter. In particular, scattered continuum photons and line photons respond differently to the angular effects, since the optical depth can be significantly different at different energies. Notably, the Fe K emission line complex becomes weaker for high viewing angles, which agrees with earlier claims of the decrease of the line equivalent width when observed at grazing angles [e.g. @bas78; @mat92; @mar00; @rey04]. The angular effects are introduced in the reprocessed radiation mainly by changing the effective optical depth $\tau_\mathrm{eff} = \tau/\mu$, where $\mu$ is the cosine of the viewing angle. The angular effects are much more extreme at energies where the photoelectric absorption dominates over the Compton scattering. At high energies ($\gtrsim 20$ keV), Compton scattering is the only source of opacity, and thus the angular effects are independent of the ionization of the gas. An important consequence of these effects is that they change the slope of the reflected spectrum as the angular-dependency acts differently at different energies.

The new reflection models predict a limb-brightening law, i.e., the reflected flux integrated over all energies becomes larger at high inclination angles, which agrees well with hydrostatic calculations of [@roz08] and [@roz11]. However, its functional form does not agree with the commonly used law $I \propto
\ln(1+1/\mu)$ [@haa93], instead, our calculations show a profile less steep. When the atmosphere is highly ionized (log $\xi > 3$), the profile becomes even flatter, approaching to the isotropic case. By carrying out calculations of a pure scattering atmosphere, we have shown that this is expected as Compton scattering is essentially isotropic for the reflected spectra of irradiated plane-parallel slabs. In the context of [*unblurred reflection*]{}, we expect that these models will be most relevant in analyzing data from astrophysical systems at either very low or very high inclination angles. A new [xillver]{} table which includes the inclination as a model parameter is provided in the proper format[^3] for its implementation via the [atable]{} model in [xspec]{}.

The table of reflection spectra mentioned above is intended to be used to model an unblurred reflection component, as is the case of the distant cold reflection observed in many AGN. However, it is important to emphasize that when relativistic effects need to be taken into account, a convolution with blurring models in the usual fashion (e.g. [relconv $\times$ xillver]{}) is in fact incorrect, as the proper link of the two models needs to be done intrinsically at each radial zone. Instead, we advise observers to implement the newly developed model [relxill]{}[^4], which properly integrates the two models. An additional version of this model ([relxill\_lp]{}), in which the illumination is prescribed under a [*lampost*]{} scenario, is also provided. The relativistic effects introduce additional directionality dependences in the final reflected spectrum. In general, deviations from the angle-averaged solution can be as large as $20\%$, and like in the unblurred reflection case, lines and continuum respond in a different way. The main effect is that, at any given inclination, a distant observer will not only see those photons that are emitted into its line of sight, but also those emitted at different directions, but for which the light bending effects are enough to modify their path so that they eventually point toward the observer (for instance, at grazing angles, an observer can even see photons emitted behind the black hole, due to the extreme space-time curvature).

Models calculated in a lampost geometry, where the illumination of the disk is mostly prescribed by the height of the source over the accretion disk, also differ significantly from their angle-averaged counterparts. We found that for either very low or very high inclinations, the discrepancies in the Fe K line band can be as large as $40\%$. It is worthwhile to mention that all these relative differences are for the reflection spectrum only. In real observations, the observed spectrum is a combination of the reflected and the direct component which illuminates the disk. In many cases the reflection spectrum may only contribute $10-20\%$ of the $2-10$ keV flux of the direct continuum, which dilutes the reflection features and consequently the deviations seen between the new and the old versions of the models. In cases where the height of the illumination source is low ($h\sim 2-10~r_g$), the distribution of emission angles that are actually observed at a given inclination becomes very flat. This degree of angular mixing is most pronounced when the disk is observed at high inclination angles. Accordingly, when relativistic effects are extreme, an observer is detecting photons from almost all angles with similar contribution, which essentially acts as an averaging of the reflected spectra over all angles.

We have also carried out a detailed error analysis of the physical parameters derived from fitting these models to simulated data. Our main goal is to show an estimate of the bias introduced by the implementation of the approximate (and commonly used) angle-averaged solution for the reflected spectra. By looking at two distinctive cases, i.e., steep and flat irradiation of the disk, we found that the Fe abundance tends to be over-predicted in almost all cases, but in particular for low inclinations. For a large spin and low ionization, the differences can be as large as a factor of two. This could explain the high abundances usually observed in the spectra of many AGN [see for example, @rey97]. Angular effects also influence the spin determination, although the effect is not dramatic for typical parameters. In particular, for large spin and low ionization the spin recovered from the fits is underestimated at low inclinations ($\theta < 50{^{\circ}}$), while for larger inclinations originally low spins ($a_*=0.0$) are overestimated. The bias on the inclination and height of the primary source is only significant in the case of low spin and large inclination angles, for which both quantities can be underestimated. Oddly, the determination of the inclination angle is mostly unaffected when the spin is large.

We have also included an illustrative example of the application of the new model to real observational data. For this, we analyzed the [*Suzaku*]{} spectrum of the Seyfert 1 galaxy Ark 120. The same fit is performed using the angle-averaged solution and the complete angle-dependent model. The differences in the recovered parameters are modest. Nevertheless, the differences in the constrain of the physical parameters is quite significant. The accuracy of the new model in the determination of both the inclination and the spin parameter $a_*$ at a 90% confidence level is superior than the old, angle-averaged model at 69% confidence level. This test case demonstrates the capabilities of the new model, which is more consistent and physical than previous versions. It also illustrates its relevance in the interpretation of observational data from accreting sources, and we expect it will have an even greater impact on analyzing higher quality data from the new generation of observatories such as [*NuSTAR*]{} and [*Astro-H*]{}.

The models presented here are a new step toward a more physical and self-consistent representation of the X-ray reflection from accretion disks. However, additional sources of systematics that could affect the retrieval of physical parameters (e.g. black hole spin) need still to be further investigated. Although recent studies suggest that for some sources the emission region is compact [e.g. @dau13], the exact geometry of the primary source of X-rays, whether it is a corona or the base of jet, is not very well understood. The assumptions made on this respect affect the way X-rays illuminate the the accretion disk. By providing a lamppost version of our reflection model we are taking the first steps to unveil the geometry of the illuminating source. However, further studies [e.g., timing measurements; see @zog10; @kar13a; @kar13b], are required to accurately determine the illumination profile and the physical origin of the primary source itself. The vertical ionization structure of the disk is treated in detail with the models presented here, however, the same ionization parameter is assumed for all radii. The radial profile of the ionization depends on both the illumination and the gas density profiles. This will be subject of future studies. Other sources of systematics are not directly related to the reflection spectrum. For instance, the presence of a warm absorber in several AGN is known to modify the prediction of the fits depending on how this component is being treated [e.g., MCG-6-30-15, @bre06]. However, observations carried out with the new observatory [*NuSTAR*]{} in combination with [*XMM-Newton*]{} show that we are currently on the edge of being able to distinguish between blurred reflection and absorption [with NGC 1365 being the best example so far, @ris13]. The improved physical description of the illuminated gas is expected to play a major role to disentangle these components by reducing the uncertainties in the X-ray reflection spectrum.

The calculations presented here were performed in the Odyssey cluster of the Research Computing Facilities of the Harvard University. JG and JEM acknowledge the support of NASA grant NNX11AD08G. JFS was supported by NASA Hubble Fellowship grant HST-HF-51315.01. CSR thanks support from NASA under grant NNX10AE41G.

[lcc|rrrr|rrrrr|r]{} & & & &\
& $a_*$ & $h$ & $\theta=20^\circ$& $\theta=40^\circ$& $\theta=60^\circ$& $\theta=80^\circ$& $\theta=20^\circ$& $\theta=40^\circ$& $\theta=60^\circ$& $\theta=80^\circ$\
$A_\mathrm{Fe}$ & 0.99 & 3 & $1.14\pm0.07$& $1.00^{+0.15}_{-0.09}$& $0.42\pm0.07$& $0.02\pm0.03$& $0.57^{+0.09}_{-0.12}$& $0.44^{+0.10}_{-0.09}$& $0.27\pm0.11$& $-0.00^{+0.10}_{-0.02}$\
$a_*$ & 0.99 & 3 & $0.00\pm0.01$& $-0.00^{+0.01}_{-0.02}$& $0.01^{+0.00}_{-0.01}$& $0.01^{+0.00}_{-0.01}$& $-0.12^{+0.05}_{-0.04}$& $-0.12^{+0.06}_{-0.05}$& $-0.02\pm0.02$& $-0.00^{+0.01}_{-0.27}$\
$\theta$ & 0.99 & 3 & $1.97^{+0.69}_{-0.67}$& $0.44^{+0.31}_{-0.36}$& $-0.39^{+0.18}_{-0.20}$& $-0.23^{+0.20}_{-0.16}$& $-1.92^{+1.51}_{-1.00}$& $-2.03^{+1.08}_{-1.16}$& $-0.09^{+0.43}_{-0.48}$& $0.00^{+0.00}_{-1.86}$\
$A_\mathrm{Fe}$ & 0.00 & 3 & $0.21\pm0.02$& $0.18^{+0.02}_{-0.03}$& $0.09\pm0.02$& $-0.16\pm0.03$& $0.16\pm0.04$& $0.13\pm0.04$& $0.08\pm0.05$& $0.07^{+0.11}_{-0.09}$\
$a_*$ & 0.00 & 3 & $-0.06^{+0.04}_{-0.03}$& $-0.06^{+0.08}_{-0.05}$& $-0.05^{+0.09}_{-0.06}$& $0.13^{+0.17}_{-0.36}$& $-0.00^{+0.12}_{-0.07}$& $-0.07^{+0.17}_{-0.15}$& $0.56^{+0.20}_{-0.57}$& $0.97^{+0.03}_{-1.09}$\
$\theta$ & 0.00 & 3 & $0.09^{+0.31}_{-0.39}$& $-0.04^{+0.14}_{-0.20}$& $0.27^{+0.20}_{-0.19}$& $-4.79^{+4.16}_{-1.13}$& $-0.34^{+1.30}_{-1.34}$& $0.11^{+0.49}_{-0.59}$& $0.81^{+0.91}_{-1.06}$& $-4.91^{+1.03}_{-2.20}$\
$A_\mathrm{Fe}$ & 0.99 & 50 & $0.25\pm0.02$& $0.14\pm0.02$& $-0.03\pm0.03$& $-0.07^{+0.03}_{-0.02}$& $0.16\pm0.03$& $0.12\pm0.03$& $-0.00^{+0.08}_{-0.01}$& $-0.06\pm0.03$\
$a_*$ & 0.99 & 50 & $-1.99^{+0.37}_{-0.01}$& $-1.99^{+1.01}_{-0.00}$& $0.01^{+0.00}_{-1.76}$& $0.01^{+0.00}_{-0.03}$& $-1.99^{+2.00}_{-0.01}$& $-1.98^{+1.99}_{-0.01}$& $0.01^{+0.00}_{-2.00}$& $0.00^{+0.01}_{-2.00}$\
$\theta$ & 0.99 & 50 & $0.29\pm0.20$& $-0.07^{+0.44}_{-0.42}$& $0.43^{+2.95}_{-1.99}$& $-21.22^{+1.60}_{-1.23}$& $0.62^{+0.79}_{-0.82}$& $-0.06^{+1.81}_{-1.33}$& $0.80^{+4.90}_{-3.74}$& $-12.55^{+12.55}_{-5.09}$\
$A_\mathrm{Fe}$ & 0.00 & 50 & $0.27\pm0.02$& $0.15\pm0.02$& $-0.03\pm0.03$& $-0.08^{+0.03}_{-0.02}$& $0.17\pm0.03$& $0.12\pm0.03$& $0.00^{+0.07}_{-0.02}$& $-0.07\pm0.03$\
$a_*$ & 0.00 & 50 & $-1.00^{+0.49}_{-0.00}$& $-1.00^{+0.90}_{-0.00}$& $1.00^{+0.01}_{-1.55}$& $1.00^{+0.00}_{-0.04}$& $-1.00^{+2.00}_{-0.00}$& $-1.00^{+2.00}_{-0.01}$& $0.77^{+0.23}_{-1.78}$& $0.99^{+0.01}_{-1.99}$\
$\theta$ & 0.00 & 50 & $0.25\pm0.20$& $0.10^{+0.46}_{-0.41}$& $-0.86^{+1.95}_{-1.78}$& $-22.10^{+2.19}_{-0.86}$& $0.72^{+0.80}_{-0.82}$& $0.03^{+1.84}_{-1.35}$& $1.20^{+5.80}_{-3.65}$& $-16.09^{+5.17}_{-3.49}$\

[lccc]{} \[tfit\] & $q_{in}=q_{out}$ & $4.4_{-1.2}^{+2.0}$ & $4.8_{-1.1}^{+1.8}$\
[relxill]{} & $a$ & $0.674_{-0.203}^{+0.118}$ & $0.655_{-0.126}^{+0.122}$\
[relxill]{} & $i$ (deg) & $45.0_{-2.7}^{+5.3}$ & $45.3_{-2.4}^{+4.8}$\
[relxill]{} & $R_{in}$ (ISCO) & $1$ & $1$\
[relxill]{} & $R_{out}$ (ISCO) & $400$ & $400$\
[relxill]{} & $z$ & $0.0327$ & $0.0327$\
[relxill]{} & $\Gamma$ & $2.17_{-0.01}^{+0.02}$ & $2.17_{-0.01}^{+0.02}$\
[relxill]{} & log $\xi$ & $1.03_{-0.26}^{+0.16}$ & $0.84_{-0.10}^{+0.24}$\
[relxill]{} & $A_{Fe}$ & $1.54_{-0.45}^{+0.44}$ & $1.78_{-0.40}^{+0.39}$\
[relxill]{} & [angleon]{} & 0 & 1\
[relxill]{} & $N (10^{-4})$ & $1.81_{-0.29}^{+0.32}$ & $2.07_{-0.2}$\
[xillver]{} & $N (10^{-4})$ & $1.55_{-0.40}^{+0.42}$ & $0.14_{-0.04}^{+0.07}$\

[^1]: See [@dau10; @dau13] for more details on the calculation and the code.

[^2]: Note that setting [angleon=0]{} is equivalent to implement a relativistically blurred reflection in the common fashion such as `relconv_lp` $\times$ `xillver`

[^3]: <http://hea-www.cfa.harvard.edu/~javier/xillver/>

[^4]: <http://www.sternwarte.uni-erlangen.de/research/relxill/>
