{"id":"d1e1e824-727d-4b11-a19b-d24eeaf30ed5","arxiv_id":"2512.18786","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The PeV gamma-ray signal from Cygnus X-3 can be roughly reproduced by a hadronic jet model in which protons scatter on gas and stellar UV photons, with enhanced wind absorption steepening the sub-PeV spectrum.","lead":"Astronomers test whether newly released LHAASO gamma-ray observations of Cygnus X-3, up to 4 PeV, match predictions that the binary system accelerates cosmic rays in its jet. They find rough agreement if the gamma rays come from protons colliding with gas and starlight, and if wind absorption shapes the spectrum.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Agreement with LHAASO hinges on an unvalidated ~2000-fold wind-photon density enhancement; without it the model's own spectrum rises too slowly.","rationale":"The reader's weakest assumption is precisely the enhanced UV photon density in the line-driven wind. As a stress-tester, I find this is indeed the single most load-bearing element of the paper's central claim. The spectral shape of the LHAASO data in Fig. 1 is the primary quantitative evidence for the model, and the model's own hadronic component from gas interactions is flat, while the pγ component rises too slowly below the threshold. The only mechanism introduced to force agreement is the enhanced absorption, based on a rough random-walk estimate whose physical application to continuum photons is questionable. The paper itself flags the limitation: 'is too slow to reproduce well the measurements.' This is not an external disagreement with consensus; it is an internally acknowledged reliance on an order-of-magnitude estimate. The model does have independent support: the threshold peak near PeV and the orbital modulation are predictive and not tuned to fit the data, so the paper merits a conditional rather than outright rejection. My recommendation is therefore to leave the reader's verdict unchanged: the paper is plausible but should be conditional on validating the wind-photon enhancement with a detailed wind model or independent measurement. A concrete Monte Carlo radiative-transfer calculation for the WR wind would settle whether a ~2000-fold enhancement actually occurs; until then, the central quantitative claim is not firmly established.","tokens_in":5234,"tokens_out":5607,"duration_ms":67177,"concrete_test":"Run a Monte Carlo radiative-transfer simulation of a Wolf–Rayet wind with Cygnus X-3 parameters (v_inf, M_dot, R_s, L_UV) as in Vilhu et al. 2021, and compute the angle-averaged UV photon density within 3 R_s as a function of radius and photon energy. If the enhancement over the free-streaming dilution factor (R_s/d)^2 is below ~100, or if the enhanced energy density is dominated by narrow line photons that contribute negligibly to γγ opacity, then τ_γγ at 1 PeV will be <1 and the 'absorbed' curve in Fig. 1 overestimates attenuation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's quantitative agreement with the LHAASO sub-PeV drop relies entirely on an ad hoc enhancement of the stellar UV photon density by a factor of order 100–1000 (adopted as ~2000) within ~3 R_s of the Wolf–Rayet star. The enhancement is justified by a random-walk estimate N ~ (v_inf/Δv)^2, with Δv ~ 1 km/s, giving sqrt(N) ~ v_inf/Δv ~ 2000. This argument is not derived from a detailed wind model and is applied to the continuum UV photon field, even though line-driven wind trapping concerns photons in spectral lines. The paper explicitly acknowledges in the Numerical results section that without this enhanced absorption the predicted flux 'is too slow to reproduce well the measurements.' Thus the claimed match to LHAASO data is not a robust prediction of the hadronic-jet model; it is a post hoc rescaling tuned to make the spectrum fit. If the true density enhancement is, say, an order of magnitude smaller, or if the enhanced photons occupy line frequencies with negligible total energy, then τ_γγ at 1 PeV would drop below ~0.1 and the fast sub-PeV rise would remain unexplained.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper is a short addendum to the authors' earlier model of Cygnus X-3 as a PeV cosmic-ray accelerator. It confronts the KL25 hadronic-jet model with the LHAASO 2025 measurement of orbitally modulated photons up to 4 PeV. The model has two photon components: hadronic interactions on stellar-wind gas, which follow the injected proton spectrum, and photo-hadronic interactions on stellar UV photons, which produce a component rising up to about 1 PeV. To reproduce the sharp sub-PeV rise in the data, the authors add gamma-gamma absorption on an enhanced UV photon density inside about 3 stellar radii; the enhancement is estimated as ~2000 from a random-walk argument for line-driven winds. They also argue that X-ray target photons are negligible and that a leptonic origin cannot reach PeV energies, and they predict an orbitally modulated neutrino flux. The agreement with LHAASO is presented by visual comparison in Fig. 1.","tokens_in":5550,"tokens_out":8615,"duration_ms":89626,"significance":"Should the interpretation survive scrutiny, it would be an important step: it identifies a microquasar as a PeV hadronic accelerator and connects the LHAASO spectrum to wind physics and a testable neutrino orbital modulation. The paper has several strengths: it is based on a published Monte Carlo framework, separates the gas and UV target components, uses the phase information of the LHAASO data qualitatively, and makes falsifiable predictions for neutrino time templates and for a declining modulation below 100 TeV. Its main weakness is that the quantitative match to the sub-PeV drop rests entirely on an absorption enhancement factor that is calibrated to the data and derived with a simplified, possibly misapplied, line-scattering argument. As presented, the result is a proof-of-concept rather than a validated model; a revision that turns the enhancement into a derived or scanned parameter with a real fit statistic is necessary.","major_comments":[{"comment":"The central agreement with LHAASO is obtained by tuning the absorption. The text states that the UV-only spectrum is 'too slow to reproduce well the measurements' and that matching the fast rise requires absorption 'increased by a factor of order 100-1000 compared to the standard treatment'; the adopted value, ~2000 from v_inf/Delta_v, is then used to display a curve that 'agrees better with the data.' Because this factor is not measured independently and no sensitivity scan is shown, the procedure is circular: the model is adjusted to the same spectral feature that it is then used to explain. Please provide a scan in the enhancement factor (including 0, 100, 1000, 2000) and a quantitative fit statistic, or derive the factor from a concrete wind model.","section":"Numerical results; Fig. 1"},{"comment":"The random-walk estimate N ~ (v_inf/Delta_v)^2 and the density enhancement sqrt(N) ~ 2000 describe photons trapped in spectral lines of a line-driven wind. It is not established that the continuum differential density n_gamma(epsilon) that enters Eq. (1) and the gamma-gamma opacity is enhanced by the same factor at the photon energies relevant for p-gamma and gamma-gamma interactions. The paper applies this single factor to the full UV target field. Please specify the frequency-dependent photon density after line transport, and identify which target energies are actually enhanced; otherwise the claim tau_gamma_gamma ~ 1 at 1 PeV is not supported.","section":"Numerical results; random-walk argument"},{"comment":"The comparison to LHAASO data is qualitative. No statistical test, residuals, or quantified uncertainty is presented, and the predicted orbital modulation in Fig. 3 is not compared with the measured modulation amplitude or phase. For a paper whose main conclusion is agreement with data, a simple chi-squared or likelihood over the binned spectrum, plus a comparison of the phase curve, should be added. This is particularly important because the data set has only a handful of energy bins.","section":"Numerical results; Fig. 1 and Fig. 3"},{"comment":"The assumption that X-ray photons are a negligible target is based on the absence of orbital modulation and on the IXPE picture. However, X-rays from the inner envelope could still contribute a non-modulated opacity in the acceleration region, and the previous KL25 model included them. Please quantify the X-ray photon density at the assumed acceleration distance and show that the p-gamma rate on X-rays is small compared to the UV rate for the adopted B=2000 G and L0=7e11 cm. As written, setting the X-ray contribution to zero is another parameter choice that selectively favors the UV interpretation.","section":"Reaction rates; neglect of X-ray targets"}],"minor_comments":[{"comment":"The notation for the angle is inconsistent: Eq. (1) uses mu_0 while Fig. 2 and the text use mu. Define the angle once and keep the notation fixed.","section":"Eq. (1) and Fig. 2"},{"comment":"The phrase 'the slope in an E2dN/dE plot' should read E^2 dN/dE.","section":"Introduction"},{"comment":"Typos: 'relaying on strong magnetisation' should be 'relying'; 'photo-hadronic interaction of stellar UV photons' should be 'interactions with stellar UV photons'.","section":"Throughout"},{"comment":"The reference to Barrios-Jimenez (2025) is a conference talk; if a published version exists, it should be cited instead.","section":"References"},{"comment":"The axis label 'R/s-1' should be typeset as 'R [s^-1]'.","section":"Fig. 2"},{"comment":"Clarify that the LHAASO data are the flaring-state, de-absorbed spectrum while the model curve is phase-averaged. If the comparison is phase-averaged, state that explicitly in the caption and in Section 'Numerical results'.","section":"Fig. 1 caption/main text"}],"recommendation":"major_revision","confidential_remarks":"For the editor: The paper is a short, phenomenon-driven addendum with a clear and testable hypothesis. The main concern is the post hoc absorption factor. I think the appropriate decision is major revision rather than rejection, because the authors can address the circularity by reformulating the absorption factor as a scanned free parameter, by adding a line-frequency-resolved estimate, and by reporting a fit statistic. If the authors are unwilling or unable to do so, the manuscript's central quantitative claim would remain unsupported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the Kachelrieß & Lammert addendum. The short version: it is a legitimate extension of their earlier KL25 model, and it does what an addendum should—confronts new data and proposes a concrete mechanism. But the agreement with LHAASO's sub-PeV drop is not a robust prediction: it rests on an order-of-magnitude estimate of photon-density enhancement in a line-driven wind that is effectively tuned to make the spectrum fall where it does. The authors say so themselves: without the enhanced absorption, the predicted flux 'is too slow to reproduce well the measurements.'\n\nWhat's genuinely good: the two-component interpretation (gas + UV) is physically reasonable; the UV-threshold rise around PeV is a natural prediction of photo-hadronic interactions; and the orbital modulation of the p-gamma channel is a real, testable signature, including for neutrinos. I also appreciate that they are explicit about what changed relative to KL25—lower B field, neglect of X-ray targets—and give a physical reason (orbital phase and IXPE) for dropping X-ray photons.\n\nThe soft spot is the absorption factor. The random-walk estimate N ~ (v_inf/Δv)^2 leads to sqrt(N) ~ 2000 when Δv ~ 1 km/s, but that argument is about line scattering, and it is being applied to the continuum UV field. That is a category slip that needs scrutiny. The factor of 100–1000 needed by the data is consistent with the estimate, but consistency is not derivation. If the real enhancement is an order of magnitude smaller, τ_γγ at 1 PeV drops below ~0.1 and the sharp rise stays unexplained. The paper would be much stronger with a wind model or at least a sensitivity scan. The comparison to data is also visual—no statistical tests—and the B field and source size remain free parameters, though the choices are not crazy.\n\nThat said, this is not a crackpot paper. It is a short, clear, honest contribution from people who know the system. The central idea is plausible and the limitations are not hidden. What it needs is a referee who can ask for a proper treatment of the wind or, failing that, an explicit acknowledgment that the absorption enhancement is a phenomenological parameter.\n\nBottom line: send it to peer review. It deserves referee time, but accept only if the absorption mechanism is either better justified or reframed as a fit parameter with a clear uncertainty.","headline":"A short, honest addendum that ties a hadronic-jet model to the new LHAASO Cygnus X-3 data, but whose agreement with the sub-PeV drop rests on an unvalidated ~2000-fold wind-photon enhancement, so the quantitative match is provisional.","tokens_in":5970,"tokens_out":3137,"would_cite":true,"duration_ms":33097,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Cygnus X-3's PeV photons trace to protons accelerated in its jet","keywords":["Cygnus X-3","PeVatron","cosmic-ray acceleration","microquasar jet","photo-hadronic interactions","gamma-ray emission","orbital modulation","line-driven stellar wind"],"falsifier":"Measure the spectrum of Cygnus X-3 between 0.1 and 1 PeV in two orbital phase bins (inferior and superior conjunction): the absorption hypothesis predicts a significant phase-dependent difference in the turnover energy, while free-streaming models predict a monotonic power law; observing no phase dependence would falsify the wind-absorption claim.","tokens_in":5125,"feed_emoji":"⚡","tokens_out":5595,"duration_ms":55879,"temperature":0.7,"pith_summary":"Using the recently measured photon flux from Cygnus X-3 that extends to 4 PeV and varies with the binary orbit, this paper argues that the source is a cosmic-ray accelerator in its jet, not just a leptonic emitter. The authors show that protons accelerated beyond tens of PeV produce the gamma rays by colliding with the companion star's UV photons and with wind gas; the UV channel naturally peaks near a PeV and is orbitally modulated because the viewing angle between the jet and the stellar light changes during the orbit. To explain the sharp drop of the flux just below 1 PeV, they propose that the line-driven wind of the Wolf-Rayet star traps and multiplies UV photons, increasing pair-production absorption. If this interpretation holds, Cygnus X-3 becomes the first binary system identified as a PeVatron.","feed_headline":"Cygnus X-3's PeV light points to jet-accelerated protons","feed_subtitle":"Hadronic collisions with the companion star's UV light explain the spectrum and predict an orbiting neutrino signal.","key_machinery":"The paper's key objects are (1) the photo-hadronic reaction rate on an anisotropic UV photon background, approximated by integrating the cross section with a fixed cosine µ0 of the angle between the proton and photon momenta—this angle changes over the orbit and produces the modulation—and (2) the enhanced pair-production absorption in the wind, modeled as a random walk of photons with a density enhancement factor N^{1/2} ~ v∞/Δv ~ 2000 in the region where the wind is still accelerating. The first element converts the spectrum into a probe of jet physics; the second explains the abrupt turnover near 1 PeV.","core_discovery":"The central claim is that the PeV photon emission observed from Cygnus X-3 is generated by hadronic interactions of ultra-relativistic protons accelerated in the jet. The paper demonstrates that photo-hadronic interactions on stellar UV photons produce a spectrum that rises toward a PeV and exhibits orbital modulation with the same phase as the GeV emission, consistent with observations. The fast drop below the PeV peak is then attributed to pair-production absorption on an enhanced UV photon density in the line-driven stellar wind, where multiple line scatterings amplify the photon density by a factor of order 2000 within a few stellar radii. With this adjustment, the predicted average spec","pith_inferences":["If this mechanism is generic, other high-mass X-ray binaries with Wolf-Rayet companions might also be PeVatrons, and could contribute to the galactic cosmic-ray flux; searching for orbitally modulated PeV emission from similar systems would test this.","The wind-absorption explanation makes a phase-dependent prediction: the depth of the sub-PeV drop should vary with orbital phase, because the line of sight crosses different amounts of the dense wind; a phase-resolved spectral analysis could confirm or refute it.","The random-walk enhancement factor is a rough estimate; a detailed radiative-transfer model of the line-driven wind could replace it with a first-principles prediction, turning the absorption hypothesis into a falsifiable quantitative model."],"forward_implications":["If correct, Cygnus X-3 is a PeVatron, establishing that microquasar jets accelerate cosmic rays to at least tens of PeV.","The neutrino flux from the same photo-hadronic interactions is predicted to share the orbital modulation, so a time-dependent search in neutrino telescopes would improve sensitivity.","The sharp sub-PeV turnover becomes a diagnostic of the wind's UV photon density and can be used to measure the structure of the line-driven wind.","In the authors' parameter regime, leptonic emission cannot reach PeV energies, so future PeV observations would strengthen the hadronic interpretation."],"fun_headline_variants":["Jet protons forge Cygnus X-3's PeV photons","PeV photons from Cygnus X-3: a jet proton story","Cygnus X-3's PeV peak: hadronic jet acceleration","How Cygnus X-3's jets make PeV gamma rays","Cygnus X-3: PeV light born from jet protons"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The linchpin is that the line-driven wind of the Wolf-Rayet star raises the UV photon density near the star by a factor of roughly 2000, making pair-production absorption strong enough to turn over the spectrum just below 1 PeV—the paper notes that without this enhancement, the predicted flux rises too slowly to match the data.","fun_headline_variants_meta":{"raw":{"variants":["Jet protons forge Cygnus X-3's PeV photons","PeV photons from Cygnus X-3: a jet proton story","Cygnus X-3's PeV peak: hadronic jet acceleration","How Cygnus X-3's jets make PeV gamma rays","Cygnus X-3: PeV light born from jet protons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000225,"raw_usage":{"total_tokens":1263,"prompt_tokens":666,"completion_tokens":597,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":410,"completion_tokens_details":{"reasoning_tokens":500}},"tokens_in":410,"tokens_out":597,"duration_ms":5663,"temperature":1.0,"reasoning_tokens":500,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T14:51:17.135525+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spectrum of Cygnus X-3 between 0.1 and 1 PeV in two orbital phase bins (inferior and superior conjunction): the absorption hypothesis predicts a significant phase-dependent difference in the turnover energy, while free-streaming models predict a monotonic power law; observing no phase dependence would falsify the wind-absorption claim.","supporting_citations":[],"review_version":1}