Source author record

H. -T. Janka

H. -T. 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

19works
6topics
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

19 published item(s)

preprint2021arXiv

Fission fragment distributions and their impact on the r-process nucleosynthesis in neutron star mergers

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. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-rich fissioning nuclei of relevance in r-process calculations. These improvements include a phenomenological modification of the scission distance and a smoothing procedure of the distribution. Such corrections lead to a much better agreement with experimental fission yields. Those yields are also used to estimate the number of neutrons emitted by the excited fragments on the basis of different neutron evaporation models. Our new fission yields are extensively compared to those predicted by the so-called GEF model. The impact of fission on the r-process nucleosynthesis in binary neutron mergers is also reanalyzed. Two scenarios are considered, the first one with low initial electron fraction is subject to intense fission recycling, in contrast to the second one which includes weak interactions on nucleons. The various regions of the nuclear chart responsible for fission recycling during the neutron irradiation as well as after freeze-out are discussed. The contribution fission processes may have to the final abundance distribution is also studied in detail in the light of newly defined quantitative indicators describing the fission recycling, the fission seeds and the fission progenitors. In particular, those allow us to estimate the contribution of fission to the final abundance distribution stemming from specific heavy nuclei. Calculations obtained with SPY and GEF fission fragment distributions are compared for both r-process scenarios.

preprint2020arXiv

Properties of gamma-ray decay lines in 3D core-collapse supernova models, with application to SN 1987A and Cas A

Comparison of theoretical line profiles to observations provides important tests for supernova explosion models. We study the shapes of radioactive decay lines predicted by current 3D core-collapse explosion simulations, and compare these to observations of SN 1987A and Cas A. Both the widths and shifts of decay lines vary by several thousand kilometers per second depending on viewing angle. The line profiles can be complex with multiple peaks. By combining observational constraints from 56Co decay lines, 44Ti decay lines, and Fe IR lines, we delineate a picture of the morphology of the explosive burning ashes in SN 1987A. For M_ZAMS=15-20 Msun progenitors exploding with ~1.5 *10^51 erg, ejecta structures suitable to reproduce the observations involve a bulk asymmetry of the 56Ni of at least ~400 km/s and a bulk velocity of at least ~1500 km/s. By adding constraints to reproduce the UVOIR bolometric light curve of SN 1987A up to 600d, an ejecta mass around 14 Msun is favoured. We also investigate whether observed decay lines can constrain the neutron star (NS) kick velocity. The model grid provides a constraint V_NS > V_redshift, and applying this to SN 1987A gives a NS kick of at least 500 km/s. For Cas A, our single model provides a satisfactory fit to the NuSTAR observations and reinforces the result that current neutrino-driven core-collapse SN models can achieve enough bulk asymmetry in the explosive burning material. Finally, we investigate the internal gamma-ray field and energy deposition, and compare the 3D models to 1D approximations.

preprint2020arXiv

The Birth Function for Black Holes and Neutron Stars in Close Binaries

The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5$\%$ to 30$\%$, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32 - 1.37 $M_\odot$ regardless of metallicity. The median black hole mass for solar metallicity is typically 8 to 9 $M_\odot$ if only initial helium cores below 40 $M_\odot$ (ZAMS mass less than 80 $M_\odot$) are counted, and 9 - 13 $M_\odot$, in most cases, if helium cores with initial masses up to 150 $M_\odot$ (ZAMS mass less than 300 $M_\odot$) contribute. As long as the mass-loss rate as a function of mass exhibits no strong non-linearities, the black hole birth function from 15 to 35 $M_\odot$ has a slope that depends mostly on the initial mass function for main sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational wave radiation. The expected rotation rates of the black holes are briefly discussed.

preprint2020arXiv

The fully developed remnant of a neutrino-driven supernova: Evolution of ejecta structure and asymmetries in SNR Cassiopeia A

Abridged. We aim at exploring to which extent the remnant keeps memory of the asymmetries that develop stochastically in the neutrino-heating layer due to hydrodynamic instabilities (e.g., convective overturn and the standing accretion shock instability) during the first second after core bounce. We coupled a 3D HD model of a neutrino-driven SN explosion with 3D MHD/HD simulations of the remnant formation. The simulations cover 2000 years of expansion and include all physical processes relevant to describe the complexities in the SN evolution and the subsequent interaction of the stellar debris with the wind of the progenitor star. The interaction of large-scale asymmetries left from the earliest phases of the explosion with the reverse shock produces, at the age of $\approx 350$~years, an ejecta structure and a remnant morphology which are remarkably similar to those observed in Cas A. Small-scale structures in the large-scale Fe-rich plumes created during the initial stages of the SN, combined with HD instabilities that develop after the passage of the reverse shock, naturally produce a pattern of ring- and crown-like structures of shocked ejecta. The consequence is a spatial inversion of the ejecta layers with Si-rich ejecta being physically interior to Fe-rich ejecta. The full-fledged remnant shows voids and cavities in the innermost unshocked ejecta resulting from the expansion of Fe-rich plumes and their inflation due to the decay of radioactive species. The asymmetric distributions of $^{44}$Ti and $^{56}$Fe and their abundance ratio are both compatible with those inferred from high-energy observations of Chandra and NuSTAR. The main asymmetries observed in the ejecta distribution of Cas A can be explained by the interaction of the reverse shock with the large-scale asymmetries that developed from stochastic processes that originate during the first seconds of the SN blast.

preprint2016arXiv

Core-Collapse Supernovae from 9 to 120 Solar Masses Based on Neutrino-powered Explosions

Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive stars with solar metallicity and masses from 9.0 to 120 solar masses. The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using several standard progenitors, and for the Crab supernova using a 9.6 solar mass progenitor. As a result of using a calibrated central engine, the final kinetic energy of the supernova is variable and sensitive to the structure of the presupernova star. Many progenitors with extended core structures do not explode, but become black holes, and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the sun provided that a reasonable contribution from Type Ia supernovae is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star IMF has a mean gravitational mass near 1.4 solar masses. The average black hole mass is about 9 solar masses if only the helium core implodes, and 14 solar masses if the entire presupernova star collapses. Only ~10% of supernovae come from stars over 20 solar masses and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.

preprint2016arXiv

Dynamics and gravitational-wave emission of neutron-star merger remnants

The coalescence of a neutron-star binary is likely to result in the formation of a neutron-star merger remnant for a large range of binary mass configurations. The massive merger remnant shows strong oscillations, which are excited by the merging process, and emits gravitational waves. Here we discuss possibilities and prospects of inferring unknown stellar properties of neutron stars by the detection of postmerger gravitational-wave emission, which thus leads to constraints of the equation of state of high-density matter. In particular, the dominant oscillation frequency of the postmerger remnant provides tight limits to neutron-star radii. We mention first steps towards a practical implementation of future gravitational-wave searches for the postmerger emission. Moreover, we outline possibilities to estimate the unknown maximum mass of nonrotating neutron stars from such types of measurements. Finally, we review the origin and scientific implications of secondary peaks in the gravitational-wave spectrum of neutron-star mergers and differences in the dynamical behavior of the postmerger remnant depending on the binary configuration. These considerations lead to a unified picture of the post-merger gravitational-wave emission and the post-merger dynamics.

preprint2015arXiv

Inferring neutron-star properties from gravitational-wave signals of binary mergers

The oscillations of a merger remnant forming after the coalescence of two neutron stars are very characteristic for the high-density equation of state. The dominant oscillation frequency occurs as a pronounced peak in the kHz range of the gravitational-wave spectrum. We describe how the dominant oscillation frequency of the remnant can be employed to infer the radii of non-rotating neutron stars.

preprint2014arXiv

Prospects For High Frequency Burst Searches Following Binary Neutron Star Coalescence With Advanced Gravitational Wave Detectors

The equation of state plays a critical role in the physics of the merger of two neutron stars. Recent numerical simulations with microphysical equation of state suggest the outcome of such events depends on the mass of the neutron stars. For less massive systems, simulations favor the formation of a hypermassive, quasi-stable neutron star, whose oscillations produce a short, high frequency burst of gravitational radiation. Its dominant frequency content is tightly correlated with the radius of the neutron star, and its measurement can be used to constrain the supranuclear equation of state. In contrast, the merger of higher mass systems results in prompt gravitational collapse to a black hole. We have developed an algorithm which combines waveform reconstruction from a morphology-independent search for gravitational wave transients with Bayesian model selection, to discriminate between post-merger scenarios and accurately measure the dominant oscillation frequency. We demonstrate the efficacy of the method using a catalogue of simulated binary merger signals in data from LIGO and Virgo, and we discuss the prospects for this analysis in advanced ground-based gravitational wave detectors. From the waveforms considered in this work and assuming an optimally oriented source, we find that the post-merger neutron star signal may be detectable by this technique to $\sim 10\text{--}25$\,Mpc. We also find that we successfully discriminate between the post-merger scenarios with $\sim 95\%$ accuracy and determine the dominant oscillation frequency of surviving post-merger neutron stars to within $\sim 10$\,Hz, averaged over all detected signals. This leads to an uncertainty in the estimated radius of a non-rotating 1.6\,M$_{\odot}$ reference neutron star of $\sim 100\,$m.

preprint2014arXiv

Revealing the high-density equation of state through binary neutron star mergers

We present a novel method for revealing the equation of state of high-density neutron star matter through gravitational waves emitted during the postmerger phase of a binary neutron star system. The method relies on a small number of detections of the peak frequency in the postmerger phase for binaries of different (relatively low) masses, in the most likely range of expected detections. From such observations, one can construct the derivative of the peak frequency versus the binary mass, in this mass range. Through a detailed study of binary neutron star mergers for a large sample of equations of state, we show that one can extrapolate the above information to the highest possible mass (the threshold mass for black hole formation in a binary neutron star merger). In turn, this allows for an empirical determination of the maximum mass of cold, nonrotating neutron stars to within 0.1 M_sun, while the corresponding radius is determined to within a few percent. Combining this with the determination of the radius of cold, nonrotating neutron stars of 1.6 M_sun (to within a few percent, as was demonstrated in Bauswein et al., PRD, 86, 063001, 2012), allows for a clear distinction of a particular candidate equation of state among a large set of other candidates. Our method is particularly appealing because it reveals simultaneously the moderate and very high-density parts of the equation of state, enabling the distinction of mass-radius relations even if they are similar at typical neutron star masses. Furthermore, our method also allows to deduce the maximum central energy density and maximum central rest-mass density of cold, nonrotating neutron stars with an accuracy of a few per cent.

preprint2013arXiv

Prompt merger collapse and the maximum mass of neutron stars

We perform hydrodynamical simulations of neutron-star mergers for a large sample of temperature-dependent, nuclear equations of state, and determine the threshold mass above which the merger remnant promptly collapses to form a black hole. We find that, depending on the equation of state, the threshold mass is larger than the maximum mass of a non-rotating star in isolation by between 30 and 70 per cent. Our simulations also show that the ratio between the threshold mass and maximum mass is tightly correlated with the compactness of the non-rotating maximum-mass configuration. We speculate on how this relation can be used to derive constraints on neutron-star properties from future observations.

preprint2013arXiv

Systematics of dynamical mass ejection, nucleosynthesis, and radioactively powered electromagnetic signals from neutron-star mergers

We investigate systematically the dynamical mass ejection, r-process nucleosynthesis, and properties of electromagnetic counterparts of neutron-star (NS) mergers in dependence on the uncertain properties of the nuclear equation of state (EoS) by employing 40 representative, microphysical high-density EoSs in relativistic, hydrodynamical simulations. The crucial parameter determining the ejecta mass is the radius R_1.35 of a 1.35 M_sun NS. NSs with smaller R_1.35 ("soft" EoS) eject systematically higher masses. These range from ~10^-3 M_sun to ~10^-2 M_sun for 1.35-1.35 M_sun binaries and from ~5*10^-3 M_sun to ~2*10^-2 M_sun for 1.2-1.5 M_sun systems (with kinetic energies between ~5*10^49 erg and 10^51 erg). Correspondingly, the bolometric peak luminosities of the optical transients of symmetric (asymmetric) mergers vary between 3*10^41 erg/s and 14*10^41 erg/s (9*10^41 erg/s and 14.5*10^41 erg/s) on timescales between ~2 h and ~12 h. If these signals with absolute bolometric magnitudes from -15.0 to -16.7 are measured, the tight correlation of their properties with those of the merging NSs might provide valuable constraints on the high-density EoS. The r-process nucleosynthesis exhibits a remarkable robustness independent of the EoS, producing a nearly solar abundance pattern above mass number 130. By the r-process content of the Galaxy and the average production per event the Galactic merger rate is limited to 4*10^-5/yr (4*10^-4/yr) for a soft (stiff) NS EoS, if NS mergers are the main source of heavy r-nuclei. The production ratio of radioactive 232Th to 238U attains a stable value of 1.64-1.67, which does not exclude NS mergers as potential sources of heavy r-material in the most metal-poor stars.

preprint2012arXiv

Equation-of-state dependence of the gravitational-wave signal from the ring-down phase of neutron-star mergers

Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.

preprint2011arXiv

The decompression of the outer neutron star crust and r-process nucleosynthesis

The rapid neutron-capture process, or r-process, is known to be fundamental for explaining the origin of approximately half of the A>60 stable nuclei observed in nature. In recent years nuclear astrophysicists have developed more and more sophisticated r-process models, by adding new astrophysical or nuclear physics ingredients to explain the solar system composition in a satisfactory way. Despite these efforts, the astrophysical site of the r-process remains unidentified. The composition of the neutron star outer crust material is investigated after the decompression that follows its possible ejection. The composition of the outer crust of a neutron star is estimated before and after decompression. Two different possible initial conditions are considered, namely an idealized crust composed of cold catalyzed matter and a crust initially in nuclear statistical equilibrium at temperatures around 10 GK. We show that in this second case before decompression and at temperatures typically corresponding to 8 GK, the Coulomb effect due to the high densities in the crust leads to an overall composition of the outer crust in neutron-rich nuclei with a mass distribution close to the solar system r-abundance distribution. Such distributions differ, however, from the solar one due to a systematic shift in the second peak to lower values. After decompression, the capture of the few neutrons per seed nucleus available in the hot outer crust leads to a final distribution of stable neutron-rich nuclei with a mass distribution of 80 < A < 140 nuclei in excellent agreement with the solar distribution, provided the outer crust is initially at temperatures around 8 GK and all layers of the outer crust are ejected. The decompression of the neutron star matter from the outer crust provides suitable conditions for a robust r-processing of the light species, i.e., r-nuclei with A < 140.

preprint2010arXiv

Nucleosynthesis-relevant conditions in neutrino-driven supernova outflows. II. The reverse shock in two-dimensional simulations

After the initiation of the explosion of core-collapse supernovae, neutrinos emitted from the nascent neutron star drive a supersonic baryonic outflow. This neutrino-driven wind interacts with the more slowly moving, earlier supernova ejecta forming a wind termination shock (or reverse shock), which changes the local wind conditions and their evolution. Important nucleosynthesis processes (alpha-process, charged-particle reactions, r-process, and vp-process) occur or might occur in this environment. The nucleosynthesis depends on the long-time evolution of density, temperature, and expansion velocity. Here we present two-dimensional hydrodynamical simulations with an approximate description of neutrino-transport effects, which for the first time follow the post-bounce accretion, onset of the explosion, wind formation, and the wind expansion through the collision with the preceding supernova ejecta. Our results demonstrate that the anisotropic ejecta distribution has a great impact on the position of the reverse shock, the wind profile, and the long-time evolution. This suggests that hydrodynamic instabilities after core bounce and the consequential asymmetries may have important effects on the nucleosynthesis-relevant conditions in the neutrino-heated baryonic mass flow from proto-neutron stars.

preprint2009arXiv

Axisymmetric General Relativistic Simulations of the Accretion-Induced Collapse of White Dwarfs

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

preprint2008arXiv

Nucleosynthesis in Early Neutrino Driven Winds

Nucleosynthesis in early neutrino winds is investigated. Presented is a brief overview of two recent problems of supernova nucleosynthesis. In the first part we investigate the effect of nuclear parameters on the synthesis of Mo92 and Mo94. Based on recent experimental results, we find that the proton rich winds of the model investigated here can not be the source of solar Mo92 and Mo94. In the second part we investigate the nucleosynthesis from neutron rich bubbles and show that they do not contribute to the overall nucleosynthesis.

preprint2006arXiv

The Physics of Core-Collapse Supernovae

Supernovae are nature's grandest explosions and an astrophysical laboratory in which unique conditions exist that are not achievable on Earth. They are also the furnaces in which most of the elements heavier than carbon have been forged. Scientists have argued for decades about the physical mechanism responsible for these explosions. It is clear that the ultimate energy source is gravity, but the relative roles of neutrinos, fluid instabilities, rotation and magnetic fields continue to be debated.

preprint2004arXiv

Relativistic outflows from remnants of compact object mergers and their viability for short gamma-ray bursts

We present the first general relativistic hydrodynamic models of the launch and evolution of relativistic jets and winds, driven by thermal energy deposition, possibly due to neutrino-antineutrino annihilation, in the close vicinity of black hole-accretion torus systems. The latter are considered to be the remnants of compact object mergers. Our two-dimensional simulations establish the link between such mergers and future observations of short gamma-ray bursts (GRBs) by the SWIFT satellite. They show that ultrarelativistic outflow with maximum terminal Lorentz factors (Gamma) around 1000 develops for polar energy deposition rates above some 1e48 erg/s per steradian, provided the merger environment has a sufficiently low baryon density. Due to the collimation by the dense accretion torus the typical semi-opening angles of the Gamma > 100 cone are 5-10 degrees, corresponding to about 0.4-1.5% of the hemisphere and apparent isotropized energies (kinetic plus internal) up to ~1e51 erg. 10-30% of the deposited energy are transferred to the outflow with Gamma > 100. Our models confirm the viability of post-merger BH-torus systems as engines of short, hard GRBs and can explain the durations of all observed short GRBs, because different propagation velocities of the front and rear ends lead to a radial stretching of the ultrarelativistic fireball before transparency is reached. The ultrarelativistic flow reveals a highly non-uniform structure with Lorentz factor variations up to factors of a few, caused by the action of Kelvin-Helmholtz instabilities that originate at the fireball-torus interface (abbreviated).

preprint1998arXiv

No Pulsar Kicks from Deformed Neutrinospheres

In a supernova core, magnetic fields cause a directional variation of the neutrino refractive index so that resonant flavor oscillations would lead to a deformation of the "neutrinosphere" for, say, tau neutrinos. The associated anisotropic neutrino emission was proposed as a possible origin of the observed pulsar proper motions. We argue that this effect was vastly overestimated because the variation of the temperature over the deformed neutrinosphere is not an adequate measure for the anisotropy of neutrino emission. The neutrino flux is generated inside the neutron star core and is transported through the atmosphere at a constant luminosity, forcing the temperature gradient in the atmosphere to adjust to the inflow of energy from below. Therefore, no emission anisotropy is caused by a deformation of the neutrinosphere to lowest order. An estimate of the higher-order corrections must take into account the modified atmospheric temperature profile in response to the deformation of the neutrinosphere and the corresponding feedback on the core. We go through this exercise in the framework of a simplified model which can be solved analytically.