---
title: Local Orbital Magnetization
url: https://www.emergentmind.com/topics/local-orbital-magnetization
type: topic
---

# Local Orbital Magnetization

Local orbital magnetization denotes a real-space-resolved orbital magnetic response whose spatial average reproduces the macroscopic orbital magnetization of an insulating or otherwise bounded electronic system. In the modern literature it appears in several closely related but non-identical forms: as a projector-based local marker built from the ground-state density matrix, as a single-site quantity decomposed into local orbital moment and local Berry-curvature terms, and as a current-generating field defined by local-flux response of the grand potential. These formulations are used to analyze boundary currents, local Chern structure, sublattice textures, orbital magnetoelectricity, and finite-size orbitronics in topological, trivial, crystalline, amorphous, fractal, and correlated systems [1301.6275][2208.11285][2512.12343][2606.29330].

## 1. Formal definitions and decompositions

A central projector-based construction begins from the occupied-state projector
$$
\hat P=\sum_{\varepsilon_n<\mu}|\phi_n\rangle\langle \phi_n|,\qquad \hat Q=1-\hat P.
$$
For a non-interacting insulator, Bianco and Resta rewrote the bulk orbital magnetization as the spatial average of a local real-space density. Their local magnetization density is
$$
\mathcal M_1(\mathbf r)=\frac{e}{\hbar c}\,\mathrm{Im}\Bigl\{
\langle \mathbf r|\hat P\,\hat x\,\hat Q\,\hat H\,\hat Q\,\hat y\,\hat P|\mathbf r\rangle
-\langle \mathbf r|\hat Q\,\hat x\,\hat P\,\hat H\,\hat P\,\hat y\,\hat Q|\mathbf r\rangle
\Bigr\},
$$
with
$$
M_1=\frac1A\int d^dr\,\mathcal M_1(\mathbf r).
$$
For a Chern insulator one adds a $\mu$-dependent Středa term through the local Chern marker
$$
\mathcal F(\mathbf r)=4\pi\,\mathrm{Im}\,\langle \mathbf r|\hat P\,\hat x\,\hat Q\,\hat y\,\hat P|\mathbf r\rangle,
$$
so that
$$
\mathcal M_2(\mathbf r)=\frac{\mu}{\phi_0}\,\mathcal F(\mathbf r),
\qquad
M(\mu)=\frac1A\int d^dr\,\bigl[\mathcal M_1(\mathbf r)+\mathcal M_2(\mathbf r)\bigr].
$$
This formulation is explicitly real-space and is intended to express orbital magnetization as a genuine bulk property of insulators [1301.6275].

A closely related decomposition, used in finite quantum anomalous Hall bilayers, writes a local marker
$$
m(\mathbf r)=m_{\mathrm{LC}}(\mathbf r)+m_{\mathrm{IC}}(\mathbf r)+m_{\mathrm{BC}}(\mathbf r),
$$
with LC, IC, and BC denoting local circulation, itinerant circulation, and Berry-curvature contributions. In that notation,
$$
m_{\mathrm{LC}}(\mathbf r)\equiv \mathrm{Im}\,\langle \mathbf r|P\,x\,Q\,H\,Q\,y\,P|\mathbf r\rangle,
$$
$$
m_{\mathrm{IC}}(\mathbf r)\equiv -\mathrm{Im}\,\langle \mathbf r|Q\,x\,P\,H\,P\,y\,Q|\mathbf r\rangle,
$$
$$
m_{\mathrm{BC}}(\mathbf r)\equiv -\frac{e}{2\pi c}\,c(\mathbf r),
\qquad
c(\mathbf r)\equiv 4\pi\,\mathrm{Im}\,\langle \mathbf r|Q\,x\,P\,y\,Q|\mathbf r\rangle.
$$
The corresponding sample-averaged terms are
$$
M_{\mathrm{LC}}=\frac{e}{2cA}\,\mathrm{Im}\,\mathrm{Tr}\{P\,x\,Q\,H\,Q\,y\,P\},
$$
$$
M_{\mathrm{IC}}=-\frac{e}{2cA}\,\mathrm{Im}\,\mathrm{Tr}\{Q\,x\,P\,H\,P\,y\,Q\},
$$
$$
M_{\mathrm{BC}}=-\frac{e}{2\pi cA}\,\mathrm{Tr}[c(\mathbf r)]=-\frac{e}{2\pi c}\,C.
$$
This three-term split is gauge invariant term by term and is particularly useful for diagnosing boundary effects in finite Chern systems [2208.11285].

A later microscopic theory defines the local orbital magnetization density from the local grand-potential density,
$$
\Xi(\mathbf r,\mu,T)=-\int_{-\infty}^{+\infty}d\varepsilon\,F(\varepsilon)\,\nu(\mathbf r,\varepsilon),
$$
with
$$
\nu(\mathbf r,\varepsilon)=\sum_n\delta(\varepsilon-\varepsilon_n)\,|\langle \mathbf r|n\rangle|^2,
\qquad
F(\varepsilon)=k_BT\ln\!\bigl[1+e^{-\beta(\varepsilon-\mu)}\bigr].
$$
The local orbital magnetization is then
$$
\mathbf M(\mathbf r)=\frac e\hbar\sum_n\Bigl[f(\varepsilon_n)\,\mathbf m_n(\mathbf r)+F(\varepsilon_n)\,\mathbf\Omega_n(\mathbf r)\Bigr],
$$
where $\mathbf m_n(\mathbf r)$ is a local orbital magnetic moment density and $\mathbf\Omega_n(\mathbf r)$ is an effective local Berry-curvature density. This formulation is explicitly state resolved and extends local orbital magnetization to single-site and sublattice texture analysis [2512.12343].

## 2. Nearsightedness, bulk character, and boundary compensation

In the projector formalism, locality is tied to the exponential decay of the one-particle density matrix in a gapped system. Both $P(\mathbf r,\mathbf r')$ and $Q(\mathbf r,\mathbf r')$ decay exponentially for $|\mathbf r-\mathbf r'|$ beyond the localization length $\xi$, so each matrix element entering $\mathcal M(\mathbf r)$ depends only on a finite neighborhood of $\mathbf r$. This is the sense in which the local magnetization density is “nearsighted.” Boundary conditions change $P(\mathbf r,\mathbf r')$ only when either argument lies within $O(\xi)$ of the boundary, and therefore a large central-region average is insensitive to whether one uses open or periodic boundaries [1301.6275].

The finite-size role of boundaries was examined numerically in a bilayer quantum anomalous Hall model with adjustable Chern number at half filling. For an $N\times N$ bilayer with open boundaries, three prescriptions were compared: a $k$-space Brillouin-zone formula, a real-space bulk average over an inner block, and a real-space whole-sample average over all sites. In every topological sector studied, $C=0,1,2$, all three methods converge to the same value up to $10^{-4}$ as $N\to 30$–$40$ [2208.11285].

The technically important subtlety is the behavior of the Berry-curvature term. For nonzero Chern number, the bulk average retains
$$
M_{\mathrm{BC}}^{(b)}=-\frac{e}{2\pi c}\,C,
$$
whereas the whole-sample average of the BC term vanishes identically because the local Chern marker in the rim carries the opposite sign and exactly cancels the bulk. Taken alone, that would appear to spoil the whole-sample average. The numerical result is that the rim simultaneously induces a nonzero boundary contribution to $M_{\mathrm{LC}}+M_{\mathrm{IC}}$ which exactly compensates the lost BC piece,
$$
A_{\mathrm{bdy}}(M_{\mathrm{LC}}+M_{\mathrm{IC}})
\to +\,M_{\mathrm{BC}}^{(b)}.
$$
Consequently,
$$
M^{(w)}=[M_{\mathrm{LC}}+M_{\mathrm{IC}}]^{(w)}+M_{\mathrm{BC}}^{(w)}=M^{(b)}.
$$
The same work showed that the slope $dM/d\mu=-(e/2\pi c)\,C$ is captured by both real-space averaging methods, and that even very small clusters, $N=5,10$, agree perfectly with the whole-sample real-space formula when compared with the first-principles definition
$$
M=-\frac{e}{2cA}\sum_{i,\mathrm{occupied}}\langle \psi_i|\,\mathbf r\times \mathbf v\,|\psi_i\rangle.
$$
This establishes that excising the rim is not required, even at quite small sizes, for the total orbital magnetization in that setting [2208.11285].

## 3. Competing local formulations and their relations

The phrase “local orbital magnetization” now covers several constructions with different defining properties.

| Formulation | Defining quantity | Distinctive statement |
|---|---|---|
| Bianco–Resta marker | $\mathcal M(\mathbf r)$ from projectors $\hat P,\hat Q$ and local Chern marker | Nearsighted bulk density, boundary-condition independent in the bulk [1301.6275] |
| Local Berry-curvature theory | $\mathbf M(\mathbf r)=\frac e\hbar\sum_n[f\,\mathbf m_n(\mathbf r)+F\,\mathbf\Omega_n(\mathbf r)]$ | Reveals topological and geometric sublattice textures; may differ from Bianco–Resta on individual sites [2512.12343] |
| Current-generating theory | $m_{\mathrm{orb}}(\mathbf r)$ from local-flux response and natural local marker $\mathcal M_i(\mathbf r)$ | Generates equilibrium current; site-position moment equals the orbital magnetic quadrupole moment [2606.29330] |

In the single-site theory, each occupied state contributes a local orbital magnetic moment density
$$
\mathbf m_n(\mathbf r)=|\langle \mathbf r|n\rangle|^2
\Bigl[-\tfrac12\,\mathrm{Im}\,\langle n|\hat{\mathbf r}\times \hat{\mathbf v}|n\rangle\Bigr]
$$
and an effective local Berry-curvature density
$$
\mathbf\Omega_n(\mathbf r)=
\mathrm{Re}\sum_{m\neq n}
\langle \mathbf r|n\rangle
\langle m|\mathbf r\rangle
\frac{\langle n|\hat{\mathbf r}\times\hat{\mathbf v}|m\rangle}{\varepsilon_n-\varepsilon_m}.
$$
The global sum of $\mathbf\Omega_n(\mathbf r)$ vanishes, $\sum_{\mathbf r}\mathbf\Omega_n(\mathbf r)=\mathbf 0$, but locally it carries topological and geometric information. For Bloch electrons this leads to a sublattice-resolved Berry curvature
$$
\mathbf\Omega_{n\mathbf k}(\mathbf r_\alpha)=
\mathbf\Omega_{n\mathbf k}^{\rm topo}(\mathbf r_\alpha)+
\mathbf\Omega_{n\mathbf k}^{\rm geom}(\mathbf r_\alpha),
$$
where the “topo” piece integrates to the usual intraband Berry curvature and the “geom” piece is purely interband and sums to zero over sublattices. Numerically, the $k$-space and $r$-space onsite magnetizations coincide, but they differ from the Bianco–Resta approach, especially in trivial insulators, because the new $\mathbf\Omega_n(\mathbf r)$ retains explicit interband and geometric content [2512.12343].

A different criterion is imposed by the current-generating theory. There the defining relation is
$$
\mathbf j(\mathbf r)=\nabla\times \mathbf m(\mathbf r),
$$
with
$$
\delta\Omega=-\int d^dr\,\mathbf m(\mathbf r)\cdot \delta\mathbf B(\mathbf r).
$$
In two dimensions the local-flux response fixes the local orbital magnetization uniquely; in three dimensions it fixes the field only up to a longitudinal ambiguity. The same work defines a coarse-grained “natural local marker” $\mathcal M_i(\mathbf r)$ which generates the current to third-derivative order and whose first spatial moment reproduces the orbital magnetic quadrupole moment of a finite sample. On the Haldane model, the cell average of the natural local marker matches the bulk magnetization exactly, while markers of Bianco–Resta, Seleznev–Vanderbilt, and Saati–Hur–Pi fail to reproduce the correct local current in the orbital quadrupole-moment state [2606.29330].

A plausible implication is that locality is not a single universal constraint. One formulation may be optimized for bulk sum rules, another for site-resolved geometric texture, and another for exact current generation. The resulting local fields can agree on macroscopic averages while differing qualitatively on individual sites.

## 4. Electric fields, adiabatic evolution, and magnetoelectric response

A semiclassical theory of adiabatically induced orbital magnetization treats a slowly varying Bloch Hamiltonian $H(\mathbf k;\mathbf r,t)$ and identifies extra current terms of the form $\nabla\times M(\mathbf r,t)$ and $\partial_t P(\mathbf r,t)$. In that framework a gauge-invariant, bulk-defined induced orbital magnetization arises only in two cases: when explicit time dependence is absent, giving the orbital magnetoelectric effect, and when one averages over a period, giving pumped orbital magnetization in insulators. The same theory also yields an electric-field-induced intrinsic orbital magnetization in two-dimensional metals and Chern insulators [2012.08750].

For static systems under a perpendicular electric field, the real-space bilayer Chern analysis introduces
$$
M=M_{\mathrm{LC}}+M_{\mathrm{IC}}+M_{\mathrm{CS}}+M_{\mathrm{BC}},
$$
with
$$
M_{\mathrm{CS}}=\frac{eE}{2cA}\,\mathrm{Im}\,\mathrm{Tr}\{P\,x\,P\,y\,P\,z-(x\leftrightarrow y)\},
\qquad
M_{\mathrm{BC}}=-\frac{eE}{cA}\,C.
$$
The local density is correspondingly
$$
m_{\mathrm{orb}}(\mathbf r)=m_{\mathrm{LC}}(\mathbf r)+m_{\mathrm{IC}}(\mathbf r)+m_{\mathrm{CS}}(\mathbf r)+m_{\mathrm{BC}}(\mathbf r).
$$
In the absence of $E$, bulk and whole-sample averages are both correct in the thermodynamic limit. When $E\neq0$ and $C\neq0$, however, the bulk-only average acquires a gauge-dependent shift under $z\to z+\Delta z$,
$$
\Delta M_{\mathrm{bulk}}=\frac{eE}{cA}\,C\,\Delta z,
$$
or, in multilayers,
$$
\Delta M_{\mathrm{bulk}}=\frac{eE}{cA}\sum_i C_i\,\Delta z_i.
$$
The whole-sample average is free of this gauge shift and is therefore the reliable real-space prescription for Chern insulators under finite electric field [2301.06282].

The same finite-field framework defines the orbital magnetoelectric polarizability
$$
\alpha=\frac{\partial M}{\partial E}\Big|_{E\to0}
=\alpha_{\mathrm{LC}}+\alpha_{\mathrm{IC}}+\alpha_{\mathrm{CS}},
$$
with
$$
\alpha_{\mathrm{CS}}=
\frac{e}{2cA}\,\mathrm{Im}\,\mathrm{Tr}\{P\,x\,P\,y\,P\,z-(x\leftrightarrow y)\}.
$$
The numerical finding is that the peaks of both the total OMP $\alpha$ and the Chern–Simons OMP $\alpha_{\mathrm{CS}}$ track the strongest field-induced change in the integrated Berry curvature,
$$
\Delta C/\Delta E=[C(E)-C(0)]/E.
$$
This supports the concrete statement that the stronger the response of Berry curvature to electric field, the stronger is the OMP and the CSOMP [2301.06282].

Adiabatic driving provides a distinct extension. In a two-band honeycomb toy model with a circular optical phonon, the pumped orbital magnetization per cycle is of order the nuclear magneton for typical phonon frequencies, and the predicted signal is comparable to pumped spin magnetization via strong Rashba spin-orbit coupling [2012.08750]. This suggests a direct link between local orbital magnetization theory and phonon angular momentum.

## 5. Model systems, textures, and finite geometries

Honeycomb two-band models provide a controlled setting for sublattice-resolved local orbital magnetization. In the topological Haldane phase with Chern number $\pm1$, one finds $M_A(\mu)\approx M_B(\mu)$ in the gap, interpreted as orbital ferromagnetism, and the slope $\partial_\mu[M_A+M_B]=eC/h$ is quantized, with each sublattice carrying half the Chern number. In the trivial Haldane phase, $M_A(\mu)=-M_B(-\mu)$ in the bands and $M_A+M_B\approx0$ in the trivial gap, giving an almost antiferromagnetic pattern, while each sublattice still shows a linear slope due to the geometric Berry curvature. In the modified Haldane model, $M_A(\mu)=M_B(-\mu)$ and the texture is orbital ferrimagnetic; the nonuniversal slopes in the trivial gap come entirely from $\mathbf\Omega^{\rm geom}$ [2512.12343].

Finite topological structures highlight the role of boundary-localized states. Quantum spin Hall nanoislands support robust orbital edge magnetism when the highest occupied Kramers doublet is singly occupied. For Dirac edge states the orbital moment grows linearly with size,
$$
M_{\max}(L)\simeq \bigl(2.5\,\mu_B/\mathrm{nm}\bigr)\,L,
$$
reaching $\sim42\,\mu_B$ at $L=18\,\mathrm{nm}$. The moment remains nearly full up to $T\sim\Delta_B/2k_B$ and is only weakly affected by Anderson disorder up to $W=0.5\,\mathrm{eV}$, edge-positional disorder, island shape, or crystallographic orientation [1509.02661].

Fractal geometries alter orbital magnetization through their boundary hierarchy. In Sierpinski carpets, higher generations create a dense set of edge states and produce a staircase profile with oscillations in magnetization as a function of chemical potential. In Sierpinski triangles, self-similarity creates distinct fractal-induced spectral gaps that appear as constant plateaus in the magnetization; these structures are strongly sensitive to edge termination. In both cases the local-marker and direct-definition methods coincide numerically, and the site-resolved marker concentrates on external and internal contours [2510.14556].

The Haldane model also serves as a benchmark for the current-generating local theory. In a topological phase with $t_1=1$, $t_2=0.5$, $\phi=\pi/2$, and $v_{SL}=0.3$, the three-point local magnetization inside a bulk unit cell forms piecewise constant domains whose jumps are fixed by bond currents, while Gaussian smearing of the natural local marker yields edge magnetization and circulating edge currents consistent with the exact bond current [2606.29330].

## 6. First-principles and correlated-material implementations

In crystalline solids, a practical first-principles route uses maximally localized Wannier functions to interpolate the modern-theory formula for orbital magnetization on very fine $k$ meshes. The gauge-invariant formulation expresses the orbital magnetization in terms of traces over small Wannier-space matrices, typically requiring only on the order of $20$ Wannier functions. Real-space decomposition can then be carried down to individual Wannier functions or unit cells, yielding a local magnetization density
$$
m_{\mathrm{loc}}(\mathbf r)=
-\frac e2\,\mathrm{Im}\sum_n
\psi_n^*(\mathbf r)\,
\bigl[(\mathbf r-\mathbf R_n)\times \hat H(\mathbf r-\mathbf R_n)\bigr]
\,\psi_n(\mathbf r),
$$
with
$$
M=\int_{\mathrm{cell}} d^3r\,m_{\mathrm{loc}}(\mathbf r).
$$
In practice, coarse meshes such as $8^3$–$10^3$ are adequate for building the Wannier representation, while the anomalous-Hall-type terms in the itinerant contribution often require fine meshes of order $100^3$ for sub-percent accuracy [1112.1938].

Applications to magnetic oxides show that the separation between local and itinerant pieces can be numerically small in the total moment even when the momentum-resolved kernel is highly structured. In insulating perovskite transition-metal oxides, the modern theory splits the orbital magnetization into a local circulation part and an itinerant circulation part, with the site-diagonal contribution
$$
\mu^0=-\mu_B\,\mathrm{Tr}\{L\,D\}.
$$
For YTiO$_3$, LaMnO$_3$, and orthorhombic and monoclinic YVO$_3$, the itinerant correction $\Delta\mathcal M$ is $10^2$–$10^4$ times smaller than $M^0$, so that $\mu^0$ captures more than $99\%$ of the net ferromagnetic orbital moment. At the same time, $\Delta\mathcal M(\mathbf k)$ follows the behavior of the Chern-invariant kernel and is suppressed only after strong Brillouin-zone cancellation [1311.4263].

Correlated materials can show the opposite trend, with dynamical correlations greatly enhancing the local-circulation channel. In layered ferromagnetic VI$_3$, the band-theory decomposition
$$
M=M_{\mathrm{LC}}+M_{\mathrm{IC}}
$$
was generalized through a DMFT Green’s-function expression. The reported values per V are

| Method | $M_{\mathrm{LC}}$ | $M_{\mathrm{IC}}$ |
|---|---:|---:|
| DFT+U | $-0.012$ | $+0.033$ |
| DMFT $(T=290\,\mathrm K)$ | $+0.208$ | $-0.148$ |
| DMFT $(T=29\,\mathrm K)$ | $+0.208$ | $-0.128$ |

The total orbital magnetization rises from $+0.021\,\mu_B$ in DFT+U to $+0.079\,\mu_B$ and $+0.080\,\mu_B$ in DMFT. The interpretation given is that dynamical correlations localize the V-$d$ electrons more strongly, boost the atomic-like circulating currents, and suppress or reverse the itinerant intercell contribution [2012.13151].

A further theoretical extension emphasizes boundary anomalies. By defining an extended velocity
$$
\hat v_{\mathrm{ext}}=\hat v+\hat v_b,
\qquad
\hat v=i[H,r],
\qquad
\hat v_b=i(H^\dagger-H)\,r,
$$
a non-Hermitian reformulation attributes additional boundary contributions to orbital magnetization, predicts an emergent covariant derivative in the one-band approximation, and identifies a many-band boundary term that can become locally giant near band crossings in the presence of Hall voltage or surface charge imbalance [1811.03090]. A plausible implication is that the treatment of boundaries remains an active fault line between bulk-only, whole-sample, and explicitly current-based local theories.

Local orbital magnetization is therefore best regarded as a family of rigorously defined real-space objects, each tied to a specific physical criterion: nearsighted bulk averaging, site-resolved Berry and geometric texture, exact current generation, or interacting-material response. Their shared macroscopic limit is well established, but their local nonequivalence has become an important part of the subject rather than a defect in it [1301.6275][2512.12343][2606.29330].

Source: https://www.emergentmind.com/topics/local-orbital-magnetization