{"id":"d8b826ed-5258-4c17-b608-badd0ee552fa","arxiv_id":"2501.03966","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"DESMO, a simulated balloon experiment, could measure the deuteron-to-proton ratio in cosmic rays from 10 to 100 GeV per nucleon with 5-8% statistical precision using multiple scattering and Cherenkov ring imaging.","lead":"This paper proposes a compact balloon-borne detector that would measure cosmic-ray deuterium up to 100 GeV per nucleon using particle scattering instead of a magnet. A generalist reader might care because the origin of deuterium in cosmic rays is an open question that affects models of galactic propagation and dark matter searches.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The D/p precision claim rests on GEANT4 template shapes that are not validated against data; deferred fragmentation systematics (2–5% each) are comparable to the 5–8% statistical errors.","rationale":"The reader's weakest-assumption analysis correctly identifies the template shapes from GEANT4 as the load-bearing element, and the paper's own deferral of atmospheric and in-detector fragmentation systematics (Section IV) supports that concern. I considered whether the Appendix B derivation of the displacement resolution is a more fundamental flaw; while the intermediate integrals in the text are indeed inconsistent (the quoted values correspond to a folded Gaussian rather than the half-range integrals written), the final formula used in the optimization matches the folded-distribution variance, and the actual detector performance is taken from the GEANT4 simulation rather than from that analytic formula. Therefore that issue is not the primary threat to the central claim. The OPERA multiple-scattering technique [19] provides independent support for the general momentum-measurement principle, and the DESMO closure test demonstrates statistical feasibility, but it does not validate the absolute template shapes. Because the claimed few-percent precision is conditional on unquantified systematics of comparable size, the conditional verdict stands; the proposed test would settle whether the precision claim survives plausible template perturbations. No reason to escalate to reject, as the detector concept is plausible and the paper is honest about the deferred uncertainties, but the precision claim cannot be accepted as stated without that validation.","tokens_in":24441,"tokens_out":8747,"duration_ms":87841,"concrete_test":"Generate a set of 'perturbed' templates by (a) repeating the GEANT4 simulation with two different hadronic physics lists (e.g., FTFP_BERT vs QGSP_BERT) and (b) reweighting the template shapes by the 2–5% atmospheric-fragmentation uncertainty quoted from BESS [41] (e.g., varying the D survival probability linearly with energy across the 20–200 GV range). Re-run the 60-day template fit to the nominal simulated data with these perturbed templates. If the reconstructed D/p shifts by more than the quoted 5–8% statistical precision in any of the six Ek bins, the 'few percent' claim is not robust and a calibration strategy (e.g., beam-test data or a dedicated atmospheric-fragmentation measurement) is required before the sensitivity projection can be accepted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Section V) is that DESMO can measure the D/p fraction with few-percent precision in 20–200 GV, based on a template fit to the inverse-displacement distribution (Figure 8). This projection is a Monte Carlo closure test: both the simulated 60-day data and the p/D templates are generated with the same GEANT4 simulation, so the fit will extract the injected fraction by construction up to statistical fluctuations. The accuracy of the extracted D/p in reality therefore hinges entirely on whether the simulated template shapes (the 1/d distributions for p and D at each reconstructed Ek bin) match the true detector response. The paper explicitly defers the two dominant sources of shape distortion: fragmentation of D and heavier nuclei in the residual atmosphere above the balloon (estimated 2–5%) and in the detector itself (estimated <5%) (Section IV). These are comparable to the quoted 5–8% statistical precision. If the templates are wrong at the few-percent level, the fitted D/p could be biased by more than the quoted precision, potentially mimicking or hiding the primary-versus-secondary spectral signature that the experiment is designed to discriminate. No beam-test or independent data validation of the GEANT4 multiple-scattering and hadronic-interaction model for this specific configuration is presented. The Appendix B analytical optimization also contains intermediate integrals that are mis-stated, though the final folded-distribution formula is plausible; this does not change the main concern.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes DESMO, a compact balloon-borne detector combining a multiple-scattering module (MSIS) with a RICH velocity detector, to measure cosmic-ray deuterons. The authors use GEANT4 simulations of the full detector, including a BDT-calibrated RICH and a template-fit analysis of simulated 60-day data, to extract the D/p ratio in six reconstructed energy bins covering roughly 10-100 GeV/nucleon. They report a statistical precision of 5-8% on D/p and argue that this would allow discrimination among secondary, mixed, and primary deuterium scenarios.","tokens_in":24760,"tokens_out":8006,"duration_ms":79101,"significance":"The proposed technique is original and the simulation study is unusually detailed, including realistic optical effects, SiPM dark counts, detector acceptance, and a minimization-based template fit. If validated experimentally, DESMO would offer a low-cost path to deuterium measurements at rigidities beyond AMS-02, potentially settling the current debate about a primary-like deuteron component. The central limitation is that the claimed precision is statistical only; the two dominant systematic uncertainties are deferred and are of comparable size to the quoted errors, and the detector response has not been validated against beam or flight data. The paper is therefore best read as a feasibility study whose quantitative physics case will stand or fall on the unvalidated template accuracy.","major_comments":[{"comment":"The quoted 5-8% precision on D/p is purely statistical. The paper itself cites 2-5% uncertainty from fragmentation in the residual atmosphere and below 5% from in-detector fragmentation, both on the same order as the statistical error. Because the p and D template shapes are generated with GEANT4 and are not validated against any beam-test or independent data, an uncontrolled few-percent template distortion could bias the fitted D/p by an amount comparable to the spectral signature that DESMO is designed to detect. Please provide an explicit systematic budget or clearly state in the abstract and conclusions that the quoted precision is statistical only.","section":"Section IV / Fig. 9"},{"comment":"The sensitivity projection is a Monte Carlo closure test: the simulated 60-day 'data' and the p/D templates are both produced by the same GEANT4 simulation. The fit is not circular in the narrow sense because the templates are statistically independent of the injected fraction, but the study cannot probe the fidelity of the detector response model. The manuscript should explicitly acknowledge this limitation and outline a calibration plan, for example beam-test measurements of the MSIS and RICH response or cross-validation with existing cosmic-ray isotope measurements, before the claimed precision is used to support scenario discrimination.","section":"Section IV, Fig. 8"},{"comment":"There is a direct numerical inconsistency in the target thickness: Section III A specifies a 14 mm thick BGO target, while Appendix B states that the final prototype has x = 9 mm. The target thickness directly controls the multiple-scattering angle and hence the D/p separation power, so the value actually used in the simulation must be identified and the discrepancy resolved.","section":"Section III A vs Appendix B"},{"comment":"The paper motivates the experiment by its ability to distinguish the three deuterium scenarios (a), (b), and (c), but the simulation only demonstrates recovery of a constant injected D/p = 2.5% across the energy range. To support the central physics claim, the authors should propagate the expected statistical (and, if available, systematic) uncertainties onto the scenario curves of Fig. 1 and show that the reconstructed D/p spectra remain distinguishable.","section":"Section II / IV, Fig. 1"},{"comment":"The analytic derivation of the average displacement uses a one-sided (half-Gaussian) integral, but the main text does not state that the measured displacement is folded or taken as an absolute value. If signed displacements are used, the mean displacement is zero and the relation d ∝ θ_MS is not valid. Please define the measured quantity precisely and confirm that the GEANT4 simulation implements the same folding procedure.","section":"Appendix B, Eqs. (B1)-(B2)"}],"minor_comments":[{"comment":"'Shade light' should be 'shed light' in both occurrences.","section":"Introduction and Section V"},{"comment":"The reference to 'Appendix III B' should be to 'Section III B'.","section":"Section III B"},{"comment":"The words 'maintainance' and 'mantainment' should be spelled 'maintenance'.","section":"Section III B"},{"comment":"'Inhelastic cross section' is a typo for 'inelastic cross section'.","section":"Section I"},{"comment":"The phrase 'x=9 mm (included in S)' is unclear; please explain what 'included in S' means and how the target thickness relates to the inter-module distance S.","section":"Appendix B"},{"comment":"'An hypothetical measurement' should be 'a hypothetical measurement'.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"This is a simulation-based design study. The central quantitative claim of 'few percent level precision' is not yet supported by an end-to-end systematic uncertainty assessment, and the internal 14 mm/9 mm target-thickness discrepancy must be corrected. If the journal is open to detector feasibility papers, a major revision addressing the major comments would be appropriate; otherwise the manuscript may be considered too preliminary for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nTwo things to know about the DESMO paper. First, the genuinely new part is the integration: multiple-scattering momentum measurement (from their earlier work and OPERA) with a RICH velocity measurement plus a BDT-based calibration, all packaged into a 90 kg balloon payload that claims D/p separation up to 100 GeV/nucl. Second, the central performance claim – a few percent precision on D/p – is statistical only, and the two dominant systematics (fragmentation in the residual atmosphere and in the detector) are deferred with estimates of 2–5% each, which is comparable to the 5–8% statistical error. That's the soft spot.\n\nThe paper does well on detail. The GEANT4 simulation covers realistic inclinations, dark counts, mirror defects, and a template-fit closure test with an injected D/p of 2.5%, which is recovered across six bins. The closure test is not circular: the templates are independent simulations and the fit is free. But that only proves the fitting machinery works, not that the templates are right. The authors are honest about this – they explicitly say the systematics are beyond the scope of the work – but the abstract and conclusion state 'few percent level precision' without that qualification, which overstates what is established.\n\nThere is also a concrete error in Appendix B. The integrals for the mean and mean-square of the positive-scattering angle are written without the factor 1/2 from the positive-deviation selection. As printed they give half the quoted values. The final conditional-moment formulas look plausible, but the derivation is mis-stated.\n\nThis is a serious concept paper, not a finished measurement. The physics motivation is solid, the design is concrete, and the simulation study is detailed enough to take seriously. The right next step is a systematic uncertainty budget for fragmentation and template validation, ideally with a test beam. That should not prevent the paper from going to peer review; it should be sent, with the expectation that the authors will need to address the systematics or scale back the headline claim.\n\nRecommendation: send it out. It earns referee time.","headline":"DESMO's integrated RICH+MSIS design is genuinely new, but the D/p precision claim is statistical-only and the deferred fragmentation systematics are comparable to the quoted errors.","tokens_in":25234,"tokens_out":5778,"would_cite":false,"duration_ms":54031,"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":"DESMO, a magnet-free balloon instrument, can measure the cosmic-ray deuteron-to-proton ratio at 20-200 GV with few-percent precision, giving a 60-day campaign enough power to decide whether deuterium is secondary or primary in origin.","keywords":["cosmic-ray deuterium","isotope identification","multiple Coulomb scattering","balloon-borne detector","RICH detector","template fit","cosmic-ray propagation","deuteron flux"],"falsifier":"A beam test of a DESMO prototype with tagged protons and deuterons of known momenta from 20 to 200 GeV/c, comparing the template-fit D/p fraction with the known mixture, would settle the claim; a bias larger than the quoted few percent, or a disagreement with AMS-02 in the overlapping 20-21 GV range beyond systematic errors, would falsify it.","tokens_in":1739,"feed_emoji":"🎈","tokens_out":1992,"duration_ms":81594,"temperature":0.7,"pith_summary":"The paper argues that a compact, magnet-free balloon instrument called DESMO can measure the cosmic-ray deuteron-to-proton ratio at high rigidities with enough precision to settle whether deuterium is purely secondary or has a primary component. Deuterium measurements currently stop at 21 GV, where the AMS-02 experiment saw more deuterons than secondary-production models predict, and the competing explanations differ mainly in how the spectrum behaves above that range. DESMO would extend the measurement to 200 GV (about 100 GeV per nucleon) by combining a multiple-scattering momentum measurement with a Cherenkov velocity measurement, and the authors' simulation shows a 60-day campaign recovering an injected 2.5% deuteron abundance with 5-8% statistical precision per energy bin. If the real detector behaves as simulated, this would discriminate between the secondary, mixed, and primary deuterium scenarios currently under debate.","feed_headline":"Magnet-free balloon detector measures cosmic deuterium to ~5 percent","feed_subtitle":"A 60-day balloon campaign would use scattering, not magnets, to test deuterium's primary or secondary origin.","key_machinery":"The Multiple Scattering Isotope Separator (MSIS) is the central device: eight identical PPT (Plane-Plane-Target) modules, each with two silicon pixel planes and a dense BGO scintillator target, measure the angular kicks a particle receives from multiple Coulomb scattering. Eight kicks per event are averaged into a displacement $d$ proportional to the scattering angle, hence inversely proportional to momentum, with resolution around 21-23%. A compact RICH with a silica aerogel radiator, spherical mirror, and SiPM readout measures particle velocity from the reconstructed Cherenkov ring, calibrated by a boosted decision tree to remove trajectory dependence; this yields energy-per-nucleon resolutions of 2-3% below 20 GeV/nucleon, degrading to about 50% at 100 GeV/nucleon. Combining momentum from scattering with velocity from the Cherenkov angle gives mass, separating deuterons from protons, and template fits of the inverse-displacement distribution extract the D/p fraction in each energy bin.","core_discovery":"The central claim is stated in Section V: DESMO can distinguish deuterons from protons and achieve a few-percent-level precision on the D/p fraction using a template-fit approach, without a magnetic spectrometer. Specifically, the authors simulate a 60-day exposure to a cosmic hydrogen flux with a deuterium abundance of 2.5% and a spectrum $\\propto E_k^{-2.7}$, reconstruct events with the MSIS and RICH detectors, and fit the distribution of inverse average displacement with proton and deuteron templates. They recover the injected D/p ratio across six energy bins spanning 10-100 GeV/nucleon, with statistical uncertainties of 5-8%, an acceptance of about 40 cm$^2$ sr, and a detector mass below 100 kg. They conclude that this is sufficient to significantly detect a deuterium abundance of 1-2% up to 100 GeV/nucleon and thereby to distinguish the three competing spectral scenarios.","pith_inferences":["The heavy reliance on simulated templates makes a beam-test calibration the natural next milestone; without it, the quoted precision rests on how well GEANT4's hadronic and scattering models reproduce real detector response at these energies.","Because the displacement resolution degrades as momentum rises and approaches the pixel spatial resolution, 200 GV is a practical ceiling for this geometry; pushing beyond would likely require finer position resolution or more stations rather than longer exposure.","If DESMO's highest bin aligns with the extrapolated AMS-02 trend, it would bridge the gap between AMS-02 and the SOKOL 1-4 TeV deuteron fraction, effectively testing whether the hard spectrum continues to TeV energies.","A precise D/p measurement would also sharpen predictions for the cosmic-ray antideuteron background in dark-matter searches, since deuteron propagation is the main physics proxy for antideuteron propagation; the paper notes this connection but does not quantify it."],"forward_implications":["A 60-day balloon campaign, split into two 30-day flights, would produce six energy bins of the D/p ratio covering about 10-100 GeV/nucleon with 5-8% statistical precision per bin if the simulated acceptance of about 40 cm$^2$ sr holds.","The measured spectral shape would discriminate among secondary-only production, secondary production with different light- and heavy-primary acceleration, and a primary-like power-law continuation of the AMS-02 trend.","A significant detection at the 1-2% deuterium abundance level up to 100 GeV/nucleon would extend the deuteron excess observed by AMS-02 and directly confront the high-energy estimates based on the SOKOL measurement.","The same MSIS technique could be adapted to other isotope ratios and, with an added calorimeter module, to antideuteron searches, both of which the paper sketches as reachable design variants."],"supporting_citations":[{"why":"Validates the multiple-scattering method for Z=1 isotopic separation that the DESMO design is built on.","marker":"[18]"},{"why":"Supplies the established technique of momentum measurement via multiple Coulomb scattering that replaces a magnetic spectrometer.","marker":"[19]"},{"why":"Provides the AMS-02 2-21 GV deuteron data whose high-energy excess motivates the experiment and defines the scenarios to be tested.","marker":"[6]"},{"why":"Shows the AMS-02 deuteron flux can be explained as secondary, the scenario DESMO must distinguish from primary production.","marker":"[11]"},{"why":"Proposes different acceleration mechanisms for light primaries as an explanation of the deuteron excess, another scenario DESMO would test.","marker":"[4]"},{"why":"Calculates an upper limit on a primary deuteron component that DESMO aims to confirm or exclude.","marker":"[15]"},{"why":"The GEANT4 simulation toolkit used to model detector response and assess DESMO performance.","marker":"[20]"},{"why":"Used to estimate the 2-5% systematic uncertainty from deuteron fragmentation in the residual atmosphere.","marker":"[41]"}],"fun_headline_variants":["Magnet-free detector simulates precise cosmic deuterium measurement","Magnet-free detector could test deuterium's origin up to 100 GeV","Scattering-based detector to measure cosmic deuterium without magnets","New detector design probes deuterium excess up to 100 GeV/nucl","Balloon detector uses multiple scattering to measure deuterium precisely"],"cache_read_input_tokens":27392,"weakest_assumption_plain":"The load-bearing premise is that the proton and deuteron templates used to fit the scattering distributions are accurate, including effects the paper did not simulate, such as fragmentation of deuterium and heavier nuclei in the atmosphere and detector, because the quoted few-percent precision on D/p comes directly from fitting those templates.","fun_headline_variants_meta":{"raw":{"variants":["Magnet-free detector simulates precise cosmic deuterium measurement","Magnet-free detector could test deuterium's origin up to 100 GeV","Scattering-based detector to measure cosmic deuterium without magnets","New detector design probes deuterium excess up to 100 GeV/nucl","Balloon detector uses multiple scattering to measure deuterium precisely"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000728,"raw_usage":{"total_tokens":3273,"prompt_tokens":967,"completion_tokens":2306,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":2232}},"tokens_in":583,"tokens_out":2306,"duration_ms":14351,"temperature":1.0,"reasoning_tokens":2232,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:42:48.484222+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A beam test of a DESMO prototype with tagged protons and deuterons of known momenta from 20 to 200 GeV/c, comparing the template-fit D/p fraction with the known mixture, would settle the claim; a bias larger than the quoted few percent, or a disagreement with AMS-02 in the overlapping 20-21 GV range beyond systematic errors, would falsify it.","supporting_citations":[{"cited_title":"Cosmic-ray deuteron excess from a pri- mary component, 2024","cited_arxiv_id":null,"evidence_quote":"Validates the multiple-scattering method for Z=1 isotopic separation that the DESMO design is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the established technique of momentum measurement via multiple Coulomb scattering that replaces a magnetic spectrometer."},{"cited_title":"Moskalenko, and Michael Unger","cited_arxiv_id":null,"evidence_quote":"Provides the AMS-02 2-21 GV deuteron data whose high-energy excess motivates the experiment and defines the scenarios to be tested."},{"cited_title":"Papini, S","cited_arxiv_id":null,"evidence_quote":"Shows the AMS-02 deuteron flux can be explained as secondary, the scenario DESMO must distinguish from primary production."},{"cited_title":"horizontal","cited_arxiv_id":null,"evidence_quote":"Proposes different acceleration mechanisms for light primaries as an explanation of the deuteron excess, another scenario DESMO would test."},{"cited_title":"Donato, N","cited_arxiv_id":null,"evidence_quote":"Calculates an upper limit on a primary deuteron component that DESMO aims to confirm or exclude."},{"cited_title":"Aguilar, L","cited_arxiv_id":null,"evidence_quote":"The GEANT4 simulation toolkit used to model detector response and assess DESMO performance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Used to estimate the 2-5% systematic uncertainty from deuteron fragmentation in the residual atmosphere."}],"review_version":1}