QCommute: a tool for symbolic computation of nested commutators in quantum many-body spin-1/2 systems

Probleemstelling

Het kwantitatief beschrijven van de dynamica van kwantum-veeldeeltjessystemen in niet-perturbatieve regimes is een uiterst moeilijke taak. Hoewel er in één dimensie aanzienlijke vooruitgang is geboekt met technieken zoals Matrix Product States (MPS), blijft het probleem in twee en vooral drie dimensies grotendeels onopgelost. Bestaande methoden (zoals variatiemethoden, Monte-Carlo en neurale netwerken) werken vaak in het Schrödinger-beeld en hebben te kampen met de exponentiële groei van de Hilbertruimte.

Methoden in het Heisenberg-beeld (zoals de recursiemethode en sparse Pauli-dynamica) zijn veelbelovend, maar hun efficiëntie hangt volledig af van de snelheid en nauwkeurigheid waarmee geneste commutatoren $[H, [H, \dots, [H, A]\dots]]$ kunnen worden berekend. Een groot nadeel van bestaande softwaretools is dat ze vaak:

1. Op eindige roosters werken (wat randeffecten introduceert).
2. Numerieke waarden voor Hamiltonian-parameters vereisen (wat betekent dat elke nieuwe parametercombinatie een nieuwe berekening vereist).
3. Beperkt zijn door het geheugengebruik bij hoge geneste orde.

Methodologie

De auteurs presenteren QCOMMUTE, een C++-tool die is ontworpen om geneste commutatoren symbolisch te berekenen voor spin-1/2 systemen op hyperkubische roosters (1D, 2D en 3D).

Kernconcepten:

• Heisenberg-beeld: De tijdsontwikkeling van een observabele $A$ wordt benaderd via een Taylor-reeks van geneste commutatoren met de Liouvilliaan $L = [H, \cdot]$.
• Pauli-strings: Operatoren worden uitgebreid in een basis van Pauli-strings (producten van Pauli-matrices op verschillende roosterplekken).
• Translatie-invariantie: Het systeem veronderstelt oneindige roosters zonder randen. Operatoren worden weergegeven in een "canonieke vorm" waarbij een lokale dichtheid wordt gedistribueerd over het hele roosters via een translatie-operator $T$. Dit zorgt voor een enorme compressie van de data.
• Symbolische Berekening: In plaats van numerieke waarden, worden de koppelingsconstanten van de Hamiltoniaan behandeld als symbolische parameters (polynomen). Dit betekent dat één enkele berekening de volledige parameter ruimte dekt.

Algoritme en Implementatie:

1. Iteratieve Berekening: Het algoritme berekent iteratief $A^{(n+1)} = [H, A^{(n)}]$.
2. Datastructuur:
   • Ruimtelijke coördinaten worden opgeslagen als template-klassen.
   • Coëfficiënten worden opgeslagen als polynomen met gehele getallen.
   • GMP-bibliotheek: Om overloop te voorkomen bij de exponentiële groei van coëfficiënten, wordt de GMP-bibliotheek gebruikt voor willekeurig nauwkeurige geheeltallige rekenkunde.
3. Geheugenbeheer (Disk Swapping): De grootste uitdaging is de exponentiële groei van het aantal termen. QCOMMUTE slaat niet alles in het RAM-geheugen op. In plaats daarvan wordt de data dynamisch opgesplitst in "chunks" en tussentijds op de schijf (in binaire .comm bestanden) geschreven. Dit maakt het mogelijk om diep geneste commutatoren te berekenen op standaard hardware.
4. Parallelisatie: De code is sterk geoptimaliseerd voor parallelle verwerking om de prestaties te maximaliseren.

Belangrijkste Bijdragen

1. Werking in de Thermodynamische Limiet: QCOMMUTE werkt direct op oneindige roosters, waardoor randeffecten volledig worden geëlimineerd. Dit is uniek vergeleken met tools die werken op eindige roosters.
2. Symbolische Parametrisatie: De tool behoudt Hamiltonian-parameters als symbolische variabelen. Dit stelt onderzoekers in staat om de dynamica voor een continuüm van parameterwaarden te analyseren zonder herhaalde berekeningen.
3. Schaalbaarheid: Door het gebruik van schijfcaching en willekeurig nauwkeurige rekenkunde, kan de tool veel dieper geneste commutatoren berekenen dan eerder mogelijk was, zelfs in 2D en 3D.
4. Open Source: De code is beschikbaar via GitLab en is geschreven in C++17, afhankelijk van de GMP-bibliotheek.

Resultaten

De auteurs testen QCOMMUTE op het Ising-model in 1D, 2D en 3D:

• Prestaties: De tool bereikt recordniveaus voor de maximale geneste orde ($n_{max}$).
  • In 1D (gemengd veld Ising-model) wordt $n_{max} = 48$ bereikt.
  • In 2D wordt $n_{max} = 23$ bereikt (waarbij eerdere methoden vaak faalden bij $n \approx 20$ door geheugenoverschrijding).
  • In 3D wordt $n_{max} = 17$ bereikt.
• Validatie: De resultaten zijn geverifieerd tegen exacte diagonalisatie (in 1D), andere softwaretools (zoals FORM-gebaseerde codes en Julia-packages), en de recursiemethode.
• Toepassing op Dynamica:
  • Quench-dynamica: De tool berekent de tijdsontwikkeling van magnetisatie na een kwantum-quench. De resultaten komen overeen met de meest nauwkeurige bestaande methoden (iPEPS, SPD).
  • Autocorrelatiefuncties: De berekende momenten worden gebruikt om Taylor-reeksen voor autocorrelatiefuncties op te stellen. Omdat dit alternerende reeksen zijn, bieden de afgeknipte reeksen strikte bovengrenzen en ondergrenzen voor de werkelijke dynamica. In veel gevallen zijn deze grenzen zo strak dat ze de exacte oplossing benaderen.

Betekenis en Toekomstperspectief

QCOMMUTE vertegenwoordigt een doorbraak in de numerieke studie van kwantum-dynamica in hoge dimensies. Door de combinatie van werken in de thermodynamische limiet, symbolische parametrisatie en efficiënt geheugenbeheer, opent het nieuwe mogelijkheden voor het bestuderen van sterk gecorreleerde systemen die voorheen ontoegankelijk waren.

Toekomstige ontwikkelingen die door de auteurs worden genoemd, zijn:

• Uitbreiding naar niet-hyperkubische roosters.
• Ondersteuning voor willekeurige beginstaten (nu beperkt tot productstaten).
• Directe implementatie van het Pauli-propagatie-algoritme binnen de package.

Samenvattend biedt QCOMMUTE een robuust, schaalbaar en nauwkeurig instrument voor theoretische fysici die zich bezighouden met de dynamica van kwantum-veeldeeltjessystemen, en het overbrugt de kloof tussen perturbatieve benaderingen en volledig numerieke methoden die lijden onder randeffecten of parameterafhankelijkheid.