Abstract:The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ together, a discrepancy between model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95%), respectively. Our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

Bayesian uncertainty quantification for nuclear matter incompressibility

Bayesian inference of nuclear incompressibility from proton elliptic flow in central Au+Au collisions at 400 MeV/nucleon

Uncertainty Quantification for Nuclear Density Functional Theory and Information Content of New Measurements

A Bayesian mixture model approach to quantifying the empirical nuclear saturation point

Constraining isovector nuclear interactions with giant resonances within a Bayesian approach

Bayesian Inference on Isospin Splitting of Nucleon Effective Mass from Giant Resonances in 208 Pb

Uncertainty quantification in $(p,n)$ reactions

Uncertainty evaluation and correlation analysis of single-particle energies in phenomenological nuclear mean field: An investigation of propagating uncertainties for independent model parameters

Effect of nuclear matter incompressibility on the 16O+208Pb system

Microscopic calculations of nuclear incompressibility from cluster radioactivity

Bayesian Inference of Fine-Features of Nuclear Equation of State from Future Neutron Star Radius Measurements to 0.1km Accuracy

Uncertainty evaluation and correlation analysis of single-particle energies in phenomenological nuclear mean field: an investigation into propagating uncertainties for independent model parameters

Effects of incompressibility $K_{0}$ in heavy-ion collisions at intermediate energies

Bayesian inference of the dense matter equation of state built upon extended Skyrme interactions

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

Quantified limits of the nuclear landscape

Nambu-Jona-Lasinio description of hadronic matter from a Bayesian approach

Structure effects on the giant monopole resonance and determinations of the nuclear incompressibility

Bayesian Survey of the Dense Matter Equation of State Built upon Skyrme Effective Interactions