{"id":"5d298103-c42e-451a-a6b0-294dda786056","arxiv_id":"1908.06902","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"Using the TWIST spectrometer, the authors measured proton and deuteron momentum spectra from muon capture on aluminum, obtaining partial yields of 0.0322 and 0.0122 per capture above 80 and 130 MeV/c, respectively.","lead":"This paper reports the first measurement of proton and deuteron momentum spectra from negative muon capture on aluminum, using the TWIST tracking spectrometer. The results provide design-critical input for the Mu2e and COMET lepton-flavor-violation experiments and benchmark nuclear models used to simulate muon capture.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Energy-loss model in the detector response matrix is the weakest link: GEANT3.21 shows a 3–4% data–MC discrepancy (Fig. 5e), and the assigned systematic rests on a single material-density variation rather than an independent validation.","rationale":"The reader and I have converged on the GEANT3.21 energy-loss model as the most load-bearing assumption. It is not a hidden flaw: the paper explicitly notes the 4% data-MC discrepancy in Sec. V D and attempts a bounding systematic. The energy-loss contribution to the partial yields is nevertheless smaller than the crosstalk systematic for protons, and the unfolding regularization and MC closure tests are plausibly adequate. The extrapolated total yields also carry an explicit model-dependence uncertainty from GEANT4-precompound/MECO averaging, which the authors separate from the main systematic. I therefore do not see a basis to change the ACCEPT verdict. The appropriate action is a targeted verification: replacing the GEANT3.21 stopping power with an independent model and re-running the unfolding would settle whether the quoted energy-loss uncertainty is realistic. If the check failed, the verdict would need to move toward CONDITIONAL or REJECT for the precision claims, but based on the current evidence I would not move the verdict pre-emptively.","tokens_in":17893,"tokens_out":10225,"duration_ms":114817,"concrete_test":"Regenerate the detector response matrix with an independent energy-loss implementation for protons and deuterons in all TWIST materials (e.g., PSTAR/SRIM range tables or GEANT4 ion energy loss) instead of the GEANT3.21 default, while keeping the geometry and all other response ingredients unchanged. Rerun the complete unfolding on the same data and compare the partial yields above 80 MeV/c (protons) and 130 MeV/c (deuterons) with Eqs. (11) and (12). If the shift from the baseline exceeds the energy-loss uncertainty row of Table I (0.00025/0.00032), the assigned systematic is underestimated and the central precision claim should be reconsidered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central measured yields and unfolded spectra inherit the accuracy of the detector response matrix, which is generated with out-of-the-box GEANT3.21 energy loss for protons/deuterons. Section V D documents that this model does not fully describe the data: for long contained tracks (Fig. 5e) the reconstructed momentum peak is 92.4 MeV/c in data versus 88.7 MeV/c in MC, about a 4% discrepancy. The systematic is estimated by changing the DME gas density by 20% and re-running the full unfolding, but this is a proxy that may not span the true model uncertainty in proton/deuteron stopping power, especially for low-energy particles near the 80/130 MeV/c thresholds where target and chamber energy loss dominate the response. If the true energy loss differs from GEANT3.21 beyond the few-percent level, the unfolded yields in Eqs. (11) and (12) could shift by more than the quoted energy-loss terms (0.00025 for protons, 0.00032 for deuterons), and the claimed precision would not hold. The extrapolated total yields would also be affected because the correction for the unmeasured sub-threshold region is normalized to the measured spectrum just above threshold.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports the first momentum spectra of protons and deuterons emitted after nuclear muon capture on aluminum, measured with a magnetic tracking spectrometer (the TWIST detector) rather than with calorimetric or emulsion techniques. The analysis selects a sample of 18M muon stops in a thin aluminum target, reconstructs positively charged tracks in downstream drift chambers, separates protons from deuterons using range-versus-momentum information, and determines the absolute normalization from DIO electrons plus the known muon-capture-to-decay ratio. The detector response is obtained from a GEANT3.21 simulation, and the true momentum spectra are extracted with a regularized maximum-likelihood unfolding in which each spectrum is represented by a cubic-spline deviation from an exponential-in-kinetic-energy form. The paper quotes partial yields of 0.0322 +/- 0.0007(stat) +/- 0.0022(syst) per capture for protons above 80 MeV/c and 0.0122 +/- 0.0009(stat) +/- 0.0006(syst) for deuterons above 130 MeV/c, plus extrapolated total yields of 0.045 +/- 0.003 and 0.018 +/- 0.002 per capture. The results are compared with the MECO parameterization, with GEANT4 precompound-model predictions, and with earlier activation measurements.","tokens_in":18141,"tokens_out":10527,"duration_ms":119333,"significance":"If the central results hold, this is a substantial new measurement: it is the first tracking-spectrometer measurement of muon-capture charged-particle emission, and it provides the most precise per-capture yields of protons and deuterons from muon capture on aluminum to date. The measurement is directly relevant to background estimates for the Mu2e and COMET lepton-flavor-violation searches and to the validation of nuclear-response models relevant to neutrino physics. The paper is careful in its treatment of the unfolding problem: the method is validated on statistically independent Monte Carlo samples, the statistical uncertainty combines data and simulation statistics via pseudo-data, and many systematic effects (crosstalk, muon stopping position, stopping-distribution cut, energy loss, heavy-particle background, and method bias) are evaluated explicitly.","major_comments":[],"minor_comments":[{"comment":"The energy-loss systematic is estimated from a 20% increase in the DME gas density, but the manuscript does not explicitly state how much this density variation shifts the reconstructed-momentum peak relative to the observed 4% data-MC discrepancy at 92.4 vs 88.7 MeV/c. Please add one or two sentences (or a small quantitative statement) showing that the density variation brackets the observed peak shift; this would make the coverage of the assigned uncertainty in Table I more transparent without requiring a new analysis.","section":"Sec. V.D, Fig. 5(e)"},{"comment":"The extrapolation to total yields normalizes GEANT4 precompound and MECO predictions to the measured spectrum just above threshold and averages the two results. It would be helpful to state explicitly that both of these are models for the sub-threshold region and that the quoted extrapolation uncertainty does not cover a common-mode error in which both models share the same incorrect low-energy shape; the current wording is clear about the procedure but could be more explicit about this residual model dependence.","section":"Sec. VI, Eqs. (13)-(14)"},{"comment":"The sentence 'the parameterization shape uncertainty is not representable in terms of the fixed set of fit parameters, and is the only contribution that is not included' is ambiguous. Please specify whether this uncertainty is included in the total-uncertainty bands of Figs. 10-11 and in the yields of Eqs. (11)-(14), or explain how it would be estimated; as written, the reader cannot tell which final quantities include it.","section":"Appendix, Table II"},{"comment":"There are several typographical errors that should be corrected: 'partices' in Sec. I, 'withing' in Sec. VI, 'uncertainy' in Sec. V.G, and 'pararameterization' in Sec. V.F. These do not affect the physics but should be cleaned up before publication.","section":"Throughout"}],"recommendation":"minor_revision","confidential_remarks":"This is a solid experimental paper. The central issue raised by the external stress test, the GEANT3 energy-loss model, is handled transparently with a dedicated systematic; my request to quantify the coverage of the observed 4% peak discrepancy is a local clarification rather than a fatal objection. The self-citation to the unfolding method (ref. [44]) is appropriate and does not create circularity, because the unfolding technique is validated on Monte Carlo and is independent of the physics result. The paper is within the scope of the journal and, with the minor clarifications above, is suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWhat you should know first: this is a genuinely new measurement—the first tracking-spectrometer study of charged particles from muon capture—and it gives the most precise proton and deuteron yields per capture on aluminum. The partial yields above 80 and 130 MeV/c are 0.0322 ± 0.0022(syst) ± 0.0007(stat) and 0.0122 ± 0.0006(syst) ± 0.0009(stat); extrapolated totals are 0.045 ± 0.003 and 0.018 ± 0.002. These numbers are direct inputs for Mu2e/COMET background simulations and useful discriminators for nuclear response models.\n\nThe analysis is careful. The DIO normalization is a clean way to convert to per-capture yields, canceling common systematics. The unfolding method from the companion paper is a sensible solution to low statistics and sharply falling spectra; regularizing only the spline deviation from an exponential tail is more principled than a blind cutoff. The systematics table is thorough—crosstalk, stopping position, stopping distribution, energy loss, heavy-particle background, and method bias are each evaluated by re-running the full unfolding chain. The crosstalk study using un-flagging fractions is especially nice.\n\nThe soft spot is the detector response model for protons and deuterons. It relies on out-of-the-box GEANT3.21 energy loss, and the authors are honest that long contained tracks show a 3–4% shift in the reconstructed momentum peak (Fig. 5e). The assigned systematic uses a 20% DME density change as a proxy, not an independent stopping-power validation. That is the weakest link. For the partial yields, the energy-loss term is small (0.00025 proton, 0.00032 deuteron), and the authors did try several variations, most of which broke agreement elsewhere, so I don't think the central numbers are in danger. But for spectra near threshold, where energy loss dominates, the model uncertainty could be larger than the quoted term. The extrapolated totals are explicitly model-dependent; the separate extrapolation uncertainty handles this, but it should be read as a shape assumption, not a measurement.\n\nThe citation to the companion paper for the unfolding method is not circular; it's a statistical technique validated on MC.\n\nBottom line: solid, significant, carefully reported measurement. It deserves a serious referee. I'd send it to review, with the request to expand the energy-loss systematic discussion and show the momentum-dependent systematic contributions for the lowest bins.\n\nSend it.","headline":"First tracking-spectrometer spectra of muon-capture protons/deuterons on Al; careful analysis, defensible if proxy-based energy-loss systematic; deserves peer review.","tokens_in":18717,"tokens_out":3291,"would_cite":true,"duration_ms":32725,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.30.Mr","29.40.Gx","23.40.-s"],"model":"deepseek-v4-flash","headline":"The paper reports the first tracking-spectrometer measurement of the momentum spectra of protons and deuterons from muon capture on aluminum, with partial yields per capture of $0.0322 \\pm 0.0022$ above 80 MeV/c and $0.0122 \\pm 0.0006$…","keywords":["muon capture","aluminum-27","proton momentum spectrum","deuteron momentum spectrum","spectrum unfolding","tracking spectrometer","lepton flavor violation backgrounds","nuclear response"],"falsifier":"Measure the range-momentum relation for protons and deuterons in dimethyl ether at 80\\textendash 250 MeV/c in a dedicated beam test, rerun the response matrix and unfolding with the measured energy loss, and see whether the partial yields move by more than the quoted systematic uncertainties; if they do, the central yields are biased.","tokens_in":17676,"feed_emoji":"⚛️","tokens_out":10848,"duration_ms":108530,"temperature":0.7,"pith_summary":"Negative muons stopped in a thin aluminum foil are captured by the nucleus about 60% of the time, releasing protons, deuterons, and heavier fragments. This paper reports the first measurement of the proton and deuteron momentum spectra from those captures using a magnetic tracking spectrometer, which reconstructs the curved trajectories rather than only recording deposited energy. The central results are partial yields per capture of $0.0322 \\pm 0.0022$ for protons above 80 MeV/c and $0.0122 \\pm 0.0006$ for deuterons above 130 MeV/c, with extrapolated total yields of $0.045 \\pm 0.003$ and $0.018 \\pm 0.002$ per capture. These are the most precise values to date, and they matter for two practical reasons: experiments that stop muons in aluminum need to know how many highly ionizing particles their trackers will see, and neutrino-nucleus calculations use the same nuclear response probed by muon capture.","feed_headline":"First tracking-spectrometer spectra of muon capture on aluminum","feed_subtitle":"The measured yields give lepton-flavor-violation searches a benchmark for the ionizing background they must survive.","key_machinery":"The argument rests on a detector-response matrix $R_{\\eta,ij}$ built from a detailed simulation: for a true proton or deuteron of momentum $p_j$, it gives the expected number of reconstructed tracks in each bin of the measured range-versus-momentum and momentum histograms. Because protons and deuterons with the same momentum stop at different depths, the track-range variable separates the two species. Absolute normalization per capture is set by counting decay-in-orbit electrons, whose spectrum and detector response are already known precisely, and scaling by the measured capture-to-decay ratio. The spectra are then extracted by a regularized maximum-likelihood unfolding, a statistical inversion that removes detector smearing; its distinctive feature is that below an unfolding cutoff the spectrum is a free cubic-spline shape, while above the cutoff it continues as an exponential in kinetic energy whose slope is fitted, with the regularization applied only to the spline part. That design removes the unphysical rising high-momentum tails that standard Tikhonov or maximum-entropy regularizations produced.","core_discovery":"The central claim is that the momentum distributions of protons and deuterons emitted after nuclear muon capture on aluminum are now measured, not inferred. Using a magnetic tracking spectrometer, the analysis finds sharply falling spectra that are approximately exponential in kinetic energy, with about 70% of the proton yield below 80 MeV/c and about 70% of the deuteron yield below 130 MeV/c. The observed partial yields above those thresholds are $0.0322 \\pm 0.0007(\\mathrm{stat}) \\pm 0.0022(\\mathrm{syst})$ per capture for protons and $0.0122 \\pm 0.0009(\\mathrm{stat}) \\pm 0.0006(\\mathrm{syst})$ for deuterons; extrapolating the measured shapes to zero momentum gives total yields $0.045 \\pm 0.003$ and $0.018 \\pm 0.002$ per capture. These totals are the most precise obtained to date, sit close to the theoretical values of 0.040 and 0.012, and bring the sum of proton and deuteron emission on aluminum into line with older radioisotope-activation sums once the unmeasured low-momentum parts are included.","pith_inferences":["A dedicated beam-test measurement of low-energy proton and deuteron energy loss in the drift-chamber gas would test whether the observed 3\\textendash 4% momentum-peak shift is an energy-loss-model bias; if it is, the deuteron yields, whose dominant systematic is energy loss, could move by roughly the quoted 0.00032.","The unfolding recipe, free spline below a cutoff plus fitted exponential above it with regularization applied only to the spline part, should transfer to other steeply falling spectra with sparse high-energy tails, such as neutrino-induced nucleon spectra, where standard regularizations produce artificial rising tails.","Reinterpreting the same data with a heavy-particle background model other than the two tested, or adding time-of-flight particle identification, would provide a sharper test of the deuteron yield, since the deuteron channel is the one most sensitive to the triton and alpha background shape."],"forward_implications":["The differential spectra, provided in bin-ready form, can be fed directly into event-generator simulations for any experiment that stops negative muons in aluminum, replacing the older ad hoc parameterization that had overestimated the proton rate.","The total yields of $0.045\\pm0.003$ and $0.018\\pm0.002$ per capture reconcile the charged-particle emission sum on aluminum with radioisotope-activation totals, since activation cannot separately count the $p+n$ and deuteron channels.","The measured proton spectrum is selective between nuclear models: it favors the precompound description, in which energy is shared before the nucleus thermalizes, over the cascade description for the high-energy tail, making these data a benchmark for muon-capture and neutrino-nucleus response models.","With better than 10% precision for protons from 100 to 190 MeV/c and better than 20% for deuterons from 145 to 250 MeV/c, the high-energy parts of both spectra are now fixed well enough for detector-rate studies in lepton-flavor-violation searches.","The reported correlation of $-0.25$ between the visible proton and deuteron yields means the two yields are not independent constraints; any future model must match both simultaneously."],"supporting_citations":[{"why":"Supplies the precisely known decay-in-orbit electron spectrum and detector response used for absolute per-capture normalization.","marker":"[11]"},{"why":"Validates the detector description, muon-beam tuning, and reconstruction software on which the analysis relies.","marker":"[10]"},{"why":"Provides the simulation's energy-loss and multiple-scattering model that generates the detector response matrix and defines the dominant systematic.","marker":"[28]"},{"why":"Introduces the regularized unfolding with exponential continuation that prevents the unphysical high-momentum tail.","marker":"[44]"},{"why":"Earlier silicon active-target measurement whose proton spectrum shape and yield are the main prior comparison for low-energy protons.","marker":"[18]"},{"why":"Supplies the parameterization used in lepton-flavor-violation design studies; the paper compares its spectra with this and shows it bounds the data.","marker":"[23]"},{"why":"Gives the theoretical proton yield of 0.040 and deuteron yield of 0.012 per capture against which the extrapolated totals are compared.","marker":"[13]"},{"why":"Provides the measured muon capture-to-decay branching ratio used to convert the decay-in-orbit normalization to per-capture yields.","marker":"[35]"},{"why":"Earlier calorimetric measurement of energetic charged particles from muon capture; comparison shows agreement in slope but a different normalization above 40 MeV.","marker":"[20]"}],"fun_headline_variants":["First tracking spectrometer muon capture spectra on aluminum","Precise proton and deuteron yields from muon capture on Al","Muon capture on Al measured with tracking spectrometer for LFV","Tracking study yields precise muon capture spectra on aluminum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on the assumption that the simulation's energy-loss model for protons, deuterons, and heavier ions reproduces reality to within the assigned few-percent systematic, because the only direct data check, a 3\\textendash 4% shift in the reconstructed momentum peak of long-range contained tracks, does not uniquely pin down the cause.","fun_headline_variants_meta":{"raw":{"variants":["First tracking spectrometer muon capture spectra on aluminum","Precise proton and deuteron yields from muon capture on Al","Muon capture on Al measured with tracking spectrometer for LFV","Tracking study yields precise muon capture spectra on aluminum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000495,"raw_usage":{"total_tokens":2496,"prompt_tokens":1082,"completion_tokens":1414,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":698,"completion_tokens_details":{"reasoning_tokens":1346}},"tokens_in":698,"tokens_out":1414,"duration_ms":12568,"temperature":1.0,"reasoning_tokens":1346,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:31:29.099881+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the range-momentum relation for protons and deuterons in dimethyl ether at 80\\textendash 250 MeV/c in a dedicated beam test, rerun the response matrix and unfolding with the measured energy loss, and see whether the partial yields move by more than the quoted systematic uncertainties; if they do, the central yields are biased.","supporting_citations":[{"cited_title":"MECO spectrum","cited_arxiv_id":null,"evidence_quote":"Supplies the precisely known decay-in-orbit electron spectrum and detector response used for absolute per-capture normalization."},{"cited_title":"un-flagging","cited_arxiv_id":null,"evidence_quote":"Validates the detector description, muon-beam tuning, and reconstruction software on which the analysis relies."},{"cited_title":"Heusser and T","cited_arxiv_id":null,"evidence_quote":"Introduces the regularized unfolding with exponential continuation that prevents the unphysical high-momentum tail."},{"cited_title":"10: Yield of protons per muon capture vs momentum","cited_arxiv_id":null,"evidence_quote":"Earlier silicon active-target measurement whose proton spectrum shape and yield are the main prior comparison for low-energy protons."},{"cited_title":"Haff and T","cited_arxiv_id":null,"evidence_quote":"Supplies the parameterization used in lepton-flavor-violation design studies; the paper compares its spectra with this and shows it bounds the data."},{"cited_title":"0009(stat) ± 0","cited_arxiv_id":null,"evidence_quote":"Gives the theoretical proton yield of 0.040 and deuteron yield of 0.012 per capture against which the extrapolated totals are compared."},{"cited_title":"Lifshitz and P","cited_arxiv_id":null,"evidence_quote":"Provides the measured muon capture-to-decay branching ratio used to convert the decay-in-orbit normalization to per-capture yields."},{"cited_title":"11: Yield of deuterons per muon capture vs momentum","cited_arxiv_id":null,"evidence_quote":"Earlier calorimetric measurement of energetic charged particles from muon capture; comparison shows agreement in slope but a different normalization above 40 MeV."}],"review_version":1}