Search for: All records

Abstract We introduceWeakhub, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino–matter, neutrino–neutrino interactions and plasma process, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations ofβ-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino–antineutrino annihilation, and nuclear de-excitation are also included. We also present the detailed derivations of weak interactions and the coupling to the two-moment based general-relativistic multigroup radiation transport in the general-relativisticmultigridnumerical (Gmunu) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernovae and binary neutron star (BNS) postmerger remnants, and predict the effects of improved opacities in comparison to conventional ones for a BNS postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino libraryNuLibin a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than ∼10% in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modeling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations. 

Abstract We present the implementation of a two-moment-based general-relativistic multigroup radiation transport module in theGeneral-relativisticmultigridnumerical (Gmunu) code. On top of solving the general-relativistic magnetohydrodynamics and the Einstein equations with conformally flat approximations, the code solves the evolution equations of the zeroth- and first-order moments of the radiations in the Eulerian-frame. An analytic closure relation is used to obtain the higher order moments and close the system. The finite-volume discretization has been adopted for the radiation moments. The advection in spatial space and frequency-space are handled explicitly. In addition, the radiation–matter interaction terms, which are very stiff in the optically thick region, are solved implicitly. The implicit–explicit Runge–Kutta schemes are adopted for time integration. We test the implementation with a number of numerical benchmarks from frequency-integrated to frequency-dependent cases. Furthermore, we also illustrate the astrophysical applications in hot neutron star and core-collapse supernovae modelings, and compare with other neutrino transport codes.