Neutron-quark stars: Discerning viable alternatives for the higher-density part of the equation of state of compact stars

TL;DR

This paper explores different models of neutron stars with quark matter, using simulations to identify gravitational wave signatures that distinguish between hadronic and quark-including equations of state, including scenarios with phase transitions.

Contribution

It introduces the concept of neutron-quark stars (NQS) and demonstrates how gravitational wave observations can differentiate between hadronic and quark matter models in neutron star mergers.

Findings

QHC models can be distinguished from hadronic models via gravitational wave signals.

The dominant postmerger frequency ($f_2$) is lower in QHC models for the same tidal deformability.

A linear correlation exists between postmerger GW energy and $f_2$ frequency, independent of quark matter presence.

Abstract

We investigate binary neutron star (BNS) mergers using general-relativistic numerical simulations with hadronic and hybrid equations of state (EOSs), incorporating the latest observations and theoretical constraints. We address two viable scenarios for the transition to quark matter: a quark-hadron crossover (QHC) or a strong first-order phase transition (1PT). To distinguish between different models, we define neutron-quark stars (NQS) as configurations where quark effects emerge at masses below the lowest observed neutron-star mass. While traditional "hybrid stars" may be distinguished by purely hadronic configurations through mass-radius measurements, the mass-radius relations of NQSs resemble those of purely hadronic models, with no sharp boundary between hadrons and quarks. The name NQS effectively captures the absence of a phase boundary between hadrons and quarks in QHC scenarios. Our results indicate that QHC models can be distinguished from hadronic ones if both the inspiral and postmerger gravitational waves (GWs) are observed. In particular, the dominant postmerger frequency ($f_2$) tends to be lower than in hadronic models with the same tidal deformability ($\Lambda$). We also present the first general-relativistic simulations of BNS mergers where the stars already contain quark matter before merging. These involve a strong first-order phase transition (1PT) at 1.8 times nuclear saturation density, followed by a stiff quark EOS. Finally, we identify a robust linear correlation between the total GW energy emitted after the merger and the $f_2$ frequency. Remarkably, this relation holds regardless of the quark presence.