The remnant of neutron star mergers is dense in neutrinos. By employing inputs from one hydrodynamical simulation of a binary neutron star merger remnant with a black hole of 3 M-circle dot in the center, dimensionless spin parameter 0.8 and an accretion torus of 0.3 M-circle dot, the neutrino emission properties are investigated as the merger remnant evolves. Initially, the local number density of (nu) over bar (e) is larger than that of nu(e) everywhere above the remnant. Then, as the torus approaches self-regulated equilibrium, the local abundance of neutrinos overcomes that of antineutrinos in a funnel around the polar region. The region where the fast pairwise flavor conversions can occur shrinks accordingly as time evolves. Still, we find that fast flavor conversions do affect most of the neutrino-driven ejecta. Assuming that fast flavor conversions lead to flavor equilibration, a significant enhancement of nuclei with mass numbers A > 130 is found as well as a change of the lanthanide mass fraction by more than a factor of a thousand. Our findings hint towards a potentially relevant role of neutrino flavor oscillations for the prediction of the kilonova (macronova) light curves and motivate further work in this direction.
Considering six general relativistic, two-dimensional (2D) supernova (SN) explosion models of progenitor stars between 8.1 and 27 M-circle dot, we systematically analyze the properties of the neutrino emission from core collapse and bounce to the post-explosion phase. The models were computed with the Vertex-CoCoNuT code, using three-flavor, energy-dependent neutrino transport in the ray-by-ray-plus approximation. Our results confirm the close similarity of the mean energies, < E >, of (nu) over bar (e) and heavy-lepton neutrinos and even their crossing during the accretion phase for stars with M >= 10 M-circle dot as observed in previous 1D and 2D simulations with state-of-the-art neutrino transport. We establish a roughly linear scaling of < E-(nu) over bare > with the proto-neutron star (PNS) mass, which holds in time as well as for different progenitors. Convection inside the PNS affects the neutrino emission on the 10%-20% level, and accretion continuing beyond the onset of the explosion prevents the abrupt drop of the neutrino luminosities seen in artificially exploded 1D models. We demonstrate that a wavelet-based time-frequency analysis of SN neutrino signals in IceCube will offer sensitive diagnostics for the SN core dynamics up to at least similar to 10 kpc distance. Strong, narrow-band signal modulations indicate quasi-periodic shock sloshing motions due to the standing accretion shock instability (SASI), and the frequency evolution of such "SASI neutrino chirps" reveals shock expansion or contraction. The onset of the explosion is accompanied by a shift of the modulation frequency below 40-50 Hz, and post-explosion, episodic accretion downflows will be signaled by activity intervals stretching over an extended frequency range in the wavelet spectrogram.
The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.
We perform hydrodynamic supernova (SN) simulations in spherical symmetry for over 100 single stars of solar metallicity to explore the progenitor-explosion and progenitor-remnant connections established by the neutrino-driven mechanism. We use an approximative treatment of neutrino transport and replace the high-density interior of the neutron star (NS) by an inner boundary condition based on an analytic proto-NS core-cooling model, whose free parameters are chosen such that explosion energy, nickel production, and energy release by the compact remnant of progenitors around 20 M-circle dot are compatible with SN 1987A. Thus, we are able to simulate the accretion phase, initiation of the explosion, subsequent neutrino-driven wind phase for 15-20 s, and the further evolution of the blast wave for hours to days until fallback is completed. Our results challenge long-standing paradigms. We find that remnant mass, launch time, and properties of the explosion depend strongly on the stellar structure and exhibit large variability even in narrow intervals of the progenitors' zero-age main-sequence mass. While all progenitors with masses below similar to 15 M-circle dot yield NSs, black hole (BH) as well as NS formation is possible for more massive stars, where partial loss of the hydrogen envelope leads to weak reverse shocks and weak fallback. Our NS baryonic masses of similar to 1.2-2.0 M-circle dot and BH masses >6 M-circle dot are compatible with a possible lack of low-mass BHs in the empirical distribution. Neutrino heating accounts for SN energies between some 10(50) erg and similar to 2 x 10(51) erg but can hardly explain more energetic explosions and nickel masses higher than 0.1-0.2 M-circle dot. These seem to require an alternative SN mechanism.
Following a simulation approach of recent publications, we explore the viability of the neutrino-heating explosion mechanism's dependence on the spatial dimension. Our results disagree with previous findings. While we also observe that two-dimensional (2D) models explode for lower driving neutrino luminosity than spherically symmetric (1D) models, we do not find that explosions in 3D occur easier and earlier than in 2D. Moreover, we find that the average entropy of matter in the gain layer hardly depends on the dimension and thus is not a good diagnostic quantity for the readiness to explode. Instead, mass, integrated entropy, total neutrino-heating rate, and non-radial kinetic energy in the gain layer are higher when models are closer to explosion. Coherent, large-scale mass motions as typically associated with the standing accretion-shock instability (SASI) are observed to be supportive for explosions because they drive strong shock expansion and thus enlarge the gain layer. While 2D models with better angular resolution clearly explode more easily, the opposite trend is seen in 3D. We interpret this as a consequence of the turbulent energy cascade, which transports energy from small to large spatial scales in 2D, thus fostering SASI activity. In contrast, the energy flow in 3D is in the opposite direction, feeding fragmentation and vortex motions on smaller scales and thus making the 3D evolution with finer grid resolution more similar to 1D. More favorable conditions for explosions in 3D may therefore be tightly linked to efficient growth of low-order SASI modes including nonaxisymmetric ones.
Self-induced flavor conversions of supernova (SN) neutrinos can strongly modify the flavor-dependent fluxes. We perform a linearized flavor stability analysis with accretion-phase matter profiles of a 15M[symbol: see text] spherically symmetric model and corresponding neutrino fluxes. We use realistic energy and angle distributions, the latter deviating strongly from quasi-isotropic emission, thus accounting for both multiangle and multienergy effects. For our matter and neutrino density profile we always find stable conditions: flavor conversions are limited to the usual Mikheyev-Smirnov-Wolfenstein effect. In this case one may distinguish the neutrino mass hierarchy in a SN neutrino signal if the mixing angle θ13 is as large as suggested by recent experiments.
Supernova theory, numerical and analytic, has made remarkable progress in the past decade. This progress was made possible by more sophisticated simulation tools, especially for neutrino transport, improved microphysics, and deeper insights into the role of hydrodynamic instabilities. Violent, large-scale nonradial mass motions are generic in supernova cores. The neutrino-heating mechanism, aided by nonradial flows, drives explosions, albeit low-energy ones, of ONeMg-core and some Fe-core progenitors. The characteristics of the neutrino emission from new-born neutron stars were revised, new features of the gravitational-wave signals were discovered, our notion of supernova nucleosynthesis was shattered, and our understanding of pulsar kicks and explosion asymmetries was significantly improved. But simulations also suggest that neutrino-powered explosions might not explain the most energetic supernovae and hypernovae, which seem to demand magnetorotational driving. Now that modeling is being advanced from two to three dimensions, more realism, new perspectives, and hopefully answers to long-standing questions are coming into reach.
We examine r-process nucleosynthesis in the neutrino-driven wind from the thick accretion disk (or “torus”) around a black hole. Such systems are expected as remnants of binary neutron star or neutron star–black hole mergers. We consider a simplified, analytic, time-dependent evolution model of a 3 M☉ central black hole surrounded by a neutrino emitting accretion torus with 90 km radius, which serves as basis for computing spherically symmetric neutrino-driven wind solutions. We find that ejecta with modest entropies (∼30 per nucleon in units of the Boltzmann constant) and moderate expansion timescales (∼100 ms) dominate in the mass outflow. The mass-integrated nucleosynthetic abundances are in good agreement with the solar system r-process abundance distribution if a minimal value of the electron fraction at the charged-particle freezeout, Ye, min ∼ 0.2, is achieved. In the case of Ye, min ∼ 0.3, the production of r-elements beyond A ∼ 130 does not reach to the third peak but could still be important for an explanation of the abundance signatures in r-process deficient stars in the early Galaxy. The total mass of the ejected r-process nuclei is estimated to be ∼1 × 10−3 M☉. If our model was representative, this demands a Galactic event rate of ∼2 × 10−4 yr−1 for black-hole-torus winds from merger remnants to be the dominant source of the r-process elements. Our result thus suggests that black-hole-torus winds from compact binary mergers have the potential to be a major, but probably not the dominant, production site of r-process elements.
Nucleosynthesis by alpha particles and heavier 4n nuclei are of great interest as they would involve nuclear cluster resonances. The role of nuclear clustering is discussed for nucleosynthesis with the Cluster Nucleosynthesis Diagram (CND) proposed before, especially those involving alpha induced reactions, based on our recent works of (alpha,p) reactions with low energy RI beams. We present experimental results that alpha resonances play a crucial role for the (alpha,p) reaction cross sections. Molecular resonances are also briefly discussed along this line for O- and C-burning.
Although the rapid neutron-capture process, or r-process, is fundamentally important for explaining the origin of approximately half of the stable nuclei with A > 60, the astrophysical site of this process has not been identified yet. Here we study r-process nucleosynthesis in material that is dynamically ejected by tidal and pressure forces during the merging of binary neutron stars (NSs) and within milliseconds afterwards. For the first time we make use of relativistic hydrodynamical simulations of such events, defining consistently the conditions that determine the nucleosynthesis, i.e., neutron enrichment, entropy, early density evolution and thus expansion timescale, and ejecta mass. We find that 10^{-3}-10^{-2} solar masses are ejected, which is enough for mergers to be the main source of heavy (A > 140) galactic r-nuclei for merger rates of some 10^{-5} per year. While asymmetric mergers eject 2-3 times more mass than symmetric ones, the exact amount depends weakly on whether the NSs have radii of ~15 km for a "stiff" nuclear equation of state (EOS) or ~12 km for a "soft" EOS. R-process nucleosynthesis during the decompression becomes largely insensitive to the detailed conditions because of efficient fission recycling, producing a composition that closely follows the solar r-abundance distribution for nuclei with mass numbers A > 140. Estimating the light curve powered by the radioactive decay heating of r-process nuclei with an approximative model, we expect high emission in the B-V-R bands for 1-2 days with potentially observable longer duration in the case of asymmetric mergers because of the larger ejecta mass.
This manuscript reviews recent progress in our understanding of the nucleosynthesis of medium and heavy elements in supernovae.Recent hydrodynamical models of core-collapse supernovae show that a large amount of proton rich matter is ejected under strong neutrino fluxes.This matter constitutes the site of the ν p-process where antineutrino absorption reactions catalyze the nucleosynthesis of nuclei with A > 64.Supernovae are also associated with the r-process responsible for the synthesis of the heaviest elements in nature.Fission during the r-process can play a major role in determining the final abundance patter and in explaining the almost universal features seen in metal-poor r-process-rich stars.
Compact object mergers have long been speculated to be a possible site of r-process nucleosynthesis.While most attention has been focused on the cold decompression of neutron star matter ejected from the merger, other sites within the merger likely contribute to its nucleosynthetic output.Here we consider hot outflows from the accretion disk that forms around the black hole following a black hole -neutron star merger.We begin with the results of a three-dimensional numerical merger model and carefully calculate the neutrino and antineutrino fluxes emitted from the accretion disk.We find that neutrino interactions on free nucleons in the outflowing material result in neutron excesses such that at least a weak r-process is produced and in some cases a main r-process as well.Additionally, we find that the weak r-process pattern calculated for certain trajectories compares favorably to the pattern observed in a weak r-process-enhanced halo star.
Based on the shell model for Gamow-Teller and the random phase approximation for forbidden transitions, we calculate cross sections for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by nuclear statistical equilibrium. The cross sections are incorporated into state-of-the-art stellar core-collapse simulations with detailed energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%.
Based on recent models of relativistic jet formation by thermal energy deposition around black hole-torus systems, the relation between the onand off-axis appearance of short, hard gamma-ray bursts (GRBs) is discussed in terms of energetics, duration, average Lorentz factor, and probability of observation, assuming that the central engines are remnants of binary neutron star (NS+NS) or neutron star-black hole (NS+BH) mergers. As a consequence of the interaction with the torus matter at the jet basis and the subsequent expansion of the jets into an extremely low-density environment, the collimated ultrarelativistic outflows possess flat core profiles with only little variation of radially-averaged properties, and are bounded by very steep lateral edges. Owing to the rapid decrease of the isotropic-equivalent energy near the jet edges, the probability of observing the lateral, lower Lorentz factor wings is significantly reduced and most short GRBs should be seen with on-axis-like properties. Taking into account cosmological and viewing angle effects, theoretical predictions are made for the short-GRB distributions with redshift z, fluence, and isotropic-equivalent energy. The observational data for short bursts with determined redshifts are found to be compatible with the predictions only if either the intrinsic GRB rate density drops rapidly at z & 1, or a large number of events at z > 1 are missed, implying that the subenergetic GRB 050509b was an extremely rare low-fluence event with detectable photon flux only because of its proximity and shortness. It appears Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1, D-85741 Garching, Germany Departamento de Astronomı́a y Astrof́ısica, Universidad de Valencia, 46100 Burjassot, Spain Istituto Naz. di Astrofisica OATs, Via Tiepolo, 11, I-34131 Trieste, Italy