{"id":"2f60d7b2-025e-4997-98ee-86e4983c9934","arxiv_id":"2504.20920","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An intranuclear cascade Monte Carlo generator reproduces measured inelastic antideuteron-nucleus cross sections above roughly 0.65 GeV/c and predicts multiplicities and stripping fractions for ALICE-like detector nuclei.","lead":"A single-author paper applies an established intranuclear cascade model to antideuteron-nucleus collisions and reports inelastic cross sections from 100 MeV/n to about 30 GeV/n, comparing them with ALICE and older Serpukhov data. The work is relevant to dark-matter searches that look for low-energy antideuterons near Earth, where detector and galactic-transmission cross sections are needed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3) normalization depends on an unreported radius parameter; the 'no additional parameters' agreement claim is not reproducible.","rationale":"Good-faith reading: the paper applies a well-established intranuclear cascade model to scarce antideuteron data, compares with old IHEP measurements and new ALICE results, and proposes a Monte Carlo generator for dark-matter searches. The weakest point is the normalization of Eq. (3): the interaction radius is identified as the only free parameter, but its value is never stated and the code is not released, so the numerical agreement claimed in the conclusion cannot be independently verified. The stripping subtraction is secondary because stripping is a physical channel that should be modeled, but the subtraction uses the model's own prediction and is therefore not an independent check. The paper itself notes several limitations, including the approximate reproduction of experimental event selection, the need to check averaging over detector nuclei, and the omission of Coulomb and diffraction-splitting contributions; these are consistent with a conditional assessment. My proposed test would determine whether the agreement survives a reasonable variation of the radius. Because the issue is an addressable reproducibility gap rather than a demonstrated internal inconsistency, the reader's CONDITIONAL verdict should stand unchanged.","tokens_in":13016,"tokens_out":6022,"duration_ms":72552,"concrete_test":"Request the numerical value of r0 (or r_int) and the impact-parameter sampling used in the code; then re-run the INC simulation for the 13.3 GeV/c targets in Table I and for the 16O/31P ALICE cases at 2 GeV/c, using r0 = 1.1, 1.2, and 1.3 fm. Compare sigma_in and sigma_in - sigma_stripping to the IHEP and ALICE data points. If the variation in the final cross sections exceeds the quoted 2% simulation statistical uncertainty or the experimental error bars, then the claimed parameter-free agreement depends on an unstated tuning and the central claim fails; if it does not, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the INC model describes antideuteron-nucleus inelastic cross sections 'without any additional parameters' rests on Eq. (3): sigma_in = pi (r_nucl + R_d + lambda/2)^2 * N_in/N_tot. The normalization is controlled by r_nucl = r0 * A^(1/3), and Sec. II explicitly calls r_int = r_nucl + lambda/2 'the only parameter of the model.' Yet no numerical value for r0 or r_int is given anywhere in the paper. Because this prefactor multiplies every quoted cross section, all comparisons with IHEP data (Figs. 6-7, Tables I-II) and with ALICE data (Figs. 8-9) scale with this unspecified radius. The situation is compounded by Sec. III, which says r_int 'can be refined' using p, pbar, and deuteron data; this raises the possibility that the radius was adjusted to earlier data before being applied to antideuterons. If so, the claim of a parameter-free prediction for antideuterons is not established; if not, the value must still be reported for reproducibility. The paper also excludes the two lowest ALICE momentum bins and subtracts a model-dependent spectator-antiproton stripping cross section whose event-selection matching is admitted to be only approximate, but the more basic issue is that the normalization parameter is absent, so the numerical agreement cannot be independently checked.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an intranuclear cascade (INC) model for inelastic antideuteron-nucleus interactions, treating the antideuteron as a loosely bound dumb-bell of an antiproton and an antineutron. The inelastic cross section is written in Eq. (3) as a geometric prefactor times the simulated ratio Nin/Ntot, and the model is compared with IHEP antideuteron data at 13.3 and 25 GeV/c, with the old Ta multiplicity data, and with ALICE low-momentum data. The authors report satisfactory agreement and state that the model works 'without any additional parameters', proposing a Monte Carlo generator for dark-matter searches.","tokens_in":13339,"tokens_out":5985,"duration_ms":61201,"significance":"If the central claim is established, the model would provide a physically motivated alternative to Glauber-based parameterizations for antideuteron propagation in dark-matter search experiments, and the treatment of peripheral one-antinucleon interactions and stripping channels is a useful contribution. The reproduction of the old Ta multiplicities after applying the experimental cuts (Table III) is a genuine success of the model and gives this work independent value beyond the cross-section comparisons. However, as presented, the validation is not complete: the radius parameter controlling the normalization is not reported, the ALICE agreement relies on excluding the two lowest momentum bins and on a model-dependent stripping subtraction, and the IHEP comparison is partly circular because the same model stripping cross section is used to correct the data and then subtracted from the model. These issues are load-bearing for the main claim, though they are potentially fixable in revision.","major_comments":[{"comment":"The normalization of the reported cross sections is not reproducible. Equation (3) is sigma_in = pi (r_nucl + R_d + lambda/2)^2 (Nin/Ntot), and Section II states that r_int = r_nucl + lambda/2, with r_nucl = r0 A^(1/3), is 'the only parameter of the model.' No numerical value of r0 or r_int is given anywhere in the paper. Because this prefactor multiplies every quoted cross section, the values in Tables I-II and the curves in Figs. 6-9 scale with this unspecified radius. Please report the value(s) used, state whether r0 was fitted to the (anti)proton/deuteron data shown in Figs. 1-5 before the antideuteron calculation, and provide an uncertainty estimate. Without this, the claim in Section V of agreement 'without any additional parameters' cannot be checked.","section":"II, Eq. (3); III"},{"comment":"The ALICE comparison is presented as the main validation for the low-momentum region, but the agreement is obtained only after excluding the first two momentum bins and subtracting the model's own spectator-antiproton stripping cross section. The paper itself notes that the experimental event selection is reproduced only 'fairly approximately' in the discussion near Fig. 4. The two excluded bins correspond to antideuteron kinetic energies of roughly 20 and 40 MeV/nucleon, below the declared validity range of 100 MeV/n, so the exclusion may be justified, but this should be stated explicitly and the comparison should be shown with all bins present. Please quantify the sensitivity of the stripping subtraction and of the normalization to model assumptions, and plot model bands rather than single curves.","section":"III, Figs. 8-9"},{"comment":"The IHEP comparison is partly circular. The text states that the experimental sigma_in values in Tables I and II 'include the calculated corrections sigma_st^barp [3]', and the plotted model quantity is sigma_in_calc - sigma_st^barp, i.e., the same model stripping cross section is used to correct the data and then subtracted from the model. For 12C at 13.3 GeV/c the tabulated experimental value equals the model value sigma_in - sigma_st by construction. The comparison should instead be made directly with the measured quantity sigma'_in = sigma_in - sigma_st^barp, with the model's stripping cross section reported separately, so that the stripping channel is tested independently.","section":"III, Tables I-II, Figs. 6-7"}],"minor_comments":[{"comment":"The heading of Section V is misspelled as 'CONCLUTION', and 'antideutron'/'deutron' are used in several places where 'antideuteron'/'deuteron' is intended.","section":"Throughout"},{"comment":"Equations (2) and (3) contain unbalanced parentheses, and 'Rd = 2 , 16 fm' should read 2.16 fm.","section":"II, Eqs. (2)-(3)"},{"comment":"The model curves in Figs. 8-9 are shown without uncertainty bands; the text gives only a 2% statistical error for the Monte Carlo, not the systematic uncertainty from the averaged-nucleus approximation or from the stripping subtraction.","section":"III, Figs. 8-9"},{"comment":"The paper acknowledges that diffraction and Coulomb splitting of the deuteron are not included at low energies; this limitation should be stated in the abstract or conclusions where the claimed energy range is quoted.","section":"III, Fig. 3"}],"recommendation":"major_revision","confidential_remarks":"I am not recommending rejection because the model's physical content is relevant and the missing radius parameter is, in principle, reportable. However, the authors should be asked to provide the numerical value and provenance of r0/r_int, to remove or substantially qualify the 'without any additional parameters' claim if the radius was adjusted to earlier data, and to restructure the IHEP comparison so that the stripping correction is not used on both sides of the comparison. The generator is mentioned but not made available; providing the code or at least a clear specification of all input parameters would strengthen the reproducibility of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a competent application of an established cascade model to a scarce-data problem, and the paper is worth a referee's time, but the central numerical claim is not reproducible as submitted because the normalization parameter is never given.\n\nWhat's new: the INC treatment of the antideuteron as a two-center cascade is not new (the same author did Ta at 12.2 GeV/c), but the paper extends it to the ALICE momentum range, to averaged detector nuclei, and to multiplicity/momentum predictions for 16O and 31P. Those predictions are the practically useful part for AMS/GAPS. The comparisons to old IHEP data and the Ta multiplicity table are genuine external checks, and the model does reproduce the high-momentum ALICE points and the Ta multiplicities once stripping is subtracted. That is real evidence and should be credited.\n\nThe soft spots are substantial. Eq. (3) multiplies every cross section by π(r_nucl + R_d + λ/2)^2, and the text calls r_int = r_nucl + λ/2 the only parameter. But no value for r0 or r_int appears anywhere. Since the formula for r_nucl uses r0 A^(1/3) and r0 is never stated, all quoted numbers scale with an unknown. Section III also says r_int can be 'refined' using p, pbar, and d data, so the 'without any additional parameters' claim in the conclusion is not established. This is a reproducibility failure, not a minor omission. Second, the ALICE agreement is obtained after dropping the two lowest-momentum bins (which are off by about a factor of 2.1) and subtracting the model's own spectator stripping cross section; the paper admits the event-selection match is only approximate. Third, the generator is described but not released, so the central deliverable for the experimental community cannot be used.\n\nNone of this makes the model silly. The physics is plausible, and the paper is honest about the stripping complication and about the low-energy failure. It just needs a major revision: report r_int and r0, show sensitivity of all cross sections to that parameter, compare with and without stripping at all ALICE bins, and make the generator accessible.\n\nWho it's for: people building antinucleus propagation into AMS/GAPS or cosmic-ray simulation. They would use the generator if it existed. As it stands, I wouldn't cite the numbers, but I would send it to peer review. A serious referee can force the parameter to be reported.","headline":"A plausible INC extension to antideuteron–nucleus scattering, but the cross-section normalization depends on an unreported radius and the generator is not released, so the headline agreement claim is not reproducible.","tokens_in":13823,"tokens_out":2502,"would_cite":false,"duration_ms":28273,"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":"An intranuclear cascade model that treats the antideuteron as a dumb-bell of two antinucleons reproduces measured inelastic antideuteron-nucleus cross sections without additional parameters.","keywords":["antideuteron","intranuclear cascade model","inelastic cross section","dark matter search","antinucleus-nucleus interaction","stripping reaction","Monte Carlo generator"],"falsifier":"A precise measurement of the inelastic antideuteron cross section on a light nucleus at momenta between about 0.65 and 4 GeV/c, with explicit identification of spectator antiprotons, would settle the claim: if the model's $\\sigma_{in}-\\sigma_{\\bar p}^{\\mathrm{st}}$ misses those data beyond the quoted simulation statistics, the geometric normalization or the stripping subtraction is wrong.","tokens_in":12806,"feed_emoji":"⚛️","tokens_out":7037,"duration_ms":67184,"temperature":0.7,"pith_summary":"The paper argues that the intranuclear cascade model, in which the antideuteron is treated as a loosely bound dumb-bell of an antiproton and an antineutron, can predict the inelastic antideuteron-nucleus cross section over about 100 MeV/n to 25–30 GeV/n for a wide range of target nuclei. Using this model, the author obtains cross sections that agree with the recent low-energy measurements once events with a spectator antiproton, the stripping channel, are subtracted and the two lowest-momentum bins are set aside. The motivation is dark-matter searches: antideuterons are a promising low-background signature of dark-matter annihilation, and computing their expected flux near Earth requires knowing how they interact with interstellar matter and with detector material. On this basis the paper delivers a Monte Carlo generator for simulating antideuteron interactions with nuclei. The central claim, stated by the author, is that without any additional parameters the model is in satisfactory agreement with the available experimental data.","feed_headline":"Cascade model matches antideuteron-nucleus data","feed_subtitle":"A dumb-bell treatment of the antideuteron predicts inelastic cross sections for dark-matter searches.","key_machinery":"The machinery is the intranuclear cascade (INC) model: a quasi-classical statistical simulation of successive hadron-nucleon collisions inside the target nucleus, with the target treated as a Fermi gas in a potential well, followed by evaporation or explosive disintegration of the residual nucleus and a local decrease of nuclear density as the cascade develops. The antideuteron is modelled as a dumb-bell of an antiproton and an antineutron at a fixed separation $l = 2R_d = 4.32\\times10^{-13}$ cm with an isotropically distributed axis; the internal momentum distribution is the square of the Fourier transform of an approximate deuteron wavefunction. The load-bearing formula is $\\sigma_{in} = \\pi(r_{\\mathrm{nucl}}+R_d+\\lambda/2)^2 N_{\\mathrm{in}}/N_{\\mathrm{tot}}$, where $r_{\\mathrm{nucl}} = r_0 A^{1/3}$ is the target radius, $\\lambda$ the de Broglie wavelength, and $N_{\\mathrm{in}}/N_{\\mathrm{tot}}$ the simulated probability of an inelastic interaction. Stripping is treated as a special case of the inelastic channel: when one antinucleon misses the nucleus the event remains in $\\sigma_{in}$, but it is subtracted when comparing with experiments that rejected spectator-antiproton events.","core_discovery":"The central claim is that the inelastic antideuteron-nucleus interaction can be simulated as two coupled intranuclear cascades initiated by the antideuteron's constituent antinucleons, and that this reproduces the measured cross sections. Taking the antideuteron as a dumb-bell with fixed separation $l = 2R_d$, the model computes the inelastic cross section from the geometric area $\\pi(r_{\\mathrm{nucl}}+R_d+\\lambda/2)^2$ times the simulated ratio of inelastic to total encounters. After subtracting the calculated antiproton-stripping cross section $\\sigma_{\\bar p}^{\\mathrm{st}}$, the channel in which one antinucleon passes through without interacting, the predicted $\\sigma_{in}$ matches the low-energy detector data except at the two lowest momentum bins, and matches the older high-momentum measurements for carbon, aluminium, copper and lead at 13.3 and 25 GeV/c, and for tantalum at 12.2 GeV/c. The paper also reports a structural result: peripheral interaction of only one antinucleon dominates, with about 72% of inelastic events on $^{16}\\mathrm{O}$ at 2 GeV/c involving a single antinucleon, while both antinucleons annihilate in only about 20% of events.","pith_inferences":["The same dumb-bell cascade construction could be extended to anti-$^3\\mathrm{He}$ and anti-$^4\\mathrm{He}$, which the paper notes are observed cosmic-ray species; treating them as rigid clusters of antinucleons would produce cross-section predictions that cosmic-ray antinuclei searches could use.","Because the paper never states the numerical value of the interaction radius $r_{\\mathrm{int}}$, a reader's first reproducibility check would be to fit $r_{\\mathrm{int}}$ on the deuteron-nucleus data shown in the paper and see whether the antideuteron predictions move into or out of agreement.","The low-momentum excess at 0.3–0.65 GeV/c is plausibly the signature of the Coulomb polarization and diffraction-splitting contributions that the paper lists but does not implement; adding those mechanisms to the cascade would be a concrete, testable next step."],"forward_implications":["For target nuclei and energies where no data exist, the model supplies explicit $\\sigma_{in}$ predictions covering $100\\ \\mathrm{MeV/n} \\le T_{\\mathrm{kin}} \\le 25$–$30\\ \\mathrm{GeV/n}$ and nuclei from $^{12}\\mathrm{C}$ upward.","The delivered Monte Carlo generator can be used to simulate antideuteron interactions in the detector material of cosmic-ray search experiments, including the multiplicity and spectra of secondary particles.","Any comparison between model and experiment must subtract the stripping channel when the experiment rejects spectator-antiproton events; the paper shows that this subtraction is essential for agreement.","Peripheral interaction of a single antinucleon is the dominant channel at the energies studied, so antideuteron interactions cannot be approximated as simple twofold annihilation; the model predicts about 72% of inelastic events on $^{16}\\mathrm{O}$ at 2 GeV/c involve only one antinucleon."],"supporting_citations":[{"why":"Provides the low-energy antideuteron inelastic cross-section data that the model is mainly compared against.","marker":"[14]"},{"why":"Supplies the measured antideuteron absorption and stripping cross sections at 13.3 GeV/c shown in Table I.","marker":"[2]"},{"why":"Supplies the measured absorption cross sections at 25 GeV/c shown in Table II.","marker":"[3]"},{"why":"Provides the antideuteron-tantalum interaction data at 12.2 GeV/c used to test exclusive features of the model.","marker":"[4]"},{"why":"Defines the intranuclear cascade model of antinucleon-nucleus inelastic interactions that the antideuteron simulation is built on.","marker":"[29]"},{"why":"Provides the deuteron-nucleus interaction treatment that motivates the dumb-bell representation of the antideuteron.","marker":"[32]"},{"why":"Supplies the compilation of (anti)proton and deuteron cross-section data used to validate the model before applying it to antideuterons.","marker":"[33]"},{"why":"Describes the Glauber-based parameterization used in the detector simulation that the model is implicitly being compared against for the low-energy data.","marker":"[25]"}],"fun_headline_variants":["Antideuteron dumb-bell model nails nuclear collision data","Two-cascade model reproduces antideuteron-nucleus cross sections","Peripheral antinucleon hits dominate antideuteron interactions","Coupled cascades simulate antideuteron-nucleus inelastic events","New generator predicts antideuteron cross sections for dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the single geometric interaction radius, whose numerical value is never stated, can be fixed once and reused for all target nuclei and energies, and that the experimental rejection of spectator-antiproton events corresponds to the model's stripping subtraction.","fun_headline_variants_meta":{"raw":{"variants":["Antideuteron dumb-bell model nails nuclear collision data","Two-cascade model reproduces antideuteron-nucleus cross sections","Peripheral antinucleon hits dominate antideuteron interactions","Coupled cascades simulate antideuteron-nucleus inelastic events","New generator predicts antideuteron cross sections for dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000726,"raw_usage":{"total_tokens":3261,"prompt_tokens":958,"completion_tokens":2303,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":2210}},"tokens_in":574,"tokens_out":2303,"duration_ms":16912,"temperature":1.0,"reasoning_tokens":2210,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:16:27.436456+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A precise measurement of the inelastic antideuteron cross section on a light nucleus at momenta between about 0.65 and 4 GeV/c, with explicit identification of spectator antiprotons, would settle the claim: if the model's $\\sigma_{in}-\\sigma_{\\bar p}^{\\mathrm{st}}$ misses those data beyond the quoted simulation statistics, the geometric normalization or the stripping subtraction is wrong.","supporting_citations":[{"cited_title":"Measurement of the low-energy antideuteron inelastic cross section","cited_arxiv_id":"2005.11122","evidence_quote":"Provides the low-energy antideuteron inelastic cross-section data that the model is mainly compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the measured antideuteron absorption and stripping cross sections at 13.3 GeV/c shown in Table I."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the measured absorption cross sections at 25 GeV/c shown in Table II."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the antideuteron-tantalum interaction data at 12.2 GeV/c used to test exclusive features of the model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the intranuclear cascade model of antinucleon-nucleus inelastic interactions that the antideuteron simulation is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the deuteron-nucleus interaction treatment that motivates the dumb-bell representation of the antideuteron."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the compilation of (anti)proton and deuteron cross-section data used to validate the model before applying it to antideuterons."},{"cited_title":"Uzhinsky, J","cited_arxiv_id":null,"evidence_quote":"Describes the Glauber-based parameterization used in the detector simulation that the model is implicitly being compared against for the low-energy data."}],"review_version":1}