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

Stephan Rosswog

Stephan 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.HEnucl-thastro-ph.SRgr-qcastro-phastro-ph.IMphysics.comp-phastro-ph.GA

Catalog footprint

What is connected

28works
8topics
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

28 published item(s)

preprint2026arXiv

Binary neutron star mergers with SPHINCS_BSSN: temperature-dependent equations of state and damping of constraint violations

Neutron star mergers hold the key to several grand challenges of contemporary (astro-)physics. In view of the upcoming next generation of ground-based detectors, it is crucial to keep improving theoretical predictions to harvest the full scientific returns from these investments. We introduce here a substantial update of our Lagrangian numerical relativity code SPHINCS_BSSN. Apart from changing our unit system, we add constraint damping terms to the BSSN spacetime evolution equations. We demonstrate that this measure reduces, without noteworthy computational cost, the Hamiltonian constraint violations by more than an order of magnitude. We further implement contributions to thermal energy and pressure that are based on Fermi liquid theory and contain a parametrization of the Dirac effective mass. These terms can be combined with any cold equation of state, and they enhance the physical realism of our simulations and introduce a physics-based concept of a temperature. In a set of merger simulations, we demonstrate good agreement with other temperature-dependent numerical relativity simulations. We find that different parametrizations of the Dirac effective mass can translate into shifts of $\sim 150$ Hz in the dominant post-merger gravitational wave peak frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlations of r-Process Elements in Very Metal-Poor Stars as Clues to their Nucleosynthesis Sites

Various nucleosynthesis studies have pointed out that the r-process elements in very metal-poor (VMP) halo stars might have different origins. By means of familiar concepts from statistics (correlations, cluster analysis, rank tests of elemental abundances), we look for causally correlated elemental abundance patterns and attempt to link them to astrophysical events. Some of these events produce the r-process elements jointly with iron, while others do not have any significant iron contribution. In the early stage of our Galaxy, at least three r-process nucleosynthesis sites have been active. The first two produce and eject iron and the majority of the lighter r-process elements. We assign them to two different types of core-collapse events, not identical to regular core-collapse supernovae (CCSNe), which produce only light trans-Fe elements. The third category is characterized by a strong r-process and responsible for the major fraction of the heavy main r-process elements without a significant co-production of Fe. It does not appear to be connected to CCSNe, in fact the Fe found in the related r-process enriched stars must come from previously occurring CCSNe. The existence of actinide boost stars indicates a further division among strong r-process sites. We assign these two strong r-process sites to neutron star mergers without fast black hole formation and to events where the ejecta are dominated by black hole accretion disk outflows. Indications from the lowest-metallicity stars hint at a connection to massive single stars (collapsars) forming black holes in the early Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A simple, entropy-based dissipation trigger for SPH

Smoothed Particle Hydrodynamics (SPH) schemes need to be enhanced by dissipation mechanisms to handle shocks. Most SPH formulations rely on artificial viscosity and while this is working well in pure shocks, attention has to be payed to avoid dissipation where it is not wanted. Commonly used approaches include limiters and time-dependent dissipation parameters. The former try to distinguish shocks from other types of flows that do not require dissipation while in the latter approach the dissipation parameters are steered by some source term ("trigger") and, if not triggered, they decay to a pre-described floor value. The commonly used source terms trigger on either compression, $-\nabla\cdot\vec{v}$, or its time derive. Here we explore a novel way to trigger SPH-dissipation: based on the entropy growth rate between two time steps we identify "troubled particles" that need to have dissipation added because they are either passing through a shock wave or become noisy. Our new scheme is implemented into the Lagrangian hydrodynamics code MAGMA2 and scrutinized in a number of shock and fluid instability tests. We find excellent results in shocks and only a moderate (and desired) switch-on in instability tests, despite our conservatively chosen trigger parameters. The new scheme is robust, trivial to implement into existing SPH codes and does not add any computational overhead.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Fate of the Merger Remnant in GW170817 and its Imprint on the Jet Structure

The first neutron star binary merger detected in gravitational waves, GW170817 and the subsequent detection of its emission across the electromagnetic spectrum showed that these systems are viable progenitors of short $γ$-ray bursts (sGRB). The afterglow signal of GW170817 has been found to be consistent with a structured GRB jet seen off-axis, requiring significant amounts of relativistic material at large angles. This trait can be attributed to the interaction of the relativistic jet with the external wind medium. Here we perform numerical simulations of relativistic jets interacting with realistic wind environments in order to explore how the properties of the wind and central engine affect the structure of successful jets. We find that the angular energy distribution of the jet depends primarily on the ratio between the lifetime of the jet and the time it takes the merger remnant to collapse. We make use of these simulations to constrain the time it took for the merger remnant in GW170817 to collapse into a black hole based on the angular structure of the jet as inferred from afterglow observations. We conclude that the lifetime of the merger remnant in GW170817 was $\approx 1-1.7$s, which, after collapse, triggered the formation of the jet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Lagrangian hydrodynamics code MAGMA2

We present the methodology and performance of the new Lagrangian hydrodynamics code MAGMA2, a Smoothed Particle Hydrodynamics code that benefits from a number of non-standard enhancements. By default it uses high-order smoothing kernels and wherever gradients are needed, they are calculated via accurate matrix inversion techniques, but a more conventional formulation with kernel gradients has also been implemented for comparison purposes. We also explore a matrix inversion formulation of SPH with a symmetrisation in the particle indices that is not frequently used. We find interesting advantages of this formulation in some of the tests, for example, a substantial reduction of surface tension effects for non-ideal particle setups and more accurate peak densities in Sedov blast waves. MAGMA2 uses artificial viscosity, but enhanced by techniques that are commonly used in finite volume schemes such as reconstruction and slope limiting. While simple to implement, this approach efficiently suppresses particle noise, but at the same time drastically reduces dissipation in locations where it is not needed and actually unwanted. We demonstrate the performance of the new code in a number of challenging benchmark tests including e.g. multi-dimensional vorticity creating Schulz-Rinne-type Riemann problems and more astrophysical tests such as a collision between two stars to demonstrate its robustness and excellent conservation properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tidal Disruptions of White Dwarfs: Theoretical Models and Observational Prospects

White dwarf stars that enter the tidal radius of black holes with masses $\lesssim$ $10^5$ M$_{\odot}$ are doomed to be ripped apart by tidal forces. Black holes in this mass range between stellar black holes and supermassive black holes have not been conclusively identified so the detection of a tidal disruption of a white dwarf would provide clear evidence for the existence of intermediate-mass black holes. In this review, we present a theoretical and observational overview of the transient events that result from the tidal disruptions of white dwarfs by intermediate-mass black holes. This includes discussion of the latest simulations and predicted properties, the results of observational searches, as well as a summary of the potential for gravitational wave emission to be detected with upcoming missions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutrino-driven winds in the aftermath of a neutron star merger: nucleosynthesis and electromagnetic transients

We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to $9 \cdot 10^{-3} M_\odot$ becomes unbound until $\sim 200~{\rm ms}$. Due to electron fractions of $Y_{\rm e} \approx 0.2 - 0.4$ mainly nuclei with mass numbers $A < 130$ are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about $4~{\rm h}$. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after $3-4~{\rm d}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic effects on tidal disruption kicks of solitary stars

Solitary stars that wander too close to their galactic centres can become tidally disrupted, if the tidal forces due to the supermassive black hole (SMBH) residing there overcome the self-gravity of the star. If the star is only partially disrupted, so that a fraction survives as a self-bound object, this remaining core will experience a net gain in specific orbital energy, which translates into a velocity "kick" of up to $\sim 10^3$ km/s. In this paper, we present the result of smoothed particle hydrodynamics (SPH) simulations of such partial disruptions, and analyse the velocity kick imparted on the surviving core. We compare $γ$ = 5/3 and $γ$ = 4/3 polytropes disrupted in both a Newtonian potential, and a generalized potential that reproduces most relativistic effects around a Schwarzschild black hole either exactly or to excellent precision. For the Newtonian case, we confirm the results of previous studies that the kick velocity of the surviving core is virtually independent of the ratio of the black hole to stellar mass, and is a function of the impact parameter $β$ alone, reaching at most the escape velocity of the original star. For a given $β$, relativistic effects become increasingly important for larger black hole masses. In particular, we find that the kick velocity increases with the black hole mass, making larger kicks more common than in the Newtonian case, as low-$β$ encounters are statistically more likely than high-$β$ encounters. The analysis of the tidal tensor for the generalized potential shows that our results are robust lower limits on the true relativistic kick velocities, and are generally in very good agreement with the exact results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Fate of the Compact Remnant in Neutron Star Mergers

Neutron star (binary neutron star and neutron star - black hole) mergers are believed to produce short-duration gamma-ray bursts. They are also believed to be the dominant source of gravitational waves to be detected by the advanced LIGO and the dominant source of the heavy r-process elements in the universe. Whether or not these mergers produce short-duration GRBs depends sensitively on the fate of the core of the remnant (whether, and how quickly, it forms a black hole). In this paper, we combine the results of merger calculations and equation of state studies to determine the fate of the cores of neutron star mergers. Using population studies, we can determine the distribution of these fates to compare to observations. We find that black hole cores form quickly only for equations of state that predict maximum non-rotating neutron star masses below 2.3-2.4 solar masses. If quick black hole formation is essential in producing gamma-ray bursts, LIGO observed rates compared to GRB rates could be used to constrain the equation of state for dense nuclear matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The interplay of disk wind and dynamical ejecta in the aftermath of neutron star - black hole mergers

We explore the evolution of the different ejecta components generated during the merger of a neutron star (NS) and a black hole (BH). Our focus is the interplay between material ejected dynamically during the merger, and the wind launched on a viscous timescale by the remnant accretion disk. These components are expected to contribute to an electromagnetic transient and to produce r-process elements, each with a different signature when considered separately. Here we introduce a two-step approach to investigate their combined evolution, using two- and three-dimensional hydrodynamic simulations. Starting from the output of a merger simulation, we identify each component in the initial condition based on its phase space distribution, and evolve the accretion disk in axisymmetry. The wind blown from this disk is injected into a three-dimensional computational domain where the dynamical ejecta is evolved. We find that the wind can suppress fallback accretion on timescales longer than ~100 ms. Due to self-similar viscous evolution, the disk accretion at late times nevertheless approaches a power-law time dependence $\propto t^{-2.2}$. This can power some late-time GRB engine activity, although the available energy is significantly less than in traditional fallback models. Inclusion of radioactive heating due to the r-process does not significantly affect the fallback accretion rate or the disk wind. We do not find any significant modification to the wind properties at large radius due to interaction with the dynamical ejecta. This is a consequence of the different expansion velocities of the two components.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Role of Fission in Neutron Star Mergers and its Impact on the r-Process Peaks

Comparing observational abundance features with nucleosynthesis predictions of stellar evolution or explosion simulations can scrutinize two aspects: (a) the conditions in the astrophysical production site and (b) the quality of the nuclear physics input utilized. We test the abundance features of r-process nucleosynthesis calculations for the dynamical ejecta of neutron star merger simulations based on three different nuclear mass models: The Finite Range Droplet Model (FRDM), the (quenched version of the) Extended Thomas Fermi Model with Strutinsky Integral (ETFSI-Q), and the Hartree-Fock-Bogoliubov (HFB) mass model. We make use of corresponding fission barrier heights and compare the impact of four different fission fragment distribution models on the final r-process abundance distribution. In particular, we explore the abundance distribution in the second r-process peak and the rare-earth sub-peak as a function of mass models and fission fragment distributions, as well as the origin of a shift in the third r-process peak position. The latter has been noticed in a number of merger nucleosynthesis predictions. We show that the shift occurs during the r-process freeze-out when neutron captures and β-decays compete and an (n,γ)-(γ,n) equilibrium is not maintained anymore. During this phase neutrons originate mainly from fission of material above A = 240. We also investigate the role of β-decay half-lives from recent theoretical advances, which lead either to a smaller amount of fissioning nuclei during freeze-out or a faster (and thus earlier) release of fission neutrons, which can (partially) prevent this shift and has an impact on the second and rare-earth peak as well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermonuclear detonations ensuing white dwarf mergers

The merger of two white dwarfs (WDs) has for many years not been considered as the favoured model for the progenitor system of type Ia supernovae (SNe Ia). But recent years have seen a change of opinion as a number of studies, both observational and theoretical, have concluded that they should contribute significantly to the observed type Ia supernova rate. In this paper, we study the ignition and propagation of detonation through post-merger remnants and we follow the resulting nucleosynthesis up to the point where a homologous expansion is reached. In our study we cover the entire range of WD masses and compositions. For the emergence of a detonation we study three different setups. The first two are guided by the merger remnants from our earlier simulations (Dan et al. 2014}, while for the third one the ignitions were set by placing hotspots with properties determined by spatially resolved calculations taken from the literature. There are some caveats to our approach which we investigate. We carefully compare the nucleosynthetic yields of successful explosions with SN Ia observations. Only three of our models are consistent with all the imposed constraints and potentially lead to a standard type Ia event. The first one, a $0.45\,M_\odot\ {\rm helium\ (He)} + 0.9 \,M_\odot$ carbon-oxygen (CO) WD system produces a sub-luminous, SN 1991bg-like event while the other two, a $0.45\,M_\odot\ {\rm He} + 1.1\,M_\odot$ oxygen-neon (ONe) WD system and a $1.05+1.05\,M_\odot$ system with two CO WDs, are good candidates for common SNe Ia.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Generalized Newtonian description of particle motion in spherically symmetric spacetimes

We present a generalized Newtonian description of particle dynamics valid for any spherically symmetric, static black hole spacetime. This approach is derived from the geodesic motion of test particles in the low-energy limit. It reproduces exactly the location of the marginally stable, marginally bound, and photon circular orbits; the radial dependence of the energy and angular momentum of circular orbits; parabolic motion; pericentre shift; and the spatial projection of general trajectories. As explicit examples of the new prescription, we apply it to the Schwarzschild, Schwarzschild-de Sitter, Reissner-Nordström, Ayón Beato-García, and Kehagias-Sfetsos spacetimes. In all of these examples, the orbital and epicyclic frequencies are reproduced to better than 10%. The resulting equations of motion can be implemented easily and efficiently within existing Newtonian frameworks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Implications of GRB 130603B and its macronova for r-process nucleosynthesis

The tentative identification of a Li-Paczynski macronova following the short GRB 130603B indicated that a few hundredths of a solar mass of neutron star matter were ejected and that this ejected mass has radioactively decayed into heavy r-process elements. If correct, this confirms long standing predictions (Eichler, Livio, Piran and Schramm 1989) that on the one hand, sGRBs are produced in compact binary mergers (CBMs) and on the other hand that these events are significant and possibly dominant sources of the heavy (A>130) r-process nuclei. Assuming that this interpretation is correct we obtain a lower limit of 0.02 m_sun on the ejected mass. Using the current estimates of the rate of sGRBs and with a beaming factor of 50, mergers associated with sGRBs can produce all the observed heavy r-process material in the Universe. We confront this conclusion with cosmochemistry and show that even though such events are rare, mixing is sufficient to account for the current homogeneous distribution of r-process material in the Galaxy. However, the appearance of significant amounts of Eu in some very low metallicity stars requires that some mergers took place very early on, namely with a very short time delay after the earliest star formation episodes. Alternatively, an additional early r-processsource may have contributed at that early stage. Finally, we note that evidence for short lived {244}^Pu in the very early solar system suggests that a merger of this kind took place within the vicinity of the solar system shortly (a few hundred million years) before its formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutrino-driven winds from neutron star merger remnants

We present a detailed, 3D hydrodynamics study of the neutrino-driven winds that emerge from the remnant of a NS merger. Our simulations are performed with the Newtonian, Eulerian code FISH, augmented by a detailed, spectral neutrino leakage scheme that accounts for heating due to neutrino absorption in optically thin conditions. Consistent with the 2D study of Dessart et al. (2009), we find that a strong baryonic wind is blown out along the original binary rotation axis within $100$ ms after the merger. We compute a lower limit on the expelled mass of $3.5 \times 10^{-3} M_{\odot}$, large enough to be relevant for heavy element nucleosynthesis. The physical properties vary significantly between different wind regions. For example, due to stronger neutrino irradiation, the polar regions show substantially larger $Y_e$ than those at lower latitudes. This has its bearings on the nucleosynthesis: the polar ejecta produce interesting r-process contributions from $A\sim 80$ to about 130, while the more neutron-rich, lower-latitude parts produce also elements up to the third r-process peak near $A\sim 195$. We also calculate the properties of electromagnetic transients that are powered by the radioactivity in the wind, in addition to the macronova transient that stems from the dynamic ejecta. The high-latitude (polar) regions produce UV/optical transients reaching luminosities up to $10^{41} {\rm erg \, s^{-1}}$, which peak around 1 day in optical and 0.3 days in bolometric luminosity. The lower-latitude regions, due to their contamination with high-opacity heavy elements, produce dimmer and more red signals, peaking after $\sim 2$ days in optical and infrared. Our numerical experiments indicate that it will be difficult to infer the collapse time-scale of the HMNS to a BH based on the wind electromagnetic transient, at least for collapse time-scales larger than the wind production time-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The long-term evolution of neutron star merger remnants - II. Radioactively powered transients

We use 3D hydrodynamic simulations of the long-term evolution of neutron star merger ejecta to predict the light curves of electromagnetic transients that are powered by the decay of freshly produced r-process nuclei. For the dynamic ejecta that are launched by tidal and hydrodynamic interaction, we adopt grey opacities of 10 cm$^2$/g, as suggested by recent studies. For our reference case of a 1.3-1.4 $M_\odot$ merger, we find a broad IR peak 2-4 d after the merger. The peak luminosity is $\approx 2\times 10^{40}$ erg/s for an average orientation, but increased by up to a factor of 4 for more favourable binary parameters and viewing angles. These signals are rather weak and hardly detectable within the large error box (~100 deg$^2$) of a gravitational wave trigger. A second electromagnetic transient results from neutrino-driven winds. These winds produce `weak' r-process material with $50 < A < 130$ and abundance patterns that vary substantially between different merger cases. For an adopted opacity of 1 cm$^2$/g, the resulting transients peak in the UV/optical about 6 h after the merger with a luminosity of $\approx 10^{41}$ erg/s (for a wind of 0.01 $M_\odot$) These signals are marginally detectable in deep follow-up searches (e.g. using Hypersuprime camera on Subaru). A subsequent detection of the weaker but longer lasting IR signal would allow an identification of the merger event. We briefly discuss the implications of our results to the recent detection of an nIR transient accompanying GRB 130603B.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An accurate Newtonian description of particle motion around a Schwarzschild black hole

A generalized Newtonian potential is derived from the geodesic motion of test particles in Schwarzschild spacetime. This potential reproduces several relativistic features with higher accuracy than commonly used pseudo-Newtonian approaches. The new potential reproduces the exact location of the marginally stable, marginally bound, and photon circular orbits, as well as the exact radial dependence of the binding energy and the angular momentum of these orbits. Moreover, it reproduces the orbital and epicyclic angular frequencies to better than 6%. In addition, the spatial projections of general trajectories coincide with their relativistic counterparts, while the time evolution of parabolic-like trajectories and the pericentre advance of elliptical-like trajectories are both reproduced exactly. We apply this approach to a standard thin accretion disc and find that the efficiency of energy extraction agrees to within 3% with the exact relativistic value, while the energy flux per unit area as a function of radius is reproduced everywhere to better than 7%. As a further astrophysical application we implement the new approach within a smoothed particle hydrodynamics code and study the tidal disruption of a main-sequence star by a supermassive black hole. The results obtained are in very good agreement with previous relativistic simulations of tidal disruptions in Schwarzschild spacetime. The equations of motion derived from this potential can be implemented easily within existing Newtonian hydrodynamics codes with hardly any additional computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How the merger of two white dwarfs depends on their mass ratio: orbital stability and detonations at contact

Despite their unique astrophysical relevance, the outcome of white dwarf binary mergers has so far only been studied for a very restricted number of systems. Here we present the results of a survey with more than two hundred simulations systematically scanning the white dwarf binary parameter space. We consider white dwarf masses ranging from 0.2 to 1.2 $M_\odot$ and account for their different chemical compositions. We find excellent agreement with the orbital evolution predicted by mass transfer stability analysis. Much of our effort in this paper is dedicated to determining which binary systems are prone to a thermonuclear explosion just prior to merger or at surface contact. We find that a large fraction of He-accreting binary systems explode: all dynamically unstable systems with accretor masses below 1.1 $M_\odot$ and donor masses above $\sim$ 0.4 $M_\odot$ are found to trigger a helium detonation at surface contact. A substantial fraction of these systems could explode at earlier times via detonations induced by instabilities in the accretion stream, as we have demonstrated in our previous work. We do not find definitive evidence for an explosion prior to merger or at surface contact in any of the studied double carbon-oxygen systems. Although we cannot exclude their occurrence if some helium is present, the available parameter space for a successful detonation in a white dwarf binary of pure carbon-oxygen composition is small. We demonstrate that a wide variety of dynamically unstable systems are viable type Ia candidates. The next decade thus holds enormous promise for the study of these events, in particular with the advent of wide-field synoptic surveys allowing a detailed characterization of their explosive properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Viscous Evolution of White Dwarf Merger Remnants

The merger of two white dwarfs (WDs) creates a differentially rotating remnant which is unstable to magnetohydrodynamic instabilities. These instabilities can lead to viscous evolution on a time-scale short compared to the thermal evolution of the remnant. We present multi-dimensional hydrodynamic simulations of the evolution of WD merger remnants under the action of an $α$-viscosity. We initialize our calculations using the output of eight WD merger simulations from Dan et al. (2011), which span a range of mass ratios and total masses. We generically find that the merger remnants evolve towards spherical states on time-scales of hours, even though a significant fraction of the mass is initially rotationally supported. The viscous evolution unbinds only a very small amount of mass $(< 10^{-5} M_\odot)$. Viscous heating causes some of the systems we study with He WD secondaries to reach conditions of nearly dynamical burning. It is thus possible that the post-merger viscous phase triggers detonation of the He envelope in some WD mergers, potentially producing a Type Ia supernova via a double detonation scenario. Our calculations provide the proper initial conditions for studying the long-term thermal evolution of WD merger remnants. This is important for understanding WD mergers as progenitors of Type Ia supernovae, neutron stars, R Coronae Borealis stars and other phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A fast recursive coordinate bisection tree for neighbour search and gravity

We introduce our new binary tree code for neighbour search and gravitational force calculations in an N-particle system. The tree is built in a "top-down" fashion by "recursive coordinate bisection" where on each tree level we split the longest side of a cell through its centre of mass. This procedure continues until the average number of particles in the lowest tree level has dropped below a prescribed value. To calculate the forces on the particles in each lowest-level cell we split the gravitational interaction into a near- and a far-field. Since our main intended applications are SPH simulations, we calculate the near-field by a direct, kernel-smoothed summation, while the far field is evaluated via a Cartesian Taylor expansion up to quadrupole order. Instead of applying the far-field approach for each particle separately, we use another Taylor expansion around the centre of mass of each lowest-level cell to determine the forces at the particle positions. Due to this "cell-cell interaction" the code performance is close to O(N) where N is the number of used particles. We describe in detail various technicalities that ensure a low memory footprint and an efficient cache use. In a set of benchmark tests we scrutinize our new tree and compare it to the "Press tree" that we have previously made ample use of. At a slightly higher force accuracy than the Press tree, our tree turns out to be substantially faster and increasingly more so for larger particle numbers. For four million particles our tree build is faster by a factor of 25 and the time for neighbour search and gravity is reduced by more than a factor of 6. In single processor tests with up to 10^8 particles we confirm experimentally that the scaling behaviour is close to O(N). The current Fortran 90 code version is OpenMP-parallel and scales excellently with the processor number (=24) of our test machine.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Prelude to a double degenerate merger: the onset of mass transfer and its impact on gravitational waves and surface detonations

We present the results of a systematic numerical study of the onset of mass transfer in double degenerate binary systems and its impact on the subsequent evolution. All investigated systems belong to the regime of direct impact, unstable mass transfer. In all of the investigated cases, even those considered unstable by conventional stability analysis, we find a long-lived mass transfer phase continuing for as many as several dozen orbital periods. This settles a recent debate sparked by a discrepancy between earlier SPH calculations that showed disruptions after a few orbital periods and newer grid-based studies in which mass transfer continued for tens of orbits. As we show that these binaries can survive at small separation for hundreds of orbital periods, their associated gravitational wave signal should be included when calculating the gravitational wave foreground (although expected to below LISA's sensitivity at these high frequencies). We also show that the inclusion of the entropy increase associated with shock-heating of the accreted material reduces the number of orbits a binary survives given the same initial conditions, although the effect is not as pronounced when using the appropriate initial conditions. The use of accurate initial conditions and a correct treatment of shock heating allows for a reliable time evolution of the temperature, density, and angular momentum, which are important when considering thermonuclear events that may occur during the mass transfer phase and/or after merger. Our treatment allows us to accurately identify when surface detonations may occur in the lead-up to the merger, as well as the properties of final merger products.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Precursors of short gamma-ray bursts

We carried out a systematic search of precursors on the sample of short GRBs observed by Swift. We found that ~8-10% of short GRBs display such early episode of emission. One burst (GRB 090510) shows two precursor events, the former ~13 s and the latter ~0.5 s before the GRB. We did not find any substantial difference between the precursor and the main GRB emission, and between short GRBs with and without precursors. We discuss possible mechanisms to reproduce the observed precursor emission within the scenario of compact object mergers. The implications of our results on quantum gravity constraints are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Going out with a bang: compact object collisions resulting from supernovae in binary systems

Binary star systems containing a neutron star or a black hole with an evolved, massive star are dynamically perturbed when the latter undergoes a supernova explosion. It is possible that the natal kick received by the newly-formed neutron star in the supernova may place the stellar remnants into a bound, highly eccentric orbit. In this case, the two compact objects can tidally interact and spiral into one another on a short timescale. The interaction with an accretion disc of supernova debris is also considered. We quantify the likelihood of such events and show that they would be expected to produce a high-energy transient, possibly a short gamma-ray burst, typically within a few days of the supernova.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Surface Detonations in Double Degenerate Binary Systems Triggered by Accretion Stream Instabilities

We present three-dimensional simulations on a new mechanism for the detonation of a sub-Chandrasekhar CO white dwarf in a dynamically unstable system where the secondary is either a pure He white dwarf or a He/CO hybrid. For dynamically unstable systems where the accretion stream directly impacts the surface of the primary, the final tens of orbits can have mass accretion rates that range from $10^{-5}$ to $10^{-3} M_{\odot}$ s$^{-1}$, leading to the rapid accumulation of helium on the surface of the primary. After $\sim 10^{-2}$ $M_{\odot}$ of helium has been accreted, the ram pressure of the hot helium torus can deflect the accretion stream such that the stream no longer directly impacts the surface. The velocity difference between the stream and the torus produces shearing which seeds large-scale Kelvin-Helmholtz instabilities along the interface between the two regions. These instabilities eventually grow into dense knots of material that periodically strike the surface of the primary, adiabatically compressing the underlying helium torus. If the temperature of the compressed material is raised above a critical temperature, the timescale for triple-$α$ reactions becomes comparable to the dynamical timescale, leading to the detonation of the primary's helium envelope. This detonation drives shockwaves into the primary which tend to concentrate at one or more focal points within the primary's CO core. If a relatively small amount of mass is raised above a critical temperature and density at these focal points, the CO core may itself be detonated.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The star ingesting luminosity of intermediate mass black holes in globular clusters

The dynamics of stars in the inner regions of nearby globular clusters (GCs) such as G1 indicate the presence of central concentrated dark masses, and one would like to know whether these are indeed intermediate mass black holes (IMBHs). As the number of surrounding stars, and their motions, are roughly known, the capture rate can be estimated; the question then arises of whether the apparent quiescence of the nuclei of these GCs is compatible with the IMBH's presence. The role of debris from disrupted stars in activating quiescent nuclei of GCs is studied here employing three-dimensional hydrodynamics simulations. It is argued that when individual stars are disrupted, the bulk of the debris would be swallowed or expelled rapidly compared with the interval between successive disruptions. Only a small fraction of the tightly bound mass is likely to be swallowed, yielding radiation with high efficiency, the surplus being expelled in a radiation-driven wind. A transient (predominantly of soft X-ray emission) signal could persist steadily with L ~ L_Edd =10^{41}(M_h/10^3M_sun) erg/s for at most tens of years; thereafter the flare would rapidly fade. While the infall rate declines as t^{-5/3}, some material may be stored for a longer time in an accretion disk. The IMBH luminosity could then remain as high as L_X > 10^{39} erg/s for several hundreds of years after disruption. In a given object, this ultraluminous X-ray activity would have a duty cycle of order 10^{-4}. Quiescent GCs, those with IMBHs now starved of fuel, should greatly outnumber active ones; observational constraints would not then be stringent until we had observed enough candidates to constitute a proper ensemble average.(abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Three-Dimensional Simulations of Tidally Disrupted Solar-Type Stars and the Observational Signatures of Shock Breakout

We describe a three-dimensional simulation of a $1 M_{\odot}$ solar-type star approaching a $10^{6} M_{\odot}$ black hole on a parabolic orbit with a pericenter distance well within the tidal radius. While falling towards the black hole, the star is not only stretched along the orbital direction but even more severely compressed at right angles to the orbit. The overbearing degree of compression achieved shortly after pericenter leads to the production of strong shocks which largely homogenize the temperature profile of the star, resulting in surface temperatures comparable to the initial temperature of the star's core. This phenomenon, which precedes the fallback accretion phase, gives rise to a unique double-peaked X-ray signature that, if detected, may be one of the few observable diagnostics of how stars behave under the influence of strong gravitational fields. If $\sim 10^{6} M_{\odot}$ black holes were prevalent in small or even dwarf galaxies, the nearest of such flares may be detectable by EXIST from no further away than the Virgo Cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

High Resolution Calculations of Merging Neutron Stars III: Gamma-Ray Bursts

Recent three dimensional, high-resolution simulations of neutron star coalescences are analysed to assess whether short gamma-ray bursts (GRBs) could originate from such encounters. The two most popular modes of energy extraction -- namely the annihilation of $ν\barν$ and magnetohydrodynamic processes- are explored in order to investigate their viability in launching GRBs. We find that $ν\barν$ annihilation can provide the necessary stresses to drive a highly relativistic expansion. But unless the outflow is beamed into less than one percent of the solid angle this mechanism may fail to explain the apparent isotropized energies implied for short GRBs at cosmological distances. We argue that the energetic, neutrino-driven wind that accompanies the merger event will have enough pressure to provide adequate collimation to the $ν\barν$-annihilation-driven jet, thereby comfortably satisfying constraints on event rate and apparent luminosity. We also assess magnetic mechanisms to transform the available energy into a GRB. If the central object does not collapse immediately into a black hole it will be convective and it is expected to act as an effective large scale dynamo amplifying the seed magnetic fields to a few times $10^{17}$ G within a small fraction of a second. The associated spindown time scale is 0.2 s, coinciding with the typical duration of a short GRB. The efficiencies of the various assessed magnetic processes are high enough to produce isotropized luminosities in excess of $10^{52}$ ergs/s even without beaming.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

An Adaptive Grid, Implicit Code for Spherically Symmetric, General Relativistic Hydrodynamics in Comoving Coordinates

We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and numerical diffusion are discussed. The time integration is performed with AGILE, an implicit solver for stiff algebrodifferential equations on a dynamical adaptive grid. We extend the adaptive grid technique, known from nonrelativistic hydrodynamics, to the general relativistic application and identify it with the concept of shift vectors in a 3+1 decomposition. The adaptive grid minimizes the number of required computational zones without compromising the resolution in physically important regions. Thus, the computational effort is greatly reduced when the zones are subject to computationally expensive additional processes, such as Boltzmann radiation transport or a nuclear reaction network. We present accurate results in the standard tests for supernova simulations: Sedov's point-blast explosion, the nonrelativistic and relativistic shock tube, the Oppenheimer-Snyder dust collapse, and homologous collapse.

0 citations0 reviewsNo code linkedOpen access