Deviations from the Isobaric Multiplet Mass Equation due to threshold states

Recent studies have completed the A=16 isospin quintets for states with spin/parity Jπ =0+ and 2+. The dependence of their masses as a function of isospin projection shows evidence for deviations from quadratic behavior indicating isospin violation beyond the expectation from two- body forces. The deviation is most pronounced for the 2+ states. Predictions from the Shell Model Embedded in the Continuum (SMEC) allow us to explain that this isospin violation is associated with a modification of the nuclear structure due to the open-quantum-system nature of the proton- rich members of the quintet. In particular, the 0+ and 2+ states in 16Ne and the 2+ state in 16F are threshold resonances located just above a proton-decay threshold where s-wave coupling to the continuum is expected. The measured deviations of these threshold states from the quadratic behavior of the remaining members of the multiplets makes it possible to obtain information on the magnitude and the energy dependence of the continuum-coupling energy correction. Continuum coupling is also indicated for the ground state of 8C, but this time through p-wave coupling.

💡 Research Summary

The paper investigates the origin of deviations from the Isobaric Multiplet Mass Equation (IMME) in the A = 16 isospin quintets for Jπ = 0⁺ and 2⁺ states. While the IMME predicts a quadratic dependence of nuclear masses on the isospin projection Tz when only two‑body forces are at work, recent high‑precision mass measurements reveal significant departures from this quadratic behavior, especially for the 2⁺ quintet. By measuring the masses of the proton‑rich members ¹⁶Ne and ¹⁶F using a ¹⁷Ne secondary beam and invariant‑mass spectroscopy, the authors obtain mass excesses that, when fitted with a pure quadratic IMME, yield poor χ² values (χ²/ν ≈ 3–4). Adding a cubic term dramatically improves the fits, producing positive d coefficients of +4.0(22) keV for the 0⁺ quintet and +8.9(31) keV for the 2⁺ quintet; the latter represents a 2.9σ deviation, indicating a genuine breakdown of the quadratic law.

The authors propose that these anomalies arise not from conventional isospin mixing with nearby T ≠ 2 states, but from coupling of the resonant states to the particle‑emission continuum. The 0⁺ and 2⁺ states in ¹⁶Ne and the 2⁺ state in ¹⁶F lie just above proton‑decay thresholds and couple predominantly to s‑wave proton channels. Such “threshold resonances” experience a strong modification of their wave functions because the open quantum‑system nature of the nucleus forces the bound‑state component to mix with the scattering continuum. This mixing lowers the energy of the resonant state by an amount termed the continuum‑coupling correlation energy, E_corr, which is negative and most pronounced when the state sits a few hundred keV above the threshold.

To quantify this effect, the authors employ the Shell Model Embedded in the Continuum (SMEC), an extension of the traditional shell model that partitions the Hilbert space into a bound‑state subspace Q₀ and a one‑particle scattering subspace Q₁. The effective Hamiltonian H(E) = H_Q₀Q₀ + W_Q₀Q₀(E) contains an energy‑dependent term W that couples Q₀ to Q₁ via the one‑nucleon Green’s function. Using Zuker‑Buck‑McGrory and Cohen‑Kurath effective interactions for A = 16 and A = 8, respectively, and a Wigner contact force for the continuum coupling, the SMEC calculations produce E_corr values that depend on the distance E_p above the relevant threshold. For the three threshold states (0⁺ in ¹⁶Ne, 2⁺ in ¹⁶Ne, and 2⁺ in ¹⁶F) the calculated E_corr lies between –0.35 and –0.40 MeV, in excellent agreement with the experimental “missing” mass inferred from the quadratic extrapolation of the neutron‑rich members. The ground state of ⁸C, which couples via a p‑wave channel, shows a much smaller correction (≈ –0.15 MeV), consistent with its weaker continuum effect.

The authors illustrate the analysis by plotting residuals of quadratic IMME fits (Fig. 2) and by comparing the experimental residuals with the SMEC‑predicted E_corr (Fig. 3). The discrepancy between the quadratic extrapolation (based on non‑threshold members) and the measured masses of the threshold states directly yields the magnitude of the continuum correction. This approach demonstrates that the positive cubic coefficients observed in the IMME fits are a natural consequence of the extra binding provided by continuum coupling for states located just above decay thresholds.

In conclusion, the work provides the first quantitative link between threshold‑state continuum coupling and IMME violations. It shows that open‑quantum‑system effects can produce sizable, localized deviations from the quadratic mass law, especially for s‑wave resonances near proton thresholds. The findings suggest that future high‑precision mass measurements of nuclei near particle emission thresholds, combined with SMEC or similar continuum‑shell‑model frameworks, will be essential for a comprehensive understanding of isospin symmetry breaking and for refining the IMME to include continuum‑induced corrections. This represents a significant step toward integrating bound‑state nuclear structure with reaction dynamics in a unified theoretical description.