BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Rosswog

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

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-phastro-ph.IMastro-ph.SRastro-ph.COnucl-thphysics.comp-phastro-ph.GAphysics.flu-dyn

Catalog footprint

What is connected

32works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

32 published item(s)

preprint2025arXiv

A thermonuclear supernova interacting with hydrogen- and helium-deficient circumstellar material. SN 2020aeuh as a SN Ia-CSM-C/O?

Identifying the progenitors of thermonuclear supernovae (Type Ia supernovae; SNe Ia) remains a key objective in contemporary astronomy. The rare subclass of SNe Ia that interacts with circumstellar material (Type Ia-CSM) allows for studies of the progenitor's environment before explosion, and generally favours single-degenerate progenitor channels. The case of SN Ia-CSM PTF11kx clearly connected thermonuclear explosions with hydrogen-rich CSM-interacting events, and the more recent SN 2020eyj connected SNe Ia with helum-rich companion progenitors. Here we present a study of SN 2020aeuh, a Type Ia-CSM with delayed interaction. We analyse photometric and spectroscopic data that monitor the evolution of SN 2020aeuh and compare its properties with those of peculiar SNe Ia and core-collapse SNe. At early times, the evolution of SN 2020aeuh resembles a slightly overluminous SN Ia. Later, the interaction-dominated spectra develop the same pseudocontinuum seen in Type Ia-CSM PTF11kx and SN 2020eyj. However, the later-time spectra of SN 2020aeuh lack hydrogen and helium narrow lines. Instead, a few narrow lines could be attributed to carbon and oxygen. We fit the pseudobolometric light curve with a CSM-interaction mode, yielding a CSM mass of 1-2 M$_{\odot}$. We propose that SN 2020aeuh was a Type Ia supernova that eventually interacted with a dense medium which was deficient in both hydrogen and helium. Whereas previous SNe Ia-CSM constitute our best evidence for nondegenerate companion progenitors, the CSM around SN 2020aeuh is more difficult to understand. We include a hydrodynamical simulation for a double-degenerate system to showcase how the dynamical evolution of such a progenitor scenario could produce the CSM observed around SN 2020aeuh. It is clear that SN 2020aeuh challenges current models for stellar evolution leading up to a SN Ia explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Simulating neutron star mergers with the Lagrangian Numerical Relativity code SPHINCS_BSSN

We present the first neutron star merger simulations performed with the newly developed Numerical Relativity code SPHINCS_BSSN. This code evolves the spacetime on a mesh using the BSSN formulation, but matter is evolved via Lagrangian particles according to a high-accuracy version of general-relativistic Smooth Particle Hydrodynamics (SPH). Our code contains a number of new methodological elements compared to other Numerical Relativity codes. The main focus here is on the new elements that were introduced to model neutron star mergers. These include a) a refinement (fixed in time) of the spacetime-mesh, b) corresponding changes in the particle--mesh mapping algorithm and c) a novel way to construct SPH initial data for binary systems via the recently developed "Artificial Pressure Method." This latter method makes use of the spectral initial data produced by the library LORENE, and is implemented in a new code called SPHINCS_ID. While our main focus is on introducing these new methodological elements and documenting the current status of SPHINCS_BSSN, we also show as a first application a set of neutron star merger simulations employing "soft" ($Γ=2.00$) and "stiff" ($Γ=2.75$) polytropic equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thinking outside the box: Numerical Relativity with particles

To date, essentially all simulation codes that solve the full set of Einstein's equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations we focus on perfectly symmetric binary systems that are irrotational and have 1.3 $M_\odot$ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle-mesh mapping) and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the $Ψ_4$-formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times $10^{-3}$ \Msun as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity ($\sim0.5 .. 0.7c$) ejecta component of $\sim 10^{-4}$ \Msun that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission and it could also potentially cause a kilonova afterglow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Localization and broadband follow-up of the gravitational-wave transient GW150914

A gravitational-wave (GW) transient was identified in data recorded by the Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) detectors on 2015 September 14. The event, initially designated G184098 and later given the name GW150914, is described in detail elsewhere. By prior arrangement, preliminary estimates of the time, significance, and sky location of the event were shared with 63 teams of observers covering radio, optical, near-infrared, X-ray, and gamma-ray wavelengths with ground- and space-based facilities. In this Letter we describe the low-latency analysis of the GW data and present the sky localization of the first observed compact binary merger. We summarize the follow-up observations reported by 25 teams via private Gamma-ray Coordinates Network circulars, giving an overview of the participating facilities, the GW sky localization coverage, the timeline and depth of the observations. As this event turned out to be a binary black hole merger, there is little expectation of a detectable electromagnetic (EM) signature. Nevertheless, this first broadband campaign to search for a counterpart of an Advanced LIGO source represents a milestone and highlights the broad capabilities of the transient astronomy community and the observing strategies that have been developed to pursue neutron star binary merger events. Detailed investigations of the EM data and results of the EM follow-up campaign are being disseminated in papers by the individual teams.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Supplement: Localization and broadband follow-up of the gravitational-wave transient GW150914

This Supplement provides supporting material for arXiv:1602.08492 . We briefly summarize past electromagnetic (EM) follow-up efforts as well as the organization and policy of the current EM follow-up program. We compare the four probability sky maps produced for the gravitational-wave transient GW150914, and provide additional details of the EM follow-up observations that were performed in the different bands.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Boosting the accuracy of SPH techniques: Newtonian and special-relativistic tests

We study the impact of different discretization choices on the accuracy of SPH and we explore them in a large number of Newtonian and special-relativistic benchmark tests. As a first improvement, we explore a gradient prescription that requires the (analytical) inversion of a small matrix. For a regular particle distribution this improves gradient accuracies by approximately ten orders of magnitude and the SPH formulations with this gradient outperform the standard approach in all benchmark tests. Second, we demonstrate that a simple change of the kernel function can substantially increase the accuracy of an SPH scheme. While the "standard" cubic spline kernel generally performs poorly, the best overall performance is found for a high-order Wendland kernel which allows for only very little velocity noise and enforces a very regular particle distribution, even in highly dynamical tests. Third, we explore new SPH volume elements that enhance the treatment of fluid instabilities and, last, but not least, we design new dissipation triggers. They switch on near shocks and in regions where the flow --without dissipation-- starts to become noisy. The resulting new SPH formulation yields excellent results even in challenging tests where standard techniques fail completely.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetohydrodynamical simulations of a tidal disruption in general relativity

We perform hydro- and magnetohydrodynamical general relativistic simulations of a tidal disruption of a $0.1\,M_\odot$ red dwarf approaching a $10^5\,M_\odot$ non-rotating massive black hole on a close (impact parameter $β=10$) elliptical (eccentricity $e=0.97$) orbit. We track the debris self-interaction, circularization, and the accompanying accretion through the black hole horizon. We find that the relativistic precession leads to the formation of a self-crossing shock. The dissipated kinetic energy heats up the incoming debris and efficiently generates a quasi-spherical outflow. The self-interaction is modulated because of the feedback exerted by the flow on itself. The debris quickly forms a thick, almost marginally bound disc that remains turbulent for many orbital periods. Initially, the accretion through the black hole horizon results from the self-interaction, while in the later stages it is dominated by the debris originally ejected in the shocked region, as it gradually falls back towards the hole. The effective viscosity in the debris disc stems from the original hydrodynamical turbulence, which dominates over the magnetic component. The radiative efficiency is very low because of low energetics of the gas crossing the horizon and large optical depth that results in photon trapping.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-messenger picture of compact binary mergers

In the last decade, enormous progress has been achieved in the understanding of the various facets of coalescing double neutron star and neutron black hole binary systems. One hopes that the mergers of such compact binaries can be routinely detected with the advanced versions of the ground-based gravitational wave detector facilities, maybe as early as in 2016. From the theoretical side, there has also been mounting evidence that compact binary mergers could be major sources of heavy elements and these ideas have gained recent observational support from the detection of an event that has been interpreted as a "macronova", an electromagnetic transient powered by freshly produced, radioactively decaying heavy elements. In addition, compact binaries are the most plausible triggers of short gamma-ray bursts (sGRBs) and the last decade has witnessed the first detection of a sGRB afterglow and subsequent observations have delivered a wealth of information on the environments in which such bursts occur. To date, compact binary mergers can naturally explain most --though not all-- of the observed sGRB properties. This article reviews major recent developments in various areas related to compact binary mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

SPH Methods in the Modelling of Compact Objects

We review the current status of compact object simulations that are based on the Smooth Particle Hydrodynamics (SPH) method. The first main part of this review is dedicated to SPH as a numerical method. We begin by discussing relevant kernel approximation techniques and discuss the performance of different kernel functions. Subsequently, we review a number of different SPH formulations of Newtonian, special- and general relativistic ideal fluid dynamics. We particularly point out recent developments that increase the accuracy of SPH with respect to commonly used techniques. The second main part of the review is dedicated to the application of SPH in compact object simulations. We discuss encounters between two white dwarfs, between two neutron stars and between a neutron star and a stellar-mass black hole. For each type of system, the main focus is on the more common, gravitational wave-driven binary mergers, but we also discuss dynamical collisions as they occur in dense stellar systems such as cores of globular clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MODA: a new algorithm to compute optical depths in multi-dimensional hydrodynamic simulations

We introduce a new algorithm for the calculation of multidimensional optical depths in approximate radiative transport schemes, equally applicable to neutrinos and photons. Motivated by (but not limited to) neutrino transport in three-dimensional simulations of core-collapse supernovae and neutron star mergers, our method makes no assumptions about the geometry of the matter distribution, apart from expecting optically transparent boundaries. Based on local information about opacities, the algorithm figures out an escape route that tends to minimize the optical depth without assuming any pre-defined paths for radiation. Its adaptivity makes it suitable for a variety of astrophysical settings with complicated geometry (e.g., core-collapse supernovae, compact binary mergers, tidal disruptions, star formation, etc.). We implement the MODA algorithm into both a Eulerian hydrodynamics code with a fixed, uniform grid and into an SPH code where we make use a tree structure that is otherwise used for searching neighbours and calculating gravity. In a series of numerical experiments, we compare the MODA results with analytically known solutions. We also use snapshots from actual 3D simulations and compare the results of MODA with those obtained with other methods such as the global and local ray-by-ray method. It turns out that MODA achieves excellent accuracy at a moderate computational cost. In an appendix we also discuss implementation details and parallelization strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The long-term evolution of neutron star merger remnants - I. The impact of r-process nucleosynthesis

We follow the longterm evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the longterm evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce "strong" r-process elements with $A > 130$ with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield "weak" r-process contributions with $50 < A < 130$ whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any $^{56}{\rm Ni}$, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the "macronova" from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Europium production: neutron star mergers versus core-collapse supernovae

We have explored the Eu production in the Milky Way by means of a very detailed chemical evolution model. In particular, we have assumed that Eu is formed in merging neutron star (or neutron star black hole) binaries as well as in type II supernovae. We have tested the effects of several important parameters influencing the production of Eu during the merging of two neutron stars, such as: i) the time scale of coalescence, ii) the Eu yields and iii) the range of initial masses for the progenitors of the neutron stars. The yields of Eu from type II supernovae are very uncertain, more than those from coalescing neutron stars, so we have explored several possibilities. We have compared our model results with the observed rate of coalescence of neutron stars, the solar Eu abundance, the [Eu/Fe] versus [Fe/H] relation in the solar vicinity and the [Eu/H] gradient along the Galactic disc. Our main results can be summarized as follows: i) neutron star mergers can be entirely responsible for the production of Eu in the Galaxy if the coalescence time scale is no longer than 1 Myr for the bulk of binary systems, the Eu yield is around $3 \times 10^{-7}$ M$_\odot$, and the mass range of progenitors of neutron stars is 9-50 M$_\odot$; ii) both type II supernovae and merging neutron stars can produce the right amount of Eu if the neutron star mergers produce $2 \times 10^{-7}$ M$_\odot$ per system and type II supernovae, with progenitors in the range 20-50 M$_\odot$, produce yields of Eu of the order of $10^{-8}-10^{-9}$ M$_\odot$; iii) either models with only neutron stars producing Eu or mixed ones can reproduce the observed Eu abundance gradient along the Galactic disc.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Light from the Cosmic Frontier: Gamma-Ray Bursts

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Electromagnetic Signals of Compact Binary Mergers

Compact binary mergers are prime sources of gravitational waves, targeted by current and next generation detectors. The question "what is the observable electromagnetic (EM) signature of a compact binary merger?" is an intriguing one with crucial consequences to the quest for gravitational waves. We present a large set of numerical simulations that focus on the electromagnetic signals that emerge from the dynamically ejected sub-relativistic material. These outflows produce on a time scale of a day macronovae - short-lived IR to UV signals powered by radioactive decay. The interaction of this outflow with the surrounding matter inevitably leads to a long-lasting remnant. The expected radio signals of these remnants last longer than a year, when the sub-relativistic ejecta dominate the emission. We discuss their detectability in 1.4 GHz and 150 MHz and compare it with the detectability of short GRBs' orphan afterglows (which are produced by a different component of this outflow). Mergers with characteristics similar to those of the Galactic neutron star binary population (similar masses and typical circum-merger Galactic disk density of ~1cm^-3) taking place at the detection horizon of advanced GW detectors (300 Mpc) yield 1.4 GHz [150 MHz] signals of ~50 [300] microJy. The signal on time scales of weeks, is dominated by the mildly and/or ultra-relativistic outflow, which is not accounted for by our simulations, and is expected to be even brighter. Upcoming all sky surveys are expected to detect a few dozen, and possibly more, merger remnants at any given time thereby providing robust lower limits to the mergers rate even before the advanced GW detectorsbecome operational. Macronovae from the same distance peak in the IR to UV range at an observed magnitude that may be as bright as 22-23 about 10 hours after the merger but dimmer, redder and longer if the opacity is larger.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Gravitational Universe

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The multi-messenger picture of compact object encounters: binary mergers versus dynamical collisions

We explore the multi-messenger signatures of encounters between two neutron stars and between a neutron star and a stellar-mass black hole. We focus on the differences between gravitational wave driven binary mergers and dynamical collisions that occur, for example, in globular clusters. For both types of encounters we compare the gravitational wave and neutrino emission properties. We also calculate fallback rates and analyze the properties of the dynamically ejected matter. Last but not least we address the electromagnetic transients that accompany each type of encounter. The canonical nsns merger case ejects more than 1% of a solar mass of extremely neutron-rich ($Y_e\sim 0.03$) material, an amount that is consistent with double neutron star mergers being a major source of r-process in the galaxy. nsbh collisions eject very large amounts of matter ($\sim 0.15$ \msun) which seriously constrains their admissible occurrence rates. The compact object collision rate must therefore be less, likely much less, than 10% of the nsns merger rate. The radioactively decaying ejecta produce optical-UV "macronova" which, for the canonical merger case, peak after $\sim 0.4$ days with a luminosity of $\sim 10^{42}$ erg/s. nsns (nsbh) collisions reach up to 3 (7) times larger peak luminosities. The dynamic ejecta deposit a kinetic energy comparable to a supernova in the ambient medium. The canonical merger case releases approximately $2 \times 10^{50}$ erg, the most extreme (but likely rare) cases deposit kinetic energies of up to $10^{52}$ erg. The deceleration of this mildly relativistic material by the ambient medium produces long lasting radio flares. A canonical ns$^2$ merger at the detection horizon of advanced LIGO/Virgo produces a radio flare that peaks on a time scale of one year with a flux of $\sim$0.1 mJy at 1.4 GHz.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The runaway instability in general relativistic accretion disks

When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The structure and fate of white dwarf merger remnants

We present a large parameter study where we investigate the structure of white dwarf (WD) merger remnants after the dynamical phase. A wide range of WD masses and compositions are explored and we also probe the effect of different initial conditions. We investigated the degree of mixing between the WDs, the conditions for detonations as well as the amount of gas ejected. We find that systems with lower mass ratios have more total angular momentum and as a result more mass is flung out in a tidal tail. Nuclear burning can affect the amount of mass ejected. Many WD binaries that contain a helium-rich WD achieve the conditions to trigger a detonation. In contrast, for carbon-oxygen transferring systems only the most massive mergers with a total mass above ~2.1 solar masses detonate. Even systems with lower mass may detonate long after the merger if the remnant remains above the Chandrasekhar mass and carbon is ignited at the centre. Finally, our findings are discussed in the context of several possible observed astrophysical events and stellar systems, such as hot subdwarfs, R Coronae Borealis stars, single massive white dwarfs, supernovae of type Ia and other transient events. A large database containing 225 white dwarf merger remnants is made available via a dedicated web page.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the astrophysical robustness of neutron star merger r-process

In this study we explore the nucleosynthesis in the dynamic ejecta of compact binary mergers. We are particularly interested in the question how sensitive the resulting abundance patterns are to the parameters of the merging system. Therefore, we systematically investigate combinations of neutron star masses in the range from 1.0 to 2.0 \Msun and, for completeness, we compare the results with those from two simulations of a neutron star black hole merger. The ejecta masses vary by a factor of five for the studied systems, but all amounts are (within the uncertainties of the merger rates) compatible with being a major source of cosmic r-process. The ejecta undergo a robust r-process nucleosynthesis which produces all the elements from the second to the third peak in close-to-solar ratios. Most strikingly, this r-process is extremely robust, all 23 investigated binary systems yield practically identical abundance patterns. This is mainly the result of the ejecta being extremely neutron rich (\ye $\approx0.04$) and the r-process path meandering along the neutron drip line so that the abundances are determined entirely by nuclear rather than by astrophysical properties. This robustness together with the ease with which both the second and third peak are reproduced make compact binary mergers the prime candidate for the source of the observed unique heavy r-process component.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the point mass approximation to calculate the gravitational wave signal from white dwarf binaries

Double white dwarf binaries in the Galaxy dominate the gravitational wave sky and would be detectable for an instrument such as LISA. Most studies have calculated the expected gravitational wave signal under the assumption that the binary white dwarf system can be represented by two point masses in orbit. We discuss the accuracy of this approximation for real astrophysical systems. For non-relativistic binaries in circular orbit the gravitational wave signal can easily be calculated. We show that for these systems the point mass approximation is completely justified when the individual stars are axisymmetric irrespective of their size. We find that the signal obtained from Smoothed-Particle Hydrodynamics simulations of tidally deformed, Roche-lobe filling white dwarfs, including one case when an accretion disc is present, is consistent with the point mass approximation. The difference is typically at the level of one per cent or less in realistic cases, yielding small errors in the inferred parameters of the binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The dynamic ejecta of compact object mergers and eccentric collisions

Compact object mergers eject neutron-rich matter in a number of ways: by the dynamical ejection mediated by gravitational torques, as neutrino-driven winds and probably also a good fraction of the resulting accretion disc finally becomes unbound by a combination of viscous and nuclear processes. If compact binary mergers produce indeed gamma-ray bursts there should also be an interaction region where an ultra-relativistic outflow interacts with the neutrino-driven wind and produces moderately relativistic ejecta. Each type of ejecta has different physical properties and therefore plays a different role for nucleosynthesis and for the electromagnetic transients that go along with compact object encounters. Here we focus on the dynamic ejecta and present results for over 30 hydrodynamical simulations of both gravitational wave-driven mergers and parabolic encounters as they may occur in globular clusters. We find that mergers eject $\sim 1$% of a solar mass of extremely neutron-rich material. The exact amount as well as the ejection velocity depends on the involved masses with asymmetric systems ejecting more material at higher velocities. This material undergoes a robust r-process and both ejecta amount and abundance pattern are consistent with neutron star mergers being a major source of the "heavy" ($A>130$) r-process isotopes. Parabolic collisions, especially those between neutron stars and black holes, eject substantially larger amounts of mass and therefore cannot occur frequently without overproducing galactic r-process matter. We also discuss the electromagnetic transients that are powered by radioactive decays within the ejecta ("Macronovae"), and the radio flares that emerge when the ejecta dissipate their large kinetic energies in the ambient medium.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Compact binary mergers: an astrophysical perspective

This paper reviews the current understanding of double neutron star and neutron star black hole binaries. It addresses mainly (nuclear) astrophysics aspects of compact binary mergers and thus complements recent reviews that have emphasized the numerical relativity viewpoint. In particular, the paper discusses different channels to release neutron-rich matter into the host galaxy, connections between compact binary mergers and short Gamma-ray bursts and accompanying electromagnetic signals.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Conservative, special-relativistic smoothed particle hydrodynamics

We present and test a new, special-relativistic formulation of Smoothed Particle Hydrodynamics (SPH). Our approach benefits from several improvements with respect to earlier relativistic SPH formulations. It is self-consistently derived from the Lagrangian of an ideal fluid and accounts for special-relativistic "grad-h terms". In our approach, we evolve the canonical momentum and the canonical energy per baryon and thus circumvent some of the problems that have plagued earlier formulations of relativistic SPH. We further use a much improved artificial viscosity prescription which uses the extreme local eigenvalues of the Euler equations and triggers selectively on a) shocks and b) velocity noise. The shock trigger accurately monitors the relative density slope and uses it to fine-tune the amount of artificial viscosity that is applied. This procedure substantially sharpens shock fronts while still avoiding post-shock noise. If not triggered, the viscosity parameter of each particle decays to zero. None of these viscosity triggers is specific to special relativity, both could also be applied in Newtonian SPH. The performance of the new scheme is explored in a large variety of benchmark tests where it delivers excellent results. Generally, the grad-h terms deliver minor, though worthwhile, improvements. The scheme performs close to perfect in supersonic advection tests, but also in strong relativistic shocks, usually considered a particular challenge for SPH, the method yields convincing results. For example, due to its perfect conservation properties, it is able to handle Lorentz-factors as large as $γ= 50 \; 000$ in the so-called wall shock test. Moreover, we find convincing results in a rarely shown, but challenging test that involves so-called relativistic simple waves and also in multi-dimensional shock tube tests.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Special-relativistic Smoothed Particle Hydrodynamics: a benchmark suite

In this paper we test a special-relativistic formulation of Smoothed Particle Hydrodynamics (SPH) that has been derived from the Lagrangian of an ideal fluid. Apart from its symmetry in the particle indices, the new formulation differs from earlier approaches in its artificial viscosity and in the use of special-relativistic ``grad-h-terms''. In this paper we benchmark the scheme in a number of demanding test problems. Maybe not too surprising for such a Lagrangian scheme, it performs close to perfectly in pure advection tests. What is more, the method produces accurate results even in highly relativistic shock problems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Collisions of white dwarfs as a new progenitor channel for type Ia supernovae

We present the results of a systematic numerical study of an alternative progenitor scenario to produce type Ia supernova explosions, which is not restricted to the ignition of a CO white dwarf near the Chandrasekhar mass. In this scenario, a shock-triggered thermonuclear explosion ensues from the collision of two white dwarfs. Consistent modeling of the gas dynamics together with nuclear reactions using both a smoothed particle and a grid-based hydrodynamics code are performed to study the viability of this alternative progenitor channel. We find that shock-triggered ignition and the synthesis of Ni are in fact a natural outcome for moderately massive white dwarf pairs colliding close to head-on. We use a multi-dimensional radiative transfer code to calculate the emergent broadband light curves and spectral time-series of these events. The synthetic spectra and lightcurves compare well with those of normal type Ia supernovae over a similar B-band decline rate and are broadly consistent with the Phillips relation, although a mild dependence on viewing angle is observed due to the asymmetry of the ejected debris. While event rates within galactic centers and globular clusters are found to be much too low to explain the bulk of the type Ia supernovae population, they may be frequent enough to make as much as a one percent contribution to the overall rate. Although these rate estimates are still subject to substantial uncertainties, they do suggest that dense stellar systems should provide upcoming supernova surveys with hundreds of such collision-induced thermonuclear explosions per year.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Tidal disruption and ignition of white dwarfs by moderately massive black holes

We present a numerical investigation of the tidal disruption of white dwarfs by moderately massive black holes, with particular reference to the centers of dwarf galaxies and globular clusters. Special attention is given to the fate of white dwarfs of all masses that approach the black hole close enough to be disrupted and severely compressed to such extent that explosive nuclear burning can be triggered. Consistent modeling of the gas dynamics together with the nuclear reactions allows for a realistic determination of the explosive energy release. In the most favorable cases, the nuclear energy release may be comparable to that of typical type Ia supernovae. Although the explosion will increase the mass fraction escaping on hyperbolic orbits, a good fraction of the debris remains to be swallowed by the hole, causing a bright soft X-ray flare lasting for about a year. Such transient signatures, if detected, would be a compelling testimony for the presence of a moderately mass black hole (below $10^5 M_\odot$).

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

From Neutron Star Binaries to Gamma-ray bursts

I summarize recent results about how a neutron star binary coalescence can produce short gamma-ray bursts (GRBs). Two possibilities are discussed: the annihilation of neutrino anti-neutrino pairs above the merged remnant and the exponential amplification of magnetic fields in the central object up to values close to equipartition. We find that the neutrino annihilation drives bipolar, relativistic outflows with Lorentz-factors large enough to circumvent the GRB 'compactness problem'. The total energy within these outflows is moderate by GRB-standards ($\sim 10^{48}-10^{49}$ ergs), but the interaction with the baryonic material blown-off by the neutrinos collimates the outflows into opening angles of typically 0.1 sterad, yielding isotropic energies close to $10^{51}$ ergs. We further want to stress the plausibility of the central object resisting the immediate collapse to a black hole. In this case the central object will --similar to a proto-neutron star-- be subject to neutrino driven convection that --together with the rapid, differential rotation-- will lead to a drastic amplification of pre-existing magnetic fields. Within fractions of a second, field strengths comparable to equipartition field strength ($> 10^{17}$ G) will be reached. These will produce large torques that will spin-down the object within about 0.2 s, and would thus naturally explain the duration of short GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

High Resolution Calculations of Merging Neutron Stars II: Neutrino Emission

The remnant resulting from the merger of two neutron stars produces neutrinos in copious amounts. In this paper we present the neutrino emission results obtained via Newtonian, high-resolution simulations of the coalescence event. These simulations use three-dimensional smoothed particle hydrodynamics together with a nuclear, temperature dependent equation of state and a multi-flavour neutrino leakage scheme. We present the details of our scheme, discuss the neutrino emission results from a neutron star coalescence and compare them to the core-collapse supernova case where neutrino emission has been studied for several decades. The average neutrino energies are similar to those in the supernova case, but contrary to the latter, the luminosities are dominated by electron-type antineutrinos which are produced in the hot, neutron-rich, thick disk of the merger remnant. The cooler parts of this disk contain substantial fractions of heavy nuclei, which, however, do not influence the overall neutrino emission results significantly. Our total neutrino luminosities from the merger event are considerably lower than those found in previous investigations. This has serious consequences for the ability to produce a gamma-ray burst via neutrino annihilation. The neutrinos are emitted preferentially along the initial binary rotation axis, an event seen ``pole-on'' would appear much brighter in neutrinos than a similar event seen ``edge-on''.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

High Resolution Calculations of Merging Neutron Stars I: Model Description and Hydrodynamic Evolution

We present the results of 3D, high-resolution calculations of neutron star mergers. Using Newtonian gravity, but adding the forces from the emission of gravitational waves, we solve the hydrodynamics equations using the smoothed particle hydrodynamics method. The microscopic properties of the neutron star matter are described in terms of a new nuclear equation of state (EOS) based on the tables of Shen et al. (1998). The effects of neutrino emission of all flavours are accounted for. The new EOS is noticeably stiffer than the Lattimer-Swesty-EOS used in previous investigations. We find in general somewhat lower temperatures than previous investigations. This allows a substantial fraction of heavy nuclei to be present even in the inner parts of the torus that forms around the central object and yields important contributions to the neutrino opacity. The central object shows a rapid differential rotation with periods of $\sim$ 2ms, which we suspect to stabilize the central object on (at least) the simulation time scale against gravitational collapse. The differential rotation could wind up initial magnetic fields to $\sim 10^17$ G.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Post Newtonian SPH

We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Coalescing Neutron Stars: A Solution to the R-Process Problem?

Recent calculations place questions on the ability of supernovae to produce r-process nuclei in the correct amounts. We present results from 3D Newtonian SPH calculations of mergers of equal mass neutron star binaries. We find the amounts of ejected mass to be very promising for an explanation of the observed r-process abundances. Preliminary abundance calculations show that practically all the material is subject to r-process conditions. The calculated abundance patterns reproduce the basic features of the solar r-process abundances very well. Especially the peak around A= 195 is easily reproduced without any tuning of the initial entropies. Thus, eventually all the observed r-process material could be explained by mass ejection during neutron star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Mass ejection in neutron star mergers

We present the results of 3D Newtonian SPH simulations of the merger of a neutron star binary. The microscopic properties of matter are described by the physical equation of state of Lattimer and Swesty (LS-EOS). To test for the robustness of our results we check the sensitivity to the approximations of our model as well as to the binary system parameters. The main and new result is that for the realistic LS-EOS, depending on the initial spin, between 4e-3 and 4e-2 solar masses of material become unbound. If, as suggested, large parts of this matter consist of r-process nuclei, neutron star mergers could account for the whole observed r-process material in the Galaxy.

0 citations0 reviewsNo code linkedOpen access