Supernova Dynamics Research Articles

Forbidden transitions in neutral- and charged-current interactions between low-energy neutrinos and argon

Background: The study of low-energy neutrinos and their interactions with atomic nuclei is crucial to several open problems in physics, including the neutrino mass hierarchy, CP-violation, candidates of Beyond Standard Model physics and supernova dynamics. Examples of experiments include CAPTAIN at SNS as well as DUNE's planned detection program of supernova neutrinos. Purpose: We present cross section calculations for quasielastic charged current and neutral current neutrinos at low energies, with a focus on $^{40}$Ar. We also take a close look at pion decay-at-rest neutrino spectra, which are used in e.g. the SNS experiment at Oakridge. Method and results: We employ a Hartree Fock + Continuum Random Phase Approximations (HF+CRPA) framework, which allows us to model the responses and include the effects of long-range correlations. It is expected to provide a good framework to calculate forbidden transitions, whose contribution which we show to be non-negligible. Conclusions: Forbidden transitions can be expected to contribute sizeably to the reaction strength at typical low-energy kinematics, such as DAR neutrinos. Modeling and Monte Carlo simulations need to take all due care to account for the influence of their contributions.

Sensitivity of nuclear statistical equilibrium to nuclear uncertainties during stellar core collapse

I have systematically investigated the equations of state (EOSs) in nuclear statistical equilibrium under thermodynamic conditions relevant for core collapse of massive stars by varying the bulk properties of nuclear matter, the mass data for neutron-rich nuclei, and the finite-temperature modifications of the nuclear model. It is found that the temperature dependence of the nuclear free energies has a significant impact on the entropy and nuclear composition, which affect the dynamics of core-collapse supernovae. There is a little influence from the bulk properties and the mass data. For all models, common nuclei that are likely to contribute to core-deleptonization are those near $Z\approx30$ and $N\approx50$. A model with a semi-empirical expression for internal degrees of freedom, however, overestimates the number densities of magic nuclei with $N\approx50$ and $82$, while a model, in which nuclear shell effects are not considered, underestimates the number densities of heavy nuclei, and especially of the magic nuclei. Other models, which include the temperature dependence of shell effects in the internal degrees of freedom and/or in the nuclear internal free energy, indicate that the difference in population between magic nuclei and non-magic nuclei disappears as the temperature increases. The construction of complete statistical EOS will require further theoretical and experimental studies of medium-mass, neutron-rich nuclei with proton numbers 25-45 and neutron numbers 40-85.

Simple method of diagnosing fast flavor conversions of supernova neutrinos

Neutrinos emitted from a supernova may undergo flavor conversions almost immediately above the core, with possible consequences for supernova dynamics and nucleosynthesis. However, the precise conditions for such fast conversions can be difficult to compute and require knowledge of the full angular distribution of the flavor-dependent neutrino fluxes, that is not available in typical supernova simulations. In this paper, we show that the overall flavor evolution is qualitatively similar to the growth of a so-called `zero mode', determined by the background matter and neutrino densities, which can be reliably predicted using only the second angular moments of the electron lepton number distribution, i.e., the difference in the angular distributions of $\nu_e$ and $\bar{\nu}_e$ fluxes. We propose that this zero mode, which neither requires computing the full Green's function nor a detailed knowledge of the angular distributions, may be useful for a preliminary diagnosis of possible fast flavor conversions in supernova simulations with modestly resolved angular distributions

R-process Nucleosynthesis from Three-dimensional Magnetorotational Core-collapse Supernovae

Abstract We investigate r-process nucleosynthesis in 3D general-relativistic magnetohydrodynamic simulations of rapidly rotating strongly magnetized core collapse. The simulations include a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to postbounce neutrino emission and absorption. We track the composition of the ejected material using the nuclear reaction network SkyNet. Our results show that the 3D dynamics of magnetorotational core-collapse supernovae (CCSN) are important for their nucleosynthetic signature. We find that production of r-process material beyond the second peak is reduced by a factor of 100 when the magnetorotational jets produced by the rapidly rotating core undergo a kink instability. Our results indicate that 3D magnetorotationally powered CCSNe are robust r-process sources only if they are obtained by the collapse of cores with unrealistically large precollapse magnetic fields of the order of 1013 G. Additionally, a comparison simulation that we restrict to axisymmetry results in overly optimistic r-process production for lower magnetic field strengths.

Neutrino signals of core-collapse supernovae in underground detectors

For a suite of fourteen core-collapse models during the dynamical first second after bounce, we calculate the detailed neutrino "light" curves expected in the underground neutrino observatories Super-Kamiokande, DUNE, JUNO, and IceCube. These results are given as a function of neutrino-oscillation modality (normal or inverted hierarchy) and progenitor mass (specifically, post-bounce accretion history), and illuminate the differences between the light curves for 1D (spherical) models that don't explode with the corresponding 2D (axisymmetric) models that do. We are able to identify clear signatures of explosion (or non-explosion), the post-bounce accretion phase, and the accretion of the silicon/oxygen interface. In addition, we are able to estimate the supernova detection ranges for various physical diagnostics and the distances out to which various temporal features embedded in the light curves might be discerned. We find that the progenitor mass density profile and supernova dynamics during the dynamical explosion stage should be identifiable for a supernova throughout most of the galaxy in all the facilities studied and that detection by any one of them, but in particular more than one in concert, will speak volumes about the internal dynamics of supernovae.

Neutrino-nucleus reactions and their role in supernova dynamics and nucleosynthesis

Neutrino reactions on nuclei play important roles for the dynamics of supernovae and their associated nucleosynthesis. This manuscript summarizes the current status in deriving the relevant cross sections for supernova neutrinos and briefly discusses a few recent advances in supernova simulations where these reactions play a role.

Investigation of Gamow Teller transition properties in [formula omitted]Ni isotopes using QRPA methods

Weak rates in nickel isotopes play an integral role in the dynamics of supernovae. Electron capture and β-decay of nickel isotopes, dictated by Gamow–Teller transitions, significantly alter the lepton fraction of the stellar matter. In this paper we calculate Gamow–Teller (GT) transitions for isotopes of nickel, Ni56–64, using QRPA methods. The GT strength distributions were calculated using four different QRPA models. Our results are also compared with previous theoretical calculations and measured strength distributions wherever available. Our investigation concluded that amongst all RPA models, the pn-QRPA(C) model best described the measured GT distributions (including total GT strength and centroid placement). It is hoped that the current investigation of GT properties would prove handy and may lead to a better understanding of the presupernova evolution of massive stars.

How to form a millisecond magnetar? Magnetic field amplification in protoneutron stars

AbstractExtremely strong magnetic fields of the order of 1015G are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

Probing nonlinear electrodynamics in slowly rotating spacetimes through neutrino astrophysics

Huge electromagnetic fields are known to be present during the late stages of the dynamics of supernovae. Thus, when dealing with electrodynamics in this context, the possibility may arise to probe nonlinear theories (generalizations of the Maxwellian electromagnetism). We firstly solve Einstein field equations minimally coupled to an arbitrary (current-free) nonlinear Lagrangian of electrodynamics (NLED) in the slow rotation regime $a\ll M$ (black hole's mass), up to first order in $a/M$. We then make use of the robust and self-contained Born-Infeld Lagrangian in order to compare and contrast the physical properties of such NLED spacetime with its Maxwellian counterpart (a slowly rotating Kerr-Newman spacetime), especially focusing on the astrophysics of both neutrino flavor oscillations ($\nu_e \rightarrow \nu_\mu, \nu_\tau$) and spin-flip ($\nu_l \rightarrow \nu_r$, "$l$" stands for "left" and "$r$" stands for "right", change of neutrino handedness) mass level-crossings, the equivalent to gyroscopic precessions. Such analysis proves that in the spacetime of a slowly rotating nonlinear charged black hole (RNCBH), intrinsically associated with the assumption the electromagnetism is nonlinear, the neutrino dynamics in core-collapse supernovae could be significantly changed. In such astrophysical environment a positive enhancement (reduction of the electron fraction $Y_e<0.5$) of the r-process may take place. Consequently, it might result in hyperluminous supernova explosions due to enlargement, in atomic number and amount, of the decaying nuclides. Finally, we envisage some physical scenarios that may lead to short-lived charged black holes with high charge-to-mass ratios (associated with unstable highly magnetized neutron stars) and ways to possibly disentangle theories of the electromagnetism from other black holes observables (by means of light polarization measurements).

Direct URCA Processes in Supernovae and Accretion Disks with Arbitrary Magnetic Field

An effect of a magnetic field of an arbitrary strength on the beta-decay and reactions related with it by the crossing symmetry (the beta-processes) in supernovae and accretion disks around black holes is analyzed. Rates of the beta-processes and the energy and momentum transfered through them to an optically transparent matter are calculated. It is shown that the macroscopic momentum transferred to the medium increases linearly with the magnetic field strength and can substantially affect the dynamics of supernovae and accretion disks especially when a matter inside is degenerate. It is also demonstrated that the rates of the beta-processes and the energy deposition in these reactions for the magnetic field strength $B \lesssim 10^{15}$ G, which is assumed to be typical in supernovae and accretion disks, are lower than in the absence of the field. This suppression is more pronounced for reactions with neutrinos.

SHOULD ONE USE THE RAY-BY-RAY APPROXIMATION IN CORE-COLLAPSE SUPERNOVA SIMULATIONS?

ABSTRACT We perform the first self-consistent, time-dependent, multi-group calculations in two dimensions (2D) to address the consequences of using the ray-by-ray+ transport simplification in core-collapse supernova simulations. Such a dimensional reduction is employed by many researchers to facilitate their resource-intensive calculations. Our new code (Fornax) implements multi-D transport, and can, by zeroing out transverse flux terms, emulate the ray-by-ray+ scheme. Using the same microphysics, initial models, resolution, and code, we compare the results of simulating 12, 15, 20, and 25 M ⊙ progenitor models using these two transport methods. Our findings call into question the wisdom of the pervasive use of the ray-by-ray+ approach. Employing it leads to maximum post-bounce/pre-explosion shock radii that are almost universally larger by tens of kilometers than those derived using the more accurate scheme, typically leaving the post-bounce matter less bound and artificially more “explodable.” In fact, for our 25 M ⊙ progenitor, the ray-by-ray+ model explodes, while the corresponding multi-D transport model does not. Therefore, in two dimensions, the combination of ray-by-ray+ with the axial sloshing hydrodynamics that is a feature of 2D supernova dynamics can result in quantitatively, and perhaps qualitatively, incorrect results.

Effect of magnetic field on beta processes in a relativistic moderately degenerate plasma

The effect of a magnetic field of arbitrary strength on the beta decay and crossing symmetric processes is analyzed. A covariant calculation technique is used to derive the expression for the squares of S-matrix elements of these reactions, which is also valid in reference frames in which the medium moves as a single whole along magnetic field lines. Simple analytic expressions obtained for the neutrino and antineutrino emissivities for a moderately degenerate plasma fully characterize the emissivity and absorbability of the studied medium. It is shown that the approximation used here is valid for core collapse supernovae and accretion disks around black holes; beta processes in these objects are predominantly neutrino reactions. The analytic expressions obtained for the emissivities can serve as a good approximation for describing the interaction of electron neutrinos and antineutrinos with the medium of the objects in question and hold for an arbitrary magnetic field strength. Due to their simplicity, these expressions can be included in the magnetohydrodynamic simulation of supernovae and accretion disks to calculate neutrino and antineutrino transport in them. The rates of beta processes and the energy and momentum emitted in them are calculated for an optically transparent matter. It is shown that the macroscopic momentum transferred in the medium increases linearly with the magnetic field strength and can substantially affect the dynamics of supernovae and accretion disks in the regions of a degenerate matter. It is also shown that the rates of beta processes and the energy emission for a magnetic field strength of B ≲ 1015 G typical of supernovae and accretion disks are lower than in the absence of field. This suppression is stronger for reactions with neutrinos.

The role of neutrino-nucleus reactions in supernova dynamics and nucleosynthesis

Neutrino reactions on nuclei play important roles for the dynamics of supernovae and their associated nucleosynthesis. This manuscript summarizes the current status in deriving the relevant cross sections for supernova neutrinos and briefly discusses a few recent advances where

Constraining the supersaturation density equation of state from core-collapse supernova simulations?

In this article the role of the supersaturation density equation of state (EOS) is explored in simulations of failed core-collapse supernova explosions. Therefore the nuclear EOS is extended via a one-parameter excluded volume description for baryons, taking into account their finite and increasing volume with increasing density in excess of saturation density. Parameters are selected such that the resulting supernova EOS represent extreme cases, with high pressure variations at supersaturation density which feature extreme stiff and soft EOS variants of the reference case, i.e. without excluded volume corrections. Unlike in the interior of neutron stars with central densities in excess of several times saturation density, central densities of core-collapse supernovae reach only slightly above saturation density. Hence, the impact of the supersaturation density EOS on the supernova dynamics as well as the neutrino signal is found to be negligible. It is mainly determined from the low- and intermediate-density domain, which is left unmodified within this generalized excluded volume approach.

THE EVOLUTION AND IMPACTS OF MAGNETOROTATIONAL INSTABILITY IN MAGNETIZED CORE-COLLAPSE SUPERNOVAE

ABSTRACT We carried out two-dimensional axisymmetric MHD simulations of core-collapse supernovae for rapidly rotating magnetized progenitors. By changing both the strength of the magnetic field and the spatial resolution, the evolution of the magnetorotational instability (MRI) and its impacts upon the dynamics are investigated. We found that the MRI greatly amplifies the seed magnetic fields in the regime where the buoyant mode, not the Alfvén mode, plays a primary role in the exponential growth phase. The MRI indeed has a powerful impact on the supernova dynamics. It makes the shock expansion faster and the explosion more energetic, with some models being accompanied by the collimated jet formations. These effects, however, are not made by the magnetic pressure except for the collimated jet formations. The angular momentum transfer induced by the MRI causes the expansion of the heating region, by which the accreting matter gain additional time to be heated by neutrinos. The MRI also drifts low-Y p matter from deep inside of the core to the heating region, which makes the net neutrino heating rate larger by the reduction of the cooling due to the electron capture. These two effects enhance the efficiency of the neutrino heating, which is found to be the key to boosting the explosion. Indeed, we found that our models explode far more weakly when the net neutrino heating is switched off. The contribution of the neutrino heating to the explosion energy could reach 60% even in the case of strongest magnetic field in the current simulations.

Expected impact from weak reactions with light nuclei in corecollapse supernova simulations

We study the role of light nuclear clusters in simulations of core-collapse supernovae. Expressions for the reaction rates are developed for a large selection of charged current absorption and scattering processes with light clusters. Medium modifications are taken into account at the mean-field level. We explore the possible impact on the supernova dynamics and the neutrino signal during the mass accretion phase prior to the possible explosion onset as well as during the subsequent protoneutron star deleptnoization after the explosion onset has been launched.

Design, characterization, and sensitivity of the supernova trigger system at Daya Bay

Providing an early warning of galactic supernova explosions from neutrino signals is important in studying supernova dynamics and neutrino physics. A dedicated supernova trigger system has been designed and installed in the data acquisition system at Daya Bay and integrated into the worldwide Supernova Early Warning System (SNEWS). Daya Bay’s unique feature of eight identically-designed detectors deployed in three separate experimental halls makes the trigger system naturally robust against cosmogenic backgrounds, enabling a prompt analysis of online triggers and a tight control of the false-alert rate. The trigger system is estimated to be fully sensitive to 1987A-type supernova bursts throughout most of the Milky Way. The significant gain in sensitivity of the eight-detector configuration over a mass-equivalent single detector is also estimated. The experience of this online trigger system is applicable to future projects with spatially distributed detectors.

FAIR - Cosmic Matter in the Laboratory

To explore cosmic matter in the laboratory - this fascinating research prospect becomes available at the Facility for Antiproton and Ion Research, FAIR. The new facility is being constructed within the next five years adjacent to the existing accelerator complex of the GSI Helmholtz Centre for Heavy Ion Research at Darmstadt/Germany, expanding the research goals and technical possibilities substantially. This includes new insights into the dynamics of supernovae depending on the properties of short-lived neutron-rich nuclei which will be investigated with intense rare isotope beams. New insights will be provided into the interior of stars by exploring dense plasmas with intense heavy-ion beams combined with a high-performance laser - or into neutron star cores by probing the highest baryon densities in relativistic nucleus-nucleus collisions at unprecedented collision rates. To the latter, the properties of hadrons play an important part which will be systematically studied by high precision hadron spectroscopy with antiproton beams at unmatched intensities. The worldwide unique accelerator and experimental facilities of FAIR will open the way for a broad spectrum of unprecedented fore-front research supplying a large variety of experiments in hadron, nuclear, atomic and plasma physics as well as biomedical and material science which will be briefly described in this article. This article is based on the FAIR Green Paper [4] and gives an update of former publications [5] - [12].

Highlight on Supernova Early Warning at Daya Bay

Providing an early warning of supernova burst neutrinos is of importance in studying both supernova dynamics and neutrino physics. The Daya Bay Reactor Neutrino Experiment, with a unique feature of multiple liquid scintillator detectors, is sensitive to the full energy spectrum of supernova burst electron-antineutrinos. By utilizing 8 Antineutrino Detectors (ADs) in the three different experimental halls which are about 1 km's apart from each other, we obtain a powerful and prompt rejection of muon spallation background than single-detector experiments with the same target volume. A dedicated trigger system embedded in the data acquisition system has been installed to allow the detection of a coincidence of neutrino signals of all ADs via an inverse beta-decay (IBD) within a 10-second window, thus providing a robust early warning of a supernova occurrence within the Milky Way. An 8-AD associated supernova trigger table has been established theoretically to tabulate the 8-AD event counts’ coincidence vs. the trigger rate. As a result, a golden trigger threshold, i.e. with a false alarm rate < 1/3-months, can be set as low as 6 candidates among the 8 detectors, leading to a 100% detection probability for all 1987A type supernova bursts at the distance to the Milky Way center and a 96% detection probability to those at the edge of the Milky Way.

Constraining the12C+12C fusion cross section for astrophysics

The 12 C+12 C reaction is one of the single most important nuclear reactions in astrophysics. It strongly influences late evolution of massive stars as well as the dynamics of type Ia supernovae and x-ray superbursts. An accurate estimation of the cross section at relevant astrophysical energies is extremely important for modeling these systems. However, the situation is complicated by the unpredictable resonance structure observed at higher energies. Two recent studies at Notre Dame have produced results which help reduce the uncertainty associated with this reaction. The first uses correlations with the isotope fusion systems, 12 C+13 C and 13 C+13 C, to establish an upper limit on the resonance strengths in 12 C+12 C. The other focuses on the specific channel 12 C+12 C→23 Mg+n and its low-energy measurement and extrapolation which is relevant to s-process nucleosynthesis. The results from each provide important constraints for astrophysical models.