Functional forms of the neutron star Equation of State (EoS) are required to extract the viable EoS band from neutron star observations. Realistic nuclear EoSs, containing deconfined quarks or hyperons, present nontrivial features in the speed of sound such as bumps, kinks and plateaus. Using modified Gaussian processes to model EoSs with nontrivial features, we show in a fully Bayesian analysis incorporating measurements from X-ray sources, gravitational wave observations, and perturbative QCD results that these features are compatible with current constraints. We find nontrivial behavior in the EoS plays a role in understanding the possible phase structure of neutron stars at densities around twice saturation density [1]. Lastly, we perform a large-scale systematic analysis of the impact of perturbative QCD constraints when they are applied beyond the maximum central densities realized in realistic neutron star EoSs.
[1] D. Mroczek et al, arXiv:2302.07978