An anisotropic functional for two-dimensional material systems

1. Enunciato del Problema

La Teoria del Funzionale Densità (DFT) è lo strumento standard per i calcoli della struttura elettronica, ma la sua accuratezza è limitata dalla scelta del potenziale di scambio-correlazione (XC).

• La Sfida dei Materiali 2D: I materiali bidimensionali (2D) (ad es. MoS₂, hBN) presentano strutture atomiche fortemente disomogenee nello spazio e uno screening dielettrico debole nel piano fuori asse. Ciò porta a uno screening dielettrico anisotropo e non locale, che differisce fondamentalmente dai materiali massivi 3D.
• Limiti dei Metodi Attuali:
  • Funzionali Semilocali (GGA/PBE): Sottostimano i gap di banda e mostrano una convessità spuria nell'energia totale in funzione dell'occupazione frazionaria (violando il teorema generalizzato di Koopmans), portando a una delocalizzazione artificiale della carica e a errori di auto-interazione.
  • Funzionali Ibridi Standard (PBE0, HSE): Progettati per materiali massivi utilizzando interazioni di Coulomb isotrope o schermate empiricamente a corto raggio. Non riescono a catturare il corretto comportamento asintotico dello screening dielettrico nei sistemi 2D, risultando in gap di banda inaccurati e allineamenti dei livelli di difetto errati.
  • Approssimazione GW: Sebbene accurata per le energie dei quasi-particelle, la GW è computazionalmente proibitiva per supercelle grandi (essenziali per gli studi sui difetti) e manca di implementazioni routinarie per le forze e le ottimizzazioni geometriche.

2. Metodologia

L'autore propone un funzionale di scambio schermato (SX) anisotropo autoconsistente all'interno del quadro della DFT di Kohn-Sham.

• Concetto Chiave: Il funzionale sostituisce l'interazione di Coulomb nuda nell'operatore di Fock con un'interazione schermata staticamente ($W$) che tiene esplicitamente conto della dipendenza dal vettore d'onda ($q$) e dell'anisotropia dello screening dielettrico 2D.
• Formulazione Matematica:
  • L'energia di scambio schermato è definita come:
    $$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'interazione schermata efficace $V_{eff}^{2D}(q)$ per un monostrato incorporato in ambienti dielettrici (substrato/copertura) è derivata da un quadro di screening macroscopico:
    $$V_{eff}^{2D}(q) = \frac{e^2 F(q)}{2\tilde{V} q \epsilon_{eff}^{2D}(q)}$$
    Dove $\epsilon_{eff}^{2D}(q)$ è una funzione dielettrica dipendente dal vettore d'onda derivata dalle costanti dielettriche dello strato e del suo ambiente.
  • Il modello incorpora uno spessore finito dello strato ($h$) e una funzione dielettrica modello per il materiale a monostrato che soddisfa i corretti limiti asintotici per $q$ piccoli e grandi.
• Implementazione: Il funzionale è stato implementato in una versione modificata del codice VASP (framework Projector-Augmented-Wave). Mantiene l'efficienza computazionale dei calcoli con funzionali ibridi (paragonabile a PBE0) incorporando al contempo la fisica dello scambio non locale.

3. Contributi Chiave

1. Primo Funzionale Anisotropo Autoconsistente: Questo è il primo funzionale di Kohn-Sham a incorporare uno screening dielettrico anisotropo, dipendente da $q$, specificamente per i materiali 2D direttamente nell'operatore di scambio.
2. Generalizzazione dello SX Massivo: Generalizza il funzionale di scambio schermato massivo di successo (Ref. 26) alla dimensionalità ridotta, recuperando il limite massivo quando lo spessore dello strato tende all'infinito.
3. Ponte tra DFT e GW: Il metodo cattura la fisica dominante dello scambio non locale dell'autoenergia GW (in particolare la componente di scambio schermato statico) a una frazione del costo computazionale.
4. Accesso a Forze e Geometria: A differenza della GW, questo funzionale fornisce l'accesso diretto alle energie totali, alle forze atomiche e alle rilassamenti strutturali, rendendolo adatto allo studio di difetti e distorsioni strutturali in grandi supercelle.

4. Risultati

Il funzionale è stato validato contro calcoli di riferimento GW₀ per un insieme diversificato di semiconduttori 2D (ad es. hBN, GaSe, MoS₂, WSe₂, Cu₂Se).

• Gap di Banda dei Quasi-Particelle:
  • Il funzionale anisotropo riproduce i gap di banda GW₀ con alta accuratezza attraverso materiali con gap di banda e caratteri di legame variabili.
  • Supera significativamente i funzionali di scambio schermato isotropi 3D, in particolare per isolanti a gap largo (hBN) e dicalcogenuri di metalli di transizione (MoS₂).
  • Lo spessore efficace dello strato ($h$) derivato dal modello si allinea con l'estensione fisica degli orbitali di valenza.
• Struttura a Bande e Fisica dei Difetti:
  • Il funzionale riproduce accuratamente la dispersione degli stati ai bordi di banda (massimo della banda di valenza e minimo della banda di conduzione), che sono critici per le energie di formazione dei difetti e la localizzazione.
  • Descrive correttamente la struttura elettronica vicino ai bordi di banda, essenziale per la modellazione degli stati di difetto.
• Linearità a Pezzi (Teorema di Koopmans):
  • I calcoli dell'energia totale in funzione dell'occupazione frazionaria per difetti sostituzionali (ad es. Ge in GaSe, C in hBN) mostrano un comportamento quasi lineare a pezzi.
  • La curvatura ($b_2$) è molto piccola (~0.03–0.08 eV), indicando una sostanziale soppressione degli errori di auto-interazione e della delocalizzazione artificiale della carica. Ciò conferma che il funzionale soddisfa approssimativamente il teorema generalizzato di Koopmans.
• Proprietà Ottiche:
  • Gli spettri di assorbimento ottico calcolati tramite TDDFT a risposta lineare utilizzando questo funzionale come kernel mostrano un eccellente accordo con i risultati GW₀ + Equazione di Bethe-Salpeter (BSE).
  • Cattura correttamente le energie di legame degli eccitoni e le forze degli oscillatori, a differenza delle approssimazioni locali adiabatiche standard (ALDA) che non riescono a descrivere le interazioni elettrone-lacuna a lungo raggio in 2D.

5. Significato

• Potere Predittivo: Il metodo fornisce un quadro senza parametri (nessun adattamento specifico del materiale) per prevedere le proprietà elettroniche, ottiche e dei difetti dei materiali 2D.
• Efficienza Computazionale: Offre un'alternativa pratica alla GW per simulazioni su larga scala (ad es. supercelle di difetti, adsorbati) dove la GW è attualmente non fattibile, mantenendo un'accuratezza paragonabile ai funzionali ibridi.
• Quadro Unificato: Permette simulazioni coerenti dall'ottimizzazione della geometria dello stato fondamentale ai calcoli della risposta ottica utilizzando un singolo funzionale motivato fisicamente.
• Prospettive Future: Sebbene attualmente focalizzato su semiconduttori non magnetici, il quadro stabilisce una base per estendere lo screening anisotropo a sistemi 2D magnetici e eterostrutture complesse.

In sintesi, questo lavoro risolve una lacuna critica nelle simulazioni di prima principio dei materiali 2D introducendo un funzionale autoconsistente ed efficiente dal punto di vista computazionale che modella correttamente lo screening dielettrico anisotropo unico di questi sistemi, consentendo così previsioni accurate di gap di banda, comportamenti dei difetti e proprietà ottiche.