Modeling X-ray reflection spectra from disk returning radiation: Application to 4U 1630-47

arXiv:2601.14860 · astro-ph.HE · Submitted 2026-01-21 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: We're looking at the paper: "Modeling X-ray reflection spectra from disk returning radiation".

Jocelyn: Here's the gist of the paper. The gist: A new table model for reflection spectroscopy that self-consistently accounts for returning radiation has been presented and applied to 4U 1630–47, revealing a rapidly spinning black hole with specific physical parameters.

Vera: Model Development The authors developed a new table model for reflection spectroscopy that "self-consistently accounts for the returning radiation" <ref:2601.14860#pg4>.

Jocelyn: To isolate this effect, they adopted the "standard disk-corona configuration but disable the corona, allowing the reflection spectrum to be produced solely by self-irradiation of the disk" <ref:2601.14860#pg3>.

Vera: The model utilizes an updated version of the ziji code and treats returning radiation as "returning thermal radiation (purple arrows in Fig.

Paper summary: Vera: Now, let’s look at what they actually found when they applied this new model to the black hole X-ray binary 4U one thousand six hundred thirty-forty-seven during its soft state <ref:2601.14860#pg4>.

Jocelyn: They used simultaneous observations from both NICER and NuSTAR, which gives them a solid set of data to test this new idea <ref:2601.14860#pg5>.

Subrahmanyan: The analysis suggests that the model fits the data quite well, which implies returning radiation from the disk could potentially explain those reflection features they observed in this source <ref:2601.14860#pg3>.

Vera: Specifically, their best-fit model, called zijiRetRad, points toward a rapidly spinning black hole with a spin parameter around zero point nine nine <ref:2601.14860#pg5>.

Jocelyn: And they also found an inclination angle estimate for the disk of about fifty-three degrees <ref:2601.14860#pg5>.

Subrahmanyan: Those specific constraints on spin and inclination are important because they tell us something concrete about how the inner accretion flow is shaped near a black hole, which is a major theoretical link <ref:2601.14860#pg2>.

Conclusion: Vera: To wrap up this paper "Modeling X-ray reflection spectra from disk returning radiation: Application to 4U one thousand six hundred thirty-forty-seven" we've seen how they use a new table model that focuses purely on the radiation bouncing back from the disk itself.

Jocelyn: Yeah, we’ve seen how this new table model fits a source called 4U one thousand six hundred thirty-forty-seven pretty well, and it seems to point toward some specific physical parameters for that system <ref:2601.14860#pg5>.

Subrahmanyan: The authors are focused on how this self-consistent accounting of returning radiation changes the standard reflection models we usually use <ref:2601.14860#pg3>.

Vera: They built this whole new table model to isolate just the effect of that returning radiation from the disk, without any corona involved <ref:2601.14860#pg3>.

Jocelyn: That’s a big move because usually, we have to invoke a corona to get that kind of hard X-ray excess you see in these spectra <ref:2601.14860#pg3>.

Vera: They found a rapidly spinning black hole, with a spin parameter around zero point nine nine, which is quite high for these types of systems <ref:2601.14860#pg5>.

Jocelyn: And they’re getting an inclination angle estimate around fifty-three degrees for the disk as well <ref:2601.14860#pg5>.

Subrahmanyan: That level of spin and that inclination constrain the geometry of the inner accretion flow in a way that links directly to how matter spirals onto the black hole <ref:2601.14860#pg2>.

Vera: The mass accretion rate they found is about fifteen percent of Eddington, which gives us a clear picture of how much fuel is actually going into the system right now <ref:2601.14860#pg5>.

Jocelyn: They also mentioned an inferred source distance around eight point two kiloparsecs, which is a bit lower than some of the literature values we’ve seen <ref:2601.14860#pg4>.

Subrahmanyan: That distance measurement ties back into how you scale everything else up when you try to connect this observation to the larger cosmological structure formation picture <ref:2601.14860#pg5>.

Vera: So, what does this mean for us looking at these X-ray binaries? It means returning radiation isn't just noise; it’s a key component we need to model accurately <ref:2601.14860#pg1>.

Jocelyn: Exactly. We need models that self-consistently handle the disk's own radiation, not just assume something else is doing the heavy lifting in those reflection spectra <ref:2601.14860#pg3>.

Subrahmanyan: If this iterative reflection mechanism is real, then future observations will need to be very precise about the energy dependence of those reflected features to confirm this physical picture <ref:2601.14860#pg8>.

Vera: The authors are pushing that returning radiation alone can account for the observed spectra, suggesting it’s a fundamental part of the high-soft state physics <ref:2601.14860#pg9>.

Jocelyn: That means we should start thinking about including disk self-irradiation as a standard part of the reflection physics right away, instead of adding it in later <ref:2601.14860#pg3>.

Subrahmanyan: If this iterative mechanism holds up across different sources, then future observations will need to be very sensitive to those spectral features that depend on the disk inclination and spin, because those are what they found constrain most strongly here <ref:2601.14860#pg8>.

Kostas Kourmpetis, Shafqat Riaz, Honghui Liu, Temurbek Mirzaev, Cosimo Bambi, Debtroy Das, Jiachen Jiang, Kostas D. Kokkotas, Andrea Santangelo

Institut für Astronomie und Astrophysik Eberhard-Karls University Tübingen Department of Theoretical Astrophysics Institute for Astronomy and Astrophysics Center for Astronomy and Astrophysics Center for Field Theory and Particle Physics Department of Physics Fudan University School of Natural Sciences and Humanities New Uzbekistan University Department of Physics University of Warwick

astro-ph.HE

Submitted: 2026-01-21

Updated: 2026-10-21

Comments: 16 pages, 10 figures, 5 tables

Journal ref: A&A, 714, A428 (2026)

DOI: 10.1051/0004-6361/202659128

Code: https://github.com/reltrans/reltrans

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 74/100

The gist: The gist: A new table model for reflection spectroscopy that self-consistently accounts for returning radiation has been presented and applied to 4U 1630–47, revealing a rapidly spinning black hole

Key concepts

Returning Radiation
This refers to thermal radiation emitted by the accretion disk that travels back toward the central black hole. The new model self-consistently accounts for this returning radiation, which is thought to be a key component in illuminating the disk and producing observed reflection features in black-hole X-ray binaries.
Self-Irradiation of the Disk
In this study, they isolated the effect of returning radiation by disabling the corona. This allowed them to model how the thermal radiation from the disk reflects off itself, meaning it is being irradiated by its own emitted light rather than external coronal emission.
Black Hole Spin Parameter (a*)
This parameter describes how fast a black hole is spinning, ranging from 0 to 1. The analysis of 4U 1630-47 yielded a very high spin value, approximately 0.99, suggesting the black hole is spinning almost at its theoretical maximum limit.

Summary

The gist: A new table model for reflection spectroscopy that self-consistently accounts for returning radiation has been presented and applied to 4U 1630–47, revealing a rapidly spinning black hole with specific physical parameters.

Model Development

The authors developed a new table model for reflection spectroscopy that self-consistently accounts for the returning radiation <ref:2601.14860#pg4>. To isolate this effect, they adopted the standard disk-corona configuration but disable the corona, allowing the reflection spectrum to be produced solely by self-irradiation of the disk <ref:2601.14860#pg3>. The model utilizes an updated version of the ziji code and treats returning radiation as returning thermal radiation (purple arrows in Fig. 1(b); 0th iteration) and higher-order returning reflection radiation up to 3 iterations <ref:2601.14860#pg4>. The observed spectral flux is computed via ray-tracing, with the observed spectral flux given by the equation (B.2) <ref:2601.14860#pg5>.

Model Parameters and Grid

The XSPEC table model includes five parameters: electron density (ne), dimensionless spin parameter (a∗), black hole mass (MBH), mass accretion rate (M˙BH), and inclination angle (i) <ref:2601.14860#pg5>. The grid for the 9,100 points was chosen based on parameter values commonly reported in the literature for the X-ray binary 4U 1630–47 during the high-soft state <ref:2601.14860#pg6>. The electron density spans from 1018 to 1022 cm−3 <ref:2601.14860#pg6>, and the mass accretion rate covers a range of (0.14 − 2.1) 10 18 g/s <ref:2601.14860#pg6>. The black hole spin parameter (a∗) ranges from approximately 0.8 up to the Thorne limit (0.9982) <ref:2601.14860#pg6>.

Spectral Analysis of 4U 1630–47

The analysis was performed using NuS T AR data from an observation taken while the source was in the soft state <ref:2601.14860#pg4>. The best-fit model, zijiRetRad, yielded a rapidly spinning black hole (a∗ ∼ 0.99), a disk inclination of i ∼ 53◦, a mass accretion rate of M˙BH ∼ 15% M˙Edd, and an electron density of ne ∼ 1021 cm−3 <ref:2601.14860#pg5>. The inferred source distance from this analysis is D ∼ 8.2 kpc, slightly below the commonly adopted value of 10 kpc in the literature <ref:2601.14860#pg4>.

Comparison with Existing Models

The performance of zijiRetRad was compared against relxillNS, a widely used model for returning radiation <ref:2601.14860#pg10>. The study found that relxillNS requires a weak corona to adequately describe the data, in contrast to zijiRetRad <ref:2601.14860#pg10>. Furthermore, the authors demonstrated that the zijiRetRad reflection spectrum is significantly harder at high energies relative to relxillNS, and this behavior arises from the iterative reflection implemented in the ziji code <ref:2601.14860#pg8>.

Conclusion

The study concludes that returning radiation alone can account for the observed spectra and should therefore be considered an important component during the high-soft state <ref:2601.14860#pg9>. The best-fit indicates a high black hole spin (a∗ ∼ 0.99), a disk inclination of i ∼ 53◦, an increased electron density (ne ∼ 1021 cm−3) and a mass accretion rate of M˙BH ∼ 15% M˙Edd <ref:2601.14860#pg5>. The inferred source distance is lower than the values reported by Kalemci et al. (2025) and Kalemci et al. (2018) when modeling the dust-scattering halo, as well as the 10 kpc distance commonly adopted in the literature <ref:2601.14860#pg8>. The authors note that including a coronal component alongside zijiRetRad does not improve the quality of the fit <ref:2601.14860#pg10>. The paper suggests that a natural next step is to include a weak corona self-consistently, in order to assess the relative contributions of coronal emission and self-irradiation <ref:2601.14860#pg9>. The inferred source distance, D ∼ 8.2 kpc, is lower than the ∼ 10 kpc value commonly adopted in the literature <ref:2601.14860#pg8>.

How it works

The model construction involves computing the blackbody and reflection spectra over a grid of parameters <ref:2601.14860#pg6>. The code treats the disk as an infinitesimally thin and oriented perpendicular to the black hole’s spin axis <ref:2601.14860#pg4>. The reflection spectrum is computed using the "reflionx code9 (Ross & Fabian 2005) <ref:2601.14860#pg4>, which includes radiation from returning thermal radiation (purple arrows in Fig. 1(b); 0th iteration), higher-order returning reflection radiation up to 3 iterations" <ref:2601.14860#pg4>. The observed spectral flux is given by the equation (B.2) <ref:2601.14860#pg5>.

Key Findings

Key findings include:

a rapidly spinning black hole (a∗ ∼ 0.99), a disk inclination of i ∼ 53◦, a mass accretion rate of M˙BH ∼ 15% M˙Edd, and an electron density of ne ∼ 1021 cm−3 to reproduce the observed reflection features

<ref:2601.14860#pg5>.

The model also yields a source distance of D ∼ 8.2 kpc, slightly below the commonly adopted value of 10 kpc in the literature

<ref:2601.14860#pg4>.

Compared to the widely used relxillNS, our model naturally produces a harder high-energy reflection spectrum, fitting the data without invoking a Comptonized component

<ref:2601.14860#pg4>.

The hard excess arises naturally from the iterative reflection included in the model, which is absent in relxillNS

<ref:2601.14860#pg8>.

Appendix Analysis

The MCMC analysis revealed strong correlation between the two parameters of the xstar model; column density and ionization parameter <ref:2601.14860#pg11>. The authors found that the mass accretion rate and the source distance are also strongly correlated, since both primarily act as scaling factors for the observed flux <ref:2601.14860#pg11>. The analysis confirms that returning radiation alone can account for the observed spectra <ref:2601.14860#pg9>. The model's normalization plays the role of the geometric prefactor, as shown by the relation N = 0.229 · 10−20 M 2(M⊙) D 2(kpc) <ref:2601.14860#pg11>.

Summary of Results

The paper presents an XSPEC table model that isolates the reflection spectrum produced by returning radiation from the accretion disk, excluding any contribution from the corona <ref:2601.14860#pg4>. When applied to 4U 1630–47 in the soft state, this model provides a good fit to the data (χ 2/ν = 1.07) <ref:2601.14860#pg10>. The best-fit indicates a high black hole spin (a∗ ∼ 0.99), a disk inclination of i ∼ 53◦, an increased electron density (ne ∼ 1021 cm−3) and a mass accretion rate of M˙BH ∼ 15% M˙Edd <ref:2601.

Improvements for AI systems

  1. Improved spectral modeling capability: The AI system can now self-consistently model X-ray reflection spectra by incorporating returning radiation, which is described as self-irradiation of the disk in a pure-disk scenario. This allows the system to produce a harder high-energy reflection spectrum naturally, bypassing the need to invoke a separate Comptonized component when fitting data from sources in the high soft state.

  2. Enhanced parameter inference for black hole properties: The system can accurately infer key astrophysical parameters like a rapidly spinning black hole (a∗ ∼ 0.99), disk inclination, mass accretion rate, and electron density by fitting observed reflection features using the zijiRetRad model. This provides a more physically constrained set of results compared to models that neglect returning radiation or use simpler power-law assumptions.

  3. Comparative modeling against existing tools: The AI system can perform direct comparisons between its self-irradiation model (zijiRetRad) and widely used frameworks like relxillNS, identifying when one requires an additional weak coronal component versus the other producing a hard excess naturally. This enables researchers to choose the most appropriate physical description for specific observational data.

  4. Distance estimation precision: The system can estimate the source distance with a derived relation of N = 0.229 · 10−20 M squared (M⊙) D squared (kpc), offering an inferred distance of D ∼ 8.2 kpc for the target source while quantifying the uncertainty associated with this geometric prefactor.

  5. Automated identification of spectral signatures: The AI can identify specific absorption residuals, such as those corresponding to disk-wind features previously identified as Fe XXV and Fe XXVI absorption lines at 6.697 keV and 6.966 keV, allowing for the automated flagging of these high-ionization features in fit results.

Sources

Related papers