FleCSPH is a smoothed particle hydrodynamics simulation tool, based on the compile-time configurable framework FleCSI. The asynchronous distributed tree topology combined with a fast multipole method allows FleCSPH to efficiently compute hydrodynamics and long range particle-particle interactions. FleCSPH provides initial data generators, particle relaxation techniques, and standard evolution drivers, which can be easily modified and extended to user-specific setups. Data input/output uses the H5part format, compatible with modern visualization software.

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle 3D Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of $^{16} \mathrm{O}$, $^{208} \mathrm{Pb}$ and $^{238} \mathrm{U}$ nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24\:fm with an average density of about $\rho = 0.05 \:\mathrm{fm}^{-3}$, proton fraction of $Y_p = 0.3$ and temperature of $T=0\:$MeV. Our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. In the present pasta calculations spin-orbit forces are not included but will be added in the future. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

We present a study of single-mode Rayleigh-Taylor instabilities (smRTI) with a modified Direct Simulation Monte Carlo (mDSMC) code in two dimensions. The mDSMC code is aimed to capture the dynamics of matter for a large range of Knudsen numbers within one approach. Our method combines the traditional Monte Carlo technique to efficiently propagate particles and the Point-of-Closest-Approach method for high spatial resolution. Simulations are performed using different particle mean-free-paths and we compare the results to linear theory predictions for the growth rate including diffusion and viscosity. We find good agreement between theoretical predictions and simulations and, at late times, observe the development of secondary instabilities, similar to hydrodynamic simulations and experiments. Large mean-free-paths favor particle diffusion, reduce the occurrence of secondary instabilities and approach the non-interacting gas limit.

We report on an ongoing project to develop a large scale Direct Simulation Monte Carlo code. The code is primarily aimed towards applications in astrophysics such as simulations of core-collapse supernovae. It has been tested on shock wave phenomena in the continuum limit and for matter out of equilibrium. In the current work we focus on the study of fluid instabilities. Like shock waves these are routinely used as test-cases for hydrodynamic codes and are discussed to play an important role in the explosion mechanism of core-collapse supernovae. As a first test we study the evolution of a single-mode Rayleigh-Taylor instability at the interface of a light and a heavy fluid in the presence of a gravitational acceleration. To suppress small-wavelength instabilities caused by the irregularity in the separation layer we use a large particle mean free path. The latter leads to the development of a diffusion layer as particles propagate from one fluid into the other. For small amplitudes, when the instability is in the linear regime, we compare its position and shape to the analytic prediction. Despite the broadening of the fluid interface we see a good agreement with the analytic solution. At later times we observe the development of a mushroom like shape caused by secondary Kelvin-Helmholtz instabilities as seen in hydrodynamic simulations and consistent with experimental observations.

Kinetic approaches are routinely employed to simulate the dynamics of systems that are too rarified to be described by the Navier-Stokes equations. However, generally they are far too computationally expensive to be applied for systems that are governed by continuum hydrodynamics. In this paper, we introduce a massively parallelized test-particle based kinetic Monte Carlo code that is capable of modeling the phase space evolution of an arbitrarily sized system that is free to move in and out of the continuum limit. Using particle mean free paths which are small with respect to the characteristic length scale of the simulated system, we retrieve continuum behavior, while non-equilibrium effects are observed when the mean free path is increased. To demonstrate the ability of our code to reproduce hydrodynamic solutions, we apply a test-suite of classic hydrodynamic shock problems. Simulations using tens of millions of test-particles are found to reproduce the analytical solutions well.

We present a review of a broad selection of nuclear matter equations of state (EOSs) applicable in core-collapse supernova studies. The large variety of nuclear matter properties, such as the symmetry energy, which are covered by these EOSs leads to distinct outcomes in supernova simulations. Many of the currently used EOS models can be ruled out by nuclear experiments, nuclear many-body calculations, and observations of neutron stars. In particular the two classical supernova EOS describe neutron matter poorly. Nevertheless, we explore their impact in supernova simulations since they are commonly used in astrophysics. They serve as extremely soft and stiff representative nuclear models. The corresponding supernova simulations represent two extreme cases, e.g., with respect to the protoneutron star (PNS) compactness and shock evolution. Moreover, in multi-dimensional supernova simulations EOS differences have a strong effect on the explosion dynamics. Because of the extreme behaviors of the classical supernova EOSs we also include DD2, a relativistic mean field EOS with density-dependent couplings, which is in satisfactory agreement with many current nuclear and observational constraints. This is the first time that DD2 is applied to supernova simulations and compared with the classical supernova EOS. We find that the overall behaviour of the latter EOS in supernova simulations lies in between the two extreme classical EOSs. As pointed out in previous studies, we confirm the impact of the symmetry energy on the electron fraction. Furthermore, we find that the symmetry energy becomes less important during the post bounce evolution, where conversely the symmetric part of the EOS becomes increasingly dominating, which is related to the high temperatures obtained. Moreover, we study the possible impact of quark matter at high densities and light nuclear clusters at low and intermediate dens

We report on the development of a test-particle based kinetic Monte Carlo code for large systems and its application to simulate matter in the continuum regime. Our code combines advantages of the Direct Simulation Monte Carlo and the Point-of-Closest-Approach methods to solve the collision integral of the Boltzmann equation. With that, we achieve a high spatial accuracy in simulations while maintaining computational feasibility when applying a large number of test-particles. The hybrid setup of our approach allows us to study systems which move in and out of the hydrodynamic regime, with low and high particle densities. To demonstrate our code's ability to reproduce hydrodynamic behavior we perform shock wave simulations and focus here on the Sedov blast wave test. The blast wave problem describes the evolution of a spherical expanding shock front and is an important verification problem for codes which are applied in astrophysical simulation, especially for approaches which aim to study core-collapse supernovae.

We study the implications on compact star properties of a soft nuclear equation of state determined from kaon production at subthreshold energies in heavy-ion collisions. On one hand, we apply these results to study radii and moments of inertia of light neutron stars. Heavy-ion data provides constraints on nuclear matter at densities relevant for those stars and, in particular, to the density dependence of the symmetry energy of nuclear matter. On the other hand, we derive a limit for the highest allowed neutron star mass of three solar masses. For that purpose, we use the information on the nucleon potential obtained from the analysis of the heavy-ion data combined with causality on the nuclear equation of state.

Measurements of kaon production at subthreshold energies in heavy-ion collisions point to a soft nuclear equation-of-state for densities up to 2-3 times nuclear matter saturation density. We apply these results to study the implications on compact star properties, especially in the context of the recent measurement of the two solar mass pulsar PSR J1614-2230. The implications are two-fold: Firstly, the heavy-ion results constrain nuclear matter at densities relevant to light neutron stars. Hence, a radius measurement could provide information about the density dependence of the symmetry energy which is a crucial quantity in nuclear physics. Secondly, the information on the nucleon potential obtained from the analysis of the heavy-ion data can be combined with restrictions from causality on the nuclear equation-of-state. From this we can derive a limit for the highest allowed compact star mass of three solar masses.

We explore heavy element nucleosynthesis in the explosion of massive stars which are triggered by a quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. The present study is based on general relativistic radiation hydrodynamics simulations with three-flavor Boltzmann neutrino transport in spherical symmetry, which utilize a quark-hadron hybrid equation of state based on the MIT bag model for strange quark matter. The quark-hadron phase transition inside the stellar core forms a shock wave propagating towards the surface of the proto-neutron star. The shock wave results in an explosion and ejects neutron-rich matter which is piled up or accreting on the proto-neutron star. Later, during the cooling phase, the proto-neutron star develops a proton-rich neutrino-driven wind. We present a detailed analysis of the nucleosynthesis outcome in both neutron-rich and proton-rich ejecta and compare our integrated nucleosynthesis with observations of metal poor stars.

We discuss the possible impact of strange quark matter on the evolution of core-collapse supernovae with emphasis on low critical densities for the quark-hadron phase transition. For such cases the hot proto-neutron star can collapse to a more compact hybrid star configuration hundreds of milliseconds after core-bounce. The collapse triggers the formation of a second shock wave. The latter leads to a successful supernova explosion and leaves an imprint on the neutrino signal. These dynamical features are discussed with respect to their compatibility with recent neutron star mass measurements which indicate a stiff high density nuclear matter equation of state.

We compare two classes of hybrid equations of state with a hadron-to-quark matter phase transition in their application to core collapse supernova simulations. The first one uses the quark bag model and describes the transition to three-flavor quark matter at low critical densities. The second one employs a Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) model with parameters describing a phase transition to two-flavor quark matter at higher critical densities. These models possess a distinctly different temperature dependence of their transition densities which turns out to be crucial for the possible appearance of quark matter in supernova cores. During the early post bounce accretion phase quark matter is found only if the phase transition takes place at sufficiently low densities as in the study based on the bag model. The increase critical density with increasing temperature, as obtained for our PNJL parametrization, prevents the formation of quark matter. The further evolution of the core collapse supernova as obtained applying the quark bag model leads to a structural reconfiguration of the central proto-neutron star where, in addition to a massive pure quark matter core, a strong hydrodynamic shock wave forms and a second neutrino burst is released during the shock propagation across the neutrinospheres. We discuss the severe constraints in the freedom of choice of quark matter models and their parametrization due to the recently observed 2 solar mass pulsar and their implications for further studies of core collapse supernovae in the QCD phase diagram.

The recent observation of the pulsar PSR J1614-2230 with a mass of 1.97 +/- 0.04 M_sun gives a strong constraint on the quark and nuclear matter equations of state (EoS). We explore the parameter ranges for a parameterized EoS for quark stars. We find that strange stars, made of absolutely stable strange quark matter, comply with the new constraint only if effects from the strong coupling constant and color-superconductivity are taken into account. Hybrid stars, compact stars with a quark matter core and an hadronic outer layer, can be as massive as 2 M_sun, but only for a significantly limited range of parameters. We demonstrate that the appearance of quark matter in massive stars depends crucially on the stiffness of the nuclear matter EoS. We show that the masses of hybrid stars stay below the ones of hadronic and pure quark stars, due to the softening of the EoS at the quark-hadron phase transition.

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition during the early post bounce phase of core-collapse supernovae. We construct a quark equation of state, based on the bag model for strange quark matter. The transition between the hadronic and the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equations of state are used in core-collapse supernova simulations, based on general relativistic radiation hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central protoneutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. A strong accretion shock forms, which propagates towards the protoneutron star surface. Due to the density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with matter outflow. This moment defines the onset of the explosion in supernova models that allow for a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock propagation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable identification for the structural reconfiguration of the stellar core. The ejected matter expands on a short timescale and remains neutron-rich. These conditions might be suitable for the production of heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven ejecta in the later cooling phase of the protoneutron star where the vp-process might occur.

The implications of a QCD phase transition at high temperatures and densities for core-collapse supernovae are discussed. For a strong first order phase transition to quark matter, various scenarios have been put forward in the literature. Here, detailed numerical simulations including neutrino transport are presented, where it is found that a second shock wave due to the QCD phase transition emerges shortly after bounce. It is demonstrated that such a supernova banging twice results in a second peak in the antineutrino spectrum. This second peak is clearly detectable in present neutrino detectors for a galactic supernova.

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

Predictions of the thermodynamic conditions for phase transitions at high baryon densities and large chemical potentials are currently uncertain and largely phenomenological. Neutrino observations of core-collapse supernovae can be used to constrain the situation. Recent simulations of stellar core collapse that include a description of quark matter predict a sharp burst of anti \nu_e several hundred milliseconds after the prompt \nu_e neutronization burst. We study the observational signatures of that anti \nu_e burst at current neutrino detectors - IceCube and Super-Kamiokande. For a Galactic core-collapse supernova, we find that signatures of the QCD phase transition can be detected, regardless of the neutrino oscillation scenario. The detection would constitute strong evidence of a phase transition in the stellar core, with implications for the equation of state at high matter density and the supernova explosion mechanism.

Explosive astrophysical systems, such as supernovae or compact star binary mergers, provide conditions where strange quark matter can appear. The high degree of isospin asymmetry and temperatures of several MeV in such systems may cause a transition to the quark phase already around saturation density. Observable signals from the appearance of quark matter can be predicted and studied in astrophysical simulations. As input in such simulations, an equation of state with an integrated quark matter phase transition for a large temperature, density and proton fraction range is required. Additionally, restrictions from heavy ion data and pulsar observation must be considered. In this work we present such an approach. We implement a quark matter phase transition in a hadronic equation of state widely used for astrophysical simulations and discuss its compatibility with heavy ion collisions and pulsar data. Furthermore, we review the recently studied implications of the QCD phase transition during the early post-bounce evolution of core-collapse supernovae and introduce the effects from strong interactions to increase the maximum mass of hybrid stars. In the MIT bag model, together with the strange quark mass and the bag constant, the strong coupling constant $\alpha_s$ provides a parameter to set the beginning and extension of the quark phase and with this the mass and radius of hybrid stars.

Some recent developments concerning the role of strange quark matter for astrophysical systems and the QCD phase transition in the early universe are addressed. Causality constraints of the soft nuclear equation of state as extracted from subthreshold kaon production in heavy-ion collisions are used to derive an upper mass limit for compact stars. The interplay between the viscosity of strange quark matter and the gravitational wave emission from rotation-powered pulsars are outlined. The flux of strange quark matter nuggets in cosmic rays is put in perspective with a detailed numerical investigation of the merger of two strange stars. Finally, we discuss a novel scenario for the QCD phase transition in the early universe, which allows for a small inflationary period due to a pronounced first order phase transition at large baryochemical potential.

We determine the Galactic production rate of strangelets as a canonical input to calculations of the measurable cosmic ray flux of strangelets by performing simulations of strange star mergers and combining the results with recent estimates of stellar binary populations. We find that the flux depends sensitively on the bag constant of the MIT bag model of QCD and disappears for high values of the bag constant and thus more compact strange stars. In the latter case strange stars could coexist with ordinary neutron stars as they are not converted by the capture of cosmic ray strangelets. An unambiguous detection of an ordinary neutron star would then not rule out the strange matter hypothesis.

The implications of the formation of strange quark matter in neutron stars and in core-collapse supernovae is discussed with special emphasis on the possibility of having a strong first order QCD phase transition at high baryon densities. If strange quark matter is formed in core-collapse supernovae shortly after the bounce, it causes the launch of a second outgoing shock which is energetic enough to lead to a explosion. A signal for the formation of strange quark matter can be read off from the neutrino spectrum, as a second peak in antineutrinos is released when the second shock runs over the neutrinosphere.

The possible role of a first order QCD phase transition at nonvanishing quark chemical potential and temperature for cold neutron stars and for supernovae is delineated. For cold neutron stars, we use the NJL model with nonvanishing color superconducting pairing gaps, which describes the phase transition to the 2SC and the CFL quark matter phases at high baryon densities. We demonstrate that these two phase transitions can both be present in the core of neutron stars and that they lead to the appearance of a third family of solution for compact stars. In particular, a core of CFL quark matter can be present in stable compact star configurations when slightly adjusting the vacuum pressure to the onset of the chiral phase transition from the hadronic model to the NJL model. We show that a strong first order phase transition can have strong impact on the dynamics of core collapse supernovae. If the QCD phase transition sets in shortly after the first bounce, a second outgoing shock wave can be generated which leads to an explosion. The presence of the QCD phase transition can be read off from the neutrino and antineutrino signal of the supernova.

We explore the implications of the QCD phase transition during the postbounce evolution of core-collapse supernovae. Using the MIT bag model for the description of quark matter and assuming small bag constants, we find that the phase transition occurs during the early postbounce accretion phase. This stage of the evolution can be simulated with general relativistic three-flavor Boltzmann neutrino transport. The phase transition produces a second shock wave that triggers a delayed supernova explosion. If such a phase transition happens in a future galactic supernova, its existence and properties should become observable as a second peak in the neutrino signal that is accompanied by significant changes in the energy of the emitted neutrinos. In contrast to the first neutronization burst, this second neutrino burst is dominated by the emission of anti-neutrinos because the electron-degeneracy is lifted when the second shock passes through the previously neutronized matter.

The effect of the QCD phase transition is studied for the mass-radius relation of compact stars and for hot and dense matter at a given proton fraction used as input in core-collapse supernova simulations. The phase transitions to the 2SC and CFL color superconducting phases lead to stable hybrid star configurations with a pure quark matter core. In supernova explosions quark matter could be easily produced due to $\beta$-equilibrium, small proton fractions and nonvanishing temperatures. A low critical density for the phase transition to quark matter is compatible with present pulsar mass measurements.

Recent astrophysical observations of neutron stars and heavy-ion data are confronted with our present understanding of the equation of state of dense hadronic matter. Emphasis is put on the possible role of the presence of hyperons in the interior of compact stars. We argue that data from low-mass pulsars provide an important cross-check between high-density astrophysics and heavy-ion physics.

We discuss the recent constraints on the nuclear equation of state from pulsar mass measurements and from subthreshold production of kaons in heavy-ion collisions. While recent pulsar data points towards a hard equation of state, the analysis of the heavy-ion data allows only for soft equations of state. We resolve the apparent contradiction by considering the different density regimes probed. We argue that future measurements of global properties of low-mass pulsars can serve as an excellent cross-check to heavy-ion data.

We discuss a pulsar acceleration mechanism based on asymmetric neutrino emission from the direct quark Urca process in the interior of proto neutron stars. The anisotropy is caused by a strong magnetic field which polarises the spin of the electrons opposite to the field direction. Due to parity violation the neutrinos and anti-neutrinos leave the star in one direction accelerating the pulsar. We calculate for varying quark chemical potentials the kick velocity in dependence of the quark phase temperature and its radius. Ignoring neutrino quark scattering we find that within a quark phase radius of 10 km and temperatures larger than 5 MeV kick velocities of 1000km s$^{-1}$ can be reached very easily. On the other hand taking into account the small neutrino mean free paths it seems impossible to reach velocities higher than 100km s$^{-1}$ even when including effects from colour superconductivity where the neutrino quark interactions are suppressed.

We discuss an acceleration mechanism for pulsars out of their supernova remnants based on asymmetric neutrino emission from quark matter in the presence of a strong magnetic field. The polarized electron spin fixes the neutrino emission from the direct quark Urca process in one direction along the magnetic field. We calculate the magnetic field strength which is required to polarize the electron spin as well as the required initial proto-neutron star temperature for a successfull acceleration mechanism. In addition we discuss the neutrino mean free paths in quark as well as in neutron matter which turn out to be very small. Consequently, the high neutrino interaction rates will wash out the asymmetry in neutrino emission. As a possible solution to this problem we take into account effects from colour superconductivity.

The puzzling phenomenon of pulsar kicks, i.e. the observed large escape velocities of pulsars out of supernova remnants, is examined for compact stars with a strange quark matter core. The direct Urca process in quark matter is studied in the presence of a strong magnetic field. Conditions for an asymmetric emission of the produced neutrinos are worked out in detail, giving constraints on the temperature, the strength of the magnetic field and the electron chemical potential in the quark matter core. In addition, the neutrino mean free paths for quark matter and a possible hadronic mantle are considered.

We discuss the possible appearance of strange exotic multi-quark states in the interior of neutron stars and signals for the existence of strange quark matter in the core of compact stars. We show how the in-medium properties of possible pentaquark states are constrained by pulsar mass measurements. The possibility of generating the observed large pulsar kick velocities by asymmetric emission of neutrinos from strange quark matter in magnetic fields is outlined.

We report on an undergraduate student project initiated in the summer semester of 2004 with the aim to establish equations of state for white dwarfs and neutron stars for computing mass-radius relations as well as corresponding maximum masses. First, white dwarfs are described by a Fermi gas model of degenerate electrons and neutrons and effects from general relativity are examined. For neutron star matter, the influence of a finite fraction of protons and electrons and of strong nucleon-nucleon interactions are studied. The nucleon--nucleon interactions are introduced within a Hartree-Fock scheme using a Skyrme-type interaction. Finally, masses and radii of neutron stars are computed for given central pressure.