Normalized to: Rosswog, S.
[1]
oai:arXiv.org:2002.01950 [pdf] - 2118898
Observational constraints on the optical and near-infrared emission from
the neutron star-black hole binary merger S190814bv
Ackley, K.;
Amati, L.;
Barbieri, C.;
Bauer, F. E.;
Benetti, S.;
Bernardini, M. G.;
Bhirombhakdi, K.;
Botticella, M. T.;
Branchesi, M.;
Brocato, E.;
Bruun, S. H.;
Bulla, M.;
Campana, S.;
Cappellaro, E.;
Castro-Tirado, A. J.;
Chambers, K. C.;
Chaty, S.;
Chen, T. -W.;
Ciolfi, R.;
Coleiro, A.;
Copperwheat, C. M.;
Covino, S.;
Cutter, R.;
D'Ammando, F.;
D'Avanzo, P.;
De Cesare, G.;
D'Elia, V.;
Della Valle, M.;
Denneau, L.;
De Pasquale, M.;
Dhillon, V. S.;
Dyer, M. J.;
Elias-Rosa, N.;
Evans, P. A.;
Eyles-Ferris, R. A. J.;
Fiore, A.;
Fraser, M.;
Fruchter, A. S.;
Fynbo, J. P. U.;
Galbany, L.;
Gall, C.;
Galloway, D. K.;
Getman, F. I.;
Ghirlanda, G.;
Gillanders, J. H.;
Gomboc, A.;
Gompertz, B. P.;
González-Fernández, C.;
González-Gaitán, S.;
Grado, A.;
Greco, G.;
Gromadzki, M.;
Groot, P. J.;
Gutiérrez, C. P.;
Heikkilä, T.;
Heintz, K. E.;
Hjorth, J.;
Hu, Y. -D.;
Huber, M. E.;
Inserra, C.;
Izzo, L.;
Japelj, J.;
Jerkstrand, A.;
Jin, Z. P.;
Jonker, P. G.;
Kankare, E.;
Kann, D. A.;
Kennedy, M.;
Kim, S.;
Klose, S.;
Kool, E. C.;
Kotak, R.;
Kuncarayakti, H.;
Lamb, G. P.;
Leloudas, G.;
Levan, A. J.;
Longo, F.;
Lowe, T. B.;
Lyman, J. D.;
Magnier, E.;
Maguire, K.;
Maiorano, E.;
Mandel, I.;
Mapelli, M.;
Mattila, S.;
McBrien, O. R.;
Melandri, A.;
Michałowski, M. J.;
Milvang-Jensen, B.;
Moran, S.;
Nicastro, L.;
Nicholl, M.;
Guelbenzu, A. Nicuesa;
Nuttal, L.;
Oates, S. R.;
O'Brien, P. T.;
Onori, F.;
Palazzi, E.;
Patricelli, B.;
Perego, A.;
Torres, M. A. P.;
Perley, D. A.;
Pian, E.;
Pignata, G.;
Piranomonte, S.;
Poshyachinda, S.;
Possenti, A.;
Pumo, M. L.;
Quirola-Vásquez, J.;
Ragosta, F.;
Ramsay, G.;
Rau, A.;
Rest, A.;
Reynolds, T. M.;
Rosetti, S. S.;
Rossi, A.;
Rosswog, S.;
Sabha, N. B.;
Carracedo, A. Sagués;
Salafia, O. S.;
Salmon, L.;
Salvaterra, R.;
Savaglio, S.;
Sbordone, L.;
Schady, P.;
Schipani, P.;
Schultz, A. S. B.;
Schweyer, T.;
Smartt, S. J.;
Smith, K. W.;
Smith, M.;
Sollerman, J.;
Srivastav, S.;
Stanway, E. R.;
Starling, R. L. C.;
Steeghs, D.;
Stratta, G.;
Stubbs, C. W.;
Tanvir, N. R.;
Testa, V.;
Thrane, E.;
Tonry, J. L.;
Turatto, M.;
Ulaczyk, K.;
van der Horst, A. J.;
Vergani, S. D.;
Walton, N. A.;
Watson, D.;
Wiersema, K.;
Wiik, K.;
Wyrzykowski, L.;
Yang, S.;
Yi, S. -X.;
Young, D. R.
Submitted: 2020-02-05, last modified: 2020-06-22
On 2019 August 14, the LIGO and Virgo interferometers detected a
high-significance event labelled S190814bv. Preliminary analysis of the GW data
suggests that the event was likely due to the merger of a compact binary system
formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave
sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out
an intensive multi-epoch, multi-instrument observational campaign to identify
the possible optical/near infrared counterpart of the event. In addition, the
ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a
search on this event. Our observations allow us to place limits on the presence
of any counterpart and discuss the implications for the kilonova (KN) possibly
generated by this NS-BH merger, and for the strategy of future searches.
Altogether, our observations allow us to exclude a KN with large ejecta mass
$M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can
exclude much smaller masses in a subsample of our observations. This disfavours
the tidal disruption of the neutron star during the merger. Despite the
sensitive instruments involved in the campaign, given the distance of S190814bv
we could not reach sufficiently deep limits to constrain a KN comparable in
luminosity to AT 2017gfo on a large fraction of the localisation probability.
This suggests that future (likely common) events at a few hundreds Mpc will be
detected only by large facilities with both high sensitivity and large field of
view. Galaxy-targeted observations can reach the needed depth over a relevant
portion of the localisation probability with a smaller investment of resources,
but the number of galaxies to be targeted in order to get a fairly complete
coverage is large, even in the case of a localisation as good as that of this
event.
[2]
oai:arXiv.org:1912.01095 [pdf] - 2112574
A simple, entropy-based dissipation trigger for SPH
Submitted: 2019-12-02, last modified: 2020-06-10
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.
[3]
oai:arXiv.org:2004.00102 [pdf] - 2074062
Axisymmetric Radiative Transfer Models of Kilonovae
Korobkin, Oleg;
Wollaeger, Ryan;
Fryer, Christopher;
Hungerford, Aimee L.;
Rosswog, Stephan;
Fontes, Christopher;
Mumpower, Matthew;
Chase, Eve;
Even, Wesley;
Miller, Jonah M.;
Misch, G. Wendell;
Lippuner, Jonas
Submitted: 2020-03-31
The detailed observations of GW170817 proved for the first time directly that
neutron star mergers are a major production site of heavy elements. The
observations could be fit by a number of simulations that qualitatively agree,
but can quantitatively differ (e.g. in total r-process mass) by an order of
magnitude. We categorize kilonova ejecta into several typical morphologies
motivated by numerical simulations, and apply a radiative transfer Monte Carlo
code to study how the geometric distribution of the ejecta shapes the emitted
radiation. We find major impacts on both spectra and light curves. The peak
bolometric luminosity can vary by two orders of magnitude and the timing of its
peak by a factor of five. These findings provide the crucial implication that
the ejecta masses inferred from observations are uncertain by at least an order
of magnitude. Mixed two-component models with lanthanide-rich ejecta are
particularly sensitive to geometric distribution. A subset of mixed models
shows very strong viewing angle dependence due to lanthanide "curtaining",
which persists even if the relative mass of lanthanide-rich component is small.
The angular dependence is weak in the rest of our models, but different
geometric combinations of the two components lead to a highly diverse set of
light curves. We identify geometry-dependent P Cygni features in late spectra
that directly map out strong lines in the simulated opacity of neodymium, which
can help to constrain the ejecta geometry and to directly probe the r-process
abundances.
[4]
oai:arXiv.org:2004.00146 [pdf] - 2074066
Tidal Disruptions of White Dwarfs: Theoretical Models and Observational
Prospects
Submitted: 2020-03-31
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.
[5]
oai:arXiv.org:1911.13093 [pdf] - 2006767
The Lagrangian hydrodynamics code MAGMA2
Submitted: 2019-11-29
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. It uses high-order smoothing
kernels and wherever gradients are needed, they are calculated via accurate
matrix inversion techniques. Our default version does not make use of any
kernel gradients, but a more conventional formulation has also been implemented
for comparison purposes. 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. complex, multi-dimensional Riemann problems and
more astrophysical tests such as a collision between two stars to demonstrate
its robustness and excellent conservation properties.
[6]
oai:arXiv.org:1908.11258 [pdf] - 2011916
Rotating neutron stars with non-barotropic thermal profile
Submitted: 2019-08-29, last modified: 2019-11-12
Neutron stars provide an excellent laboratory for physics under the most
extreme conditions. Up to now, models of axisymmetric, stationary,
differentially rotating neutron stars were constructed under the strong
assumption of barotropicity, where a one-to-one relation between all
thermodynamic quantities exists. This implies that the specific angular
momentum of a matter element depends only on its angular velocity. The physical
conditions in the early stages of neutron stars, however, are determined by
their violent birth processes, typically a supernova or in some cases the
merger of two neutron stars, and detailed numerical models show that the
resulting stars are by no means barotropic. Here, we construct models for
stationary, differentially rotating, non-barotropic neutron stars, where the
equation of state and the specific angular momentum depend on more than one
independent variable. We show that the potential formulation of the
relativistic Euler equation can be extended to the non-barotropic case, which,
to the best of our knowledge, is a new result even for the Newtonian case. We
implement the new method into the XNS code and construct equilibrium
configurations for non-barotropic equations of state. We scrutinize the
resulting configurations by evolving them dynamically with the numerical
relativity code BAM, thereby demonstrating that the new method indeed produces
stationary, differentially rotating, non-barotropic neutron star
configurations.
[7]
oai:arXiv.org:1906.11494 [pdf] - 1994076
A multi-dimensional implementation of the Advanced Spectral neutrino
Leakage scheme
Submitted: 2019-06-27, last modified: 2019-10-15
We present a new, multi-dimensional implementation of the Advanced Spectral
Leakage (ASL) scheme with the purpose of modelling neutrino-matter interactions
in neutron star mergers. A major challenge is the neutrino absorption in the
semi-transparent regime, which is responsible for driving winds from the merger
remnant. The composition of such winds is crucial in the understanding of the
electromagnetic emission in the recently observed macronova following GW170817.
Compared to the original version, we introduce an optical-depth-dependent flux
factor to model the average angle of neutrino propagation, and a modulation
that accounts for flux anisotropies in non-spherical geometries. We scrutinise
our approach by first comparing the new scheme against the original one for a
spherically symmetric core-collapse supernova snapshot, both in 1D and in 3D,
and additionally against a two-moment (M1) scheme as implemented in 1D into the
code GR1D. The luminosities and mean energies agree to a few percents in most
tests. Finally, for the case of a binary merger remnant snapshot we compare the
new ASL scheme with the M1 scheme that is implemented in the Eulerian adaptive
mesh refinement code FLASH. We find that the neutrino absorption distribution
in the semi-transparent regime is overall well reproduced. Both approaches
agree to within $\lesssim 15\%$ for the average energies and to better than
$\sim 35 \%$ in the total luminosities.
[8]
oai:arXiv.org:1905.02159 [pdf] - 1966751
Short GRB 160821B: a reverse shock, a refreshed shock, and a
well-sampled kilonova
Lamb, G. P.;
Tanvir, N. R.;
Levan, A. J.;
Postigo, A. de Ugarte;
Kawaguchi, K.;
Corsi, A.;
Evans, P. A.;
Gompertz, B.;
Malesani, D. B.;
Page, K. L.;
Wiersema, K.;
Rosswog, S.;
Shibata, M.;
Tanaka, M.;
van der Horst, A. J.;
Cano, Z.;
Fynbo, J. P. U.;
Fruchter, A. S.;
Greiner, J.;
Heintz, K.;
Higgins, A.;
Hjorth, J.;
Izzo, L.;
Jakobsson, P.;
Kann, D. A.;
O'Brien, P. T.;
Perley, D. A.;
Pian, E.;
Pugliese, G.;
Starling, R. L. C.;
Thöne, C. C.;
Watson, D.;
Wijers, R. A. M. J.;
Xu, D.
Submitted: 2019-05-06, last modified: 2019-08-05
We report our identification of the optical afterglow and host galaxy of the
short-duration gamma-ray burst GRB 160821B. The spectroscopic redshift of the
host is $z=0.162$, making it one of the lowest redshift sGRBs identified by
Swift. Our intensive follow-up campaign using a range of ground-based
facilities as well as HST, XMM and Swift, shows evidence for a late-time excess
of optical and near-infrared emission in addition to a complex afterglow. The
afterglow light-curve at X-ray frequencies reveals a narrow jet,
$\theta_j\sim1.9^{+0.10}_{-0.03}$ deg, that is refreshed at $>1$ day post-burst
by a slower outflow with significantly more energy than the initial outflow
that produced the main GRB. Observations of the 5 GHz radio afterglow shows a
reverse shock into a mildly magnetised shell. The optical and near-infrared
excess is fainter than AT2017gfo associated with GW170817, and is well
explained by a kilonova with dynamic ejecta mass $M_{\rm
dyn}=(1.0\pm0.6)\times10^{-3}$ M$_{\odot}$ and a secular (postmerger) ejecta
mass with $M_{\rm pm}=(1.0\pm0.6)\times10^{-2}$ M$_\odot$, consistent with a
binary neutron star merger resulting in a short-lived massive neutron star.
This optical and near-infrared dataset provides the best-sampled kilonova
light-curve without a gravitational wave trigger to date.
[9]
oai:arXiv.org:1907.10091 [pdf] - 1925167
An unusual transient following the short GRB 071227
Eyles, R. A. J.;
O'Brien, P. T.;
Wiersema, K.;
Starling, R. L. C.;
Gompertz, B. P.;
Lamb, G. P.;
Lyman, J. D.;
Levan, A. J.;
Rosswog, S.;
Tanvir, N. R.
Submitted: 2019-07-23
We present X-ray and optical observations of the short duration gamma-ray
burst GRB 071227 and its host at $z=0.381$, obtained using \textit{Swift},
Gemini South and the Very Large Telescope. We identify a short-lived and
moderately bright optical transient, with flux significantly in excess of that
expected from a simple extrapolation of the X-ray spectrum at 0.2-0.3 days
after burst. We fit the SED with afterglow models allowing for high extinction
and thermal emission models that approximate a kilonova to assess the excess'
origins. While some kilonova contribution is plausible, it is not favoured due
to the low temperature and high luminosity required, implying superluminal
expansion and a large ejecta mass of $\sim 0.1$ M$_{\odot}$. We find, instead,
that the transient is broadly consistent with power-law spectra with additional
dust extinction of $E(B-V)\sim0.4$ mag, although a possibly thermal excess
remains in the \textit{z}-band. We investigate the host, a spiral galaxy with
an edge-on orientation, resolving its spectrum along its major axis to
construct the galaxy rotation curve and analyse the star formation and chemical
properties. The integrated host emission shows evidence for high extinction,
consistent with the afterglow findings. The metallicity and extinction are
consistent with previous studies of this host and indicate the galaxy is a
typical, but dusty, late-type SGRB host.
[10]
oai:arXiv.org:1903.09147 [pdf] - 1909276
Tidal disruptions by rotating black holes: effects of spin and impact
parameter
Submitted: 2019-03-21, last modified: 2019-07-02
We present the results of relativistic smoothed particle hydrodynamics
simulations of tidal disruptions of stars by rotating supermassive black holes,
for a wide range of impact parameters and black hole spins. For deep
encounters, we find that: relativistic precession creates debris geometries
impossible to obtain with the Newtonian equations; part of the fluid can be
launched on plunging orbits, reducing the fallback rate and the mass of the
resulting accretion disc; multiple squeezings and bounces at periapsis may
generate distinctive X-ray signatures resulting from the associated shock
breakout; disruptions can occur inside the marginally bound radius, if the
angular momentum spread launches part of the debris on non-plunging orbits.
Perhaps surprisingly, we also find relativistic effects important in partial
disruptions, where the balance between self-gravity and tidal forces is so
precarious that otherwise minor relativistic effects can have decisive
consequences on the stellar fate. In between, where the star is fully disrupted
but relativistic effects are mild, the difference resides in a gentler rise of
the fallback rate, a later and smaller peak, and longer return times. However,
relativistic precession always causes thicker debris streams, both in the bound
part (speeding up circularization) and in the unbound part (accelerating and
enhancing the production of separate transients). We discuss various properties
of the disruption (compression at periapsis, shape and spread of the energy
distribution) and potential observables (peak fallback rate, times of rise and
decay, duration of super-Eddington fallback) as a function of the impact
parameter and the black hole spin.
[11]
oai:arXiv.org:1811.05506 [pdf] - 1896980
A luminosity distribution for kilonovae based on short gamma-ray burst
afterglows
Ascenzi, Stefano;
Coughlin, Michael W.;
Dietrich, Tim;
Foley, Ryan J.;
Ramirez-Ruiz, Enrico;
Piranomonte, Silvia;
Mockler, Brenna;
Murguia-Berthier, Ariadna;
Fryer, Chris L.;
Lloyd-Ronning, Nicole M.;
Rosswog, Stephan
Submitted: 2018-11-13, last modified: 2019-06-09
The combined detection of a gravitational-wave signal, kilonova, and short
gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the
field of multi-messenger astronomy. But even before GW170817, there have been a
number of sGRBs with possible associated kilonova detections. In this work, we
re-examine these "historical" sGRB afterglows with a combination of
state-of-the-art afterglow and kilonova models. This allows us to include
optical/near-infrared synchrotron emission produced by the sGRB as well as
ultraviolet/optical/near-infrared emission powered by the radioactive decay of
$r$-process elements (i.e., the kilonova). Fitting the lightcurves, we derive
the velocity and the mass distribution as well as the composition of the
ejected material. The posteriors on kilonova parameters obtained from the fit
were turned into distributions for the peak magnitude of the kilonova emission
in different bands and the time at which this peak occurs. From the sGRB with
an associated kilonova, we found that the peak magnitude in H bands falls in
the range [-16.2, -13.1] ($95\%$ of confidence) and occurs within $0.8-3.6\,\rm
days$ after the sGRB prompt emission. In g band instead we obtain a peak
magnitude in range [-16.8, -12.3] occurring within the first $18\,\rm hr$ after
the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo,
kilonova candidates GRB130603B, GRB050709 and GRB060614 (with the possible
inclusion of GRB150101B) and the upper limits from all the other sGRBs not
associated with any kilonova detection we obtain for the first time a kilonova
luminosity function in different bands.
[12]
oai:arXiv.org:1903.09224 [pdf] - 1854203
The Yet-Unobserved Multi-Messenger Gravitational-Wave Universe
Kalogera, Vassiliki;
Bizouard, Marrie-Anne;
Burrows, Adam;
Janka, Thomas;
Kotake, Kei;
Messer, Bronson;
Mezzacappa, Tony;
Mueller, Bernhard;
Mueller, Ewald;
Papa, Maria Alessandra;
Reddy, Sanjay;
Rosswog, Stephan
Submitted: 2019-03-21
Observations with next-generation ground-based detectors further enhanced
with multi-messenger (electromagnetic and neutrino) detections will allow us to
probe new extreme astrophysics. Target sources included: core-collapse
supernovae, continuous emission from isolated or accreting neutron stars, and
bursts from magnetars and other pulsars.
[13]
oai:arXiv.org:1903.04553 [pdf] - 1847028
Gravity and Light: Combining Gravitational Wave and Electromagnetic
Observations in the 2020s
Foley, R. J.;
Alexander, K. D.;
Andreoni, I.;
Arcavi, I.;
Auchettl, K.;
Barnes, J.;
Baym, G.;
Bellm, E. C.;
Beloborodov, A. M.;
Blagorodnova, N.;
Blakeslee, J. P.;
Brady, P. R.;
Branchesi, M.;
Brown, J. S.;
Butler, N.;
Cantiello, M.;
Chornock, R.;
Cook, D. O.;
Cooke, J.;
Coppejans, D. L.;
Corsi, A.;
Couch, S. M.;
Coughlin, M. W.;
Coulter, D. A.;
Cowperthwaite, P. S.;
Dietrich, T.;
Dimitriadis, G.;
Drout, M. R.;
Elias, J. H.;
Farr, B.;
Fernandez, R.;
Filippenko, A. V.;
Fong, W.;
Fragos, T.;
Frail, D. A.;
Freedman, W. L.;
Fryer, C. L.;
Golkhou, V. Z.;
Hiramatsu, D.;
Hjorth, J.;
Horesh, A.;
Hosseinzadeh, G.;
Hotokezaka, K.;
Howell, D. A.;
Hung, T.;
Jones, D. O.;
Kalogera, V.;
Kasen, D.;
Kerzendorf, W. E.;
Kilpatrick, C. D.;
Kirshner, R. P.;
Krisciunas, K.;
Lattimer, J. M.;
Lazzati, D.;
Levan, A. J.;
MacFadyen, A. I.;
Maeda, K.;
Mandel, I.;
Mandel, K. S.;
Margalit, B.;
Margutti, R.;
McIver, J.;
Metzger, B. D.;
Mooley, K.;
Moriya, T.;
Murguia-Berthier, A.;
Narayan, G.;
Nicholl, M.;
Nissanke, S.;
Nomoto, K.;
O'Meara, J. M.;
O'Shaughnessy, R.;
O'Connor, E.;
Palmese, A.;
Pan, Y. -C.;
Pankow, C.;
Paterson, K.;
Perley, D. A.;
Perna, R.;
Piro, A. L.;
Pritchard, T. A.;
Quataert, E.;
Radice, D.;
Ramirez-Ruiz, E.;
Reddy, S.;
Rest, A.;
Riess, A. G.;
Rodriguez, C. L.;
Rojas-Bravo, C.;
Rossi, E. M.;
Rosswog, S.;
Ruiz, M.;
Shapiro, S. L.;
Shoemaker, D. H.;
Siebert, M. R.;
Siegel, D. M.;
Siellez, K.;
Smith, N.;
Soares-Santos, M.;
Suntzeff, N. B.;
Surman, R.;
Tanaka, M.;
Tanvir, N. R.;
Terreran, G.;
Valenti, S.;
Villar, V. A.;
Wang, L.;
Webb, S. A.;
Wheeler, J. C.;
Williams, P. K. G.;
Woosley, S.;
Zaldarriaga, M.;
Zevin, M.
Submitted: 2019-03-11
As of today, we have directly detected exactly one source in both
gravitational waves (GWs) and electromagnetic (EM) radiation, the binary
neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and
the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect
hundreds of events, including new classes of events such as
neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly
something completely unexpected. As we build this sample, we will explore
exotic astrophysical topics ranging from nucleosynthesis, stellar evolution,
general relativity, high-energy astrophysics, nuclear matter, to cosmology. The
discovery potential is extraordinary, and investments in this area will yield
major scientific breakthroughs. Here we outline some of the most exciting
scientific questions that can be answered by combining GW and EM observations.
[14]
oai:arXiv.org:1812.10492 [pdf] - 1833225
Serendipitous Discoveries of Kilonovae in the LSST Main Survey:
Maximising Detections of Sub-Threshold Gravitational Wave Events
Submitted: 2018-12-26, last modified: 2019-02-18
We investigate the ability of the Large Synoptic Survey Telescope (LSST) to
discover kilonovae (kNe) from binary neutron star (BNS) and neutron star-black
hole (NSBH) mergers, focusing on serendipitous detections in the Wide-Fast-Deep
(WFD) survey. We simulate observations of kNe with proposed LSST survey
strategies, paying particular attention to cadence choices that are compatible
with the broader LSST cosmology programme. We find that if all kNe are
identical to GW170817, the baseline survey strategy will yield 58 kNe over the
survey lifetime. If we instead assume a representative population model of BNS
kNe, we expect to detect only 27 kNe. However, we find the choice of survey
strategy significantly impacts these numbers and can increase them to 254 kNe
and 82 kNe over the survey lifetime, respectively. This improvement arises from
an increased cadence of observations between different filters with respect to
the baseline. We then consider the ability of the Advanced LIGO/Virgo (ALV)
detector network to detect these BNS mergers. If the optimal survey strategy is
adopted, 202 of the GW170817-like kNe and 56 of the BNS population model kNe
are detected with LSST but are below the threshold for detection by the ALV
network. This represents, for both models, an increase by a factor greater than
4.5 in the number of detected sub-threshold events over the baseline survey
strategy. Such a population of sub-threshold events would provide an
opportunity to conduct electromagnetic-triggered searches for signals in
gravitational-wave detector data and assess selection effects in measurements
of the Hubble constant from standard sirens, e.g., related to viewing angle
effects.
[15]
oai:arXiv.org:1902.02915 [pdf] - 1828935
Catching Element Formation In The Act
Fryer, Chris L.;
Timmes, Frank;
Hungerford, Aimee L.;
Couture, Aaron;
Adams, Fred;
Aoki, Wako;
Arcones, Almudena;
Arnett, David;
Auchettl, Katie;
Avila, Melina;
Badenes, Carles;
Baron, Eddie;
Bauswein, Andreas;
Beacom, John;
Blackmon, Jeff;
Blondin, Stephane;
Bloser, Peter;
Boggs, Steve;
Boss, Alan;
Brandt, Terri;
Bravo, Eduardo;
Brown, Ed;
Brown, Peter;
Budtz-Jorgensen, Steve Bruenn. Carl;
Burns, Eric;
Calder, Alan;
Caputo, Regina;
Champagne, Art;
Chevalier, Roger;
Chieffi, Alessandro;
Chipps, Kelly;
Cinabro, David;
Clarkson, Ondrea;
Clayton, Don;
Coc, Alain;
Connolly, Devin;
Conroy, Charlie;
Cote, Benoit;
Couch, Sean;
Dauphas, Nicolas;
deBoer, Richard James;
Deibel, Catherine;
Denisenkov, Pavel;
Desch, Steve;
Dessart, Luc;
Diehl, Roland;
Doherty, Carolyn;
Dominguez, Inma;
Dong, Subo;
Dwarkadas, Vikram;
Fan, Doreen;
Fields, Brian;
Fields, Carl;
Filippenko, Alex;
Fisher, Robert;
Foucart, Francois;
Fransson, Claes;
Frohlich, Carla;
Fuller, George;
Gibson, Brad;
Giryanskaya, Viktoriya;
Gorres, Joachim;
Goriely, Stephane;
Grebenev, Sergei;
Grefenstette, Brian;
Grohs, Evan;
Guillochon, James;
Harpole, Alice;
Harris, Chelsea;
Harris, J. Austin;
Harrison, Fiona;
Hartmann, Dieter;
Hashimoto, Masa-aki;
Heger, Alexander;
Hernanz, Margarita;
Herwig, Falk;
Hirschi, Raphael;
Hix, Raphael William;
Hoflich, Peter;
Hoffman, Robert;
Holcomb, Cole;
Hsiao, Eric;
Iliadis, Christian;
Janiuk, Agnieszka;
Janka, Thomas;
Jerkstrand, Anders;
Johns, Lucas;
Jones, Samuel;
Jose, Jordi;
Kajino, Toshitaka;
Karakas, Amanda;
Karpov, Platon;
Kasen, Dan;
Kierans, Carolyn;
Kippen, Marc;
Korobkin, Oleg;
Kobayashi, Chiaki;
Kozma, Cecilia;
Krot, Saha;
Kumar, Pawan;
Kuvvetli, Irfan;
Laird, Alison;
Laming, Martin;
Larsson, Josefin;
Lattanzio, John;
Lattimer, James;
Leising, Mark;
Lennarz, Annika;
Lentz, Eric;
Limongi, Marco;
Lippuner, Jonas;
Livne, Eli;
Lloyd-Ronning, Nicole;
Longland, Richard;
Lopez, Laura A.;
Lugaro, Maria;
Lutovinov, Alexander;
Madsen, Kristin;
Malone, Chris;
Matteucci, Francesca;
McEnery, Julie;
Meisel, Zach;
Messer, Bronson;
Metzger, Brian;
Meyer, Bradley;
Meynet, Georges;
Mezzacappa, Anthony;
Miller, Jonah;
Miller, Richard;
Milne, Peter;
Misch, Wendell;
Mitchell, Lee;
Mosta, Philipp;
Motizuki, Yuko;
Muller, Bernhard;
Mumpower, Matthew;
Murphy, Jeremiah;
Nagataki, Shigehiro;
Nakar, Ehud;
Nomoto, Ken'ichi;
Nugent, Peter;
Nunes, Filomena;
O'Shea, Brian;
Oberlack, Uwe;
Pain, Steven;
Parker, Lucas;
Perego, Albino;
Pignatari, Marco;
Pinedo, Gabriel Martinez;
Plewa, Tomasz;
Poznanski, Dovi;
Priedhorsky, William;
Pritychenko, Boris;
Radice, David;
Ramirez-Ruiz, Enrico;
Rauscher, Thomas;
Reddy, Sanjay;
Rehm, Ernst;
Reifarth, Rene;
Richman, Debra;
Ricker, Paul;
Rijal, Nabin;
Roberts, Luke;
Ropke, Friedrich;
Rosswog, Stephan;
Ruiter, Ashley J.;
Ruiz, Chris;
Savin, Daniel Wolf;
Schatz, Hendrik;
Schneider, Dieter;
Schwab, Josiah;
Seitenzahl, Ivo;
Shen, Ken;
Siegert, Thomas;
Sim, Stuart;
Smith, David;
Smith, Karl;
Smith, Michael;
Sollerman, Jesper;
Sprouse, Trevor;
Spyrou, Artemis;
Starrfield, Sumner;
Steiner, Andrew;
Strong, Andrew W.;
Sukhbold, Tuguldur;
Suntzeff, Nick;
Surman, Rebecca;
Tanimori, Toru;
The, Lih-Sin;
Thielemann, Friedrich-Karl;
Tolstov, Alexey;
Tominaga, Nozomu;
Tomsick, John;
Townsley, Dean;
Tsintari, Pelagia;
Tsygankov, Sergey;
Vartanyan, David;
Venters, Tonia;
Vestrand, Tom;
Vink, Jacco;
Waldman, Roni;
Wang, Lifang;
Wang, Xilu;
Warren, MacKenzie;
West, Christopher;
Wheeler, J. Craig;
Wiescher, Michael;
Winkler, Christoph;
Winter, Lisa;
Wolf, Bill;
Woolf, Richard;
Woosley, Stan;
Wu, Jin;
Wrede, Chris;
Yamada, Shoichi;
Young, Patrick;
Zegers, Remco;
Zingale, Michael;
Zwart, Simon Portegies
Submitted: 2019-02-07
Gamma-ray astronomy explores the most energetic photons in nature to address
some of the most pressing puzzles in contemporary astrophysics. It encompasses
a wide range of objects and phenomena: stars, supernovae, novae, neutron stars,
stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays
and relativistic-particle acceleration, and the evolution of galaxies. MeV
gamma-rays provide a unique probe of nuclear processes in astronomy, directly
measuring radioactive decay, nuclear de-excitation, and positron annihilation.
The substantial information carried by gamma-ray photons allows us to see
deeper into these objects, the bulk of the power is often emitted at gamma-ray
energies, and radioactivity provides a natural physical clock that adds unique
information. New science will be driven by time-domain population studies at
gamma-ray energies. This science is enabled by next-generation gamma-ray
instruments with one to two orders of magnitude better sensitivity, larger sky
coverage, and faster cadence than all previous gamma-ray instruments. This
transformative capability permits: (a) the accurate identification of the
gamma-ray emitting objects and correlations with observations taken at other
wavelengths and with other messengers; (b) construction of new gamma-ray maps
of the Milky Way and other nearby galaxies where extended regions are
distinguished from point sources; and (c) considerable serendipitous science of
scarce events -- nearby neutron star mergers, for example. Advances in
technology push the performance of new gamma-ray instruments to address a wide
set of astrophysical questions.
[16]
oai:arXiv.org:1806.05195 [pdf] - 1916735
Black holes, gravitational waves and fundamental physics: a roadmap
Barack, Leor;
Cardoso, Vitor;
Nissanke, Samaya;
Sotiriou, Thomas P.;
Askar, Abbas;
Belczynski, Krzysztof;
Bertone, Gianfranco;
Bon, Edi;
Blas, Diego;
Brito, Richard;
Bulik, Tomasz;
Burrage, Clare;
Byrnes, Christian T.;
Caprini, Chiara;
Chernyakova, Masha;
Chrusciel, Piotr;
Colpi, Monica;
Ferrari, Valeria;
Gaggero, Daniele;
Gair, Jonathan;
Garcia-Bellido, Juan;
Hassan, S. F.;
Heisenberg, Lavinia;
Hendry, Martin;
Heng, Ik Siong;
Herdeiro, Carlos;
Hinderer, Tanja;
Horesh, Assaf;
Kavanagh, Bradley J.;
Kocsis, Bence;
Kramer, Michael;
Tiec, Alexandre Le;
Mingarelli, Chiara;
Nardini, Germano;
Nelemans, Gijs;
Palenzuela, Carlos;
Pani, Paolo;
Perego, Albino;
Porter, Edward K.;
Rossi, Elena M.;
Schmidt, Patricia;
Sesana, Alberto;
Sperhake, Ulrich;
Stamerra, Antonio;
Stein, Leo C.;
Tamanini, Nicola;
Tauris, Thomas M.;
Urena-Lopez, L. Arturo;
Vincent, Frederic;
Volonteri, Marta;
Wardell, Barry;
Wex, Norbert;
Yagi, Kent;
Abdelsalhin, Tiziano;
Aloy, Miguel Angel;
Amaro-Seoane, Pau;
Annulli, Lorenzo;
Arca-Sedda, Manuel;
Bah, Ibrahima;
Barausse, Enrico;
Barakovic, Elvis;
Benkel, Robert;
Bennett, Charles L.;
Bernard, Laura;
Bernuzzi, Sebastiano;
Berry, Christopher P. L.;
Berti, Emanuele;
Bezares, Miguel;
Blanco-Pillado, Jose Juan;
Blazquez-Salcedo, Jose Luis;
Bonetti, Matteo;
Boskovic, Mateja;
Bosnjak, Zeljka;
Bricman, Katja;
Bruegmann, Bernd;
Capelo, Pedro R.;
Carloni, Sante;
Cerda-Duran, Pablo;
Charmousis, Christos;
Chaty, Sylvain;
Clerici, Aurora;
Coates, Andrew;
Colleoni, Marta;
Collodel, Lucas G.;
Compere, Geoffrey;
Cook, William;
Cordero-Carrion, Isabel;
Correia, Miguel;
de la Cruz-Dombriz, Alvaro;
Czinner, Viktor G.;
Destounis, Kyriakos;
Dialektopoulos, Kostas;
Doneva, Daniela;
Dotti, Massimo;
Drew, Amelia;
Eckner, Christopher;
Edholm, James;
Emparan, Roberto;
Erdem, Recai;
Ferreira, Miguel;
Ferreira, Pedro G.;
Finch, Andrew;
Font, Jose A.;
Franchini, Nicola;
Fransen, Kwinten;
Gal'tsov, Dmitry;
Ganguly, Apratim;
Gerosa, Davide;
Glampedakis, Kostas;
Gomboc, Andreja;
Goobar, Ariel;
Gualtieri, Leonardo;
Guendelman, Eduardo;
Haardt, Francesco;
Harmark, Troels;
Hejda, Filip;
Hertog, Thomas;
Hopper, Seth;
Husa, Sascha;
Ihanec, Nada;
Ikeda, Taishi;
Jaodand, Amruta;
Jimenez-Forteza, Philippe Jetzer Xisco;
Kamionkowski, Marc;
Kaplan, David E.;
Kazantzidis, Stelios;
Kimura, Masashi;
Kobayashi, Shiho;
Kokkotas, Kostas;
Krolik, Julian;
Kunz, Jutta;
Lammerzahl, Claus;
Lasky, Paul;
Lemos, Jose P. S.;
Said, Jackson Levi;
Liberati, Stefano;
Lopes, Jorge;
Luna, Raimon;
Ma, Yin-Zhe;
Maggio, Elisa;
Montero, Marina Martinez;
Maselli, Andrea;
Mayer, Lucio;
Mazumdar, Anupam;
Messenger, Christopher;
Menard, Brice;
Minamitsuji, Masato;
Moore, Christopher J.;
Mota, David;
Nampalliwar, Sourabh;
Nerozzi, Andrea;
Nichols, David;
Nissimov, Emil;
Obergaulinger, Martin;
Obers, Niels A.;
Oliveri, Roberto;
Pappas, George;
Pasic, Vedad;
Peiris, Hiranya;
Petrushevska, Tanja;
Pollney, Denis;
Pratten, Geraint;
Rakic, Nemanja;
Racz, Istvan;
Radia, Miren;
Ramazanouglu, Fethi M.;
Ramos-Buades, Antoni;
Raposo, Guilherme;
Rosca-Mead, Roxana;
Rogatko, Marek;
Rosinska, Dorota;
Rosswog, Stephan;
Morales, Ester Ruiz;
Sakellariadou, Mairi;
Sanchis-Gual, Nicolas;
Salafia, Om Sharan;
Samajdar, Anuradha;
Sintes, Alicia;
Smole, Majda;
Sopuerta, Carlos;
Souza-Lima, Rafael;
Stalevski, Marko;
Stergioulas, Nikolaos;
Stevens, Chris;
Tamfal, Tomas;
Torres-Forne, Alejandro;
Tsygankov, Sergey;
Unluturk, Kivanc;
Valiante, Rosa;
van de Meent, Maarten;
Velhinho, Jose;
Verbin, Yosef;
Vercnocke, Bert;
Vernieri, Daniele;
Vicente, Rodrigo;
Vitagliano, Vincenzo;
Weltman, Amanda;
Whiting, Bernard;
Williamson, Andrew;
Witek, Helvi;
Wojnar, Aneta;
Yakut, Kadri;
Yan, Haopeng;
Yazadjiev, Stoycho;
Zaharijas, Gabrijela;
Zilhao, Miguel
Submitted: 2018-06-13, last modified: 2019-02-01
The grand challenges of contemporary fundamental physics---dark matter, dark
energy, vacuum energy, inflation and early universe cosmology, singularities
and the hierarchy problem---all involve gravity as a key component. And of all
gravitational phenomena, black holes stand out in their elegant simplicity,
while harbouring some of the most remarkable predictions of General Relativity:
event horizons, singularities and ergoregions. The hitherto invisible landscape
of the gravitational Universe is being unveiled before our eyes: the historical
direct detection of gravitational waves by the LIGO-Virgo collaboration marks
the dawn of a new era of scientific exploration. Gravitational-wave astronomy
will allow us to test models of black hole formation, growth and evolution, as
well as models of gravitational-wave generation and propagation. It will
provide evidence for event horizons and ergoregions, test the theory of General
Relativity itself, and may reveal the existence of new fundamental fields. The
synthesis of these results has the potential to radically reshape our
understanding of the cosmos and of the laws of Nature. The purpose of this work
is to present a concise, yet comprehensive overview of the state of the art in
the relevant fields of research, summarize important open problems, and lay out
a roadmap for future progress.
[17]
oai:arXiv.org:1812.08708 [pdf] - 1818826
Spitzer Mid-Infrared Detections of Neutron Star Merger GW170817 Suggests
Synthesis of the Heaviest Elements
Kasliwal, Mansi M.;
Kasen, Daniel;
Lau, Ryan M.;
Perley, Daniel A.;
Rosswog, Stephan;
Ofek, Eran O.;
Hotokezaka, Kenta;
Chary, Ranga-Ram;
Sollerman, Jesper;
Goobar, Ariel;
Kaplan, David L.
Submitted: 2018-12-20
We report our Spitzer Space Telescope observations and detections of the
binary neutron star merger GW170817. At 4.5um, GW170817 is detected at 21.9 mag
AB at +43 days and 23.9 mag AB at +74 days after merger. At 3.6um, GW170817 is
not detected to a limit of 23.2 mag AB at +43 days and 23.1 mag AB at +74 days.
Our detections constitute the latest and reddest constraints on the
kilonova/macronova emission and composition of heavy elements. The 4.5um
luminosity at this late phase cannot be explained by elements exclusively from
the first abundance peak of the r-process. Moreover, the steep decline in the
Spitzer band, with a power-law index of 3.4 +/- 0.2, can be explained by a few
of the heaviest isotopes in the third abundance peak with half-life around 14
days dominating the luminosity (e.g. 140Ba, 143Pr, 147Nd, 156Eu, 191Os, 223Ra,
225Ra, 233Pa, 234Th) or a model with lower deposition efficiency. This data
offers evidence that the heaviest elements in the second and third r-process
abundance peak were indeed synthesized. Our conclusion is verified by both
analytics and network simulations and robust despite intricacies and
uncertainties in the nuclear physics. Future observations with Spitzer and
James Webb Space Telescope will further illuminate the relative abundance of
the synthesized heavy elements.
[18]
oai:arXiv.org:1801.02669 [pdf] - 1753218
The optical afterglow of the short gamma-ray burst associated with
GW170817
Lyman, J. D.;
Lamb, G. P.;
Levan, A. J.;
Mandel, I.;
Tanvir, N. R.;
Kobayashi, S.;
Gompertz, B.;
Hjorth, J.;
Fruchter, A. S.;
Kangas, T.;
Steeghs, D.;
Steele, I. A.;
Cano, Z.;
Copperwheat, C.;
Evans, P. A.;
Fynbo, J. P. U.;
Gall, C.;
Im, M.;
Izzo, L.;
Jakobsson, P.;
Milvang-Jensen, B.;
O'Brien, P.;
Osborne, J. P.;
Palazzi, E.;
Perley, D. A.;
Pian, E.;
Rosswog, S.;
Rowlinson, A.;
Schulze, S.;
Stanway, E. R.;
Sutton, P.;
Thöne, C. C.;
Postigo, A. de Ugarte;
Watson, D. J.;
Wiersema, K.;
Wijers, R. A. M. J.
Submitted: 2018-01-08, last modified: 2018-09-19
The binary neutron star merger GW170817 was the first multi-messenger event
observed in both gravitational and electromagnetic waves. The electromagnetic
signal began approximately 2 seconds post-merger with a weak, short burst of
gamma-rays, which was followed over the next hours and days by the ultraviolet,
optical and near-infrared emission from a radioactively- powered kilonova.
Later, non-thermal rising X-ray and radio emission was observed. The low
luminosity of the gamma-rays and the rising non-thermal flux from the source at
late times could indicate that we are outside the opening angle of the beamed
relativistic jet. Alternatively, the emission could be arising from a cocoon of
material formed from the interaction between a jet and the merger ejecta. Here
we present late-time optical detections and deep near-infrared limits on the
emission from GW170817 at 110 days post-merger. Our new observations are at
odds with expectations of late-time emission from kilonova models, being too
bright and blue. Instead, the emission arises from the interaction between the
relativistic ejecta of GW170817 and the interstellar medium. We show that this
emission matches the expectations of a Gaussian structured relativistic jet,
which would have launched a high luminosity short GRB to an aligned observer.
However, other jet structure or cocoon models can also match current data - the
future evolution of the afterglow will directly distinguish the origin of the
emission.
[19]
oai:arXiv.org:1806.07775 [pdf] - 1731108
Disk formation in the collapse of supramassive neutron stars
Submitted: 2018-06-20, last modified: 2018-08-10
Short gamma-ray bursts (sGRBs) show a large diversity in their properties.
This suggests that the observed phenomenon can be caused by different "central
engines" or that the engine produces a variety of outcomes depending on its
parameters, or possibly both. The most popular engine scenario, the merger of
two neutron stars, has received support from the recent Fermi and INTEGRAL
detection of a burst of gamma rays (GRB170817A) following the neutron star
merger GW170817, but at the moment it is not clear how peculiar this event
potentially was. Several sGRBs engine models involve the collapse of a
supramassive neutron star that produces a black hole plus an accretion disk. We
study this scenario for a variety of equations of states both via angular
momentum considerations based on equilibrium models and via fully dynamical
Numerical Relativity simulations. We obtain a broader range of disk forming
configurations than earlier studies but we agree with the latter that none of
these configurations is likely to produce a phenomenon that would be classified
as an sGRB.
[20]
oai:arXiv.org:1710.05445 [pdf] - 1728811
The first direct double neutron star merger detection: implications for
cosmic nucleosynthesis
Submitted: 2017-10-16, last modified: 2018-04-04
The astrophysical r-process site where about half of the elements heavier
than iron are produced has been a puzzle for several decades. Here we discuss
the role of neutron star mergers (NSMs) in the light of the first direct
detection of such an event in both gravitational (GW) and electromagnetic (EM)
waves. We analyse bolometric and NIR lightcurves of the first detected double
neutron star merger and compare them to nuclear reaction network-based
macronova models. The slope of the bolometric lightcurve is consistent with the
radioactive decay of neutron star ejecta with $Y_e \lesssim 0.3$ (but not
larger), which provides strong evidence for an r-process origin of the
electromagnetic emission. This rules out in particular "nickel winds" as major
source of the emission. We find that the NIR lightcurves can be well fitted
either with or without lanthanide-rich ejecta. Our limits on the ejecta mass
together with estimated rates directly confirm earlier purely theoretical or
indirect observational conclusions that double neutron star mergers are indeed
a major site of cosmic nucleosynthesis. If the ejecta mass was {\em typical},
NSMs can easily produce {\em all} of the estimated Galactic r-process matter,
and --depending on the real rate-- potentially even more. This could be a hint
that the event ejected a particularly large amount of mass, maybe due to a
substantial difference between the component masses. This would be compatible
with the mass limits obtained from the GW-observation. The recent observations
suggests that NSMs are responsible for a broad range of r-process nuclei and
that they are at least a major, but likely the dominant r-process site in the
Universe.
[21]
oai:arXiv.org:1710.05455 [pdf] - 1589771
The Emergence of a Lanthanide-Rich Kilonova Following the Merger of Two
Neutron Stars
Tanvir, N. R.;
Levan, A. J.;
Gonzalez-Fernandez, C.;
Korobkin, O.;
Mandel, I.;
Rosswog, S.;
Hjorth, J.;
D'Avanzo, P.;
Fruchter, A. S.;
Fryer, C. L.;
Kangas, T.;
Milvang-Jensen, B.;
Rosetti, S.;
Steeghs, D.;
Wollaeger, R. T.;
Cano, Z.;
Copperwheat, C. M.;
Covino, S.;
D'Elia, V.;
Postigo, A. de Ugarte;
Evans, P. A.;
Even, W. P.;
Fairhurst, S.;
Jaimes, R. Figuera;
Fontes, C. J.;
Fujii, Y. I.;
Fynbo, J. P. U.;
Gompertz, B. P.;
Greiner, J.;
Hodosan, G.;
Irwin, M. J.;
Jakobsson, P.;
Jorgensen, U. G.;
Kann, D. A.;
Lyman, J. D.;
Malesani, D.;
McMahon, R. G.;
Melandri, A.;
O'Brien, P. T.;
Osborne, J. P.;
Palazzi, E.;
Perley, D. A.;
Pian, E.;
Piranomonte, S.;
Rabus, M.;
Rol, E.;
Rowlinson, A.;
Schulze, S.;
Sutton, P.;
Thoene, C. C.;
Ulaczyk, K.;
Watson, D.;
Wiersema, K.;
Wijers, R. A. M. J.
Submitted: 2017-10-16
We report the discovery and monitoring of the near-infrared counterpart
(AT2017gfo) of a binary neutron-star merger event detected as a gravitational
wave source by Advanced LIGO/Virgo (GW170817) and as a short gamma-ray burst by
Fermi/GBM and Integral/SPI-ACS (GRB170817A). The evolution of the transient
light is consistent with predictions for the behaviour of a
"kilonova/macronova", powered by the radioactive decay of massive neutron-rich
nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In
particular, evidence for this scenario is found from broad features seen in
Hubble Space Telescope infrared spectroscopy, similar to those predicted for
lanthanide dominated ejecta, and the much slower evolution in the near-infrared
Ks-band compared to the optical. This indicates that the late-time light is
dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to
the 3rd r-process peak (atomic masses A~195). This discovery confirms that
neutron-star mergers produce kilo-/macronovae and that they are at least a
major - if not the dominant - site of rapid neutron capture nucleosynthesis in
the universe.
[22]
oai:arXiv.org:1710.05444 [pdf] - 1589762
The environment of the binary neutron star merger GW170817
Levan, A. J.;
Lyman, J. D.;
Tanvir, N. R.;
Hjorth, J.;
Mandel, I.;
Stanway, E. R.;
Steeghs, D.;
Fruchter, A. S.;
Troja, E.;
Schrøder, S. L;
Wiersema, K.;
Bruun, S. H.;
Cano, Z.;
Cenko, S. B.;
Postigo, A de Ugarte;
Evans, P.;
Fairhurst, S.;
Fox, O. D.;
Fynbo, J. P. U.;
Gompertz, B.;
Greiner, J.;
Im, M.;
Izzo, L.;
Jakobsson, P.;
Kangas, T.;
Khandrika, H. G.;
Lien, A. Y.;
Malesani, D.;
O'Brien, P.;
Osborne, J. P.;
Plazzi, E.;
Pian, E.;
Perley, D. A.;
Rosswog, S.;
Ryan, R. E.;
Schulze, S.;
Sutton, P.;
Thöne, C. C.;
Watson, D. J.;
Wijers, R. A. M. J.
Submitted: 2017-10-16
We present Hubble Space Telescope and Chandra imaging, combined with Very
Large Telescope MUSE integral field spectroscopy of the counterpart and host
galaxy of the first binary neutron star merger detected via gravitational wave
emission by LIGO & Virgo, GW170817. The host galaxy, NGC 4993, is an S0 galaxy
at z=0.009783. There is evidence for large, face-on spiral shells in continuum
imaging, and edge-on spiral features visible in nebular emission lines. This
suggests that NGC 4993 has undergone a relatively recent (<1 Gyr) ``dry''
merger. This merger may provide the fuel for a weak active nucleus seen in
Chandra imaging. At the location of the counterpart, HST imaging implies there
is no globular or young stellar cluster, with a limit of a few thousand solar
masses for any young system. The population in the vicinity is predominantly
old with <1% of any light arising from a population with ages <500 Myr. Both
the host galaxy properties and those of the transient location are consistent
with the distributions seen for short-duration gamma-ray bursts, although the
source position lies well within the effective radius (r_e ~ 3 kpc), providing
an r_e-normalized offset that is closer than ~90% of short GRBs. For the long
delay time implied by the stellar population, this suggests that the kick
velocity was significantly less than the galaxy escape velocity. We do not see
any narrow host galaxy interstellar medium features within the counterpart
spectrum, implying low extinction, and that the binary may lie in front of the
bulk of the host galaxy.
[23]
oai:arXiv.org:1710.05437 [pdf] - 1589755
Swift and NuSTAR observations of GW170817: detection of a blue kilonova
Evans, P. A.;
Cenko, S. B.;
Kennea, J. A.;
Emery, S. W. K.;
Kuin, N. P. M.;
Korobkin, O.;
Wollaeger, R. T.;
Fryer, C. L.;
Madsen, K. K.;
Harrison, F. A.;
Xu, Y.;
Nakar, E.;
Hotokezaka, K.;
Lien, A.;
Campana, S.;
Oates, S. R.;
Troja, E.;
Breeveld, A. A.;
Marshall, F. E.;
Barthelmy, S. D.;
Beardmore, A. P.;
Burrows, D. N.;
Cusumano, G.;
D'Ai, A.;
D'Avanzo, P.;
D'Elia, V.;
de Pasquale, M.;
Even, W. P.;
Fontes, C. J.;
Forster, K.;
Garcia, J.;
Giommi, P.;
Grefenstette, B.;
Gronwall, C.;
Hartmann, D. H.;
Heida, M.;
Hungerford, A. L.;
Kasliwal, M. M.;
Krimm, H. A.;
Levan, A. J.;
Malesani, D.;
Melandri, A.;
Miyasaka, H.;
Nousek, J. A.;
O'Brien, P. T.;
Osborne, J. P.;
Pagani, C.;
Page, K. L.;
Palmer, D. M.;
Perri, M.;
Pike, S.;
Racusin, J. L.;
Rosswog, S.;
Siegel, M. H.;
Sakamoto, T.;
Sbarufatti, B.;
Tagliaferri, G.;
Tanvir, N. R.;
Tohuvavohu, A.
Submitted: 2017-10-16
With the first direct detection of merging black holes in 2015, the era of
gravitational wave (GW) astrophysics began. A complete picture of compact
object mergers, however, requires the detection of an electromagnetic (EM)
counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the
Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the
binary neutron star merger GW170817. The bright, rapidly fading ultraviolet
emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow
with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray
limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from
the orbital rotation axis, which avoids both obscuration from the heaviest
elements in the orbital plane and a direct view of any ultra-relativistic,
highly collimated ejecta (a gamma-ray burst afterglow).
[24]
oai:arXiv.org:1710.05436 [pdf] - 1589754
Illuminating Gravitational Waves: A Concordant Picture of Photons from a
Neutron Star Merger
Kasliwal, M. M.;
Nakar, E.;
Singer, L. P.;
Kaplan, D. L.;
Cook, D. O.;
Van Sistine, A.;
Lau, R. M.;
Fremling, C.;
Gottlieb, O.;
Jencson, J. E.;
Adams, S. M.;
Feindt, U.;
Hotokezaka, K.;
Ghosh, S.;
Perley, D. A.;
Yu, P. -C.;
Piran, T.;
Allison, J. R.;
Anupama, G. C.;
Balasubramanian, A.;
Bannister, K. W;
Bally, J.;
Barnes, J.;
Barway, S.;
Bellm, E.;
Bhalerao, V.;
Bhattacharya, D.;
Blagorodnova, N.;
Bloom, J. S.;
Brady, P. R.;
Cannella, C.;
Chatterjee, D.;
Cenko, S. B.;
Cobb, B. E.;
Copperwheat, C.;
Corsi, A.;
De, K.;
Dobie, D.;
Emery, S. W. K.;
Evans, P. A.;
Fox, O. D.;
Frail, D. A.;
Frohmaier, C.;
Goobar, A.;
Hallinan, G.;
Harrison, F.;
Helou, G.;
Hinderer, T.;
Ho, A. Y. Q.;
Horesh, A.;
Ip, W. -H.;
Itoh, R.;
Kasen, D.;
Kim, H.;
Kuin, N. P. M.;
Kupfer, T.;
Lynch, C.;
Madsen, K.;
Mazzali, P. A.;
Miller, A. A.;
Mooley, K.;
Murphy, T.;
Ngeow, C. -C.;
Nichols, D.;
Nissanke, S.;
Nugent, P.;
Ofek, E. O.;
Qi, H.;
Quimby, R. M.;
Rosswog, S.;
Rusu, F.;
Sadler, E. M.;
Schmidt, P.;
Sollerman, J.;
Steele, I.;
Williamson, A. R.;
Xu, Y.;
Yan, L.;
Yatsu, Y.;
Zhang, C.;
Zhao, W.
Submitted: 2017-10-16
Merging neutron stars offer an exquisite laboratory for simultaneously
studying strong-field gravity and matter in extreme environments. We establish
the physical association of an electromagnetic counterpart EM170817 to
gravitational waves (GW170817) detected from merging neutron stars. By
synthesizing a panchromatic dataset, we demonstrate that merging neutron stars
are a long-sought production site forging heavy elements by r-process
nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to
classical short gamma-ray bursts with ultra-relativistic jets. Instead, we
suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the
jet elegantly explains the low-luminosity gamma-rays, the high-luminosity
ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit
that all merging neutron stars may lead to a wide-angle cocoon breakout;
sometimes accompanied by a successful jet and sometimes a choked jet.
[25]
oai:arXiv.org:1710.05863 [pdf] - 1589851
Lanthanides or dust in kilonovae: lessons learned from GW170817
Submitted: 2017-10-16
The unprecedented optical and near-infrared lightcurves of the first
electromagnetic counterpart to a gravitational wave source, GW170817, a binary
neutron star merger, exhibited a strong evolution from blue to near-infrared (a
so-called 'kilonova' or 'macronova'). The emerging near-infrared component is
widely attributed to the formation of r-process elements which provide the
opacity to shift the blue light into the near infrared. An alternative scenario
is that the light from the blue component gets extinguished by dust formed by
the kilonova and subsequently is re-emitted at near-infrared wavelengths. We
here test this hypothesis using the lightcurves of AT2017gfo, the kilonova
accompanying GW170817. We find that of order 10$^{-5}$ $M_\odot$ of carbon is
required to reproduce the optical/near-infrared lightcurves as the kilonova
fades. This putative dust cools from $\sim$ 2000 K at $\sim$ 4 d after the
event to $\sim$ 1500 K over the course of the following week, thus requiring
dust with a high condensation temperature, such as carbon. We contrast this
with the nucleosynthetic yields predicted by a range of kilonova wind models.
These suggest that at most 10$^{-9}$ $M_\odot$ of carbon is formed. Moreover,
the decay in the inferred dust temperature is slower than that expected in
kilonova models. We therefore conclude that in current models of the blue
component of the kilonova, the near-infrared component in the kilonova
accompanying GW170817 is unlikely to be due to dust.
[27]
oai:arXiv.org:1508.02399 [pdf] - 1579685
Optical Thermonuclear Transients From Tidal Compression of White Dwarfs
as Tracers of the Low End of the Massive Black Hole Mass Function
Submitted: 2015-08-10, last modified: 2017-06-14
In this paper, we model the observable signatures of tidal disruptions of
white dwarf (WD) stars by massive black holes (MBHs) of moderate mass, $\approx
10^3 - 10^5 M_\odot$. When the WD passes deep enough within the MBH's tidal
field, these signatures include thermonuclear transients from burning during
maximum compression. We combine a hydrodynamic simulation that includes nuclear
burning of the disruption of a $0.6 M_\odot$ C/O WD with a Monte Carlo
radiative transfer calculation to synthesize the properties of a representative
transient. The transient's emission emerges in the optical, with lightcurves
and spectra reminiscent of type I SNe. The properties are strongly
viewing-angle dependent, and key spectral signatures are $\approx 10,000$ km
s$^{-1}$ Doppler shifts due to the orbital motion of the unbound ejecta.
Disruptions of He WDs likely produce large quantities of intermediate-mass
elements, offering a possible production mechanism for Ca-rich transients.
Accompanying multiwavelength transients are fueled by accretion and arise from
the nascent accretion disk and relativistic jet. If MBHs of moderate mass exist
with number densities similar to those of supermassive MBHs, both high energy
wide-field monitors and upcoming optical surveys should detect tens to hundreds
of WD tidal disruptions per year. The current best strategy for their detection
may therefore be deep optical follow up of high-energy transients of
unusually-long duration. The detection rate or the non-detection of these
transients by current and upcoming surveys can thus be used to place meaningful
constraints on the extrapolation of the MBH mass function to moderate masses.
[28]
oai:arXiv.org:1701.00303 [pdf] - 1581043
Tidal disruptions by rotating black holes: relativistic hydrodynamics
with Newtonian codes
Submitted: 2017-01-01, last modified: 2017-05-24
We propose an approximate approach for studying the relativistic regime of
stellar tidal disruptions by rotating massive black holes. It combines an exact
relativistic description of the hydrodynamical evolution of a test fluid in a
fixed curved spacetime with a Newtonian treatment of the fluid's self-gravity.
Explicit expressions for the equations of motion are derived for Kerr spacetime
using two different coordinate systems. We implement the new methodology within
an existing Newtonian Smoothed Particle Hydrodynamics code and show that
including the additional physics involves very little extra computational cost.
We carefully explore the validity of the novel approach by first testing its
ability to recover geodesic motion, and then by comparing the outcome of tidal
disruption simulations against previous relativistic studies. We further
compare simulations in Boyer--Lindquist and Kerr--Schild coordinates and
conclude that our approach allows accurate simulation even of tidal disruption
events where the star penetrates deeply inside the tidal radius of a rotating
black hole. Finally, we use the new method to study the effect of the black
hole spin on the morphology and fallback rate of the debris streams resulting
from tidal disruptions, finding that while the spin has little effect on the
fallback rate, it does imprint heavily on the stream morphology, and can even
be a determining factor in the survival or disruption of the star itself. Our
methodology is discussed in detail as a reference for future astrophysical
applications.
[29]
oai:arXiv.org:1705.07084 [pdf] - 1674697
Impact of ejecta morphology and composition on the electromagnetic
signatures of neutron star mergers
Wollaeger, Ryan T.;
Korobkin, Oleg;
Fontes, Christopher J.;
Rosswog, Stephan K.;
Even, Wesley P.;
Fryer, Christopher L.;
Sollerman, Jesper;
Hungerford, Aimee L.;
van Rossum, Daniel R.;
Wollaber, Allan B.
Submitted: 2017-05-19
The electromagnetic transients accompanying compact binary mergers
($\gamma$-ray bursts, afterglows and 'macronovae') are crucial to pinpoint the
sky location of gravitational wave sources. Macronovae are caused by the
radioactivity from freshly synthesised heavy elements, e.g. from dynamic ejecta
and various types of winds. We study macronova signatures by using
multi-dimensional radiative transfer calculations. We employ the radiative
transfer code SuperNu and state-of-the-art LTE opacities for a few
representative elements from the wind and dynamical ejecta (Cr, Pd, Se, Te, Br,
Zr, Sm, Ce, Nd, U) to calculate synthetic light curves and spectra for a range
of ejecta morphologies. The radioactive power of the resulting macronova is
calculated with the detailed input of decay products. We assess the detection
prospects for our most complex models, based on the portion of viewing angles
that are sufficiently bright, at different cosmological redshifts ($z$). The
brighter emission from the wind is unobscured by the lanthanides (or actinides)
in some of the models, permitting non-zero detection probabilities for
redshifts up to $z=0.07$. We also find the nuclear mass model and the resulting
radioactive heating rate are crucial for the detectability. While for the most
pessimistic heating rate (from the FRDM model) no reasonable increase in the
ejecta mass or velocity, or wind mass or velocity, can possibly make the light
curves agree with the observed nIR excess after GRB130603B, a more optimistic
heating rate (from the Duflo-Zuker model) leads to good agreement. We conclude
that future reliable macronova observations would constrain nuclear heating
rates, and consequently help constrain nuclear mass models.
[30]
oai:arXiv.org:1702.02990 [pdf] - 1535649
A line-smeared treatment of opacities for the spectra and light curves
from macronovae
Submitted: 2017-02-09
Gravitational wave observations need accompanying electromagnetic signals to
accurately determine the sky positions of the sources. The ejecta of neutron
star mergers are expected to produce such electromagnetic transients, called
macronovae. Characteristics of the ejecta include large velocity gradients and
the presence of heavy $r$-process elements, which pose significant challenges
to the accurate calculation of radiative opacities and radiation transport. For
example, these opacities include a dense forest of bound-bound features arising
from near-neutral lanthanide and actinide elements. Here we investigate the use
of fine-structure, line-smeared opacities that preserve the integral of the
opacity over frequency, which is an alternative approach to the
expansion-opacity formalism. The use of area-preserving line profiles produces
frequency-dependent opacities that are one to two orders of magnitude greater
than those obtained with the use of expansion opacities in the Sobolev
approximation. Opacities are calculated for four $r$-process elements (cerium,
neodymium, samarium and uranium) using fully and semi-relativistic methods, as
well as different amounts of configuration interaction in the atomic structure,
in order to test the sensitivity of the emission to the underlying atomic
physics. We determine the effect of these new opacities on simulated spectra
and broad-band light curves by applying a multi-dimensional ray-trace method to
the ejecta predicted from 3-dimensional merger calculations. With our
substantially larger opacities the simulations yield slightly lower
luminosities that peak in the mid-IR. These results suggest that those
radioactively powered transients that are related to the very heaviest
$r$-process material are more difficult to observe than previously believed.
[31]
oai:arXiv.org:1611.09822 [pdf] - 1563970
Detectability of compact binary merger macronovae
Submitted: 2016-11-29
We study the optical and near-infrared luminosities and detectability of
radioactively powered electromagnetic transients ('macronovae') occuring in the
aftermath of binary neutron star and neutron star black hole mergers. We
explore the transients that result from the dynamic ejecta and those from
different types of wind outflows. Based on full nuclear network simulations we
calculate the resulting light curves in different wavelength bands. We
scrutinize the robustness of the results by comparing a) two different nuclear
reaction networks and b) two macronova models. We explore in particular how
sensitive the results are to the production of alpha-decaying trans-lead
nuclei. We compare two frequently used mass models: the Finite-Range Droplet
Model (FRDM) and the nuclear mass model of Duflo and Zuker (DZ31). We find that
the abundance of alpha-decaying trans-lead nuclei has a significant impact on
the observability of the resulting macronovae. For example, the DZ31 model
yields considerably larger abundances resulting in larger heating rates and
thermalization efficiencies and therefore predicts substantially brighter
macronova transients. We find that the dynamic ejecta from NSNS models can
reach peak K-band magnitudes in excess of -15 while those from NSBH cases can
reach beyond -16. Similar values can be reached by some of our wind models.
Several of our models (both wind and dynamic ejecta) yield properties that are
similar to the transient that was observed in the aftermath of the short GRB
130603B. We further explore the expected macronova detection frequencies for
current and future instruments such as VISTA, ZTF and LSST.
[32]
oai:arXiv.org:1609.04828 [pdf] - 1531614
The Properties of Short gamma-ray burst Jets Triggered by neutron star
mergers
Submitted: 2016-09-15
The most popular model for short gamma-ray bursts (sGRBs) involves the
coalescence of binary neutron stars. Because the progenitor is actually hidden
from view, we must consider under which circumstances such merging systems are
capable of producing a successful sGRB. Soon after coalescence, winds are
launched from the merger remnant. In this paper, we use realistic wind profiles
derived from global merger simulations in order to investigate the interaction
of sGRB jets with these winds using numerical simulations. We analyze the
conditions for which these axisymmetric winds permit relativistic jets to
breakout and produce a sGRB. We find that jets with luminosities comparable to
those observed in sGRBs are only successful when their half-opening angles are
below ~20{\deg}. This jet collimation mechanism leads to a simple physical
interpretation of the luminosities and opening angles inferred for sGRBs. If
wide, low luminosity jets are observed, they might be indicative of a different
progenitor avenue such as the merger of a neutron star with a black hole. We
also use the observed durations of sGRB to place constraints on the lifetime of
the wind phase, which is determined by the time it takes the jet to breakout.
In all cases we find that the derived limits argue against completely stable
remnants for binary neutron star mergers that produce sGRBs.
[35]
oai:arXiv.org:1512.04865 [pdf] - 1378744
Magnetohydrodynamical simulations of a tidal disruption in general
relativity
Submitted: 2015-12-15
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
$\beta=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.
[36]
oai:arXiv.org:1406.4224 [pdf] - 1292115
SPH Methods in the Modelling of Compact Objects
Submitted: 2014-06-16, last modified: 2015-10-13
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.
[37]
oai:arXiv.org:1508.02402 [pdf] - 1304052
Thermonuclear detonations ensuing white dwarf mergers
Submitted: 2015-08-10, last modified: 2015-10-01
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.
[38]
oai:arXiv.org:1506.05048 [pdf] - 1298131
Neutrino-driven winds in the aftermath of a neutron star merger:
nucleosynthesis and electromagnetic transients
Submitted: 2015-06-16, last modified: 2015-09-24
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}$.
[39]
oai:arXiv.org:1411.0974 [pdf] - 1258679
The Role of Fission in Neutron Star Mergers and its Impact on the
r-Process Peaks
Eichler, Marius;
Arcones, Almudena;
Kelic, Alexandra;
Korobkin, Oleg;
Langanke, Karlheinz;
Marketin, Tomislav;
Martinez-Pinedo, Gabriel;
Panov, Igor V.;
Rauscher, Thomas;
Rosswog, Stephan;
Winteler, Christian;
Zinner, Nikolaj T.;
Thielemann, Friedrich-Karl
Submitted: 2014-11-04, last modified: 2015-05-29
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
{\beta}-decays compete and an (n,{\gamma})-({\gamma},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 {\beta}-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.
[40]
oai:arXiv.org:1504.07605 [pdf] - 1292178
The Fate of the Compact Remnant in Neutron Star Mergers
Submitted: 2015-04-28
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.
[41]
oai:arXiv.org:1502.02039 [pdf] - 968972
Relativistic effects on tidal disruption kicks of solitary stars
Submitted: 2015-02-06, last modified: 2015-02-17
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 $\gamma$ = 5/3 and $\gamma$ = 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 $\beta$ alone, reaching at most the escape
velocity of the original star. For a given $\beta$, 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-$\beta$ encounters are
statistically more likely than high-$\beta$ 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.
[42]
oai:arXiv.org:1405.6034 [pdf] - 1209705
Boosting the accuracy of SPH techniques: Newtonian and
special-relativistic tests
Submitted: 2014-05-23, last modified: 2015-02-16
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.
[43]
oai:arXiv.org:1501.02081 [pdf] - 1223937
Multi-messenger picture of compact binary mergers
Submitted: 2015-01-09, last modified: 2015-02-11
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.
[44]
oai:arXiv.org:1412.5588 [pdf] - 1223695
The interplay of disk wind and dynamical ejecta in the aftermath of
neutron star - black hole mergers
Submitted: 2014-12-17, last modified: 2015-02-03
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.
[45]
oai:arXiv.org:1405.6730 [pdf] - 1209746
Neutrino-driven winds from neutron star merger remnants
Submitted: 2014-05-26
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.
[46]
oai:arXiv.org:1403.1297 [pdf] - 862557
MODA: a new algorithm to compute optical depths in multi-dimensional
hydrodynamic simulations
Submitted: 2014-03-05
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.
[47]
oai:arXiv.org:1402.1171 [pdf] - 781114
Generalized Newtonian description of particle motion in spherically
symmetric spacetimes
Submitted: 2014-02-05
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\"om, Ay\'on Beato-Garc\'ia, 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.
[48]
oai:arXiv.org:1307.2939 [pdf] - 779576
The long-term evolution of neutron star merger remnants - I. The impact
of r-process nucleosynthesis
Submitted: 2013-07-10, last modified: 2014-01-16
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.
[49]
oai:arXiv.org:1307.2943 [pdf] - 779577
The long-term evolution of neutron star merger remnants - II.
Radioactively powered transients
Submitted: 2013-07-10, last modified: 2014-01-16
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.
[50]
oai:arXiv.org:1401.2166 [pdf] - 771931
Implications of GRB 130603B and its macronova for r-process
nucleosynthesis
Submitted: 2014-01-09
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.
[51]
oai:arXiv.org:1311.6980 [pdf] - 1201897
Europium production: neutron star mergers versus core-collapse
supernovae
Submitted: 2013-11-27
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.
[52]
oai:arXiv.org:1308.1667 [pdf] - 1173261
The structure and fate of white dwarf merger remnants
Submitted: 2013-08-07
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.
[53]
oai:arXiv.org:1306.5259 [pdf] - 683968
Light from the Cosmic Frontier: Gamma-Ray Bursts
Amati, L.;
Atteia, J. -L.;
Balazs, L.;
Basa, S.;
Tjus, J. Becker;
Bersier, D. F.;
Boer, M.;
Campana, S.;
Ciardi, B.;
Covino, S.;
Daigne, F.;
Feroci, M.;
Ferrara, A.;
Frontera, F.;
Fynbo, J. P. U.;
Ghirlanda, G.;
Ghisellini, G.;
Glover, S.;
Greiner, J.;
Gotz, D.;
Hanlon, L.;
Hjorth, J.;
Hudec, R.;
Katz, U.;
Khochfar, S.;
Klessen, R.;
Kowalski, M.;
Levan, A. J.;
McBreen, S.;
Mesinger, A.;
Mochkovitch, R.;
O'Brien, P.;
Osborne, J. P.;
Petitjean, P.;
Reimer, O.;
Resconi, E.;
Rosswog, S.;
Ryde, F.;
Salvaterra, R.;
Savaglio, S.;
Schneider, R.;
Tagliaferri, G.;
Tanvir, N. R.;
van der Horst, A.
Submitted: 2013-06-21
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.
[54]
oai:arXiv.org:1303.4068 [pdf] - 680253
An accurate Newtonian description of particle motion around a
Schwarzschild black hole
Submitted: 2013-03-17, last modified: 2013-06-14
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.
[55]
oai:arXiv.org:1305.5720 [pdf] - 667680
The Gravitational Universe
Consortium, The eLISA;
:;
Seoane, P. Amaro;
Aoudia, S.;
Audley, H.;
Auger, G.;
Babak, S.;
Baker, J.;
Barausse, E.;
Barke, S.;
Bassan, M.;
Beckmann, V.;
Benacquista, M.;
Bender, P. L.;
Berti, E.;
Binétruy, P.;
Bogenstahl, J.;
Bonvin, C.;
Bortoluzzi, D.;
Brause, N. C.;
Brossard, J.;
Buchman, S.;
Bykov, I.;
Camp, J.;
Caprini, C.;
Cavalleri, A.;
Cerdonio, M.;
Ciani, G.;
Colpi, M.;
Congedo, G.;
Conklin, J.;
Cornish, N.;
Danzmann, K.;
de Vine, G.;
DeBra, D.;
Freitag, M. Dewi;
Di Fiore, L.;
Aguilo, M. Diaz;
Diepholz, I.;
Dolesi, R.;
Dotti, M.;
Barranco, G. Fernández;
Ferraioli, L.;
Ferroni, V.;
Finetti, N.;
Fitzsimons, E.;
Gair, J.;
Galeazzi, F.;
Garcia, A.;
Gerberding, O.;
Gesa, L.;
Giardini, D.;
Gibert, F.;
Grimani, C.;
Groot, P.;
Cervantes, F. Guzman;
Haiman, Z.;
Halloin, H.;
Heinzel, G.;
Hewitson, M.;
Hogan, C.;
Holz, D.;
Hornstrup, A.;
Hoyland, D.;
Hoyle, C. D.;
Hueller, M.;
Hughes, S.;
Jetzer, P.;
Kalogera, V.;
Karnesis, N.;
Kilic, M.;
Killow, C.;
Klipstein, W.;
Kochkina, E.;
Korsakova, N.;
Krolak, A.;
Larson, S.;
Lieser, M.;
Littenberg, T.;
Livas, J.;
Lloro, I.;
Mance, D.;
Madau, P.;
Maghami, P.;
Mahrdt, C.;
Marsh, T.;
Mateos, I.;
Mayer, L.;
McClelland, D.;
McKenzie, K.;
McWilliams, S.;
Merkowitz, S.;
Miller, C.;
Mitryk, S.;
Moerschell, J.;
Mohanty, S.;
Monsky, A.;
Mueller, G.;
Müller, V.;
Nelemans, G.;
Nicolodi, D.;
Nissanke, S.;
Nofrarias, M.;
Numata, K.;
Ohme, F.;
Otto, M.;
Perreur-Lloyd, M.;
Petiteau, A.;
Phinney, E. S.;
Plagnol, E.;
Pollack, S.;
Porter, E.;
Prat, P.;
Preston, A.;
Prince, T.;
Reiche, J.;
Richstone, D.;
Robertson, D.;
Rossi, E. M.;
Rosswog, S.;
Rubbo, L.;
Ruiter, A.;
Sanjuan, J.;
Sathyaprakash, B. S.;
Schlamminger, S.;
Schutz, B.;
Schütze, D.;
Sesana, A.;
Shaddock, D.;
Shah, S.;
Sheard, B.;
Sopuerta, C. F.;
Spector, A.;
Spero, R.;
Stanga, R.;
Stebbins, R.;
Stede, G.;
Steier, F.;
Sumner, T.;
Sun, K. -X.;
Sutton, A.;
Tanaka, T.;
Tanner, D.;
Thorpe, I.;
Tröbs, M.;
Tinto, M.;
Tu, H. -B.;
Vallisneri, M.;
Vetrugno, D.;
Vitale, S.;
Volonteri, M.;
Wand, V.;
Wang, Y.;
Wanner, G.;
Ward, H.;
Ware, B.;
Wass, P.;
Weber, W. J.;
Yu, Y.;
Yunes, N.;
Zweifel, P.
Submitted: 2013-05-24
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.
[56]
oai:arXiv.org:1210.1214 [pdf] - 657312
The runaway instability in general relativistic accretion disks
Submitted: 2012-10-03, last modified: 2013-04-27
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.
[57]
oai:arXiv.org:1204.6240 [pdf] - 1118331
The multi-messenger picture of compact object encounters: binary mergers
versus dynamical collisions
Submitted: 2012-04-27, last modified: 2013-01-09
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.
[58]
oai:arXiv.org:1204.6242 [pdf] - 1118332
The Electromagnetic Signals of Compact Binary Mergers
Submitted: 2012-04-27, last modified: 2013-01-07
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.
[59]
oai:arXiv.org:1210.6549 [pdf] - 1152407
The dynamic ejecta of compact object mergers and eccentric collisions
Submitted: 2012-10-24, last modified: 2012-11-23
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.
[60]
oai:arXiv.org:1207.0512 [pdf] - 587233
The Viscous Evolution of White Dwarf Merger Remnants
Submitted: 2012-07-02, last modified: 2012-11-06
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
$\alpha$-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.
[61]
oai:arXiv.org:1206.2379 [pdf] - 1124046
On the astrophysical robustness of neutron star merger r-process
Submitted: 2012-06-11, last modified: 2012-08-04
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.
[62]
oai:arXiv.org:1206.0744 [pdf] - 1123879
On the point mass approximation to calculate the gravitational wave
signal from white dwarf binaries
Submitted: 2012-06-04
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.
[63]
oai:arXiv.org:1201.2406 [pdf] - 1092885
How the merger of two white dwarfs depends on their mass ratio: orbital
stability and detonations at contact
Submitted: 2012-01-11
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.
[64]
oai:arXiv.org:1108.0028 [pdf] - 442271
A fast recursive coordinate bisection tree for neighbour search and
gravity
Submitted: 2011-07-29
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.
[65]
oai:arXiv.org:1101.5132 [pdf] - 1051705
Prelude to a double degenerate merger: the onset of mass transfer and
its impact on gravitational waves and surface detonations
Submitted: 2011-01-26
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.
[66]
oai:arXiv.org:1012.0912 [pdf] - 955682
Compact binary mergers: an astrophysical perspective
Submitted: 2010-12-04, last modified: 2010-12-22
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.
[67]
oai:arXiv.org:1009.1385 [pdf] - 1034776
Precursors of short gamma-ray bursts
Submitted: 2010-09-07
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.
[68]
oai:arXiv.org:0907.4890 [pdf] - 1003140
Conservative, special-relativistic smoothed particle hydrodynamics
Submitted: 2009-07-28, last modified: 2010-08-19
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 $\gamma= 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.
[69]
oai:arXiv.org:1005.1679 [pdf] - 163302
Special-relativistic Smoothed Particle Hydrodynamics: a benchmark suite
Submitted: 2010-05-10
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.
[70]
oai:arXiv.org:0911.0416 [pdf] - 1018187
Surface Detonations in Double Degenerate Binary Systems Triggered by
Accretion Stream Instabilities
Submitted: 2009-11-03
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-$\alpha$ 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.
[71]
oai:arXiv.org:0907.3196 [pdf] - 1002959
Collisions of white dwarfs as a new progenitor channel for type Ia
supernovae
Submitted: 2009-07-18, last modified: 2009-10-01
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.
[72]
oai:arXiv.org:0909.3632 [pdf] - 1017608
Going out with a bang: compact object collisions resulting from
supernovae in binary systems
Submitted: 2009-09-21
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.
[73]
oai:arXiv.org:0811.1370 [pdf] - 1001060
Three-Dimensional Simulations of Tidally Disrupted Solar-Type Stars and
the Observational Signatures of Shock Breakout
Submitted: 2008-11-10, last modified: 2009-09-17
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.
[74]
oai:arXiv.org:0903.5075 [pdf] - 22864
Astrophysical Smooth Particle Hydrodynamics
Submitted: 2009-03-29, last modified: 2009-08-05
The paper presents a detailed review of the smooth particle hydrodynamics
(SPH) method with particular focus on its astrophysical applications. We start
by introducing the basic ideas and concepts and thereby outline all ingredients
that are necessary for a practical implementation of the method in a working
SPH code. Much of SPH's success relies on its excellent conservation properties
and therefore the numerical conservation of physical invariants receives much
attention throughout this review. The self-consistent derivation of the SPH
equations from the Lagrangian of an ideal fluid is the common theme of the
remainder of the text. We derive a modern, Newtonian SPH formulation from the
Lagrangian of an ideal fluid. It accounts for changes of the local resolution
lengths which result in corrective, so-called "grad-h-terms". We extend this
strategy to special relativity for which we derive the corresponding grad-h
equation set. The variational approach is further applied to a
general-relativistic fluid evolving in a fixed, curved background space-time.
Particular care is taken to explicitely derive all relevant equations in a
coherent way.
[75]
oai:arXiv.org:0808.3847 [pdf] - 315041
The star ingesting luminosity of intermediate mass black holes in
globular clusters
Submitted: 2008-08-28, last modified: 2009-04-10
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)
[76]
oai:arXiv.org:0808.2143 [pdf] - 900366
Tidal disruption and ignition of white dwarfs by moderately massive
black holes
Submitted: 2008-08-15, last modified: 2008-11-17
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$).
[77]
oai:arXiv.org:0811.2129 [pdf] - 18504
Tidally-induced thermonuclear Supernovae
Submitted: 2008-11-13
We discuss the results of 3D simulations of tidal disruptions of white dwarfs
by moderate-mass black holes as they may exist in the cores of globular
clusters or dwarf galaxies. Our simulations follow self-consistently the
hydrodynamic and nuclear evolution from the initial parabolic orbit over the
disruption to the build-up of an accretion disk around the black hole. For
strong enough encounters (pericentre distances smaller than about 1/3 of the
tidal radius) the tidal compression is reversed by a shock and finally results
in a thermonuclear explosion. These explosions are not restricted to progenitor
masses close to the Chandrasekhar limit, we find exploding examples throughout
the whole white dwarf mass range. There is, however, a restriction on the
masses of the involved black holes: black holes more massive than $2\times
10^5$ M$_\odot$ swallow a typical 0.6 M$_\odot$ dwarf before their tidal forces
can overwhelm the star's self-gravity. Therefore, this mechanism is
characteristic for black holes of moderate masses. The material that remains
bound to the black hole settles into an accretion disk and produces an X-ray
flare close to the Eddington limit of $L_{\rm Edd} \simeq 10^{41} {\rm erg/s}
M_{\rm bh}/1000 M$_\odot$), typically lasting for a few months. The combination
of a peculiar thermonuclear supernova together with an X-ray flare thus
whistle-blows the existence of such moderate-mass black holes. The next
generation of wide field space-based instruments should be able to detect such
events.
[78]
oai:arXiv.org:0811.1517 [pdf] - 18370
Mass transfer dynamics in double degenerate binary systems
Submitted: 2008-11-10
We present a numerical study of the mass transfer dynamics prior to the
gravitational wave-driven merger of a double white dwarf system. Recently,
there has been some discussion about the dynamics of these last stages,
different methods seemed to provide qualitatively different results. While
earlier SPH simulations indicated a very quick disruption of the binary on
roughly the orbital time scale, more recent grid-based calculations find
long-lived mass transfer for many orbital periods. Here we demonstrate how
sensitive the dynamics of this last stage is to the exact initial conditions.
We show that, after a careful preparation of the initial conditions, the
reportedly short-lived systems undergo mass transfer for many dozens of orbits.
The reported numbers of orbits are resolution-biased and therefore represent
only lower limits to what is realized in nature. Nevertheless, the study shows
convincingly the convergence of different methods to very similar results.
[79]
oai:arXiv.org:0806.4380 [pdf] - 13951
Neutrino signatures and the neutrino-driven wind in Binary Neutron Star
Mergers
Submitted: 2008-06-26
We present VULCAN/2D multi-group flux-limited-diffusion radiation
hydrodynamics simulations of binary neutron star (BNS) mergers, using the Shen
equation of state, covering ~100 ms, and starting from azimuthal-averaged 2D
slices obtained from 3D SPH simulations of Rosswog & Price for 1.4 Msun
(baryonic) neutron stars with no initial spins, co-rotating spins, and
counter-rotating spins. Snapshots are post-processed at 10 ms intervals with a
multi-angle neutrino-transport solver. We find polar-enhanced neutrino
luminosities, dominated by $\bar{\nu}_e$ and ``$\nu_\mu$'' neutrinos at peak,
although $\nu_e$ emission may be stronger at late times. We obtain typical peak
neutrino energies for $\nu_e$, $\bar{\nu}_e$, and ``$\nu_\mu$'' of ~12, ~16,
and ~22 MeV. The super-massive neutron star (SMNS) formed from the merger has a
cooling timescale of ~1 s. Charge-current neutrino reactions lead to the
formation of a thermally-driven bipolar wind with <$\dot{M}$> ~10$^{-3}$
Msun/s, baryon-loading the polar regions, and preventing any production of a
GRB prior to black-hole formation. The large budget of rotational free energy
suggests magneto-rotational effects could produce a much greater polar mass
loss. We estimate that ~10$^{-4}$ Msun of material with electron fraction in
the range 0.1-0.2 become unbound during this SMNS phase as a result of neutrino
heating. We present a new formalism to compute the $\nu_i\bar{\nu}_i$
annihilation rate based on moments of the neutrino specific intensity computed
with our multi-angle solver. Cumulative annihilation rates, which decay as
$t^{-1.8}$, decrease over our 100 ms window from a few 10$^{50}$ to ~10$^{49}$
erg/s, equivalent to a few 10$^{54}$ to ~10$^{53}$ $e^-e^+$ pairs per second.
[80]
oai:arXiv.org:0802.0418 [pdf] - 9729
3D meshfree magnetohydrodynamics
Submitted: 2008-02-04
We describe a new method to include magnetic fields into smooth particle
hydrodynamics. The derivation of the self-gravitating hydrodynamics equations
from a variational principle is discussed in some detail. The non-dissipative
magnetic field evolution is instantiated by advecting so-called Euler
potentials. This approach enforces the crucial
$\nabla\cdot\vec{B}=0$-constraint by construction. These recent developments
are implemented in our three-dimensional, self-gravitating magnetohydrodynamics
code MAGMA. A suite of tests is presented that demonstrates the superiority of
this new approach in comparison to previous implementations.
[81]
oai:arXiv.org:0801.1582 [pdf] - 8865
Simulating Black Hole White Dwarf Encounters
Submitted: 2008-01-10, last modified: 2008-01-17
The existence of supermassive black holes lurking in the centers of galaxies
and of stellar binary systems containing a black hole with a few solar masses
has been established beyond reasonable doubt. The idea that black holes of
intermediate masses ($\sim 1000$ \msun) may exist in globular star clusters has
gained credence over recent years but no conclusive evidence has been
established yet. An attractive feature of this hypothesis is the potential to
not only disrupt solar-type stars but also compact white dwarf stars. In close
encounters the white dwarfs can be sufficiently compressed to thermonuclearly
explode. The detection of an underluminous thermonuclear explosion accompanied
by a soft, transient X-ray signal would be compelling evidence for the presence
of intermediate mass black holes in stellar clusters. In this paper we focus on
the numerical techniques used to simulate the entire disruption process from
the initial parabolic orbit, over the nuclear energy release during tidal
compression, the subsequent ejection of freshly synthesized material and the
formation process of an accretion disk around the black hole.
[82]
oai:arXiv.org:0712.2513 [pdf] - 8167
Atypical Thermonuclear Supernovae from Tidally Crushed White Dwarfs
Submitted: 2007-12-15
Suggestive evidence has accumulated that intermediate mass black holes (IMBH)
exist in some globular clusters. As stars diffuse in the cluster, some will
inevitable wander sufficiently close to the hole that they suffer tidal
disruption. An attractive feature of the IMBH hypothesis is its potential to
disrupt not only solar-type stars but also compact white dwarf stars. Attention
is given to the fate of white dwarfs that approach the hole close enough to be
disrupted and compressed to such extent that explosive nuclear burning may be
triggered. Precise modeling of the dynamics of the encounter coupled with a
nuclear network allow for a realistic determination of the explosive energy
release, and it is argued that ignition is a natural outcome for white dwarfs
of all varieties passing well within the tidal radius. Although event rates are
estimated to be significantly less than the rate of normal Type Ia supernovae,
such encounters may be frequent enough in globular clusters harboring an IMBH
to warrant a search for this new class of supernova.
[83]
oai:arXiv.org:0705.1441 [pdf] - 1145
MAGMA: a 3D, Lagrangian magnetohydrodynamics code for merger
applications
Submitted: 2007-05-10
We present a new, completely Lagrangian magnetohydrodynamics code that is
based on the SPH method. The equations of self-gravitating hydrodynamics are
derived self-consistently from a Lagrangian and account for variable smoothing
length (``grad-h''-) terms in both the hydrodynamic and the gravitational
acceleration equations. The evolution of the magnetic field is formulated in
terms of so-called Euler potentials which are advected with the fluid and thus
guarantee the MHD flux-freezing condition. This formulation is equivalent to a
vector potential approach and therefore fulfills the
$\vec{\nabla}\cdot\vec{B}=0$-constraint by construction. Extensive tests in
one, two and three dimensions are presented. The tests demonstrate the
excellent conservation properties of the code and show the clear superiority of
the Euler potentials over earlier magnetic SPH formulations.
[84]
oai:arXiv.org:0704.0297 [pdf] - 1944028
Remnant evolution after a carbon-oxygen white dwarf merger
Submitted: 2007-04-02, last modified: 2007-04-04
We systematically explore the evolution of the merger of two carbon-oxygen
(CO) white dwarfs. The dynamical evolution of a 0.9 Msun + 0.6 Msun CO white
dwarf merger is followed by a three-dimensional SPH simulation. We use an
elaborate prescription in which artificial viscosity is essentially absent,
unless a shock is detected, and a much larger number of SPH particles than
earlier calculations. Based on this simulation, we suggest that the central
region of the merger remnant can, once it has reached quasi-static equilibrium,
be approximated as a differentially rotating CO star, which consists of a
slowly rotating cold core and a rapidly rotating hot envelope surrounded by a
centrifugally supported disc. We construct a model of the CO remnant that
mimics the results of the SPH simulation using a one-dimensional hydrodynamic
stellar evolution code and then follow its secular evolution. The stellar
evolution models indicate that the growth of the cold core is controlled by
neutrino cooling at the interface between the core and the hot envelope, and
that carbon ignition in the envelope can be avoided despite high effective
accretion rates. This result suggests that the assumption of forced accretion
of cold matter that was adopted in previous studies of the evolution of double
CO white dwarf merger remnants may not be appropriate. Our results imply that
at least some products of double CO white dwarfs merger may be considered good
candidates for the progenitors of Type Ia supernovae. In this case, the
characteristic time delay between the initial dynamical merger and the eventual
explosion would be ~10^5 yr. (Abridged).
[85]
oai:arXiv.org:astro-ph/0611440 [pdf] - 86826
Fallback accretion in the aftermath of a compact binary merger
Submitted: 2006-11-14, last modified: 2007-01-16
Recent observations of long and short gamma-ray bursts have revealed a
puzzling X-ray activity that in some cases continues for hours after the burst.
It is difficult to reconcile such time scales with the viscous time scales that
an accretion disk can plausibly provide. Here I discuss the accretion activity
expected from the material that is launched into eccentric, but gravitationally
bound orbits during a compact binary merger coalescence. From a simple
analytical model the time scales and accretion luminosities that result from
fallback in the aftermath of a compact binary merger are derived. For the
considered mass range, double neutron star binaries are relatively homogeneous
in their fallback luminosities. Neutron star black hole systems show a larger
spread in their fallback behaviour. While the model is too simple to make
predictions about the detailed time structure of the fallback, it makes
reasonable predictions about the gross properties of the fallback. About one
hour after the coalescence the fallback accretion luminosity can still be as
large as $\sim 10^{45}$ erg/s, a fraction of which will be transformed into
X-rays. Large-scale amplitude variations in the X-ray luminosities can
plausibly be caused by gravitational fragmentation, which for the
high-eccentricity fallback should occur more easily than in an accretion disk.
[86]
oai:arXiv.org:astro-ph/0612572 [pdf] - 87910
Last moments in the life of a compact binary system: gravitational
waves, gamma-ray bursts and magnetar formation
Submitted: 2006-12-20
The first detections of afterglows from short gamma-ray bursts (GRBs) have
confirmed the previous suspicion that they are triggered by a different central
engine than long bursts. In particular, the recent detections of short GRBs in
galaxies without star formation lends support to the idea that an old stellar
population is involved. Most prominent are mergers of either double neutron
stars or of a neutron star with a stellar-mass black hole companion. Since the
final identification of the central engine will only come from an integral view
of several properties, we review the observable signatures that can be expected
from both double neutron stars and neutron star black hole systems. We discuss
the gravitational wave emission, the structure of the neutrino-cooled accretion
disks, the resulting neutrino signal and possible mechanisms to launch a GRB.
In addition, we address the speculative idea that in some cases a magnetar-like
object may be the final outcome of a double neutron star merger. We also
discuss possibilities to explain the late-time X-ray activity that has been
observed in several bursts.
[87]
oai:arXiv.org:astro-ph/0603845 [pdf] - 81054
Producing ultra-strong magnetic fields in neutron star mergers
Submitted: 2006-03-31
We report an extremely rapid mechanism for magnetic field amplification
during the merger of a binary neutron star system. This has implications for
the production of the short class of Gamma-Ray Bursts, which recent
observations suggest may originate in such mergers. In detailed
magnetohydrodynamic simulations of the merger process, the fields are amplified
via Kelvin-Helmholtz instabilities beyond magnetar field strength and may
therefore represent the strongest magnetic fields in the Universe. The
amplification occurs in the shear layer which forms between the neutron stars
and on a time scale of only 1 millisecond, i.e. long before the remnant can
collapse into a black hole.
[88]
oai:arXiv.org:astro-ph/0508138 [pdf] - 75027
Mergers of neutron star black hole binaries with small mass ratios:
nucleosynthesis, gamma-ray bursts and electromagnetic transients
Submitted: 2005-08-04
I discuss simulations of the coalescence of black hole neutron star binary
systems with black hole masses between 14 and 20 \msun. The calculations use a
three-dimensional smoothed particle hydrodynamics code, a
temperature-dependent, nuclear equation of state and a multi-flavor neutrino
scheme. General relativistic effects are mimicked using the \Pacz-Wiita
pseudo-potential and gravitational radiation reaction forces. Opposite to
previous, purely Newtonian calculations, in none of the explored cases episodic
mass transfer occurs. The neutron star is always completely disrupted after
most of its mass has been transferred directly into the hole. For black hole
masses between 14 and 16 \Msun an accretion disk forms, large parts of it,
however, are inside the last stable orbit and therefore falling with large
radial velocities into the hole. These disks are (opposite to the neutron star
merger case) thin and -apart from a spiral shock- essentially cold. For higher
mass black holes ($M_{\rm BH} \ge 18$ \msun) almost the complete neutron star
disappears in the hole without forming an accretion disk. In these cases the
surviving material is spun up by tidal torques and ejected as a half-ring of
neutron-rich matter. None of the investigated systems is a promising GRB
central engine. We find between 0.01 and 0.2 \Msun of the neutron star to be
dynamically ejected. Like in a type Ia supernova, the radioactive decay of this
material will power a light curve with a peak luminosity of a few times
$10^{44}$ erg/s. The maximum will be reached about three days after the
coalescence and will be mainly visible in the optical/near infrared band. The
coalescence itself may produce a precursor pulse with a thermal spectrum of
$\sim 10$ ms duration.
[89]
oai:arXiv.org:astro-ph/0505007 [pdf] - 72777
On the viability of neutron star black hole binaries as central engines
of gamma-ray bursts
Submitted: 2005-04-30
I discuss three-dimensional SPH simulations of neutron star black hole
encounters. The calculations are performed using a nuclear equation of state
and a multi-flavor neutrino treatment, general relativistic effects are
mimicked using the Paczynski-Wiita pseudo-potential and gravitational radiation
reaction forces. Most of the explored mass range (14 to 20 \msun) has not been
considered before in numerical simulations. The neutron star is always
disrupted during the first approach after most of its mass has been transferred
directly into the hole. In none of the analyzed cases episodic mass transfer is
found. If an accretion disk forms at all, it is relatively cold ($T < 2.5$
MeV), of only moderate density ($\rho \sim 10^{10}$ g/ccm) and substantially
consumed by the hole on the simulation time scale of a few milliseconds. None
of the investigated systems of this study yields conditions that are promising
to launch a GRB. While we cannot completely exclude that a subset of neutron
star black hole binaries, maybe black holes with low masses and very large
initial spins, can produce a GRB, this seems to happen -if at all- only in a
restricted region of the available parameter space. I argue that the difficulty
to form promising disks together with the absence of any observed neutron star
black hole binary may mean that they are insignificant as central engines of
the observed, short-hard GRBs and that the vast majority of the latter ones is
caused by double neutron star coalescences.
[90]
oai:arXiv.org:astro-ph/0504368 [pdf] - 260587
From Neutron Star Binaries to Gamma-ray bursts
Submitted: 2005-04-17
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.
[91]
oai:arXiv.org:astro-ph/0403500 [pdf] - 63691
Accretion dynamics in neutron star black hole binaries
Submitted: 2004-03-20
We perform three-dimensional, Newtonian hydrodynamic simulations with a
nuclear equation of state to investigate the accretion dynamics in neutron star
black hole systems. We find as a general result that non-spinning donor stars
yield larger circularization radii than corotating donors. Therefore, the
matter from a neutron star without spin will more likely settle into an
accretion disk outside the Schwarzschild radius. With the used stiff equation
of state we find it hard to form an accretion disk that is promising to launch
a gamma-ray burst. In all relevant cases the core of the neutron star survives
and keeps orbiting the black hole as a mini neutron star for the rest of the
simulation time (up to several hundred dynamical neutron star times scales).
The existence of this mini neutron star leaves a clear imprint on the
gravitational wave signal which thus can be used to probe the physics at
supra-nuclear densities.
[92]
oai:arXiv.org:astro-ph/0401022 [pdf] - 61892
Gamma-ray Bursts and their Central Engines
Submitted: 2004-01-04, last modified: 2004-01-09
Gamma-ray bursts are the most luminous and probably the most relativistic
events in the universe. The last few years have seen a tremendous increase in
our knowledge of these events, but the source of the bursts still remains
elusive. I will summarise recent progress in this field with special emphasis
on our understanding of the possible progenitor systems.
[93]
oai:arXiv.org:astro-ph/0312191 [pdf] - 61404
Dynamos, Super-pulsars and Gamma-ray bursts
Submitted: 2003-12-07, last modified: 2003-12-10
The remnant of a neutron star binary coalescence is expected to be
temporarily stabilised against gravitational collapse by its differential
rotation. We explore the possibility of dynamo activity in this remnant and
assess the potential for powering a short-duration gamma-ray burst (GRB). We
analyse our three-dimensional hydrodynamic simulations of neutron star mergers
with respect to the flow pattern inside the remnant. If the central, newly
formed super-massive neutron star remains stable for a good fraction of a
second an efficient low-Rossby number $\alpha-\Omega$-dynamo will amplify the
initial seed magnetic fields exponentially. We expect that values close to
equipartition field strength will be reached within several tens of
milliseconds. Such a super-pulsar could power a GRB via a relativistic wind,
with an associated spin-down time scale close to the typical duration of a
short GRB. Similar mechanisms are expected to be operational in the surrounding
torus formed from neutron star debris.
[94]
oai:arXiv.org:astro-ph/0309201 [pdf] - 59028
On the Central Engine of short Gamma-ray Bursts
Submitted: 2003-09-06, last modified: 2003-09-16
We assess two mechanisms to produce a short GRB from a neutron star
coalescence: the annihilation of neutrino and anti-neutrino pairs and magnetic
processes within the remnant. We find that the neutrino-driven wind from the
remnant can collimate the jets produced via neutrino annihilation into a small
solid angle. This mechanism satisfies the observational constraints on the
apparent luminosity and on the merger event rate. Moreover, the supermassive
neutron star that is formed in the center of the remnant is expected to act as
a convective dynamo reaching equipartition field strengths within fractions of
a second. The associated spin-down time scale is similar to the typical
duration of short GRBs.
[95]
oai:arXiv.org:astro-ph/0306418 [pdf] - 1233164
High Resolution Calculations of Merging Neutron Stars III: Gamma-Ray
Bursts
Submitted: 2003-06-20, last modified: 2003-07-15
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 $\nu\bar{\nu}$ and magnetohydrodynamic processes-
are explored in order to investigate their viability in launching GRBs. We find
that $\nu\bar{\nu}$ 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
$\nu\bar{\nu}$-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.
[96]
oai:arXiv.org:astro-ph/0306172 [pdf] - 57253
On the diversity of short GRBs
Submitted: 2003-06-09
Hydrodynamical simulations of the last inspiral stages and the final
coalescence of a double neutron star system are used to investigate the power
of the neutrino-driven wind, the energy and momentum of the fireball produced
via $\nu \bar{\nu}$-annihilation, and the intensity and character of their
interaction. It is argued that the outflow that derives from the debris will
have enough pressure to collimate the relativistic fireball that it surrounds.
The low luminosity relativistic jet will then appear brighter to an observer
within the beam although most of the energy of the event is in the unseen, less
collimated and slower wind. This model leads to a simple physical
interpretation of the isotropic luminosities implied for short GRBs at
cosmological distances. A wide variety of burst phenomenology could be
attributable to the dependence of the neutrino luminosity on the initial mass
of the double NS binary.
[97]
oai:arXiv.org:astro-ph/0302301 [pdf] - 386165
High Resolution Calculations of Merging Neutron Stars II: Neutrino
Emission
Submitted: 2003-02-15
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''.
[98]
oai:arXiv.org:astro-ph/0207576 [pdf] - 50719
Jets, winds and bursts from coalescing neutron stars
Submitted: 2002-07-26
Recent high-resolution calculations of neutron star coalescences are used to
investigate whether nu bar{nu} annihilation can provide sufficient energy to
power gamma-ray bursts, especially those belonging to the short duration
category. Late time slices of the simulations, where the neutrino emission has
reached a maximum, stationary level are analyzed to address this question. We
find that nu bar{nu} annihilation can provide the necessary driving stresses
that lead to relativistic jet expansion. Maximum Lorentz factors of the order
of 15 are found along the binary rotation axis, but larger values are expected
to arise from higher numerical resolution. Yet the accompanying neutrino-driven
wind must be absent from the axis when the burst occurs or prohibitive baryon
loading may occur. We argue that even under the most favorable conditions, nu
bar{nu} annihilation is unlikely to power a burst in excess of ~ 10^{48} erg.
Unless the emission is beamed into less than one percent of the solid angle,
which we argue is improbable if it is collimated by gas pressure, this may fail
to satisfy the apparent isotropic energies inferred at cosmological distances.
This mechanism may nonetheless drive a precursor fireball, thereby evacuating a
cavity into which a later magnetic driven jet could expand. A large range of
time delays between the merger and the black hole formation are to be expected.
If the magnetic driven jet occurs after the black hole has formed, a time span
as long as weeks could pass between the neutrino powered precursor and the
magnetic driven GRB.
[99]
oai:arXiv.org:astro-ph/0204494 [pdf] - 49039
Neutron Star Binaries as Central Engines of GRBs
Submitted: 2002-04-29
We describe the results high resolution, hydrodynamic calculations of neutron
star mergers. The model makes use of a new, nuclear equation of state, accounts
for multi-flavour neutrino emission and solves the equations of hydrodynamics
using the smoothed particle hydrodynamics method with more than $10^6$
particles. The merger leaves behind a strongly differentially rotating central
object of $\sim 2.5$ M$_{\odot}$ together with a distribution of hot debris
material. For the most realistic case of initial neutron star spins, no sign of
a collapse to a black hole can be seen. We argue that the differential rotation
stabilizes the central object for $\sim 10^2$ s and leads to superstrong
magnetic fields. We find the neutrino emission from the hot debris around the
freshly-formed, supermassive neutron star to be substantially lower than
predicted previously. Therefore the annihilation of neutrino anti-neutrino
pairs will have difficulties to power very energetic bursts ($\gg 10^{49}$
erg).
[100]
oai:arXiv.org:astro-ph/0106539 [pdf] - 163474
An Adaptive Grid, Implicit Code for Spherically Symmetric, General
Relativistic Hydrodynamics in Comoving Coordinates
Submitted: 2001-06-28, last modified: 2002-04-22
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.
[101]
oai:arXiv.org:astro-ph/0204358 [pdf] - 48903
Gamma-ray bursts, supernova kicks, and gravitational radiation
Submitted: 2002-04-22
We suggest that the collapsing core of a massive rotating star may fragment
to produce two or more compact objects. Their coalescence under gravitational
radiation gives the resulting black hole or neutron star a significant kick
velocity, which may explain those observed in pulsars. A gamma-ray burst can
result only when this kick is small. Thus only a small fraction of
core-collapse supernovae produce gamma-ray bursts. The burst may be delayed
significantly (hours - days) after the supernova, as suggested by recent
observations. If this picture is correct, core-collapse supernovae should be
significant sources of gravitational radiation with a chirp signal similar to a
coalescing binary.
[102]
oai:arXiv.org:astro-ph/0202453 [pdf] - 47947
Stellar Sources of the Interstellar Medium
Thielemann, F. -K.;
Argast, D.;
Brachwitz, F.;
Martinez-Pinedo, G.;
Oechslin, R.;
Rauscher, T.;
Hix, W. R.;
Liebendoerfer, M.;
Mezzacappa, A.;
Hoeflich, P.;
Iwamoto, K.;
Nomoto, K.;
Schatz, H.;
Wiescher, M. C.;
Kratz, K. -L.;
Pfeiffer, B.;
Rosswog, S.
Submitted: 2002-02-25
With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li,
stars act as sources for the composition of the interstellar medium. Cosmic
rays are related to the latter and very probably due to acceleration of the
mixed interstellar medium by shock waves from supernova remnants. Thus, the
understanding of the abundance evolution in the interstellar medium and
especially the enrichment of heavy elements, as a function of space and time,
is essential. It reflects the history of star formation and the lifetimes of
the diverse contributing stellar objects. Therefore, the understanding of the
endpoints of stellar evolution is essential as well. These are mainly planetary
nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars,
but also events in binary systems can contribute, like e.g. supernovae of type
Ia, novae and possibly X-ray bursts and neutron star or neutron star - black
hole mergers. Despite many efforts, a full and self-consistent understanding of
supernovae (the main contributors to nucleosynthesis in galaxies) is not
existing, yet. Their fingerprints, however, seen either in spectra,
lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can
help to constrain the composition of their ejecta and related model
uncertainties.
[103]
oai:arXiv.org:astro-ph/0110180 [pdf] - 1354497
High Resolution Calculations of Merging Neutron Stars I: Model
Description and Hydrodynamic Evolution
Submitted: 2001-10-08, last modified: 2002-02-23
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.
[104]
oai:arXiv.org:gr-qc/0111005 [pdf] - 110871
Conformally Flat Smoothed Particle Hydrodynamics: Application to Neutron
Star Mergers
Submitted: 2001-11-02, last modified: 2002-02-07
We present a new 3D SPH code which solves the general relativistic field +
hydrodynamics equations in the conformally flat approximation. Several test
cases are considered to test different aspects of the code. We finally apply
then the code to the coalescence of a neutron star binary system. The neutron
stars are modeled by a polytropic equation of state (EoS) with adiabatic
indices $\Gamma=2.0$, $\Gamma=2.6$ and $\Gamma=3.0$. We calculate the
gravitational wave signals, luminosities and frequency spectra by employing the
quadrupole approximation for emission and back reaction in the slow motion
limit. In addition, we consider the amount of ejected mass.
[105]
oai:arXiv.org:astro-ph/0112168 [pdf] - 46541
The Supernova-GRB Connection
Submitted: 2001-12-06, last modified: 2002-01-24
We discuss the possible connection between supernova explosions (SN) and
gamma-ray bursters (GRB) from the perspective of our current understanding of
SN physics. Recent evidence strongly suggests that the explosion mechanism of
core collapse SN is intrinsically aspherical. Typically, a neutron star is
formed. However, the observed properties of the expanding SN envelopes remnants
make these objects very unlikely candidates for GRBs. Most candidates for a
GRB/SN connection seem to require the prompt or delayed formation of a black
hole. These include the collapse of very massive stars (e.g. hypernovae) and
'classical' SNe with a significant fallback of material over time scales of
hours to days, resulting in the collapse of the neutron star to a black hole.
We suggest the merger of a neutron star with a white dwarf as a subclass of
thermonuclear SNe and a potential candidate for a SN/GRB connection.
[106]
oai:arXiv.org:astro-ph/9910154 [pdf] - 387118
Post Newtonian SPH
Submitted: 1999-10-08, last modified: 2001-04-15
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.
[107]
oai:arXiv.org:astro-ph/0101476 [pdf] - 40614
Element Synthesis in Stars
Thielemann, F. -K.;
Brachwitz, F.;
Freiburghaus, C.;
Kolbe, E.;
Martinez-Pinedo, G.;
Rauscher, T.;
Rembges, F.;
Hix, W. R.;
Liebendoerfer, M.;
Mezzacappa, A.;
Kratz, K. -L.;
Pfeiffer, B.;
Langanke, K.;
Nomoto, K.;
Rosswog, S.;
Schatz, H.;
Wiescher, M.
Submitted: 2001-01-26, last modified: 2001-02-12
Except for 1H, 2H, 3He, 4He, and 7Li, originating from the Big Bang, all
heavier elements are made in stellar evolution and stellar explosions. Nuclear
physics, and in many cases nuclear structure far from stability, enters in a
crucial way. Therefore, we examine in this review the role of nuclear physics
in astrophysics in general and in particular how it affects stellar events and
the resulting nucleosynthesis. Stellar modeling addresses four major aspects:
1. energy generation and nucleosynthesis, 2. energy transport via conduction,
radiation or possibly convection, 3. hydrodynamics/hydrostatics, and finally 4.
thermodynamic properties of the matter involved. Nuclear Physics enters via
nuclear reaction cross sections and nuclear structure (affecting the
composition changes and nuclear energy generation), neutrino-nucleon and
neutrino-nucleus cross sections (affecting neutrino opacities and transport),
and e.g. the equation of state at and beyond nuclear densities which creates a
relation between the nuclear many-body problem and the hydrodynamic response
like pressure and entropy. In the following we review these four topics by
highlighting the role and impact of nuclear physics in each of these aspects of
stellar modeling. The main emphasis is put on the connection to element
synthesis.
[108]
oai:arXiv.org:astro-ph/0012046 [pdf] - 39635
Nucleosynthesis Calculations for the Ejecta of Neutron Star Coalescences
Submitted: 2000-12-02, last modified: 2000-12-09
We present the results of fully dynamical r-process network calculations for
the ejecta of neutron star mergers (NSMs). The late stages of the inspiral and
the final violent coalescence of a neutron star binary have been calculated in
detail using a 3D hydrodynamics code (Newtonian gravity plus backreaction
forces emerging from the emission of gravitational waves) and a realistic
nuclear equation of state. The found trajectories for the ejecta serve as input
for dynamical r-process calculations where all relevant nuclear reactions
(including beta-decays depositing nuclear energy in the expanding material) are
followed. We find that all the ejected material undergoes r-process. For an
initial Ye close to 0.1 the abundance distributions reproduce very accurately
the solar r-process pattern for nuclei with A above 130. For lighter nuclei
strongly underabundant (as compared to solar) distributions are encountered. We
show that this behaviour is consistent with the latest observations of very
old, metal-poor stars, despite simplistic arguments that have recently been
raised against the possibility of NSM as possible sources of Galactic r-process
material.
[109]
oai:arXiv.org:astro-ph/0005550 [pdf] - 36321
Merging neutron stars: asymmetric systems
Submitted: 2000-05-29
Results of 3D, Newtonian hydrodynamic calculations of the last stages of the
inspiral and the final coalescence of neutron star binary systems are
described. The focus is on sligthly asymmetric systems with the asymmetry
stemming either from different spins or different masses of the binary
components. Implications for nucleosynthesis and gamma ray bursts are
discussed. Test calculations for the used viscosity scheme are provided in the
appendix.
[110]
oai:arXiv.org:astro-ph/9811367 [pdf] - 401144
Mass ejection in neutron star mergers
Submitted: 1998-11-23
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.
[111]
oai:arXiv.org:astro-ph/9804332 [pdf] - 1469635
Coalescing Neutron Stars: A Solution to the R-Process Problem?
Submitted: 1998-04-30, last modified: 1998-05-11
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.
