An anisotropic functional for two-dimensional material systems

1. Énoncé du problème

La théorie de la fonctionnelle de la densité (DFT) est l'outil standard pour les calculs de structure électronique, mais sa précision est limitée par le choix du potentiel d'échange-corrélation (XC).

• Le défi des matériaux 2D : Les matériaux bidimensionnels (2D) (par exemple, MoS₂, hBN) présentent des structures atomiques fortement inhomogènes dans l'espace et un écrantage diélectrique faible hors du plan. Cela conduit à un écrantage diélectrique anisotrope et non local, qui diffère fondamentalement des matériaux 3D en volume.
• Limites des méthodes actuelles :
  • Fonctionnelles semi-locales (GGA/PBE) : Elles sous-estiment les gaps de bande et présentent une convexité spurieuse dans l'énergie totale en fonction de l'occupation fractionnaire (violant le théorème de Koopmans généralisé), conduisant à une délocalisation artificielle de la charge et à des erreurs d'auto-interaction.
  • Fonctionnelles hybrides standards (PBE0, HSE) : Conçues pour les matériaux en volume en utilisant des interactions de Coulomb isotropes ou écrantées empiriquement à courte portée. Elles échouent à capturer le comportement asymptotique correct de l'écrantage diélectrique dans les systèmes 2D, résultant en des gaps de bande et des alignements de niveaux de défauts inexacts.
  • Approximation GW : Bien que précise pour les énergies de quasi-particules, GW est prohibitivement coûteuse en calcul pour les grandes supercellules (essentielles pour les études de défauts) et manque d'implémentations courantes pour les forces et les optimisations de géométrie.

2. Méthodologie

L'auteur propose une fonctionnelle d'échange écranté (SX) anisotrope auto-cohérente dans le cadre de la DFT de Kohn-Sham.

• Concept central : La fonctionnelle remplace l'interaction de Coulomb nue dans l'opérateur de Fock par une interaction écrantée statique ($W$) qui prend explicitement en compte la dépendance au vecteur d'onde ($q$) et l'anisotropie de l'écrantage diélectrique 2D.
• Formulation mathématique :
  • L'énergie d'échange écranté est définie comme :
    $$E_{SX} = -\frac{1}{2} \sum_{i,j}^{occ} \int dr \int dr' \psi_i^*(r)\psi_j(r) W(r, r') \psi_j^*(r')\psi_i(r')$$
  • L'interaction écrantée effective $V_{eff}^{2D}(q)$ pour une monocouche intégrée dans des environnements diélectriques (substrat/couverture) est dérivée d'un cadre d'écrantage macroscopique :
    $$V_{eff}^{2D}(q) = \frac{e^2 F(q)}{2\tilde{V} q \epsilon_{eff}^{2D}(q)}$$
    Où $\epsilon_{eff}^{2D}(q)$ est une fonction diélectrique dépendante du vecteur d'onde dérivée des constantes diélectriques de la couche et de son environnement.
  • Le modèle intègre une épaisseur de couche finie ($h$) et une fonction diélectrique modèle pour le matériau monocouche qui satisfait les limites asymptotiques correctes aux petits et grands $q$.
• Implémentation : La fonctionnelle a été implémentée dans une version modifiée du code VASP (cadre des ondes augmentées par projecteurs). Elle conserve l'efficacité computationnelle des calculs de fonctionnelles hybrides (comparable à PBE0) tout en incorporant la physique de l'échange non local.

3. Contributions clés

1. Première fonctionnelle anisotrope auto-cohérente : Il s'agit de la première fonctionnelle de Kohn-Sham intégrant un écrantage diélectrique anisotrope dépendant de $q$ spécifiquement pour les matériaux 2D directement dans l'opérateur d'échange.
2. Généralisation du SX en volume : Elle généralise la fonctionnelle d'échange écranté en volume réussie (Réf. 26) à une dimensionnalité réduite, retrouvant la limite du volume lorsque l'épaisseur de la couche tend vers l'infini.
3. Pont entre DFT et GW : La méthode capture la physique dominante de l'échange non local de l'énergie propre GW (spécifiquement le composant d'échange écranté statique) à une fraction du coût computationnel.
4. Accès aux forces et à la géométrie : Contrairement à GW, cette fonctionnelle fournit un accès direct aux énergies totales, aux forces atomiques et aux relaxations structurales, la rendant adaptée à l'étude des défauts et des distorsions structurales dans de grandes supercellules.

4. Résultats

La fonctionnelle a été validée par rapport à des calculs de référence GW₀ pour un ensemble diversifié de semi-conducteurs 2D (par exemple, hBN, GaSe, MoS₂, WSe₂, Cu₂Se).

• Gaps de bande de quasi-particules :
  • La fonctionnelle anisotrope reproduit les gaps de bande GW₀ avec une grande précision pour des matériaux aux gaps de bande et caractères de liaison variés.
  • Elle surpasse nettement les fonctionnelles d'échange écranté 3D isotropes, en particulier pour les isolants à large gap (hBN) et les dichalcogénures de métaux de transition (MoS₂).
  • L'épaisseur de couche effective ($h$) dérivée du modèle s'aligne avec l'étendue physique des orbitales de valence.
• Structure de bande et physique des défauts :
  • La fonctionnelle reproduit avec précision la dispersion des états aux bords de bande (maximum de la bande de valence et minimum de la bande de conduction), qui sont critiques pour les énergies de formation des défauts et la localisation.
  • Elle décrit correctement la structure électronique près des bords de bande, essentiel pour la modélisation des états de défauts.
• Linéarité par morceaux (Théorème de Koopmans) :
  • Les calculs d'énergie totale en fonction de l'occupation fractionnaire pour des défauts de substitution (par exemple, Ge dans GaSe, C dans hBN) montrent un comportement quasi linéaire par morceaux.
  • La courbure ($b_2$) est très faible (~0,03–0,08 eV), indiquant une suppression substantielle des erreurs d'auto-interaction et de la délocalisation artificielle de la charge. Cela confirme que la fonctionnelle satisfait approximativement le théorème de Koopmans généralisé.
• Propriétés optiques :
  • Les spectres d'absorption optique calculés via la TDDFT à réponse linéaire utilisant cette fonctionnelle comme noyau montrent un excellent accord avec les résultats GW₀ + Équation de Bethe-Salpeter (BSE).
  • Elle capture correctement les énergies de liaison des excitons et les forces d'oscillateur, contrairement aux approximations locales adiabatiques standard (ALDA) qui échouent à décrire les interactions électron-trou à longue portée en 2D.

5. Importance

• Pouvoir prédictif : La méthode fournit un cadre sans paramètres (sans ajustement spécifique au matériau) pour prédire les propriétés électroniques, optiques et des défauts des matériaux 2D.
• Efficacité computationnelle : Elle offre une alternative pratique à GW pour les simulations à grande échelle (par exemple, supercellules de défauts, adsorbats) où GW est actuellement irréalisable, tout en maintenant une précision comparable aux fonctionnelles hybrides.
• Cadre unifié : Elle permet des simulations cohérentes, de l'optimisation de la géométrie de l'état fondamental aux calculs de réponse optique, en utilisant une seule fonctionnelle motivée physiquement.
• Perspectives futures : Bien que actuellement centrée sur les semi-conducteurs non magnétiques, le cadre établit une base pour étendre l'écrantage anisotrope aux systèmes 2D magnétiques et aux hétérostructures complexes.

En résumé, ce travail résout une lacune critique dans les simulations ab initio des matériaux 2D en introduisant une fonctionnelle auto-cohérente et efficace en calcul qui modélise correctement l'écrantage diélectrique anisotrope unique de ces systèmes, permettant ainsi des prédictions précises des gaps de bande, des comportements des défauts et des propriétés optiques.