Theoretical study on ambient pressure superconductivity in La$_3$Ni$_2$O$_7$ thin films : structural analysis, model construction, and robustness of $s\pm$-wave pairing

Questo studio teorico dimostra che, sebbene la struttura elettronica dei film sottili di La$_3$Ni$_2$O$_7$ a pressione ambiente vari in base alla struttura cristallina e ai dettagli computazionali, la robustezza dell'accoppiamento $s\pm$-wave mediato da fluttuazioni di spin a energia finita persista, sebbene la riduzione osservata della temperatura critica rispetto al bulk pressurizzato sia meglio spiegata da modelli con piccolo hopping interstrato derivati dalle strutture reticolari sperimentali.

L'essenziale

Il Quadro Generale: Un Superconduttore a Pressione Ambiente

Immagina un materiale chiamato La₃Ni₂O₇ (un tipo di cristallo a base di nichel) che può condurre elettricità senza resistenza (superconduttività). Gli scienziati hanno scoperto recentemente che se si comprime questo materiale con una pressione enorme, diventa un superconduttore a circa 80 Kelvin (molto freddo, ma caldo per i superconduttori).

Recentemente, i ricercatori hanno trovato un modo per ottenere questo risultato senza comprimerlo. Hanno fatto crescere il materiale come un film sottilissimo su un tipo specifico di "pavimento" cristallino (un substrato). Il pavimento era leggermente più piccolo del film, il che lo schiacciava dai lati, mimando l'effetto dell'alta pressione. Questo film è diventato superconduttivo a circa 40 Kelvin.

La Domanda: Perché il film sottile funziona a una temperatura più bassa (40 K) rispetto al materiale massiccio compresso (80 K)? E qual è la "ricetta" esatta all'interno del materiale che fa fluire l'elettricità senza resistenza?

L'Approccio degli Scienziati: Costruire un Modello Digitale

Gli autori di questo documento non hanno solo indovinato; hanno costruito una simulazione al computer dettagliata. Pensateci come a un motore per videogiochi in cui hanno cercato di ricreare la fisica di questo materiale da zero.

1. La Plantilla (Struttura): Hanno esaminato la "plantilla" degli atomi. Hanno provato due plantille diverse:
   • La Plantilla Teorica: Ciò che i loro calcoli al computer dicevano che gli atomi dovrebbero essere.
   • La Plantilla Sperimentale: Ciò che gli scienziati hanno effettivamente misurato in laboratorio di recente.
2. Il Motore (FLEX): Hanno utilizzato un motore matematico complesso chiamato FLEX (Approssimazione di Scambio delle Fluttuazioni). Immaginate questo motore come un simulatore meteorologico super-preciso. Invece di prevedere la pioggia, prevede come gli elettroni (le minuscole particelle che trasportano l'elettricità) danzano e interagiscono tra loro. Tiene conto di ogni possibile mossa che gli elettroni possono fare, non solo di quelle ovvie.

Risultati Chiave: La "Danza" degli Elettroni

1. Il Mistero della "Borsa γ"

Nel mondo di questi materiali, esiste una forma specifica della folla di elettroni chiamata borsa γ. Alcuni scienziati pensavano che questa borsa fosse essenziale per la superconduttività; altri pensavano che non importasse.

• Il Verdetto del Documento: Gli autori hanno scoperto che se questa "borsa γ" esiste o scompare dipende interamente da quale plantilla si utilizza (teorica vs sperimentale) e da come si sintonizza la matematica.
• L'Analogia: È come guardare una folla attraverso occhiali di colori diversi. Attraverso una lente, vedi un gruppo distinto di persone (la borsa); attraverso un'altra, si fondono con il resto.
• Il Risultato: Sorprendentemente, non importava. Che la borsa ci fosse o meno, la superconduttività rimaneva forte. La "danza" degli elettroni era abbastanza robusta da gestire questi cambiamenti strutturali.

2. La "Colla" che Tiene Insieme

Come si accoppiano gli elettroni per condurre elettricità? Di solito, hanno bisogno di una "colla".

• Il Verdetto del Documento: La colla qui sono le fluttuazioni di spin. Immaginate gli elettroni come ballerini che ruotano su se stessi. A volte, vacillano o fluttuano nel loro spin. Questi vacillamenti agiscono come un ritmo che aiuta i ballerini ad accoppiarsi.
• La Svolta: Il documento sostiene che questo "ritmo" proviene da vacillamenti ad alta energia, non solo dai movimenti lenti e ovvi vicino alla superficie della folla di elettroni. Poiché la colla si basa su questi vacillamenti ad alta energia, la superconduttività è molto stabile e non si rompe facilmente se la forma della folla di elettroni cambia leggermente.

3. Perché il Film è più Freddo (40 K) rispetto al Massiccio (80 K)?

Questo era il più grande enigma. Il film sottile diventa superconduttivo alla metà della temperatura del materiale massiccio pressurizzato.

• Il Verdetto del Documento: La differenza si riduce a un numero specifico: $|t_\perp|$.
• L'Analogia: Pensate al materiale come a un edificio a due piani dove gli elettroni possono saltare tra i piani.
  • Nel massiccio pressurizzato, i piani sono perfettamente allineati e il salto tra di loro è forte e facile (Alto $|t_\perp|$). Questo crea una pista da ballo molto efficiente (80 K).
  • Nel film sottile, le misurazioni sperimentali mostrano che i piani sono leggermente disallineati o il salto è più debole (Basso $|t_\perp|$).
• La Conclusione: Quando gli autori hanno utilizzato la "Plantilla Sperimentale" (che mostrava questo salto più debole) nella loro simulazione, la temperatura superconduttiva è scesa esattamente ai 40 K osservati. Quando hanno utilizzato la "Plantilla Teorica" (che prevedeva un salto più forte), la temperatura è rimasta alta a 80 K.
• Il Messaggio Chiave: Il motivo per cui il film è "più debole" è probabilmente perché la struttura fisica reale del film ha una connessione più debole tra i suoi strati rispetto a quanto previsto dalla teoria.

Riepilogo in Pillole

Gli scienziati hanno costruito una simulazione high-tech per capire perché un nuovo film superconduttore funziona a pressione ambiente. Hanno scoperto che:

1. Il meccanismo di accoppiamento è resistente: Gli elettroni si accoppiano utilizzando "vacillamenti" ad alta energia (fluttuazioni di spin), rendendo la superconduttività molto robusta contro piccoli cambiamenti nella forma del materiale.
2. La "borsa γ" non importa: Che esista o meno una specifica forma di elettrone non cambia il risultato.
3. Il calo di temperatura è strutturale: Il film raggiunge solo 40 K (invece di 80 K) perché la distanza fisica reale tra gli strati atomici nel film è leggermente diversa da quanto previsto dalla teoria, rendendo il "salto" tra gli strati più debole.

Il documento dice essenzialmente: "Conosciamo la ricetta per la superconduttività in questo materiale. Il motivo per cui il film è leggermente meno efficiente del blocco pressurizzato è semplicemente perché gli strati del film non sono collegati perfettamente come pensavamo che sarebbero stati".

Sommario tecnico

1. Enunciato del Problema

La scoperta della superconduttività ad alta temperatura ($T_c \approx 80$ K) nel nickelato bilayer pressurizzato $La_3Ni_2O_7$ ha suscitato un intenso interesse. Recentemente, è stata osservata la superconduttività a pressione ambiente ($T_c \approx 40$ K) in film sottili di $La_3Ni_2O_7$ cresciuti su substrati come $SrLaAlO_4$ (SLAO). Tuttavia, rimane una sfida teorica significativa:

• La discrepanza di $T_c$: La $T_c$ nei film sottili è circa la metà di quella del sistema bulk pressurizzato.
• Ambiguità strutturale: Esiste una grande deviazione tra le strutture cristalline determinate dall'ottimizzazione della Teoria del Funzionale Densità (DFT) e quelle determinate sperimentalmente. Nello specifico, la struttura sperimentale mostra un integrale di hopping interlayer ($|t_\perp| \approx 0.4$ eV) significativamente più piccolo rispetto alle previsioni DFT ($|t_\perp| \approx 0.6$ eV).
• Controversia sulla superficie di Fermi: Studi sperimentali e teorici non concordano sull'esistenza della "$\gamma$-tasca" (una superficie di Fermi originata dalla parte superiore della banda di legame $d_{3z^2-r^2}$). La sua presenza o assenza è dibattuta per quanto riguarda il suo impatto sulla simmetria di accoppiamento.
• Obiettivo: Gli autori mirano a riprodurre teoricamente la superconduttività a pressione ambiente nei film sottili, determinare l'origine della $T_c$ ridotta e investigare la robustezza della simmetria di accoppiamento rispetto alle variazioni strutturali ed elettroniche.

2. Metodologia

Lo studio adotta un approccio multistadio che combina calcoli da primi principi con la teoria delle perturbazioni a molti corpi:

• Ottimizzazione Strutturale:

  • Sono stati eseguiti calcoli da primi principi utilizzando VASP con il funzionale PBEsol.
  • Le costanti reticolari nel piano ($a$) sono state fissate per corrispondere a tre substrati sperimentali: LSAT, LAO e SLAO.
  • La costante reticolare fuori dal piano ($c$) e le posizioni atomiche sono state ottimizzate per minimizzare l'energia totale.
  • Sono state testate due ipotesi di simmetria: la $I4/mmm$ (tetragonale) suggerita sperimentalmente e la $Amam$ (ortorombica) ottimizzata teoricamente.
  • È stato costruito un modello aggiuntivo utilizzando la struttura determinata sperimentalmente dalla Ref. [24] (Yue et al.), che presenta una distanza Ni-Ni distinta.

• Costruzione del Modello:

  • Funzioni di Wannier: Sono state estratte funzioni di Wannier massimamente localizzate utilizzando WANNIER90, focalizzandosi sugli orbitali Ni-$3d_{x^2-y^2}$ e $3d_{3z^2-r^2}$.
  • Hamiltoniana: È stato costruito un modello tight-binding bilayer a due orbitali.
  • Parametri di Interazione: La repulsione Coulombiana on-site ($U$), la repulsione interorbitale ($U'$), l'accoppiamento di Hund ($J$) e l'hopping di coppia ($J'$) sono stati adottati da precedenti calcoli RPA vincolati ($U \approx 3.6-3.8$ eV).
  • Struttura di Banda: I calcoli sono stati eseguiti con e senza la correzione +U (formulazione di Dudarev, $U=3$ eV) per valutare la sensibilità.

• Analisi della Superconduttività (FLEX):

  • L'approssimazione Fluctuation Exchange (FLEX) è stata applicata al modello di Hubbard a due orbitali.
  • Questo metodo calcola in modo autoconsistente l'autoenergia e l'interazione di accoppiamento mediate dalle fluttuazioni di spin, tenendo conto delle dipendenze complete da momento e frequenza della funzione di Green.
  • L'equazione di Eliashberg linearizzata è stata risolta per ottenere l'autovalore $\lambda$ (un proxy per $T_c$) e la funzione di gap superconduttivo $\Delta(k)$.

3. Contributi Chiave e Risultati

A. Tendenza Strutturale ed Elettronica

• Angolo di Legame: Al diminuire della costante reticolare nel piano $a$ (aumento della deformazione compressiva), l'angolo di legame interlayer Ni-O-Ni ($\alpha$) si avvicina a $180^\circ$, coerentemente con la transizione verso la simmetria $I4/mmm$. Tuttavia, anche per il substrato SLAO, l'ottimizzazione DFT non raggiunge esattamente $180^\circ$ (la simmetria $Amam$ persiste leggermente), sebbene la differenza di energia tra $Amam$e$I4/mmm$ sia trascurabile.
• Occupazione Orbitale: La riduzione di $a$ aumenta l'altezza dell'ossigeno apicale e lo sfasamento energetico ($\Delta E$) tra gli orbitali $d_{x^2-y^2}$ e $d_{3z^2-r^2}$. Questo sposta gli elettroni nell'orbitale $d_{3z^2-r^2}$, portandolo più vicino alla metà riempimento (occupazione $\approx 0.95$ per SLAO).
• Integrali di Hopping:
  • L'hopping interlayer $|t_\perp|$ aumenta al diminuire di $a$ nelle strutture ottimizzate DFT (raggiungendo $\sim 0.6$ eV).
  • Crucialmente, la struttura determinata sperimentalmente produce un $|t_\perp|$ molto più piccolo, $\approx 0.45$ eV, a causa di una maggiore distanza Ni-Ni.

B. Robustezza dell'Accoppiamento $s_{\pm}$-Wave

• Stabilità della Simmetria: Lo studio rileva che la simmetria di accoppiamento $s_{\pm}$-wave è robusta. Rimane la soluzione dominante indipendentemente da:
  • La specifica simmetria cristallina ($Amam$ vs $I4/mmm$).
  • La presenza o assenza della $\gamma$-tasca (che dipende dalla correzione $U$ e dal riempimento della banda).
  • I dettagli specifici della topologia della superficie di Fermi.
• Caratteristiche della Funzione di Gap:
  • Nella rappresentazione orbitale, la funzione di gap è dominata dall'accoppiamento interlayer degli orbitali $d_{3z^2-r^2}$.
  • Questo gap interlayer è quasi indipendente dal momento (accoppiamento intra-cellula unitaria).
  • La funzione di gap nella rappresentazione a bande riflette il peso orbitale: è grande dove è presente il carattere $d_{3z^2-r^2}$ e cambia segno tra le bande di legame e anti-legame.

C. Origine della $T_c$ Ridotta nei Film Sottili

• Ruolo di $t_\perp$: La ragione principale della $T_c$ dimezzata nei film sottili (rispetto al bulk pressurizzato) è identificata come la riduzione dell'integrale di hopping interlayer $|t_\perp|$.
  • I modelli basati su strutture ottimizzate DFT (grande $|t_\perp| \approx 0.6$ eV) prevedono un alto $\lambda$ simile al bulk pressurizzato.
  • I modelli basati sulla struttura determinata sperimentalmente (piccolo $|t_\perp| \approx 0.45$ eV) producono un $\lambda$ significativamente soppresso, coerente con la $T_c \approx 40$ K osservata.
• Meccanismo: La soppressione della $T_c$ con la riduzione di $|t_\perp|$ è attribuita alla natura della "colla" di accoppiamento. L'accoppiamento è mediato da fluttuazioni di spin a energia finita (piuttosto che dal nesting della superficie di Fermi a bassa energia).
  • Le fluttuazioni di spin ad alta energia (che guidano l'accoppiamento) sono soppresse quando $|t_\perp|$ è ridotto.
  • Al contrario, le fluttuazioni a bassa energia (associate al nesting della superficie di Fermi) sono potenziate dalla riduzione di $|t_\perp|$, ma queste sono meno efficaci per l'accoppiamento $s_{\pm}$ in questo sistema.

4. Significato e Conclusione

• Risoluzione dell'Enigma $T_c$: Il documento fornisce una spiegazione teorica coerente per la $T_c$ ridotta nei film sottili a pressione ambiente: la struttura cristallina sperimentale possiede intrinsecamente un hopping interlayer ($t_\perp$) più piccolo rispetto a quello previsto dall'ottimizzazione DFT standard, il quale sopprime direttamente la forza di accoppiamento mediata dalle fluttuazioni di spin a energia finita.
• Validazione del Meccanismo di Accoppiamento: Lo studio rafforza la visione secondo cui la superconduttività in $La_3Ni_2O_7$ è guidata da fluttuazioni di spin a energia finita che mediano un accoppiamento interlayer $d_{3z^2-r^2}$. Questo meccanismo è robusto rispetto alle variazioni nella topologia della superficie di Fermi (come la presenza della $\gamma$-tasca) e ai dettagli della simmetria strutturale.
• Coerenza Teorica: I risultati si allineano con approcci non perturbativi (cDMFT, DCA, VMC), confermando che il sistema può essere visto come un modello di Hubbard bilayer (formato dagli orbitali $d_{3z^2-r^2}$) immerso in un ambiente multiorbitale.
• Direzioni Future: La grande discrepanza tra le strutture cristalline previste dalla DFT e quelle osservate sperimentalmente (in particolare la distanza Ni-Ni) rimane una questione aperta. Gli autori suggeriscono che ordini competitivi (ad esempio, onde di densità) non catturati dallo schema FLEX corrente potrebbero svolgere un ruolo anche nella stabilità e nelle proprietà del materiale reale.

In sintesi, il documento stabilisce che mentre la simmetria di accoppiamento ($s_{\pm}$) è robusta, l'entità della $T_c$ è altamente sensibile all'hopping interlayer $t_\perp$, che è dettato dalla precisa struttura atomica del film sottile.