Bayesian Quantification of Observability and Equation of State of Twin Stars

TL;DR

This paper uses Bayesian methods with neutron star radius data to quantify the likelihood and properties of twin stars, revealing their robustness and potential to constrain the neutron star equation of state.

Contribution

It introduces a Bayesian framework with a flexible meta-model to analyze the observability and EOS of twin stars using current observational data.

Findings

12-18% of EOSs produce twin stars

Twin star scenarios are observable with current measurement accuracy

EOS parameters depend on twin star category and presence

Abstract

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equation of state (EOS) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a meta-model with 9-parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12-18\% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. We show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared $c_{\rm qm}^2$ has two peaks, one below and another above the conformal limit $c_{\rm qm}^2=1/3$ predicted by perturbative QCD.