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

Hans-Thomas Janka

Hans-Thomas 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-ph.HEhep-phnucl-thhep-exgr-qcnucl-exastro-ph.COastro-ph.GAastro-ph.IMphysics.comp-phphysics.ins-det

Catalog footprint

What is connected

50works
12topics
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

50 published item(s)

preprint2022arXiv

Dynamical ejecta of neutron star mergers with nucleonic weak processes II: Kilonova emission

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fast Neutrino Conversion in Hydrodynamic Simulations of Neutrino-Cooled Accretion Disks

The outflows from neutrino-cooled black-hole (BH) accretion disks formed in neutron-star mergers or cores of collapsing stars are expected to be neutron-rich enough to explain a large fraction of elements created by the rapid neutron-capture (r-) process, but their precise chemical composition remains elusive. Here, we investigate the role of fast neutrino flavor conversion, motivated by the findings of our post-processing analysis that shows evidence of electron-neutrino lepton-number (ELN) crossings deep inside the disk, hence suggesting possibly non-trivial effects due to neutrino flavor mixing. We implement a parametric, dynamically self-consistent treatment of fast conversion in time-dependent simulations and examine the impact on the disk and its outflows. By activating the, otherwise inefficient, emission of heavy-lepton neutrinos, fast conversions enhance the disk cooling rates and reduce the absorption rates of electron-type neutrinos, causing a reduction of the electron fraction in the disk by 0.03-0.06 and in the ejected material by 0.01-0.03. The r-process yields are enhanced by typically no more than a factor of two, rendering the overall impact of fast conversions modest. The kilonova is prolonged as a net result of increased lanthanide opacities and enhanced radioactive heating rates. We observe only mild sensitivity to the disk mass, the condition for the onset of flavor conversion, and to the considered cases of flavor mixing. Remarkably, parametric models of flavor mixing that conserve the lepton numbers per family result in an overall smaller impact than models invoking three-flavor equipartition, often assumed in previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Low-Energy Supernovae Severely Constrain Radiative Particle Decays

The hot and dense core formed in the collapse of a massive star is a powerful source of hypothetical feebly-interacting particles such as sterile neutrinos, dark photons, axion-like particles (ALPs), and others. Radiative decays such as $a\to2γ$ deposit this energy in the surrounding material if the mean free path is less than the radius of the progenitor star. For the first time, we use a supernova (SN) population with particularly low explosion energies as the most sensitive calorimeters to constrain this possibility. These SNe are observationally identified as low-luminosity events with low ejecta velocities and low masses of ejected $^{56}$Ni. Their low energies limit the energy deposition from particle decays to less than about 0.1 B, where $1~{\rm B~(bethe)}=10^{51}~{\rm erg}$. For 1-500 MeV-mass ALPs, this generic argument excludes ALP-photon couplings $G_{aγγ}$ in the $10^{-10}$-$10^{-8}~{\rm GeV}^{-1}$ range.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Low-luminosity type IIP supernovae: SN 2005cs and SN 2020cxd as very low-energy iron core-collapse explosions

SN 2020cxd is a representative of the family of low-energy, underluminous Type IIP supernovae (SNe), whose observations and analysis were recently reported by Yang et al. (2021). Here we re-evaluate the observational data for the diagnostic SN properties by employing the hydrodynamic explosion model of a 9 MSun red supergiant progenitor with an iron core and a pre-collapse mass of 8.75 Msun. The explosion of the star was obtained by the neutrino-driven mechanism in a fully self-consistent simulation in three dimensions (3D). Multi-band light curves and photospheric velocities for the plateau phase are computed with the one-dimensional radiation-hydrodynamics code STELLA, applied to the spherically averaged 3D explosion model as well as spherisized radial profiles in different directions of the 3D model. We find that the overall evolution of the bolometric light curve, duration of the plateau phase, and basic properties of the multi-band emission can be well reproduced by our SN model with its explosion energy of only 0.7x10^50 erg and an ejecta mass of 7.4 Msun. These values are considerably lower than the previously reported numbers, but they are compatible with those needed to explain the fundamental observational properties of the prototype low-luminosity SN 2005cs. Because of the good compatibility of our photospheric velocities with line velocities determined for SN 2005cs, we conclude that the line velocities of SN 2020cxd are probably overestimated by up to a factor of about 3. The evolution of the line velocities of SN 2005cs compared to photospheric velocities in different explosion directions might point to intrinsic asymmetries in the SN ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Pulsational pair-instability supernovae: gravitational collapse, black-hole formation, and beyond

We investigate the final collapse of rotating and non-rotating pulsational pair-instability supernova progenitors with zero-age-main-sequence masses of 60, 80, and 115$\mathrm{M}_\odot$ and iron cores between 2.37$\mathrm{M}_\odot$ and 2.72$\mathrm{M}_\odot$ by 2D hydrodynamics simulations. Using the general relativistic NADA-FLD code with energy-dependent three-flavor neutrino transport by flux-limited diffusion allows us to follow the evolution beyond the moment when the transiently forming neutron star (NS) collapses to a black hole (BH), which happens within 350$-$580 ms after bounce in all cases. Because of high neutrino luminosities and mean energies, neutrino heating leads to shock revival within $\lesssim$250 ms post bounce in all cases except the rapidly rotating 60$\mathrm{M}_\odot$ model. In the latter case, centrifugal effects support a 10% higher NS mass but reduce the radiated neutrino luminosities and mean energies by $\sim$20% and $\sim$10%, respectively, and the neutrino-heating rate by roughly a factor of two compared to the non-rotating counterpart. After BH formation, the neutrino luminosities drop steeply but continue on a 1$-$2 orders of magnitude lower level for several 100 ms because of aspherical accretion of neutrino and shock-heated matter, before the ultimately spherical collapse of the outer progenitor shells suppresses the neutrino emission to negligible values. In all shock-reviving models BH accretion swallows the entire neutrino-heated matter and the explosion energies decrease from maxima around 1.5$\times$10$^{51}$erg to zero within a few seconds latest. Nevertheless, the shock or a sonic pulse moves outward and may trigger mass loss, which we estimate by long-time simulations with the PROMETHEUS code. We also provide gravitational-wave signals.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Synthetic observables for electron-capture supernovae and low-mass core collapse supernovae

Stars in the mass range from 8 to 10 solar masses are expected to produce one of two types of supernovae (SNe), either electron-capture supernovae (ECSNe) or core-collapse supernovae (CCSNe), depending on their previous evolution. Either of the associated progenitors retain extended and massive hydrogen-rich envelopes, the observables of these SNe are, therefore, expected to be similar. In this study we explore the differences in these two types of SNe. Specifically, we investigate three different progenitor models: a solar-metallicity ECSN progenitor with an initial mass of 8.8 solar masses, a zero-metallicity progenitor with 9.6 solar masses, and a solar-metallicity progenitor with 9 solar masses, carrying out radiative transfer simulations for these progenitors. We present the resulting light curves for these models. The models exhibit very low photospheric velocity variations of about 2000 km/s, therefore, this may serve as a convenient indicator of low-mass SNe. The ECSN has very unique light curves in broad bands, especially the U band, and does not resemble any currently observed SN. This ECSN progenitor being part of a binary will lose its envelope for which reason the light curve becomes short and undetectable. The SN from the 9.6 solar masses progenitor exhibits also quite an unusual light curve, explained by the absence of metals in the initial composition. The artificially iron polluted 9.6 solar masses model demonstrates light curves closer to normal SNe IIP. The SN from the 9 solar masses progenitor remains the best candidate for so-called low-luminosity SNe IIP like SN 1999br and SN 2005cs.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The infancy of core-collapse supernova remnants

We present 3D hydrodynamic simulations of neutrino-driven supernovae (SNe) with the PROMETHEUS-HOTB code, evolving the asymmetrically expanding ejecta from shock breakout until they reach the homologous expansion phase after roughly one year. Our calculations continue the simulations for two red supergiant (RSG) and two blue supergiant (BSG) progenitors by Wongwathanarat et al., who investigated the growth of explosion asymmetries produced by hydrodynamic instabilities during the first second of the explosion and their later fragmentation by Rayleigh-Taylor instabilities. We focus on the late time acceleration and inflation of the ejecta caused by the heating due to the radioactive decay of $^{56}$Ni to $^{56}$Fe and by a new outward-moving shock, which forms when the reverse shock from the He/H-shell interface compresses the central part of the ejecta. The mean velocities of the iron-rich ejecta increase between 100 km/s and 350 km/s ($\sim$8-30\%), and the fastest one percent of the iron accelerates by up to $\sim$1000 km/s ($\sim$20-25\%). This 'Ni-bubble effect', known from 1D models, accelerates the bulk of the nickel in our 3D models and causes an inflation of the initially overdense Ni-rich clumps, which leads to underdense, extended fingers, enveloped by overdense skins of compressed surrounding matter. We also provide volume and surface filling factors as well as a spherical harmonics analysis to characterize the spectrum of Ni-clump sizes quantitatively. Three of our four models give volume filling factors larger than 0.3, consistent with what is suggested for SN 1987A by observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetar formation through a convective dynamo in protoneutron stars

The release of spin-down energy by a magnetar is a promising scenario to power several classes of extreme explosive transients. However, it lacks a firm basis because magnetar formation still represents a theoretical challenge. Using the first three-dimensional simulations of a convective dynamo based on a protoneutron star interior model, we demonstrate that the required dipolar magnetic field can be consistently generated for sufficiently fast rotation rates. The dynamo instability saturates in the magnetostrophic regime with the magnetic energy exceeding the kinetic energy by a factor of up to 10. Our results are compatible with the observational constraints on galactic magnetar field strength and provide strong theoretical support for millisecond protomagnetar models of gamma-ray burst and superluminous supernova central engines.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutrino emission characteristics of black hole formation in three-dimensional simulations of stellar collapse

Neutrinos are unique probes of core-collapse supernova dynamics, especially in the case of black hole (BH) forming stellar collapses, where the electromagnetic emission may be faint or absent. By investigating two 3D hydrodynamical simulations of BH-forming stellar collapses of mass 40 and 75 M_sun, we identify the physical processes preceding BH formation through neutrinos, and forecast the neutrino signal expected in the existing IceCube and Super-Kamiokande detectors, as well as in the future generation DUNE facility. Prior to the abrupt termination of the neutrino signal corresponding to BH formation, both models develop episodes of strong and long-lasting activity by the spiral standing accretion shock instability (SASI). We find that the spiral SASI peak in the Fourier power spectrum of the neutrino event rate will be distinguishable at 3 sigma above the detector noise for distances up to O(30) kpc in the most optimistic scenario, with IceCube having the highest sensitivity. Interestingly, given the long duration of the spiral SASI episodes, the spectrograms of the expected neutrino event rate carry clear signs of the evolution of the blue spiral SASI frequency as a function of time, as the shock radius and post-shock fluid velocity evolve. Due to the high accretion luminosity and its large-amplitude SASI-induced modulations, any contribution from asymmetric (dipolar or quadrupolar) neutrino emission associated with the lepton emission self-sustained asymmetry (LESA) is far subdominant in the neutrino signal.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

NS 1987A in SN 1987A

The possible detection of a compact object in the remnant of SN 1987A presents an unprecedented opportunity to follow its early evolution. The suspected detection stems from an excess of infrared emission from a dust blob near the compact object's predicted position. The infrared excess could be due to the decay of isotopes like 44Ti, accretion luminosity from a neutron star or black hole, magnetospheric emission or a wind originating from the spindown of a pulsar, or thermal emission from an embedded, cooling neutron star (NS 1987A). It is shown that the last possibility is the most plausible as the other explanations are disfavored by other observations and/or require fine-tuning of parameters. Not only are there indications the dust blob overlaps the predicted location of a kicked compact remnant, but its excess luminosity also matches the expected thermal power of a 30 year old neutron star. Furthermore, models of cooling neutron stars within the Minimal Cooling paradigm readily fit both NS 1987A and Cas A, the next-youngest known neutron star. If correct, a long heat transport timescale in the crust and a large effective stellar temperature are favored, implying relatively limited crustal n-1S0 superfluidity and an envelope with a thick layer of light elements, respectively. If the locations don't overlap, then pulsar spindown or accretion might be more likely, but the pulsar's period and magnetic field or the accretion rate must be rather finely tuned. In this case, NS 1987A may have enhanced cooling and/or a heavy-element envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Impact of Nucleon-Nucleon Bremsstrahlung Rates Beyond One-Pion Exchange

Neutrino-pair production and annihilation through nucleon-nucleon bremsstrahlung is included in current supernova simulations by rates that are based on the one-pion-exchange approximation. Here we explore the consequences of bremsstrahlung rates based on a modern nuclear interactions for proto-neutron star cooling and the corresponding neutrino emission. We find that despite a reduction of the bremsstrahlung emission by a factor of 2-5 in the neutrinospheric region, models with the improved treatment exhibit only $\lesssim$5% changes of the neutrino luminosities and an increase of $\lesssim$0.7 MeV of the average energies of the radiated neutrino spectra, with the largest effects for the antineutrinos of all flavors and at late times. Overall, the proto-neutron star cooling evolution is slowed down modestly by $\lesssim$0.5-1 s.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Progenitor-dependent Explosion Dynamics in Self-consistent, Axisymmetric Simulations of Neutrino-driven Core-collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11-28 solar masses, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si-O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si-O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection time scales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supernova Neutrinos: Production, Oscillations and Detection

Neutrinos play a crucial role in the collapse and explosion of massive stars, governing the infall dynamics of the stellar core, triggering and fueling the explosion and driving the cooling and deleptonization of the newly formed neutron star. Due to their role neutrinos carry information from the heart of the explosion and, due to their weakly interacting nature, offer the only direct probe of the dynamics and thermodynamics at the center of a supernova. In this paper, we review the present status of modelling the neutrino physics and signal formation in collapsing and exploding stars. We assess the capability of current and planned large underground neutrino detectors to yield faithful information of the time and flavor dependent neutrino signal from a future Galactic supernova. We show how the observable neutrino burst would provide a benchmark for fundamental supernova physics with unprecedented richness of detail. Exploiting the treasure of the measured neutrino events requires a careful discrimination of source-generated properties from signal features that originate on the way to the detector. As for the latter, we discuss self-induced flavor conversions associated with neutrino-neutrino interactions that occur in the deepest stellar regions; matter effects that modify the pattern of flavor conversions in the dynamical stellar envelope; neutrino-oscillation signatures that result from structural features associated with the shock-wave propagation as well as turbulent mass motions in post-shock layers. Finally, we highlight our current understanding of the formation of the diffuse supernova neutrino background and we analyse the perspectives for a detection of this relic signal that integrates the contributions from all past core-collapse supernovae in the Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Last Minutes of Oxygen Shell Burning in a Massive Star

We present the first 3D simulation of the last minutes of oxygen shell burning in an 18 solar mass supernova progenitor up to the onset of core collapse. A moving inner boundary is used to accurately model the contraction of the silicon and iron core according to a 1D stellar evolution model with a self-consistent treatment of core deleptonization and nuclear quasi-equilibrium. The simulation covers the full solid angle to allow the emergence of large-scale convective modes. Due to core contraction and the concomitant acceleration of nuclear burning, the convective Mach number increases to ~0.1 at collapse, and an l=2 mode emerges shortly before the end of the simulation. Aside from a growth of the oxygen shell from 0.51 to 0.56 solar masses due to entrainment from the carbon shell, the convective flow is reasonably well described by mixing length theory, and the dominant scales are compatible with estimates from linear stability analysis. We deduce that artificial changes in the physics, such as accelerated core contraction, can have precarious consequences for the state of convection at collapse. We argue that scaling laws for the convective velocities and eddy sizes furnish good estimates for the state of shell convection at collapse and develop a simple analytic theory for the impact of convective seed perturbations on shock revival in the ensuing supernova. We predict a reduction of the critical luminosity for explosion by 12--24% due to seed asphericities for our 3D progenitor model relative to the case without large seed perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Exploring properties of high-density matter through remnants of neutron-star mergers

Remnants of neutron-star mergers are essentially massive, hot, differentially rotating neutron stars, which are initially strongly oscillating. They represent a unique probe for high-density matter because the oscillations are detectable via gravitational-wave measurements and are strongly dependent on the equation of state. The impact of the equation of state is apparent in the frequency of the dominant oscillation mode of the remnant. For a fixed total binary mass a tight relation between the dominant postmerger frequency and the radii of nonrotating neutron stars exists. Inferring observationally the dominant postmerger frequency thus determines neutron star radii with high accuracy of the order of a few hundred meters. By considering symmetric and asymmetric binaries of the same chirp mass, we show that the knowledge of the binary mass ratio is not critical for this kind of radius measurements. We summarize different possibilities to deduce the maximum mass of nonrotating neutron stars. We clarify the nature of the three most prominent features of the postmerger gravitational-wave spectrum and argue that the merger remnant can be considered to be a single, isolated, self-gravitating object that can be described by concepts of asteroseismology. The understanding of the different mechanisms shaping the gravitational-wave signal yields a physically motivated analytic model of the gravitational-wave emission, which may form the basis for template-based gravitational-wave data analysis. We explore the observational consequences of a scenario of two families of compact stars including hadronic and quark stars. We find that this scenario leaves a distinctive imprint on the postmerger gravitational-wave signal. In particular, a strong discontinuity in the dominant postmerger frequency as function of the total mass will be a strong indication for two families of compact stars. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impact of weak interactions of free nucleons on the r-process in dynamical ejecta from neutron-star mergers

We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino viscosity and drag: impact on the magnetorotational instability in protoneutron stars

The magnetorotational instability (MRI) is a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. The diffusion of neutrinos trapped in the PNS induces a transport of momentum, which can be modelled as a viscosity on length-scales longer than the neutrino mean free path. This neutrino-viscosity can slow down the growth of MRI modes to such an extent that a minimum initial magnetic field strength of $\gtrsim 10^{12}\, {\rm G}$ is needed for the MRI to grow on a sufficiently short time-scale to potentially affect the explosion. It is uncertain whether the magnetic field of fast rotating progenitor cores is strong enough to yield such an initial magnetic field in PNS. At MRI wavelengths shorter than the neutrino mean free path, on the other hand, neutrino radiation does not act as a viscosity but rather induces a drag on the velocity with a damping rate independent of the wavelength. We perform a linear analysis of the MRI in this regime, and apply our analytical results to the PNS structure from a one-dimensional numerical simulation. We show that in the outer layers of the PNS, the MRI can grow from weak magnetic fields at wavelengths shorter than the neutrino mean free path, while deeper in the PNS MRI growth takes place in the viscous regime and requires a minimum magnetic field strength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino-driven explosion of a 20 solar-mass star in three dimensions enabled by strange-quark contributions to neutrino-nucleon scattering

Interactions with neutrons and protons play a crucial role for the neutrino opacity of matter in the supernova core. Their current implementation in many simulation codes, however, is rather schematic and ignores not only modifications for the correlated nuclear medium of the nascent neutron star, but also free-space corrections from nucleon recoil, weak magnetism or strange quarks, which can easily add up to changes of several 10% for neutrino energies in the spectral peak. In the Garching supernova simulations with the Prometheus-Vertex code, such sophistications have been included for a long time except for the strange-quark contributions to the nucleon spin, which affect neutral-current neutrino scattering. We demonstrate on the basis of a 20 M_sun progenitor star that a moderate strangeness-dependent contribution of g_a^s = -0.2 to the axial-vector coupling constant g_a = 1.26 can turn an unsuccessful three-dimensional (3D) model into a successful explosion. Such a modification is in the direction of current experimental results and reduces the neutral-current scattering opacity of neutrons, which dominate in the medium around and above the neutrinosphere. This leads to increased luminosities and mean energies of all neutrino species and strengthens the neutrino-energy deposition in the heating layer. Higher nonradial kinetic energy in the gain layer signals enhanced buoyancy activity that enables the onset of the explosion at ~300 ms after bounce, in contrast to the model with vanishing strangeness contributions to neutrino-nucleon scattering. Our results demonstrate the close proximity to explosion of the previously published, unsuccessful 3D models of the Garching group.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino-driven supernova of a low-mass iron-core progenitor boosted by three-dimensional turbulent convection

We present the first successful simulation of a neutrino-driven supernova explosion in three dimensions (3D), using the Prometheus-Vertex code with an axis-free Yin-Yang grid and a sophisticated treatment of three-flavor, energy-dependent neutrino transport. The progenitor is a nonrotating, zero-metallicity 9.6 Msun star with an iron core. While in spherical symmetry outward shock acceleration sets in later than 300 ms after bounce, a successful explosion starts at ~130 ms postbounce in two dimensions (2D). The 3D model explodes at about the same time but with faster shock expansion than in 2D and a more quickly increasing and roughly 10 percent higher explosion energy of >10^50 erg. The more favorable explosion conditions in 3D are explained by lower temperatures and thus reduced neutrino emission in the cooling layer below the gain radius. This moves the gain radius inward and leads to a bigger mass in the gain layer, whose larger recombination energy boosts the explosion energy in 3D. These differences are caused by less coherent, less massive, and less rapid convective downdrafts associated with postshock convection in 3D. The less violent impact of these accretion downflows in the cooling layer produces less shock heating and therefore diminishes energy losses by neutrino emission. We thus have, for the first time, identified a reduced mass accretion rate, lower infall velocities, and a smaller surface filling factor of convective downdrafts as consequences of 3D postshock turbulence that facilitate neutrino-driven explosions and strengthen them compared to the 2D case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the nuclear robustness of the r process in neutron-star mergers

We have performed r-process calculations for matter ejected dynamically in neutron star mergers based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic simulation. Our calculations consider an extended nuclear network, including spontaneous, $β$- and neutron-induced fission and adopting fission yield distributions from the ABLA code. We have studied the sensitivity of the r-process abundances to nuclear masses by using different models. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the abundance yields at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number $N=130$. The remaining trajectories, which contribute 10% by mass to the total integrated abundances, follow such a fast expansion that the r process does not use all the neutrons. This also leads to a larger variation of abundances among trajectories as fission does not dominate the r-process dynamics. The total integrated abundances are dominated by contributions from the slow abundances and hence reproduce the general features of the observed r-process abundances. We find that at timescales of weeks relevant for kilonova light curve calculations, the abundance of actinides is larger than the one of lanthanides. Hence actinides can be even more important than lanthanides to determine the photon opacities under kilonova conditions. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Supernova deleptonization asymmetry: Impact on self-induced flavor conversion

During the accretion phase of a core-collapse supernova (SN), the deleptonization flux has recently been found to develop a global dipole pattern (LESA---Lepton Emission Self-sustained Asymmetry). The $ν_e$ minus $\barν_e$ flux essentially vanishes in one direction, potentially facilitating self-induced flavor conversion. On the other hand, below the stalled shock wave, self-induced flavor conversion is typically suppressed by multi-angle matter effects, preventing any impact of flavor conversion on SN explosion dynamics. In a schematic model of SN neutrino fluxes, we study the impact of modified $\barν_e$-$ν_e$ flux asymmetries on collective flavor conversion. In the parameter space consisting of matter density and effective neutrino density, the region of instability with regard to self-induced flavor conversion is much larger for a vanishing lepton number flux, yet this modification does not intersect a realistic SN profile. Therefore, it appears that, even in the presence of LESA, self-induced flavor conversion remains suppressed below the shock front.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterizing SASI- and Convection-Dominated Core-Collapse Supernova Explosions in Two Dimensions

The success of the neutrino mechanism of core-collapse supernovae relies on the supporting action of two hydrodynamic instabilities: neutrino-driven convection and the Standing Accretion Shock Instability (SASI). Depending on the structure of the stellar progenitor, each of these instabilities can dominate the evolution of the gain region prior to the onset of explosion, with implications for the ensuing asymmetries. Here we examine the flow dynamics in the neighborhood of explosion by means of parametric two-dimensional, time-dependent hydrodynamic simulations for which the linear stability properties are well understood. We find that systems for which the convection parameter is sub-critical (SASI-dominated) develop explosions once large-scale, high-entropy bubbles are able to survive for several SASI oscillation cycles. These long-lived structures are seeded by the SASI during shock expansions. Finite-amplitude initial perturbations do not alter this outcome qualitatively, though they can lead to significant differences in explosion times. Supercritical systems (convection-dominated) also explode by developing large-scale bubbles, though the formation of these structures is due to buoyant activity. Non-exploding systems achieve a quasi-steady state in which the time-averaged flow adjusts itself to be convectively sub-critical. We characterize the turbulent flow using a spherical Fourier-Bessel decomposition, identifying the relevant scalings and connecting temporal and spatial components. Finally, we verify the applicability of these principles on the general relativistic, radiation-hydrodynamic simulations of Mueller, Janka, & Heger (2012), and discuss implications for the three-dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutrino emission characteristics and detection opportunities based on three-dimensional supernova simulations

The neutrino emission characteristics of the first full-scale three-dimensional supernova simulations with sophisticated three-flavor neutrino transport for three models with masses 11.2, 20 and 27 M_sun are evaluated in detail. All the studied progenitors show the expected hydrodynamical instabilities in the form of large-scale convective overturn. In addition, the recently identified LESA phenomenon (lepton-number emission self-sustained asymmetry) is generic for all our cases. Pronounced SASI (standing accretion-shock instability) activity appears in the 20 and 27 M_sun cases, partly in the form of a spiral mode, inducing large but direction and flavor-dependent modulations of neutrino emission. These modulations can be clearly identified in the existing IceCube and future Hyper-Kamiokande detectors, depending on distance and detector location relative to the main SASI sloshing direction.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self-sustained asymmetry of lepton-number emission: A new phenomenon during the supernova shock-accretion phase in three dimensions

During the stalled-shock phase of our 3D hydrodynamical core-collapse simulations with energy-dependent, 3-flavor neutrino transport, the lepton-number flux (nue minus antinue) emerges predominantly in one hemisphere. This novel, spherical-symmetry breaking neutrino-hydrodynamical instability is termed LESA for "Lepton-number Emission Self-sustained Asymmetry." While the individual nue and antinue fluxes show a pronounced dipole pattern, the heavy-flavor neutrino fluxes and the overall luminosity are almost spherically symmetric. Initially, LESA seems to develop stochastically from convective fluctuations, it exists for hundreds of milliseconds or more, and it persists during violent shock sloshing associated with the standing accretion shock instability. The nue minus antinue flux asymmetry originates mainly below the neutrinosphere in a region of pronounced proto-neutron star (PNS) convection, which is stronger in the hemisphere of enhanced lepton-number flux. On this side of the PNS, the mass-accretion rate of lepton-rich matter is larger, amplifying the lepton-emission asymmetry, because the spherical stellar infall deflects on a dipolar deformation of the stalled shock. The increased shock radius in the hemisphere of less mass accretion and minimal lepton-number flux (antinue flux maximum) is sustained by stronger convection on this side, which is boosted by stronger neutrino heating because the average antinue energy is higher than the average nue energy. Asymmetric heating thus supports the global deformation despite extremely nonstationary convective overturn behind the shock. While these different elements of LESA form a consistent picture, a full understanding remains elusive at present. There may be important implications for neutrino-flavor oscillations, the neutron-to-proton ratio in the neutrino-heated supernova ejecta, and neutron-star kicks, which remain to be explored.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Towards Petaflops Capability of the VERTEX Supernova Code

The VERTEX code is employed for multi-dimensional neutrino-radiation hydrodynamics simulations of core-collapse supernova explosions from first principles. The code is considered state-of-the-art in supernova research and it has been used for modeling for more than a decade, resulting in numerous scientific publications. The computational performance of the code, which is currently deployed on several high-performance computing (HPC) systems up to the Tier-0 class (e.g. in the framework of the European PRACE initiative and the German GAUSS program), however, has so far not been extensively documented. This paper presents a high-level overview of the relevant algorithms and parallelization strategies and outlines the technical challenges and achievements encountered along the evolution of the code from the gigaflops scale with the first, serial simulations in 2000, up to almost petaflops capabilities, as demonstrated lately on the SuperMUC system of the Leibniz Supercomputing Centre (LRZ). In particular, we shall document the parallel scalability and computational efficiency of VERTEX at the large scale and on the major, contemporary HPC platforms. We will outline upcoming scientific requirements and discuss the resulting challenges for the future development and operation of the code.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Turbulence patterns and neutrino flavor transitions in high-resolution supernova models

During the shock-wave propagation in a core-collapse supernova (SN), matter turbulence may affect neutrino flavor conversion probabilities. Such effects have been usually studied by adding parametrized small-scale random fluctuations (with arbitrary amplitude) on top of coarse, spherically symmetric matter density profiles. Recently, however, two-dimensional (2D) SN models have reached a space resolution high enough to directly trace anisotropic density profiles, down to scales smaller than the typical neutrino oscillation length. In this context, we analyze the statistical properties of a large set of SN matter density profiles obtained in a high-resolution 2D simulation, focusing on a post-bounce time (2 s) suited to study shock-wave effects on neutrino propagation on scales as small as O(100) km and possibly below. We clearly find the imprint of a broken (Kolmogorov-Kraichnan) power-law structure, as generically expected in 2D turbulence spectra. We then compute the flavor evolution of SN neutrinos along representative realizations of the turbulent matter density profiles, and observe no or modest damping of the neutrino crossing probabilities on their way through the shock wave. In order to check the effect of possibly unresolved fluctuations at scales below O(100) km, we also apply a randomization procedure anchored to the power spectrum calculated from the simulation, and find consistent results within \pm 1 sigma fluctuations. These results show the importance of anchoring turbulence effects on SN neutrinos to realistic, fine-grained SN models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Multi-Dimensional General Relativistic Neutrino Hydrodynamics Code of Core-Collapse Supernovae III. Gravitational Wave Signals from Supernova Explosion Models

We present a detailed theoretical analysis of the gravitational-wave (GW) signal of the post-bounce evolution of core-collapse supernovae (SNe), employing for the first time relativistic, two-dimensional (2D) explosion models with multi-group, three-flavor neutrino transport based on the ray-by-ray-plus approximation. The waveforms reflect the accelerated mass motions associated with the characteristic evolutionary stages that were also identified in previous works: A quasi-periodic modulation by prompt postshock convection is followed by a phase of relative quiescence before growing amplitudes signal violent hydrodynamical activity due to convection and the standing accretion shock instability during the accretion period of the stalled shock. Finally, a high-frequency, low-amplitude variation from proto-neutron star (PNS) convection below the neutrinosphere appears superimposed on the low-frequency trend associated with the aspherical expansion of the SN shock after the onset of the explosion. Relativistic effects in combination with detailed neutrino transport are shown to be essential for quantitative predictions of the GW frequency evolution and energy spectrum, because they determine the structure of the PNS surface layer and its characteristic g-mode frequency. Burst-like high-frequency activity phases, correlated with sudden luminosity increase and spectral hardening of electron (anti-)neutrino emission for some 10ms, are discovered as new features after the onset of the explosion. They correspond to intermittent episodes of anisotropic accretion by the PNS in the case of fallback SNe. We find stronger signals for more massive progenitors with large accretion rates. The typical frequencies are higher for massive PNSs, though the time-integrated spectrum also strongly depends on the model dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron-capture supernovae as origin of 48Ca

We report that electron-capture supernovae (ECSNe), arising from collapsing oxygen-neon-magnesium cores, are a possible source of 48Ca, whose origin has remained a long-standing puzzle. Our two-dimensional, self-consistent explosion model of an ECSN predicts ejection of neutron-rich matter with electron fractions Ye = 0.40-0.42 and relatively low entropies, s = 13-15 kB per nucleon (kB is the Boltzmann constant). Post-processing nucleosynthesis calculations result in appreciable production of 48Ca in such neutron-rich and low-entropy matter during the quasi-nuclear equilibrium and subsequent freezeout phases. The amount of ejected 48Ca can account for that in the solar inventory when we consider possible uncertainties in the entropies or ejecta-mass distribution. ECSNe could thus be a site of 48Ca production in addition to a hypothetical, rare class of high-density Type Ia supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron-capture supernovae as sources of 60Fe

We investigate the nucleosynthesis of the radionuclide 60Fe in electron-capture supernovae (ECSNe). The nucleosynthetic results are based on a self-consistent, two-dimensional simulation of an ECSN as well as models in which the densities are systematically increased by some factors (low-entropy models). 60Fe is found to be appreciably made in neutron-rich ejecta during the nuclear quasi-equilibrium phase with greater amounts being produced in the lower-entropy models. Our results, combining them with the yields of core-collapse supernovae (CCSNe) in the literature, suggest that ECSNe account for at least 4-30% of live 60Fe in the Milky Way. ECSNe co-produce neutron-rich isotopes, 48Ca, 50Ti, 54Cr, some light trans-iron elements, and possibly weak r-process elements including some radionuclides such as 93Zr, 99Tc, and 107Pd, whose association with 60Fe might have been imprinted in primitive meteorites or in the deep ocean crust on the Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High-resolution supernova neutrino spectra represented by a simple fit

To study the capabilities of supernova neutrino detectors, the instantaneous spectra are often represented by a quasi-thermal distribution of the form f(E) = E^alpha e^{-(alpha+1)E/E_{av}} where E_{av} is the average energy and alpha a numerical parameter. Based on a spherically symmetric supernova model with full Boltzmann neutrino transport we have, at a few representative post-bounce times, re-converged the models with vastly increased energy resolution to test the fit quality. For our examples, the spectra are well represented by such a fit in the sense that the counting rates for a broad range of target nuclei, sensitive to different parts of the spectrum, are reproduced very well. Therefore, the mean energy and root-mean-square energy of numerical spectra hold enough information to provide the correct alpha and to forecast the response of multi-channel supernova neutrino detection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutrino signature of supernova hydrodynamical instabilities in three dimensions

The first full-scale three-dimensional (3D) core-collapse supernova (SN) simulations with sophisticated neutrino transport show pronounced effects of the standing accretion shock instability (SASI) for two high-mass progenitors (20 and 27 M_sun). In a low-mass progenitor (11.2 M_sun), large-scale convection is the dominant nonradial hydrodynamic instability in the postshock accretion layer. The SASI-associated modulation of the neutrino signal (80 Hz in our two examples) will be clearly detectable in IceCube or the future Hyper-Kamiokande detector, depending on progenitor properties, distance, and observer location relative to the main SASI sloshing direction. The neutrino signal from the next galactic SN can therefore diagnose the nature of the hydrodynamic instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fast time variations of supernova neutrino signals from 3-dimensional models

We study supernova neutrino flux variations in the IceCube detector, using 3D models based on a simplified neutrino transport scheme. The hemispherically integrated neutrino emission shows significantly smaller variations compared with our previous study of 2D models, largely because of the reduced SASI activity in this set of 3D models which we interpret as a pessimistic extreme. For the studied cases, intrinsic flux variations up to about 100 Hz frequencies could still be detected in a supernova closer than about 2 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Impact of eV-mass sterile neutrinos on neutrino-driven supernova outflows

Motivated by recent hints for sterile neutrinos from the reactor anomaly, we study active-sterile conversions in a three-flavor scenario (2 active + 1 sterile families) for three different representative times during the neutrino-cooling evolution of the proto-neutron star born in an electron-capture supernova. In our "early model" (0.5 s post bounce), the nu_e-nu_s MSW effect driven by Delta m^2=2.35 eV^2 is dominated by ordinary matter and leads to a complete nu_e-nu_s swap with little or no trace of collective flavor oscillations. In our "intermediate" (2.9 s p.b.) and "late models" (6.5 s p.b.), neutrinos themselves significantly modify the nu_e-nu_s matter effect, and, in particular in the late model, nu-nu refraction strongly reduces the matter effect, largely suppressing the overall nu_e-nu_s MSW conversion. This phenomenon has not been reported in previous studies of active-sterile supernova neutrino oscillations. We always include the feedback effect on the electron fraction Y_e due to neutrino oscillations. In all examples, Y_e is reduced and therefore the presence of sterile neutrinos can affect the conditions for heavy-element formation in the supernova ejecta, even if probably not enabling the r-process in the investigated outflows of an electron-capture supernova. The impact of neutrino-neutrino refraction is strong but complicated, leaving open the possibility that with a more complete treatment, or for other supernova models, active-sterile neutrino oscillations could generate conditions suitable for the r-process.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Is Strong SASI Activity the Key to Successful Neutrino-Driven Supernova Explosions?

Following a simulation approach of recent publications we explore the viability of the neutrino-heating explosion mechanism in dependence on the spatial dimension. Our results disagree with previous findings. While we also observe that two-dimensional (2D) models explode for lower driving neutrino luminosity than spherically symmetric (1D) models, we do not find that explosions in 3D occur easier and earlier than in 2D. Moreover, we find that the average entropy of matter in the gain layer hardly depends on the dimension and thus is no good diagnostic quantity for the readiness to explode. Instead, mass, integrated entropy, total neutrino-heating rate, and nonradial kinetic energy in the gain layer are higher when models are closer to explosion. Coherent, large-scale mass motions as typically associated with the standing accretion-shock instability (SASI) are observed to be supportive for explosions because they drive strong shock expansion and thus enlarge the gain layer. While 2D models with better angular resolution explode clearly more easily, the opposite trend is seen in 3D. We interpret this as a consequence of the turbulent energy cascade, which transports energy from small to large spatial scales in 2D, thus fostering SASI activity. In contrast, the energy flow in 3D is in the opposite direction, feeding fragmentation and vortex motions on smaller scales and thus making the 3D evolution with finer grid resolution more similar to 1D. More favorable conditions for explosions in 3D may therefore be tightly linked to efficient growth of low-order SASI modes including nonaxisymmetric ones.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing Lorentz Violation in Neutrino Propagation from a Core-Collapse Supernova

Supernova explosions provide the most sensitive probes of neutrino propagation, such as the possibility that neutrino velocities might be affected by the foamy structure of space-time thought to be generated by quantum-gravitational (QG) effects. Recent two-dimensional simulations of the neutrino emissions from core-collapse supernovae suggest that they might exhibit variations in time on the scale of a few milliseconds. We analyze simulations of such neutrino emissions using a wavelet technique, and consider the limits that might be set on a linear or quadratic violation of Lorentz invariance in the group velocities of neutrinos of different energies, v/c = [1 \pm (E/M_{nuLV1})] or [1 \pm (E/M_{\nuLV2})^2], if variations on such short time scales were to be observed, where the mass scales M_{nuLVi} might appear in models of quantum gravity. We find prospective sensitivities to M_{nuLV1} ~ 2 X 10^{13} GeV and M_{nuLV2} ~ 10^6 GeV at the 95% confidence level, up to two orders of magnitude beyond estimates made using previous one-dimensional simulations of core-collapse supernovae. We also analyze the prospective sensitivities to scenarios in which the propagation times of neutrinos of fixed energies are subject to stochastic fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing the neutrino mass hierarchy with the rise time of a supernova burst

The rise time of a Galactic supernova (SN) bar-nue lightcurve, observable at a high-statistics experiment such as the IceCube Cherenkov detector, can provide a diagnostic tool for the neutrino mass hierarchy at "large" 1-3 leptonic mixing angle theta_13. Thanks to the combination of matter suppression of collective effects at early postbounce times on one hand and the presence of the ordinary Mikheyev-Smirnov-Wolfenstein effect in the outer layers of the SN on the other hand, a sufficiently fast rise time on O(100) ms scale is indicative of an inverted mass hierarchy. We investigate results from an extensive set of stellar core-collapse simulations, providing a first exploration of the astrophysical robustness of these features. We find that for all the models analyzed (sharing the same weak interaction microphysics) the rise times for the same hierarchy are similar not only qualitatively, but also quantitatively, with the signals for the two classes of hierarchies significantly separated. We show via Monte Carlo simulations that the two cases should be distinguishable at IceCube for SNe at a typical Galactic distance 99% of the times. Finally, a preliminary survey seems to show that the faster rise time for inverted hierarchy as compared to normal hierarchy is a qualitatively robust feature predicted by several simulation groups. Since the viability of this signature ultimately depends on the quantitative assessment of theoretical/numerical uncertainties, our results motivate an extensive campaign of comparison of different code predictions at early accretion times with implementation of microphysics of comparable sophistication, including effects such like nucleon recoils in weak interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Prospective Constraints on Neutrino Masses from a Core-Collapse Supernova

We discuss the prospects for improved upper limits on neutrino masses that may be provided by a core-collapse supernova explosion in our galaxy, if it exhibits time variations in the neutrino emissions on the scale of a few milliseconds as suggested by recent two-dimensional simulations. Analyzing simulations of such neutrino emissions using the wavelet technique adopted in [1], we find that an upper limit m_nu ~ 0.14 eV could be established at the 95% confidence level if the time variations in emissions were to be preserved during neutrino propagation to the Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Relativistic collapse and explosion of rotating supermassive stars with thermonuclear effects

We present results of general relativistic simulations of collapsing supermassive stars with and without rotation using the two-dimensional general relativistic numerical code Nada, which solves the Einstein equations written in the BSSN formalism and the general relativistic hydrodynamics equations with high resolution shock capturing schemes. These numerical simulations use an equation of state which includes effects of gas pressure, and in a tabulated form those associated with radiation and the electron-positron pairs. We also take into account the effect of thermonuclear energy released by hydrogen and helium burning. We find that objects with a mass of 5x10^{5} solar mass and an initial metallicity greater than Z_{CNO}~0.007 do explode if non-rotating, while the threshold metallicity for an explosion is reduced to Z_{CNO}~0.001 for objects uniformly rotating. The critical initial metallicity for a thermonuclear explosion increases for stars with mass ~10^{6} solar mass. For those stars that do not explode we follow the evolution beyond the phase of black hole formation. We compute the neutrino energy loss rates due to several processes that may be relevant during the gravitational collapse of these objects. The peak luminosities of neutrinos and antineutrinos of all flavors for models collapsing to a BH are ~10^{55} erg/s. The total radiated energy in neutrinos varies between ~10^{56} ergs for models collapsing to a BH, and ~10^{45}-10^{46} ergs for models exploding.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Supernova neutrino halo and the suppression of self-induced flavor conversion

Neutrinos streaming from a supernova (SN) core occasionally scatter in the envelope, producing a small "neutrino halo" with a much broader angle distribution than the primary flux originating directly from the core. Cherry et al. (2012) have recently pointed out that, during the accretion phase, the halo actually dominates neutrino-neutrino refraction at distances exceeding some 100 km. However, the multiangle matter effect (which increases if the angle distribution is broader) still appears to suppress self-induced flavor conversion during the accretion phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Suppression of Self-Induced Flavor Conversion in the Supernova Accretion Phase

Self-induced flavor conversions of supernova (SN) neutrinos can strongly modify the flavor dependent fluxes. We perform a linearized flavor stability analysis with accretion-phase matter profiles of a 15 M_sun spherically symmetric model and corresponding neutrino fluxes. We use realistic energy and angle distributions, the latter deviating strongly from quasi-isotropic emission, thus accounting for both multi-angle and multi-energy effects. For our matter and neutrino density profile we always find stable conditions: flavor conversions are limited to the usual MSW effect. In this case one may distinguish the neutrino mass hierarchy in a SN neutrino signal if the mixing angle theta_13 is as large as suggested by recent experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The next-generation liquid-scintillator neutrino observatory LENA

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Core-Collapse Supernovae: Explosion Dynamics, Neutrinos and Gravitational Waves

The quest for the supernova explosion mechanism has been one of the outstanding challenges in computational astrophysics for several decades. Simulations have now progressed to a stage at which the solution appears close and neutrino and gravitational wave signals from self-consistent explosion models are becoming available. Here we focus one of the recent advances in supernova modeling, the inclusion of general relativity in multi-dimensional neutrino hydrodynamics simulations, and present the latest simulation results for an 11.2 and a 15 solar mass progenitor. We also mention 3D effects as another aspect in supernova physics awaiting further, more thorough investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

General Relativistic Explosion Models of Core-Collapse Supernovae

We present results from the first generation of multi-dimensional general relativistic neutrino hydrodynamics simulations of core-collapse supernovae. A comparison with models computed using either the purely Newtonian approximation or the "effective gravitational potential" approach reveals appreciable quantitative differences in the heating conditions and the gravitational wave spectra. Our results underscore the important role of general relativity in the supernova problem (which appears to be on par with other important factors such as the dimensionality and the equation of state) both for our understanding of the explosion dynamics as well as for predictions of observable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational waves and nonaxisymmetric oscillation modes in mergers of compact object binaries

We study the excitation of nonaxisymmetric modes in the post-merger phase of binary compact object mergers and the associated gravitational wave emission. Our analysis is based on general-relativistic simulations, in the spatial conformal flatness approximation, using smoothed-particle-hydrodynamics for the evolution of matter, and we use a set of equal and unequal mass models, described by two nonzero-temperature hadronic equations of state and by one strange star equation of state. Through Fourier transforms of the evolution of matter variables, we can identify a number of oscillation modes, as well as several nonlinear components (combination frequencies). We focus on the dominant m=2 mode, which forms a triplet with two nonlinear components that are the result of coupling to the quasiradial mode. A corresponding triplet of frequencies is identified in the gravitational wave spectrum, when the individual masses of the compact objects are in the most likely range of 1.2 to 1.35 $M_\odot$. We can thus associate, through direct analysis of the dynamics of the fluid, a specific frequency peak in the gravitational wave spectrum with the nonlinear component resulting from the difference between the m=2 mode and the quasiradial mode. Once such observations become available, both the m=2 and quasiradial mode frequencies could be extracted, allowing for the application of gravitational-wave asteroseismology to the post-merger remnant and leading to tight constraints on the equation of state of high-density matter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetic field amplification in collapsing, non-rotating stellar cores

Context. The influence of magnetic fields on stellar core collapse and explosion is not well explored. It depends on the possibility to amplify the pre-collapse fields. Without rotation this can happen by compression, convection, the standing accretion shock instability, and the accumulation and growth of Alfven waves in the accretion flows. Aims. We investigate such amplification mechanisms of the magnetic field during the collapse and post-bounce evolution of the core of a non-rotating 15 solar mass star with varied initial field strengths, taking into account the microphysical equation of state and neutrino physics crucial for supernova cores. Methods. We perform simulations of ideal MHD with neutrino transport in axisymmetry. The transport of electron neutrinos and antineutrinos is treated with a new scheme that solves the energy-dependent set of radiation energy and momentum equations in 2d by using an analytic closure relation. Results. The magnetic field undergoes amplification by turbulent flows. We find indications for amplification by interacting waves in accretion streams. The fields can reach up to equipartition with the velocity field. Very high magnetic field strengths require very strong pre-collapse fields and are able to shape the post-bounce flow, leading to a pattern dominated by low-order multipoles. Such models are closest to a successful explosion. Conclusions. Magnetic fields can build up to interesting strengths even in non-rotating collapsing stellar cores. Starting with fields in the pre-collapse core as predicted by present stellar evolution models, typical neutron star fields emerge, whereas stronger progenitor fields lead to fields of magnetar strength. Only in the latter case the fields have dynamical effects on the flows in the supernova core. However, in none of our 2D simulations, we find an explosion until 500 ms post-bounce.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The r-process in the neutrino-driven wind from a black-hole torus

We examine r-process nucleosynthesis in the neutrino-driven wind from the thick accretion disk (or "torus") around a black hole. Such systems are expected as emnants of binary neutron star or neutron star -- black hole mergers. We consider a simplified, analytic, time-dependent evolution model of a 3M_sun central black hole surrounded by a neutrino emitting accretion torus with 90km radius, which serves as basis for computing spherically symmetric neutrino-driven wind solutions. We find that ejecta with modest entropies (~30 per nucleon in units of the Boltzmann constant) and moderate expansion timescales (~100ms) dominate in the mass outflow. The mass-integrated nucleosynthetic abundances are in good agreement with the solar system r-process abundance distribution if a minimal value of the electron fraction at the charged-particle freezeout, Ye,min~0.2, is achieved. In the case of Ye,min~0.3, the production of r-elements beyond A~130 does not reach to the third peak but could be still important for an explanation of the abundance signatures in r-process deficient stars in the early Galaxy. The total mass of the ejected r-process nuclei is estimated to be ~1x10^-3 M_sun. If our model was representative, this demands a Galactic event rate of ~2x10^-4 yr^-1 for black-hole-torus winds from merger remnants to be the dominant source of the r-process elements. Our result thus suggests that black-hole-torus winds from compact binary mergers have the potential to be a major but probably not the dominant production site of r-process elements.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Electron-capture supernovae as the origin of elements beyond iron

We examine electron-capture supernovae (ECSNe) as sources of elements heavier than iron in the solar system and in Galactic halo stars. Nucleosynthesis calculations are performed on the basis of thermodynamic histories of mass elements from a fully self-consistent, two-dimensional (2D) hydrodynamic explosion model of an ECSN. We find that neutron-rich convective lumps with an electron fraction down to Ye,min=0.40, which are absent in the one-dimensional (1D) counterpart, allow for interesting production of elements between the iron group and N=50 nuclei (from Zn to Zr, with little Ga) in nuclear (quasi-)equilibrium. Our models yield very good agreement with the Ge, Sr, Y, and Zr abundances of r-process deficient Galactic halo stars and constrain the occurrence of ECSNe to ~4% of all stellar core-collapse events. If tiny amounts of additional material with slightly lower Ye,min down to ~0.30-0.35 were also ejected - which presently cannot be excluded because of the limitations of resolution and two-dimensionality of the model -, a weak r-process can yield elements beyond N=50 up to Pd, Ag, and Cd as observed in the r-process deficient stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fast time variations of supernova neutrino fluxes and their detectability

In the delayed explosion scenario of core-collapse supernovae (SNe), the accretion phase shows pronounced convective overturns and a low-multipole hydrodynamic instability, the standing accretion shock instability (SASI). These effects imprint detectable fast time variations on the emerging neutrino flux. Among existing detectors, IceCube is best suited to this task, providing an event rate of ~1000 events per ms during the accretion phase for a fiducial SN distance of 10 kpc, comparable to what could be achieved with a megaton water Cherenkov detector. If the SASI activity lasts for several hundred ms, a Fourier component with an amplitude of 1% of the average signal clearly sticks out from the shot noise. We analyze in detail the output of axially symmetric hydrodynamical simulations that predict much larger amplitudes up to frequencies of a few hundred Hz. If these models are roughly representative for realistic SNe, fast time variations of the neutrino signal are easily detectable in IceCube or future megaton-class instruments. We also discuss the information that could be deduced from such a measurement about the physics in the SN core and the explosion mechanism of the SN.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The r-Process in Black Hole Winds

All the current r-process scenarios relevant to core-collapse supernovae are facing severe difficulties. In particular, recent core-collapse simulations with neutrino transport show no sign of a neutron-rich wind from the proto-neutron star. In this paper, we discuss nucleosynthesis of the r-process in an alternative astrophysical site, "black hole winds", which are the neutrino-driven outflow from the accretion torus around a black hole. This condition is assumed to be realized in double neutron star mergers, neutron star - black hole mergers, or hypernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Uncertainties in the nu p-process: supernova dynamics versus nuclear physics

We examine how the uncertainties involved in supernova dynamics as well as in nuclear data inputs affect the nup-process in the neutrino-driven winds. For the supernova dynamics, we find that the wind termination by the preceding dense ejecta shell, as well as the electron fraction (Y_{e, 3}; at 3 10^9 K) play a crucial role. A wind termination within the temperature range of (1.5-3) 10^9 K greatly enhances the efficiency of the nu p-process. This implies that the early wind phase, when the innermost layer of the preceding supernova ejecta is still 200-1000 km from the center, is most relevant to the nup-process. The outflows with Y_{e, 3} = 0.52-0.60 result in the production of the p-nuclei up to A=108 with interesting amounts. Furthermore, the p-nuclei up to A=152 can be produced if Y_{e, 3} = 0.65 is achieved. For the nuclear data inputs, we test the sensitivity to the rates relevant to the breakout from the pp-chain region (A < 12), to the (n, p) rates on heavy nuclei, and to the nuclear masses along the nup-process pathway. We find that a small variation of the rates of triple-alpha and of the (n, p) reaction on 56Ni leads to a substantial change in the p-nuclei production. We also find that 96Pd (N=50) on the nup-process path plays a role as a second seed nucleus for the production of heavier p-nuclei. The uncertainty in the nuclear mass of 82Zr can lead to a factor of two reduction in the abundance of the p-isotope 84Sr.

0 citations0 reviewsNo code linkedOpen access