The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity.

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

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.

With the same method as used previously, we investigate neutrino-driven explosions of a larger sample of blue supergiant models. The larger sample includes three new presupernova stars. The results are compared with light-curve observations of the peculiar type IIP SN 1987A. The explosions were modeled in 3D with the neutrino-hydrodynamics code PROMETHEUS-HOTB, and light-curve calculations were performed in spherical symmetry with the radiation-hydrodynamics code CRAB. Our results confirm the basic findings of the previous work: 3D neutrino-driven explosions with SN 1987A-like energies synthesize an amount of Ni-56 that is consistent with the radioactive tail of the light curve. Moreover, the models mix hydrogen inward to minimum velocities below 400 km/s as required by spectral observations. Hydrodynamic simulations with the new progenitor models, which possess smaller radii than the older ones, show much better agreement between calculated and observed light curves in the initial luminosity peak and during the first 20 days. A set of explosions with similar energies demonstrated that a high growth factor of Rayleigh-Taylor instabilities at the (C+O)/He composition interface combined with a weak interaction of fast Rayleigh-Taylor plumes, where the reverse shock occurs below the He/H interface, provides a sufficient condition for efficient outward mixing of Ni-56 into the hydrogen envelope. This condition is realized to the required extent only in one of the older stellar models, which yielded a maximum velocity of around 3000 km/s for the bulk of ejected Ni-56, but failed to reproduce the helium-core mass of 6 Msun inferred from the absolute luminosity of the presupernova star. We conclude that none of the single-star progenitor models proposed for SN 1987A to date satisfies all constraints set by observations. (Abridged)

We develop a new method, which is based on the optimal redshift weighting scheme, to extract the maximal tomographic information of baryonic acoustic oscillations (BAO) and redshift space distortions (RSD) from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 quasar (DR14Q) survey. We validate our method using the EZ mocks, and apply our pipeline to the eBOSS DR14Q sample in the redshift range of $0.8<z<2.2$. We report a joint measurement of $f\sigma_8$ and two-dimensional BAO parameters $D_{\rm A}$ and $H$ at four effective redshifts of $z_{\rm eff}=0.98, 1.23, 1.52$ and $1.94$, and provide the full data covariance matrix. Using our measurement combined with BOSS DR12, MGS and 6dFGS BAO measurements, we find that the existence of dark energy is supported by observations at a $7.4\sigma$ significance level. Combining our measurement with BOSS DR12 and Planck observations, we constrain the gravitational growth index to be $\gamma=0.580\pm0.082$, which is fully consistent with the prediction of general relativity. This paper is part of a set that analyses the eBOSS DR14 quasar sample.

We develop a practical methodology to remove modes from a galaxy survey power spectrum that are associated with systematic errors. We apply this to the BOSS CMASS sample, to see if it removes the excess power previously observed beyond the best-fit $\Lambda$CDM model on very large scales. We consider several possible sources of data contamination, and check whether they affect the number of targets that can be observed and the power spectrum measurements. We describe a general framework for how such knowledge can be transformed into template fields. Mode subtraction can then be used to remove these systematic contaminants at least as well as applying corrective weighting to the observed galaxies, but benefits from giving an unbiased power. Even after applying templates for all known systematics, we find a large-scale power excess, but this is reduced compared with that observed using the weights provided by the BOSS team. This excess is at much larger scales than the BAO scale and does not affect the main results of BOSS. However, it will be important for the measurement of a scale-dependent bias due to primordial non-Gaussianity. The excess is beyond that allowed by any simple model of non-Gaussianity matching Planck data, and is not matched in other surveys. We show that this power excess can further be reduced but is still present using "phenomenological" templates, designed to consider further potentially unknown sources of systematic contamination. As all discrepant angular modes can be removed using "phenomenological" templates, the potentially remaining contaminant acts radially.

We present a new method for solving the three-dimensional gravitational potential of a density field on the Yin-Yang grid. Our algorithm is based on a multipole decomposition and completely symmetric with respect to the two Yin-Yang grid patches. It is particularly efficient on distributed-memory machines with a large number of compute tasks, because the amount of data being explicitly communicated is minimized. All operations are performed on the original grid without the need for interpolating data onto an auxiliary spherical mesh.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We measure the growth rate and its evolution using the anisotropic clustering of the extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample, which includes $148\,659$ quasars covering the wide redshift range of $0.8 < z < 2.2$ and a sky area of $2112.90$ $\rm deg^2$. To optimise measurements we deploy a redshift-dependent weighting scheme, which allows us to avoid binning, and perform the data analysis consistently including the redshift evolution across the sample. We perform the analysis in Fourier space, and use the redshift evolving power spectrum multipoles to measure the redshift space distortion parameter $f\sigma_8$ and parameters controlling the anisotropic projection of the cosmological perturbations. We measure $f \sigma_8(z=1.52)=0.43 \pm 0.05 $ and $df\sigma_8/dz (z=1.52)= - 0.16 \pm 0.08$, consistent with the expectation for a $\Lambda$CDM cosmology as constrained by the Planck experiment.

We present an anisotropic analysis of Baryon Acoustic Oscillation (BAO) signal from the SDSS-IV extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 (DR14) quasar sample. The sample consists of 147,000 quasars distributed over a redshift range of $0.8 < z < 2.2$. We apply the redshift weights technique to the clustering of quasars in this sample and achieve a 4.6 per cent measurement of the angular distance measurement $D_M$ at $z = 2.2$ and Hubble parameter $H$ at $z=0.8$. We parameterize the distance-redshift relation, relative to a fiducial model, as a quadratic expansion. The coefficients of this expansion are used to reconstruct the distance-redshift relation and obtain distance and Hubble parameter measurements at all redshifts within the redshift range of the sample. Reporting the result at two characteristic redshifts, we determine $D_M(z=1) = 3405\pm305 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=1) = 120.7\pm 7.3 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$ and $D_M(z=2) = 5325\pm249 \ (r_{\rm d} / r_{\rm d, fid}) \ {\rm Mpc}$, $H(z=2) = 189.9\pm 32.9 \ (r_{\rm d,fid} / r_{\rm d}) \ {\rm km} \ {\rm s}^{-1}{\rm Mpc}^{-1}$. These measurements are highly correlated. We assess the outlook of BAO analysis from the final quasar sample by testing the method on a set of mocks that mimic the noise level in the final sample. We demonstrate on these mocks that redshift weighting shrinks the measurement error by over 25 per cent on average. We conclude redshift weighting can bring us closer to the cosmological goal of the final quasar sample.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 14 sample includes 80,118 Luminous Red Galaxies. By combining these galaxies with the high-redshift tail of the BOSS galaxy sample, we form a sample of LRGs at an effective redshift $z=0.72$, covering an effective volume of 0.9~Gpc$^3$. We introduce new techniques to account for spurious fluctuations caused by targeting and by redshift failures which were validated on a set of mock catalogs. This analysis is sufficient to provide a $2.6$\% measurement of spherically averaged BAO, $D_V(z=0.72) = 2353^{+63}_{-61} (r_d/r_{d,\rm{fid}}) h^{-1}$Mpc, at 2.8$\sigma$ of significance. Together with the recent quasar-based BAO measurement at $z=1.5$, and forthcoming Emission Line Galaxy-based measurements, this measurement demonstrates that eBOSS is fulfilling its remit of extending the range of redshifts covered by such measurements, laying the ground work for forthcoming surveys such as the Dark Energy Spectroscopic Survey and Euclid.

The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is one of the first of a new generation of galaxy redshift surveys that will cover a large range in redshift with sufficient resolution to measure the baryon acoustic oscillations (BAO) signal. For surveys covering a large redshift range we can no longer ignore cosmological evolution, meaning that either the redshift shells analysed have to be significantly narrower than the survey, or we have to allow for the averaging over evolving quantities. Both of these have the potential to remove signal: analysing small volumes increases the size of the Fourier window function, reducing the large-scale information, while averaging over evolving quantities can, if not performed carefully, remove differential information. It will be important to measure cosmological evolution from these surveys to explore and discriminate between models. We apply a method to optimally extract this differential information to mock catalogues designed to mimic the eBOSS quasar sample. By applying a set of weights to extract redshift space distortion measurements as a function of redshift, we demonstrate an analysis that does not invoke the problems discussed above. We show that our estimator gives unbiased constraints.

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$\sigma$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $\Lambda$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous cosmological tests to be performed with eBOSS quasar samples. We combine our result with previous, independent, BAO distance measurements to construct an updated BAO distance-ladder. Using these BAO data alone and marginalizing over the length of the standard ruler, we find $\Omega_{\Lambda} > 0$ at 6.6$\sigma$ significance when testing a $\Lambda$CDM model with free curvature.

Type II-plateau supernovae (SNe IIP) are the most numerous subclass of core-collapse SNe originating from massive stars. In the framework of the neutrino-driven explosion mechanism, we study the SN outburst properties for a red supergiant progenitor model and compare the corresponding light curves with observations of the ordinary Type IIP SN 1999em. Three-dimensional (3D) simulations of (parametrically triggered) neutrino-driven explosions are performed with the (explicit, finite-volume, Eulerian, multifluid hydrodynamics) code PROMETHEUS, using a presupernova model of a 15 Msun star as initial data. At approaching homologous expansion, the hydrodynamical and composition variables of the 3D models are mapped to a spherically symmetric configuration, and the simulations are continued with the (implicit, Lagrangian radiation-hydrodynamics) code CRAB to follow the blast-wave evolution during the SN outburst. Our 3D neutrino-driven explosion model with an explosion energy of about 0.5x10^51 erg produces Ni-56 in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail of SN 1999em. The considered presupernova model, 3D explosion simulations, and light-curve calculations can explain the basic observational features of SN 1999em, except for those connected to the presupernova structure of the outer stellar layers. Our 3D simulations show that the distribution of Ni-rich matter in velocity space is asymmetric with a strong dipole component that is consistent with the observations of SN 1999em. The monotonic luminosity decline from the plateau to the radioactive tail in ordinary SNe IIP is a manifestation of the intense turbulent mixing at the He/H composition interface.

Extremely strong magnetic fields of the order of $10^{15}\,{\rm G}$ are required to explain the properties of magnetars, the most magnetic neutron stars. Such a strong magnetic field is expected to play an important role for the dynamics of core-collapse supernovae, and in the presence of rapid rotation may power superluminous supernovae and hypernovae associated to long gamma-ray bursts. The origin of these strong magnetic fields remains, however, obscure and most likely requires an amplification over many orders of magnitude in the protoneutron star. One of the most promising agents is the magnetorotational instability (MRI), which can in principle amplify exponentially fast a weak initial magnetic field to a dynamically relevant strength. We describe our current understanding of the MRI in protoneutron stars and show recent results on its dependence on physical conditions specific to protoneutron stars such as neutrino radiation, strong buoyancy effects and large magnetic Prandtl number.

The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive $^{44}$Ti with high-energy satellite telescopes (Nuclear Spectroscopic Telescope Array [NuSTAR], INTEGRAL) have measured gamma-ray line details, which provide direct evidence of large-scale explosion asymmetries in SN 1987A and in Cassiopeia A (Cas A) even by mapping of the spatial brightness distribution (NuSTAR). Here we discuss a 3D simulation of a neutrino-driven explosion, using a parameterized neutrino engine, whose $^{44}$Ti distribution is mostly concentrated in one hemisphere pointing opposite to the neutron star (NS) kick velocity. Both exhibit intriguing resemblance to the observed morphology of the Cas A remnant, although neither the progenitor nor the explosion was fine-tuned for a perfect match. Our results demonstrate that the asymmetries observed in this remnant can, in principle, be accounted for by a neutrino-driven explosion, and that the high $^{44}$Ti abundance in Cas A may be explained without invoking rapid rotation or a jet-driven explosion, because neutrino-driven explosions generically eject large amounts of high-entropy matter. The recoil acceleration of the NS is connected to mass ejection asymmetries and is opposite to the direction of the stronger explosion, fully compatible with the gravitational tugboat mechanism. Our results also imply that Cas A and SN 1987A could possess similarly "one-sided" Ti and Fe asymmetries, with the difference that Cas A is viewed from a direction with large inclination angle to the NS motion, whereas the NS in SN 1987A should have a dominant velocity component pointing toward us.

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.

Galaxy clustering data from current and upcoming large scale structure surveys can provide strong constraints on primordial non-Gaussianity through the scale-dependent halo bias. To fully exploit the information from galaxy surveys, optimal analysis methods need to be developed and applied to the data. Since the halo bias is sensitive to local non-Gaussianity predominately at large scales, the volume of a given survey is crucial. Consequently, for such analyses we do not want to split into redshift bins, which would lead to information loss due to edge effects, but instead analyse the full sample. We present an optimal technique to directly constrain local non-Gaussianity parametrised by $f_\mathrm{NL}^\mathrm{loc}$, from galaxy clustering by applying redshift weights to the galaxies. We derive a set of weights to optimally measure the amplitude of local non-Gaussianity, $f_\mathrm{NL}^\mathrm{loc}$, discuss the redshift weighted power spectrum estimators, outline the implementation procedure and test our weighting scheme against Lognormal catalogs for two different surveys: the quasar sample of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and the emission line galaxy sample of the Dark Energy Spectroscopic Instrument (DESI) survey. We find an improvement of 30 percent for eBOSS and 6 percent for DESI compared to the standard Feldman, Kaiser $\&$ Peacock weights, although these predictions are sensitive to the bias model assumed.

We use baryon acoustic oscillation and redshift space distortion from the completed Baryon Oscillation Spectroscopic Survey, corresponding to data release 12 of the Sloan Digital Sky Survey, combined sample analysis in combination with cosmic microwave background, supernova and redshift space distortion measurements from additional spectroscopic surveys to test deviations from general relativity. We present constraints on several phenomenological models of modified gravity: First, we parametrise the growth of structure using the growth index $\gamma$, finding $\gamma=0.566\pm0.058$ (68\% C.L.). Second, we modify the relation of the two Newtonian potentials by introducing two additional parameters, $G_M$ and $G_L$. In this approach, $G_M$ refers to modifications of the growth of structure whereas $G_L$ to modification of the lensing potential. We consider a power law to model the redshift dependency of $G_M$ and $G_L$ as well as binning in redshift space, introducing four additional degrees of freedom, $G_M(z<0.5)$, $G_M(z>0.5)$, $G_L(z<0.5)$, $G_L(z>0.5)$. At 68\% C.L. we measure $G_M=0.980 \pm 0.096$ and $G_L=1.082 \pm 0.060$ for a linear model, $G_M= 1.01\pm 0.36$ and $G_L=1.31 \pm 0.19$ for a cubic model as well as $G_M(z<0.5)=1.26 \pm 0.32$, $G_M(z>0.5)=0.986\pm0.022$, $G_L(z<0.5)=1.067\pm0.058$ and $G_L(z>0.5)=1.037\pm0.029$. Thirdly, we investigate general scalar tensor theories of gravity, finding the model to be mostly unconstrained by current data. Assuming a one-parameter $f(R)$ model we can constrain $B_0< 7.7 \times 10^{-5}$ (95\% C.L). For all models we considered we find good agreement with general relativity.

During the common envelope (CE) phase, a giant star in a binary system overflows its Roche lobe and unstable mass transfer leads to a spiral-in of the companion, resulting in a close binary system or in a merger of the stellar cores. Dynamo processes during the CE phase have been proposed as a mechanism to generate magnetic fields that are important for forming magnetic white dwarfs (MWDs) and for shaping planetary nebulae. Here, we present the first magnetohydrodynamics simulations of the dynamical spiral-in during a CE phase. We find that magnetic fields are strongly amplified in the accretion stream around the $1M_\odot$ companion as it spirals into the envelope of a $2M_\odot$ RG. This leads to field strengths of 10 to 100 kG throughout the envelope after 120 d. The magnetic field amplification is consistent with being driven by the magnetorotational instability. The field strengths reached in our simulation make the magnetic field interesting for diagnostic purposes, but they are dynamically irrelevant. They are also too small to explain the formation of the highest fields found in MWDs, but may be relevant for luminous red novae, and detecting magnetic fields in these events would support the scenario as proposed here.

The well-studied type IIP SN 1987A, produced by the explosion of a blue supergiant (BSG) star, is a touchstone for massive-star evolution, simulations of neutrino-driven explosions, and modeling of light curves and spectra. In the framework of the neutrino-driven mechanism, we study the dependence of explosion properties on the structure of four different BSGs and compare the corresponding light curves with observations of SN 1987A. We perform 3D simulations with the PROMETHEUS code until about one day and map the results to the 1D code CRAB for the light curve calculations. All of our 3D models with explosion energies compatible with SN 1987A produce 56Ni in rough agreement with the amount deduced from fitting the radioactively powered light-curve tail. One of the progenitors yields maximum velocities of ~3000 km/s for the bulk of ejected 56Ni, consistent with observations. In all of our models inward mixing of hydrogen during the 3D evolution leads to minimum H-velocities below 100 km/s, in good agreement with spectral observations. The considered BSG models, 3D explosion simulations, and light-curve calculations can thus explain basic observational features of SN 1987A. However, all progenitors have too large pre-SN radii to reproduce the narrow initial luminosity peak, and the structure of their outer layers is not suitable to match the observed light curve during the first 30-40 days. Only one stellar model has a structure of the He core and the He/H composition interface that enables sufficient outward mixing of 56Ni and inward mixing of hydrogen to produce a good match of the dome-like shape of the observed light-curve maximum. But this model falls short of the He-core mass of 6 Msun inferred from the absolute luminosity of the pre-SN star. The lack of an adequate pre-SN model for SN 1987A is a pressing challenge for the theory of massive-star evolution. (Abridged)

The magneto-rotational instability (MRI) is considered to be a promising mechanism to amplify the magnetic field in fast rotating protoneutron stars. In contrast to accretion disks, radial buoyancy driven by entropy and lepton fraction gradients is expected to have a dynamical role as important as rotation and shear. We investigate the poorly known impact of buoyancy on the non-linear phase of the MRI, by means of three dimensional numerical simulations of a local model in the equatorial plane of a protoneutron star. The use of the Boussinesq approximation allows us to utilise a shearing box model with clean shearing periodic boundary conditions, while taking into account the buoyancy driven by radial entropy and composition gradients. We find significantly stronger turbulence and magnetic fields in buoyantly unstable flows. On the other hand, buoyancy has only a limited impact on the strength of turbulence and magnetic field amplification for buoyantly stable flows in the presence of a realistic thermal diffusion. The properties of the turbulence are, however, significantly affected in the latter case. In particular, the toroidal components of the magnetic field and of the velocity become even more dominant with respect to the poloidal ones. Furthermore, we observed in the regime of stable buoyancy the formation of long lived coherent structures such as channel flows and zonal flows. Overall, our results support the ability of the MRI to amplify the magnetic field significantly even in stably stratified regions of protoneutron stars.

We present 3D simulations of core-collapse supernovae from blast-wave initiation by the neutrino-driven mechanism to shock breakout from the stellar surface, considering two 15 Msun red supergiants (RSG) and two blue supergiants (BSG) of 15 Msun and 20 Msun. We demonstrate that the metal-rich ejecta in homologous expansion still carry fingerprints of asymmetries at the beginning of the explosion, but the final metal distribution is massively affected by the detailed progenitor structure. The most extended and fastest metal fingers and clumps are correlated with the biggest and fastest-rising plumes of neutrino-heated matter, because these plumes most effectively seed the growth of Rayleigh-Taylor (RT) instabilities at the C+O/He and He/H composition-shell interfaces after the passage of the SN shock. The extent of radial mixing, global asymmetry of the metal-rich ejecta, RT-induced fragmentation of initial plumes to smaller-scale fingers, and maximal Ni and minimal H velocities do not only depend on the initial asphericity and explosion energy (which determine the shock and initial Ni velocities) but also on the density profiles and widths of C+O core and He shell and on the density gradient at the He/H transition, which lead to unsteady shock propagation and the formation of reverse shocks. Both RSG explosions retain a great global metal asymmetry with pronounced clumpiness and substructure, deep penetration of Ni fingers into the H-envelope (with maximum velocities of 4000-5000 km/s for an explosion energy around 1.5 bethe) and efficient inward H-mixing. While the 15 Msun BSG shares these properties (maximum Ni speeds up to ~3500 km/s), the 20 Msun BSG develops a much more roundish geometry without pronounced metal fingers (maximum Ni velocities only ~2200 km/s) because of reverse-shock deceleration and insufficient time for strong RT growth and fragmentation at the He/H interface.

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

We present the mean pairwise momentum of clusters, as observed through the kinematic Sunvaev-Zel'dovich (kSZ) effect, as a novel probe of massive neutrinos. We find that kSZ measurements with current and upcoming surveys will provide complementary constraints on the sum of neutrino masses from large scale structure and will improve on Planck satellite measurements of the primordial cosmic microwave background (CMB) and CMB lensing. Central to the constraints is a distinctive scale dependency of the kSZ neutrino signature on the mean pairwise momentum of clusters that we do not expect to be mirrored in systematic effects that change the overall amplitude of the signal, like the cluster optical depth. Assuming a minimal $\Lambda$CDM cosmology including massive neutrinos with Planck primordial CMB priors combined with conservative kSZ specifications, we forecast $68%$ upper limits on the neutrino mass sum of $290$ meV, $220$ meV, $96$ meV for "Stage II" (ACTPol + BOSS), "Stage III" (Advanced ACTPol + BOSS), and "Stage IV" (CMB-S4 + DESI) surveys respectively, compared to the Planck alone forecast of 540 meV. These forecasts include the ability to simultaneously constrain the neutrino mass sum and the mass-averaged optical depth of the cluster sample in each redshift bin. If the averaged optical depth of clusters can be measured with few percent accuracy and a lower limiting mass is assumed, the projected kSZ constraints improve further to $120$ meV, $90$ meV and $33$ meV (Stage II, III and IV). These forecasts represent a conservative estimate of neutrino constraints using cross-correlations of arcminute-resolution CMB measurements and spectroscopic galaxy surveys. More information relevant for neutrino constraints is available from these surveys, such as galaxy clustering, weak lensing, and CMB temperature and polarization.

We calculate the constraints on dark energy and cosmic modifications to gravity achievable with upcoming cosmic microwave background (CMB) surveys sensitive to the Sunyaev-Zeldovich (SZ) effects. The analysis focuses on using the mean pairwise velocity of clusters as observed through the kinematic SZ effect (kSZ), an approach based on the same methods used for the first detection of the kSZ effect, and includes a detailed derivation and discussion of this statistic's covariance under a variety of different survey assumptions. The potential of current, Stage II, and upcoming, Stage III and Stage IV, CMB observations are considered, in combination with contemporaneous spectroscopic and photometric galaxy observations. A detailed assessment is made of the sensitivity to the assumed statistical and systematic uncertainties in the optical depth determination, the magnitude and uncertainty in the minimum detectable mass, and the importance of pairwise velocity correlations at small separations, where non-linear effects can start to arise. In combination with Stage III constraints on the expansion history, such as those projected by the Dark Energy Task Force, we forecast 5\% and 2\% for fractional errors on the growth factor, $\gamma$, for Stage III and Stage IV surveys respectively, and 2\% constraints on the growth rate, $f_g$, for a Stage IV survey for $0.2<z<0.6$. The results suggest that kSZ measurements of cluster peculiar velocities, obtained from cross-correlation with upcoming spectroscopic galaxy surveys, could provide robust tests of dark energy and theories of gravity on cosmic scales.

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

We consider cosmological constraints arising from the background expansion history on the ef- fective field theory of cosmic acceleration, a theoretical framework that allows for a unified way to classify both models of dark energy and modified gravity within the linear regime. In the Einstein frame, the most general action for the background can be written in terms of a canonical scalar field which is non-minimally coupled to matter. The leading corrections to the action are expressible through a quartic kinetic term, and scalar couplings to a Gauss-Bonnet curvature term and the Einstein tensor. We determine the implications of the terms in this general action for the predicted expansion history in the context of dynamical attractors. We find that each modifies the matter dominated and/or accelerative eras in ways that allow us to place cosmological constraints on them. We present current constraints on the effective action using the latest Type Ia supernovae, Cosmic Microwave Background, and Baryonic Acoustic Oscillation data. This includes finding that the scalar field EFT with a coupled Gauss-Bonnet term and the data are significantly discrepant.

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

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.

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

Mixing processes in stars driven by composition gradients as a result of the Rayleigh-Taylor instability are not anticipated. They are supported only by hydrodynamic studies of stellar convection. We find that such mixing occurs below the bottom edge of convection zones in our multidimensional hydrodynamic shell convection models. It operates at interfaces created by off-center nuclear burning, where less dense gas with higher mean molecular weight is located above denser gas with a lower mean molecular weight. We discuss the mixing under various conditions with hydrodynamic convection models based on stellar evolutionary calculations of the core helium flash in a 1.25 Msun star, the core carbon flash in a 9.3 Msun star, and of oxygen burning shell in a star with a mass of 23 Msun. We simulate the hydrodynamic behavior of shell convection during various phases of stellar evolution with the Eulerian hydrodynamics code HERAKLES in two and three spatial dimensions. Initial models for this purpose are obtained by state-of-the-art stellar evolutionary codes, namely GARSTEC, STAREVOL, and TYCHO for the core helium flash, core carbon flash, and oxygen shell burning, respectively. Most of our analysis is performed for two-dimensional hydrodynamic models of shell convection during the core helium flash at its peak covering approximately 250 convective turnover timescales. The mixing manifests itself in the form of overdense and cold fingers enriched with matter of higher mean molecular weight, originating from density fluctuations at the lower boundary of the convective shell, and "shooting" down into the core. They result from the Rayleigh-Taylor instability at the lower convection zone boundary due to a negative mean molecular weight gradient. They do not appear when the mean molecular weight gradient is positive.

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

Degenerate ignition of helium in low-mass stars at the end of the red giant branch phase leads to dynamic convection in their helium cores. One-dimensional (1D) stellar modeling of this intrinsically multi-dimensional dynamic event is likely to be inadequate. Previous hydrodynamic simulations imply that the single convection zone in the helium core of metal-rich Pop I stars grows during the flash on a dynamic timescale. This may lead to hydrogen injection into the core, and a double convection zone structure as known from one-dimensional core helium flash simulations of low-mass Pop III stars. We perform hydrodynamic simulations of the core helium flash in two and three dimensions to better constrain the nature of these events. To this end we study the hydrodynamics of convection within the helium cores of a 1.25 \Msun metal-rich Pop I star (Z=0.02), and a 0.85 \Msun metal-free Pop III star (Z=0) near the peak of the flash. These models possess single and double convection zones, respectively. We use 1D stellar models of the core helium flash computed with state-of-the-art stellar evolution codes as initial models for our multidimensional hydrodynamic study, and simulate the evolution of these models with the Riemann solver based hydrodynamics code Herakles which integrates the Euler equations coupled with source terms corresponding to gravity and nuclear burning. The hydrodynamic simulation of the Pop I model involving a single convection zone covers 27 hours of stellar evolution, while the first hydrodynamic simulations of a double convection zone, in the Pop III model, span 1.8 hours of stellar life. We find differences between the predictions of mixing length theory and our hydrodynamic simulations. The simulation of the single convection zone in the Pop I model shows a strong growth of the size of the convection zone due to turbulent entrainment. Hence we predict that for the Pop I model a hydrogen injection phase (i.e. hydrogen injection into the helium core) will commence after about 23 days, which should eventually lead to a double convection zone structure known from 1D stellar modeling of low-mass Pop III stars. Our two and three-dimensional hydrodynamic simulations of the double (Pop III) convection zone model show that the velocity field in the convection zones is different from that predicted by stellar evolutionary calculations. The simulations suggest that the double convection zone decays quickly, the flow eventually being dominated by internal gravity waves.

We describe and discuss hydrodynamic simulations of the core helium flash using an initial model of a 1.25 M_sol star with a metallicity of 0.02 near at its peak. Past research concerned with the dynamics of the core helium flash is inconclusive. Its results range from a confirmation of the standard picture, where the star remains in hydrostatic equilibrium during the flash (Deupree 1996), to a disruption or a significant mass loss of the star (Edwards 1969; Cole & Deupree 1980). However, the most recent multidimensional hydrodynamic study (Dearborn 2006) suggests a quiescent behavior of the core helium flash and seems to rule out an explosive scenario. Here we present partial results of a new comprehensive study of the core helium flash, which seem to confirm this qualitative behavior and give a better insight into operation of the convection zone powered by helium burning during the flash. The hydrodynamic evolution is followed on a computational grid in spherical coordinates using our new version of the multi-dimensional hydrodynamic code HERAKLES, which is based on a direct Eulerian implementation of the piecewise parabolic method.

We present 2D hydrodynamic simulations of the long-time accretion phase of a 15 solar mass star after core bounce and before the launch of a supernova explosion. Our simulations are performed with the Prometheus-Vertex code, employing multi-flavor, energy-dependent neutrino transport and an effective relativistic gravitational potential. Testing the influence of a stiff and a soft equation of state for hot neutron star matter, we find that the non-radial mass motions in the supernova core due to the standing accretion shock instability (SASI) and convection impose a time variability on the neutrino and gravitational-wave signals. These variations have larger amplitudes as well as higher frequencies in the case of a more compact nascent neutron star. After the prompt shock-breakout burst of electron neutrinos, a more compact accreting remnant radiates neutrinos with higher luminosities and larger mean energies. The observable neutrino emission in the direction of SASI shock oscillations exhibits a modulation of several 10% in the luminosities and ~1 MeV in the mean energies with most power at typical SASI frequencies of 20-100 Hz. At times later than 50-100 ms after bounce the gravitational-wave amplitude is dominated by the growing low-frequency (<200 Hz) signal associated with anisotropic neutrino emission. A high-frequency wave signal is caused by nonradial gas flows in the outer neutron star layers, which are stirred by anisotropic accretion from the SASI and convective regions. The gravitational-wave power then peaks at about 300-800 Hz with distinctively higher spectral frequencies originating from the more compact and more rapidly contracting neutron star. The detectability of the SASI effects in the neutrino and gravitational-wave signals is briefly discussed. (abridged)

We study turbulent convection during the core helium flash close to its peak by comparing the results of two and three-dimensional hydrodynamic simulations. We use a multidimensional Eulerian hydrodynamics code based on state-of-the-art numerical techniques to simulate the evolution of the helium core of a $1.25 M_{\odot}$ Pop I star. Our three-dimensional hydrodynamic simulations of the evolution of a star during the peak of the core helium flash do not show any explosive behavior. The convective flow patterns developing in the three-dimensional models are structurally different from those of the corresponding two-dimensional models, and the typical convective velocities are smaller than those found in their two-dimensional counterparts. Three-dimensional models also tend to agree better with the predictions of mixing length theory. Our hydrodynamic simulations show the presence of turbulent entrainment that results in a growth of the convection zone on a dynamic time scale. Contrary to mixing length theory, the outer part of the convection zone is characterized by a sub-adiabatic temperature gradient.

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 investigate the hydrodynamics of the core helium flash near its peak. Past research concerned with the dynamics of this event is inconclusive. However, the most recent multidimensional hydrodynamic studies suggest a quiescent behavior and seem to rule out an explosive scenario. Previous work indicated, that depending on initial conditions, employed turbulence models, grid resolution, and dimensionality of the simulation, the core helium flash leads either to the disruption of a low-mass star or to a quiescent quasi-hydrostatic evolution. We try to clarify this issue by simulating the evolution with advanced numerical methods and detailed microphysics. Assuming spherical or axial symmetry, we simulate the evolution of the helium core of a $1.25 M_{\odot}$ star with a metallicity Z=0.02 during the core helium flash at its peak with a grid-based hydrodynamics code. We find that the core helium flash neither rips the star apart, nor that it significantly alters its structure, as convection plays a crucial role in keeping the star in hydrostatic equilibrium. In addition, our simulations show the presence of overshooting, which implies new predictions concerning mixing of chemical species in red giants.

We perform general relativistic simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state as well as an approximate description of deleptonization. We show that for a wide variety of rotation rates and profiles the gravitational wave burst signals from the core bounce are of a generic type, known as Type I in the literature. In our systematic study, using both general relativity and Newtonian gravity, we identify and individually quantify the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the equation of state, and deleptonization. Such a generic type of signal templates will likely facilitate a more efficient search in current and future gravitational wave detectors of both interferometric and resonant type.

We present detailed results from performing general relativistic (GR) simulations of stellar core collapse to a proto-neutron star, using a microphysical equation of state (EoS) as well as an approximate description of deleptonization during the collapse phase. We show that for a wide variety of rotation rates and profiles the gravitational wave (GW) burst signals from the core bounce are of a generic type, already known as Type I in the literature. In addition, for most models the characteristic frequency of the GW burst signal lies in a narrow range around approximately 718 Hz. In our systematic study, using both GR and Newtonian gravity, we identify, individually quantify, and discuss in detail the micro- and macrophysical mechanisms leading to this result, i.e. the effects of rotation, the EoS, and deleptonization. We also discuss the detectability prospects of such GW burst signals by GW detectors, and infer that such a generic type of signal templates will likely facilitate a more efficient search in current and future detectors of both interferometric and resonant type.

We investigate the role of rotational instabilities in the context of black hole formation in relativistic stars. In addition to the standard scenario - an axially symmetric dynamical instability forming a horizon at the star's center - the recently found low-$T/|W|$ instabilities are shown to lead to fragmentation and off-center horizon formation in differentially rotating stars. This process might be an alternative pathway to produce SMBHs from supermassive stars with inefficient angular momentum transport.

(Abridged): We study the collision of magnetized irregularities (shells) in relativistic outflows in order to explain the origin of the generic phenomenology observed in the non-thermal emission of both blazars and gamma-ray bursts. We focus on the influence of the magnetic field on the collision dynamics, and we further investigate how the properties of the observed radiation depend on the strength of the initial magnetic field and on the initial internal energy density of the flow. The collisions of magnetized shells and the radiation resulting from these collisions are calculated using the 1D relativistic magnetohydrodynamics code MRGENESIS. The interaction of the shells with the external medium prior to their collision is also determined using an exact solver for the corresponding 1D relativistic magnetohydrodynamic Riemann problem. Our simulations show that two magnetization parameters - the ratio of magnetic energy density and thermal energy density, \alpha_B, and the ratio of magnetic energy density and mass-energy density, \sigma - play an important role in the pre-collision phase, while the dynamics of the collision and the properties of the light curves depend mostly on the magnetization parameter \sigma. The interaction of the shells with the external medium changes the flow properties at their edges prior to the collision. For sufficiently dense shells moving at large Lorentz factors (\simgt 25) these properties depend only on the magnetization parameter \sigma. Internal shocks in GRBs may reach maximum efficiencies of conversion of kinetic into thermal energy between 6% and 10%, while in case of blazars, the maximum efficiencies are \sim 2%.

We investigate the energy-momentum deposition due to neutrino-antineutrino annihilation in the vicinity of axisymmetric, accreting black holes (BHs) by numerically ray-tracing neutrino trajectories in a Kerr space-time. Hyperaccreting stellar-mass BHs are widely considered as energy sources that can drive ultrarelativistic outflows with the potential to produce gamma-ray bursts. In contrast to earlier works, we provide an extensive and detailed parameter study of the influence of general relativistic (GR) effects and of different neutrinosphere geometries. These include idealized thin disks, tori, and spheres, or are constructed as non-selfgravitating equilibrium matter distributions for varied BH rotation. Considering isothermal neutrinospheres with the same temperature and surface area, we confirm previous results that compared to Newtonian calculations, GR effects increase the annihilation rate measured by an observer at infinity by a factor of 2 when the neutrinosphere is a disk. However, in case of a torus and a sphere the influence of GR effects is globally only ~25%, although locally it can be significantly larger. Independent of whether GR effects are included, disks yield the highest energy deposition rates, followed by tori and, with the lowest rates, spheres. For disks and tori, increasing the angular momentum of the BH from 0 to 1 enhances the energy deposition rate measured by an observer at infinity by roughly a factor of 2 due to the shrinking inner radius of the neutrinosphere. (abridged)

In a recent publication, we have demonstrated that differentially rotating stars admit new channels of black hole formation via fragmentation instabilities. Since a higher order instability of this kind could potentially transform a differentially rotating supermassive star into a multiple black hole system embedded in a massive accretion disk, we investigate the dependence of the instability on parameters of the equilibrium model. We find that many of the models constructed exhibit non-axisymmetric instabilities with corotation points, even for low values of T/|W|, which lead to a fission of the stars into one, two or three fragments, depending on the initial perturbation. At least in the models selected here, an m=1 mode becomes unstable at lower values of T/|W|, which would seem to favor a scenario where one black hole with a massive accretion disk forms. In this case, we have gained evidence that low values of compactness of the initial model can lead to a stabilization of the resulting fragment, thus preventing black hole formation in this scenario.

We continue our investigations of the magnetorotational collapse of stellar cores discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV potential used in our simulations catches several features of fully relativistic simulations quite well. It is able to correctly reproduce the behavior of models which show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. If this is not the case, the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential allowing for a more efficient amplification of the magnetic field by differential rotation. Sufficiently strong magnetic fields brake down the core's rotation and trigger a contraction phase to higher densities. Several models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce does not suffice to eject any mass, and the core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Instead it belongs to the signal type IV introduced by us in the first paper of this series. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and > 100 kHz denotes the generation of a jet-like outflow. [Abstract abbreviated]

By means of two-dimensional (2D) simulations we study hydrodynamic instabilities during the first seconds of neutrino-driven supernova explosions, using a PPM hydrodynamics code, supplemented with a gray, non-equilibrium approximation of radial neutrino transport. We consider three 15 solar mass progenitors with different structures and one rotating model, in which we replace the dense core of the newly formed neutron star (NS) by a contracting inner grid boundary, and trigger neutrino-driven explosions by systematically varying the neutrino fluxes emitted at this boundary. Confirming more idealized studies as well as supernova simulations with spectral transport, we find that random seed perturbations can grow by hydrodynamic instabilities to a globally asymmetric mass distribution, leading to a dominance of dipole (l=1) and quadrupole (l=2) modes in the explosion ejecta. Anisotropic gravitational and hydrodynamic forces are found to accelerate the NS on a timescale of 2-3 seconds. Since the explosion anisotropies develop chaotically, the magnitude of the corresponding kick varies stochastically in response to small differences in the fluid flow. Our more than 70 models separate into two groups, one with high and the other with low NS velocities and accelerations after 1s of post-bounce evolution, depending on whether the l=1 mode is dominant in the ejecta or not. This leads to a bimodality of the distribution when the NS velocities are extrapolated to their terminal values. The fast group has an average velocity of about 500 km/s and peak values in excess of 1000 km/s. Establishing a link to the measured distribution of pulsar velocities, however, requires a much larger set of calculations and ultimately 3D modeling. (abridged)

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l=2 and l=1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**51 erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event. (abbreviated)

We have performed a comprehensive parameter study of the collapse of rotating, strongly magnetized stellar cores in axisymmetry to determine their gravitational wave signature based on the Einstein quadrupole formula. We use a Newtonian explicit magnetohydrodynamic Eulerian code based on the relaxing-TVD method for the solution of the ideal MHD equations, and apply the constraint-transport method to guarantee a divergence--free evolution of the magnetic field. We neglect effects due to neutrino transport and employ a simplified equation of state. The pre--collapse initial models are polytropes in rotational equilibrium with a prescribed degree of differential rotation and rotational energy (~ 1 % of the gravitational energy). The initial magnetic fields are purely poloidal the field strength ranging from 10^10 G to 10^13 G. The evolution of the core, whose collapse is initiated by reducing the gas pressure by a prescribed amount, is followed until a few ten milliseconds past core bounce. The initial magnetic fields are amplified mainly by the differential rotation of the core giving rise to a strong toroidal field component. The poloidal field component grows by compression during collapse, but does not change significantly after core bounce if (abbreviated)

Based on new developments in the understanding of supernovae (SNe) as gravitational-wave (GW) sources we estimate the GW background from all cosmic SNe. For a broad range of frequencies around 1 Hz, this background is crudely comparable to the GW background expected from standard inflationary models. While our estimate remains uncertain within several orders of magnitude, the SN GW background may become detectable by second-generation space-based interferometers such as the proposed Big-Bang Observatory (BBO). By the same token, the SN GWs may become a foreground for searches of the inflationary GWs, in particular for sub-Hz frequencies where the SN background is Gaussian and where the BBO will be most sensitive. SN simulations lasting far beyond the usual cutoff of about 1 second are needed for more robust predictions in the sub-Hz frequency band. An even larger GW background can arise from a hypothetical early population of massive stars, although their GW source strength as well as their abundance are currently poorly understood.

We investigate the possibility to approximate relativistic effects in hydrodynamical simulations of stellar core collapse and post-bounce evolution by using a modified gravitational potential in an otherwise standard Newtonian hydrodynamic code. Different modifications of the effective relativistic potential introduced by Rampp & Janka (2002) are discussed. Corresponding hydrostatic solutions are compared with solutions of the TOV equations, and hydrodynamic simulations with two different codes are compared with fully relativistic results. One code is applied for one- and two-dimensional calculations with a simple equation of state, and employs either the modified effective relativistic potential in a Newtonian framework or solves the general relativistic field equations under the assumption of the conformal flatness condition (CFC) for the three-metric. The second code allows for full-scale supernova runs including a microphysical equation of state and neutrino transport based on the solution of the Boltzmann equation and its moments equations. We present prescriptions for the effective relativistic potential for self-gravitating fluids to be used in Newtonian codes, which produce excellent agreement with fully relativistic solutions in spherical symmetry, leading to significant improvements compared to previously published approximations. Moreover, they also approximate qualitatively well relativistic solutions for models with rotation.

We present results of a detailed numerical study and theoretical analysis of the dynamics of internal shocks in relativistic jets and the non-thermal flares associated with these shocks. In our model internal shocks result from collisions of density inhomogeneities (shells) in relativistic jet flows. We find that the merged shell resulting from the inelastic collision of shells has a complicated internal structure due to the non-linear dynamics of the interaction. Furthermore, the instantaneous efficiency for converting kinetic energy into thermal energy is found to be almost twice as high as theoretically expected during the period of significant emission. The Lorentz factors of the internal shocks are correlated with the initial inertial masses of the shells. Because of the complexity of the non-linear evolution the merged shell becomes very inhomogeneous and simple one-zone models are inadequate to extract physical parameters of the emitting region from the resulting light curves. In order to improve on these one-zone approximations, we propose a novel way of analyzing the space-time properties of the emission. Based on these properties we construct an analytic model of non-thermal flares which can be used to constrain some (unobservable) physical parameters of the internal shocks. These are the ratio of the Lorentz factors between the forward and the reverse shock (caused by the shell collision), and the shell crossing times of these shocks. The analytic model is validated by applying it to the synthetic light curves computed from our models. It can equally well be applied to observations.

We present a new 3D general relativistic hydrodynamics code for simulations of stellar core collapse to a neutron star, as well as pulsations and instabilities of rotating relativistic stars. It uses spectral methods for solving the metric equations, assuming the conformal flatness approximation for the three-metric. The matter equations are solved by high-resolution shock-capturing schemes. We demonstrate that the combination of a finite difference grid and a spectral grid can be successfully accomplished. This "Mariage des Maillages" (French for grid wedding) approach results in high accuracy of the metric solver and allows for fully 3D applications using computationally affordable resources, and ensures long term numerical stability of the evolution. We compare our new approach to two other, finite difference based, methods to solve the metric equations. A variety of tests in 2D and 3D is presented, involving highly perturbed neutron star spacetimes and (axisymmetric) stellar core collapse, demonstrating the ability to handle spacetimes with and without symmetries in strong gravity. These tests are also employed to assess gravitational waveform extraction, which is based on the quadrupole formula.

We investigate new paths to black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to follow the evolution to the formation of an apparent horizon centered on the fragment. This is the first study of the one-armed instability in full general relativity.

Core collapse supernovae are a promising source of detectable gravitational waves. Most of the existing (multidimensional) numerical simulations of core collapse in general relativity have been done using approximations of the Einstein field equations. As recently shown by Dimmelmeier et al (2002a,b), one of the most interesting such approximation is the so-called conformal flatness condition (CFC) of Isenberg, Wilson and Mathews. Building on this previous work we present here new results from numerical simulations of relativistic rotational core collapse in axisymmetry, aiming at improving the dynamics and the gravitational waveforms. The computer code used for these simulations evolves the coupled system of metric and fluid equations using the 3+1 formalism, specialized to a new framework for the gravitational field equations which we call CFC+. In this approach we add new degrees of freedom to the original CFC equations, which extend them by terms of second post-Newtonian order. The corrections for CFC+ are computed solving a system of elliptic linear equations. The new formalism is assessed with time evolutions of both rotating neutron stars in equilibrium and gravitational core collapse of rotating polytropes. Gravitational wave signals for a comprehensive sample of collapse models are extracted using either the quadrupole formula or directly from the metric. We discuss our results on the dynamics and the gravitational wave emission through a detailed comparison between CFC and CFC+ simulations. The main conclusion is that, for the neutron star spacetimes analyzed in the present work, no significant differences are found among CFC, CFC+, and full general relativity, which highlights the suitability of the former.

Two- and three-dimensional simulations demonstrate that hydrodynamic instabilities can lead to low-mode (l=1,2) asymmetries of the fluid flow in the neutrino-heated layer behind the supernova shock. This provides a natural explanation for aspherical mass ejection and for pulsar recoil velocities even in excess of 1000 km/s. We propose that the bimodality of the pulsar velocity distribution might be a consequence of a dominant l=1 mode in case of the fast component, while higher-mode anisotropy characterizes the postshock flow and SN ejecta during the birth of the slow neutron stars. We argue that the observed large asymmetries of supernovae and the measured high velocities of young pulsars therefore do not imply rapid rotation of the iron core of the progenitor star, nor do they require strong magnetic fields to play a crucial role in the explosion. Anisotropic neutrino emission from accretion contributes to the neutron star acceleration on a minor level, and pulsar kicks do not make a good case for non-standard neutrino physics in the nascent neutron star.

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

We present the results of relativistic hydrodynamic simulations of the collision of two dense shells in a uniform external medium, as envisaged in the internal shock model for BL Lac jets. The non-thermal radiation produced by highly energetic electrons injected at the relativistic shocks is computed following their temporal and spatial evolution. The acceleration of electrons at the relativistic shocks is parametrized using two different models and the corresponding X-ray light curves are computed. We find that the interaction time scale of the two shells is influenced by an interaction with the external medium. For the chosen parameter sets, the efficiency of the collision in converting dissipated kinetic energy into the observed X-ray radiation is of the order of one percent.

Assuming that the neutrino luminosity from the neutron star core is sufficiently high to drive supernova explosions by the neutrino-heating mechanism, we show that low-mode (l = 1, 2) convection can develop from random seed perturbations behind the shock. A slow onset of the explosion is crucial, requiring the core luminosity to vary slowly with time, in contrast to the burst-like exponential decay assumed in previous work. Gravitational and hydrodynamic forces by the globally asymmetric supernova ejecta were found to accelerate the remnant neutron star on a timescale of more than a second to velocities above 500 km/s, in agreement with observed pulsar proper motions.

We have computed the gravitational wave signal from supernova core collapse using the presently most realistic input physics available. We start from state-of-the-art progenitor models of rotating and non-rotating massive stars, and simulate the dynamics of their core collapse by integrating the equations of axisymmetric hydrodynamics together with the Boltzmann equation for the neutrino transport including an elaborate description of neutrino interactions, and a realistic equation of state. We compute the quadrupole wave amplitudes, the Fourier wave spectra, the amount of energy radiated in form of gravitational waves, and the S/N ratios for the LIGO and the tuned Advanced LIGO interferometers resulting both from non-radial mass motion and anisotropic neutrino emission. The simulations demonstrate that the dominant contribution to the gravitational wave signal is produced by neutrino-driven convection behind the supernova shock. For stellar cores rotating at the extreme of current stellar evolution predictions, the core-bounce signal is detectable with advanced LIGO up to a distance of 5kpc, whereas the signal from post-shock convection is observable up to a distance of about 100kpc. If the core is non-rotating its gravitational wave emission can be measured up to a distance of 15kpc, while the signal from the Ledoux convection in the deleptonizing, nascent neutron star can be detected up to a distance of 10kpc. Both kinds of signals are generically produced by convection in any core collapse supernova.

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

We investigate explosive nuclear burning in core collapse supernovae by coupling a tracer particle method to one and two-dimensional Eulerian hydrodynamic calculations. Adopting the most recent experimental and theoretical nuclear data, we compute the nucleosynthetic yields for 15 Msun stars with solar metallicity, by post-processing the temperature and density history of advected tracer particles. We compare our results to 1D calculations published in the literature.

We present the results of a three-dimensional, relativistic, hydrodynamic simulation of a precessing jet into which a compact blob of matter is injected. A comparison of synthetic radio maps computed from the hydrodynamic model, taking into account the appropriate light travel time delays, with those obtained from observations of actual superluminal sources shows that the variability of the jet emission is the result of a complex combination of phase motions, viewing angle selection effects, and non-linear interactions between perturbations and the underlying jet and/or the external medium. These results question the hydrodynamic properties inferred from observed apparent motions and radio structures, and reveal that shock-in-jet models may be overly simplistic.

We present results from axisymmetric stellar core collapse simulations in general relativity. Our hydrodynamics code has proved robust and accurate enough to allow for a detailed analysis of the global dynamics of the collapse. Contrary to traditional approaches based on the 3+1 formulation of the gravitational field equations, our framework uses a foliation based on a family of outgoing light cones, emanating from a regular center, and terminating at future null infinity. Such a coordinate system is well adapted to the study of interesting dynamical spacetimes in relativistic astrophysics such as stellar core collapse and neutron star formation. Perhaps most importantly this procedure allows for the unambiguous extraction of gravitational waves at future null infinity without any approximation, along with the commonly used quadrupole formalism for the gravitational wave extraction. Our results concerning the gravitational wave signals show noticeable disagreement when those are extracted by computing the Bondi news at future null infinity on the one hand and by using the quadrupole formula on the other hand. We have strong indication that for our setup the quadrupole formula on the null cone does not lead to physical gravitational wave signals. The Bondi gravitational wave signals extracted at infinity show typical oscillation frequencies of about 0.5 kHz.

We study the influence of the matter content of extragalactic jets on their morphology, dynamics and emission properties. For this purpose we consider jets of extremely different compositions including pure leptonic and baryonic plasmas. Our work is based on two-dimensional relativistic hydrodynamic simulations of the long-term evolution of powerful extragalactic jets propagating into a homogeneous environment. The equation of state used in the simulations accounts for an arbitrary mixture of electrons, protons and electron-positron pairs. Using the hydrodynamic models we have also computed synthetic radio maps and the thermal Bremsstrahlung X-ray emission from their cavities. Although there is a difference of about three orders of magnitude in the temperatures of the cavities inflated by the simulated jets, we find that both the morphology and the dynamic behaviour are almost independent on the assumed composition of the jets. Their evolution proceeds in two distinct epochs. During the first one multidimensional effects are unimportant and the jets propagate ballistically. The second epoch starts when the first larger vortices are produced near the jet head causing the beam cross section to increase and the jet to decelerate. The evolution of the cocoon and cavity is in agreement with a simple theoretical model. The beam velocities are relativistic ($\Gamma \simeq 4$) at kiloparsec scales supporting the idea that the X-ray emission of several extragalactic jets may be due to relativistically boosted CMB photons. The radio emission of all models is dominated by the contribution of the hot spots. All models exhibit a depression in the X-rays surface brightness of the cavity interior in agreement with recent observations.

We have performed hydrodynamic simulations of relativistic rotational supernova core collapse in axisymmetry and have computed the gravitational radiation emitted by such an event. Details of the methodology and of the numerical code have been given in an accompanying paper. We have simulated the evolution of 26 models in both Newtonian and relativistic gravity. Our simulations show that the three different types of rotational supernova core collapse and gravitational waveforms identified in previous Newtonian simulations (regular collapse, multiple bounce collapse, and rapid collapse) are also present in relativistic gravity. However, rotational core collapse with multiple bounces is only possible in a much narrower parameter range in relativistic gravity. The relativistic models cover almost the same range of gravitational wave amplitudes and frequencies as the corresponding Newtonian ones. For a given model, relativistic gravity can cause a large increase of the characteristic signal frequency of up to a factor of five, which may have important consequences for the signal detection. The gravitational wave signals obtained in our study are within the sensitivity range of the first generation laser interferometer detectors if the source is located within the Local Group.

We describe an axisymmetric general relativistic code for rotational core collapse. The code evolves the coupled system of metric and fluid equations using the ADM 3+1 formalism and a conformally flat metric approximation of the Einstein equations. The relativistic hydrodynamics equations are formulated as a first-order flux-conservative hyperbolic system and are integrated using high-resolution shock-capturing schemes based on Riemann solvers. We assess the quality of the conformally flat metric approximation for relativistic core collapse and present a comprehensive set of tests which the code successfully passed. The tests include relativistic shock tubes, the preservation of the rotation profile and of the equilibrium of rapidly and differentially rotating neutron stars (approximated as rotating polytropes), spherical relativistic core collapse, and the conservation of rest-mass and angular momentum in dynamic spacetimes. The application of the code to relativistic rotational core collapse, with emphasis on the gravitational waveform signature, is presented in an accompanying paper.

We present new numerical algorithms for the coupled Einstein-perfect fluid system in axisymmetry. Our framework uses a foliation based on a family of light cones, emanating from a regular center, and terminating at future null infinity. This coordinate system is well adapted to the study of the dynamical spacetimes associated with isolated relativistic compact objects such as neutron stars. In particular, the approach allows the unambiguous extraction of gravitational waves at future null infinity and avoids spurious outer boundary reflections. The code can accurately maintain long-term stability of polytropic equilibrium models of relativistic stars. We demonstrate global energy conservation in a strongly perturbed neutron star spacetime, for which the total energy radiated away by gravitational waves corresponds to a significant fraction of the Bondi mass. As a first application we present results in the study of pulsations of axisymmetric relativistic stars, extracting the frequencies of the different fluid modes in fully relativistic evolutions of the Einstein-perfect fluid system and making a first comparison between the gravitational news function and the predicted wave using the approximations of the quadrupole formula.

We use hydrodynamic simulations with adaptive grid refinement to study the dependence of hot gas flows in X-ray luminous giant elliptical galaxies on the efficiency of heat supply to the gas. We consider a number of potential heating mechanisms including Type Ia supernovae and sporadic nuclear activity of a central supermassive black hole. As a starting point for this research we use an equilibrium hydrostatic recycling model (Kritsuk 1996). We show that a compact cooling inflow develops, if the heating is slightly insufficient to counterbalance radiative cooling of the hot gas in the central few kiloparsecs. An excessive heating in the centre, instead, drives a convectively unstable outflow. We model the onset of the instability and a quasi-steady convective regime in the core of the galaxy in two-dimensions assuming axial symmetry. Provided the power of net energy supply in the core is not too high, the convection remains subsonic. The convective pattern is dominated by buoyancy driven large-scale mushroom-like structures. Unlike in the case of a cooling inflow, the X-ray surface brightness of an (on average) isentropic convective core does not display a sharp maximum at the centre. A hybrid model, which combines a subsonic peripheral cooling inflow with an inner convective core, appears to be stable. We also discuss observational implications of these results.

We present the results of numerical simulations of the spherically symmetric gravitational collapse of supermassive stars (SMS). The collapse is studied using a general relativistic hydrodynamics code. The coupled system of Einstein and fluid equations is solved employing observer time coordinates, by foliating the spacetime by means of outgoing null hypersurfaces. The code contains an equation of state which includes effects due to radiation, electrons and baryons, and detailed microphysics to account for electron-positron pairs. In addition energy losses by thermal neutrino emission are included. We are able to follow the collapse of SMS from the onset of instability up to the point of black hole formation. Several SMS with masses in the range $5\times 10^5 M_{\odot}- 10^9 M_{\odot}$ are simulated. In all models an apparent horizon forms initially, enclosing the innermost 25% of the stellar mass. From the computed neutrino luminosities, estimates of the energy deposition by $\nu\bar{\nu}$-annihilation are obtained. Only a small fraction of this energy is deposited near the surface of the star, where, as proposed recently by Fuller & Shi (1998), it could cause the ultrarelativistic flow believed to be responsible for $\gamma$-ray bursts. Our simulations show that for collapsing SMS with masses larger than $5\times 10^5 M_{\odot}$ the energy deposition is at least two orders of magnitude too small to explain the energetics of observed long-duration bursts at cosmological redshifts. In addition, in the absence of rotational effects the energy is deposited in a region containing most of the stellar mass. Therefore relativistic ejection of matter is impossible.

We present results from simulations of axisymmetric relativistic rotational core collapse. The general relativistic hydrodynamic equations are formulated in flux-conservative form and solved using a high-resolution shock-capturing scheme. The Einstein equations are approximated with a conformally flat 3-metric. We use the quadrupole formula to extract waveforms of the gravitational radiation emitted during the collapse. A comparison of our results with those of Newtonian simulations shows that the wave amplitudes agree within 30%. Surprisingly, in some cases, relativistic effects actually diminish the amplitude of the gravitational wave signal. We further find that the parameter range of models suffering multiple coherent bounces due to centrifugal forces is considerably smaller than in Newtonian simulations.

We present tests and results of a new axisymmetric, fully general relativistic code capable of solving the coupled Einstein-matter system for a perfect fluid matter field. Our implementation is based on the Bondi metric, by which the spacetime is foliated with a family of outgoing light cones. We use high-resolution shock-capturing schemes to solve the fluid equations. The code can accurately maintain long-term stability of a spherically symmetric, relativistic, polytropic equilibrium model of a neutron star. In axisymmetry, we demonstrate global energy conservation of a perturbed neutron star in a compactified spacetime, for which the total energy radiated away by gravitational waves corresponds to a significant fraction of the Bondi mass.

After decades of one-dimensional nucleosynthesis calculations, the growth of computational resources has meanwhile reached a level, which for the first time allows astrophysicists to consider performing routinely realistic multidimensional nucleosynthesis calculations in explosive and, to some extent, also in non-explosive environments. In the present contribution we attempt to give a short overview of the physical and numerical problems which are encountered in these simulations. In addition, we assess the accuracy that can be currently achieved in the computation of nucleosynthetic yields, using multidimensional simulations of core collapse supernovae as an example.

Implementation details and test cases of a newly developed hydrodynamic code, AMRA, are presented. The numerical scheme exploits the adaptive mesh refinement technique coupled to modern high-resolution schemes which are suitable for relativistic and non-relativistic flows. Various physical processes are incorporated using the operator splitting approach, and include self-gravity, nuclear burning, physical viscosity, implicit and explicit schemes for conductive transport, simplified photoionization, and radiative losses from an optically thin plasma. Several aspects related to the accuracy and stability of the scheme are discussed in the context of hydrodynamic and astrophysical flows.

Motivated by observations of supernova SN 1987A, various authors have simulated Rayleigh-Taylor (RT) instabilities in the envelopes of core collapse supernovae (for a review, see Mueller 1998). The non-radial motion found in these simulations qualitatively agreed with observations in SN 1987A, but failed to explain the extent of mixing of newly synthesized 56Ni quantitatively. Here we present results of a 2D hydrodynamic simulation which re-addresses this failure and covers the entire evolution of the first 5 hours after core bounce.

New two-dimensional, high-resolution calculations of a core collapse supernova in a 15 Msol blue supergiant are presented, which cover the entire evolution from shock revival until the first few hours of the explosion. Explosive nucleosynthesis, its dependence upon convective mixing during the first second of the evolution and the growth of Rayleigh-Taylor instabilities at the composition interfaces of the progenitor star are all modeled consistently and allow for a comparison with observational data. We confirm our earlier findings, that the perturbations induced by neutrino driven convection are sufficiently strong to seed large-scale Rayleigh-Taylor mixing and to destroy the onion-shell structure of the stellar He-core. As in our earlier calculations, the strong deceleration of the nickel clumps in the layers adjacent to the He/H interface suggests that the high velocities of iron-group elements observed in SN 1987A cannot be explained on the basis of currently favored progenitor models. Possible solutions to this dilemma and the implications of the mixing for type Ib explosions are briefly discussed.

We have generalised the exact solution of the Riemann problem in special relativistic hydrodynamics for arbitrary tangential flow velocities. The solution is obtained by solving the jump conditions across shocks plus an ordinary differential equation arising from the self-similarity condition along rarefaction waves, in a similar way as in purely normal flow. The dependence of the solution on the tangential velocities is analysed, and the impact of this result on the development of multidimensional relativistic hydrodynamic codes (of Godunov type) is discussed.

Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ cone around the rotation axis. The jet flow is strongly beamed ($\la$ few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor ($R_{\ast} = 2.98 \times 10^{10} $cm) intact. At breakout the maximum Lorentz factor of the jet flow is 33. After breakout the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then being constant $\rho_{\rm ext} = 10^{-5}$ gcm$^{-3}$. Outside the star the flow begins to expand also laterally ($v \sim c$), but the beam remains very well collimated. At a distance of $2.54 R_{\ast}$, where the simulation ends, the Lorentz factor has increased to 44.

High-resolution two-dimensional simulations were performed for the first five minutes of the evolution of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout the inner half of the helium core by Rayleigh-Taylor instabilities at the Ni+Si/O and C+O/He interfaces, seeded by convective overturn during the early stages of the explosion. Fast moving nickel mushrooms with velocities up to about 4000 km/s are observed. This offers a natural explanation for the mixing required in light curve and spectral synthesis studies of Type Ib explosions. A continuation of the calculations to later times, however, indicates that the iron velocities observed in SN 1987 A cannot be reproduced because of a strong deceleration of the clumps in the dense shell left behind by the shock at the He/H interface.

We present results of high-resolution two-dimensional simulations which follow the first five minutes of a core collapse supernova explosion in a 15 solar mass blue supergiant progenitor. The computations start shortly after core bounce and include neutrino-matter interactions by using a light-bulb approximation for the neutrinos, and a treatment of the nucleosynthesis due to explosive silicon and oxygen burning. We find that newly formed iron-group elements are distributed throughout a significant fraction of the stellar helium core by the concerted action of convective and Rayleigh-Taylor instabilities. Fast moving nickel mushrooms with velocities up to 4000 km/s are observed. A continuation of the calculations to later times, however, indicates, that the iron velocities observed in SN 1987 A cannot be reproduced due to a strong deceleration of the clumps during their interaction with the dense shell left behind by the shock at the He/H interface. Therefore, we cannot confirm the claim that convective "premixing" of the nickel in the early phases of the explosion solves the problem of the high iron velocities.

We present the first radio emission simulations from high resolution three dimensional relativistic hydrodynamic jets, which allow for a study of the observational implications of the interaction between the jet and external medium. This interaction gives rise to a stratification of the jet where a fast spine is surrounded by a slow high energy shear layer. The stratification, and in particular the large specific internal energy and slow flow in the shear layer largely determines the emission from the jet. If the magnetic field in the shear layer becomes helical (e.g., resulting from an initial toroidal field and an aligned field component generated by shear) the emission shows a cross section asymmetry, in which either the top or the bottom of the jet dominates the emission. This, as well as limb or spine brightening, is a function of the viewing angle and flow velocity, and the top/bottom jet emission predominance can be reversed if the jet changes direction with respect to the observer, or presents a change in velocity. The asymmetry is more prominent in the polarized flux, because of field cancellation (or amplification) along the line of sight. Recent observations of jet cross section emission asymmetries in the blazar 1055+018 can be explained assuming the existence of a shear layer with a helical magnetic field.

We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r $\approx 3 10^{10} $cm) intact. At break-out the maximum Lorentz factor of the jet flow is about 33. Simulations have been performed with the GENESIS multi-dimensional relativistic hydrodynamic code.

The multidimensional relativistic hydrodynamical code GENESIS has been used to obtain first results of {\it 3D} simulations of relativistic jets. We have studied the influence of a slight perturbation of the injection velocity field on the morphodynamics of otherwise axisymmetric relativistic jets.

This review is concerned with a discussion of numerical methods for the solution of the equations of special relativistic hydrodynamics (SRHD). Particular emphasis is put on a comprehensive review of the application of high-resolution shock-capturing methods in SRHD. Results obtained with different numerical SRHD methods are compared, and two astrophysical applications of SRHD flows are discussed. An evaluation of the various numerical methods is given and future developments are analyzed.

Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to one in every computational zone throughout the simulation is fulfilled by renormalizing the mass fractions during the advection step. The proposed modification is appropriate for any difference scheme written in conservation form. Unlike other commonly used methods it does not violate the conservative character of the advection method. A new steepening mechanism, which is based on modification of interpolation profiles, is used to reduce numerical diffusion across composition discontinuities. Additional procedures are described, which are necessary to enforce monotonicity. Several numerical experiments are presented which demonstrate the capability of our Consistent Multi-fluid Advection (CMA) method in case of smooth and discontinuous distributions of fluid phases and under different hydrodynamic conditions. It is shown that due to the reduced diffusivity of the proposed scheme the abundance of some heavy elements obtained from hydrodynamic simulations of type II supernova explosions can change by a factor of a few in the most extreme cases.

In recent years the technique of Adaptive Mesh Refinement (AMR) has been successfully applied to several astrophysical problems. This method should also allow a consistent multi-dimensional modelling of supernova explosions. In this contribution we address some of the computational difficulties encountered when trying to apply the AMR method to explosive nucleosynthesis and supernova envelope ejection.

Using numerical techniques we studied the global stability of cooling flows in giant elliptical galaxies. As an initial equilibrium state we choose the hydrostatic gas recycling model (Kritsuk 1996). Non-equilibrium radiative cooling, stellar mass loss, heating by type Ia supernovae, distributed mass deposition, and thermal conductivity are included. Although the recycling model reproduces the basic X-ray observables, it appears to be unstable with respect to the development of inflow or outflow. In spherically symmetry the inflows are subject to a central cooling catastrophe, while the outflows saturate in a form of a subsonic galactic wind. Two-dimensional axisymmetric random velocity perturbations of the equilibrium model trigger the onset of a cooling catastrophe, which develops in an essentially non-spherical way. The simulations show a patchy pattern of mass deposition and the formation of hollow gas jets, which penetrate through the outflow down to the galaxy core. The X-ray observables of such a hybrid gas flow mimic those of the equilibrium recycling model, but the gas temperature exhibits a central depression. The mass deposition rate M_dot consists of two contributions of similar size: (i) a hydrostatic one resembling that of the equilibrium model, and (ii) a dynamical one which is related to the jets and is more concentrated to the centre. For a model galaxy, like NGC 4472, our 2D simulations predict M_dot = 2 M_sun/yr within the cooling radius for the advanced non-linear stage of the instability. We discuss the implications of these results to H_alpha nebulae and star formation in cooling flow galaxies and emphasize the need for high-resolution 3D simulations.

We report on the first three-dimensional hydrodynamic simulations of secular and dynamical non-axisymmetric instabilities in collapsing, rapidly rotating stellar cores which extend well beyond core bounce. The resulting gravitational radiation has been calculated using the quadrupole approximation. We find that secular instabilities do not occur during the simulated time interval of several 10 ms. Models which become dynamically unstable during core collapse show a strong nonlinear growth of non-axisymmetric instabilities. Both random and coherent large scale initial perturbations eventually give rise to a dominant bar-like deformation ($\exp(\pm im\phi$) with $m=2$). In spite of the pronounced tri-axial deformation of certain parts of the core no considerable enhancement of the gravitational radiation is found. This is due to the fact that rapidly rotating cores re-expand after core bounce on a dynamical time scale before non-axisymmetric instabilities enter the nonlinear regime. Hence, when the core becomes tri-axial, it is no longer very compact.

We present simulations of relativistic jets propagating in a nonuniform medium. Specifically, we study the bending of jets propagating obliquely to the vector of the density-gradient. Our results are applied to the NGC 4258, where such a medium is assumed to be provided by the atmosphere of the sub-parsec accretion disk tilted with respect to the original direction of the jet propagation. As a result, the jet is bent on a scale comparable to the density scaleheight of the disk atmosphere. The magnitude of the bending effect is found to be largest for light jets with low Lorentz factors. The predicted direction of bending is consistent with the observations.

One- and two-dimensional hydrodynamical simulations of neutrino-driven supernova explosions are discussed. The simulations cover the phase between the stagnation of the prompt shock and about one second after core bounce. Systematic variation of the neutrino fluxes from the neutrino sphere shows that the explosion energy, explosion time scale, initial mass of the protoneutron star, and explosive nucleosynthesis of iron-group elements depend sensitively on the strength of the neutrino heating during the first few 100 ms after shock formation. Convective overturn in the neutrino-heated region behind the shock is a crucial help for the explosion only in a narrow window of neutrino luminosities. Here powerful explosions can be obtained only in the multi-dimensional case. For higher core-neutrino fluxes also spherically symmetrical models yield energetic explosions, while for lower luminosities even with convection no strong explosions occur.

We present results concerning the internal structure and kinematics of disk galaxies formed in cosmologically motivated simulations. The calculations include dark matter, gas dynamics, radiative cooling, star formation, supernova feedback and metal enrichment. The initial model is a rigidly rotating overdense sphere with a mass of about 8 10^11 Msol which is perturbed by small scale fluctuations according to a biased CDM power spectrum. Converging, Jeans unstable and rapidly cooling regions are allowed to form stars. Via supernovae, metal enriched gas is returned to the interstellar medium. {}From these initial conditions a galaxy forms which shows the main properties of spiral galaxies: a rotationally supported exponential disk which consists of young stars with about solar metallicity, a slowly rotating halo of old metal poor stars, a bulge of old metal rich stars and a slowly rotating extended halo of dark matter. Bulge, stellar and dark halo are supported by an anisotropic velocity dispersion and have a de Vaucouleurs surface density profile. The flattening of the dark and stellar halo is too large to be explained by rotation only. Whether the flattening of the bulge is caused by an anisotropic velocity dispersion or by its rotation cannot be answered, because of the limited numerical resolution due to gravitational softening. The velocity dispersion and the thickness of the stellar disk increase with the age of the stars. Considering only the young stellar component, the disk is cold (sigma=20 km/sec) and thin (z <1 kpc). The dynamical formation process ends after about 4\,Gyr, when

An efficient algorithm for solving Poisson's equation in two and three spatial dimensions is discussed. The algorithm, which is described in detail, is based on the integral form of Poisson's equation and utilizes spherical coordinates and an expansion into spherical harmonics. The solver can be applied to and works well for all problems for which the use of spherical coordinates is appropriate. We also briefly discuss the implementation of the algorithm into hydrodynamic codes which are based on a conservative finite--difference scheme.

We present first results concerning the metallicities and stellar populations of galaxies formed in a cosmologically motivated simulation. The calculations include dark matter, gas dynamics, radiation processes, star formation, supernovae feedback, and metal enrichment. A rotating, overdense sphere with a mass of $8\,10^{11}\,$\Msol\ serves as initial model. Converging and Jeans unstable regions are allowed to form stars, which get their metallicity from the gas they are formed from. Via supernovae, metal enriched gas is given back to the interstellar medium. The forming galaxy shows the main properties of spiral galaxies: A metal rich bulge, a metal poor stellar halo and a disk of nearly solar composition. Halo and bulge consist predominantly of old stars ($>10.5\,$Gyrs). The disk has a metallicity gradient of $\dd(\log Z) / \dd r = - 0.05\,$kpc$^{-1}$, whereas the halo shows none. The models also exhibit a correlation between the metallicity of Pop~II stars and the power of small scale fluctuations. The stars of the bulge form from gas which is initially located in the largest maxima of the primordial density fluctuations, whereas the halo stars originate from gas accumulated in less pronounced maxima.

The formation of galaxies in a CDM cosmogony is investigated by following the evolution of dark and baryonic matter and of the frequency- dependent spatially averaged radiation field. The gas is allowed to form stars which are independently treated as a fourth collisionless component. By supernovae, energy and metal enriched gas is returned to the interstellar medium. Two extreme cases, the formation of an isolated field galaxy and the merging of two gas rich proto spirals are considered the latter being used as a model for galaxy formation in a cluster. The isolated galaxy is assumed to interact with the surrounding matter via tidal fields only, resulting in a rotation with a spin parameter of about 0.08. The forming field galaxy shows the main properties of spiral galaxies: A metal rich, rotatianally flattened bulge, a metal poor anisotrope stellar halo and a exponetial disk of nearly solar composition. The spiral is embedded in a triaxial halo of dark matter. The stars of the bulge form from gas which is initially located in the largest maxima of the primordial density fluctuations, whereas the halo stars originate from gas accumulated in less pronounced maxima. Merging proto spirals evolve to objects which can be interpreted as luminous ellipticals . The metal rich anisotrope elliptical rotates only slowly. Depending on the geometry of the merger event sometimes counter rotating gas is observed.