Magnetoelectric Multipole Formalism

• Magnetoelectric multipole formalism is a symmetry-based framework that defines coupled electric and magnetic multipole moments using tensor decompositions.
• It underpins analysis of spontaneous multipolar ordering and non-linear magnetoelectric responses in systems breaking inversion and time-reversal symmetry.
• The approach connects microscopic Hamiltonians with macroscopic observables, aiding experimental probes via neutron scattering and transport measurements.

A magnetoelectric multipole formalism provides a rigorous microscopic and symmetry-based framework for describing coupled magnetic and electric multipole moments, focusing on the operators, tensor structures, selection rules, and physical consequences of moments that govern the linear and non-linear magnetoelectric responses of complex materials. These multipoles—monopole (often called "axion"), toroidal dipole, quadrupole, and higher-rank moments—are crucial in systems that break both spatial inversion ($\mathcal{P}$) and time-reversal ($\mathcal{T}$) symmetries, leading to emergent collective behaviors such as spontaneous multipolar order, unusual magnetoelectric coupling, topological responses, and symmetry-protected surface phenomena.

1. Fundamental Magnetoelectric Multipole Operators and Tensor Decomposition

The magnetoelectric multipole tensor is most naturally defined as the rank-2 moment of the magnetization density: $M_{ij} = \int r_i\, \mu_j(\mathbf{r})\, d^3r$ where $\mu_j(\mathbf{r})$ is the local magnetization. This tensor can be decomposed into irreducible symmetry components:

• Monopole (pseudoscalar, axion): $a = \frac{1}{3} M_{kk}$, odd under $\mathcal{P}$ and $\mathcal{T}$.
• Toroidal dipole (axial vector): $t_k = \frac{1}{2}\varepsilon_{kij} M_{ij}$, also odd under both $\mathcal{P}$ and $\mathcal{T}$.
• Quadrupole (symmetric, traceless tensor):

$$q_{ij} = \frac{1}{2}(M_{ij} + M_{ji}) - \frac{1}{3}\delta_{ij}M_{kk}$$

This carries the five independent components of a second-rank irreducible tensor and is also odd under both symmetries.

Explicit operator representations link these moments to spin and position: $$\mathbf T = \frac{1}{2}\sum_i \mathbf r_i \times \mathbf s_i, \quad Q_{\alpha\beta} = \sum_i \left(r_{i\alpha} s_{i\beta} + r_{i\beta} s_{i\alpha}\right)$$ Higher-order multipoles, e.g., octupoles, involve moments such as $$\mathcal{M}_{ijk} = \int \mu_i(\mathbf{r}) r_j r_k\, d^3r$$, whose decomposition yields the rank-3 irreducible representation consistent with SO(3) tensor algebra (Urru et al., 2022, Urru et al., 2022, Thöle et al., 2018).

2. Microscopic Lattice Models and Spontaneous Multipole Ordering

A canonical microscopic route starts from Hamiltonians with explicit local parity mixing, such as the periodic Anderson model with sublattice-dependent hybridization: $$\mathcal H_1 = \sum_{l,\mathbf k,\sigma}V_l(\mathbf k)\bigl(c^\dagger_{l\mathbf k\sigma}f_{l\mathbf k\sigma}+\mathrm{H.c.}\bigr) +\sum_{l,\mathbf k}\mathbf g^{cf}_l(\mathbf k)\cdot\mathbf s^{cf}_l(\mathbf k)$$ On a zig-zag chain, local parity mixing yields $$\mathbf{g}_l^{cf}(\mathbf{k}) = (-1)^l g \hat{\mathbf{z}} \sin k_x$$ distinguishing between A/B sublattices (Hayami et al., 2015).

Applying a Schrieffer-Wolff transformation yields an extended Kondo lattice model with antisymmetric exchange interactions: $$\mathcal{H}_{\rm ex-KLM} = \cdots + \sum_l D_l\text{ --exchange} + G\text{ --exchange terms}$$ where $D_l = (-1)^l \sqrt{J G}$ encodes sublattice-odd exchange, directly linked to local parity mixing.

At half-filling, variational, simulated annealing, and Monte Carlo simulations demonstrate a stable $z$-UD (up/down) ordering on the zigzag chain, carrying both odd-parity toroidal ($T_x$) and quadrupolar ($Q_{yz}$) order: $$T_x\sim\sum_l(-1)^l S^z_l,\qquad Q_{yz}\sim\sum_l(-1)^l y_l S^z_l$$ This state breaks inversion but preserves combined $\mathcal{PT}$, and the toroidal component $T_x \neq 0$ implies a robust spontaneous odd-parity multipolar phase (Hayami et al., 2015).

3. Linear Magnetoelectric Response and Tensor Structure

The magnetoelectric effect is formally captured by the linear-response tensor $\alpha_{ij}$ coupling electric field to magnetization: $m^z_{\rm s} = \alpha_{zx} E_x$ with the Kubo formula: $$\alpha_{zx} = \lim_{\omega\to0}\frac{1}{i\omega}\frac{1}{V_s}\sum_{m,n,\mathbf k} \frac{f(\varepsilon_{n\mathbf k})-f(\varepsilon_{m\mathbf k})}{\varepsilon_{n\mathbf k}-\varepsilon_{m\mathbf k}}\sigma^z_{\mathbf k}{}^{nm}J_{x,\mathbf k}{}^{mn}\frac{1}{\varepsilon_{n\mathbf k}-\varepsilon_{m\mathbf k}+i0^+}$$ In the $z$-UD multipole background, antisymmetric exchange ($D_l$) is directly tied to nonzero off-diagonal $\alpha_{zx}$, evidence of the fundamental role of the toroidal multipole in generating magnetoelectricity when both $\mathcal{P}$ and $\mathcal{T}$ are locally broken (Hayami et al., 2015, Thöle et al., 2018).

4. Emergence, Band Structure, and Experimental Signatures

The presence of spontaneous odd-parity multipole order induces characteristic signatures:

• Band-structure deformation: including a band-bottom shift proportional to $\sqrt{J G} \sin k_x$
• Phase diagram: Multipole-ordered phases appear over a wide range near half filling ($n_e=1$), with finite $T_c$ confirmed numerically
• Response functions: Direct correspondence between the emergence of $T_x$ and finite $\alpha_{zx}$; changes in magnetic or transport properties act as indirect probes.

In experimental contexts, the consequences include nontrivial electromagnetic responses such as nonlinear polarization, band splittings, and field-orientation-dependent magnetization. Neutron and x-ray scattering, particularly with polarization analysis, couple directly to the spatial structure of ME multipoles (Urru et al., 2022, Lovesey, 2014).

5. Symmetry, Selection Rules, and Higher-Rank Generalizations

Symmetry dictates the allowed magnetoelectric multipoles and their physical observability:

• Only sites lacking spatial inversion and time-reversal may host nonzero $a$, $\mathbf t$, or $q_{ij}$.
• Odd-rank tensors (e.g., octupoles) enter higher-order (second-order) ME effects, with selection rules requiring ferroic alignment of local octupoles for macroscopic response (Urru et al., 2022).
• Multipole orders in real or momentum space distinctly pattern spontaneous ME responses and can be classified in tensor tables according to system topology and texture, as in skyrmion-hosting materials (Bhowal et al., 2022).

The formalism extends naturally to finite systems (clusters, heterostructures), surfaces, and interfaces, with boundary conditions and surface multipolization paralleling those for classical polarization in ferroelectrics (Spaldin, 2020).

6. Principal Theoretical Implications and Outlook

The formalism developed from local parity mixing to sublattice-odd antisymmetric exchange derives direct links between microscopic symmetry breaking and emergent macroscopic ME effects. It establishes a unifying operator/tensor language—grounded in quantum mechanics and symmetry representation theory—for describing and classifying all symmetry-allowed cross-correlated ME phases. The framework is directly extendable to multi-band, spatially inhomogenous, time-dependent, and driven systems, relevant to quantum materials, engineered heterostructures, and device applications.

Multipole-based formulations are critical for bridging ab initio electronic structure, effective model Hamiltonians, and phenomenological response tensors in magnetic, multiferroic, and low-dimensional quantum matter. As experimental techniques and first-principles modeling converge on probing "hidden order" parameters, the rigorous identification and control of magnetoelectric multipoles is central to future progress in fundamental and applied quantum materials research (Hayami et al., 2015, Lovesey, 2014, Thöle et al., 2018, Urru et al., 2022).