Source author record

H. -Thomas Janka

H. -Thomas Janka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2026arXiv

Gravitational-Wave Signatures of Nonstandard Neutrino Properties in Collapsing Stellar Cores

We present a novel multi-messenger approach for probing nonstandard neutrino properties through the detection of gravitational waves (GWs) from collapsing stellar cores and associated supernova explosions. We show that neutrino flavor conversion inside the proto-neutron star (PNS), motivated by physics Beyond the Standard Model (BSM), can significantly boost PNS convection. This effect leads to large-amplitude GW emission over a wide frequency range during an otherwise relatively quiescent GW phase shortly after core bounce. Such a signal provides a promising new avenue for exploring nonstandard neutrino phenomena and other BSM physics impacting PNS convection.

preprint2021arXiv

Supernova Fallback as Origin of Neutron Star Spins and Spin-kick Alignment

Natal kicks and spins are characteristic properties of neutron stars (NSs) and black holes (BHs). Both offer valuable clues to dynamical processes during stellar core collapse and explosion. Moreover, they influence the evolution of stellar multiple systems and the gravitational-wave signals from their inspiral and merger. Observational evidence of possibly generic spin-kick alignment has been interpreted as indication that NS spins are either induced with the NS kicks or inherited from progenitor rotation, which thus might play a dynamically important role during stellar collapse. Current three-dimensional supernova simulations suggest that NS kicks are transferred in the first seconds of the explosion, mainly by anisotropic mass ejection and, on a secondary level, anisotropic neutrino emission. In contrast, the NS spins are only determined minutes to hours later by angular momentum associated with fallback of matter that does not become gravitationally unbound in the supernova. Here, we propose a novel scenario to explain spin-kick alignment as a consequence of tangential vortex flows in the fallback matter that is accreted mostly from the direction of the NS's motion. For this effect the initial NS kick is crucial, because it produces a growing offset of the NS away from the explosion center, thus promoting onesided accretion. In this new scenario conclusions based on traditional concepts are reversed. For example, pre-kick NS spins are not required, and rapid progenitor-core rotation can hamper spin-kick alignment. We also discuss implications for natal BH kicks and the possibility of tossing the BH's spin axis during its formation.

preprint2020arXiv

Fast Neutrino Flavor Instability in the Neutron-star Convection Layer of Three-dimensional Supernova Models

Neutrinos from a supernova (SN) might undergo fast flavor conversions near the collapsed stellar core. We perform a detailed study of this intriguing possibility, analyzing time-dependent state-of-the-art 3D SN models of 9 and 20 Msun. Both models were computed with multi-D three-flavor neutrino transport based on a two-moment solver, and both exhibit the presence of the lepton-number emission self-sustained asymmetry (LESA). The transport solution does not provide the angular distributions of the neutrino fluxes, which are crucial to track the fast flavor instability. To overcome this limitation, we use a recently proposed approach based on the angular moments of the energy-integrated electron lepton-number distribution. With this method we find the possibility of fast neutrino flavor instability at radii <~20 km, which is well interior to the neutrinosphere. Our results confirm recent observations in a 2D SN model and in 2D/3D models with fixed matter background, which were computed with Boltzmann neutrino transport. However, the flavor unstable locations are not isolated points as discussed previously, but thin skins surrounding volumes where electron antineutrinos are more abundant than electron neutrinos. These volumes grow with time and appear first in the convective layer of the proto-neutron star (PNS), where a decreasing electron fraction (Ye) and high temperatures favor the occurrence of regions with negative neutrino chemical potential. Since Ye remains higher in the LESA dipole direction, where convective lepton-number transport out from the nonconvective PNS core slows down the deleptonization, flavor unstable conditions become more widespread in the opposite hemisphere. This interesting phenomenon deserves further investigation, since its impact on SN modeling and possible consequences for SN dynamics and neutrino observations are presently unclear. (abridged)

preprint2020arXiv

How much H and He is "hidden" in SNe Ib/c? -- II. Intermediate-mass objects: a 22 M$_{\odot}$ progenitor case study

Stripped envelope supernovae are a sub-class of core collapse supernovae showing several stages of H/He shell stripping that determines the class: H-free/He-poor SNe are classified as Type Ic, H-poor/He-rich are Type Ib, and H/He-rich are Type IIb. Stripping H/He with only stellar wind requires significantly higher mass loss rates than observed while binary-involved mass transfer may usually not strip enough to produce H/He free SNe. Type Ib/c SNe are sometimes found to include weak H/He transient lines as a product of a trace amount of H/He left over from stripping processes. The extent and mass of the H/He required to produce these lines is not well known. In this work, a 22 M$_{\odot}$ progenitor model is stripped of the H/He shells in five steps prior to collapse and then exploded at four explosion energies. Requiring both optical and NIR He I lines for helium identification does not allow much He mass to be hidden in SE--SNE. Increasing the mass of He above the CO core delays the visibility of O I 7774 in early spectra. Our SN Ib-like models are capable of reproducing the spectral evolution of a set of observed SNe with reasonable estimated $E_\mathrm{k}$ accuracy. Our SN\,IIb-like models can partially reproduce low energy observed SN IIb, but we find no observed comparison for the SN IIb-like models with high $E_\mathrm{k}$.

preprint2020arXiv

Resolution Study for Three-dimensional Supernova Simulations with the Prometheus-Vertex Code

We present a carefully designed, systematic study of the angular resolution dependence of simulations with the Prometheus-Vertex neutrino-hydrodynamics code. Employing a simplified neutrino heating-cooling scheme in the Prometheus hydrodynamics module allows us to sample the angular resolution between 4 degrees and 0.5 degrees. With a newly-implemented static mesh refinement (SMR) technique on the Yin-Yang grid, the angular coordinates can be refined in concentric shells, compensating for the diverging structure of the spherical grid. In contrast to previous studies with Prometheus and other codes, we find that higher angular resolution and therefore lower numerical viscosity provides more favorable explosion conditions and faster shock expansion. We discuss the possible reasons for the discrepant results. The overall dynamics seem to converge at a resolution of about 1 degree. Applying the SMR setup to marginally exploding progenitors is disadvantageous for the shock expansion, however, because kinetic energy of downflows is dissipated to internal energy at resolution interfaces, leading to a loss of turbulent pressure support and a steeper temperature gradient. We also present a way to estimate the numerical viscosity on grounds of the measured turbulent kinetic-energy spectrum, leading to smaller values that are better compatible with the flow behavior witnessed in our simulations than results following calculations in previous literature. Interestingly, the numerical Reynolds numbers in the turbulent, neutrino-heated postshock layer (some 10 to several 100) are in the ballpark of expected neutrino-drag effects on the relevant length scales in the turbulent postshock layer. We provide a formal derivation and quantitative assessment of the neutrino drag terms in an appendix.

preprint2019arXiv

J0453+1559: a neutron star-white dwarf binary from a thermonuclear electron-capture supernova?

The compact binary radio pulsar system J0453+1559 (Martinez et al. 2015) consists of a recycled pulsar as primary component of $1.559(5)\,M_\odot$ and an unseen companion star of $1.174(4)\,M_\odot$. Because of the relatively large orbital eccentricity of $e = 0.1125$, it was argued that the companion is a neutron star, making it the neutron star with the lowest accurately determined mass to date. However, a direct observational determination of the nature of the companion is currently not feasible. Moreover, state-of-the-art stellar evolution and supernova modeling are contradictive concerning the possibility to produce such a low-mass neutron star remnant. Here we challenge the neutron star interpretation by reasoning that the lower-mass component could instead be a white dwarf born in a thermonuclear electron-capture supernova (tECSN) event, in which oxygen-neon deflagration in the degenerate stellar core of an ultra-stripped progenitor ejects several $0.1\,M_\odot$ of matter and leaves a bound ONeFe white dwarf as the second-formed compact remnant. We determine the ejecta mass and remnant kick needed in this scenario to explain the properties of PSR~J0453+1559 by a neutron star--white dwarf system. More work on tECSNe is needed to assess the viability of this scenario.

preprint2019arXiv

The Explosion of Helium Stars Evolved With Mass Loss

Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5 to 40\,\Msun, which correspond to main sequence masses of about 13 to 90\,\Msun. Common Type Ib and Ic supernovae result from stars whose final masses are approximately 2.5 to 5.6\,\Msun. For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10\,\Msun. The median neutron star mass in binaries is 1.35--1.38\,\Msun \ and the median black hole mass is between 9 and 11\,\Msun. Even though black holes less massive than 5 \Msun\ are rare, they are predicted down to the maximum neutron star mass. There is no empty ``gap'', only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common Type Ib and Ic supernovae. Even with a very liberal, but physically plausible increase in $^{56}$Ni production, the highest energy models are fainter, at peak, than 10$^{42.6}$\,erg\,s$^{-1}$, and very few approach that limit. The median peak luminosity ranges from 10$^{42.0}$ to 10$^{42.3}$\,erg\,s$^{-1}$. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of Type I supernovae at both high and low mass are explored.

preprint2016arXiv

Erratum: Progenitor-explosion connection and remnant birth masses for neutrino-driven supernovae of iron-core progenitors (2012, ApJ, 757, 69)

An erroneous interpretation of the hydrodynamical results led to an incorrect determination of the fallback masses in Ugliano et al. (2012), which also (on a smaller level) affects the neutron star masses provided in that paper. This problem was already addressed and corrected in the follow-up works by Ertl et al. (2015) and Sukhbold et al. (2015). Therefore, the reader is advised to use the new data of the latter two publications. In the remaining text of this Erratum we present the differences of the old and new fallback results in detail and explain the origin of the mistake in the original analysis by Ugliano et al. (2012).

preprint2016arXiv

Neutron-star merger ejecta as obstacles to neutrino-powered jets of gamma-ray bursts

We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.

preprint2016arXiv

Physics of Core-Collapse Supernovae in Three Dimensions: a Sneak Preview

Nonspherical mass motions are a generic feature of core-collapse supernovae, and hydrodynamic instabilities play a crucial role for the explosion mechanism. First successful neutrino-driven explosions could be obtained with self-consistent, first-principle simulations in three spatial dimensions (3D). But 3D models tend to be less prone to explosion than corresponding axisymmetric (2D) ones. This has been explained by 3D turbulence leading to energy cascading from large to small spatial scales, inversely to the 2D case, thus disfavoring the growth of buoyant plumes on the largest scales. Unless the inertia to explode simply reflects a lack of sufficient resolution in relevant regions, it suggests that some important aspect may still be missing for robust and sufficiently energetic neutrino-powered explosions. Such deficits could be associated with progenitor properties like rotation, magnetic fields or pre-collapse perturbations, or with microphysics that could lead to an enhancement of neutrino heating behind the shock. 3D simulations have also revealed new phenomena that are not present in 2D, for example spiral modes of the standing accretion shock instability (SASI) and a stunning dipolar lepton-emission self-sustained asymmetry (LESA). Both impose time- and direction-dependent variations on the detectable neutrino signal. The understanding of these effects and of their consequences is still in its infancy.

preprint2015arXiv

A new multidimensional, energy-dependent two-moment transport code for neutrino-hydrodynamics

We present the new code ALCAR developed to model multidimensional, multi energy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution equations of the 0th- and 1st-order angular moments of the specific intensity, supplemented by an algebraic relation for the 2nd-moment tensor to close the system. The scheme takes into account frame-dependent effects of order O(v/c) as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem setups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.

preprint2015arXiv

Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

preprint2015arXiv

Impact of neutrino flavor oscillations on the neutrino-driven wind nucleosynthesis of an electron-capture supernova

Neutrino oscillations, especially to light sterile states, can affect the nucleosynthesis yields because of their possible feedback effect on the electron fraction (Ye). For the first time, we perform nucleosynthesis calculations for neutrino-driven wind trajectories from the neutrino-cooling phase of an 8.8 Msun electron-capture supernova, whose hydrodynamic evolution was computed in spherical symmetry with sophisticated neutrino transport and whose Ye evolution was post-processed by including neutrino oscillations both between active and active-sterile flavors. We also take into account the alpha-effect as well as weak magnetism and recoil corrections in the neutrino absorption and emission processes. We observe effects on the Ye evolution which depend in a subtle way on the relative radial positions of the sterile MSW resonances, of collective flavor transformations, and on the formation of alpha particles. For the adopted supernova progenitor, we find that neutrino oscillations, also to a sterile state with eV-mass, do not significantly affect the element formation and in particular cannot make the post-explosion wind outflow neutron rich enough to activate a strong r-process. Our conclusions become even more robust when, in order to mimic equation-of-state dependent corrections due to nucleon potential effects in the dense-medium neutrino opacities, six cases with reduced Ye in the wind are considered. In these cases, despite the conversion of active neutrinos to sterile neutrinos, Ye increases or is not significantly lowered compared to the values obtained without oscillations and active flavor transformations. This is a consequence of a complicated interplay between sterile-neutrino production, neutrino-neutrino interactions, and alpha-effect.

preprint2015arXiv

Nucleosynthesis in dynamical and torus ejecta of compact binary mergers

We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

preprint2015arXiv

Supernova 1987A: neutrino-driven explosions in three dimensions and light curves

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

preprint2015arXiv

Three-Dimensional Simulations of Core-Collapse Supernovae: From Shock Revival to Shock Breakout

We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.

preprint2014arXiv

Nucleosynthesis constraints on the neutron star-black hole merger rate

We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolved with theoretical predictions of the binary population. This allows us to determine a strict upper limit of the average NS-BH merger rate of ~6*10^-5 per year. We quantify the uncertainties of this estimate to be within factors of a few mostly because of the unknown BH spin distribution of such systems, the uncertain equation of state of NS matter, and possible errors in the Galactic content of r-process material. Our approach implies a correlation between the merger rates of NS-BH binaries and of double NS systems. Predictions of the detection rate of gravitational-wave signals from such compact-object binaries by Advanced LIGO and Advanced Virgo on the optimistic side are incompatible with the constraints set by our analysis.

preprint2013arXiv

Natal Kicks of Stellar-Mass Black Holes by Asymmetric Mass Ejection in Fallback Supernovae

Integrating trajectories of low-mass X-ray binaries containing black holes within the Galactic potential, Repetto, Davies & Sigurdsson recently showed that the large distances of some systems above the Galactic plane can only be explained if black holes receive appreciable natal kicks. Surprisingly, they found that the distribution of black hole kick velocities (rather than that of the momenta) should be similar to that of neutron stars. Here I argue that this result can be understood if neutron star and black hole kicks are a consequence of large-scale asymmetries created in the supernova ejecta by the explosion mechanism. The corresponding anisotropic gravitational attraction of the asymmetrically expelled matter does not only accelerate new-born neutron stars by the "gravitational tug-boat mechanism". It can also lead to delayed black-hole formation by asymmetric fallback of the slowest parts of the initial ejecta onto the transiently existing neutron star, in course of which the momentum of the black hole can grow with the fallback mass. Black hole kick velocities will therefore not be reduced by the ratio of neutron star to black hole mass as would be expected for kicks caused by anisotropic neutrino emission of the nascent neutron star.

preprint2013arXiv

New fission fragment distributions and r-process origin of the rare-earth elements

Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A > 140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties, but calculations sensitively depend on nuclear physics. In particular the fission fragment yields determine the creation of 110 < A < 170 nuclei. Here we apply a new scission-point model, called SPY, to derive the fission fragment distribution (FFD) of all relevant neutron-rich, fissioning nuclei. The model predicts a doubly asymmetric FFD in the abundant A ~ 278 mass region that is responsible for the final recycling of the fissioning material. Using ejecta conditions based on relativistic NS merger calculations we show that this specific FFD leads to a production of the A ~ 165 rare-earth peak that is nicely compatible with the abundance patterns in the Sun and metal-poor stars. This new finding further strengthens the case of NS mergers as possible dominant origin of r-nuclei with A > 140.

preprint2012arXiv

Core-Collapse Supernovae: Reflections and Directions

Core-collapse supernovae are among the most fascinating phenomena in astrophysics and provide a formidable challenge for theoretical investigation. They mark the spectacular end of the lives of massive stars and, in an explosive eruption, release as much energy as the sun produces during its whole life. A better understanding of the astrophysical role of supernovae as birth sites of neutron stars, black holes, and heavy chemical elements, and more reliable predictions of the observable signals from stellar death events are tightly linked to the solution of the long-standing puzzle how collapsing stars achieve to explode. In this article our current knowledge of the processes that contribute to the success of the explosion mechanism are concisely reviewed. After a short overview of the sequence of stages of stellar core-collapse events, the general properties of the progenitor-dependent neutrino emission will be briefly described. Applying sophisticated neutrino transport in axisymmetric (2D) simulations with general relativity as well as in simulations with an approximate treatment of relativistic effects, we could find successful neutrino-driven explosions for a growing set of progenitor stars. First results of three-dimensional (3D) models have been obtained, and magnetohydrodynamic simulations demonstrate that strong initial magnetic fields in the pre-collapse core can foster the onset of neutrino-powered supernova explosions even in nonrotating stars. These results are discussed in the context of the present controversy about the value of 2D simulations for exploring the supernova mechanism in realistic 3D environments, and they are interpreted against the background of the current disagreement on the question whether the standing accretion shock instability (SASI) or neutrino-driven convection is the crucial agency that supports the onset of the explosion.

preprint2012arXiv

Explosion Mechanisms of Core-Collapse Supernovae

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.

preprint2012arXiv

Progenitor-Explosion Connection and Remnant Birth Masses for Neutrino-Driven Supernovae of Iron-Core Progenitors

We perform hydrodynamic supernova 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 solar masses are compatible with Supernova 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 about 15 solar masses 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 ~1.2-2.0 solar masses and BH masses >6 solar masses 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 about 2*10^{51} erg, but can hardly explain more energetic explosions and nickel masses higher than 0.1-0.2 solar masses. These seem to require an alternative SN mechanism.

preprint2011arXiv

R-Process Nucleosynthesis in Dynamically Ejected Matter of Neutron Star Mergers

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.

preprint1995arXiv

The First Second of a Type-II Supernova: Convection, Accretion, and Shock Propagation

One- and two-dimensional hydrodynamical simulations of neutrino-driven supernova explosions are discussed. The simulations cover the phase between the stagnation of the prompt shock and about one second after core bounce. Systematic variation of the neutrino fluxes from the neutrino sphere shows that the explosion energy, explosion time scale, initial mass of the protoneutron star, and explosive nucleosynthesis of iron-group elements depend sensitively on the strength of the neutrino heating during the first few 100 ms after shock formation. Convective overturn in the neutrino-heated region behind the shock is a crucial help for the explosion only in a narrow window of neutrino luminosities. Here powerful explosions can be obtained only in the multi-dimensional case. For higher core-neutrino fluxes also spherically symmetrical models yield energetic explosions, while for lower luminosities even with convection no strong explosions occur.

preprint1995arXiv

WHEN DO SUPERNOVA NEUTRINOS OF DIFFERENT FLAVORS HAVE SIMILAR LUMINOSITIES BUT DIFFERENT SPECTRA ?

Muon and tau neutrinos ($ν_x$) interact with protoneutron star matter only via neutral currents and exchange energy with the stellar gas predominantly by neutrino-electron scattering and neutrino-pair processes. In contrast, electron neutrinos and antineutrinos ($ν_e$ and $\barν_e$) are frequently absorbed and produced in charged-current mediated reactions with nucleons. Therefore the emergent $ν_e$ and $\barν_e$ originate from layers with lower temperatures further out in the star and are emitted with much lower characteristic spectral temperatures. In addition, a major contribution to the $ν_e$ and $\barν_e$ opacities is due to absorptions, while the opacity of $ν_x$ is strongly dominated by scattering reactions with nucleons and nuclei in which the $ν_x$ energy is (essentially) conserved. Therefore the $ν_x$ distribution is nearly isotropic when $ν_x$ decouple energetically and their outward diffusion is slowed down. In a generalized form to include this effect, the Stefan-Boltzmann Radiation Law can account for both the facts that $ν_e$ ($\barν_e$) and $ν_x$ emerge from the star with similar luminosities but with very different characteristic spectral temperatures. Simple analytical expressions to estimate the effect are given. If, as recently argued, even at densities significantly below nuclear matter density neutral-current scatterings were associated with considerable energy transfer between neutrino and target particle, one might expect spectral temperatures of $ν_x$ much closer to those of $ν_e$ and $\barν_e$. This is of relevance for the detection of neutrino signals from supernovae.