Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

We introduce a new discrete-ordinate scheme for solving the general relativistic Boltzmann transport equation in the context of core-collapse supernovae. Our algorithm avoids the need to spell out the complicated advection terms in energy and angle that arise when the transport equation is formulated in spherical polar coordinates, in the comoving frame, or in a general relativistic spacetime. We instead approach the problem by calculating the advection of neutrinos across momentum space using an intuitive particle-like approach that has excellent conservation properties and fully accounts for Lorentz boosts, general relativistic effects, and grid geometry terms. In order to avoid the need for a global implicit solution, time integration is performed using a locally implicit Lax-Wendroff scheme that correctly reproduces the diffusion limit. This will facilitate the use of our method on massively parallel distributed-memory architectures. We have verified the accuracy and stability of our scheme with a suite of test problems in spherical symmetry and axisymmetry. To demonstrate that the new algorithm works stably in core-collapse supernova simulations, we have coupled it to the general relativistic hydrodynamics code CoCoNuT and present a first demonstration run of a $20M_\odot$ progenitor with a reduced set of neutrino opacities.

Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spin is imparted by partial fallback in an asymmetric explosion. Black hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the initial explosion energy barely exceeds the envelope binding energy, stronger fallback results, and the remnant kick and spin remain small. If the explosion energy is high with respect to the envelope binding energy, there is little fallback with a small effect on the remnant kick, but the spin-up by fallback can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we find that the neutron star (NS) is spun up to millisecond periods. The high specific angular momentum of the fallback material can also lead to disk formation around black holes. Fallback may thus be a pathway towards millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray bursts that does not require rapid progenitor rotation. Within our small set of simulations, none reproduced the peculiar layered fallback necessary to explain the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP) stars. Models with different explosion energy and different realisations of asymmetries may, however, be compatible with CEMP abundance patterns.

Observations with next-generation ground-based detectors further enhanced with multi-messenger (electromagnetic and neutrino) detections will allow us to probe new extreme astrophysics. Target sources included: core-collapse supernovae, continuous emission from isolated or accreting neutron stars, and bursts from magnetars and other pulsars.

Samarasinha & Mueller (2013) related changes of cometary rotation to other physical parameters for four Jupiter family comets defining a parameter $X$, which is approximately constant within a factor of two irrespective of the active fraction of a comet. Two additional comets are added to this sample in this paper and the claim of a nearly constant parameter $X$ for these six comets is confirmed, albeit with a larger scatter. Taking the geometric mean of $X$ for all the comets above excluding 2P/Encke (as $X$ for each comet was determined with respect to that of 2P/Encke), the expected changes in the rotation periods for a sample of 24 periodic comets are derived. We identify comets from this sample that are most likely to show observationally detectable changes in their rotation periods. Using this sample and including the six comets used to determine $X$, we find a correlation between the parameter $\zeta$ (i.e. the total water production per unit surface area per orbit approximated by that inside of 4 au) and the perihelion distance $q$; specifically we derive $\zeta$ $\propto$ $q^{-0.8}$ and provide a theoretical basis for this in Appendix A. This relationship between $\zeta$ and $q$ enables ready comparisons of activity due to insolation between comets. Additionally, a relationship between the nuclear radius $R$ and the rotation period $P$ is found. Specifically, we find that on average smaller nuclei have smaller rotation periods compared to the rotation periods of larger nuclei. This is consistent with expectations for rotational evolution and spin-up of comet nuclei, providing strong observational evidence for sublimation-driven rotational changes in comets.

We investigate the relationship between explosion energy and nucleosynthesis in Population III supernovae and provide nucleosynthetic results for the explosions of stars with progenitor masses of $15\,\mathrm{M}_\odot$, $20\,\mathrm{M}_\odot$, $30\,\mathrm{M}_\odot$, $40\,\mathrm{M}_\odot$, $60\,\mathrm{M}_\odot$, and $80\,\mathrm{M}_\odot$, and explosion energies between approximately $10^{50}$ erg and $10^{53}$ erg. We find that the typical abundance pattern observed in metal-poor stars are best matched by supernovae with progenitor mass in the range $15\,\mathrm{M}_\odot$ - $30\,\mathrm{M}_\odot$, and explosion energy of $\sim$ ($5$ - $10$) $\times$ $10^{51}$ erg. In these models, a reverse shock caused by jumps in density between shells of different composition serves to decrease synthesis of chromium and manganese, which is favourable to matching the observed abundances in metal-poor stars. Spherically symmetric explosions of our models with progenitor mass $\geq$ $40\,\mathrm{M}_\odot$ do not provide yields that are compatible with the iron-peak abundances that are typically observed in metal-poor stars, however, by approximating the yields that we might expect from these models in highly aspherical explosions, we find indications that explosions of stars $40\,\mathrm{M}_\odot$ - $80\,\mathrm{M}_\odot$ with bipolar jets may be good candidates for the enrichment sources of metal-poor stars with enhanced carbon abundances.

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum.

We present gravitational wave (GW) signal predictions from four 3D multi-group neutrino hydrodynamics simulations of core-collapse supernovae of progenitors with 11.2 Msun, 20 Msun, and 27 Msun. GW emission in the pre-explosion phase strongly depends on whether the post-shock flow is dominated by the standing accretion shock instability (SASI) or convection and differs considerably from 2D models. SASI activity produces a strong signal component below 250 Hz through asymmetric mass motions in the gain layer and a non-resonant coupling to the proto-neutron star (PNS). Both convection- and SASI-dominated models show GW emission above 250 Hz, but with considerably lower amplitudes than in 2D. This is due to a different excitation mechanism for high-frequency l=2 motions in the PNS surface, which are predominantly excited by PNS convection in 3D. Resonant excitation of high-frequency surface g-modes in 3D by mass motions in the gain layer is suppressed compared to 2D because of smaller downflow velocities and a lack of high-frequency variability in the downflows. In the exploding 20 Msun model, shock revival results in enhanced low-frequency emission due to a change of the preferred scale of the convective eddies in the PNS convection zone. Estimates of the expected excess power in two frequency bands suggests that second-generation detectors will only be able to detect very nearby events, but that third-generation detectors could distinguish SASI- and convection-dominated models at distances of ~10 kpc.

We present a search for an electromagnetic counterpart of the gravitational wave source GW151226. Using the Pan-STARRS1 telescope we mapped out 290 square degrees in the optical i_ps filter starting 11.5hr after the LIGO information release and lasting for a further 28 days. The first observations started 49.5hr after the time of the GW151226 detection. We typically reached sensitivity limits of i_ps = 20.3-20.8 and covered 26.5% of the LIGO probability skymap. We supplemented this with ATLAS survey data, reaching 31% of the probability region to shallower depths of m~19. We found 49 extragalactic transients (that are not obviously AGN), including a faint transient in a galaxy at 7Mpc (a luminous blue variable outburst) plus a rapidly decaying M-dwarf flare. Spectral classification of 20 other transient events showed them all to be supernovae. We found an unusual transient, PS15dpn, with an explosion date temporally coincident with GW151226 which evolved into a type Ibn supernova. The redshift of the transient is secure at z=0.1747 +/- 0.0001 and we find it unlikely to be linked, since the luminosity distance has a negligible probability of being consistent with that of GW151226. In the 290 square degrees surveyed we therefore do not find a likely counterpart. However we show that our survey strategy would be sensitive to NS-NS mergers producing kilonovae at D < 100 Mpc which is promising for future LIGO/Virgo searches.

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 $\nu_e$ minus $\bar\nu_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\nu_e$-$\nu_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.

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.

We present an analysis of Kitt Peak National Observatory and Lowell Observatory observations of comet 103P/Hartley 2 obtained from August through December 2010. The results are then compared with contemporaneous observations made by the EPOXI spacecraft. Each ground-based dataset has previously been investigated individually; the combined dataset has complementary coverage that reduces the time between observing runs and allows us to determine additional apparent periods at intermediate times. We compare CN coma morphology between ground-based datasets, making nine new measurements of apparent periods. The first five are consistent with the roughly linearly increasing apparent period during the apparition found by previous authors. The final four suggest that the change in apparent period slowed or stopped by late November. We also measure an inner coma lightcurve in both CN and R-band ground-based images, finding a single-peaked lightcurve which repeats in phase with the coma morphology. The apparent period from the lightcurve had significantly larger uncertainties than from the coma morphology, but varied over the apparition in a similar manner. Our ground-based lightcurve aligns with the published EPOXI lightcurve, indicating that the lightcurve represents changing activity rather than viewing geometry of structures in the coma. The EPOXI lightcurve can best be phased by a triple-peaked period near 54-55 hr that increases from October to November. This phasing reveals that the spacing between maxima is not constant, and that the overall lightcurve shape evolves from one triple-peaked cycle to the next. These behaviors suggest that much of the scatter in apparent periods derived from ground-based datasets acquired at similar epochs are likely due to limited sampling of the data.

Lightcurve observations of asteroids and bare cometary nuclei are the most widely used observational tool to derive the rotational parameters. Therefore, an in-depth understanding of how component periods of dynamically excited non-principal axis (NPA) rotators manifest in lightcurves is a crucial step in this process. We investigated this with the help of numerically generated lightcurves of NPA rotators with component periods known a priori. The component periods of NPA rotation were defined in terms of two widely used yet complementary conventions. We derive the relationships correlating the component rotation periods in the two conventions. These relationships were then used to interpret the periodicity signatures present in the simulated lightcurves and rationalize them in either convention.

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.

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.

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.

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.

We introduce a parameter, X, to predict the changes in the rotational period of a comet in terms of the rotational period itself, the nuclear radius, and the orbital characteristics. We show that X should be a constant if the bulk densities and shapes of nuclei are nearly identical and the activity patterns are similar for all comets. For four nuclei for which rotational changes are well documented, despite the nearly factor 30 variation observed among the effective active fractions of these comets, X is constant to within a factor two. We present an analysis for the sungrazing comet C/2012 S1 (ISON) to explore what rotational changes it could undergo during the upcoming perihelion passage where its perihelion distance will be ~2.7 solar radii. When close to the sun, barring a catastrophic disruption of the nucleus, the activity of ISON will be sufficiently strong to put the nucleus into a non-principal-axis rotational state and observable changes to the rotational period should also occur. Additional causes for rotational state changes near perihelion for ISON are tidal torques caused by the sun and the significant mass loss due to a number of mechanisms resulting in alterations to the moments of inertia of the nucleus.

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.

The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.

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.

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.

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.

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

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

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

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.

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.

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.

The CN coma structure of the EPOXI mission target, comet 103P/Hartley 2, was observed during twenty nights from September to December 2010. These CN images probe the rotational state of the comet's nucleus and provide a ground-based observational context to complement the EPOXI observations. A dynamically excited cometary nucleus with a changing rotational rate is observed, a characteristic not seen in any comet in the past. The lack of rotational damping during the four-month observing interval places constraints on the interior structure of the nucleus.

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

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.

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

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

We present the generalization of a recently introduced modified gravitational potential for self-gravitating fluids. The use of this potential allows for an accurate approximation of general relativistic effects in an otherwise Newtonian hydrodynamics code also in cases of rapid rotation. We test this approach in numerical simulations of astrophysical scenarios related to compact stars, like supernova core collapse with both a simplified and detailed microphysical description of matter, and rotating neutron stars in equilibrium. We assess the quality of the new potential, and demonstrate that it provides a significant improvement compared to previous formulations for such potentials. Newtonian simulations of compact objects employing such an effective relativistic potential predict inaccurate pulsation frequencies despite the excellent agreement of the collapse dynamics and structure of the compact objects with general relativistic results. We analyze and discuss the reason for this behavior.

We discuss the flavor conversion of neutrinos from core collapse supernovae that have oxygen-neon-magnesium (ONeMg) cores. Using the numerically calculated evolution of the star up to 650 ms post bounce, we find that, for the normal mass hierarchy, the electron neutrino flux in a detector shows signatures of two typical features of an ONeMg-core supernova: a sharp step in the density profile at the base of the He shell and a faster shock wave propagation compared to iron core supernovae. Before the shock hits the density step (t ~ 150 ms), the survival probability of electron neutrinos is about 0.68, in contrast to values of 0.32 or less for an iron core supernova. The passage of the shock through the step and its subsequent propagation cause a decrease of the survival probability and a decrease of the amplitude of oscillations in the Earth, reflecting the transition to a more adiabatic propagation inside the star. These changes affect the lower energy neutrinos first; they are faster and more sizable for larger theta_13. They are unique of ONeMg-core supernovae, and give the possibility to test the speed of the shock wave. The time modulation of the Earth effect and its negative sign at the neutronization peak are the most robust signatures in a detector.

It has been recently proposed that the shocked surface layers of exploding O-Ne-Mg cores provide the conditions for r-process nucleosynthesis, because their rapid expansion and high entropies enable heavy r-process isotopes to form even in an environment with very low initial neutron excess of the matter. We show here that the most sophisticated available hydrodynamic simulations (in spherical and axial symmetry) do not support this new r-process scenario because they fail to provide the necessary conditions of temperature, entropy, and expansion timescale by significant factors. This suggests that, either the formation of r-process elements works differently than suggested by Ning et al. (2007, NQM07), or that some essential core properties with influence on the explosion dynamics might be different from those predicted by Nomoto's progenitor model.

We report here on recent progress in understanding the birth conditions of neutron stars and the way how supernovae explode. More sophisticated numerical models have led to the discovery of new phenomena in the supernova core, for example a generic hydrodynamic instability of the stagnant supernova shock against low-mode nonradial deformation and the excitation of gravity-wave activity in the surface and core of the nascent neutron star. Both can have supportive or decisive influence on the inauguration of the explosion, the former by improving the conditions for energy deposition by neutrino heating in the postshock gas, the latter by supplying the developing blast with a flux of acoustic power that adds to the energy transfer by neutrinos. While recent two-dimensional models suggest that the neutrino-driven mechanism may be viable for stars from about 8 solar masses to at least 15 solar masses, acoustic energy input has been advocated as an alternative if neutrino heating fails. Magnetohydrodynamic effects constitute another way to trigger explosions in connection with the collapse of sufficiently rapidly rotating stellar cores, perhaps linked to the birth of magnetars. The global explosion asymmetries seen in the recent simulations offer an explanation of even the highest measured kick velocities of young neutron stars.

Based on the shell model for Gamow-Teller and the Random Phase Approximation for forbidden transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. The rates have been incorporated into state-of-the-art stellar core-collapse simulations with detailed energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%.

Advances in our understanding and the modeling of stellar core-collapse and supernova explosions over the past 15 years are reviewed, concentrating on the evolution of hydrodynamical simulations, the description of weak interactions and nuclear equation of state effects, and new insights into the nucleosynthesis occurring in the early phases of the explosion, in particular the neutrino-p process. The latter is enabled by the proton-richness of the early ejecta, which was discovered because of significant progress has been made in the treatment of neutrino transport and weak interactions. This progress has led to a new generation of sophisticated Newtonian and relativistic hydrodynamics simulations in spherical symmetry. Based on these, it is now clear that the prompt bounce-shock mechanism is not the driver of supernova explosions, and that the delayed neutrino-heating mechanism can produce explosions without the aid of multi-dimensional processes only if the progenitor star has an ONeMg core inside a very dilute He-core, i.e., has a mass in the 8--10 solar mass range. Hydrodynamic instabilities of various kinds have indeed been recognized to occur in the supernova core and to be of potential importance for the explosion. Neutrino-driven explosions, however, have been seen in two-dimensional simulations with sophisticated neutrino transport so far only when the star has a small iron core and low density in the surrounding shells as being found in stars near 10--11 solar masses. The explosion mechanism of more massive progenitors is still a puzzle. It might involve effects of three-dimensional hydrodynamics or might point to the relevance of rapid rotation and magnetohydrodynamics, or to still incompletely explored properties of neutrinos and the high-density equation of state.

In a fully relativistic approach, a RLSM description of nuclei colliding at ultra-relativistic energies can be formulated within the framework of a classical transport theory. The valence quarks of the nucleons are described through collections of classical point-like particles moving in the continuum. They are coupled to soft gluon fields which are described through the Yang Mills equations on a gauge lattice. In a first step, we focus on the range of low-$p_t$ interactions. Results from numerical model simulations of pure gluonic nucleus-nucleus collisions on SU(2) gauge lattices in 3+1 dimensions are presented. They show an effect which we call the glue burst.