Saved in:
Bibliographic Details
Main Authors: Alzás, Pablo Barham, Radev, Radi
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2404.00180
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910392493539328
author Alzás, Pablo Barham
Radev, Radi
author_facet Alzás, Pablo Barham
Radev, Radi
contents 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. Automatic differentiation frameworks like PyTorch or JAX, widely used in the field of Machine Learning, provide us the tools to compute exact gradients of the simulation outputs with respect to the model parameters. Such differentiable simulators can be applied in simulator tuning to match observed data and forward modeling to efficiently infer the impact of parameter distributions to the physics output, or for parameter inference acting as a likelihood estimator. This paradigm can pave the way to a fully differentiable analysis of the whole simulation chain or direct integration with preexisting Machine Learning tools. As a proof of concept, in this work we implement a fully differentiable nuclear deexcitation simulation based on a simplified version of the Hauser-Feshbach statistical emission model.
format Preprint
id arxiv_https___arxiv_org_abs_2404_00180
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Differentiable nuclear deexcitation simulation for low energy neutrino physics
Alzás, Pablo Barham
Radev, Radi
High Energy Physics - Experiment
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. Automatic differentiation frameworks like PyTorch or JAX, widely used in the field of Machine Learning, provide us the tools to compute exact gradients of the simulation outputs with respect to the model parameters. Such differentiable simulators can be applied in simulator tuning to match observed data and forward modeling to efficiently infer the impact of parameter distributions to the physics output, or for parameter inference acting as a likelihood estimator. This paradigm can pave the way to a fully differentiable analysis of the whole simulation chain or direct integration with preexisting Machine Learning tools. As a proof of concept, in this work we implement a fully differentiable nuclear deexcitation simulation based on a simplified version of the Hauser-Feshbach statistical emission model.
title Differentiable nuclear deexcitation simulation for low energy neutrino physics
topic High Energy Physics - Experiment
url https://arxiv.org/abs/2404.00180