🔬 materials science

The Critical Role of Itinerant Contributions to Orbital Angular Momentum Relaxation and Dynamics

This study resolves the discrepancy between theoretical predictions and experimental observations of orbital angular momentum (OAM) relaxation by demonstrating, through first-principles simulations of MoS2 and silicene, that itinerant OAM contributions are essential for accurately describing long diffusion lengths, as atom-centered models fail to capture these effects and significantly underestimate OAM lifetimes.

Original authors: Andrew C. Grieder, Luis M. Canonico, Frederico Simões, Aron W. Cummings, Yuan Ping

Published 2026-09-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Andrew C. Grieder, Luis M. Canonico, Frederico Simões, Aron W. Cummings, Yuan Ping

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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

In the quest to build faster, more efficient electronics, scientists are looking beyond the simple flow of electric charge. For decades, the industry has relied on manipulating the charge of electrons, but this approach generates heat and wastes energy. A newer field, known as orbitronics, proposes using a different property of the electron: its orbital angular momentum. Imagine an electron not just as a tiny particle moving through a wire, but as a spinning top that also orbits around the nucleus of an atom. This orbital motion carries a specific type of momentum that can be used to store information or control magnetic fields with far less energy loss than traditional methods. However, for this technology to work, scientists must understand how long this orbital motion lasts before it fades away. If the momentum disappears too quickly, the information is lost. Recent experiments have shown that electrons can maintain this orbital state for surprisingly long distances, yet older computer models predicted they would vanish almost instantly. This contradiction has left researchers unsure of what actually controls the lifespan of these orbital states.

A team of researchers has now resolved this mystery by developing a new way to simulate how electrons behave in solid materials. They focused on two specific materials: molybdenum disulfide, a compound with strong internal magnetic interactions, and silicene, a form of silicon that behaves very differently. Using advanced computer simulations that track the movement of electrons and their interactions with the vibrating atoms of the material, the team discovered that the answer lies in how we define the orbital momentum itself. For years, scientists have used a simplified model that treats the orbital motion as if it were locked strictly to the individual atoms, like a planet orbiting a single star. The researchers found that this "atom-centered" view is incomplete. It misses a crucial component where the orbital motion is tied to the overall shape and structure of the electron's path through the entire crystal. This "itinerant" contribution, which moves freely with the electron, turns out to be the dominant factor in how long the orbital state survives.

When the team applied their new framework to molybdenum disulfide, they observed a complex, two-step process for how the orbital momentum relaxes. Initially, the electrons, which were excited by a pulse of light, rapidly redistributed their energy between different regions of the material's internal structure. This first stage happened in a fraction of a picosecond, causing the orbital signal to drop quickly. However, this initial drop was not the end of the story. The remaining orbital momentum then decayed much more slowly, over a period of about 6.1 picoseconds, in a way that was tightly linked to the electron's spin. The simulations revealed that the old, atom-centered models predicted the orbital state would vanish in less than three femtoseconds, a time so short it is effectively instantaneous. This massive discrepancy occurred because the old models only saw the fast, chaotic wobbling of the electron around its specific atom, driven by the electric fields of the crystal. They completely missed the smoother, longer-lasting motion of the electron as it traveled through the material.

In the case of silicene, where the internal magnetic interactions are weak, the researchers found that the spin and orbital motions could be controlled independently. By applying an electric field, they could speed up or slow down the relaxation of the electron's spin, yet the orbital lifetime remained stubbornly unchanged. This confirmed that the two properties are governed by different rules. The simulations showed that the atom-centered component of the orbital motion still wobbled and decayed rapidly, driven by the crystal's internal structure. In contrast, the itinerant component, which represents the electron's journey through the material, persisted for tens of picoseconds. The key finding was that the itinerant orbital momentum is so much more stable that it dictates the overall behavior of the system. In these materials, the orbital lifetime is determined not by how the electron spins around a single atom, but by how it moves between different valleys in the material's energy landscape.

The researchers concluded that the wide range of experimental results seen in the past, where some measurements showed long orbital distances and others showed very short ones, likely stems from which part of the orbital motion is being measured. The atom-centered models, which have been the standard for theoretical predictions, consistently underestimate the lifespan of orbital states by at least a factor of ten. By ignoring the itinerant contribution, these models miss the very mechanism that allows orbital angular momentum to survive long enough to be useful. The study demonstrates that to accurately describe how electrons move and relax in modern materials, scientists must look beyond the individual atoms and consider the electron's journey through the entire crystal structure. This shift in perspective is essential for designing the next generation of low-energy electronic devices, ensuring that the theoretical tools used to predict their behavior match the reality of the materials they are built from.

Technical Summary: The Critical Role of Itinerant Contributions to Orbital Angular Momentum Relaxation and Dynamics

Problem Statement
Orbital angular momentum (OAM) is recognized as a promising degree of freedom for low-dissipation transport and magnetization control. However, a significant discrepancy exists between experimental measurements and theoretical predictions regarding OAM relaxation mechanisms. Experimentally, orbital diffusion lengths vary widely, ranging from ~70–80 nm in some systems to merely 2–9 nm in others. Theoretically, first-principles quantum transport simulations based on the atom-centered approximation (ACA) predict orbital relaxation lengths of less than 1 nm, whereas nonlocal transport calculations in multi-orbital lattices suggest responses persisting over hundreds or thousands of lattice constants.

The authors posit that this conflict arises not from material differences or methodological flaws, but from the underlying definition of OAM. Existing ACA-based approaches treat OAM as a purely local quantity projected onto atomic orbitals, neglecting itinerant contributions associated with the geometric structure of Bloch states (Berry-phase contributions). Conversely, some transport simulations implicitly capture these itinerant effects by inferring transport from orbital character rather than direct OAM evaluation. The paper aims to resolve this by developing a framework that treats local and itinerant OAM on equal footing while incorporating quantum scattering and relaxation dynamics.

Methodology
The study employs two complementary first-principles approaches to disentangle spin and orbital relaxation channels in monolayer MoS $_2$ (strong spin-orbit coupling, SOC) and silicene (weak SOC):

First-Principles Lindbladian Density-Matrix Dynamics (FPDMD):
- This method computes OAM using a Berry-phase formalism derived from $k$ -space derivatives of Bloch states, consistent with the modern theory of polarization. This formulation inherently includes itinerant contributions.
- The dynamics are governed by the Lindbladian equation of motion for the one-particle electron density matrix, explicitly including electron-phonon scattering and spin-orbit coupling.
- Observables (spin and OAM) are calculated as the trace of the operator with the time-evolving density matrix, and lifetimes are extracted by fitting the decay to exponential functions.
Real-Space Tight-Binding (TB) Quantum Transport:
- This approach utilizes a linear-scaling method based on Chebyshev polynomial expansion of spectral functions.
- Crucially, this method allows for the explicit separation of OAM components: the atom-centered component ( $L^{\text{ACA}}$ ) and the full itinerant component including Berry-phase contributions ( $L^{\text{Berry}}$ ).
- The TB Hamiltonians are parameterized using first-principles data (e.g., via PAOFLOW for MoS $_2$ and DFT band structures for silicene). Disorder (Anderson) is introduced to match scattering rates observed in FPDMD.

Key Results

MoS $_2$ (Strong SOC Regime):
- Multi-timescale Dynamics: Orbital relaxation exhibits two distinct timescales. An initial fast decay (0.25 ps) corresponds to intervalley redistribution of valley-polarized populations, where carriers scatter into opposite valleys with opposite orbital but identical spin polarization. A slower decay (6.1 ps) follows, where spin and orbital degrees of freedom are locked together by strong SOC.
- Itinerant Dominance: When comparing components, the ACA-only orbital lifetime is calculated to be < 3 fs, which is three to four orders of magnitude shorter than the full FPDMD result. This indicates that in strong-SOC systems, the relaxation is governed by the itinerant OAM, which the ACA omits.
Silicene (Weak SOC Regime):
- Decoupled Dynamics: Spin and orbital dynamics decouple. While an external electric field dramatically enhances spin relaxation (via field-induced splitting and precession), orbital lifetimes remain largely independent of the field.
- ACA vs. Itinerant Behavior: The ACA component exhibits rapid oscillations (frequency $\omega_{p_z} \approx 0.7 \times \Delta_{p_z}/\hbar$ ) driven by crystal-field splitting, decaying on a sub-picosecond scale. In contrast, the itinerant OAM shows minimal oscillation and decays on a timescale of tens of picoseconds.
- Lifetime Discrepancy: The ACA orbital lifetime is at least one order of magnitude shorter than the itinerant lifetime across all electric fields. The ACA decay is driven by dephasing of coherent crystal-field precession, whereas the itinerant decay is scattering-limited.
Relaxation Mechanisms:
- ACA Scaling: The ACA orbital relaxation rate scales approximately linearly with the momentum relaxation rate ( $\tau_p^{-1}$ ), mimicking the Elliott-Yafet (EY) mechanism but driven by crystal-field precession interruption.
- Itinerant Scaling: The itinerant relaxation is determined by the distribution of orbital texture in momentum space. In silicene, intervalley scattering dominates, connecting pockets of opposite OAM. While the relative contribution of intervalley vs. intravalley scattering shifts with the electric field, the total lifetime remains nearly field-independent.

Significance and Claims
The paper claims that the wide range of predicted and measured orbital diffusion lengths originates from the definition of OAM itself. By demonstrating that ACA-based models yield lifetimes at least an order of magnitude shorter than those including itinerant contributions, the authors conclude that:

Itinerant OAM is the governing quantity: In the systems studied (MoS $_2$ and silicene), the itinerant component, associated with the geometric structure of Bloch states, dictates the actual orbital relaxation and diffusion.
Limitations of Local Models: Atom-centered approximations fail to capture the long-lived nature of orbital dynamics because they miss the itinerant contributions and are dominated by ultrafast crystal-field precession.
General Applicability: The findings suggest that theoretical studies of orbital dynamics and the interpretation of experimental data must include itinerant OAM contributions. This is asserted as essential for 2D materials like MoS $_2$ and silicene, and likely applicable to the broader family of 2D materials.
Material Dependence: In metals, the balance between atomic and itinerant channels is expected to depend on orbital character (localized vs. delocalized) and sample quality (disorder), potentially explaining the variability in experimental diffusion lengths observed in transition metals.

The study establishes a unified microscopic picture where orbital relaxation is not a single mechanism but a competition between local crystal-field effects and itinerant geometric contributions, with the latter often dominating the observable transport properties.

More like this