No Arabic abstract
Neutrino-matter interactions play a key role in binary neutron star mergers. Thermodynamics conditions at the surfaces where neutrinos decouple from matter influence neutrino spectra, ultimately affecting the evolution of the remnant and the properties of the ejecta. In this work, we post-process results of general relativistic merger simulations employing microphysical equations of state and approximate neutrino transport to investigate the thermodynamics conditions at which weak and thermal equilibrium freezes out (equilibrium surfaces), as well as conditions at which the transition between diffusion and free-streaming regime occurs (diffusion surfaces). We find that the rest mass density and the neutrino energy are the most relevant quantities in determining the location of the decoupling surfaces. For mean energy neutrinos ($langle E_{{ u}_e} rangle approx 9~{rm MeV}$, $langle E_{{bar{ u}}_e} rangle approx 15~{rm MeV}$, $langle E_{{ u}_{mu,tau}} rangle approx 25~{rm MeV}$), diffusion surfaces are located around $10^{11}{rm g~cm^{-3}}$ for all neutrino species, while equilibrium surfaces for heavy flavor neutrinos are significantly deeper (several $10^{12}{rm g~cm^{-3}}$) than the ones of $bar{ u}_e$ and $ u_e$ ($gtrsim 10^{11}{rm g~cm^{-3}}$). The resulting decoupling temperatures are in good agreement with the average neutrino energies ($langle E_{ u} rangle sim 3.15~T$), with the softer equation of state characterized by systematically larger decoupling temperatures ($Delta T lesssim 1~{rm MeV}$). Neutrinos streaming at infinity with different energies come from very different regions of the remnant. The presence of a massive NS or of a BH in the remnant influences the neutrino thermalization process.
In order to extract maximal information from neutron-star merger signals, both gravitational and electromagnetic, we need to ensure that our theoretical models/numerical simulations faithfully represent the extreme physics involved. This involves a range of issues, with the finite temperature effects regulating many of the relevant phenomena. As a step towards understanding these issues, we explore the conditions for $beta$-equilibrium in neutron star matter for the densities and temperatures reached in a binary neutron star merger. Using the results from our out-of-equilibrium merger simulation, we consider how different notions of equilibrium may affect the merger dynamics, raising issues that arise when attempting to account for these conditions in future simulations. These issues are both computational and conceptual. We show that the effects lead to, in our case, a softening of the equation of state in some density regions, and to composition changes that affect processes that rely on deviation from equilibrium, such as bulk viscosity, both in terms of the magnitude and the equilibration timescales inherent to the relevant set of reactions. We also demonstrate that it is difficult to determine exactly which equilibrium conditions are relevant in which regions of the matter due to the dependence on neutrino absorption, further complicating the calculation of the reactions that work to restore the matter to equilibrium.
The detections of gravitational waves (GWs) from binary neutron star (BNS) systems and neutron star--black hole (NSBH) systems provide new insights into dense matter properties in extreme conditions and associated high-energy astrophysical processes. However, currently information about NS equation of state (EoS) is extracted with very limited precision. Meanwhile, the fruitful results from the serendipitous discovery of the $gamma$-ray burst alongside GW170817 show the necessity of early warning alerts. Accurate measurements of the matter effects and sky location could be achieved by joint GW detection from space and ground. In our work, based on two example cases, GW170817 and GW200105, we use the Fisher information matrix analysis to investigate the multiband synergy between the space-borne decihertz GW detectors and the ground-based Einstein Telescope (ET). We specially focus on the parameters pertaining to spin-induced quadrupole moment, tidal deformability, and sky localization. We demonstrate that, (i) only with the help of multiband observations can we constrain the quadrupole parameter; and (ii) with the inclusion of decihertz GW detectors, the errors of tidal deformability would be a few times smaller, indicating that many more EoSs could be excluded; (iii) with the inclusion of ET, the sky localization improves by about an order of magnitude. Furthermore, we have systematically compared the different limits from four planned decihertz detectors and adopting two widely used waveform models.
With the first observation of a binary neutron star merger through gravitational waves and light GW170817, compact binary mergers have now taken the center stage in nuclear astrophysics. They are thought to be one of the main astrophysical sites of production of r-process elements, and merger observations have become a fundamental tool to constrain the properties of matter. Here, we review our current understanding of the dynamics of neutron star mergers, in general, and of GW170817 in particular. We discuss the physical processes governing the inspiral, merger, and postmerger evolution, and we highlight the connections between these processes, the dynamics, and the multimessenger observables. Finally, we discuss open questions and issues in the field and the need to address them through a combination of better theoretical models and new observations.
We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_odot-10M_odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperature nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($chi_{rm BH}gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($Tlesssim 1,{rm MeV}$), unbound, neutron-rich material ($M_{rm ej}sim 0.05M_odot-0.20M_odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{rm max}sim 15,{rm MeV}$) with typical neutrino luminosity $L_ usim 10^{53},{rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $sim 10,{rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{rm BH}sim 7M_odot-10M_odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.
(abridged) We investigate the quark deconfinement phase transition in the context of binary neutron star (BNS) mergers. We employ a new finite-temperature composition-dependent equation of state (EOS) with a first order phase transition between hadrons and deconfined quarks to perform numerical relativity simulations of BNS mergers. 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 (BH) 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 quasi-universal relations for the postmerger GW peak frequency. We also study the impact of the phase transition on dynamical ejecta, remnant accretion disk masses, r-process nucleosynthetic yields and associated electromagnetic (EM) counterparts. While there are differences in the EM 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. An exception is the non-thermal 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.