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.

The improved saturation model in nuclei

Color dipole cross-section from unifying the color dipole picture and improved saturation models

Suppression of diffraction in deep-inelastic scattering on nuclei and dynamical mechanism of leading twist nuclear shadowing

The dipole cross section by the unintegrated gluon distribution at small $x$

Reduced cross-section in Electron-Ion Colliders at small $x$

CGC/saturation approach: Impact-parameter dependent model in next-to-leading order of perturbative QCD and combined HERA data

Dynamical model of antishadowing of the nuclear gluon distribution

Shadowing and antishadowing in the rescaling model

A simple model for the charm structure function of nuclei

Effects of heavy quarks in the dipole cross section with the KLN model

Perturbative Gluon Shadowing in Heavy Nuclei

Shadowing of Gluons in Perturbative Qcd: A Comparison of Different Models

Single Inclusive Hadron Production at the LHC with an Improved Impact Parameter Dependent Saturation Model

Color Transparency and Saturation in QCD

The behavior of the heavy-quarks structure functions at small-$x$

Small x, Saturation, and Diffraction in Collisions with electrons, protons, and nuclei

Improved Single Particle Potential for Transport Model Simulations of Nuclear Reactions Induced by Rare Isotope Beams

The Color Dipole Picture of low-x DIS: Model-Independent and Model-Dependent Results

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

Forward dijet production at the LHC within an impact parameter dependent TMD approach