Dynamical Ejection Research Articles

Numerical relativity simulations of the neutron star merger GW190425: microphysics and mass ratio effects

ABSTRACT GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

Kinematics of the Central Stars Powering Bowshock Nebulae and the Large Multiplicity Fraction of Runaway OB Stars

OB stars powering stellar bowshock nebulae (SBNe) have been presumed to have large peculiar velocities. We measured peculiar velocities of SBN central stars to assess their kinematics relative to the general O-star population using Gaia EDR3 data for 267 SBN central stars and a sample of 455 Galactic O stars to derive projected velocities v 2D. For a subset of each sample, we obtained new optical spectroscopy to measure radial velocities and identify multiple-star systems. We find a minimum multiplicity fraction of 36% ± 6% among SBN central stars, consistent with >28% among runaway Galactic O stars. The large multiplicity fraction among runaways implicates very efficient dynamical ejection rather than binary-supernova origins. The median v 2D of SBN central stars is v 2D = 14.6 km s−1, larger than the median v 2D = 11.4 km s−1 for non-bowshock O stars. Central stars of SBNe have a runaway (v 2D > 25 km s−1) fraction of 24%, consistent with the % for control-sample O stars. Most (76%) SBNe central stars are not runaways. Our analysis of alignment (ΔPA) between the nebular morphological and v 2D kinematic position angles reveals two populations: a highly aligned (σ PA = 25°) population that includes stars with the largest v 2D (31% of the sample) and a random (nonaligned) population (69%). SBNe that lie within or near H ii regions comprise a larger fraction of this latter component than SBNe in isolated environments, implicating localized ISM flows as a factor shaping their orientations and morphologies. We outline a new conceptual approach to computing the solar local standard of rest motion, yielding [U ⊙, V ⊙, W ⊙] = [5.5, 7.5,4.5] km s−1.

The viability of the optically thick wind model in accreting white dwarfs at high accretion rates

ABSTRACT The single degenerate (SD) model is one of the leading models for the progenitors of Type Ia supernovae (SNe Ia). The optically thick wind (OTW) model is considered to be essential for the SD model. In this article, we aim to constrain the condition for the occurrence of the OTW. Using the stellar evolution code, Modules for Experiments in Stellar Astrophysics, we carried out a series of one-dimensional hydrodynamical simulations with different accretion rates and different white dwarf (WD) masses to simulate the dynamical ejection process on the accreting WDs. We find that, by neglecting interaction with the accreted material, we can get a similar result to the previous works using steady-state method, i.e. the OTW occurs. However, if we consider the interaction between the accreted and the outflowing materials, whether the OTW can occur or not becomes strongly dependent on the accretion rate. For a given WD mass, the higher the accretion rate, the less likely the OTW can occur. In addition, the WD masses can also affect the occurrence of the OTW, i.e. the higher the WD mass, the easier the OTW can occur. As a result, the OTW in novae may occur because the accretion rate is low enough, while in the SD progenitor system of SNe Ia, a high accretion rate is likely to prevent the occurrence of the OTW.

General-relativistic neutrino-radiation magnetohydrodynamic simulation of seconds-long black hole-neutron star mergers

Seconds-long numerical-relativity simulations for black hole-neutron star mergers are performed for the first time to obtain a self-consistent picture of the merger and post-merger evolution processes. To investigate the case that tidal disruption takes place, we choose the initial mass of the black hole to be $5.4M_\odot$ or $8.1M_\odot$ with the dimensionless spin of 0.75. The neutron-star mass is fixed to be $1.35M_\odot$. We find that after the tidal disruption, dynamical mass ejection takes place spending $\lesssim 10$ ms together with the formation of a massive accretion disk. Subsequently, the magnetic field in the disk is amplified by the magnetic winding and magnetorotational instability, establishing a turbulent state and inducing the angular momentum transport. The post-merger mass ejection by the magnetically-induced viscous effect sets in at $\sim 300$-500\,ms after the tidal disruption, at which the neutrino luminosity drops below $\sim 10^{51.5}\,{\rm erg/s}$, and continues for several hundreds ms. A magnetosphere near the rotational axis of the black hole is developed after the matter and magnetic flux fall into the black hole from the accretion disk, and high-intensity Poynting flux generation sets in at a few hundreds ms after the tidal disruption. The intensity of the Poynting flux becomes low after the significant post-merger mass ejection, because the opening angle of the magnetosphere increases. The lifetime for the stage with the strong Poynting flux is $1$-2 s, which agrees with the typical duration of short-hard gamma-ray bursts.

The origin of early-type runaway stars from open clusters

Runaway stars are ejected from their place of birth in the Galactic disk, with some young B-type runaways found several tens of kiloparsecs from the plane traveling at speeds beyond the escape velocity, which calls for violent ejection processes. Young open clusters are a likely place of origin, and ejection may be either through N-body interactions or in binary supernova explosions. The most energetic events may require dynamical interaction with massive black holes. The excellent quality of Gaia astrometry opens up the path to study the kinematics of young runaway stars to such a high precision that the place of origin in open stellar clusters can be identified uniquely even when the star is a few kiloparsecs away. We developed an efficient minimization method to calculate whether two or more objects may come from the same place, which we tested against samples of Orion runaways. Our fitting procedure was then used to calculate trajectories for known runaway stars, which have previously been studied from HIPPARCOS astrometry as well as known open clusters. For runaways in our sample we used Gaia data and updated radial velocities, and found that only half of the sample could be classified as runaways. The other half of the sample moves so slowly (< 30 km s−1) that they have to be considered as walkaway stars. Most of the latter stars turned out to be binaries. We identified parent clusters for runaways based on their trajectories. We then used cluster age and flight time of the stars to investigate whether the ejection was likely due to a binary supernova or due to a dynamical ejection. In particular we show that the classical runaways AE Aurigae and μ Columbae might not have originated together, with μ Columbae having an earlier ejection from Collinder 69, a cluster near the ONC. The second sample investigated comprises a set of distant runaway B stars in the halo which have been studied carefully by quantitative spectral analyses. We are able to identify candidate parent clusters for at least four stars including the hyper-runaway candidate HIP 60350. The ejection events had to be very violent, ejecting stars at velocities as large as 150–400 km s−1.

On the Jet–Ejecta Interaction in 3D GRMHD Simulations of a Binary Neutron Star Merger Aftermath

Short γ-ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta’s isotropic equivalent mass along the pole at the time of the BH formation is ≲10−4 M ⊙, setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳1015 G; and (iii) if the jets are weak (≲1050 erg), the ejecta’s isotropic equivalent mass along the pole must be small (≲10−2 M ⊙). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ∼1.5).

Thinking Outside the Box: Numerical Relativity with Particles

The observation of gravitational waves from compact objects has now become an active part of observational astronomy. For a sound interpretation, one needs to compare such observations against detailed Numerical Relativity simulations, which are essential tools to explore the dynamics and physics of compact binary mergers. To date, essentially all simulation codes that solve the full set of Einstein’s equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations, we focus on perfectly symmetric binary systems that are irrotational and have 1.3 M⊙ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle–mesh mapping), and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the Ψ4 formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times 10−3M⊙ as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity (∼0.5…0.7c) ejecta component of ∼10−4M⊙ that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission, and it could also potentially cause a kilonova afterglow.

Dynamical Evolution of Ejecta from the DART Impact on Dimorphos

The DART spacecraft will impact Dimorphos (the secondary body of the Didymos binary asteroid) to test the kinetic impactor deflection method against possibly hazardous near-Earth asteroids. The DART impact ejecta plume, and possibly the impact crater, will be imaged by the LICIACube spacecraft, hosted as a piggyback and released by DART just before the impact, and then, several years later, by the Hera probe. To exploit the wealth of data obtained and understand the physics of the whole impact experiment, it is of paramount importance to properly model the dynamics of the binary system pre- and postimpact and the dynamics of the particles ejected from the impact crater. A model was developed to simulate the evolution of the ejecta particles created during the impact in order to first interpret the LICIACube images and then test the survival of particles on long intervals of time that might be detected by the Hera mission either as individual bodies or as parts of rings. The dynamical evolution of the particles is simulated over different timescales to highlight the most important perturbations and their relative importance. The ejecta dynamics turns out to be highly chaotic due to repeated close encounters with the two asteroids. However, we find that some ejecta survive in the binary orbital environment for timescales comparable to the Hera arrival time. The effects of the particles reimpacting against either one of the components is also analyzed to estimate the amount of momentum transfer to the target bodies.

Impact of extreme spins and mass ratios on the post-merger observables of high-mass binary neutron stars

ABSTRACTThe gravitational-wave events GW170817 and GW190425 have led to a number of important insights on the equation of state of dense matter and the properties of neutron stars, such as their radii and the maximum mass. Some of these conclusions have been drawn on the basis of numerical-relativity simulations of binary neutron-star mergers with vanishing initial spins. While this may be a reasonable assumption in equal-mass systems, it may be violated in the presence of large mass asymmetries accompanied by the presence of high spins. To quantify the impact of high spins on multimessenger gravitational-wave events, we have carried out a series of high-mass binary neutron-star mergers with a highly spinning primary star and large mass asymmetries that have been modelled self-consistently using two temperature-dependent equations of state. We show that, when compared with equal-mass, irrotational binaries, these systems can lead to significant differences in the remnant lifetime, in the dynamical ejecta, in the remnant disc masses, in the secular ejecta, and on the bulk kilonova properties. These differences could be exploited to remove the degeneracy between low- and high-spin priors in the detection of gravitational waves from binary neutron-star mergers.

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACTWe present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

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−6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r-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.

Jet launching from merging magnetized binary neutron stars with realistic equations of state

We perform general relativistic, magnetohydrodynamic (GRMHD) simulations of binary neutron stars in quasi-circular orbit that merge and undergo delayed or prompt collapse to a black hole (BH). The stars are irrotational and modeled using an SLy or an H4 nuclear equation of state. To assess the impact of the initial magnetic field configuration on jet launching, we endow the stars with a purely poloidal magnetic field that is initially unimportant dynamically and is either confined to the stellar interior or extends from the interior into the exterior as in typical pulsars. Consistent with our previous results, we find that only the BH + disk remnants originating from binaries that form hypermassive neutron stars (HMNSs) and undergo delayed collapse can drive magnetically-powered jets. We find that the closer the total mass of the binary is to the threshold value for prompt collapse, the shorter is the time delay between the gravitational wave peak amplitude and jet launching. This time delay also strongly depends on the initial magnetic field configuration. We also find that seed magnetic fields confined to the stellar interior can launch a jet over $\sim 25\,\rm ms$ later than those with pulsar-like magnetic fields. The lifetime of the jet [$\Delta t\lesssim 150\,\rm ms$] and its outgoing Poynting luminosity [$L_{\rm EM}\sim 10^{52\pm 1}\rm erg/s$] are consistent with typical short gamma-ray burst central engine lifetimes, as well as with the Blandford--Znajek mechanism for launching jets and their associated Poynting luminosities. Our numerical results also suggest that the dynamical ejection of matter can be enhanced by the magnetic field. Therefore, GRMHD studies are required to fully understand kilonova signals from GW170818-like events.

Mapping dynamical ejecta and disk masses from numerical relativity simulations of neutron star mergers

We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions, that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear EOS lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.

Impact of the nuclear symmetry energy on the post-merger phase of a binary neutron star coalescence

The nuclear symmetry energy plays a key role in determining the equation of state of dense, neutron-rich matter, which governs the properties of both terrestrial nuclear matter as well as astrophysical neutron stars. A recent measurement of the neutron skin thickness from the PREX Collaboration has lead to new constraints on the slope of the nuclear symmetry energy, $L$, which can be directly compared to inferences from gravitational wave observations of the first binary neutron star merger inspiral, GW170817. In this paper, we explore a new regime for potentially constraining the slope, $L$, of the nuclear symmetry energy with future gravitational wave events: the post-merger phase of a binary neutron star coalescence. In particular, we go beyond the inspiral phase, where imprints of the slope parameter $L$ may be inferred from measurements of the tidal deformability, to consider imprints on the post-merger dynamics, gravitational wave emission, and dynamical mass ejection. To this end, we perform a set of targeted neutron star merger simulations in full general relativity using new finite-temperature equations of state, which systematically vary $L$ while keeping the magnitude of the symmetry energy at the saturation density, $S$, fixed. We find that the post-merger dynamics and gravitational wave emission are mostly insensitive to the slope of the nuclear symmetry energy. In contrast, we find that dynamical mass ejection contains a weak imprint of $L$, with large values of $L$ leading to systematically enhanced ejecta.

Dynamical ejecta of neutron star mergers with nucleonic weak processes I: nucleosynthesis

ABSTRACT 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 synthesized 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.

Dynamical ejecta of neutron star mergers with nucleonic weak processes – II: kilonova emission

ABSTRACT 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 $\sim 10\,$ ms 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 multidimensional 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.5\,$ d in the near-infrared frequency regime with luminosities between $3\!-\!7\times 10^{40}\,$ erg s−1 and at photospheric temperatures of $2.2\!-\!2.8\times 10^3\,$ K. 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 modelling. 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.

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.6c 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.

Signatures of deconfined quark phases in binary neutron star mergers

We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We treat hadronic matter using a Brueckner-Hartree-Fock quantum many-body approach and modern two-body and three-body nuclear interactions derived within chiral effective field theory. Quark matter is modeled using an extended version of the bag model. We combine these approaches to construct a new finite-temperature composition-dependent equation of state (EOS) with a first-order phase transition between hadrons and deconfined quarks. We perform numerical relativity simulations of BNS mergers with this new EOS and compare results obtained with or without the deconfinment phase transition. We find that deconfined quark production in a neutron star merger results from matter crossing the phase boundary over a wide range of temperatures and densities. The softening of the EOS due to the phase transition causes the merger remnants to be more compact and to collapse to a black hole at earlier times. The phase transition is imprinted on the postmerger gravitational wave (GW) signal duration, amplitude, and peak frequency. However, this imprint is only detectable for binaries with sufficiently long-lived remnants. Moreover, the phase transition does not result in significant deviations from quasiuniversal relations for the postmerger GW peak frequency. Consequently, the postmerger GW peak frequency alone is not sufficient to conclusively exclude or confirm the presence of a phase transition in a BNS merger. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, $r$-process nucleosynthetic yields and associated electromagnetic counterparts. While there are differences in the electromagnetic counterparts and nucleosynthesis yields between the purely hadronic models and the models with phase transitions, these can be primarily ascribed to the difference in remnant collapse time between the two, so they are degenerate with other effects. An exception is the nonthermal afterglow caused by the interaction of the fastest component of the dynamical ejecta and the interstellar medium, which is systematically boosted in the binaries with phase transition as a consequence of the more violent merger they experience.

A Global Numerical Model of the Prompt Emission in Short Gamma-ray Bursts

Abstract We present the first global model of prompt emission from a short gamma-ray burst (GRB) that consistently describes the evolution of the central black hole (BH) torus system, the propagation of the jet through multicomponent merger ejecta, the transition into free expansion, and the photospheric emission from the relativistic jet. To this end, we perform a special relativistic neutrino-hydrodynamics simulation of a viscous BH-torus system, which is formed about 500 ms after the merger and is surrounded by dynamical ejecta as well as neutron star winds, along with a jet that is injected in the vicinity of the central BH. In a postprocessing step, we compute the photospheric emission using a relativistic Monte Carlo radiative transfer code. It is found that the wind from the torus leaves a strong imprint on the jet as well as on the emission, causing narrow collimation and rapid time variability. The dependence of the emission on viewing angle gives rise to correlations among the spectral peak energy, E p , isotropic energy, E iso, and peak luminosity, L p , which may provide natural explanations for the Amati and Yonetoku relations. We also find that the degree of polarization is small for emission from the jet core (≲2%), while it tends to increase with viewing angle outside the core and can become as high as ∼10%–40% for energies larger than the peak energy. Finally, the comparison of our model with GRB 170817A strongly disfavors the photospheric emission scenario and therefore supports alternative scenarios, such as cocoon shock breakout.

Gamma-Ray Emission Produced by r-process Elements from Neutron Star Mergers

Abstract The observation of a radioactively powered kilonova AT 2017gfo associated with the gravitational wave event GW170817 from a binary neutron star merger proves that these events are ideal sites for the production of heavy r-process elements. The gamma-ray photons produced by the radioactive decay of heavy elements are unique probes for the detailed nuclide compositions. Based on the detailed r-process nucleosynthesis calculations and considering radiative transport calculations for the gamma rays in different shells, we study the gamma-ray emission in a merger ejecta on a timescale of a few days. It is found that the total gamma-ray energy generation rate evolution is roughly depicted as E ̇ ∝ t − 1.3 . For the dynamical ejecta with a low electron fraction (Y e ≲ 0.20), the dominant contributors of gamma-ray energy are the nuclides around the second r-process peak (A ∼ 130) and the decay chain of 132Te (t 1/2 = 3.21 days) → 132I (t 1/2 = 0.10 days) → 132Xe produces gamma-ray lines at 228, 668, and 773 keV. For the case of a wind ejecta with Y e ≳ 0.30, the dominant contributors of gamma-ray energy are the nuclides around the first r-process peak (A ∼ 80) and the decay chain of 72Zn (t 1/2 = 1.93 days) → 72Ga (t 1/2 = 0.59 days) → 72Ge produces gamma-ray lines at 145, 834, 2202, and 2508 keV. The peak fluxes of these lines are 10−9 ∼ 10−7 ph cm−2 s−1, which are marginally detectable with the next-generation MeV gamma-ray detector ETCC if the source is at a distance of 40 Mpc.