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.

Degeneracy in the Inference of Phase Transitions in the Neutron Star Equation of State from Gravitational Wave Data

Gravitational wave (GW) detections of binary neutron star inspirals will be crucial for constraining the dense matter equation of state (EOS). We demonstrate a new degeneracy in the mapping from tidal deformability data to the EOS, which occurs for models with strong phase transitions. We find that there exists a new family of EOS with phase transitions that set in at different densities and that predict neutron star radii that differ by up to ∼500  m but that produce nearly identical tidal deformabilities for all neutron star masses. Next-generation GW detectors and advances in nuclear theory may be needed to resolve this degeneracy.

A physics-based neural network reconstruction of the dense matter equation of state from neutron star observables

We introduce a novel technique that utilizes a physics-driven deep learning method to reconstruct the dense matter equation of state from neutron star observables, particularly the masses and radii. The proposed framework involves two neural networks: one to optimize the EoS using Automatic Differentiation in the unsupervised learning scheme; and a pre-trained network to solve the Tolman–Oppenheimer–Volkoff (TOV) equations. The gradient-based optimization process incorporates a Bayesian picture into the proposed framework. The reconstructed EoS is proven to be consistent with the results from conventional methods. Furthermore, the resulting tidal deformation is in agreement with the limits obtained from the gravitational wave event, GW170817.

Trace Anomaly as Signature of Conformality in Neutron Stars

We discuss an interpretation that a peak in the sound velocity in neutron star matter, as suggested by the observational data, signifies strongly coupled conformal matter. The normalized trace anomaly is a dimensionless measure of conformality leading to the derivative and the nonderivative contributions to the sound velocity. We find that the peak in the sound velocity is attributed to the derivative contribution from the trace anomaly that steeply approaches the conformal limit. Smooth continuity to the behavior of high-density QCD implies that the matter part of the trace anomaly may be positive definite. We discuss a possible implication of the positivity condition of the trace anomaly on the M-R relation of the neutron stars.

A hybrid GNA instability

A semi-analytic admixed model formalism to study the stability effects of the inner crust regions against the local collective perturbations in non-rotating neutron stars is proposed. It consists of the viscoelastic heavy neutron-rich nuclei, superfluid neutrons, and degenerate quantum electrons. A normal spherical mode analysis yields a generalized linear dispersion relation multiparametrically mimicking the inner crust features of neutron stars. A hybrid gravito-nucleo-acoustic (GNA) instability mode is found to be excited. It is demonstrated that the electron density and the inner crust curvature act as its accelerating and antidispersive agents. In contrast, the heavy neutron-rich nucleus and neutron densities act as decelerating factors. The heavy nucleus density, electron density, and geometric curvature act as its destabilizers. It is only the neutron density that acts as the GNA stabilizing agent. The heavy neutron-rich nucleus and neutron densities are found to act as dispersive broadening factors to it. The high-[Formula: see text] regions are the more unstable spectral windows indicating that the GNA mode plays a dominant role in the inner crust zone towards the local stability. Its fair reliability is indicated in light of the recent astronomic observed scenarios. It could be useful to explore acoustic mode signatures in non-rotating neutron stars and similar other compact astroobjects.

Compact Binary Coalescences: Astrophysical Processes and Lessons Learned

On 11 February 2016, the LIGO and Virgo scientific collaborations announced the first direct detection of gravitational waves, a signal caught by the LIGO interferometers on 14 September 2015, and produced by the coalescence of two stellar-mass black holes. The discovery represented the beginning of an entirely new way to investigate the Universe. The latest gravitational-wave catalog by LIGO, Virgo and KAGRA brings the total number of gravitational-wave events to 90, and the count is expected to significantly increase in the next years, when additional ground-based and space-born interferometers will be operational. From the theoretical point of view, we have only fuzzy ideas about where the detected events came from, and the answers to most of the five Ws and How for the astrophysics of compact binary coalescences are still unknown. In this work, we review our current knowledge and uncertainties on the astrophysical processes behind merging compact-object binaries. Furthermore, we discuss the astrophysical lessons learned through the latest gravitational-wave detections, paying specific attention to the theoretical challenges coming from exceptional events (e.g., GW190521 and GW190814).

Thinking Outside the Box: Numerical Relativity with Particles

The observation of gravitational waves from compact objects has now become an active part of observational astronomy. For a sound interpretation, one needs to compare such observations against detailed Numerical Relativity simulations, which are essential tools to explore the dynamics and physics of compact binary mergers. To date, essentially all simulation codes that solve the full set of Einstein’s equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations, we focus on perfectly symmetric binary systems that are irrotational and have 1.3 M⊙ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle–mesh mapping), and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the Ψ4 formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times 10−3M⊙ as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity (∼0.5…0.7c) ejecta component of ∼10−4M⊙ that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission, and it could also potentially cause a kilonova afterglow.

Constraining neutron-star matter with microscopic and macroscopic collisions

Interpreting high-energy, astrophysical phenomena, such as supernova explosions or neutron-star collisions, requires a robust understanding of matter at supranuclear densities. However, our knowledge about dense matter explored in the cores of neutron stars remains limited. Fortunately, dense matter is not probed only in astrophysical observations, but also in terrestrial heavy-ion collision experiments. Here we use Bayesian inference to combine data from astrophysical multi-messenger observations of neutron stars1-9 and from heavy-ion collisions of gold nuclei at relativistic energies10,11 with microscopic nuclear theory calculations12-17 to improve our understanding of dense matter. We find that the inclusion of heavy-ion collision data indicates an increase in the pressure in dense matter relative to previous analyses, shifting neutron-star radii towards larger values, consistent with recent observations by the Neutron Star Interior Composition Explorer mission5-8,18. Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and provide complementary information on nuclear matter at intermediate densities. This work combines nuclear theory, nuclear experiment and astrophysical observations, and shows how joint analyses can shed light on the properties of neutron-rich supranuclear matter over the density range probed in neutron stars.

How Perturbative QCD Constrains the Equation of State at Neutron-Star Densities

We demonstrate in a general and analytic way how high-density information about the equation of state (EOS) of strongly interacting matter obtained using perturbative quantum chromodynamics constrains the same EOS at densities reachable in physical neutron stars. Our approach is based on utilizing the full information of the thermodynamic potentials at the high-density limit together with thermodynamic stability and causality. This requires considering the pressure as a function of chemical potential p(μ) instead of the commonly used pressure as a function of energy density p(ε). The results can be used to propagate the perturbative quantum chromodynamics calculations reliable around 40n_{s} to lower densities in the most conservative way possible. We constrain the EOS starting from only a few times the nuclear saturation density n≳2.2n_{s}, and at n=5n_{s} we exclude at least 65% of otherwise allowed area in the ε-p plane. This provides information complementary to astrophysical observations that should be taken into account in any complete statistical inference study of the EOS. These purely theoretical results are independent of astrophysical neutron-star input, and hence, they can also be used to test theories of modified gravity and beyond the standard model physics in neutron stars.

Dark Matter Effects on the Compact Star Properties

The neutron star properties are generally determined by the equation of state of β-equilibrated dense matter. In this work, we consider the interaction of fermionic dark matter (DM) particles with nucleons via Higgs exchange and investigate the effect on the neutron star properties with the relativistic mean-field model equation of state coupled with DM. We deduce that DM significantly affects the neutron star properties, such as considerably reducing the maximum mass of the star, which depends on the percentage of the DM considered inside the neutron star. The tidal Love numbers both for electric and magnetic cases and surficial Love numbers are also studied for DM admixed NS. We observed that the magnitude of tidal and surficial Love numbers increases with a greater DM percentage. Further, we present post-Newtonian tidal corrections to gravitational waves decreased by increasing the DM percentage. The DM effect on the GW signal is significant during the late inspiral and merger stages of binary evolution for GW frequencies >500 Hz.

Thoughts about the utility of perturbative QCD in the cores of neutron stars – contribution to a roundtable discussion on neutron stars and QCD

In this contribution, I discuss the utility that perturbative QCD offers in studying the matter in the cores of neutron stars. I discuss the reasons why perturbative QCD can constrain the equation of state at densities far below the densities where we can perform controlled calculations. I discuss how perturbative QCD can inform nuclear modelling of neutron stars and how it influences equation-of-state inference. And finally, I discuss the implications to the QCD phase diagram and argue that interesting features in the equation of state revealed by the QCD input may be used to argue for the existence of quark-matter cores in most massive neutron stars.

Uniqueness of the Inflationary Higgs Scalar for Neutron Stars and Failure of Non-Inflationary Approximations

Neutron stars are perfect candidates to investigate the effects of a modified gravity theory, since the curvature effects are significant and more importantly, potentially testable. In most cases studied in the literature in the context of massive scalar-tensor theories, inflationary models were examined. The most important of scalar-tensor models is the Higgs model, which, depending on the values of the scalar field, can be approximated by different scalar potentials, one of which is the inflationary. Since it is not certain how large the values of the scalar field will be at the near vicinity and inside a neutron star, in this work we will answer the question, which potential form of the Higgs model is more appropriate in order for it to describe consistently a static neutron star. As we will show numerically, the non-inflationary Higgs potential, which is valid for certain values of the scalar field in the Jordan frame, leads to extremely large maximum neutron star masses; however, the model is not self-consistent, because the scalar field approximation used for the derivation of the potential, is violated both at the center and at the surface of the star. These results shows the uniqueness of the inflationary Higgs potential, since it is the only approximation for the Higgs model, that provides self-consistent results.

Phase Conversions in Neutron Stars: Implications for Stellar Stability and Gravitational Wave Astrophysics

We review the properties of hybrid stars with a quark matter core and a hadronic mantle, focusing on the role of key micro-physical properties such as the quark/hadron surface and curvature tensions and the conversion speed at the interface between both phases. We summarize the results of works that have determined the surface and curvature tensions from microscopic calculations. If these quantities are large enough, mixed phases are energetically suppressed and the quark core would be separated from the hadronic mantle by a sharp interface. If the conversion speed at the interface is slow, a new class of dynamically stable hybrid objects is possible. Densities tens of times larger than the nuclear saturation density can be attained at the center of these objects. We discuss possible formation mechanisms for the new class of hybrid stars and smoking guns for their observational identification.

Hyperons in Neutron Stars

I review the issues related to the appearance of hyperons in neutron star matter, focusing in particular on the problem of the maximum mass supported by hyperonic equations of state. I discuss the general mechanism that leads to the formation of hyperons in the core of neutron stars and I review the main techniques and many-body methods used to construct an appropriate equation of state to describe the strongly interacting system of hadrons hosted in the core of neutron stars. I outline the consequences on the structure and internal composition of neutron stars and also discuss the possible signatures of the presence of hyperons in astrophysical dynamical systems like supernova explosions and binary neutron star mergers. Finally, I briefly report about the possible important role played by hyperons in the transport properties of neutron star matter and on the consequences of neutron star cooling and gravitational wave instabilities induced by the presence of hyperons.

High-Order Multipole and Binary Love Number Universal Relations

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λl, for several high-order multipoles (l=5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M⊙ NS, R1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, RM. We find that there is an optimal value of M for every M such that the uncertainty in the estimate of RM is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.

Nuclear astrophysics in our time: supernovae, neutron stars and binary neutron star mergers

This volume focuses on topical issues in nuclear astrophysics, particularly supernovae, neutron stars, and binary neutron star mergers. An overview of those topics is presented in this editorial. The lasting impact of the Saha equation on astrophysics over the past hundred years is highlighted here.

Neutron stars harboring a primordial black hole: Maximum survival time

We explore in general relativity the survival time of neutron stars that host an endoparasitic, possibly primordial, black hole at their center. Corresponding to the minimum steady-state Bondi accretion rate for adiabatic flow that we found earlier for stiff nuclear equations of state (EOSs), we derive analytically the maximum survival time after which the entire star will be consumed by the black hole. We also show that this maximum survival time depends only weakly on the stiffness for polytropic EOSs with Γ ≥ 5/3, so that this survival time assumes a nearly universal value that depends on the initial black-hole mass alone. Establishing such a value is important for constraining the contribution of primordial black holes in the mass range 10-16 M ⊙ ≲ M ≲ 10-10 M ⊙ to the dark-matter content of the Universe.

Multimessenger Binary Mergers Containing Neutron Stars: Gravitational Waves, Jets, and γ-Ray Bursts

Neutron stars (NSs) are extraordinary not only because they are the densest form of matter in the visible Universe but also because they can generate magnetic fields ten orders of magnitude larger than those currently constructed on earth. The combination of extreme gravity with the enormous electromagnetic (EM) fields gives rise to spectacular phenomena like those observed on August 2017 with the merger of a binary neutron star system, an event that generated a gravitational wave (GW) signal, a short γ -ray burst (sGRB), and a kilonova. This event serves as the highlight so far of the era of multimessenger astronomy. In this review, we present the current state of our theoretical understanding of compact binary mergers containing NSs as gleaned from the latest general relativistic magnetohydrodynamic simulations. Such mergers can lead to events like the one on August 2017, GW170817, and its EM counterparts, GRB 170817 and AT 2017gfo. In addition to exploring the GW emission from binary black hole-neutron star and neutron star-neutron star mergers, we also focus on their counterpart EM signals. In particular, we are interested in identifying the conditions under which a relativistic jet can be launched following these mergers. Such a jet is an essential feature of most sGRB models and provides the main conduit of energy from the central object to the outer radiation regions. Jet properties, including their lifetimes and Poynting luminosities, the effects of the initial magnetic field geometries and spins of the coalescing NSs, as well as their governing equation of state, are discussed. Lastly, we present our current understanding of how the Blandford-Znajek mechanism arises from merger remnants as the trigger for launching jets, if, when and how a horizon is necessary for this mechanism, and the possibility that it can turn on in magnetized neutron ergostars, which contain ergoregions, but no horizons.

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.

Neutron Star Equation of State in Light of GW190814

The observation of gravitational waves from an asymmetric binary opens the possibility for heavy neutron stars, but these pose challenges to models of the neutron star equation of state. We construct heavy neutron stars by introducing nontrivial structure in the speed of sound sourced by deconfined QCD matter, which cannot be well recovered by spectral representations. Their moment of inertia, Love number, and quadrupole moment are very small, so a tenfold increase in sensitivity may be needed to test this possibility with gravitational waves, which is feasible with third generation detectors.

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.

Black hole-neutron star coalescence: Effects of the neutron star spin on jet launching and dynamical ejecta mass

Black hole-neutron star (BHNS) mergers are thought to be sources of gravitational waves (GWs) with coincident electromagnetic (EM) counterparts. To further probe whether these systems are viable progenitors of short gamma-ray bursts (SGRBs) and kilonovas, and how one may use (the lack of) EM counterparts associated with LIGO/Virgo candidate BHNS GW events to sharpen parameter estimation, we study the impact of neutron star spin in BHNS mergers. Using dynamical spacetime magnetohydrodynamic simulations of BHNSs initially on a quasicircular orbit, we survey configurations that differ in the BH spin (a BH/M BH = 0 and 0.75), the NS spin (a NS/M NS = -0.17, 0, 0.23, and 0.33), and the binary mass ratio (q = M BH:M NS = 3:1 and 5:1). The general trend we find is that increasing the NS prograde spin increases both the rest mass of the accretion disk onto the remnant black hole, and the rest mass of dynamically ejected matter. By a time Δt ~ 3500-5500M ~ 88-138(M NS/1.4 M ⊙) ms after the peak gravitational-wave amplitude, a magnetically driven jet is launched only for q = 3:1 regardless of the initial NS spin. The lifetime of the jets [Δt ~ 0.5-0.8(M NS/1.4 M ⊙) s] and their outgoing Poynting luminosity [L Poyn ~ 1051.5±0.5 erg/s] are consistent with typical SGRBs' luminosities and expectations from the Blandford-Znajek mechanism. By the time we terminate our simulations, we do not observe either an outflow or a large-scale magnetic-field collimation for the other systems we consider. The mass range of dynamically ejected matter is 10-4.5-10-2(M NS/1.4 M ⊙) M ⊙, which can power kilonovas with peak bolometric luminosities L knova ~ 1040-1041.4 erg/s with rise times ≲6.5 h and potentially detectable by the LSST.

Transport in Strongly Coupled Quark Matter

Motivated by the possible presence of deconfined quark matter in neutron stars and their mergers and the important role of transport phenomena in these systems, we perform the first-ever systematic study of different viscosities and conductivities of dense quark matter using the gauge/gravity duality. Using the V-QCD model, we arrive at results that are in qualitative disagreement with the predictions of perturbation theory, which highlights the differing transport properties of the system at weak and strong coupling and calls for caution in the use of the perturbative results in neutron star applications.

Constraint on Hybrid Stars with Gravitational Wave Events

Motivated by the recent discoveries of compact objects from LIGO/Virgo observations, we study the possibility of identifying some of these objects as compact stars made of dark matter called dark stars, or the mix of dark and nuclear matters called hybrid stars. In particular, in GW190814, a new compact object with 2.6 M⊙ is reported. This could be the lightest black hole, the heaviest neutron star, and a dark or hybrid star. In this work, we extend the discussion on the interpretations of the recent LIGO/Virgo events as hybrid stars made of various self-interacting dark matter (SIDM) in the isotropic limit. We pay particular attention to the saddle instability of the hybrid stars which will constrain the possible SIDM models.

Probing Crust Meltdown in Inspiraling Binary Neutron Stars

Thanks to recent measurements of tidal deformability and radius, the nuclear equation of state and structure of neutron stars are now better understood. Here, we show that through resonant tidal excitations in a binary inspiral, the neutron crust generically undergoes elastic-to-plastic transition, which leads to crust heating and eventually meltdown. This process could induce ∼O(0.1) phase shift in the gravitational waveform. Detecting the timing and induced phase shift of this crust meltdown will shed light on the crust structure, such as the core-crust transition density, which previous measurements are insensitive to. A direct search using GW170817 data has not found this signal, possibly due to limited signal-to-noise ratio. We predict that such a signal may be observable with Advanced LIGO Plus and more likely with third-generation gravitational-wave detectors such as the Einstein Telescope and Cosmic Explorer.

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.

A lower bound on the maximum mass if the secondary in GW190814 was once a rapidly spinning neutron star

The recent detection of GW190814 featured the merger of a binary with a primary having a mass of $\sim 23\, \mathrm{ M}_{\odot }$ and a secondary with a mass of $\sim 2.6\, \mathrm{ 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 although the secondary in GW190814 is most likely a black hole at merger, it 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_{-0.05}^{+0.08} \lesssim \chi \lesssim 0.68_{-0.05}^{+0.11}$, of the secondary – a quantity not constrained so far by the detection – and a novel strict lower bound on the maximum mass, $M_{_{\mathrm{TOV}}}\gt 2.08^{+0.04}_{-0.04}\, \, \mathrm{ M}_{\odot }$ and an optimal bound of $M_{_{\mathrm{TOV}}}\gt 2.15^{+0.04}_{-0.04}\, \, \mathrm{ M}_{\odot }$, of non-rotating 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 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.

Nuclear Pairing Gaps and Neutron Star Cooling

We study the cooling of isolated neutron stars with particular regard to the importance of nuclear pairing gaps. A microscopic nuclear equation of state derived in the Brueckner-Hartree-Fock approach is used together with compatible neutron and proton pairing gaps. We then study the effect of modifying the gaps on the final deduced neutron star mass distributions. We find that a consistent description of all current cooling data can be achieved and a reasonable neutron star mass distribution can be predicted employing the (slightly reduced by about 40%) proton 1S0 Bardeen-Cooper-Schrieffer (BCS) gaps and no neutron 3P2 pairing.

General Relativistic Mean-Field Dynamo Model for Proto-Neutron Stars

Neutron stars, and magnetars in particular, are known to host the strongest magnetic fields in the Universe. The origin of these strong fields is a matter of controversy. In this preliminary work, via numerical simulations, we study, for the first time in non-ideal general relativistic magnetohydrodynamic (GRMHD) regime, the growth of the magnetic field due to the action of the mean-field dynamo due to sub-scale, unresolved turbulence. The dynamo process, combined with the differential rotation of the (proto-)star, is able to produce an exponential growth of any initial magnetic seed field up to the values required to explain the observations. By varying the dynamo coefficient we obtain different growth rates. We find a quasi-linear dependence of the growth rates on the intensity of the dynamo. Furthermore, the time interval in which exponential growth occurs and the growth rates also seems to depend on the initial configuration of the magnetic field.

Postmerger Gravitational-Wave Signatures of Phase Transitions in Binary Mergers

With the first detection of gravitational waves from a binary system of neutron stars GW170817, a new window was opened to study the properties of matter at and above nuclear-saturation density. Reaching densities a few times that of nuclear matter and temperatures up to 100 MeV, such mergers also represent potential sites for a phase transition (PT) from confined hadronic matter to deconfined quark matter. While the lack of a postmerger signal in GW170817 has prevented us from assessing experimentally this scenario, two theoretical studies have explored the postmerger gravitational-wave signatures of PTs in mergers of a binary system of neutron stars. We here extend and complete the picture by presenting a novel signature of the occurrence of a PT. More specifically, using fully general-relativistic hydrodynamic simulations and employing a suitably constructed equation of state that includes a PT, we present the occurrence of a "delayed PT," i.e., a PT that develops only some time after the merger and produces a metastable object with a quark-matter core, i.e., a hypermassive hybrid star. Because in this scenario, the postmerger signal exhibits two distinct fundamental gravitational-wave frequencies-before and after the PT-the associated signature promises to be the strongest and cleanest among those considered so far, and one of the best signatures of the production of quark matter in the present Universe.