Superheavy Element Physics

• Superheavy element physics is the study of nuclei beyond Z=104, characterized by extreme Coulomb repulsion, unique shell effects, and relativistic atomic behavior.
• Researchers employ heavy-ion fusion-evaporation and multi-nucleon transfer methods alongside advanced models to predict synthesis cross sections and decay mechanisms.
• Insights from this field enhance understanding of atomic structure, astrophysical r-process nucleosynthesis, and exotic nuclear configurations like bubble structures.

Superheavy element (SHE) physics encompasses the synthesis, structure, stability, and properties of nuclei beyond the actinide series, typically defined as those with atomic number $Z>104$. This regime probes the limits of nuclear binding under extreme Coulomb repulsion, the emergence of new shell closures (“islands of stability”), high-$Z$ atomic and chemical phenomena governed by strong relativistic effects, and the mechanisms by which these nuclei form both in the laboratory and astrophysical contexts. Experimental access to SHEs is primarily via complete-fusion reactions in heavy-ion collisions, with successful synthesis and identification dictated by cross sections in the picobarn to femtobarn range and by the identification of characteristic decay chains. Theoretical approaches integrate advanced nuclear mean-field models, macroscopic-microscopic shell corrections, high-precision atomic structure calculations, and dynamical models for reaction pathways. The field sits at the intersection of nuclear structure, atomic theory, and astrophysics.

1. Synthesis and Reaction Mechanisms

The primary route to laboratory SHE synthesis is low-energy heavy-ion fusion-evaporation, where a projectile and a target combine to form a compound nucleus (CN) which cools via neutron evaporation, producing an evaporation residue (ER). The ER cross section is decomposed as:

$$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$

• $\sigma_{\rm cap}(E_{\rm cm}, J)$: capture cross section for forming the dinuclear system.
• $P_{\rm CN}(E_{\rm cm}, J)$: probability of evolving past quasifission to a fully equilibrated CN.
• $W_{\rm sur}(E_{\rm cm}, J)$: survival probability against fission for the hot CN de-exciting by neutron emission.

The interplay between the entrance-channel properties—charge/mass asymmetry, mean fissility ($\chi_m$), and nuclear deformation—and the dynamics of multinucleon transfer, shell effects, and orientation angles controls the probability of forming and observing an evaporation residue. Systematic analysis has demonstrated sharply peaked excitation functions: optimal beam energies must be within 3–5 MeV of a narrow maximum for measurable cross sections (e.g. $\sim$12.3 fb for $Z=119$ with $^{51}$V+$Z$0Cm at $Z$1 MeV and optimal $Z$2–$Z$3) (Nasirov et al., 2023, Manjunatha et al., 2021, Manjunatha et al., 2021).

Table 1: Sample cross sections, $Z$4 synthesis (DNS model, $Z$5V+$Z$6Cm) (Nasirov et al., 2023)

Channel	$Z$7 (MeV)	$Z$8 (fb)	Optimal orientation $Z$9 (deg)
4n	232	12.3	60–70
3n	225	6.6	60–70

Further increases in $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$0 or the use of heavier projectiles (e.g. 50Ti, 54Cr) on actinide targets result in drastic cross section suppression due to enhanced quasifission probabilities and reduced survival (III et al., 2012, Nasirov et al., 2013). The use of $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$1Ca has proven particularly favorable, but new superheavy elements will require projectiles such as Ti or Cr, with cross sections in the 10–100 fb range for $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$2 (III et al., 2012).

Alternative routes, such as multi-nucleon transfer and neutron capture processes, have been proposed for populating neutron-rich SHEs, but are experimentally challenging due to small yields and short half-lives (Zagrebaev et al., 2012).

2. Shell Structure, Deformation, and the Island of Stability

The stability of SHEs is dictated by quantal shell corrections superimposed on the macroscopic trends of the liquid-drop model. Closed (magic) shells provide localized stabilization against spontaneous fission. Predicted major closures include $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$3 and $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$4, with the expected “island of stability” centered in this vicinity (Ramirez et al., 2014, Hessberger, 2021, Stone et al., 2019).

Penning-trap mass spectrometry and high-precision binding-energy measurements directly pin down the strength and location of neutron shell gaps: for example, a deformed neutron shell gap of $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$5 MeV at $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$6 has been established in nobelium, providing a benchmark for global models (Ramirez et al., 2014). Self-consistent mean-field models (e.g. QMC, Skyrme-HFB, RMF) predict pronounced deformed gaps at $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$7 and $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$8 in $$\sigma_{\rm ER}(E_{\rm cm}) = \sum_{J=0}^{J_{\rm max}} (2J{+}1)\; \sigma_{\rm cap}(E_{\rm cm}, J) P_{\rm CN}(E_{\rm cm}, J) W_{\rm sur}(E_{\rm cm}, J)$$9–$\sigma_{\rm cap}(E_{\rm cm}, J)$0 nuclei, but shape transitions (prolate–oblate–spherical) at $\sigma_{\rm cap}(E_{\rm cm}, J)$1 and the possible absence or smearing of the $\sigma_{\rm cap}(E_{\rm cm}, J)$2 closure in higher $\sigma_{\rm cap}(E_{\rm cm}, J)$3 isotopes (Stone et al., 2019). Bulk observables (e.g. $\sigma_{\rm cap}(E_{\rm cm}, J)$4, $\sigma_{\rm cap}(E_{\rm cm}, J)$5) only weakly reflect sub-MeV shell gaps, so single-particle spectra and neutron pairing energies are critical for gap identification.

Isotope shift measurements and extraction of nuclear radii in superheavy atoms further constrain shell-model and mean-field predictions and provide guidance for future searches (Lackenby et al., 2018).

3. Alpha Decay, Spontaneous Fission, and Decay Chains

The dominant decay mode for most known SHEs is $\sigma_{\rm cap}(E_{\rm cm}, J)$6 emission, with spontaneous fission branching ratios growing for more neutron-deficient or higher-$\sigma_{\rm cap}(E_{\rm cm}, J)$7 isotopes. Systematic measurements of $\sigma_{\rm cap}(E_{\rm cm}, J)$8 values and half-lives display signatures of shell stabilization (suppressed $\sigma_{\rm cap}(E_{\rm cm}, J)$9, increased $P_{\rm CN}(E_{\rm cm}, J)$0) at known and predicted magic numbers ($P_{\rm CN}(E_{\rm cm}, J)$1, $P_{\rm CN}(E_{\rm cm}, J)$2), but clear evidence for the $P_{\rm CN}(E_{\rm cm}, J)$3 closure remains elusive (Hessberger, 2021, Kumar, 2011).

Theoretical $P_{\rm CN}(E_{\rm cm}, J)$4-decay half-lives are calculated via barrier-penetration models which account for nuclear deformation, pairing, and preformation probabilities. A generic expression is

$P_{\rm CN}(E_{\rm cm}, J)$5

where $P_{\rm CN}(E_{\rm cm}, J)$6 is the cluster preformation probability, $P_{\rm CN}(E_{\rm cm}, J)$7 the assault frequency, and $P_{\rm CN}(E_{\rm cm}, J)$8 the WKB penetrability. Isospin-dependent radii and improved proximity potentials increase the accuracy of models such as the Isospin Cluster Model, reproducing decay properties and highlighting enhanced stability near $P_{\rm CN}(E_{\rm cm}, J)$9 (Kumar, 2011).

Spontaneous fission lifetimes vary over many orders of magnitude within the superheavy landscape, complicating identification of decay chains. Accurate theoretical predictions require microscopic treatments of collective inertia, shell effects, and pairing correlations (Hessberger, 2021). Ambiguities in chain assignments and the presence of isomerism or alternative (e.g. EC) decay paths remain major sources of uncertainty.

4. Atomic Physics of Superheavy Elements

High-$W_{\rm sur}(E_{\rm cm}, J)$0 atoms experience strong relativistic and QED effects, radically altering their electronic structure. Solving the Dirac–Coulomb–Breit–QED Hamiltonian is essential for predicting atomic spectra, ground-state configurations, and parity nonconservation phenomena (Smits et al., 2023, Lackenby et al., 2018).

Configuration-interaction plus perturbation-theory (CIPT) and coupled-cluster approaches allow explicit calculation of level energies, transition amplitudes, and isotope shifts, exemplified by work on dubnium ($W_{\rm sur}(E_{\rm cm}, J)$1) and other 6d elements. Key features include:

• Substantial level mixing from strong spin-orbit splitting and relativistic contraction.
• QED radiative corrections (self-energy, vacuum polarization) introducing $W_{\rm sur}(E_{\rm cm}, J)$2 shifts in binding energies and line positions.
• Critical nuclear charge, $W_{\rm sur}(E_{\rm cm}, J)$3, beyond which 1s levels dive into the negative-energy continuum, leading to Gamow resonances and new theoretical challenges for multi-electron systems (Smits et al., 2023).

Atomic structure calculations guide experimental laser spectroscopy campaigns, interpretation of chemistry, and extraction of nuclear properties (e.g. charge radii via isotope shifts).

5. Astrophysical Context and the r-Process

The synthesis of superheavy nuclides in astrophysical environments occurs, if at all, in extremely neutron-rich, high-entropy, rapid-capture (r-process) conditions. Neutron-star mergers, magnetorotational supernovae, and collapsars are candidate sites. Key nuclear physics inputs are nuclear masses, separation energies, β-decay and fission rates, and fission fragment distributions (Holmbeck et al., 2023, Petermann et al., 2012).

Dynamical network calculations demonstrate that depending on fission barrier heights, the r-process may reach $W_{\rm sur}(E_{\rm cm}, J)$4 ($W_{\rm sur}(E_{\rm cm}, J)$5) during freeze-out, with survival limited by fission recycling and short spontaneous fission half-lives ($W_{\rm sur}(E_{\rm cm}, J)$6 hours to days). The actual extent is highly sensitive to uncertainties in mass models (FRDM, ETFSI, HFB), barrier predictions, and astrophysical parameters (Holmbeck et al., 2023, Petermann et al., 2012). Observational constraints from meteoritic actinides, metal-poor stars, and kilonova light curves provide indirect access to the existence and properties of the heaviest r-process nuclei.

6. Exotic Structures and Macroscopic Properties

Predictions for exotic density distributions, such as hollow (“bubble” or “fullerene-like”) configurations, appear for doubly magic or near-magic nuclei ($W_{\rm sur}(E_{\rm cm}, J)$7, $W_{\rm sur}(E_{\rm cm}, J)$8) in both macroscopic liquid-drop and microscopic soliton approaches. These models yield central density depletion, altered giant resonance spectra, and prospective metastable minima stabilized by shell corrections (Misicu et al., 2018).

The bulk mass density of SHE matter is predicted via relativistic Thomas–Fermi theory to reach 36–68 g/cm$W_{\rm sur}(E_{\rm cm}, J)$9 at $\chi_m$0=164, far exceeding terrestrial elements (LaForge et al., 2023). The presence of “shared” conduction electrons and strong relativistic contraction of atomic radii affect the material properties. Compressed, superheavy atoms, potentially forming in neutron-star crusts or ultradense astrophysical environments, exhibit nuclear surface and Coulomb energies closely tied to electron penetration and background densities, influencing stability against fission-like deformations (Rueda et al., 2017, LaForge et al., 2023).

7. Prospects, Uncertainties, and Future Directions

The search for new SHEs beyond $\chi_m$1 is constrained by minuscule cross sections, target material availability, and requirements for long beam times and sensitive detection schemes. Accurate predictions of optimal beam energies, fusion probabilities, and survival factors are critical (Manjunatha et al., 2021, Nasirov et al., 2023). Multi-nucleon transfer and neutron-capture approaches remain underexplored but may provide access to more neutron-rich (and possibly more stable) superheavies (Zagrebaev et al., 2012).

On the theory side, resolving the location and magnitude of major shell closures ($\chi_m$2 vs $\chi_m$3, $\chi_m$4 vs $\chi_m$5), improving microscopic treatments of fission barriers and pairing, and achieving tighter connections between atomic, nuclear, and astrophysical observations are current objectives. The possible discovery of fullerene-like superheavy nuclei, extreme-density materials, and observational evidence for superheavy nucleosynthesis in the cosmos represent the next frontiers for superheavy element physics.