Search

Recent neutron star surface observations corroborate a long-standing theory that neutron stars may be accelerated over extended periods after their birth. We analyze how these prolonged rocket-like accelerations, combined with rapid birth kicks, impact binary orbits. We find that even a small contribution of rocket kicks combined with instantaneous natal kicks can allow binaries to reach period--eccentricity combinations unattainable in standard binary evolution models. We propose these kick + rocket combinations as a new channel to form wide low-eccentricity neutron star binaries such as Gaia NS1, as well as inducing stellar mergers months to years after a supernova to cause peculiar high-energy transients., Comment: 11 pages, 7 figures, published in ApJL

In the standard model of core-collapse supernova (CCSN), all neutrinos are assumed to be in pure flavor eigenstates in CCSN cores, but the assumption becomes invalid if neutrino distributions are unstable to flavor conversions. In this paper, we present a study of the occurrences of two representative neutrino-flavor instabilities, fast- and collisional flavor instabilities, in the cooling phase of proto-neutron star (PNS) from 1- to 50 seconds. We follow the long-term evolution of a PNS under spherically symmetric and quasi-static approximations, in which the matter profile is determined by solving the Tolman-Oppenheimer-Volkoff equation with neutrino feedback under the treatment of multi-group flux limited diffusion. For the stability analysis of neutrino flavor conversions, we recompute neutrino distributions using Monte Carlo transport in order to obtain the full angular distribution needed to compute the dispersion relations. We find no signs of flavor conversions in our models; the physical reason is thoroughly investigated. We also argue that the negative conclusion in flavor conversions could be changed qualitatively if multi-dimensional effects are included, as similar to cases in the earlier phase of CCSN., Comment: 14 pages, 9 figures. Submitted to ApJ

Neutrinos are known to undergo flavor conversion among their three flavors. In the theoretical modeling of core-collapse supernova (CCSN), there has been a great deal of attention to recent discoveries of a new type of neutrino flavor conversions, namely collisional flavor instability (CFI), in which the instability is induced by the flavor-dependent decoherence due to the disparity of neutrino-matter interactions among flavors. In this paper, we study how the appearance of on-shell muons and associated neutrino-matter interactions can impact CFIs based on linear stability analysis of flavor conversions. Some striking results emerge from the present study. First, we analytically show that breaking beta- and pair equilibrium is a necessary condition to trigger CFIs. This also indicates that CFIs with on-shell muons could appear in $e \tau$ and $\mu \tau$ neutrino mixing sectors in very high-density region ($\gtrsim 10^{13} {\rm g/cm^{3}}$), exhibiting a possibility of large impacts of CFIs on CCSN. Second, resonance-like CFIs, having a much higher growth rate than normal CFIs, can be triggered by muons. The resonance point of CFIs is different between $e \tau$ and $\mu \tau$ sectors; the former (latter) occurs at $\mu_{e (\mu)} = \mu_{n} - \mu_{p}$, where $\mu_{i}$ denotes the chemical potential of $i$ constitute ($n$ and $p$ represent neutrons and protons, respectively). Our result suggests that the non-linear evolution of CFI with on-shell muons would induce flavor conversions with the complex interplay among all three different neutrino-mixing sectors.

Multi-dimensional nature of core-collapse supernova (CCSN) leads to asymmetric matter ejection and neutrino emission, that potentially accounts for the origin of neutron star (NS) kick. Asymmetric neutrino radiation fields are, in general, accompanied by large-scale inhomogeneous fluid distributions, in particular for electron-fraction ($Y_e$) distributions. Recently, it has also been revealed that lower $Y_e$ environments in proto-neutron star envelope can offer preferable conditions for collective neutrino oscillations. In this paper, we show that a dipole asymmetry of fast neutrino-flavor conversion (FFC), one of the collective neutrino oscillation modes, can power a NS kick, and that it would generate a characteristic correlation between asymmetric distributions of heavy elements in the ejecta and the direction of NS kick. We strengthen our argument for the FFC-driven NS kick mechanism by performing axisymmetric neutrino transport simulations with full Boltzmann neutrino transport. We show that this mechanism can generate sufficient linear momentum of neutrinos to account for typical proper motions of NS. Although more detailed studies are necessary, the present study opens a new channel to give a natal NS kick., Comment: Accepted to PRD

In dense neutrino environments like core-collapse supernovae (CCSNe) and neutron star mergers (NSMs), neutrinos can undergo fast flavor conversions (FFC) when their angular distribution of neutrino electron lepton number ($\nu$ELN) crosses zero along some directions. While previous studies have demonstrated the detection of axisymmetric $\nu$ELN crossings in these extreme environments, non-axisymmetric crossings have remained elusive, mostly due to the absence of models for their angular distributions. In this study, we present a pioneering analysis of the detection of non-axisymmetric $\nu$ELN crossings using machine learning (ML) techniques. Our ML models are trained on data from two CCSN simulations, one with rotation and one without, where non-axisymmetric features in neutrino angular distributions play a crucial role. We demonstrate that our ML models achieve detection accuracies exceeding 90\%. This is an important improvement, especially considering that a significant portion of $\nu$ELN crossings in these models eluded detection by earlier methods., Comment: 10 pages, 2 figures

We present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in quantum kinetic equation by Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that non-linear phases of collective neutrino oscillations settle into quasi-steady structures. We explicitly provide basic equations of the BGK subgrid model for both multi-angle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within $\sim 20 \%$. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model., Comment: Accepted to PRD

We assert that non-linear features of fast neutrino-flavor conversion (FFC) can be qualitatively different between core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs). This argument arises from recent global FFC simulations in BNSM, in which fast flavor swap (FFS) emerges in very narrow spatial regions, whereas neutrinos in CCSN tend to evolve towards flavor equipartition. In this {\it Letter}, we provide the physical mechanism of FFS based on a colliding neutrino beam model. Neutrinos/antineutrinos can undergo FFS when they propagate in ambient neutrino gas that propagates in the opposite direction and also has the opposite sign of ELN-XLN, where ELN and XLN denote electron- and heavy-leptonic neutrino number, respectively. Such environments can be naturally realized in BNSMs, whereas they are unlikely in CCSNe unless the neutrino sphere is strongly deformed aspherically. Our study exhibits the diversity of non-linear dynamics of FFC., Comment: Accepted to PRD. 7 pages, 4 figures

We present a comprehensive study on the occurrences of the collisional flavor instability (CFI) and the fast flavor instability (FFI) of neutrinos based on a two-dimensional (2D) core-collapse supernova (CCSN) simulation performed with a Boltzmann radiation hydrodynamics code. We find that CFI occurs in a region with the baryon-mass density of $10^{10}\lesssim \rho \lesssim 10^{12}\,\mathrm{g}\,\mathrm{cm}^{-3}$, which is similar to the previous results in one-dimensional (1D) CCSN models. In contrast to 1D, however, the CFI region varies with time vigorously in the 2D model, whereas it had a quiescent structure in 1D. This is attributed to the fact that the turbulent flows advected from a gain region account for the temporal variations. Another noticeable difference from the 1D models is the appearance of resonance-like CFI where number densities of $\nu_e$, $\bar\nu_e$ nearly coincide each other. The CFI growth rate there is enhanced and can reach $\sim10^{-3}\,\mathrm{cm}^{-1}$. As for FFI, on the other hand, it appears in three different regions; (1) the region overlapped with the resonance-like CFI, (2) neutrino decoupling regions where $\bar{\nu}_e$'s are strongly emitted, and (3) optically thin regions where neutral-current scatterings dominate over charged-current reactions. Although overall properties for FFI are consistent with previous studies, we find that the number of electron-neutrinos lepton number crossing (ELN crossing) temporary becomes multiple, which can be assessed accurately only by multi-angle treatments in neutrino transport. We find that the growth rate of FFI is always higher than CFI if both of them occur, which suggests that the former is dominant for the linear evolution., Comment: published in PRD

The sensitivity of current and future neutrino detectors like Super-Kamiokande (SK), JUNO, Hyper-Kamiokande (HK), and DUNE is expected to allow for the detection of the diffuse supernova neutrino background (DSNB). However, the DSNB model ingredients like the core-collapse supernova (CCSN) rate, neutrino emission spectra, and the fraction of failed supernovae are not precisely known. We quantify the uncertainty on each of these ingredients by (i) compiling a large database of recent star formation rate density measurements, (ii) combining neutrino emission from long-term axisymmetric CCSNe simulations and strategies for estimating the emission from the protoneutron star cooling phase, and (iii) assuming different models of failed supernovae. Finally, we calculate the fluxes and event rates at multiple experiments and perform a simplified statistical estimate of the time required to significantly detect the DSNB at SK with the gadolinium upgrade and JUNO. Our fiducial model predicts a flux of $5.1\pm0.4^{+0.0+0.5}_{-2.0-2.7}\,{\rm cm^2~s^{-1}}$ at SK employing Gd-tagging, or $3.6\pm0.3^{+0.0+0.8}_{-1.6-1.9}$ events per year, where the errors represent our uncertainty from star formation rate density measurements, uncertainty in neutrino emission, and uncertainty in the failed-supernova scenario. In this fiducial calculation, we could see a $3\sigma$ detection by $\sim2030$ with SK-Gd and a $5\sigma$ detection by $\sim2035$ with a joint SK-Gd/JUNO analysis, but background reduction remains crucial., Comment: 19 pages, 9 figures, 3+2 tables. Comments welcome

Neutrinos are known to undergo flavor conversion processes among the three flavors. The fast flavor conversion (FFC) has been the central piece of flavor conversions taking place in core-collapse supernovae (CCSNe) due to its shorter timescale to the completion of flavor conversion compared to other types of flavor conversion. Although the ordinary collisions between neutrinos and matter were once thought to decohere neutrinos and thus damp flavor conversions, it was recently realized that they can also induce the flavor conversion. The linear analysis showed that the so-called collisional flavor instability or CFI occurs in the absence of FFC. In this paper, we investigate if CFI takes place in of the post-bounce core of CCSNe, using the results of spherically symmetric Boltzmann simulations of CCSNe for four progenitor models with different masses. We also provide a necessary (but not sufficient) condition of matter properties for the occurrence of CFI in optically thick and semi-transparent regions; baryon mass density ($\rho$), electron fraction ($Y_e$), and the degeneracy of electron-type neutrinos ($\eta_{\nu_e}$) need to be $10^{10} {\rm g/cm^3} \lesssim \rho \lesssim 10^{12} {\rm g/cm^3}$, $Y_e\lesssim 0.4$, and $\eta_{\nu_e} \lesssim 0.5$, respectively. This condition allows us to easily locate the place of possible CFI occurence without detailed stability analyses, which is useful for analyzing CFI in CCSN models phenomenologically, Comment: accepted in PRD

Neutrinos in dense environments like core-collapse supernovae (CCSNe) and neutron star mergers (NSMs) can undergo fast flavor conversions (FFCs) once the angular distribution of neutrino lepton number crosses zero along a certain direction. Recent advancements have demonstrated the effectiveness of machine learning (ML) in detecting these crossings. In this study, we enhance prior research in two significant ways. Firstly, we utilize realistic data from CCSN simulations, where neutrino transport is solved using the full Boltzmann equation. We evaluate the ML methods' adaptability in a real-world context, enhancing their robustness. In particular, we demonstrate that when working with artificial data, simpler models outperform their more complex counterparts, a noteworthy illustration of the bias-variance tradeoff in the context of ML. We also explore methods to improve artificial datasets for ML training. In addition, we extend our ML techniques to detect the crossings in the heavy-leptonic channels, accommodating scenarios where $\nu_x$ and $\bar\nu_x$ may differ. Our research highlights the extensive versatility and effectiveness of ML techniques, presenting an unparalleled opportunity to evaluate the occurrence of FFCs in CCSN and NSM simulations., Comment: 12 pages, 3 figures. v2: discussions on the performance of ML in the neutrino decoupling region added, matches published version

Neutrinos play pivotal roles in determining fluid dynamics, nucleosynthesis, and their observables in core-collapse supernova (CCSN) and binary neutron star merger (BNSM). In this paper, we present a novel phenomenon, collisional flavor swap, in which neutrino-matter interactions trigger the complete interchange of neutrino spectra between two different flavors, aided by neutrino self-interactions. We find that a necessary condition to trigger the collisional swap is occurrences of resonance-like collisional flavor instability. In cases where neutrino self-interactions substantially dominate over the collision rate, the collisional swap occurs in the entire neutrino energy spectrum, while intriguing energy dependent features can emerge after the completion of flavor swap. Since flavor swaps correspond to the most extreme case in flavor conversions, they have a great potential to affect CCSN and BNSM phenomena., Comment: 8 pages, 4 figures, accepted to PRD

One of the active debates in core-collapse supernova (CCSN) theory is how significantly neutrino flavor conversions induced by neutrino-neutrino self-interactions change the conventional picture of CCSN dynamics. Recent studies have indicated that strong flavor conversions can occur inside neutrino spheres where neutrinos are tightly coupled to matter. These flavor conversions are associated with either collisional instability or fast neutrino-flavor conversion (FFC) or both. The impact of these flavor conversions on CCSN dynamics is, however, still highly uncertain due to the lack of global simulations of quantum kinetic neutrino transport with appropriate microphysical inputs. Given fluid profiles from a recent CCSN model at three different time snapshots in the early post-bounce phase, we perform global quantum kinetic simulations in spherical symmetry with an essential set of microphysics. We find that strong flavor conversions occur in optically thick regions, resulting in a substantial change of neutrino radiation field. The neutrino heating in the gain region is smaller than the case with no flavor conversions, whereas the neutrino cooling in the optically thick region is commonly enhanced. Based on the neutrino data obtained from our multi-angle neutrino transport simulations, we also assess some representative classical closure relations by applying them to diagonal components of density matrix of neutrinos. We find that Eddington tensors can be well approximated by these closure relations except for the region where flavor conversions occur vividly. We also analyze the neutrino signal by carrying out detector simulations for Super-Kamiokande, DUNE, and JUNO. We propose a useful strategy to identify the sign of flavor conversions in neutrino signal, that can be easily implemented in real data analyses of CCSN neutrinos., Comment: Accepted to PRD

Recent multi-dimensional (multi-D) core-collapse supernova (CCSN) simulations characterize gravitational waves (GWs) and neutrino signals, offering insight into universal properties of CCSN independent of progenitor. Neutrino analysis in real observations, however, will be complicated due to the ambiguity of self-induced neutrino flavor conversion (NFC), which poses an obstacle to extracting detailed physical information. In this paper, we propose a novel approach to place a constraint on NFC from observed quantities of GWs and neutrinos based on correlation analysis from recent, detailed multi-D CCSN simulations. The proposed method can be used even in cases with low significance - or no detection of GWs. We also discuss how we can utilize electro-magnetic observations to complement the proposed method. Although our proposed method has uncertainties associated with CCSN modeling, the present result will serve as a base for more detailed studies. Reducing the systematic errors involved in CCSN models is a key to success in this multi-messenger analysis that needs to be done in collaboration with different theoretical groups., Comment: Adding references, 9 pages, 4 figures, accepted to PRD

Binary neutron star merger (BNSM) offers an environment where fast neutrino-flavor conversion (FFC) can vividly occur, that potentially leads to a considerable change of neutrino radiation field. In this Letter, we investigate global features of FFC by general relativistic quantum kinetic neutrino transport simulations in spatial axisymmetry. Our result suggests that global advection of neutrinos plays a crucial role in FFC dynamics. Although flavor conversions occur ubiquitously in the early phase, they can be active only in a narrow region in the late phase. This region includes an ELN-XLN Zero Surface (EXZS), corresponding to a surface where electron-neutrinos lepton number (ELN) equals to heavy-leptonic one (XLN). The EXZS is not stationary, but dynamically evolve in a time scale of global advection. We also find that neutrinos can undergo a flavor swap when they pass through the EXZS, resulting in qualitatively different neutrino radiation fields between both sides of EXZS. Our result suggests that EXZS is one of the key ingredients to characterize FFC in BNSM., Comment: Accepted to PRD

We investigate the late-time neutrino emission powered by fallback mass accretion onto proto-neutron star (PNS), using neutrino radiation-hydrodynamic simulations with full Boltzmann neutrino transport. We follow the time evolution of accretion flow onto PNS until the system reaches a steady state. A standing shock wave is commonly formed in the accretion flow, whereas the shock radius varies depending on mass accretion rate and PNS mass. A sharp increase in temperature emerges in the vicinity of PNS ($\sim 10$ km), which characterizes neutrino emission. Both neutrino luminosity and average energy become higher with increasing mass accretion rate and PNS mass. The mean energy of emitted neutrinos is in the range of $10\lesssim\epsilon\lesssim20\,\mathrm{MeV}$, which is higher than that estimated from PNS cooling models ($\lesssim10\,\mathrm{MeV}$). Assuming a distance to core-collapse supernova of $10\,\mathrm{kpc}$, we quantify neutrino event rates for Super-Kamiokande (Super-K) and 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 adiabatic Mikheev-Smirnov-Wolfenstein model, and show that the normal- and inverted mass hierarchy offer the large number of neutrino detection in Super-K and DUNE, respectively. Hence the simultaneous observation with Super-K and DUNE of the fallback neutrinos will provide a strong constraint on neutrino mass hierarchy., Comment: published in ApJ

The question of what ingredients characterize the quasi-steady state of fast neutrino-flavor conversion (FFC) is one of the long-standing riddles in neutrino oscillation. Addressing this issue is necessary for accurate modeling of neutrino transport in core-collapse supernova and binary neutron star merger. Recent numerical simulations of FFC have shown, however, that the quasi-steady state is sensitively dependent on boundary conditions in space, and the physical reason for the dependence is not clear at present. In this study, we provide a physical interpretation of this issue based on arguments with stability and conservation laws. The stability can be determined by the disappearance of ELN(electron neutrino-lepton number)-XLN(heavy-leptonic one) angular crossings, and we also highlight two conserved quantities characterizing the quasi-steady state of FFC: (1) lepton number conservation along each neutrino trajectory and (2) conservation law associated with angular moments, depending on boundary conditions, for each flavor of neutrinos. We present an analytic prescription that matches the results of the nonlinear simulations presented in this work. This study represents a major step forward a unified picture determining asymptotic states of FFCs., Comment: 10 pages, 5 figures. Accepted in PRD

Fast neutrino flavor conversions (FFCs) and collisional flavor instabilities (CFIs) potentially affect the dynamics of core-collapse supernovae (CCSNe) and binary neutron star mergers (BNSMs). Under the assumption of homogeneous neutrinos, we investigate effects of neutrino emission and absorption (EA) by matters through both single and multi-energy numerical simulations with physically motivated setup. In our models, FFCs dominate over CFIs in the early phase, while EA secularly and significantly give impacts on flavor conversions. They facilitate angular swaps, or the full exchange between electron neutrinos ($\nu_e$) and heavy-leptonic neutrinos ($\nu_x$). As a result, the number density of $\nu_x$ becomes more abundant than the case without EA, despite the fact that the isotropization by EA terminates the FFCs earlier. In the later phase, the system approaches new asymptotic states characterized by EA and CFIs, in which rich energy-dependent structures also emerge. Multi-energy effects sustain FFCs and the time evolution of the flavor conversion becomes energy dependent, which are essentially in line with effects of the isoenergetic scattering studied in our previous paper. We also find that $\nu_x$ in the high-energy region convert into $\nu_e$ via flavor conversions and then they are absorbed through charged current reactions, exhibiting the possibility of new path of heating matters., Comment: 18 pages, 17 figures, accepted to PRD

One of the greatest uncertainties in any modeling of inner engine of core-collapse supernova (CCSN) is neutrino flavor conversions driven by neutrino self-interactions. We carry out large-scale numerical simulations of multi-energy, multi-angle, three-flavor framework, and general relativistic quantum kinetic neutrino transport in spherical symmetry with an essential set of neutrino-matter interactions under a realistic fluid profile of CCSN. Our result suggests that the neutrino heating in the gain region is reduced by $\sim 40\%$ due to fast neutrino-flavor conversion (FFC). We also find that the total luminosity of neutrinos is enhanced by $\sim 30 \%$, for which the substantial increase of heavy-leptonic neutrinos by FFCs are mainly responsible. This study provides evidence that FFC has a significant impact on the delayed neutrino-heating mechanism., Comment: Accepted to PRL

Neutrino-neutrino forward scatterings potentially induce collective neutrino oscillation in dense neutrino gases in astrophysical sites such as core-collapse supernovae (CCSN) and binary neutron star mergers (BNSM). In this paper, we present a detailed study of fast neutrino-flavor conversion (FFC), paying special attention to asymptotic states, by linear stability analysis and local simulations with a periodic boundary condition. We find that asymptotic states can be characterized by two key properties of FFC: (1) the conservation of lepton number for each flavor of neutrinos and (2) the disappearance of ELN(electron neutrino-lepton number)-XLN(heavy-leptonic one) angular crossings in the spatial- or time-averaged distributions. The system which initially has the positive (negative) ELN-XLN density reaches a flavor equipartition in the negative (positive) ELN-XLN angular directions, and the other part compensates it to preserve the conservation laws. These properties of FFCs offer an approximate scheme determining the survival probability of neutrinos in asymptotic states without solving quantum kinetic equations. We also demonstrate that the total amount of flavor conversion can vary with species-dependent neutrino distributions for identical ELN-XLN ones. Our results suggest that even shallow or narrow ELN angular crossings have the ability to drive large flavor conversion, exhibiting the need for including the effects of FFCs in the modeling of CCSN and BNSM., Comment: 11 pages, 6 figures, Accepted to PRD

We present a systematic study of fast neutrino-flavor conversion (FFC) with both small-scale and large-scale numerical simulations in spherical symmetry. We find that FFCs can, in general, reach a quasi-steady state, and these features in the non-linear phase are not characterized by the growth rate of FFC instability but rather angular structures of electron neutrino lepton number (ELN) and heavy one (XLN). Our result suggests that neutrinos can almost reach a flavor equipartition even in cases with low growth rate of instability (e.g., shallow ELN crossing) and narrow angular regions (in momentum space) where flavor conversions occur vigorously. This exhibits that ELN and XLN angular distributions can not provide a sufficient information to determine total amount of flavor conversion in neutrinos and antineutrinos of all flavors. Based on the results of our numerical simulations, we provide a new approximate scheme of FFC that is designed so that one can easily incorporate effects of FFCs in existing classical neutrino transport codes for the study of core-collapse supernova (CCSN) and binary neutron star merger (BNSM). The scheme has an ability to capture key features of quasi-steady state of FFCs without solving quantum kinetic neutrino transport, which will serve to facilitate access to FFCs for CCSN and BNSM theorists., Comment: Accepted to PRD

Herein, we review the nuclear equations of state (EOSs) %for core-collapse supernova simulations and the constituent nuclei of core-collapse supernovae (CCSNe) and their roles in CCSN simulations. Various nuclei such as deuterons, iron, and extremely neutron-rich nuclei compose in the central engines of CCSNe. The center of a collapsing core is dominated by neutron-rich heavy nuclei prior to the occurrence of core bounce. Their weak interactions significantly affect the neutrino emission and the size of the produced proto-neutron star. After a core bounce, heavy nuclei are dissolved to protons, neutrons, and light nuclei between the expanding shock wave and the newly formed neutron star (NS). Some of the key components in determining the shock-wave dynamics and supernova explosion of outer envelopes are neutrino interactions of nucleons and light nuclei such as deuterons. An EOS provides the relations between thermodynamical properties and the nuclear composition, and is needed to simulate this explosion. Further investigations on uniform and non-uniform nuclear matter are needed to improve the understanding of the mechanism of CCSNe and the properties of supernova nuclei. The knowledge of the EOS for uniform nuclear matter is being continually improved by a combination of microscopic calculations, terrestrial experiments, and NS observations. With reference to various nuclear experiments and current theories, the finite temperature effects on heavy nuclei, formation of light nuclei in dilute nuclear matter, and transition to uniform nuclear matter should be improved in the model of the EOS for non-uniform nuclear matter., Comment: 57 pages, 23 figures, Accepted by Prog. Part. Nucl. Phys

Fast neutrino-flavor conversion (FFC) is a possible game-changing ingredient in core-collapse supernova (CCSN) theory. In this paper, we examine the impact of FFC on explosive nucleosynthesis by including the effects of FFC in conjunction with asymmetric neutrino emission into nucleosynthetic computations in a parametric way. We find that the ejecta compositions are not appreciably affected by FFC for elements lighter than Co, while the compositions are influenced by FFC for the heavier elements. We also find that the role of FFC varies depending on the asymmetric degree of neutrino emission ($m_{\rm asy}$) and the degree of neutrino-flavor mixing. The impact of FFC is not monotonic to $m_{\rm asy}$; The change in the ejecta composition increases for higher $m_{\rm asy}$ up to $\sim 10\%$ compared with that without FFC, whereas FFC has little effect on the nucleosynthesis in very large asymmetric neutrino emission ($\gtrsim 30 \%$). Our results suggest that FFC facilitates the production of neutron-rich ejecta in most cases, although it makes the ejecta more proton-rich if anti-neutrino conversion is more vigorous than that of neutrino. The key ingredient accounting for this trend is neutrino absorption, whose effects on nucleosynthesis can be quantified by simple diagnostics., Comment: 7 pages, 6 figures, accepted for the publication in MNRAS

Upcoming LIGO/Virgo/KAGRA (LVK) observing runs are expected to detect a variety of inspiralling gravitational-wave (GW) events, that come from black-hole and neutron-star binary mergers. Detection of non-inspiral GW sources is also anticipated. We report the discovery of a new class of non-inspiral GW sources - the end states of massive stars - that can produce the brightest simulated stochastic GW burst signal in LVK bands known to date, and could be detectable in the LVK run A+. Some dying massive stars launch bipolar relativistic jets, which inflate a turbulent energetic bubble - cocoon - inside of the star. We simulate such a system using state-of-the-art 3D general-relativistic magnetohydrodynamic simulations and show that these cocoons emit quasi-isotropic GW emission in the LVK band, $\sim 10-100$ Hz, over a characteristic jet activity timescale, $\sim 10-100$ s. Our first-principles simulations show that jets exhibit a wobbling behavior, in which case cocoon-powered GWs might be detected already in LVK run A+, but it is more likely that these GWs will be detected by the third generation GW detectors with estimated rate of $ \sim 10 $ events/year. The detection rate drops to $ \sim 1\% $ of that value if all jets were to feature a traditional axisymmetric structure instead of a wobble. Accompanied by electromagnetic emission from the energetic core-collapse supernova and the cocoon, we predict that collapsars are powerful multi-messenger events., Comment: For movies of the simulation, see https://oregottlieb.com/gw.html

We study the influence of hot and dense matter in core-collapse supernovae by adopting up-to-date nuclear equation of state (EOS) based on the microscopic nuclear many-body frameworks. We explore effects of EOS based on the Dirac Brueckner Hartree-Fock theory through comparisons with those based on the variational method. We also examine effects of the differences in the composition of nuclei and nucleons by using the same EOS by the variational method but employing two different treatments in computations of nuclear abundances. We perform numerical simulations of core-collapse supernovae adopting the three EOSs. We also perform numerical simulations of the long-term evolution over 70 s of the proto-neutron star cooling. We show that impacts by different modeling of composition are remarkable as in those by different treatments of uniform matter in the gravitational collapse, bounce, and shock propagation. The cooling of proto-neutron star and the resulting neutrino emission are also affected by the compositional difference even if the same treatment in computing uniform matter of EOS., Comment: 33 pages, 23 figures, 1 table, accepted for publication in Progress of Theoretical and Experimental Physics (PTEP)

Neutrino self-interactions in a dense neutrino gas can induce collective neutrino flavor conversions. Fast neutrino flavor conversions (FFCs), one of the collective neutrino conversion modes, potentially change the dynamics and observables in core-collapse supernovae and binary neutron star mergers. In cases without neutrino-matter interactions (or collisions), FFCs are essentially energy-independent, and therefore the single energy treatment has been used in previous studies. However, neutrino-matter collisions in general depend on neutrino energy, suggesting that energy-dependent features may emerge in FFCs with collisions. In this paper, we perform dynamical simulations of FFCs with iso-energetic scatterings (emulating nucleon scatterings) under multi-energy treatment. We find that cancellation between in- and out-scatterings happens in high energy region, which effectively reduces the number of collisions and then affects the FFC dynamics. In fact, the lifetime of FFCs is extended compared to the single-energy case, leading to large flavor conversions. Our result suggests that the multi-energy treatment is mandatory to gauge the sensitivity of FFCs to collisions. We also provide a useful quantity to measure the importance of multi-energy effects of collisions on FFCs., Comment: 13 pages, 7 figures, accepted to PRD

Several recent numerical studies have examined the effects of neutrino neutral-current scattering on fast flavor conversion (FFC). Those studies are in apparent conflict, with some finding enhancement of flavor conversion and others finding suppression. The ones that report enhancement all use homogeneous models, and we discuss in detail the self-consistency issues that they face as a result. We reproduce scattering-enhanced FFC in our own homogeneous calculations, showing that it occurs for both neutral- and charged-current scattering and that it is due to isotropization of the angular distributions over time. Because this is the very feature of the calculations that is not self-consistent, we conclude that the enhancement effect is of unclear astrophysical relevance and may not occur in natural environments., Comment: 11 pages, 2 figures

Despite the theoretical indication that fast neutrino-flavor conversion (FFC) ubiquitously occurs in core-collapse supernova and binary neutron star merger, the lack of global simulations has been the greatest obstacle to study their astrophysical consequences. In this {\it Letter}, we present large-scale ($50 {\rm km}$) simulations of FFC in spherical symmetry by using a novel approach. We effectively rescale the oscillation scale of FFC by reducing the number of injected neutrinos in the simulation box, and then extrapolate back to the case of the target density of neutrinos with a convergence study. We find that FFC in all models achieves quasi-steady state in the non-linear regime, and its saturation property of FFC is universal. We also find that temporal- and spatial variations of FFC are smeared out at large radii due to phase cancellation through neutrino self-interactions. Finally, we provide a new diagnostic quantity, ELN-XLN angular crossing, to assess the non-linear saturation of FFC., Comment: Accepted for publication in PRL

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., Comment: Accepted to PRD

We investigate the protoneutron star (PNS) convection using our newly developed general relativistic Boltzmann neutrino radiation-hydrodynamics code. This is a pilot study for more comprehensive investigations later. As such, we take a snapshot of a PNS at 2.3 seconds after bounce from a 1D PNS cooling calculation and run our simulation for $\sim160\,\mathrm{ms}$ in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial condition is linearly unstable to convection. We find in our 2D simulation that convection is instigated there indeed and expands inward after being full-fledged. The convection is then settled to a quasi-steady state in $\sim100\,\mathrm{ms}$, being sustained by the negative $Y_e$ gradient, which is in turn maintained by neutrino emissions. It enhances the luminosities and mean energies of all species of neutrinos compared to 1D. Taking advantage of the Boltzmann solver, we analyze the possible occurrence of the neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in the regions, where $Y_e$ is lower, and that the growth rate can be as high as $\sim 10^{-1}\,{\mathrm{cm}^{-1}}$., Comment: accepted to ApJ

In this paper we present a novel method to estimate the time evolution of proto-neutron star (PNS) structure from the neutrino signal in core-collapse supernovae (CCSN). Employing recent results of multi-dimensional CCSN simulations, we delve into a relation between total emitted neutrino energy (TONE) and PNS mass/radius, and we find that they are strongly correlated with each other. We fit the relation by simple polynomial functions connecting TONE to PNS mass and radius as a function of time. By combining another fitting function representing the correlation between TONE and cumulative number of event at each neutrino observatory, PNS mass and radius can be retrieved from purely observed neutrino data. We demonstrate retrievals of PNS mass and radius from mock data of neutrino signal, and we assess the capability of our proposed method. While underlining the limitations of the method, we also discuss the importance of the joint analysis with gravitational wave signal. This would reduce uncertainties of parameter estimations in our method, and may narrow down the possible neutrino oscillation model. The proposed method is a very easy and inexpensive computation, which will be useful in real data analysis of CCSN neutrino signal., Comment: Accepted to MNRAS

There is mounting evidence that neutrinos undergo fast flavor conversion (FFC) in core-collapse supernova (CCSN). In this letter, 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 ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$). This has a significant impact on the angular distribution of neutrinos in momentum space, in which $\nu_{\rm e}$ tends to be more isotropic than $\bar{\nu}_{\rm 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., Comment: 8 pages, 5 figures, accepted for publication in the ApJ

Using full Boltzmann neutrino transport, we performed two-dimensional (2D) core-collapse supernova simulations in axisymmetry for two progenitor models with 11.2M and 15.0M both rotational and non-rotational. 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 (LB) 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., Comment: 28pages, 22figures, accepted for publication in the ApJ

Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings is a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the post-shock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. We also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star (PNS) convection, asymmetric neutrino emission, neutrino absorptions and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models., Comment: Accepted to PRD

Neutrinos have an unique quantum feature as flavor conversions. Recent studies suggested that collective neutrino oscillations play important roles in high-energy astrophysical phenomena. Quantum kinetic equation (QKE) is capable of describing the neutrino flavor conversion, transport and matter collision self-consistently. However, we have experienced many technical difficulties in their numerical implementation. In this paper, we present a new QKE solver based on Monte Carlo (MC) approach. This is an upgraded version of our classical MC neutrino transport solver; in essence, a flavor degree of freedom including mixing state is added into each MC particle. This extension requires updating numerical treatments of collision terms, in particular for scattering processes. We deal with the technical problem by generating a new MC particle at each scattering event. To reduce statistical noise inherent in MC methods, we develop the effective mean free path method. This suppresses a sudden change of flavor state due to collisions without increasing the number of MC particles. We present a suite of code tests to validate these new modules with comparing to the results reported in previous studies. Our QKE-MC solver is developed with fundamentally different philosophy and design from other deterministic- and mesh methods, suggesting that it will be complementary to others, and potentially provide new insights into physical processes of neutrino dynamics., Comment: 22 pages, 21 figures, Submitted to ApJS, we made minor revisions for V2 and V3

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., Comment: Minor update to match version published in PRD

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., Comment: Minor corrections were made in Figs. (Accepted to PRD)

With the recognition that fast flavor instabilities likely affect supernova and neutron-star-merger neutrinos, using simulation data to pin down when and where the instabilities occur has become a high priority. The effort faces an interesting problem. Fast instabilities are related to neutrino angular crossings, but simulations often employ moment methods, sacrificing momentum-space angular resolution in order to allocate resources elsewhere. How can limited angular information be used most productively? The main aims here are to sharpen this question and examine some of the available answers. A recently proposed method of searching for angular crossings is scrutinized, the limitations of moment closures are highlighted, and two ways of reconstructing angular distributions solely from the flux factors (based respectively on maximum-entropy and sharp-decoupling assumptions) are compared. In (semi)transparent regions, the standard closure prescriptions likely miss some crossings that should be there and introduce others that should not., Comment: 9 pages, 5 figures. reference added

We study theoretical neutrino signals from core-collapse supernova (CCSN) computed using axisymmetric CCSN simulations that cover the post-bounce phase up to $\sim 4$~s. We provide basic quantities of the neutrino signals such as event rates, energy spectra, and cumulative number of events at some terrestrial neutrino detectors, and then discuss some new features in the late phase that emerge in our models. Contrary to popular belief, neutrino emissions in the late phase are not always steady, but rather have temporal fluctuations, the vigor of which hinges on the CCSN model and neutrino flavor. We find that such temporal variations are not primarily driven by proto-neutron star (PNS) convection, but by fallback accretion in exploding models. We assess the detectability of these temporal variations, and find that IceCube is the most promising detector with which to resolve them. We also update fitting formulae first proposed in our previous paper for which the total neutrino energy (TONE) emitted at the CCSN source is estimated from the cumulative number of events in each detector. This will be a powerful technique with which to analyze real observations, particularly for low-statistics data., Comment: Matching the publish version

We investigate the impact of asymmetric neutrino-emissions on explosive nucleosynthesis in core-collapse supernovae (CCSNe) of progenitors with a mass range of 9.5 to 25$M_{\odot}$. We perform axisymmetric, hydrodynamic simulations of the CCSN explosion with a simplified neutrino-transport, in which anti-correlated dipolar emissions of $\nu_{\rm e}$ and ${\bar \nu}_{\rm e}$ are imposed. We then evaluate abundances and masses of the CCSN ejecta in a post-processing manner. We find that the asymmetric $\nu$-emission leads to the abundant ejection of $p$- and $n$-rich matter in the high-$\nu_{\rm e}$ and -${\bar \nu}_{\rm e}$ hemispheres, respectively. It substantially affects the abundances of the ejecta for elements heavier than Ni regardless of progenitors, although those elements lighter than Ca are less sensitive. Based on these results, we calculate the IMF-averaged abundances of the CCSN ejecta with taking into account the contribution from Type Ia SNe. For $m_{\rm asy} = 10/3\%$ and $10\%$, where $m_{\rm asy}$ denotes the asymmetric degree of the dipole components in the neutrino emissions, the averaged abundances for elements lighter than Y are comparable to those of the solar abundances, whereas those of elements heavier than Ge are overproduced in the case with $m_{\rm asy} \ge 30\%$. Our result also suggests that the effect of the asymmetric neutrino emissions is imprinted in the difference of abundance ratio of [Ni/Fe] and [Zn/Fe] between the high-$\nu_{\rm e}$ and -${\bar \nu}_{\rm e}$ hemispheres, indicating that the future spectroscopic X-ray observations of a CCSN remnant will bring evidence of the asymmetric neutrino emissions if exist., Comment: 14 pages, 10 figures, Accepted for the publication in MNRAS

We present a comprehensive study of neutrino shock acceleration in core-collapse supernova (CCSN). The leading players are heavy leptonic neutrinos, $\nu_{\mu}$ and $\nu_{\tau}$; the former and latter potentially gain the energy up to $\sim 100$ MeV and $\sim 200$ MeV, respectively, through the shock acceleration. Demonstrating the neutrino shock acceleration by Monte Carlo neutrino transport, we make a statement that it commonly occurs in the early post bounce phase ($\lesssim 50$ ms after bounce) for all massive stellar collapse experiencing nuclear bounce and would reoccur in the late phase ($\gtrsim 100$ ms) for failed CCSNe. This opens up a new possibility to detect high energy neutrinos by terrestrial detectors from Galactic CCSNe; hence, we estimate the event counts for Hyper(Super)-Kamiokande, DUNE, and JUNO. We find that the event count with the energy of $\gtrsim 80$ MeV is a few orders of magnitude higher than that of the thermal neutrinos regardless of the detectors, and muon production may also happen in these detectors by $\nu_{\mu}$ with the energy of $\gtrsim 100$ MeV. The neutrino signals provide a precious information on deciphering the inner dynamics of CCSN and placing a constraint on the physics of neutrino oscillation; indeed, the detection of the high energy neutrinos through charged current reaction channels will be a smoking gun evidence of neutrino flavor conversion., Comment: Accepted to MNRAS

We develop a neutrino transfer code for core-collapse simulations, that directly solves the multidimensional Boltzmann equations in full general relativity. We employ the discrete ordinate method, which discretizes the six dimensional phase space. The code is an extension of our special relativistic code coupled to a Newtonian hydrodynamics code, which is currently employed for core-collapse supernova simulations. In order to demonstrate our code's capability to treat general relativistic effects, we conduct some tests: we first compute the free streaming of neutrinos in the Schwarzschild and Kerr spacetimes and compare the results with the geodesic curves; in the Schwarzschild case we deploy not only a 1-dimensional grid in space under spherical symmetry but also a 2-dimensional spatial mesh under axisymmetry in order to assess the capability of the code to compute the spatial advection of neutrinos; secondly, we calculate the neutrino transport in a fixed matter background, which is taken from a core-collapse supernova simulation with our general relativistic but spherically symmetric Boltzmann-hydrodynamics code, to obtain a steady neutrino distribution; the results are compared with those given by the latter code., Comment: 23 pages, 19 figures, accepted version

Mounting evidence indicates that neutrinos likely undergo fast flavor conversion (FFC) in at least some core-collapse supernovae. Outcomes of FFC, however, remain highly uncertain. Here we study the cascade of flavor-field power from large angular scales in momentum space down to small ones, showing that FFC enhances this process and thereby hastens relaxation. Cascade also poses a computational challenge, which is present even if the flavor field is stable: When power reaches the smallest angular scale of the calculation, error from truncating the angular-moment expansion propagates back to larger scales, to disastrous effect on the overall evolution. Essentially the same issue has prompted extensive work in the context of plasma kinetics. This link suggests new approaches to averting spurious evolution, a problem that presently puts severe limitations on the feasibility of realistic oscillation calculations., Comment: 15 pages, 8 figures

We present a new method by which to retrieve energy spectrum for all flavor of neutrinos from core-collapse supernova (CCSN). In the retrieval process, we do not assume any analytic formulae to express the energy spectrum of neutrinos but rather take a direct way of spectrum reconstruction from the observed data; the Singular Value Decomposition algorithm with a newly developed adaptive energy-gridding technique is adopted. We employ three independent reaction channels having different flavor sensitivity to neutrinos. Two reaction channels, inverse beta decay on proton and elastic scattering on electrons, from a water Cherenkov detector such as Super-Kamiokande (SK) and Hyper-Kamiokande (HK), and a charged current reaction channel with Argon from the Deep Underground Neutrino Experiment (DUNE) are adopted. Given neutrino oscillation models, we iteratively search the neutrino energy spectra at the CCSN source until they provide the consistent event counts in the three reaction channels. We test the capability of our method by demonstrating the spectrum retrieval to a theoretical neutrino data computed by our recent three-dimensional CCSN simulation. Although the energy spectrum with either electron-type or electron-type anti-neutrinos at the CCSN source has relatively large error compared to that of other species, the joint analysis with HK + DUNE or SK + DUNE will provide precise energy spectrum of all flavors of neutrinos at the source. Finally, we discuss perspectives for improvements of our method by using neutrino data of other detectors., Comment: Accepted to MNRAS

Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($\gtrsim 50 \%$), but that the time-integrated signal has a relatively modest asymmetry ($\lesssim 20 \%$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $\sim$20\% if the distance to the CCSN is $\lesssim 6$ kpc., Comment: Accepted to MNRAS

We examine the sensitivity of neutrino emissions to stellar evolution models for a 15$M_\odot$ progenitor, paying particular attention to a phase prior to the collapse. We demonstrate that the number luminosities in both electron-type neutrinos ($\nu_e$) and their anti-partners ($\bar{\nu}_e$) differ by more than an order of magnitude by changing spatial resolutions and nuclear network sizes on stellar evolution models. We also develop a phenomenological model to capture the essential trend of the diversity, in which neutrino luminosities are expressed as a function of central density, temperature and electron fraction. In the analysis, we show that neutrino luminosity can be well characterized by these central quantities. This analysis also reveals that the most influential quantity to the time evolution of $\nu_e$ luminosity is matter density, while it is temperature for $\bar{\nu}_e$. These qualitative trends will be useful and applicable to constrain the physical state of progenitors at the final stages of stellar evolution from future neutrino observations, although more detailed systematic studies including various mass progenitors are required to assess the applicability., Comment: 13 pages, 12 figures, submitted to MNRAS

We report on the core-collapse supernova simulation we conducted for a 11.2 M 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., Comment: 29 pages, 22 figures, submitted to ApJ

Using the Boltzmann-radiation-hydrodynamics code, which solves the Boltzmann equation for neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer and Swesty (LS) and Furusawa and Shen (FS). We extend the simulation time of the LS model and conduct thorough investigations, though our previous paper briefly reported some of the results. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The protoneutron star seen in the FS model is more compact than that in the LS model because the existence of multinuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the motion of matter determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semitransparent region., Comment: 30 pages, 24 figures, 1 explosion, published in ApJ

The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods., Comment: 22 pages, 11 figures, submitted to ApJS

This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical Fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behavior emerges. We find that the vigor of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS $-$ heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of $\nu_{\mu}$, $\nu_{\tau}$, $\bar{\nu}_{\mu}$, and $\bar{\nu}_{\tau}$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields., Comment: Accepted to MNRAS