Neutrino Radiation Transport and Other Topics in High Energy Density Astrophysics

Abstract

Neutron star mergers and the collapse of massive stars result in some of the universe’s most violet explosions. However, the detailed mechanisms behind all of these astrophysical explosions remain elusive. Their strongly nonlinear and complicated nature makes them difficult and expensive to simulate, and the properties of matter in these extreme conditions are poorly constrained. I use a variety of computational tools to understand the detailed mechanisms behind both types of events. I describe my relativistic time-independent multidimensional Monte Carlo neutrino radiation transport code Sedonu that provides an accurate account of the neutrino radiation fields and the interaction with neutrinos and background fluid. Though Sedonu calculations are time-independent, I demonstrate their utility in dynamical general relativistic variable Eddington tensor radiation hydrodynamics simulations. I apply Sedonu to simulations of accretion disks following neutron star mergers to demonstrate that more realistic disk cooling and neutrino-driven mass ejection rates are larger than is predicted using approximate transport methods. I also reinforce that neutrino pair annihilation from these disk configurations is unlikely to be able to energize a gamma-ray burst jet. I subject Sedonu to the first thorough comparison of Boltzmann neutrino radiation transport methods in multiple spatial dimensions in the context of core-collapse supernovae. The comparisons with the other highly accurate discrete ordinates-based transport scheme show remarkably similar results, verifying the accuracy of both methods and underscoring the importance of numerical fidelity. I perform the first broad parameter study on how different descriptions of dense nuclear matter and star rotation rates influence the dynamics of, and hence gravitational waves from, the bounce and early post-bounce phase of rapidly rotating core collapse supernovae. Using the results of 1824 two-dimensional general relativistic core-collapse simulations, I demonstrate that the equation of state is unlikely to be constrained by LIGO observations. I show that the effect of the equation of state on the gravitational wave frequency can be described by a single universal relation. Finally, I use results of three-dimensional general relativistic magnetohydrodynamics simulations of rapidly rotating core collapse to demonstrate that the polar magnetic structures that form are destroyed by a magnetohydrodynamic kink instability.