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 mixture model approach to quantifying the empirical nuclear saturation point

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

Bayesian Survey of the Dense Matter Equation of State built upon Skyrme effective interactions

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

QCD Equation of State of Dense Nuclear Matter from a Bayesian Analysis of Heavy-Ion Collision Data

From chiral EFT to perturbative QCD: a Bayesian model mixing approach to symmetric nuclear matter

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

Bayesian inference of the dense-matter equation of state of neutron stars with antikaon condensation

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

Bayesian uncertainty quantification for nuclear matter incompressibility

Simultaneously Constraining the Neutron Star Equation of State and Mass Distribution through Multimessenger Observations and Nuclear Benchmarks

Nucleon-quark mixed matter and neutron-star equation of state

Bayesian Inference of dense matter equation of state of neutron star with antikaon condensation

Quantum van der Waals theory meets quarkyonic matter

Quantified limits of the nuclear landscape

Nuclear incompressibility and speed of sound in uniform matter and finite nuclei

Non-Parametric Multi-Messenger Constraints on the Equation of State of Cold Dense Nuclear Matter

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