{"id":"21fefd59-4b6b-4e8b-9976-55e83b92000d","arxiv_id":"2507.21048","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Super-accreting X-ray binaries with jet and wind power above 1e39 erg/s can plausibly accelerate protons past several PeV and may feed the galactic cosmic-ray population above the knee.","lead":"This paper argues that X-ray binaries accreting near or above the Eddington limit can accelerate cosmic-ray protons to energies above several PeV. It proposes that a handful of such 'super-PeVatrons' could supply the galactic cosmic rays above the knee and predicts gamma-ray halos around them.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper explicitly assumes rather than demonstrates near-Hillas acceleration; the V4641 Sgr energy budget hints that the CR-production efficiency may be ~10^-2, which, if confirmed, would invalidate the population claim.","rationale":"The reader's weakest assumption identifies the same structural weakness: the Hillas limit is necessary but not sufficient, and the paper's central statement that these systems 'can accelerate protons to energies above several PeV' relies on the unproven assumption that acceleration operates near the theoretical maximum. I find no internal algebraic error in Eq. (1) or in the population counting; the concern is about the physical sufficiency of the argument. The paper is transparent about the assumption, which is why the result is a plausibility argument rather than a demonstration. The additional order-of-magnitude estimate from V4641 Sgr strengthens the concern: the published UHE bubble energy, when divided by the available kinetic energy over the inferred active period, suggests a CR-production efficiency well below the 10% assumed in the population calculation. This does not change the verdict because the reader already assigned CONDITIONAL with the same caveat; it confirms that the conditionality is justified. The proposed check — converting observed UHE gamma-ray data into a direct measurement of ϵ_CR — would settle whether the population claim is viable or should be downgraded.","tokens_in":18790,"tokens_out":10987,"duration_ms":139276,"concrete_test":"Perform a hadronic fit to the LHAASO/H.E.S.S. UHE gamma-ray spectra of V4641 Sgr and SS 433 that directly yields the total proton energy W_p(>E0). Combine this with the jet kinetic luminosity and the nebula age inferred from Eqs. (A1)–(A4) to compute ϵ_CR = W_p/(L_K t_active). If the inferred ϵ_CR is below ~0.1 by an order of magnitude or more, the population estimate in Sec. 3 and the associated claim that 'a few to tens of such systems could supply the galactic cosmic-ray flux above the knee' fail; the paper should then be weakened to an individual-source plausibility claim. If, instead, ϵ_CR is consistent with ~0.1, the population claim survives this check.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (1) is a necessary upper bound: E_max ∝ (σ L_K β)^{1/2}. The paper converts this bound into an achieved energy through two untested assumptions: (i) in Sec. 2.1, that each source 'operate[s] at or near its maximum capability,' and (ii) in Sec. 3, that the outflow kinetic power is comparable to the X-ray luminosity. Neither is demonstrated. If real acceleration reaches only a small fraction of the Hillas voltage — as is typical for many astrophysical accelerators — the fiducial source with L_K=10^39 erg/s, β=0.1, σ=0.1 no longer reaches multi-PeV energies, and the individual claim weakens. The population claim is even more sensitive. The estimate N≈2.5–21 in Sec. 3 assumes ϵ_CR≈0.1; if ϵ_CR≈10^-2, the required N becomes 250–2100, exceeding the adopted ULX population estimate N∼10. The paper cites V4641 Sgr as an example, but the quoted W_p≈10^50 erg (Alfaro et al. 2024), combined with L_K≈10^41 erg/s over an active time of ~10^4 yr, implies ϵ_CR≈10^-3–10^-2 — an order of magnitude below the assumed value. This is only an indicative estimate, but it shows the efficiency assumption is load-bearing and presently unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper argues that Galactic super-accreting X-ray binaries (including ULXs and microquasars in extreme states) are viable super-PeVatron candidates. It derives a Hillas-type maximum energy from the kinetic luminosity, magnetization, and velocity of trans-relativistic jets/winds (Eq. 1), compares several candidate source classes, estimates that a few to a few tens of such sources could supply the Galactic cosmic-ray flux above the knee (Section 3), and presents a forward model for UHE gamma-ray halos and nebulae (Section 4 and Appendix A). The central claim is that systems with kinetic luminosity ≥10^39 erg/s can accelerate protons to several PeV and above.","tokens_in":19110,"tokens_out":8689,"duration_ms":87060,"significance":"If the two load-bearing assumptions—near-Hillas acceleration efficiency and outflow kinetic power comparable to X-ray luminosity—are correct, the proposal is timely and physically interesting, and the paper provides a useful comparative survey of Galactic super-PeVatron candidates as well as a concrete, observationally testable halo-emission model. The halo flux formula (Eq. 4) and the nebular dynamics in Appendix A are constructive steps that can be matched against LHAASO, CTAO, SWGO, and future neutrino observations. However, the central claim is currently conditional on assumptions that are not demonstrated, and one internal example (V4641 Sgr) suggests an efficiency far below the adopted value.","major_comments":[{"comment":"Eq. (1) is an upper limit from the Hillas criterion, but the abstract and Section 5 convert it into the claim that super-accreting XRBs 'can accelerate protons to energies above several PeV' by taking the position that the accelerators operate at or near their maximum capability. This assumption is not tested; for a generic accelerator with efficiency well below the Hillas voltage, the fiducial source with L_K=10^39 erg/s, β=0.1, σ=0.1 would not reach multi-PeV energies. The paper should either provide evidence for near-Hillas acceleration in these systems or explicitly restrict the claim to an upper limit on source requirements.","section":"Section 2.1, Eq. (1) and following paragraph"},{"comment":"The assumption that outflow kinetic power is comparable to X-ray luminosity is load-bearing for both the individual Emax values in Table 1 and the population estimate. The paper does not justify this comparability; radiative efficiencies, beaming, and the distinction between transient and persistent states can make the X-ray luminosity either larger or smaller than the mechanical power. A sensitivity analysis varying L_K/L_X over a decade or two would clarify how robust the claimed multi-PeV capability is.","section":"Section 3, paragraph beginning 'For super-accreting XRBs identified in X-rays'"},{"comment":"Using the numbers quoted in the paper, the hadronic energy W_p≈10^50 erg (Alfaro et al. 2024) and the active time t≈10 kyr estimated from Eqs. (A1)-(A2) with L_j=10^41 erg/s imply a CR production efficiency ϵ_CR=W_p/(L_K t)≈3×10^-3, two orders of magnitude below the ϵ_CR≈0.1 adopted later for the population estimate. This is an internal consistency problem for the claim that a handful of such sources can supply the CR flux above the knee.","section":"Section 3, V4641 Sgr example"},{"comment":"The estimate N≈(2.5−21) ϵ_CR,-1^-1 ϵ_DC^-1 (τ_10/100)^-1 assumes ϵ_CR≈0.1. With ϵ_CR≈10^-2, the required number becomes 250–2100, and with ϵ_CR≈3×10^-3 it becomes roughly 10^3–10^4, far exceeding the adopted ULX population of N∼10. Because no independent constraint on ϵ_CR is provided, the population claim is not robust, and the summary statement that around ten sources may be sufficient should be explicitly qualified by the efficiency assumption.","section":"Section 3, final paragraph"}],"minor_comments":[{"comment":"The running header and title contain 'Super-PeV atron' due to an apparent line-break artifact; this should be corrected to 'Super-PeVatron'.","section":"Title / running header"},{"comment":"The coefficient in Eq. (2) is rounded from the exact inversion of Eq. (1) (which gives about 8.2×10^37 erg/s for Emax=10 PeV and the quoted fiducial parameters); the rounding is fine but a remark would avoid confusion.","section":"Section 2.1, Eq. (2)"},{"comment":"The flux expression omits CMB absorption, which is discussed later in the section; for sources at distances ≳8 kpc and photon energies above ≳PeV, the predicted halo flux should be attenuated before comparing with LHAASO or CTAO sensitivities.","section":"Section 4, Eq. (4)"},{"comment":"The paper states that the resulting photons have energies extending up to Eγ≈0.1Ep, but does not discuss the corresponding spectral shape or the influence of the diffusion index δ on the observable spectrum; a brief elaboration would make the predicted halo signal more directly testable.","section":"Section 4, text after Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is potentially suitable for ApJL, and the comparative survey plus the halo model are useful contributions. However, the central claim currently rests on two unverified assumptions, and the V4641 Sgr energy-budget argument in Section 3 directly points to an efficiency ϵ_CR≈10^-3, which would invalidate the population estimate as stated. I recommend major revision rather than rejection because the framework is transparent and the assumptions could be addressed with a more careful energy-budget analysis or by reframing the results as conditional upper limits."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper is a clean and transparent application of the Hillas criterion to super-accreting X-ray binaries. The opening claim that the authors 'demonstrate' proton acceleration above several PeV overstates what the paper actually delivers: they re-express the Hillas bound in terms of Poynting flux and kinetic luminosity, then explicitly assume the accelerators run near the theoretical maximum. That is a plausibility argument, not a demonstration. But the paper never hides this, which I respect.\n\nWhat is genuinely useful: Eq. (1) and Eq. (2) give a tidy form of the Hillas limit that will be handy for order-of-magnitude work. The comparative survey of YSCs, SNRs, PWNe, the Galactic Centre, and XRBs in Figure 1 and Table 1 is informative, and the individual source parameters are gathered from the literature in a compact way. The halo flux forecast and the population estimate are forward models with explicit parameters, so they are testable as stated. The authors also clearly distinguish their approach from Peretti et al. and Ohira, and they cite the prior work on microquasars and ULXs as PeVatrons fairly.\n\nWhere the soft spots are: everything hinges on two unverified assumptions. First, that acceleration operates at or near the Hillas maximum. Second, that the outflow kinetic power is comparable to the X-ray luminosity. The paper states both explicitly, so this is not a hidden flaw, but it means the 'super-PeVatron' conclusion is conditional. The V4641 Sgr energy budget is a useful sanity check: the quoted W_p ~ 10^50 erg and an active time of ~10^4 yr imply a CR production efficiency of roughly 10^-3 to 10^-2, an order of magnitude below the assumed 0.1. If that efficiency applies broadly, the required source count rises from N ~ 2.5-21 to N ~ 25-210 (the stress-test's 250-2100 is off by a factor of ten). That still exceeds the adopted ULX population of N ~ 10, though it might fit a starburst episode. So the population claim is sensitive, but not dead on arrival.\n\nThe halo flux estimates carry wide uncertainty because they depend on the diffusion coefficient, target gas density, source spectrum, and duty cycle. The authors acknowledge these choices, and the formulas are clearly presented, so a reader can adjust them. No algebraic errors jumped out at me, and the citation pattern is honest.\n\nWho should read this: anyone working on galactic cosmic-ray origins above the knee or on UHE gamma-ray sources will find the checklist and the source table worth having, even if they remain sceptical of the central claim. It deserves a real peer-review process, not a desk rejection. My suggestion to the editor: send it out, but ask the referee to push the authors to soften the 'demonstrate' language and to add an explicit discussion of how the inferred efficiency from V4641 Sgr affects the population claim.","headline":"Useful, honest Hillas-criterion case for super-accreting XRBs as super-PeVatrons, but the central claim is a plausibility argument, not a demonstration.","tokens_in":19683,"tokens_out":3024,"would_cite":true,"duration_ms":34943,"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":"Super-accreting X-ray binaries, whose jets and winds carry kinetic luminosity above $10^{39}\\ \\mathrm{erg\\,s^{-1}}$, can accelerate protons past several PeV, and a handful of them could supply the galactic cosmic-ray flux above the knee.","keywords":["cosmic rays","PeVatrons","super-accreting X-ray binaries","ultraluminous X-ray sources","Hillas criterion","ultra-high-energy gamma-ray halos","microquasars","galactic cosmic-ray knee"],"falsifier":"A concrete test is to measure the gamma-ray spectrum of a super-accreting X-ray binary halo, such as the V4641 Sgr bubble, above 1 PeV: a cutoff well below the predicted several-PeV maximum, combined with a confirmed kinetic luminosity around $10^{40}$ erg/s, would show the accelerators fall far short of the Hillas maximum and overturn the central claim.","tokens_in":18555,"feed_emoji":"🌌","tokens_out":14819,"duration_ms":142533,"temperature":0.7,"pith_summary":"The paper argues that super-accreting X-ray binaries — compact binary systems accreting near or above the Eddington limit, including ultraluminous X-ray sources — can accelerate protons to energies above several PeV, which would make them galactic 'super-PeVatrons.' The case rests on the kinetic power of their trans-relativistic jets and winds, which can exceed $10^{39}\\ \\mathrm{erg\\,s^{-1}}$, and on a Hillas-type maximum-energy formula that converts that power into a maximum particle rigidity. If the argument is right, a few to a few tens of such systems (roughly 3–20, depending on escape index, conversion efficiency, and duty cycle) could maintain the observed cosmic-ray flux at and above the knee, an energy range that supernova remnants struggle to reach. The paper also predicts that these systems should appear as ultra-high-energy gamma-ray halos and large-scale nebulae, giving observers concrete signatures to look for.","feed_headline":"X-ray binaries may push protons past the cosmic-ray knee","feed_subtitle":"A few super-accreting X-ray binaries with powerful outflows could supply the cosmic rays that cross the PeV knee.","key_machinery":"The machinery is the super-accreting X-ray binary outflow — a collimated jet ($\\tilde\\omega=1$) or quasi-spherical wind ($\\tilde\\omega=4$) — characterized by kinetic luminosity $L_K$, velocity $\\beta c$, and magnetization $\\sigma$. The load-bearing identity is the maximum-energy formula $E_{\\max} = 35Z\\sigma_{-1}^{1/2}(L_{K,39}\\beta)^{1/2}\\tilde\\omega^{-1/2}\\ \\mathrm{PeV}$, obtained by relating the outflow's Poynting flux to the electric potential across the acceleration zone; it turns observable outflow power into a maximum proton rigidity. The second mechanism is the population conversion: the ULX luminosity function $\\mathrm{d}N/\\mathrm{d}L\\propto L^{-1.6}$ with cutoff at $10^{41}\\ \\mathrm{erg\\,s^{-1}}$, together with population synthesis giving $N\\sim 10$ systems, yields the count $N\\approx (2.5\\text{–}21)\\epsilon_{CR,-1}^{-1}\\epsilon_{DC}^{-1}(\\tau_{10}/100)^{-1}$ needed to supply the knee flux. Nebula and halo dynamics (jet-cocoon and wind-termination-shock scalings) convert these predictions into observable sizes and fluxes.","core_discovery":"The central claim is that galactic super-accreting X-ray binaries, with outflow kinetic luminosity $L_K \\gtrsim 10^{39}\\ \\mathrm{erg\\,s^{-1}}$, are viable super-PeVatrons: protons can be accelerated above several PeV and, in flaring states, to rigidities of order $10^2$ PV. The authors derive this from the maximum-energy bound $E_{\\max} = 35\\,Z\\,\\sigma_{-1}^{1/2}(L_{K,39}\\beta)^{1/2}\\tilde\\omega^{-1/2}\\ \\mathrm{PeV}$, where $\\sigma$ is the outflow magnetization, $\\beta c$ its bulk velocity, and $\\tilde\\omega$ a geometry factor ($1$ for a jet, $4$ for a wind). Inverting this relation yields a minimum kinetic luminosity $L_K \\gtrsim 10^{38}(E_{\\max}/10\\ \\mathrm{PeV})^2\\tilde\\omega\\beta^{-1}\\sigma_{-1}^{-1}\\ \\mathrm{erg\\,s^{-1}}$, which super-Eddington outflows satisfy. Combining population synthesis ($N\\sim 10$ ULX-type systems in a Milky Way-like galaxy) with a luminosity function $\\mathrm{d}N/\\mathrm{d}L\\propto L^{-1.6}$, the paper estimates that about 3–20 sources, assuming $\\sim 10\\%$ conversion of outflow power into cosmic rays, can account for the galactic cosmic-ray flux above the knee. It further predicts PeV gamma-ray halos with fluxes near $10^{-13}\\,\\epsilon_{CR,-1}L_{K,39}(\\theta/1^\\circ)^2 n_t D_{30}^{-1}\\ \\mathrm{erg\\,cm^{-2}\\,s^{-1}}$ and comparable neutrino fluxes from the same halos.","pith_inferences":["A corollary the paper leaves implicit: the flaring-state power of sources like V4641 Sgr ($L_K\\sim 10^{40}$–$10^{41}\\ \\mathrm{erg\\,s^{-1}}$) is so far above the $10^{39}$ erg/s threshold that even a factor-of-ten shortfall from the theoretical maximum would still permit multi-PeV protons, so the near-maximum assumption could be relaxed for the most luminous flares.","The duty-cycle factor $\\epsilon_{DC}$ is likely the controlling unknown: if super-accreting phases are short and rare, the current X-ray census would undercount the relevant population, and nebular ages rather than flaring rates would be the better population measure.","The predicted MeV synchrotron cascade from pair production of absorbed PeV gamma rays offers a cheaper multi-wavelength test of hadronic halos than neutrinos alone, and could be searched for in hard X-ray data around known microquasars.","Sub-EeV cosmic-ray composition measurements could test the XRB origin directly: the super-solar $\\alpha$-process abundances measured in SS 433 and V4641 Sgr jets imply a distinctive heavy-element signature that would be hard to reproduce with supernova-remnant sources."],"forward_implications":["If the claim holds, a population of roughly 3–20 super-accreting X-ray binaries can explain the galactic cosmic-ray flux at and above the knee, with the exact number set by the escape-time index, the cosmic-ray conversion efficiency, and the duty cycle.","The outflows should inflate nebulae tens to hundreds of parsecs across and produce PeV gamma-ray halos with angular sizes $\\theta_h\\approx 0.7^\\circ D_{30}^{1/2}E^{\\delta/2}_{p,\\mathrm{PeV}}t_3^{1/2}(d/10\\ \\mathrm{kpc})^{-1}$, within reach of current and next-generation observatories.","PeV halos larger than the roughly 100 pc synchrotron and inverse-Compton cooling distance of PeV electrons would point to a hadronic, not leptonic, origin for the emission.","A comparable flux of neutrinos at $E_\\nu\\approx 0.05E_p$ should accompany the halos, offering an independent detection channel.","If the accelerated spectrum extends with index $s=2$ down to GeV energies, these systems could contribute to the whole galactic cosmic-ray population from GeV to tens of PeV, not just the above-knee component."],"supporting_citations":[{"why":"Supplies the maximum-energy (Hillas) bound from which Eq. (1) for $E_{\\max}$ is built.","marker":"A. M. Hillas 1984"},{"why":"Provides the classical-electrodynamics constraints that frame the maximum-energy estimate in terms of Poynting flux and potential across the acceleration zone.","marker":"F. A. Aharonian et al. 2002"},{"why":"Establishes that super-Eddington outflows carry kinetic luminosity near $L_{\\rm Edd}\\gtrsim 10^{39}$ erg/s and can inflate nebulae of tens to hundreds of parsecs.","marker":"A. King et al. 2023"},{"why":"Global radiation-MHD simulations of super-Eddington accretion disks that justify powerful trans-relativistic winds.","marker":"Y.-F. Jiang et al. 2014"},{"why":"Population synthesis giving about 10 ULX systems in a Milky Way-like galaxy, the basis of the source-count estimate.","marker":"G. Wiktorowicz et al. 2017"},{"why":"Population synthesis of black hole X-ray binaries supporting the same $N\\sim 10$ abundance.","marker":"Y. Shao & X.-D. Li 2020"},{"why":"Provides the observed high-mass X-ray binary luminosity function $\\mathrm{d}N/\\mathrm{d}L\\propto L^{-1.6}$ used to compute the average kinetic power.","marker":"S. Mineo et al. 2012"},{"why":"Ultra-high-energy gamma-ray bubble around V4641 Sgr used as the worked example; its hadronic energy budget motivates the multi-PeV requirement.","marker":"R. Alfaro et al. 2024"},{"why":"Ultra-high-energy gamma-ray detection from V4641 Sgr that anchors the example and motivates super-accreting XRBs as multi-PeV accelerators.","marker":"LHAASO Collaboration 2024b"},{"why":"TeV detection of SS 433 showing acceleration close to the maximum limit and supporting the spectral index $s=2$ adopted for the population estimate.","marker":"H.E.S.S. Collaboration et al. 2024"}],"fun_headline_variants":["Super-accreting X-ray binaries may be super-PeVatrons","PeV cosmic rays from super-Eddington X-ray binaries?","X-ray binary outflows could accelerate protons past the PeV knee","Galactic X-ray binaries as candidate ultra-high-energy cosmic-ray sources","Do super-accreting X-ray binaries power the cosmic-ray spectrum?"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument relies on super-accreting X-ray binaries accelerating particles at or near their theoretical maximum energy, and on outflow kinetic power being comparable to X-ray luminosity; if real acceleration is orders of magnitude less efficient, the multi-PeV conclusion fails even with correct source parameters.","fun_headline_variants_meta":{"raw":{"variants":["Super-accreting X-ray binaries may be super-PeVatrons","PeV cosmic rays from super-Eddington X-ray binaries?","X-ray binary outflows could accelerate protons past the PeV knee","Galactic X-ray binaries as candidate ultra-high-energy cosmic-ray sources","Do super-accreting X-ray binaries power the cosmic-ray spectrum?"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000392,"raw_usage":{"total_tokens":2127,"prompt_tokens":1077,"completion_tokens":1050,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":693,"completion_tokens_details":{"reasoning_tokens":958}},"tokens_in":693,"tokens_out":1050,"duration_ms":10900,"temperature":1.0,"reasoning_tokens":958,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T13:00:48.358024+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test is to measure the gamma-ray spectrum of a super-accreting X-ray binary halo, such as the V4641 Sgr bubble, above 1 PeV: a cutoff well below the predicted several-PeV maximum, combined with a confirmed kinetic luminosity around $10^{40}$ erg/s, would show the accelerators fall far short of the Hillas maximum and overturn the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the maximum-energy (Hillas) bound from which Eq. (1) for $E_{\\max}$ is built."},{"cited_title":"2023, Ultraluminous X-ray sources, NewAR, 96, 101672, doi: 10.1016/j.newar.2022.101672","cited_arxiv_id":null,"evidence_quote":"Establishes that super-Eddington outflows carry kinetic luminosity near $L_{\\rm Edd}\\gtrsim 10^{39}$ erg/s and can inflate nebulae of tens to hundreds of parsecs."}],"review_version":1}