Differentiable nuclear deexcitation simulation for low energy neutrino physics

Pablo Barham Alzás,Radi Radev
2024-03-30
Abstract:Neutrino-nucleus interactions play an important role in present and future neutrino experiments. The accurate simulation of these interactions at low energies ($<$100 MeV) is crucial for the detection and study of supernova, solar and atmospheric neutrinos. In particular, the reconstruction of the incoming neutrino properties depends on the ability to measure the products from the deexcitation of the final state nucleus after the initial neutrino-nucleus scattering reaction. A realistic nuclear deexcitation model that can correctly manage the theoretical uncertainties in the process is key to determine the response of a detector to low energy neutrinos and estimate the overall systematic uncertainties that affect it.
High Energy Physics - Experiment
What problem does this paper attempt to address?
The problem addressed in this paper is to improve the simulation of nuclear de-excitation in low-energy neutrino physics. Simulation of the interaction between low-energy neutrinos (below 100 MeV) and atomic nuclei is crucial for detecting and studying supernovae, the sun, and atmospheric neutrinos. Existing simulation methods typically divide the process into two steps: neutrino-nucleon scattering and de-excitation of the final-state nucleus. However, this separate simulation approach neglects the possibility that the de-excitation process may depend on the initial-state properties of the neutrinos. The paper proposes using automatic differentiation frameworks (such as PyTorch or JAX) to build a differentiable simulation model for nuclear de-excitation, based on a simplified version of the Hauser-Feshbach statistical emission model. This allows for simulation optimization to match observed data, estimation of parameter distributions on the physical outputs, and parameter inference. The authors implement a simplified differentiable simulation of nuclear de-excitation as a proof of concept, demonstrating its vectorized execution and gradient checking on GPUs. The results of the paper show that this approach is feasible and has fast and robust sampling and gradient estimation. Although there are some instabilities, especially with small sample sizes, the variance of the estimation loss can be reduced by increasing the sample size. Future work includes improving the Hauser-Feshbach model, enhancing the stability and speed of the gradient descent algorithm, and ultimately integrating this model into a fully differentiable simulation chain from generation to reconstruction.