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

Deze theoretische studie toont aan dat, hoewel de elektronische structuur van La$_3$Ni$_2$O$_7$-dunne films bij omgevingsdruk varieert met kristalstructuur en computationele details, de robuustheid van $s\pm$-golfkoppeling gemedieerd door spinfluctuaties met eindige energie behouden blijft, alhoewel de waargenomen verlaging van de kritische temperatuur ten opzichte van onder druk staand bulkmateriaal het beste wordt verklaard door modellen met kleine interlaag-hopping die zijn afgeleid van experimentele roosterstructuren.

De kern

Het Grote Geheel: Een Supergeleider bij Omgevingsdruk

Stel je een materiaal voor genaamd La₃Ni₂O₇ (een soort nikkelgebaseerd kristal) dat elektriciteit kan geleiden zonder weerstand (supergeleiding). Wetenschappers ontdekten recentelijk dat als je dit materiaal met enorme druk samendrukt, het een supergeleider wordt bij ongeveer 80 Kelvin (zeer koud, maar warm voor supergeleiders).

Recentelijk vonden onderzoekers een manier om dit te laten gebeuren zonder het te comprimeren. Ze kweekten het materiaal als een zeer dunne film op een specifiek type kristal "vloer" (een substraat). De vloer was iets kleiner dan de film, waardoor de film van de zijkanten werd samengedrukt, wat het effect van hoge druk nabootste. Deze film werd supergeleidend bij ongeveer 40 Kelvin.

De Vraag: Waarom werkt de dunne film bij een lagere temperatuur (40 K) dan het samengeperste bulkmateriaal (80 K)? En wat is het exacte "recept" binnenin het materiaal dat zorgt voor de weerstandloze stroom van elektriciteit?

De Aanpak van de Wetenschappers: Een Digitaal Model Bouwen

De auteurs van dit artikel gokten niet zomaar; ze bouwden een gedetailleerde computersimulatie. Denk hierbij aan een videospel-engine waarmee ze de fysica van dit materiaal vanaf nul probeerden na te bootsen.

1. De Blauwdruk (Structuur): Ze keken naar de "blauwdruk" van de atomen. Ze probeerden twee verschillende blauwdrukken:
   • De Theoretische Blauwdruk: Wat hun computerberekeningen zeiden dat de atomen zouden moeten zijn.
   • De Experimentele Blauwdruk: Wat wetenschappers recentelijk daadwerkelijk in het lab hebben gemeten.
2. De Engine (FLEX): Ze gebruikten een complexe wiskundige engine genaamd FLEX (Fluctuation Exchange Approximation). Stel je deze engine voor als een super-nauwkeurige weersimulator. In plaats van regen te voorspellen, voorspelt het hoe elektronen (de tiny deeltjes die elektriciteit dragen) dansen en met elkaar interageren. Het houdt rekening met elke mogelijke beweging die elektronen kunnen maken, niet alleen de voor de hand liggende.

Belangrijkste Bevindingen: De "Dans" van Elektronen

1. Het "γ-Pocket" Mysterie

In de wereld van deze materialen is er een specifieke vorm van de elektronenmenigte genaamd de γ-pocket. Sommige wetenschappers dachten dat deze pocket essentieel was voor supergeleiding; anderen dachten dat het er niet toe deed.

• Het Oordeel van het Artikel: De auteurs vonden dat of deze "γ-pocket" bestaat of verdwijnt, volledig afhangt van welke blauwdruk je gebruikt (theoretisch versus experimenteel) en hoe je de wiskunde bijstelt.
• De Analogie: Het is alsof je door verschillende gekleurde brillen naar een menigte kijkt. Door één lens zie je een duidelijke groep mensen (de pocket); door een andere lens vloeien ze samen.
• Het Resultaat: Verrassend genoeg maakte het niet uit. Of de pocket er nu was of niet, de supergeleiding bleef sterk. De "dans" van de elektronen was robuust genoeg om deze structurele veranderingen aan te kunnen.

2. De "Lijm" die het Bij elkaar Houdt

Hoe paren de elektronen zich om elektriciteit te geleiden? Normaal gesproken hebben ze een "lijm" nodig.

• Het Oordeel van het Artikel: De lijm hier is spinfluctuaties. Stel je de elektronen voor als dansers die ronddraaien. Soms waggelen of fluctueren ze in hun spin. Deze waggels fungeren als een ritmisch beat dat de dansers helpt paren.
• De Twist: Het artikel betoogt dat deze "beat" komt van hoge-energie waggels, niet alleen van de trage, voor de hand liggende bewegingen nabij het oppervlak van de elektronenmenigte. Omdat de lijm gebaseerd is op deze hoge-energie waggels, is de supergeleiding zeer stabiel en breekt deze niet gemakkelijk als de vorm van de elektronenmenigte licht verandert.

3. Waarom is de Film Koeler (40 K) dan het Bulk (80 K)?

Dit was het grootste raadsel. De dunne film is supergeleidend bij de helft van de temperatuur van het samengeperste bulkmateriaal.

• Het Oordeel van het Artikel: Het verschil komt neer op één specifiek getal: $|t_\perp|$.
• De Analogie: Denk aan het materiaal als een twee verdiepingen tellend gebouw waar elektronen tussen de verdiepingen kunnen springen.
  • In het samengeperste bulk zijn de verdiepingen perfect uitgelijnd, en is de sprong ertussen sterk en makkelijk (Hoge $|t_\perp|$). Dit zorgt voor een zeer efficiënte dansvloer (80 K).
  • In de dunne film tonen de experimentele metingen aan dat de verdiepingen licht uitgelijnd zijn of de sprong zwakker is (Lage $|t_\perp|$).
• De Conclusie: Toen de auteurs de "Experimentele Blauwdruk" (die deze zwakkere sprong liet zien) in hun simulatie gebruikten, daalde de supergeleidende temperatuur precies naar de waargenomen 40 K. Toen ze de "Theoretische Blauwdruk" gebruikten (die een sterkere sprong voorspelde), bleef de temperatuur hoog op 80 K.
• De Kernboodschap: De reden dat de film "zwakker" is, ligt waarschijnlijk aan het feit dat de daadwerkelijke fysische structuur van de film een zwakkere verbinding tussen zijn lagen heeft dan de theorie voorspelde.

Samenvatting in het Kort

De wetenschappers bouwden een high-tech simulatie om te begrijpen waarom een nieuwe supergeleidende film werkt bij omgevingsdruk. Ze ontdekten dat:

1. Het koppelingsmechanisme stevig is: De elektronen paren zich met behulp van hoge-energie "waggels" (spinfluctuaties), waardoor de supergeleiding zeer robuust is tegen kleine veranderingen in de vorm van het materiaal.
2. De "γ-pocket" niet uitmaakt: Of een specifieke elektronenvorm nu bestaat of niet, verandert het resultaat niet.
3. De temperatuurdaling structureel is: De film bereikt slechts 40 K (in plaats van 80 K) omdat de daadwerkelijke fysieke afstand tussen atoomlagen in de film iets anders is dan wat de theorie voorspelde, waardoor de "sprong" tussen de lagen zwakker is.

Het artikel zegt in feite: "We kennen het recept voor supergeleiding in dit materiaal. De reden dat de film iets minder efficiënt is dan het samengeperste blok, is simpelweg omdat de lagen van de film niet helemaal zo perfect verbonden zijn als we dachten dat ze zouden zijn."

Technische samenvatting

1. Probleemstelling

De ontdekking van supergeleiding bij hoge temperatuur ($T_c \approx 80$ K) in onder druk gezette bilagen nikkelaat $La_3Ni_2O_7$ heeft intense belangstelling gewekt. Onlangs werd supergeleiding bij omgevingsdruk ($T_c \approx 40$ K) waargenomen in dunne films van $La_3Ni_2O_7$ die op substraten zoals $SrLaAlO_4$ (SLAO) zijn gekweekt. Er blijft echter een aanzienlijke theoretische uitdaging bestaan:

• De $T_c$-discrepantie: De $T_c$ in dunne films is ongeveer de helft van die van het onder druk gezette bulk-systeem.
• Structurele ambiguïteit: Er is een grote afwijking tussen kristalstructuren die zijn bepaald via optimalisatie met Dichtefunctietheorie (DFT) en die experimenteel zijn bepaald. Specifiek toont de experimentele structuur een aanzienlijk kleinere interlaag hopping-integraal ($|t_\perp| \approx 0,4$ eV) in vergelijking met DFT-predicties ($|t_\perp| \approx 0,6$ eV).
• Fermi-oppervlak controverse: Experimentele en theoretische studies zijn het oneens over het bestaan van de "$\gamma$-pocket" (een Fermi-oppervlak dat voortkomt uit de top van de $d_{3z^2-r^2}$-bindingsband). De aanwezigheid of afwezigheid ervan wordt besproken met betrekking tot de impact op de paringsymmetrie.
• Doel: De auteurs streven ernaar om theoretisch de supergeleiding bij omgevingsdruk in dunne films na te bootsen, de oorsprong van de verlaagde $T_c$ te bepalen en de robuustheid van de paringsymmetrie te onderzoeken tegenover structurele en elektronische variaties.

2. Methodologie

De studie hanteert een meervoudige aanpak die berekeningen uit eerste principes combineert met veeldeeltjes-stoornistheorie:

• Structurele optimalisatie:

  • Berekeningen uit eerste principes werden uitgevoerd met VASP en de PBEsol-functionaal.
  • De in-vlak roosterconstanten ($a$) werden vastgezet om overeen te komen met drie experimentele substraten: LSAT, LAO en SLAO.
  • De uit-vlak roosterconstante ($c$) en atoomposities werden geoptimaliseerd om de totale energie te minimaliseren.
  • Twee symmetrie-aannames werden getest: de experimenteel voorgestelde $I4/mmm$ (tetragonaal) en de theoretisch geoptimaliseerde $Amam$ (orthorombisch).
  • Een extra model werd geconstrueerd met behulp van de experimenteel bepaalde structuur uit Ref. [24] (Yue et al.), die een onderscheidende Ni-Ni-afstand vertoont.

• Modelconstructie:

  • Wannier-functies: Maximaal gelokaliseerde Wannier-functies werden geëxtraheerd met WANNIER90, met focus op de Ni-$3d_{x^2-y^2}$ en $3d_{3z^2-r^2}$-orbitalen.
  • Hamiltoniaan: Een tweedimensionaal bilagen tight-binding-model werd geconstrueerd.
  • Interactieparameters: On-site Coulomb-afstoting ($U$), interorbitale afstoting ($U'$), Hund-koppeling ($J$) en paar-hopping ($J'$) werden overgenomen uit eerdere berekeningen met beperkte RPA ($U \approx 3,6-3,8$ eV).
  • Bandstructuur: Berekeningen werden uitgevoerd met en zonder de +U-correctie (Dudarev's formulering, $U=3$ eV) om de gevoeligheid te beoordelen.

• Analyse van supergeleiding (FLEX):

  • De Fluctuation Exchange (FLEX)-benadering werd toegepast op het tweedimensionale Hubbard-model.
  • Deze methode berekent zelfconsistent de zelfenergie en de paringsinteractie die wordt bemiddeld door spinfluctuaties, waarbij rekening wordt gehouden met de volledige impuls- en frequentie-afhankelijkheid van de Green-functie.
  • De gelijneerde Eliashvergelijking werd opgelost om de eigenwaarde $\lambda$ (een proxy voor $T_c$) en de supergeleidende gapfunctie $\Delta(k)$ te verkrijgen.

3. Belangrijkste bijdragen en resultaten

A. Structurele en elektronische trends

• Bindingshoek: Naarmate de in-vlak roosterconstante $a$ afneemt (toenemende compressieve spanning), nadert de interlaag Ni-O-Ni-bindinghoek ($\alpha$) $180^\circ$, in overeenstemming met de overgang naar $I4/mmm$-symmetrie. Echter, zelfs voor het SLAO-substraat bereikt DFT-optimalisatie niet exact $180^\circ$ (de $Amam$-symmetrie blijft lichtjes bestaan), hoewel het energieverschil tussen $Amam$ en $I4/mmm$ verwaarloosbaar is.
• Orbitale bezetting: Het verkleinen van $a$ verhoogt de hoogte van de apicale zuurstof en het niveauverschil ($\Delta E$) tussen de $d_{x^2-y^2}$ en $d_{3z^2-r^2}$-orbitalen. Dit verschuift elektronen naar de $d_{3z^2-r^2}$-orbitaal, waardoor deze dichter bij halfvulling komt (bezetting $\approx 0,95$ voor SLAO).
• Hopping-integralen:
  • De interlaag hopping $|t_\perp|$ neemt toe naarmate $a$ afneemt in DFT-geoptimaliseerde structuren (bereikend $\sim 0,6$ eV).
  • Cruciaal is dat de experimenteel bepaalde structuur een veel kleinere $|t_\perp| \approx 0,45$ eV oplevert vanwege een grotere Ni-Ni-afstand.

B. Robuustheid van $s_{\pm}$-golf paring

• Symmetrie-stabiliteit: De studie concludeert dat de $s_{\pm}$-golf paringsymmetrie robuust is. Het blijft de dominante oplossing ongeacht:
  • De specifieke kristalsymmetrie ($Amam$ versus $I4/mmm$).
  • De aanwezigheid of afwezigheid van de $\gamma$-pocket (die afhangt van de $U$-correctie en bandvulling).
  • De specifieke details van de topologie van het Fermi-oppervlak.
• Kenmerken van de gapfunctie:
  • In de orbitale representatie wordt de gapfunctie gedomineerd door de interlaagparing van $d_{3z^2-r^2}$-orbitalen.
  • Deze interlaaggap is bijna impuls-onafhankelijk (paring binnen de eenheidscel).
  • De gapfunctie in de bandrepresentatie weerspiegelt de orbitale gewicht: deze is groot waar $d_{3z^2-r^2}$-karakter aanwezig is en wisselt van teken tussen bindings- en antibindingsbanden.

C. Oorsprong van de verlaagde $T_c$ in dunne films

• Rol van $t_\perp$: De primaire reden voor de gehalveerde $T_c$ in dunne films (in vergelijking met onder druk gezet bulk) wordt geïdentificeerd als de reductie van de interlaag hopping-integraal $|t_\perp|$.
  • Modellen gebaseerd op DFT-geoptimaliseerde structuren (grote $|t_\perp| \approx 0,6$ eV) voorspellen een hoge $\lambda$ die vergelijkbaar is met het onder druk gezette bulk.
  • Modellen gebaseerd op de experimenteel bepaalde structuur (kleine $|t_\perp| \approx 0,45$ eV) leveren een aanzienlijk onderdrukte $\lambda$ op, in overeenstemming met de waargenomen $T_c \approx 40$ K.
• Mechanisme: De onderdrukking van $T_c$ bij verminderde $|t_\perp|$ wordt toegeschreven aan de aard van de paringslijm. De paring wordt bemiddeld door spinfluctuaties met eindige energie (in plaats van lage-energie nesten van het Fermi-oppervlak).
  • Spinfluctuaties met hoge energie (die de paring aandrijven) worden onderdrukt wanneer $|t_\perp|$ wordt verkleind.
  • Omgekeerd worden lage-energie fluctuaties (geassocieerd met nesten van het Fermi-oppervlak) versterkt door verminderde $|t_\perp|$, maar deze zijn minder effectief voor $s_{\pm}$-paring in dit systeem.

4. Betekenis en conclusie

• Oplossing van de $T_c$-puzzel: Het artikel biedt een samenhangende theoretische verklaring voor de verlaagde $T_c$ in dunne films bij omgevingsdruk: de experimentele kristalstructuur bezit inherent een kleinere interlaag hopping ($t_\perp$) dan voorspeld door standaard DFT-optimalisatie, wat de paringssterkte die wordt bemiddeld door spinfluctuaties met eindige energie direct onderdrukt.
• Validatie van het paringsmechanisme: De studie versterkt het inzicht dat supergeleiding in $La_3Ni_2O_7$ wordt aangedreven door spinfluctuaties met eindige energie die een interlaag $d_{3z^2-r^2}$-paring bemiddelen. Dit mechanisme is robuust tegenover variaties in de topologie van het Fermi-oppervlak (zoals de aanwezigheid van de $\gamma$-pocket) en details van de structurele symmetrie.
• Theoretische consistentie: De resultaten sluiten aan bij niet-perturbatieve benaderingen (cDMFT, DCA, VMC), wat bevestigt dat het systeem kan worden beschouwd als een bilagen Hubbard-model (gevormd door $d_{3z^2-r^2}$-orbitalen) ingebed in een multiorbitale omgeving.
• Toekomstige richtingen: De grote discrepantie tussen DFT-voorspelde en experimenteel waargenomen kristalstructuren (specifiek de Ni-Ni-afstand) blijft een open vraag. De auteurs suggereren dat concurrerende ordeningen (bijvoorbeeld dichtheidsgolven) die niet worden vastgelegd in het huidige FLEX-schema, ook een rol kunnen spelen in de stabiliteit en eigenschappen van het feitelijke materiaal.

Kortom, het artikel stelt vast dat hoewel de paringsymmetrie ($s_{\pm}$) robuust is, de grootte van $T_c$ zeer gevoelig is voor de interlaag hopping $t_\perp$, die wordt bepaald door de precieze atoomstructuur van de dunne film.