Source author record

Irina Sagert

Irina Sagert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Realistic Kilonova Up Close

Neutron star mergers are cosmic catastrophes that produce some of the most energetic observed phenomena: short gamma-ray bursts, gravitational wave signals, and kilonovae. The latter are optical transients, powered by radioactive nuclides which are synthesized when the neutron-rich ejecta of a disrupted neutron star undergoes decompression. We model this decompression phase using data from simulations of post-merger accretion disk winds. We use smoothed particle hydrodynamics with realistic nuclear heating to model the expansion over multiple scales, from initially several thousand km to billions of km. We then render a realistic image of a kilonova ejecta as it would appear for a nearby observer. This is the first time such a visualization is performed using input from state-of-the-art accretion disk simulations, nuclear physics and atomic physics. The volume rendering in our model computes an opacity transfer function on the basis of the physical opacity, varying significantly with the inhomogeneity of the neutron richness in the ejecta. Other physical quantities such as temperature or electron fraction can be visualized using an independent color transfer function. We discuss several difficulties with the ParaView application that we encountered during the visualization process, and give descriptions of our solutions and workarounds which could be used for future improvements.

preprint2015arXiv

Knudsen number dependence of 2D single-mode Rayleigh-Taylor fluid instabilities

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.

preprint2014arXiv

Hydrodynamic Shock Wave Studies within a Kinetic Monte Carlo Approach

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.

preprint2014arXiv

Kinetic Simulations of Rayleigh-Taylor Instabilities

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.

preprint2014arXiv

Symmetry energy impact in simulations of core-collapse supernovae

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

preprint2013arXiv

Building a Hydrodynamics Code with Kinetic Theory

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.

preprint2013arXiv

Parallelization of Kinetic Theory Simulations

Numerical studies of shock waves in large scale systems via kinetic simulations with millions of particles are too computationally demanding to be processed in serial. In this work we focus on optimizing the parallel performance of a kinetic Monte Carlo code for astrophysical simulations such as core-collapse supernovae. Our goal is to attain a flexible program that scales well with the architecture of modern supercomputers. This approach requires a hybrid model of programming that combines a message passing interface (MPI) with a multithreading model (OpenMP) in C++. We report on our approach to implement the hybrid design into the kinetic code and show first results which demonstrate a significant gain in performance when many processors are applied.

preprint2012arXiv

Nucleosynthesis in core-collapse supernova explosions triggered by a quark-hadron phase transition

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.

preprint2012arXiv

Soft nuclear equation-of-state from heavy-ion data and implications for compact stars

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.

preprint2011arXiv

Quark Matter In Massive Compact Stars

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.

preprint2010arXiv

Can a supernova bang twice?

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.

preprint2010arXiv

Strangeness in Astrophysics and Cosmology

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.

preprint2009arXiv

The Strange Prospects for Astrophysics

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.