No Arabic abstract
We study the long-term evolution of ejecta formed in a binary neutron star (BNS) merger that results in a long-lived remnant NS by performing a hydrodynamics simulation with the outflow data of a numerical relativity simulation as the initial condition. At the homologously expanding phase, the total ejecta mass reaches $approx0.1,M_odot$ with an average velocity of $approx0.1,c$ and lanthanide fraction of $approx 0.005$. We further perform the radiative transfer simulation employing the obtained ejecta profile. We find that, contrary to a naive expectation from the large ejecta mass and low lanthanide fraction, the optical emission is not as bright as that in GW170817/AT2017gfo, while the infrared emission can be brighter. This light curve property is attributed to preferential diffusion of photons toward the equatorial direction due to the prolate ejecta morphology, large opacity contribution of Zr, Y, and lanthanides, and low specific heating rate of the ejecta. Our results suggest that these light curve features could be used as an indicator for the presence of a long-lived remnant NS. We also found that the bright optical emission broadly consistent with GW170817/AT2017gfo is realized for the case that the high-velocity ejecta components in the polar region are suppressed. These results suggest that the remnant in GW170817/AT2017gfo is unlikely to be a long-lived NS, but might have collapsed to a black hole within ${cal O}(0.1)$ s.
Two neutron stars merge somewhere in the Universe approximately every 10 seconds, creating violent explosions observable in gravitational waves and across the electromagnetic spectrum. The transformative coincident gravitational-wave and electromagnetic observations of the binary neutron star merger GW170817 gave invaluable insights into these cataclysmic collisions, probing bulk nuclear matter at supranuclear densities, the jet structure of gamma-ray bursts, the speed of gravity, and the cosmological evolution of the local Universe, among other things. Despite the wealth of information, it is still unclear when the remnant of GW170817 collapsed to form a black hole. Evidence from other short gamma-ray bursts indicates a large fraction of mergers may form long-lived neutron stars. We review what is known observationally and theoretically about binary neutron star post-merger remnants. From a theoretical perspective, we review our understanding of the evolution of short- and long-lived merger remnants, including fluid, magnetic-field, and temperature evolution. These considerations impact prospects of detection of gravitational waves from either short- or long-lived neutron star remnants which potentially allows for new probes into the hot nuclear equation of state in conditions inaccessible in terrestrial experiments. We also review prospects for determining post-merger physics from current and future electromagnetic observations, including kilonovae and late-time x-ray and radio afterglow observations.
We present twelve new simulations of unequal mass neutron star mergers. The simulations were preformed with the SpEC code, and utilize nuclear-theory based equations of state and a two-moment gray neutrino transport scheme with an improved energy estimate based on evolving the number density. We model the neutron stars with the SFHo, LS220 and DD2 equations of state (EOS) and we study the neutrino and matter emission of all twelve models to search for robust trends between binary parameters and emission characteristics. We find that the total mass of the dynamical ejecta exceeds $0.01M_odot$ only for SFHo with weak dependence on the mass-ratio across all models. We find that the ejecta have a broad electron fraction ($Y_e$) distribution ($approx 0.06-0.48$), with mean $0.2$. $Y_e$ increases with neutrino irradiation over time, but decreases with increasing binary asymmetry. We also find that the models have ejecta with a broad asymptotic velocity distribution ($approx 0.05-0.7c$). The average velocity lies in the range $0.2c - 0.3c$ and decreases with binary asymmetry. Furthermore, we find that disk mass increases with binary asymmetry and stiffness of the EOS. The $Y_e$ of the disk increases with softness of the EOS. The strongest neutrino emission occurs for the models with soft EOS. For (anti) electron neutrinos we find no significant dependence of the magnitude or angular distribution or neutrino luminosity with mass-ratio. The heavier neutrino species have a luminosity dependence on mass-ratio but an angular distribution which does not change with mass-ratio.
We construct a set of binary evolutionary sequences for systems composed by a normal, solar composition, donor star together with a neutron star. We consider a variety of masses for each star as well as for the initial orbital period corresponding to systems that evolve to ultra-compact or millisecond pulsar-helium white dwarf pairs. Specifically, we select a set of donor star masses of 0.50, 0.65, 0.80, 1.00, 1.25, 1.50, 1.75, 2.00, 2.25, 2.50, 3.00, and 3.50 solar masses, whereas for the accreting neutron star we consider initial masses values of 0.8, 1.0, 1.2, and 1.4 solar masses. The considered initial orbital period interval ranges from 0.5 to 12 days. It is found that the evolution of systems, with fixed initial values for the orbital period and the mass of the normal donor star, heavily depends upon the mass of the neutron star. In some cases, varying the initial value of the neutron star mass, we obtain evolved configurations ranging from ultra-compact to widely separated objects. We also analyse the dependence of the final orbital period with the mass of the white dwarf. In agreement with previous expectations, our calculations show that the final orbital period-white dwarf mass relation is fairly insensitive to the initial neutron star mass value. A new period-mass relation based on our own calculations is proposed, which is in good agreement with period-mass relations available in the literature. As consequence of considering a set of values for the initial neutron star mass, these models allow finding different plausible initial configurations (donor and neutron star masses and orbital period interval) for some of the best observed binary systems of the kind we are interested in here. We apply our calculations to analyse the case of PSR J0437-4715.
Weak reactions are critical for the neutron richness of the matter dynamically ejected after the merger of two neutron stars. The neutron richness, defined by the electron fraction (Ye), determines which heavy elements are produced by the r-process and thus directly impacts the kilonova light curve. In this work, we have performed a systematic and detailed post-processing study of the impact of weak reactions on the distribution of the electron fraction and of the entropy on the dynamic ejecta obtained from an equal mass neutron star binary merger simulated in full general relativity and with microscopic equation of state. Previous investigations indicated that shocks increase Ye, however our results show that shocks can also decrease Ye, depending on their thermodynamical conditions. Moreover, we have found that neutrino absorption are key and need to be considered in future simulations. We also demonstrated that the angular dependence of the neutrino luminosity and the spatial distribution of the ejecta can lead to significant difference in the electron fraction distribution. In addition to the detailed study of the Ye evolution and its dependences, we have performed nucleosynthesis calculations. They clearly point to the necessity of improving the neutrino treatment in current simulations to be able to predict the contribution of neutron star mergers to the chemical history of the universe and to reliable calculate their kilonova light curves.
With the epochal first detection of gravitational waves from a binary neutron star (NS) merger with the GW170817 event, and its direct confirmation that NS-NS mergers are significant sources of the of the r-process nucleosynthesis of heavy elements, an immense new arena for prompt EM (X-rays through IR and radio) studies of fundamental physics has been opened. Over the next decade, GW observatories will expand in scale and sensitivity so the need for facilities that can provide prompt, high sensitivity, broad-band EM followup becomes more urgent. NS-NS or NS-black hole (BH) mergers will be instantly recognized (and announced) by the LIGO-international collaboration. LSST will be a prime resource for rapid tiling of what will usually be large (~10-100 degree squared) error boxes. X-ray through IR Telescopes in space with (nearly) full-sky access that can rapidly image and tile are crucial for providing the earliest imaging and spectroscopic studies of the kilonova emission immediately following NS-NS mergers. The Time-domain Spectroscopic Observatory (TSO) is a proposed Probe-class 1.3 m telescope at L2, with imaging and spectroscopy (R = 200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky. TSO nUV-mid-IR spectra will enable new constraints on NS structure and nucleosynthesis.