skip to main content


Title: Resolving the Fastest Ejecta from Binary Neutron Star Mergers: Implications for Electromagnetic Counterparts
Abstract We examine the effect of spatial resolution on initial mass ejection in grid-based hydrodynamic simulations of binary neutron star mergers. The subset of the dynamical ejecta with velocities greater than ∼0.6 c can generate an ultraviolet precursor to the kilonova on approximately hour timescales and contribute to a years long nonthermal afterglow. Previous work has found differing amounts of this fast ejecta, by one to two orders of magnitude, when using particle-based or grid-based hydrodynamic methods. Here, we carry out a numerical experiment that models the merger as an axisymmetric collision in a corotating frame, accounting for Newtonian self-gravity, inertial forces, and gravitational wave losses. The lower computational cost allows us to reach spatial resolutions as high as 4 m, or ∼3 × 10 −4 of the stellar radius. We find that fast ejecta production converges to within 10% for a cell size of 20 m. This suggests that fast ejecta quantities found in existing grid-based merger simulations are unlikely to increase to the level needed to match particle-based results upon further resolution increases. The resulting neutron-powered precursors are in principle detectable out to distances ≲200 Mpc with upcoming facilities.We also find that head-on collisions at the freefall speed, relevant for eccentric mergers, yield fast and slow ejecta quantities of order 10 −2 M ⊙ , with a kilonova signature distinct from that of quasi-circular mergers.  more » « less
Award ID(s):
2002577
NSF-PAR ID:
10329669
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
921
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
161
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We study the production of very light elements (Z< 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10−6Mof hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6Mof helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavyr-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.

     
    more » « less
  2. ABSTRACT

    Rapid identification of the optical counterparts of neutron star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and what kinds of photometric observations are needed to achieve this goal? NS merger kilonova exhibits low ejecta mass (∼5 × 10−2 M⊙) and a rapidly evolving photospheric radius (with a velocity ∼0.2c). As a consequence, these sources display rapid optical-flux evolution. Indeed, selection based on fast flux variations is commonly used for young supernovae and NS mergers. In this study, we leverage the best currently available flux-limited transient survey – the Zwicky Transient Facility Bright Transient Survey – to extend and quantify this approach. We focus on selecting transients detected in a 3-day cadence survey and observed at a one-day cadence. We explore their distribution in the phase space defined by g–r, $\dot{g}$, and $\dot{r}$. Our analysis demonstrates that for a significant portion of the time during the first week, the kilonova AT 2017gfo stands out in this phase space. It is important to note that this investigation is subject to various biases and challenges; nevertheless, it suggests that certain photometric observations can be leveraged to identify transients with the highest probability of being fast-evolving events. We also find that a large fraction (≈75 per cent) of the transient candidates with $\vert\dot{g}\vert>0.7$ mag d−1, are cataclysmic variables or active galactic nuclei with radio counterparts.

     
    more » « less
  3. Abstract

    Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resultingr-process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several$$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$1053ergs-1) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

     
    more » « less
  4. ABSTRACT We present a systematic numerical relativity study of the impact of different physics input and grid resolution in binary neutron star mergers. We compare simulations employing a neutrino leakage scheme, leakage plus M0 scheme, the M1 transport scheme, and pure hydrodynamics. Additionally, we examine the effect of a sub-grid scheme for turbulent viscosity. We find that the overall dynamics and thermodynamics of the remnant core are robust, implying that the maximum remnant density could be inferred from gravitational wave observations. Black hole collapse instead depends significantly on viscosity and grid resolution. Differently from recent work, we identify possible signatures of neutrino effects in the gravitational waves only at the highest resolutions considered; new high-resolution simulations will be thus required to build accurate gravitational wave templates to observe these effects. Different neutrino transport schemes impact significantly mass, geometry, and composition of the remnant’s disc and ejecta; M1 simulations show systematically larger proton fractions, reaching maximum values larger than 0.4. r-process nucleosynthesis yields reflect the different ejecta compositions; they are in agreement and reproduce residual solar abundances only if M0 or M1 neutrino transport schemes are adopted. We compute kilonova light curves using spherically-symmetric radiation-hydrodynamics evolutions up to 15 d post-merger, finding that they are mostly sensitive to the ejecta mass and electron fraction; accounting for multiple ejecta components appears necessary for reliable light curve predictions. We conclude that advanced neutrino schemes and resolutions higher than current standards are essential for robust long-term evolutions and detailed astrophysical predictions. 
    more » « less
  5. Abstract

    The Wide-Field Infrared Transient Explorer (WINTER) is a new 1 deg2seeing-limited time-domain survey instrument designed for dedicated near-infrared follow-up of kilonovae from binary neutron star (BNS) and neutron star–black hole mergers. WINTER will observe in the near-infraredY,J, and short-Hbands (0.9–1.7μm, toJAB= 21 mag) on a dedicated 1 m telescope at Palomar Observatory. To date, most prompt kilonova follow-up has been in optical wavelengths; however, near-infrared emission fades more slowly and depends less on geometry and viewing angle than optical emission. We present an end-to-end simulation of a follow-up campaign during the fourth observing run (O4) of the LIGO, Virgo, and KAGRA interferometers, including simulating 625 BNS mergers, their detection in gravitational waves, low-latency and full parameter estimation skymaps, and a suite of kilonova lightcurves from two different model grids. We predict up to five new kilonovae independently discovered by WINTER during O4, given a realistic BNS merger rate. Using a larger grid of kilonova parameters, we find that kilonova emission is ≈2 times longer lived and red kilonovae are detected ≈1.5 times further in the infrared than in the optical. For 90% localization areas smaller than 150 (450) deg2, WINTER will be sensitive to more than 10% of the kilonova model grid out to 350 (200) Mpc. We develop a generalized toolkit to create an optimal BNS follow-up strategy with any electromagnetic telescope and present WINTER’s observing strategy with this framework. This toolkit, all simulated gravitational-wave events, and skymaps are made available for use by the community.

     
    more » « less