---
title: Static Hexadecapole Deformation in Nuclear Shapes
url: https://www.emergentmind.com/topics/static-hexadecapole-deformation
type: topic
---

# Static Hexadecapole Deformation in Nuclear Shapes

Static hexadecapole deformation is the intrinsic, equilibrium component of the $\lambda=4$ degree of freedom in nuclear shape, conventionally denoted by $\beta_4$. In the standard surface expansion for an axially symmetric nucleus, it enters through the $Y_{40}$ term in
$$
R(\theta,\phi)=R_0\left[1+\sum_\lambda \beta_\lambda Y_{\lambda 0}(\theta,\phi)\right],
$$
or, when only the leading even multipoles are retained,
$$
R(\theta)=R_0\big[1+\beta_2Y_{20}(\theta)+\beta_4Y_{40}(\theta)+\cdots\big].
$$
Here $\beta_2$ controls the dominant quadrupole elongation or flattening, whereas $\beta_4$ modifies finer surface curvature, including sharpening or flattening of the ends and equatorial pinching, typically on top of a quadrupole background [1811.12756]. In contemporary usage, “static” distinguishes a nonzero equilibrium $\beta_4$ at the ground-state or mean-field minimum from dynamic hexadecapole correlations associated with vibrations, softness, or collective fluctuations [2312.01593][2401.06117].

## 1. Formal definition and shape parameterizations

Static hexadecapole deformation is defined through the nonzero axial $\lambda=4$ component of the intrinsic nuclear shape. In the axial limit, the relevant spherical harmonics are
$$
Y_{20}(\theta)=\sqrt{\frac{5}{16\pi}}(3\cos^2\theta-1),\qquad
Y_{40}(\theta)=\frac{3}{16\sqrt{\pi}}[35\cos^4\theta-30\cos^2\theta+3],
$$
with equivalent normalizations and notational variants used across reaction, mean-field, and heavy-ion initial-state studies [1811.12756][2603.24088][2312.01593]. For axial shapes in the notation of the spherical multipole expansion,
$$
R(\theta,\phi)=R_0\left[1+\sum_{\lambda\mu}\alpha_{\lambda\mu}Y_{\lambda\mu}(\theta,\phi)\right],\qquad
\beta_\lambda=\sqrt{4\pi}\,\alpha_{\lambda 0},
$$
so $\beta_4$ is the axial projection of the $\lambda=4$ distortion [2401.06117].

The deformation parameter may also be defined from intrinsic moments. In self-consistent mean-field and related Gogny or relativistic calculations, the axial moments $Q_{20}$ and $Q_{40}$ are mapped to dimensionless deformations by conventions of the form
$$
\beta_\lambda=\frac{4\pi}{3AR^\lambda}\langle \hat Q_{\lambda 0}\rangle,\qquad R=1.2A^{1/3}\ \mathrm{fm},
$$
or equivalent variants with $R_0=1.2A^{1/3}$ [2312.01593][2502.04985][2508.05268][2403.03393]. The microscopic hexadecapole moment is correspondingly
$$
Q_{40}=\int \rho(r)\,r^4Y_{40}(\theta,\phi)\,d^3r
$$
in the axial case [2312.01593].

Several works emphasize that deformation parameters depend on the representation. In deformed Woods–Saxon densities,
$$
\rho(r,\theta)=\frac{\rho_0}{1+\exp((r-R(\theta))/a)},
$$
the surface parameters $\beta_\lambda^{WS}$ are not identical to the volume multipole moments $\beta_{\lambda m}$ extracted from the density. For $^{238}$U, the mapping is nonlinear, and the volume quadrupole moment receives a contribution from the surface hexadecapole term:
$$
\beta_{20}=\left(\frac{R_d^2}{R_0^2}\right)\left[\beta^{WS}_{20}+\frac{2}{7}\sqrt{\frac{5}{\pi}}(\beta^{WS}_{20})^2+\frac{12}{7\sqrt{\pi}}\beta^{WS}_{20}\beta^{WS}_{40}\right].
$$
This distinction is central in relativistic heavy-ion implementations of deformed nuclei [2302.13617].

A separate but related convention arises in the Pt–Hg–Pb macroscopic–microscopic study using a rapidly converging Fourier shape parametrization. There the coordinate $q_4$ is identified as the hexadecapole-like degree of freedom, but no explicit closed-form conversion $q_4\to\beta_4$ is given [2001.09997]. This is a reminder that “static hexadecapole deformation” is model-independent as a concept, whereas the numerical deformation parameter depends on the chosen shape representation.

## 2. Static deformation, softness, and collective diagnostics

A nonzero $\beta_4$ at an energy minimum is the direct signature of static hexadecapole deformation in self-consistent or macroscopic–microscopic energy surfaces. In the axially deformed Gd isotopes studied with constrained relativistic mean field mapped onto the $sdg$ interacting-boson model, the global minima occur at
$(\beta_2^{\min},\beta_4^{\min})\approx (-0.05,0)$, $(0.15,0.05)$, $(0.2,0.1)$, $(0.3,0.15)$, $(0.3,0.15)$, $(0.35,0.15)$, and $(0.35,0.15)$ for $^{148-160}$Gd, respectively, so static hexadecapole deformation is present from $^{150}$Gd onward, whereas $^{148}$Gd is soft in $\beta_4$ but has $\beta_4^{\min}=0$ [2312.01593].

The distinction between static deformation and softness is treated systematically in spherical HFBCS+QRPA studies. There, the spherical solution is diagnosed as unstable in a given multipolarity when any QRPA eigenfrequency becomes imaginary,
$$
\omega_\nu^2<0,
$$
which corresponds, in a harmonic picture,
$$
E(\beta_\lambda)\approx E_0+\frac{1}{2}C_\lambda\beta_\lambda^2,
$$
to negative curvature $C_\lambda<0$ [2401.06117]. In nuclei without collapse, hexadecapole softness is quantified by the inverse-energy-weighted sum rule
$$
m_{-1}(F_4)=\sum_\nu \frac{|\langle 0|\hat F_4|\nu\rangle|^2}{\omega_\nu},
$$
and by the polarizability
$$
\mathcal C_4=\frac{2m_{-1}(4)}{A},
$$
with larger $\mathcal C_4$ indicating weaker stiffness against $\lambda=4$ distortion [2401.06117]. That framework identifies $\lambda=4$ collapse mainly around neodymium and polonium, and stresses that the method does not determine the sign of $\beta_4$ [2401.06117].

Beyond-mean-field quadrupole–hexadecapole coupling has been quantified with two-dimensional GCM in both actinides and rare earths. In Ra–Pu isotopes, constrained Gogny-HFB plus 2D-GCM finds large positive static $\beta_4$ around $^{238}$U at the HFB and 2D-GCM levels, but also a dynamically stable region with weak negative $\beta_4$ just below $N=184$ [2502.04985]. In Yb, Hf, W, and Os, HFB curvature analysis and 2D-GCM collective wave functions show that quadrupole and hexadecapole degrees of freedom are interwoven up to approximately $A=184$–$188$, and that a square-like region with $\beta_4<0$ persists below $N=126$ after zero-point fluctuations are included [2508.05268].

These studies also connect static $\beta_4$ to spectroscopy. In the $sdg$-IBM and mapped bosonic descriptions, explicit inclusion of the $g$ boson improves high-spin yrast states near shell closure and produces $K^\pi=4^+$ bands with strong $E4$ transitions in deformed nuclei [2312.01593][2403.03393]. A plausible implication is that static hexadecapole deformation is best viewed not as an isolated parameter but as one component of a coupled even-multipole geometry, especially in regions where $\beta_2$ and $\beta_4$ are strongly correlated.

## 3. Experimental determination and extraction strategies

Backward-angle quasi-elastic scattering near the Coulomb barrier has become a high-sensitivity probe of static $\beta_4$ in light nuclei. In this method, the quasi-elastic barrier distribution is obtained from
$$
D_{\mathrm{qel}}(E_{\mathrm{eff}})=-\frac{d}{dE_{\mathrm{eff}}}\left[\frac{\sigma_{\mathrm{qel}}(E_{\mathrm{eff}})}{\sigma_R(E_{\mathrm{eff}})}\right],
$$
with angle-dependent effective energy mapped through
$$
E_{\mathrm{eff}}=\frac{2E_{\mathrm{c.m.}}}{1+\cosec(\theta_{\mathrm{c.m.}}/2)}.
$$
Coupled-channels analyses with modified CCFULL then compare measured excitation functions and barrier distributions to rotor-plus-phonon calculations [1811.12756].

For $^{24}$Mg, quasi-elastic scattering on $^{90}$Zr together with Bayesian analysis yielded
$\beta_2=+0.43\pm0.02$ and $\beta_4=-0.11\pm0.02$ at 95% confidence, with a moderate anticorrelation between the parameters and a clearly identified negative hexadecapole deformation [1811.12756]. For $^{28}$Si+$^{90}$Zr, the analogous analysis found
$\beta_2=-0.38\pm0.01$ and $\beta_4=+0.03\pm0.01$, with the oblate solution decisively favored over a prolate alternative [2303.12495]. Earlier fusion-barrier-distribution analysis of $^{28}$Si+$^{92}$Zr had already shown that the $2^+$ reorientation term,
$$
O_{22}=\langle Y_{20}|\beta_2R_PY_{20}+\beta_4R_PY_{40}|Y_{20}\rangle,
$$
is strongly sensitive to the sign and value of $\beta_4$, and that for $\beta_2=-0.407$ and $\beta_4=+0.25$ the quadrupole and hexadecapole contributions nearly cancel [1805.03395].

Inverse-kinematics inelastic proton scattering provides another route. For $^{74,76}$Kr, coupled-channels fits to the $4_1^+$ cross sections gave two possible $\beta_4$ solutions because the measured cross sections are insensitive to the sign of $\beta_4$:
for $^{76}$Kr, $\beta_4=+0.201\pm0.009\text{ (stat.)}\pm0.016\text{ (sys.)}$ or $\beta_4=-0.127\pm0.009\text{ (stat.)}\pm0.022\text{ (sys.)}$;
for $^{74}$Kr, $\beta_4=+0.23\pm0.02\text{ (stat.)}\pm0.02\text{ (sys.)}$ or $\beta_4=-0.17\pm0.02\text{ (stat.)}\pm0.02\text{ (sys.)}$. Comparison to non-relativistic and relativistic EDF calculations favored the large positive solutions and linked them to well-deformed prolate configurations [2304.14246].

A concise set of representative nucleus-specific values illustrates the diversity of extracted or predicted static $\beta_4$:

| Nucleus/system | Static hexadecapole result | Source |
|---|---:|---|
| $^{24}$Mg | $\beta_4=-0.11\pm0.02$ | [1811.12756] |
| $^{28}$Si | $\beta_4=+0.03\pm0.01$ | [2303.12495] |
| $^{76}$Kr | $\beta_4=+0.201\pm0.009\pm0.016$ or $-0.127\pm0.009\pm0.022$ | [2304.14246] |
| $^{74}$Kr | $\beta_4=+0.23\pm0.02\pm0.02$ or $-0.17\pm0.02\pm0.02$ | [2304.14246] |
| $^{238}$U | HFB $\beta_4=0.16$; 2D-GCM ground state $\beta_4=0.16$ | [2502.04985] |

A recurring methodological point is that barrier distributions or one-step $\lambda=4$ cross sections are often more discriminating than raw excitation functions. In the $^{24}$Mg and $^{28}$Si quasi-elastic studies, the barrier distribution was described as much more sensitive to structural couplings than the excitation function itself [1811.12756][2303.12495].

## 4. Nuclear-structure systematics across mass regions

Static hexadecapole deformation is not confined to one part of the nuclide chart. In the $sd$ shell, quasi-elastic analyses establish opposite-sign examples: $^{24}$Mg is strongly prolate with negative $\beta_4$, whereas $^{28}$Si is oblate with small positive $\beta_4$ [1811.12756][2303.12495]. These two cases are frequently treated as benchmarks for the sign sensitivity of barrier distributions.

In rare-earth nuclei near $N=90$, mapped $sdg$-IBM and relativistic mean-field studies report nonzero equilibrium $\beta_4$ that increases with neutron number. The mean-field minima reach approximately $\beta_4^{\min}\approx0.25$ in lighter rare earths such as Nd and Sm around $N\approx90$–$94$, and $\beta_4^{\min}\approx0.15$ in Gd, Dy, and Er for $N\ge 90$ [2403.03393]. The same region is also highlighted in spherical QRPA as one where large $\mathcal C_4$ and, in some nuclei, $\lambda=4$ collapse occur, especially around neodymium [2401.06117].

In Gd isotopes, constrained relativistic calculations explicitly show the onset of static $\beta_4$ in the ground-state minimum from $^{150}$Gd onward, with $\beta_4^{\min}=0.05$ in $^{150}$Gd, $0.10$ in $^{152}$Gd, and $0.15$ in $^{154-160}$Gd [2312.01593]. In a broader rare-earth set, Gogny HFB and 2D-GCM find a structural evolution from diamond-like shapes with $\beta_4>0$ in lighter isotopes to square-like shapes with $\beta_4<0$ below $N=126$; for the 2D-GCM ground states, the square-like region satisfies
$$
-0.07\le \beta_{4,2D\text{-}GCM}^{\sigma=1}\le -0.01
$$
for selected Yb, Hf, W, and Os isotopes [2508.05268].

In the A$\approx180$ region, macroscopic–microscopic Woods–Saxon plus HFBC cranking calculations identify an “island of negative axial $\beta_4$” in $^{180-184}$Yb, $^{182-186}$Hf, and $^{184-188}$W, with equilibrium values spanning roughly $-0.07$ to $-0.09$ and with observable consequences for moments of inertia [2601.10052]. In neutron-rich Zr, Skyrme-HFB predicts that $\beta_4$ is suddenly enhanced at $N=60$, remains sizable and positive across $N=60$–$74$, shows a kink around $N\approx70$, and collapses at $N=76$ when the shape changes to oblate [2304.06238]. The microscopic driver is identified as occupation of the intruder Nilsson orbits $[550]1/2$ and then $[530]1/2$ [2304.06238].

Actinides display another characteristic pattern. Gogny HFB and 2D-GCM calculations for Ra, Th, U, and Pu find HFB ground-state $\beta_4$ values around $^{238}$U of $0.09$, $0.13$, $0.16$, and $0.17$ for $^{238}$Ra, $^{238}$Th, $^{238}$U, and $^{238}$Pu, respectively, with 2D-GCM values remaining close to the HFB results [2502.04985]. With increasing mass number, $\beta_4$ decreases toward zero and becomes weakly negative just below $N=184$, yet those small negative values remain dynamically stable under quadrupole–hexadecapole configuration mixing [2502.04985].

These regional patterns show that both positive and negative static $\beta_4$ occur, often in correlation with shell structure, intruder occupation, and the underlying quadrupole background. This suggests that there is no universal sign rule for $\beta_4$; its sign is nucleus-specific and strongly model- and region-dependent.

## 5. Spectroscopy, rotational dynamics, and reaction observables

Static hexadecapole deformation affects spectra, transition rates, and rotational response. In the $sdg$-IBM description of axially deformed Gd isotopes, inclusion of the $g$ boson does not qualitatively alter most low-spin low-lying states, but it lowers calculated excitation energies of ground-band states with $I^\pi\ge 6^+$ in nuclei with $N=84$ and 86 and produces a distinct $K^\pi=4^+$ band in strongly deformed nuclei [2312.01593]. For $^{154}$Gd, the calculated band built on the $4_3^+$ state is identified as $K^\pi=4^+$, and the predicted
$$
B(E4;4^+_{K=4^+}\to 0^+_1)\approx 93\ \text{W.u.}
$$
in $sdg$-IBM contrasts with approximately $1.3$ W.u. in $sd$-IBM, while both are constrained to reproduce
$$
B(E4;4^+_1\to0^+_1)=38\pm3\ \text{W.u.}
$$
in $^{154}$Gd [2312.01593].

Near $N=90$, the mapped $sdg$-IBM finds that $g$ bosons improve the description of $J^\pi\ge6^+$ yrast energies in nuclei with $N=84$ and 86 and increase quadrupole transition strengths between yrast states in the well-deformed $N=90$ and 92 nuclei [2403.03393]. The same work reports reduced $E4$ matrix elements for Gd isotopes, including for $^{160}$Gd an $sdg$-IBM value of $0.35$ e·b$^2$, matching the experimental $0.35^{+0.09}_{-0.07}$ e·b$^2$ [2403.03393].

Rotational observables in the A$\approx180$ region show that negative $\beta_4$ lowers the high-frequency moments of inertia while leaving low-spin moments of inertia largely unaffected. The HFBC and rigid-body calculations display similar trends: at normal deformation, axial $\beta_4<0$ reduces the moment of inertia, whereas $\beta_4>0$ increases it [2601.10052]. The single-particle interpretation given there emphasizes enhanced mixing of $\Delta l=\Delta j=4$ partners and stronger shell gaps near the Fermi surface when $\beta_4\neq0$ [2601.10052].

Reaction observables are often especially sensitive. In $^{28}$Si+$^{92}$Zr fusion, the barrier distribution
$$
D_{\mathrm{fus}}(E)=\frac{d^2[E\sigma_{\mathrm{fus}}(E)]}{dE^2}
$$
changes markedly when $\beta_4$ is switched from $+0.25$ to $-0.25$, despite identical $|\beta_4|$ [1805.03395]. In neutron-rich Zr isotopes, deformation-induced changes in radii and surface diffuseness generate an approximately $100$ mb enhancement of total reaction cross sections relative to spherical-constrained calculations across $N=60$–$74$ [2304.06238].

A common misconception is that $\beta_4$ is only a minor correction to $\beta_2$. The cited studies do not support that simplification. In some nuclei it is small, as in $^{28}$Si; in others it changes the barrier distribution qualitatively, generates $K^\pi=4^+$ collectivity, modifies high-spin rotational spacings, or contributes correlation energy comparable to the quadrupole correlation energy itself [2303.12495][2312.01593][2508.05268].

## 6. High-energy collisions, machine-learning identifiability, and unresolved issues

Static hexadecapole deformation has also entered relativistic heavy-ion phenomenology. For $^{238}$U, hydrodynamic studies argue that previous implementations conflated surface and volume deformations. Skyrme-HFB fits to microscopic densities give, for a representative BSkG2 parametrization,
$$
\beta_{20}=+0.280,\qquad \beta_{40}=+0.153,\qquad
\beta^{WS}_{20}=+0.247,\qquad \beta^{WS}_{40}=+0.081,
$$
so the realistic Woods–Saxon surface quadrupole is significantly smaller than the volume quadrupole because of the nonzero surface hexadecapole [2302.13617]. Correcting this mapping restores agreement between IP-Glasma+MUSIC+UrQMD simulations and RHIC data for central U+U collisions [2302.13617].

A more targeted proposal uses the nonlinear response coefficient
$$
\chi_{4,22}\equiv \frac{v_4\{\Psi_2\}}{\langle v_2^4\rangle^{1/2}}
= \frac{\langle\!\langle 3\rangle\!\rangle_{2,2,-4}}{\langle v_2^4\rangle}
$$
in ultra-central U+U versus Au+Au collisions. The relative difference
$$
R(X)=\frac{2(X_{UU}-X_{AuAu})}{X_{UU}+X_{AuAu}}
$$
is reported to be nearly zero and flat in centrality when $\beta_{4,U}=0$, insensitive to $\beta_2$, and clearly nonzero when $\beta_{4,U}>0$ [2402.16550]. This suggests a route to constraining $\beta_4$ of $^{238}$U that is complementary to low-energy electromagnetic data, where the $\beta_4$ effect is overwhelmed by large $\beta_2$ [2402.16550].

Machine-learning analyses of heavy-ion initial conditions reach a related conclusion. For deformed Woods–Saxon configurations of $^{238}$U sampled on a $21\times21$ grid in
$\beta_2\in[-0.5,0.5]$ and $\beta_4\in[-0.2,0.2]$, permutation-invariant point-cloud networks recover $\beta_4$ with test $R^2=0.5825$ for a single configuration and $0.9600$ for $N=20$ aggregated configurations [2603.24088]. In TRENTo entropy-density images, $\beta_4$ is much less identifiable in single events: at 0–10% centrality the regression test scores are $0.0871$ for $M=1$, $0.3682$ for $M=10$, $0.5604$ for $M=50$, and $0.6197$ for $M=100$, while SBI posterior means give $0.1034$, $0.4048$, $0.6075$, and $0.6566$ for the same bag sizes [2603.24088]. The work states that multi-event averaging is essential and that $\beta_4$ remains intrinsically harder than $\beta_2$ to extract [2603.24088].

Several unresolved issues recur across the literature. One is sign ambiguity: in $^{74,76}$Kr the reaction cross sections alone admit positive and negative $\beta_4$ solutions, requiring EDF input to choose between them [2304.14246]. Another is model dependence: quasi-elastic extractions depend on coupled-channels truncations and assumptions such as $\beta_N=\beta_C$ [1811.12756][2303.12495], QRPA softness maps do not determine the sign of $\beta_4$ [2401.06117], and beyond-mean-field actinide calculations remain restricted to axial symmetry without octupole coupling [2502.04985]. A further point is representational ambiguity: surface $\beta^{WS}_4$, volume moments, transition-derived deformations, and intrinsic mean-field $\beta_4$ are related but not interchangeable [2302.13617].

Taken together, these studies define static hexadecapole deformation as a measurable and theoretically consequential component of nuclear structure rather than a peripheral correction. Its manifestation spans precision low-energy reaction analyses, collective spectroscopy, mean-field topology, rotational dynamics, and relativistic heavy-ion observables, while its quantitative determination remains contingent on the chosen deformation convention, the degree of collective-mode coupling retained, and the probe used to isolate the $\lambda=4$ degree of freedom.

Source: https://www.emergentmind.com/topics/static-hexadecapole-deformation