Quantum chemistry based on classical mechanics inspired by simulated bifurcation

Deze paper introduceert de Simulated Bifurcation-based CI (SBCI)-algoritme, een op klassieke mechanica gebaseerde methode die de rekenkosten voor nauwkeurige quantumchemische berekeningen verlaagt zonder in te boeten aan precisie.

De kern

De Digitale Gids voor Moleculaire Puzzels: Een Nieuwe Weg naar Snellere Berekeningen

Stel je voor dat je een gigantische, ingewikkelde puzzel moet oplossen. De stukjes zijn atomen en elektronen, en de puzzel is een molecule (zoals water of stikstof). Om te begrijpen hoe deze moleculen werken, moeten we de "energie" van deze puzzel berekenen. Hoe lager de energie, hoe stabieler de molecule.

In de chemie noemen we dit een Configuration Interaction (CI) berekening. Het probleem? Voor grote moleculen zijn er zoveel mogelijke manieren om de elektronen te rangschikken, dat het aantal combinaties explosief groeit. Het is alsof je probeert elke mogelijke combinatie van een slot met 100 cijfers te proberen.

Het Oude Probleem: De "Davidson" Methode

Voor decennia hebben wetenschappers een standaardmethode gebruikt, de Davidson-methode.

• De Analogie: Stel je voor dat je in een donker, enorm kasteel (het molecule) op zoek bent naar de laagste verdieping (de laagste energie). De Davidson-methode is als een team van verkenners dat systematisch elke kamer afloopt, maar ze moeten constant terugkeren naar de ingang om te checken welke kamers ze al hebben bezocht. Ze dragen een zware rugzak met alle eerdere routes bij zich.
• Het Nadeel: Naarmate het kasteel groter wordt, wordt de rugzak zo zwaar dat het team langzaam wordt en veel ruimte in beslag neemt. Het kost enorm veel tijd en computergeheugen.

De Nieuwe Oplossing: SBCI (De "Simulated Bifurcation" Methode)

De auteurs van dit paper, Fumihiko Aiga en Hayato Goto, hebben een slimme nieuwe methode bedacht die SBCI heet. Ze zijn geïnspireerd door een techniek die oorspronkelijk werd ontwikkeld voor het oplossen van complexe logistieke problemen (zoals het optimaliseren van verkeerslichten of bezorgroutes).

• De Analogie: In plaats van een team dat systematisch kamers afloopt, sturen ze bij SBCI een energieke, springerige bal het kasteel in.
  • Deze bal rolt niet zomaar; hij volgt de wetten van de klassieke mechanica (zoals een balletje dat over heuvels en in dalen rolt).
  • De "heuvels" in dit landschap zijn de hoge energiewaarden (wat we niet willen), en de "dalen" zijn de lage energiewaarden (wat we wel willen).
  • De bal heeft een eigen momentum (snelheid). Als hij in een dal rolt, schiet hij er soms overheen, maar door slimme regels (die ze "bifurcatie" noemen) wordt hij steeds vaker naar het diepste dal getrokken.

Wat maakt SBCI zo speciaal?

1. Geen zware rugzak:
   De Davidson-methode moet onthouden welke routes ze al hebben genomen. SBCI heeft dat niet nodig. De "bal" (de berekening) onthoudt alleen waar hij nu is en hoe snel hij gaat.

   • Vergelijking: Davidson moet een kaart van het hele kasteel bijhouden. SBCI kijkt gewoon waar hij nu staat en waar de grond naar beneden loopt. Dit bespaart enorm veel computergeheugen.

2. Sneller vinden van het diepste punt:
   Omdat de bal momentum heeft, kan hij sneller over kleine hindernissen springen en direct naar het diepste dal stuiven.

   • Vergelijking: Waar de Davidson-verkenners elke kamer één voor één aflopen, springt de SBCI-bal als een acrobaat direct naar de beste plek. Dit bespaart tijd.

3. Twee versies voor verschillende taken:

   • SBCI1: Zoekt één oplossing tegelijk (bijvoorbeeld de grondtoestand van een molecule). Dit is als het vinden van het laagste punt in een berglandschap.
   • SBCI2: Zoekt twee oplossingen tegelijk (bijvoorbeeld de grondtoestand én de eerste aangeslagen toestand). Dit is als twee ballen die parallel rollen, maar elkaar niet blokkeren. Dit werkt heel goed als de "dalen" dicht bij elkaar liggen (nearly degenerate states), wat voor de oude methode vaak een probleem was.

De Resultaten: Sneller, Lichter, Net zo Accuraat

De auteurs hebben deze nieuwe methode getest op verschillende moleculen (zoals stikstof, water en koolstof).

• Snelheid: De berekeningen waren overal sneller dan met de oude Davidson-methode. Soms wel 2 tot 3 keer sneller.
• Geheugen: De computer had minder geheugen nodig om dezelfde berekening te doen.
• Nauwkeurigheid: Het belangrijkste: de uitkomsten waren exact hetzelfde als de oude, betrouwbare methode. Ze hebben geen kwaliteit ingeleverd voor snelheid.

Waarom is dit belangrijk?

Vandaag de dag zijn supercomputers nodig om complexe moleculen te simuleren, wat essentieel is voor het ontwikkelen van nieuwe medicijnen, batterijen of materialen.

• Met SBCI kunnen wetenschappers grotere en complexere moleculen bestuderen zonder dat ze duizelingwekkende rekenkracht nodig hebben.
• Het is alsof je van een oude, zware fiets (Davidson) overstapt op een snelle, lichte racefiets (SBCI) die precies dezelfde route aflegt, maar veel sneller en met minder inspanning.

Kortom: Dit paper introduceert een slimme, klassieke mechanische "dans" die computers helpt om moleculaire puzzels veel efficiënter op te lossen, zonder de nauwkeurigheid te verliezen. Het is een grote stap voorwaarts in de digitale chemie.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "Quantum chemistry based on classical mechanics inspired by simulated bifurcation" in het Nederlands.

Titel: Kwantumchemie gebaseerd op klassieke mechanica, geïnspireerd door gesimuleerde bifurcatie

Auteurs: Fumihiko Aiga en Hayato Goto (Toshiba Corporation en RIKEN Center for Quantum Computing)

---

1. Het Probleem

Binnen de kwantumchemie is het berekenen van de elektronische structuur van moleculen cruciaal voor het begrijpen van moleculaire eigenschappen. De Configuratie-Interactie (CI) methode, en specifiek de Full Configuration Interaction (FCI), biedt exacte oplossingen binnen een gegeven basisset. Echter, het aantal configuraties neemt combinatorisch toe met het aantal elektronen en orbitalen, wat leidt tot een enorme rekenlast.

Voor grote systemen is FCI op klassieke computers vaak onhaalbaar. Hoewel kwantumcomputers veelbelovend zijn, zijn ze nog niet operationeel op de vereiste schaal. Op klassieke computers wordt momenteel de Davidson-methode veel gebruikt voor het vinden van eigenwaarden en eigenvectoren in CI-berekeningen. Deze methode heeft echter beperkingen:

• Het vereist het opslaan van de geschiedenis van residu-vector voor alle niet-geconvergeerde toestanden, wat leidt tot hoge geheugengebruik.
• Het kost veel tijd voor inwendige producten en matrix-productconstructies, vooral bij het berekenen van meerdere aangeslagen toestanden of bijna-ontartde toestanden.

Er is een dringende behoefte aan algoritmen die op klassieke computers draaien, maar met lagere rekenkosten (tijd en geheugen) en vergelijkbare nauwkeurigheid als de Davidson-methode.

2. Methodologie: SBCI (Simulated Bifurcation-based CI)

De auteurs introduceren een nieuw algoritme genaamd SBCI, dat is afgeleid van een kwantum-geïnspireerd algoritme voor combinatorische optimalisatie, bekend als Simulated Bifurcation (SB). In plaats van het eigenwaardeprobleem om te zetten in een Ising-probleem (zoals bij Ising-machines), mapt SBCI het optimalisatieprobleem van de CI-coëfficiënten direct af op klassieke mechanische dynamica.

Kernprincipes van SBCI:

• Klassieke Hamiltoniaan: De CI-coëfficiënten worden behandeld als posities ($\mathbf{x}$) van een klassiek systeem, en er worden geconjugeerde impulsen ($\mathbf{y}$) geïntroduceerd. Een klassieke Hamiltoniaan wordt geconstrueerd waarbij de potentiële energie evenredig is met de Rayleigh-kwotiënt (de energie die geminimaliseerd moet worden).
• Symplectische Euler-methode: De bewegingsvergelijkingen van Hamilton worden opgelost met behulp van de symplectische Euler-methode (een numeriek integratieschema).
• Variatiele parameterbepaling: In tegenstelling tot de originele SB-algoritmen die adiabatische evolutie gebruiken, bepaalt SBCI de parameters ($b_t, c_t$) variatieel. Dit betekent dat ze zo worden gekozen dat de geüpdatete vector evenwijdig is aan de laagste (of opvolgende) eigenvector in de onderruimte die wordt opgespannen door de huidige vector, de impuls en het residu.
• Adaptieve Restart: Een uniek kenmerk is dat de trial golffunctie niet genormaliseerd is tijdens de iteratie. De grootte van de vector en de parameters worden gebruikt om adaptieve restarts te triggeren, wat de convergentie versnelt.

Twee Varianten:

1. SBCI1: Update één staat tegelijkertijd. Toestanden worden sequentieel van laag naar hoog bepaald. Dit verlaagt de rekenkosten door minder vectoren op te slaan.
2. SBCI2: Update twee toestanden ($\alpha$ en $\alpha+1$) gelijktijdig. Dit is vooral effectief voor bijna-ontartde toestanden, waar sequentiële updates (SBCI1) trager kunnen convergeren.

3. Belangrijkste Bijdragen

• Nieuw Algoritme: De ontwikkeling van SBCI1 en SBCI2 als efficiënte alternatieven voor de Davidson-methode voor het oplossen van CI-eigenwaardeproblemen.
• Implementatie: De algoritmen zijn geïmplementeerd in PySCF (Python-based Simulations of Chemistry Framework), een veelgebruikt open-source pakket.
• Theoretische Innovatie: Het toepassen van principes uit de gesimuleerde bifurcatie (oorspronkelijk voor combinatorische optimalisatie) op het probleem van elektronische structuurberekening via klassieke mechanica, zonder gebruik te maken van kwantumhardware.
• Efficiëntie: Het verminderen van het geheugengebruik door slechts een beperkt aantal vectoren (2 voor SBCI1, 4 voor SBCI2) op te slaan tijdens het update-proces, in tegenstelling tot de Davidson-methode die de geschiedenis van alle residu-vectoren bewaart.

4. Resultaten

De auteurs hebben SBCI getest op diverse moleculaire systemen, waaronder $N_2$, $CN$, $H_2O$, $HF$, $BH$, en $C_2$, voor zowel grondtoestanden als aangeslagen toestanden.

• Nauwkeurigheid: De berekende energieën van SBCI1 en SBCI2 vertonen dezelfde nauwkeurigheid als de Davidson-methode in PySCF en zijn consistent met waarden uit de literatuur. De verschillen in energie liggen vaak in de orde van $10^{-9}$tot$10^{-10}$ Hartree.
• Rekentijd:
  • Voor grondtoestanden ($N_2$ en $CN$) was de uitvoeringstijd van SBCI1 korter dan die van de Davidson-methode bij alle geteste bindingslengten. Het voordeel nam toe bij langere bindingslengten (waar de elektronische correlatie sterker is).
  • Voor aangeslagen toestanden was SBCI2 over het algemeen de snelste methode, gevolgd door SBCI1, en vervolgens de Davidson-methode.
• Geheugengebruik: Zowel SBCI1 als SBCI2 verbruiken aanzienlijk minder geheugen dan de Davidson-methode, wat essentieel is voor grotere systemen.
• Convergentie: De adaptieve restart-mechanismen bleken effectief om stagnatie te voorkomen en snellere convergentie te bereiken.

5. Betekenis en Toekomstperspectief

Deze studie toont aan dat het mogelijk is om hoge-precisie elektronische structuurberekeningen (FCI) uit te voeren op klassieke computers met aanzienlijk lagere kosten dan de huidige industriestandaard (Davidson), zonder in te leveren op betrouwbaarheid.

• Praktische Toepassing: SBCI kan direct worden gebruikt in bestaande workflows (zoals PySCF) om grotere moleculaire systemen te bestuderen of om berekeningen sneller uit te voeren.
• Schalbaarheid: Omdat SBCI minder vectoren opslaat, wordt het geacht zeer goed te schalen naar multi-node omgevingen (via MPI), wat de weg vrijmaakt voor nog grobere berekeningen.
• Breder Toepassingsgebied: Hoewel deze studie zich richtte op CI, kan het algoritme ook worden toegepast op het optimaliseren van CI-coëfficiënten in multiconfiguratie self-consistent field (MCSCF) berekeningen, zoals CASSCF.

Kortom, SBCI belooft een nieuwe standaard te worden voor CI-berekeningen, waarbij klassieke mechanica wordt gebruikt om de complexiteit van kwantumchemische problemen op een efficiëntere manier aan te pakken.