Orbital-Moment Octupole in Crystals

• Orbital-moment octupole is a rank-3 magnetic multipole defined by gauge-invariant expressions using quantum-geometric quantities, such as Berry curvature and quantum metric, in periodic crystals.
• The formulation bridges local time-reversal-odd octupolar operators and transport coefficients, enabling higher-rank Hall responses under spatially nonuniform electric fields with precise symmetry constraints.
• Experimental realizations in altermagnets, magnetic multilayers, and Mott systems provide platforms for tuning and detecting octupole contributions via advanced transport and spectroscopic techniques.

Orbital-moment octupole denotes a rank-3 magnetic multipole associated with orbital degrees of freedom. In current usage, the term covers several closely related but not identical constructions: a bulk orbital magnetic octupole tensor $M_{ijk}$ in periodic crystals; local time-reversal-odd octupolar operators such as $T_i^y$ or $T_{xyz}$ in orbital and spin-orbit-entangled manifolds; and octupolar contributions to anomalous Hall transport extracted from symmetry-resolved angular decompositions. A central recent development is the formulation of a gauge-invariant bulk expression for the orbital magnetic octupole in crystalline solids and its connection to a higher-rank Hall response induced by spatially nonuniform electric fields, including cases in which symmetry forbids the conventional anomalous Hall effect against uniform electric fields (Sato et al., 30 Dec 2025).

1. Terminology and scope

In the periodic-crystal setting, the orbital magnetic octupole is defined thermodynamically as the response of the free-energy density to the second spatial derivative of a magnetic field. This places it in the same hierarchy as dipole and quadrupole magnetic moments, but at the next rank. The resulting object is a rank-3 axial tensor $M_{ijk}$ (Sato et al., 30 Dec 2025).

In local orbital models, the same octupolar designation is often attached to a time-reversal-odd operator acting within a restricted orbital manifold. In the spin-less doubly degenerate $e_g$ model, the site-octupole operator is identified as

$$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$

whose eigenstates are the complex orbitals

$$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$

and which transforms with $A_{2g}$ octupolar symmetry (Nasu et al., 2010). In the $J=2$ non-Kramers-doublet setting, the fully symmetrized octupole is

$$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$

and projection into the low-energy doublet gives

$T_i^y$0

(Paramekanti et al., 2019).

A further symmetry-based usage appears in dipole-octupole doublets on the triangular lattice. There, $T_i^y$1 is invariant under the two-fold rotation $T_i^y$2, odd under time reversal, and is identified as the component of a rank-3 octupolar moment, whereas $T_i^y$3 transform as the two components of an ordinary magnetic dipole (Li et al., 2016). In a surface-state $T_i^y$4 model on the simple cubic lattice, the single-site octupole operator is written as

$T_i^y$5

with $T_i^y$6 mixing the two $T_i^y$7 orbitals with an imaginary coefficient (Kubo, 2024).

A common source of confusion is that not every octupole-related transport theory starts from the same microscopic object. In the multilayer anomalous-Hall literature, the “octupole term” is a rank-4 tensor $T_i^y$8 entering the magnetization-direction expansion

$T_i^y$9

where $T_{xyz}$0 is symmetric and traceless in its last three indices. That work explicitly does not write down an explicit second-quantized or Bloch-basis operator $T_{xyz}$1, and instead parameterizes the Hall conductivity by the coefficients $T_{xyz}$2 (Niu et al., 26 Aug 2025). The bulk-crystal tensor $T_{xyz}$3 and the transport-expansion tensor $T_{xyz}$4 therefore belong to different but related frameworks.

2. Gauge-invariant bulk formulation in periodic crystals

For a periodic crystal at finite chemical potential $T_{xyz}$5 and temperature $T_{xyz}$6, the orbital magnetic octupole is given by the gauge-invariant bulk expression

$T_{xyz}$7

with

$T_{xyz}$8

The band-resolved tensors $T_{xyz}$9, $M_{ijk}$0, and $M_{ijk}$1 are fully gauge-invariant rank-3 axial tensors built out of the Berry curvature $M_{ijk}$2, quantum metric $M_{ijk}$3, orbital moment $M_{ijk}$4, and higher derivatives of the Hamiltonian (Sato et al., 30 Dec 2025).

The geometric building blocks are

$M_{ijk}$5

$M_{ijk}$6

together with the covariant derivative

$M_{ijk}$7

An explicit example is

$M_{ijk}$8

The expressions for $M_{ijk}$9 and $e_g$0 are likewise completely gauge-invariant.

The gauge structure is central. Under a local phase redefinition $e_g$1, all terms remain unchanged. The position operator never appears explicitly; only derivatives with respect to $e_g$2 enter. Gauge invariance is guaranteed by using projectors $e_g$3 and covariant derivatives (Sato et al., 30 Dec 2025).

3. Thermodynamic definition and response-theory derivation

The thermodynamic starting point is the free-energy density $e_g$4 in a slowly varying magnetic field $e_g$5. Expanding $e_g$6 in gradients of $e_g$7 introduces the dipole, quadrupole, and octupole magnetic moments. By definition,

$e_g$8

This identifies the orbital magnetic octupole as the coefficient conjugate to the second spatial derivative of the magnetic field (Sato et al., 30 Dec 2025).

A Maxwell relation makes the calculation algebraically simpler: $e_g$9 The derivation then proceeds through linear response for the density variation $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$0. One computes the density-current correlation $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$1 up to $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$2, expands $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$3, and uses

$$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$4

to isolate the contribution proportional to $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$5. Integrating the $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$6-derivative and reorganizing the result yields the final gauge-invariant bulk formula (Sato et al., 30 Dec 2025).

This framework clarifies the status of $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$7. It is not introduced as a phenomenological fitting parameter but as a thermodynamic quantity derivable from Bloch-state response theory. The formulation also makes explicit that the bulk theory relies on Bloch periodicity on the Brillouin zone and on a representation in terms of quantum-geometric objects rather than real-space position operators. A plausible implication is that the formulation places the orbital magnetic octupole on the same conceptual footing as other modern-band-theory observables built from Berry curvature, quantum metric, and orbital moments.

4. Symmetry constraints and model realizations in crystalline solids

A minimal realization of a nonzero orbital magnetic octupole is provided by the two-sublattice, collinear antiferromagnet described as a “$$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$8-wave altermagnet” on a tetragonal lattice: $$\hat O_i \equiv T_i^y =\frac1{2i}\bigl(c_{iu}^\dagger c_{iv}-c_{iv}^\dagger c_{iu}\bigr),$$9 Here $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$0 act in sublattice space, $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$1 in spin space, and

$$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$2

$$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$3

The term $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$4 breaks time-reversal but preserves mirror symmetries $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$5 and $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$6. As a result, the only nonzero octupole components are

$$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$7

Numerical evaluation shows that both vanish for $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$8 and grow roughly linearly with spin-orbit coupling $$|\,\pm\>\;=\;\tfrac1{\sqrt2}\bigl(d_u\pm i\,d_v\bigr),$$9, which highlights the relativistic origin of the orbital contribution (Sato et al., 30 Dec 2025).

A distinct realization arises in $A_{2g}$0-symmetric magnetic multilayers, where the anomalous Hall response is decomposed into dipole and octupole harmonics. For in-plane magnetization $A_{2g}$1 and current $A_{2g}$2 along $A_{2g}$3,

$A_{2g}$4

$A_{2g}$5

The first two coefficients are

$A_{2g}$6

Here $A_{2g}$7 is the conventional dipole-dominated anomalous Hall effect with a small octupole correction, whereas $A_{2g}$8 is exclusively the third-order octupole contribution. In $A_{2g}$9, only combinations such as

$J=2$0

can be nonzero, and they reduce to two independent parameters in the Hall geometry (Niu et al., 26 Aug 2025).

The multilayer calculations use a fully-relativistic, spin-orbit-coupled exact muffin-tin orbitals scattering-wave-function method for alternating bcc-Fe and fcc-Ag layers, sandwiched between semi-infinite Ag leads, with the Ag lattice rotated $J=2$1 to obtain $J=2$2 in-plane symmetry. Charge currents between atoms $J=2$3 are evaluated from

$J=2$4

and a dense $J=2$5 $J=2$6 mesh in the 2D Brillouin zone ensures convergence of the transport integrals. In $J=2$7 multilayers, the pure-octupole term $J=2$8 reaches up to $J=2$9 in the $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$0 ferromagnet, while in a simple Fe$$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$1Ag interface it is only $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$2. The relative octupole strength $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$3 changes from $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$4 to $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$5 to $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$6 as $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$7 varies from $$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$8, demonstrating tunability by layer number, and the “$$T_{xyz}\equiv \overline{J_xJ_yJ_z} =\frac1{6}\sum_{\pi\in S_3}J_{\pi(x)}J_{\pi(y)}J_{\pi(z)},$$9-AFM” configuration reverses the sign of both $T_i^y$00 and $T_i^y$01 (Niu et al., 26 Aug 2025).

5. Higher-rank Hall response and experimental access

The principal transport consequence of the bulk orbital magnetic octupole is a higher-rank Hall response driven by a spatially nonuniform electric field. Since a spatially nonuniform electric field can be written as a gradient of a local chemical potential, $T_i^y$02, the octupole induces a transverse current under $T_i^y$03. The corresponding “octupolar anomalous Hall conductivity” is

$T_i^y$04

In real space, the leading linear response is

$T_i^y$05

This relation is the higher-rank analogue of the ordinary Středa formula $T_i^y$06 (Sato et al., 30 Dec 2025).

This framework has a direct symmetry implication for altermagnets. Mirror or rotation symmetries can forbid the conventional anomalous Hall effect $T_i^y$07 against uniform electric fields while still allowing the octupolar conductivity $T_i^y$08. In that situation, the first nonzero Hall response under a spatially varying $T_i^y$09 is purely octupolar. The proposed detection route is to engineer a controlled electric-field gradient, for example via split-gate electrodes or near-field optical excitation, and measure the resulting transverse current. A static alternative is the quadrupolar magnetoelectric effect $T_i^y$10 with

$T_i^y$11

so that applying a magnetic-field gradient generates an electric polarization with the symmetry of $T_i^y$12. Because $T_i^y$13 is linear in $T_i^y$14 within an insulating gap, one expects a strictly linear-in-$T_i^y$15 plateau in the octupolar response versus doping (Sato et al., 30 Dec 2025).

The multilayer literature provides a complementary experimental strategy based on angle-resolved transport. Because $T_i^y$16 is a pure octupole contribution, measuring the $T_i^y$17 and $T_i^y$18 components of the Hall current completely isolates the octupole piece. This method explicitly incorporates discrete crystal symmetries and is stated to enable the investigation of octupole contribution in non-periodic systems, particularly at interfaces and surfaces (Niu et al., 26 Aug 2025).

A recurrent misconception is that anomalous Hall phenomena in magnetic conductors are exhausted by the dipole moment or net magnetization. The multilayer results explicitly state that even the conventional contribution arises not only from the dipole moment, and the crystalline-octupole theory shows that a Hall response can persist in a higher-rank form when the uniform-field anomalous Hall channel is symmetry-forbidden (Niu et al., 26 Aug 2025, Sato et al., 30 Dec 2025).

6. Correlated-orbital, surface, and collective manifestations

In Mott-insulating orbital models, the orbital-moment octupole appears as an ordered degree of freedom rather than as a band-structure response coefficient. In the doubly degenerate $T_i^y$19 model, fourth-order ring exchange generates the plaquette interaction

$T_i^y$20

with $T_i^y$21. Because $T_i^y$22 naturally involves the imaginary combination $T_i^y$23, it directly couples to the time-reversal-odd operator $T_i^y$24, identified as the magnetic octupole. The reported phase structure includes quadrupole order up to $T_i^y$25, onset of genuine $T_i^y$26 order at $T_i^y$27, and collapse of both quadrupole and octupole moments by $T_i^y$28 because of strong fluctuations (Nasu et al., 2010).

In $T_i^y$29 double perovskites, the octupole resides in a low-lying non-Kramers doublet of a $T_i^y$30 ion. The doublet carries purely quadrupolar and octupolar moments and no matrix elements of $T_i^y$31 within it. A second-order perturbative mechanism from Heisenberg exchange and orbital repulsion produces the ferro-octupolar coupling

$T_i^y$32

The same work identifies experimental signatures in $T_i^y$33SR, Raman scattering, X-ray diffraction, inelastic neutron scattering, and a gapped, dispersive magnetic exciton that can condense into conventional type-I antiferromagnetic order when $T_i^y$34 (Paramekanti et al., 2019).

In dipole-octupole doublets on the triangular lattice, octupolar order may be hidden from conventional magnetization measurements. The nearest-neighbor pseudospin Hamiltonian is

$T_i^y$35

After an $T_i^y$36 rotation in pseudospin space, $T_i^y$37 remains purely octupolar while $T_i^y$38 form the dipole doublet that couples linearly to an external field. The model supports ferro-octupolar order with

$T_i^y$39

and antiferro-octupolar three-sublattice order. In the ferro-octupolar phase, the octupolar susceptibility $T_i^y$40 diverges at $T_i^y$41, whereas the dipolar susceptibilities remain finite. The predicted experimental consequences include finite Kerr rotation below the ordering temperature and sharp octupolar-wave modes in inelastic neutron scattering (Li et al., 2016).

Surface-localized orbital octupoles also occur in a purely orbital tight-binding setting. For the simple-cubic $T_i^y$42 model with finite thickness, (110) and (111) surfaces host electronic states localized around the surfaces, and the (111) surface supports a perfectly flat band when the bulk band projected onto the surface Brillouin zone is gapped. The surface-state spinor is proportional to the $T_i^y$43 direction, $T_i^y$44, and the expectation value

$T_i^y$45

shows that the surface-state manifold carries a quantized octupole moment density of unit magnitude on the top surface and the opposite sign on the bottom (Kubo, 2024).

A different use of orbital moment appears in the octupole collective Hamiltonian. In the intrinsic-frame description of octupole vibrations, the operators

$T_i^y$46

satisfy the usual angular-momentum algebra and are identified as an intrinsic orbital moment carried by the octupole vector vibrations $T_i^y$47. This establishes a collective-mechanical notion of orbital moment that is distinct from the Bloch-band and local-orbital octupoles discussed above (Rohozinski et al., 2018).

Taken together, these results show that orbital-moment octupoles are not restricted to a single formalism. They appear as gauge-invariant bulk tensors in periodic crystals, as local order parameters in orbital and spin-orbit-entangled Mott systems, as symmetry-resolved transport coefficients in interfacial anomalous Hall effects, as quantized surface-state moments in $T_i^y$48 bands, and as intrinsic orbital moments of octupole vibrations in collective Hamiltonians. This suggests that the unifying content of the term is not a unique operator but the organization of time-reversal-odd, higher-rank orbital structure and its electromagnetic or dynamical consequences.