Development of ab initio Hubbard parameter calculation schemes in the k-point sampling real-time TDDFT program in CP2K

De auteurs hebben in het CP2K-programma nieuwe *ab initio*-methoden voor het berekenen van Hubbard-parameters geïmplementeerd, waaronder een uitgebreide lineaire-reactiebenadering voor energie-afhankelijke parameters, en vergelijken deze met het ACBN0-schema om hun respectievelijke toepassingsgebieden in statische en dynamische berekeningen te verduidelijken.

De kern

Stel je voor dat je een enorme, ingewikkelde stad probeert te begrijpen. In deze stad wonen miljarden kleine deeltjes: elektronen. De wetenschappers die deze stad bestuderen, gebruiken een soort "stadsplanner" genaamd DFT (Density Functional Theory). Deze planner is geweldig voor de meeste gebouwen en straten, maar hij heeft een groot probleem: hij faalt volledig als hij naar de wijken met sterke spanningen kijkt.

In deze specifieke wijken (zoals in materialen met overgangsmetalen, bijvoorbeeld roest of koper) duwen de elektronen elkaar zo hard weg dat de standaardplanner het niet meer snapt. Hij maakt fouten, alsof hij denkt dat twee buren die elkaar haten, toch hand in hand wandelen.

Om dit op te lossen, hebben wetenschappers een extra tooltje uitgevonden: de Hubbard-parameter. Je kunt dit zien als een speciale "afstandshouder". Als je deze afstandshouder in de planner zet, zegt hij: "Oké, deze elektronen mogen niet te dicht bij elkaar komen, ze moeten op een veilige afstand blijven."

Het probleem is echter: Hoe groot moet die afstandshouder precies zijn?
Vroeger moesten onderzoekers dit raden of het afstemmen op experimenten, wat vaak leidde tot een "one-size-fits-all" oplossing die niet voor elke situatie werkte.

In dit paper hebben Kota Hanasaki en Sandra Luber twee nieuwe, slimme manieren bedacht om die afstandshouder automatisch en precies te berekenen, direct vanuit de natuurwetten (zodat je niet hoeft te gissen). Ze hebben dit ingebouwd in een krachtige software genaamd CP2K.

Hier is hoe hun twee nieuwe methoden werken, vertaald in alledaagse termen:

1. De "Real-time Camera" Methode (ACBN0)

Stel je voor dat je een film maakt van de elektronen die door de stad rennen.

• Hoe het werkt: Deze methode kijkt naar de "foto's" van de elektronen op elk moment. Als de elektronen beginnen te dansen (bijvoorbeeld door een laserflits), past deze methode de afstandshouder direct aan op basis van wat ze op dat moment zien.
• De kracht: Het is supersnel en werkt perfect voor films (dynamische situaties). Het is alsof je een slimme camera hebt die de afstand tussen mensen automatisch aanpast terwijl ze rennen.
• Het nadeel: De formule die ze gebruiken is een beetje een "zwarte doos". Het werkt in de praktijk, maar de theoretische wiskunde erachter is niet volledig verklaard. Het is alsof je een auto hebt die perfect rijdt, maar niemand weet precies hoe de motor onder de motorkap werkt.

2. De "Echo-test" Methode (Lineair Respons)

Stel je voor dat je in een grote grot staat en je roept "Hallo!". De echo die je terugkrijgt, vertelt je iets over de grootte en vorm van de grot.

• Hoe het werkt: Deze methode geeft een heel klein, zacht duwtje aan de elektronen (een "echo-test") en luistert heel nauwkeurig naar hoe ze reageren. Uit die reactie kunnen ze precies afleiden hoe sterk de elektronen elkaar afstoten.
• De kracht: Dit is gebaseerd op stevige, klassieke natuurkunde. Het is alsof je een wetenschappelijke meting doet in plaats van te gokken. Bovendien hebben ze deze methode nu zo ver ontwikkeld dat ze niet alleen de "standaard-echo" kunnen meten, maar ook hoe de echo verandert als je het geluid sneller of trager maakt (energie-afhankelijkheid).
• Het nadeel: Het is rekenkundig zwaar. Het is alsof je in plaats van één keer te roepen, duizenden keren moet roepen in verschillende richtingen om het perfecte plaatje te krijgen. Het kost dus veel meer tijd en rekenkracht.

Wat hebben ze ontdekt?

De auteurs hebben getest of deze methoden werken op verschillende materialen (zoals roest en koper-oxide).

• Resultaat: Beide methoden geven goede resultaten voor statische situaties (een stilstaande foto van de stad). Ze voorspellen de eigenschappen van het materiaal vrijwel even goed.
• Het verschil:
  • Wil je kijken hoe een materiaal reageert op een flitsend laserlicht (zoals in een futuristische computer of een snelle schakelaar)? Gebruik dan de "Real-time Camera" (ACBN0). Die is snel en flexibel.
  • Wil je de fundamentele, exacte natuurwetten begrijpen en weten hoe de elektronen reageren op verschillende energie-niveaus? Gebruik dan de "Echo-test". Die is theoretisch sterker, maar kost meer tijd.

Waarom is dit belangrijk?

Vroeger moesten wetenschappers de "afstandshouder" (de Hubbard-parameter) handmatig instellen, wat vaak leidde tot fouten in het ontwerp van nieuwe materialen. Met deze nieuwe tools kunnen ze nu:

1. Nieuwe materialen ontwerpen voor betere batterijen, snellere computers en efficiëntere zonnepanelen.
2. Simulaties draaien die eerder onmogelijk waren, zoals het simuleren van hoe materialen reageren op extreem snelle lichtflitsen (attoseconden).

Kortom: Ze hebben twee nieuwe, slimme gereedschappen in de gereedschapskist van de materiaalkundige gelegd. De ene is een snelle, slimme camera voor beweging, en de andere is een nauwkeurige meetlat voor de diepere waarheid. Samen maken ze het mogelijk om de complexe wereld van elektronen beter te begrijpen dan ooit tevoren.

Technische samenvatting

Titel: Ontwikkeling van ab initio berekeningsschema's voor Hubbard-parameters in het k-punt-sampling real-time TDDFT-programma in CP2K

Auteurs: Kota Hanasaki en Sandra Luber (Universiteit van Zürich)

---

1. Het Probleem

De dichtheidsfunctionaaltheorie (DFT) is een krachtig hulpmiddel voor de berekening van materiaaleigenschappen, maar faalt vaak bij materialen met sterke elektron-elektroncorrelaties, zoals overgangsmetaaloxiden (met gelokaliseerde d- of f-orbitalen). Dit komt door zelfinteractie-fouten in lokale en semi-lokale uitwisselings-correlatie (xc) functionalen.

• DFT+U: Een veelgebruikte oplossing is DFT+U, dat een expliciete on-site interactie introduceert via de Hubbard-parameter ($U$) en de Hund-koppeling ($J$).
• De Uitdaging: De waarden van $U$ en $J$ zijn cruciaal. Empirische aanpassing is onbetrouwbaar omdat de optimale waarden afhangen van de chemische omgeving en de te modelleren fysische grootheid.
• Beperkingen van bestaande methoden:
  • cRPA (Constrained Random Phase Approximation): Kan energievafhankelijke parameters berekenen, maar neigt tot overschatting van afscherming in materialen met sterke hybridisatie.
  • ACBN0: Een "mean-field"-benadering die goed werkt voor dynamische simulaties, maar whose formules niet strikt uit de standaard veel-deeltjestheorie zijn afgeleid, waardoor het gedrag in dynamische situaties moeilijk te analyseren is.
  • Lineaire Respons: Bestaande lineaire respons-methoden zijn voornamelijk beperkt tot statische (energie-onafhankelijke) parameters.

Het doel van dit paper is het implementeren en uitbreiden van ab initio methoden voor het berekenen van Hubbard-parameters binnen het CP2K-framework, specifiek gericht op real-time TDDFT (RT-TDDFT) met k-punt-sampling, inclusief een nieuwe methode voor energievafhankelijke parameters.

---

2. Methodologie

De auteurs hebben twee verschillende schema's geïmplementeerd en uitgebreid in CP2K:

A. ACBN0 (Quasi-hybrid functional)

• Principe: Berekent $U$ en $J$ puur op basis van de statische correlaties in de renormaliseerde dichtheidsmatrix en de Kohn-Sham orbitalen.
• Implementatie: De auteurs hebben de projectieschema's voor de projectie-operator $P^{(I,A)}$ aangepast voor k-punt-sampling in periodieke systemen.
• Voordeel: Omdat de berekening puur op de momentane dichtheidsmatrix gebaseerd is, is deze methode direct toepasbaar op real-time simulaties (adiabatisch) om tijdsafhankelijke Hubbard-parameters te verkrijgen.

B. Minimum-Tracking Lineaire Respons (MT-LR)

• Principe: Gebaseerd op de respons van het systeem op een externe verstoring in de geprojecteerde subspace. De Hubbard-parameters worden bepaald door de afgeleide van het Hartree-xc-potentieel ten opzichte van de bezettingsgetallen.
• Implementatie: Uitgebreid van de eerdere $\Gamma$-punt-implementatie naar k-punt-sampling.
  • Een kleine verstoring ($\lambda$) wordt toegepast op specifieke orbitalen.
  • De respons van de bezetting en het potentieel wordt gemeten.
  • $U$ en $J$ worden berekend via numerieke differentiatie (verschil tussen verstoorde en onverstoorde toestanden).

C. Nieuwe Uitbreiding: Energie-afhankelijke Hubbard-parameters

• Theorie: De auteurs hebben de lineaire respons-theorie uitgebreid naar de tijdsdomein (en via Fourier-transformatie naar het frequentiedomein).
• Formulering: Door een impulsverstorende operator toe te passen in RT-TDDFT en de respons te volgen, kunnen ze de responskernen $U(\omega)$ en $J(\omega)$ berekenen.
• Voordeel: Deze methode houdt rekening met de xc-functional van DFT (in tegenstelling tot cRPA) en kan dynamische afscherming en energievafhankelijkheid modelleren.

---

3. Belangrijkste Resultaten

Statische Berekeningen

De auteurs hebben de methoden getest op zes overgangsmetaalverbindingen (ScN, TiO2, MnO, NiO, ZnO, CuO).

• Nauwkeurigheid: Zowel ACBN0 als de MT-LR methode leveren bandgaps en magnetische momenten op die redelijk overeenkomen met experimentele waarden.
• Verschil in Parameters: Er is geen duidelijke overeenstemming tussen de berekende Hubbard-parameters ($U_{eff}$) van de twee methoden. De waarden verschillen aanzienlijk (bijv. voor NiO: ACBN0 geeft ~3.6 eV, MT-LR geeft ~5.4 eV).
• Conclusie: Geen van beide methoden overtreft de andere duidelijk in nauwkeurigheid voor statische eigenschappen; de keuze hangt af van het specifieke systeem.
• Kosten: De MT-LR methode is computatieel veel duurder vanwege de noodzaak van supercellen en iteratieve zelfconsistentie.

Dynamische Berekeningen (RT-TDDFT)

• ACBN0 in sterke velden: Simulaties van NiO onder een sterke laserpuls tonen een daling van $U_{eff}$ tijdens de excitatie. Dit komt overeen met eerdere studies (Cazali et al.) en wordt toegeschreven aan elektronische excitatie en versterkte afscherming. De amplitude van de daling (~140 meV) is redelijk in lijn met experimenten.
• Energie-afhankelijkheid (MT-LR): De nieuwe methode voor $U(\omega)$ toont dat de parameters bij $\omega=0$ overeenkomen met de statische waarden. Bij hoge frequenties ($\omega \to \infty$) convergeren sommige parameters correct naar de onafgeschermde Coulomb-interactie, terwijl andere (zoals bij Mn) numerieke problemen vertonen die nog verder onderzocht moeten worden.

---

4. Bijdragen en Innovaties

1. Implementatie in CP2K: Eerste implementatie van zowel ACBN0 als de MT-LR methode in een k-punt-sampling RT-TDDFT-framework.
2. Nieuwe Theorie voor $U(\omega)$: Een nieuwe, lineaire respons-gebaseerde methode om energievafhankelijke Hubbard-parameters te berekenen die de xc-functional van de gebruikte DFT-benadering meeneemt. Dit is een alternatief voor cRPA.
3. Vergelijking van Methodes: Een grondige vergelijking tussen ACBN0 (mean-field, dynamisch toepasbaar maar theoretisch minder strikt) en MT-LR (theoretisch robuust, maar rekenintensief).
4. Validatie: Uitgebreide validatie op een reeks TMO/N-materialen, inclusief analyse van bandstructuren en magnetische momenten.

---

5. Significatie en Toekomstperspectief

• Dynamische Materialen: De methode maakt het mogelijk om niet-evenwichtsdynamica (zoals laser-gedreven processen) in sterk gecorreleerde materialen nauwkeuriger te modelleren door tijds- en energievafhankelijke Hubbard-parameters te gebruiken.
• Theoretische Robuustheid: De nieuwe lineaire respons-aanpak voor $U(\omega)$ biedt een theoretisch onderbouwde manier om afschermingseffecten te berekenen die afhankelijk zijn van de gekozen xc-functional, wat cruciaal is voor DFT+U toepassingen.
• Toekomstige Toepassingen: De resultaten kunnen leiden tot verbeterde DFT+U($\omega$) schema's (zoals voorgesteld door Vanzini en Marzari) voor het modelleren van complexe materialen met sterke hybridisatie.
• Beperkingen: De hoge rekenkosten van de lineaire respons-methode en de numerieke instabiliteit bij hoge frequenties zijn huidige beperkingen die verdere optimalisatie vereisen.

Kortom, dit paper levert een essentiële stap voorwaarts in de berekening van dynamische elektronische eigenschappen van gecorreleerde materialen door geavanceerde ab initio Hubbard-parameters te integreren in moderne RT-TDDFT-simulaties.