Microscopic origin of Rashba coupling from first principles: Layer-resolved orbital asymmetry in transition metal dichalcogenides

This paper provides a first-principles microscopic explanation for Rashba spin splitting in transition metal dichalcogenides, revealing that the effect arises from a competition between internal polarization and interlayer hybridization and can be characterized by a new orbital polarization imbalance order parameter.

The Gist

The Tale of the Unbalanced Seesaw: Understanding Spin in 2D Materials

Imagine you are looking at a perfectly balanced seesaw in a playground. If everything is symmetrical—the board is flat, the weight is centered, and there’s no wind—the seesaw stays perfectly level. In the world of physics, this is like a material with perfect symmetry: the electrons (the tiny particles that carry electricity and "spin") behave in a very predictable, balanced way.

But what happens if you tilt the seesaw? Or what if one side of the playground is built on a hill? Suddenly, everything changes.

This scientific paper explores a phenomenon called Rashba coupling in ultra-thin materials called Transition Metal Dichalcogenides (TMDs). To understand their discovery, let’s use three simple metaphors.

---

1. The "Spin" of the Electron: The Tiny Spinning Top

Every electron acts like a tiny, microscopic spinning top. This "spin" is a fundamental property. In a perfectly symmetrical material, these tiny tops spin in a way that is balanced and uniform.

However, when you break the symmetry of the material (by stacking layers unevenly or applying an electric field), you create a "slope." This slope forces the electrons to spin in specific directions depending on how they are moving. This "forced spin direction" is what scientists call Rashba coupling. It is the "holy grail" for spintronics—a future technology where we use the spin of an electron, rather than just its charge, to make computers much faster and more energy-efficient.

2. The "Orbital Imbalance": The Uneven Room

The researchers wanted to know why this happens at a microscopic level. They discovered it isn't just about the material being "tilted"; it’s about how the "furniture" (the electron clouds, or orbitals) is arranged in the room.

Think of an atom like a room. In a perfectly symmetrical material, the furniture is spread evenly across the floor. But in these TMD materials, the researchers found an "Orbital Polarization Imbalance."

Imagine a room where all the heavy chairs are pushed toward the left wall. Even if the floor is level, the distribution of weight is lopsided. This "lopsidedness" in the electron clouds is what actually triggers the Rashba effect. The researchers even created a mathematical "scorecard" (an order parameter) to measure exactly how lopsided these electron clouds are.

3. The "Competition": The Tug-of-War

One of the most surprising findings in the paper is that when you stack two layers of these materials together (a bilayer), the Rashba effect actually gets weaker in some cases, even though you’d expect more layers to mean more effect.

Why? The researchers describe this as a microscopic tug-of-war.

• On one side, you have the Atomic Strength: The heavy atoms want to create a strong spin effect.
• On the other side, you have Polarizability: How easily the electron clouds can be "pushed" or "stretched" by an electric field.

In certain materials, as the atoms get heavier, they get "stiffer" and harder to move. It’s like trying to tilt a heavy lead seesaw versus a light wooden one. Even though the lead seesaw is more powerful, it’s so heavy and stubborn that it doesn't respond to your touch as easily. This competition determines whether the material will be a good "highway" for spin-based technology.

---

Why does this matter?

Right now, our computers get hot and waste energy because we are just pushing "charge" around. If we can master the Rashba effect, we can control the spin of electrons with much less effort.

By understanding the "lopsidedness" of these tiny electron clouds, these scientists have provided a "blueprint." They’ve shown us exactly how to stack and tune these 2D materials to create the next generation of ultra-fast, cool-running electronic devices.

Technical Summary

Technical Summary: Microscopic Origin of Rashba Coupling in Transition Metal Dichalcogenides (TMDs)

1. Problem Statement

While it is well-established that breaking inversion ($I$) and mirror ($m_z$) symmetries in two-dimensional (2D) materials induces Rashba spin splitting, the microscopic origin and quantitative characterization of this effect in transition metal dichalcogenides (TMDs) have remained incomplete. Specifically, there has been a lack of understanding regarding:

• Why certain bilayer (BL) configurations exhibit significantly smaller Rashba splitting than monolayers (ML) despite having built-in electric dipoles.
• The fundamental competition between atomic spin-orbit coupling (SOC) and orbital hybridization that dictates the magnitude of the Rashba coefficient.
• How to bridge the gap between macroscopic observables (energy-momentum extrema) and atomic-scale orbital asymmetry.

2. Methodology

The authors employ a multi-scale theoretical approach to resolve the Rashba effect from first principles:

• Ab Initio Calculations: Density Functional Theory (DFT) using the Generalized Gradient Approximation (PBE) and DFT-D3-BJ for van der Waals interactions. They utilized fully-relativistic norm-conserving pseudopotentials within the Quantum ESPRESSO package.
• Wannierization: The DFT Kohn-Sham states were transformed into a Maximally Localized Wannier Function (MLWF) basis using wannier90. This allowed for an orbital-resolved Hamiltonian focusing on transition metal $d$-orbitals ($d_{z^2}, d_{xz}, d_{yz}$) and chalcogen $p$-orbitals ($p_x, p_y, p_z$).
• Perturbative Microscopic Model: A tight-binding framework was developed to treat atomic SOC and orbital hybridization perturbatively. This model specifically examines the mixing of $m_l = 0$ states (like $p_z, d_{z^2}$) with $m_l = \pm 1$ states (like $p_{x,y}, d_{xz,yz}$) mediated by momentum-dependent hopping.
• Order Parameter Definition: They introduced the orbital polarization imbalance ($\xi_n$) and orbital polarizability ($\chi_n^{orb}$) to quantify uniaxial charge asymmetry and its response to external fields.

3. Key Contributions

• New Descriptors: The introduction of $\xi_n$ (orbital polarization imbalance) as a central order parameter that captures the asymmetry of valence bands and determines the spin ordering of Rashba-split states.
• Decomposition of the Rashba Effect: They decoupled the Rashba response into two components: the Rashba coefficient ($\lambda_R^n$), representing the response to an external field, and the intrinsic orbital field ($E_n^0$), representing the built-in asymmetry.
• First Microscopic Model for $\Gamma$-point Rashba in TMDs: A theoretical derivation showing that Rashba splitting arises from two primary pathways:
  1. $|p_{zi}\rangle \xrightarrow{\text{hopping}} |d_{xz}\rangle \xrightarrow{\text{SOC}} |d_{z^2}\rangle$
  2. $|d_{z^2}\rangle \xrightarrow{\text{hopping}} |p_{xi}\rangle \xrightarrow{\text{SOC}} |p_{zi}\rangle$
• Identification of Mediator Orbitals: The study identifies the in-plane $d_{xz/yz}$ and $p_{x/y}$ orbitals as the essential mediators for the Rashba process.

4. Key Results

• Bilayer vs. Monolayer Paradox: The study explains why the Anti-bonding Valence Band (AVB) in bilayers shows reduced Rashba splitting compared to monolayers. This is due to a competition between internal polarization (which promotes splitting) and interlayer hybridization (which can suppress it).
• Material Trends:
  • Increasing the chalcogen mass (e.g., S $\to$ Se $\to$ Te) monotonically increases the Rashba coefficient due to increased atomic SOC and higher orbital polarizability.
  • The influence of the transition metal (Mo vs. W) is surprisingly weak because the increase in atomic SOC from W is offset by a decrease in orbital polarizability.
• Spin Ordering: The sign of the spin splitting (and thus the spin texture) is directly governed by the sign of the orbital polarization imbalance $\xi_n$.
• Orbital Polarizability: The BVB (Bonding Valence Band) exhibits a much larger Rashba response than the AVB because it possesses higher in-plane $d$-orbital weight, allowing for more efficient hybridization.

5. Significance

This work provides a unified, quantitative, and qualitative framework for understanding Rashba physics in TMDs. By moving beyond simple energy-momentum mapping to an orbital-resolved description, the authors enable:

• Predictive Engineering: The ability to tailor spin-orbit effects by manipulating stacking configurations (e.g., $R$-type vs. $H$-type) and chalcogen selection.
• Spintronic Design: Insights into how external gating can control spin textures, which is critical for low-power, non-volatile spintronic devices.
• Generalizability: The proposed descriptors ($\xi_n$ and $\chi_n^{orb}$) are applicable to more complex 2D systems, such as Janus monolayers and heterobilayers, providing a roadmap for the design of next-generation van der Waals heterostructures.