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

Riassunto Tecnico: Origine microscopica del coupling Rashba nei TMD

1. Il Problema (Problem Statement)

Sebbene sia noto che la rottura delle simmetrie di inversione ($I$) e di specchiatura ($m_z$) nei materiali bidimensionali (2D) generi lo splitting di spin di tipo Rashba, la sua origine microscopica e la sua caratterizzazione quantitativa nei dicalcogenuri dei metalli di transizione (TMD) rimangono parzialmente oscure. In particolare, non era chiaro perché in alcune strutture bilayer (BL) lo splitting di Rashba fosse significativamente inferiore rispetto ai monolayer (ML) sottoposti a campi elettrici simili, né come l'asimmetria orbitale interna influenzasse tale fenomeno.

2. Metodologia (Methodology)

Gli autori utilizzano un approccio combinato che integra calcoli ab initio e modellazione teorica:

• Calcoli DFT (Density Functional Theory): Utilizzando il pacchetto Quantum ESPRESSO con approssimazione GGA-PBE e correzioni di van der Waals (DFT-D3-BJ), sono stati studiati sistemi di MoS₂, MoSe₂, MoTe₂, WS₂, WSe₂ e WTe₂ in fase 1H, sia come monolayer che come bilayer con stacking di tipo $R_M X$ (che rompe sia $I$ che $m_z$).
• Wannierizzazione: Attraverso il pacchetto wannier90, gli stati di Kohn-Sham sono stati proiettati in funzioni di Wannier per ottenere un Hamiltoniano risolto per orbitali (in particolare gli orbitali $d$ del metallo e $p$ del calcogeno).
• Modello Perturbativo: È stato sviluppato un modello di tight-binding che tratta l'accoppiamento spin-orbita (SOC) e l'ibridazione orbitale in modo perturbativo per identificare i processi atomistici responsabili dello splitting.
• Parametri di Ordine: Viene introdotto il parametro di squilibrio della polarizzazione orbitale ($\xi_n$) per quantificare l'asimmetria della distribuzione di carica rispetto al piano del metallo.

3. Contributi Chiave (Key Contributions)

• Definizione di nuovi descrittori: Gli autori introducono il coefficiente di Rashba $\lambda_R^n$ (risposta al campo esterno) e il campo orbitale intrinseco $E_n^0$ (che cattura l'asimmetria built-in).
• Parametro d'ordine $\xi_n$: L'introduzione della polarizzazione orbitale come descrittore per governare l'ordine di spin e la magnitudo dello splitting.
• Modello Microscopico: La formulazione di un modello che spiega lo splitting attraverso due percorsi di ibridazione orbitale mediati dal momento ($k$):
  1. $|p_z\rangle \xrightarrow{\text{hopping}} |d_{xz}\rangle \xrightarrow{\text{SOC}} |d_{z^2}\rangle$
  2. $|d_{z^2}\rangle \xrightarrow{\text{hopping}} |p_x\rangle \xrightarrow{\text{SOC}} |p_z\rangle$
• Identificazione dei mediatori: La dimostrazione che gli orbitali $d_{xz}/d_{yz}$ del metallo e $p_{x}/p_{y}$ del calcogeno sono i mediatori essenziali del processo Rashba.

4. Risultati Principali (Results)

• Competizione tra polarizzazione e ibridazione: Lo splitting di Rashba nei bilayer è ridotto rispetto ai monolayer a causa di una competizione tra la polarizzazione interna e l'ibridazione interstrato.
• Effetto del Calcogeno: La magnitudo di $\lambda_R$ aumenta monotonicamente con la massa del calcogeno (da S a Te). Questo è dovuto alla diminuzione dell'ionicità (maggiore polarizzabilità) e all'aumento del SOC atomico del calcogeno.
• Effetto del Metallo di Transizione: L'influenza del metallo (Mo vs W) è debole poiché l'aumento del SOC atomico del tungsteno (W) è compensato da una riduzione della polarizzabilità orbitale.
• Comportamento dei Band (AVB vs BVB): Nei bilayer, la banda anti-legante (AVB) mostra uno splitting molto più piccolo rispetto alla banda legante (BVB). Ciò è spiegato dalla maggiore densità di orbitali $d$ in-plane nella BVB, che facilita l'ibridazione.
• Correlazione con $\xi_n$: Lo squilibrio della polarizzazione orbitale $\xi_n$ predice con precisione il campo elettrico critico al quale lo splitting di Rashba scompare.

5. Significato e Implicazioni (Significance)

Il lavoro fornisce una base unificata e quantitativa per confrontare l'accoppiamento Rashba in tutta la famiglia dei TMD. La capacità di comprendere come la struttura atomica e lo stacking influenzino lo spin permette un controllo predittivo dello spin tramite ingegneria dei difetti, stacking (es. Janus layers o etero-strutture) e gating elettrico. Questo è fondamentale per lo sviluppo di dispositivi spintronici a basso consumo e per l'integrazione di funzionalità spin-carica-ottica nei materiali 2D.