Reuse & Permissions

It is not necessary to obtain permission to reuse this article or its components as it is available under the terms of the Creative Commons Attribution 4.0 International license. This license permits unrestricted use, distribution, and reproduction in any medium, provided attribution to the author(s) and the published article's title, journal citation, and DOI are maintained. Please note that some figures may have been included with permission from other third parties. It is your responsibility to obtain the proper permission from the rights holder directly for these figures.

Export citation

Export citation

Choose format for download:

Download Citation
  • Open Access

Emulator-assisted nuclear density-functional-theory inference and its consequences for the structure of neutron stars

Pietro Klausner1,2,3,*, Marco Antonelli1,†, Gianluca Colò2,3,‡, Francesca Gulminelli1,4,§, Xavier Roca-Maza2,3,5,6,∥, and Enrico Vigezzi3,¶

  • *Contact author: pietro.klausner@unimi.it
  • Contact author: antonelli@lpccaen.in2p3.fr
  • Contact author: gianluca.colo@mi.infn.it
  • §Contact author: gulminelli@lpccaen.in2p3.fr
  • Contact author: xavier.roca.maza@fqa.ub.es
  • Contact author: enrico.vigezzi@mi.infn.it

Phys. Rev. C 114, 034320 – Published 16 September, 2026

DOI: https://doi.org/10.1103/p6xt-x2zp

Abstract

Nuclear density functional theory provides a unified description of finite nuclei and bulk nuclear matter, and is widely used to model the neutron star (NS) equation of state. However, extrapolations to suprasaturation densities require a quantified treatment of uncertainties arising from parameter estimation and functional choices. We present an updated Bayesian inference of a Skyrme energy density functional augmented by a flexible metamodel density dependence at high density. Nuclear observables are computed using a Gaussian emulator of the publicly available Milano hfbcs-qrpa code, enabling efficient exploration of a high-dimensional parameter space. Relative to previous analyses, we extend the calibration set with isospin-sensitive data, including masses and charge radii along selected Ca and Sn isotopic chains, and updated constraints from giant monopole resonances. The resulting posteriors are further constrained by ab initio neutron-matter calculations and astrophysical observations, including recent NICER measurements, yielding consistent crust and core properties of catalyzed NS compatible with current constraints. Bulk nuclear-matter parameters are well approximated by a multivariate Gaussian with covariance matrix provided for direct reuse, while several finite-nucleus parameters exhibit pronounced non-Gaussianity.

View figure in article

Physics Subject Headings (PhySH)

Article Text

Supplemental Material

References (70)

  1. N. Schunck, Energy Density Functional Methods for Atomic Nuclei (IOP Publishing, Bristol, 2019), pp. 2053–2563.
  2. G. Grams, N. N. Shchechilin, A. Sánchez-Fernández, W. Ryssens, N. Chamel, and S. Goriely, Skyrme–Hartree–Fock–Bogoliubov mass models on a 3D mesh: IV. Improved description of the isospin dependence of pairing, Eur. Phys. J. A 61, 35 (2025).
  3. G. Grams, W. Ryssens, A. Sánchez-Fernández, N. N. Shchechilin, L. G.-M. Zaragoza, P. Demol, N. Chamel, S. Goriely, and M. Bender, Skyrme-Hartree-Fock-Bogoliubov mass models on a 3D mesh V. The N2LO extension of the Skyrme EDF, Phys. Lett. B 879, 140590 (2026).
  4. L. Scurto, H. Pais, and F. Gulminelli, General predictions of neutron star properties using unified relativistic mean-field equations of state, Phys. Rev. D 109, 103015 (2024).
  5. P. Char and C. Mondal, Exploring the limits of nucleonic metamodeling using different relativistic density functionals, Phys. Rev. D 111, 103024 (2025).
  6. H. Dinh Thi, C. Mondal, and F. Gulminelli, The nuclear matter density functional under the nucleonic hypothesis, Universe 7, 373 (2021).
  7. M. V. Beznogov and A. R. Raduta, Bayesian inference of the dense matter equation of state built upon extended Skyrme interactions, Phys. Rev. C 110, 035805 (2024).
  8. X. Roca-Maza and N. Paar, Nuclear equation of state from ground and collective excited state properties of nuclei, Prog. Part. Nucl. Phys. 101, 96 (2018).
  9. R. Kumar, V. Dexheimer, J. Jahan, J. Noronha, J. Noronha-Hostler, C. Ratti, N. Yunes, A. R. Nava Acuna, M. Alford, M. H. Anik, D. Chatterjee, K. Chatziioannou, H.-Y. Chen, A. Clevinger, C. Conde, N. Cruz-Camacho, T. Dore, C. Drischler, H. Elfner, R. Essick, et al., Theoretical and experimental constraints for the equation of state of dense and hot matter, Living Rev. Relativity 27, 3 (2024).
  10. E. Yüksel, T. Marketin, and N. Paar, Optimizing the relativistic energy density functional with nuclear ground state and collective excitation properties Phys. Rev. C 99, 034318 (2019).
  11. Z. Zhang and L.-W. Chen, Bayesian inference of the symmetry energy and the neutron skin in 48Ca and 208Pb from CREX and PREX-2, Phys. Rev. C 108, 024317 (2023).
  12. J. Xu and P. Papakonstantinou, Bayesian inference of finite-nuclei observables based on the KIDS model Phys. Rev. C 105, 044305 (2022).
  13. J. Erler, N. Birge, M. Kortelainen, W. Nazarewicz, E. Olsen, A. M. Perhac, and M. Stoitsov, The limits of the nuclear landscape, Nature (London) 486, 509 (2012).
  14. P. Klausner, G. Colò, X. Roca-Maza, and E. Vigezzi, Impact of ground-state properties and collective excitations on the Skyrme ansatz: A Bayesian study, Phys. Rev. C 111, 014311 (2025).
  15. P. Klausner, M. Antonelli, and F. Gulminelli, Properties of the neutron star crust informed by nuclear structure data, Phys. Rev. C 113, 025808 (2026).
  16. C. E. Rasmussen and C. K. I. Williams, Gaussian Processes for Machine Learning (MIT Press, Cambridge, 2006).
  17. MADAI collaboration (Models and D. A. Initiative), http://madai.phy.duke.edu/.
  18. J. Margueron, R. Hoffmann Casali, and F. Gulminelli, Equation of state for dense nucleonic matter from metamodeling. I. Foundational aspects, Phys. Rev. C 97, 025805 (2018).
  19. P. Landry and R. Essick, Nonparametric inference of the neutron star equation of state from gravitational wave observations, Phys. Rev. D 99, 084049 (2019).
  20. R. Essick, P. Landry, and D. E. Holz, Nonparametric inference of neutron star composition, equation of state, and maximum mass with GW170817, Phys. Rev. D 101, 063007 (2020).
  21. J. S. Read, B. D. Lackey, B. J. Owen, and J. L. Friedman, Constraints on a phenomenologically parametrized neutron-star equation of state, Phys. Rev. D 79, 124032 (2009).
  22. L. Suleiman, M. Fortin, J. L. Zdunik, and C. M. C. Providência, Polytropic fits of modern and unified equations of state, Phys. Rev. C 106, 035805 (2022).
  23. E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and R. Schaeffer, A Skyrme parametrization from subnuclear to neutron star densities, Nucl. Phys. A 627, 710 (1997).
  24. L.-W. Chen, B.-J. Cai, C. M. Ko, B.-A. Li, C. Shen, and J. Xu, Higher-order effects on the incompressibility of isospin asymmetric nuclear matter, Phys. Rev. C 80, 014322 (2009).
  25. L.-W. Chen, C. M. Ko, B.-A. Li, and J. Xu, Density slope of the nuclear symmetry energy from the neutron skin thickness of heavy nuclei Phys. Rev. C 82, 024321 (2010).
  26. P. Ring, P. Schuck, and M. R. Strayer, The nuclear many-body problem, Phys. Today 36(7), 70 (1983).
  27. W. Satuła and R. Wyss, Competition between T=0 and T=1 pairing in proton-rich nuclei, Phys. Lett. B 393, 1 (1997).
  28. W. Satuła, D. Dean, J. Gary, S. Mizutori, and W. Nazarewicz, On the origin of the Wigner energy, Phys. Lett. B 407, 103 (1997).
  29. W. Satuła and R. Wyss, A number projected model with generalized pairing interaction, Nucl. Phys. A 676, 120 (2000).
  30. M. Wang, W. Huang, F. Kondev, G. Audi, and S. Naimi, The AME 2020 atomic mass evaluation (II). Tables, graphs and references, Chin. Phys. C 45, 030003 (2021).
  31. I. Angeli and K. Marinova, Table of experimental nuclear ground state charge radii: An update, At. Data Nucl. Data Tables 99, 69 (2013).
  32. W. Satuła, J. Dobaczewski, and W. Nazarewicz, Odd-even staggering of nuclear masses: Pairing or shape effect? Phys. Rev. Lett. 81, 3599 (1998).
  33. Y. K. Gupta, K. B. Howard, U. Garg, J. T. Matta, M. Şenyiğit, M. Itoh, S. Ando, T. Aoki, A. Uchiyama, S. Adachi, M. Fujiwara, C. Iwamoto, A. Tamii, H. Akimune, C. Kadono, Y. Matsuda, T. Nakahara, T. Furuno, T. Kawabata, M. Tsumura, et al., Isoscalar giant monopole, dipole, and quadrupole resonances in Zr90,92 and Mo92, Phys. Rev. C 97, 064323 (2018).
  34. Y. Gupta, U. Garg, K. Howard, J. Matta, M. Şenyiğit, M. Itoh, S. Ando, T. Aoki, A. Uchiyama, S. Adachi, M. Fujiwara, C. Iwamoto, A. Tamii, H. Akimune, C. Kadono, Y. Matsuda, T. Nakahara, T. Furuno, T. Kawabata, M. Tsumura, et al., Are there nuclear structure effects on the isoscalar giant monopole resonance and nuclear incompressibility near A90, Phys. Lett. B 760, 482 (2016).
  35. J. A. Melendez, R. J. Furnstahl, D. R. Phillips, M. T. Pratola, and S. Wesolowski, Quantifying correlated truncation errors in effective field theory, Phys. Rev. C 100, 044001 (2019).
  36. G. Colò, L. Cao, N. Van Giai, and L. Capelli, Self-consistent RPA calculations with Skyrme-type interactions: The skyrme_rpa program, Comput. Phys. Commun. 184, 142 (2013).
  37. G. Colò and X. Roca-Maza, User guide for the hfbcs-qrpa(v1) code, arXiv:2102.06562.
  38. M. Kortelainen, T. Lesinski, J. Moré, W. Nazarewicz, J. Sarich, N. Schunck, M. V. Stoitsov, and S. Wild, Nuclear energy density optimization Phys. Rev. C 82, 024313 (2010).
  39. P. Möller, W. D. Myers, H. Sagawa, and S. Yoshida, New finite-range droplet mass model and equation-of-state parameters Phys. Rev. Lett. 108, 052501 (2012).
  40. S. Shlomo, V. M. Kolomietz, and G. Colò, Deducing the nuclear-matter incompressibility coefficient from data on isoscalar compression modes, Eur. Phys. J. A 30, 23 (2006).
  41. E. Khan, J. Margueron, and I. Vidaña, Constraining the nuclear equation of state at subsaturation densities Phys. Rev. Lett. 109, 092501 (2012).
  42. B.-A. Li, B.-J. Cai, W.-J. Xie, and N.-B. Zhang, Progress in constraining nuclear symmetry energy using neutron star observables since GW170817, Universe 7, 182 (2021).
  43. P.-G. Reinhard, X. Roca-Maza, and W. Nazarewicz, Combined theoretical analysis of the parity-violating asymmetry for Ca48 and Pb208, Phys. Rev. Lett. 129, 232501 (2022).
  44. X. Roca-Maza, M. Brenna, G. Colò, M. Centelles, X. Viñas, B. K. Agrawal, N. Paar, D. Vretenar, and J. Piekarewicz, Electric dipole polarizability in Pb208: Insights from the droplet model, Phys. Rev. C 88, 024316 (2013).
  45. X. Roca-Maza, X. Viñas, M. Centelles, B. K. Agrawal, G. Colò, N. Paar, J. Piekarewicz, and D. Vretenar, Neutron skin thickness from the measured electric dipole polarizability in Ni68, Sn120, and Pb208, Phys. Rev. C 92, 064304 (2015).
  46. T.-G. Yue, Z. Zhang, and L.-W. Chen, PREX and CREX: Evidence for strong isovector spin-orbit interaction, Sci. Bull. 71, 1270 (2026).
  47. T. Zhao, Z. Lin, B. Kumar, A. W. Steiner, and M. Prakash, Characterizing the nuclear models informed by PREX and CREX: A view from Bayesian inference, Phys. Rev. Res. 7, 043335 (2025).
  48. M. Qiu, T.-G. Yue, Z. Zhang, and L.-W. Chen, A relativistic mechanism for the enhanced isovector spin-orbit interaction suggested by parity-violating electron scattering experiments, Phys. Lett. B 879, 140648 (2026).
  49. A. Kunjipurayil, J. Piekarewicz, and M. Salinas, Role of the isovector spin-orbit potential in mitigating the CREX-PREX dilemma, Phys. Rev. C 112, 014310 (2025).
  50. M. Salinas and J. Piekarewicz, Impact of tensor couplings with scalar mixing on covariant energy density functionals, Phys. Rev. C 109, 045807 (2024).
  51. P. Papakonstantinou, CREX and PREX-II reconciled within energy-density functional theory, arXiv:2602.22722.
  52. G. Montefusco, M. Antonelli, and F. Gulminelli, Frozen and β-equilibrated f and p modes of cold neutron stars: Nuclear metamodel predictions A& A 694, A150 (2025).
  53. A. Saffer, E. Fonseca, S. Ransom, I. Stairs, R. Lynch, D. Good, K. W. Masui, J. W. McKee, B. W. Meyers, S. S. Patil, and C. M. Tan, A lower mass estimate for PSR J0348+0432 based on CHIME/Pulsar precision timing, Astrophys. J. Lett. 983, L20 (2025).
  54. J. Antoniadis, P. C. C. Freire, N. Wex, T. M. Tauris, R. S. Lynch, M. H. van Kerkwijk, M. Kramer, C. Bassa, V. S. Dhillon, T. Driebe, J. W. T. Hessels, V. M. Kaspi, V. I. Kondratiev, N. Langer, T. R. Marsh, M. A. McLaughlin, T. T. Pennucci, S. M. Ransom, I. H. Stairs, J. van Leeuwen, et al., A massive pulsar in a compact relativistic binary, Science 340, 1233232 (2013).
  55. E. Fonseca, H. T. Cromartie, T. T. Pennucci, P. S. Ray, A. Y. Kirichenko, S. M. Ransom, P. B. Demorest, I. H. Stairs, Z. Arzoumanian, L. Guillemot, A. Parthasarathy, M. Kerr, I. Cognard, P. T. Baker, H. Blumer, P. R. Brook, M. DeCesar, T. Dolch, F. A. Dong, E. C. Ferrara, et al., Refined mass and geometric measurements of the high-mass PSR J0740+6620, Astrophys. J. Lett. 915, L12 (2021).
  56. J. M. Pearson, N. Chamel, A. Y. Potekhin, A. F. Fantina, C. Ducoin, A. K. Dutta, and S. Goriely, Unified equations of state for cold non-accreting neutron stars with Brussels–Montreal functionals – I. Role of symmetry energy Mon. Not. R. Astron. Soc. 481, 2994 (2018).
  57. X. Roca-Maza and J. Piekarewicz, Impact of the symmetry energy on the outer crust of nonaccreting neutron stars, Phys. Rev. C 78, 025807 (2008).
  58. H. Sotani, Shear oscillations in neutron stars and the nuclear symmetry energy Phys. Rev. D 109, 023030 (2024).
  59. N. A. Zemlyakov and A. I. Chugunov, Neutron star inner crust: Reduction of shear modulus by nuclei finite size effect, Mon. Not. R. Astron. Soc. 518, 3813 (2023).
  60. N. A. Zemlyakov and A. I. Chugunov, Constraining the shear modulus of a polycrystalline neutron star crust: Hashin-Shtrikman variational approach, Phys. Rev. D 112, 043032 (2025).
  61. M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Arzoumanian, K. C. Gendreau, S. Guillot, W. C. G. Ho, J. M. Lattimer, M. Loewenstein, S. M. Morsink, P. S. Ray, M. T. Wolff, C. L. Baker, T. Cazeau, S. Manthripragada, C. B. Markwardt, T. Okajima, S. Pollard, I. Cognard, et al., The radius of PSR J0740+6620 from NICER and XMM-Newton data, Astrophys. J. Lett. 918, L28 (2021).
  62. D. Choudhury, T. Salmi, S. Vinciguerra, T. E. Riley, Y. Kini, A. L. Watts, B. Dorsman, S. Bogdanov, S. Guillot, P. S. Ray, D. J. Reardon, R. A. Remillard, A. V. Bilous, D. Huppenkothen, J. M. Lattimer, N. Rutherford, Z. Arzoumanian, K. C. Gendreau, S. M. Morsink, and W. C. G. Ho, A NICER view of the nearest and brightest millisecond pulsar: PSR J0437–4715, Astrophys. J. Lett. 971, L20 (2024).
  63. S. Vinciguerra, T. Salmi, A. L. Watts, D. Choudhury, T. E. Riley, P. S. Ray, S. Bogdanov, Y. Kini, S. Guillot, D. Chakrabarty, W. C. G. Ho, D. Huppenkothen, S. M. Morsink, Z. Wadiasingh, and M. T. Wolff, An updated mass–radius analysis of the 2017–2018 NICER data set of PSR J0030+0451, Astrophys. J. 961, 62 (2024).
  64. L. Mauviard, S. Guillot, T. Salmi, D. Choudhury, B. Dorsman, D. González-Caniulef, M. Hoogkamer, D. Huppenkothen, C. Kazantsev, Y. Kini, J.-F. Olive, P. Stammler, A. L. Watts, M. Mendes, N. Rutherford, A. Schwenk, I. Svensson, S. Bogdanov, M. Kerr, P. S. Ray, et al., A NICER view of the 1.4M edge-on pulsar PSR J0614-3329, Astrophys. J. 995, 60 (2025).
  65. See Supplemental Material at http://link.aps.org/supplemental/10.1103/p6xt-x2zp for files containing numerical values of Tables 5 and 6.
  66. N. Rutherford, M. Mendes, I. Svensson, A. Schwenk, A. L. Watts, K. Hebeler, J. Keller, C. Prescod-Weinstein, D. Choudhury, G. Raaijmakers, T. Salmi, P. Timmerman, S. Vinciguerra, S. Guillot, and J. M. Lattimer, Constraining the dense matter equation of state with new NICER mass–radius measurements and new chiral effective field theory inputs, Astrophys. J. Lett. 971, L19 (2024).
  67. H. Gil, Y. Oh, C. H. Hyun, and P. Papakonstantinou, Skyrme-type nuclear force for the KIDS energy density functional, New Phys.: Sae Mulli 67, 456 (2017).
  68. H. Gil, P. Papakonstantinou, C. H. Hyun, T. S. Park, and Y. Oh, Acta Phys. Pol. B 48, 305 (2017).
  69. P. Papakonstantinou, T.-S. Park, Y. Lim, and C. H. Hyun, Density dependence of the nuclear energy-density functional Phys. Rev. C 97, 014312 (2018).
  70. S. A. Fayans, Towards a universal nuclear density functional, Pis'ma Zh. Éksp. Teor. Fiz. 68, 161 (1998) [JETP Lett. 68, 169 (1998)].

Outline

Information

Sign In to Your Journals Account

Filter

Filter

Article Lookup

Enter a citation