Representing Equations of State with Strong First-Order Phase Transitions

Parametric representations of the high-density nuclear equation of state are used in constructing models for interpreting the astrophysical observations of neutron stars. This study explores how accurately equations of state with strong first-order phase transitions can be represented using spectral or piecewise analytic methods that assume no a priori knowledge of the location or the strength of the phase transition. The model equations of state used in this study have phase transitions strong enough to induce a gravitational instability that terminates the sequence of stable neutron stars. These equations of state also admit a second sequence of stable stars with core matter that has undergone this strong first-order phase transition (possibly driven by quark deconfinement). These results indicate that spectral representations generally achieve somewhat higher accuracy than piecewise analytic representations having the same number of parameters. Both types of representation show power-law convergence at approximately the same rate.