Flux variability of the "10 keV feature" of 4U 0115+63
Pith reviewed 2026-05-10 16:51 UTC · model grok-4.3
The pith
The 10 keV feature in 4U 0115+63 varies independently of the continuum and cyclotron lines as flux changes during QPOs.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Flux-resolved spectroscopy shows that the 10 keV feature exhibits an anticorrelation between its flux relative to the continuum and the continuum flux level itself, while its centroid energy remains constant within a given observation. No analogous flux-dependent deviation appears at the energies of the harmonic cyclotron resonant scattering features. These observations strengthen earlier indications that the 10 keV feature is formed by a mechanism distinct from the continuum emission, even though both components originate in the same physical environment near the neutron star.
What carries the argument
The ratio of the flux in an added Gaussian 10 keV emission component to the flux of the cutoff power-law continuum, extracted from spectra binned by count rate during the ~500 s QPOs.
If this is right
- The 10 keV feature cannot be dismissed as an artifact of continuum modeling or as a simple blend with the cyclotron harmonics.
- Its formation must respond to changes in local accretion conditions on the 500 s timescale without altering the cyclotron line energies or depths.
- Spectral models of accreting pulsars should retain the 10 keV component as an independent term rather than absorbing it into the continuum parameters.
- Quasiperiodic flux swings offer a practical way to separate spectral components without assuming that global parameters remain fixed across an entire observation.
Where Pith is reading between the lines
- If the anticorrelation persists across multiple outbursts, it could constrain the radial or vertical location where the 10 keV photons are produced relative to the cyclotron scattering region.
- Similar flux-resolved tests on other sources that display a 10 keV feature would show whether the independence is generic or specific to 4U 0115+63.
- The constancy of the feature's energy while its strength changes argues against explanations that tie it directly to a fixed cyclotron energy or to a simple temperature shift in the emitting plasma.
Load-bearing premise
The chosen spectral model isolates the 10 keV feature without leaving residual curvature from the continuum or cyclotron lines that could mimic the observed changes in flux ratio.
What would settle it
Re-fitting the same flux-binned spectra with a different continuum shape or without an explicit 10 keV Gaussian component removes the apparent anticorrelation between the 8-12 keV excess and overall flux.
Figures
read the original abstract
X-ray spectra of accretion-powered X-ray pulsars can often be described using a power-law continuum with a high-energy cutoff, which might be further modified by additional spectral components. The Be X-ray binary system 4U 0115+63 is well known for having one of the highest numbers of detected harmonics of its cyclotron resonant scattering features (CRSFs), a pronounced spectral component known as the ''10 keV feature,'' and quasiperiodic oscillations (QPOs) with a period of about 500 s during outbursts. The changes in count rate by a factor of two during the approximately 500 s QPOs allow us to probe the variation in the spectral components with flux. We study the ''10 keV feature'' in emission, aiming to disentangle it from the broadband continuum and CRSFs and investigate its origin. We focus on the flux-dependent behavior of the CRSF and its harmonics, and particularly the contribution of the ''10 keV feature,'' as seen in the flux-resolved analysis of two NuSTAR observations of the 2015 outburst. Comparing the flux-resolved spectra of a given observation with the respective total dataset revealed a distinct change in overall spectral shape at the position of the ''10 keV feature'' but no comparable deviation at the energies of the harmonic CRSFs. The change associated with the ''10 keV feature'' does not seem to involve its centroid energy, which remains constant within a given observation. We find indications for an anticorrelation between the continuum flux and the ratio of the ''10 keV feature'' flux to the continuum flux within each observation. The analysis strengthens previous claims that the ''10 keV feature'' shows some independence from the remaining features. This result supports the interpretation that the ''10 keV feature'' has a different formation mechanism than the continuum emission, although its origin lies within the same physical environment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes two NuSTAR observations of the 2015 outburst of the Be/X-ray binary 4U 0115+63, exploiting count-rate variations by a factor of ~2 during ~500 s QPOs to perform flux-resolved spectroscopy. It reports a distinct change in spectral shape at the position of the '10 keV feature' (modeled as a separate component) when comparing flux-binned spectra to the time-averaged spectrum, with no analogous change at the energies of the CRSFs and their harmonics. The analysis finds indications of an anticorrelation between the continuum flux and the ratio of the 10 keV feature flux to the continuum flux within each observation, while the feature's centroid energy remains constant. This is interpreted as evidence that the 10 keV feature has a different formation mechanism from the continuum and CRSFs, although it originates in the same physical environment.
Significance. If the reported anticorrelation is robust, the work strengthens prior suggestions of independence for the 10 keV feature and supplies a concrete, observationally testable distinction between its behavior and that of the CRSFs. The direct use of flux segments from public NuSTAR data, without reliance on fitted parameters that could introduce circularity, is a methodological strength that makes the result falsifiable and potentially reproducible by other groups.
major comments (2)
- [Flux-resolved spectral analysis] The central claim that the 10 keV feature exhibits flux-independent behavior (distinct from the continuum and CRSFs) rests on the assumption that the adopted spectral model (power-law plus high-energy cutoff, CRSFs, and isolated 10 keV component) fully accounts for the broadband shape. The manuscript provides no quantitative fit statistics, residual plots, or alternative-model comparisons for the flux-resolved bins, leaving open the possibility that unmodeled continuum curvature or flux-dependent absorption produces the reported shape change at 10 keV.
- [Results and interpretation] The anticorrelation between continuum flux and the 10 keV feature-to-continuum flux ratio is described only as 'indications' without reported correlation coefficients, uncertainties, or Monte Carlo tests of parameter covariances. This weakens the load-bearing contrast drawn with the CRSFs, for which no equivalent flux-ratio variations are quantified.
minor comments (2)
- The abstract and title consistently place '10 keV feature' in quotes; the main text should adopt the same notation for clarity when referring to the modeled component.
- [Introduction] The description of the QPO period as 'about 500 s' could be made more precise by citing the exact value or range used to define the flux segments.
Simulated Author's Rebuttal
We thank the referee for their careful and constructive review of our manuscript. Their comments highlight important aspects of the analysis that we will address to improve clarity and robustness. We respond to each major comment below.
read point-by-point responses
-
Referee: [Flux-resolved spectral analysis] The central claim that the 10 keV feature exhibits flux-independent behavior (distinct from the continuum and CRSFs) rests on the assumption that the adopted spectral model (power-law plus high-energy cutoff, CRSFs, and isolated 10 keV component) fully accounts for the broadband shape. The manuscript provides no quantitative fit statistics, residual plots, or alternative-model comparisons for the flux-resolved bins, leaving open the possibility that unmodeled continuum curvature or flux-dependent absorption produces the reported shape change at 10 keV.
Authors: We agree that quantitative fit statistics and residual plots for the individual flux bins would strengthen the presentation. In the revised manuscript we will add a table listing the reduced chi-squared and degrees of freedom for each flux-resolved spectrum, together with residual plots (data-to-model ratios) in an appendix. The model itself was selected because it provides an acceptable description of the time-averaged spectrum, which is consistent with the extensive prior literature on 4U 0115+63; we will add a short paragraph explaining why flux-dependent absorption is disfavored by the absence of corresponding changes at the CRSF energies. revision: yes
-
Referee: [Results and interpretation] The anticorrelation between continuum flux and the 10 keV feature-to-continuum flux ratio is described only as 'indications' without reported correlation coefficients, uncertainties, or Monte Carlo tests of parameter covariances. This weakens the load-bearing contrast drawn with the CRSFs, for which no equivalent flux-ratio variations are quantified.
Authors: We accept that the qualitative phrasing 'indications' is insufficient. In the revision we will report the Spearman rank correlation coefficient and its uncertainty for the 10 keV feature-to-continuum ratio versus continuum flux in each observation. We will also describe a Monte Carlo test in which the flux ratios are randomly reassigned to assess the probability of obtaining the observed anticorrelation by chance. For the CRSFs we will explicitly compute and report the corresponding correlation coefficients, which remain consistent with zero, thereby quantifying the contrast. revision: yes
Circularity Check
No significant circularity: empirical flux-resolved spectral analysis
full rationale
The paper reports direct comparisons of flux-binned NuSTAR spectra for 4U 0115+63, identifying shape changes at the 10 keV feature and an anticorrelation between continuum flux and the feature-to-continuum flux ratio. These results follow from standard data reduction and model fitting applied to observed count-rate segments; no equation or parameter is defined in terms of the target anticorrelation, no prediction is statistically forced by a prior fit to the same data, and no load-bearing step reduces to a self-citation or imported uniqueness theorem. The central claim remains an empirical observation rather than a tautological renaming or self-referential derivation.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
G., Meszaros, P., & Bussard, R
Alexander, S. G., Meszaros, P., & Bussard, R. W. 1989, ApJ, 342, 928
work page 1989
-
[2]
Arnaud, K., Gordon, C., & Dorman, B. 2022, Xspec – An X-Ray Spec- tral Fitting Package, https://heasarc.gsfc.nasa.gov/xanadu/xspec/ manual/XspecManual.html, accessed on 2022-09-21
work page 2022
- [3]
-
[4]
Bachetti, M., Harrison, F. A., Cook, R., et al. 2015, ApJ, 800, 109
work page 2015
- [5]
-
[6]
Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Demleitner, M., & Andrae, R. 2021, AJ, 161, 147
work page 2021
-
[7]
Barthelmy, S. D., Barbier, L. M., Cummings, J. R., et al. 2005, Space Sci. Rev., 120, 143
work page 2005
-
[8]
Becker, P. A. & Wolff, M. T. 2007, ApJ, 654, 435
work page 2007
-
[9]
Becker, P. A. & Wolff, M. T. 2022, ApJ, 939, 67
work page 2022
-
[10]
Berger, K. 2022, MSc thesis, Friedrich-Alexander-Universität Erlangen-Nürnberg Bissinger, né Kühnel, M., Kreykenbohm, I., Ferrigno, C., et al. 2020, A&A, 634, A99
work page 2022
-
[11]
Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2005, Space Sci. Rev., 120, 165
work page 2005
-
[12]
Coburn, W., Heindl, W. A., Rothschild, R. E., et al. 2002, ApJ, 580, 394
work page 2002
-
[13]
W., Li, F., Mayer, W., & Rappaport, S
Cominsky, L., Clark, G. W., Li, F., Mayer, W., & Rappaport, S. 1978, Nature, 273, 367
work page 1978
-
[14]
M., Grinberg, V ., Fürst, F., et al
Diez, C. M., Grinberg, V ., Fürst, F., et al. 2022, A&A, 660, A19
work page 2022
-
[15]
Z., Wang, W., Zhang, P., et al
Ding, Y . Z., Wang, W., Zhang, P., et al. 2021, MNRAS, 503, 6045
work page 2021
-
[16]
Dugair, M. R., Jaisawal, G. K., Naik, S., & Jaaffrey, S. N. A. 2013, MNRAS, 434, 2458
work page 2013
-
[17]
Farinelli, R., Ferrigno, C., Bozzo, E., & Becker, P. A. 2016, A&A, 591, A29
work page 2016
-
[18]
A., Segreto, A., Mineo, T., & Santangelo, A
Ferrigno, C., Becker, P. A., Segreto, A., Mineo, T., & Santangelo, A. 2009, A&A, 498, 825
work page 2009
-
[19]
2020, A&A, 642, A73 Fürst, F., Kreykenbohm, I., Pottschmidt, K., et al
Ferrigno, C., Bozzo, E., & Romano, P. 2020, A&A, 642, A73 Fürst, F., Kreykenbohm, I., Pottschmidt, K., et al. 2010, A&A, 519, A37 Fürst, F., Pottschmidt, K., Wilms, J., et al. 2014, ApJ, 780, 133 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1
work page 2020
- [20]
- [21]
-
[22]
Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103
work page 2013
-
[23]
Heindl, W. A., Coburn, W., Gruber, D. E., et al. 2000, BAAS, 32, 1230
work page 2000
-
[24]
Heindl, W. A., Coburn, W., Gruber, D. E., et al. 1999, ApJ, 521, L49
work page 1999
-
[25]
Houck, J. C. & Denicola, L. A. 2000, in Astronomical Society of the Pacific Con- ference Series, V ol. 216, Astronomical Data Analysis Software and Systems IX, ed. N. Manset, C. Veillet, & D. Crabtree, 591
work page 2000
-
[26]
C., Bhattacharya, D., & Seetha, S
Iyer, N., Mukherjee, D., Dewangan, G. C., Bhattacharya, D., & Seetha, S. 2015, MNRAS, 454, 741
work page 2015
- [27]
-
[28]
Kortright, J. B. & Thompson, A. C. 2000, X-Ray Data Booklet – Section 1.2 X- ray emission energies, https://xdb.lbl.gov/Section1/Sec_1-2.html, accessed on 2021-06-16
work page 2000
- [29]
- [30]
- [31]
- [32]
- [33]
- [34]
-
[35]
Madsen, K. K., Beardmore, A., & Page, K. 2020, X-calibration of NuSTAR and Swift,url: https://iachec.org/wp-content/presentations/2020/ Xcal_swift_nustar.pdf
work page 2020
-
[36]
K., Burwitz, V ., Forster, K., et al
Madsen, K. K., Burwitz, V ., Forster, K., et al. 2021, arXiv:2111.01613
- [37]
-
[38]
Manikantan, H., Paul, B., Sharma, R., Pradhan, P., & Rana, V . 2024, MNRAS, 531, 530
work page 2024
-
[39]
Martin, R. G. & Charles, P. A. 2024, MNRAS, 528, L59
work page 2024
-
[40]
Martin, R. G., Nixon, C., Armitage, P. J., Lubow, S. H., & Price, D. J. 2014, ApJ, 790, L34 Mészáros, P. 1992, High-energy radiation from magnetized neutron stars. (Chicago: Univ. Chicago Press)
work page 2014
-
[41]
Miller, G., Wasserman, I., & Salpeter, E. E. 1989, ApJ, 346, 405
work page 1989
-
[42]
1991, ApJ, 383, 784 Müller, S., Ferrigno, C., Kühnel, M., et al
Miyamoto, S., Kimura, K., Kitamoto, S., Dotani, T., & Ebisawa, K. 1991, ApJ, 383, 784 Müller, S., Ferrigno, C., Kühnel, M., et al. 2013, A&A, 551, A6
work page 1991
- [43]
-
[44]
Nakajima, M., Mihara, T., Makishima, K., & Niko, H. 2006, ApJ, 646, 1125
work page 2006
- [45]
- [46]
-
[47]
Nelson, R. W., Wang, J. C. L., Salpeter, E. E., & Wasserman, I. 1995, ApJ, 438, L99
work page 1995
-
[48]
T., Hayasaki, K., & Moritani, Y
Okazaki, A. T., Hayasaki, K., & Moritani, Y . 2013, PASJ, 65, 41
work page 2013
-
[49]
Pottschmidt, K., Kreykenbohm, I., Wilms, J., et al. 2005, ApJ, 634, L97
work page 2005
-
[50]
Rappaport, S., Clark, G. W., Cominsky, L., Joss, P. C., & Li, F. 1978, ApJ, 224, L1 Rouco Escorial, A., Bak Nielsen, A. S., Wijnands, R., et al. 2017, MNRAS, 472, 1802
work page 1978
- [51]
- [52]
-
[53]
1999, ApJ, 523, L85 Schönherr, G., Wilms, J., Kretschmar, P., et al
Santangelo, A., Segreto, A., Giarrusso, S., et al. 1999, ApJ, 523, L85 Schönherr, G., Wilms, J., Kretschmar, P., et al. 2007, A&A, 472, 353
work page 1999
-
[54]
W., Schönherr, G., Falkner, S., et al
Schwarm, F. W., Schönherr, G., Falkner, S., et al. 2017, A&A, 597, A3
work page 2017
-
[55]
2023, PhD thesis, Friedrich-Alexander-Universität Erlangen- Nürnberg
Sokolova-Lapa, E. 2023, PhD thesis, Friedrich-Alexander-Universität Erlangen- Nürnberg
work page 2023
-
[56]
Sokolova-Lapa, E., Stierhof, J., Dauser, T., & Wilms, J. 2023, A&A, 674, L2
work page 2023
-
[57]
Staubert, R., Trümper, J., Kendziorra, E., et al. 2019, A&A, 622, A61
work page 2019
-
[58]
Stierhof, J. J. R., Sokolova-Lapa, E., Berger, K., et al. 2025, A&A, 698, A308
work page 2025
-
[59]
Thalhammer, P., Bissinger, M., Ballhausen, R., et al. 2021, A&A, 656, A105
work page 2021
-
[60]
Troja, E. 2020, The Neil Gehrels Swift Observatory, Technical Handbook, Version 17.0, https://swift.gsfc.nasa.gov/proposals/tech_appd/ swiftta_v17.pdf, accessed on 2023-08-17
work page 2020
-
[61]
Tsygankov, S. S., Lutovinov, A. A., Doroshenko, V ., et al. 2016, A&A, 593, A16
work page 2016
-
[62]
J., Roche, P., Negueruela, I., et al
Unger, S. J., Roche, P., Negueruela, I., et al. 1998, A&A, 336, 960
work page 1998
-
[63]
Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487
work page 1996
-
[64]
Vybornov, V ., Klochkov, D., Gornostaev, M., et al. 2017, A&A, 601, A126
work page 2017
- [65]
- [66]
-
[67]
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 Article number, page 12 of 15 Katrin Berger et al.: Flux variability of the “10 keV feature” of 4U 0115+63 Appendix A: Additional figures (N1) In the main body of the paper we showed the results from the flux-resolved analysis of N2, see Fig. 7. The residuals of the analysis of observation N1 are show...
work page 2000
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.