Pith. sign in

REVIEW 2 major objections 3 minor 3 references

The same NuSTAR/Suzaku reflection spectra of Cyg X-1 are fit equally well by a disk-like corona that changes the inferred black hole spin from ~0.99 to ~0.73–0.82 and the disk inclination from ~70° to ~30°.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 14:08 UTC pith:HDW64D2T

load-bearing objection A careful model-comparison study with a real systematic result, but the 'statistically similar' framing doesn't survive its own chi2 numbers. the 2 major comments →

arxiv 2512.21230 v2 pith:HDW64D2T submitted 2025-12-24 astro-ph.HE

Assessing systematic uncertainties from spectral re-analysis of Cyg X-1 with different coronal geometries

Abdurakhmon Nosirov , Jiachen Jiang , Cosimo Bambi , John A. Tomsick This is my paper
classification astro-ph.HE
keywords X-ray reflection spectroscopycoronal geometryblack hole spinCyg X-1accretion disk inclinationsystematic uncertaintiesdisk-like coronalamppost model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper re-analyzes the same X-ray spectra of the black hole binary Cyg X-1 under three different assumptions about the corona—the hot electron cloud that produces the X-rays—and asks how much the inferred black hole properties change. It finds that the geometry assumption is the largest source of systematic error: a horizontally extended, disk-shaped corona fits the data almost as well as the phenomenological broken power-law emissivity model and the compact lamppost geometry, but gives a disk inclination of about 30°, consistent with the binary orbit, instead of ~70°, and a spin of about 0.73–0.82 instead of ~0.99. The extended corona also removes the need for an unphysical negative outer emissivity index that previous fits required, by attributing the outer-disk illumination to a distant reflector. The paper concludes that a disk-like corona is physically well-motivated and statistically competitive, though spectral data alone cannot decide between the geometries.

Core claim

On its own terms, the paper's central claim is that the geometric modelling of the corona—not the spectral data themselves—is what pushes previous reflection fits of Cyg X-1 to extreme parameters. Replacing the broken power-law emissivity profile with a disk-like corona and allowing the disk density to be a free parameter produces a statistically comparable fit (Δχ²≈11) while lowering the inferred inclination from ~70° to ~30° and the spin from ~0.99 to a*=0.73^{+0.21}_{-0.07}. The disk density converges to n_e≈10^20 cm^-3, and the physically problematic negative outer emissivity index disappears. The paper presents this as evidence that systematic uncertainties from assumed geometry can exc

What carries the argument

The central object is the disk-like corona model (relconvdisk_nk): an infinitesimally thin, stationary disk of isotropic point sources at height h above the accretion disk, with sources placed at equal radial intervals so the intensity scales as 1/R, and photon paths ray-traced in Kerr spacetime. It computes the disk emissivity profile that is then convolved with a reflection spectrum, allowing the inner/outer coronal radii, height, spin, and inclination to be fit. Its key role is to replace the arbitrary broken power-law emissivity with a physically derived illumination pattern; in the best fit it flattens the emissivity and shifts the narrow Fe Kα line to a distant reflector, which is what

Load-bearing premise

The load-bearing assumption is that the real corona is close to an infinitesimally thin, stationary disk of isotropic point sources with no scattering or re-interaction of reflected photons inside the corona; if the corona has vertical thickness, bulk outflow, or reprocesses reflected photons, the inferred low inclination and density could be biased.

What would settle it

Fit the same NuSTAR/Suzaku data with a disk-like corona that includes a mildly relativistic outflow (β≈0.4) and Compton scattering of reflected photons in the corona; if the best-fit inclination rises above about 45° with a comparable fit statistic, the 30° inclination is an artifact of the static, optically thin geometry assumed here.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Spin and inclination measurements from reflection spectroscopy carry systematic uncertainties of order Δa* ≈ 0.2 and Δi ≈ 40° depending on the assumed corona shape, so single-geometry spin estimates should be quoted with model-uncertainty caveats.
  • A disk-like corona is consistent with the low binary inclination (~27°) of Cyg X-1, implying the inner disk need not be misaligned with the orbital plane.
  • The disk density in Cyg X-1 is near 10^20 cm^-3, so reflection models with a fixed density of 10^15 cm^-3 bias inferred iron abundance and ionization parameters.
  • Lamppost and disk-like geometries produce statistically indistinguishable fits, so polarimetric (spectro-polarimetric) data will be needed to discriminate between geometries.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the same geometry sensitivity applies to other black hole binaries and active galactic nuclei, many published spin constraints may be dominated by the assumed corona geometry rather than by data quality.
  • The paper's low-inclination solution implies that the high X-ray polarization degree of Cyg X-1 must be produced either by a relativistic outflow (β≥0.4) or by a different emitting geometry; measuring the polarization angle rotation across the Fe K line could test this.
  • Freeing the disk electron density in future fits of other sources may systematically reduce iron-abundance outliers, since density and abundance are partially degenerate.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper re-analyzes NuSTAR and Suzaku spectra of Cyg X-1 with three coronal-geometry assumptions: a phenomenological broken power-law emissivity profile (Model 0), an extended disk-like corona (Models 2–3), and a compact lamppost corona (Models 4–5). The authors report spin, inclination, and disk density for each geometry and argue that the disk-like corona gives a fit statistically similar to the others while yielding more physically reasonable parameters, notably an inclination near 30° and a disk density near 10^20 cm^-3. The main scientific claim is that coronal geometry is a major source of systematic uncertainty in reflection-based spin and inclination measurements.

Significance. If the statistical comparison were correct, this would be a useful addition to the ongoing discussion of coronal geometry in black hole X-ray binaries, especially given the recent IXPE polarimetric results motivating extended coronae. The paper carefully documents the fitting procedure, provides MCMC corner plots for key models, and reproduces earlier results. However, the central statistical claim is not supported by the reported chi-square values: the lamppost models are strongly preferred over the disk-like corona model by Δχ² ≈ 56 with equal degrees of freedom. The systematic-uncertainty point about geometry-dependent spin/inclination is still valuable, but it needs reframing: the disk-like corona is a worse-fitting alternative, not a statistically comparable one. The modeling limitations (stationary, thin, non-scattering corona) and the need for an outflow to reconcile polarization are already acknowledged, further weakening the claim that the disk-like corona is equally supported.

major comments (2)
  1. [§4.3.2, Tables 3–4; Abstract] The claim that the disk-like corona (Model 3, n_e free) is statistically indistinguishable from the lamppost models is contradicted by the reported fits. Model 3 has χ²/ν = 1296.5/1093; Model 4 has χ²/ν = 1240.3/1093 (Δχ² = 56.2, same dof); Model 5 has χ²/ν = 1249.3/1094 (Δχ² = 47.2, one fewer dof). These are large differences by standard likelihood-ratio or information-criterion standards, and they cannot be called 'statistically similar' or 'comparable.' The abstract's 'statistically similar' and Section 4.3.2's 'unable to statistically distinguish' should be corrected. The results should be presented as showing that the lamppost geometry is clearly preferred spectrally, while the disk-like corona gives a physically motivated but significantly worse fit.
  2. [§5, Table 4] The statement that the lamppost geometry gives 'a Δχ² improvement of 45 and one additional free parameter compared to Model 0' is inaccurate. Model 4 has χ²/ν = 1240.3/1093 and Model 0 has χ²/ν = 1285.3/1095, so Model 4 has two additional free parameters (ν decreases by 2), not one. Moreover, Model 4 includes a different reflection model with returning radiation and an ionization gradient, so the comparison is not a simple one-parameter geometry change. This point matters because the paper uses this comparison to rank the geometries; the number of added parameters should be stated correctly, or the non-nested nature of the models should be acknowledged explicitly.
minor comments (3)
  1. [Abstract and §4.2.3] The comparison between Model 3 (n_e free) and Model 0 is described as having 'relatively small' Δχ² ≈ 11, but Model 3 has two additional free parameters and a worse χ². A Δχ² of 11.2 for two extra parameters is not 'relatively small' in a standard likelihood-ratio sense (p ≈ 0.004). Please quantify or rephrase.
  2. [General] There are several typographical issues: 'exsitence' (Introduction), 'absorbtion' (Tables 2–4), 'above the the equatorial plane' (Section 3), and 'broken solid lines' in the Figure 3 caption. These should be corrected.
  3. [§4.3.2] The sentence 'Based on the spectral analysis alone, we are unable to statistically distinguish whether the lamppost geometry or the extended coronal geometry provides a superior description of the data' is at odds with the Δχ² values and should be revised.

Circularity Check

0 steps flagged

No significant circularity: the analysis is an empirical spectral-fit comparison; the reported parameters are fitted outputs, and the extended-corona self-citation is not used to derive the central result.

full rationale

The paper's central claim is a model-comparison result obtained by fitting NuSTAR and Suzaku spectra in XSPEC. The quoted spin, inclination, and disk density values are fitted parameters, not inputs used to construct the models, so no step reduces to its own inputs by construction. The disk-like corona model is adopted from S. Riaz et al. (2022), a self-citation with a coauthor overlap (C. Bambi), but the paper does not invoke that work as a uniqueness theorem or as a substitute for the fit; instead, the comparison is performed here and includes independent lamppost and broken power-law models. The model's simplifying assumptions are explicitly stated in Section 3 and Footnote 6, and the limitations are acknowledged in Section 5.1 and the Conclusion. I also note the internal statistical tension between the abstract's 'statistically similar' and the reported chi2 values (e.g., Model 3: 1296.5/1093 vs Model 4: 1240.3/1093), and the paper itself concedes 'the slight statistical inferiority' of the disk-like model. That is a consistency/correctness concern, not a circularity, and it does not change the finding that no fitted quantity is renamed as a prediction or derived from a self-citation chain.

Axiom & Free-Parameter Ledger

8 free parameters · 7 axioms · 0 invented entities

The paper is a model-comparison fitting exercise. Its key outputs (spin, inclination, density, emissivity indices) are free parameters fitted to data, while the extended corona geometry and reflection tables are imported from prior literature, including the authors' own relconvdisk_nk model. The main unmeasured assumptions are the corona's geometric idealization (stationary, thin, no scattering), the fixed Rout=24 rg, and the use of the binary inclination as the physical plausibility benchmark.

free parameters (8)
  • Black hole spin a* = 0.989, 0.986, 0.82, 0.73, 0.875, 0.92 (Models 0,1,2,3,4,5)
    Central inferred parameter; shifts by ~0.3 across geometries, driving the systematic-uncertainty claim.
  • Disk inclination i = 69.8, 69.1, 27.0, 30.6, 40.6, 38.9 deg
    Key parameter; the disk-like corona gives ~30 deg, matching binary orbit, while broken power-law gives ~70 deg.
  • Corona height h = >6.7, <4.5, <5.3, <2.1, 2.6 rg (Models 2,3,3,4,5)
    Free parameter in extended and lamppost models; often pegged at boundaries.
  • Inner corona radius R_in = <5, <10, <6 rg (Models 2,3,3)
    Free parameter for disk-like corona with lower limit 0.5 rg; only upper limits obtained.
  • Disk electron density log ne = 15 fixed; 20.2 (Model 3 free); 18.34, 18.4 (Models 4,5)
    Freeing density changes iron abundance and ionization; 20.2 is claimed as a measurement.
  • Emissivity indices qin, qout, rbreak = qin>9.3, qout=-6.96, rbreak=31 rg (Model 0)
    Phenomenological broken power-law parameters; negative qout is unphysical and motivates distant reflector/extended corona.
  • Iron abundance Fe/solar = 3.0, 3.7, >4.94, 2.0, 2.14, 2.6 (Models 0-5)
    Correlated with density and ionization; values shift substantially when ne is free.
  • Continuum/absorber parameters (Gamma, Ecut, kTe, kTin, NH, norms, distant-reflector xi/norm) = listed in Tables 2-4 per model
    Standard spectral fit parameters; they co-vary with geometry and affect the reported uncertainties.
axioms (7)
  • standard math Kerr spacetime and relativistic ray tracing for photon paths to the disk
    All relativistic blurring models (relconv, relxill, relconvdisk_nk) assume general relativity; used throughout Section 3 and 4.
  • domain assumption Rest-frame reflection tables (reflionx, reflionx_HD, xillver) correctly describe disk reflection physics
    The fitted spectra are convolved with precomputed reflection tables, whose atomic physics and densities are accepted as input (Section 4).
  • domain assumption Disk-like corona is infinitesimally thin, stationary, emits isotropically from point sources at equal radial intervals, and ignores coronal scattering/re-interaction
    Stated in Section 3 and footnote 6; this geometry is the defining feature of relconvdisk_nk and drives the low-inclination result.
  • ad hoc to paper Outer radius of the disk-like corona fixed at 24 rg
    Footnote 10 of Section 4.2.1 fixes Rout=24 rg due to model reliability limits; freeing it changes chi^2 by ~7 and may affect parameter constraints.
  • ad hoc to paper Disk inner radius at ISCO, outer radius 400 rg, redshift z=0
    Stated in Section 4.1.1 as modeling choices for the accretion disk geometry.
  • domain assumption Ionized wind absorption modeled with XSTAR grid at n=1e12 cm^-3, vturb=300 km/s
    Section 4.1.1 describes the XSTAR absorption grid; the assumed turbulence and density affect the continuum shape.
  • domain assumption Binary orbital inclination 27.1 deg is the appropriate prior for judging 'physically reasonable' inner disk inclination
    Section 5.1 uses the orbital inclination as the reference value; this is an external constraint, not derived from the X-ray spectra.

pith-pipeline@v1.3.0-alltime-deepseek · 20806 in / 12959 out tokens · 131058 ms · 2026-08-03T14:08:48.280428+00:00 · methodology

0 comments
read the original abstract

In this work, we carry out a new spectral reanalysis of NuSTAR and Suzaku observations of the disk reflection spectra in the stellar-mass black hole X-ray binary Cyg~X-1. We compare three types of models: a broken power-law disk emissivity profile with no assumption about the coronal shape used in the previous work of the same observations, a compact lamppost corona, and an extended disk-like corona motivated by recent X-ray polarization results. Our goal is to measure the systematic uncertainties caused by the assumed geometry, with a focus on key parameters such as the black hole spin and the inclination of the inner accretion disk. We find that the disk-like corona gives a fit that is statistically similar to the broken power-law and lamppost models, but it leads to more physically reasonable results, such as a lower inclination angle of about $30^{\circ}$. By using a variable disk density model, we measure the disk density to be $n_{\rm e}\approx10^{20}$\,cm$^{-3}$, which is similar to earlier results. While the extended corona model infers a wider allowed parameter space for black hole spin and the inner radius of the disk-shaped coronal region, this reflects the additional physical freedom of the model. Even so, the disk-like corona remains a strong and physically well-motivated candidate for explaining the X-ray emission from Cyg~X-1.

Figures

Figures reproduced from arXiv: 2512.21230 by Abdurakhmon Nosirov, Cosimo Bambi, Jiachen Jiang, John A. Tomsick.

Figure 1
Figure 1. Figure 1: Cartoon of the astrophysical system. The corona is described by an infinitesimal thin disk of inner radius of Rin, outer radius of Rout at height h above the the equatorial plane. The corona has a central hole with radius Rin. In the model, Rin is a free parameter with a lower limit of 0.5 rg, where rg = M is the grav￾itational radius of the black hole. This central gap is included for two computational re… view at source ↗
Figure 2
Figure 2. Figure 2: In the upper panel, the black solid line corresponds to the total model and other lines cor￾respond to its components. In the lower panel, ra￾tio plots of Model 0 are given and red, blue, green, orange, brown, and magenta correspond to NuS￾TAR’s FPMA, FPMB and Suzaku’s XIS0, XIS1, PIN, GSO data respectively. Analysis of the best-fit parameters in [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Error bars represent the ratio plots (data/model) for the models in the denominator for each case (see upper left corner of each plot). The black broken solid lines show the ratios of different models, as indicated in the upper right corner of each plot. This allows us to clearly assess how each model in the numerator improves the fit of the model in the denominator. Since significant changes are observed … view at source ↗
Figure 4
Figure 4. Figure 4: One-dimensional confidence intervals for the spin parameter across the different models considered in this work. The green dashed line in￾dicates the 90% confidence level (∆χ 2 = 2.706). statistical difference between these two models is relatively small, with ∆χ 2 ≈ 11. This dif￾ference could be further reduced through addi￾tional parameter exploration, as discussed ear￾lier, for example by freeing the ou… view at source ↗
Figure 6
Figure 6. Figure 6: Emissivity profiles for Models: 0, 3, 4 and 5 together with shaded 90% confidence inter￾vals. See the text for more information. additional distant reflector. This result is con￾sistent with the known binary orbital inclina￾tion (27.1 ◦ ±0.8 ◦ ; J. A. Orosz et al. 2011). Fur￾thermore, recent radio jet analysis of this source by S. Prabu et al. (2025) yields consistent re￾sults with very small misalignment … view at source ↗
Figure 7
Figure 7. Figure 7: Ratio plots (data/model) for all model combinations are provided for the 1-200 keV energy range, with corresponding titles. The colored error bars—red, blue, green, orange, brown, and magenta—represent data from NuSTAR’s FPMA, FPMB, and Suzaku’s XIS0, XIS1, PIN, and GSO detectors, respectively. A thick gray line marks the iron Kα line at 6.4 keV, with a 0.1 keV width [PITH_FULL_IMAGE:figures/full_fig_p017… view at source ↗
Figure 8
Figure 8. Figure 8: Corner plot of MCMC result of Model 0. angle of the accretion disk), h (height of the corona), kTe (electron temperature in the corona), log ne (electron density in logarithmic scale), and Γ (photon index). Acknowledgments – J.J. and A.N are supported by the Warwick-Fudan Joint Seed Grant. This work was supported by the National Natural Science Foundation of China (NSFC), Grant No. W2531002 [PITH_FULL_IMA… view at source ↗
Figure 9
Figure 9. Figure 9: Corner plot of MCMC result of Model 3. Here ξ and ξdist are the ionisation parameters of main and distant reflectors, respectively. REFERENCES Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 Baker, F. J. E., & Young, A. J. 2025, MNRAS, doi: 10.1093/mnras/staf1770 Bambi, C. 2024, arXiv… view at source ↗
Figure 10
Figure 10. Figure 10: Corner plot of MCMC result of Model 4. Here log ξ and log ξdist are ionisation parameters of main and distant reflectors in a logarithmic scale, respectively. Bambi, C., Brenneman, L. W., Dauser, T., et al. 2021, SSRv, 217, 65, doi: 10.1007/s11214-021-00841-8 Basak, R., Zdziarski, A. A., Parker, M., & Islam, N. 2017, MNRAS, 472, 4220, doi: 10.1093/mnras/stx2283 Bolton, C. T. 1972, nature, 235, 271, doi: 1… view at source ↗
Figure 11
Figure 11. Figure 11: Corner plot of MCMC result of Model 5. Here log ξ and log ξdist are ionisation parameters of main and distant reflectors in a logarithmic scale, respectively. Dauser, T., Garcia, J., Wilms, J., et al. 2013, MNRAS, 430, 1694, doi: 10.1093/mnras/sts710 Dauser, T., Wilms, J., Reynolds, C. S., & Brenneman, L. W. 2010, MNRAS, 409, 1534, doi: 10.1111/j.1365-2966.2010.17393.x Dove, J. B., Wilms, J., Maisack, M.,… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 1 linked inside Pith

  1. [1]

    Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes, 17 Baker, F. J. E., & Young, A. J. 2025, MNRAS, doi: 10.1093/mnras/staf1770 Bambi, C. 2024, arXiv e-prints, arXiv:2408.12262, doi: 10.48550/arXiv.2408.12262 20 0.645+0.023 0.023 0.528 0.53...

    Pith/arXiv arXiv doi:10.48550/arxiv.2408.12262 1996

  2. [4]

    Bambi, C., Brenneman, L

    Here logξand logξ dist are ionisation parameters of main and distant reflectors in a logarithmic scale, respectively. Bambi, C., Brenneman, L. W., Dauser, T., et al. 2021, SSRv, 217, 65, doi: 10.1007/s11214-021-00841-8 Basak, R., Zdziarski, A. A., Parker, M., & Islam, N. 2017, MNRAS, 472, 4220, doi: 10.1093/mnras/stx2283 Bolton, C. T. 1972, nature, 235, 2...

    arXiv 2021

  3. [5]

    Dauser, T., Garcia, J., Wilms, J., et al

    Here logξand logξ dist are ionisation parameters of main and distant reflectors in a logarithmic scale, respectively. Dauser, T., Garcia, J., Wilms, J., et al. 2013, MNRAS, 430, 1694, doi: 10.1093/mnras/sts710 Dauser, T., Wilms, J., Reynolds, C. S., & Brenneman, L. W. 2010, MNRAS, 409, 1534, doi: 10.1111/j.1365-2966.2010.17393.x Dove, J. B., Wilms, J., Ma...

    arXiv 2013