BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

H. -Th. Janka

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

ResearcherUnclaimed source record
astro-ph.SRastro-phastro-ph.HEgr-qchep-phhep-exnucl-thnucl-exphysics.ins-det

Catalog footprint

What is connected

31works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2022arXiv

Low-Energy Physics in Neutrino LArTPCs

In this white paper, we outline some of the scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) detectors. Key takeaways are summarized as follows. 1) LArTPCs have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below the few tens of MeV range. 2) Low-energy signatures are an integral part of GeV-scale accelerator neutrino interaction final states, and their reconstruction can enhance the oscillation physics sensitivities of LArTPC experiments. 3) BSM signals from accelerator and natural sources also generate diverse signatures in the low-energy range, and reconstruction of these signatures can increase the breadth of BSM scenarios accessible in LArTPC-based searches. 4) Neutrino interaction cross sections and other nuclear physics processes in argon relevant to sub-hundred-MeV LArTPC signatures are poorly understood. Improved theory and experimental measurements are needed. Pion decay-at-rest sources and charged particle and neutron test beams are ideal facilities for experimentally improving this understanding. 5) There are specific calibration needs in the low-energy range, as well as specific needs for control and understanding of radiological and cosmogenic backgrounds. 6) Novel ideas for future LArTPC technology that enhance low-energy capabilities should be explored. These include novel charge enhancement and readout systems, enhanced photon detection, low radioactivity argon, and xenon doping. 7) Low-energy signatures, whether steady-state or part of a supernova burst or larger GeV-scale event topology, have specific triggering, DAQ and reconstruction requirements that must be addressed outside the scope of conventional GeV-scale data collection and analysis pathways.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Three-dimensional Models of Core-collapse Supernovae From Low-mass Progenitors With Implications for Crab

We present 3D full-sphere supernova simulations of non-rotating low-mass (~9 Msun) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions [~(0.5-1.0)x 10^{50} erg] of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen-neon-magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the Vertex-Prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the Prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less 2nd dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of <~5x10^{-3} Msun have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly-born NS. The LMCCSNe model with less 2nd dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km/s and a NS spin period of ~30 ms, both not largely different from those of the Crab pulsar at birth.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A two-parameter criterion for classifying the explodability of massive stars by the neutrino-driven mechanism

Thus far, judging the fate of a massive star (either a neutron star (NS) or a black hole) solely by its structure prior to core collapse has been ambiguous. Our work and previous attempts find a non-monotonic variation of successful and failed supernovae with zero-age main-sequence mass, for which no single structural parameter can serve as a good predictive measure. However, we identify two parameters computed from the pre-collapse structure of the progenitor, which in combination allow for a clear separation of exploding and non-exploding cases with only few exceptions (~1-2.5%) in our set of 621 investigated stellar models. One parameter is M4, defining the normalized enclosed mass for a dimensionless entropy per nucleon of s=4, and the other is mu4 = d(m/M_sun)/d(r/1000 km) at s=4, being the normalized mass-derivative at this location. The two parameters mu4 and M4*mu4 can be directly linked to the mass-infall rate, Mdot, of the collapsing star and the electron-type neutrino luminosity of the accreting proto-NS, L_nue ~ M_ns*Mdot, which play a crucial role in the "critical luminosity" concept for the theoretical description of neutrino-driven explosions as runaway phenomenon of the stalled accretion shock. All models were evolved employing the approach of Ugliano et al. for simulating neutrino-driven explosions in spherical symmetry. The neutrino emission of the accretion layer is approximated by a gray transport solver, while the uncertain neutrino emission of the 1.1 M_sun proto-NS core is parametrized by an analytic model. The free parameters connected to the core-boundary prescription are calibrated to reproduce the observables of Supernova 1987A for five different progenitor models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Solar r-process-constrained actinide production in neutrino-driven winds of supernovae

Long-lived radioactive nuclei play an important role as nucleo-cosmochronometers and as cosmic tracers of nucleosynthetic source activity. In particular nuclei in the actinide region like thorium, uranium, and plutonium can testify to the enrichment of an environment by the still enigmatic astrophysical sources that are responsible for the production of neutron-rich nuclei by the rapid neutron-capture process (r-process). Supernovae and merging neutron-star (NS) or NS-black hole binaries are considered as most likely sources of the r-nuclei. But arguments in favour of one or the other or both are indirect and make use of assumptions; they are based on theoretical models with remaining simplifications and shortcomings. An unambiguous observational determination of a production event is still missing. In order to facilitate searches in this direction, e.g.\ by looking for radioactive tracers in stellar envelopes, the interstellar medium or terrestrial reservoirs, we provide improved theoretical estimates and corresponding uncertainty ranges for the actinide production ($^{232}$Th, $^{235,236,238}$U, $^{237}$Np, $^{244}$Pu, and $^{247}$Cm) in neutrino-driven winds of core-collapse supernovae. Since state-of-the-art supernova models do not yield r-process viable conditions ---but still lack, for example, the effects of strong magnetic fields--- we base our investigation on a simple analytical, Newtonian, adiabatic and steady-state wind model and consider the superposition of a large number of contributing components, whose nucleosynthesis-relevant parameters (mass weight, entropy, expansion time scale, and neutron excess) are constrained by the assumption that the integrated wind nucleosynthesis closely reproduces the solar system distribution of r-process elements. We also test the influence of uncertain nuclear physics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-Radial Instabilities and Progenitor Asphericities in Core-Collapse Supernovae

Since core-collapse supernova simulations still struggle to produce robust neutrino-driven explosions in 3D, it has been proposed that asphericities caused by convection in the progenitor might facilitate shock revival by boosting the activity of non-radial hydrodynamic instabilities in the post-shock region. We investigate this scenario in depth using 42 relativistic 2D simulations with multi-group neutrino transport to examine the effects of velocity and density perturbations in the progenitor for different perturbation geometries that obey fundamental physical constraints (like the anelastic condition). As a framework for analysing our results, we introduce semi-empirical scaling laws relating neutrino heating, average turbulent velocities in the gain region, and the shock deformation in the saturation limit of non-radial instabilities. The squared turbulent Mach number, <Ma^2>, reflects the violence of aspherical motions in the gain layer, and explosive runaway occurs for <Ma^2>~0.3, corresponding to a reduction of the critical neutrino luminosity by ~25% compared to 1D. In the light of this theory, progenitor asphericities aid shock revival mainly by creating anisotropic mass flux onto the shock: Differential infall efficiently converts velocity perturbations in the progenitor into density perturbations (Delta rho/rho) at the shock of the order of the initial convective Mach number Ma. The anisotropic mass flux and ram pressure deform the shock and thereby amplify post-shock turbulence. Large-scale (l=2,l=1) modes prove most conducive to shock revival, whereas small-scale perturbations require unrealistically high convective Mach numbers. Initial density perturbations in the progenitor are only of order Ma^2 and therefore play a subdominant role.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A New Multi-Dimensional General Relativistic Neutrino Hydrodynamics Code for Core-Collapse Supernovae IV. The Neutrino Signal

Considering general relativistic, two-dimensional (2D) supernova (SN) explosion models of progenitor stars between 8.1 and 27 solar masses, 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 of electron antineutrinos and heavy-lepton neutrinos and even their crossing during the accretion phase for stars with M>10 M_sun as observed in previous 1D and 2D simulations with state-of-the-art neutrino transport. We establish a roughly linear scaling of the electron antineutrino mean energy 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 ~10kpc 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-50Hz, and post-explosion, episodic accretion downflows will be signaled by activity intervals stretching over an extended frequency range in the wavelet spectrogram.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Production of Nickel-56 in black hole-neutron star merger accretion disk outflows

The likely outcome of a compact object merger event is a central black hole surrounded by a rapidly accreting torus of debris. This disk of debris is a rich source of element synthesis, the outcome of which is needed to predict electromagnetic counterparts of individual events and to understand the contribution of mergers to galactic chemical evolution. Here we study disk outflow nucleosynthesis in the context of a two-dimensional, time-dependent black hole-neutron star merger accretion disk model. We use two time snapshots from this model to examine the impact of the evolution of the neutrino fluxes from the disk on the element synthesis. While the neutrino fluxes from the early-time disk snapshot appear to favor neutron-rich outflows, by the late-time snapshot the situation is reversed. As a result we find copious production of Nickel-56 in the outflows.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SASI Activity in Three-Dimensional Neutrino-Hydrodynamics Simulations of Supernova Cores

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three-dimensional neutrino-driven supernovae: Neutron star kicks, spins, and asymmetric ejection of nucleosynthesis products

We present 3D simulations of supernova (SN) explosions of nonrotating stars, triggered by the neutrino-heating mechanism with a suitable choice of the core-neutrino luminosity. Our results show that asymmetric mass ejection caused by hydrodynamic instabilities can accelerate the neutron star (NS) up to recoil velocities of more than 700 km/s by the "gravitational tug-boat mechanism", which is enough to explain most observed pulsar velocities. The associated NS spin periods are about 100 ms to 8 s without any correlation between spin and kick magnitudes or directions. This suggests that faster spins and a possible spin-kick alignment might require angular momentum in the progenitor core prior to collapse. Our simulations for the first time demonstrate a clear correlation between the size of the NS kick and anisotropic ejection of heavy elements created by explosive burning behind the shock. In the case of large NS kicks the explosion is significantly stronger opposite to the kick vector. Therefore the bulk of the Fe-group elements, in particular nickel, is ejected mostly in large clumps against the kick direction. This contrasts with the case of low recoil velocity, where the Ni-rich lumps are more isotropically distributed. Intermediate-mass nuclei heavier than Si (like Ca and Ti) also exhibit a significant enhancement in the hemisphere opposite to the direction of fast NS motion, while the distribution of C, O, and Ne is not affected, and that of Mg only marginally. Mapping the spatial distribution of the heavy elements in SN remnants with identified pulsar motion may offer an important diagnostic test of the kick mechanism. Different from kick scenarios based on anisotropic neutrino emission, our hydrodynamical acceleration model predicts enhanced ejection of Fe-group elements and of their nuclear precursors in the direction opposite to the NS recoil. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A New Multi-Dimensional General Relativistic Neutrino Hydrodynamics Code for Core-Collapse Supernovae II. Relativistic Explosion Models of Core-Collapse Supernovae

We present the first two-dimensional general relativistic (GR) simulations of stellar core collapse and explosion with the CoCoNuT hydrodynamics code in combination with the VERTEX solver for energy-dependent, three-flavor neutrino transport, using the extended conformal flatness condition for approximating the spacetime metric and a ray-by-ray-plus ansatz to tackle the multi-dimensionality of the transport. For both of the investigated 11.2 and 15 solar mass progenitors we obtain successful, though seemingly marginal, neutrino-driven supernova explosions. This outcome and the time evolution of the models basically agree with results previously obtained with the PROMETHEUS hydro solver including an approximative treatment of relativistic effects by a modified Newtonian potential. However, GR models exhibit subtle differences in the neutrinospheric conditions compared to Newtonian and pseudo-Newtonian simulations. These differences lead to significantly higher luminosities and mean energies of the radiated electron neutrinos and antineutrinos and therefore to larger energy-deposition rates and heating efficiencies in the gain layer with favorable consequences for strong non-radial mass motions and ultimately for an explosion. Moreover, energy transfer to the stellar medium around the neutrinospheres through nucleon recoil in scattering reactions of heavy-lepton neutrinos also enhances the mentioned effects. Together with previous pseudo-Newtonian models the presented relativistic calculations suggest that the treatment of gravity and energy-exchanging neutrino interactions can make differences of even 50-100% in some quantities and is likely to contribute to a finally successful explosion mechanism on no minor level than hydrodynamical differences between different dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

New Two-Dimensional Models of Supernova Explosions by the Neutrino-Heating Mechanism: Evidence for Different Instability Regimes in Collapsing Stellar Cores

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible candidates, convection and the so-called standing accretion shock instability (SASI), have been proposed for this role. In this paper, we discuss new successful simulations of supernova explosions that shed light on the relative importance of these two instabilities. While convection has so far been observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, we here present the first such simulation in which the SASI grows faster while the development of convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an explosion model of a 27 solar mass progenitor, which is contrasted with a convectively-dominated model of an 8.1 solar mass progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with previous 11.2 and 15 solar mass explosion models. We analyze the conditions discriminating between the two different regimes, showing that a high mass-accretion rate and a short advection time-scale are conducive for strong SASI activity. We also briefly discuss some important factors for capturing the SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state leading to a compact proto-neutron star, and the neutrino treatment. Finally, we evaluate possible implications of our findings for 2D and 3D supernova simulations. Our results show that a better understanding of the SASI and convection in the non-linear regime is required.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Measuring neutron-star properties via gravitational waves from binary mergers

We demonstrate by a large set of merger simulations for symmetric binary neutron stars (NSs) that there is a tight correlation between the frequency peak of the postmerger gravitational-wave (GW) emission and the physical properties of the nuclear equation of state (EoS), e.g. expressed by the radius of the maximum-mass Tolman-Oppenheimer-Volkhoff configuration. Therefore, a single measurement of the peak frequency of the postmerger GW signal will constrain the NS EoS significantly. For plausible optimistic merger-rate estimates a corresponding detection with Advanced LIGO is likely to happen within an operation time of roughly a year.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Neutrino Signal of Electron-Capture Supernovae from Core Collapse to Cooling

An 8.8 solar mass electron-capture supernova (SN) was simulated in spherical symmetry consistently from collapse through explosion to nearly complete deleptonization of the forming neutron star. The evolution time of about 9 s is short because of nucleon-nucleon correlations in the neutrino opacities. After a brief phase of accretion-enhanced luminosities (~200 ms), luminosity equipartition among all species becomes almost perfect and the spectra of electron antineutrinos and muon/tau antineutrinos very similar. We discuss consequences for the neutrino-driven wind as a nucleosynthesis site and for flavor oscillations of SN neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Parametrized 3D models of neutrino-driven supernova explosions: Neutrino emission asymmetries and gravitational-wave signals

Time-dependent and direction-dependent neutrino and gravitational-wave (GW) signatures are presented for a set of 3D hydrodynamic models of parametrized, neutrino-driven supernova explosions of non-rotating 15 and 20 solar mass stars. We employ an approximate treatment of neutrino transport. Due to the excision of the high-density core of the proto-neutron star and the use of an axis-free overset grid, the models can be followed from the post-bounce accretion phase for more than one second without imposing any symmetry restrictions. GW and neutrino emission exhibit the generic time-dependent features known from 2D models. Non-radial hydrodynamic mass motions in the accretion layer and their interaction with the outer layers of the proto-neutron star together with anisotropic neutrino emission give rise to a GW signal with an amplitude of ~5-20 cm and frequencies 100--500 Hz. The GW emission from mass motions reaches a maximum before the explosion sets in. Afterwards the GW signal exhibits a low-frequency modulation, in some cases describing a quasi-monotonic growth, associated with the non-spherical expansion of the explosion shock wave and the large-scale anisotropy of the escaping neutrino flow. Variations of the mass-quadrupole moment due to convective activity inside the nascent neutron star contribute a high-frequency component to the GW signal during the post-explosion phase. The GW signals exhibit strong variability between the two polarizations, different explosion simulations and different observer directions, and does not possess any template character. The neutrino emission properties show fluctuations over the neutron star surface on spatial and temporal scales that reflect the different types of non-spherical mass motions. The modulation amplitudes of the measurable neutrino luminosities and mean energies are significantly smaller than predicted by 2D simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. I. Method and code tests in spherical symmetry

We present a new general relativistic (GR) code for hydrodynamic supernova simulations with neutrino transport in spherical and azimuthal symmetry (1D/2D). The code is a combination of the CoCoNuT hydro module, which is a Riemann-solver based, high-resolution shock-capturing method, and the three-flavor, energy-dependent neutrino transport scheme VERTEX. VERTEX integrates the neutrino moment equations with a variable Eddington factor closure computed from a model Boltzmann equation and uses the ray-by-ray plus approximation in 2D, assuming the neutrino distribution to be axially symmetric around the radial direction, and thus the neutrino flux to be radial. Our spacetime treatment employs the ADM 3+1 formalism with the conformal flatness condition for the spatial three-metric. This approach is exact in 1D and has been shown to yield very accurate results also for rotational stellar collapse. We introduce new formulations of the energy equation to improve total energy conservation in relativistic and Newtonian hydro simulations with Eulerian finite-volume codes. Moreover, a modified version of the VERTEX scheme is developed that simultaneously conserves energy and lepton number with better accuracy and higher numerical stability. To verify our code, we conduct a series of tests, including a detailed comparison with published 1D results for stellar core collapse. Long-time simulations of proto-neutron star cooling over several seconds both demonstrate the robustness of the new CoCoNuT-VERTEX code and show the approximate treatment of GR effects by means of an effective gravitational potential as in PROMETHEUS-VERTEX to be remarkably accurate in 1D. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Discriminating Strange Star Mergers from Neutron Star Mergers by Gravitational-Wave Measurements

We perform three-dimensional relativistic hydrodynamical simulations of the coalescence of strange stars (SSs) and explore the possibility to decide on the strange matter hypothesis by means of gravitational-wave (GW) measurements. Selfbinding of strange quark matter (SQM) and the generally more compact stars yield features that clearly distinguish SS from neutron star (NS) mergers, e.g. hampering tidal disruption during the plunge of quark stars. Furthermore, instead of forming dilute halos around the remnant as in the case of NS mergers, the coalescence of SSs results in a differentially rotating hypermassive object with a sharp surface layer surrounded by a geometrically thin, clumpy high-density SQM disk. We also investigate the importance of including non-zero temperature equations of state (EoSs) in NS and SS merger simulations. In both cases we find a crucial sensitivity of the dynamics and outcome of the coalescence to thermal effects, which, e.g., determine the outer remnant structure and the delay time of the dense remnant core to black hole collapse. For comparing and classifying the GW signals, we use a number of characteristic quantities like the maximum frequency during inspiral or the dominant frequency of oscillations of the postmerger remnant. In general, these frequencies are higher for SS mergers. If not, additional features of the GW luminosity spectrum may help to discriminate coalescence events of the different types. Future GW measurements may thus help to decide on the existence of SQM stars. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hydrodynamical Neutron Star Kicks in Three Dimensions

Using three-dimensional (3D) simulations of neutrino-powered supernova explosions we show that the hydrodynamical kick scenario proposed by Scheck et al. on the basis of two-dimensional (2D) models can yield large neutron star (NS) recoil velocities also in 3D. Although the shock stays relatively spherical, standing accretion-shock and convective instabilities lead to a globally asymmetric mass and energy distribution in the postshock layer. An anisotropic momentum distribution of the ejecta is built up only after the explosion sets in. Total momentum conservation implies the acceleration of the NS on a timescale of 1-3 seconds, mediated mainly by long-lasting, asymmetric accretion downdrafts and the anisotropic gravitational pull of large inhomogeneities in the ejecta. In a limited set of 15 solar-mass models with an explosion energy of about 10^51 erg this stochastic mechanism is found to produce kicks from <100 km/s to >500 km/s, and >1000 km/s seem possible. Strong rotational flows around the accreting NS do not develop in our collapsing, non-rotating progenitors. The NS spins therefore remain low with estimated periods of about 500-1000 ms and no alignment with the kicks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polytropic neutron star - black hole merger simulations with a Paczynski-Wiita potential

Context: Mergers of neutron stars (NS) and black holes (BH) are among the strongest sources of gravitational waves and are potential central engines for short gamma-ray bursts. Aims: We aim to compare the general relativistic (GR) results by other groups with Newtonian calculations of models with equivalent parameters. We vary the mass ratios between NS and BH and the compactness of the NS. The mass of the NS is 1.4 M_sol. We compare the dynamics in the parameter-space regions where the NS is expected to reach the innermost stable circular orbit (ISCO) before being tidally disrupted (mass shedding, MS) and vice versa. Methods: The hydrodynamics is evolved by a Newtonian PPM scheme with four levels of nested grids. We use a polytropic EoS (Gamma=2), as was done in the GR simulations. However, instead of full GR we use a Newtonian potential supplemented by a Paczynski-Wiita-Artemova potential for the BH, both disregarding and including rotation of the BH. Results: If the NS is compact (C=0.18) it is accreted by the BH more quickly, and only a small amount of mass remains outside the BH. If the mass ratio is small (Q=2 or 3) or the NS is less compact (C=0.16 or less) the NS is tidally torn apart before being accreted. Although most of the mass is absorbed by the BH, some 0.1 M_sol remain in a tidal arm. For small mass ratios the tidal arm can wrap around the BH to form a thick disk. When including the effects of BH spin-up or spin-down by the accreted matter, more mass remains in the surroundings (0.2-0.3 M_sol). Conclusions: Although details and quantitative results differ, the general trends of our Newtonian calculations are similar to the GR calculations. A clear delimiting line that separates ISCO from the MS cases is not found. Inclusion of BH rotation as well as sufficient numerical resolution are extremely important.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing Approximations of Thermal Effects in Neutron Star Merger Simulations

We perform three-dimensional relativistic hydrodynamical calculations of neutron star mergers to assess the reliability of an approximate treatment of thermal effects in such simulations by combining an ideal-gas component with zero-temperature, micro-physical equations of state. To this end we compare the results of simulations that make this approximation to the outcome of models with a consistent treatment of thermal effects in the equation of state. In particular we focus on the implications for observable consequences of merger events like the gravitational-wave signal. It is found that the characteristic gravitational-wave oscillation frequencies of the post-merger remnant differ by about 50 to 250 Hz (corresponding to frequency shifts of 2 to 8 per cent) depending on the equation of state and the choice of the characteristic index of the ideal-gas component. In addition, the delay time to black hole collapse of the merger remnant as well as the amount of matter remaining outside the black hole after its formation are sensitive to the description of thermal effects.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Three-Dimensional Simulations of Mixing Instabilities in Supernova Explosions

We present the first three-dimensional (3D) simulations of the large-scale mixing that takes place in the shock-heated stellar layers ejected in the explosion of a 15.5 solar-mass blue supergiant star. The outgoing supernova shock is followed from its launch by neutrino heating until it breaks out from the stellar surface more than two hours after the core collapse. Violent convective overturn in the post-shock layer causes the explosion to start with significant asphericity, which triggers the growth of Rayleigh-Taylor (RT) instabilities at the composition interfaces of the exploding star. Deep inward mixing of hydrogen (H) is found as well as fast-moving, metal-rich clumps penetrating with high velocities far into the H-envelope of the star as observed, e.g., in the case of SN 1987A. Also individual clumps containing a sizeable fraction of the ejected iron-group elements (up to several 0.001 solar masses) are obtained in some models. The metal core of the progenitor is partially turned over with Ni-dominated fingers overtaking oxygen-rich bullets and both Ni and O moving well ahead of the material from the carbon layer. Comparing with corresponding 2D (axially symmetric) calculations, we determine the growth of the RT fingers to be faster, the deceleration of the dense metal-carrying clumps in the He and H layers to be reduced, the asymptotic clump velocities in the H-shell to be higher (up to ~4500 km/s for the considered progenitor and an explosion energy of 10^{51} ergs, instead of <2000 km/s in 2D), and the outward radial mixing of heavy elements and inward mixing of hydrogen to be more efficient in 3D than in 2D. We present a simple argument that explains these results as a consequence of the different action of drag forces on moving objects in the two geometries. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nucleosynthesis in Electron Capture Supernovae of AGB Stars

We examine nucleosynthesis in the electron capture supernovae of progenitor AGB stars with an O-Ne-Mg core (with the initial stellar mass of 8.8 M_\odot). Thermodynamic trajectories for the first 810 ms after core bounce are taken from a recent state-of-the-art hydrodynamic simulation. The presented nucleosynthesis results are characterized by a number of distinct features that are not shared with those of other supernovae from the collapse of stars with iron core (with initial stellar masses of more than 10 M_\odot). First is the small amount of 56Ni (= 0.002-0.004 M_\odot) in the ejecta, which can be an explanation for observed properties of faint supernovae such as SNe 2008S and 1997D. In addition, the large Ni/Fe ratio is in reasonable agreement with the spectroscopic result of the Crab nebula (the relic of SN 1054). Second is the large production of 64Zn, 70Ge, light p-nuclei (74Se, 78Kr, 84Sr, and 92Mo), and in particular, 90Zr, which originates from the low Y_e (= 0.46-0.49, the number of electrons per nucleon) ejecta. We find, however, that only a 1-2% increase of the minimum Y_e moderates the overproduction of 90Zr. In contrast, the production of 64Zn is fairly robust against a small variation of Y_e. This provides the upper limit of the occurrence of this type of events to be about 30% of all core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

r-Process Nucleosynthesis in Hot Accretion Disk Flows from Black Hole - Neutron Star Mergers

We consider hot accretion disk outflows from black hole - neutron star mergers in the context of the nucleosynthesis they produce. We begin with a three dimensional numerical model of a black hole - neutron star merger and calculate the neutrino and antineutrino fluxes emitted from the resulting accretion disk. We then follow the element synthesis in material outflowing the disk along parameterized trajectories. We find that at least a weak r-process is produced, and in some cases a main r-process as well. The neutron-rich conditions required for this production of r-process nuclei stem directly from the interactions of the neutrinos emitted by the disk with the free neutrons and protons in the outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Effects of Inelastic Neutrino-Nucleus Scattering on Supernova Dynamics and Radiated Neutrino Spectra

Based on the shell model for Gamow-Teller and the Random Phase Approximation for forbidden transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. The rates have been 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%.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Explosions of O-Ne-Mg Cores, the Crab Supernova, and Subluminous Type II-P Supernovae

We present results of simulations of stellar collapse and explosions in spherical symmetry for progenitor stars in the 8-10 solar mass range with an O-Ne-Mg core. The simulations were continued until nearly one second after core bounce and were performed with the Prometheus/Vertex code with a variable Eddington factor solver for the neutrino transport, including a state-of-the-art treatment of neutrino-matter interactions. Particular effort was made to implement nuclear burning and electron capture rates with sufficient accuracy to ensure a smooth continuation, without transients, from the progenitor evolution to core collapse. Using two different nuclear equations of state (EoSs), a soft version of the Lattimer & Swesty EoS and the significantly stiffer Wolff & Hillebrandt EoS, we found no prompt explosions, but instead delayed explosions, powered by neutrino heating and the neutrino-driven baryonic wind which sets in about 200 ms after bounce. The models eject little nickel (< 0.015 solar masses), explode with an energy of about or slightly more than 10**50 erg, and leave behind neutron stars (NSs) with a baryonic mass near 1.36 solar masses. Different from previous models of such explosions, the ejecta during the first second have a proton-to-baryon ratio of Ye > 0.46, which suggests a chemical composition that is not in conflict with galactic abundances. No low-entropy matter with Ye << 0.5 is ejected. This excludes such explosions as sites of a low-entropy r-process. The low explosion energy and nucleosynthetic implications are compatible with the observed properties of the Crab supernova, and the small nickel mass supports the possibility that our models explain some subluminous Type II-P supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Exploiting the neutronization burst of a galactic supernova

One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt neutronization burst, i.e. the first $\sim 25$ milliseconds after bounce when the SN emits with very high luminosity mainly $ν_e$ neutrinos. We examine the dependence of this burst on variations in the input of current SN models and find that recent improvements of the electron capture rates as well as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle $θ_{13}$, where the peak is absent. The robustness of the predicted total event number in the neutronization burst makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the likely case that the SN is optically obscured.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Supernova Simulations with Boltzmann Neutrino Transport: A Comparison of Methods

Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, Agile-Boltztran of the Oak Ridge-Basel group and Vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 solar mass star with Newtonian gravity and a 15 solar mass star with relativistic gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Non-spherical Core Collapse Supernovae I. Neutrino-Driven Convection, Rayleigh-Taylor Instabilities, and the Formation and Propagation of Metal Clumps

Two-dimensional simulations of a Type II and a Type Ib-like supernova explosion are presented that encompass shock revival by neutrino heating, neutrino-driven convection, explosive nucleosynthesis, the growth of Rayleigh-Taylor instabilities, and the propagation of newly formed metal clumps through the exploding star. In both cases we find very high Ni56 velocities of 17000 km/s shortly after shock-revival, and a complete fragmentation of the progenitor's metal core within the first few minutes after core bounce, due to the growth of Rayleigh-Taylor instabilities at the Si/O and (C+O)/He composition interfaces. This leads to the formation of high-velocity, metal-rich clumps which eventually decouple from the flow and move ballistically through the ejecta. Maximum final metal velocities of 3500-5500 km/s and 1200 km/s are obtained for the Type Ib model and the Type II model, respectively. The low maximum metal velocities in the Type II model, which are significantly smaller than those observed in SN 1987A, are due to the massive hydrogen envelope of the progenitor. The envelope forces the supernova shock to decelerate strongly, leaving behind a reverse shock below the He/H interface, which interacts with the clumps and slows them down significantly. This reverse shock is absent in the Type Ib-like model. The latter is in fairly good agreement with observations of Type Ib supernovae. In addition, in this case the pattern of clump formation in the ejecta is correlated with the convective pattern prevailing during the shock-revival phase. This might be used to deduce observational constraints for the dynamics during this early phase of the evolution, and the role of neutrino heating in initiating the explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Supernova Explosions and Neutron Star Formation

The current picture of the collapse and explosion of massive stars and the formation of neutron stars is reviewed. According to the favored scenario, however by no means proven and undisputed, neutrinos deposit the energy of the explosion in the stellar medium which surrounds the nascent neutron star. Observations, in particular of Supernova~1987A, suggest that mixing processes play an important role in the expanding star, and multi-dimensional simulations show that these are linked to convective instabilities in the immediate vicinity of the neutron star. Convectively enhanced energy transport inside the neutron star can have important consequences for the neutrino emission and thus the neutrino-heating mechanism. This also holds for a suppression of the neutrino interactions at nuclear densities. Multi-dimensional hydrodynamics, general relativity, and a better understanding of the neutrino interactions in neutron star matter may be crucial to resolve the problem that state-of-the-art spherical models do not yield explosions even with a very accurate treatment of neutrino transport by solving the Boltzmann equation.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Can neutrinos from neutron star mergers power gamma-ray bursts?

We perform three-dimensional hydrodynamical simulations of the coalescence of binary neutron stars. We include the emission and backreaction of gravitational waves into the Newtonian ``Piecewise Parabolic Method''. The use of the physical equation of state of Lattimer & Swesty (1991) allows us to take into account the production of neutrinos. We evaluate our models for the efficiency of neutrino-antineutrino annihilation in the surroundings of the coalescing neutron stars. The corresponding energy deposition prior to and during merging turns out to be 2--3 orders of magnitude too small to power a typical gamma-ray burst with an energy output of about 10^{51}/(4 pi) erg/sterad at cosmological distances. Analytical estimates of the subsequent evolution of the disk which possibly surrounds the central black hole show that even under the most favorable conditions the energy in a fireball of electron-positron pairs and photons that is created by neutrino-antineutrino annihilation falls short by at least an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

COALESCING NEUTRON STARS AS GAMMA RAY BURSTERS ?

We investigate the dynamics and evolution of coalescing neutron stars. The three-dimensional Newtonian equations of hydrodynamics are integrated by the `Piecewise Parabolic Method' However, we do include the effects of the emission of gravitational waves on the hydrodynamics. The properties of neutron star matter are described by the equation of state of Lattimer & Swesty. In addition to the fundamental hydrodynamic quantities, density, momentum, and energy, we follow the time evolution of the electron density in the stellar gas. Energy losses and changes of the electron abundance due to the emission of neutrinos are taken into account by an elaborate ``neutrino leakage scheme'', which employs a careful calculation of the lepton number and energy source terms of all neutrino types. The grid is Cartesian and equidistant with a resolution of 64**3 or 128**3, which allows us to calculate the self-gravity via fast Fourier transforms.

0 citations0 reviewsNo code linkedOpen access