Simulations Of Binary Neutron Star Mergers Research Articles

What Can We Learn about the Unstable Equation-of-state Branch from Neutron Star Mergers?

The equation of state (EOS) of dense strongly interacting matter can be probed by astrophysical observations of neutron stars (NS), such as X-ray detections of pulsars or the measurement of the tidal deformability of NSs during the inspiral stage of NS mergers. These observations constrain the EOS at most up to the density of the maximum-mass configuration, n TOV, which is the highest density that can be explored by stable NSs for a given EOS. However, under the right circumstances, binary neutron star (BNS) mergers can create a postmerger remnant that explores densities above n TOV. In this work, we explore whether the EOS above n TOV can be measured from gravitational-wave or electromagnetic observations of the postmerger remnant. We perform a total of 25 numerical-relativity simulations of BNS mergers for a range of EOSs and find no case in which different descriptions of the matter above n TOV have a detectable impact on postmerger observables. Hence, we conclude that the EOS above n TOV can likely not be probed through BNS merger observations for the current and next generation of detectors.

Searching for Gravitational-wave Counterparts Using the Transiting Exoplanet Survey Satellite

In 2017, the LIGO and Virgo gravitational-wave (GW) detectors, in conjunction with electromagnetic (EM) astronomers, observed the first GW multimessenger astrophysical event, the binary neutron star (BNS) merger GW170817. This marked the beginning of a new era in multimessenger astrophysics. To discover further GW multimessenger events, we explore the synergies between the Transiting Exoplanet Survey Satellite (TESS) and GW observations triggered by the LIGO–Virgo–KAGRA Collaboration (LVK) detector network. TESS's extremely wide field of view (∼2300 deg2) means that it could overlap with large swaths of GW localizations, which often span hundreds of square degrees or more. In this work, we use a recently developed transient detection pipeline to search TESS data collected during the LVK’s third observing run, O3, for any EM counterparts. We find no obvious counterparts brighter than about 17th magnitude in the TESS bandpass. Additionally, we present end-to-end simulations of BNS mergers, including their detection in GWs and simulations of light curves, to identify TESS's kilonova discovery potential for the LVK's next observing run (O4). In the most optimistic case, TESS will observe up to one GW-found BNS merger counterpart per year. However, TESS may also find up to five kilonovae that did not trigger the LVK network, emphasizing that EM-triggered GW searches may play a key role in future kilonova detections. We also discuss how TESS can help place limits on EM emission from binary black hole mergers and rapidly exclude large sky areas for poorly localized GW events.

Muons in the aftermath of neutron star mergers and their impact on trapped neutrinos

Context. In the upcoming years, present and next-generation gravitational wave observatories will detect a larger number of binary neutron star (BNS) mergers with increasing accuracy. In this context, improving BNS merger numerical simulations is crucial to correctly interpret the data and constrain the equation of state (EOS) of neutron stars (NSs). Aims. State-of-the-art simulations of BNS mergers do not include muons. However, muons are known to be relevant in the microphysics of cold NSs and are expected to have a significant role in mergers, where the typical thermodynamic conditions favour their production. Our work is aimed at investigating the impact of muons on the merger remnant. Methods. We post-process the outcome of four numerical relativity simulations of BNS mergers performed with three different baryonic EOSs and two mass ratios considering the first 15 milliseconds after merger. We compute the abundance of muons in the remnant and analyse how muons affect the trapped neutrino component and the fluid pressure. Results. We find that depending on the baryonic EOS, the net fraction of muons is between 30% and 70% the net fraction of electrons. Muons change the flavour hierarchy of trapped (anti-)neutrinos such that deep inside the remnant, muon anti-neutrinos are the most abundant, followed by electron anti-neutrinos. Finally, muons and trapped neutrinos modify the neutron-to-proton ratio, affecting the remnant pressure by up to 7% when compared with calculations neglecting them. Conclusions. This work demonstrates that muons have a non-negligible effect on the outcome of BNS merger simulations, and they should be included to improve the accuracy of a simulation.

3D radiative transfer kilonova modelling for binary neutron star merger simulations

ABSTRACTThe 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 three-dimensional (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 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 that along equatorial lines of sight. Instead of a distinct emission peak, our bolometric light curve exhibits a monotonic decline, characterized 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\,$ d the emission is strongly enhanced and dominated by the secular ejecta, leading to the shift of the shoulder from ∼1–2 to 6–10 d.

GRMHD Simulations of Neutron-star Mergers with Weak Interactions: r-process Nucleosynthesis and Electromagnetic Signatures of Dynamical Ejecta

Fast neutron-rich material ejected dynamically over ≲10 ms during the merger of a binary neutron star (BNS) can give rise to distinctive electromagnetic counterparts to the system’s gravitational-wave emission that serve as a “smoking gun” to distinguish between a BNS and an NS–black hole merger. We present novel ab initio modeling of the kilonova precursor and kilonova afterglow based on 3D general-relativistic magnetohydrodynamic simulations of BNS mergers with nuclear, tabulated, finite-temperature equations of state (EOSs), weak interactions, and approximate neutrino transport. We analyze dynamical mass ejection from 1.35–1.35 M ⊙ binaries, consistent with properties of the first observed BNS merger GW170817, using three nuclear EOSs that span the range of allowed compactness of 1.35 M ⊙-neutron stars. Nuclear reaction network calculations yield a robust second-to-third-peak r-process. We find few ×10−6 M ⊙ of fast (v > 0.6c) ejecta that give rise to broadband synchrotron emission on ∼years timescales, consistent with tentative evidence for excess X-ray/radio emission following GW170817. We find ≈2 × 10−5 M ⊙ of free neutrons that power a kilonova precursor on ≲ hours timescale. A boost in early UV/optical brightness by a factor of a few due to previously neglected relativistic effects, with enhancements up to ≲10 hr post-merger, is promising for future detection with UV/optical telescopes like Swift or ULTRASAT. We find that a recently predicted opacity boost due to highly ionized lanthanides at ≳70,000 K is unlikely to affect the early kilonova based on the obtained ejecta structures. Azimuthal inhomogeneities in dynamical ejecta composition for soft EOSs found here (“lanthanide/actinide pockets”) may have observable consequences for both early kilonova and late-time nebular emission.

Modelling kilonova afterglows: Effects of the thermal electron population and interaction with GRB outflows

AbstractGiven an increasing number of gamma-ray bursts accompanied by potential kilonovae, there is a growing importance to advance modelling of kilonova afterglows. In this work, we investigate how the presence of two electron populations that follow a Maxwellian (thermal) and a power-law (non-thermal) distribution affect kilonova afterglow light curves. We employ semi-analytic afterglow model, PyBlastAfterglow. We consider kilonova ejecta profiles from ab-initio numerical relativity binary neutron star merger simulations, targeted to GW170817. We do not perform model selection. We find that the emission from thermal electrons dominates at early times. If the interstellar medium density is high (${\simeq }0.1\, \, \text{cm}^{-3}$), it adds an early time peak to the light curve. As ejecta decelerates, the spectral and temporal indexes change in a characteristic way that, if observed, can be used to reconstruct the ejecta velocity distribution. For the low interstellar medium density, inferred for GRB 170817A, the emission from the non-thermal electron population generally dominates. We also assess how kilonova afterglow light curves change if the interstellar medium has been partially removed and pre-accelerated by laterally expanding gamma-ray burst ejecta. For the latter, we consider properties informed by observations of GRB170817A. We find that the main effect is the emission suppression at early time ${\lesssim }10^{3}\,$ days, and at its maximum it reaches ${\sim }40{{\ \rm per\ cent}}$ when the fast tail of the kilonova ejecta moves subsonically through the wake of laterally spreading gamma-ray burst ejecta. The subsequent rebrightening, when these ejecta break through and shocks form, is very mild (${\lesssim }10{{\ \rm per\ cent}}$) and may not be observable.

Binary neutron star merger simulations with neutrino transport and turbulent viscosity: impact of different schemes and grid resolution

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.

Binary neutron star mergers as a probe of quark-hadron crossover equations of state

It is anticipated that the gravitational radiation detected in future gravitational wave (GW) detectors from binary neutron star (NS) mergers can probe the high-density equation of state (EOS). We perform the first simulations of binary NS mergers which adopt various parametrizations of the quark-hadron crossover (QHC) EOS. These are constructed from combinations of a hadronic EOS ($n_{b} ~5~n_0$), where $n_{b}$ and $n_0$ are the baryon number density and the nuclear saturation density, respectively. At the crossover densities ($2~ n_0 < n_{b} < 5~ n_0$) the QHC EOSs continuously soften, while remaining stiffer than hadronic and first-order phase transition EOSs, achieving the stiffness of strongly correlated quark matter. This enhanced stiffness leads to significantly longer lifetimes of the postmerger NS than that for a pure hadronic EOS. We find a dual nature of these EOSs such that their maximum chirp GW frequencies $f_{max}$ fall into the category of a soft EOS while the dominant peak frequencies ($f_{peak}$) of the postmerger stage fall in between that of a soft and stiff hadronic EOS. An observation of this kind of dual nature in the characteristic GW frequencies will provide crucial evidence for the existence of strongly interacting quark matter at the crossover densities for QCD.

Equation of state of neutron star matter and its warm extension with an interacting hadron resonance gas

We propose an interpolating equation of state that satisfies phenomenologically established boundary conditions in two extreme regimes at high temperature and low baryon density and at low temperature and high baryon density. We confirm that the hadron resonance gas model with the Carnahan-Starling excluded volume effect can reasonably fit the empirical equation of state at high density up to several times the normal nuclear density. We identify the onsets of strange particles and quantify the strangeness contents in dense matter. We finally discuss the finite temperature effects and estimate the thermal index Γth as a function of the baryon density, which should be a crucial input for the core-collapse supernova and the binary neutron star merger simulations.

Finite-temperature effects in dynamical spacetime binary neutron star merger simulations: validation of the parametric approach

ABSTRACT Parametric equations of state (EoSs) provide an important tool for systematically studying EoS effects in neutron star merger simulations. In this work, we perform a numerical validation of the M*-framework for parametrically calculating finite-temperature EoS tables. The framework, introduced by Raithel et al., provides a model for generically extending any cold, β-equilibrium EoS to finite temperatures and arbitrary electron fractions. In this work, we perform numerical evolutions of a binary neutron star merger with the SFHo finite-temperature EoS, as well as with the M*-approximation of this same EoS, where the approximation uses the zero-temperature, β-equilibrium slice of SFHo and replaces the finite-temperature and composition-dependent parts with the M*-model. We find that the approximate version of the EoS is able to accurately recreate the temperature and thermal pressure profiles of the binary neutron star remnant, when compared to the results found using the full version of SFHo. We additionally find that the merger dynamics and gravitational wave signals agree well between both cases, with differences of $\lesssim 1\!-\!2\,{\textrm{per cent}}$ introduced into the post-merger gravitational wave peak frequencies by the approximations of the EoS. We conclude that the M*-framework can be reliably used to probe neutron star merger properties in numerical simulations.

Improving the convergence order of binary neutron star merger simulations in the Baumgarte- Shapiro-Shibata-Nakamura formulation

High-accuracy numerical relativity simulations of binary neutron star mergers are a necessary ingredient for constructing gravitational waveform templates to analyze and interpret observations of compact object mergers. Numerical convergence in the post-merger phase of such simulations is challenging to achieve with many modern codes. In this paper, we study two ways of improving the convergence properties of binary neutron star merger simulations within the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's equations. We show that discontinuities in a particular constraint damping scheme in this formulation can destroy the post-merger convergence of the simulation. A continuous prescription, in contrast, ensures convergence until late times. We additionally study the impact of the equation of state parametrization on the pre- and post-merger convergence properties of the simulations. In particular, we compare results for a piecewise polytropic parametrization, which is commonly used in merger simulations but suffers unphysical discontinuities in the sound speed, with results using a "generalized" piecewise polytropic parametrization, which was designed to ensure both continuity and differentiability of the equation of state. We report on the differences in the gravitational waves and any spurious pre-merger heating, depending on which equation of state parametrization is used.

Neutrino emission from binary neutron star mergers: characterising light curves and mean energies

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 resulting r-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}}) 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.

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.

Universality of the Turbulent Magnetic Field in Hypermassive Neutron Stars Produced by Binary Mergers

Abstract The detection of a binary neutron star (BNS) merger in 2017 through both gravitational waves and electromagnetic emission opened a new era of multimessenger astronomy. The understanding of the magnetic field amplification triggered by the Kelvin–Helmholtz instability during the merger is still a numerically unresolved problem because of the relevant small scales involved. One of the uncertainties comes from the simplifications usually assumed in the initial magnetic topology of merging neutron stars. We perform high-resolution, convergent large-eddy simulations of BNS mergers, following the newly formed remnant for up to 30 ms. Here we specifically focus on the comparison between simulations with different initial magnetic configurations, going beyond the widespread-used aligned dipole confined within each star. The results obtained show that the initial topology is quickly forgotten, in a timescale of a few milliseconds after the merger. Moreover, at the end of the simulations, the average intensity (B ∼ 1016 G) and the spectral distribution of magnetic energy over spatial scales barely depend on the initial configuration. This is expected due to the small-scale efficient dynamo involved, and thus it holds as long as (i) the initial large-scale magnetic field is not unrealistically high (as often imposed in mergers studies), and (ii) the turbulent instability is numerically (at least partially) resolved, so that the amplified magnetic energy is distributed across a wide range of scales and becomes orders of magnitude larger than the initial one.

Binary neutron star mergers of quark matter based nuclear equations of state

With observations of gravitational wave signals from binary neutron star mergers (BNSM) by LIGO-Virgo-KAGRA (LVK) Collaboration and NICER, the nuclear equation of state (EOS) is becoming increasingly testable by complementary numerical simulations. Numerous simulations currently explore the EOS at different density regimes for the constituent neutron stars specifically narrowing the uncertainty in the sub-nuclear densities. In this paper we summarize the three-dimensional general relativistic-hydrodynamics based simulations of BNSMs for EOSs with a specific emphasis on the quark matter EOS at the highest densities.

Models of binary neutron star remnants with tabulated equations of state

ABSTRACT The emergence of novel differential rotation laws that can reproduce the rotational profile of binary neutron star merger remnants has opened the way for the construction of equilibrium models with properties that resemble those of remnants in numerical simulations. We construct models of merger remnants, using a recently introduced 4-parameter differential rotation law and three tabulated, zero-temperature equations of state. The models have angular momenta that are determined by empirical relations, constructed through numerical simulations. After a systematic exploration of the parameter space of merger remnant equilibrium sequences, which includes the determination of turning points along constant angular momentum sequences, we find that a particular rotation law can reproduce the threshold mass to prompt collapse to a black hole with a relative difference of only $\sim 1{{\ \rm per\ cent}}$ with respect to numerical simulations, in all cases considered. Furthermore, our results indicate a possible correlation between the compactness of equilibrium models of remnants at the threshold mass and the compactness of maximum-mass non-rotating models. Another key prediction of binary neutron star merger simulations is a relatively slowly rotating inner region, where the angular velocity Ω (as measured by an observer at infinity) is mostly due to the frame dragging angular velocity ω. In our investigation of the parameter space of the adopted differential rotation law, we naturally find quasi-spherical (Type A) remnant models with this property. Our investigation clarifies the impact of the differential rotation law and of the equation of state on key properties of binary neutron star remnants and lays the groundwork for including thermal effects in future studies.

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.

Realistic finite-temperature effects in neutron star merger simulations

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.

Dynamical ejecta synchrotron emission as a possible contributor to the changing behaviour of GRB170817A afterglow

ABSTRACT Over the past 3 yr, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equations of state (EOSs) and high-mass ratio models with stiff EOSs are disfavoured as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.

Binary neutron star merger simulations with hot microscopic equations of state

We perform binary neutron star merger simulations using a newly derived set of finite-temperature equations of state in the Brueckner-Hartree-Fock approach. We point out the important and opposite roles of finite temperature and rotation for stellar stability and systematically investigate the gravitational-wave properties, matter distribution, and ejecta properties in the postmerger phase for the different cases. The validity of several universal relations is also examined and the most suitable EOSs are identified.