Published September 2, 2025 | Published
Journal Article Open

Nonparametric extensions of nuclear equations of state: Probing the breakdown scale of relativistic mean-field theory

  • 1. ROR icon California Institute of Technology
  • 2. ROR icon Washington University in St. Louis
  • 3. ROR icon University of Southampton
  • 4. ROR icon Canadian Institute for Theoretical Astrophysics
  • 5. ROR icon University of Toronto

Abstract

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. We further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

Copyright and License

© 2025 American Physical Society.

Acknowledgement

We thank Thibeau Wouters for helpful comments on this manuscript. L. B. and A. H. thank Mark Alford for useful conversations. I. L. thanks Philippe Landry for insightful discussions. I. L. and K. C. acknowledge support from the Department of Energy under Award No. DE-SC0023101, and by a grant from the Simons Foundation (MP-SCMPS-00001470). L. B. and A. H. are partly supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, under Award No. DE-FG02-05ER41375. A. H. furthermore acknowledges financial support by the UKRI under the Horizon Europe Guarantee project EP/Z000939/1. R. E. is supported by the Natural Sciences and Engineering Research Council of Canada (NSERC) through a Discovery Grant (RGPIN-2023-03346). The authors are grateful for computational resources provided by the LIGO Laboratory and supported by National Science Foundation Grants No. PHY-0757058 and No. PHY-0823459. This material is based upon work supported by NSF’s LIGO Laboratory which is a major facility fully funded by the National Science Foundation.

Data Availability

The data that support the findings of this article are openly available [104].

Files

9kh9-xfpd.pdf
Files (10.5 MB)
Name Size Download all
md5:a64134b143c30f90b81a8ce3e8dfdd88
10.5 MB Preview Download

Additional details

Created:
September 4, 2025
Modified:
September 4, 2025