Nuclide Shift in Nuclear & Atomic Observables

• Nuclide shift is a collective term for displacement phenomena in atomic energies and nuclear decay driven by changes in nuclide mass and charge distribution.
• Theoretical frameworks such as relativistic Hartree-Fock and CI+all-order methods quantify shifts via mass and field contributions, revealing measurable effects in heavy elements.
• Observable nuclide shifts impact applications ranging from astrophysical isotope identification to detector calibration, reflecting underlying nuclear structure and reaction dynamics.

In the literature considered here, “nuclide shift” functions as an umbrella designation for several distinct but related displacement phenomena produced by changing nuclide identity or by revising a nuclide’s nuclear parameters. Its most established meaning is the isotope shift of atomic energies, transition frequencies, wavelengths, or total electron binding energies, where the underlying drivers are changes in nuclear mass and charge distribution. Closely related usages occur in nuclear decay, where a revised $Q$-value shifts phase space and predicted half-lives, and in reaction or measurement contexts, where nuclide production distributions or detector peak positions are displaced under specific dynamical conditions (Dzuba et al., 29 Jul 2025, Fink et al., 2011, Feng, 2023, Bazlov et al., 2023).

1. Formal definitions and principal observables

For atomic spectroscopy, the isotope-shift problem is conventionally decomposed into a mass shift and a field shift. In the broad formulation used for the total electron binding energy,

$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$

where $F$ is the field-shift coefficient, $G$ is a small quadratic correction, and $K_{\rm NMS}$ and $K_{\rm SMS}$ are the normal and specific mass-shift coefficients. For spectral lines, the same physics is usually written as a transition shift,

$$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$

These formulas encode two complementary sensitivities: nuclear recoil and finite nuclear size (Dzuba et al., 29 Jul 2025, Nazé et al., 2014).

In heavy systems, the field shift typically dominates. For the total electron binding energy with $\Delta A=2$, the cited calculations show that the normal mass shift dominates strongly at low $Z$, that mass shift and field shift are comparable around $Z\approx 38$, and that the field shift dominates above that point; representative ratios given are $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$0 for Ne, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$1 for Sr, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$2 for Xe, and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$3 for Og (Dzuba et al., 29 Jul 2025). In superheavy-atom work, the same dependence is often compressed into the interpolation law

$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$4

which is motivated by the scaling $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$5 and is used when explicit radius differences are unavailable (Dzuba et al., 2017).

A central misconception in older practice is that the field shift is exhausted by $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$6. The cited work consistently shows that this is only a leading approximation. For heavy and deformed nuclei, higher radial moments, deformation, surface diffuseness, and even central depression can contribute at measurable levels (Papoulia et al., 2016, Flambaum et al., 2019, Allehabi et al., 2020).

2. Computational frameworks for nuclide-shift theory

The modern theory of nuclide shifts is methodologically heterogeneous because the relevant observable may be a total binding energy, a line shift, or a highly charged-ion transition. For total electron binding energies up to $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$7, relativistic Hartree-Fock calculations including the Breit interaction are used; field-shift coefficients are extracted by varying the nuclear charge radius and fitting the energy variation with a parabola in $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$8 (Dzuba et al., 29 Jul 2025). For heavy multivalent systems such as $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$9, a finite-field method is combined with CI+all-order many-body theory: the Hamiltonian is modified as $F$0, and the field-shift constant is obtained from

$F$1

This framework was introduced specifically to handle heavy atoms and ions with several valence electrons and strong configuration mixing (Okhapkin et al., 2015).

For isoelectronic sequences, the relevant electronic factors are computed with multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction. In the Be-, B-, C-, and N-like sequences, these calculations supply $F$2, $F$3, and the electronic density at the nucleus, thereby allowing direct construction of level and transition isotope shifts from nuclear inputs (Nazé et al., 2014). For highly charged Li-like ions, the treatment is more specialized: the nuclear recoil contribution is split into $F$4, $F$5, $F$6, and $F$7, and evaluated with a hybrid perturbative plus CI-DFS strategy, while the field shift is computed with CI-DFS including electron correlation, Breit, and QED corrections (Zubova et al., 2014).

A further methodological development is the explicit coupling of nuclear and atomic structure theory. Realistic nuclear charge distributions from Skyrme-Hartree-Fock-Bogoliubov calculations are combined with multiconfiguration Dirac-Hartree-Fock atomic calculations to produce field shifts that depend on deformation and diffuseness, rather than on a schematic Fermi model alone (Papoulia et al., 2016). In the nobelium case, covariant density functional theory charge densities were interfaced with CI+MBPT atomic calculations, and the resulting isotope shifts were used as a benchmark for competing nuclear models (Allehabi et al., 2020).

3. Sensitivity to nuclear size, deformation, and shell structure

The field shift is a probe of the nuclear charge distribution inside the volume sampled by the electronic density. In reformulated form, the first-order field shift can be expanded as

$F$8

so that $F$9, $G$0, and $G$1 enter beyond the usual rms-radius term. The cited analysis shows that nuclear deformation and variations in diffuseness give measurable contributions, and that with a suitably chosen orthogonal basis $G$2, two sufficiently independent transitions can in principle be used to extract $G$3 and $G$4 (Papoulia et al., 2016).

The same point appears in the nobelium literature in a different parameterization. There, the isotope shift is fitted as

$G$5

and for neighboring isotopes reduced to

$G$6

This demonstrates that quadrupole deformation $G$7 can enter explicitly. The work further argues that earlier interpretation of the measured $G$8No shift purely in terms of $G$9 should be amended once deformation is taken into account, and that at least two transitions are required to disentangle $K_{\rm NMS}$0 from $K_{\rm NMS}$1 (Allehabi et al., 2020).

Superheavy systems amplify these effects. In calculations for E120$K_{\rm NMS}$2, a central depletion of charge density changes the isotope shift by about $K_{\rm NMS}$3 of the reference shift, while quadrupole deformation with $K_{\rm NMS}$4 produces an effect of order $K_{\rm NMS}$5 for $K_{\rm NMS}$6 states, roughly $K_{\rm NMS}$7 of the reference isotope shift. For the $K_{\rm NMS}$8 transition, fitting a single size parameter leaves a residual mismatch of about $K_{\rm NMS}$9–$K_{\rm SMS}$0, which the authors regard as potentially detectable. The same study finds that the relativistic $K_{\rm SMS}$1 law is extremely stable for spherical nuclei, with only about $K_{\rm SMS}$2 variation in $K_{\rm SMS}$3, but fluctuates by several percent for deformed nuclei (Flambaum et al., 2019).

Nuclide shifts along isotope chains can also reflect shell-structure physics at the level of nuclear wave functions. In Pb nuclei, the isotope shift

$K_{\rm SMS}$4

has a pronounced kink at $K_{\rm SMS}$5. The cited Hartree-Fock-Bogolyubov calculations show that a density-dependent three-nucleon spin-orbit term can reproduce this kink without requiring near-degeneracy of the $K_{\rm SMS}$6 and $K_{\rm SMS}$7 neutron levels. The mechanism is wave-function reshaping: $K_{\rm SMS}$8 states shrink, $K_{\rm SMS}$9 states broaden, and the mean-square radius of the $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$0 state increases by about $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$1 when switching from M3Y-P6 to M3Y-P6a (Nakada et al., 2014).

A further nuance concerns King-plot interpretations. One analysis of nobelium states that deformation can induce King-plot nonlinearity even without new particles (Allehabi et al., 2020), whereas the E120$$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$2 study finds that the deformations and central depressions examined there do not break King-plot linearity at a detectable level (Flambaum et al., 2019). Taken together, these results suggest that nuclear-structure-induced nonlinearity is system dependent rather than universal.

4. Heavy, superheavy, and astrophysical applications

Heavy atoms and ions provide some of the most striking nuclide-shift phenomena because relativistic and finite-size effects are simultaneously large. In $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$3, the 402.0 nm resonance line has a positive isotope shift of $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$4 GHz relative to $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$5, whereas the 399.6 nm line shows the unexpected negative value $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$6 GHz. The explanation is a corrected classification of the $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$7 level: instead of being predominantly $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$8, the CI+all-order calculation finds it to be $$\delta \nu_k^{A,A'} = \left(\frac{M'-M}{MM'}\right)\Delta \widetilde K_{\rm MS} + F_k\,\delta\langle r^2\rangle^{A,A'}.$$9, $\Delta A=2$0, $\Delta A=2$1, and only $\Delta A=2$2. The resulting strong $\Delta A=2$3 and $\Delta A=2$4-like admixture gives a field-shift constant $\Delta A=2$5 and explains the negative transition isotope shift (Okhapkin et al., 2015).

For total electron binding energies, the heavy-element trend is steep. Closed-shell neutral atoms from Ne to Og and a hypothetical $\Delta A=2$6 have tabulated field-shift coefficients that increase from $\Delta A=2$7 a.u./fm$\Delta A=2$8 for Cd and $\Delta A=2$9 for Xe to $Z$0 for No, $Z$1 for Og, and $Z$2 for $Z$3. A simple interpolation law $Z$4 reproduces calculated field shifts between neighboring closed shells to better than about $Z$5, with the effective exponent growing from about $Z$6 near $Z$7 to about $Z$8 near $Z$9. Because deep inner shells dominate, neutral and singly charged systems differ by only a few percent, mainly when an outer $Z\approx 38$0 electron is removed (Dzuba et al., 29 Jul 2025).

Superheavy-isotope searches in astrophysical spectra use isotope shifts as a translation rule between laboratory isotopes and neutron-richer isotopes near the proposed $Z\approx 38$1 island of stability. The strongest optical transitions in No, Lr, Nh, Fl, and element 120 were analyzed with the approximation

$Z\approx 38$2

For $Z\approx 38$3No, using $Z\approx 38$4 from CI+MBPT shifts the laboratory line at $Z\approx 38$5 to $Z\approx 38$6 for the hypothetical $Z\approx 38$7No; using the experimental $Z\approx 38$8 gives $Z\approx 38$9. The element-120 transition $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$00 is predicted to have particularly large isotope-shift coefficients, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$01 analytically and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$02 from many-body theory (Dzuba et al., 2017).

Nobelium also illustrates how atomic shifts can test nuclear theory directly. For the measured $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$03 isotope shift between $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$04No and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$05No, the experimental value is $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$06. Comparison with CDFT-based charge densities was used to discriminate among functionals, to argue that deformation contributes to the interpretation, and to predict $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$07No–$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$08No shifts of roughly $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$09 to $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$10 for four transitions (Allehabi et al., 2020).

5. Decay-energy and reaction-distribution shifts

A distinct but related use of shift language appears when improved nuclide characterization alters decay observables. For double-beta decay of $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$11, a direct Penning-trap measurement with ISOLTRAP gave

$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$12

compared with the AME2003 value $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$13. The result is therefore almost $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$14 keV higher and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$15 times more precise. Since the phase space grows strongly with $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$16, the revised value increases the phase-space factors from $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$17 and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$18 at $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$19 keV to $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$20 and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$21 at $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$22 keV. Under the single-state-dominance hypothesis, the expected two-neutrino half-life was reevaluated as $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$23 yr (Fink et al., 2011).

In multinucleon transfer, “shift” can designate motion of the fragment distribution through nuclide space. The DNS model study of $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$24Xe+$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$25Pb and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$26Ni+$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$27Pt incorporates deuteron, triton, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$28He, and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$29 transfer directly into the master equations. The inclusion of cluster transfer is found to be favorable for fragment formation with increasing the transferring nucleons and to lead to a broad mass distribution. For $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$30Xe+$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$31Pb at $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$32 MeV, the isotopic cross sections of W, Os, Rn, and Fr are described as nicely consistent with the Argonne data, and new neutron-rich isotopes of W and Os are predicted with cross sections above $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$33 nb. In the $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$34Ni+$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$35Pt and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$36Ni+$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$37Pt systems, neutron-rich isotopic maxima are shifted away from $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$38-stability, with reported maxima of $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$39 mb for $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$40Pt, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$41 mb for $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$42Ir, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$43 mb for $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$44Os, and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$45 mb for $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$46Re (Feng, 2023).

An analogous redistribution appears in true ternary fission of $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$47Cf. In the almost sequential collinear picture, the yield

$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$48

is reshaped because the Coulomb field of the first emitted outer fragment lowers the second pre-scission barrier. In the $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$49 example, the presence of the heavy third fragment makes the residual barrier shallower by about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$50 MeV, and the configuration $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$51 lies about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$52 MeV below $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$53. The resulting probability of about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$54 yields heavy clusters such as $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$55Ni, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$56Ge, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$57Se, and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$58Kr together with fragments in the $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$59–$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$60 region, while the middle fragment is predicted to have very small velocity and therefore often evade detection (Tashkhodjaev et al., 2015).

6. Instrumental response and shifting experimental frontiers

Not every reported “shift” is intrinsic to a nuclide; some are detector-response effects caused by nuclide irradiation. Under prolonged irradiation by fission products from $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$61Cf, Si(Li) p-i-n, Si surface-barrier, and planar p$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$62n detectors all show a linear shift of the heavy-fragment and light-fragment peaks toward lower visible energies as dose increases. The dependence of peak position on exposure is described as well represented by linear functions with negative slopes, and the light-fragment peak shifts about twice as fast as the heavy-fragment peak across all investigated detectors and conditions. The shift rate depends strongly on detector type and electric-field strength, but not on irradiation temperature; the inferred operational lifetime for use in a neutron calibration source is $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$63 to $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$64 years depending on detector choice, with the planar detector identified as the most radiation-hard among those tested (Bazlov et al., 2023).

At a broader historical scale, the frontier of known nuclides has itself shifted across the nuclear chart as instrumentation and production methods changed. One review states that presently about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$65 different nuclei are known and about another $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$66–$$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$67 are predicted to exist. It also gives several more specific counts under different conventions: $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$68 nuclides in the 2012 comprehensive overview, $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$69 observed by the end of 2011, and $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$70 presently reported in the published literature. With about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$71 nuclei calculated to be bound against neutron or proton emission, and after subtracting regions considered experimentally inaccessible, the review estimates that about $$\Delta E_{\rm IS} = F\,\delta\langle r^2\rangle + G\,\delta\langle r^2\rangle^2 + \left(\frac{m_e}{M_1}-\frac{m_e}{M_2}\right)\left(K_{\rm NMS}+K_{\rm SMS}\right),$$72 nuclides remain to be discovered. The historical sequence of discovery moved from stable nuclides found by mass spectroscopy, to neutron-deficient and transuranium systems produced by fusion-evaporation, and then to neutron-rich systems reached by fragmentation and fission; in this sense, a “shift” of the discovery frontier is a recurrent structural feature of nuclear science rather than a single observable (Thoennessen, 2013).

Taken together, these usages show that nuclide shift is not a single invariant quantity but a family of shift phenomena linking nuclear identity to observables in atomic structure, nuclear decay, reaction dynamics, detector response, and the evolving map of the nuclide chart. The common principle is that changing the nuclide, or changing what is known about a nuclide, propagates into measurable displacements whose interpretation requires both accurate many-body theory and careful experimental definition.