Massive stars end their life in an explosion event with kinetic energies of the order 1 Bethe. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large entropies per baryon of about 100 and the fast expansion timescale, the conditions found in our model are unlikely to favor...

By performing axisymmetric hydrodynamic simulations of core-collapse supernovae with spectral neutrino transport based on the isotropic diffusion source approximation scheme, we support the assumption that the neutrino-heating mechanism aided by the standing accretion shock instability and convection can initiate an explosion of a 13 $M_{\odot}$ star. Our results show that bipolar explosions are more likely to be associated with models which include rotation. We point out that models, which form a north-south symmetric bipolar explosion, can lead to larger explosion energies than for the corresponding unipolar explosions.

We have performed a set of 11 three-dimensional magnetohydrodynamical core collapse supernova simulations in order to investigate the dependencies of the gravitational wave signal on the progenitor's initial conditions. We study the effects of the initial central angular velocity and different variants of neutrino transport. Our models are started up from a 15 solar mass progenitor and incorporate an effective general relativistic gravitational potential and a finite temperature nuclear equation of state. Furthermore, the electron flavour neutrino transport is tracked by efficient algorithms for the radiative transfer of massless fermions. We find that non- and slowly rotating models show gravitational wave emission due to prompt- and lepton driven convection that reveals details about the hydrodynamical state of the fluid inside the protoneutron stars. Furthermore we show that protoneutron stars can become dynamically unstable to rotational instabilities at T/|W| values as low as ~2 % at core bounce. We point out that the inclusion of deleptonization during the postbounce phase is very important for the quantitative GW prediction, as it enhances the absolute values of the gravitational wave trains up to a factor of ten with respect to a lepton-conserving treatment.

We present the gravitational wave (GW) analysis of an extensive series of 3D MHD core-collapse simulations. Our 25 models are launched from a 15 solar mass progenitor, a spherically symmetric effective general relativistic potential, the Lattimer-Swesty or the Shen equation of state (EoS), and a neutrino parametrisation scheme which is accurate until about 5ms postbounce. For 3 representative models, we also include long-term neutrino physics by means of a leakage scheme. Non- or only slowly rotating models show GW emission due to prompt and proto-neutron star convection, allowing the distinction between the two different nuclear EoS. For moderately or fast rotation rates models, we find, in agreement with recent results, only a type I GW signature at core bounce. Models which are set up with an initial central angular velocity of >~ 2pi rad/s emit GWs due to the low T/|W| dynamical instability during the postbounce phase. Weak B-fields do not notably influence the dynamical evolution of the core and thus the GW emission. However, for strong initial poloidal B-fields (~1e12 G),flux-freezing and field winding leads to conditions where P_{mag}/P_{mat} ~ 1, causing the onset of a jet-like supernova explosion and hence the emission of a type IV GW signal. In contradiction to axisymmetric simulations, we find evidence that nonaxisymmetric fluid modes can counteract or even suppress jet formation for models with strong initial toroidal B-fields. We point out that the inclusion of the deleptonisation during the postbounce phase is an indispensable issue for the quantitative prediction of GWs from core-collapse supernovae, as it can alter the GW amplitude up to a factor of 10 compared to a pure hydrodynamical treatment.

FISH is a fast and simple ideal magneto-hydrodynamics code that scales to ~10 000 processes for a Cartesian computational domain of ~1000^3 cells. The simplicity of FISH has been achieved by the rigorous application of the operator splitting technique, while second order accuracy is maintained by the symmetric ordering of the operators. Between directional sweeps, the three-dimensional data is rotated in memory so that the sweep is always performed in a cache-efficient way along the direction of contiguous memory. Hence, the code only requires a one-dimensional description of the conservation equations to be solved. This approach also enable an elegant novel parallelisation of the code that is based on persistent communications with MPI for cubic domain decomposition on machines with distributed memory. This scheme is then combined with an additional OpenMP parallelisation of different sweeps that can take advantage of clusters of shared memory. We document the detailed implementation of a second order TVD advection scheme based on flux reconstruction. The magnetic fields are evolved by a constrained transport scheme. We show that the subtraction of a simple estimate of the hydrostatic gradient from the total gradients can significantly reduce the dissipation of the advection scheme in simulations of gravitationally bound hydrostatic objects. Through its simplicity and efficiency, FISH is as well-suited for hydrodynamics classes as for large-scale astrophysical simulations on high-performance computer clusters. In preparation for the release of a public version, we demonstrate the performance of FISH in a suite of astrophysically orientated test cases.

We discuss the formation of stellar mass black holes via protoneutron star (PNS) collapse. In the absence of an earlier explosion, the PNS collapses to a black hole due to the continued mass accretion onto the PNS. We present an analysis of the emitted neutrino spectra of all three flavors during the PNS contraction. Special attention is given to the physical conditions which depend on the input physics, e.g. the equation of state (EoS) and the progenitor model. The PNSs are modeled as the central object in core collapse simulations using general relativistic three-flavor Boltzmann neutrino transport in spherical symmetry. The simulations are launched from several massive progenitors of 40 and 50 solar mass. We analyze the electron-neutrino luminosity dependencies and construct a simple approximation for the electron-neutrino luminosity, which depends only on the physical conditions at the electron-neutrinosphere. In addition, we analyze different mu/tau-neutrino pair-reactions separately and compare the differences during the post-bounce phase of failed core collapse supernova explosions of massive progenitors. We also investigate the connection between the increasing mu/tau-neutrino luminosity and the PNS contraction during the accretion phase before black hole formation. Comparing the different post bounce phase of the progenitor models under investigation, we find large differences in the emitted neutrino spectra. These differences and the analysis of the electron-neutrino luminosity indicate a strong progenitor model dependency of the emitted neutrino signal.

Astrophysical observations originate from matter that interacts with radiation or transported particles. We develop a pragmatic approximation in order to enable multi-dimensional simulations with basic spectral radiative transfer when the computational resources are not sufficient to solve the complete Boltzmann transport equation. The distribution function of the transported particles is decomposed into trapped and streaming particle components. Their separate evolution equations are coupled by a source term that converts trapped particles into streaming particles. We determine this source term by requiring the correct diffusion limit. For a smooth transition to the free streaming regime, this 'diffusion source' is limited by the matter emissivity. The resulting streaming particle emission rates are integrated over space to obtain the streaming particle flux. A geometric estimate of the flux factor is used to convert the particle flux to the streaming particle density. The efficiency of the scheme results from the freedom to use different approximations for each particle component. In supernovae, reactions with trapped particles on fast time scales establish equilibria that reduce the number of primitive variables required to evolve the trapped particle component. On the other hand, a stationary-state approximation facilitates the treatment of the streaming particle component. Different approximations may apply in applications to stellar atmospheres, star formation, or cosmological radiative transfer. We compare the isotropic diffusion source approximation with Boltzmann neutrino transport of electron flavour neutrinos in spherically symmetric supernova models and find good agreement. An extension of the scheme to the multi-dimensional case is also discussed.

We present the gravitational wave analysis from rotating (model s15g) and nearly non-rotating (model s15h) 3D MHD core collapse supernova simulations at bounce and the first couple of ten milliseconds afterwards. The simulations are launched from 15M_{\odot} progenitor models stemming from stellar evolution calculations. Gravity is implemented by a spherically symmetric effective general relativistic potential. The input physics uses the Lattimer-Swesty equation of state for hot, dense matter and a neutrino parametrisation scheme that is accurate until the first few ms after bounce. The 3D simulations allow us to study features already known from 2D simulations as well as nonaxisymmetric effects. In agreement with recent results we find only type I gravitational wave signals at core bounce. In the later stage of the simulations, one of our models (s15g) shows nonaxisymmetric gravitational wave emission caused by a low T/|W| dynamical instability, while the other model radiates gravitational waves due to a convective instability in the protoneutron star. The total energy released in gravitational waves within the considered time intervals is 1.52\times10^{-7}M_{\odot} (s15g) and 4.72\times10^{-10}M_{\odot} (s15h). Both core collapse simulations indicate that corresponding events in our Galaxy would be detectable either by the LIGO or Advanced LIGO detector.

Few years ago, Boltzmann neutrino transport led to a new and reliable generation of spherically symmetric models of stellar core collapse and postbounce evolution. After the failure to prove the principles of the supernova explosion mechanism, these sophisticated models continue to illuminate the close interaction between high-density matter under extreme conditions and the transport of leptons and energy in general relativistically curved space-time. We emphasize that very different input physics is likely to be relevant for the different evolutionary phases, e.g. nuclear structure for weak rates in collapse, the equation of state of bulk nuclear matter during bounce, multidimensional plasma dynamics in the postbounce evolution, and neutrino cross sections in the explosive nucleosynthesis. We illustrate the complexity of the dynamics using preliminary 3D MHD high-resolution simulations based on parameterized deleptonization. With established spherically symmetric models we show that typical features of the different phases are reflected in the predicted neutrino signal and that a consistent neutrino flux leads to electron fractions larger than 0.5 in neutrino-driven supernova ejecta.

We present the results of a series of calculations studying the collapse of molecular cloud cores performed using a three-dimensional smoothed particle hydr odynamics code with radiative transfer in the flux-limited diffusion approximation. The opacities and specific heat capacities are identical for each calculation. However, we find that the temperature evolution during the simulations varies significantly when starting from different initial conditions. Even spherically-symmetric clouds with different initial densities show markedly different development. We conclude that simple barotropic equations of state like those used in some previous calculations provide at best a crude approximation to the thermal behaviour of the gas. Radiative transfer is necessary to obtain accurate temperatures.

We describe a new, faster implicit algorithm for solving the radiation hydrodynamics equations in the flux-limited diffusion approximation for smoothed particle hydrodynamics. This improves on the method elucidated in Whitehouse & Bate by using a Gauss-Seidel iterative method rather than iterating over the exchange of energy between pairs of particles. The new algorithm is typically many thousands of times faster than the old one, which will enable more complex problems to be solved. The new algorithm is tested using the same tests performed by Turner & Stone for ZEUS-2D, and repeated by Whitehouse & Bate.

We describe the implementation and testing of a smoothed particle hydrodynamics (SPH) code that solves the equations of radiation hydrodynamics in the flux-limited diffusion (FLD) approximation. The SPH equations of radiation hydrodynamics for an explicit integration scheme are derived and tested. We also discuss the implementation of an implicit numerical scheme for solving the radiation equations that allows the system to be evolved using timesteps much longer than the explicit radiation timestep. The code is tested on a variety of one-dimensional radiation hydrodynamics problems including radiation propagating in an optically thin medium, optically thick and thin shocks, subcritical and supercritical radiating shocks, and a radiation dominated shock. Many of the tests were also performed by Turner and Stone (2001) to test their implementation of a FLD module for the ZEUS-2D code. The SPH code performs at least as well as the ZEUS-2D code in these tests.

The geodesic equations resulting from the Schwarzschild gravitational metric element are solved exactly including the contribution from the Cosmological constant. The exact solution is given by genus 2 Siegelsche modular forms. For zero cosmological constant the hyperelliptic curve degenerates into an elliptic curve and the resulting geodesic is solved by the Weierstra$\ss$ Jacobi modular form. The solution is applied to the precise calculation of the perihelion precession of the orbit of planet Mercury around the Sun.

Strong field (exact) solutions of the gravitational field equations of General Relativity in the presence of a Cosmological Constant are investigated. In particular, a full exact solution is derived within the inhomogeneous Szekeres-Szafron family of space-time line element with a nonzero Cosmological Constant. The resulting solution connects, in an intrinsic way, General Relativity with the theory of modular forms and elliptic curves. The homogeneous FLRW limit of the above space-time elements is recovered and we solve exactly the resulting Friedmann Robertson field equation with the appropriate matter density for generic values of the Cosmological Constant %Lambda and curvature constant K. A formal expression for the Hubble constant is derived. The cosmological implications of the resulting non-linear solutions are systematically investigated. Two particularly interesting solutions i) the case of a flat universe K=0, Lambda not= 0 and ii) a case with all three cosmological parameters non-zero, are described by elliptic curves with the property of complex multiplication and absolute modular invariant j=0 and 1728, respectively. The possibility that all non-linear solutions of General Relativity are expressed in terms of theta functions associated with Riemann-surfaces is discussed.

This paper describes how the non-gravitational contribution to Galactic Velocity Rotation Curves can be explained in terms of a negative Cosmological Constant ($\Lambda$). It will be shown that the Cosmological Constant leads to a velocity contribution proportional to the radii, at large radii, and depending on the mass of the galaxy. This explanation contrasts with the usual interpretation that this effect is due to Dark Matter halos. The velocity rotation curve for the galaxy NGC 3198 will be analysed in detail, while several other galaxies will be studied superficially. The Cosmological Constant derived experimentally from the NGC 3198 data was found to be:$|\Lambda|_{Exp}= 5.0\times 10^{-56} cm^{-2}$. This compares favourably with the theoretical value obtained from the Large Number Hypothesis of: $|\Lambda|_{Theory}=2.1\times 10^{-56}cm^{-2}$. The Extended LNH is then used to define other cosmological parameters: gravitational modification constant, energy density, and the Cosmological Constant in terms of a fundamental length. A speculative theory for the evolution of the Universe is outlined where it is shown how the Universe can be defined, in any particular era, by two parameters: the fundamental length and the energy density of the vacuum for that epoch. The theory is applied to the time evolution of the universe where a possible explanation for the $\rho_{Planck}/\rho_{\Lambda}^{QH} \approx 10^{120}$ problem is proposed. The nature of the ''vacuum'' is reviewed along with a speculative approach for calculating the Cosmological Constant via formal M-theory.The experimentally derived results presented in this paper support a decelerating Universe, in contrast with recent indicationsfrom Type Ia Supernovae experiments, for an accelerating Universe.

In this paper a new theory of Dark Matter is proposed. Experimental analysis of several Galaxies show how the non-gravitational contribution to galactic Velocity Rotation Curves can be interpreted as that due to the Cosmological Constant $\Lambda$. The experimentally determined values for $\Lambda$ are found to be consistent with those expected from Cosmological Constraints. The Cosmological Constant is interpreted as leading to a constant energy density which in turn can be used to partly address the energy deficit problem (Dark Energy) of the Universe. The work presented here leads to the conclusion that the Cosmological Constant is negative and that the universe is de-accelerating. This is in clear contradiction to the Type Ia Supernovae results which support an accelerating universe.