⚛️ nuclear theory

A common finite-density charge-symmetry-breaking response in mirror displacement energies and charge radii

This study demonstrates that a unified charge-symmetry-breaking functional can simultaneously describe mirror displacement energies and charge-radius differences, revealing a common finite-density response that constrains effective coupling combinations and highlights the necessity of quantifying surface-gradient-sensitive corrections before using mirror charge radii as precise probes of neutron skins.

Original authors: Myeong-Hwan Mun, Kyoungsu Heo, Jubin Park, Myung-Ki Cheoun, H. Sagawa

Published 2026-07-22

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Myeong-Hwan Mun, Kyoungsu Heo, Jubin Park, Myung-Ki Cheoun, H. Sagawa

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Atomic Mirror and the Invisible Tug-of-War

Imagine the atomic nucleus as a bustling city built from two types of citizens: protons, who carry a positive electric charge, and neutrons, who are electrically neutral. In the world of nuclear physics, there's a fascinating concept called "mirror nuclei." These are pairs of atomic twins where the roles are swapped: if one twin has 10 protons and 14 neutrons, its mirror partner has 14 protons and 10 neutrons. In a perfect, charge-symmetric world, these twins would be identical in every way except for their names. However, reality is messier. The protons repel each other because they are all positively charged, like magnets with the same pole facing each other. This repulsion, known as the Coulomb force, stretches the proton-rich twin slightly differently than the neutron-rich one.

But there's a deeper mystery. Even after we account for that electric repulsion, the twins still don't behave exactly as physics textbooks predict. They have slightly different sizes and binding energies. This leftover difference is a clue that the strong nuclear force—the glue holding the city together—has a tiny, subtle bias against protons and neutrons swapping places. Physicists call this "charge-symmetry breaking." Understanding this bias is crucial because it acts like a ruler for measuring "neutron skins," the fuzzy layers of extra neutrons that coat heavy atoms. If we don't understand the bias, our ruler is crooked, and we can't accurately measure the size of the universe's most extreme objects, like neutron stars.

The Paper's Story: Tuning the Invisible Ruler

In this study, the authors act like master mechanics trying to tune a very sensitive instrument. They are using a mathematical framework called the "Skyrme energy-density functional," which is essentially a sophisticated recipe for calculating how protons and neutrons behave inside a nucleus. They started with a standard recipe that only included the electric repulsion (the Coulomb force) and found that it didn't quite match the real-world measurements of mirror nuclei. The "mirror displacement energy" (MDE)—the energy cost to swap protons and neutrons—was off, and the difference in their charge radii (sizes) was also slightly wrong.

To fix this, the researchers introduced a new ingredient to their recipe: a specific type of "charge-symmetry breaking" (CSB) force. Think of this force as a tiny, invisible hand that nudges protons and neutrons differently. They tested two different versions of this hand: one that acts uniformly throughout the entire nucleus (the "volume" term) and one that acts more strongly at the surface, like a crust (the "surface-gradient" term).

The team focused on four specific pairs of mirror nuclei that serve as their "anchors": Argon-34/Sulfur-34, Calcium-36/Sulfur-36, Calcium-38/Argon-38, and Nickel-54/Iron-54. By carefully adjusting the strength of their new CSB force to match the energy differences (MDEs) of these four pairs, they discovered something interesting. The data didn't tell them exactly how strong the "volume" force was versus the "surface" force. Instead, it told them that these two forces work together in a very specific, nearly locked ratio. It's as if they found a single "master dial" that controls the combination of both forces, rather than two independent knobs.

Once they calibrated this master dial using the energy data, they tested it on the size of the nuclei. They found that for one of their recipes (called SLy4), the calibration worked beautifully, reducing the errors in the predicted sizes to a tiny 0.0031 fm (femtometers). However, for the second recipe (SkM*), while the "type" of response was the same, the actual size corrections were larger and less perfect. This suggests that while the pattern of how the force works is robust, the exact amount of correction depends on the specific mathematical recipe used.

The paper then used this calibrated "master dial" to make predictions for four other proton-rich mirror pairs that haven't been measured yet: Titanium-40, Titanium-42, Chromium-46, and Iron-50. They predicted the charge radii for these atoms, finding that the differences between their two recipes were small but noticeable, ranging up to 0.021 fm.

Crucially, the authors point out that not all mirror nuclei follow this same rule. They checked a pair called Sulfur-30/Silicon-30 and found it behaved differently, with a much weaker surface response. This tells us that we can't just apply a single rule to every atom; the "response class" depends on the specific shell structure of the nucleus.

In the end, the paper concludes that before we can use mirror nuclei as clean tools to measure neutron skins or the properties of the nuclear equation of state, we must first quantify this specific, finite-density charge-symmetry breaking response. The authors have successfully mapped out a "response class" that works for a broad group of nuclei, but they warn that the exact size corrections remain dependent on the theoretical model chosen. It's a significant step forward in calibrating our nuclear rulers, ensuring that when we look at the smallest mirrors in the universe, we see a reflection that is as clear and true as possible.

Technical Summary: A Common Finite-Density Charge-Symmetry-Breaking Response in Mirror Displacement Energies and Charge Radii

Problem Statement
Mirror nuclei serve as differential probes of isospin symmetry, yet extracting clean information regarding neutron-skin thickness and the symmetry energy slope ( $L$ ) from mirror charge-radius differences ( $\Delta R^{\text{mirr}}_{\text{ch}}$ ) is complicated by nuclear isospin-symmetry breaking (ISB) beyond the Coulomb interaction. Traditionally, Mirror Displacement Energies (MDEs) and $\Delta R^{\text{mirr}}_{\text{ch}}$ are analyzed separately: MDEs constrain the isospin-breaking interaction, while radius differences are interpreted as proxies for neutron skins. This separation ignores a critical consistency requirement: the same Charge-Symmetry Breaking (CSB) response that accounts for the Coulomb-subtracted MDE residual must also generate a specific spatial response in charge radii and neutron skins. If MDE-calibrated CSB corrections significantly alter $\Delta R^{\text{mirr}}_{\text{ch}}$ at current experimental precision, a Coulomb-only analysis introduces a systematic, pair-dependent bias. Furthermore, the Coulomb-subtracted MDE residual is not a model-independent measure of bare CSB interaction but a differential quantity dependent on the chosen many-body framework.

Methodology
The authors formulate a direct mass–radius consistency test within a single self-consistent Energy-Density Functional (EDF) framework.

Baselines: Calculations utilize modified HFBTHO v4.0 solvers with two Skyrme EDFs: SLy4 and SkM*. A "Coulomb-only" baseline is established for each, subtracting Coulomb effects from experimental MDEs and $\Delta R^{\text{mirr}}_{\text{ch}}$ to define residuals.
CSB Functional: The study supplements the baseline with a class-III CSB functional comprising two terms:
1. A volume term: $H^{(0)}_{\text{CSB}} = \frac{1}{2}t^{\text{III}}_0 (\rho_n^2 - \rho_p^2)$ .
2. A surface-gradient term: $H^{(\Delta)}_{\text{CSB}} = C^{\text{III}}_\Delta (\rho_0 \Delta \rho_1 + \rho_1 \Delta \rho_0)$ .
Linear Response Analysis: Near the Coulomb-only solution, the MDE response is treated as locally linear. The authors define a response coefficient $\lambda_B$ for each mirror pair, relating the volume and surface-gradient couplings.
Fitting Strategy:
- MDE-only fit: Determines the effective coupling combination $t^{\text{III}}_{\text{eff}} = t^{\text{III}}_0 - \lambda_{\text{cl}} C^{\text{III}}_\Delta$ that minimizes MDE residuals.
- Joint fit: Uses measured anchor charge radii to select a specific point along the MDE-compatible "valley" of couplings, without introducing radius-specific parameters.
Anchors and Targets: Four measured pairs ( $^{34}\text{Ar}$ – $^{34}\text{S}$ , $^{36}\text{Ca}$ – $^{36}\text{S}$ , $^{38}\text{Ca}$ – $^{38}\text{Ar}$ , $^{54}\text{Ni}$ – $^{54}\text{Fe}$ ) serve as anchors. Predictions are made for unmeasured proton-rich partners ( $^{40}\text{Ti}$ , $^{42}\text{Ti}$ , $^{46}\text{Cr}$ , $^{50}\text{Fe}$ ) based on the jointly calibrated response.

Key Contributions and Results

Common Response Class: The four anchor pairs form a compact "surface-gradient-like" response class in both EDFs. The effective coupling ratio is determined to be $\lambda^{\text{SLy4}}_{\text{cl}} = 1.32 \pm 0.04 \, \text{fm}^{-2}$ and $\lambda^{\text{SkM*}}_{\text{cl}} = 1.20 \pm 0.05 \, \text{fm}^{-2}$ . This indicates that MDEs primarily constrain the effective combination $t^{\text{III}}_0 - \lambda_{\text{cl}} C^{\text{III}}_\Delta$ rather than the individual couplings.
Mass–Radius Consistency:
- SLy4: The joint volume-plus-surface-gradient fit successfully reduces residuals for both MDEs (RMS $\approx$ 0.0529 MeV) and $\Delta R^{\text{mirr}}_{\text{ch}}$ (RMS $\approx$ 0.0031 fm).
- SkM:* While SkM* recovers the same response class, it exhibits a larger tension between mass and radius, resulting in larger $\Delta R^{\text{mirr}}_{\text{ch}}$ residuals (RMS $\approx$ 0.0136 fm) compared to SLy4.
EDF Dependence: The spatial propagation of the MDE-selected CSB direction is EDF-dependent. While the response class (the ratio of volume to surface sensitivity) is robust across baselines, the absolute radius corrections vary. The calibrated response shifts proton-rich mirror skins by $\sim 10^{-2}$ fm, a magnitude comparable to current experimental precision.
Predictions: The study provides out-of-sample predictions for the charge radii of $^{40}\text{Ti}$ , $^{42}\text{Ti}$ , $^{46}\text{Cr}$ , and $^{50}\text{Fe}$ . The two-EDF spread in these predictions reaches up to 0.021 fm (notably for $^{46}\text{Cr}$ ), highlighting model dependence. The $^{42}\text{Ti}$ prediction is the most stable, with a spread of only 0.0020 fm.
Regime Validity: The analysis confirms that the response ratio is not universal across all mirror nuclei. Systems like $^{30}\text{S}$ – $^{30}\text{Si}$ and $^{32}\text{Ar}$ – $^{32}\text{Si}$ fall outside the compact response regime defined by the anchors, suggesting that extrapolation must be restricted to nuclei with similar radial and shell structures.

Significance and Claims
The paper asserts that the MDE-selected and radius-constrained finite-density isospin-breaking response must be quantified within the same framework before $\Delta R^{\text{mirr}}_{\text{ch}}$ can be reliably used as a clean probe for neutron skins or the symmetry energy.

Robustness vs. Precision: The agreement between SLy4 and SkM* regarding the response-class assignment is more robust and meaningful than their agreement on absolute radius corrections, which remain EDF-dependent.
Systematic Bias: A Coulomb-only analysis carries a systematic bias because the CSB corrections required to fit MDEs inevitably alter charge radii.
Limitations: The authors explicitly state that the extracted couplings are effective, baseline-dependent parameters, not model-independent measurements of a bare CSB interaction. The difference between SLy4 and SkM* is presented as a limited diagnostic of model dependence, not a statistically complete uncertainty budget. Future precision constraints require a broader ensemble of EDFs and systematic studies of deformation, continuum, and shell effects.

The Atomic Mirror and the Invisible Tug-of-War

The Paper's Story: Tuning the Invisible Ruler

Technical Summary: A Common Finite-Density Charge-Symmetry-Breaking Response in Mirror Displacement Energies and Charge Radii

More like this