LUNAR: a Monte Carlo generator for bound-nucleon decay in liquid argon
Abstract
The search for nucleon decay in liquid-argon time-projection chambers requires a quantitative description of how the bound nuclear environment reshapes the decay-product kinematics. We present LUNAR, a fast, openly available Monte Carlo generator dedicated to two-body decays of protons and neutrons bound in argon-40, the target nucleus of the DUNE far detector. The parent nucleon is drawn from a selectable nuclear ground state -- ten momentum distributions ranging from mean-field Fermi gases to argon spectral functions -- and bound off the mass shell by one of three removal-energy prescriptions, including the momentum-dependent optical potential of Juszczak \textit{et al}. The two-body decay is performed off-shell and boosted to the laboratory frame, and the daughter meson is then propagated out of the nucleus by a semi-classical intranuclear cascade with an optional formation zone for the freshly produced meson. We use the generator to separate the distinct roles of Fermi motion and binding in shaping the observable meson spectrum, to quantify final-state interactions channel by channel, and to translate present Super-Kamiokande limits into expected DUNE event yields for the full set of standard decay modes. Final-state interactions leave the supersymmetry-favored p→ K+ν̅ signal essentially intact while roughly halving the pion, η, and antikaon rates -- an effect that dominates over the ±10% spread induced by the choice of nuclear model. The code is released to the community as a lightweight, extensible tool for signal efficiency and systematics studies.
Versions
Related
- Wigner distribution modeling of nuclear state in NuWro Monte Carlo event generator
- CGC-py: A Monte Carlo Event Generator for Gluon Saturation Physics
- Uncertainty in AI driven Monte Carlo simulations
- Unbiased Diffusion Monte Carlo for non local operators
- A control variates technique for the null-collision Monte Carlo algorithms in radiative transfer
BibTeX
@misc{nowak2026lunar,
title = {{LUNAR: a Monte Carlo generator for bound-nucleon decay in liquid argon}},
author = {Jaroslaw Nowak},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2606.30872},
eprint = {2606.30872},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2606.30872},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - LUNAR: a Monte Carlo generator for bound-nucleon decay in liquid argon
AU - Jaroslaw Nowak
PY - 2026
DA - 2026/09/11
PB - arXiv
DO - 10.48550/arxiv.2606.30872
UR - https://arxiv.org/abs/2606.30872
AB - The search for nucleon decay in liquid-argon time-projection chambers requires a quantitative description of how the bound nuclear environment reshapes the decay-product kinematics. We present LUNAR, a fast, openly available Monte Carlo generator dedicated to two-body decays of protons and neutrons bound in argon-40, the target nucleus of the DUNE far detector. The parent nucleon is drawn from a selectable nuclear ground state -- ten momentum distributions ranging from mean-field Fermi gases to argon spectral functions -- and bound off the mass shell by one of three removal-energy prescriptions, including the momentum-dependent optical potential of Juszczak \textit{et al}. The two-body decay is performed off-shell and boosted to the laboratory frame, and the daughter meson is then propagated out of the nucleus by a semi-classical intranuclear cascade with an optional formation zone for the freshly produced meson. We use the generator to separate the distinct roles of Fermi motion and binding in shaping the observable meson spectrum, to quantify final-state interactions channel by channel, and to translate present Super-Kamiokande limits into expected DUNE event yields for the full set of standard decay modes. Final-state interactions leave the supersymmetry-favored p→ K^+ν̅ signal essentially intact while roughly halving the pion, η, and antikaon rates -- an effect that dominates over the ±10% spread induced by the choice of nuclear model. The code is released to the community as a lightweight, extensible tool for signal efficiency and systematics studies.
ER -