Proto-neutron Star Research Articles

Ab-initio General-relativistic Neutrino-radiation Hydrodynamics Simulations of Long-lived Neutron Star Merger Remnants to Neutrino Cooling Timescales

We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼1014.5 g cm−3, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T/∣W∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Constraining the long-lived supramassive neutron stars by magnetar boosted kilonovae

ABSTRACT Kilonovae are optical transients following the merger of neutron star binaries, which are powered by the r-process heating of merger ejecta. However, if a merger remnant is a long-lived supramassive neutron star supported by its uniform rotation, it will inject energy into the ejecta through spin-down power. The energy injection can boost the peak luminosity of a kilonova by many orders of magnitudes, thus significantly increasing the detectable volume. Therefore, even if such events are only a small fraction of the kilonova population, they could dominate the detection rates. However, after many years of optical sky surveys, no such event has been confirmed. In this work, we build a boosted kilonova model with rich physical details, including the description of the evolution and stability of a proto neutron star, and the energy absorption through X-ray photoionization. We simulate the observation prospects and find the only way to match the absence of detection is to limit the energy injection by the newly born magnetar to only a small fraction of the neutron star rotational energy, thus they should collapse soon after the merger. Our result indicates that most supramassive neutron stars resulting from binary neutron star mergers are short lived and they are likely to be rare in the Universe.

Anomaly induced cooling of neutron stars: a Standard Model contribution

Young neutron stars cool via the emission of neutrinos from their core. A precise understanding of all the different processes producing neutrinos in the hot and degenerate matter is essential for assessing the cooling rate of such stars. The main Standard Model processes contributing to this effect are ν bremsstrahlung, mURCA among others. In this paper, we investigate another Standard Model process initiated by the Wess-Zumino-Witten term, leading to the emission of neutrino pairs via Nγ → Nνν̅. We find that for proto-neutron stars, such processes with degenerate neutrons can be comparable and even dominate over the typical and well-known cooling mechanisms.

The early evolution of magnetar rotation – II. Rapidly rotating magnetars: implications for gamma-ray bursts and superluminous supernovae

ABSTRACT Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and superluminous supernovae (SLSNe). Using a suite of two-dimensional magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytical estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength B0, PNS angular velocity Ω⋆, PNS radius R⋆, and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods P⋆ ≲ 4 ms and polar magnetic field strength B0 ≳ 1015 G can release 1050 to 5 × 1051 erg of energy during the first ∼2 s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of B0 ≲ 5 × 1014 G do not release a large fraction of their rotational kinetic energy during the cooling phase and, hence, are not likely to power GRBs. Although we cannot simulate to times greater than ∼3–5 s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on time-scales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase.

Failed supernova simulations beyond black hole formation

ABSTRACT We present an axisymmetric failed supernova simulation beyond black hole formation, for the first time with numerical relativity and two-moment multi-energy neutrino transport. To ensure stable numerical evolution, we use an excision method for neutrino radiation hydrodynamics within the inner part of black hole domain. We demonstrate that our excision method is capable of stably evolving the radiation hydrodynamics in dynamical black hole space–time. As a remarkable signature of the final moment of proto-neutron star (PNS), we find the emergence of high-energy neutrinos. Those high-energy neutrinos are associated with the PNS shock surface being swallowed by the central black hole and could be a possible observable of failed supernovae.

Three-dimensional General-relativistic Simulations of Neutrino-driven Winds from Magnetized Proto–Neutron Stars

Formed in the aftermath of a core-collapse supernova or neutron star merger, a hot proto–neutron star (PNS) launches an outflow driven by neutrino heating lasting for up to tens of seconds. Though such winds are considered potential sites for the nucleosynthesis of heavy elements via the rapid neutron capture process (r-process), previous work has shown that unmagnetized PNS winds fail to achieve the necessary combination of high entropy and/or short dynamical timescale in the seed nucleus formation region. We present three-dimensional general-relativistic magnetohydrodynamical simulations of PNS winds which include the effects of a dynamically strong (B ≳ 1015 G) dipole magnetic field. After initializing the magnetic field, the wind quickly develops a helmet-streamer configuration, characterized by outflows along open polar magnetic field lines and a “closed” zone of trapped plasma at lower latitudes. Neutrino heating within the closed zone causes the thermal pressure of the trapped material to rise in time compared to the polar outflow regions, ultimately leading to the expulsion of this matter from the closed zone on a timescale of ∼60 ms, consistent with the predictions of Thompson. The high entropies of these transient ejecta are still growing at the end of our simulations and are sufficient to enable a successful second-peak r-process in at least a modest ≳1% of the equatorial wind ejecta.

The proto-neutron star inner crust in a multi-component plasma approach

Context. Proto-neutron stars are born hot, with temperatures exceeding a few times 1010 K. In these conditions, the crust of the proto-neutron star is expected to be made of a Coulomb liquid and composed of an ensemble of different nuclear species. Aims. In this work, we perform a study of the beta-equilibrated proto-neutron-star crust in the liquid phase in a self-consistent multi-component approach. This also allows us to perform a consistent calculation of the impurity parameter, which is often taken as a free parameter in cooling simulations. Methods. To this aim, we developed a self-consistent multi-component approach at finite temperature using a compressible liquid-drop description of the ions, with surface parameters adjusted to reproduce experimental masses. The treatment of the ion centre-of-mass motion was included through a translational free-energy term accounting for in-medium effects. The results of the self-consistent calculations of the multi-component plasma are systematically compared with those performed in a perturbative treatment as well as in the one-component plasma approximation. Results. We show that the inclusion of non-linear mixing terms arising from the ion centre-of-mass motion leads to a breakdown of the ensemble equivalence between the one-component and multi-component approach. Our findings also illustrate that the abundance of light nuclei becomes important and eventually dominates the whole distribution at higher density and temperature in the crust. This is reflected in the impurity parameter, which, in turn, may have a potential impact on neutron-star cooling. For practical application to astrophysical simulations, we also provide a fitting formula for the impurity parameter in the proto-neutron-star inner crust. Conclusions. Our results obtained within a self-consistent multi-component approach show important differences in the prediction of the proto-neutron-star composition with respect to those obtained with a one-component approximation or a perturbative multi-component approximation, particularly in the deeper region of the crust. This highlights the importance of a full, self-consistent multi-component plasma calculation for reliable predictions of the proto-neutron-star crust composition.

Gravitational Wave Eigenfrequencies from Neutrino-driven Core-collapse Supernovae

Core-collapse supernovae (CCSNe) are predicted to produce gravitational waves (GWs) that may be detectable by Advanced LIGO/Virgo. These GW signals carry information from the heart of these cataclysmic events, where matter reaches nuclear densities. Recent studies have shown that it may be possible to infer the properties of the proto-neutron star (PNS) via GWs generated by hydrodynamic perturbations of the PNS. However, we lack a comprehensive understanding of how these relationships may change with the properties of CCSNe. In this work, we build a self-consistent suite of over 1000 exploding CCSNe from a grid of progenitor masses and metallicities combined with six different nuclear equations of state (EOS). Performing a linear perturbation analysis on each model, we compute the resonant GW frequencies of the PNS, and we motivate a time-agnostic method for identifying characteristic frequencies of the dominant GW emission. From this, we identify two characteristic frequencies, of the early- and late-time signal, that measure the surface gravity of the cold remnant neutron star, and simultaneously constrain the hot nuclear EOS. However, we find that the details of the CCSN model, such as the treatment of gravity or the neutrino transport, and whether it explodes, noticeably change the magnitude and evolution of the PNS eigenfrequencies.

Energy and angle dependence of neutrino scattering rates in proto–neutron star and supernova matter within Skyrme RPA

Energy and angle dependence of neutrino scattering rates in proto–neutron star and supernova matter within Skyrme RPA

Observing Supernova Neutrino Light Curves with Super-Kamiokande. IV. Development of SPECIAL BLEND: A New Public Analysis Code for Supernova Neutrinos

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.

Numerical simulations of the Tayler–Spruit dynamo in proto-magnetars

ABSTRACT The Tayler–Spruit dynamo is one of the most promising mechanisms proposed to explain angular momentum transport during stellar evolution. Its development in proto-neutron stars spun-up by supernova fallback has also been put forward as a scenario to explain the formation of very magnetized neutron stars called magnetars. Using three-dimensional direct numerical simulations, we model the proto-neutron star interior as a stably stratified spherical Couette flow with the outer sphere that rotates faster than the inner one. We report the existence of two subcritical dynamo branches driven by the Tayler instability. They differ by their equatorial symmetry (dipolar or hemispherical) and the magnetic field scaling, which is in agreement with different theoretical predictions (by Fuller and Spruit, respectively). The magnetic dipole of the dipolar branch is found to reach intensities compatible with observational constraints on magnetars.

Numerical simulations of the random angular momentum in convection – II. Delayed explosions of red supergiants following ‘failed’ supernovae

ABSTRACT When collapse of the iron core in a massive red or yellow supergiant does not lead to an energetic supernova, a significant fraction of the convective hydrogen envelope will fall in towards the black hole formed from the collapsing core. The random velocity field in the convective envelope results in finite specific angular momentum in each infalling shell. Using 3D hydrodynamical simulations, we follow the infall of this material to small radii, resolving the circularization radii of the flow. We show that infall of the convective envelope leads to nearly complete envelope ejection in a ≳1048 erg explosion with outflow speeds of ≳200 km s−1. The light curve of such an explosion would show a characteristic, red plateau as the ejecta cools and a hydrogen recombination front recedes through the expanding ejecta. Adopting supernova IIp scalings, the event would have a plateau luminosity of ≳1040 erg s−1 and a duration of several hundreds of days. These events would appear quite similar to luminous red novae with red or yellow supergiant progenitors; some luminous red novae may, in fact, be signposts of black hole formation. The mechanism studied here produces more energetic explosions than the weak shock generated from radiation of neutrino energy during the protoneutron star phase. Because we cannot simulate all the way to the horizon, our results are likely lower limits on the energy and luminosity of transients produced during the collapse of a red or yellow supergiant to form a black hole.

3D evolution of neutron star magnetic fields from a realistic core-collapse turbulent topology

ABSTRACT We perform the first 3D fully coupled magneto-thermal simulations of neutron stars (including the most realistic background structure and microphysical ingredients so far) applied to a very complex initial magnetic field topology in the crust, similar to what was recently obtained by proto-neutron stars dynamo simulations. In such configurations, most of the energy is stored in the toroidal field, while the dipolar component is a few per cent of the mean magnetic field. This initial feature is maintained during the long-term evolution (∼106 yr), since the Hall term favours a direct cascade (compensating for Ohmic dissipation) rather than a strong inverse cascade, for such an initial field topology. The surface dipolar component, responsible for the dominant electromagnetic spin-down torque, does not show any increase in time, when starting from this complex initial topology. This is in contrast to the timing properties of young pulsars and magnetars which point to higher values of the surface dipolar fields. A possibility is that the deep-seated magnetic field (currents in the core) is able to self-organize in large scales (during the collapse or in the early life of a neutron star). Alternatively, the dipolar field might be lower than is usually thought, with magnetosphere substantially contributing to the observed high spin-down, via e.g. strong winds or strong coronal magnetic loops, which can also provide a natural explanation to the tiny surface hotspots inferred from X-ray data.

Measuring the properties of f−mode oscillations of a protoneutron star by third-generation gravitational-wave detectors

Core-collapse supernovae are among the astrophysical sources of gravitational waves that could be detected by third-generation gravitational-wave detectors. Here, we analyze the gravitational-wave strain signals from two- and three-dimensional simulations of core-collapse supernovae generated using the code fornax. A subset of the two-dimensional simulations has nonzero core rotation at the core bounce. A dominant source of time changing quadrupole moment is the $l=2$ fundamental mode ($f\text{\ensuremath{-}}\mathrm{mode}$) oscillation of the protoneutron star. From the time-frequency spectrogram of the gravitational-wave strain we see that, starting $\ensuremath{\sim}400\text{ }\text{ }\mathrm{ms}$ after the core bounce, most of the power lies within a narrow track that represents the frequency evolution of the $f\text{\ensuremath{-}}\mathrm{mode}$ oscillations. The $f\text{\ensuremath{-}}\mathrm{mode}$ frequencies obtained from linear perturbation analysis of the angle-averaged profile of the protoneutron star corroborate what we observe in the spectrograms of the gravitational-wave signal. We explore the measurability of the $f\text{\ensuremath{-}}\mathrm{mode}$ frequency evolution of a protoneutron star for a supernova signal observed in the third-generation gravitational-wave detectors. Measurement of the frequency evolution can reveal information about the masses, radii, and densities of the protoneutron stars. We find that if the third-generation detectors observe a supernova within 10 kpc, then we can measure these frequencies to within 5 Hz rms error. We can also measure the energy emitted in the fundamental $f\text{\ensuremath{-}}\mathrm{mode}$ using the spectrogram data of the strain signal. We find that the energy in the $f\text{\ensuremath{-}}\mathrm{mode}$ can be measured to within 20% error for signals observed by Cosmic Explorer using simulations with successful explosion, assuming source distances within 10 kpc.

Effect of stellar rotation on the development of post-shock instabilities during core-collapse supernovae

Context. The growth of hydrodynamical instabilities is key to triggering a core-collapse supernova explosion during the phase of stalled accretion shock, immediately after the birth of a proto-neutron star (PNS). Stellar rotation is known to affect the standing accretion shock instability (SASI) even for small rotation rates, but its effect on the onset of neutrino-driven convection is still poorly known. Aims. We assess the effect of stellar rotation on SASI when neutrino heating is taken into account as well as the effect of rotation on neutrino-driven convection. The interplay of rotation with these two instabilities affects the frequency of the mode m = 2, which can be detected with gravitational waves at the onset of a supernova explosion. Methods. We used a linear stability analysis to study the dynamics of the accreting gas in the equatorial plane between the surface of the PNS and the stationary shock. We explored rotation effects on the relative strength of SASI and convection by considering a large range of specific angular momenta and neutrino luminosities. Results. The nature of the dominant non-axisymmetric instability developing in the equatorial post-shock region depends on both the convection parameter, χ, and the rotation rate. Equatorial convective modes with χ ≳ 5 are hampered by differential rotation. At smaller χ, however, mixed SASI-convective modes with a large angular scale, m = 1, 2, 3, can take advantage of rotation and become dominant for relatively low rotation rates, at which centrifugal effects are small. For rotation rates exceeding ∼30% of the Keplerian rotation at the PNS surface, a new instability regime is characterised by a frequency that, when measured in units of the post-shock velocity and radius, vsh/rsh, is nearly independent of the convection parameter, χ. A strong prograde m = 2 spiral dominates over a large parameter range and is favorable to the production of gravitational waves. In this regime, a simple linear relation exists between the oscillation frequency of the dominant mode and the specific angular momentum of the accreted gas. Conclusions. Three different regimes of post-shock instabilities can be distinguished depending on the rotation rate. For low rotation rates (less than 10% of the Keplerian rotation at the PNS surface), differential rotation has a linear destabilising effect on SASI and a quadratic stabilising or destabilising effect on the purely convective equatorial modes depending on their azimuthal wavenumber. Intermediate rotation rates (10% to 30% of the Keplerian rotation) lead to the emergence of mixed SASI-convection-rotation modes that involve large angular scales. Finally, strong rotation erases the influence of the buoyancy and heating rate on the instability. This independence allows for a reduction in the parameter space, which can be helpful for gravitational wave analysis.

Fast neutrino flavor conversion in core-collapse supernovae: A parametric study in 1D models

We explore the impact of small-scale flavor conversions of neutrinos, the so-called fast flavor conversions (FFCs), on the dynamical evolution and neutrino emission of core-collapse supernovae (CCSNe). In order to do that, we implement FFCs in the spherically symmetric (1D) CCSN simulations of a 20 solar-mass progenitor model parametrically, assuming that FFCs happen at densities lower than a systematically varied threshold value and lead to an immediate flavor equilibrium consistent with lepton number conservation. We find that besides hardening the electron neutrino and antineutrino spectra, which helps the expansion of the shock by enhanced postshock heating, FFCs can cause significant, nontrivial modifications of the energy transport in the SN environment via increasing the heavy-lepton neutrino luminosities. In our non-exploding models this results in extra cooling of the layers around the neutrinospheres, which triggers a faster contraction of the proto-neutron star and hence, in our 1D models, hampers the CCSN explosion. Although our study is limited by the 1D nature of our simulations, it provides valuable insights into how neutrino flavor conversions in the deepest CCSN regions can impact the neutrino release and the corresponding response of the stellar medium.

Analytic approach to axion-like-particle emission in core-collapse supernovae

We investigate the impact of a presumed axion-like-particle (ALP) emission in a core-collapse supernova explosion on neutrino luminosities and mean energies employing a relatively simple analytic description. We compute the nuclear Bremsstrahlung and Primakoff axion luminosities as functions of the protoneutron star (PNS) parameters and discuss how the ALP luminosities compete with the neutrino emission, modifying the total PNS thermal energy dissipation. Our results are publicly available in the python package ARtiSANS, which can be used to compute the neutrino and axion observables for different choices of parameters.

Three approaches for the classification of protoneutron star oscillation modes

ABSTRACT The future detection of gravitational waves (GWs) from a Galactic core-collapse supernova will provide information on the physics inside protoneutron stars (PNS). In this work, we apply three different classification methods for the PNS non-radial oscillation modes: Cowling classification, Generalized Cowling Nomenclature (GCN), and a classification based on modal properties (CBMP). Using PNS models from 3D simulations of core-collapse supernovae, we find that in the early stages of the PNS evolution, typically 0.4 s after the bounce, the Cowling classification is inconsistent, but the GCN and the CBMP provide complementary information that helps to understand the evolution of the modes. In the GCN, we note several avoided crossings as the mode frequencies evolve at early times, while the CBMP tracks the modes across the avoided crossings. We verify that the strongest emission of GWs by the PNS corresponds to the f mode in the GCN, indicating that the mode trapping region alternates between the core and the envelope at each avoided crossing. At later times, approximately 0.4 s after the bounce, the three classification methods present a similar description of the mode spectrum. We use our results to test universal relations for the PNS modes according to their classification and find that the behaviour of the universal relations for f and p modes is remarkably simple in the CBMP.

Spontaneous scalarization as a new core-collapse supernova mechanism and its multimessenger signals

We perform multi-dimensional core-collapse supernova (CCSN) simulations in a massive scalar-tensor theory for the first time with a realistic equation of state and multi-energy neutrino radiation. Among the set of our models varying the scalar mass and the coupling strength between the scalar and gravitational fields, a particular model allows for recurrent spontaneous scalarizations (SSs) in the proto-neutron star (PNS). Each SS induces the PNS collapse and subsequent bounce, from which devastating shock waves emanate and eject the PNS envelope. The explosion energy can easily exceed $\mathcal O(10^{51})$ erg. This study reveals new aspects of SS as the explosion mechanism of CCSNe. We also discuss its characteristic multi-messenger signals: neutrinos and gravitational waves.