---
title: 'Toroidal Nuclei: Dynamics and Structure'
url: https://www.emergentmind.com/topics/toroidal-nuclei
type: topic
---

# Toroidal Nuclei: Dynamics and Structure

Toroidal nuclei denote two related but distinct nuclear phenomena. In one usage, dominant in contemporary dipole spectroscopy, the term refers to excited states whose **transition current** forms a torus-like vortical flow: the toroidal dipole mode or toroidal dipole resonance. In this dynamical sense, the equilibrium density is not donut-shaped; what is toroidal is the current field in an excited \(1^-\) state. In a second usage, the term denotes nuclei whose **intrinsic density** itself acquires a ring-like topology, as in toroidal high-spin isomers, toroidal cluster states, or hyperheavy mean-field solutions. Modern work therefore treats toroidicity both as a property of nuclear current and as a property of nuclear density, with different operators, observables, and stability criteria attached to each case [2510.09157] [1902.10108].

## 1. Terminological scope and geometric content

In nuclear-structure physics, “toroidal nucleus” most often means a nucleus in a toroidal **excitation**, not a static torus. The review of the toroidal dipole mode makes this point explicitly: the phrase usually refers to states in which the transition current forms a torus-like flow pattern, rather than to a ground-state density shaped like a torus [2510.09157]. By contrast, studies of \(^{12}\)C, high-spin light nuclei, superheavy nuclei, and hyperheavy nuclei treat toroidicity as an intrinsic property of the density itself, characterized by a major radius \(R\), a minor radius \(d\), and a central void [1911.11918] [1701.06327].

| Sense of “toroidal nucleus” | Hallmark | Representative cases |
|---|---|---|
| Dynamical toroidicity | Toroidal transition current in an excited state | TDM/TDR in \(^{58}\)Ni, \(^{24}\)Mg, \(^{170}\)Yb |
| Static toroidicity | Toroidal intrinsic density distribution | \(^{12}\)C toroidal states, light high-spin isomers, superheavy and hyperheavy toroids |

The distinction is not merely semantic. Dynamical toroidicity is diagnosed through current transition densities, toroidal operators, and transverse electron-scattering form factors. Static toroidicity is diagnosed through intrinsic density distributions, shell structure at fixed \(R/d\), deformation-energy surfaces, and stability against fission or multifragmentation. A persistent misconception is to treat these two literatures as interchangeable. They are contiguous, but not identical: one concerns vortical nuclear response, the other equilibrium or metastable ring-shaped matter distributions.

## 2. Vortical electric modes and the toroidal dipole operator

The theoretical starting point for dynamical toroidicity is the multipole decomposition of the nuclear transition current \(\delta \mathbf{j}(\mathbf r)\). In the Chandrasekhar–Moffatt or Debye-potential form,
\[
\delta\mathbf j(\mathbf r)
=
\nabla \phi(\mathbf r)
+
\nabla\times\big(\mathbf r\,\psi(\mathbf r)\big)
+
\nabla\times\nabla\times\big(\mathbf r\,\chi(\mathbf r)\big),
\]
the gradient term is irrotational and underlies conventional electric multipoles, the single-curl term is vortical magnetic, and the double-curl term is vortical electric; its long-wavelength part generates the toroidal current [2510.09157]. This decomposition is the formal reason toroidal multipoles are treated as a “third family” of electromagnetic modes, distinct from ordinary electric and magnetic multipoles.

For isoscalar dipole motion, toroidal and compressional operators separate the vortical and irrotational parts of the E1 response. A standard current-space form of the toroidal operator is
\[
\hat M_{\text{tor}(E1\mu)}
=
-\frac{i}{2\sqrt{3}c}
\int d^3r\,
\hat{\mathbf j}_c(\mathbf r)\cdot
\left[
\frac{\sqrt{2}}{5}r^2\mathbf Y_{12\mu}
+
\bigl(r^2-\langle r^2\rangle_0\bigr)\mathbf Y_{10\mu}
\right],
\]
while the compressional partner depends on the different linear combination of the same \(j_{10}\) and \(j_{12}\) current components, or equivalently on a density operator weighted by \(r^3Y_{1\mu}\) after the continuity equation is used [1602.03326]. In this formulation, toroidal motion is explicitly tied to the curl of the current, compressional motion to its divergence.

The continuity equation,
\[
\frac{\partial \rho(\mathbf r,t)}{\partial t}
+
\nabla\cdot \mathbf j(\mathbf r,t)=0,
\]
is central to the distinction. Ordinary electric dipole motion is largely longitudinal and continuity-equation constrained. Toroidal motion is approximately transverse, with
\[
\nabla\cdot \mathbf j_{\text{toroidal}}(\mathbf r)\approx 0,
\]
so it contributes only weakly to density-based observables in the long-wavelength limit [2510.09157]. This is why toroidal E1 modes can carry strong current structure while remaining weak in ordinary \(B(E1)\).

This current-based definition also sharpened the debate over “nuclear vorticity.” A criterion based solely on the \(j_{12}\) or \(j_+\) current component was shown to be inadequate because both toroidal and compressional modes contain \(j_{10}\) and \(j_{12}\). The hydrodynamically consistent quantity is the transverse current or, operationally, the toroidal strength function itself [1602.03326]. In practical self-consistent Skyrme RPA and QRPA calculations, this requires retaining current-dependent and time-odd terms in the functional; the formalism of full Skyrme RPA was developed precisely to calculate electric, magnetic, vortical, toroidal, and compression transitions on the same footing [1510.01248].

## 3. Low-energy E1 response, pygmy strength, and the \(^{58}\)Ni evidence

A major development in the field is the reinterpretation of low-energy dipole strength. In systematic Skyrme-QRPA studies of \(^{40,48}\)Ca, \(^{58,72}\)Ni, \(^{90,100}\)Zr, and \(^{100,120,132}\)Sn, the lower part of the region often labeled the pygmy dipole resonance was found to be basically **isoscalar vortical toroidal motion with a minor irrotational fraction**, independently of whether obvious PDR strength exists in the standard E1 channel [1903.01348]. In \(^{132}\)Sn, for example, the \(6\text{–}10\) MeV window simultaneously shows PDR-like transition densities and clear toroidal current fields, implying that the familiar neutron-skin-against-core picture is, at minimum, incomplete [1602.03326].

The dynamical picture is therefore mixed. Transition densities in neutron-rich nuclei may indeed display a neutron-dominated surface hump, but current transition densities reveal that the underlying flow is predominantly toroidal. This does not invalidate earlier \(B(E1)\)-based PDR analyses; rather, it separates observables. Density-sensitive probes emphasize the irrotational fraction, while current-sensitive probes emphasize the toroidal one.

The clearest experimental case is \(^{58}\)Ni. Long-standing TU Darmstadt \((e,e')\) data had shown at least six dipole states in the \(5\text{–}10\) MeV region with unusually steep transverse form factors at backward angles. These states were first labeled M1, then reclassified as E1 after polarized-photon measurements, leaving the transverse response unexplained within a conventional irrotational picture [2510.09157]. The combined analysis of high-resolution photon, proton, and electron scattering identified low-lying \(1^-\) candidates for toroidal dipole excitation at \(6.03\), \(8.24\), and \(8.87\) MeV, in correspondence with QRPA toroidal states at \(6.18\), \(8.26\), and \(8.95\) MeV [2310.04736].

Electron scattering is decisive because, in plane-wave Born approximation,
\[
\frac{d^2\sigma}{d\Omega\,dE'}
=
\sigma_{\text{Mott}}
\left[
v_L |F_L(q)|^2 + v_T |F_T(q)|^2
\right],
\]
and at very large angles the longitudinal piece becomes negligible, leaving the cross section directly sensitive to the transverse current form factor [2510.09157]. In \(^{58}\)Ni, the large-angle \((e,e')\) slopes are reproduced only when the states are assigned a strong toroidal component. The same work proposed an additional discriminator: for toroidal states the relative sign
\[
\Pi=\mathrm{sign}(F^C/F^T)
\]
is \(+1\), whereas for GDR-like and compression-like states it is \(-1\), so \((e,e'\gamma)\) interference measurements could separate vortical from irrotational E1 motion [2310.04736].

The consequence is broader than one nucleus. The \(^{58}\)Ni case shows that weakly collective, low-lying electric dipole states can exhibit current patterns more commonly associated with magnetic or transverse motion. It also suggests that a fraction of historical “anomalous M1” assignments in backward-angle electron scattering may require reinterpretation in terms of E1 toroidal dynamics.

## 4. Deformation, anomalous \(K\)-splitting, and individual toroidal states

Axial deformation reorganizes toroidal response in characteristic ways. In prolate \(^{170}\)Yb, Skyrme-RPA calculations showed that the low-energy toroidal region lies at \(5\text{–}20\) MeV and exhibits an **anomalous branch ordering**:
\[
E_{\text{tor}}(\mu=1) < E_{\text{tor}}(\mu=0),
\]
opposite to the normal prolate ordering of the GDR and compression mode, for which \(E(\mu=0) < E(\mu=1)\) [1311.4366]. The effect appears already at the unperturbed two-quasiparticle level and was later emphasized as a robust deformation fingerprint of the toroidal dipole resonance in prolate nuclei [1602.03326].

In light strongly deformed nuclei, deformation can isolate **individual toroidal states** rather than merely reshuffle a broad resonance. The clearest example is \(^{24}\)Mg, where axial QRPA with SLy6 predicts the lowest \(I^\pi K=1^-1\) excitation at \(E=7.92\) MeV to be a vortical toroidal state. Its current transition density forms a vortex–antivortex realization of Hill’s spherical vortex in strong axial confinement, and the result persists for SLy6, SVbas, and SkM* [1711.08953]. In the same nucleus, the nearby \(K=0\) state at \(9.56\) MeV is compression-dominated, making \(^{24}\)Mg a textbook case of side-by-side vortical and irrotational low-energy E1 motion.

The detailed spectroscopy of light deformed nuclei also showed that the lowest toroidal state is **not universal**. A comparative QRPA study of \(^{24}\)Mg and \(^{20}\)Ne found that the lowest toroidal \(K=1\) state is a peculiarity of \(^{24}\)Mg. In \(^{20}\)Ne, the toroidal state appears at \(10.11\) MeV, whereas the \(K=0\) compression state lies lower, at \(7.91\) MeV [1809.01097]. The difference was traced to deformation-driven rearrangements of specific Nilsson orbitals near the Fermi surface. This suggests that strong prolate deformation is necessary but not sufficient; the detailed single-particle spectrum matters.

These light-nucleus results also connect toroidicity to clustering. In \(^{24}\)Mg, the toroidal \(7.92\) MeV state and the compression \(9.56\) MeV state both lie near the \(\alpha\)-particle threshold \(S_\alpha=9.3\) MeV, and the density/current patterns were interpreted as intertwined with cluster structure [1711.08953]. This does not imply that toroidal states are reducible to cluster states, but it does indicate that in light nuclei the two descriptions can overlap in the same energy domain.

## 5. Static toroidal densities, shell structure, and high-spin toroidal isomers

Static toroidicity enters nuclear structure in several forms. In \(^{12}\)C, planar intrinsic \(K=0\) and \(K=I\) states were constructed both from rotated triangular \(3\alpha\) cluster configurations and from toroidal shell-model Slater determinants. In that framework, the ground state was found to contain a toroidal core, and a toroidal mean-field configuration with
\[
R=1.16~\text{fm},\qquad d=1.37~\text{fm},\qquad a_2=-0.007
\]
occurs at the Hoyle energy \(E_x=7.65\) MeV, although it is not a local minimum in \(R\) within a single-determinant treatment [1911.11918]. The significance of this result is not a definitive toroidal assignment of the Hoyle state, but the demonstration that toroidal density profiles and \(3\alpha\) generator-coordinate constructions can be placed within one framework.

A different route to static toroidicity is high spin. Self-consistent Skyrme-Hartree-Fock calculations predicted even-even toroidal high-spin isomers in light nuclei with \(28\le A\le 52\), including \(N\neq Z\) systems. Explicit examples are \(^{36}\mathrm S(I=74\hbar)\) and \(^{40}\mathrm Ar(I=80,102\hbar)\), which fall on the same regular multi-particle–multi-hole patterns as the \(N=Z\) toroidal high-spin isomers [1412.0050]. In these systems, aligned particle–hole excitations create large \(I_z\), and the toroidal density is stabilized over a finite angular-momentum window.

The shell-model underpinning of such states was generalized to the intermediate-mass region by a toroidal single-particle potential study. That work found large toroidal shell gaps at various nucleon numbers and aspect ratios \(R/d\), showed how Bohr–Mottelson spin-aligning particle–hole excitations can generate toroidal high-spin isomers, and introduced the possibility of **toroidal vortex nuclei**, in which particle–hole excitations change a nucleon’s vorticity quantum number [1807.11646]. This suggests that static toroidal densities support their own shell structure, high-spin spectroscopy, and even intrinsic vorticity labels.

In superheavy nuclei, the same logic reappears at much larger scales. For \(^{304}{120}_{184}\), strongly deformed oblate configurations bifurcate into a toroidal branch, and cranked Skyrme-Hartree-Fock predicts toroidal high-spin isomers at
\[
I_z=81\hbar \quad \text{and}\quad I_z=208\hbar,
\]
with \(Q_{20}\) near \(-300\) b [1705.01408]. Yet zero-spin calculations for \(Z=120\) isotopes and \(N=184\) isotones showed that, although toroidal solutions become lower in energy than biconcave discs beyond sufficiently negative \(Q_{20}\), the toroidal branch itself still lacks a local minimum at those proton and neutron numbers [1701.06327]. High spin therefore acts as a stabilizer where shell structure and Coulomb effects alone are not yet sufficient.

## 6. Hyperheavy toroids, shell stabilization, and the extended nuclear landscape

The most systematic static-toroidicity results come from covariant density functional theory in the hyperheavy region. Axial RHB calculations predict that beyond about \(Z\sim 130\), the energetically favored axial shapes change qualitatively from ellipsoidal-like to toroidal ones, and in the \(Z=140\text{–}180\) region the lowest-energy axial solutions carry very large negative \(\beta_2\) and ring-shaped densities [1902.10108]. The physical driver is the Coulomb-energy gain from spreading charge into a torus: the reduction in \(E_{\text{Coul}}\) between spherical and toroidal configurations rises from about \(346\) MeV in \(^{208}\)Pb to about \(721\) MeV in \(^{354}_{134}\), \(874\) MeV in \(^{466}_{156}\), and \(1126\) MeV in \(^{426}_{176}\mathrm X\) [1902.10108].

Subsequent CDFT systematics refined this picture. Toroidal hyperheavy even-even nuclei were found to be **stable with respect to breathing deformations**; the most compact fat tori lie near \(Z\approx 136\), \(N\approx 206\), whereas thin toroidal nuclei become dominant with increasing proton number and toward proton and neutron drip lines [2012.13799]. The lowest axial toroidal solutions are characterized either by large shell gaps or by low single-particle level density near the Fermi level in at least one subsystem, and the \(N=210\) toroidal shell gap was identified as an important stabilizer of fat toroidal nuclei [2012.13799].

The stability problem, however, is not closed. The same hyperheavy studies emphasize that many toroidal nuclei are expected to be unstable toward multifragmentation or sausage-type distortions, even when they are stable against breathing motion [1902.10108]. In parallel, spherical shell closures at \(Z=154,186\) and \(N=228,308,406\) define competing islands of spherical hyperheavy stability [2012.13799]. The global landscape is therefore a competition between Coulomb-driven toroidal axial minima and shell-driven spherical minima, with the final answer depending on the balance among shell structure, triaxiality, and non-axial instabilities.

A recent extension in EQMD has added a morphological classification of **bubble-like** nuclei using the dimensionless parameters \(BHTU\). In that scheme, \(B=2\) defines “toroidal bubble nuclei,” identified from two inflection points in the radial density profile; these nuclei emerge for \(Z\approx 25\) and become prevalent in heavy systems [2605.24676]. This is not the same concept as the macroscopic toroidal mean-field shapes of hyperheavy CDFT, but it shows that toroidal-like hollows now enter nuclear taxonomy through both self-consistent mean-field and cluster-based descriptions.

Taken together, the literature establishes toroidicity as a genuine organizing principle of nuclear structure rather than an isolated curiosity. In the dynamical sector, it identifies a vortical electric mode beyond the continuity-equation description of ordinary E1 motion and now supported experimentally in \(^{58}\)Ni. In the static sector, it links ring-shaped densities to cluster geometry, aligned high-spin configurations, shell structure in intermediate-mass nuclei, and the Coulomb-dominated shape evolution of superheavy and hyperheavy matter. The term “toroidal nucleus” therefore names not one object but a family of nuclear states in which either the **current** or the **density** acquires toroidal topology.

Source: https://www.emergentmind.com/topics/toroidal-nuclei