Proto-neutron Star Convection Research Articles

Supernova Explosions of the Lowest-mass Massive Star Progenitors

We here focus on the behavior of supernovae that technically explode in 1D (spherical symmetry). When simulated in 3D, however, the outcomes of representative progenitors of this class are quite different in almost all relevant quantities. In 3D, the explosion energies can be 2 to 10 times higher, and there are correspondingly large differences in the 56Ni yields. These differences between the 3D and 1D simulations reflect in part the relative delay to explosion of the latter and in the former the presence of protoneutron star convection that boosts the driving neutrino luminosities by as much as ∼50% at later times. In addition, we find that the ejecta in 3D models are more neutron-rich, resulting in significant weak r-process and 48Ca yields. Furthermore, we find that in 3D the core is an interesting, though subdominant, source of acoustic power. In summary, we find that though a model might be found theoretically to explode in 1D, one must perform supernova simulations in 3D to capture most of the associated observables. The differences between 1D and 3D models are just too large to ignore.

Proto-neutron star convection and the neutrino-driven wind: implications for the νp-process

ABSTRACT Recent studies of the neutrino-driven wind from proto-neutron stars have indicated that the wind is likely proton rich for much of its lifetime, and the high flux of neutrinos can induce νp-process nucleosynthesis allowing for the formation of heavy elements. It has also been shown that gravito-acoustic waves, generated by convection within the proto-neutron star, can significantly alter the dynamics and nucleosynthesis in the wind. Therefore, we present a study of the effects of convection-driven waves on the nucleosynthesis in proton-rich neutrino-driven winds, focusing on the νp-process. We find that wave effects can strongly impact νp-process nucleosynthesis even at wave luminosities a factor of 10−5 smaller than the total neutrino luminosity. The momentum flux of the waves accelerates the wind, reducing the net neutrino heating and the persistent neutron abundance created by p($\bar{\nu }_\mathrm{ e},e^+$), which impedes νp-process nucleosynthesis. However, this effect is generally counteracted by the effects of waves on seed nucleus formation, as the acceleration of the wind and the heating that occurs as these waves shock both favour a more α-rich environment with very little heavy seed nucleus formation. Overall, higher wave luminosities correlate (albeit non-monotonically) with heavier element νp-processing, up to A ≈ 200 in some cases. At very high wave luminosities (≳10−3Lν), early shock heating by the waves disrupts α recombination, and drives a suppressed, fast-outflow r-process proceeding up to A ≈ 200. This occurs despite an assumed neutrino spectrum that predicts a proton-rich wind with equilibrium Ye = 0.6.

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25M⊙

We present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and $25{M}_{\ensuremath{\bigodot}}$, are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model is distinct from the temporal evolution of our solar metallicity $15{M}_{\ensuremath{\bigodot}}$ progenitor model and our zero metallicity $25{M}_{\ensuremath{\bigodot}}$ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period $\ensuremath{\sim}75\text{ }\text{ }\mathrm{ms}$ after bounce, whereas in the latter two cases high-frequency emission does not commence until $\ensuremath{\sim}125\text{ }\text{ }\mathrm{ms}$ after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below $\ensuremath{\sim}250\text{ }\text{ }\mathrm{Hz}$, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at $\ensuremath{\sim}125$, $\ensuremath{\sim}500$, and $\ensuremath{\sim}250\text{ }\text{ }\mathrm{ms}$ after bounce in the 9.6, 15, and $25{M}_{\ensuremath{\bigodot}}$ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below $\ensuremath{\sim}10\text{ }\text{ }\mathrm{Hz}$. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5--10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today's detectors. Finally, in our $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model, we see very high-frequency gravitational radiation, extending up to $\ensuremath{\sim}2000\text{ }\text{ }\mathrm{Hz}$. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

Proto-neutron star convection and the neutrino-driven wind: implications for the r-process

ABSTRACT The neutrino-driven wind from proto-neutron stars is a proposed site for r-process nucleosynthesis, although most previous work has found that a wind heated only by neutrinos cannot produce the third r-process peak. However, several groups have noted that introducing a secondary heating source within the wind can change the hydrodynamic conditions sufficiently for a strong r-process to proceed. One possible secondary heating source is gravito-acoustic waves, generated by convection inside the proto-neutron star. As these waves propagate into the wind, they can both accelerate the wind and shock and deposit energy into the wind. Additionally, the acceleration of the wind by these waves can reduce the total number of neutrino captures and thereby reduce the final electron fraction of the wind. In neutron rich conditions, all of these effects can make conditions more favourable for r-process nucleosynthesis. Here, we present a systematic investigation of the impact of these convection-generated gravito-acoustic waves within the wind on potential nucleosynthesis. We find that wave effects in the wind can generate conditions favourable for a strong r-process, even when the energy flux in the waves is a factor of 10−4 smaller than the total neutrino energy flux and the wind is marginally neutron rich. Nevertheless, this depends strongly on the radius at which the waves become non-linear and form shocks. We also find that both entropy production after shock formation and the acceleration of the wind due to stresses produced by the waves prior to shock formation impact the structure and nucleosynthesis of these winds.

Protoneutron Star Convection Simulated with a New General Relativistic Boltzmann Neutrino Radiation Hydrodynamics Code

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.

Effect of the Nuclear Equation of State on Relativistic Turbulence-induced Core-collapse Supernovae

Abstract The nuclear equation of state (EOS) is an important component in the evolution of core-collapse supernovae. In this paper we make a survey of various EOSs in the literature and analyze their effect on spherical core-collapse models in which the effects of three-dimensional turbulence is modeled by a general relativistic formulation of Supernova Turbulence In Reduced-dimensionality (STIR). We show that the viability of the explosion is quite EOS dependent and that it best correlates with the early-time interior entropy density of the proto–neutron star. We check that this result is not progenitor dependent, although the lowest-mass progenitors show different explosion properties, due to the different pre-collapse nuclear composition. Larger central entropies also induce more vigorous proto–neutron star convection in our one-dimensional turbulence model, as well as a wider convective layer.

Convection and Dynamo in Newly Born Neutron Stars

Abstract To study properties of magnetohydrodynamic (MHD) convection and resultant dynamo activities in proto-neutron stars (PNSs), we construct a “PNS in a box” simulation model and solve the compressible MHD equation coupled with a nuclear equation of state (EOS) and simplified leptonic transport. As a demonstration, we apply it to two types of PNS model with different internal structures: a fully convective model and a spherical-shell convection model. By varying the spin rate of the models, the rotational dependence of convection and the dynamo that operate inside the PNS is investigated. We find that, as a consequence of turbulent transport by rotating stratified convection, large-scale structures of flow and thermodynamic fields are developed in all models. Depending on the spin rate and the depth of the convection zone, various profiles of the large-scale structures are obtained, which can be physically understood as steady-state solutions to the “mean-field” equation of motion. Additionally to those hydrodynamic structures, a large-scale magnetic component of  ( 10 15 ) G is also spontaneously organized in disordered tangled magnetic fields in all models. The higher the spin rate, the stronger the large-scale magnetic component grows. Intriguingly, as an overall trend, the fully convective models have a stronger large-scale magnetic component than that in the spherical-shell convection models. The deeper the convection zone extends, the larger the size of the convective eddies becomes. As a result, rotationally constrained convection seems to be more easily achieved in the fully convective model, resulting in a higher efficiency of the large-scale dynamo there. To gain a better understanding of the origin of the diversity of a neutron star’s magnetic field, we need to study the PNS dynamo in a wider parameter range.

Equation-of-state Dependence of Gravitational Waves in Core-collapse Supernovae

Abstract Gravitational waves (GWs) provide unobscured insight into the birthplace of neutron stars and black holes in core-collapse supernovae (CCSNe). The nuclear equation of state (EOS) describing these dense environments is yet uncertain, and variations in its prescription affect the proto−neutron star (PNS) and the post-bounce dynamics in CCSN simulations, subsequently impacting the GW emission. We perform axisymmetric simulations of CCSNe with Skyrme-type EOSs to study how the GW signal and PNS convection zone are impacted by two experimentally accessible EOS parameters, (1) the effective mass of nucleons, m ⋆, which is crucial in setting the thermal dependence of the EOS, and (2) the isoscalar incompressibility modulus, K sat. While K sat shows little impact, the peak frequency of the GWs has a strong effective mass dependence due to faster contraction of the PNS for higher values of m ⋆ owing to a decreased thermal pressure. These more compact PNSs also exhibit more neutrino heating, which drives earlier explosions and correlates with the GW amplitude via accretion plumes striking the PNS, exciting the oscillations. We investigate the spatial origin of the GWs and show the agreement between a frequency-radial distribution of the GW emission and a perturbation analysis. We do not rule out overshoot from below via PNS convection as another moderately strong excitation mechanism in our simulations. We also study the combined effect of effective mass and rotation. In all our simulations we find evidence for a power gap near ∼1250 Hz; we investigate its origin and report its EOS dependence.

Gravitational wave signature of proto-neutron star convection: I. MHD numerical simulations

ABSTRACT Gravitational waves provide a unique and powerful opportunity to constrain the dynamics in the interior of proto-neutron stars during core collapse supernovae. Convective motions play an important role in generating neutron stars magnetic fields, which could explain magnetar formation in the presence of fast rotation. We compute the gravitational wave emission from proto-neutron star convection and its associated dynamo, by post-processing three-dimensional MHD simulations of a model restricted to the convective zone in the anelastic approximation. We consider two different proto-neutron star structures representative of early times (with a convective layer) and late times (when the star is almost entirely convective). In the slow rotation regime, the gravitational wave emission follows a broad spectrum peaking at about three times the turnover frequency. In this regime, the inclusion of magnetic fields slightly decreases the amplitude without changing the spectrum significantly compared to a non-magnetized simulation. Fast rotation changes both the amplitude and spectrum dramatically. The amplitude is increased by a factor of up to a few thousands. The spectrum is characterized by several peaks associated with inertial modes, whose frequency scales with the rotation frequency. Using simple physical arguments, we derive scalings that reproduce quantitatively several aspects of these numerical results. We also observe an excess of low-frequency gravitational waves, which appears at the transition to a strong field dynamo characterized by a strong axisymmetric toroidal magnetic field. This signature of dynamo action could be used to constrain the dynamo efficiency in a proto-neutron star with future gravitational wave detections.

Where, when, and why: Occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

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.

Supernova neutrino signals based on long-term axisymmetric simulations

ABSTRACT We study theoretical neutrino signals from core-collapse supernova (CCSN) computed using axisymmetric CCSN simulations that cover the post-bounce phase up to ∼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 vigour of which hinges on the CCSN model and neutrino flavour. We find that such temporal variations are not primarily driven by proto-neutron star 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 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 analyse real observations, particularly for low-statistics data.

Magnetorotational Explosion of a Massive Star Supported by Neutrino Heating in General Relativistic Three-dimensional Simulations

Abstract We present results of three-dimensional (3D), radiation-magnetohydrodynamics (MHD) simulations of core-collapse supernovae in full general relativity (GR) with spectral neutrino transport. In order to study the effects of the progenitor’s rotation and magnetic fields, we compute three models, where the precollapse rotation rate and magnetic fields are included parametrically to a 20 M ⊙ star. While we find no shock revival in our two nonmagnetized models during our simulation times (∼500 ms after bounce), the magnetorotational (MR) driven shock expansion immediately initiates after bounce in our rapidly rotating and strongly magnetized model. We show that the expansion of the MR-driven flows toward the polar directions is predominantly driven by the magnetic pressure, whereas the shock expansion toward the equatorial direction is supported by neutrino heating. Our detailed analysis indicates that the growth of the so-called kink instability may hinder the collimation of jets, resulting in the formation of broader outflows. Furthermore, we find a dipole emission of lepton number, only in the MR explosion model, whose asymmetry is consistent with the explosion morphology. Although it is similar to the lepton number emission self-sustained asymmetry (LESA), our analysis shows that the dipole emission occurs not from the proto–neutron star convection zone but from above the neutrino sphere, indicating that it is not associated with the LESA. We also report several unique neutrino signatures, which are significantly dependent on both the time and the viewing angle, if observed, possibly providing rich information regarding the onset of the MR-driven explosion.

A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations

ABSTRACT 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 behaviour emerges. We find that the vigour 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 νμ, ντ, $\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.

Wave heating from proto-neutron star convection and the core-collapse supernova explosion mechanism

ABSTRACT Our understanding of the core-collapse supernova explosion mechanism is incomplete. While the favoured scenario is delayed revival of the stalled shock by neutrino heating, it is difficult to reliably compute explosion outcomes and energies, which depend sensitively on the complex radiation hydrodynamics of the post-shock region. The dynamics of the (non-)explosion depend sensitively on how energy is transported from inside and near the proto-neutron star (PNS) to material just behind the supernova shock. Although most of the PNS energy is lost in the form of neutrinos, hydrodynamic and hydromagnetic waves can also carry energy from the PNS to the shock. We show that gravity waves excited by core PNS convection can couple with outgoing acoustic waves that present an appreciable source of energy and pressure in the post-shock region. Using one-dimensional simulations, we estimate the gravity wave energy flux excited by PNS convection and the fraction of this energy transmitted upwards to the post-shock region as acoustic waves. We find wave energy fluxes near $10^{51}\, \mathrm{erg}\, \mathrm{s}^{-1}\,$ are likely to persist for $\sim \! 1\, \mathrm{s}$ post-bounce. The wave pressure on the shock may exceed $10{{\ \rm per\ cent}}$ of the thermal pressure, potentially contributing to shock revival and, subsequently, a successful and energetic explosion. We also discuss how future simulations can better capture the effects of waves, and more accurately quantify wave heating rates.

Gravitational waves from 3D core-collapse supernova models: The impact of moderate progenitor rotation

We present predictions for the gravitational-wave (GW) emission of three-dimensional supernova (SN) simulations performed for a 15 solar-mass progenitor with the Prometheus-Vertex code using energy-dependent, three-flavor neutrino transport. The progenitor adopted from stellar evolution calculations including magnetic fields had a fairly low specific angular momentum (j_Fe <~ 10^{15} cm^2/s) in the iron core (central angular velocity ~0.2 rad/s), which we compared to simulations without rotation and with artificially enhanced rotation (j_Fe <~ 2*10^{16} cm^2/s; central angular velocity ~0.5 rad/s). Our results confirm that the time-domain GW signals of SNe are stochastic, but possess deterministic components with characteristic patterns at low frequencies (<~200 Hz), caused by mass motions due to the standing accretion shock instability (SASI), and at high frequencies, associated with gravity-mode oscillations in the surface layer of the proto-neutron star (PNS). Non-radial mass motions in the post-shock layer as well as PNS convection are important triggers of GW emission, whose amplitude scales with the power of the hydrodynamic flows. There is no monotonic increase of the GW amplitude with rotation, but a clear correlation with the strength of SASI activity. Our slowly rotating model is a fainter GW emitter than the non-rotating model because of weaker SASI activity and damped convection in the post-shock layer and PNS. In contrast, the faster rotating model exhibits a powerful SASI spiral mode during its transition to explosion, producing the highest GW amplitudes with a distinctive drift of the low-frequency emission peak from ~80-100 Hz to ~40-50 Hz. This migration signifies shock expansion, whereas non-exploding models are discriminated by the opposite trend.

The impact of different neutrino transport methods on multidimensional core-collapse supernova simulations

Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the FLASH simulation framework. We compare a two-moment M1 scheme with an analytic closure (M1), the isotropic diffusion source approximation (IDSA), and the advanced spectral leakage method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a ∼10%–20% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations.

Electron-capture and Low-mass Iron-core-collapse Supernovae: New Neutrino-radiation-hydrodynamics Simulations

Abstract We present new 1D (spherical) and 2D (axisymmetric) simulations of electron-capture (EC) and low-mass iron-core-collapse supernovae (SN). We consider six progenitor models: the ECSN progenitor from Nomoto; two ECSN-like low-mass low-metallicity iron-core progenitors from A. Heger (2016, private communication); and the 9, 10, and 11 (zero-age main-sequence) progenitors from Sukhbold et al. We confirm that the ECSN and ESCN-like progenitors explode easily even in 1D with explosion energies of up to a 0.15 Bethes ( ), and are a viable mechanism for the production of very-low-mass neutron stars. However, the 9, 10, and 11 progenitors do not explode in 1D and are not even necessarily easier to explode than higher-mass progenitor stars in 2D. We study the effect of perturbations and of changes to the microphysics and we find that relatively small changes can result in qualitatively different outcomes, even in 1D, for models sufficiently close to the explosion threshold. Finally, we revisit the impact of convection below the protoneutron star (PNS) surface. We analyze 1D and 2D evolutions of PNSs subject to the same boundary conditions. We find that the impact of PNS convection has been underestimated in previous studies and could result in an increase of the neutrino luminosity by up to factors of two.

Gravitational wave signals from 3D neutrino hydrodynamics simulations of core-collapse supernovae

We present gravitational wave (GW) signal predictions from four 3D multi-group neutrino hydrodynamics simulations of core-collapse supernovae of progenitors with 11.2 Msun, 20 Msun, and 27 Msun. GW emission in the pre-explosion phase strongly depends on whether the post-shock flow is dominated by the standing accretion shock instability (SASI) or convection and differs considerably from 2D models. SASI activity produces a strong signal component below 250 Hz through asymmetric mass motions in the gain layer and a non-resonant coupling to the proto-neutron star (PNS). Both convection- and SASI-dominated models show GW emission above 250 Hz, but with considerably lower amplitudes than in 2D. This is due to a different excitation mechanism for high-frequency l=2 motions in the PNS surface, which are predominantly excited by PNS convection in 3D. Resonant excitation of high-frequency surface g-modes in 3D by mass motions in the gain layer is suppressed compared to 2D because of smaller downflow velocities and a lack of high-frequency variability in the downflows. In the exploding 20 Msun model, shock revival results in enhanced low-frequency emission due to a change of the preferred scale of the convective eddies in the PNS convection zone. Estimates of the expected excess power in two frequency bands suggests that second-generation detectors will only be able to detect very nearby events, but that third-generation detectors could distinguish SASI- and convection-dominated models at distances of ~10 kpc.

SELF-SUSTAINED ASYMMETRY OF LEPTON-NUMBER EMISSION: A NEW PHENOMENON DURING THE SUPERNOVA SHOCK-ACCRETION PHASE IN THREE DIMENSIONS

During the stalled-shock phase of our 3D hydrodynamical core-collapse simulations with energy-dependent, 3-flavor neutrino transport, the lepton-number flux (nue minus antinue) emerges predominantly in one hemisphere. This novel, spherical-symmetry breaking neutrino-hydrodynamical instability is termed LESA for "Lepton-number Emission Self-sustained Asymmetry." While the individual nue and antinue fluxes show a pronounced dipole pattern, the heavy-flavor neutrino fluxes and the overall luminosity are almost spherically symmetric. Initially, LESA seems to develop stochastically from convective fluctuations, it exists for hundreds of milliseconds or more, and it persists during violent shock sloshing associated with the standing accretion shock instability. The nue minus antinue flux asymmetry originates mainly below the neutrinosphere in a region of pronounced proto-neutron star (PNS) convection, which is stronger in the hemisphere of enhanced lepton-number flux. On this side of the PNS, the mass-accretion rate of lepton-rich matter is larger, amplifying the lepton-emission asymmetry, because the spherical stellar infall deflects on a dipolar deformation of the stalled shock. The increased shock radius in the hemisphere of less mass accretion and minimal lepton-number flux (antinue flux maximum) is sustained by stronger convection on this side, which is boosted by stronger neutrino heating because the average antinue energy is higher than the average nue energy. Asymmetric heating thus supports the global deformation despite extremely nonstationary convective overturn behind the shock. While these different elements of LESA form a consistent picture, a full understanding remains elusive at present. There may be important implications for neutrino-flavor oscillations, the neutron-to-proton ratio in the neutrino-heated supernova ejecta, and neutron-star kicks, which remain to be explored.

A NEW MULTI-DIMENSIONAL GENERAL RELATIVISTIC NEUTRINO HYDRODYNAMICS CODE OF CORE-COLLAPSE SUPERNOVAE. III. GRAVITATIONAL WAVE SIGNALS FROM SUPERNOVA EXPLOSION MODELS

We present a detailed theoretical analysis of the gravitational-wave (GW) signal of the post-bounce evolution of core-collapse supernovae (SNe), employing for the first time relativistic, two-dimensional (2D) explosion models with multi-group, three-flavor neutrino transport based on the ray-by-ray-plus approximation. The waveforms reflect the accelerated mass motions associated with the characteristic evolutionary stages that were also identified in previous works: A quasi-periodic modulation by prompt postshock convection is followed by a phase of relative quiescence before growing amplitudes signal violent hydrodynamical activity due to convection and the standing accretion shock instability during the accretion period of the stalled shock. Finally, a high-frequency, low-amplitude variation from proto-neutron star (PNS) convection below the neutrinosphere appears superimposed on the low-frequency trend associated with the aspherical expansion of the SN shock after the onset of the explosion. Relativistic effects in combination with detailed neutrino transport are shown to be essential for quantitative predictions of the GW frequency evolution and energy spectrum, because they determine the structure of the PNS surface layer and its characteristic g-mode frequency. Burst-like high-frequency activity phases, correlated with sudden luminosity increase and spectral hardening of electron (anti-)neutrino emission for some 10ms, are discovered as new features after the onset of the explosion. They correspond to intermittent episodes of anisotropic accretion by the PNS in the case of fallback SNe. We find stronger signals for more massive progenitors with large accretion rates. The typical frequencies are higher for massive PNSs, though the time-integrated spectrum also strongly depends on the model dynamics.