Full Boltzmann Neutrino Transport Research Articles

Detectability of Late-time Supernova Neutrinos with Fallback Accretion onto Protoneutron Star

We investigate the late-time neutrino emission powered by fallback mass accretion onto a protoneutron star (PNS), using neutrino radiation-hydrodynamic simulations with full Boltzmann neutrino transport. We follow the time evolution of the accretion flow onto the PNS until the system reaches a quasi-steady state. A standing shock wave is commonly formed in the accretion flow, whereas the shock radius varies depending on the mass accretion rate and the PNS mass. A sharp increase in temperature emerges in the vicinity of the PNS (∼10 km), which characterizes neutrino emission. Both the neutrino luminosity and the average energy become higher with increasing mass accretion rate and PNS mass. The mean energy of the emitted neutrinos is in the range of 10 ≲ ϵ ≲ 20 MeV, which is higher than that estimated from PNS cooling models (≲10 MeV). Assuming a distance to core-collapse supernova of 10 kpc, we quantify neutrino event rates for Super-Kamiokande (Super-K) and Deep Underground Neutrino Experiment (DUNE). The estimated detection rates are well above the background, and their energy-dependent features are qualitatively different from those expected from PNS cooling models. Another notable feature is that the neutrino emission is strongly flavor dependent, exhibiting that the neutrino event rate hinges on the neutrino oscillation model. We estimate them in the case with the adiabatic Mikheev–Smirnov–Wolfenstein model, and show that the normal and inverted mass hierarchy offer a large number of neutrino detections in Super-K and DUNE, respectively. Hence the simultaneous observation with Super-K and DUNE of fallback neutrinos will provide a strong constraint on the neutrino mass hierarchy.

General-relativistic quantum-kinetics neutrino transport

We developed a new general-relativistic quantum-kinetics neutrino transport code, GRQKNT, for numerical studies of quantum kinetics of non-equilibrium neutrinos in six-dimensional phase space. This code is intended for use in both local and global simulations of neutrino transport in core-collapse supernova and binary neutron star merger. It has been widely recognized that global simulations of collective neutrino oscillations, in particular fast neutrino-flavor conversions (FFC), require unfeasible computational resources due to large disparity of scales between flavor conversion and astrophysical phenomena. We propose a novel approach to tackle the issue. This paper is devoted to describe the philosophy, design, and numerical implementation of GRQKNT with a number of tests ensuring correct implementation of each module. The code is based on a discrete-ordinate Sn method, finite-difference realization of mean-field quantum kinetic equation. The transport equation is solved based on a conservative formalism, and we use a fifth-order weighted essentially non-oscillatory scheme with fourth-order TVD Runge-Kutta time-integration. The transport module is designed to work with arbitrary spacetimes and currently three different stationary spacetimes (flat spacetime, Schwarzschild black hole, and Kerr black hole) are implemented. The collision term including neutrino emission, absorption, and momentum-exchanged scatterings are also implemented into our code. The oscillation Hamiltonian consists of vacuum, matter, and self-interactions. Both two- and three neutrino-flavor scenarios can be applied. Fluid-velocity dependences in transport-, collision-, and oscillation modules, are also treated self-consistently by using two-energy-grid technique, that has been already established in our another code with full Boltzmann neutrino transport.

Principal-axis Analysis of the Eddington Tensor for the Early Post-bounce Phase of Rotational Core-collapse Supernovae

Using full Boltzmann neutrino transport, we performed 2D core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2 and 15.0 M ⊙, both rotational and nonrotational. We employed the results obtained in the early post-bounce phase (t ≲ 20 ms) to assess performance under rapid rotation of some closure relations commonly employed in the truncated moment method. We first made a comparison in 1D under spherical symmetry, though, of the Eddington factor p defined in the fluid rest frame (FR). We confirmed that the maximum entropy closure for the Fermionic distribution (MEFD) performs better than others near the proto–neutron star surface, where p < 1/3 occurs, but does not work well even in 1D when the phase-space occupancy satisfies e < 0.5 together with p < 1/3, the condition known to be not represented by MEFD. For the 2D models with the rapid rotation, we employed the principal-axis analysis of the Eddington tensor. We paid particular attention to the direction of the longest principal axis. We observed in FR that it is aligned neither with the radial direction nor with the neutrino flux in 2D, particularly so in convective and/or rapidly rotating regions, the fact not accommodated in the moment method. We repeated the same analysis in the laboratory frame and found again that the direction of the longest principal axis is not well reproduced by MEFD because the interpolation between the optically thick and thin limits is not very accurate in this frame.

Prospects of Fast Flavor Neutrino Conversion in Rotating Core-collapse Supernovae

Abstract There is mounting evidence that neutrinos undergo fast flavor conversion (FFC) in core-collapse supernova (CCSN). In this paper, we investigate the roles of stellar rotation on the occurrence of FFC by carrying out axisymmetric CCSN simulations with full Boltzmann neutrino transport. Our result suggests that electron neutrino lepton number (ELN) angular crossings, which are the necessary and sufficient condition to trigger FFC, preferably occur in the equatorial region for rotating CCSNe. By scrutinizing the neutrino–matter interaction and neutrino radiation field, we find some pieces of evidence that the stellar rotation facilitates the occurrence of FFC. The low-electron-fraction region in the post-shock layer expands by centrifugal force, enhancing the disparity of neutrino absorption between electron-type neutrinos (ν e) and their anti-particles ( ν ¯ e ). This has a significant impact on the angular distribution of neutrinos in momentum space, in which ν e tends to be more isotropic than ν ¯ e ; consequently, ELN crossings emerge. The ELN crossing found in this study is clearly associated with rotation, which motivates further investigation on how the subsequent FFC influences explosion dynamics, nucleosynthesis, and neutrino signals in rotating CCSNe.

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernova (CCSN): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron neutrinos lepton number (ELN) angular crossings in momentum space is a good indicator of occurrences of the flavor conversions. However, many multi-dimensional (multi-D) CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits the access to the angular distributions of neutrinos, i.e., inhibits ELN-crossing searches. In this paper, we develop a new method of ELN-crossing search in these CCSN simulations. The required data is the zero-th and first angular moments of neutrinos and matter profile, all of which are available in CCSN models with two-moment method. One of the novelties of our new method is to use a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches only with the two angular moments. We assess the capability of the method by carrying out a detailed comparison to results of full Boltzmann neutrino transport in 1D and 2D CCSN models. We find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has a benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models employed two-moment neutrino transport.

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernova (CCSN), since they substantially reduce the computational burden inherent in the multi-dimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zero-th and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. We find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\kappa$). This issue originates from the insensitiveness of incoming neutrinos to $\kappa$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of ELN-crossing searches with two-moment neutrino transport. This analysis is complementary to our another paper that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\kappa \gtrsim 0.5$. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

Simulations of the Early Postbounce Phase of Core-collapse Supernovae in Three-dimensional Space with Full Boltzmann Neutrino Transport

Abstract We report on the core-collapse supernova simulation we conducted for a 11.2M ⊙ progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation-hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ∼10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation, whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ∼0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and to possibly improve the closure relation there.

Fast-pairwise Collective Neutrino Oscillations Associated with Asymmetric Neutrino Emissions in Core-collapse Supernovae

Abstract We present a linear stability analysis of the fast-pairwise neutrino flavor conversion based on a result of our latest axisymmetric core-collapse supernova (CCSN) simulation with full Boltzmann neutrino transport. In the CCSN simulation, coherent asymmetric neutrino emissions of electron-type neutrinos (ν e) and their antiparticles ( ), in which the asymmetries of ν e and are anticorrelated with each other, occur at almost the same time as the onset of aspherical shock expansion. We find that the asymmetric neutrino emissions play a crucial role on occurrences of fast flavor conversions. The linear analysis shows that unstable modes appear in both pre- and post-shock flows; for the latter, they appear only in the hemisphere of higher emissions (the same hemisphere with stronger shock expansion). We analyze the characteristics of electron–lepton number (ELN) crossing in depth by closely inspecting the angular distributions of neutrinos in momentum space. The ELN crossing happens in various ways, and the property depends on the radius: in the vicinity of neutron star, (ν e) dominates over ν e ( ) in the forward (backward) direction; at the larger radius, the ELN crossing occurs in the opposite way. We also find that the non-radial ELN crossing occurs at the boundary between no ELN crossing and the radial one, which is an effect of genuine multi-dimensional transport. Our findings indicate that the collective neutrino oscillation may occur more commonly in CCSNe and suggest that the CCSN community needs to accommodate these oscillations self-consistently in the modeling of CCSNe.

Three-dimensional Boltzmann-hydro Code for Core-collapse in Massive Stars. III. A New Method for Momentum Feedback from Neutrino to Matter

Abstract We present a new method for neutrino-matter coupling in multi-dimensional radiation-hydrodynamic simulations of core-collapse supernovae (CCSNe) with the full Boltzmann neutrino transport. This development is motivated by the fact that accurate conservation of momentum is required for reliable numerical modelings of CCSN dynamics including a recoil of proto-neutron stars (PNSs). The new method is built on a hybrid approach in which we use the energy-momentum tensor of neutrinos to compute the momentum feedback from neutrino to matter in the optically thick region while we employ the collision term in the optically thin region. In this method we utilize a general relativistic description of radiation-hydrodynamics with angular moments, which allows us to evaluate the momentum feedback from neutrino to matter without inconsistency with our Boltzmann solver. We demonstrate that the new method substantially improves the accuracy of linear momentum conservation in our CCSN simulations under reasonable angular resolutions in momentum space, alleviating the difficulty in giving the diffusion limit precisely with the discrete ordinate (S n ) method. It is the first ever demonstration that the PNS kick can be handled directly and properly in multi-dimensional radiation-hydrodynamic simulations with the full Boltzmann neutrino transport.

Radiation-Hydrodynamic Simulations of Core-collapse Supernovae with 6 Dimensional Boltzmann Neutrino Transport

Core-collapse supernovae (CCSNe) are intrinsically multi-scale, multi-physics and multi-dimensional phenomena. Because of the enormous complexity, the first-principles numerical simulations under realistic input physics are strongly required to uncover the explosion mechanism, predict observational signals (neutrinos, gravitational waves and electromagnetic waves) and prove physical state in extremely hot and dense matter of supernova core. We have tackled the development of multi-dimensional radiation-hydrodynamic code with full Boltzmann neutrino transport and performed several scientific CCSNe simulations in the last few years. In this article, we report the recent progress of our CCSNe numerical simulations with the most up-to-date equation-of-state (EOS) and nuclear weak interactions. We also present preliminary results of non-rotating CCSNe simulations in spatial axisymmetry.

Comparing Treatments of Weak Reactions with Nuclei in Simulations of Core-collapse Supernovae

Abstract We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae, paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ∼10% and in average neutrino energies by 0.25–1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

Simulations of Core-collapse Supernovae in Spatial Axisymmetry with Full Boltzmann Neutrino Transport

Abstract We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent five-dimensional (two in space and three in momentum space) problem. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2 M ☉, employing different nuclear equations of state (EOSs): Lattimer and Swesty’s EOS with the incompressibility of K = 220 MeV (LS EOS) and Furusawa’s EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS, the shock wave reaches ∼700 km at 300 ms after bounce and is still expanding, whereas in the FS EOS it stalled at ∼200 km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire postbounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find nonaxisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the rθ component reaches ∼10% of the dominant rr component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the θθ component, in the semitransparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

A NEW MULTI-ENERGY NEUTRINO RADIATION-HYDRODYNAMICS CODE IN FULL GENERAL RELATIVITY AND ITS APPLICATION TO THE GRAVITATIONAL COLLAPSE OF MASSIVE STARS

ABSTRACT We present a new multi-dimensional radiation-hydrodynamics code for massive stellar core-collapse in full general relativity (GR). Employing an M1 analytical closure scheme, we solve spectral neutrino transport of the radiation energy and momentum based on a truncated moment formalism. Regarding neutrino opacities, we take into account a baseline set in state-of-the-art simulations, in which inelastic neutrino–electron scattering, thermal neutrino production via pair annihilation, and nucleon–nucleon bremsstrahlung are included. While the Einstein field equations and the spatial advection terms in the radiation-hydrodynamics equations are evolved explicitly, the source terms due to neutrino–matter interactions and energy shift in the radiation moment equations are integrated implicitly by an iteration method. To verify our code, we first perform a series of standard radiation tests with analytical solutions that include the check of gravitational redshift and Doppler shift. A good agreement in these tests supports the reliability of the GR multi-energy neutrino transport scheme. We then conduct several test simulations of core-collapse, bounce, and shock stall of a 15 M &CircleDot; ?> star in the Cartesian coordinates and make a detailed comparison with published results. Our code performs quite well to reproduce the results of full Boltzmann neutrino transport especially before bounce. In the postbounce phase, our code basically performs well, however, there are several differences that are most likely to come from the insufficient spatial resolution in our current 3D-GR models. For clarifying the resolution dependence and extending the code comparison in the late postbounce phase, we discuss that next-generation Exaflops class supercomputers are needed at least.

AN OPEN-SOURCE NEUTRINO RADIATION HYDRODYNAMICS CODE FOR CORE-COLLAPSE SUPERNOVAE

We present an open-source update to the spherically symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D. The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general-relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al. and Cardall et al. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations is the neutrino–matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress toward robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino-transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general-relativistic transport code to model failed CCSNe by evolving a 40-solar-mass progenitor to the onset of collapse to a black hole.

High-resolution supernova neutrino spectra represented by a simple fit

To study the capabilities of supernova neutrino detectors, the instantaneous spectra are often represented by a quasi-thermal distribution of the form f(E) = E^alpha e^{-(alpha+1)E/E_{av}} where E_{av} is the average energy and alpha a numerical parameter. Based on a spherically symmetric supernova model with full Boltzmann neutrino transport we have, at a few representative post-bounce times, re-converged the models with vastly increased energy resolution to test the fit quality. For our examples, the spectra are well represented by such a fit in the sense that the counting rates for a broad range of target nuclei, sensitive to different parts of the spectrum, are reproduced very well. Therefore, the mean energy and root-mean-square energy of numerical spectra hold enough information to provide the correct alpha and to forecast the response of multi-channel supernova neutrino detection.