CREX and PREX-II reconciled within energy-density functional theory

By relaxing the implicit surface-bulk coupling constraint in energy-density functional theory, this study demonstrates that new functionals can simultaneously reconcile the conflicting CREX and PREX-II neutron-skin measurements with other nuclear and astrophysical observables without requiring extreme values of the symmetry-energy slope parameter.

The Gist

The Great Neutron Skin Puzzle

Imagine the atomic nucleus as a crowded dance floor. In some heavy atoms, like Lead-208, there are more "dancers" called neutrons than "dancers" called protons. Because they don't like to be squeezed too tightly by the protons, the extra neutrons tend to hang out on the very edge of the dance floor, forming a fuzzy outer layer. Physicists call this the "neutron skin."

For a long time, scientists had a rule of thumb: The stiffer the dance floor (the nuclear force), the thicker the fuzzy skin.

• Stiff floor (Hard EoS): Neutrons push out hard, creating a thick skin.
• Soft floor (Soft EoS): Neutrons stay closer to the center, creating a thin skin.

The Conflict: Two Experiments, Two Answers

Recently, two major experiments tried to measure the thickness of this skin, but they got results that seemed to break the rules:

1. The CREX Experiment (measuring Calcium-48): Found a very thin skin. This suggested the nuclear floor is soft.
2. The PREX-II Experiment (measuring Lead-208): Found a very thick skin. This suggested the nuclear floor is stiff.

It was like measuring the height of two different people and concluding that one is a dwarf while the other is a giant, even though standard physics said they should be similar sizes based on their weight. This created a huge headache for nuclear physicists: How can the same laws of physics produce a thin skin for one atom and a thick skin for another?

The Old Way vs. The New Idea

The Old Way (The "One-Size-Fits-All" Recipe):
Standard theories assumed that the rules governing the center of the atom (where it's dense and crowded) were exactly the same as the rules governing the edge (where it's sparse and dilute). It was like assuming a recipe for a dense chocolate cake works perfectly if you just spread it out thin on a plate. The theory said, "If we know how the center behaves, we can mathematically guess how the edge behaves."

This assumption created a "tension." If you tweaked the math to make Lead's skin thick, Calcium's skin became too thick. If you fixed Calcium, Lead's skin became too thin. You couldn't satisfy both.

The New Idea (The "Special Edge" Recipe):
The author of this paper, Panagiota Papakonstantinou, suggests that the old recipe was missing a key ingredient.

Imagine the edge of the atom not as a thin version of the center, but as a completely different environment.

• The Center: A dense, uniform liquid (like a packed subway car).
• The Edge: A sparse, gas-like cloud (like people standing far apart in a park).

The paper argues that the physics of the "park" (the dilute edge) doesn't have to follow the exact same rules as the "subway car" (the dense center). In the real world, matter behaves differently when it's spread out; particles might clump together or act like a gas rather than a liquid.

The Solution: Decoupling the Two Worlds

The author introduced a new mathematical "switch" (a special term in the equations) that only turns on at the very edge of the atom.

• Before: The math forced the edge to be a smooth, predictable extension of the center.
• After: The math allows the edge to have its own unique personality.

By "decoupling" the edge from the center, the author was able to tweak the physics of the surface without breaking the physics of the core.

The Result: Everyone Wins

With this new approach, the author found a set of rules that works for everything:

1. Calcium-48: Gets its thin skin (satisfying CREX).
2. Lead-208: Gets its thick skin (satisfying PREX-II).
3. Neutron Stars: The math still correctly predicts how massive neutron stars (which are essentially giant atomic nuclei) should behave.
4. Electric Polarizability: The atoms still wiggle and react to electric fields exactly as observed in experiments.

The Big Picture Analogy

Think of the nucleus like a city.

• The Old Theory assumed that the rules for building a skyscraper in the city center were exactly the same as the rules for building a tent in the suburbs. If you changed the zoning laws to make skyscrapers taller, the tents would automatically get taller too. This didn't match reality.
• The New Theory realizes that the suburbs (the nuclear surface) have their own unique zoning laws. You can build a massive skyscraper in the center (Lead-208) and a small tent in the suburbs (Calcium-48) without them fighting each other.

Why This Matters

This paper solves a major mystery in nuclear physics. It shows that we don't need to invent "extreme" or "weird" physics to explain the data. Instead, we just needed to stop assuming that the edge of an atom behaves exactly like its center.

It also gives us a better understanding of neutron stars. Since these stars are made of the same stuff as atomic nuclei, understanding how the "skin" of an atom works helps us predict how big and heavy these cosmic giants can be before they collapse.

In short: The author fixed the math by realizing that the "edge" of an atom is a special place with its own rules, allowing us to finally agree on the size of the neutron skin for both Calcium and Lead.

Technical Summary

1. The Problem: The CREX-PREX-II Tension

The paper addresses a significant contradiction in nuclear physics arising from two recent high-precision experiments:

• CREX (Calcium-48): Measured a thin neutron skin ($R_{np} \approx 0.121$ fm), favoring "soft" equations of state (EoS) with a low symmetry-energy slope parameter ($L$).
• PREX-II (Lead-208): Measured a thick neutron skin ($R_{np} \approx 0.283$ fm), favoring "stiff" EoS with a high $L$.

Standard Energy-Density Functional (EDF) theories predict a strong, linear correlation between neutron-skin thickness and $L$. Consequently, standard models cannot simultaneously reproduce both measurements without violating other constraints, such as the electric dipole polarizability ($a_D$) of $^{208}$Pb or the mass-radius relations of neutron stars. Previous attempts to resolve this via Bayesian analysis, isovector spin-orbit forces, or alpha clustering have either failed or required unphysical parameter values.

2. Methodology: Decoupling Surface and Bulk Physics

The author proposes that the tension arises from an implicit constraint in standard EDFs: the assumption that the density dependence of the functional at the dilute nuclear surface is a smooth extrapolation of the functional near saturation density (uniform matter).

Key Methodological Innovations:

• The KIDS Framework: The study utilizes the KIDS (Korea-IBS-Daegu-SKKU) framework, which translates EoS parameters into a non-relativistic EDF resembling a Skyrme potential but with extended density dependence.
• Introduction of a Dilute-Matter Term ($V_d$): To break the surface-bulk coupling, a new term is added to the pseudo-potential, active only at low densities (surface region):
  $$V_d = \frac{1}{6}(t_d + y_d P_\sigma)\rho_a d e^{-b_d \rho^2} \delta(\vec{r}_1 - \vec{r}_2)$$
  This term includes an exponential factor ensuring it vanishes at densities around $\rho_0/10$, leaving the high-density behavior (neutron star interiors) and saturation properties ($\rho_0, E_0, K_0, J, L$) untouched.
• Monte Carlo Sampling: The author employs the Metropolis-Hastings Monte Carlo (MHMC) method to sample the parameter space of the new term ($t_d, y_d, a_d, b_d$) alongside standard gradient and spin-orbit coefficients ($C_{12}, W_{02}$).
• Constraints: The sampling is driven by a log-posterior distribution maximizing agreement with:
  1. PREX-II ($^{208}$Pb skin) and CREX ($^{48}$Ca skin) data.
  2. Global nuclear properties (binding energies and charge radii of 19 nuclei, quantified by Average Deviation Per Datum, ADPD $\leq 1.5\%$).
  3. Symmetry energy at low density ($S(\rho_d) \approx 5 \pm 5$ MeV).
  4. Post-hoc validation against electric dipole polarizabilities ($a_D$) of $^{208}$Pb and $^{48}$Ca using Random Phase Approximation (RPA).

3. Key Contributions

• Relaxation of Surface-Bulk Coupling: The paper demonstrates that the neutron skin thickness is not solely determined by the symmetry energy slope $L$ at saturation. By allowing independent control over the dilute-density sector (nuclear surface), the strict correlation between $L$ and $R_{np}$ is broadened.
• Simultaneous Reconciliation: The study successfully identifies EDFs that reproduce the neutron skins of both $^{48}$Ca and $^{208}$Pb, their dipole polarizabilities, and realistic neutron-star mass-radius relations, all while maintaining standard saturation properties.
• Identification of a New Degree of Freedom: The work highlights that the low-density behavior of nuclear matter is currently underconstrained by existing data and represents a critical degree of freedom in EDF theory.

4. Results

• Success of the Dilute Term: When $V_d$ is set to zero (standard EDF behavior), the distributions for neutron skins are too narrow to accommodate both CREX and PREX-II simultaneously. Introducing $V_d$ broadens the distributions significantly, allowing thousands of parameter sets to satisfy both skin constraints.
• Validation with Dipole Polarizability:
  • EDFs based on a stiff EoS (EoS(65)) were rejected because they predicted dipole polarizabilities ($a_D$) far exceeding experimental values.
  • EDFs based on APR (Akmal-Pandharipande-Ravenhall) and QMC (Quantum Monte Carlo) EoSs, modified with $V_d$, successfully reproduced all four observables ($R_{np}$ for both nuclei and $a_D$ for both nuclei).
  • Specific successful models (e.g., APR(B), APR(C), QMC(A)) were identified.
• Physical Consistency:
  • The successful models preserve realistic single-particle spectra (no unphysical orbital ordering) and reproduce the neutron skin of $^{90}$Zr.
  • The neutron-star mass-radius relations remain consistent with dense matter constraints because the modification vanishes at high densities.
  • The symmetry energy $S(\rho)$ shows a strong suppression at low densities ($\rho < \rho_0/10$), effectively absorbing surface effects that standard gradient terms cannot capture.

5. Significance

This paper resolves the apparent contradiction between CREX and PREX-II without requiring extreme or unphysical values for the symmetry-energy slope parameter $L$.

• Theoretical Implication: It challenges the standard assumption that the nuclear surface is merely a smooth extrapolation of uniform matter. Instead, it suggests that the "dilute" regime on the nuclear surface involves complex physics (potentially related to clustering or gas-like states) that must be treated independently in EDFs.
• Astrophysical Impact: The reconciliation is achieved without compromising the description of neutron stars, maintaining consistency between terrestrial nuclear experiments and multimessenger astrophysics.
• Future Directions: The results suggest that future EDF optimizations must explicitly address the low-density sector, as it is a major source of uncertainty in predicting neutron-skin thicknesses and symmetry energy properties.