The recent detection of GW190814 featured the merger of a binary with a primary having a mass of $\sim 23\,M_{\odot}$ and a secondary with a mass of $\sim 2.6\,M_{\odot}$. While the primary was most likely a black hole, the secondary could be interpreted as either the lightest black hole or the most massive neutron star ever observed, but also as the indication of a novel class of exotic compact objects. We here argue that the secondary in GW190814 needs not be an ab-initio black hole nor an exotic object; rather, based on our current understanding of the nuclear-matter equation of state, it can be a rapidly rotating neutron star that collapsed to a rotating black hole at some point before merger. Using universal relations connecting the masses and spins of uniformly rotating neutron stars, we estimate the spin, $0.49 \lesssim \chi \lesssim 0.68$, of the secondary -- a quantity not constrained so far by the detection -- and a novel strict lower bound on the maximum mass, $M_{\rm TOV} > 2.08^{+0.04}_{-0.04}\, \,M_{\odot}$, of nonrotating neutron stars, consistent with recent observations of a very massive pulsar. The new lower bound also remains valid even in the less likely scenario in which the secondary neutron star never collapsed to a black hole.

We investigate the impact of using high-order numerical methods to study the merger of magnetised neutron stars with finite-temperature microphysics and neutrino cooling in full general relativity. By implementing a fourth-order accurate conservative finite-difference scheme we model the inspiral together with the early post-merger and highlight the differences to traditional second-order approaches at the various stages of the simulation. We find that even for finite-temperature equations of state, convergence orders higher than second order can be achieved in the inspiral and post-merger for the gravitational-wave phase. We further demonstrate that the second-order scheme overestimates the amount of proton-rich shock-heated ejecta, which can have an impact on the modelling of the dynamical part of the kilonova emission. Finally, we show that already at low resolution the growth rate of the magnetic energy is consistently resolved by using a fourth-order scheme.

We study in detail the nuclear aspects of a neutron-star merger in which deconfinement to quark matter takes place. For this purpose, we make use of the Chiral Mean Field (CMF) model, an effective relativistic model that includes self-consistent chiral symmetry restoration and deconfinement to quark matter and, for this reason, predicts the existence of different degrees of freedom depending on the local density/chemical potential and temperature. We then use the out-of-chemical-equilibrium finite-temperature CMF equation of state in full general-relativistic simulations to analyze which regions of different QCD phase diagrams are probed and which conditions, such as strangeness and entropy, are generated when a strong first-order phase transition appears. We also investigate the amount of electrons present in different stages of the merger and discuss how far from chemical equilibrium they can be and, finally, draw some comparisons with matter created in supernova explosions and heavy-ion collisions.

Following the detection of GW170817 and the accompanying kilonova AT2017gfo, it has become crucial to model and understand the various channels through which mass is ejected in neutron-star binary mergers. We discuss the impact that high stellar spins prior to merger have on the ejection of mass focussing, in particular, on the dynamically ejected mass by performing general-relativistic magnetohydrodynamic simulations employing finite-temperature equations of state and neutrino-cooling effects. Using eight different models with dimensionless spins ranging from $\chi\simeq-0.14$ to $\chi\simeq0.29$ we discuss how the presence of different spins affects the angular distribution and composition of the ejected matter. Most importantly, we find that the dynamical component of the ejected mass can be strongly suppressed in the case of high spins aligned with the orbital angular momentum. In this case, in fact, the merger remnant has an excess angular momentum yielding a more extended and "colder" object, with reduced ability to shed mass dynamically. We discuss how this result impacts the analysis of the recent merger event GW170817 and its kilonova afterglow.

In this work, we discuss the dense matter equation of state (EOS) for the extreme range of conditions encountered in neutron stars and their mergers. The calculation of the properties of such an EOS involves modeling different degrees of freedom (such as nuclei, nucleons, hyperons, and quarks), taking into account different symmetries, and including finite density and temperature effects in a thermodynamically consistent manner. We begin by addressing subnuclear matter consisting of nucleons and a small admixture of light nuclei in the context of the excluded volume approach. We then turn our attention to supranuclear homogeneous matter as described by the Chiral Mean Field (CMF) formalism. Finally, we present results from realistic neutron-star-merger simulations performed using the CMF model that predict signatures for deconfinement to quark matter in gravitational wave signals.

Merging binaries of neutron stars are not only strong sources of gravitational waves, but also have the potential of revealing states of matter at densities and temperatures not accessible in laboratories. A crucial and long-standing question in this context is whether quarks are deconfined as a result of the dramatic increase in density and temperature following the merger. We present the first fully general-relativistic simulations of merging neutron stars including quarks at finite temperatures that can be switched off consistently in the equation of state. Within our approach, we can determine clearly what signatures a quark-hadron phase transition would leave in the gravitational-wave signal. In particular, we show that if the conditions are met for a phase transition to take place at several times nuclear saturation density, they would lead to a post-merger signal considerably different from the one expected from the inspiral, that can only probe the hadronic part of the equations of state, and to an anticipated collapse of the merged object. We also show that the phase transition leads to a very hot and dense quark core that, when it collapses to a black hole, produces a ringdown signal different from the hadronic one. Finally, in analogy with what is done in heavy-ion collisions, we use the evolution of the temperature and density in the merger remnant to illustrate the properties of the phase transition in a QCD phase diagram.

With the recent advent of multi-messenger gravitational-wave astronomy and in anticipation of more sensitive, next-generation gravitational-wave detectors, we investigate the dynamics, gravitational-wave emission, and nucleosynthetic yields of numerous eccentric binary neutron-star mergers having different equations of state. For each equation of state we vary the orbital properties around the threshold of immediate merger, as well as the binary mass ratio. In addition to a study of the gravitational-wave emission including $f$-mode oscillations before and after merger, we couple the dynamical ejecta output from the simulations to the nuclear-reaction network code \texttt{SkyNet} to compute nucleosynthetic yields and compare to the corresponding results in the case of a quasi-circular merger. We find that the amount and velocity of dynamically ejected material is always much larger than in the quasi-circular case, reaching maximal values of $M_{\rm ej, max} \sim 0.1 \, M_{\odot}$ and $v_{\rm max}/c \sim 0.75$. At the same time, the properties of this material are rather insensitive to the details of the orbit, such as pericenter distance or post-encounter apoastron distance. Furthermore, while the composition of the ejected matter depends on the orbital parameters and on the equation of state, the relative nucleosynthetic yields do not, thus indicating that kilonova signatures could provide information on the orbital properties of dynamically captured neutron-star binaries.