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

Cette étude théorique démontre que, bien que la structure électronique des films minces de La$_3$Ni$_2$O$_7$ à pression ambiante varie selon la structure cristalline et les détails computationnels, la robustesse de l'appariement de type onde $s\pm$ médié par des fluctuations de spin à énergie finie persiste, bien que la réduction observée de la température critique par rapport au volume sous pression s'explique au mieux par des modèles présentant un faible saut intercouche dérivé des structures de réseau expérimentales.

L'essentiel

La Grande Image : Un supraconducteur à pression ambiante

Imaginez un matériau appelé La₃Ni₂O₇ (un type de cristal à base de nickel) capable de conduire l'électricité sans aucune résistance (supraconductivité). Des scientifiques ont récemment découvert que si vous serrez ce matériau avec une pression massive, il devient supraconducteur à environ 80 Kelvin (très froid, mais chaud pour les supraconducteurs).

Récemment, des chercheurs ont trouvé un moyen de faire cela sans le serrer. Ils ont fait pousser le matériau sous forme d'un film très mince sur un type spécifique de « sol » cristallin (un substrat). Le sol était légèrement plus petit que le film, ce qui a serré le film sur les côtés, imitant l'effet d'une haute pression. Ce film est devenu supraconducteur à environ 40 Kelvin.

La Question : Pourquoi le film mince fonctionne-t-il à une température plus basse (40 K) que le matériau massif comprimé (80 K) ? Et quelle est la « recette » exacte à l'intérieur du matériau qui permet au courant de circuler sans résistance ?

L'Approche des Scientifiques : Construire un Modèle Numérique

Les auteurs de ce papier n'ont pas seulement deviné ; ils ont construit une simulation informatique détaillée. Imaginez cela comme un moteur de jeu vidéo où ils ont tenté de recréer la physique de ce matériau à partir de zéro.

1. Le Plan (Structure) : Ils ont examiné le « plan » des atomes. Ils ont essayé deux plans différents :
   • Le Plan Théorique : Ce que leurs calculs informatiques disaient que les atomes devraient être.
   • Le Plan Expérimental : Ce que les scientifiques ont réellement mesuré en laboratoire récemment.
2. Le Moteur (FLEX) : Ils ont utilisé un moteur mathématique complexe appelé FLEX (Approximation d'Échange de Fluctuations). Imaginez ce moteur comme un simulateur de météo ultra-précis. Au lieu de prédire la pluie, il prédit comment les électrons (les minuscules particules transportant l'électricité) dansent et interagissent entre eux. Il prend en compte chaque mouvement possible que les électrons peuvent faire, pas seulement les plus évidents.

Résultats Clés : La « Danse » des Électrons

1. Le Mystère de la « Poche γ »

Dans le monde de ces matériaux, il existe une forme spécifique de la foule d'électrons appelée la poche γ. Certains scientifiques pensaient que cette poche était essentielle à la supraconductivité ; d'autres pensaient qu'elle n'avait pas d'importance.

• Le Verdict du Papier : Les auteurs ont découvert que la présence ou la disparition de cette « poche γ » dépend entièrement du plan que vous utilisez (théorique vs expérimental) et de la façon dont vous ajustez les mathématiques.
• L'Analogie : C'est comme regarder une foule à travers des lunettes de couleurs différentes. À travers une lentille, vous voyez un groupe distinct de personnes (la poche) ; à travers une autre, ils se fondent dans le décor.
• Le Résultat : Étonnamment, cela n'avait pas d'importance. Que la poche soit là ou non, la supraconductivité restait forte. La « danse » des électrons était suffisamment robuste pour supporter ces changements structurels.

2. La « Colle » qui Tout Maintient

Comment les électrons s'apparient-ils pour conduire l'électricité ? Habituellement, ils ont besoin d'une « colle ».

• Le Verdict du Papier : La colle ici est constituée de fluctuations de spin. Imaginez les électrons comme des danseurs qui tournent sur eux-mêmes. Parfois, ils vacillent ou fluctuent dans leur spin. Ces vacillements agissent comme un rythme qui aide les danseurs à s'apparier.
• La Surprise : Le papier soutient que ce « rythme » provient de vacillements à haute énergie, et non pas seulement des mouvements lents et évidents près de la surface de la foule d'électrons. Parce que la colle est basée sur ces vacillements à haute énergie, la supraconductivité est très stable et ne se brise pas facilement si la forme de la foule d'électrons change légèrement.

3. Pourquoi le Film est-il plus Froid (40 K) que le Massif (80 K) ?

C'était le plus grand mystère. Le film mince est supraconducteur à la moitié de la température du matériau massif sous pression.

• Le Verdict du Papier : La différence se résume à un nombre spécifique : $|t_\perp|$.
• L'Analogie : Imaginez le matériau comme un bâtiment à deux étages où les électrons peuvent sauter entre les étages.
  • Dans le massif sous pression, les étages sont parfaitement alignés, et le saut entre eux est fort et facile (Haut $|t_\perp|$). Cela crée une piste de danse très efficace (80 K).
  • Dans le film mince, les mesures expérimentales montrent que les étages sont légèrement désalignés ou que le saut est plus faible (Bas $|t_\perp|$).
• La Conclusion : Lorsque les auteurs ont utilisé le « Plan Expérimental » (qui montrait ce saut plus faible) dans leur simulation, la température supraconductrice est tombée exactement aux 40 K observés. Lorsqu'ils ont utilisé le « Plan Théorique » (qui prédisait un saut plus fort), la température est restée élevée à 80 K.
• L'Essentiel : La raison pour laquelle le film est « plus faible » est probablement que la structure physique réelle du film a une connexion plus faible entre ses couches que ce que la théorie prévoyait.

Résumé en Bref

Les scientifiques ont construit une simulation de haute technologie pour comprendre pourquoi un nouveau film supraconducteur fonctionne à pression ambiante. Ils ont découvert que :

1. Le mécanisme d'appariement est robuste : Les électrons s'apparient en utilisant des « vacillements » à haute énergie (fluctuations de spin), rendant la supraconductivité très résistante aux petits changements de forme du matériau.
2. La « poche γ » n'a pas d'importance : Que la présence d'une forme spécifique d'électrons existe ou non ne change pas le résultat.
3. La chute de température est structurelle : Le film n'atteint que 40 K (au lieu de 80 K) parce que la distance physique réelle entre les couches atomiques dans le film est légèrement différente de ce que la théorie prévoyait, rendant le « saut » entre les couches plus faible.

Le papier dit essentiellement : « Nous connaissons la recette de la supraconductivité dans ce matériau. La raison pour laquelle le film est légèrement moins efficace que le bloc sous pression est simplement que les couches du film ne sont pas tout à fait aussi parfaitement connectées que nous le pensions. »

Résumé technique

1. Énoncé du problème

La découverte de la supraconductivité à haute température ($T_c \approx 80$ K) dans le nickelate bicouche $La_3Ni_2O_7$ sous pression a suscité un vif intérêt. Récemment, une supraconductivité à pression ambiante ($T_c \approx 40$ K) a été observée dans des films minces de $La_3Ni_2O_7$ déposés sur des substrats tels que $SrLaAlO_4$ (SLAO). Cependant, un défi théorique majeur subsiste :

• La divergence de $T_c$ : La $T_c$ dans les films minces est environ la moitié de celle du système massif sous pression.
• L'ambiguïté structurale : Il existe une grande divergence entre les structures cristallines déterminées par l'optimisation de la théorie de la fonctionnelle de la densité (DFT) et celles déterminées expérimentalement. Plus précisément, la structure expérimentale montre un intégrale de saut intercouche ($|t_\perp| \approx 0,4$ eV) significativement plus faible que les prédictions DFT ($|t_\perp| \approx 0,6$ eV).
• La controverse sur la surface de Fermi : Les études expérimentales et théoriques ne s'accordent pas sur l'existence de la « poche $\gamma$ » (une surface de Fermi provenant du sommet de la bande de liaison $d_{3z^2-r^2}$). Sa présence ou son absence fait débat quant à son impact sur la symétrie de l'appariement.
• Objectif : Les auteurs visent à reproduire théoriquement la supraconductivité à pression ambiante dans les films minces, à déterminer l'origine de la réduction de $T_c$, et à étudier la robustesse de la symétrie de l'appariement face aux variations structurelles et électroniques.

2. Méthodologie

L'étude emploie une approche multi-étapes combinant des calculs ab initio avec la théorie de la perturbation à N corps :

• Optimisation structurale :

  • Des calculs ab initio ont été effectués en utilisant VASP avec le fonctionnel PBEsol.
  • Les constantes de réseau dans le plan ($a$) ont été fixées pour correspondre à trois substrats expérimentaux : LSAT, LAO et SLAO.
  • La constante de réseau hors plan ($c$) et les positions atomiques ont été optimisées pour minimiser l'énergie totale.
  • Deux hypothèses de symétrie ont été testées : la $I4/mmm$ (tétragonale) suggérée expérimentalement et la $Amam$ (orthorhombique) optimisée théoriquement.
  • Un modèle supplémentaire a été construit en utilisant la structure déterminée expérimentalement de la référence [24] (Yue et al.), qui présente une distance Ni-Ni distincte.

• Construction du modèle :

  • Fonctions de Wannier : Des fonctions de Wannier maximalement localisées ont été extraites en utilisant WANNIER90, en se concentrant sur les orbitales Ni-$3d_{x^2-y^2}$ et $3d_{3z^2-r^2}$.
  • Hamiltonien : Un modèle de liaison forte bicouche à deux orbitales a été construit.
  • Paramètres d'interaction : La répulsion de Coulomb sur site ($U$), la répulsion interorbitale ($U'$), le couplage de Hund ($J$) et le saut de paires ($J'$) ont été adoptés à partir de calculs antérieurs RPA contraints ($U \approx 3,6-3,8$ eV).
  • Structure de bande : Des calculs ont été effectués avec et sans la correction +U (formulation de Dudarev, $U=3$ eV) pour évaluer la sensibilité.

• Analyse de la supraconductivité (FLEX) :

  • L'approximation Fluctuation Exchange (FLEX) a été appliquée au modèle de Hubbard à deux orbitales.
  • Cette méthode calcule de manière auto-cohérente l'énergie propre et l'interaction d'appariement médiée par les fluctuations de spin, en tenant compte des dépendances complètes en impulsion et en fréquence de la fonction de Green.
  • L'équation d'Eliashberg linéarisée a été résolue pour obtenir la valeur propre $\lambda$ (un proxy pour $T_c$) et la fonction de gap supraconducteur $\Delta(k)$.

3. Contributions et résultats clés

A. Tendances structurelles et électroniques

• Angle de liaison : Lorsque la constante de réseau dans le plan $a$ diminue (augmentation de la contrainte de compression), l'angle de liaison intercouche Ni-O-Ni ($\alpha$) s'approche de $180^\circ$, ce qui est cohérent avec la transition vers la symétrie $I4/mmm$. Cependant, même pour le substrat SLAO, l'optimisation DFT n'atteint pas exactement $180^\circ$ (la symétrie $Amam$ persiste légèrement), bien que la différence d'énergie entre $Amam$ et $I4/mmm$ soit négligeable.
• Occupation orbitale : La réduction de $a$ augmente la hauteur de l'oxygène apical et le décalage de niveau ($\Delta E$) entre les orbitales $d_{x^2-y^2}$ et $d_{3z^2-r^2}$. Cela déplace les électrons vers l'orbitale $d_{3z^2-r^2}$, la rapprochant de la demi-remplissage (occupation $\approx 0,95$ pour SLAO).
• Intégrales de saut :
  • Le saut intercouche $|t_\perp|$ augmente lorsque $a$ diminue dans les structures optimisées par DFT (atteignant $\sim 0,6$ eV).
  • Crucialement, la structure déterminée expérimentalement produit un $|t_\perp|$ beaucoup plus faible, $\approx 0,45$ eV, en raison d'une distance Ni-Ni plus grande.

B. Robustesse de l'appariement en onde $s_{\pm}$

• Stabilité de la symétrie : L'étude révèle que la symétrie d'appariement en onde $s_{\pm}$ est robuste. Elle reste la solution dominante indépendamment de :
  • La symétrie cristalline spécifique ($Amam$ vs $I4/mmm$).
  • La présence ou l'absence de la poche $\gamma$ (qui dépend de la correction $U$ et du remplissage de bande).
  • Les détails spécifiques de la topologie de la surface de Fermi.
• Caractéristiques de la fonction de gap :
  • Dans la représentation orbitale, la fonction de gap est dominée par l'appariement intercouche des orbitales $d_{3z^2-r^2}$.
  • Ce gap intercouche est quasi indépendant de l'impulsion (appariement intra-maille).
  • La fonction de gap dans la représentation de bande reflète le poids orbital : elle est grande là où le caractère $d_{3z^2-r^2}$ est présent et change de signe entre les bandes de liaison et antiliantes.

C. Origine de la $T_c$ réduite dans les films minces

• Rôle de $t_\perp$ : La raison principale de la $T_c$ divisée par deux dans les films minces (par rapport au massif sous pression) est identifiée comme étant la réduction de l'intégrale de saut intercouche $|t_\perp|$.
  • Les modèles basés sur des structures optimisées par DFT (grand $|t_\perp| \approx 0,6$ eV) prédisent un $\lambda$ élevé, similaire au massif sous pression.
  • Les modèles basés sur la structure déterminée expérimentalement (petit $|t_\perp| \approx 0,45$ eV) produisent un $\lambda$ considérablement supprimé, cohérent avec la $T_c$ observée d'environ 40 K.
• Mécanisme : La suppression de $T_c$ avec la réduction de $|t_\perp|$ est attribuée à la nature de la « colle » d'appariement. L'appariement est médié par des fluctuations de spin à énergie finie (plutôt que par le nesting de surface de Fermi à basse énergie).
  • Les fluctuations de spin à haute énergie (qui pilotent l'appariement) sont supprimées lorsque $|t_\perp|$ est réduit.
  • À l'inverse, les fluctuations à basse énergie (associées au nesting de surface de Fermi) sont renforcées par la réduction de $|t_\perp|$, mais celles-ci sont moins efficaces pour l'appariement $s_{\pm}$ dans ce système.

4. Importance et conclusion

• Résolution de l'énigme de la $T_c$ : L'article fournit une explication théorique cohérente de la $T_c$ réduite dans les films minces à pression ambiante : la structure cristalline expérimentale possède intrinsèquement un saut intercouche ($t_\perp$) plus faible que celui prédit par l'optimisation DFT standard, ce qui supprime directement la force d'appariement médiée par les fluctuations de spin à énergie finie.
• Validation du mécanisme d'appariement : L'étude renforce la vision selon laquelle la supraconductivité de $La_3Ni_2O_7$ est pilotée par des fluctuations de spin à énergie finie médiatisant un appariement intercouche en $d_{3z^2-r^2}$. Ce mécanisme est robuste face aux variations de la topologie de la surface de Fermi (telles que la présence de la poche $\gamma$) et aux détails de la symétrie structurale.
• Cohérence théorique : Les résultats s'alignent avec des approches non perturbatives (cDMFT, DCA, VMC), confirmant que le système peut être considéré comme un modèle de Hubbard bicouche (formé par les orbitales $d_{3z^2-r^2}$) intégré dans un environnement multi-orbitales.
• Perspectives futures : La grande divergence entre les structures cristallines prédites par DFT et celles observées expérimentalement (spécifiquement la distance Ni-Ni) reste une question ouverte. Les auteurs suggèrent que des ordres compétitifs (par exemple, les ondes de densité) non capturés par le schéma FLEX actuel pourraient également jouer un rôle dans la stabilité et les propriétés du matériau réel.

En résumé, l'article établit que si la symétrie d'appariement ($s_{\pm}$) est robuste, l'ampleur de $T_c$ est hautement sensible au saut intercouche $t_\perp$, lequel est dicté par la structure atomique précise du film mince.