Toward constraining QCD phase transitions in neutron star interiors: Bayesian inference with a Tolman-Oppenheimer-Volkof linear response analysis

The potential hadron-to-quark phase transition in neutron stars has not been fully understood as the property of cold, dense, and strongly interacting matter cannot be theoretically described by the first-principle perturbative calculations, nor have they been systematically measured through terrestrial low-to-intermediate energy heavy-ion experiments. Given the Tolman-Oppenheimer-Volkoff (TOV) equations, the equation of state (EoS) of the neutron star (NS) matter can be constrained by the observations of NS mass, radius, and tidal deformability. However, large observational uncertainties and the limited number of observations currently make it challenging to strictly reconstruct the EoS, especially to identify interesting features such as a strong first-order phase transition. In this work, we study the dependency of reconstruction quality of the phase transition on the number of NS observations of mass and radius as well as their uncertainty, based on a fiducial EoS. We conquer this challenging problem by constructing a neural network, which allows one to parametrize the EoS with minimum model-dependency, and by devising an algorithm of parameter optimization based on the analytical linear response analysis of the TOV equations. This work may pave the way for the understanding of the phase transition features in NSs using future x-ray and gravitational wave measurements.