Extracting the Neutron Star Equation of State from Theoretical and Observational Constraints Using Modified Gaussian Processes

Functional forms of the neutron star Equation of State (EoS) are required to extract the viable EoS band from neutron star observations. Typically, one of three methods are used-- spectral functions, piecewise polytropes, or gaussian processes. However, realistic nuclear EoS, containing deconfined quarks or hyperons, present nontrivial features in the speed of sound such as bumps, kinks and plateaus. These features in the speed of sound are difficult to capture using previous methods [1]. We modify gaussian processes by introducing spikes and plateaus in the speed of sound and check how those features impact the posterior distribution obtained using constraints from NICER, gravitational-wave observations, and perturbative QCD. We find that introducing modifications to the EoS plays a role in understanding the possible phase structure of neutron stars at densities around twice saturation density. 

[1] Tan et al, Phys.Rev.Lett. 125 (2020) 26, 261104