🔬 materials science

Real-space first-principles approach to orbitronic phenomena in metallic multilayers

This paper presents a scalable, real-space first-principles method based on RS-LMTO-ASA and Chebyshev polynomial expansion to simulate orbital and spin transport phenomena in complex metallic multilayers, successfully demonstrating substantial orbital accumulation in centrosymmetric systems while naturally accounting for disorder and interface effects.

Original authors: Ramon Cardias, Hugo U. R. Strand, Anders Bergman, A. B. Klautau, Tatiana G. Rappoport

Published 2026-06-19

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Ramon Cardias, Hugo U. R. Strand, Anders Bergman, A. B. Klautau, Tatiana G. Rappoport

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you have a massive, complex city made of atoms. In this city, electrons are the citizens moving around. Usually, scientists study how these citizens move by looking at the city from a distance, using a map that assumes the city repeats itself perfectly forever (like a wallpaper pattern). But real materials, like the layers of metal in a computer chip, aren't perfect wallpapers. They have edges, rough surfaces, and different neighborhoods (layers) that don't match up.

This paper introduces a new way to study these electrons: a "real-space" method. Instead of looking at the perfect map, the researchers built a tool that lets them walk through the city street-by-street, atom-by-atom, to see exactly what's happening at the edges and interfaces.

Here is a breakdown of their work using simple analogies:

1. The New Tool: A "Street-Level" Microscope

The researchers developed a computer method called RS-LMTO-ASA. Think of this as a high-tech, street-level microscope.

The Old Way: Previous methods were like trying to understand a city by only looking at a satellite photo of a perfect, repeating grid. If you wanted to see what happened at a specific messy intersection (an interface between two metals), you had to guess or build a simplified model.
The New Way: This new method looks directly at the atoms in their actual positions. It doesn't care if the city is messy, if there are potholes (disorder), or if the buildings are different sizes (roughness). It calculates the electron behavior right where the atoms are.

2. The "Orbitronic" Phenomenon: Spinning Top vs. Dancing Top

The paper focuses on orbitronics. To understand this, imagine an electron as a dancer.

Spin: The dancer spinning in place on one foot. (This is what "spintronics" studies).
Orbit: The dancer running in a circle around the stage. (This is "orbitronics").

The researchers wanted to see what happens when you push these dancers with an electric field (like a wind blowing across the stage).

The Hall Effect: When you push the dancers, they don't just move forward; they drift sideways.
The Orbital Hall Effect: The researchers found that the "dancing" (orbital motion) creates a massive sideways flow, sometimes even bigger than the "spinning" (spin) flow.

3. The Experiment: Testing the City Blocks

The team tested their new microscope on different "neighborhoods" (metals like Copper, Titanium, Platinum, Tungsten, and Nickel).

The Bulk Test: First, they looked at a solid block of metal. They confirmed that their new tool sees the same things other scientists have seen with older tools: heavy metals like Platinum and Tungsten are great at creating these sideways "orbital currents."
The Layer Test: Then, they looked at sandwiches of metals (multilayers). They stacked a magnetic metal (like Nickel) on top of a non-magnetic one.

4. The Big Surprise: The Edge is Different from the Middle

Here is the most interesting finding, explained with a crowd analogy:
Imagine a stadium full of people (electrons).

In the middle of the stadium (the bulk): The crowd moves in a predictable pattern.
At the very edge (the surface): The crowd behaves differently because they can't move past the wall.

The researchers found that in some metals (specifically Tungsten), the "sideways flow" at the edge of the metal was actually pointing in the opposite direction compared to the flow in the middle of the metal.

Why? It's because the "seats" available for the electrons at the edge are different from the seats in the middle. The edge has a different "density of people" (density of states) right where the action happens. This proves you can't just guess what happens at the edge by looking at the middle; you have to look at the edge directly.

5. Magnetic Neighbors

When they put a magnetic metal (Nickel) next to a non-magnetic one, the magnetic "neighbor" changed the behavior of the electrons. It created new types of flows that wouldn't exist in a non-magnetic world. Their tool successfully mapped out exactly how these flows built up layer-by-layer at the interface.

Summary

The paper doesn't claim to have built a new device yet. Instead, it claims to have built a better, more flexible ruler for measuring how electrons move in complex, messy, real-world materials.

What they did: Created a computer method that calculates electron movement in real space, handling disorder and interfaces naturally.
What they found: Confirmed that "orbital" currents are huge in metals, and that the behavior at the surface of a metal can be completely different (even opposite) from the behavior in the center.
Why it matters: It gives scientists a way to simulate complex, real-life metal layers without having to simplify them, which is crucial for designing future electronic devices that rely on these subtle electron movements.

Technical Summary: Real-space first-principles approach to orbitronic phenomena in metallic multilayers

Problem Statement
The field of orbitronics seeks to exploit the orbital degree of freedom of electrons for information processing and spintronic applications. While significant progress has been made in understanding orbitronic effects in bulk quantum materials and 2D systems, transition metals (TMs) remain a central platform for both fundamental research and technological applications due to their complex $d$ -band structures and strong hybridization. However, standard theoretical approaches often rely on momentum-space density functional theory (DFT) calculations followed by the construction of tight-binding Hamiltonians via Wannier functions or pseudo-atomic orbital projections. These methods face challenges when applied to complex heterostructures, disordered systems, or interfaces where translational symmetry is broken. Specifically, there is a need for a scalable, first-principles framework that can directly compute orbital (spin) Hall transport and accumulation in real space, naturally accommodating finite-size effects, interface roughness, and structural imperfections without intermediate approximations.

Methodology
The authors develop a real-space first-principles method based on the Real-Space Linear Muffin-Tin Orbital method within the Atomic Sphere Approximation (RS-LMTO-ASA). This approach retains the full real-space DFT Hamiltonian, avoiding the need for intermediate projection schemes. The methodology consists of two main components:

Electronic Structure Calculation: The RS-LMTO-ASA method is used to obtain self-consistent electronic structures and real-space Hamiltonians. To solve the eigenvalue problem and calculate the Local Density of States (LDOS) efficiently, the authors employ a Chebyshev polynomial-based recursion scheme. This replaces the traditional Lanczos algorithm, offering superior numerical stability, particularly for systems with substantial $s$ - and $p$ -orbital contributions. The LDOS is computed via a Chebyshev expansion of the Green's functions, utilizing the Kernel Polynomial Method (KPM) with a Jackson kernel to ensure convergence and smooth spectral functions.
Linear Response Transport: Transport quantities, such as orbital and spin Hall conductivities and non-equilibrium accumulations, are computed directly in real space using the Kubo-Bastin formalism. The authors employ a Chebyshev polynomial expansion of the Green's functions to evaluate the Kubo-Bastin formula. This allows for the systematic calculation of linear response coefficients (conductivities and magnetoelectric susceptibilities) without relying on stochastic trace evaluations. Instead, observables are computed locally at representative sites, exploiting translational symmetry in periodic directions while fully incorporating contributions from all other sites via Chebyshev recursion.

Key Contributions

Development of a Scalable Framework: The paper presents a methodology that scales linearly with the number of nonequivalent atoms, enabling the simulation of large-scale heterostructures and multilayers with complex geometries.
Direct Real-Space Transport: By integrating the Chebyshev expansion with the RS-LMTO-ASA Hamiltonian, the method computes transport properties directly in real space, eliminating the need for Wannier projections and allowing for the natural inclusion of open boundary conditions and disorder.
Layer-Resolved Analysis: The approach provides atomic-resolution access to the spatial profiles of spin and orbital accumulation across interfaces and within finite slabs, distinguishing between bulk and surface contributions.

Results
The method was applied to a series of transition-metal-based heterostructures and elemental systems:

Bulk Transition Metals: Calculations for fcc Pt, fcc Ti, fcc Cu, fcc Ni, and bcc W reproduced established trends in orbital Hall conductivity (OHC) and spin Hall conductivity (SHC). The results confirmed that $d$ -orbitals dominate the response near the Fermi level in heavy transition metals like Pt and W, while showing good agreement with previous momentum-space studies (e.g., Salemi et al., Go et al.).
Ferromagnetic Ni: In ferromagnetic Ni, the method successfully captured time-reversal symmetry breaking, revealing finite values for both time-reversal even and odd components of the conductivity tensor (e.g., $\sigma^S_{xz}$ and $\sigma^S_{zx}$ ), consistent with theoretical predictions of anomalous spin and orbital Hall effects.
Multilayer Accumulations: In 18-layer slabs of Cu, Ti, Pt, W, and Ni, the study demonstrated that substantial orbital and spin accumulation occurs near surfaces, decaying toward the slab center. Notably, in Tungsten (W), the sign of the surface orbital accumulation differed from the bulk OHC sign. This discrepancy was attributed to significant differences between the surface and bulk density of states (DOS) near the Fermi level, highlighting the critical role of surface electronic structure.
Heterostructures: Simulations of bilayers consisting of two ferromagnetic layers (Fe, Co, Ni) adjacent to nine layers of a non-magnetic TM (Ti or W) revealed that the largest accumulations occur near the interface. Fe-based heterostructures exhibited the strongest orbital accumulation, particularly in Fe/Ti systems. The study also noted that while accumulation is a proxy for torque, the actual torque is governed by the divergence of the angular momentum flux, a distinction the authors note requires future flux calculations.

Significance
The paper claims that the developed real-space methodology provides a scalable and flexible framework for realistic simulations of orbital transport phenomena. Its primary significance lies in its ability to model systems where spatial inhomogeneity and finite-size effects are crucial, such as interfaces and disordered multilayers, which are difficult to treat with standard momentum-space approaches. The authors demonstrate that substantial orbital (spin) accumulation can arise even in centrosymmetric systems, driven by band-structure asymmetries and interfacial scattering. Furthermore, the work establishes that the relationship between bulk conductivities and surface accumulations is non-trivial and heavily influenced by surface electronic structure. The methodology is presented as a powerful tool for investigating orbitronic phenomena in complex heterostructures, offering direct access to local responses and paving the way for studying systems with reduced symmetry and disorder relevant to next-generation spintronic and orbitronic devices.