Neutrino interactions are essential for an accurate understanding of the binary neutron star merger process. In this article, we extend the code infrastructure of the well-established numerical-relativity code BAM that until recently neglected neutrino-driven interactions. In fact, while previous work allowed already the usage of nuclear-tabulated equations of state and employing a neutrino leakage scheme, we are moving forward by implementing a first-order multipolar radiation transport scheme (M1) for the advection of neutrinos. After testing our implementation on a set of standard scenarios, we apply it to the evolution of four low-mass binary systems, and we perform an analysis of ejecta properties. We also show that our new ejecta analysis infrastructure is able to provide numerical relativity-informed inputs for the codes POSSIS and Skynet, for the computation of kilonova lightcurves and nucleosynthesis yields, respectively.
@article{arxiv.2307.04572,
title = {M1 neutrino transport within the numerical-relativistic code BAM with application to low mass binary neutron star mergers},
author = {Federico Schianchi and Henrique Gieg and Vsevolod Nedora and Anna Neuweiler and Maximiliano Ujevic and Mattia Bulla and Tim Dietrich},
journal= {arXiv preprint arXiv:2307.04572},
year = {2023}
}