Long-term Monte Carlo-based neutrino-radiation hydrodynamics simulations for a black hole-torus system

Kyohei Kawaguchi,Sho Fujibayashi,Masaru Shibata
2024-11-29
Abstract:We present our new general relativistic Monte Carlo (MC)-based neutrino radiation hydrodynamics code designed to solve axisymmetric systems with several improvements. The main improvements are as follows: (i) the development of an extended version of the implicit MC method for multi-species radiation fields; (ii) modeling of neutrino pair process rates based on a new numerically efficient and asymptotically correct fitting function for the kernel function; (iii) the implementation of new numerical limiters on the radiation-matter interaction to ensure a stable and physically correct evolution of the system. We apply our code to a black hole (BH)-torus system with a BH mass of $3\,M_\odot$, BH dimmensionless spin of 0.8, and a torus mass of $0.1\,M_\odot$, which mimics a post-merger remnant of a binary neutron star merger in the case that the massive neutron star collapses to a BH within a short time scale ($\sim10\,{\rm ms}$). We follow the evolution of the BH-torus system up to more than $1\,{\rm s}$ with our MC-based radiation viscous-hydrodynamics code that dynamically takes into account non-thermal pair annihilation. We find that the system evolution and the various key quantities, such as neutrino luminosity, ejecta mass, torus $Y_e$, and pair annihilation luminosity, are broadly in agreement with the results of the previous studies. We also find that the $\nu_e{\bar \nu}_e$ pair annihilation can launch a relativistic outflow for a time scale of $\sim 0.1\,{\rm s}$, and it can be energetic enough to explain some of short-hard gamma-ray bursts and the precursors. Finally, we calculate the indicators of the fast flavor instability directly from the obtained neutrino distribution functions, which indicate that the instability can occur particularly near the equatorial region of the torus.
High Energy Astrophysical Phenomena,General Relativity and Quantum Cosmology
What problem does this paper attempt to address?
This paper aims to solve the problem of long - term Monte Carlo (MC) - based neutrino radiation hydrodynamics simulations in black - hole - accretion - disk systems. Specifically, the author has developed a new general - relativistic MC neutrino radiation hydrodynamics code to address several improvements in axisymmetric systems, including: 1. **Extended Implicit MC Method**: An extended version of the implicit MC method has been developed for multi - species radiation fields. 2. **Modeling of Neutrino Pair - Process Rates**: The neutrino pair - process rates are modeled based on a new numerically efficient and asymptotically correct fitting function. 3. **New Numerical Limiters for Radiation - Matter Interactions**: New numerical limiters are implemented to ensure the stable and physically correct evolution of the system. The author applies this code to study a black - hole - accretion - disk system, which consists of a black hole with a mass of 3 solar masses and a dimensionless spin of 0.8, and an accretion disk with a mass of 0.1 solar masses. This system simulates the black - hole - accretion - disk system formed after the merger of two neutron stars. Through the simulation, the author tracks the evolution of the system for more than 1 second and takes into account the dynamic influence of non - thermal pair annihilation. The study finds that the system evolution and various key parameters (such as neutrino luminosity, ejecta mass, accretion - disk electron fraction, and pair - annihilation luminosity) are roughly consistent with the results of previous studies. In addition, the author also finds that neutrino - antineutrino pair annihilation can initiate a relativistic outflow within a time scale of about 0.1 seconds, and its energy is sufficient to explain some short - hard gamma - ray bursts and their precursors. Finally, the author directly calculates the indicators of fast - flavor instability from the obtained neutrino distribution function, indicating that this instability is particularly likely to occur in the equatorial region of the accretion disk. Overall, the main objective of this paper is to more accurately simulate the neutrino radiation process in black - hole - accretion - disk systems through the improved MC method, so as to better understand the physical phenomena in these systems, especially the influence of neutrino pair annihilation on the relativistic outflow and the potential location of fast - flavor instability.