Low-temperature transport in high-conductivity correlated metals: a density-functional plus dynamical mean-field study of cubic perovskites

1. Problématique et Contexte

Le transport électronique dans les matériaux corrélés, en particulier à basse température, pose un défi majeur pour la modélisation théorique. Bien que les méthodes basées sur la théorie de la perturbation de la fonctionnelle de la densité (DFT) aient considérablement amélioré notre compréhension de la diffusion électron-phonon, les méthodes pour capturer avec précision la diffusion électron-électron (el-el) restent moins développées.

Le cas spécifique des métaux à haute conductivité et modérément corrélés (comme les oxydes de pérovskite $ABO_3$) est particulièrement difficile car :

• Les taux de diffusion à basse température sont extrêmement faibles (de l'ordre du meV).
• Cela nécessite une précision numérique exceptionnelle de la structure électronique à basse énergie.
• La résistivité $\rho(T)$ dans ces matériaux suit souvent une loi de puissance $T^2$ (comportement de liquide de Fermi), mais extraire le coefficient $A$ de manière quantitative à partir de calculs ab initio est complexe en raison des erreurs d'extrapolation et de continuation analytique.

L'objectif de ce travail est de démontrer une méthodologie robuste permettant de calculer la contribution locale de la diffusion électron-électron à la résistivité en courant continu (DC) dans ces matériaux, en utilisant une combinaison avancée de DFT et de la théorie de la fonctionnelle de densité dynamique (DMFT).

2. Méthodologie

Les auteurs ont développé un cadre computationnel sophistiqué combinant plusieurs techniques de pointe pour surmonter les limitations numériques habituelles :

• Structure de base (DFT + Wannier) :

  • Calculs DFT effectués avec VASP pour les structures cubiques ($Pm\bar{3}m$) de quatre oxydes : $SrVO_3$, $SrMoO_3$, $PbMoO_3$ et $SrRuO_3$.
  • Utilisation de fonctions de Wannier maximales localisées (MLWF) pour interpoler la structure de bande sur des maillages réciproques denses, permettant un calcul précis des intégrales de la zone de Brillouin (BZ).

• Intégration de la Zone de Brillouin (BZ) Adaptative :

  • Pour les matériaux à faible diffusion, les intégrales de BZ convergent très lentement avec les grilles uniformes. Les auteurs utilisent une méthode d'intégration adaptative itérée (IAI) (implémentée dans le package AutoBZ.jl). Cette méthode réduit la complexité de calcul et automatise la convergence vers une précision numérique spécifiée par l'utilisateur.

• Résolution du Problème d'Impureté (DMFT) :

  • Utilisation de deux solveurs d'impureté complémentaires pour valider les résultats :
    1. QMC (Quantum Monte Carlo) : Basé sur l'expansion de l'hybridation en temps imaginaire (CT-QMC). Nécessite une continuation analytique (AC) vers l'axe réel.
    2. NRG (Groupe de Renormalisation Numérique) : Opère directement sur l'axe des fréquences réelles et à température nulle, évitant l'incertitude de la continuation analytique, mais limité en complexité du système.
  • Accord "poignée de main" (handshake) : Les résultats des deux solveurs sont comparés pour garantir la fiabilité des self-énergies $\Sigma(\omega)$.

• Extraction des Taux de Diffusion (Stratégie de Mise à l'Échelle) :

  • Au lieu d'extrapoler directement les données à basse fréquence (source d'erreurs majeures), les auteurs exploitent la mise à l'échelle de Fermi-Liquide (FL).
  • Ils ajustent la partie imaginaire de la self-énergie à la forme théorique $\text{Im}\Sigma(\omega) \propto C(\omega^2 + \pi^2 T^2)$ dans une région de fréquence plus élevée où les données sont stables, puis extrapolent le coefficient $C$ vers les basses fréquences.
  • Cela permet d'obtenir le coefficient de résistivité $A$ via la relation $\rho \approx A T^2$, où $A = 24C / \Phi(\epsilon_F)$, $\Phi$ étant la fonction de transport.

3. Résultats Principaux

L'étude porte sur quatre matériaux aux propriétés de transport distinctes :

1. $SrVO_3$ et $SrMoO_3$ (Métaux à haute conductivité, $d^1$ et $d^2$) :

   • Ces matériaux présentent un comportement de liquide de Fermi clair à basse température.
   • Les calculs montrent une convergence quantitative de la contribution électronique à la résistivité.
   • Le coefficient $A$ calculé est cohérent avec les mesures expérimentales, confirmant que la diffusion électron-électron est bien décrite par la théorie FL, bien que la diffusion électron-phonon domine à température ambiante.
   • $SrMoO_3$ est confirmé comme ayant la résistivité la plus faible parmi les pérovskites, mais ce n'est pas dû à une absence de corrélations électroniques.

2. $PbMoO_3$ (Isolant ou métal à haute résistivité) :

   • Bien que formellement similaire à $SrMoO_3$ ($d^2$), il présente une résistivité trois ordres de grandeur plus élevée et une dépendance en température sous-linéaire.
   • Les calculs DFT+DMFT montrent que la contribution électronique locale seule ne suffit pas à expliquer la résistivité élevée observée expérimentalement, suggérant un rôle crucial de la structure cristalline (distorsions) ou d'autres mécanismes non locaux, bien que le modèle cubique utilisé ici montre un comportement FL similaire à $SrMoO_3$.

3. $SrRuO_3$ (Système fortement corrélé, $d^4$) :

   • Ce matériau présente un comportement de non-liquide de Fermi (non-FL) à basse température expérimentalement.
   • Les calculs dans la phase paramagnétique cubique montrent un comportement FL à très basse température ($T < 4$ K), mais avec un coefficient $A$ très élevé (3 ordres de grandeur supérieur à $SrVO_3$).
   • L'étude met en évidence la sensibilité de la transition ferromagnétique aux paramètres d'interaction (Hund's coupling $J$). Les résultats dans la phase ferromagnétique (calculés avec NRG) montrent une résistivité plus faible, en accord qualitatif avec les tendances expérimentales, bien que la comparaison quantitative soit difficile en raison de la dispersion des données expérimentales.

4. Contributions Clés et Innovations

• Précision Numérique : Démonstration que l'utilisation de l'intégration adaptative de la zone de Brillouin et de l'analyse de mise à l'échelle FL permet de surmonter les problèmes de convergence dans les métaux à haute conductivité.
• Validation Croisée : Établissement d'une procédure robuste validant les résultats QMC (avec continuation analytique) par des calculs NRG directs, validant ainsi la fiabilité des taux de diffusion extraits.
• Méthodologie de Extraction de $C$ : Proposition d'une méthode pour extraire le coefficient de diffusion $C$ en évitant les extrapolations instables à $\omega \to 0$, en utilisant plutôt la région de fréquence où la loi de mise à l'échelle FL est vérifiée.
• Analyse Comparative : Mise en évidence que la conductivité exceptionnelle de $SrMoO_3$ n'est pas due à une faiblesse des corrélations électroniques, mais implique nécessairement une interaction électron-phonon très efficace.

5. Signification et Perspectives

Ce travail fournit un outil prédictif robuste pour identifier et comprendre les matériaux à haute conductivité pour des applications technologiques (ex: interconnexions électroniques). Il démontre que la combinaison DFT+DMFT, lorsqu'elle est couplée à des techniques numériques avancées, peut fournir une compréhension fondamentale des phénomènes de transport dans les métaux corrélés.

Les auteurs soulignent que cette méthodologie ouvre la voie à l'étude d'autres propriétés de transport (effet Seebeck, conductivité de Hall, magnétorésistance) et à l'investigation future des effets de corrélations non locales et des corrections de vertex, qui sont actuellement négligés dans l'approximation DMFT à site unique mais qui pourraient être cruciaux pour des systèmes de basse dimensionnalité.