Investigating Twin Star Equation of States in Light of Recent Astrophysical Observations

TL;DR

This study explores the conditions and observational constraints for twin star configurations in hybrid equations of state, analyzing how phase transition parameters influence maximum mass and structural features of such stars.

Contribution

It provides a comprehensive parameter space analysis for twin star equations of state, incorporating recent astrophysical observations to constrain phase transition properties.

Findings

Maximum twin star mass is 2.05 solar masses.

Allowed strongest density discontinuity is 7.76 times nuclear saturation density.

Upper bound on transition density is 4.03 times nuclear saturation density.

Abstract

Twin stars are predicted to exist in nature if the hadron-to-quark phase transition is strong enough to form a new branch of hybrid stars, separated from the branch of neutron stars. We adopt an agnostic approach, using transition energy density, transition pressure, the discontinuity strength, and a constant speed of sound for quark matter as our parameter space to construct a large possibility of hybrid equations of state, and thereby encapsulating a comprehensive picture of the twin star scenario. First, we report the complete conditions on our parameter space imposed by the general relativistic hydrostatic equilibrium solutions. For a fixed transition energy density and speed of sound for quark matter, we define distinct ranges of transition pressures based on the allowed strengths of discontinuity. Below a maximum transition pressure, a range of discontinuity exists that increases as the transition pressure decreases. Thereby, we identify the loci of the limits on discontinuities as the `witch-hat' curves. Based on the causality limit, the witch-hat curves can be punctured or incomplete. Strong constraints on this picture are drawn from the inferences from GW170817 and the NICER measurements. We computed the maximum mass for twin stars to be $2.05~M_\odot$, the allowed strongest discontinuity in rest-mass density to be $7.76\rho_\mathrm{sat}$, and the upper bound on transition rest-mass density to be $4.03\rho_\mathrm{sat}$. Subsequently, we compute the implications of the stiffness of the quark matter equation of state on this picture. Different confidence levels for observational inferences are considered to assess the extent of inclusion (and rejection) of hybrid equations of state and, consequently, their effects on the limits of the maximum mass of twin stars and phase transition properties.