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.

The first binary neutron-star merger event, GW170817, and its bright electromagnetic counterpart have provided a remarkable amount of information. By contrast, the second event, GW190425, with $M_{\rm tot}=3.4^{+0.3}_{-0.1}\,M_{\odot}$ and the lack of an electromagnetic counterpart, has hardly improved our understanding of neutron-star physics. While GW190425 is compatible with a scenario in which the merger has lead to a prompt collapse to a black hole and little ejected matter to power a counterpart, determining the mass ratio and the effective spin $\tilde{\chi}$ of the binary remains difficult. This is because gravitational waveforms cannot yet well constrain the component spins of the binary. However, since the mass of GW190425 is significantly larger than the maximum mass for nonrotating neutron stars, $M_{_{\rm TOV}}$, the mass ratio $q$ cannot be too small, as the heavier star would not be gravitationally stable. Making use of universal relations and a large number of equations of state, we provide limits in the $(\tilde{\chi},q)$ plane for GW190425, namely: $q_{\rm min} \geq 0.38$ and $ \tilde{\chi}_{\rm max}\leq 0.20$, assuming $M_\mathrm{tot} \simeq 3.4\, M_\odot$. Finally, we show how future observations of high-mass binaries can provide a lower bound on $M_{_{\rm TOV}}$.

The detection of gravitational waves from neutron star merger events has opened up a new field of multi-messenger astronomy linking gravitational waves events to short-gamma ray bursts and kilonova afterglows. A further - yet to be discovered - electromagnetic counterpart is precursor emission produced by the non-trivial interaction of the magnetospheres of the two neutron stars prior to merger. By performing special-relativistic force-free simulations of orbiting neutron stars we discuss the effect of different magnetic field orientations and show how the emission can be significantly enhanced by differential motion present in the binary, either due to stellar spins or misaligned stellar magnetospheres. We find that the built-up of twist in the magnetic flux tube connecting the two stars can lead to the repeated emission of powerful flares for a variety of orbital configurations. We also discuss potential coherent radio emission mechanisms in the flaring process.

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.

Exploiting a very large library of physically plausible equations of state (EOSs) containing more than $10^{7}$ members and yielding more than $10^{9}$ stellar models, we conduct a survey of the impact that a neutron-star radius measurement via electromagnetic observations can have on the EOS of nuclear matter. Such measurements are soon to be expected from the ongoing \textit{Neutron Star Interior Composition Explorer} mission and will complement the constraints on the EOS from gravitational-wave detections. Thanks to the large statistical range of our EOS library, we can obtain a first quantitative estimate of the commonly made assumption that the high-density part of the EOS is best constrained when measuring the radius of the most massive, albeit rare, neutron stars with masses $M\gtrsim2.1\,M_\odot$. At the same time, we find that radius measurements of neutron stars with masses $M\simeq1.7-1.85\,M_\odot$ can provide the strongest constraints on the low-density part of the EOS. Finally, we quantify how radius measurements by future missions can further improve our understanding of the EOS of matter at nuclear densities.

Recent advances in black hole astrophysics, particularly the first visual evidence of a supermassive black hole at the center of the galaxy M87 by the Event Horizon Telescope (EHT), and the detection of an orbiting "hot spot" nearby the event horizon of Sgr A* in the Galactic center by the Gravity Collaboration, require the development of novel numerical methods to understand the underlying plasma microphysics. Non-thermal emission related to such hot spots is conjectured to originate from plasmoids that form due to magnetic reconnection in thin current layers in the innermost accretion zone. Resistivity plays a crucial role in current sheet formation, magnetic reconnection, and plasmoid growth in black hole accretion disks and jets. We included resistivity in the three-dimensional general-relativistic magnetohydrodynamics (GRMHD) code BHAC and present the implementation of an Implicit-Explicit scheme to treat the stiff resistive source terms of the GRMHD equations. The algorithm is tested in combination with adaptive mesh refinement to resolve the resistive scales and a constrained transport method to keep the magnetic field solenoidal. Several novel methods for primitive variable recovery, a key part in relativistic magnetohydrodynamics codes, are presented and compared for accuracy, robustness, and efficiency. We propose a new inversion strategy that allows for resistive-GRMHD simulations of low gas-to-magnetic pressure ratio and highly magnetized regimes as applicable for black hole accretion disks, jets, and neutron star magnetospheres. We apply the new scheme to study the effect of resistivity on accreting black holes, accounting for dissipative effects as reconnection.

Worldwide very long baseline radio interferometry arrays are expected to obtain horizon-scale images of supermassive black hole candidates as well as of relativistic jets in several nearby active galactic nuclei. This motivates the development of models for magnetohydrodynamic flows in strong gravitational fields. The Black Hole Accretion Code (BHAC) intends to aid with the modelling of such sources by means of general relativistic magnetohydrodynamical (GRMHD) simulations in arbitrary stationary spacetimes. New additions were required to guarantee an accurate evolution of the magnetic field when small and large scales are captured simultaneously. We discuss the adaptive mesh refinement (AMR) techniques employed in BHAC, essential to keep several problems computationally tractable, as well as staggered-mesh-based constrained transport (CT) algorithms to preserve the divergence-free constraint of the magnetic field, including a general class of prolongation operators for face-allocated variables compatible with them. Through several standard tests, we show that the choice of divergence-control method can produce qualitative differences in simulations of scientifically relevant accretion problems. We demonstrate the ability of AMR to reduce the computational costs of accretion simulations while sufficiently resolving turbulence from the magnetorotational instability. In particular, we describe a simulation of an accreting Kerr black hole in Cartesian coordinates using AMR to follow the propagation of a relativistic jet while self-consistently including the jet engine, a problem set up-for which the new AMR implementation is particularly advantageous. The CT methods and AMR strategies discussed here are being employed in the simulations performed with BHAC used in the generation of theoretical models for the Event Horizon Telescope Collaboration.

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.

It has been suggested that a non-repeating fast radio burst (FRB) represents the final signal of a magnetized neutron star collapsing to a black hole. In this model, a supramassive neutron star supported by rapid rotation, will collapse to a black hole several thousand to million years after its birth as a result of spin down. The collapse violently snaps the magnetic-field lines anchored on the stellar surface, thus producing an electromagnetic pulse that will propagate outwards and accelerate electrons producing a massive radio burst, i.e. a "blitzar". We present a systematic study of the gravitational collapse of rotating and magnetised neutron stars with special attention to far-field evolution at late times after the collapse. By considering a series of neutron stars with rotation ranging from zero to millisecond periods and different magnetic-field strengths, we show that the blitzar emission is very robust and always characterised by a series sub-millisecond pulses decaying exponentially in amplitude. The luminosity and energy released when the magnetosphere is destroyed are well reproduced by a simple expression in terms of the stellar magnetic field and radius. Finally, we assess the occurrence of pair production during a blitzar scenario, concluding that for typical magnetic-field strengths of $10^{12}\,{\rm G}$ and spin frequencies of a few Hz, pair production is suppressed. Overall, the very good match between the results of the simulations and the luminosities normally observed for FRBs lends credibility to the blitzar model as a simple and yet plausible explanation for the phenomenology of non-repeating FRBs.

We explore in a parameterized manner a very large range of physically plausible equations of state (EOSs) for compact stars for matter that is either purely hadronic or that exhibits a phase transition. In particular, we produce two classes of EOSs with and without phase transitions, each containing one million EOSs. We then impose constraints on the maximum mass, ($M < 2.16 M_{\odot}$), and on the dimensionless tidal deformability ($\tilde{\Lambda} <800$) deduced from GW170817, together with recent suggestions of lower limits on $\tilde{\Lambda}$. Exploiting more than $10^9$ equilibrium models for each class of EOSs, we produce distribution functions of all the stellar properties and determine, among other quantities, the radius that is statistically most probable for any value of the stellar mass. In this way, we deduce that the radius of a purely hadronic neutron star with a representative mass of $1.4\,M_{\odot}$ is constrained to be $12.00\!<\!R_{1.4}/{\rm km}\!<\!13.45$ at a $2$-$\sigma$ confidence level, with a most likely value of $\bar{R}_{1.4}=12.39\,{\rm km}$; similarly, the smallest dimensionless tidal deformability is $\tilde{\Lambda}_{1.4}\!>\!375$, again at a $2$-$\sigma$ level. On the other hand, because EOSs with a phase transition allow for very compact stars on the so-called `twin-star' branch, small radii are possible with such EOSs although not probable, i.e. $8.53\!<\!R_{1.4}/{\rm km}\!<\!13.74$ and $\bar{R}_{1.4}=13.06\,{\rm km}$ at a $2$-$\sigma$ level, with $\tilde{\Lambda}_{1.4}\!>\!35.5$ at a $3$-$\sigma$ level. Finally, since these EOSs exhibit upper limits on $\tilde{\Lambda}$, the detection of a binary with total mass of $3.4\,M_{\odot}$ and $\tilde{\Lambda}_{1.7}\!>\!461$ can rule out twin-star solutions.

Combining the GW observations of merging systems of binary neutron stars and quasi-universal relations, we set constraints on the maximum mass that can be attained by nonrotating stellar models of neutron stars. More specifically, exploiting the recent observation of the GW event GW 170817 and drawing from basic arguments on kilonova modeling of GRB 170817A, together with the quasi-universal relation between the maximum mass of nonrotating stellar models $M_{\rm TOV}$ and the maximum mass supported through uniform rotation $M_{\rm max}=\left(1.20^{+0.02}_{-0.05}\right) M_{\rm TOV}$ we set limits for the maximum mass to be $ 2.01^{+0.04}_{-0.04}\leq M_{\rm TOV}/M_{\odot}\lesssim 2.16^{+0.17}_{-0.15}$, where the lower limit in this range comes from pulsar observations. Our estimate, which follows a very simple line of arguments and does not rely on the modeling of the electromagnetic signal in terms of numerical simulations, can be further refined as new detections become available. We briefly discuss the impact that our conclusions have on the equation of state of nuclear matter.

The stability properties of rotating relativistic stars against prompt gravitational collapse to a black hole are rather well understood for uniformly rotating models. This is not the case for differentially rotating neutron stars, which are expected to be produced in catastrophic events such as the merger of binary system of neutron stars or the collapse of a massive stellar core. We consider sequences of differentially rotating equilibrium models using the $j$-constant law and by combining them with their dynamical evolution, we show that a sufficient stability criterion for differentially rotating neutron stars exists similar to the one of their uniformly rotating counterparts. Namely: along a sequence of constant angular momentum, a dynamical instability sets in for central rest-mass densities slightly below the one of the equilibrium solution at the turning point. In addition, following Breu & Rezzolla (2016), we show that "quasi-universal" relations can be found when calculating the turning-point mass. In turn, this allows us to compute the maximum mass allowed by differential rotation, $M_{\rm max,dr}$, in terms of the maximum mass of the nonrotating configuration, $M_{_{\rm TOV}}$, finding that $M_{\rm max, dr} \simeq \left(1.54 \pm 0.05\right) M_{_{\rm TOV}}$ for all the equations of state we have considered.

We present the first systematic study of the gravitational collapse of rotating and magnetised neutron stars to charged and rotating (Kerr-Newman) black holes. In particular, we consider the collapse of magnetised and rotating neutron stars assuming that no pair-creation takes place and that the charge density in the magnetosphere is so low that the stellar exterior can be described as an electrovacuum. Under these assumptions, which are rather reasonable for a pulsar that has crossed the 'death line', we show that when the star is rotating, it acquires a net initial electrical charge, which is then trapped inside the apparent horizon of the newly formed back hole. We analyse a number of different quantities to validate that the black hole produced is indeed a Kerr-Newman one and show that, in the absence of rotation or magnetic field, the end result of the collapse is a Schwarzschild or Kerr black hole, respectively.