---
title: Dipolar-Octupolar Spin Hamiltonian Overview
url: https://www.emergentmind.com/topics/dipolar-octupolar-spin-hamiltonian
type: topic
---

# Dipolar-Octupolar Spin Hamiltonian Overview

The dipolar–octupolar spin Hamiltonian is the effective nearest-neighbor pseudospin description for rare-earth pyrochlores whose crystal-field ground state is a dipolar–octupolar Kramers doublet. Its defining feature is that the three pseudospin components do not transform uniformly under the local \(D_{3d}\) symmetry: two components transform as magnetic dipoles, while the third transforms as a magnetic octupole. On the pyrochlore lattice of corner-sharing tetrahedra, this symmetry structure permits an unusually economical exchange Hamiltonian, yet supports a wide range of phases, including all-in–all-out order, moment fragmentation, and several \(U(1)\) quantum spin ice regimes with either \(0\)-flux or \(\pi\)-flux backgrounds [2201.00828], [2002.05608], [2312.03106].

## 1. Local doublets, local frames, and multipolar content

In Ce-, Nd-, and Sm-based pyrochlores, the low-energy magnetic degree of freedom is a pseudospin-\(1/2\) defined within a crystal-field Kramers doublet on each rare-earth site. The natural local quantization frame is tied to the four sublattices of the pyrochlore lattice, with local Ising axes along the four \([111]\)-type directions, for example
\[
\hat z_0=\frac{(1,1,1)}{\sqrt3},\quad
\hat z_1=\frac{(1,-1,-1)}{\sqrt3},\quad
\hat z_2=\frac{(-1,1,-1)}{\sqrt3},\quad
\hat z_3=\frac{(-1,-1,1)}{\sqrt3}.
\]
These axes define the site-dependent local frames in which the effective Hamiltonian is written [2201.00828], [2601.03202].

A central property of the dipolar–octupolar doublet is the symmetry distinction among pseudospin components. In one standard convention, used in several Ce-based analyses, \(S^x\) and \(S^z\) transform as magnetic dipoles, while \(S^y\) transforms as a rank-3 magnetic octupole; all three are odd under time reversal [2201.00828], [2312.03106]. In other conventions, especially in some Ce\(_2\)Hf\(_2\)O\(_7\) and Ce\(_2\)Sn\(_2\)O\(_7\) analyses, \(s^z\) is taken as the dipolar component and \(s^{x,y}\) as octupolar components [2305.08261], [2211.15140]. This reflects the use of different rotated local bases rather than different underlying symmetry content.

The magnetic-field and neutron-selection rules follow directly from this multipolar structure. Dipolar components couple linearly to a uniform magnetic field, while the octupolar component does not couple linearly and is typically inefficiently probed by standard magnetic neutron scattering. In Ce\(_2\)Zr\(_2\)O\(_7\), the effective Zeeman term is dominated by the local \(\hat z_i\) projection, and \(g_x\) is strongly suppressed, so weak fields and neutrons couple predominantly to \(S^z\) [2201.00828]. In Ce\(_2\)Hf\(_2\)O\(_7\) and hydrothermal Ce\(_2\)Sn\(_2\)O\(_7\), the working approximation is likewise that only the dipolar \(s^z\) or \(\tau^z\) component couples linearly to the field [2305.08261], [2211.15140].

## 2. Symmetry-allowed Hamiltonians and rotated-basis forms

The most general nearest-neighbor bilinear Hamiltonian for a dipolar–octupolar doublet on the pyrochlore lattice contains three diagonal exchanges and a single symmetry-allowed off-diagonal \(x\)–\(z\) term,
\[
H_l=\sum_{\langle i,j\rangle}
\Big[
J_{xx} S_i^x S_j^x
+J_{yy} S_i^y S_j^y
+J_{zz} S_i^z S_j^z
+J_{xz}(S_i^x S_j^z+S_i^z S_j^x)
\Big].
\]
The absence of \(J_{xy}\) and \(J_{yz}\) terms follows from the fact that the octupolar component transforms in a different irreducible representation from the two dipolar components [2312.03106], [2311.04269].

A uniform rotation about the local \(y\)-axis diagonalizes this Hamiltonian. Defining
\[
\tau^x=S^x\cos\theta-S^z\sin\theta,\qquad
\tau^y=S^y,\qquad
\tau^z=S^x\sin\theta+S^z\cos\theta,
\]
with
\[
\tan(2\theta)=\frac{2J_{xz}}{J_{zz}-J_{xx}},
\]
one obtains the diagonal XYZ form
\[
H_l=\sum_{\langle i,j\rangle}\sum_{a\in\{x,y,z\}}\mathcal J_a\,\tau_i^a\tau_j^a,
\]
where
\[
\mathcal J_y=J_{yy},
\]
\[
\mathcal J_x=
\frac{J_{zz}+J_{xx}}{2}
-\frac{1}{2}\sqrt{(J_{zz}-J_{xx})^2+4J_{xz}^2},
\]
\[
\mathcal J_z=
\frac{J_{zz}+J_{xx}}{2}
+\frac{1}{2}\sqrt{(J_{zz}-J_{xx})^2+4J_{xz}^2}.
\]
In this rotated basis, the nearest-neighbor exchanges are bond-independent, and the bond-dependent complex phase factors familiar from conventional anisotropic pyrochlore Hamiltonians are not needed for the DO doublet [2312.03106], [2002.05608].

Several papers further specialize this Hamiltonian. For Ce\(_2\)Zr\(_2\)O\(_7\), analyses often set \(J_{xz}\approx 0\), leading to a minimal nearest-neighbor XYZ model
\[
H_{\mathrm{XYZ}}=\sum_{\langle i,j\rangle}
\big(
\mathcal J_x S_i^xS_j^x+\mathcal J_y S_i^yS_j^y+\mathcal J_z S_i^zS_j^z
\big),
\]
which is the form used in exact diagonalization, molecular dynamics, and pseudofermion FRG studies [2201.00828], [2311.04269]. In octupolar-spin-ice treatments, an equivalent DO–XXZ form is often written with a dominant \(J_{yy}\) term and transverse couplings \(J_\pm\) and \(J_{\pm\pm}\) [2406.01472], [2301.05240].

## 3. Quantum spin ice regimes and emergent gauge structure

Near an Ising limit, the dipolar–octupolar Hamiltonian supports several \(U(1)\) quantum spin ice regimes. The basic distinction is whether the dominant Ising component lies in a dipolar or octupolar channel. In the notation of the rotated XYZ model, octupolar QSI arises when \(\mathcal J_y>0\) is the dominant antiferromagnetic coupling, while dipolar QSI arises when \(\mathcal J_z>0\) or \(\mathcal J_x>0\) is dominant [2201.00828].

The low-energy description is an emergent lattice \(U(1)\) gauge theory,
\[
H_{\rm eff}=U\,\mathbf E^2+K\cos(\nabla\times\mathbf A),
\]
where the Ising component maps to an emergent electric field and the transverse components generate gauge dynamics. The sign of the effective ring-exchange term \(K\) determines the background flux sector: unfrustrated transverse interactions favor \(0\)-flux, while frustrated transverse interactions favor \(\pi\)-flux [2201.00828]. Accordingly, four QSI regimes appear in the nearest-neighbor DO model: \(0\)-O-QSI, \(\pi\)-O-QSI, \(0\)-D-QSI, and \(\pi\)-D-QSI [2201.00828], [2002.05608].

In octupolar formulations with dominant \(J_{yy}\), the same distinction is հաճախ encoded by the sign of \(J_\pm\): within the GMFT used for DO pyrochlores, \(J_\pm/J_{yy}>0\) corresponds to \(0\)-flux QSI and \(J_\pm/J_{yy}<0\) to \(\pi\)-flux QSI [2406.01472], [2406.18650]. In Ce\(_2\)Hf\(_2\)O\(_7\), the fitted parameter sets yield \(J_\pm<0\), which places the system in the \(\pi\)-flux octupolar QSI sector; the ring-exchange scale \(J_{\rm ring}\) is then of order a few \(10^{-3}\) meV [2305.08261].

The emergent excitations depend on this gauge structure. The standard \(U(1)\) description contains gapless emergent photons and gapped spinons. In DO systems, however, the neutron cross section is dominated by the dipolar channel, so the octupolar Ising variable can be nearly invisible even when it is the field that enforces the ice rule. This asymmetry between the internal gauge variable and the observable dipolar channel is one of the characteristic features of the dipolar–octupolar problem [2201.00828], [2301.05240].

## 4. Neutron selection rules and structure-factor diagnostics

Because neutrons couple to magnetic dipoles, the measured response is controlled by the dipolar sector of the DO doublet. The basic dynamical correlator is
\[
S^{\alpha\beta}(\mathbf q,\omega)=
\sum_{i,j}e^{i\mathbf q\cdot(\mathbf r_i-\mathbf r_j)}
\int dt\,e^{i\omega t}
\langle S_i^\alpha(0)S_j^\beta(t)\rangle.
\]
For Ce\(_2\)Zr\(_2\)O\(_7\), where \(g_x\approx 0\), the equal-time neutron structure factor is dominated by \(S^z\) correlations,
\[
S(\mathbf q)=\frac{1}{N}\sum_{i,j}
\left[
\hat z_i\cdot\hat z_j
-\frac{(\hat z_i\cdot\mathbf q)(\hat z_j\cdot\mathbf q)}{q^2}
\right]
e^{-i\mathbf q\cdot(\mathbf R_i-\mathbf R_j)}
\langle S_i^zS_j^z\rangle.
\]
This polarization factor is the standard local-frame projection for magnetic neutron scattering [2201.00828].

Classical and quantum calculations produce sharply different signatures. In classical spin ice or classical molecular dynamics, the Ising-sector correlations show sharp pinch points characteristic of divergence-free constraints. In exact diagonalization, these pinch points are substantially smeared by quantum fluctuations of the \(U(1)\) gauge field, even in unfrustrated regimes [2201.00828]. In the \(\pi\)-flux octupolar QSI regime, the equal-time structure factor displays high-intensity peaks at \([001]\) and \([003]\) together with rod-like motifs along \([00l]\); these peaks are absent in classical MD and in the other QSI regimes studied in the same work [2201.00828].

The dynamical structure factor provides an even sharper distinction. Exact diagonalization for \(\pi\)-O-QSI finds a pronounced high-intensity feature at the \(X=[001]\) point and a gapped spectrum with \(\Delta/\mathcal J_y\approx 0.6\) for a representative parameter set, suggestive of gapped spinon excitations; by contrast, \(0\)-O-QSI and the dipolar-QSI regimes show suppressed intensity at \(X\) and much weaker momentum dependence [2201.00828]. A later GMFT analysis of octupolar QSI predicted that \(\pi\)-O-QSI should exhibit a broad continuum with three distinctive peaks in the momentum-integrated spectrum, arising from two mostly flat spinon bands [2301.05240].

Polarized neutron channels also encode the flux sector. PFFRG calculations for the nearest-neighbor DO XYZ model showed that the \(\pi\)-flux QSI regime produces pronounced rod-like features in the spin-flip channel and intensity minima at the centers of the fcc-Brillouin-zone hexagons in the non-spin-flip channel, consistent with polarized-neutron observations on Ce\(_2\)Zr\(_2\)O\(_7\) [2311.04269]. This point is important because it rules out the common simplification that all rod-like diffuse scattering in a DO pyrochlore is equivalent: the detailed SF/NSF modulation depends on both the flux background and the role of quantum fluctuations.

## 5. Material realizations and competing interpretations

Ce\(_2\)Zr\(_2\)O\(_7\) is the most extensively analyzed DO-QSI candidate. Exact diagonalization and molecular-dynamics comparisons concluded that the quantum structure factor of the \(\pi\)-flux octupolar QSI regime is most compatible with neutron data, and that nearest-neighbor interactions suffice once quantum fluctuations are included, without requiring sizable further-neighbor couplings [2201.00828]. PFFRG studies place the experimentally relevant parameter sets inside the \(\pi\)-flux QSI region of the nearest-neighbor \(J_y\)-dominated phase diagram [2311.04269]. Field studies further interpret Ce\(_2\)Zr\(_2\)O\(_7\) as an octupolar \(U(1)\) QSL that undergoes an Anderson–Higgs transition in sufficiently large fields, with octupolar magnons that are invisible to neutrons but contribute to the heat capacity [2209.04590].

Ce\(_2\)Hf\(_2\)O\(_7\) was analyzed by neutron scattering and finite-temperature Lanczos fits, which found best-fit parameter sets with dominant octupolar exchange, significant but smaller dipolar exchange, and small negative \(J_{xz}\). The resulting phase identification is an octupolar QSI in the \(\pi\)-flux sector, with a low-energy spinon continuum centered at \(\Delta=0.024\pm0.002\) meV and a quasi-elastic low-\(Q\) dipolar signal. The work described the material as a “quantum multipolar liquid” because both dipolar and octupolar correlations are experimentally visible [2305.08261].

Ce\(_2\)Sn\(_2\)O\(_7\) remains more contentious. Earlier powder neutron data on solid-state samples showed high-\(|Q|\) diffuse scattering associated with magnetic octupoles and motivated an octupolar-QSI interpretation. By contrast, a later study of hydrothermally grown single crystals and powders fitted heat capacity, susceptibility, and diffuse neutron data to the nearest-neighbor DO XYZ model and placed the material inside the dipolar all-in–all-out Néel phase, with a finite-temperature proximate dipolar spin-ice regime above an expected low-temperature transition [2211.15140]. This contrast is a materials issue rather than a formal ambiguity of the Hamiltonian.

Ordered DO pyrochlores show a different use of the same framework. In Nd\(_2\)Hf\(_2\)O\(_7\), the rotated DO XYZ Hamiltonian with
\[
\widetilde J_{\widetilde x}=0.106(5)\ {\rm meV},\quad
\widetilde J_{\widetilde y}=0.008(5)\ {\rm meV},\quad
\widetilde J_{\widetilde z}=-0.057(4)\ {\rm meV}
\]
accounts for all-in–all-out order together with a flat, gapped pinch-point mode at \(0.094\) meV and dispersive half-moon features, which were interpreted as signatures of magnetic fragmentation [2208.12369]. In Sm\(_2\)Ti\(_2\)O\(_7\), inelastic neutron scattering identified a DO Ising doublet, neutron diffraction found all-in–all-out order below \(T_N=0.35\) K with ordered moment \(0.44(7)\,\mu_B\), and \(\mu\)SR showed persistent low-energy spin dynamics, making the material a candidate for moment fragmentation physics [1805.09472].

## 6. Magnetic fields, advanced probes, and current theoretical extensions

Magnetic-field response is highly anisotropic in DO pyrochlores because the field couples strongly only to the dipolar component projected onto the local \(\hat z_i\) axes. For fields along \([110]\), \([111]\), and \([001]\), the sublattice-dependent projections differ qualitatively, leading to distinct phase diagrams and structure-factor evolution [2406.18650], [2502.14067]. In GMFT calculations for Ce\(_2\)Zr\(_2\)O\(_7\)-like parameters, the \(\pi\)-flux QSI is destabilized by comparatively small critical fields, while Ce\(_2\)Hf\(_2\)O\(_7\)-like parameters yield a more robust \(\pi\)-flux QSI and, near the Ising limit under \([110]\) field, even allow staggered-flux states such as \((0,\pi,\pi,0)\) [2406.18650], [2502.14067].

A recent neutron study under \([111]\) field introduced a same-temperature high-field subtraction protocol to separate photon and spinon sectors in Ce\(_2\)Zr\(_2\)O\(_7\). Weak fields of about \(0.15\) T suppress the low-energy photon weight while leaving the higher-energy spinon continuum robust, albeit hardened; GMFT and exact diagonalization were used to interpret this as spectroscopic demarcation of emergent photons and spinons in a \(\pi\)-flux QSI [2601.03202]. This is a notable development because it uses the selective Zeeman coupling of the DO Hamiltonian as a spectroscopic control parameter.

Nonlinear spectroscopy has also been proposed as a probe of the DO Hamiltonian. A two-dimensional coherent spectroscopy study of octupolar QSI showed that, in the intermediate-temperature window
\[
|J_{\pm}|^3/J_{yy}^2,\ |J_{\pm\pm}|^3/J_{yy}^2 \ll T \ll |J_{\pm}|,
\]
spinons remain quantum coherent but move in an incoherent ice background, producing a broad response. At lower temperature, a sharp rephasing signal emerges and distinguishes \(0\)-flux from \(\pi\)-flux QSI [2406.01472]. This suggests that the hierarchy of scales built into the DO Hamiltonian can be accessed by probes other than neutron scattering.

Theoretical generalizations now extend beyond bosonic spinons. A fermionic parton mean-field treatment classified 12 fully symmetric uniform \(U(1)\) QSLs for DO pyrochlores, including four “monopole-flux” states. Several of these states show linear low-temperature specific heat, while others show \(C\sim T^\alpha\) with \(\alpha\) close to \(1\); the work connected these results to the metallic specific-heat response in Nd\(_2\)ScNbO\(_7\) [2312.03106]. A plausible implication is that the dipolar–octupolar Hamiltonian is not tied to a single parton language: its restricted symmetry structure admits both gauge-theoretic and fermionic-spinon descriptions, with experimentally distinguishable consequences.

Across these developments, the central lesson remains stable. The dipolar–octupolar spin Hamiltonian is minimal at the level of symmetry, but not at the level of phenomenology. Its single octupolar channel, its basis-dependent dipole–octupole decomposition, and its selective coupling to fields and neutrons together generate a phase structure and spectroscopy that differ qualitatively from conventional pyrochlore anisotropic exchange models [2002.05608], [2201.00828], [2406.18650].

Source: https://www.emergentmind.com/topics/dipolar-octupolar-spin-hamiltonian