We study the impact of a small-scale dynamo in core-collapse supernovae using a 3D neutrino magnetohydrodynamics simulation of a $15 M_\odot$ progenitor. The weak seed field is amplified exponentially in the gain region once neutrino-driven convection develops, and remains dominated by small-scale structures. About $250\, \mathrm{ms}$ after bounce, the field energy in the gain region reaches $\mathord{\sim} 50\%$ of kinetic equipartition. This supports the development of a neutrino-driven explosion with modest global anisotropy, which does not occur in a corresponding model without magnetic fields. Our results suggest that magnetic fields may play a beneficial subsidiary role in neutrino-driven supernovae even without rapid progenitor rotation. Further investigation into the nature of magnetohydrodynamic turbulence in the supernova core is required.

Based on recent results from three-dimensional supernova simulations and semi-analytical parametrised models, we develop analytical prescriptions for the dependence of the mass of neutron stars and black holes and the natal kicks, if any, on the pre-supernova carbon-oxygen core and helium shell masses. Our recipes are probabilistic rather than deterministic in order to account for the intrinsic stochasticity of stellar evolution and supernovae. We anticipate that these recipes will be particularly useful for rapid population synthesis, and we illustrate their application to distributions of remnant masses and kicks for a population of single stars.

Multi-dimensional fluid flow plays a paramount role in the explosions of massive stars as core-collapse supernovae. In recent years, three-dimensional (3D) simulations of these phenomena have matured significantly. Considerable progress has been made towards identifying the ingredients for shock revival by the neutrino-driven mechanism, and successful explosions have already been obtained in a number of self-consistent 3D models. These advances also bring new challenges, however. Prompted by a need for increased physical realism and meaningful model validation, supernova theory is now moving towards a more integrated view that connects multi-dimensional phenomena in the late convective burning stages prior to collapse, the explosion engine, and mixing instabilities in the supernova envelope. Here we review our current understanding of multi-D fluid flow in core-collapse supernovae and their progenitors. We start by outlining specific challenges faced by hydrodynamic simulations of core-collapse supernovae and of the late convective burning stages. We then discuss recent advances and open questions in theory and simulations.

We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor stars. This requires a hypernova occurence rate that is large (50%) in the early Universe, decreasing throughout evolution to a value that is within present day observational constraints (>~ 1%). A large hypernova occurence rate is beneficial to matching the high [Zn/Fe] observed in the most metal-poor stars, although including hypernovae with progenitor mass >= 60 solar masses is detrimental to matching the observed [(Mn,Co)/Fe] evolution at low [Fe/H]. A significant contribution from HNe seems to be critical for producing supersolar [(Co,Zn)/Fe] at low metallicity, though more work will need to be done in order to match the most extreme values. We also emphasise the need to update models for the enrichment sources at higher metallicity, as the satisfactory recovery of the solar values of [(Cr,Mn,Co,Zn)/Fe] still presents a challenge.

It has been suggested based on analytic theory that even in non-rotating supernova progenitors stochastic spin-up by internal gravity waves (IGWs) during the late burning stages can impart enough angular momentum to the core to result in neutron star birth spin periods below 100ms, and a relatively firm upper limit of 500ms for the spin period. We here investigate this process using a 3D simulation of oxygen shell burning in a $3M_\odot$ He star. Our model indicates that stochastic spin-up by IGWs is less efficient than previously thought. We find that the stochastic angular momentum flux carried by waves excited at the shell boundary is significantly smaller for a given convective luminosity and turnover time than would be expected from simple dimensional analysis. This can be explained by noting that the waves launched by overshooting convective plumes contain modes of opposite angular wave number with similar amplitudes, so that the net angular momentum of excited wave packets almost cancels. We find that the wave-mediated angular momentum flux from the oxygen shell follows a random walk, but again dimensional analysis overestimates the random walk amplitudes since the correlation time is only a fraction of the convective turnover time. Extrapolating our findings over the entire life time of the last burning stages prior to collapse, we predict that the core angular momentum from stochastic spin-up would translate into long birth spin periods of several seconds for low-mass progenitors and no less than 100ms even for high-mass progenitors.

We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200 Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$ Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma & Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers.

We present three-dimensional simulations of the core-collapse of massive rotating and non-rotating progenitors performed with the general relativistic neutrino hydrodynamics code CoCoNuT-FMT and analyse their explosion properties and gravitational-wave signals. The progenitor models include Wolf-Rayet stars with initial helium star masses of $39\,M_{\odot}$ and $20\,M_{\odot}$, and an $18\,M_{\odot}$ red supergiant. The $39\,M_{\odot}$ model is a rapid rotator, whereas the two other progenitors are non-rotating. Both Wolf-Rayet models produce healthy neutrino-driven explosions, whereas the red supergiant model fails to explode. By the end of the simulations, the explosion energies have already reached $1.1\times 10^{51}\,\mathrm{erg}$ and $0.6\times 10^{51}\,\mathrm{erg}$ for the $39\,M_{\odot}$ and $20\,M_{\odot}$ model, respectively. The explosions produce neutron stars of relatively high mass, but with modest kicks. Due to the alignment of the bipolar explosion geometry with the rotation axis, there is a relatively small misalignment of $30^\circ$ between the spin and the kick in the $39\,M_{\odot}$ model. In terms of gravitational-wave signals, the massive and rapidly rotating $39\,M_{\odot}$ progenitor stands out by large gravitational-wave amplitudes that would make it detectable out to almost 2 Mpc by the Einstein Telescope. For this model, we find that rotation significantly changes the dependence of the characteristic gravitational-wave frequency of the f-mode on the proto-neutron star parameters compared to the non-rotating case. The other two progenitors have considerably smaller detection distances, despite significant low-frequency emission in the most sensitive frequency band of current gravitational-wave detectors due to the standing accretion shock instability in the $18\,M_{\odot}$ model.

We present a seven-minute long $4\pi$-3D simulation of a shell merger event in a non-rotating $18.88\, M_\odot$ supernova progenitor before the onset of gravitational collapse. The key motivation is to capture the large-scale mixing and asymmetries in the wake of the shell merger before collapse using a self-consistent approach. The $4\pi$ geometry is crucial as it allows us to follow the growth and evolution of convective modes on the largest possible scales. We find significant differences between the kinematic, thermodynamic and chemical evolution of the 3D and the 1D model. The 3D model shows vigorous convection leading to more efficient mixing of nuclear species. In the 3D case the entire oxygen shell attains convective Mach numbers of $\mathord{\approx}\, 0.1$, whereas in the 1D model, the convective velocities are much lower and there is negligible overshooting across convective boundaries. In the 3D case, the convective eddies entrain nuclear species from the neon (and carbon) layers into the deeper part of the oxygen burning shell, where they burn and power a violent convection phase with outflows. This is a prototypical model of a convective-reactive system. Due to the strong convection and the resulting efficient mixing, the interface between the neon layer and the silicon-enriched oxygen layer disappears during the evolution, and silicon is mixed far out into merged oxygen/neon shell. Neon entrained inwards by convective downdrafts burns, resulting in lower neon mass in the 3D model compared to the 1D model at time of collapse. In addition, the 3D model develops remarkable large-scale, large-amplitude asymmetries, which may have important implications for the impending gravitational collapse and the subsequent explosion.

We study the effect of tidal forcing on gravitational wave signals from tidally relaxed white dwarf pairs in the LISA, DECIGO and BBO frequency band ($0.1-100\,{\rm mHz}$). We show that for stars not in hydrostatic equilibrium (in their own rotating frames), tidal forcing will result in energy and angular momentum exchange between the orbit and the stars, thereby deforming the orbit and producing gravitational wave power in harmonics not excited in perfectly circular synchronous binaries. This effect is not present in the usual orbit-averaged treatment of the equilibrium tide, and is analogous to transit timing variations in multiplanet systems. It should be present for all LISA white dwarf pairs since gravitational waves carry away angular momentum faster than tidal torques can act to synchronize the spins, and when mass transfer occurs as it does for at least eight LISA verification binaries. With the strain amplitudes of the excited harmonics depending directly on the density profiles of the stars, gravitational wave astronomy offers the possibility of studying the internal structure of white dwarfs, complimenting information obtained from asteroseismology of pulsating white dwarfs. Since the vast majority of white-dwarf pairs in this frequency band are expected to be in the quasi-circular state, we focus here on these binaries, providing general analytic expressions for the dependence of the induced eccentricity and strain amplitudes on the stellar apsidal motion constants and their radius and mass ratios. Tidal dissipation and gravitation wave damping will affect the results presented here and will be considered elsewhere.

Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Coherent THz radiation allows direct measurement of the complex dielectric function (refractive index) of astrophysically relevant ice species in the THz range. Methods: The time-domain waveforms and the frequency-domain spectra of reference samples of CO ice, deposited at a temperature of 28.5 K and annealed to 33 K at different thicknesses, have been recorded. A new algorithm is developed to reconstruct the real and imaginary parts of the refractive index from the time-domain THz data. Results: The complex refractive index in the wavelength range of 1 mm - 150 ${\mu}$m (0.3 - 2.0 THz) has been determined for the studied ice samples, and compared with available data found in the literature. Conclusions: The developed algorithm of reconstructing the real and imaginary parts of the refractive index from the time-domain THz data enables, for the first time, the determination of optical properties of astrophysical ice analogs without using the Kramers-Kronig relations. The obtained data provide a benchmark to interpret the observational data from current ground based facilities as well as future space telescope missions, and have been used to estimate the opacities of the dust grains in presence of CO ice mantles.

It has recently been proposed that global or local turbulence models can be used to simulate core-collapse supernova explosions in spherical symmetry (1D) more consistently than with traditional approaches for parameterised 1D models. However, a closer analysis of the proposed schemes reveals important consistency problems. Most notably, they systematically violate energy conservation as they do not balance buoyant energy generation with terms that reduce potential energy, thus failing to account for the physical source of energy that buoyant convection feeds on. We also point out other non-trivial consistency requirements for viable turbulence models. The model of Kuhfuss (1986) proves more consistent than the newly proposed approaches for supernovae, but still cannot account naturally for all the relevant physics for predicting explosion properties. We perform numerical simulations for a $20 M_\odot$ progenitor to further illustrate problems of 1D turbulence models. If the buoyant driving term is formulated in a conservative manner, the explosion energy of $\mathord{\sim}2\times 10^{51}\, \mathrm{erg}$ for the corresponding non-conservative turbulence model is reduced to $<10^{48}\,\mathrm{erg}$ even though the shock expands continuously. This demonstrates that the conservation problem cannot be ignored. Although plausible energies can be reached using an energy-conserving model when turbulent viscosity is included, it is doubtful whether the energy budget of the explosion is regulated by the same mechanism as in multi-dimensional models. We conclude that 1D turbulence models based on a spherical Reynolds decomposition cannot provide a more consistent approach to supernova explosion and remnant properties than other phenomenological approaches before some fundamental problems are addressed.

Understanding gravitational wave emission from core-collapse supernovae will be essential for their detection with current and future gravitational wave detectors. This requires a sample of waveforms from modern 3D supernova simulations reaching well into the explosion phase, where gravitational wave emission is expected to peak. However, recent waveforms from 3D simulations with multi-group neutrino transport do not reach far into the explosion phase, and some are still obtained from non-exploding models. We therefore calculate waveforms up to 0.9\,s after bounce using the neutrino hydrodynamics code \textsc{CoCoNuT-FMT}. We consider two models with low and normal explosion energy, namely explosions of an ultra-stripped progenitor with an initial helium star mass of $3.5\,M_{\odot}$, and of an $18\,M_{\odot}$ single star. Both models show gravitational wave emission from the excitation of surface g-modes in the proto-neutron star with frequencies between $\mathord{\sim}800\,\mathrm{Hz}$ and 1000\,Hz at peak emission. The peak amplitudes are about $6\, \mathrm{cm}$ and $10\, \mathrm{cm}$, respectively, which is somewhat higher than in most recent 3D models of the pre-explosion or early explosion phase. Using a Bayesian analysis, we determine the maximum detection distances for our models in simulated Advanced LIGO, Advanced Virgo, and Einstein Telescope design sensitivity noise. The more energetic $18 M_\odot$ explosion will be detectable to about $17.5 \,\mathrm{kpc}$ by the LIGO/Virgo network and to about $180\, \mathrm{kpc}$ with the Einstein Telescope.

With myriads of detection events from a prospective Galactic core-collapse supernova, current and future neutrino detectors will be able to sample detailed, time-dependent neutrino fluxes and spectra. This offers enormous possibilities for inferring supernova physics from the various phases of the neutrino signal from the neutronization burst through the accretion and early explosion phase to the cooling phase. The signal will constrain the time evolution of bulk parameters of the young proto-neutron star like its mass and radius as well as the structure of the progenitor, probe multi-dimensional phenomena in the supernova core, and constrain thedynamics of the early explosion phase. Aside from further astrophysical implications, supernova neutrinos may also shed further light on the properties of matter at supranuclear densities and on open problems in particle physics.

State-of-the-art numerical simulations of core-collapse supernovae reveal that the main source of gravitational waves is the excitation of proto-neutron star modes during post-bounce evolution. In this work we derive universal relations that relate the frequencies of the most common oscillation modes observed, i.e. g-modes, p-modes and the f-mode, with fundamental properties of the system, such as the surface gravity of the proto-neutron star or the mean density in the region enclosed by the shock. These relations are independent of the equation of state, the neutrino treatment, and the progenitor mass and hence can be used to build methods to infer proto-neutron star properties from gravitational-wave observations alone. We outline how these measurements could be done and the constraints that could be placed on the proto-neutron star properties.

The solution of the Poisson equation is a ubiquitous problem in computational astrophysics. Most notably, the treatment of self-gravitating flows involves the Poisson equation for the gravitational field. In hydrodynamics codes using spherical polar grids, one often resorts to a truncated spherical harmonics expansion for an approximate solution. Here we present a non-iterative method that is similar in spirit, but uses the full set of eigenfunctions of the discretized Laplacian to obtain an exact solution of the discretized Poisson equation. This allows the solver to handle density distributions for which the truncated multipole expansion fails, such as off-center point masses. In three dimensions, the operation count of the new method is competitive with a naive implementation of the truncated spherical harmonics expansion with $N_\ell \approx 15$ multipoles. We also discuss the parallel implementation of the algorithm. The serial code and a template for the parallel solver are made publicly available.

We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.

Aims. In this paper we investigate the detectability of the molecular oxygen in icy dust grain mantles towards astronomical objects. Methods. We present a systematic set of experiments with O2-H2O ice mixtures designed to disentangle how the molecular ratio affects the O2 signature in the mid- and near-infrared spectral regions. All the experiments were conducted in a closed-cycle helium cryostat coupled to a Fourier transform infrared spectrometer. The ice mixtures comprise varying thicknesses from 8 $\times$ 10$^{-3}$ to 3 $\mu$m. The absorption spectra of the O2-H2O mixtures are also compared to the one of pure water. In addition, the possibility to detect the O2 in icy bodies and in the interstellar medium is discussed. Results. We are able to see the O2 feature at 1551 cm$^{-1}$ even for the most diluted mixture of H2O : O2 = 9 : 1, comparable to a ratio of O2/H2O = 10 % which has already been detected in situ in the coma of the comet 67P/Churyumov-Gerasimenko. We provide an estimate for the detection of O2 with the future mission of the James Webb Space Telescope (JWST).

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions.

The recent discoveries of many double neutron star systems and their detection as LIGO-Virgo merger events call for a detailed understanding of their origin. Explosions of ultra-stripped stars in binary systems have been shown to play a key role in this context and have also generated interest as a potential explanation for rapidly evolving hydrogen-free transients. Here we present the first attempt to model such explosions based on binary evolution calculations that follow the mass transfer to the companion to obtain a consistent core-envelope structure as needed for reliable predictions of the supernova transient. We simulate the explosion in 2D and 3D, and confirm the modest explosion energies ~10^50erg and small kick velocities reported earlier in 2D models based on bare carbon-oxygen cores. The spin-up of the neutron star by asymmetric accretion is small in 3D with no indication of spin-kick alignment. Simulations up to shock breakout show the mixing of sizeable amounts of iron group material into the helium envelope. In view of recent ideas for a mixing-length treatment (MLT) of Rayleigh-Taylor instabilities in supernovae, we perform a detailed analysis of the mixing, which reveals evidence for buoyancy-drag balance, but otherwise does not support the MLT approximation. The mixing may have implications for the spectroscopic signatures of ultra-stripped supernovae that need to be investigated in the future. Our stellar evolution calculation also predicts presupernova mass loss due to an off-centre silicon deflagration flash, which suggests that supernovae from extremely stripped cores may show signs of interactions with circumstellar material.

We examine the nucleosynthesis in the innermost, neutrino-processed ejecta (a few $10^{-3}\,M_\odot$) of self-consistent, two-dimensional explosion models of core-collapse supernovae for six progenitor stars with different initial masses. Three models have initial masses near the low-mass end of the supernova range, $8.8\,M_\odot$ (e8.8; electron-capture supernova), $9.6\,M_\odot$ (z9.6), and $8.1\,M_\odot$ (u8.1), with initial metallicities of 1, 0, and $10^{-4}$ times the solar metallicity, respectively. The other three are solar-metallicity models with initial masses of $11.2\,M_\odot$ (s11), $15\,M_\odot$ (s15), and $27\,M_\odot$ (s27). The low-mass models e8.8, z9.6, and u8.1 exhibit high production factors (nucleosynthetic abundances relative to the solar ones) of 100--200 for light trans-iron elements from Zn to Zr. This is associated with appreciable ejection of neutron-rich matter in these models. Remarkably, the nucleosynthetic outcomes for progenitors e8.8 and z9.6 are almost identical, including interesting productions of $^{48}$Ca and $^{60}$Fe, irrespective of their quite different (O-Ne-Mg and Fe) cores prior to collapse. In the more massive models s11, s15, and s27, several proton-rich isotopes of light trans-iron elements, including the $p$-isotope $^{92}$Mo (for s27) are made, up to production factors of $\sim$30. Both electron-capture and core-collapse supernovae near the low-mass end can therefore be dominant contributors to the Galactic inventory of light trans-iron elements from Zn to Zr and probably $^{48}$Ca and live $^{60}$Fe. The innermost ejecta of more massive supernovae may have only sub-dominant contributions to the chemical enrichment of the Galaxy except for $^{92}$Mo.

The evolutionary pathways of core-collapse supernova progenitors at the low-mass end of the spectrum are beset with major uncertainties. In recent years, a variety of evolutionary channels has been discovered in addition to the classical electron capture supernova channel of super-AGB stars. The few available progenitor models at the low-mass end have been studied with great success in supernova simulations as the peculiar density structure makes for robust neutrino-driven explosions in this mass range. Detailed nucleosynthesis calculations have been conducted both for models of electron capture supernovae and low-mass iron core supernovae and revealed an interesting production of the lighter trans-iron elements (such as Zn, Sr, Y, Zr) as well as rare isotopes like Ca-48 and Fe-60. We stress the need to explore the low-mass end of the supernova spectrum further and link various observables to understand the diversity of explosions in this regime.

Fallback in core-collapse supernovae is considered a major ingredient for explaining abundance anomalies in metal-poor stars and the natal kicks and spins of black holes (BHs). We present a first 3D simulation of BH formation and fallback in an "aborted" neutrino-driven explosion of a $40$ solar mass zero-metallicity progenitor from collapse to shock breakout. We follow the phase up to BH formation using the relativistic CoCoNuT-FMT code. For the subsequent evolution to shock breakout we apply the moving-mesh code Arepo to core-collapse supernovae for the first time. Our simulation shows that despite early BH formation, neutrino-heated bubbles can survive for tens of seconds before being accreted, leaving them sufficient time to transfer part of their energy to sustain the shock wave as is propagates through the envelope. Although the initial net energy ($\sim 2$ Bethe) of the neutrino-heated ejecta barely equals the binding energy of the envelope, $11\,\mathrm{M}_\odot$ of hydrogen are still expelled with an energy of $0.23$ Bethe. We find no significant mixing and only a modest BH kick and spin, but speculate that stronger effects could occur for slightly more energetic explosions or progenitors with less tightly bound envelopes.

Recent three-dimensional simulations have suggested that convective seed perturbations from shell burning can play an important role in triggering neutrino-driven supernova explosions. Since isolated simulations cannot determine whether this perturbation-aided mechanism is of general relevance across the progenitor mass range, we here investigate the pertinent properties of convective oxygen and silicon burning shells in a broad range of presupernova stellar evolution models. We find that conditions for perturbation-aided explosions are most favourable in the extended oxygen shells of progenitors between about 16 and 26 solar masses, which exhibit large-scale convective overturn with high convective Mach numbers. Although the highest convective Mach numbers of up to 0.3 are reached in the oxygen shells of low-mass progenitors, convection is typically dominated by small-scale modes in these shells, which implies a more modest role of initial perturbations in the explosion mechanism. Convective silicon burning rarely provides the high Mach numbers and large-scale perturbations required for perturbation-aided explosions. We also find that about $40\%$ of progenitors between 16 and 26 solar masses exhibit simultaneous oxygen and neon burning in the same convection zone as a result of a shell merger shortly before collapse.

We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D.

Gravitational waves are a potential direct probe for the multi-dimensional flow during the first second of core-collapse supernova explosions. Here we outline the structure of the predicted gravitational wave signal from neutrino-driven supernovae of non-rotating progenitors from recent 2D and 3D simulations. We sketch some quantitative dependencies that govern the amplitudes of this signal and its evolution in the time-frequency domain.

The explosion mechanism of core-collapse supernovae is a long-standing problem in stellar astrophysics. We briefly outline the main contenders for a solution and review recent efforts to model core-collapse supernova explosions by means of multi-dimensional simulations. We discuss several suggestions for solving the problem of missing or delayed neutrino-driven explosions in three-dimensional supernova models, including -- among others -- variations in the microphysics and large seed perturbations in convective burning shells. Focusing on the neutrino-driven mechanism, we summarise currents efforts to predict supernova explosion and remnant properties based on first-principle models and on more phenomenological approaches.

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.

Models of core-collapse supernova explosions powered by the neutrino-driven mechanism have matured considerable in recent years. Explosions at the low-mass end of the progenitor spectrum can routinely be simulated in 1D, 2D, and 3D and allow us to study supernova nucleosynthesis based on first-principle models. Results of nucleosynthesis calculations indicate that supernovae of the lowest masses could be important contributors of some lighter n-rich elements beyond iron. The explosion mechanism of more massive stars is still under investigation, although first 3D models of neutrino-driven explosions employing multi-group neutrino transport have recently become available. Together with earlier 2D models and more simplified 3D simulations, these have elucidated the interplay between neutrino heating and hydrodynamic instabilities in the post-shock region that is essential for shock revival. However, some physical ingredients may still need to be added or improved before simulations can robustly explain supernova explosions over a wide mass range. We explore possible issues that may affect the accuracy of supernova simulations, and review some of the ideas that have recently been explored as avenues to robust explosions, including uncertainties in the neutrino rates, rapid rotation, and an external forcing of non-radial fluid motions by strong seed perturbations from convective shell burning. The perturbation-aided neutrino-driven mechanism and the implications of recent 3D simulations of shell burning in supernova progenitors are discussed in detail. The efficacy of the perturbation-aided mechanism is illustrated by the first successful multi-group neutrino hydrodynamics simulation of an 18 solar mass progenitor with 3D initial conditions. We conclude with speculations about the potential impact of 3D effects on the structure of massive stars through convective boundary mixing. (abridged)

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.

We present a new approach to understand the landscape of supernova explosion energies, ejected nickel masses, and neutron star birth masses. In contrast to other recent parametric approaches, our model predicts the properties of neutrino-driven explosions based on the pre-collapse stellar structure without the need for hydrodynamic simulations. The model is based on physically motivated scaling laws and simple differential equations describing the shock propagation, the contraction of the neutron star, the neutrino emission, the heating conditions, and the explosion energetics. Using model parameters compatible with multi-D simulations and a fine grid of thousands of supernova progenitors, we obtain a variegated landscape of neutron star and black hole formation similar to other parameterised approaches and find good agreement with semi-empirical measures for the "explodability" of massive stars. Our predicted explosion properties largely conform to observed correlations between the nickel mass and explosion energy. Accounting for the coexistence of outflows and downflows during the explosion phase, we naturally obtain a positive correlation between explosion energy and ejecta mass. These correlations are relatively robust against parameter variations, but our results suggest that there is considerable leeway in parametric models to widen or narrow the mass ranges for black hole and neutron star formation and to scale explosion energies up or down. Our model is currently limited to an all-or-nothing treatment of fallback and there remain some minor discrepancies between model predictions and observational constraints.

We study the growth of the explosion energy after shock revival in neutrino-driven explosions in two and three dimensions (2D/3D) using multi-group neutrino hydrodynamics simulations of an $11.2 M_\odot$ star. The 3D model shows a faster and steadier growth of the explosion energy and already shows signs of subsiding accretion after one second. By contrast, the growth of the explosion energy in 2D is unsteady, and accretion lasts for several seconds as confirmed by additional long-time simulations of stars of similar masses. Appreciable explosion energies can still be reached, albeit at the expense of rather high neutron star masses. In 2D, the binding energy at the gain radius is larger because the strong excitation of downward-propagating g-modes removes energy from the freshly accreted material in the downflows. Consequently, the mass outflow rate is considerably lower in 2D than in 3D. This is only partially compensated by additional heating by outward-propagating acoustic waves in 2D. Moreover, the mass outflow rate in 2D is reduced because much of the neutrino energy deposition occurs in downflows or bubbles confined by secondary shocks without driving outflows. Episodic constriction of outflows and vertical mixing of colder shocked material and hot, neutrino-heated ejecta due to Rayleigh-Taylor instability further hamper the growth of the explosion energy in 2D. Further simulations will be necessary to determine whether these effects are generic over a wider range of supernova progenitors.

Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $\eta \approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation.

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.

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.

The usage of the anthropic principle in modern cosmology is reviewed. It is argued that its recent use to explain the observedvalues of cosmological parameters as most probable values for an ensemble of universes, is not justified. However, the anthropic principle can be invoked to argue that a vast number of universes must exist, in which the cosmological parameters and constants of nature take on different values. This argument lends support to the hypothesis of eternal inflation and suggests that at least some of the parameters of the Standard Model cannot be quantitatively derived from an underlying theory.