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

Sintesi Tecnica: Approccio ai primi principi nello spazio reale per i fenomeni orbitronici in multistrati metallici

Problematica
Il campo dell'orbitronica mira a sfruttare il grado di libertà orbitale degli elettroni per l'elaborazione delle informazioni e le applicazioni spintroniche. Sebbene siano stati fatti progressi significativi nella comprensione degli effetti orbitronici in materiali quantistici bulk e sistemi 2D, i metalli di transizione (TM) rimangono una piattaforma centrale sia per la ricerca fondamentale che per le applicazioni tecnologiche, grazie alle loro complesse strutture a banda $d$ e alla forte ibridazione. Tuttavia, gli approoli teorici standard si affidano spesso a calcoli di teoria del funzionale della densità (DFT) nello spazio dei momenti, seguiti dalla costruzione di Hamiltoniane di tight-binding tramite funzioni di Wannier o proiezioni di orbitali pseudo-atomici. Questi metodi affrontano sfide quando applicati a eterostrutture complesse, sistemi disordinati o interfacce in cui la simmetria traslazionale è interrotta. Nello specifico, è necessario un quadro scalabile, basato sui primi principi, che possa calcolare direttamente il trasporto e l'accumulo orbitale (di spin) nello spazio reale, accogliendo naturalmente effetti di dimensione finita, rugosità delle interfacce e imperfezioni strutturali senza approssimazioni intermedie.

Metodologia
Gli autori sviluppano un metodo ai primi principi nello spazio reale basato sul metodo Real-Space Linear Muffin-Tin Orbital nell'approssimazione Atomic Sphere (RS-LMTO-ASA). Questo approccio mantiene l'intera Hamiltoniana DFT nello spazio reale, evitando la necessità di schemi di proiezione intermedi. La metodologia consiste in due componenti principali:

1. Calcolo della struttura elettronica: Il metodo RS-LMTO-ASA viene utilizzato per ottenere strutture elettroniche autoself-consistent e Hamiltoniane nello spazio reale. Per risolvere il problema degli autovalori e calcolare la Densità Locale degli Stati (LDOS) in modo efficiente, gli autori impiegano uno schema di ricorsione basato sui polinomi di Chebyshev, che offre una stabilità numerica superiore rispetto al tradizionale algoritmo di Lanczos, particolarmente per sistemi con sostanziali contributi di orbitali $s$ e $p$. La LDOS è calcolata tramite un'espansione di Chebyshev delle funzioni di Green, utilizzando il Metodo del Polinomio Kernel (KPM) con un kernel di Jackson per garantire la convergenza e funzioni spettrali lisce.
2. Trasporto a risposta lineare: Le quantità di trasporto, come le conducibilità Hall orbitali e di spin e gli accumuli fuori equilibrio, sono calcolate direttamente nello spazio reale utilizzando il formalismo di Kubo-Bastin. Gli autori impiegano un'espansione in polinomi di Chebyshev delle funzioni di Green per valutare la formula di Kubo-Bastin. Ciò consente il calcolo sistematico dei coefficienti di risposta lineare (conducibilità e suscettibilità magnetoelettrica) senza fare affidamento su valutazioni stocastiche della traccia. Invece, gli osservabili sono calcolati localmente in siti rappresentativi, sfruttando la simmetria traslazionale nelle direzioni periodiche pur incorporando pienamente i contributi di tutti gli altri siti tramite la ricorsione di Chebyshev.

Contributi Chiave

• Sviluppo di un quadro scalabile: Il documento presenta una metodologia che scala linearmente con il numero di atomi non equivalenti, consentendo la simulazione di eterostrutture su larga scala e multistrati con geometrie complesse.
• Trasporto diretto nello spazio reale: Integrando l'espansione di Chebyshev con l'Hamiltoniana RS-LMTO-ASA, il metodo calcola le proprietà di trasporto direttamente nello spazio reale, eliminando la necessità di proiezioni di Wannier e permettendo l'inclusione naturale di condizioni al contorno aperte e disordine.
• Analisi risolta per strati: L'approccio fornisce un accesso alla risoluzione atomica dei profili spaziali di accumulo di spin e orbitale attraverso le interfacce e all'interno di lastre finite, distinguendo tra contributi bulk e di superficie.

Risultati
Il metodo è stato applicato a una serie di eterostrutture basate su metalli di transizione e sistemi elementari:

• Metalli di transizione bulk: I calcoli per fcc Pt, fcc Ti, fcc Cu, fcc Ni e bcc W hanno riprodotto le tendenze stabilite nella conducibilità Hall orbitale (OHC) e nella conducità Hall di spin (SHC). I risultati hanno confermato che gli orbitali $d$ dominano la risposta vicino al livello di Fermi nei metalli di transizione pesanti come Pt e W, mostrando al contempo un buon accordo con studi precedenti nello spazio dei momenti (es. Salemi et al., Go et al.).
• Nichel ferromagnetico (Ni): Nel Ni ferromagnetico, il metodo ha catturato con successo la rottura della simmetria per inversione temporale, rivelando valori finiti sia per le componenti della tensore di conducibilità pari che dispari rispetto alla inversione temporale (es. $\sigma^S_{xz}$ e $\sigma^S_{zx}$), coerentemente con le previsioni teoriche degli effetti Hall di spin e orbitali anomali.
• Accumuli in multistrati: In lastre di 18 strati di Cu, Ti, Pt, W e Ni, lo studio ha dimostrato che un sostanziale accumulo orbitale e di spin avviene vicino alle superfici, decadendo verso il centro della lastra. Notevolmente, nel Tungsteno (W), il segno dell'accumulo orbitale superficiale differiva dal segno dell'OHC bulk. Questa discrepanza è stata attribuita a differenze significative tra la densità degli stati (DOS) superficiale e quella bulk vicino al livello di Fermi, evidenziando il ruolo critico della struttura elettronica superficiale.
• Eterostrutture: Le simulazioni di bilayer consistenti in due strati ferromagnetici (Fe, Co, Ni) adiacenti a nove strati di un metallo di transizione non magnetico (Ti o W) hanno rivelato che i maggiori accumuli avvengono vicino all'interfaccia. Le eterostrutture a base di Fe hanno mostrato l'accumulo orbitale più forte, particolarmente nei sistemi Fe/Ti. Lo studio ha inoltre notato che, sebbene l'accumulo sia un proxy per la coppia (torque), la coppia effettiva è governata dalla divergenza del flusso di momento angolare, una distinzione che gli autori osservano richiedere futuri calcoli del flusso.

Significatività
Il documento sostiene che la metodologia sviluppata nello spazio reale fornisce un quadro scalabile e flessibile per simulazioni realistiche di fenomeni di trasporto orbitale. La sua importanza primaria risiede nella capacità di modellare sistemi in cui l'inomogeneità spaziale e gli effetti di dimensione finita sono cruciali, come nelle interfacce e nei multistrati disordinati, che sono difficili da trattare con gli approcci standard nello spazio dei momenti. Gli autori dimostrano che un sostanziale accumulo orbitale (di spin) può derivare anche in sistemi centrosimmetrici, guidato da asimmetrie della struttura a bande e dallo scattering d'interfaccia. Inoltre, il lavoro stabilisce che la relazione tra conducibilità bulk e accumuli superficiali è non banale ed è fortemente influenzata dalla struttura elettronica superficiale. La metodologia è presentata come uno strumento potente per investigare i fenomeni orbitronici in eterostrutture complesse, offrendo un accesso diretto alle risposte locali e aprendo la strada allo studio di sistemi con simmetria ridotta e disordine, rilevanti per i futi dispositivi spintronici e orbitronici di prossima generazione.