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

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 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)$ from Skyrme and Relativistic Mean Field models together, a discrepancy between these 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 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 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 symmetry energy and its slope parameter 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.

close to the published version; extended analysis and minor changes; 31 pages, 14 figures, 5 tables

References in corpus (90)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Equation of state and neutron star properties constrained by nuclear physics and observation
• Skyrme Interaction and Nuclear Matter Constraints
• Accurate Determination of the Neutron Skin Thickness of ${}^{208}$Pb through Parity-Violation in Electron Scattering
• Improved nuclear matter calculations from chiral low-momentum interactions
• Implications of PREX-II on the equation of state of neutron-rich matter
• Relativistic Mean-Field Hadronic Models under Nuclear Matter Constraints
• Accurate nuclear radii and binding energies from a chiral interaction
• Nuclear effective field theory: status and perspectives
• Nuclear energy density optimization: Large deformations
• The maximum mass and radius of neutron stars and the nuclear symmetry energy
• Neutron matter at next-to-next-to-next-to-leading order in chiral effective field theory
• Constraints on neutron star radii based on chiral effective field theory interactions
• Chiral interactions up to next-to-next-to-next-to-leading order and nuclear saturation
• Constraining Neutron-Star Matter with Microscopic and Macroscopic Collisions
• How well do we know the neutron-matter equation of state at the densities inside neutron stars? A Bayesian approach with correlated uncertainties
• Precision Determination of the Neutral Weak Form Factor of ${}^{48}$Ca
• Building relativistic mean field models for finite nuclei and neutron stars
• A Guided Tour of Ab Initio Nuclear Many-Body Theory
• The equation of state for dense nucleonic matter from a metamodeling. I. Foundational aspects
• Ab initio predictions link the neutron skin of ${}^{208}$Pb to nuclear forces
• Nuclear Equation of State from ground and collective excited state properties of nuclei
• Ab initio limits of atomic nuclei
• Nuclear energy density optimization: Shell structure
• Saturation with chiral interactions and consequences for finite nuclei
• Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ${}^{208}$Pb with Minimal Modeling Assumptions
• Neutron star tidal deformabilities constrained by nuclear theory and experiment
• Dense Nuclear Matter Equation of State from Heavy-Ion Collisions
• Searching for isovector signatures in the neutron-rich oxygen and calcium isotopes
• Asymmetric nuclear matter based on chiral two- and three-nucleon interactions
• Chiral Effective Field Theory and the High-Density Nuclear Equation of State
• Incompressibility of neutron-rich matter
• Quantifying Correlated Truncation Errors in Effective Field Theory
• Accurate bulk properties of nuclei from $A=2$ to $\infty$ from potentials with $\Delta$ isobars
• Uncertainty analysis and order-by-order optimization of chiral nuclear interactions
• Quantifying uncertainties and correlations in the nuclear-matter equation of state
• Three-Nucleon Forces: Implementation and Applications to Atomic Nuclei and Dense Matter
• Information content of the parity-violating asymmetry in ${}^{208}$Pb
• Uncertainty Quantification for Nuclear Density Functional Theory and Information Content of New Measurements
• Equation of state of nuclear and neutron matter at third-order in perturbation theory from chiral EFT
• Constraints on Nuclear Symmetry Energy Parameters
• From finite nuclei to the nuclear liquid drop: leptodermous expansion based on the self-consistent mean-field theory
• Exploring sd-shell nuclei from two- and three-nucleon interactions with realistic saturation properties
• Towards order-by-order calculations of the nuclear and neutron matter equations of state in chiral effective field theory
• Family of Chiral Two- plus Three-Nucleon Interactions for Accurate Nuclear Structure Studies
• Accurate nuclear masses from a three parameter Kohn-Sham DFT approach (BCPM)
• Delta isobars and nuclear saturation
• Energy density functional for nuclei and neutron stars
• Quantified limits of the nuclear landscape
• Combined theoretical analysis of the parity-violating asymmetry for ${}^{48}$Ca and ${}^{208}Pb$
• Get on the BAND Wagon: A Bayesian Framework for Quantifying Model Uncertainties in Nuclear Dynamics
• Nuclear and neutron-star matter from local chiral interactions
• New equations of state constrained by nuclear physics, observations, and QCD calculations of high-density nuclear matter
• A Detailed Examination of Astrophysical Constraints on the Symmetry Energy and the Neutron Skin of ${}^{208}$Pb with Minimal Modeling Assumptions
• Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter
• Relativistic mean-field study of neutron-rich nuclei
• Isovector properties of the Gogny interaction
• A Minimal Nuclear Energy Density Functional
• What is ab initio in nuclear theory?
• Beyond the proton drip line: Bayesian analysis of proton-emitting nuclei
• Probing chiral interactions up to next-to-next-to-next-to-leading order in medium-mass nuclei
• Implications of PREX-2 on the electric dipole polarizability of neutron rich nuclei
• Constraints on Skyrme Equations of State from Properties of Doubly Magic Nuclei and Ab-Initio Calculations of Low-Density Neutron Matter
• New Ideas in Constraining Nuclear Forces
• Bayesian inference of nuclear symmetry energy from measured and imagined neutron skin thickness in ${}^{116,118,120,122,124,130,132}$Sn, ${}^{208}$Pb, and ${}^{48}$Ca
• Model reduction methods for nuclear emulators
• Nuclear charge and neutron radii and nuclear matter: trend analysis
• Eigenvector Continuation and Projection-Based Emulators
• Insights into nuclear saturation density from parity violating electron scattering
• Large and massive neutron stars: Implications for the sound speed in dense QCD
• Power counting in chiral effective field theory and nuclear binding
• Bayesian refinement of covariant energy density functionals
• Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter
• The first five years of gravitational-wave astrophysics
• What is ab initio?
• Symmetry energy and neutron star properties constrained by chiral effective field theory calculations
• Emulating \emph{ab initio} computations of infinite nucleonic matter
• Nuclear Structure at the Crossroads
• Applications of reduced basis methods to the nuclear single particle spectrum
• Applying the Density Matrix Expansion with Coordinate-Space Chiral Interactions
• Nuclear-matter saturation and symmetry energy within $\Delta$--full chiral effective field theory
• Constraining the nonanalytic terms in the isospin-asymmetry expansion of the nuclear equation of state
• Optimization and Supervised Machine Learning Methods for Fitting Numerical Physics Models without Derivatives
• Bayesian Analysis of $\chi$EFT at Leading Order in a Modified Weinberg Power Counting Approach
• Energy systematics of heavy nuclei -- mean field models in comparison
• Tabulated Equations of State From Models Informed by Chiral Effective Field Theory
• Calibration of Energy Density Functionals with Deformed Nuclei
• Neutron Skins: Weak Elastic Scattering and Neutron Stars
• Optimized nuclear energy density functionals including long-range pion contributions
• Comparing different density-matrix expansions for long-range pion exchange

Cited by in corpus (9)

• Diagrammatic ab initio methods for infinite nuclear matter with modern chiral interactions
• Emulators for Scarce and Noisy Data: Application to Auxiliary-Field Diffusion Monte Carlo for Neutron Matter
• Effect of Magnetic Fields on Urca Rates in Neutron Star Mergers
• The nucleardatapy toolkit for simple access to experimental nuclear data, astrophysical observations, and theoretical predictions
• Microscopic constraints for the equation of state and structure of neutron stars: a Bayesian model mixing framework
• The maximum mass and rotational kinetic energy of rapidly rotating neutron stars
• Structure of Anisotropic Magnetized Neutron Stars in f(R,T) Gravity with Realistic Equation of State
• Toward a Unified Understanding of the Dense Matter Equation of State
• Renormalization-Group Invariant Parity-Doublet Model for Nuclear and Neutron-Star Matter