Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85 million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017  s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

A Modern View of the Equation of State in Nuclear and Neutron Star Matter

Background: We analyze several constraints on the nuclear equation of state (EOS) currently available from neutron star (NS) observations and laboratory experiments and study the existence of possible correlations among properties of nuclear matter at saturation density with NS observables. Methods: We use a set of different models that include several phenomenological EOSs based on Skyrme and relativistic mean field models as well as microscopic calculations based on different many-body approaches, i.e., the (Dirac–)Brueckner–Hartree–Fock theories, Quantum Monte Carlo techniques, and the variational method. Results: We find that almost all the models considered are compatible with the laboratory constraints of the nuclear matter properties as well as with the largest NS mass observed up to now, 2.14−0.09+0.10M⊙ for the object PSR J0740+6620, and with the upper limit of the maximum mass of about 2.3–2.5M⊙ deduced from the analysis of the GW170817 NS merger event. Conclusion: Our study shows that whereas no correlation exists between the tidal deformability and the value of the nuclear symmetry energy at saturation for any value of the NS mass, very weak correlations seem to exist with the derivative of the nuclear symmetry energy and with the nuclear incompressibility.

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

Multimessenger constraints on the neutron-star equation of state and the Hubble constant

Observations of neutron-star mergers with distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus and to measure the expansion rate of the Universe as quantified by the Hubble constant. We performed a joint analysis of the gravitational-wave event GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave event GW190425, both originating from neutron-star mergers. We combined these with previous measurements of pulsars using x-ray and radio observations, and nuclear-theory computations using chiral effective field theory, to constrain the neutron-star equation of state. We found that the radius of a 1.4-solar mass neutron star is [Formula: see text] km at 90% confidence and the Hubble constant is [Formula: see text] at 1σ uncertainty.

Neutron star merger remnants

Binary neutron star mergers observations are a unique way to constrain fundamental physics and astrophysics at the extreme. The interpretation of gravitational-wave events and their electromagnetic counterparts crucially relies on general-relativistic models of the merger remnants. Quantitative models can be obtained only by means of numerical relativity simulations in 3 + 1 dimensions including detailed input physics for the nuclear matter, electromagnetic and weak interactions. This review summarizes the current understanding of merger remnants focusing on some of the aspects that are relevant for multimessenger observations.

Equation of State Constraints from the Threshold Binary Mass for Prompt Collapse of Neutron Star Mergers

Using hydrodynamical simulations for a large set of high-density matter equations of state (EOSs), we systematically determine the threshold mass M_{thres} for prompt black-hole formation in equal-mass and asymmetric neutron star (NS) mergers. We devise the so far most direct, general, and accurate method to determine the unknown maximum mass of nonrotating NSs from merger observations revealing M_{thres}. Considering hybrid EOSs with hadron-quark phase transition, we identify a new, observable signature of quark matter in NS mergers. Furthermore, our findings have direct applications in gravitational wave searches, kilonova interpretations, and multimessenger constraints on NS properties.

The Equation of State of Nuclear Matter: From Finite Nuclei to Neutron Stars

Background. We investigate possible correlations between neutron star observables and properties of atomic nuclei. In particular, we explore how the tidal deformability of a 1.4 solar mass neutron star, M1.4, and the neutron-skin thickness of 48Ca and 208Pb are related to the stellar radius and the stiffness of the symmetry energy. Methods. We examine a large set of nuclear equations of state based on phenomenological models (Skyrme, NLWM, DDM) and ab initio theoretical methods (BBG, Dirac–Brueckner, Variational, Quantum Monte Carlo). Results: We find strong correlations between tidal deformability and NS radius, whereas a weaker correlation does exist with the stiffness of the symmetry energy. Regarding the neutron-skin thickness, weak correlations appear both with the stiffness of the symmetry energy, and the radius of a M1.4. Our results show that whereas the considered EoS are compatible with the largest masses observed up to now, only five microscopic models and four Skyrme forces are simultaneously compatible with the present constraints on L and the PREX experimental data on the 208Pb neutron-skin thickness. We find that all the NLWM and DDM models and the majority of the Skyrme forces are excluded by these two experimental constraints, and that the analysis of the data collected by the NICER mission excludes most of the NLWM considered. Conclusion. The tidal deformability of a M1.4 and the neutron-skin thickness of atomic nuclei show some degree of correlation with nuclear and astrophysical observables, which however depends on the ensemble of adopted EoS.

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multi-messenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative, impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

Multimessenger Bayesian parameter inference of a binary neutron star merger

The combined detection of a binary neutron star merger in both gravitational waves (GWs) and electromagnetic (EM) radiation spanning the entire spectrum – GW170817/AT2017gfo/GRB170817A – marks a breakthrough in the field of multimessenger astronomy. Between the plethora of modelling and observations, the rich synergy that exists among the available data sets creates a unique opportunity to constrain the binary parameters, the equation of state of supranuclear density matter, and the physical processes at work during the kilonova and gamma-ray burst. We report, for the first time, Bayesian parameter estimation combining information from GW170817, AT2017gfo, and GRB170817 to obtain truly multimessenger constraints on the tidal deformability $\tilde{\Lambda } \in [302,860]$, total binary mass M ∈ [2.722, 2.751] M⊙, the radius of a 1.4 solar mass neutron star $R \in [11.3,13.5] \,\,\rm km$ (with additional $0.2\ \rm km$ systematic uncertainty), and an upper bound on the mass ratio of q ≤ 1.27, all at 90 per cent confidence. Our joint novel analysis uses new phenomenological descriptions of the dynamical ejecta, debris disc mass, and remnant black hole properties, all derived from a large suite of numerical relativity simulations.