Characterizing the Nuclear Models Informed by PREX and CREX: a View from Bayesian Inference

New measurements of the weak charge density distributions of $^{48}$Ca and $^{208}$Pb challenge existing nuclear models. In the post-PREX-CREX era, it is unclear if current models can simultaneously describe weak charge distributions along with accurate measurements of binding energy and charge radii. In this letter, we explore the parameter space of relativistic and non-relativistic models to study the differences between the form factors of the electric and weak charge distributions, $\Delta F=F_{ch}-F_{w}$, in $^{48}$Ca and $^{208}$Pb. We show, for the first time, the parts of the mean-field models that are the most important in determining the relative magnitude of the neutron skin in lead and calcium nuclei. We carefully disentangle the tension between the PREX/CREX constraints and the ability of the RMF and Skyrme models to accurately describe binding energies and charge radii. We find that the nuclear symmetry energy coefficient $S_\mathrm{V}$ and the isovector spin-orbit coefficient $b'_4$ play different roles in determining $\Delta F$ of $^{48}$Ca and $^{208}$Pb. Consequently, adjusting $S_\mathrm{V}$ or $b'_4$ shifts predicted $\Delta F$ values toward or away from PREX/CREX measurements. Additionally, $S_\mathrm{V}$ and the slope L are marginally correlated given the constraints of our Bayesian inference, allowing us to infer them separately from PREX and CREX data.