Search

Multi-messenger astrophysics has produced a wealth of data with much more to come in the future. This enormous data set will reveal new insights into the physics of core-collapse supernovae, neutron star mergers, and many other objects where it is actually possible, if not probable, that new physics is in operation. To tease out different possibilities, we will need to analyze signals from photons, neutrinos, gravitational waves, and chemical elements. This task is made all the more difficult when it is necessary to evolve the neutrino component of the radiation field and associated quantum-mechanical property of flavor in order to model the astrophysical system of interest—a numerical challenge that has not been addressed to this day. In this work, we take a step in this direction by adopting the technique of angular-integrated moments with a truncated tower of dynamical equations and a closure, convolving the flavor-transformation with spatial transport to evolve the neutrino radiation quantum field. We show that moments capture the dynamical features of fast flavor instabilities in a variety of systems, although our technique is by no means a universal blueprint for solving fast flavor transformation. To evaluate the effectiveness of our moment results, we compare to a more precise particle-in-cell method. Based on our results, we propose areas for improvement and application to complementary techniques in the future. [ABSTRACT FROM AUTHOR]

Core-collapse supernovae (SNe) are one of the most powerful cosmic sources of neutrinos, with energies of several MeV. The emission of neutrinos and antineutrinos of all flavors carries away the gravitational binding energy of the compact remnant and drives its evolution from the hot initial to the cold final states. Detecting these neutrinos from Earth and analyzing the emitted signals present a unique opportunity to explore the neutrino mass ordering problem. This research outlines the detection of neutrinos from SNe and their relevance in understanding the neutrino mass ordering. The focus is on developing a model-independent analysis strategy, achieved by comparing distinct detection channels in large underground detectors. The objective is to identify potential indicators of mass ordering within the neutrino sector. Additionally, a thorough statistical analysis is performed on the anticipated neutrino signals for both mass orderings. Despite uncertainties in supernova explosion parameters, an exploration of the parameter space reveals an extensive array of models with significant sensitivity to differentiate between mass orderings. The assessment of various observables and their combinations underscores the potential of forthcoming supernova observations in addressing the neutrino mass ordering problem. [ABSTRACT FROM AUTHOR]

Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars or black holes) and the synthesis and injection of many heavy elements into their host galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of one-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infrared light curves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine. Join us in improving and evolving this document through active community engagement. Instructions are provided at this link: https://github.com/clfryer/MM-SNe. [ABSTRACT FROM AUTHOR]

Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost region of the supernova, where neutrinos are expected to interact, allows us to examine the presence of neutrino interactions. Here we demonstrate that the amount of Mn produced without neutrino nucleosynthesis processes (i.e., the ν - and νp -processes) is too small to explain the Mn/Fe mass ratio we measure (0.14%–0.67%). This result supports the operation of significant neutrino interactions in the Cassiopeia A supernova. If the observed Mn/Fe mass ratio purely reflects the production at the innermost region of the supernova, this would be the first robust confirmation of neutrino–matter interactions in an individual supernova. We further show that the Mn/Fe mass ratio has the potential to constrain supernova neutrino parameters (i.e., total neutrino luminosity, neutrino temperature). Future spatially resolved, high-resolution X-ray spectroscopy will allow us to investigate the details of neutrino–supernova astrophysics through its signatures in elemental composition not only in Cassiopeia A but also in other remnants. [ABSTRACT FROM AUTHOR]

The Spallation Neutron Source (SNS) at Oak Ridge National Laboratory provides an intense, high-quality source of neutrinos from pion decay at rest. This source was recently used for the first measurements of coherent elastic neutrino–nucleus scattering (CEvNS) by the COHERENT Collaboration, which resulted in new constraints of physics beyond the Standard Model. The SNS neutrino source will enable further CEvNS measurements, exploration of inelastic neutrino–nucleus interactions of particular relevance for understanding supernova neutrinos, and searches for accelerator-produced sub-GeV dark matter. Taking advantage of this unique facility, COHERENT's suite of detectors in Neutrino Alley at the SNS is accumulating more data to address a broad physics program at the intersection of particle physics, nuclear physics, and astrophysics. This review describes COHERENT's first two CEvNS measurements, their interpretation, and the potential of a future physics program at the SNS. [ABSTRACT FROM AUTHOR]

Among multimessenger observations of the next Galactic core-collapse supernova, Super-Kamiokande (SK) plays a critical role in detecting the emitted supernova neutrinos, determining the direction to the supernova (SN), and notifying the astronomical community of these observations in advance of the optical signal. In 2022, SK has increased the gadolinium dissolved in its water target (SK-Gd) and has achieved a Gd concentration of 0.033%, resulting in enhanced neutron detection capability, which in turn enables more accurate determination of the supernova direction. Accordingly, SK-Gd’s real-time supernova monitoring system has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on 2021 December 13, and is available through GCN Notices. When the monitoring system detects an SN-like burst of events, SK_SN Notice will automatically distribute an alarm with the reconstructed direction to the supernova candidate within a few minutes. In this paper, we present a systematic study of SK-Gd’s response to a simulated Galactic SN. Assuming a supernova situated at 10 kpc, neutrino fluxes from six supernova models are used to characterize SK-Gd’s pointing accuracy using the same tools as the online monitoring system. The pointing accuracy is found to vary from 3° to 7° depending on the models. However, if the supernova is closer than 10 kpc, SK_SN Notice can issue an alarm with three-degree accuracy, which will benefit follow-up observations by optical telescopes with large fields of view.

When a star undergoes core collapse, a vast amount of energy is released in a ∼10 s long burst of neutrinos of all species. Inverse beta decay in the star’s hydrogen envelope causes an electromagnetic cascade that ultimately results in a flare of gamma rays—an “echo” of the neutrino burst—at the characteristic energy of 0.511 MeV. We study the phenomenology and detectability of this flare. Its luminosity curve is characterized by a fast, seconds-long rise and an equally fast decline, with a minute- or hour-long plateau in between. For a near-Earth star (distance D ≲ 1 kpc) the echo will be observable at near future gamma-ray telescopes with an effective area of 10 ^3 cm ^2 or larger. Its observation will inform us on the envelope size and composition. In conjunction with the direct detection of the neutrino burst, it will also give information on the neutrino emission away from the line of sight and will enable tests of neutrino propagation effects between the stellar surface and Earth.

In high-energy astrophysical processes involving compact objects, such as core-collapse supernovae or binary neutron star mergers, neutrinos play an important role in the synthesis of nuclides. Neutrinos in these environments can experience collective flavor oscillations driven by neutrino–neutrino interactions, including coherent forward scattering and incoherent (collisional) effects. Recently, there has been interest in exploring potential novel behaviors in collective oscillations of neutrinos by going beyond the one-particle effective or “mean-field” treatments. Here, we seek to explore implications of collective neutrino oscillations, in the mean-field treatment and beyond, for the nucleosynthesis yields in supernova environments with different astrophysical conditions and neutrino inputs. We find that collective oscillations can impact the operation of the ν p -process and r -process nucleosynthesis in supernovae. The potential impact is particularly strong in high-entropy, proton-rich conditions, where we find that neutrino interactions can nudge an initial ν p -process neutron-rich, resulting in a unique combination of proton-rich low-mass nuclei as well as neutron-rich high-mass nuclei. We describe this neutrino-induced neutron-capture process as the “ ν i -process.” In addition, nontrivial quantum correlations among neutrinos, if present significantly, could lead to different nuclide yields compared to the corresponding mean-field oscillation treatments, by virtue of modifying the evolution of the relevant one-body neutrino observables.

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

If a Galactic core-collapse supernova explosion occurs in the future, it will be critical to rapidly alert the community to the direction of the supernova by utilizing neutrino signals in order to enable the initiation of follow-up optical observations. In addition, there is anticipation that observation of the diffuse supernova neutrino background will yield discoveries in the near future, given that experimental upper limits are approaching theoretical predictions. We have developed a new supernova event simulator for water Cherenkov neutrino detectors, such as the highly sensitive Super-Kamiokande. This simulator calculates the neutrino interaction in water for two simulation purposes, individual core-collapse supernova bursts and diffuse supernova neutrino background. Based on this simulator, we can evaluate the precision in determining the location of supernovae and estimate the expected number of events related to the diffuse supernova neutrino background in Super-Kamiokande. In this paper, we describe the basic structure of the simulator and its demonstration.

The occurrence of a first-order hadron–quark matter phase transition at high baryon densities is investigated in astrophysical simulations of core-collapse supernovae, to decipher yet incompletely understood properties of the dense matter equation of state (EOS) using neutrinos from such cosmic events. It is found that the emission of a nonstandard second neutrino burst, dominated by electron antineutrinos, is not only a measurable signal for the appearance of deconfined quark matter but also reveals information about the state of matter at extreme conditions encountered at the supernova (SN) interior. To this end, a large set of spherically symmetric SN models is investigated, studying the dependence on the EOS and the stellar progenitor. General relativistic neutrino-radiation hydrodynamics is employed featuring three-flavor Boltzmann neutrino transport and a microscopic hadron-quark hybrid matter EOS class. Therefore, the DD2 relativistic mean-field hadronic model is employed, and several variations of it, and the string-flip model for the description of deconfined quark matter. The resulting hybrid model covers a representative range of onset densities for the phase transition and latent heats. This facilitates the direct connection between intrinsic signatures of the neutrino signal and properties of the EOS. In particular, a set of linear relations has been found empirically. These potentially provide a constraint for the onset density of a possible QCD phase transition from the future neutrino observation of the next galactic core-collapse SN, if a millisecond electron anti-neutrino burst is present around or less than 1 s.

Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s ^−1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s ^−1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M _⊙ , is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M _⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

Multi-messenger astrophysics has produced a wealth of data with much more to come in the future. This enormous data set will reveal new insights into the physics of core-collapse supernovae, neutron star mergers, and many other objects where it is actually possible, if not probable, that new physics is in operation. To tease out different possibilities, we will need to analyze signals from photons, neutrinos, gravitational waves, and chemical elements. This task is made all the more difficult when it is necessary to evolve the neutrino component of the radiation field and associated quantum-mechanical property of flavor in order to model the astrophysical system of interest—a numerical challenge that has not been addressed to this day. In this work, we take a step in this direction by adopting the technique of angular-integrated moments with a truncated tower of dynamical equations and a closure, convolving the flavor-transformation with spatial transport to evolve the neutrino radiation quantum field. We show that moments capture the dynamical features of fast flavor instabilities in a variety of systems, although our technique is by no means a universal blueprint for solving fast flavor transformation. To evaluate the effectiveness of our moment results, we compare to a more precise particle-in-cell method. Based on our results, we propose areas for improvement and application to complementary techniques in the future.

The IceCube Neutrino Observatory has been continuously taking data to search for ${ \mathcal O }(0.5\mbox{--}10)$ s long neutrino bursts since 2007. Even if a Galactic core-collapse supernova is optically obscured or collapses to a black hole instead of exploding, it will be detectable via the ${ \mathcal O }(10)$ MeV neutrino burst emitted during the collapse. We discuss a search for such events covering the time between 2008 April 17 and 2019 December 31. Considering the average data taking and analysis uptime of 91.7% after all selection cuts, this is equivalent to 10.735 yr of continuous data taking. In order to test the most conservative neutrino production scenario, the selection cuts were optimized for a model based on an 8.8 solar mass progenitor collapsing to an O–Ne–Mg core. Conservative assumptions on the effects of neutrino oscillations in the exploding star were made. The final selection cut was set to ensure that the probability to detect such a supernova within the Milky Way exceeds 99%. No such neutrino burst was found in the data after performing a blind analysis. Hence, a 90% C.L. upper limit on the rate of core-collapse supernovae out to distances of ≈25 kpc was determined to be 0.23 yr ^−1 . For the more distant Magellanic Clouds, only high neutrino luminosity supernovae will be detectable by IceCube, unless external information on the burst time is available. We determined a model-independent limit by parameterizing the dependence on the neutrino luminosity and the energy spectrum.

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.

We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/Si–O interface, the star will likely explode. We develop a quantitative criterion by using ∼1300 1D simulations where ν -driven turbulence is included via time-dependent mixing-length theory. This criterion correctly identifies the outcome of the supernova more than 90% of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that ν -driven turbulence plays in the explosion by comparing our results to previous works. [ABSTRACT FROM AUTHOR]

The stellar initial mass function (IMF) is expressed by ϕ (m) ∝ m − α with the slope α, and known as a poorly constrained but very important function in studies of star and galaxy formation. There are no sensible observational constraints on the IMF slopes beyond the Milky Way and nearby galaxies. Here we combine two sets of observational results, (1) cosmic densities of core-collapse supernova (CCSN) explosion rates and (2) cosmic far-UV radiation (and infrared reradiation) densities, which are sensitive to massive (≃8–50 M ⊙) and moderately massive (≃2.5–7 M ⊙) stars, respectively, and constrain the IMF slope at m > 1 M ⊙ with a freedom of redshift evolution. Although no redshift evolution is identified beyond the uncertainties, we find that the cosmic-average IMF slope at z = 0 is α = 1.8–3.2 at the 95% confidence level that is comparable with the Salpeter IMF, α = 2.35, which marks the first constraint on the cosmic-average IMF. We show a forecast for the Nancy Grace Roman Space Telescope supernova survey that will provide significantly strong constraints on the IMF slope with δ α ≃ 0.5 over z = 0–2. Moreover, as for an independent IMF probe instead of (1), we suggest to use diffuse supernovae neutrino background (DSNB), relic neutrinos from CCSNe. We expect that the Hyper-Kamiokande neutrino observations over 20 yr will improve the constraints on the IMF slope and the redshift evolution significantly better than those obtained today, if the systematic uncertainties of DSNB production physics are reduced in the future numerical simulations. [ABSTRACT FROM AUTHOR]

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova (CCSN) will be preceded by an increasing flux of neutrinos and antineutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent CCSN, hours earlier than the detection of the neutrinos from core collapse. Electron antineutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximizing detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hr before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, α-Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hr before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code F ornax. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better. [ABSTRACT FROM AUTHOR]

In this paper, we study the evolution of random hypermagnetic fields (HMFs) in the symmetric phase of the early universe before the electroweak phase transition. The behavior of HMFs is driven by the analog of the chiral magnetic effect accounting for the asymmetries of leptons and Higgs bosons. These asymmetries are also dynamical variables of the model and evolve together with HMFs. Moreover, we account for the contribution of the hyper-MHD turbulence in the effective diffusion coefficient and the α -dynamo parameter. The realistic spectrum of seed HMFs consists of two branches: Batchelor and Kolmogorov ones. The impact of HMFs on the production of relic gravitational waves (GWs) and the baryon asymmetry of the universe (BAU), as well as flavor oscillations of supernova neutrinos in the stochastic GWs generated is considered. We establish the constraint on the strength of the seed HMF comparing the spectral density of produced GWs with the observations of the LIGO-Virgo-KAGRA collaborations. The stronger upper bound on the seed HMF is obtained from the condition of not exceeding the observed value of BAU. [ABSTRACT FROM AUTHOR]

Core-collapse supernovae (SNe) are one of the most powerful cosmic sources of neutrinos, with energies of several MeV. The emission of neutrinos and antineutrinos of all flavors carries away the gravitational binding energy of the compact remnant and drives its evolution from the hot initial to the cold final states. Detecting these neutrinos from Earth and analyzing the emitted signals present a unique opportunity to explore the neutrino mass ordering problem. This research outlines the detection of neutrinos from SNe and their relevance in understanding the neutrino mass ordering. The focus is on developing a model-independent analysis strategy, achieved by comparing distinct detection channels in large underground detectors. The objective is to identify potential indicators of mass ordering within the neutrino sector. Additionally, a thorough statistical analysis is performed on the anticipated neutrino signals for both mass orderings. Despite uncertainties in supernova explosion parameters, an exploration of the parameter space reveals an extensive array of models with significant sensitivity to differentiate between mass orderings. The assessment of various observables and their combinations underscores the potential of forthcoming supernova observations in addressing the neutrino mass ordering problem.

Growing evidence from multiwavelength observations of extragalactic supernovae (SNe) has established the presence of dense circumstellar material in Type II SNe. Interaction between the SN ejecta and the circumstellar material should lead to diffusive shock acceleration of cosmic rays and associated high-energy emission. Observation of high-energy neutrinos along with the MeV neutrinos from SNe will provide unprecedented opportunities to understand unanswered questions in cosmic-ray and neutrino physics. We show that current and future neutrino detectors can identify high-energy neutrinos from an extragalactic SN in the neighborhood of the Milky Way. We present the prospects for detecting high-energy neutrino minibursts from SNe in known local galaxies, and demonstrate how the network of multiple high-energy neutrino detectors will extend the horizon for the identification of high-energy SN neutrinos. We also discuss high-energy neutrino emission from SN 2023ixf.

Supernova (SN) 2023ixf in M101 is the closest SN explosion observed in the last decade. Therefore, it is a suitable test bed to study the role of jets in powering the SN ejecta. With this aim, we explored the idea that high-energy neutrinos could be produced during the interaction between the jets and the intense radiation field produced in the SN explosion and eventually be observed by the IceCube neutrino telescope. The lack of detection of such neutrinos has significantly constrained both the fraction of stellar collapses that produce jets and/or the theoretical models for neutrino production. Finally, we investigated the possibility of detecting low-energy neutrinos from SN 2023ixf with the Super- and Hyper-Kamiokande experiments, obtaining, in both cases, subthreshold estimates.

We report the first search result for the flux of astrophysical electron antineutrinos for energies ${ \mathcal O }(10)\,\mathrm{MeV}$ in the gadolinium-loaded Super-Kamiokande (SK) detector. In 2020 June, gadolinium was introduced to the ultrapure water of the SK detector in order to detect neutrons more efficiently. In this new experimental phase, SK-Gd, we can search for electron antineutrinos via inverse beta decay with efficient background rejection thanks to the high efficiency of the neutron tagging technique. In this paper, we report the result for the initial stage of SK-Gd, during 2020 August 26, and 2022 June 1 with a 22.5 × 552 kton · day exposure at 0.01% Gd mass concentration. No significant excess over the expected background in the observed events is found for the neutrino energies below 31.3 MeV. Thus, the flux upper limits are placed at the 90% confidence level. The limits and sensitivities are already comparable with the previous SK result with pure water (22.5 × 2970 kton · day) owing to the enhanced neutron tagging. Operation with Gd increased to 0.03% started in 2022 June.

Core-collapse supernova explosions play a wide role in astrophysics by producing compact remnants (neutron stars or black holes) and the synthesis and injection of many heavy elements into their host galaxy. Because they are produced in some of the most extreme conditions in the universe, they can also probe physics in extreme conditions (matter at nuclear densities and extreme temperatures and magnetic fields). To quantify the impact of supernovae on both fundamental physics and our understanding of the universe, we must leverage a broad set of observables of this engine. In this paper, we study a subset of these probes using a suite of one-dimensional, parameterized mixing models: ejecta remnants from supernovae, ultraviolet, optical and infrared light curves, and transient gamma-ray emission. We review the other diagnostics and show how the different probes tie together to provide a more clear picture of the supernova engine. Join us in improving and evolving this document through active community engagement. Instructions are provided at this link: https://github.com/clfryer/MM-SNe .

Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost region of the supernova, where neutrinos are expected to interact, allows us to examine the presence of neutrino interactions. Here we demonstrate that the amount of Mn produced without neutrino nucleosynthesis processes (i.e., the ν - and νp -processes) is too small to explain the Mn/Fe mass ratio we measure (0.14%–0.67%). This result supports the operation of significant neutrino interactions in the Cassiopeia A supernova. If the observed Mn/Fe mass ratio purely reflects the production at the innermost region of the supernova, this would be the first robust confirmation of neutrino–matter interactions in an individual supernova. We further show that the Mn/Fe mass ratio has the potential to constrain supernova neutrino parameters (i.e., total neutrino luminosity, neutrino temperature). Future spatially resolved, high-resolution X-ray spectroscopy will allow us to investigate the details of neutrino–supernova astrophysics through its signatures in elemental composition not only in Cassiopeia A but also in other remnants.

Supernova neutrinos are invaluable signals that offer information about the interior of supernovae. Because a nearby supernova can occur at any time, preparing for future supernova neutrino observation is an urgent task. For the prompt analysis of supernova neutrinos, we have developed a new analysis code, the “Supernova Parameter Estimation Code based on Insight on Analytic Late-time Burst Light curve at Earth Neutrino Detector” ( SPECIAL BLEND ). This code estimates the parameters of supernovae based on an analytic model of supernova neutrinos from the proto-neutron star cooling phase. For easy availability to the community, this code is public and easily runs in web environments. SPECIAL BLEND can estimate the parameters better than the analysis pipeline we developed in a previous paper. By using SPECIAL BLEND , we can estimate the supernova parameters within 10% precision up to ∼20 and ∼60 kpc (Large Magellanic Cloud contained) with Super-Kamiokande and Hyper-Kamiokande, respectively.

Fluxes of the diffuse supernova neutrino background (DSNB) are calculated based on a new modeling of galactic chemical evolution, where a variable stellar initial mass function (IMF), depending on the galaxy type, is introduced and black hole (BH) formation from the failed supernova is considered for progenitors heavier than 18 M _⊙ . The flux calculations are performed for different combinations of the star formation rate, nuclear equation of state, and neutrino mass hierarchy, to examine the systematic effects from these factors. In any case, our new model predicts the enhanced DSNB ${\bar{\nu }}_{e}$ flux at E _ν ≳ 30 MeV and E _ν ≲ 10 MeV, due to more frequent BH formation and a larger core-collapse rate at high redshifts in early-type galaxies, respectively. Event rate spectra of the DSNB ${\bar{\nu }}_{e}$ at a detector from the new model are shown, and the detectability at water-based Cherenkov detectors, Super-Kamiokande with a gadolinium dissolution and Hyper-Kamiokande, is discussed. In order to investigate the impacts of the assumptions in the new model, we prepare alternative models, based on different IMF forms and treatments of BH formation, and estimate the discrimination capabilities between the new and alternative models at these detectors.

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

We present a simple criterion to predict the explodability of massive stars based on the density and entropy profiles before collapse. If a pronounced density jump is present near the Si/Si–O interface, the star will likely explode. We develop a quantitative criterion by using ∼1300 1D simulations where ν -driven turbulence is included via time-dependent mixing-length theory. This criterion correctly identifies the outcome of the supernova more than 90% of the time. We also find no difference in how this criterion performs on two different sets of progenitors, evolved using two different stellar evolution codes: FRANEC and KEPLER. The explodability as a function of mass of the two sets of progenitors is very different, showing: (i) that uncertainties in the stellar evolution prescriptions influence the predictions of supernova explosions; (ii) the most important properties of the pre-collapse progenitor that influence the explodability are its density and entropy profiles. We highlight the importance that ν -driven turbulence plays in the explosion by comparing our results to previous works.

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code F ornax . We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

We investigate 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 s after bounce from a 1D PNS cooling calculation and run our simulation for ∼160 ms in 2D under axisymmetry. The original PNS cooling calculation neglected convection entirely and the initial conditions were 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 then settled to a quasi-steady state after ∼100 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 analyse the possible occurrence of neutrino fast flavor conversion (FFC). We found that FFC is likely to occur in regions where Y _e is lower, and that the growth rate can be as high as ∼10 ^−1 cm ^−1 .

Large neutrino detectors like IceCube monitor for core-collapse supernovae using low-energy (MeV) neutrinos with a detection reach from a supernova neutrino burst to the Magellanic Cloud. However, some models predict the emission of high-energy neutrinos of GeV–TeV from core-collapse supernovae through the interaction of ejecta with circumstellar material with energies of TeV–PeV produced through choked jets. In this paper, we explore the detection horizon of IceCube for core-collapse supernovae using high-energy neutrinos from these models. We examine the potential of two high-energy neutrino data samples from IceCube, one that performs best in the northern sky and one that has better sensitivity in the southern sky. We demonstrate that, by using high-energy neutrinos from core-collapse supernovae, the detection reach can be extended to the megaparsec range, far beyond what is accessible through low-energy neutrinos. Looking ahead to IceCube-Gen2, this reach will be extended considerably.

The stellar initial mass function (IMF) is expressed by ϕ ( m ) ∝ m ^− ^α with the slope α , and known as a poorly constrained but very important function in studies of star and galaxy formation. There are no sensible observational constraints on the IMF slopes beyond the Milky Way and nearby galaxies. Here we combine two sets of observational results, (1) cosmic densities of core-collapse supernova (CCSN) explosion rates and (2) cosmic far-UV radiation (and infrared reradiation) densities, which are sensitive to massive (≃8–50 M _⊙ ) and moderately massive (≃2.5–7 M _⊙ ) stars, respectively, and constrain the IMF slope at m > 1 M _⊙ with a freedom of redshift evolution. Although no redshift evolution is identified beyond the uncertainties, we find that the cosmic-average IMF slope at z = 0 is α = 1.8–3.2 at the 95% confidence level that is comparable with the Salpeter IMF, α = 2.35, which marks the first constraint on the cosmic-average IMF. We show a forecast for the Nancy Grace Roman Space Telescope supernova survey that will provide significantly strong constraints on the IMF slope with δ α ≃ 0.5 over z = 0–2. Moreover, as for an independent IMF probe instead of (1), we suggest to use diffuse supernovae neutrino background (DSNB), relic neutrinos from CCSNe. We expect that the Hyper-Kamiokande neutrino observations over 20 yr will improve the constraints on the IMF slope and the redshift evolution significantly better than those obtained today, if the systematic uncertainties of DSNB production physics are reduced in the future numerical simulations.

Dark Matter particles can be detected directly via their elastic scattering with nuclei. Next generation experiments can eventually find physical evidences about dark matter candidates. With this motivation in mind, we have calculated the expected signals of dark matter particles in xenon detectors. The calculations were performed by considering different masses and parameters within the minimal supersymmetric standard model. Since the detectors can also detect neutrinos, we have analyzed the supernova neutrino signal including a sterile neutrino in the formalism. Using this 3 + 1 scheme, we make predictions for both the normal and inverse mass hierarchy. In order to perform a study of the response of planned direct-detection experiments, to be located in ANDES (Agua Negra Deep Experimental Site), we have calculated the neutrino contributions to the background by taken into account reactor's neutrinos and geoneutrinos at the site of the lab. As a test detector, we take a Xenon1T-like array. [ABSTRACT FROM AUTHOR]

Current neutrino detectors will observe hundreds to thousands of neutrinos from Galactic supernovae, and future detectors will increase this yield by an order of magnitude or more. With such a data set comes the potential for a huge increase in our understanding of the explosions of massive stars, nuclear physics under extreme conditions, and the properties of the neutrino. However, there is currently a large gap between supernova simulations and the corresponding signals in neutrino detectors, which will make any comparison between theory and observation very difficult. SNEWPY is an open-source software package that bridges this gap. The SNEWPY code can interface with supernova simulation data to generate from the model either a time series of neutrino spectral fluences at Earth, or the total time-integrated spectral fluence. Data from several hundred simulations of core-collapse, thermonuclear, and pair-instability supernovae is included in the package. This output may then be used by an event generator such as sntools or an event rate calculator such as the SuperNova Observatories with General Long Baseline Experiment Simulator (SNOwGLoBES). Additional routines in the SNEWPY package automate the processing of the generated data through the SNOwGLoBES software and collate its output into the observable channels of each detector. In this paper we describe the contents of the package, the physics behind SNEWPY, the organization of the code, and provide examples of how to make use of its capabilities.

We briefly review the recent developments in neutrino physics and astrophysics which have import for frontline research in nuclear physics. These developments, we argue, tie nuclear physics to exciting developments in observational cosmology and astrophysics in new ways. Moreover, the behavior of neutrinos in dense matter is itself a fundamental problem in many-body quantum mechanics, in some ways akin to well-known issues in nuclear matter and nuclei, and in some ways radically different, especially because of nonlinearity and quantum de-coherence. The self-interacting neutrino gas is the only many body system driven by the weak interactions. © 2013 Elsevier B.V. All rights reserved.

Neutrino-nucleus reaction cross sections on 18O are evaluated by shell-model calculations and compared with those on 16O. Important contributions from Gamow-Teller transitions are noticed for 18O (ν e , e−) 18F in contrary to the case for 16O, where spin-dipole transitions are dominant contributions. Calculated cross sections for 18O (ν e , e−) 18F are shown to be larger than for 16O at low neutrino energies below 20 MeV in natural water with the 0.205% admixture of 18O due to the lower threshold energy (1.66 MeV) for 18O than that for 16O (15.42 MeV). The resulting electron spectra, that is, the cross sections as functions of emitted electron energy T e , are also shown to be quite different, reflecting the different threshold energies. The electron spectra from (ν e , e−) reactions on 18O and 16O in water Cherenkov detectors for supernova neutrino detection are investigated for both the cases with and without the neutrino oscillation and compared with those of the neutrino-electron scattering. It has been shown that the contribution from 18O (0.205% mixture) enhances the rates from 16O by 60% for the case without the oscillation and by 20-30% for the case with the oscillation below T e =20 MeV. For the case with the neutrino oscillation, the event rates for 18O and 16O become comparable to those of the neutrino-electron scattering. However, their rates at low energy (T e < 20 MeV) are much lower than those of the neutrino-electron scattering, which is important for the pointing accuracy to the supernova direction. [ABSTRACT FROM AUTHOR]

The next core-collapse supernova in the Milky Way or its satellites will represent a once-in-a-generation opportunity to obtain detailed information about the explosion of a star and provide significant scientific insight for a variety of fields because of the extreme conditions found within. Supernovae in our galaxy are not only rare on a human timescale but also happen at unscheduled times, so it is crucial to be ready and use all available instruments to capture all possible information from the event. The first indication of a potential stellar explosion will be the arrival of a bright burst of neutrinos. Its observation by multiple detectors worldwide can provide an early warning for the subsequent electromagnetic fireworks, as well as signal to other detectors with significant backgrounds so they can store their recent data. The supernova early warning system (SNEWS) has been operating as a simple coincidence between neutrino experiments in automated mode since 2005. In the current era of multi-messenger astronomy there are new opportunities for SNEWS to optimize sensitivity to science from the next galactic supernova beyond the simple early alert. This document is the product of a workshop in June 2019 towards design of SNEWS 2.0, an upgraded SNEWS with enhanced capabilities exploiting the unique advantages of prompt neutrino detection to maximize the science gained from such a valuable event.

Abstract Neutrino–nuclear responses associated with astro-neutrinos, single beta decays and double beta decays are crucial in studies of neutrino properties of interest for astro-particle physics. The present report reviews briefly recent studies of the neutrino–nuclear responses from both experimental and theoretical points of view in order to obtain a consistent understanding of the many facets of the neutrino–nuclear responses. Subjects discussed in this review include (i) experimental studies of neutrino–nuclear responses by means of single beta decays, charge-exchange nuclear reactions, muon- photon- and neutrino–nuclear reactions, and nucleon-transfer reactions, (ii) implications of and discussions on neutrino–nuclear responses for single beta decays, for astro-neutrinos, and for astro-neutrino nucleosynthesis, (iii) theoretical aspects of neutrino–nuclear responses for beta and double beta decays, for nuclear muon capture and for neutrino–nucleus scattering, and (iv) critical discussions on nucleonic and non-nucleonic spin–isospin correlations and renormalization (quenching or enhancement) effects on the axial weak coupling. Remarks are given on perspectives of experimental and theoretical studies of the neutrino–nuclear responses and on future experiments of double beta decays. [ABSTRACT FROM AUTHOR]

Abstract Ton-scale direct dark matter search experiments should be sensitive to neutrino-induced recoil events from either 8B solar neutrinos or the brief but intense flux from a core collapse supernova in the Milky Way. These low-threshold detectors are sensitive to the very low recoil energies, of order 10 keV, deposited via coherent elastic scatters between supernova neutrinos and target nuclei. Large superheated fluid detectors like PICO-500, a bubble chamber to be initially filled with an active target of 1 t of C 3 F 8 , should see multiple-bubble events from CE ν NS if the detector is live during a neutrino burst from a supernova at a distance up to 10 kpc. This paper discusses conditions under which bubble chambers could be used as an independent measurement in the event of a supernova similar to SN 1987A, with particular sensitivity to the currently less-constrained heavy-lepton ν x channel. [ABSTRACT FROM AUTHOR]