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.

We report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner-Bethe-Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne $v_{18}$ interaction and two of its simplified versions, and four of the new Norfolk $\Delta$-full chiral effective field theory potentials. The results up to twice nuclear matter saturation density show some divergence among the methods, but improved agreement compared to earlier work. We find that the potentials fit to higher-energy nucleon-nucleon scattering data exhibit a much smaller spread of energies.

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 adopt two- and three-body nuclear forces derived at the next-to-next-to-leading-order (N2LO) in the framework of effective chiral perturbation theory (ChPT) to calculate the equation of state (EOS) of $\beta$-stable neutron star matter using the Brueckner--Hartree--Fock many-body approach. We use the recent optimized chiral two-body nuclear interaction at N2LO derived by \cite{ekstrom1} and two different parametrizations of the three-body N2LO interaction: the first one is fixed to reproduce the saturation point of symmetric nuclear matter while the second one is fixed to reproduce the binding energies of light atomic nuclei. We show that in the second case the properties of nuclear matter are not well determined whereas in the first case various empirical nuclear matter properties around the saturation density are well reproduced. We also calculate the nuclear symmetry energy $E_{sym}$ as a function of the nucleonic density and compare our results with the empirical constraints obtained using the excitation energies of isobaric analog states in nuclei and the experimental data on the neutron skin thickness of heavy nuclei. We next calculate various neutron star properties and in particular the mass-radius and mass-central density relations. We find that the adopted interactions based on a fully microscopic framework, are able to provide an EOS which is consistent with the present data of measured neutron star masses and in particular with the mass $M=2.01\pm0.04 M_\odot$ of the neutron star in PSR J0348+0432. We finally consider the possible presence of hyperons in the stellar core and we find a softening of the EOS and a substantial reduction of the stellar maximum mass in agreement with similar calculations present in the literature.

We report a new microscopic equation of state (EOS) of dense symmetric nuclear matter, pure neutron matter, and asymmetric and $\beta$-stable nuclear matter at zero temperature using recent realistic two-body and three-body nuclear interactions derived in the framework of chiral perturbation theory (ChPT) and including the $\Delta(1232)$ isobar intermediate state. This EOS is provided in tabular form and in parametrized form ready for use in numerical general relativity simulations of binary neutron star merging. Here we use our new EOS for $\beta$-stable nuclear matter to compute various structural properties of non-rotating neutron stars.The EOS is derived using the Brueckner--Bethe--Goldstone quantum many-body theory in the Brueckner--Hartree--Fock approximation. Neutron star properties are next computed solving numerically the Tolman--Oppenheimer--Volkov structure equations. Our EOS models are able to reproduce the empirical saturation point of symmetric nuclear matter, the symmetry energy $E_{sym}$, and its slope parameter $L$ at the empirical saturation density $n_{0}$. In addition, our EOS models are compatible with experimental data from collisions between heavy nuclei at energies ranging from a few tens of MeV up to several hundreds of MeV per nucleon. These experiments provide a selective test for constraining the nuclear EOS up to $\sim 4 n_0$. Our EOS models are consistent with present measured neutron star masses and particularly with the mass $M = 2.01 \pm 0.04 \, M_{\odot}$ of the neutron stars in PSR~J0348+0432.

It is usually thought that a single equation of state (EoS) model "correctly" represents cores of all compact stars. Here we emphasize that two families of compact stars, viz., neutron stars and strange stars, can coexist in nature, and that neutron stars can get converted to strange stars through the nucleation process of quark matter in the stellar center. From our fully general relativistic numerical computations of the structures of fast-spinning compact stars, known as millisecond pulsars, we find that such a stellar conversion causes a simultaneous spin-up and decrease in gravitational mass of these stars. This is a new type of millisecond pulsar evolution through a new mechanism, which gives rise to relatively lower mass compact stars with higher spin rates. This could have implication for the observed mass and spin distributions of millisecond pulsars. Such a stellar conversion can also rescue some massive, spin-supported millisecond pulsars from collapsing into black holes. Besides, we extend the concept of critical mass $M_{\rm cr}$ for the neutron star sequence (Berezhiani et al. 2003; Bombaci et al. 2004) to the case of fast-spinning neutron stars, and point out that neutron star EoS models cannot be ruled out by the stellar mass measurement alone. Finally, we emphasize the additional complexity for constraining EoS models, for example, by stellar radius measurements using X-ray observations, if two families of compact stars coexist.

The recent fast growth of a population of millisecond pulsars with precisely measured mass provides an excellent opportunity to characterize these compact stars at an unprecedented level. This is because the stellar parameter values can be accurately computed for known mass and spin rate and an assumed equation of state (EoS) model. For each of the 16 such pulsars and for a set of EoS models from nucleonic, hyperonic, strange quark matter and hybrid classes, we numerically compute fast spinning stable stellar parameter values considering the full effect of general relativity. This first detailed catalogue of the computed parameter values of observed millisecond pulsars provides a testbed to probe the physics of compact stars, including their formation, evolution and EoS. We estimate uncertainties on these computed values from the uncertainty of the measured mass, which could be useful to quantitatively constrain EoS models. We note that the largest value of the central density $\rho_{\rm c}$ in our catalogue is $\sim 5.8$ times the nuclear saturation density $\rho_{\rm sat}$, which is much less than the expected maximum value $13 \rho_{\rm sat}$. We argue that the $\rho_{\rm c}$-values of at most a small fraction of compact stars could be much larger than $5.8 \rho_{\rm sat}$. Besides, we find that the constraints on EoS models from accurate radius measurements could be significantly biased for some of our pulsars, if stellar $spinning$ configurations are not used to compute the theoretical radius values.

For a set of equation of state (EoS) models involving interacting strange quark matter, characterized by an effective bag constant (B_eff) and a perturbative QCD corrections term (a_4), we construct fully general relativistic equilibrium sequences of rapidly spinning strange stars for the first time. Computation of such sequences is important to study millisecond pulsars and other fast spinning compact stars. Our EoS models can support a gravitational mass (M_G) and a spin frequency at least up to approximately 3.0 solar mass and approximately 1250 Hz respectively, and hence are fully consistent with measured M_G and spin frequency values. This paper reports the effects of B_eff and a_4 on measurable compact star properties, which could be useful to find possible ways to constrain these fundamental quark matter parameters, within the ambit of our EoS models. We confirm that a lower B_eff allows a higher mass. Besides, for known M_G and spin frequency, measurable parameters, such as stellar radius, radius-to-mass ratio and moment of inertia, increase with the decrease of B_eff. Our calculations also show that a_4 significantly affects the stellar rest mass and the total stellar binding energy. As a result, a_4 can have signatures in evolutions of both accreting and non-accreting compact stars, and the observed distribution of stellar mass and spin and other source parameters. Finally, we compute the parameter values of two important pulsars, PSR J1614-2230 and PSR J1748-2446ad, which may have implications to probe their evolutionary histories, and for constraining EoS models.

A phase of strong interacting matter with deconfined quarks is expected in the core of massive neutron stars. We investigate the quark deconfinement phase transition in cold (T = 0) and hot beta-stable hadronic matter. Assuming a first order phase transition, we calculate and compare the nucleation rate and the nucleation time due to quantum and thermal nucleation mechanisms. We show that above a threshold value of the central pressure a pure hadronic star (HS) (i.e. a compact star with no fraction of deconfined quark matter) is metastable to the conversion to a quark star (QS) (i.e. a hybrid star or a strange star). This process liberates an enormous amount of energy, of the order of 10^{53}~erg, which causes a powerful neutrino burst, likely accompanied by intense gravitational waves emission, and possibly by a second delayed (with respect to the supernova explosion forming the HS) explosion which could be the energy source of a powerful gamma-ray burst (GRB). This stellar conversion process populates the QS branch of compact stars, thus one has in the Universe two coexisting families of compact stars: pure hadronic stars and quark stars. We introduce the concept of critical mass M_{cr} for cold HSs and proto-hadronic stars (PHSs), and the concept of limiting conversion temperature for PHSs. We show that PHSs with a mass M < M_{cr} could survive the early stages of their evolution without decaying to QSs. Finally, we discuss the possible evolutionary paths of proto-hadronic stars.

We calculate the energy per particle of symmetric nuclear matter and pure neutron matter using the microscopic many-body Brueckner-Hartree-Fock (BHF) approach and employing the Argonne V18 (AV18) nucleon-nucleon (NN) potential supplemented with two different three-nucleon force models recently constructed to reproduce the binding energy of $^3$H, $^3$He and $^4$He nuclei as well as the neutron-deuteron doublet scattering length. We find that none of these new three-nucleon force models is able to reproduce simultaneously the empirical saturation point of symmetric nuclear matter and the properties of three- and four-nucleon systems.

We study the consequences of quark matter nucleation in cold hadronic matter employing three relativistic-mean-field (RMF) models to describe the hadronic phase and the Nambu-Jona-Lasinio (NJL) model for the quark one. We explore the effect of a vector interaction in the NJL Lagrangian and of a phenomenological bag constant on neutron stars metastability. We delineate the region of parameters of the quark phase that allow for the formation of stable hybrid stars with mass compatible with the almost $2 \ M_\odot$ pulsars PSR J1614-2230 ($1.97 \pm 0.04 M_\odot$) and PSR J0348+0432 ($2.01 \pm 0.04 M_\odot$). It is shown, however, that not all hybrid star configurations with $\sim 2 \,M_\odot$ are populated after nucleation.

A phase of strong interacting matter with deconfined quarks is expected in the core of massive neutron stars. In this article, we perform a study of the hadron-quark phase transition in cold (T = 0) neutron star matter and we calculate various structural properties of hybrid stars. For the quark phase, we make use of an equation of state (EOS) derived with the field correlator method (FCM) recently extended to the case of nonzero baryon density. For the hadronic phase, we consider both pure nucleonic and hyperonic matter, and we derive the corresponding EOS within a relativistic mean field approach. We make use of measured neutron star masses, and particularly the mass $M = 1.97 \pm 0.04 \, M_\odot$ of PSR J1614 -2230 to constrain the values of the gluon condensate $G_2$, which is one of the EOS parameters within the FCM. We find that the values of $G_2$ extracted from the mass measurement of PSR J1614 -2230 are consistent with the values of the same quantity derived within the FCM from recent lattice QCD calculations of the deconfinement transition temperature at zero baryon chemical potential. The FCM thus provides a powerful tool to link numerical calculations of QCD on a space-time lattice with measured neutron star masses.

We study the hadron-quark phase transition in neutron star matter and the structural properties of hybrid stars using an equation of state (EOS) for the quark phase derived with the field correlator method (FCM). We make use of the measured neutron star masses, and particularly the mass of PSR J1614-2230, to constrain the values of the gluon condensate $G_2$ which is one of the EOS parameter within the FCM. We find that the values of $G_2$ extracted from the mass measurement of PSR J1614-2230 are fully consistent with the values of the same quantity derived, within the FCM, from recent lattice quantum chromodynamics (QCD) calculations of the deconfinement transition temperature at zero baryon chemical potential. The FCM thus provides a powerful tool to link numerical calculations of QCD on a space-time lattice with neutron stars physics.

The nucleation process of quark matter in cold (T = 0) stellar matter is investigated using the microscopic Brueckner-Hartree-Fock approach to describe the hadronic phase, and the MIT bag model, the Nambu-Jona-Lasinio, and the Chromo Dielectric models to describe the deconfined phase of quark matter. The consequences of the nucleation process for neutron star physics are outlined. Hyperonic stars are metastable only for some of the quark matter equations of state considered. The effect of an hyperonic three body force on the metastability of compact stars is estimated, and it is shown that, except for the Nambu-Jona-Lasinio model and the MIT bag model with a large bag pressure, the other models predict the formation of hybrid stars with a maximum mass not larger than \sim 1.62 M\odot .

(Abridged) A phase of strong interacting matter with deconfined quarks is expected in the core of massive neutron stars. If this deconfinement phase transition is of the first order then it will be triggered by the nucleation of a critical size drop of the stable quark phase in the metastable hadronic phase. Within these circumstances it has been shown that cold pure hadronic compact stars above a threshold value of their gravitational mass are metastable with respect to the "decay" to quark stars (compact stars made at least in part of quark matter). This stellar conversion process liberates a huge amount of energy, and it could be the energy source of some of the long GRBs. The main goal of the present work is to establish whether a newborn hadronic star (proto-hadronic star) could survive the early stages of its evolution without "decaying" to a quark star. To this aim, we study the nucleation process of quark matter in hot beta-stable hadronic matter, with and without trapped neutrinos. We calculate and compare the nucleation rate and the nucleation time due to thermal and quantum nucleation mechanisms. We compute the crossover temperature above which thermal nucleation dominates the finite temperature quantum nucleation mechanism. We next discuss the consequences of quark matter nucleation for the physics and the evolution of proto-neutron stars. We introduce the new concept of limiting conversion temperature and critical mass M_cr for proto-hadronic stars, and we show that proto-hadronic stars with a mass M < M_cr could survive the early stages of their evolution without decaying to a quark star. We extend the concept of maximum mass of a "neutron star" with respect to the classical one introduced by Oppenheimer & Volkoff to account for the existence of two distinct families of compact stars (hadronic stars and quark stars) as predicted by the present scenario.

We study the quark deconfinement phase transition in hot $\beta$-stable hadronic matter. Assuming a first order phase transition, we calculate the enthalpy per baryon of the hadron-quark phase transition. We calculate and compare the nucleation rate and the nucleation time due to thermal and quantum nucleation mechanisms. We compute the crossover temperature above which thermal nucleation dominates the finite temperature quantum nucleation mechanism. We next discuss the consequences for the physics of proto-neutron stars. We introduce the concept of limiting conversion temperature and critical mass $M_{cr}$ for proto-hadronic stars, and we show that proto-hadronic stars with a mass $M < M_{cr}$ could survive the early stages of their evolution without decaying to a quark star.