Reaching the thermodynamic limit of periodic CCSD cohesive energies and band gaps

Technische Samenvatting: Het bereiken van de thermodynamische limiet van periodieke CCSD-cohesie-energieën en bandkloven

Probleemstelling
Periodieke gekoppelde-cluster-theorie met enkelvoudige en dubbele excitaties (CCSD) is een niet-perturbatieve, grootte-extensieve methode die in staat is tot hoogwaardige voorspellingen van de elektronische structuren van vaste stoffen. De toepassing ervan wordt echter ernstig beperkt door de computationele kosten. De opslagvereisten schalen als $O(n^4 N_k^3)$ en de executietijd als $O(n^6 N_k^4)$, waarbij $n$ het aantal atomen per eenheidscel is en $N_k$ het aantal $k$-punten. Als gevolg hiervan waren eerdere studies doorgaans beperkt tot $k$-puntsroosters die niet groter waren dan $4^3$. Deze beperking voorkomt een betrouwbare eliminatie van eindige-omvangfouten (finite-size errors), die langzamer afnemen voor bandkloven ($N_k^{-1/3}$) vergeleken met grondtoestandsenergieën ($N_k^{-1}$). Zonder toegang tot voldoende dichte $k$-puntssampling zijn extrapolaties naar de thermodynamische limiet (TDL) onbetrouwbaar, met name voor indirecte bandkloven en materialen waarbij voorspellingen met kleine roosters buiten het asymptotische convergentieregime vallen.

Methodologie
De auteurs hebben een nieuwe software-implementatie met gedistribueerd geheugen ontwikkeld en toegepast voor periodieke CCSD en equation-of-motion CCSD (EOM-CCSD), gebaseerd op het PySCF-framework. Belangrijke methodologische kenmerken zijn:

• Parallellisatie: De code maakt gebruik van ruimtegroep-symmetrie en maakt gebruik van parallellisatie met zowel gedistribueerd als gedeeld geheugen, waardoor deze efficiënt draait op maximaal 12 nodes (96 kernen per node, 1,5 TB geheugen per node).
• Sampling-dichtheid: Deze implementatie maakt berekeningen mogelijk met tot $6^3 = 216$ $k$-punten in een double-zeta (DZ) basis, $5^3$ in triple-zeta (TZ), en $4^3$ in quadruple-zeta (QZ) bases.
• Extrapolatieschema's: De studie evalueert meerdere extrapolatiemodellen voor TDL-convergentie:
  • Model A: Single-term ($N_k^{-1/3}$).
  • Model AB: Two-term ($N_k^{-1/3} + N_k^{-2/3}$).
  • Model ABC: Three-term ($N_k^{-1/3} + N_k^{-2/3} + N_k^{-1}$).
  • GW-EOM Schema: Een methode die in eerdere literatuur is voorgesteld (Ref. 19), die uitgaat van proportionaliteit tussen de expansiecoëfficiënten van EOM-CCSD en $G_0W_0@HF$ bandkloven om convergentie te versnellen.
• Basissets en Potentialen: Berekeningen maken gebruik van correlation-consistent basissets (cc-pVXZ) en Gaussian-type orbitalen (GTH) pseudopotentialen die geoptimaliseerd zijn voor Hartree-Fock-theorie. Zowel all-electron (AE) met frozen-core (FC) benaderingen als pseudopotential (PP) berekeningen zijn uitgevoerd.
• Indirecte Kloof Strategie: Voor indirecte bandkloven stellen de auteurs een samengestelde aanpak voor. In plaats van de indirecte kloof direct te extrapoleren (wat gevoelig is voor $k$-mesh verschuivingen), extrapoleren zij de directe kloof bij de valentiebandmaximum (VBM) en de geleidingsband offset tussen de VBM en de geleidingsbandminimum (CBM) afzonderlijk, en tellen deze vervolgens op.

Belangrijkste Bijdragen en Resultaten

1. Software-implementatie: De auteurs rapporteren de grootste periodieke CCSD-berekeningen tot nu toe gebaseerd op canonieke (gedelokaliseerde) orbitalen, waarmee betrouwbare TDL-extrapolatie is bereikt voor acht eenvoudige halfgeleiders en isolatoren (MgO, LiCl, LiF, LiH, BN, BP, Si, C).
2. Cohesie-energieën:
   • Cohesie-energieën werden geconvergeerd binnen 0,1 eV.
   • CCSD-resultaten vertonen gemiddelde fouten van 0,1–0,2 eV vergeleken met experimentele waarden gecorrigeerd voor nulpuntsenergie, waarbij ze de cohesie-energie doorgaans onderschatten.
   • Pseudopotential-fouten bleken klein (<0,05 eV voor de meeste vaste stoffen) en worden primair bepaald op Hartree-Fock niveau.
   • De CCSD-nauwkeurigheid is vergelijkbaar met of beter dan gangbare dichtheidsfunctionalen (PBE, PBEsol, SCAN).
3. Bandkloven en Extrapolatie-analyse:
   • Convergentie: Bandkloven zijn moeilijker te convergeren dan cohesie-energieën. De studie stelt vast dat $N_k=2^3$ onvoldoende is om hogere-orde termen in extrapolatiemodellen te negeren, terwijl $N_k=3^3$ over het algemeen wel voldoende is.
   • Prestaties van Modellen: Het ABC-extrapolatiemodel met $k$-meshes van $3^3$ tot $6^3$ dient als de geconvergeerde referentie. De single-term (Model A) extrapolatie onderschat de bandkloven consequent met 0,2–0,5 eV wanneer deze beperkt is tot kleinere meshes.
   • GW-EOM Schema: Het GW-EOM-234 schema (het fitten van coëfficiënten van $2^3$–$4^3$ meshes) levert erratische resultaten op. Het presteert goed voor sommige materialen (MgO, LiCl, LiF), maar overschat de kloven aanzienlijk voor andere (BN, BP, Si, C) wanneer de proportionaliteitsaanname tussen GW en CCSD faalt.
   • Indirecte Kloven: De samengestelde aanpak (directe kloof + offset) mitigeert succesvol de gevoeligheid van indirecte kloven voor $k$-mesh verschuivingen, wat geconvergeerde waarden oplevert met onzekerheden rond de 0,1 eV.
4. Benchmark Waarden: Het artikel levert definitieve benchmark-getallen voor CCSD cohesie-energieën en bandkloven.
   • Nauwkeurigheid van de Bandkloof: EOM-CCSD voorspelt bandkloven met een gemiddelde absolute fout (MAE) van ~0,4 eV vergeleken met experimentele waarden gecorrigeerd voor zero-point renormalization (ZPR). De theorie overschat de kloof doorgaans.
   • Uitschieters: Significante afwijkingen (0,8–1,0 eV) worden waargenomen voor rochtselt-structuren (MgO, LiH, LiF, LiCl). Analyse van het single-excitation karakter ($n_{1}^{IP}$) suggereert dat deze materialen een lagere single-excitation gewichten hebben (93–94%) vergeleken met anderen (>95%), wat aangeeft dat de inclusie van triple excitaties (bijv. CCSD(T)) noodzakelijk kan zijn voor een hogere nauwkeurigheid in deze specifieke systemen.
5. Toepassing op TiO$_2$: De methode is toegepast op rutiel TiO$_2$, waarbij een TDL bandkloof van 4,17 eV werd voorspeld. Dit overschat de ZPR-gecorrigeerde experimentele waarde met ~0,3 eV, wat consistent is met de algemene trend in andere vaste stoffen, hoewel de auteurs opmerken dat de aanwezigheid van overgangsmetalen en anisotrope meshes extra onzekerheid introduceert.

Significantie
Het artikel beweert de eerste betrouwbare, geconvergeerde CCSD-benchmarks te leveren voor de cohesie-energieën en bandkloven van eenvoudige vaste stoffen, waarbij zowel de eindige-omvang- als de basisset-onvolledigheidsfouten effectief zijn geëlimineerd. Door aan te tonen dat standaard extrapolatieschema's (zoals Model A) bandkloven systematisch onderschatten en dat het GW-EOM schema niet universeel robuust is, vestigt dit werk een nieuwe standaard voor het evalueren van de nauwkeurigheid van many-body theorieën in vaste stoffen. De resultaten suggereren dat hoewel CCSD zeer accuraat is voor cohesie-energieën, het een systematische overschatting van bandkloven behoudt (~0,4 eV), wat waarschijnlijk voortkomt uit een onvoldoende behandeling van elektroncorrelatie (specifiek triple excitaties) of elektron-fonon interacties in ionische systemen. De ontwikkelde software en data dienen als fundament voor toekomstige hoogwaardige berekeningen in de vaste stof chemie.