Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

En utilisant des calculs DFT+DMFT auto-cohérents en charge, cette étude examine comment la structure cristalline, les champs de ligands et les forces de corrélation variables (y compris les corrélations 2p de l'oxygène) influencent la structure électronique et les fonctions spectrales du NiO et du CoO paramagnétiques.

L'essentiel

Imaginez un monde composé de petits danseurs énergiques (les électrons) se déplaçant sur une scène. Dans la plupart des matériaux, ces danseurs se déplacent librement, comme une foule à un concert. Mais dans des matériaux spéciaux appelés oxydes de métaux de transition (comme l'oxyde de nickel et l'oxyde de cobalt), les danseurs sont si nombreux et si sensibles les uns aux autres qu'ils commencent à agir comme un groupe soudé. Si l'un bouge, tous les autres réagissent instantanément. C'est ce qu'on appelle la « forte corrélation », et cela fait que ces matériaux se comportent de manière délicate et fascinante — agissant parfois comme des isolants (bloquant l'électricité) au lieu de conducteurs.

L'article que vous avez fourni est une plongée approfondie pour comprendre exactement comment ces danseurs se comportent dans deux matériaux spécifiques : l'oxyde de nickel (NiO) et l'oxyde de cobalt (CoO). Les chercheurs voulaient déterminer pourquoi ces matériaux agissent ainsi et comment prédire leur comportement avec précision en utilisant des simulations informatiques.

Voici une décomposition de leurs découvertes utilisant des analogies simples :

1. Le Problème : La « Mauvaise Carte »

Les scientifiques ont une méthode standard pour cartographier ces matériaux appelée DFT (Théorie de la Fonctionnelle de la Densité). Considérez la DFT comme un GPS qui fonctionne généralement très bien pour conduire sur des autoroutes dégagées. Cependant, pour ces matériaux « dansants » et bondés, le GPS est cassé. Il prédit que les danseurs peuvent se déplacer librement (métallique), mais en réalité, ils sont bloqués sur place (isolant). La carte est fausse car elle ignore le fait que les danseurs se cognent constamment les uns contre les autres et réagissent mutuellement.

2. La Solution : Un Meilleur GPS avec « Auto-Correction »

Pour corriger cela, les auteurs ont utilisé une méthode ultra-avancée appelée DFT+DMFT.

• DMFT (Théorie du Champ Moyen Dynamique) : Imaginez ajouter un « chorégraphe » à la simulation qui observe chaque danseur individuellement et prend en compte comment ils réagissent à leurs voisins en temps réel. Cela résout le problème « métal contre isolant ».
• SIC (Correction d'Auto-Interaction) : Les chercheurs ont réalisé que les danseurs ne sont pas seulement les atomes métalliques (Nickel ou Cobalt) ; les atomes d'oxygène qui les entourent font aussi partie de la danse. Ils ont ajouté une règle spéciale d'« auto-correction » pour les danseurs d'oxygène afin de s'assurer que la simulation ne comptait pas deux fois leur énergie. Cette combinaison (DFT+DMFT+SIC) a créé une carte qui correspondait presque parfaitement aux expériences réelles.

3. Les Deux Pistes de Danse : Octaédrique vs Tétraédrique

Les chercheurs ont testé ces matériaux sur deux « pistes de danse » différentes (structures cristallines) :

• Sel gemme (RS) : Les danseurs sont disposés dans un octaèdre (comme une boule à 6 pointes). C'est la forme stable et naturelle de ces matériaux.
• Blende de zinc (ZB) : Les danseurs sont disposés dans un tétraèdre (comme une pyramide à 4 pointes). C'est une forme instable et artificielle que les chercheurs ont créée juste pour voir comment les danseurs réagiraient à un environnement différent.

La Découverte : Tout comme un danseur peut se sentir différent dans un cercle par rapport à un carré, les électrons divisent leurs niveaux d'énergie différemment selon la forme de la piste de danse.

• Sur la piste Octaédrique (RS), les niveaux d'énergie se divisent d'une certaine manière.
• Sur la piste Tétraédrique (ZB), les niveaux s'inversent.
  L'étude a confirmé que la piste Octaédrique est plus stable car les danseurs (électrons) peuvent s'installer dans des endroits à plus basse énergie, plus confortables. La piste Tétraédrique force trop de danseurs dans des endroits « anti-liants » (inconfortables, à haute énergie), rendant la structure instable. Cela explique pourquoi vous trouvez de l'oxyde de nickel naturellement sous forme octaédrique mais pas sous forme tétraédrique.

4. Le « Satellite » et le « Gap »

L'un des objectifs principaux était de mesurer le gap de bande — le saut d'énergie requis pour mettre un danseur en mouvement.

• L'Expérience : Les expériences réelles ont montré un gap d'environ 5 à 6 électron-volts (eV).
• L'Ancienne Simulation : Sans les corrections spéciales, l'ordinateur prédisait un minuscule gap (environ 2,5–3 eV), ce qui était faux.
• La Nouvelle Simulation : En utilisant l'« Auto-Correction » (SIC) pour les atomes d'oxygène, le modèle informatique des chercheurs a prédit un gap de 5,1 eV. C'est une correspondance parfaite avec le monde réel !

Ils ont également examiné un « pic satellite » (une bosse de haute énergie dans les données). Ils ont constaté que, bien que leur modèle ait amélioré le gap principal, le pic satellite restait un peu difficile à définir parfaitement, mais l'image globale était beaucoup plus claire qu'auparavant.

5. Nickel vs Cobalt : La Différence « D'un Pas »

Le nickel et le cobalt sont voisins dans le tableau périodique. Le nickel possède un électron de plus que le cobalt.

• Oxyde de Nickel : Possède un nombre spécifique de danseurs qui crée un état « Zhang-Rice » très stable et à haute énergie (une paire spéciale de danseurs fortement liés).
• Oxyde de Cobalt : Possède un danseur de moins. Cela crée un « trou » (un espace vide). Les chercheurs ont découvert que ce danseur manquant fait que le matériau au cobalt se comporte légèrement différemment, créant un pic plus net et plus intense juste au bord de la bande d'énergie. C'est comme retirer une personne d'une pièce bondée ; les personnes restantes se déplacent et réagissent plus intensément à l'espace vide.

Résumé

En bref, cet article traite de la construction d'une simulation informatique parfaitement précise de la façon dont les électrons dansent dans les oxydes de nickel et de cobalt.

1. Ils ont prouvé que vous devez prendre en compte les atomes d'oxygène (pas seulement le métal) pour obtenir la bonne réponse.
2. Ils ont montré que la forme du cristal (Octaédrique vs Tétraédrique) modifie la façon dont les électrons divisent leur énergie, expliquant pourquoi certaines formes sont stables et d'autres non.
3. Ils ont obtenu une correspondance quasi parfaite avec les expériences réelles, prouvant que leur méthode est un outil fiable pour comprendre ces matériaux complexes.

Les auteurs concluent que cette méthode est un outil puissant pour les scientifiques qui souhaitent concevoir de meilleurs matériaux pour des applications telles que les catalyseurs (pour décomposer l'eau ou produire du carburant) et le stockage d'énergie, car comprendre la « danse » des électrons est la clé pour débloquer leur potentiel.

Résumé technique

Résumé technique : Influence du champ de ligands et des corrélations sur la structure électronique de NiO et CoO à partir de calculs DFT+DMFT

Énoncé du problème
Les oxydes de métaux de transition (OMT) tels que NiO et CoO présentent des comportements électroniques complexes pilotés par l'interplay des corrélations électroniques, des champs de ligands, de la structure cristalline et de l'hybridation. Alors que la théorie de la fonctionnelle de la densité (DFT) standard, utilisant les approximations du gradient local ou généralisé (LDA/GGA), échoue à reproduire les états fondamentaux isolants et les spectres d'excitation corrects de ces systèmes fortement corrélés, des méthodes plus avancées sont requises. Plus précisément, il est nécessaire de comprendre systématiquement comment la structure cristalline (type sel de roche vs blende de zinc), la composition (Ni vs Co) et le traitement des corrélations dans les orbitales $3d$ du métal de transition (MT) et $2p$ de l'oxygène influencent les fonctions spectrales et les largeurs de bande interdite. Des études antérieures se sont souvent concentrées sur des paramètres d'interaction spécifiques ou ont négligé les corrélations dans les ligands d'oxygène, qui sont cruciaux pour décrire avec précision les isolants à transfert de charge.

Méthodologie
Les auteurs emploient une combinaison auto-cohérente en charge de la théorie de la fonctionnelle de la densité et de la théorie du champ moyen dynamique (DFT+DMFT), augmentée d'un schéma de correction d'auto-interaction (SIC) pour les orbitales de l'oxygène (DFT+sicDMFT).

• Cadre de calcul : Le composant DFT utilise le code de pseudopotentiels à base mixte (MBPP) avec des pseudopotentiels à conservation de norme. Le sous-espace corrélé est défini à l'aide d'orbitales de type atomique $d$ centrées sur les ions MT, décomposées dans une base de Kohn-Sham tronquée.
• Traitement des corrélations : Le problème d'impureté DMFT est résolu à l'aide du solveur Monte Carlo quantique en temps continu (QMC) CTHYB (paquet TRIQS). L'interaction électron-électron au sein du sous-espace $3d$ du MT est modélisée par un Hamiltonien de Slater invariant par rotation avec un paramètre de Hubbard $U$ et un échange de Hund $J$.
  • Pour NiO : $U = 10$ eV, $J = 1$ eV.
  • Pour CoO : $U = 9$ eV (testé également avec $U=6$ eV), $J = 1$ eV (testé également avec $J=0.9$ eV).
• Corrélations des ligands : Pour rendre compte des corrélations dans les orbitales O($2p$), qui sont hybridées avec les états du MT, les auteurs appliquent un schéma SIC aux pseudopotentiels de l'oxygène. Les paramètres sont fixés à $\alpha=0.8$, avec une correction complète ($w_l=1.0$) pour O($2s$) et une correction à 80 % ($w_l=0.8$) pour O($2p$).
• Systèmes étudiés : Le NiO et le CoO paramagnétiques sont analysés dans les structures cubiques sel de roche (RS) et blende de zinc (ZB). Alors que RS est la phase stable, ZB est étudié comme système modèle pour isoler les effets du champ de ligands (coordination octaédrique vs tétraédrique).
• Analyse : Les fonctions spectrales en fréquence réelle sont dérivées par continuation analytique (méthode du maximum d'entropie). L'étude se concentre sur les fonctions spectrales totales, projetées par élément (MT et O) et résolues par orbitale locale, en analysant spécifiquement les centres de gravité (COG) et les caractéristiques des bandes interdites.

Résultats clés

1. Rôle du SIC et de l'intensité de l'interaction ($U$) :

   • Pour RS-CoO, l'augmentation de $U$ de 6 eV à 9 eV déplace le pic de basse énergie $d^7\underline{L}$ (où $\underline{L}$ désigne un trou de ligand) de $-1.0$ eV à $-1.6$ eV.
   • L'application du O($2p$)-SIC déplace ce pic dans la direction opposée vers $-0.6$ eV. La combinaison de $U=9$ eV et du O($2p$)-SIC produit une forme spectrale et une position de pic ($-3.1$ eV pour la bande O($2p$)) qui s'alignent bien avec les données expérimentales UPS.
   • Crucialement, une large bande interdite significative ne s'ouvre que pour des valeurs de $U$ suffisamment grandes. L'inclusion du O($2p$)-SIC déplace la bande de Hubbard supérieure inoccupée (UHB) vers des énergies plus élevées, résultant en une bande interdite calculée de 5,1 eV pour RS-CoO. Cette valeur tombe dans la plage expérimentale de 5–6 eV, constituant une amélioration marquée par rapport aux études DFT+DMFT précédentes sans O($2p$)-SIC, qui rapportaient des largeurs de bande aussi faibles que 2,5–3,8 eV.

2. Effets structuraux et du champ de ligands (RS vs ZB) :

   • Les calculs confirment les prédictions de la théorie du champ cristallin (CF). Dans RS (octaédrique), les états $t_{2g}$ sont plus bas en énergie que les $e_g$. Dans ZB (tétraédrique), la hiérarchie s'inverse ($e$ plus bas que $t_2$).
   • L'analyse de l'occupation via les fonctions spectrales intégrées confirme des configurations à haut spin : $t_{2g}^6 e_g^2$ pour RS-NiO et $t_{2g}^5 e_g^2$ pour RS-CoO. Pour les structures ZB, les occupations sont $e^4 t_2^4$ (NiO) et $e^4 t_2^3$ (CoO).
   • Les structures ZB sont trouvées moins stables en raison de la population d'états antiliants. Dans ZB, les orbitales $t_2$ sont antiliantes ; la présence de 3 (CoO) et 4 (NiO) électrons dans ces états déstabilise la structure par rapport à la phase RS.

3. Caractéristiques spectrales et écrantage :

   • Les spectres de CoO présentent des pics plus nets et un poids spectral accru à basse énergie par rapport à NiO, attribué au trou supplémentaire dans la coquille $d$.
   • Un pic net près du bord de la bande de valence dans CoO est interprété comme un état $d^7\underline{Z}$ (où $\underline{Z}$ désigne un trou de Zhang-Rice), analogue au comportement $d^8\underline{Z}$ discuté dans le NiO dopé au Li.
   • L'énergie de transfert de charge ($\Delta$), définie comme la différence d'énergie entre les états O($2p$) occupés et les états MT($3d$) inoccupés, est calculée entre 4,3 et 4,9 eV selon la structure et l'élément, montrant une hiérarchie cohérente avec les tendances expérimentales (NiO < CoO).

Signification et affirmations
L'article affirme que l'approche DFT+sicDMFT fournit un cadre robuste pour décrire la structure électronique des OMT fortement corrélés. En incluant explicitement les effets de corrélation dans les orbitales $2p$ de l'oxygène via le SIC, la méthode reproduit avec succès les spectres d'excitation expérimentaux et les valeurs de bande interdite pour le NiO et le CoO paramagnétiques, que les approches de champ moyen statique (comme DFT+U) ou la DFT+DMFT standard échouent souvent à capturer avec précision.

L'étude démontre que :

• L'interplay entre le paramètre de Hubbard $U$ et les corrélations des ligands est crucial pour un accord quantitatif avec l'expérience.
• La méthodologie peut systématiquement démêler les effets de la structure cristalline (champ de ligands) et de la composition sur les fonctions spectrales.
• Les caractéristiques spectrales calculées, telles que les largeurs de bande interdite et les énergies de transfert de charge, servent de descripteurs fiables qui pourraient soutenir l'optimisation des OMT pour des applications en catalyse et conversion d'énergie, à condition que la structure électronique soit modélisée avec précision.

Les auteurs adoptent une stance modeste, notant que bien que le ZB-NiO ne soit pas une phase stable en volume, son étude en tant que système modèle permet une comparaison systématique des effets du champ de ligands. Le travail s'appuie sur des résultats antérieurs pour NiO afin de valider l'approche pour CoO et d'étendre l'analyse aux variations structurelles.