We present new numerical relativity results of neutron star mergers with chirp mass $1.188M_\odot$ and mass ratios $q=1.67$ and $q=1.8$ using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model derived from ab-initio calculations based on the Brueckner-Hartree-Fock approach. We report for the first time evidence for accretion-induced prompt collapse in high-mass-ratio mergers, in which the tidal disruption of the companion and its accretion onto the primary star determine prompt black hole formation. As a result of the tidal disruption, an accretion disc of neutron-rich and cold matter forms with baryon masses ${\sim}0.15M_\odot$, and it is significantly heavier than the remnant discs in equal-masses prompt collapse mergers. Massive dynamical ejecta of order ${\sim}0.01M_\odot$ also originate from the tidal disruption. They are neutron rich and expand from the orbital plane with a crescent-like geometry. Consequently, bright, red and temporally extended kilonova emission is predicted from these mergers. Our results show that prompt black hole mergers can power bright electromagnetic counterparts for high-mass-ratio binaries, and that the binary mass ratio can be in principle constrained from multimessenger observations.

The AT2017gfo kilonova counterpart of the binary neutron star merger event GW170817 was characterized by an early-time bright peak in optical and UV bands. Such blue kilonova is commonly interpreted as a signature of weak $r$-process nucleosynthesis in a fast expanding wind whose origin is currently debated. Numerical-relativity simulations with microphysical equations of state, approximate neutrino transport, and turbulent viscosity reveal a new mechanism that can power the blue kilonova. Spiral density waves in the remnant generate a characteristic wind of mass ${\sim}10^{-2}~M_{\odot}$ and velocity ${\sim}0.2$c. The ejected material has electron fraction mostly distributed above $0.25$ being partially reprocessed by hydrodynamic shocks in the expanding arms. The combination of dynamical ejecta and spiral-wave wind can account for solar system abundances of $r$-process elements and early-time observed light curves.

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_{{\nu}_e} \rangle \approx 9~{\rm MeV}$, $\langle E_{{\bar{\nu}}_e} \rangle \approx 15~{\rm MeV}$, $\langle E_{{\nu}_{\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{\nu}_e$ and $\nu_e$ ($\gtrsim 10^{11}{\rm g~cm^{-3}}$). The resulting decoupling temperatures are in good agreement with the average neutrino energies ($\langle E_{\nu} \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.

We present fully general relativistic simulations of binary neutron star mergers, employing a new zero- temperature chiral effective field theory equation of state, the BL EOS. We offer a comparison with respect to the older GM3 EOS, which is based on standard relativistic mean field theory, and separately determine the impact of the mass. We provide a detailed analysis of the dynamics, with focus on the post-merger phase. For all models, we extract the gravitational wave strain and the post-merger frequency spectrum. Further, we determine the amount, velocity, and polar distribution of ejected matter, and provide estimates for the resulting kilonova signals. We also study the evolution of the disk while it is interacting with the hypermassive remnant, and dis- cuss the merits of different disk mass definitions applicable before collapse, with regard to the mass remaining after BH formation. Finally, we investigate the radial mass distribution and rotation profile of the remnants, which validate previous results and also corroborate a recently proposed stability criterion.

We present general relativistic numerical simulations of binary neutron star (BNS) mergers with different initial spin configurations. We focus on models with stars of mass 1.4 M_sol each, which employ the equation of state (EOS) by Shen, Horowitz, and Teige, and which result in stable NSs as merger remnants. For comparison, we consider two irrotational equal mass (M=1.35 M_sol) and unequal mass (M=1.29,1.42 M_sol) BNS models using the APR4 EOS, which result in a supramassive merger remnant. We present visualizations of the fluid flow and temperature distribution and find a strong impact of the spin on vortex structure and nonaxisymmetric deformation. We compute the radial mass distribution and the rotation profile in the equatorial plane using recently developed measures independent of spatial gauge, revealing slowly rotating cores that can be well approximated by the cores of spherical stars. We also study the influence of the spin on the inspiral phase and the gravitational wave (GW) signal. Using a newly developed analysis method, we further show that gravitational waveforms from BNS mergers can exhibit one or more phase jumps after merger, which occur together with minima of the strain amplitude. We provide a natural explanation in terms of the remnant's quadrupole moment, and show that cancellation effects due to phase jumps can have a strong impact on the GW power spectrum. Finally, we discuss the impact of the spin on the amount of ejected matter.

Merging binary neutron stars (BNSs) represent the ultimate targets for multimessenger astronomy, being among the most promising sources of gravitational waves (GWs), and, at the same time, likely accompanied by a variety of electromagnetic counterparts across the entire spectrum, possibly including short gamma-ray bursts (SGRBs) and kilonova/macronova transients. Numerical relativity simulations play a central role in the study of these events. In particular, given the importance of magnetic fields, various aspects of this investigation require general relativistic magnetohydrodynamics (GRMHD). So far, most GRMHD simulations focused the attention on BNS mergers leading to the formation of a hypermassive NS, which, in turn, collapses within few tens of ms into a black hole surrounded by an accretion disk. However, recent observations suggest that a significant fraction of these systems could form a long-lived NS remnant, which will either collapse on much longer timescales or remain indefinitely stable. Despite the profound implications for the evolution and the emission properties of the system, a detailed investigation of this alternative evolution channel is still missing. Here, we follow this direction and present a first detailed GRMHD study of BNS mergers forming a long-lived NS. We consider magnetized binaries with different mass ratios and equations of state and analyze the structure of the NS remnants, the rotation profiles, the accretion disks, the evolution and amplification of magnetic fields, and the ejection of matter. Moreover, we discuss the connection with the central engine of SGRBs and provide order-of-magnitude estimates for the kilonova/macronova signal. Finally, we study the GW emission, with particular attention to the post-merger phase.

We present fully GRMHD simulations of the merger of binary neutron star (BNS) systems. We consider BNSs producing a hypermassive neutron star (HMNS) that collapses to a spinning black hole (BH) surrounded by a magnetized accretion disk in a few tens of ms. We investigate whether such systems may launch relativistic jets and power short gamma-ray bursts. We study the effects of different equations of state (EOSs), different mass ratios, and different magnetic field orientations. For all cases, we present a detailed investigation of the matter dynamics and of the magnetic field evolution, with particular attention to its global structure and possible emission of relativistic jets. The main result of this work is that we found the formation of an organized magnetic field structure. This happens independently of EOS, mass ratio, and initial magnetic field orientation. We also show that those models that produce a longer-lived HMNS lead to a stronger magnetic field before collapse to BH. Such larger fields make it possible, for at least one of our models, to resolve the MRI and hence further amplify the magnetic field. However, by the end of our simulations, we do not observe a magnetically dominated funnel and hence neither a relativistic outflow. With respect to the recent simulations of Ruiz et al 2016, we evolve models with lower and more realistic initial magnetic field strengths and, because of computational reasons, we do not evolve the accretion disk for the long timescales that seem to be required in order to see a relativistic outflow. Since all our models produce a similar ordered magnetic field structure, we expect that the results found in Ruiz et al 2016, where they only considered an equal-mass system with an ideal fluid EOS, should be general and, at least from a qualitative point of view, independent from mass-ratio, magnetic field orientation, and EOS.

We present new results of fully general relativistic magnetohydrodynamic simulations of binary neutron star (BNS) mergers performed with the Whisky code. All the models use a piecewise polytropic approximation of the APR4 equation of state for cold matter, together with a 'hybrid' part to incorporate thermal effects during the evolution. We consider both equal and unequal-mass models, with total masses such that either a supramassive NS or a black hole is formed after merger. Each model is evolved with and without a magnetic field initially confined to the stellar interior. We present the different gravitational wave (GW) signals as well as a detailed description of the matter dynamics (magnetic field evolution, ejected mass, post-merger remnant/disk properties). Our simulations provide new insights into BNS mergers, the associated GW emission and the possible connection with the engine of short gamma-ray bursts (both in the 'standard' and in the 'time-reversal' scenarios) and other electromagnetic counterparts.