REVIEW 4 major objections 5 minor 1 cited by
X-BPASS : Self-consistent modelling of stellar populations and their associated X-ray Binary emission in a binary stellar evolution framework
T0 review · 4 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Stellar evolution models now predict galaxies' X-ray binaries.
desk verdict Self-consistent BPASS+XRB spectral synthesis is a real step forward, but the SMC validation is partly circular and the model ships no code. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central machinery is the per-binary accretion disc model attached to the BPASS stellar evolution grid. For every interacting binary, the donor's mass-loss rate — via stellar wind, Roche-lobe overflow, or a Be-star decretion disc — sets the accretion rate; a uniform 20 percent duty cycle shortens each accretion episode; and the disc's temperature profile, computed in 100 annuli from standard thin-disc equations and truncated at an inner radius that is a free parameter, determines the emitted spectrum. Super-Eddington black hole accretion is treated with a slim-disc luminosity formula, and summing these spectra across the population in time bins matched to the stellar synthesis produces the combined stellar plus X-ray spectral energy distribution.
What would settle it
If a large, unbiased monitoring programme of individual X-ray binaries in the Small Magellanic Cloud showed that the fraction of time sources spend in outburst varies strongly with donor type or companion mass, the blanket 20 percent duty cycle would be invalid and the predicted number counts and flux normalization would shift. Alternatively, finding far more M51 sources above $10^{40}$ erg s$^{-1}$ than the Eddington-limited models produce would require raising the neutron star accretion cap or revising the assumed inner truncation radii.
Extended reading notes
Core claim
The central discovery is that the X-ray emission of a stellar population can be made a direct prediction of the binary stellar evolution models already used to synthesize its optical and ultraviolet light. By taking each compact remnant's mass accretion rate from the BPASS binary grid, imposing a 20 percent duty cycle, and computing disc spectra from a thin-disc temperature profile cut at physically motivated inner truncation radii — one white dwarf radius, eight Schwarzschild radii for neutron stars, three Schwarzschild radii for black holes — the models match the observed X-ray binary number evolution in the Small Magellanic Cloud and the X-ray flux evolution of M51. The same models predict that X-ray binaries contribute substantially to He II ionizing photons at ages beyond a few tens of millions of years, with the contribution comparable to the stellar one, while the hydrogen-ionizing photon budget remains dominated by stars.
Load-bearing premise
The uniform 20 percent duty cycle applied to every wind-fed and Roche-lobe overflow X-ray binary is the load-bearing assumption: the physics is poorly understood, no reliable method exists to assign a duty cycle from binary parameters, and this single factor directly scales both the number of X-ray binaries above the luminosity threshold and the total X-ray flux.
Editorial extensions
If this is right
- X-ray binary emission becomes a prediction of the same models that produce the stellar spectral energy distribution, so a galaxy's ultraviolet, optical, and X-ray light can be described by one consistent stellar population.
- For ages beyond roughly ten million years, X-ray binary photons are a non-negligible source of helium-ionizing radiation, with the X-ray contribution comparable to the stellar contribution at late ages.
- The predicted X-ray flux per unit star formation depends strongly on the duration of the star formation episode, so using X-ray luminosity as a star formation rate indicator requires a known star formation history.
- X-ray binaries alone are unlikely to account for the most extreme observed He II to H-beta ratios; the models indicate that additional sources or super-Eddington neutron star accretion would be needed to close the gap.
Reading between the lines
- The uniform 20 percent duty cycle is an acknowledged placeholder, and the same SMC and M51 data could be used to map how the duty cycle varies with donor type or companion mass rather than assuming a single value.
- Because the models already track compact remnant masses and binary orbital evolution, the same grid could connect X-ray binary luminosity with gravitational-wave merger rates predicted from the same stellar population.
- The predicted late-age supply of helium-ionizing photons suggests that quiescent or old stellar populations, even those too faint to detect in the ultraviolet, might show observable helium-recombination signatures; deep spectroscopy targeting such populations would be a direct test.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents X-BPASS, an extension of the BPASS binary population and spectral synthesis suite that computes accretion-disc X-ray emission for every accreting compact remnant in the BPASS binary grid and adds this emission to the stellar SEDs of the same population. The model applies standard thin-disc and slim-disc physics, Eddington limits, radiative efficiencies, and a duty cycle, with separate treatments for wind-fed, Roche-lobe overflow, Be, and super-Eddington black-hole systems. The authors validate the model against the SMC HMXB census of Antoniou et al. (2019), the M51 X-ray flux-age relation of Lehmer et al. (2017), and the GOODS relation of Gilbertson et al. (2022), and then use the models to predict the contribution of XRBs to He II and O VI ionizing photons. They conclude that XRBs can contribute to He II nebular emission without overproducing hydrogen-ionizing photons and that the same models reproduce the observed X-ray number and flux evolution of stellar populations.
Significance. If the framework is sound, X-BPASS would be a significant community asset: it provides a self-consistent prediction of stellar and XRB emission from the same evolutionary models, which is exactly what is needed for SED fitting, X-ray scaling relations, and nebular He II studies at high redshift. The construction is largely transparent and physical: the disc temperature profiles, Eddington limits, and spectral scaling follow standard accretion theory, and the inner-disc radii are referenced to observational constraints such as Aquila X-1. The paper makes falsifiable predictions, including the He II/H-beta versus L_X/SFR plane, O VI photon production, and the metallicity dependence of the X-ray luminosity function. The honest and detailed caveats section is a strength. The main weaknesses are the partly circular SMC validation from the calibrated Be normalization, the unconstrained global duty cycle, and the acknowledged but unquantified excess in the M51/GOODS flux comparisons at intermediate ages; these are fixable with additional tests and reframed claims.
major comments (4)
- [§2.4.1 and §3.1, Eq. (5), Fig. 1] The Be XRB accretion-rate normalization in Eq. (5) is explicitly described as 'selected to broadly reproduce the observed number of Be XRBs in the Small Magellanic Cloud as outlined in Antoniou et al. (2019)', and §3.1 then compares the full model prediction to that same Antoniou et al. (2019) HMXB census in Fig. 1. Because Be systems make up a substantial fraction of the SMC HMXB population, the agreement in Fig. 1 is partly imposed by construction and cannot be presented as an independent validation of the Be treatment. The abstract's 'reproduces ... validating our models' and conclusion 4 overstate the evidential weight of the SMC comparison. Please quantify the fraction of predicted SMC HMXBs that are Be systems, show the comparison with the Be component removed, and either calibrate Eq. (5) on an independent sample or explicitly re-label the SMC comparison as a consistency check of a calibrated model.
- [§2.3, §§3.1-3.2] The uniform 20 per cent duty cycle is an unconstrained global input that multiplies both the number of systems above the luminosity threshold and the integrated X-ray flux, and the paper states that there is no reliable method to assign a duty cycle from binary parameters. No sensitivity study is presented, so it is not possible to tell whether the SMC counts and M51/GOODS flux agreement are robust to the duty cycle, which Sidoli & Paizis (2018) find varies between roughly 10 and 55 per cent by HMXB type. In addition, the implementation is not equivalent across accretion channels: wind-fed and RLOF systems are active for only 20 per cent of each timestep, while Be systems have their peak luminosity multiplied by 0.2 and are assumed active over the full timestep, which is not the same for snapshot number counts. Please provide duty-cycle sensitivity tests and justify the Be prescription.
- [§5.4 and Figs. 7-9] Section 5.4 candidly states that the model X-ray flux shows an excess in the region 8.0 < log(age/yr) < 9.0, attributed to overproduction of low-mass black holes, yet §3.2 presents the M51 comparison in Figs. 7 and 8 as validation, and the abstract claims the models reproduce the observed X-ray flux evolution of M51. The same excess is visible in the GOODS comparison in Fig. 9. The validation claim needs to be quantified: state the size of the excess relative to the observational uncertainties, or restrict the claim to the age ranges where the model is actually within uncertainties.
- [§2.5, Table 1] The inner disc radius is described as 'a free parameter which encompasses much of the uncertainty in our disc model', and the disc luminosity scales linearly with 1/R_inner. The central validation claim is that observationally motivated values of R_inner reproduce the SMC and M51 data, but no sensitivity analysis is shown for the chosen values (R_inner = 1 R_wd, 8 R_s, 3 R_s). Please show how the predicted counts and fluxes vary when R_inner is moved within the plausible range from the literature (for example tau = 5-10 for neutron stars), to demonstrate that the agreement does not simply follow from the particular choice of this parameter.
minor comments (5)
- [Eq. (1)] The typesetting of Eq. (1) is garbled in the manuscript ('𝐿disc= G𝑀CR⁄𝑀CR 2𝑅inner'); please ensure it reads L_disc = G M_CR Mdot_CR / (2 R_inner).
- [Throughout] There are several typographical errors: 'fuducial' appears in the captions of Figs. 4-6 and in Section 3.3, and 'viscious decretion disc' appears in Section 2.4.1; these should be corrected to 'fiducial' and 'viscous'.
- [§3.2-3.3] The text discusses the M51 comparison in the 2-7 keV band and the GOODS comparison in the 2-10 keV band, but it does not explicitly state whether the luminosity functions in Fig. 10 and the flux comparisons in Figs. 7-9 are corrected for this bandpass difference; please clarify the energy bands used in each figure.
- [Data Availability] The Data Availability statement lists the BPASS models but not the X-BPASS post-processing code; if the code is to be released, please provide a link or state its availability.
- [§4.2] The conversion factors used to estimate line fluxes from photon production rates are taken from 'table A1 of Eldridge & Stanway (2022)', but this table is not reproduced or summarised; a short description of the conversion would make the calculation more transparent.
Circularity Check
The SMC number-count validation is partly circular because the Be XRB accretion-rate normalization in Eq. 5 is tuned to the same Antoniou et al. (2019) census that is later used as an independent validation in Section 3.1.
-
fitted input called prediction
[Section 2.4.1 (Eq. 5); Section 3.1 (Figure 1)]
"For Be XRB systems we use a modified version of the accretion rate formulae of Liu et al. (2023) ... replace it with a unique mass accretion rate for each system equal to the wind-loss rate of the donor reduced by a factor of 1e-5. The reduction factor is selected to broadly reproduce the observed number of Be XRBs in the Small Magellanic Cloud as outlined in Antoniou et al. (2019). ... We validate our model construction through comparison with observational data. Antoniou et al. (2019) compiled a census of HMXBs, in the Small Magellanic Cloud."
Equation 5 contains a freely chosen factor of 1e-5, and the paper states explicitly that this factor is selected to reproduce the observed number of Be XRBs in the SMC from Antoniou et al. (2019). Section 3.1 then uses the same Antoniou et al. census as the validation target for the predicted HMXB number evolution in Figure 1. Because the fitted Be normalization directly contributes to the early-time HMXB counts, the agreement shown in Figure 1 is not an independent confirmation of the Be treatment; part of the match is imposed by construction. The paper does not remove or separately report the calibrated Be contribution when assessing the SMC agreement, so the SMC validation is partially circular.
full rationale
The paper's load-bearing validation is partly circular at one specific point: the Be XRB accretion-rate factor in Eq. 5 is fitted to the Antoniou et al. (2019) SMC census, and that same census is then presented as an independent validation of the predicted SMC HMXB numbers in Section 3.1 and Figure 1. This is a genuine fitted-input-called-prediction step. However, the circularity is partial rather than total: the M51 comparison with Lehmer et al. (2017) and the GOODS comparison with Gilbertson et al. (2022) are external datasets not used in the calibration, and the accretion-disc spectral calculation, inner-radius choices, and Eddington-limited luminosity treatment are based on independent literature prescriptions rather than on the validation data. The duty-cycle assumption is uncertain but is not circular, since it is taken from independent observational estimates rather than fitted to the validation targets. The paper also explicitly acknowledges the Be normalization choice, making the step transparent rather than hidden. No self-citation chain is load-bearing, and no uniqueness theorem or ansatz-smuggling pattern is present. Overall, the central claim that the models are 'validated' by the SMC comparison is overstated because one component of that comparison was tuned to it, but the independent M51 and GOODS checks mean the derivation is not wholly circular.
Assumptions & free parameters
free parameters (4)
- R_inner (inner disc truncation radius) per remnant type =
BH: 3 R_s; NS: 8 R_s; WD: 1 R_wd
- Be XRB accretion reduction factor =
1e-5
- Duty cycle =
0.2 (20 percent) for wind-fed, RLOF, and Be systems
- Radiative efficiency for WD and NS accretors =
0.1
assumptions (5)
- domain assumption BPASS v2.2.1 binary stellar evolution models provide accurate remnant masses, mass-transfer histories, and binary survival probabilities.
- domain assumption Accretion onto compact remnants forms a steady-state, geometrically thin blackbody disc whose inner radius is R_inner and whose emission is entirely thermal.
- domain assumption WD and NS accretion is capped at the Eddington rate in the fiducial models; BH accretion is uncapped.
- ad hoc to paper A uniform 20 percent duty cycle applies to all wind-fed, RLOF, and Be XRBs, with the active interval placed randomly within each timestep.
- ad hoc to paper Be donors can be identified from stellar structure (no He core, temperature 10-30 kK, H fraction > 0.4, mass > 6 Msun, initial mass < 30 Msun) rather than from rotation.
Cite this review
Pith. "Pith review of X-BPASS : Self-consistent modelling of stellar populations and their associated X-ray Binary emission in a binary stellar evolution framework." pith.science (2026). https://pith.science/paper/IO4UCQML
@misc{pith2026250818628,
author = {Pith},
title = {Pith review of: X-BPASS : Self-consistent modelling of stellar populations and their associated X-ray Binary emission in a binary stellar evolution framework},
year = {2026},
howpublished = {\url{https://pith.science/paper/IO4UCQML}},
note = {Machine review of arXiv:2508.18628}
}
read the original abstract
X-ray binaries play a significant role in the thermal and ionization history of galaxies. Their X-ray luminosity can shed light on galactic star formation rates and histories. Compact objects are also crucial in the evolution of gravitational wave progenitors. Here we present the results from our work to extend the binary population and spectral synthesis (BPASS) code suite to incorporate X-ray emission onto compact remnants in binary systems. We self-consistently model the accretion disc for each interacting binary system in a grid of stellar evolution models and then combine these to obtain the total X-ray spectra for stellar populations over a range of ages and metallicities. Crucially, these are estimated using the same stellar models as those used for modelling the stellar spectral energy distribution. We utilise first principle equations to calculate the X-ray binary (XRB) evolution, luminosity and spectral energy densities of individual accreting compact objects. Population synthesis using observationally motivated values for R_inner (the accretion disc inner truncation radius) reproduces the observed X-ray number evolution in the Small Magellanic Cloud and the inferred X-ray flux evolution for M51, validating our models. Using these models, we explore the implications of a self-consistent stellar and XRB emission population synthesis for ionizing photon production, the XRB dependence on metallicity and, for XRBs as a potential source of nebular He II emission seen in the spectra of high redshift galaxies. We conclude that XRBs contribute towards powering nebular He II emission without causing significant overestimates of hydrogen ionization.
Figures
Figures from the paper (13 more)
Forward citations
Cited by 1 Pith paper
-
Constructing X-ray Spectral Models of Galaxies: Varying Contributions from X-ray Binary Populations with Host Galaxy Properties
The authors present and validate a stochastic sampling method that generates X-ray spectral models for X-ray binary populations from galaxy properties and delivers a main-sequence spectral library as a function of ste...
Reference graph
Works this paper leans on
-
[1]
Antoniou V., et al., 2019, ApJ, 887, 20 AvdanS.,AkyuzA.,AcarS.,AkkayaOralhanI.,AllakS.,AksakerN.,2023, MNRAS, 519, 4826 Belczynski K., Kalogera V., Rasio F. A., Taam R. E., Zezas A., Bulik T., Maccarone T. J., Ivanova N., 2008, ApJS, 174, 223 Belczynski K., Wiktorowicz G., Fryer C. L., Holz D. E., Kalogera V., 2012, ApJ, 757, 91 BergD.A.,ChisholmJ.,ErbD.K...
work page 2019
-
[2]
Douna, V. M. Pellizza, L. J. Mirabel, I. F. Pedrosa, S. E. 2015, A&A, 579, A44 Dray L. M., 2006, Monthly Notices of the Royal Astronomical Society, 370, 2079 Dziełak M. A., Zdziarski A. A., Szanecki M., De Marco B., Niedźwiecki A., MarkowitzA.,2019,MonthlyNoticesoftheRoyalAstronomicalSociety, 485, 3845 Eldridge J. J., Stanway E. R., 2022, ARA&A, 60, 455 E...
work page Pith review arXiv 2015
-
[3]
Decretion disc evolution and neutron star accretion in short-period eccentric Be/X-ray binaries
Ranalli, P. Comastri, A. Setti, G. 2003, A&A, 399, 39 Rast R. G., Jones C. E., Suffak M. W., Carciofi A. C., 2025, Decretion disc evolution and neutron star accretion in short-period eccentric Be/X- ray binaries (arXiv:2502.04705), https://arxiv.org/abs/2502. 04705 Rocha K. A., et al., 2024, The Astrophysical Journal, 971, 133 Rosdahl J., et al., 2018, MN...
work page Pith review arXiv 2015
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.