PREX, CREX, and MREX: Probing Neutron Skins

• The experiments measure neutron skin thickness via helicity-dependent cross-section asymmetry, achieving sub-0.05 fm precision in heavy nuclei.
• They employ parity-violating electron scattering and Bayesian inference to extract the weak charge form factor and neutron density profile with minimal model dependence.
• Results from these initiatives provide crucial constraints on the nuclear symmetry energy, impacting neutron star modeling and energy density functional calibrations.

Parity-violating electron scattering (PVES) has become the preeminent experimental tool for model-independent measurements of the neutron density distribution in heavy nuclei. The PREX (Lead Radius Experiment), CREX (Calcium Radius Experiment), and the planned MREX (Mainz Radius EXperiment) form a sequence of high-precision initiatives that collectively address fundamental issues in nuclear structure, the density dependence of the symmetry energy, and the consistency between nuclear and astrophysical constraints.

1. Experimental Objectives, Technique, and Infrastructure

The principal aim of the PREX, CREX, and MREX programs is to extract the neutron skin thickness, defined as $$\Delta r_{np} = \langle r_n^2 \rangle^{1/2} - \langle r_p^2 \rangle^{1/2}$$, with sub-0.05 fm precision for heavy neutron-rich nuclei such as $^{208}$Pb and $^{48}$Ca. The significance of $\Delta r_{np}$ lies in its sensitivity to the nuclear symmetry energy and its slope parameter $L$ at nuclear saturation density $\rho_0$—a central quantity for both nuclear structure and neutron star modeling (Xu et al., 2020).

These experiments employ parity-violating elastic electron scattering at forward angles, using high-current, highly polarized, longitudinal electron beams incident on isotopically pure targets. The observable is the helicity-dependent cross-section asymmetry:

$$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$

which—dominated by weak neutral current (Z$^0$) exchange with the neutron distribution—directly probes the weak charge form factor $F_W(q)$. By comparing $F_W(q)$ to the precisely known charge form factor $^{208}$0 (from conventional electron scattering), the neutron distribution can be extracted with minimal model dependence (Zhao et al., 2024, Davoudiasl et al., 17 Dec 2025).

• PREX targets $^{208}$1Pb at JLab, measuring $^{208}$2 at $^{208}$3 fm$^{208}$4.
• CREX focuses on $^{208}$5Ca at JLab, at $^{208}$6 fm$^{208}$7. The lighter mass and amenability to ab initio calculations provide a critical cross-check for nuclear models.
• MREX, planned at Mainz/MESA, aims for sub-0.03 fm precision in both nuclei and possibly additional targets.

Current PVES facilities achieve statistical and systematic error control at the permille to few-percent level, with experimental precision continually improving due to advances in beam stability, detector technology, and background suppression (Reed et al., 2020).

2. Theoretical Framework, Model Independence, and Extraction Formalism

PVES observables are rigorously interpretable via the Standard Model. The leading-order Born formula for the parity-violating asymmetry relates $^{208}$8 to nuclear structure quantities as:

$^{208}$9

where $^{48}$0 is the Fermi constant, $^{48}$1 the fine-structure constant, and $^{48}$2 the weak charge form factor (Davoudiasl et al., 17 Dec 2025).

The extraction of $^{48}$3 uses distorted wave Born approximation (DWBA) calculations that incorporate strong coulomb distortions, but the method remains exceptionally model-independent compared to hadronic probes. Single-parameter (rms radius) or two-parameter Fermi (half-density radius $^{48}$4, surface diffuseness $^{48}$5) shapes for $^{48}$6 are fit to the data. Statistical consistency and systematic robustness are ensured by exploiting Markov-Chain Monte Carlo methods and Bayesian parameter estimation frameworks (Xu et al., 2020, Zhao et al., 2024, Zhang et al., 2022).

A recent refinement is the use of the difference $^{48}$7, which is measured directly and constrains both the neutron skin and the surface structure. In particular, with measurements at multiple momentum transfers (as foreseen for MREX), both $^{48}$8 and $^{48}$9 can be determined, fully specifying the neutron density profile (Reed et al., 2020).

3. Bayesian Inference, Symmetry Energy Constraints, and Posterior Structure

PREX, CREX, and MREX data are interpreted within Bayesian statistical frameworks using Skyrme-Hartree-Fock, relativistic mean field (RMF), and energy density functional (EDF) approaches. The focus is on the symmetry energy $\Delta r_{np}$0 and particularly its slope $\Delta r_{np}$1 at saturation density:

$\Delta r_{np}$2

Experimental likelihoods incorporate the measured $\Delta r_{np}$3 (or, equivalently, $\Delta r_{np}$4), with uncertainties dominated by beam statistics and detector systematics. Priors on nuclear matter parameters are typically chosen broad and uniform, reflecting the goal of unbiased inference. Posterior distributions for $\Delta r_{np}$5 and related parameters are obtained by marginalization over nuisance quantities (Xu et al., 2020, Zhao et al., 2024).

Numerically, the PREX and CREX data pull the inferred posteriors in different directions:

• PREX favors $\Delta r_{np}$6 fm, implying $\Delta r_{np}$7 MeV within most functionals.
• CREX yields $\Delta r_{np}$8 fm, implying $\Delta r_{np}$9 MeV (Reed et al., 2023, Kumar et al., 2023, Zhang et al., 2022).

The conflicting inferences highlight the need for model extensions, notably in handling the surface symmetry energy and higher-order empirical parameters that can decouple low- and high-density regimes (Mondal et al., 2022).

4. Model Tensions, Isovector Spin-Orbit Effects, and Resolution Strategies

A dominant theoretical challenge is the so-called CREX–PREX dilemma: canonical EDFs predict a strong correlation between $L$0 in $L$1Pb and $L$2Ca, yet experiment finds a large skin in Pb and a much smaller one in Ca. This suggests missing physics in the standard isovector channel.

Recent studies identify the isovector spin-orbit interaction, parameterized (in Skyrme or RMF models) by $L$3, as a powerful lever arm. Enhancement of the isovector spin-orbit term can selectively reduce $L$4 in $L$5Ca without strongly affecting $L$6Pb, matching the reported data. However, the magnitude of the required enhancement (e.g., $L$7 in RMF) is unphysically large and destroys agreement with empirical single-particle spectra and shell closures (Kunjipurayil et al., 10 Mar 2025, Zhao et al., 2024).

A more physically consistent approach is simultaneous refitting of central, spin-orbit, and tensor terms under combined experimental constraints. Linear "pocket formulas" express $L$8 and $L$9 as separate functions of $\rho_0$0 and $\rho_0$1, facilitating the efficient planning of future experiments such as MREX (Zhao et al., 2024).

5. MREX and the Next Generation of Parity-Violating Measurements

MREX at MESA will extend the capacity of PVES programs by enabling a second measurement at larger $\rho_0$2 for both $\rho_0$3Pb ($\rho_0$4 fm$\rho_0$5) and $\rho_0$6Ca ($\rho_0$7 fm$\rho_0$8). This provides heightened sensitivity to the surface diffuseness parameter $\rho_0$9 of the weak charge distribution. Simulation studies indicate that modest beam times ($$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$050–100 days) suffice to extract $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$1 to a few percent precision, provided detector and beam systematic controls (Reed et al., 2020).

By jointly fitting data at two $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$2-values, both the half-density radius $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$3 and diffuseness $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$4 are constrained, yielding an essentially model-independent neutron density profile. This fully determines $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$5 and $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$6 around and below $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$7, with minimal sensitivity to modeling assumptions. MREX's results will permit rigorous global Bayesian updating, integrating nuclear, astrophysical, and ab initio data to map the symmetry energy (and its density derivatives) (Xu et al., 2020, Davoudiasl et al., 17 Dec 2025).

6. Impact on Nuclear Matter Properties, Astrophysical EOS, and Future Outlook

The results from PREX, CREX, and anticipated MREX analyses feed directly into calibrations of nuclear density functionals, with broad consequences for neutron-rich nuclei and neutron star structure. Posterior distributions on $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$8, $$A_{PV} = \frac{\sigma_R - \sigma_L}{\sigma_R + \sigma_L}\ ,$$9, and other EOS parameters impact predictions for neutron star radii $^0$0, tidal deformabilities $^0$1, and core–crust transition densities (Kumar et al., 2023, Reed et al., 2023, Mondal et al., 2022).

A key finding is that models accommodating both PREX and CREX generally require either (a) extreme surface or higher-order couplings, or (b) a very stiff EOS at supra-nuclear densities—leading to neutron star radii and $^0$2 values exceeding current multimessenger bounds. This suggests that either the standard nuclear structure paradigm requires revision or the data reflect a confluence of "soft" and "stiff" symmetry energy behavior at different densities, or both (Reed et al., 2023, Zhao et al., 2024).

Planned and proposed future experiments, including multi-$^0$3 PVES (e.g., double-q CREX or extended MREX) and continuous-$^0$4 mapping at an Electron-Ion Collider, will further decouple and diagnose the various empirical parameters underlying the isovector sector and the symmetry energy (Davoudiasl et al., 17 Dec 2025, Reed et al., 2020). Systematic inclusion of QED radiative corrections, especially photon vacuum polarization, will be necessary as experimental precision approaches the sub-percent level (Reed et al., 4 Jan 2026).

---

Summary Table: Key Targets and Precision Goals for PREX, CREX, and MREX

Experiment	Target	$^0$5 (fm$^0$6)	Reported/Goal $^0$7 (fm)
PREX-II	$^0$8Pb	0.398	$^0$9 (goal), $F_W(q)$0 (reported)
CREX	$F_W(q)$1Ca	0.873	$F_W(q)$2 (goal), $F_W(q)$3 (reported)
MREX	$F_W(q)$4Pb	0.76	$F_W(q)$5 (goal, projected)
MREX	$F_W(q)$6Ca	1.28	$F_W(q)$7 (goal, projected)

---

The PREX, CREX, and MREX program, interpreted using Bayesian theory with advanced EDFs, is driving a paradigm shift in the precision mapping of the nuclear symmetry energy and confronting long-standing theoretical assumptions about the density dependence of the isovector channel in nuclear matter. The results inform both terrestrial and astrophysical observables and promise sharper convergence between nuclear structure and neutron star modeling as experimental, theoretical, and computational methodology advances further (Xu et al., 2020, Zhao et al., 2024, Kunjipurayil et al., 10 Mar 2025, Mondal et al., 2022, Reed et al., 2023).