Modern Approach to Orbital Hall Effect Based on Wannier Picture of Solids

This paper introduces a new Wannier-function-based approach grounded in the modern theory of orbital magnetization to accurately quantify orbital angular momentum by capturing both local and itinerant contributions, thereby revealing significant non-local corrections to the orbital Hall conductivity that are missed by traditional atom-centered approximations.

The Gist

The Big Picture: The "Orbital" Revolution

Imagine you have a tiny, spinning top. In the world of atoms, electrons are like those tops. They have two main ways of moving:

1. Spin: Spinning on their own axis (like a top spinning in place).
2. Orbit: Moving around the nucleus (like the Earth orbiting the Sun).

For decades, scientists focused almost entirely on the Spin. It's the "classic" way we make hard drives and magnetic memory work. But recently, scientists realized that the Orbit is actually a much more powerful, untapped resource. If you can control the "orbit" of electrons, you can create electricity and magnetic effects that are orders of magnitude stronger than what we get from spin alone. This new field is called Orbitronics.

The Problem: The "Blind Map"

To use this power, scientists need to measure something called the Orbital Hall Conductivity (OHC). Think of this as a map that tells you how much "orbital traffic" flows through a material when you apply electricity.

However, measuring this map has been incredibly difficult.

• The Old Way (The "Atom-Centered" Approximation): Imagine trying to understand the traffic in a massive city by only looking at the front door of every single house. You see what's happening right at the doorstep (local effects), but you miss the highways, the bridges, and the traffic flowing between neighborhoods.
• The Flaw: This old method assumes electrons only care about the atom they are currently sitting on. It ignores the fact that electrons are "itinerant"—they roam freely and interact with their neighbors. Because of this, the old maps were often wrong, missing huge amounts of the actual traffic.

The Solution: The "Wannier" GPS

The authors of this paper introduced a new, high-tech GPS system based on something called Wannier functions.

• The Analogy: Instead of just looking at the front door of a house, imagine you have a drone that can see the entire neighborhood, the roads connecting the houses, and even the traffic flowing between them.
• How it works: They used a mathematical framework (based on the "Modern Theory of Orbital Magnetization") to create a complete picture. This new approach accounts for two types of traffic:
  1. Local Circulation: The electron spinning in place (the front door traffic).
  2. Itinerant Circulation: The electron moving around the whole crystal (the highway traffic).

The Discovery: The "Hidden" Traffic

When the team applied this new GPS to real materials (like Platinum, Iron, and Tungsten), they found something shocking:

1. The Old Maps Were Wrong: In many materials, the "local" traffic (what the old method saw) was actually being canceled out by the "itinerant" traffic (the new stuff they found).
2. The Sign Flip: In some cases, like Platinum, the old method said the traffic was flowing one way (positive), but the new method showed it was actually flowing the opposite way (negative) and was much stronger.
3. The "J-Matrix" Effect: They discovered that the "highway" part of the traffic is incredibly sensitive to the specific shape of the material's electronic structure. It's not just a generic background noise; it's a complex, dynamic flow that changes depending on the material's details.

Why This Matters for the Future

Think of this like the difference between building a car based on a sketch of a wheel versus building it based on a full 3D simulation of the engine, the transmission, and the road.

• Better Devices: If we want to build next-generation computers or sensors that use "orbital currents" instead of just spin, we need accurate maps. If we use the old, incomplete maps, we might build a device that doesn't work or works in the opposite way we expect.
• Designing Materials: This new method allows scientists to predict exactly how a material will behave before they even build it. They can now say, "If we tweak the crystal structure here, we can flip the direction of the current," which is crucial for designing efficient electronics.

Summary in One Sentence

The paper introduces a new, super-accurate mathematical tool that stops scientists from "looking at the trees and missing the forest" when studying electron orbits, revealing that the "roaming" electrons contribute just as much (or more) to magnetic effects as the "stationary" ones, completely changing how we predict and design future electronic devices.

Technical Summary

1. Problem Statement

The field of orbital dynamics and orbitronics has gained significant traction due to the discovery that orbital angular momentum (OAM) can generate currents orders of magnitude larger than spin currents in light materials. However, a rigorous theoretical framework for quantifying the Orbital Hall Conductivity (OHC) has been missing.

• The Core Challenge: In solid-state theory, calculating the OAM operator using Bloch states is mathematically difficult because the position operator ($\hat{r}$) is unbounded.
• Limitations of Current Methods: The standard approach, the Atom-Centered Approximation (ACA), assumes OAM is localized to atomic sites. While efficient, ACA:
  • Captures only local contributions.
  • Fails to account for itinerant (non-local) contributions arising from surface effects or the circulation of Wannier functions across the unit cell.
  • Deviates from the rigorous modern theory of orbital magnetization (OM).
  • Cannot distinguish between local and itinerant contributions to orbital dynamics.
• Gap: Previous attempts to quantify non-local OAM failed to correctly incorporate surface contributions or lacked gauge invariance, leading to inaccurate predictions of the OHC.

2. Methodology

The authors propose a new framework based on the Wannier function picture of solids, extending the modern theory of orbital magnetization to derive a complete, gauge-covariant OAM operator.

• Wannier Basis: Instead of Bloch states, the method utilizes maximally localized Wannier functions, which are well-defined in real space and handle the unbounded position operator naturally.
• Extension of Modern OM Theory:
  • The authors separate OAM into Local-Circulation (LC) (self-rotation of Wannier functions) and Itinerant-Circulation (IC) (circulation of surface Wannier functions).
  • They derive the matrix elements of the OAM operator ($L^H_{nm}$) in the Hamiltonian gauge by carefully treating the position operator in bulk and surface regions.
• Gauge Transformation and Covariant Derivatives:
  • The formalism involves switching between the Wannier gauge and the Hamiltonian gauge.
  • To ensure numerical stability and gauge covariance (independence of the specific gauge choice), the authors introduce a covariant derivative ($D_k$). This replaces standard $k$-derivatives in the OAM expression.
  • This approach eliminates numerical instabilities (poles) that typically arise at band crossings when using Berry-phase quantities directly.
• Dynamic Gauge for OHC:
  • To calculate the Orbital Hall Conductivity, a dynamic gauge is introduced to account for time-dependent perturbations (electric fields).
  • Crucially, the authors identify an additional term ($\delta \tilde{L}_\alpha$) in the orbital current change. This term arises from the non-locality of the modern OAM operator and was missing in previous ACA-based calculations.

3. Key Contributions

1. Derivation of the Modern OAM Operator: A rigorous expression for the full OAM operator in the Wannier basis that includes both local and itinerant contributions and is fully gauge-covariant.
2. Identification of Non-Local Corrections: The demonstration that the OAM operator contains non-local terms (involving the $J$-matrix and covariant derivatives) that significantly alter the calculated OHC compared to the ACA.
3. Separation of LC and IC Contributions: A method to explicitly decompose the OHC into Local-Circulation and Itinerant-Circulation components, revealing their distinct behaviors and signs.
4. First-Principles Implementation: The development of a computational workflow using ORBITRANS (based on Wannier90 and FLEUR) to calculate OHC in complex materials without empirical fitting.

4. Results

The authors applied their formalism to various transition metals (Pt, W, Fe, V, Cr, Cu, Ti) and MoS$_2$. Key findings include:

• Discrepancy with ACA: In FCC Platinum (Pt), the total OHC calculated via the modern approach differs drastically from the ACA prediction.
  • Sign Reversal: The LC and IC contributions are of similar magnitude but opposite signs. In Pt, the IC contribution dominates, resulting in a total OHC with the opposite sign compared to the ACA result.
  • Band Filling Dependence: The non-local contributions show a highly non-trivial dependence on band filling, contradicting the common wisdom that OHE is a robust, material-independent property.
• Material-Specific Trends:
  • In many materials (e.g., bcc Fe, hcp Ti), the itinerant (IC) contribution is larger than the local (LC) contribution and often opposes it.
  • The $O(J^2)$ term (part of the non-local expansion) dominates the OHC in bcc Fe but vanishes in monolayer MoS$_2$.
• Local vs. Non-Local: When terms containing the $J$-matrix (non-local) are removed, the modern OHC closely matches the ACA result, confirming that ACA captures local effects but misses critical non-local physics.
• K-Space Sensitivity: The contributions to OHC vary significantly across the Brillouin zone. Different regions of $k$-space contribute differently to LC and IC, suggesting that experimental probes (like torque measurements) might be sensitive to specific parts of the Fermi surface.

5. Significance

• Theoretical Advancement: This work resolves the long-standing difficulty of defining the OAM operator in solids by providing a mathematically rigorous, gauge-invariant, and numerically stable framework.
• Re-evaluation of Orbitronics: The findings suggest that previous estimates of the orbital Hall effect (based on ACA) may be qualitatively incorrect regarding both magnitude and sign in many materials. This necessitates a re-evaluation of experimental data interpreted through the lens of the OHE.
• Device Design Implications: The strong dependence of OHC on the itinerant contribution and specific $k$-space regions implies that crystallographic engineering and interface control can be used to selectively enhance or suppress specific orbital currents. This is crucial for designing next-generation orbitronic devices and understanding orbital torque mechanisms in magnetic heterostructures.
• Experimental Guidance: The paper highlights that different experimental techniques (e.g., XMCD vs. torque measurements) may probe different components (local vs. itinerant) of the OAM, explaining potential discrepancies between theory and experiment.

In summary, the paper establishes that non-local effects are not negligible corrections but fundamental components of the orbital Hall effect, fundamentally changing our understanding of orbital transport in solids.