An official website of the United States government
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−6M⊙of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6M⊙of 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
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
- PAR ID:
- 10329669
- 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
-
-
null (Ed.)ABSTRACT Neutron star mergers produce a substantial amount of fast-moving ejecta, expanding outwardly for years after the merger. The interaction of these ejecta with the surrounding medium may produce a weak isotropic radio remnant, detectable in relatively nearby events. We use late-time radio observations of short duration gamma-ray bursts (sGRBs) to constrain this model. Two samples of events were studied: four sGRBs that are possibly in the local (<200 Mpc) Universe were selected to constrain the remnant non-thermal emission from the sub-relativistic ejecta, whereas 17 sGRBs at cosmological distances were used to constrain the presence of a proto-magnetar central engine, possibly re-energizing the merger ejecta. We consider the case of GRB 170817A/GW170817 and find that in this case the early radio emission may be quenched by the jet blast-wave. In all cases, for ejecta mass range of $${M}_{\rm {ej}}\lesssim 10^{-2}\, (5\times 10^{-2})\, \mathrm{M}_\odot$$, we can rule out very energetic merger ejecta $${E}_{\rm {ej}}\gtrsim 5\times 10^{52}\, (10^{53})\, \rm erg$$, thus excluding the presence of a powerful magnetar as a merger remnant.more » « less
-
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
-
Kilonovae, one source of electromagnetic emission associated with neutron star mergers, are powered by the decay of radioactive isotopes in the neutron-rich merger ejecta. Models for kilonova emission consistent with the electromagnetic counterpart to GW170817 predict characteristic abundance patterns, determined by the relative balance of different types of material in the outflow. Assuming that the observed source is prototypical, this inferred abundance pattern in turn must matchr-process abundances deduced by other means, such as what is observed in the solar system. We report on analysis comparing the input mass-weighted elemental compositions adopted in our radiative transfer simulations to the mass fractions of elements in the Sun, as a practical prototype for the potentially universal abundance signature from neutron star mergers. We characterize the extent to which our parameter inference results depend on our assumed composition for the dynamical and wind ejecta and examine how the new results compare to previous work. We find that a dynamical ejecta composition calculated using the FRDM2012 nuclear mass and FRLDM fission models with extremely neutron-rich ejecta (Ye= 0.035) along with moderately neutron-rich (Ye= 0.27) wind ejecta composition yields a wind-to-dynamical mass ratio ofMw/Md= 0.47, which best matches the observed AT2017gfo kilonova light curves while also producing the best-matching abundance of neutron capture elements in the solar system, though, allowing for systematics, the ratio may be as high as of order unity.more » « less
-
null (Ed.)Binary neutron star mergers provide a unique probe of the dense-matter equation of state (EoS) across a wide range of parameter space, from the zero-temperature EoS during the inspiral to the high-temperature EoS following the merger. In this paper, we implement a new model for calculating parametrized finite-temperature EoS effects into numerical relativity simulations. This "M* model" is based on a two-parameter approximation of the particle effective mass and includes the leading-order effects of degeneracy in the thermal pressure and energy. We test our numerical implementation by performing evolutions of rotating single stars with zero- and non-zero temperature gradients, as well as evolutions of binary neutron star mergers. We find that our new finite-temperature EoS implementation can support stable stars over many dynamical timescales. We also perform a first parameter study to explore the role of the M* parameters in binary neutron star merger simulations. All simulations start from identical initial data with identical cold EoSs, and differ only in the thermal part of the EoS. We find that both the thermal profile of the remnant and the post-merger gravitational wave signal depend on the choice of M* parameters, but that the total merger ejecta depends only weakly on the finite-temperature part of the EoS across a wide range of parameters. Our simulations provide a first step toward understanding how the finite-temperature properties of dense matter may affect future observations of binary neutron star mergers.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/ac1f20
