We describe an unambiguous gravitational-wave signature to identify the occurrence of a strong phase transition from hadronic matter to deconfined quark matter in neutron star mergers. Such a phase transition leads to a strong softening of the equation of state and hence to more compact merger remnants compared to purely hadronic models. If a phase transition takes place during merging, this results in a characteristic increase of the dominant postmerger gravitational-wave frequency relative to the tidal deformability characterizing the inspiral phase. By comparing results from different purely hadronic and hybrid models we show that a strong phase transition can be identified from a single, simultaneous measurement of pre- and postmerger gravitational waves. Furthermore, we present new results for hybrid star mergers, which contain quark matter already during the inspiral stage. Also for these systems we find that the postmerger GW frequency is increased compared to purely hadronic models. We thus conclude that also hybrid star mergers with an onset of the hadron-quark phase transition at relatively low densities may lead to the very same characteristic signature of quark deconfinement in the postmerger GW signal as systems undergoing the phase transition during merging.

We study the possible occurrence of the hadron-quark phase transition (PT) during the merging of neutron star binaries by hydrodynamical simulations employing a set of temperature dependent hybrid equations of state (EoSs). Following previous work we describe an unambiguous and measurable signature of deconfined quark matter in the gravitational-wave (GW) signal of neutron star binary mergers including equal-mass and unequal-mass systems of different total binary mass. The softening of the EoS by the PT at higher densities, i.e. after merging, leads to a characteristic increase of the dominant postmerger GW frequency f_peak relative to the tidal deformability Lambda inferred during the premerger inspiral phase. Hence, measuring such an increase of the postmerger frequency provides evidence for the presence of a strong PT. If the postmerger frequency and the tidal deformability are compatible with results from purely baryonic EoS models yielding very tight relations between f_peak and Lambda, a strong PT can be excluded up to a certain density. We find tight correlations of f_peak and Lambda with the maximum density during the early postmerger remnant evolution. These GW observables thus inform about the density regime which is probed by the remnant and its GW emission. Exploiting such relations we devise a directly applicable, concrete procedure to constrain the onset density of the QCD PT from future GW measurements. We point out two interesting scenarios: if no indications for a PT are inferred from a GW detection, our procedure yields a lower limit on the onset density of the hadron quark PT. On the contrary, if a merger event reveals evidence for the occurrence of deconfined quark matter, the inferred GW parameters set an upper limit on the PT onset density. (abridged)

Using hydrodynamical simulations for a large set of high-density matter equations of state (EoSs) we systematically determine the threshold mass M_thres for prompt black-hole formation in equal-mass and asymmetric neutron star (NS) mergers. We devise the so far most direct, general and accurate method to determine the unknown maximum mass of nonrotating NSs from merger observations revealing M_thres. Considering hybrid EoSs with hadron-quark phase transition, we identify a new, observable signature of quark matter in NS mergers. Furthermore, our findings have direct applications in gravitational wave searches, kilonova interpretations and multi-messenger constraints on NS properties.

Half of all the elements in the universe heavier than iron were created by rapid neutron capture. The theory for this astrophysical `$r$-process' was worked out six decades ago and requires an enormous neutron flux to make the bulk of these elements. Where this happens is still debated. A key piece of missing evidence is the identification of freshly-synthesised $r$-process elements in an astrophysical site. Current models and circumstantial evidence point to neutron star mergers as a probable $r$-process site, with the optical/infrared `kilonova' emerging in the days after the merger a likely place to detect the spectral signatures of newly-created neutron-capture elements. The kilonova, AT2017gfo, emerging from the gravitational-wave--discovered neutron star merger, GW170817, was the first kilonova where detailed spectra were recorded. When these spectra were first reported it was argued that they were broadly consonant with an outflow of radioactive heavy elements, however, there was no robust identification of any element. Here we report the identification of the neutron-capture element strontium in a re-analysis of these spectra. The detection of a neutron-capture element associated with the collision of two extreme-density stars establishes the origin of $r$-process elements in neutron star mergers, and demonstrates that neutron stars contain neutron-rich matter.

We construct new, multivariate empirical relations for measuring neutron star radii and tidal deformabilities from the dominant gravitational wave frequency in the post-merger phase of binary neutron star mergers. The relations determine neutron star radii and tidal deformabilities for specific neutron star masses with consistent accuracy and depend only on two observables: the post-merger peak frequency $f_{\rm peak}$ and the chirp mass $M_{\rm chirp}$. The former could be measured with good accuracy from gravitational waves emitted in the post-merger phase using next-generation detectors, whereas the latter is already obtained with good accuracy from the inspiral phase with present-day detectors. Our main data set consists of a gravitational wave catalogue obtained with CFC/SPH simulations. We also extract the $f_{\rm peak}$ frequency from the publicly available CoRe data set, obtained through grid-based GRHD simulations and find good agreement between the extracted frequencies of the two data sets. As a result, we can construct empirical relations for the combined data sets. Furthermore, we investigate empirical relations for two secondary peaks, $f_{2-0}$ and $f_{\rm spiral}$, and show that these relations are distinct in the whole parameter space, in agreement with a previously introduced spectral classification scheme. Finally, we show that the spectral classification scheme can be reproduced using machine-learning techniques.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

We describe a multi-messenger interpretation of GW170817, which yields a robust lower limit on NS radii. This excludes NSs with radii smaller than about 10.7 km and thus rules out very soft nuclear matter. We stress the potential of this type of constraints when future detections become available. A very similar argumentation may yield an upper bound on the maximum mass of nonrotating NSs. We also discuss simulations of NS mergers, which undergo a first-order phase transition to quark matter. We point out a different dynamical behavior. Considering the gravitational-wave signal, we identify an unambiguous signature of the QCD phase transition in NS mergers. The occurrence of quark matter through a strong first-order phase transition during merging leads to a characteristic shift of the dominant postmerger frequency. The frequency shift is indicative for a phase transition if it is compared to the postmerger frequency which is expected for purely hadronic EoS models. A very strong deviation of several 100 Hz is observed for hybrid EoSs in an otherwise tight relation between the tidal deformability and the postmerger frequency. We address the potential impact of a first-order phase transition on the electromagnetic counterpart of NS mergers. Our simulations suggest that there would be no significant qualitative differences between a system undergoing a phase transition to quark matter and purely hadronic mergers. The quantitative differences are within the spread which is found between different hadronic EoS models. This implies on the one hand that GW170817 is compatible with a possible transition to quark matter. On the other hand these considerations show that it may not be easy to identify quantitative differences between purely hadronic mergers and events in which quark matter occurs considering solely their electromagnetic counterpart or their nucleosynthesis products. (abridged)

We present a new, computationally efficient, energy-integrated approximation for neutrino effects in hot and dense astrophysical environments such as supernova cores and compact binary mergers and their remnants. Our new method, termed ILEAS for Improved Leakage-Equilibration-Absorption Scheme, improves the lepton-number and energy losses of traditional leakage descriptions by a novel prescription of the diffusion time-scale based on a detailed energy integral of the flux-limited diffusion equation. The leakage module is supplemented by a neutrino-equilibration treatment that ensures the proper evolution of the total lepton number and medium plus neutrino energies as well as neutrino-pressure effects in the neutrino-trapping domain. Moreover, we employ a simple and straightforwardly applicable ray-tracing algorithm for including re-absorption of escaping neutrinos especially in the decoupling layer and during the transition to semi-transparent conditions. ILEAS is implemented on a three-dimensional (3D) Cartesian grid with a minimum of free and potentially case-dependent parameters and exploits the basic physics constraints that should be fulfilled in the neutrino-opaque and free-streaming limits. We discuss a suite of tests for stationary and time-dependent proto-neutron star models and post-merger black-hole-torus configurations, for which 3D ILEAS results are demonstrated to agree with energy-dependent 1D and 2D two-moment (M1) neutrino transport on the level of 10--15 percent in basic neutrino properties. This also holds for the radial profiles of the neutrino luminosities and of the electron fraction. Even neutrino absorption maps around torus-like neutrino sources are qualitatively similar without any fine-tuning, confirming that ILEAS can satisfactorily reproduce local losses and re-absorption of neutrinos as found in sophisticated transport calculations.

The recent detection of a neutron star binary through gravitational waves, GW170817, has offered another source of information about the properties of cold supranuclear matter. Information from the signal emitted before the neutron stars merged has been used to study the equation of state of these bodies, however, any complementary information included in the signal emitted after the merger has been lost in the detector noise. In this paper we investigate the prospects of studying GW170817-like post-merger signals with future gravitational-wave detectors. We first compute the expected properties of the possible GW170817 post-merger signal using information from pre-merger analyses. We then quantify the required improvement in detector sensitivity in order to extract key features of the post-merger signal. We find that if we observe a signal of similar strength to GW170817 when the aLIGO detectors have been improved by $\sim 2-3$ times over their design sensitivity in the kHz regime, we will be able to extract the dominant frequency component of the post-merger. With further improvements and next-generation detectors we will also be able to extract subdominant frequencies. We conclude that post-merger signals could be brought within our reach in the coming years given planned detector upgrades, such as A+, Voyager, and the next-generation detectors.

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.

We identify an observable imprint of a first-order hadron-quark phase transition at supranuclear densities on the gravitational-wave (GW) emission of neutron star mergers. Specifically, we show that the dominant postmerger GW frequency f_peak may exhibit a significant deviation from an empirical relation between f_peak and the tidal deformability if a strong first-order phase transition leads to the formation of a gravitationally stable extended quark matter core in the postmerger remnant. A comparison of the GW signatures from a large, representative sample of microphysical, purely hadronic equations of state indicates that this imprint is only observed in those systems which undergo a strong first-order phase transition. Such a shift of the dominant postmerger GW frequency can be revealed by future GW observations, which would provide evidence for the existence of a strong first-order phase transition in the interior of neutron stars.

The physics of neutron stars can be studied with gravitational waves emitted from coalescing binary systems. Tidal effects become significant during the last few orbits and can be visible in the gravitational-wave spectrum above 500 Hz. After the merger, the neutron star remnant oscillates at frequencies above 1 kHz and can collapse into a black hole. Gravitational-wave detectors with a sensitivity of ~10^{-24} strain/sqHz at 2-4 kHz can observe these oscillations from a source which is ~100 Mpc away. The current observatories, such as LIGO and Virgo, are limited by shot noise at high frequencies and have a sensitivity of > 2 * 10^{-23} strain/sqHz at 3 kHz. In this paper, we propose an optical configuration of gravitational-wave detectors which can be set up in present facilities using the current interferometer topology. This scheme has a potential to reach 7 * 10^{-25} strain/sqHz at 2.5 kHz without compromising the detector sensitivity to black hole binaries. We argue that the proposed instruments have a potential to detect similar amount of post-merger neutron star oscillations as the next generation detectors, such as Cosmic Explorer and Einstein Telescope. We also optimise the arm length of the future detectors for neutron star physics and find that the optimal arm length is ~20 km. These instruments have the potential to observe neutron star post-merger oscillations at a rate of ~30 events per year with a signal-to-noise ratio of 5 or more.

This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.

We present a robust method to characterize the gravitational wave emission from the remnant of a neutron star coalescence. Our approach makes only minimal assumptions about the morphology of the signal and provides a full posterior probability distribution of the underlying waveform. We apply our method on simulated data from a network of advanced ground-based detectors and demonstrate the gravitational wave signal reconstruction. We study the reconstruction quality for different binary configurations and equations of state for the colliding neutron stars. We show how our method can be used to constrain the yet-uncertain equation of state of neutron star matter. The constraints on the equation of state we derive are complimentary to measurements of the tidal deformation of the colliding neutron stars during the late inspiral phase. In the case of a non-detection of a post-merger signal following a binary neutron star inspiral we show that we can place upper limits on the energy emitted.

We introduce a new, powerful method to constrain properties of neutron stars (NSs). We show that the total mass of GW170817 provides a reliable constraint on the stellar radius if the merger did not result in a prompt collapse as suggested by the interpretation of associated electromagnetic emission. The radius R_1.6 of nonrotating NSs with a mass of 1.6 M_sun can be constrained to be larger than 10.68_{-0.04}^{+0.15} km, and the radius R_max of the nonrotating maximum mass configuration must be larger than 9.60_{-0.03}^{+0.14} km. We point out that detections of future events will further improve these constraints. Moreover, we show that a future event with a signature of a prompt collapse of the merger remnant will establish even stronger constraints on the NS radius from above and the maximum mass M_max of NSs from above. These constraints are particularly robust because they only require a measurement of the chirp mass and a distinction between prompt and delayed collapse of the merger remnant, which may be inferred from the electromagnetic signal or even from the presence/absence of a ringdown gravitational-wave (GW) signal. This prospect strengthens the case of our novel method of constraining NS properties, which is directly applicable to future GW events with accompanying electromagnetic counterpart observations. We emphasize that this procedure is a new way of constraining NS radii from GW detections independent of existing efforts to infer radius information from the late inspiral phase or postmerger oscillations, and it does not require particularly loud GW events.

A binary neutron star merger produces a rapidly and differentially rotating compact remnant whose lifespan heavily affects the electromagnetic and gravitational emissions. Its stability depends on both the equation of state (EOS) and the rotation law and it is usually investigated through numerical simulations. Nevertheless, by means of a sufficient criterion for secular instability, equilibrium sequences can be used as a computationally inexpensive way to estimate the onset of dynamical instability, which, in general, is close to the secular one. This method works well for uniform rotation and relies on the location of turning points: stellar models that are stationary points in a sequence of equilibrium solutions with constant rest mass or angular momentum. Here, we investigate differentially rotating models (using a large number of equations of state and different rotation laws) and find that several universal relations between properly scaled gravitational mass, rest mass and angular momentum of the turning-point models that are valid for uniform rotation, are insensitive to the degree of differential rotation, to high accuracy.

The threshold mass for prompt collapse in binary neutron star mergers was empirically found to depend on the stellar properties of the maximum-mass non-rotating neutron star model. Here we present a semi-analytic derivation of this empirical relation which suggests that it is rather insensitive to thermal effects, to deviations from axisymmetry and to the exact rotation law in merger remnants. We utilize axisymmetric, cold equilibrium models with differential rotation and determine the threshold mass for collapse from the comparison between an empirical relation that describes the angular momentum in the remnant for a given total binary mass and the sequence of rotating equilibrium models at the threshold to collapse (the latter assumed to be near the turning point of fixed-angular-momentum sequences). In spite of the various simplifying assumptions, the empirical relation for prompt collapse is reproduced with good accuracy, which demonstrates its robustness. We discuss implications of our methodology and results for understanding other empirical relations satisfied by neutron-star merger remnants that have been discovered by numerical simulations and that play a key role in constraining the high-density equation of state through gravitational-wave observations.

The merger of two neutron stars may give birth to a long-lived hypermassive neutron star. If it harbours a strong magnetic field of magnetar strength, its spin-down could explain several features of short gamma-ray burst afterglows. The magnetorotational instability (MRI) has been proposed as a mechanism to amplify the magnetic field to the required strength. Previous studies have, however, neglected neutrinos, which may have an important impact on the MRI by inducing a viscosity and drag. We investigate the impact of these neutrino effects on the linear growth of the MRI by applying a local stability analysis to snapshots of a neutron star merger simulation. We find that neutrinos have a significant impact inside the hypermassive neutron star, but have at most a marginal effect in the torus surrounding it. Inside the hypermassive neutron star, the MRI grows in different regimes depending on the radius and on the initial magnetic field strength. For magnetic fields weaker than $10^{13}-10^{14}\,{\rm G}$, the growth rate of the MRI is significantly reduced due to the presence of neutrinos. We conclude that neutrinos should be taken into account when studying the growth of the MRI from realistic initial magnetic fields. Current numerical simulations, which neglect neutrino viscosity, are only consistent, i.e. in the adopted ideal regime, if they start from artificially strong initial magnetic fields above $\sim10^{14}\,{\rm G}$. One should be careful when extrapolating these results to lower initial magnetic fields, where the MRI growth is strongly affected by neutrino viscosity or drag.

The coalescence of a neutron-star binary is likely to result in the formation of a neutron-star merger remnant for a large range of binary mass configurations. The massive merger remnant shows strong oscillations, which are excited by the merging process, and emits gravitational waves. Here we discuss possibilities and prospects of inferring unknown stellar properties of neutron stars by the detection of postmerger gravitational-wave emission, which thus leads to constraints of the equation of state of high-density matter. In particular, the dominant oscillation frequency of the postmerger remnant provides tight limits to neutron-star radii. We mention first steps towards a practical implementation of future gravitational-wave searches for the postmerger emission. Moreover, we outline possibilities to estimate the unknown maximum mass of nonrotating neutron stars from such types of measurements. Finally, we review the origin and scientific implications of secondary peaks in the gravitational-wave spectrum of neutron-star mergers and differences in the dynamical behavior of the postmerger remnant depending on the binary configuration. These considerations lead to a unified picture of the post-merger gravitational-wave emission and the post-merger dynamics.

We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star (NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.

Remnants of neutron-star mergers are essentially massive, hot, differentially rotating neutron stars, which are initially strongly oscillating. They represent a unique probe for high-density matter because the oscillations are detectable via gravitational-wave measurements and are strongly dependent on the equation of state. The impact of the equation of state is apparent in the frequency of the dominant oscillation mode of the remnant. For a fixed total binary mass a tight relation between the dominant postmerger frequency and the radii of nonrotating neutron stars exists. Inferring observationally the dominant postmerger frequency thus determines neutron star radii with high accuracy of the order of a few hundred meters. By considering symmetric and asymmetric binaries of the same chirp mass, we show that the knowledge of the binary mass ratio is not critical for this kind of radius measurements. We summarize different possibilities to deduce the maximum mass of nonrotating neutron stars. We clarify the nature of the three most prominent features of the postmerger gravitational-wave spectrum and argue that the merger remnant can be considered to be a single, isolated, self-gravitating object that can be described by concepts of asteroseismology. The understanding of the different mechanisms shaping the gravitational-wave signal yields a physically motivated analytic model of the gravitational-wave emission, which may form the basis for template-based gravitational-wave data analysis. We explore the observational consequences of a scenario of two families of compact stars including hadronic and quark stars. We find that this scenario leaves a distinctive imprint on the postmerger gravitational-wave signal. In particular, a strong discontinuity in the dominant postmerger frequency as function of the total mass will be a strong indication for two families of compact stars. (abridged)

We present a comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

We have performed r-process calculations for matter ejected dynamically in neutron star mergers based on a complete set of trajectories from a three-dimensional relativistic smoothed particle hydrodynamic simulation. Our calculations consider an extended nuclear network, including spontaneous, $\beta$- and neutron-induced fission and adopting fission yield distributions from the ABLA code. We have studied the sensitivity of the r-process abundances to nuclear masses by using different models. Most of the trajectories, corresponding to 90% of the ejected mass, follow a relatively slow expansion allowing for all neutrons to be captured. The resulting abundances are very similar to each other and reproduce the general features of the observed r-process abundance (the second and third peaks, the rare-earth peak and the lead peak) for all mass models as they are mainly determined by the fission yields. We find distinct differences in the abundance yields at and just above the third peak, which can be traced back to different predictions of neutron separation energies for r-process nuclei around neutron number $N=130$. The remaining trajectories, which contribute 10% by mass to the total integrated abundances, follow such a fast expansion that the r process does not use all the neutrons. This also leads to a larger variation of abundances among trajectories as fission does not dominate the r-process dynamics. The total integrated abundances are dominated by contributions from the slow abundances and hence reproduce the general features of the observed r-process abundances. We find that at timescales of weeks relevant for kilonova light curve calculations, the abundance of actinides is larger than the one of lanthanides. Hence actinides can be even more important than lanthanides to determine the photon opacities under kilonova conditions. (Abridged)

We present an effective, low-dimensionality frequency-domain template for the gravitational wave signal from the stellar remnants from binary neutron star coalescence. A principal component decomposition of a suite of numerical simulations of binary neutron star mergers is used to construct orthogonal basis functions for the amplitude and phase spectra of the waveforms for a variety of neutron star equations of state and binary mass configurations. We review the phenomenology of late merger / post-merger gravitational wave emission in binary neutron star coalescence and demonstrate how an understanding of the dynamics during and after the merger leads to the construction of a universal spectrum. We also provide a discussion of the prospects for detecting the post-merger signal in future gravitational wave detectors as a potential contribution to the science case for third generation instruments. The template derived in our analysis achieves $>90\%$ match across a wide variety of merger waveforms and strain sensitivity spectra for current and potential gravitational wave detectors. A Fisher matrix analysis yields a preliminary estimate of the typical uncertainty in the determination of the dominant post-merger oscillation frequency $f_{\mathrm{peak}}$ as $\delta f_{\mathrm{peak}} \sim 50$Hz. Using recently derived correlations between $f_{\mathrm{peak}}$ and the neutron star radii, this suggests potential constraints on the radius of a fiducial neutron star of $\sim 220$\,m. Such measurements would only be possible for nearby ($\sim 30$Mpc) sources with advanced LIGO but become more feasible for planned upgrades to advanced LIGO and other future instruments, leading to constraints on the high density neutron star equation of state which are independent and complementary to those inferred from the pre-merger inspiral gravitational wave signal.

We introduce a classification scheme of the post-merger dynamics and gravitational-wave emission in binary neutron star mergers, after identifying a new mechanism by which a secondary peak in the gravitational-wave spectrum is produced. It is caused by a spiral deformation, the pattern of which rotates slower with respect to the double-core structure in center of the remnant. This secondary peak is typically well separated in frequency from the secondary peak produced by a nonlinear interaction between a quadrupole and a quasi-radial oscillation. The new mechanism allows for an explanation of low-frequency modulations seen in a number of physical characteristics of the remnant, such as the central lapse function, the maximum density and the separation between the two cores. We find empirical relations for both types of secondary peaks between their gravitational-wave frequency and the compactness of nonrotating individual neutron stars, that exist for fixed total binary masses. These findings are derived for equal-mass binaries without intrinsic neutron-star spin analyzing hydrodynamical simulations without magnetic field effects. Our classification scheme may form the basis for the construction of detailed gravitational-wave templates of the post-merger phase. We find that the quasi-radial oscillation frequency of the remnant decreases with the total binary mass. For a given merger event our classification scheme may allow to determine the proximity of the measured total binary mass to the threshold mass for prompt black hole formation, which can, in turn, yield an estimate of the maximum neutron-star mass.

We investigate beta-interactions of free nucleons and their impact on the electron fraction (Y_e) and r-process nucleosynthesis in ejecta characteristic of binary neutron star mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model to represent the density-temperature evolutions in our parametric study. In the high-density environment, positron captures decrease the neutron richness at the high temperatures predicted by the hydrodynamic simulation. Circumventing the complexities of modelling three-dimensional neutrino transport, (anti)neutrino captures are parameterized in terms of prescribed neutrino luminosities and mean energies, guided by published results and assumed as constant in time. Depending sensitively on the adopted neutrino-antineutrino luminosity ratio, neutrino processes increase Y_e to values between 0.25 and 0.40, still allowing for a successful r-process compatible with the observed solar abundance distribution and a significant fraction of the ejecta consisting of r-process nuclei. If the electron neutrino luminosities and mean energies are relatively large compared to the antineutrino properties, the mean Y_e might reach values >0.40 so that neutrino captures seriously compromise the success of the r-process. In this case, the r-abundances remain compatible with the solar distribution, but the total amount of ejected r-material is reduced to a few percent, because the production of iron-peak elements is favored. Proper neutrino physics, in particular also neutrino absorption, have to be included in BNSM simulations before final conclusions can be drawn concerning r-processing in this environment and concerning observational consequences like kilonovae, whose peak brightness and color temperature are sensitive to the composition-dependent opacity of the ejecta.

The oscillations of a merger remnant forming after the coalescence of two neutron stars are very characteristic for the high-density equation of state. The dominant oscillation frequency occurs as a pronounced peak in the kHz range of the gravitational-wave spectrum. We describe how the dominant oscillation frequency of the remnant can be employed to infer the radii of non-rotating neutron stars.

We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black-hole (BH)-torus systems. The evolution of the BH-accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS-NS) and NS-BH merger models, producing ~3-6 Msun BHs with rotation parameters of A~0.8 and tori of 0.03-0.3 Msun. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH-torus systems. While typically ~20-25% of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ~1%, but neutrino heating enhances the viscous ejecta significantly. Since BH-torus ejecta possess a wide distribution of electron fractions (0.1-0.6) and entropies, they produce heavy elements from A~80 up to the actinides, with relative contributions of A>130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH-torus phases can reproduce the solar abundances amazingly well for A>90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40<Z<56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS-BH mergers, might explain the composition of heavy-element deficient stars.

The merger of binary neutron stars (NSs) ejects a small quantity of neutron rich matter, the radioactive decay of which powers a day to week long thermal transient known as a kilonova. Most of the ejecta remains sufficiently dense during its expansion that all neutrons are captured into nuclei during the r-process. However, recent general relativistic merger simulations by Bauswein and collaborators show that a small fraction of the ejected mass (a few per cent, or ~1e-4 Msun) expands sufficiently rapidly for most neutrons to avoid capture. This matter originates from the shocked-heated interface between the merging NSs. Here we show that the beta-decay of these free neutrons in the outermost ejecta powers a `precursor' to the main kilonova emission, which peaks on a timescale of a few hours following merger at U-band magnitude ~22 (for an assumed distance of 200 Mpc). The high luminosity and blue colors of the neutron precursor render it a potentially important counterpart to the gravitational wave source, that may encode valuable information on the properties of the merging binary (e.g. NS-NS versus NS-black hole) and the NS equation of state. Future work is necessary to assess the robustness of the fast moving ejecta and the survival of free neutrons in the face of neutrino absorptions, although the precursor properties are robust to a moderate amount of leptonization. Our results provide additional motivation for short latency gravitational wave triggers and rapid follow-up searches with sensitive ground based telescopes.

We derive constraints on the time-averaged event rate of neutron star-black hole (NS-BH) mergers by using estimates of the population-integrated production of heavy rapid neutron-capture (r-process) elements with nuclear mass numbers A > 140 by such events in comparison to the Galactic repository of these chemical species. Our estimates are based on relativistic hydrodynamical simulations convolved with theoretical predictions of the binary population. This allows us to determine a strict upper limit of the average NS-BH merger rate of ~6*10^-5 per year. We quantify the uncertainties of this estimate to be within factors of a few mostly because of the unknown BH spin distribution of such systems, the uncertain equation of state of NS matter, and possible errors in the Galactic content of r-process material. Our approach implies a correlation between the merger rates of NS-BH binaries and of double NS systems. Predictions of the detection rate of gravitational-wave signals from such compact-object binaries by Advanced LIGO and Advanced Virgo on the optimistic side are incompatible with the constraints set by our analysis.

The equation of state plays a critical role in the physics of the merger of two neutron stars. Recent numerical simulations with microphysical equation of state suggest the outcome of such events depends on the mass of the neutron stars. For less massive systems, simulations favor the formation of a hypermassive, quasi-stable neutron star, whose oscillations produce a short, high frequency burst of gravitational radiation. Its dominant frequency content is tightly correlated with the radius of the neutron star, and its measurement can be used to constrain the supranuclear equation of state. In contrast, the merger of higher mass systems results in prompt gravitational collapse to a black hole. We have developed an algorithm which combines waveform reconstruction from a morphology-independent search for gravitational wave transients with Bayesian model selection, to discriminate between post-merger scenarios and accurately measure the dominant oscillation frequency. We demonstrate the efficacy of the method using a catalogue of simulated binary merger signals in data from LIGO and Virgo, and we discuss the prospects for this analysis in advanced ground-based gravitational wave detectors. From the waveforms considered in this work and assuming an optimally oriented source, we find that the post-merger neutron star signal may be detectable by this technique to $\sim 10\text{--}25$\,Mpc. We also find that we successfully discriminate between the post-merger scenarios with $\sim 95\%$ accuracy and determine the dominant oscillation frequency of surviving post-merger neutron stars to within $\sim 10$\,Hz, averaged over all detected signals. This leads to an uncertainty in the estimated radius of a non-rotating 1.6\,M$_{\odot}$ reference neutron star of $\sim 100\,$m.

We present a novel method for revealing the equation of state of high-density neutron star matter through gravitational waves emitted during the postmerger phase of a binary neutron star system. The method relies on a small number of detections of the peak frequency in the postmerger phase for binaries of different (relatively low) masses, in the most likely range of expected detections. From such observations, one can construct the derivative of the peak frequency versus the binary mass, in this mass range. Through a detailed study of binary neutron star mergers for a large sample of equations of state, we show that one can extrapolate the above information to the highest possible mass (the threshold mass for black hole formation in a binary neutron star merger). In turn, this allows for an empirical determination of the maximum mass of cold, nonrotating neutron stars to within 0.1 M_sun, while the corresponding radius is determined to within a few percent. Combining this with the determination of the radius of cold, nonrotating neutron stars of 1.6 M_sun (to within a few percent, as was demonstrated in Bauswein et al., PRD, 86, 063001, 2012), allows for a clear distinction of a particular candidate equation of state among a large set of other candidates. Our method is particularly appealing because it reveals simultaneously the moderate and very high-density parts of the equation of state, enabling the distinction of mass-radius relations even if they are similar at typical neutron star masses. Furthermore, our method also allows to deduce the maximum central energy density and maximum central rest-mass density of cold, nonrotating neutron stars with an accuracy of a few per cent.

Neutron star (NS) merger ejecta offer a viable site for the production of heavy r-process elements with nuclear mass numbers A > 140. The crucial role of fission recycling is responsible for the robustness of this site against many astrophysical uncertainties, but calculations sensitively depend on nuclear physics. In particular the fission fragment yields determine the creation of 110 < A < 170 nuclei. Here we apply a new scission-point model, called SPY, to derive the fission fragment distribution (FFD) of all relevant neutron-rich, fissioning nuclei. The model predicts a doubly asymmetric FFD in the abundant A ~ 278 mass region that is responsible for the final recycling of the fissioning material. Using ejecta conditions based on relativistic NS merger calculations we show that this specific FFD leads to a production of the A ~ 165 rare-earth peak that is nicely compatible with the abundance patterns in the Sun and metal-poor stars. This new finding further strengthens the case of NS mergers as possible dominant origin of r-nuclei with A > 140.

We perform hydrodynamical simulations of neutron-star mergers for a large sample of temperature-dependent, nuclear equations of state, and determine the threshold mass above which the merger remnant promptly collapses to form a black hole. We find that, depending on the equation of state, the threshold mass is larger than the maximum mass of a non-rotating star in isolation by between 30 and 70 per cent. Our simulations also show that the ratio between the threshold mass and maximum mass is tightly correlated with the compactness of the non-rotating maximum-mass configuration. We speculate on how this relation can be used to derive constraints on neutron-star properties from future observations.

We investigate systematically the dynamical mass ejection, r-process nucleosynthesis, and properties of electromagnetic counterparts of neutron-star (NS) mergers in dependence on the uncertain properties of the nuclear equation of state (EoS) by employing 40 representative, microphysical high-density EoSs in relativistic, hydrodynamical simulations. The crucial parameter determining the ejecta mass is the radius R_1.35 of a 1.35 M_sun NS. NSs with smaller R_1.35 ("soft" EoS) eject systematically higher masses. These range from ~10^-3 M_sun to ~10^-2 M_sun for 1.35-1.35 M_sun binaries and from ~5*10^-3 M_sun to ~2*10^-2 M_sun for 1.2-1.5 M_sun systems (with kinetic energies between ~5*10^49 erg and 10^51 erg). Correspondingly, the bolometric peak luminosities of the optical transients of symmetric (asymmetric) mergers vary between 3*10^41 erg/s and 14*10^41 erg/s (9*10^41 erg/s and 14.5*10^41 erg/s) on timescales between ~2 h and ~12 h. If these signals with absolute bolometric magnitudes from -15.0 to -16.7 are measured, the tight correlation of their properties with those of the merging NSs might provide valuable constraints on the high-density EoS. The r-process nucleosynthesis exhibits a remarkable robustness independent of the EoS, producing a nearly solar abundance pattern above mass number 130. By the r-process content of the Galaxy and the average production per event the Galactic merger rate is limited to 4*10^-5/yr (4*10^-4/yr) for a soft (stiff) NS EoS, if NS mergers are the main source of heavy r-nuclei. The production ratio of radioactive 232Th to 238U attains a stable value of 1.64-1.67, which does not exclude NS mergers as potential sources of heavy r-material in the most metal-poor stars.

Neutron-star (NS) merger simulations are conducted for 38 representative microphysical descriptions of high-density matter in order to explore the equation-of-state dependence of the postmerger ring-down phase. The formation of a deformed, oscillating, differentially rotating very massive NS is the typical outcome of the coalescence of two stars with 1.35 $M_{\odot}$ for most candidate EoSs. The oscillations of this object imprint a pronounced peak in the gravitational-wave (GW) spectra, which is used to characterize the emission for a given model. The peak frequency of this postmerger GW signal correlates very well with the radii of nonrotating NSs, and thus allows to constrain the high-density EoS by a GW detection. In the case of 1.35-1.35 $M_{\odot}$ mergers the peak frequency scales particularly well with the radius of a NS with 1.6 $M_{\odot}$, where the maximum deviation from this correlation is only 60 meters for fully microphysical EoSs which are compatible with NS observations. Combined with the uncertainty in the determination of the peak frequency it appears likely that a GW detection can measure the radius of a 1.6 $M_{\odot}$ NS with an accuracy of about 100 to 200 meters. We also uncover relations of the peak frequency with the radii of nonrotating NSs with 1.35 $M_{\odot}$ or 1.8 $M_{\odot}$, with the radius or the central energy density of the maximum-mass Tolman-Oppenheimer-Volkoff configuration, and with the pressure or sound speed at a fiducial rest-mass density of about twice nuclear saturation density. Furthermore, it is found that a determination of the dominant postmerger GW frequency can provide an upper limit for the maximum mass of nonrotating NSs. The prospects for a detection of the postmerger GW signal and a determination of the dominant GW frequency are estimated to be in the range of 0.015 to 1.2 events per year with the upcoming Advanced LIGO detector.

We demonstrate by a large set of merger simulations for symmetric binary neutron stars (NSs) that there is a tight correlation between the frequency peak of the postmerger gravitational-wave (GW) emission and the physical properties of the nuclear equation of state (EoS), e.g. expressed by the radius of the maximum-mass Tolman-Oppenheimer-Volkhoff configuration. Therefore, a single measurement of the peak frequency of the postmerger GW signal will constrain the NS EoS significantly. For plausible optimistic merger-rate estimates a corresponding detection with Advanced LIGO is likely to happen within an operation time of roughly a year.

Although the rapid neutron-capture process, or r-process, is fundamentally important for explaining the origin of approximately half of the stable nuclei with A > 60, the astrophysical site of this process has not been identified yet. Here we study r-process nucleosynthesis in material that is dynamically ejected by tidal and pressure forces during the merging of binary neutron stars (NSs) and within milliseconds afterwards. For the first time we make use of relativistic hydrodynamical simulations of such events, defining consistently the conditions that determine the nucleosynthesis, i.e., neutron enrichment, entropy, early density evolution and thus expansion timescale, and ejecta mass. We find that 10^{-3}-10^{-2} solar masses are ejected, which is enough for mergers to be the main source of heavy (A > 140) galactic r-nuclei for merger rates of some 10^{-5} per year. While asymmetric mergers eject 2-3 times more mass than symmetric ones, the exact amount depends weakly on whether the NSs have radii of ~15 km for a "stiff" nuclear equation of state (EOS) or ~12 km for a "soft" EOS. R-process nucleosynthesis during the decompression becomes largely insensitive to the detailed conditions because of efficient fission recycling, producing a composition that closely follows the solar r-abundance distribution for nuclei with mass numbers A > 140. Estimating the light curve powered by the radioactive decay heating of r-process nuclei with an approximative model, we expect high emission in the B-V-R bands for 1-2 days with potentially observable longer duration in the case of asymmetric mergers because of the larger ejecta mass.

We study the excitation of nonaxisymmetric modes in the post-merger phase of binary compact object mergers and the associated gravitational wave emission. Our analysis is based on general-relativistic simulations, in the spatial conformal flatness approximation, using smoothed-particle-hydrodynamics for the evolution of matter, and we use a set of equal and unequal mass models, described by two nonzero-temperature hadronic equations of state and by one strange star equation of state. Through Fourier transforms of the evolution of matter variables, we can identify a number of oscillation modes, as well as several nonlinear components (combination frequencies). We focus on the dominant m=2 mode, which forms a triplet with two nonlinear components that are the result of coupling to the quasiradial mode. A corresponding triplet of frequencies is identified in the gravitational wave spectrum, when the individual masses of the compact objects are in the most likely range of 1.2 to 1.35 $M_\odot$. We can thus associate, through direct analysis of the dynamics of the fluid, a specific frequency peak in the gravitational wave spectrum with the nonlinear component resulting from the difference between the m=2 mode and the quasiradial mode. Once such observations become available, both the m=2 and quasiradial mode frequencies could be extracted, allowing for the application of gravitational-wave asteroseismology to the post-merger remnant and leading to tight constraints on the equation of state of high-density matter.

We perform three-dimensional relativistic hydrodynamical calculations of neutron star mergers to assess the reliability of an approximate treatment of thermal effects in such simulations by combining an ideal-gas component with zero-temperature, micro-physical equations of state. To this end we compare the results of simulations that make this approximation to the outcome of models with a consistent treatment of thermal effects in the equation of state. In particular we focus on the implications for observable consequences of merger events like the gravitational-wave signal. It is found that the characteristic gravitational-wave oscillation frequencies of the post-merger remnant differ by about 50 to 250 Hz (corresponding to frequency shifts of 2 to 8 per cent) depending on the equation of state and the choice of the characteristic index of the ideal-gas component. In addition, the delay time to black hole collapse of the merger remnant as well as the amount of matter remaining outside the black hole after its formation are sensitive to the description of thermal effects.

We perform three-dimensional relativistic hydrodynamical simulations of the coalescence of strange stars (SSs) and explore the possibility to decide on the strange matter hypothesis by means of gravitational-wave (GW) measurements. Selfbinding of strange quark matter (SQM) and the generally more compact stars yield features that clearly distinguish SS from neutron star (NS) mergers, e.g. hampering tidal disruption during the plunge of quark stars. Furthermore, instead of forming dilute halos around the remnant as in the case of NS mergers, the coalescence of SSs results in a differentially rotating hypermassive object with a sharp surface layer surrounded by a geometrically thin, clumpy high-density SQM disk. We also investigate the importance of including non-zero temperature equations of state (EoSs) in NS and SS merger simulations. In both cases we find a crucial sensitivity of the dynamics and outcome of the coalescence to thermal effects, which, e.g., determine the outer remnant structure and the delay time of the dense remnant core to black hole collapse. For comparing and classifying the GW signals, we use a number of characteristic quantities like the maximum frequency during inspiral or the dominant frequency of oscillations of the postmerger remnant. In general, these frequencies are higher for SS mergers. If not, additional features of the GW luminosity spectrum may help to discriminate coalescence events of the different types. Future GW measurements may thus help to decide on the existence of SQM stars. (abridged)

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