In search of the electron-phonon contribution to total energy

1. Problématique et Contexte

L'énergie totale est une grandeur fondamentale en physique du solide et en chimie quantique, déterminant la stabilité des matériaux, les transitions de phase et la formation de défauts. La méthode standard, la théorie de la fonctionnelle de la densité (DFT), repose sur l'approximation de Born-Oppenheimer (BO), qui découple les degrés de liberté électroniques et nucléaires.

Dans cette approximation, l'énergie totale est généralement calculée comme la somme de l'énergie électronique statique (noyaux fixes) et de la contribution vibrationnelle (phonons) obtenue via la statistique de Bose-Einstein ou l'approximation quasi-harmonique. Cependant, cette séparation est approximative. L'article s'interroge sur l'existence et la nature des termes d'ordre supérieur dans le développement perturbatif par rapport à la masse nucléaire ($M_0$), spécifiquement ceux qui ne relèvent ni de l'énergie électronique pure ni de la dynamique phononique standard.

Un point de confusion majeur réside dans l'interprétation récente (Allen, 2020) suggérant que la somme des renormalisations des valeurs propres (zero-point renormalization) des états occupés, calculée via la théorie Allen-Heine-Cardona (AHC), constituerait une contribution « électron-phonon » à l'énergie totale. Les auteurs remettent en cause cette interprétation.

2. Méthodologie

Les auteurs développent une analyse rigoureuse basée sur une expansion perturbative du Hamiltonien total en puissances du paramètre petit $\lambda \propto M_0^{-1/4}$, où $M_0$ est une masse nucléaire caractéristique.

• Développement formel : Ils partent de l'équation de Schrödinger exacte et utilisent une factorisation exacte de la fonction d'onde (exact factorization) pour dériver une équation effective pour la dynamique nucléaire. Cette équation inclut un opérateur de self-énergie ($\hat{\Sigma}$) qui condense tous les degrés de liberté électroniques.
• Analyse de l'échelle de masse : En développant les termes jusqu'à l'ordre 6 en $\lambda$, ils classent les contributions énergétiques selon leur dépendance en masse ($M_0^{-n/2}$).
  • Ordre 0 : Énergie électronique BO.
  • Ordre 2 ($M_0^{-1/2}$) : Contribution phononique harmonique.
  • Ordre 4 ($M_0^{-1}$) : Contributions anharmoniques, potentiel vectoriel (Berry), et contribution électron-phonon ($E_{elph}$).
  • Ordre 6 ($M_0^{-3/2}$) : Renormalisation de la masse électronique (inertie électronique), termes mixtes, etc.
• Clarification théorique : Ils démontrent mathématiquement que la quantité proposée par Allen ($E_{Allen}$), qui est la somme des renormalisations des valeurs propres, correspond en réalité à une partie de la contribution phononique d'ordre 2 (énergie potentielle des phonons) et non à une nouvelle contribution d'ordre 4.
• Implémentation numérique : Une nouvelle formule pour la contribution électron-phonon d'ordre 4 ($E_{elph}$) est dérivée et implémentée dans le code de calcul Abinit. Cette formule fait intervenir les éléments de matrice électron-phonon et une somme sur les états non occupés (ou une formulation sans somme sur les états via la méthode de Sternheimer).
• Validation : Les résultats sont validés par comparaison avec des calculs de différences finies (phonons gelés) et par des tests de cohérence de taille (size consistency) entre maille primitive et supermaille.

3. Contributions Clés

1. Formulation exacte de l'énergie totale : Dérivation d'une équation exacte (Eq. 31) pour l'énergie totale dans la base des fonctions d'onde BO, incluant les corrections au-delà de l'approximation BO standard via un opérateur de self-énergie.
2. Classification des termes : Identification et classification complète de tous les termes jusqu'à l'ordre 6 dans l'expansion en masse, distinguant clairement les contributions anharmoniques, vectorielles (Berry), et électron-phonon.
3. Correction de l'interprétation d'Allen : Démonstration que la somme des renormalisations des valeurs propres (AHC) n'est pas une contribution séparée à l'énergie totale, mais fait partie intégrante de l'énergie phononique d'ordre 2. Ajouter cette quantité à l'énergie totale reviendrait à un double comptage.
4. Définition et calcul de $E_{elph}$ : Identification de la véritable contribution électron-phonon à l'énergie totale, apparaissant au quatrième ordre ($M_0^{-1}$). Cette contribution est liée à l'énergie de Born et aux dérivées des fonctions d'onde électroniques par rapport aux déplacements nucléaires.
5. Masse inertielle électronique : Estimation de la contribution de la renormalisation de la masse des électrons de cœur (ordre 6), montrant qu'elle est négligeable pour la stabilité relative des phases.

4. Résultats

Les auteurs appliquent leur méthode à deux allotropes du carbone : le diamant cubique (stable) et le diamant hexagonal (lonsdaléite, métastable).

• Validation à $\Gamma$ : Pour le diamant, les calculs de la renormalisation des valeurs propres et des fréquences phonons par la théorie des perturbations (DFPT) concordent avec les différences finies à moins de 0,5 %.
• Cohérence de taille : Les résultats sont indépendants de la taille de la maille (primitive vs supermaille), validant la formulation théorique.
• Valeurs énergétiques :
  • La contribution phononique d'ordre 2 ($E_{ph}$) est d'environ 358 meV par cellule (2 atomes).
  • La contribution électron-phonon d'ordre 4 ($E_{elph}$) est calculée pour la première fois : elle vaut 3,64 meV pour le diamant et 7,61 meV pour la lonsdaléite (par cellule de 2 atomes).
  • La contribution de la masse inertielle électronique ($E_{elm}$) est très faible (-0,05 meV).
• Stabilité relative :
  • L'énergie BO pure indique que le diamant est plus stable que la lonsdaléite de ~54 meV/(2 atomes).
  • L'inclusion des corrections phononiques et électron-phonon réduit cet écart à ~52,7 meV/(2 atomes). Bien que la correction soit petite, elle est non négligeable et déplace légèrement la stabilité relative.
• Expansion du réseau : L'inclusion de ces termes modifie les paramètres de réseau à 0 K (expansion du réseau due au point zéro), augmentant le paramètre de maille du diamant de 0,4 %.

5. Signification et Impact

Ce travail apporte une clarification théorique cruciale pour le calcul précis des énergies totales en physique du solide :

• Correction d'erreur conceptuelle : Il empêche l'utilisation incorrecte des renormalisations de bande (AHC) comme une contribution supplémentaire à l'énergie totale, évitant ainsi des erreurs systématiques dans les prédictions de stabilité de phase.
• Précision accrue : Il fournit une méthode rigoureuse pour inclure la contribution électron-phonon d'ordre 4, essentielle pour les systèmes où les différences d'énergie sont infimes (polymorphes, défauts, configurations magnétiques, réactions chimiques).
• Outil pratique : L'implémentation dans Abinit rend cette correction accessible aux chercheurs, permettant des prédictions de structures cristallines plus proches des mesures expérimentales à basse température.
• Fondement théorique : La dérivation complète jusqu'à l'ordre 6 et l'analyse de l'échelle de masse offrent un cadre solide pour les développements futurs au-delà de l'approximation de Born-Oppenheimer.

En résumé, l'article établit que si la contribution électron-phonon directe à l'énergie totale est petite (de l'ordre de quelques meV par atome), elle est physiquement réelle, distincte des effets phononiques standards, et nécessaire pour une description quantique complète et précise de la stabilité des matériaux.