{"id":"7332b97a-6e12-4e69-b17e-2ac25fd64f93","arxiv_id":"2508.18628","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A self-consistent stellar population and X-ray binary spectral synthesis model, X-BPASS, reproduces observed X-ray binary counts and flux trends and predicts that X-ray binaries contribute to nebular He II emission without overproducing hydrogen ionizing photons.","lead":"X-BPASS is an upgraded computer model that adds X-ray light from binary star systems to the BPASS stellar population code, using the same stellar models that produce the starlight. It roughly reproduces observed X-ray binary numbers in the Small Magellanic Cloud and X-ray flux over time in M51, and suggests these systems contribute to helium ionization without overproducing hydrogen ionization.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SMC validation is partly circular: the Be XRB normalization in Eq. 5 is tuned to the same Antoniou et al. (2019) data later used as validation in Fig. 1.","rationale":"The central claim is that X-BPASS reproduces the observed SMC X-ray number evolution and M51 flux evolution, validating the model. The reader identified the uniform 20 percent duty cycle as the weakest assumption. That is a reasonable concern because the duty cycle scales all predicted numbers and fluxes. However, the more load-bearing issue is the explicit calibration of the Be accretion-rate normalization (Eq. 5) to the SMC Be XRB population. The same SMC dataset is then used in Fig. 1 to demonstrate agreement, so the SMC leg of the validation is partly circular. This does not mean the model is wrong, but it weakens the 'validation' claim in the abstract and Section 3.1. The M51 and GOODS comparisons provide independent support, although they share the Be normalization and other tuned parameters, so they are not fully orthogonal. The uniform duty cycle is also uncertain, but it is at least observationally motivated and not fitted to the SMC validation sample; the Be normalization is explicitly selected to match that sample. A simple test removing the Be component from the SMC comparison would settle whether the agreement is predictive. This concern supports the reader's CONDITIONAL verdict rather than changing it, because the issue is addressable with an independent check or a reframed validation claim.","tokens_in":26805,"tokens_out":7588,"duration_ms":73731,"concrete_test":"Recompute Fig. 1 with the Be XRB component removed and compare only wind-fed/RLOF systems to the Antoniou et al. (2019) HMXB counts; if the agreement degrades substantially, the claimed SMC number evolution is largely an artifact of the Be normalization fitted to the same data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 2.4.1 (Eq. 5) sets the Be XRB accretion rate to the donor wind-loss rate times a factor of 1e-5, stating explicitly that this normalization 'is selected to broadly reproduce the observed number of Be XRBs in the Small Magellanic Cloud' from Antoniou et al. (2019). Section 3.1 then uses that same Antoniou et al. census to validate the model's HMXB number evolution (Fig. 1). Since the SMC is known to host a substantial Be XRB population, the agreement in Fig. 1 is not an independent confirmation of the Be treatment: it is partly imposed by construction. The abstract's claim that the models are 'validated' by the SMC comparison therefore overstates the evidential weight of that comparison. The M51 flux comparison is more independent, but the Be normalization also contributes to the total X-ray flux at early times, so the model's overall calibration is entangled with the validation dataset. A direct test that removes the calibrated Be component would establish how much of the SMC agreement is genuinely predictive rather than a consequence of the fitting procedure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":27043,"tokens_out":10618,"duration_ms":99827,"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":[{"comment":"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.","section":"§2.4.1 and §3.1, Eq. (5), Fig. 1"},{"comment":"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.","section":"§2.3, §§3.1-3.2"},{"comment":"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.","section":"§5.4 and Figs. 7-9"},{"comment":"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.","section":"§2.5, Table 1"}],"minor_comments":[{"comment":"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).","section":"Eq. (1)"},{"comment":"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'.","section":"Throughout"},{"comment":"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.","section":"§3.2-3.3"},{"comment":"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.","section":"Data Availability"},{"comment":"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.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely to be a useful contribution after revision. The main risk is overclaiming validation: the SMC comparison is partly circular because the Be normalization in Eq. (5) is calibrated to the same Antoniou et al. (2019) data, and the duty cycle and inner-disc radii are unconstrained global inputs that directly affect the headline comparisons. I would ask the authors to add sensitivity tests, quantify the Be fraction in the SMC comparison, and temper the 'validating' language. No concerns about citation patterns or scope; the manuscript fits well within the journal's remit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a useful and mostly sound extension of BPASS that couples binary stellar population synthesis to accretion disc SEDs for XRBs. The new piece is real: for the first time you get a self-consistent stellar + XRB spectral synthesis from the same stellar evolution grid, with predictions for Q(HI), Q(HeII), OVI, and HeII/Hbeta. The comparisons to M51, GOODS, and the Mineo XLF are in the right ballpark, and the authors are appropriately cautious about the HeII conclusions.\n\nThe soft spots are not fatal but they matter. Most important: the Be XRB normalization in Eq. 5 is tuned to reproduce the SMC number counts from Antoniou et al. (2019), and the same dataset is then used as validation in Fig. 1. That's a circularity, and the abstract's 'validating our models' overstates the evidential weight. The M51 and GOODS comparisons are more independent and help, but the Be component also contributes at early times, so the entanglement is real. A clean test would remove the calibrated Be contribution and see what remains.\n\nThere are also free parameters: R_inner per remnant type, the duty cycle (uniform 20%), the radiative efficiency, plus the Eddington cap on NSs. The authors disclose these, and some are anchored to observations, but the combination of several adjustable parameters means the good agreement is not a sharp test of the underlying physics. The lack of any public code or data products is a genuine limitation: the paper describes a new tool but doesn't ship it. And the pass on Lecroq et al. (2024) is a missed opportunity – they acknowledge the overlap but provide no quantitative comparison, which weakens the claim of novelty.\n\nOverall the central framework holds up. The accretion disc physics is standard, the BPASS machinery is mature, and the paper is honest about its caveats and about the fact that XRBs alone don't explain the most extreme HeII. I'd send this to a serious referee; with a code release, an SMC check that excludes the fitted Be component, and a direct comparison to Lecroq et al., it would be a solid paper. Without those, it's still citable but needs the caveats taken seriously.","headline":"Self-consistent BPASS+XRB spectral synthesis is a real step forward, but the SMC validation is partly circular and the model ships no code.","tokens_in":27613,"tokens_out":1998,"would_cite":true,"duration_ms":18182,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Stellar evolution models now predict galaxies' X-ray binaries.","keywords":["X-ray binaries","binary stellar evolution","population synthesis","accretion discs","BPASS","He II ionizing photons","Small Magellanic Cloud","X-ray luminosity function"],"falsifier":"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.","tokens_in":26504,"feed_emoji":"🔭","tokens_out":11809,"duration_ms":94677,"temperature":0.7,"pith_summary":"This paper extends the BPASS binary stellar evolution and spectral synthesis suite so that the X-ray light from accreting compact remnants is computed for each interacting binary using the same stellar models that produce the population's starlight. It claims that with observationally motivated values for the inner radius of the accretion disc, the resulting X-ray binary populations reproduce the number of X-ray binaries per stellar age bin seen in the Small Magellanic Cloud and the age evolution of X-ray flux in M51. At the same time, the models predict that X-ray binaries supply a significant share of helium-ionizing photons at late ages without overproducing hydrogen-ionizing photons, helping to explain the excess He II nebular emission in high-redshift galaxies.","feed_headline":"Stellar evolution models now predict galaxies' X-ray binaries","feed_subtitle":"BPASS models match SMC counts and M51 flux, linking X-ray binaries to helium-ionizing photons.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the SMC census of high-mass X-ray binaries against which the model's X-ray binary number evolution is validated.","marker":"Antoniou et al. (2019)"},{"why":"Provides the M51 X-ray flux versus stellar age data that the model's integrated flux evolution is compared with.","marker":"Lehmer et al. (2017)"},{"why":"Describes the BPASS binary stellar evolution and population synthesis grid on which the accretion post-processing is built.","marker":"Eldridge et al. (2017)"},{"why":"Gives the Be X-ray binary accretion prescription that the paper modifies to assign a mass accretion rate to each Be system.","marker":"Liu et al. (2023)"},{"why":"Supplies the radiative efficiency formula for black holes and the Eddington luminosity and mass accretion rate prescriptions.","marker":"Podsiadlowski et al. (2003)"},{"why":"Provides the slim disc luminosity formula used for super-Eddington black hole accretion.","marker":"King et al. (2023)"},{"why":"Supplies the spectral hardening factor applied to disc annuli with temperatures above roughly 10,000 K.","marker":"Shimura & Takahara (1995)"},{"why":"Provides the classical wind-fed accretion rate formula used for non-Be wind-fed systems.","marker":"Frank et al. (2002)"},{"why":"Defines the observed He II/H-beta versus L_X/SFR diagnostic space that the models' ionization predictions are tested against.","marker":"Senchyna et al. (2020a)"},{"why":"Supplies the observed duty cycle estimates for different high-mass X-ray binary classes that motivate the paper's adopted average duty cycle.","marker":"Sidoli & Paizis (2018)"}],"fun_headline_variants":["X-ray binaries now predicted by stellar evolution models","BPASS links X-ray binaries to helium ionization","Stellar models tie X-ray binaries to galaxy spectra","Self-consistent X-ray binaries from binary evolution","X-ray binary emission emerges from BPASS models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["X-ray binaries now predicted by stellar evolution models","BPASS links X-ray binaries to helium ionization","Stellar models tie X-ray binaries to galaxy spectra","Self-consistent X-ray binaries from binary evolution","X-ray binary emission emerges from BPASS models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00063,"raw_usage":{"total_tokens":2939,"prompt_tokens":1005,"completion_tokens":1934,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":621,"completion_tokens_details":{"reasoning_tokens":1862}},"tokens_in":621,"tokens_out":1934,"duration_ms":12482,"temperature":1.0,"reasoning_tokens":1862,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:54:56.310421+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"A., Taam R","cited_arxiv_id":null,"evidence_quote":"Supplies the SMC census of high-mass X-ray binaries against which the model's X-ray binary number evolution is validated."}],"review_version":2}