---
title: Quadrupole–Hexadecapole Configuration Mixing
url: https://www.emergentmind.com/topics/quadrupole-hexadecapole-configuration-mixing
type: topic
---

# Quadrupole–Hexadecapole Configuration Mixing

Quadrupole–hexadecapole configuration mixing denotes the nonseparable coupling between rank-2 and rank-4 multipole degrees of freedom, most commonly the axial mass moments $Q_{20}$ and $Q_{40}$ or the corresponding deformation parameters $\beta_2$ and $\beta_4$. In nuclear structure it appears as a tilted collective valley in the $(\beta_2,\beta_4)$ plane, as off-diagonal couplings in a two-dimensional Generator Coordinate Method (GCM), or as admixture between $d$- and $g$-boson configurations in $sdg$ interacting-boson models. In deformed-density parameterizations it also appears as a mismatch between volume multipoles and Woods–Saxon surface parameters, so that a surface hexadecapole term contributes to the volume quadrupole moment. In correlated electron systems, an analogous phenomenon occurs when magnetic quadrupole and hexadecapole order parameters belong to the same irreducible representation and therefore admix once symmetry is lowered. Across these settings, the central theme is that quadrupole and hexadecapole modes are often interwoven rather than independent, with consequences for correlation energies, excitation spectra, transition strengths, rotational response, flow observables, and magnetoelectric transport [2308.03650] [2302.13617] [1705.10349] [2502.04985].

## 1. Operator content and symmetry conditions

The standard intrinsic mass multipole operators are
\[
\hat Q_{20}=\sum_{i=1}^A r_i^2\,Y_{20}(\hat r_i),\qquad
\hat Q_{40}=\sum_{i=1}^A r_i^4\,Y_{40}(\hat r_i),
\]
with associated deformation coordinates
\[
\beta_\ell=\frac{\sqrt{4\pi(2\ell+1)}}{3\,R_0^\ell\,A}\,
\langle\varphi|\hat Q_{\ell 0}|\varphi\rangle,\qquad
R_0=1.2\,A^{1/3}\,\mathrm{fm},\qquad \ell=2,4.
\]
Equivalent definitions are used in Gogny-HFB, GCM, and boson-mapping studies, sometimes with $Q_{20}=z^2-\tfrac12(x^2+y^2)$ and $Q_{40}=z^4-3r_\perp^2 z^2-\tfrac38 r_\perp^4$ in Cartesian form [2308.03650] [2508.05268].

A recurrent microscopic signature of mixing is that the minimum-energy path in $(\beta_2,\beta_4)$ space is tilted rather than aligned with either axis. In $^{154}$Sm, the bottom of the prolate valley follows roughly $\beta_4\approx1.8\,\beta_2-0.4$, so quadrupole and hexadecapole deformations vary in lockstep; in actinides, the principal valley direction in $^{238}$U is reported as $\beta_4=a_A\beta_2+b_A$ with $a_A\simeq2.50$ and $b_A\simeq-0.52$ [2308.03650] [2502.04985]. Such relations are the geometric expression of a nonzero mixed curvature in the collective potential.

In magnetic systems the same logic is formulated group-theoretically. For Ba$_{1-x}$K$_x$Mn$_2$As$_2$, the crystal space group is $I4/mmm$ with point group $D_{4h}$, while G-type antiferromagnetism breaks $\mathcal T$ and global inversion $\mathcal P$ but preserves $\mathcal{PT}$, reducing the point symmetry to $D_{2d}$. Under $D_{4h}\to D_{2d}$, only the $B_{1u}$ irrep reduces to the totally symmetric $A_1$ of $D_{2d}$, so the antiferromagnetic order parameter must transform as $B_{1u}$. The lowest-rank magnetic multipoles in this irrep are the magnetic quadrupole $M_{22}^+$ and the magnetic hexadecapole $M_{42}^+$, which therefore can admix on symmetry grounds [1705.10349].

## 2. Collective and microscopic formulations of the mixing problem

The most explicit formulation in nuclear structure is the two-dimensional GCM built from Hartree–Fock–Bogoliubov intrinsic vacua constrained simultaneously in $\beta_2$ and $\beta_4$:
\[
|\Psi^J_\sigma\rangle=
\int d\beta_2\,d\beta_4\;
f^J_\sigma(\beta_2,\beta_4)\,
\hat P^J_{00}\,
|\Phi(\beta_2,\beta_4)\rangle ,
\]
with weights determined by the Hill–Wheeler–Griffin equation
\[
\int d\beta'_2\,d\beta'_4\,
\Bigl[\mathcal H^J(\beta,\beta')-E^J_\sigma\mathcal N^J(\beta,\beta')\Bigr]
f^J_\sigma(\beta')=0.
\]
The normalized amplitudes $G^J_\sigma(\beta_2,\beta_4)$ then provide a probability distribution in the $(\beta_2,\beta_4)$ plane [2308.03650]. Closely related Gogny-EDF studies of rare-earth and actinide nuclei use the same structure without full angular-momentum projection, again interpreting $|G^\sigma(\beta_2,\beta_4)|^2$ as the collective probability density [2508.05268] [2502.04985].

Near a minimum, the microscopic origin of mixing is often summarized by a quadratic expansion of the collective potential,
\[
V(\beta_2,\beta_4)\simeq
V_0+\tfrac12 C_{22}q_2^2 + C_{24}q_2q_4 + \tfrac12 C_{44}q_4^2,
\]
where $q_\lambda=\beta_\lambda-\beta_\lambda^0$. The cross-term $C_{24}$ is the static quadrupole–hexadecapole coupling. In the adiabatic description of actinides there is also an off-diagonal inertia $B^{24}$, so coupling enters both the potential and kinetic sectors of the collective Hamiltonian [2502.04985].

Not every study constructs an explicit off-diagonal collective kernel. In the macroscopic–microscopic Woods–Saxon plus HFBC treatment of the $A\approx180$ region, the coupling is introduced directly through the simultaneous $\beta_2$ and $\beta_4$ dependence of the deformed surface in the mean-field Hamiltonian. The work explicitly states that it does not build a GCM-type mixing Hamiltonian in collective coordinates, but instead explores the coupling through total-energy-surface scans and the induced single-particle admixtures at a fixed equilibrium shape [2601.10052].

In odd-parity magnetic metals, a Landau description is natural. Writing $Q\equiv\langle M_{22}^+\rangle$ and $H\equiv\langle M_{42}^+\rangle$, the free energy is
\[
F[Q,H]=\tfrac12\alpha_QQ^2+\tfrac12\alpha_HH^2+\beta_{QH}QH+\cdots.
\]
In the insulating limit, $H$ condenses at $T_N$ and the order is pure hexadecapolar; in the metallic state, LS coupling generates $\beta_{QH}\neq0$, so a nonzero $H$ generically induces $Q\propto(-\beta_{QH}/\alpha_Q)H$ [1705.10349]. A plausible implication is that “configuration mixing” across subfields refers less to a unique formalism than to a common structural feature: symmetry-allowed off-diagonal coupling between rank-2 and rank-4 sectors.

## 3. Beyond-mean-field nuclear structure: correlation energies, wave functions, and shape coexistence

The axial HFB+GCM study of even-even Sm and Gd isotopes finds strong coupling between the quadrupole and hexadecapole degrees of freedom. For nuclei around $N\approx90$, the static HFB surface shows a well-developed prolate valley with positive $\beta_4$, a much higher-lying oblate pocket with $\beta_4<0$, and near-spherical minima at the semi-magic nucleus $^{190}$Gd. In $^{154}$Sm, the ground-state collective wave function is a tilted two-dimensional Gaussian elongated along the valley direction, the first excited $0^+$ is a one-phonon excitation along that same mixed $\beta_2$–$\beta_4$ direction, and the second $0^+$ is predominantly one-phonon along the perpendicular direction [2308.03650].

The same study reports that the gain from the $\beta_2$–$\beta_4$ space is $E_c\approx1.1\text{–}1.4$ MeV in Sm/Gd, about twice the $\beta_2$-only result, implying a non-negligible hexadecapole contribution of about $0.6$ MeV. The first and second $0^+$ excitations appear at $E\approx2\text{–}4$ MeV in most cases, while in isotopes such as $^{180}$Gd and lighter Sm/Gd near $N=88$ the oblate pocket comes down, so the first excited $0^+$ often localizes in the oblate minimum at $E_1<1$ MeV with $\beta_2<0$ and $\beta_4\simeq0$ [2308.03650]. This distinguishes genuine mixed-mode vibrational states from low-lying shape-coexisting configurations.

For Yb, Hf, W, and Os isotopes in the interval $170\le A\le202$, a 2D-GCM analysis finds that ground and excited states of the lighter isotopes are associated with diamond-like shapes, while a region of square-like shapes occurs below the neutron shell closure $N=126$. The quadrupole and hexadecapole degrees of freedom are reported to be interwoven up to about $A=184$–$188$, after which the collective wave functions align predominantly along the $\beta_2$ axis and the motion becomes effectively quadrupole-dominated. The total correlation energy $E_{\rm corr}^{2D}$ ranges from about $1.15$ MeV up to about $1.73$ MeV, with an additional hexadecapole gain $\delta E=E_{\rm corr}^{2D}-E_{\rm corr}^{\beta_2}\approx0.66\text{–}0.76$ MeV [2508.05268].

Actinide calculations based on Gogny HFB and GCM reach a parallel conclusion. For selected Ra, Th, U, and Pu isotopes, static hexadecapole deformations are sizable near $^{238}$U, and a region with small negative hexadecapole deformation just below $N=184$ remains stable once zero-point quadrupole–hexadecapole fluctuations are included. The quoted correlation energies are $E_{\rm corr,2D}\simeq1.42\text{–}1.72$ MeV and $E_{\rm corr,1D}^{Q2}\simeq0.45\text{–}0.98$ MeV, so the extra gain from including $\beta_4$ is $\Delta E\simeq0.74\text{–}0.97$ MeV. For $^{238}$U specifically, the reported values are $E_{\rm corr,1D}^{Q2}\simeq0.90$ MeV, $E_{\rm corr,2D}\simeq1.62$ MeV, and therefore $\Delta E\simeq0.72$ MeV [2502.04985]. The data are consistent with the repeated observation that hexadecapole correlations contribute an amount of binding comparable to the quadrupole correlation energy itself.

A more mean-field-oriented treatment of $^{180-184}$Yb, $^{182-186}$Hf, and neighboring nuclei also emphasizes the structural role of $\beta_4$, even though it does not implement explicit GCM mixing. In $^{184}$Hf the global minimum is found at $(\beta_2\approx0.237,\gamma\approx0,\beta_4\approx-0.087)$, whereas setting $\beta_4=0$ shifts the quadrupole minimum and raises the total energy by $1$–$3$ MeV. At finite $\beta_4$ a pronounced shell gap at $Z=72$, $N=112$ opens, and the proton $h_{11/2}$ orbital with $\Omega=7/2$ shows increased $\Delta\ell=\Delta j=4$ admixture, with a $j_{15/2}$ component of about $4.7\%$ instead of about $2.7\%$ at $\beta_4=0$ [2601.10052]. This suggests that even without collective superposition, quadrupole–hexadecapole coupling can be read directly from deformation-induced orbital mixing.

## 4. Boson realizations and spectroscopic fingerprints

In interacting-boson descriptions, quadrupole–hexadecapole mixing is represented by the coexistence of $s$ ($0^+$), $d$ ($2^+$), and $g$ ($4^+$) bosons. A widely used Hamiltonian is
\[
\hat H_{sdg}=
\epsilon_d\,\hat n_d+\epsilon_g\,\hat n_g
+\kappa\,\hat Q^{(2)}\!\cdot\!\hat Q^{(2)}
+\kappa(1-\chi^2)\,\hat Q^{(4)}\!\cdot\!\hat Q^{(4)},
\]
with
\[
\hat Q^{(4)}=(s^\dagger\times\tilde g+g^\dagger\times\tilde s)^{(4)},
\]
and a quadrupole operator containing the mixed $(d^\dagger\times\tilde g+g^\dagger\times\tilde d)^{(2)}$ terms [2411.11331] [2403.03393]. These are the explicit off-diagonal pieces that connect configurations differing by one $d$ and one $g$ boson and thereby mix quadrupole and hexadecapole structures.

In neutron-rich Sm and Gd isotopes, Gogny-D1S mapped $sdg$-IBM calculations show that the largest differences between $sd$-IBM and $sdg$-IBM occur in transitional nuclei around $N\approx90$. In Sm with $N\le92$, ground-state-band states carry $\langle \hat n_g\rangle\ge1$, while for $N\ge94$ they have $\langle \hat n_g\rangle\lesssim0.5$. The excited $0_2^+$ and $2_2^+$ states have $30$–$50\%$ $g$-boson fractions for $N\approx90$, but only $5$–$10\%$ for $N\ge96$. The mapped $sdg$-IBM improves the description of high-spin yrast states in lighter Sm and Gd and enhances non-yrast $B(E4)$ strengths, whereas for $N\ge94$ the main residual sensitivity is in the $E0$ monopole systematics [2411.11331].

A related mapped $sdg$-IBM study of Nd, Sm, Gd, Dy, and Er near $N=90$ concludes that the inclusion of the $g$ boson is necessary to improve the $J^\pi\ge6^+$ yrast energies in nuclei with $N=84$ and $86$, near the neutron shell closure. In the well-deformed $N=90$ and $92$ systems, the model increases the quadrupole transition strengths between yrast states, in better agreement with experiment, while the monopole strengths do not differ significantly from those of the simpler $sd$ model [2403.03393].

For axially deformed Gd isotopes, the explicit inclusion of the hexadecapole degree of freedom is reported not to affect most low-spin and low-lying states qualitatively, but to improve the description of high-spin states of ground-state bands in nearly spherical vibrational nuclei and to generate $K^\pi=4^+$ bands with strong $E4$ transitions in strongly deformed nuclei. Illustrative values include downward shifts of roughly $100$ keV, $200$ keV, and $350$ keV for the $6_1^+$, $8_1^+$, and $10_1^+$ levels in $^{148}$Gd, and for $^{154}$Gd a $K^\pi=4^+$ bandhead at about $1.65$ MeV, close to the experimental $1.646$ MeV, together with a predicted $B(E4)\approx93$ W.u. as opposed to $1.3$ W.u. in the $sd$-IBM [2312.01593].

Mixed-symmetry hexadecapole states provide a particularly transparent example. In $^{96}$Ru, the $sdg$-IBM-2 yields for the $4_2^+$ state an approximate composition of $43\%$ $s$ bosons, $34\%$ $d$ bosons, and $23\%$ $g$ bosons. The observed strong $M1$ transition $B(M1;4_2^+\to4_1^+)=0.90(18)\,\mu_N^2$ is analyzed as $83\%$ arising from the $g$-boson term and $17\%$ from the $d$-boson term, while the moderate branch $B(E2;4_2^+\to2_1^+)=1.52(19)$ W.u. is tied to the $d$-boson admixture. The model also predicts $B(E4;4_2^+\to0_1^+)=0.55$ W.u. [1502.06409]. In this setting, configuration mixing is directly observable in the decomposition of the wave function and in the relative strength of $M1$, $E2$, and $E4$ channels.

## 5. Surface–volume mixing and the heavy-ion-collision realization

A distinct but conceptually related form of quadrupole–hexadecapole mixing arises when one compares volume multipole moments with surface deformation parameters in Woods–Saxon densities. The volume deformations
\[
\beta_{\ell m}=\frac{4\pi}{3R_0^\ell A}\int d^3r\,\rho(r)\,r^\ell\,Y_{\ell m}(\theta,\phi)
\]
characterize the density throughout the nucleus, whereas the Woods–Saxon surface parameters $\beta_{\ell m}^{WS}$ enter the angle-dependent radius
\[
R(\theta,\phi)=R_d\Bigl(1+\sum_{\ell,m}\beta_{\ell m}^{WS}Y_{\ell m}(\theta,\phi)\Bigr).
\]
These two notions are not identical for a nucleus that carries both quadrupole and hexadecapole deformation [2302.13617].

For an axially symmetric nucleus with only $\beta_{20}^{WS}$ and $\beta_{40}^{WS}$ nonzero, expanding the quadrupole volume moment in the sharp-surface limit gives
\[
\beta_2 \simeq \beta_2^{WS}
+\frac{2}{7}\sqrt{\frac{5}{\pi}}(\beta_2^{WS})^2
+\frac{12}{7\sqrt{\pi}}\beta_2^{WS}\beta_4^{WS}
+\cdots.
\]
The last term is the explicit quadrupole–hexadecapole mixing term: a surface hexadecapole distortion contributes to the volume quadrupole moment [2302.13617]. One common misconception is therefore that low-energy tabulated $\beta_2$ values can be inserted directly as Woods–Saxon surface parameters in high-energy collision simulations.

The case of $^{238}$U is the canonical demonstration. Using 21 Skyrme parameter sets in 3D Skyrme-HFB, the quoted ensemble averages are $\langle\beta_2\rangle=0.280\pm0.005$ and $\langle\beta_4\rangle=0.153\pm0.014$ for the volume multipoles, with corresponding best-fit Woods–Saxon parameters $\beta_2^{WS}\simeq0.247\pm0.005$ and $\beta_4^{WS}\simeq0.081\pm0.012$ [2302.13617]. In other words, reproducing a true volume $\beta_2\approx0.28$ requires a smaller surface quadrupole once the hexadecapole deformation is included.

This correction alters the initial conditions used in hydrodynamic simulations of U+U collisions. The older practice of taking $\beta_2^{WS}(^{238}\mathrm U)=0.28$ and neglecting $\beta_4^{WS}$ led to an overestimate of elliptic flow. Using the ratio
\[
r_{Au,U}^2\equiv
\langle v_2^2\rangle_{U+U}/\langle v_2^2\rangle_{Au+Au},
\]
the reported estimate is $r^2\approx1.78$ for the old choice, compared with the STAR measurement $1.49\pm0.05$. Replacing it by $\beta_2^{WS}(^{238}\mathrm U)\approx0.25$ and $\beta_4^{WS}\approx0.08$ yields $r^2\approx1.55\pm0.10$, and full IP-Glasma + MUSIC + UrQMD calculations recover the observed reduction of $v_2$ in central U+U collisions [2302.13617]. The same work emphasizes that $\beta_4^{WS}$ has virtually no direct effect on $v_2$ fluctuations in central events; its importance is indirect, through the redefinition of the effective surface quadrupole.

## 6. Magnetic quadrupole–hexadecapole mixing in odd-parity antiferromagnets

In BaMn$_2$As$_2$ and Ba$_{1-x}$K$_x$Mn$_2$As$_2$, the relevant order is not a mass deformation but an odd-parity magnetic multipole. The magnetic quadrupole and hexadecapole operators in the $B_{1u}$ representation are given in cubic-harmonic form as
\[
M_{22}^+=\mu_B\sqrt3\,(x\,s_x-y\,s_y),
\]
and
\[
M_{42}^+=\mu_B\Bigl[
3\sqrt5\,z(x^2-y^2)s_z
+\frac{\sqrt5}{2}(7z^2-r^2)(x\,s_x-y\,s_y)
-\frac{\sqrt5}{2}(x^2-y^2)(x\,s_x+y\,s_y)
\Bigr].
\]
The ground state of semiconducting BaMn$_2$As$_2$ is identified as a magnetic hexadecapole ordered state, and microscopic analysis shows ferroic ordering of a leading magnetic hexadecapole moment together with an admixed magnetic quadrupole moment [1705.10349].

In the metallic regime the admixture becomes explicit. LS coupling induces a nonzero $\langle M_{22}^+\rangle$ while the leading order remains $\langle M_{42}^+\rangle$, and the effective single-band Hamiltonian is built on two Mn sublattices with odd-parity antisymmetric spin–orbit coupling and opposite antiferromagnetic molecular fields on the two sublattices [1705.10349]. In Landau language, once the hexadecapole condenses, the linear coupling term $\beta_{QH}QH$ generically induces the quadrupole component as well.

The mixed order produces characteristic electromagnetic and transport responses. By symmetry, the uniform magnetoelectric tensor has only $\alpha_{xx}=-\alpha_{yy}\neq0$ in the hexadecapole state; in the insulating case the reported magnitude is $\alpha\sim10^{-4}$, enhanced in the metal. The antiferromagnetic Edelstein response is expressed as
\[
M_\mu^{AF}=\alpha_{\mu\nu}^{AF}E_\nu,\qquad
\alpha_{xy}^{AF}=\alpha_{yx}^{AF}\neq0,
\]
with an intraband contribution proportional to the relaxation time $\tau$. A further consequence is current-induced nematic order: tetrahedral Fermi-surface warping proportional to $k_xk_yk_z$ implies that a current along the $z$ axis induces $xy$-plane nematicity, a $J_z\to\phi_{xy}$ coupling of $B_{2g}$ symmetry, described as magnetopiezoelectricity [1705.10349].

The proposed experimental signatures are correspondingly diverse: neutron scattering should detect the $q=0$ spin-density anisotropy associated with the primary hexadecapole form factor, with the small quadrupole admixture shifting intensity between $(h,k,0)$ reflections; an out-of-plane current should drive an in-plane shear strain observable by X-ray or ultrasonic probes; $\alpha_{xx}=-\alpha_{yy}$ may be accessed by polarized optics or dielectric tuning under magnetic field; and a basal-plane current should generate a staggered sublattice magnetization that may be detectable by XMCD or second-harmonic magneto-optics [1705.10349]. In this condensed-matter realization, quadrupole–hexadecapole mixing is not a secondary analogy to nuclear collectivity but a symmetry-equivalent instance of the same general principle: rank-2 and rank-4 multipoles of the same symmetry class do not remain independent once the microscopic Hamiltonian permits their coupling.

Source: https://www.emergentmind.com/topics/quadrupole-hexadecapole-configuration-mixing