{"id":"69fa42d7-bbae-4f1a-937a-7877b8bc0604","arxiv_id":"1909.01277","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Using the Cannon-Ball model's original priors, the author computes a positron source spectrum whose shape matches AMS data, with an absolute normalization 23% higher than observed.","lead":"This paper claims that a 20-year-old model of gamma-ray bursts, the Cannon-Ball model, predicts the cosmic-ray positron excess measured by AMS without fitting any parameters. If true, the long-standing positron excess would not need dark matter or pulsars as an explanation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Wind attenuation A(l)=exp[-(l_tr/l)^2] contradicts the stated n∝l^-2 profile; the correct form is exp[-l_tr/l], so the predicted flux and the 1.23 normalization rest on an incorrect attenuation law.","rationale":"The reader's weakest assumption identifies exactly the same load-bearing flaw, so my agreement is 'agree'. The stress-test confirms that the wind attenuation formula is not merely an unvalidated prior but contradicts the paper's own stated density profile n∝l^-2. This is decisive for the central claim because N→e(γ) determines both the source normalization and the weighting over Lorentz factors; the correct exponential form changes the integral by roughly two orders of magnitude at typical γ, so the excellent agreement claimed in the abstract and Sec. V cannot survive as stated. I also noted the dimensional inconsistency in Eq.(3) (the erfc term lacks l_tr), which independently prevents verification of the printed formula. I considered whether the CB model could still work with a recalibrated normalization, but the paper's explicit claim is that the original unmodified priors and no parameter fitting reproduce the data; the attenuation error invalidates that claim. The idea remains interesting and the propagation framework is transparent, but the derivation gap is load-bearing. The reader's REJECT verdict is therefore appropriate, and my analysis does not move it.","tokens_in":13328,"tokens_out":14926,"duration_ms":128682,"concrete_test":"Recompute the source flux with the physically correct attenuation A(l)=exp(-l_tr/l) while keeping all other priors (Table I) and the propagation treatment unchanged. Specifically, replace I in Eq.(3) with I_correct(L)=∫_0^L exp(-a/x) dx = L e^{-a/L} + a Ei(-a/L), where a=σ_T Σ/m_p and L=√3γR_pp, propagate through Eqs.(4)-(11), and compare the resulting E^3 dΦ/dE to the AMS data in Fig. 2. If the predicted-to-observed normalization changes by more than a factor of 5 from the reported 1.23, the consistency check is an artifact of the incorrect attenuation law.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III.C states that for a wind with n_e ∝ l^-2, the survival probability of a positron produced at distance l from the supernova is A(l)=exp[-(l_tr/l)^2], with l_tr=σ_T Σ/m_p. This is internally inconsistent with the paper's own density profile. The optical depth from l outward is τ(l)=σ_T ∫_l^∞ n(l') dl' = σ_T Σ/(m_p l)=l_tr/l, so A(l) must be exp(-l_tr/l). The squared exponent would require τ ∝ l^-2, i.e. n ∝ l^-3. Since A(l) multiplies the integrand in Eq.(2) to give Eq.(3), N→e(γ) and the γ-weighting in Fig. 3, and hence the entire source flux Eq.(11), are computed with the wrong attenuation law. A numerical estimate at γ=10^3 (l_tr≈4×10^15 cm, l_max≈1.7×10^15 cm) gives I_correct=∫_0^L e^{-a/l} dl ≈ 5×10^13 cm while Eq.(3)'s exponent gives I_sq ≈ 5×10^11 cm, a factor of approximately 100. Thus the claimed '1.23 times larger' normalization and the 'no fitting' agreement do not test the CB model as stated. In addition, Eq.(3) as printed has a dimensional mismatch: the erfc term is missing the factor l_tr, so the formula cannot be evaluated as written. This is an internal derivation error, not a disagreement with external consensus; a corrected calculation might still be viable, but the paper's central claim is not supported by the calculation as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the AMS-02 positron excess originates from positrons produced in proton-proton collisions between relativistic cannonballs of the Cannon-Ball GRB model and the wind of the exploding star. The authors derive a source spectrum using model priors fixed in earlier CB-model publications, with no parameters fitted to the AMS data, add a diffuse background taken from Lipari's fit to earlier AMS data, and compare the sum with the observed AMS-02 spectrum. They report a 1.23 times normalization overestimate, which they interpret as a satisfactory consistency check.","tokens_in":13736,"tokens_out":12214,"duration_ms":115480,"significance":"If the calculation were correct, this would be a noteworthy result: a pre-existing, non-dark-matter, non-pulsar astrophysical model predicting both the shape and the normalization of the AMS positron flux without free parameters. The paper's transparency in writing out the derivation is a strength, since it allows the main error to be identified precisely. However, the error in the wind-attenuation formula undermines the central quantitative claim, so the significance is conditional on a corrected calculation.","major_comments":[{"comment":"The survival probability A(l)=exp[-(l_tr/l)^2] is inconsistent with the stated wind density profile n(l)=Sigma/(m_p l^2). The opacity from a production point l to infinity is tau(l)=sigma_T integral_l^infinity n(l') dl' = l_tr/l, so the correct survival probability is exp(-l_tr/l); the squared exponent would correspond to n proportional to l^-3. Since A(l) multiplies the integrand in Eq. (2) to produce N->e(gamma) in Eq. (3), and this quantity enters the source spectrum in Eqs. (4), (6), and the final flux in Eq. (11), the central claims about the shape and the 1.23 normalization are not supported by the calculation as presented. The numerical impact is parameter-dependent; for gamma approximately 1000, with l_tr approximately 4e15 cm and l_max approximately 1.7e15 cm, the two integrals differ by close to an order of magnitude. The authors should repeat the calculation with the attenuation law A(l)=exp(-l_tr/l) and check whether the agreement with AMS persists.","section":"Section III.C, Eq. (3)"},{"comment":"The expression for I in Eq. (3) has a dimensional inconsistency: the first term, l_max exp[-(l_tr/l_max)^2], has dimensions of length, while the second term, minus sqrt(pi) erfc[l_tr/l_max], is dimensionless. The correct expression for the integral of exp[-(l_tr/l)^2] contains an additional factor l_tr multiplying the erfc term. As printed, Eq. (3) cannot be evaluated in physical units, which compounds the error identified in the previous comment.","section":"Eq. (3)"},{"comment":"The adopted diffuse background is Lipari's parametrization of earlier AMS data, and the source term is then defined as the residual of the same spectrum. The comparison in Fig. 2 is therefore partly circular: the model is effectively tested against a leftover after a fit to the total data has been subtracted. To substantiate the 'no fitting' claim, the authors should show that their conclusions are insensitive to plausible alternatives for the diffuse term, for example a GALPROP/DRAGON secondary spectrum or a different analytic background. Without such a test, the 1.23 normalization and the shape agreement are less decisive than stated.","section":"Sections I and V, Fig. 2"}],"minor_comments":[{"comment":"The statement that an increase by 40% (or 19%) of the wind surface density Sigma would reduce the cited normalization by a factor of two (or 1/1.23) is confusing, since a 19% increase producing a reduction by a factor 1/1.23 is a much weaker effect than a 40% increase producing a factor of two; please clarify the intended functional dependence.","section":"Section V, last paragraph"},{"comment":"The unit of Sigma appears as '10 16 g/cm' without a superscript; it should read 10^16 g cm^{-1}. In addition, the caption of Fig. 2 omits the symbol lambda in the label for the blue curve.","section":"Table I and Fig. 2 caption"},{"comment":"The text introducing A(l) states that l_tr is the distance at which the remaining optical depth of the wind is unity; with the correct attenuation law A(l)=exp(-l_tr/l) this definition is accurate, but it is inconsistent with the squared exponent used in Eq. (3). Please harmonize the definition with the final formula.","section":"Section III.C"},{"comment":"Reference [28] contains a typo ('or As- trophys. J.'), and the arXiv identifier for reference [8] should be checked for accuracy.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The derivation error in Eq. (3) is severe and must be addressed by recomputing the source flux with the correct attenuation law. The authors should also provide a sensitivity analysis of the adopted diffuse background, since the current comparison is partly circular. I would not consider acceptance of the manuscript in its present form; the issues are significant enough to require a major revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: De Rújula applies the Cannon-Ball model to the AMS positron excess and claims a parameter-free prediction, but the wind-attenuation expression in Eq. (3) is inconsistent with the stated n∝l^-2 wind profile. The correct survival factor is exp(-l_tr/l), not exp[-(l_tr/l)^2]; the printed integral also has a dimensional mismatch in the erfc term. This is a load-bearing error: at γ≈10^3 it changes the source normalization by roughly two orders of magnitude, so the alleged '1.23 times larger' agreement is an artifact.\n\nWhat's good: the paper is a clearly written, honest attempt to explain the positron excess without dark matter or pulsars. It gives a real no-parameter calculation (modulo a few priors) with a detailed treatment of ICS energy losses and escape, and it is transparent about the arbitrary normalization in Fig. 2 and the fitted Lipari diffuse term. The sensitivity to λ and to the LF distribution is shown. Credit where due: this is a serious test of an unfashionable model, not a toy.\n\nThe soft spots: the central derivation error; the diffuse term is fitted to the same AMS data so the comparison is not fully independent; and the 'no fitting' claim is softened by the choice U_in=10 U_loc and the (largely cosmetic) multiplicity correction. The novelty is also limited—Dado & Dar 2016 already fitted the positron excess with CB parameters; the new step is the no-parameter version, which is exactly what fails.\n\nIf the attenuation formula were corrected, the flux would be far above the observed one, and the shape would shift. A factor ~3 change in Σ might bring the normalization down, but that is a fit, and the paper does not provide such an alternative. The citation pattern is fine.\n\nRecommendation: send to peer review—this deserves a serious referee who can check the integral—but expect rejection or major revision. The idea remains interesting, but the calculation as presented does not support the conclusion.","headline":"Cannon-Ball explanation of the AMS positron excess, but the wind-attenuation formula is internally inconsistent and the no-parameter claim collapses.","tokens_in":14255,"tokens_out":7716,"would_cite":false,"duration_ms":69614,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.70.Sa","14.60.Cd","97.60.Bw","96.60.tk"],"model":"deepseek-v4-flash","headline":"A supernova cannonball gamma-ray-burst model reproduces the AMS positron spectrum without any fitted parameters.","keywords":["cosmic-ray positrons","AMS-02","Cannon-Ball model","gamma-ray bursts","positron excess","supernova winds","primary cosmic-ray sources","parameter-free prediction"],"falsifier":"Recompute Eq. (3) with the optical depth that follows from the paper's own wind density profile $n(l) = \\Sigma/(m_p l^2)$, namely $\\tau(l) = \\sigma_T \\Sigma/(m_p l)$, so that $A(l) = \\exp(-l_{\\rm tr}/l)$ instead of $\\exp(-(l_{\\rm tr}/l)^2)$; if the resulting source normalization moves away from the observed value by more than the stated uncertainties, the no-parameter match rests on the unstated attenuation ansatz.","tokens_in":13097,"feed_emoji":"🛰️","tokens_out":11077,"duration_ms":103373,"temperature":0.7,"pith_summary":"The paper argues that the rise in cosmic-ray positrons measured by AMS—the excess over secondary production—does not require dark matter or pulsars. The excess, it claims, is made at the source: relativistic plasma \"cannonballs\" ejected by stripped-envelope supernovae collide with the pre-supernova wind, and pion and kaon decays in those collisions produce positrons. Using the model's original parameter values and no fitting, the predicted source spectrum matches the AMS data in shape and, within a factor 1.23, in normalization. If this is right, the positron excess is a direct view of the same engine that generates gamma-ray bursts and the knees in the cosmic-ray spectra.","feed_headline":"Supernova cannonballs explain AMS positron excess","feed_subtitle":"A model of supernova cannonballs matches AMS data with zero fitted parameters—no dark matter or pulsars needed.","key_machinery":"The central object is the Cannon-Ball (CB) model of gamma-ray bursts: a supernova ejects a compact blob of ordinary plasma with baryon number $\\sim 10^{50}$ at Lorentz factor $\\sim 10^3$. The load-bearing step is the wind-attenuation factor $A(l) = \\exp[-(l_{\\rm tr}/l)^2]$ of Eq. (3), which gives the probability that a positron born at distance $l$ from the supernova escapes the wind without being reabsorbed, with $l_{\\rm tr} = \\sigma_T \\Sigma / m_p$ the distance at which the wind's optical depth is unity. This factor controls the number of surviving positrons $N_{\\to e}(\\gamma)$, and through Eq. (6) it up-weights higher Lorentz factors, because faster cannonballs keep colliding with the wind at larger distances where the wind is thinner. That weighting, together with the model's gamma distribution $D(y) = \\exp[-((y-y_0)/\\sigma)^2]$ with $y_0 = 6.3$ and $\\sigma = 0.5$, sets the shape and peak of the predicted source spectrum.","core_discovery":"The central claim is that the AMS positron excess is a primary source effect of the Cannon-Ball model of gamma-ray bursts, not a dark-matter or pulsar signal. A cannonball with baryon number $N_B \\sim 10^{50}$ and Lorentz factor $\\gamma\\sim 10^3$ crosses the dense wind of its parent supernova; $pp$ and $pn$ collisions feed the chain $pp$ (or $pn$) $\\to \\pi$ or $K \\to \\mu \\to e^+$, and the positrons that survive wind attenuation are then propagated with a lifetime that combines galactic escape, inverse-Compton losses on starlight, and losses on the cosmic background, far-infrared, and magnetic-field photon baths. The author adopts a single parametrization for the diffuse secondary background and derives the source term from Eq. (11); the predicted source normalization is 1.23 times the observed one. The paper presents this factor as a satisfactory consistency check rather than a discrepancy, noting that a 19% change in the assumed wind surface density would bring the two into exact agreement.","pith_inferences":["A testable extension the paper does not pursue: the same $pp$-collision machinery predicts a corresponding positron energy distribution at the source, so a simultaneous parameter-free fit to both the AMS positron and electron spectra would broaden the claim beyond positrons alone.","The attenuation functional form is the fragile point; recomputing Eq. (3) with the optical depth implied by the paper's own $n(l) \\propto l^{-2}$ profile, namely $A(l) = \\exp(-l_{\\rm tr}/l)$, would show whether the normalization match survives an alternative and arguably more natural ansatz.","If future high-statistics data resolve the arrival direction of the positron excess, the wind-attenuation weighting predicts a small anisotropy toward the inner Galaxy, whereas dark-matter annihilation would be roughly isotropic; this could separate the two origins empirically."],"forward_implications":["If the CB source is correct, the positron excess and the proton, helium, iron, and $e^+e^-$ knees all share one origin, since the same Lorentz-factor distribution $D(\\gamma)$ is used for each.","The AMS normalization then becomes a counting measurement of galactic cannonball-launching supernovae, fixing their rate at about one per century with ten cannonballs per event.","The positrons should be accompanied by prompt hard photons, but the paper estimates only about $5 \\times 10^{-6}$ of a GRB's gamma-ray count, so no observable hard-photon counterpart is expected.","Because all inputs are fixed priors, any future deviation of the measured positron flux from the predicted shape would directly constrain the wind surface density, the starlight energy density, or the CB Lorentz-factor distribution.","The factor-1.23 normalization gap provides a sharp check: a 19% increase in the wind surface density removes it entirely, so better wind-density measurements would decide whether the match is meaningful."],"supporting_citations":[{"why":"provides the AMS-02 positron flux that is the target the model must reproduce.","marker":"[1]"},{"why":"supplies the adopted parametrization of the diffuse secondary positron background, Eq. (1).","marker":"[3]"},{"why":"defines the Cannon-Ball model, its priors (gamma, N_B, wind surface density), and the galactic GRB/SN rate.","marker":"[6]"},{"why":"gives the Lorentz-factor distribution D(gamma) of Eq. (5), calibrated to the proton, He, Fe, and e+e- knees.","marker":"[8]"},{"why":"supplies the positron energy distribution F(x, E_p) from pp/pn collisions used in Eq. (4).","marker":"[16]"},{"why":"provides the independent calculation of F(x, E_p) that the paper relies on at TeV energies.","marker":"[17]"},{"why":"the AMS B/C measurement used to justify the high-rigidity confinement exponent beta_conf = 0.6.","marker":"[25]"},{"why":"models the interstellar starlight energy density used for the U_in = 10 U_loc choice in Eq. (8).","marker":"[29]"}],"fun_headline_variants":["Cannonball source explains AMS positron excess with zero fitted parameters","AMS positron excess from cannonball primaries, no dark matter needed","Cannonballs, not pulsars or dark matter, produce AMS positrons","Zero-parameter cannonball model reproduces AMS positron flux"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the assumed attenuation formula $A(l) = \\exp[-(l_{\\rm tr}/l)^2]$ for how likely a positron is to escape the supernova wind; the paper's own wind density profile would suggest a different decay law, and the source flux depends critically on this choice.","fun_headline_variants_meta":{"raw":{"variants":["Cannonball source explains AMS positron excess with zero fitted parameters","AMS positron excess from cannonball primaries, no dark matter needed","Cannonballs, not pulsars or dark matter, produce AMS positrons","Zero-parameter cannonball model reproduces AMS positron flux"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000569,"raw_usage":{"total_tokens":2671,"prompt_tokens":898,"completion_tokens":1773,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":514,"completion_tokens_details":{"reasoning_tokens":1693}},"tokens_in":514,"tokens_out":1773,"duration_ms":12226,"temperature":1.0,"reasoning_tokens":1693,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:23:09.768636+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute Eq. (3) with the optical depth that follows from the paper's own wind density profile $n(l) = \\Sigma/(m_p l^2)$, namely $\\tau(l) = \\sigma_T \\Sigma/(m_p l)$, so that $A(l) = \\exp(-l_{\\rm tr}/l)$ instead of $\\exp(-(l_{\\rm tr}/l)^2)$; if the resulting source normalization moves away from the observed value by more than the stated uncertainties, the no-parameter match rests on the unstated attenuation ansatz.","supporting_citations":[{"cited_title":"Aguilar el al., Phys","cited_arxiv_id":null,"evidence_quote":"provides the AMS-02 positron flux that is the target the model must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the positron energy distribution F(x, E_p) from pp/pn collisions used in Eq. (4)."},{"cited_title":"Moskalenko & A.W","cited_arxiv_id":null,"evidence_quote":"provides the independent calculation of F(x, E_p) that the paper relies on at TeV energies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"the AMS B/C measurement used to justify the high-rigidity confinement exponent beta_conf = 0.6."}],"review_version":1}