Search

The time until black hole formation in a binary neutron-star (NS) merger contains invaluable information about the nuclear equation of state (EoS) but has thus far been difficult to measure. We propose a new way to constrain the merger remnant's NS lifetime, which is based on the tendency of the NS remnant neutrino-driven winds to enrich the ejected material with helium. Based on the He I $\lambda 1083.3$ nm line, we show that the feature around 800-1200 nm in AT2017gfo at 4.4 days seems inconsistent with a helium mass fraction of $X_{\mathrm{He}} \gtrsim 0.05$ in the polar ejecta. Recent neutrino-hydrodynamic simulations of merger remnants are only compatible with this limit if the NS remnant collapses within 20-30 ms. Such a short lifetime implies that the total binary mass of GW170817, $M_\mathrm{\rm tot}$, lay close to the threshold binary mass for direct gravitational collapse, $M_\mathrm{thres}$, for which we estimate $M_{\mathrm{thres}}\lesssim 2.93 M_\odot$. This upper bound on $M_\mathrm{thres}$ yields upper limits on the radii and maximum mass of cold, non-rotating NSs, which rule out simultaneously large values for both quantities. In combination with causality arguments, this result implies a maximum NS mass of $M_\mathrm{max}\lesssim2.3 M_\odot$. The combination of all limits constrains the radii of 1.6 M$_\odot$ NSs to about 12$\pm$1 km for $M_\mathrm{max}$ = 2.0 M$_\odot$ and 11.5$\pm$1 km for $M_\mathrm{max}$ = 2.15 M$_\odot$. This $\sim2$ km allowable range then tightens significantly for $M_\mathrm{max}$ above $\approx2.15$ M$_\odot$. This rules out a significant number of current EoS models. The short NS lifetime also implies that a black-hole torus, not a highly magnetized NS, was the central engine powering the relativistic jet of GRB170817A. Our work motivates future developments... [abridged], Comment: 30 pages, 16 figures, 2 tables, submitted to PRX

We perform an extensive study of equation of state (EoS) models featuring a phase transition from hadronic to deconfined quark matter in neutron star merger simulations. We employ three different hadronic EoSs, a constant speed of sound parameterization for the quark phase and a Maxwell construction to generate a large sample of hybrid EoS models. We systematically vary the onset density and density jump of the phase transition as well as the quark matter stiffness and simulate binary neutron star mergers to infer how the properties of the phase transition affect the gravitational-wave signal. In total we simulate mergers with 245 different hybrid EoS models. In particular, we explore in which scenarios a phase transition would be detectable by a characteristically increased postmerger gravitational-wave frequency compared to an estimate from the inspiral signal assuming a purely hadronic EoS. We find that the density jump at the transition (latent heat) has the largest impact on the gravitational-wave frequencies, while the influence of the stiffness of quark matter is smaller. We quantify which range of phase transition properties would be compatible with a certain magnitude or absence of the gravitational-wave postmerger frequency shift. By means of these dependencies, a future detection will thus directly yield constraints on the allowed features of the hadron-quark phase transition., Comment: 15 pages, 12 figures

We present a time-domain model for the gravitational waves emitted by equal-mass binary neutron star merger remnants for a fixed equation of state. We construct a large set of numerical relativity simulations for a single equation of state consistent with current constraints, totaling 157 equal-mass binary neutron star merger configurations. The gravitational-wave model is constructed using the supervised learning method of K-nearest neighbor regression. As a first step toward developing a general model with supervised learning methods that accounts for the dependencies on equation of state and the binary masses of the system, we explore the impact of the size of the dataset on the model. We assess the accuracy of the model for a varied dataset size and number density in total binary mass. Specifically, we consider five training sets of $\{ 20,40, 60, 80, 100\}$ simulations uniformly distributed in total binary mass. We evaluate the resulting models in terms of faithfulness using a test set of 30 additional simulations that are not used during training and which are equidistantly spaced in total binary mass. The models achieve faithfulness with maximum values in the range of $0.980$ to $0.995$. We assess our models simulating signals observed by the three-detector network of Advanced LIGO-Virgo. We find that all models with training sets of size equal to or larger than $40$ achieve an unbiased measurement of the main gravitational-wave frequency. We confirm that our results do not depend qualitatively on the choice of the (fixed) equation of state. We conclude that training sets, with a minimum size of $40$ simulations, or a number density of approximately $11$ simulations per $0.1\,M_\odot$ of total binary mass, suffice for the construction of faithful templates for the post-merger signal for a single equation of state and equal-mass binaries (abbreviated)., Comment: 20 pages, 12 figures, submitted to Phys. Rev. D

Recent analysis of the kilonova, AT2017gfo, has indicated that this event was highly spherical. This may challenge hydrodynamics simulations of binary neutron star mergers, which usually predict a range of asymmetries, and radiative transfer simulations show a strong direction dependence. Here we investigate whether the synthetic spectra from a 3D kilonova simulation of asymmetric ejecta from a hydrodynamical merger simulation can be compatible with the observational constraints suggesting a high degree of sphericity in AT2017gfo. Specifically, we determine whether fitting a simple P-Cygni line profile model leads to a value for the photospheric velocity that is consistent with the value obtained from the expanding photosphere method. We would infer that our kilonova simulation is highly spherical at early times, when the spectra resemble a blackbody distribution. The two independently inferred photospheric velocities can be very similar, implying a high degree of sphericity, which can be as spherical as inferred for AT2017gfo, demonstrating that the photosphere can appear spherical even for asymmetrical ejecta. The last-interaction velocities of radiation escaping the simulation show a high degree of sphericity, supporting the inferred symmetry of the photosphere. We find that when the synthetic spectra resemble a blackbody the expanding photosphere method can be used to obtain an accurate luminosity distance (within 4-7 per cent)., Comment: Accepted by MNRAS

We provide the first comprehensive study of hyperons in neutron star mergers and quantify their specific impact. We discuss the thermal behavior of hyperonic equations of state~(EoSs) as a distinguishing feature from purely nucleonic models in the remnants of binary mergers using a large set of numerical simulations. Finite temperature enhances the production of hyperons, which leads to a reduced pressure as highly degenerate nucleons are depopulated. This results in a characteristic increase of the dominant postmerger gravitational-wave frequency by up to $\sim150$~Hz compared to purely nucleonic EoS models. By our comparative approach we can directly link this effect to the occurrence of hyperons. Although this feature is generally weak, it is in principle measurable if the EoS and stellar parameters of cold neutron stars are sufficiently well determined. Considering that the mass-radius relations of purely nucleonic and hyperonic EoSs may be indistinguishable and the overall challenge to infer the presence of hyperons in neutron stars, these findings are important as a new route to answer the outstanding question about hyperonic degrees of freedom in high-density matter., Comment: 9 pages, 7 figures, 2 tables; published in PRD

We present three-dimensional radiative transfer calculations for the ejecta from a neutron star merger that include line-by-line opacities for tens of millions of bound-bound transitions, composition from an r-process nuclear network, and time-dependent thermalization of decay products from individual $\alpha$ and $\beta^-$ decay reactions. In contrast to expansion opacities and other wavelength-binned treatments, a line-by-line treatment enables us include fluorescence effects and associate spectral features with the emitting and absorbing lines of individual elements. We find variations in the synthetic observables with both the polar and azimuthal viewing angles. The spectra exhibit blended features with strong interactions by Ce III, Sr II, Y II, and Zr II that vary with time and viewing direction. We demonstrate the importance of wavelength-calibration of atomic data using a model with calibrated Sr, Y, and Zr data, and find major differences in the resulting spectra, including a better agreement with AT2017gfo. The synthetic spectra for near-polar inclination show a feature at around 8000 A, similar to AT2017gfo. However, they evolve on a more rapid timescale, likely due to the low ejecta mass (0.005 M$_\odot$) as we take into account only the early ejecta. The comparatively featureless spectra for equatorial observers gives a tentative prediction that future observations of edge-on kilonovae will appear substantially different from AT2017gfo. We also show that 1D models obtained by spherically averaging the 3D ejecta lead to dramatically different direction-integrated luminosities and spectra compared to full 3D calculations., Comment: 12 pages, 5 figures. Accepted by ApJL

Binary neutron star mergers are the first confirmed site of rapid neutron capture (r-process) element nucleosynthesis. The kilonova AT2017gfo is the only electromagnetic counterpart of a neutron star merger spectroscopically observed. We analyse the entire spectral sequence of AT2017gfo (from merger to +10.4 days) and identify seven emission-like features. We confirm that the prominent 1.08um feature can be explained by the Sr II near-infrared triplet evolving from a P-Cygni profile through to pure emission. We calculate the expected strength of the [Sr II] doublet and show that its absence requires highly clumped ejecta. Near-infrared features at 1.58 and 2.07um emerge after three days and become more prominent as the spectra evolve. We model these as optically thick P-Cygni profiles and alternatively as pure emission features (with FWHM = 35600 +/- 6600 km/s), and favour the latter interpretation. The profile of the strong 2.07um emission feature is best reproduced with two lines, centred at 2.059 and 2.135um. We search for candidate ions for all prominent features in the spectra. Strong, permitted transitions of La III, Ce III, Gd III, Ra II and Ac I are plausible candidates for the emission features. If any of these features are produced by intrinsically weak, forbidden transitions, we highlight candidate ions spanning the three r-process peaks. The second r-process peak elements Te and I have plausible matches to multiple features. We highlight the need for more detailed and quantitative atomic line transition data., Comment: Main text: 19 pages, 7 figures, 6 tables. Appendices: 10 pages, 3 tables. Submitted to MNRAS. Comments welcome

While gravitational wave (GW) standard sirens from neutron star (NS) mergers have been proposed to offer good measurements of the Hubble constant, we show in this paper how a variation of the expanding photosphere method (EPM) or spectral-fitting expanding atmosphere method, applied to the kilonovae (KNe) associated with the mergers, can provide an independent distance measurement to individual mergers that is potentially accurate to within a few percent. There are four reasons why the KN-EPM overcomes the major uncertainties commonly associated with this method in supernovae: 1) the early continuum is very well-reproduced by a blackbody spectrum, 2) the dilution effect from electron scattering opacity is likely negligible, 3) the explosion times are exactly known due to the GW detection, and 4) the ejecta geometry is, at least in some cases, highly spherical and can be constrained from line-shape analysis. We provide an analysis of the early VLT/X-shooter spectra AT2017gfo showing how the luminosity distance can be determined, and find a luminosity distance of $D_L = 44.5\pm0.8$ Mpc in agreement with, but more precise than, previous methods. We investigate the dominant systematic uncertainties, but our simple framework, which assumes a blackbody photosphere, does not account for the full time-dependent three-dimensional radiative transfer effects, so this distance should be treated as preliminary. The luminosity distance corresponds to an estimated Hubble constant of $H_0 = 67.0\pm 3.6$ km s$^{-1}$ Mpc$^{-1}$, where the dominant uncertainty is due to the modelling of the host peculiar velocity. We also estimate the expected constraints on $H_0$ from future KN-EPM-analysis with the upcoming O4 and O5 runs of the LIGO collaboration GW-detectors, where five to ten similar KNe would yield 1\% precision cosmological constraints., Comment: A&A, In press

We study the importance of the thermal behavior of the hadron-quark phase transition in neutron star (NS) mergers. To this end, we devise a new scheme approximating thermal effects to supplement any cold, barotropic hybrid equation of state (EoS) model, i.e. two-phase EoS constructions with a hadronic regime and a phase of deconfined quark matter. The consideration of temperature-dependent phase boundaries turns out to be critical for a quantitative description of quark matter effects in NS mergers, since the coexistence phase can introduce a strong softening of the EoS at finite temperature, which is even more significant than the change of the EoS by the phase transition at T=0. We validate our approach by comparing to existing fully temperature-dependent EoS models and find a very good quantitative agreement of postmerger gravitational-wave (GW) features. Simulations with the commonly-used thermal ideal-gas approach exhibit sizable differences compared to full hybrid models implying that its use in NS merger simulations with quark matter is problematic. Our new scheme provides the means to isolate thermal effects of quark matter from the properties of the cold hybrid EoS and thus allows an assessment of the thermal behavior alone. We show that different shapes of the phase boundaries at finite temperature can have a large impact on the postmerger dynamics and GW signal for the same cold hybrid model. This finding demonstrates that postmerger GW emission contains important complementary information compared to properties extracted from cold stars. We also show by concrete examples that it is even possible for quark matter to only occur and thus be detectable in finite-temperature systems like merger remnants but not in cold NSs (abbreviated)., Comment: 28 pages, 16 figures, revised version, published by Phys. Rev. D

We are witnessing the dawn of gravitational wave (GW) astronomy. With currently available detectors, observations are restricted to GW frequencies in the range between ${\sim} 10\,\mathrm{Hz}$ and $10\,\mathrm{kHz}$, which covers the signals from mergers of compact objects. The launch of the space observatory LISA will open up a new frequency band for the detection of stellar interactions at lower frequencies. In this work, we predict the shape and strength of the GW signals associated with common-envelope interaction and merger events in binary stars, and we discuss their detectability. Previous studies estimated these characteristics based on semi-analytical models. In contrast, we used detailed three-dimensional magnetohydrodynamic simulations to compute the GW signals. We show that for the studied models, the dynamical phase of common-envelope events and mergers between main-sequence stars lies outside of the detectability band of the LISA mission. We find, however, that the final stages of common-envelope interactions leading to mergers of the stellar cores fall into the frequency band in which the sensitivity of LISA peaks, making them promising candidates for detection. These detections can constrain the enigmatic common-envelope dynamics. Furthermore, future decihertz observatories such as DECIGO or BBO would also be able to observe this final stage and the post-merger signal, through which we might be able to detect the formation of Thorne-\.Zytkow objects., Comment: Accepted for publication in A&A, 12 pages, 8 figures

We study the impact of pions in simulations of neutron star mergers and explore the impact on gravitational-wave observables. We model charged and neutral pions as a non-interacting Boson gas with a chosen, constant effective mass. We add the contributions of pions, which can occur as a condensate or as a thermal population, to existing temperature and composition dependent equations of state. Compared to the models without pions, the presence of a pion condensate decreases the characteristic properties of cold, non-rotating neutron stars such as the maximum mass, the radius and the tidal deformability. We conduct relativistic hydrodynamical simulations of neutron star mergers for these modified equations of state models and compare to the original models, which ignore pions. Generally, the inclusion of pions leads to a softening of the equation of state, which is more pronounced for smaller effective pion masses. We find a shift of the dominant postmerger gravitational-wave frequency by up to 150~Hz to higher frequencies and a reduction of the threshold binary mass for prompt black-hole formation by up to 0.07~$M_\odot$. We evaluate empirical relations between the threshold mass or the dominant postmerger gravitational-wave frequency and stellar parameters of nonrotating neutron stars. These relations are constructed to extract these stellar properties from merger observations and are built based on large sets of equation of state models which do not include pions. Comparing to our calculations with pions, we find that these empirical relations remain valid to good accuracy, which justifies their use although they neglect a possible impact of pions. We also address the mass ejection by neutron star mergers and observe a moderate enhancement of the ejecta mass by a few ten per cent. (abridged), Comment: 27 pages, 24 figures, accepted for publication in Phys. Rev. D

We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ a combination of hydrodynamics, nucleosynthesis, and radiative-transfer tools to achieve a consistent end-to-end modeling of the system and its observables. We adopt a novel version of the Shakura-Sunyaev scheme allowing to vary the approximate turbulent viscosity inside the NS remnant independently of the surrounding disk. We find that asymmetric progenitors lead to shorter remnant lifetimes and enhanced ejecta masses, although the viscosity affects the absolute values of these characteristics. The integrated production of lanthanides and heavier elements in such binary systems is sub-solar, suggesting that the considered scenarios contribute in a sub-dominant fashion to r-process enrichment. One reason is that BH-tori formed after delayed collapse exhibit less neutron-rich conditions than typically found, and often assumed in previous BH-torus models, for early BH formation. The outflows in our models feature strong anisotropy as a result of the lanthanide-poor polar neutrino-driven wind pushing aside lanthanide-rich dynamical ejecta. Considering the complexity of the models, the estimated kilonova light curves show promising agreement with AT2017gfo after times of several days, while the remaining inconsistencies at early times could possibly be overcome in binary configurations with a more dominant neutrino-driven wind relative to the dynamical ejecta., Comment: 16 pages, 9 figures, 1 table, accepted to ApJL

The mergers of neutron stars expel a heavy-element enriched fireball which can be observed as a kilonova. The kilonova's geometry is a key diagnostic of the merger and is dictated by the properties of ultra-dense matter and the energetics of the collapse to a black hole. Current hydrodynamical merger models typically show aspherical ejecta. Previously, Sr$^+$ was identified in the spectrum of the the only well-studied kilonova AT2017gfo, associated with the gravitational wave event GW170817. Here we combine the strong Sr$^+$ P Cygni absorption-emission spectral feature and the blackbody nature of kilonova spectrum, to determine that the kilonova is highly spherical at early epochs. Line shape analysis combined with the known inclination angle of the source also shows the same sphericity independently. We conclude that energy injection by radioactive decay is insufficient to make the ejecta spherical. A magnetar wind or jet from the black hole disk could inject enough energy to induce a more spherical distribution in the overall ejecta, however an additional process seems necessary to make the element distribution uniform

Neutron stars (NSs) and black holes (BHs) are born when the final collapse of the stellar core terminates the lives of stars more massive than about 9 Msun. This can trigger the powerful ejection of a large fraction of the star's material in a core-collapse supernova (CCSN), whose extreme luminosity is energized by the decay of radioactive isotopes such as 56Ni and 56Co. When evolving in close binary systems, the compact relics of such infernal catastrophes spiral towards each other on orbits gradually decaying by gravitational-wave emission. Ultimately, the violent collision of the two components forms a more massive, rapidly spinning remnant, again accompanied by the ejection of considerable amounts of matter. These merger events can be observed by high-energy bursts of gamma rays with afterglows and electromagnetic transients called kilonovae, which radiate the energy released in radioactive decays of freshly assembled rapid neutron-capture elements. By means of their mass ejection and the nuclear and neutrino reactions taking place in the ejecta, both CCSNe and compact object mergers (COMs) are prominent sites of heavy-element nucleosynthesis and play a central role in the cosmic cycle of matter and the chemical enrichment history of galaxies. The nuclear equation of state (EoS) of NS matter, from neutron-rich to proton-dominated conditions and with temperatures ranging from about zero to ~100 MeV, is a crucial ingredient in these astrophysical phenomena. It determines their dynamical processes, their remnant properties even at the level of deciding between NS or BH, and the properties of the associated emission of neutrinos, whose interactions govern the thermodynamic conditions and the neutron-to-proton ratio for nucleosynthesis reactions in the innermost ejecta. This chapter discusses corresponding EoS dependent effects of relevance in CCSNe as well as COMs. (slightly abridged), Comment: 96 pages, 23 figures; figure citations added; accepted for the Handbook of Nuclear Physics, eds. Isao Tanihata, Hiroshi Toki, Toshitaka Kajino

We present a relativistic density functional approach to color superconducting quark matter that mimics quark confinement by a fast growth of the quasiparticle selfenergy in the confining region. The approach is shown to be equivalent to a chiral model of quark matter with medium dependent couplings. While the (pseudo)scalar sector of the model is fitted to the vacuum phenomenology of quantum chromodynamics, the strength of interaction in the vector and diquark channels is varied in order to provide the best agreement with the observational constraints on the mass-radius relation and tidal deformability of neutron stars modelled with our approach. In order to recover the conformal behavior of quark matter at asymptotically high densities we introduce a medium dependence of the vector and diquark couplings motivated by the nonperturbative gluon exchange. Our analysis signals that the onset of deconfinement to color superconducting quark matter is likely to occur in neutron stars with masses below 1.0 $M_\odot$., Comment: 8 pages, 4 figuures, Contribution to Proceedings of Quark Confinement and the Hadron Spectrum XV, August 1-6, 2022, Stavanger, Norway

Binary neutron star (BNS) post-merger gravitational-wave emission can occur in the aftermath of a BNS merger -- provided the system avoids prompt collapse to a black hole -- as a quasistable hypermassive remnant experiences quadrupolar oscillations and non-axisymmetric deformations. The post-merger gravitational-wave spectrum possesses a characteristic peak frequency that has been shown to be dependent on the binary chirp mass and the neutron star equation of state (EoS), rendering post-merger gravitational waves a powerful tool for constraining neutron star composition. Unfortunately, the BNS post-merger signal is emitted at high ($\gtrsim 1.5$ kHz) frequencies, where ground-based gravitational wave detectors suffer from reduced sensitivity. It is therefore unlikely that post-merger signals will be detected with sufficient signal-to-noise ratio (SNR) until the advent of next-generation detectors. However, by employing empirical relations derived from numerical relativity simulations, we can combine information across an ensemble of BNS mergers, allowing us to obtain EoS constraints with many low-SNR signals. We present a hierarchical Bayesian method for deriving constraints on $R_{1.6}$, the radius of a 1.6$\mathrm{M_{\odot}}$ neutron star, through an ensemble analysis of binary neutron star mergers. We apply this method to simulations of the next two LIGO-Virgo-KAGRA observing runs, O4 and O5, as well as an extended 4-year run at A+ sensitivity, demonstrating the potential of our approach to yield EoS information from the post-merger signal with current-generation detectors. The A+ 4-year scenario is predicted to improve the constraint on $R_{1.6}$ from the currently available multimessenger-based 95\% credible interval (C.I.) uncertainty of $R_{1.6}=12.07^{+0.98}_{-0.77}$ km to $R_{1.6}=11.91^{+0.80}_{-0.56}$ km, a 22% reduction of the 95% C.I. width., Comment: 22 pages, 15 figures, 4 tables

The detection of GW170817 and the accompanying electromagnetic counterpart, AT2017gfo, have provided an important set of observational constraints for theoretical models of neutron star mergers, nucleosynthesis, and radiative transfer for kilonovae. We apply the 3D Monte Carlo radiative transfer code ARTIS to produce synthetic light curves of the dynamical ejecta from a neutron star merger, which has been modelled with 3D smooth-particle hydrodynamics (SPH) and included neutrino interactions. Nucleosynthesis calculations provide the energy released from radioactive decays of r-process nuclei, and radiation transport is performed using grey opacities given as functions of the electron fraction. We present line-of-sight dependent bolometric light curves, and find the emission along polar lines of sight to be up to a factor of ~2 brighter than along equatorial lines of sight. Instead of a distinct emission peak, our bolometric light curve exhibits a monotonic decline, characterised by a shoulder at the time when the bulk ejecta becomes optically thin. We show approximate band light curves based on radiation temperatures and compare these to the observations of AT2017gfo. We find that the rapidly declining temperatures lead to a blue to red colour evolution similar to that shown by AT2017gfo. We also investigate the impact of an additional, spherically symmetric secular ejecta component, and we find that the early light curve remains nearly unaffected, while after about 1 day the emission is strongly enhanced and dominated by the secular ejecta, leading to the shift of the shoulder from 1-2 to 6-10 days., Comment: Accepted by MNRAS

Based on a recently developed relativistic density functional approach to color-superconducting quark matter and a novel quark-hadron transition construction which phenomenologically accounts for the effects of inhomogeneous pasta phases and quark-hadron continuity, we construct a class of hybrid equations of state applicable at the regimes typical for compact star astrophysics and heavy ion collisions. We outline that early quark deconfinement is a notable consequence of strong diquark pairing providing a good agreement with the observational data and driving the trajectories of the matter evolution during the supernovae explosions toward the regimes typical for the compact star mergers and heavy-ion collisions., Comment: Presented at the $29^{\rm th}$ Conference ''Quark Matter 2022'' on ultrarelativistic nucleus-nucleus collisions, April 4-10, 2022, Krak\'ow, Poland; minor typo corrections in v2

We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated static neutron stars using a fixed metric or a dynamical spacetime. In both tests the frequencies of the radial oscillation mode match those of independent calculations. We run the first moving-mesh simulation of a neutron star merger. The simulation includes a scheme to adaptively refine or derefine cells and thereby adjusting the local resolution dynamically. The general dynamics are in agreement with independent smoothed particle hydrodynamics and static-mesh simulations of neutron star mergers. Coarsely comparing, we find that dynamical features like the post-merger double-core structure or the quasi-radial oscillation mode persist on longer time scales, possibly reflecting a low numerical diffusivity of our method. Similarly, the post-merger gravitational wave emission shows the same features as observed in simulations with other codes. In particular, the main frequency of the post-merger phase is found to be in good agreement with independent results for the same binary system, while, in comparison, the amplitude of the post-merger gravitational wave signal falls off slower, i.e. the post-merger oscillations are less damped. The successful implementation of general relativistic hydrodynamics in the moving-mesh AREPO code, including a dynamical spacetime evolution, provides a fundamentally new tool to simulate general relativistic problems in astrophysics., Comment: 24 pages, 18 figures, accepted for publication in MNRAS

Theoretically predicted yields of elements created by the rapid neutron capture (r-) process carry potentially large uncertainties associated with incomplete knowledge of nuclear properties and approximative hydrodynamical modelling of the matter ejection processes. We present an in-depth study of the nuclear uncertainties by varying theoretical nuclear input models that describe the experimentally unknown neutron-rich nuclei. This includes two frameworks for calculating the radiative neutron capture rates and 14 different models for nuclear masses, $\beta$-decay rates and fission properties. Our r-process nuclear network calculations are based on detailed hydrodynamical simulations of dynamically ejected material from NS-NS or NS-BH binary mergers plus the secular ejecta from BH-torus systems. The impact of nuclear uncertainties on the r-process abundance distribution and the early radioactive heating rate is found to be modest (within a factor of $\sim20$ for individual $A>90$ abundances and a factor of 2 for the heating rate). However, the impact on the late-time heating rate is more significant and depends strongly on the contribution from fission. We witness significantly larger sensitivity to the nuclear physics input if only a single trajectory is used compared to considering ensembles of $\sim$200-300 trajectories, and the quantitative effects of the nuclear uncertainties strongly depend on the adopted conditions for the individual trajectory. We use the predicted Th/U ratio to estimate the cosmochronometric age of six metal-poor stars and find the impact of the nuclear uncertainties to be up to 2 Gyr., Comment: 28 pages, 24 figures, accepted to MNRAS

Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities., Comment: 96 pages. Submitted to Journal of Physics G

We outline a theoretical approach supporting strong phase transitions from normal nuclear matter to the deconfined quark-gluon plasma, in the equation of state (EOS) for compact star matter. Implications of this hypothesis are discussed for astrophysical applications. Special emphasis is devoted to potentially detectable signatures, which can be directly related with the occurrence of a sufficiently strong phase transition. Therefore, simulations of core-collapse supernovae and binary compact star mergers are considered, including the subsequent emission of gravitational waves and, in the case of supernova, in addition the neutrinos play the role of messengers., Comment: 39 pages, 13 figures, contribution submitted on 20.02.2022 to the Book "Astrophysics in the XXI Century with Compact Stars", edited by Cesar Augusto Zen Vasconcellos and Fridolin Weber, World Scientific

The gravitational wave signal emitted during the coalescence of two neutron stars carries information about the stars' internal structure. During the long inspiral phase the main matter observable is the tidal interaction between the binary components, an effect that can be parametrically modeled with compact-binary solutions to General Relativity. After the binary merger the main observable is frequency modes of the remnant, most commonly giving rise to a short-duration signal accessible only through numerical simulations. The complicated morphology and the decreasing detector sensitivity in the relevant frequencies currently hinder detection of the post-merger signal and motivate separate analyses for the pre-merger and post-merger data. However, planned and ongoing detector improvements could soon put the post-merger signal within reach. In this study we target the whole pre-merger and post-merger signal without an artificial separation at the binary merger. We construct a hybrid analysis that models the inspiral with templates based on analytical calculations and calibrated to numerical relativity and the post-merger signal with a flexible morphology-independent analysis. Applying this analysis to GW170817 we find, as expected, that the post-merger signal remains undetected. We further study simulated signals and find that we can reconstruct the full signal and simultaneously estimate both the pre-merger tidal deformation and the post-merger signal frequency content. Our analysis allows us to study neutron star physics using all the data available and directly test the pre-merger and post-merger signal for consistency thus probing effects such as the onset of the hadron-quark phase transition., Comment: 18 pages, 17 figures, final published version

The kilonova (KN) associated with the binary neutron star (BNS) merger GW170817 is the only known electromagnetic counterpart to a gravitational wave source. Here we produce a sequence of radiative transfer models (using $\textsc{tardis}$) with updated atomic data, and compare them to accurately calibrated spectra. We use element compositions from nuclear network calculations based on a realistic hydrodynamical simulation of a BNS merger. We show that the blue spectrum at +1.4 days after merger requires a nucleosynthetic trajectory with a high electron fraction. Our best-fitting model is composed entirely of first $r$-process peak elements (Sr & Zr) and the strong absorption feature is reproduced well by Sr$\,\textsc{ii}$ absorption. At this epoch, we set an upper limit on the lanthanide mass fraction of $X_{\textsc{ln}} \lesssim 5 \times 10^{-3}$. In contrast, all subsequent spectra from $+2.4 - 6.4$ days require the presence of a modest amount of lanthanide material ($X_{\textsc{ln}} \simeq 0.05^{+0.05}_{-0.02}$), produced by a trajectory with $Y_{\rm e} = 0.29$. This produces lanthanide-induced line blanketing below 6000$\,$\AA, and sufficient light $r$-process elements to explain the persistent strong feature at $\sim 0.7 - 1.0 \, \micron$ (Sr$\,\textsc{ii}$). The composition gives good matches to the observed data, indicating that the strong blue flux deficit results in the near-infrared (NIR) excess. The disjoint in composition between the first epoch and all others indicates either ejecta stratification, or the presence of two distinct components of material. This further supports the `two-component' kilonova model, and constrains the element composition from nucleosynthetic trajectories. The major uncertainties lie in availability of atomic data and the ionisation state of the expanding material., Comment: 20 pages, 14 figures and 3 tables. Accepted for publication in MNRAS

We present a new analytic model describing gravitational wave emission in the post-merger phase of binary neutron star mergers. The model is described by a number of physical parameters that are related to various oscillation modes, quasi-linear combination tones or non-linear features that appear in the post-merger phase. The time evolution of the main post-merger frequency peak is taken into account and it is described by a two-segment linear expression. The effectiveness of the model, in terms of the fitting factor or, equivalently, the reduction in the detection rate, is evaluated along a sequence of equal-mass simulations of varying mass. We find that all parameters of the analytic model correlate with the total binary mass of the system. For high masses, we identify new spectral features originating from the non-linear coupling between the quasi-radial oscillation and the antipodal tidal deformation, the inclusion of which significantly improves the fitting factors achieved by the model. We can thus model the post-merger gravitational-wave emission with an analytic model that achieves high fitting factors for a wide range of total binary masses. Our model can be used for the detection and parameter estimation of the post-merger phase in upcoming searches with upgraded second-generation detectors, such as aLIGO+ and aVirgo+, with future, third-generation detectors (Einstein Telescope and Cosmic Explorer) or with dedicated, high-frequency detectors., Comment: 35 pages, 33 figures, accepted for publication in Phys. Rev. D

Astroparticle physics is undergoing a profound transformation, due to a series of extraordinary new results, such as the discovery of high-energy cosmic neutrinos with IceCube, the direct detection of gravitational waves with LIGO and Virgo, and many others. This white paper is the result of a collaborative effort that involved hundreds of theoretical astroparticle physicists and cosmologists, under the coordination of the European Consortium for Astroparticle Theory (EuCAPT). Addressed to the whole astroparticle physics community, it explores upcoming theoretical opportunities and challenges for our field of research, with particular emphasis on the possible synergies among different subfields, and the prospects for solving the most fundamental open questions with multi-messenger observations., Comment: White paper of the European Consortium for Astroparticle Theory (EuCAPT). 135 authors, 400 endorsers, 133 pages, 1382 references

The majority of existing results for the kilonova (or macronova) emission from material ejected during a neutron-star (NS) merger is based on (quasi-)one-zone models or manually constructed toy-model ejecta configurations. In this study we present a kilonova analysis of the material ejected during the first ~10ms of a NS merger, called dynamical ejecta, using directly the outflow trajectories from general relativistic smoothed-particle hydrodynamics simulations including a sophisticated neutrino treatment and the corresponding nucleosynthesis results, which have been presented in Part I of this study. We employ a multi-dimensional two-moment radiation transport scheme with approximate M1 closure to evolve the photon field and use a heuristic prescription for the opacities found by calibration with atomic-physics based reference results. We find that the photosphere is generically ellipsoidal but augmented with small-scale structure and produces emission that is about 1.5-3 times stronger towards the pole than the equator. The kilonova typically peaks after 0.7-1.5days in the near-infrared frequency regime with luminosities between 3-7x10^40erg/s and at photospheric temperatures of 2.2-2.8x10^3K. A softer equation of state or higher binary-mass asymmetry leads to a longer and brighter signal. Significant variations of the light curve are also obtained for models with artificially modified electron fractions, emphasizing the importance of a reliable neutrino-transport modeling. None of the models investigated here, which only consider dynamical ejecta, produces a transient as bright as AT2017gfo. The near-infrared peak of our models is incompatible with the early blue component of AT2017gfo., Comment: 23 pages, 15 figures, 1 table, accepted to MNRAS; updated for correct version of Fig. A1, left panel, and corrected Eqs. 7 and 9

We present a coherent study of the impact of neutrino interactions on the r-process element nucleosynthesis and the heating rate produced by the radioactive elements synthesised in the dynamical ejecta of neutron star-neutron star (NS-NS) mergers. We have studied the material ejected from four NS-NS merger systems based on hydrodynamical simulations which handle neutrino effects in an elaborate way by including neutrino equilibration with matter in optically thick regions and re-absorption in optically thin regions. We find that the neutron richness of the dynamical ejecta is significantly affected by the neutrinos emitted by the post-merger remnant, in particular when compared to a case neglecting all neutrino interactions. Our nucleosynthesis results show that a solar-like distribution of r- process elements with mass numbers $A \gtrsim 90$ is produced, including a significant enrichment in Sr and a reduced production of actinides compared to simulations without inclusion of the nucleonic weak processes. The composition of the dynamically ejected matter as well as the corresponding rate of radioactive decay heating are found to be rather independent of the system mass asymmetry and the adopted equation of state. This approximate degeneracy in abundance pattern and heating rates can be favourable for extracting the ejecta properties from kilonova observations, at least if the dynamical component dominates the overall ejecta. Part II of this work will study the light curve produced by the dynamical ejecta of our four NS merger models., Comment: 18 pages, 21 figures, 2 Tables, accepted by MNRAS 2021-11-19

The very first detection of gravitational waves from a neutron star binary merger, GW170817, exceeded all expectations. The event was relatively nearby, which may point to a relatively high merger rate. It was possible to extract finite-size effects from the gravitational-wave signal, which constrains the nuclear equation of state. Also, an electromagnetic counterpart was detected at many wavebands from radio to gamma rays marking the begin of a new multi-messenger era involving gravitational waves. We describe how multi-messenger observations of GW170817 are employed to constrain the nuclear equation of state. Combining the information from the optical emission and the mass measurement through gravitational waves leads to a lower limit on neutron star radii. According to this conservative analysis, which employs a minimum set of assumptions, the radii of neutron stars with typical masses should be larger than about 10.7~km. This implies a lower limit on the tidal deformability of about 210, while much stronger lower bounds are not supported by the data of GW170817. The multi-messenger interpretation of GW170817 rules out very soft nuclear matter and complements the upper bounds on NS radii which are derived from the measurement of finite-size effects during the pre-merger phase. We highlight the future potential of multi-messenger observations and of GW measurements of the postmerger phase for constraining the nuclear equation of state. Finally, we propose an observing strategy to maximize the scientific yield of future multi-messenger observations., Comment: 26 pages, 8 figures, contribution to "Special Issue on Nuclear astrophysics in the new era of multi-messenger astronomy", Annals of Physics

We relate the fundamental quadrupolar fluid mode of isolated non-rotating NSs and the dominant oscillation frequency of neutron star merger remnants. Both frequencies individually are known to correlate with certain stellar parameters like radii or the tidal deformability, which we further investigate by constructing fit formulae and quantifying the scatter of the data points from those relations. Furthermore, we compare how individual data points deviate from the corresponding fit to all data points. Considering this point-to-point scatter we uncover a striking similarity between the frequency deviations of perturbative data for isolated NSs and of oscillation frequencies of rapidly rotating, hot, massive merger remnants. The correspondence of frequency deviations in these very different stellar systems points to an underlying mechanism and EoS information being encoded in the frequency deviation. We trace the frequency scatter back to deviations of the tidal Love number from an average tidal Love number for a given stellar compactness. Our results thus indicate a possibility to break the degeneracy between NS radii, tidal deformability and tidal Love number. We also relate frequency deviations to the derivative of the tidal deformability with respect to mass. Our findings generally highlight a possibility to improve GW asteroseismology relations where the systematic behavior of frequency deviations is employed to reduce the scatter in such relationships and consequently increase the measurement accuracy. In addition, we relate the f-mode frequency of static stars and the dominant GW frequency of merger remnants. We find an analytic mapping to connect the masses of both stellar systems, which yields particularly accurate mass-independent relations between both frequencies and between the postmerger frequency and the tidal deformability., Comment: 25 pages, 16 figures, accepted for publication in Physical Review D

Black-hole (BH) accretion disks formed in compact-object mergers or collapsars may be major sites of the rapid-neutron-capture (r-)process, but the conditions determining the electron fraction (Y_e) remain uncertain given the complexity of neutrino transfer and angular-momentum transport. After discussing relevant weak-interaction regimes, we study the role of neutrino absorption for shaping Y_e using an extensive set of simulations performed with two-moment neutrino transport and again without neutrino absorption. We vary the torus mass, BH mass and spin, and examine the impact of rest-mass and weak-magnetism corrections in the neutrino rates. We also test the dependence on the angular-momentum transport treatment by comparing axisymmetric models using the standard alpha-viscosity with viscous models assuming constant viscous length scales (l_t) and three-dimensional magnetohydrodynamic (MHD) simulations. Finally, we discuss the nucleosynthesis yields and basic kilonova properties. We find that absorption pushes Y_e towards ~0.5 outside the torus, while inside increasing the equilibrium value Y_e^eq by ~0.05-0.2. Correspondingly, a substantial ejecta fraction is pushed above Y_e=0.25, leading to a reduced lanthanide fraction and a brighter, earlier, and bluer kilonova than without absorption. More compact tori with higher neutrino optical depth, tau, tend to have lower Y_e^eq up to tau~1-10, above which absorption becomes strong enough to reverse this trend. Disk ejecta are less (more) neutron-rich when employing an l_t=const. viscosity (MHD treatment). The solar-like abundance pattern found for our MHD model marginally supports collapsar disks as major r-process sites, although a strong r-process may be limited to phases of high mass-infall rates, Mdot>~ 2 x 10^(-2) Msun/s., Comment: 37 pages, 22 figures, 5 tables, accepted to MNRAS, minor changes compared to previous version

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. Here, we introduce new improvements to our scission-point model, called SPY, to derive the fission fragment distribution for all neutron-rich fissioning nuclei of relevance in r-process calculations. These improvements include a phenomenological modification of the scission distance and a smoothing procedure of the distribution. Such corrections lead to a much better agreement with experimental fission yields. Those yields are also used to estimate the number of neutrons emitted by the excited fragments on the basis of different neutron evaporation models. Our new fission yields are extensively compared to those predicted by the so-called GEF model. The impact of fission on the r-process nucleosynthesis in binary neutron mergers is also reanalyzed. Two scenarios are considered, the first one with low initial electron fraction is subject to intense fission recycling, in contrast to the second one which includes weak interactions on nucleons. The various regions of the nuclear chart responsible for fission recycling during the neutron irradiation as well as after freeze-out are discussed. The contribution fission processes may have to the final abundance distribution is also studied in detail in the light of newly defined quantitative indicators describing the fission recycling, the fission seeds and the fission progenitors. In particular, those allow us to estimate the contribution of fission to the final abundance distribution stemming from specific heavy nuclei. Calculations obtained with SPY and GEF fission fragment distributions are compared for both r-process scenarios.

We estimate the long-lasting gravitational wave (GW) emission of compact dark objects following a binary neutron-star (NS) merger. We consider compact dark objects, which initially reside in the centers of NSs and which may consist of self-interacting dark matter (DM). By approximating the compact dark objects as test particles, we model the merging of NS binaries hosting DM components with three-dimensional relativistic simulations. Our simulation results suggest that the DM components remain gravitationally bound and orbit inside the merger remnant with orbital separations of typically a few km. The subsequent orbital motion of the DM components generates a GW signal with frequencies in the range of a few kHz. When considering a range of different binary masses and high-density equations of state (EoS), we find that the GW frequency of the orbiting DM components scales with the compactness of NSs. Similarly, we find relations between the DM GW frequency and the dominant postmerger GW frequency of the stellar fluid or the tidal deformability, which quantifies EoS effects during the binary inspiral. Hence, a measurement of these quantities can be used to specify the frequency range of the GW emission by DM. Under the assumption that GW back reaction is the only relevant dissipative process, the GW signal may last between seconds and years depending on the mass of the DM component. We estimate the detectability of the GW signals and find that DM components in NS mergers may only be detectable with existing and projected GW instruments if the dark objects are as massive as about 0.01 to 0.1 M_sun. We emphasize that the GW emission is limited by the lifetime of the remnant. A forming black hole will immediately swallow the DM objects because their orbits are smaller than the innermost stable circular orbit of the black hole., Comment: 17 pages, 17 figures, accepted for publication in Phys. Rev. D

The first direct measurement of gravitational waves by the LIGO and Virgo collaborations has opened up new avenues to explore our Universe. This white paper outlines the challenges and gains expected in gravitational wave searches at frequencies above the LIGO/Virgo band, with a particular focus on Ultra High-Frequency Gravitational Waves (UHF-GWs), covering the MHz to GHz range. The absence of known astrophysical sources in this frequency range provides a unique opportunity to discover physics beyond the Standard Model operating both in the early and late Universe, and we highlight some of the most promising gravitational sources. We review several detector concepts which have been proposed to take up this challenge, and compare their expected sensitivity with the signal strength predicted in various models. This report is the summary of the workshop "Challenges and opportunities of high-frequency gravitational wave detection" held at ICTP Trieste, Italy in October 2019, that set up the stage for the recently launched Ultra-High-Frequency Gravitational Wave (UHF-GW) initiative., Comment: Published in Living Reviews in Relativity

This study addresses the collapse behavior of neutron star (NS) mergers expressed through the binary threshold mass M_thr for prompt black hole (BH) formation, which we determine by relativistic hydrodynamical simulations for 40 equation of state (EoS) models. M_thr can be well described by various fit formulae involving stellar parameters of nonrotating NSs. Using these relations we compute which constraints on NS radii and the tidal deformability are set by current and future merger detections revealing information about the merger product. We systematically investigate the impact of the binary mass ratio q=M_1/M_2 and assemble different fits, which make different assumptions about a-priori knowlegde. We find fit formulae for M_thr including an explicit q dependence, which are valid in a broad range of 0.7<=q<=1 and which are nearly as tight as relations for fixed mass ratios. For most EoS models except extreme cases M_thr of asymmetric mergers is equal or smaller than the one of equal-mass binaries. The impact of the binary mass asymmetry on M_thr becomes stronger with more extreme mass ratios, while M_thr is approximately constant for small deviations from q=1. We describe that a phase transition to deconfined quark matter can leave a characteristic imprint on the collapse behavior. The presence of quark matter can reduce the stability of the remnant and thus M_thr relative to a purely hadronic reference model. Comparing the threshold mass and the tidal deformability Lambda_thr of a system with M_thr can yield peculiar combinations of those two quantities, where M_thr is particularly small in relation to Lambda_thr. Hence, a combined measurement of both quantities can indicate the onset of quark deconfinement. We point out new univariate relations between M_thr and stellar properties of high-mass NSs, which can be employed for direct EoS constraints or consistency checks. (abridged), Comment: 35 pages, 13 figures, 14 tables, accepted for publication in Phys. Rev. D

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., Comment: 11 pages, 5 figures, submitted to EPJ Special Topic: Strong Correlations in Dense Matter Physics

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), Comment: 35 pages, 24 figures; accepted for publication in Phys. Rev. D

Finite-size effects on the gravitational wave signal from a neutron star merger typically manifest at high frequencies where detector sensitivity decreases. Proposed sensitivity improvements can give us access both to stronger signals and to a myriad of weak signals from cosmological distances. The latter will outnumber the former and the relevant part of signal will be redshifted towards the detector most sensitive band. We study the redshift dependence of information about neutron star matter and find that single-scale properties, such as the star radius or the post-merger frequency, are better measured from the distant weak sources from $z\sim 1$., Comment: 8 pages, 4 figures

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., Comment: Accepted for publication in Phys. Rev. Lett.; 14 pages, 5 figures

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., Comment: 26 pages, 16 figures

We study the partial fragmentation of a strange quark star into strangelets during the process of merger of two strange quark stars. We discuss the fate of the fragments considering their possible evaporation into nucleons. We show that only a rather small amount of large size strangelets, ejected from the spiral arms in the post-merger, survives a total evaporation into nucleons. In this way we demonstrate that: 1) the density of strangelets in the galaxy is too low to trigger the conversion of all neutron stars into strange quark stars and it allows the co-existence of both types of compact objects; 2) the probability of direct detection of a strangelet is negligible and therefore its nondetection is compatible with the strange quark matter hypothesis; 3) most of the matter ejected during and after the merger of two strange quark stars evaporates into nucleons and therefore it can generate a kilonova-like signal., Comment: 10 pages, 7 figures, 1 table. Significantly expanded version in print on Physical Review D

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., Comment: 7 pages, 3 figures