Monte Carlo Neutrino Closures in 3D GRMHD Simulations of Core-Collapse Supernovae and Neutron Star Mergers
Three dimensional (3D) core collapse supernova and neutron star merger simulations suffer from approximate treatments of neutrino transport that cripple their reliability and realism. I propose to use Blue Waters to carry out direct Monte Carlo simulations of the neutrino transport problem to calculate the true closure relation that will perfect the two-moment neutrino transport in 3D GRMHD simulations of these systems. These 3D simulations will be the most realistic to date and will be the first to have concrete error estimates on the accuracy of the neutrino transport method. I will also provide the first open-source benchmark with which all neutrino transport methods of any number of dimensions can compare.