First-principles study of bulk stacking, $J_{\rm eff}$ picture, magnetic Hamiltonian, $g$ factors, and structural distortions of $α$-RuCl$_3$

Sintesi Tecnica: Studio di Primo Principio su Impilamento Bulk, Immagine $J_{\text{eff}}$, Hamiltoniana Magnetica, Fattori $g$ e Distorsioni Strutturali di $\alpha$-RuCl$_3$

Problematica
$\alpha$-RuCl$_3$ è un candidato primario per la realizzazione dello stato di liquido di spin quantistico di Kitaev, tuttavia le sue proprietà magnetiche fondamentali rimangono oggetto di dibattito. Le questioni chiave ancora irrisolte includono:

1. Struttura Bulk: Sebbene studi recenti suggeriscano che la fase a bassa temperatura appartenga al gruppo spaziale $R\bar{3}$ (con una transizione a $C2/m$ a $\sim$150 K), non vi era precedente studio teorico che confermasse la struttura $R\bar{3}$ come la fase energeticamente favorita a bassa temperatura.
2. Struttura Elettronica: La validità dell'immagine $J_{\text{eff}} = 1/2$ per $\alpha$-RuCl$_3$ richiede una verifica, in particolare riguardo all'orientamento dell'asse di quantizzazione del momento angolare rispetto al vettore Néel.
3. Hamiltoniana Magnetica: Studi precedenti riportano discrepanze significative nei parametri di scambio magnetico (ad esempio, il segno dello scambio di Kitaev a primo vicino $K_1$) e spesso trascurano le interazioni al secondo vicino.
4. Anisotropia del Fattore $g$: I rapporti sperimentali sull'anisotropia del fattore $g$ ($g_{XY}$ vs. $g_Z$) sono contraddittori, variando da anisotropie piccole a grandi. Inoltre, i metodi teorici per calcolare questi fattori (specificamente la proiezione dell'orbitale atomico) potrebbero essere insufficienti.
5. Distorsioni Strutturali: L'impatto di specifiche distorsioni nell'ottaedro di RuCl$_6$ (distorsione trigonale vs. rotazione relativa dei triangoli di Cl) sull'anisotropia magnetica non è stato isolato sistematicamente.

Metodologia
Gli autori hanno impiegato la teoria del funzionale della densità con vincoli (cDFT) all'interno del framework Quantum ESPRESSO utilizzando l'approssimazione di densità locale (LDA) e i metodi LDA+U invarianti per rotazione.

• Rilassamento Strutturale: I calcoli sono stati eseguiti sia per le configurazioni di impilamento $C2/m$ che $R\bar{3}$ sotto diverse condizioni di rilassamento (struttura sperimentale, cella fissa e rilassamento a cella variabile) per determinare la stabilità energetica.
• Analisi $J_{\text{eff}}$: La struttura elettronica è stata analizzata proiettando le funzioni d'onda sugli stati di base $J_{\text{eff}}$. Crucialmente, l'asse di quantizzazione del momento angolare è stato allineato parallelamente al vettore Néel per testare l'immagine $J_{\text{eff}}$.
• Costruzione dell'Hamiltoniana Magnetica: Gli autori hanno adattato le energie totali da varie configurazioni magnetiche (zigzag, Néel, stripy, ferromagnetica e una specifica configurazione 1-per-4) a tre hamiltoniane modello:
  1. Modello convenzionale Kitaev-Heisenberg (KH).
  2. Modello Spin-ice (SI).
  3. Un nuovo modello Anisotropo $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) che include interazioni anisotrope al secondo e terzo vicino.
• Calcolo del Fattore $g$: Sono stati confrontati due metodi:
  1. Proiezione dell'Orbitale Atomico: Proiezione dei momenti sul manifold di Hubbard.
  2. Inclinazione del Momento (Moment Canting) con Interpolazione Wannier: Inclinazione del momento $J_{\text{eff}}$ di $1^\circ$ e calcolo della variazione del momento magnetico totale tramite interpolazione Wannier traslazionale equivariante (WannierBerri).
• Analisi delle Distorsioni: Gli effetti di due specifiche distorsioni dall'ottaedro perfetto sono stati isolati: la distorsione trigonale (cambio dell'angolo Ru-Cl-Z) e la rotazione relativa (rotazione dei triangoli di Cl superiore e inferiore nel piano XY).

Contributi Chiave e Risultati

1. Conferma dell'Impilamento $R\bar{3}$:
   Lo studio dimostra che la struttura $R\bar{3}$ è inferiore in energia rispetto alla struttura $C2/m$ di circa 5–6 meV/Ru attraverso diversi valori di Hubbard $U$ e tipi di rilassamento. Ciò fornisce il primo supporto teorico per la fase $R\bar{3}$ osservata sperimentalmente a bassa temperatura.

2. Immagine $J_{\text{eff}}$ Raffinata:
   Si dimostra che gli stati del minimo della banda di conduzione sono quasi puramente di carattere $J_{\text{eff}} = 1/2$ e $m_{\text{eff}} = -1/2$, ma solo quando l'asse di quantizzazione del momento angolare è impostato parallelamente al vettore Néel. Questa prospettiva, precedentemente trascurata, è essenziale per identificare correttamente il carattere orbitale degli stati. L'entità del momento $J_{\text{eff}}$ si è rivelata robusta tra le diverse configurazioni, sebbene leggermente inferiore all'ideale $1/2$ a causa dell'ibridazione e delle distorsioni.

3. Hamiltoniana Magnetica e Parametri di Scambio:

   • Il modello convenzionale KH e il modello Spin-ice non riescono ad aiutare accuratamente a riprodurre le energie totali DFT, in particolare per la configurazione Néel.
   • Il modello proposto A-$J_1J_2J_3$ (e una variante semplificata A-$J_1J_2$-I-$J_3$) riproduce accuratamente le energie DFT.
   • Le interazioni al secondo vicino sono critiche. Le componenti anisotrope dello scambio al secondo vicino ($K_2, \Gamma_2, \Gamma'_2$) sono non trascurabili e, in magnitudo, superiori ai termini anisotropi al primo vicino ($\Gamma'_1$) spesso considerati nei modelli precedenti. L'inclusione di questi termini ha ridotto l'errore quadratico medio (RMS) nell'energia totale di un fattore di 3,6 rispetto al modello convenzionale.

4. Anisotropia del Fattore $g$ e Metodo di Calcolo:

   • Il metodo dell'inclinazione del momento combinato con l'interpolazione Wannier fornisce fattori $g$ di $g_{XY} \approx 2,22$ e $g_Z \approx 1,84$. Questi valori indicano una piccola anisotropia, allineandosi con i recenti rapporti sperimentali (ad es. Loidl et al.) e contraddicendo le precedenti affermazioni di grande anisotropia.
   • Il metodo di proiezione dell'orbitale atomico sottostima significativamente l'entità del fattore $g$ ($g_{XY} \approx 1,54$), mancando una parte non trascurabile del momento magnetico totale. Gli autori concludono che il calcolo accurato del magnetismo orbitale richiede l'approccio di interpolazione Wannier.

5. Ruolo delle Distorsioni Strutturali:
   Isolando le distorsioni, lo studio trova che la rotazione relativa dei triangoli di Cl (all'interno del piano XY) è il fattore dominante che determina l'anisotropia magnetica e la corretta direzione del momento dello stato fondamentale ($\theta \sim 60^\circ$). La distorsione trigonale da sola non riesce a riprodurre lo stato fondamentale sperimentale. Questa distorsione di rotazione, presente nel materiale bulk senza perturbazioni esterne, è stata trascurata nei modelli teorici precedenti.

Significato
L'articolo sostiene di fornire una comprensione completa, basata sui primi principi, di $\alpha$-RuCl$_3$ risolvendo le lunghe controversie nelle sue descrizioni strutturali e magnetiche. In particolare:

• Valida teoricamente la struttura a bassa temperatura $R\bar{3}$.
• Stabilisce la necessità di allineare l'asse di quantizzazione con il vettore Néel per interpretare correttamente l'immagine $J_{\text{eff}}$.
• Dimostra che la modellazione magnetica accurata richiede l'inclusione di interazioni di scambio anisotrope al secondo vicino.
• Corregge la metodologia per il calcolo dei fattori $g$, mostrando che la proiezione dell'orbitale atomico è insufficiente e che la vera anisotropia è piccola.
• Identifica la rotazione relativa dei triangoli di Cl come un parametro strutturale critico, precedentemente trascurato, che governa l'anisotropia magnetica, suggerendo una nuova via per l'ingegneria del magnetismo nei materiali di Kitaev.

Gli autori concludono che il loro lavoro offre una prospettiva differente sul magnetismo di $\alpha$-RuCl$_3$ e suggeriscono che i futuri modelli teorici e gli esperimenti ad alta pressione debbano incorporare queste specifiche distorsioni strutturali e termini di scambio per un'interpretazione accurata.