---
title: 'REFLIONX: X-ray Reflection Model'
url: https://www.emergentmind.com/topics/reflionx
type: topic
---

# REFLIONX: X-ray Reflection Model

Searching arXiv for recent and relevant papers on REFLIONX and its use in reflection and reverberation modeling.
REFLIONX is a non-relativistic X-ray reflection model for optically thick, ionized accretion-disc material illuminated by an external continuum. In the literature summarized here, it is consistently treated as a rest-frame reflection code that computes the reflected continuum and line emission from a constant-density slab, and is then either convolved with a relativistic transfer kernel for spectroscopy or embedded in reverberation frameworks for timing analyses. Across active galactic nuclei, black-hole X-ray binaries, and neutron-star comparisons, REFLIONX functions both as a widely used baseline model and as a source of identifiable systematics relative to XILLVER-based and RELXILL-based alternatives, especially in the soft X-ray band and in high-density regimes [2401.06976].

## 1. Definition and physical scope

REFLIONX, in the usage described across these studies, computes the spectrum produced when a power-law X-ray continuum illuminates a constant-density, optically thick slab representing the disc surface. It includes continuum reflection, photoelectric absorption edges, fluorescent emission lines, and Compton scattering, and it is commonly described as a model for ionized reflection from a semi-infinite slab [1909.01897]. In the AGN reverberation literature, it is the reflection engine underlying `kynrefrev`, where it supplies the rest-frame reflection spectrum subsequently processed by relativistic transfer functions and timing calculations [2401.06976].

The model is parameterized by the ionization parameter $\xi$, iron abundance $A_{\rm Fe}$, the photon index $\Gamma$ of the illuminating continuum, and a normalization; in several applications the high-energy cutoff is fixed at $300\,\mathrm{keV}$ or represented through a tied Comptonization continuum [1909.01897]. The ionization parameter is written in the supplied studies as either
\[
\xi = \frac{L}{n r^2}
\]
or
\[
\xi = \frac{4\pi F_{\rm X}}{n_e},
\]
with the notation varying by paper but the physical meaning unchanged: $\xi$ measures the ratio of ionizing flux to gas density [1909.01897] [2305.05914] [2407.07362].

Several papers emphasize two structural assumptions. First, REFLIONX is based on a constant-density slab rather than a vertically stratified atmosphere [1909.01897]. Second, it is angle-averaged rather than angle-resolved, so inclination dependence is typically introduced by an external relativistic convolution such as `kerrconv`, `relconv`, or `relconv_lp` [1909.01897] [2007.15914]. This differs from XILLVER-based frameworks, which are described as angle-dependent and built on a richer atomic database [2007.15914].

## 2. Implementation in spectral and reverberation models

In relativistic spectroscopy, REFLIONX is not itself a relativistic model. It is convolved with transfer kernels that impose Doppler shifts, gravitational redshift, and disc geometry. In the Mrk 478 multi-epoch analysis, the exact XSPEC construction is
\[
\texttt{tbabs} \times \texttt{const\_1} \times \big[\texttt{zgauss} + (\texttt{cflux} \times \texttt{cutoffpl}) + (\texttt{const\_2} \times \texttt{cflux} \times \texttt{kerrconv} \times \texttt{reflionx})\big],
\]
where `cutoffpl` represents the primary continuum, `kerrconv` supplies relativistic blurring, and `const_2` encodes the reflection fraction [1909.01897]. In black-hole binary work, analogous constructions use `relconv*reflionx_hd`, `relconv*reflionx_hdv2`, or `relconv_lp*reflionx_HD_nthcomp_v2`, tying the illuminating continuum parameters between the direct and reflected components [2407.07362] [2109.05380].

In reverberation studies, REFLIONX enters through public lamp-post timing codes. In `kynrefrev`, the workflow described in the data is: assume a lamp-post corona at height $h$ above a Novikov–Thorne disc, compute the rest-frame reflection spectrum with REFLIONX, apply relativistic ray-tracing and time delays through the KYN transfer function, integrate over soft and hard bands, and compute lag-frequency spectra from the Fourier cross-spectrum [2401.06976]. In the general-relativistic spectral–timing model of Chainakun and collaborators, REFLIONX is used locally on each disc zone after computing the incident flux and local ionization, and the blurred reflected spectra are propagated to the observer to generate both mean spectra and lag-energy spectra [1605.01300].

The timing formalism uses
\[
C(f) = S^*(f)\,H(f),
\qquad
\tau(f) = \frac{\arg C(f)}{2\pi f},
\]
with dilution handled through the reflection fraction in each band [2401.06976]. Because REFLIONX fixes the shape and relative strength of the reflected spectrum in the soft and hard bands, it directly affects the inferred reverberation lag amplitudes and the mapping between observed delays and geometry [2401.06976] [1506.04524].

## 3. REFLIONX in comparison with XILLVER and RELXILL

A recurring theme is that REFLIONX and XILLVER-based models can fit similar datasets while implying different physical interpretations. The most systematic comparison is the AGN reverberation study of `kynrefrev` versus `kynxilrev`, where the only intended difference is the underlying reflection code: `kynrefrev` uses `reflionx.mod`, whereas `kynxilrev` uses `xillverD-5.fits` [2401.06976]. That work attributes their timing differences chiefly to soft-band spectral differences: XILLVER is described as more absorbed below $\lesssim 1$ keV, implying a lower soft-band reflection fraction and stronger lag dilution for the same geometry [2401.06976].

The same issue appears in simultaneous spectral–timing analyses of Mrk 335. There, REFLIONX and XILLVER both fit the Fe K band and hard reverberation lags well, but they extrapolate differently below 2 keV. The authors state that the soft-excess mismatch between `revb ⊗ reflionx` and `revb ⊗ xillver` is attributable to differences in atomic data and low-energy absorption, and that these systematic differences exceed the statistical uncertainties [1506.04524]. In that study, for a constant-density disc, REFLIONX and XILLVER yield markedly different combinations of black-hole mass and source height for otherwise similar-quality fits, reinforcing the claim that inferred geometry can be model-dependent [1506.04524].

Spectral fitting studies reach comparable conclusions. In Mrk 478, REFLIONX and RELXILL both fit the X-ray data statistically well, but REFLIONX implies a low-ionization, power-law-dominated solution with $R_{0.1-100}\sim 0.7$–$1.3$, whereas RELXILL implies a highly ionized, reflection-dominated solution with $R_{0.1-100}\sim 2.0$–$3.5$ [1909.01897]. The authors judge the REFLIONX interpretation to be more consistent with the source’s normal $\alpha_{\rm ox}$, showing that the two models encode different disc ionization regimes and different balances between direct and reflected emission [1909.01897].

The contrast persists in broader benchmarking. In the study of GRS 1915+105 and MCG–6–30–15, RELCONV$\times$REFLIONX, RELCONV$\times$XILLVER, and RELXILL yield similar parameters for the simpler stellar-mass case, but differ significantly for the AGN case, where REFLIONX tends to prefer lower spin and lower inclination and gives worse overall fits [2007.15914]. In the AGN-host alignment study, however, simulations show that the inclination difference induced by switching between RELXILL and REFLIONX is usually only a few degrees, with mean rms deviations of $2.22^\circ$, $1.50^\circ$, and $1.45^\circ$ for input inclinations of $30^\circ$, $50^\circ$, and $70^\circ$ respectively [1601.03090]. This suggests that model dependence is observable but context dependent: it is small for some inclination measurements, yet substantial for ionization, reflection fraction, mass, height, and soft-band interpretation.

A recent broadband Seyfert analysis sharpens this point further. For PG 1426+015, high-density REFLIONX-based models fit the soft excess, Fe K complex, and Compton hump without ad hoc modification, whereas XILLVER-based models fail unless the O VIII feature is artificially suppressed [2509.13411]. The authors attribute the divergence to a stronger soft continuum and weaker O VIII line in REFLIONX, partly linked to different resonant-line transfer treatments [2509.13411]. This suggests that, in bare Seyferts with strong soft excesses, the choice between REFLIONX and XILLVER is not merely technical but can determine whether a reflection-only interpretation remains viable.

## 4. High-density extensions and the iron-abundance problem

A major development in REFLIONX usage is the introduction of high-density variants such as `reflionx_hd`, `reflionx_hdv2`, and `reflionx_HD_nthcomp_v2`, motivated by the inadequacy of standard low-density reflection grids for black-hole binaries [2305.05914] [2407.07362]. These models extend the electron density well above the canonical $10^{15}\,\mathrm{cm^{-3}}$ used in standard reflection tables and permit explicit tests of the degeneracy between iron abundance and density [2305.05914].

In MAXI J0637-430, standard `relxillCp` requires $A_{\rm Fe} \approx 10\,A_{\rm Fe,\odot}$, even when density is allowed to rise to $10^{20}\,\mathrm{cm^{-3}}$ within the high-density RELXILL grid [2305.05914]. Replacing that with `relconv*reflionx_hd`, fixing iron abundance to solar, and allowing the density to vary yields `logN = 22.0_{-0.6}^*`, corresponding to $n_e \sim 10^{22}\,\mathrm{cm^{-3}}$, with a good fit quality $\chi^2/\nu = 1181.94/1065$ and no need for extreme iron enrichment [2305.05914]. The authors interpret this as evidence that the supersolar iron abundance inferred from low-density models is an artifact of the density assumption, while the black-hole spin remains robust across model choices [2305.05914].

The same pattern appears in XTE J2012+381. Low-density RELXILL fits produce $A_{\rm Fe}>8.39\,A_{\rm Fe,\odot}$, and even high-density `relxillCp` at $10^{20}\,\mathrm{cm^{-3}}$ still demands $A_{\rm Fe}>8.8$ [2407.07362]. In contrast, `relconv*reflionx_HD_nthcomp_v2` yields $n_e > 4.39\times10^{21}\,\mathrm{cm^{-3}}$ with $A_{\rm Fe}=2.51_{-1.46}^{+1.48}\,A_{\rm Fe,\odot}$ in one configuration, and $n_e > 3.11\times10^{21}\,\mathrm{cm^{-3}}$ with $A_{\rm Fe}=2.12_{-0.50}^{+1.01}\,A_{\rm Fe,\odot}$ in another [2407.07362]. The paper therefore uses REFLIONX primarily as a density diagnostic rather than as the primary spin-measurement tool [2407.07362].

In MAXI J1348-630, REFLIONX-based high-density models again reduce the required iron abundance to around solar while returning densities of $10^{20.3-21.4}\,\mathrm{cm^{-3}}$ [2109.05380]. The preferred `reflionx_hdv2` implementation yields $\log n_{\rm disc}=20.3$–$20.7$ across most epochs, with $A_{\rm Fe}\sim0.7$–$1.4$ except for one higher value in E3, and tracks a physically plausible evolution of inner radius and ionization through the outburst [2109.05380]. In Cygnus X-1, a modified variable-density `reflionx_hd` gives the especially clear result
\[
n_e = (3.98^{+0.12}_{-0.25}) \times 10^{20}\ {\rm cm^{-3}},
\]
while removing the need for iron abundances of $>9.96$ or $10.6^{+1.6}_{-0.9}$ times solar required by standard RELXILL and REFLIONX fits [1801.07267].

These studies converge on a common inference. High density changes the reflected continuum through enhanced free–free heating and associated soft excess, modifies line and edge strengths, and can mimic the effects otherwise ascribed to extreme $A_{\rm Fe}$ [2305.05914] [2109.05380] [1801.07267]. A plausible implication is that low-density reflection fits in bright black-hole binaries can systematically bias iron abundance and, through parameter degeneracies, bias inner radius and inclination as well.

## 5. Reverberation, dilution, and model-dependent geometry

REFLIONX has a distinct role in reverberation analyses because timing constraints are sensitive not only to light-travel time but also to the spectral partition between direct and reflected flux. In the comparative KYN study, both REFLIONX-based `kynrefrev` and XILLVER-based `kynxilrev` reproduce the observed lag–mass relation by construction, but the inferred coronal properties differ [2401.06976]. For the REFLIONX-based fits, the lag–height relation is given as
\[
\log \tau = -2.05(\pm 0.39) + 0.35 (\pm 0.03)\, h,
\]
with $R^2=0.53$, while the preferred source heights cluster around $h\sim5$–$15\,r_{\rm g}$ with a typical value near $10\,r_{\rm g}$ [2401.06976]. REFLIONX also yields a strong correlation between source height and luminosity, with Spearman $r_s(h,L)=0.80$, a trend absent in the XILLVER-based model [2401.06976].

The authors interpret the main difference as a dilution effect. Because XILLVER is more absorbed in the $0.3$–$0.8$ keV band, soft reverberation lags are more diluted, forcing the XILLVER-based timing model toward lower heights and lower black-hole masses to match the same lag-frequency observables [2401.06976]. REFLIONX, having less soft-band absorption in this comparison, can accommodate larger intrinsic heights while remaining consistent with the observed lag amplitudes and phase-wrapping frequencies [2401.06976].

The simultaneous spectral–timing modeling of Mrk 335 reaches a related conclusion. There, REFLIONX and XILLVER fit the hard band comparably, but the soft band and therefore the Fe L lag interpretation diverge [1506.04524]. More generally, the AGN reverberation model of 2016 uses REFLIONX to compute the local reflection spectrum on each disc element, then derives the energy-dependent response function $\psi(E,t)$ and reflected response fraction $R(E)$, which together determine the lag-energy spectrum including dilution [1605.01300]. The model shows that ionization gradients encoded through REFLIONX can naturally produce lag-energy dips around 3 keV and 7–10 keV if the inner disc is highly ionized and the source height exceeds $5\,r_g$ or the disc becomes colder at larger radii [1605.01300].

This establishes a technical point of broader importance. In reverberation work, REFLIONX is not merely a spectral ingredient. It sets the band-dependent reflection fractions that enter the diluted light curves and therefore the mapping from timing data to source geometry [2401.06976] [1605.01300]. Consequently, model-dependent differences in low-energy opacity and line strengths can propagate into fitted masses, heights, and scaling relations even when the relativistic transfer formalism is identical.

## 6. Limitations, revisions, and broader model context

The supplied literature repeatedly identifies REFLIONX as foundational but simplified. Its principal limitations are the constant-density slab assumption, angle-averaged output, and older or more limited atomic data relative to XILLVER-based families [2007.15914] [2512.12728]. In the DAO benchmark study, REFLIONX is contrasted with DAO and XILLVER as one of the standard non-relativistic reflection models, but it is noted to rely on a Gaussian-approximated Compton redistribution function and a custom atomic dataset rather than the newer XSTAR-based infrastructure used in XILLVER and DAO [2512.12728]. The same paper states that REFLIONX is angle-averaged and that classic tables are less flexible with respect to arbitrary illuminating spectra [2512.12728].

At the same time, REFLIONX has been modified in several directions. High-density extensions have already been noted. Another important modification appears in the returning-radiation study, where the underlying REFLIONX code, rather than its public power-law tables, is used with the actual incident spectrum at each radius: the sum of lamp-post illumination and returning radiation [2401.05582]. That work shows that standard relativistic models modified only to read REFLIONX tables, while still assuming simple power-law illumination, can misestimate spin, height, ionization, and iron abundance when returning radiation is strong, especially for very high spin and low coronal height [2401.05582]. The central lesson is that REFLIONX as a microphysical solver is more general than many public tables imply, but relativistic fitting frameworks may still inherit biases if they simplify the local illumination spectrum [2401.05582].

The neutron-star reflection paper adds another contextual refinement. It treats REFLIONX as a workhorse model originally designed for black-hole systems illuminated by cutoff power laws, and argues that neutron-star soft states require blackbody-illuminated reflection instead [2111.12838]. New models such as `xillverNS` and `relxillNS` are presented there as successors for that regime, while REFLIONX-based blackbody variants remain useful comparators [2111.12838]. This suggests that the domain of validity of REFLIONX depends not only on geometry and density but also on whether the assumed form of the illuminating continuum matches the astrophysical source.

Across the supplied studies, the most stable conclusion is not that REFLIONX is invalid, but that it is model-defining: it encodes specific assumptions about density, angular structure, incident spectrum, and atomic microphysics that materially affect inferred quantities. It remains one of the most widely used reflection models in relativistic X-ray spectroscopy, yet the modern literature treats its outputs as subject to cross-checking against alternative reflection families and to reinterpretation when high-density, soft-band, or returning-radiation physics is important [2401.06976] [1909.01897] [2007.15914] [2512.12728].

Source: https://www.emergentmind.com/topics/reflionx