Jiarui Chen, Zuhua Ji

Understanding the equation of state (EOS) of neutron stars (NSs) is a fundamental challenge in astrophysics and nuclear physics. A first-order phase transition (FOPT) at high densities could lead to the formation of a quark core, significantly affecting NS properties. This review explores observational and theoretical constraints on such transitions using multi-messenger astrophysics. X-ray observations, including mass-radius measurements from NICER and spectral features like quasi-periodic oscillations (QPOs) and cyclotron resonance scattering features (CRSFs), provide indirect evidence of EOS modifications. Gravitational wave detections, particularly from binary NS mergers such as GW170817, constrain tidal deformability and post-merger oscillations, which may carry signatures of phase transitions. Pulsar timing offers additional constraints through measurements of mass, spin evolution, and glitches, with millisecond pulsars exceeding twice the solar mass posing challenges to purely hadronic EOSs. Theoretical models and numerical simulations predict that an FOPT could impact gravitational wave signals, twin-star configurations, and NS cooling. Future advancements, including next-generation gravitational wave detectors, high-precision X-ray telescopes, and improved theoretical modeling, will enhance our ability to probe phase transitions in NSs. A combination of these approaches will provide crucial insights into the existence and properties of deconfined quark matter in NS interiors.