Efficient computation of real-time correlators using Pauli Propagation

Technische Samenvatting: Efficiënte Berekening van Real-Time Correlatoren middels Pauli-propagatie

Probleemstelling
Het berekenen van real-time twee-punts tijdgeordende correlatiefuncties, $C(t) = \text{Tr}[\rho A(t)B]$, vormt een fundamentele uitdaging in de veel-deeltjesfysica, en dient als de brug tussen microscopische Hamiltoniaanse details en macroscopische experimentele grootheden (bijv. neutronenverstrooiing, conductiviteit). Traditionele methoden kampen met aanzienlijke knelpunten: exacte diagonalisatie wordt beperkt door geheugen, real-time Quantum Monte Carlo lijdt onder het tekenprobleem, en tensornetwerk-methoden (zoals TEBD of TDVP) hebben moeite met het schalen naar meer dan één dimensie vanwege de snelle groei van de bond-dimensie tijdens de tijdsevolutie. Hoewel Pauli-propagatie (PP) een veelbelovend alternatief biedt door operatoren direct in de Heisenberg-representatie te evolueren, is het inherent beperkt door de exponentiële proliferatie van Pauli-strings over de tijd, wat afkaping (truncation) noodzakelijk maakt, wat weer bias introduceert en het toegankelijke tijdsvenster beperkt.

Methodologie
De auteurs stellen een hybride computationele strategie voor die korte-tijd Pauli-propagatie combineert met geavanceerde signaalverwerking en grondtoestands-benaderingstechnieken om correlatiefuncties buiten het directe simulatievenster uit te breiden. De methodologie bestaat uit drie primaire componenten:

1. Afgekapte Pauli-propagatie: De tijdsgeëvolueerde operator $A(t)$ wordt geëxpandeerd in een basis van Pauli-strings, $A(t) = \sum_i a_i(t) P_i$. Om de exponentiële groei van termen te beheersen, gebruiken de auteurs twee afkapstrategieën:

   • Coëfficiënt-afkap ($\epsilon$): Het verwijderen van strings waar de coëfficiënt $|a_i| < \epsilon$ is.
   • Gewichts-afkap ($w$): Het verwijderen van strings waar het Pauli-gewicht (aantal niet-identiteit sites) groter is dan $w$.
     De evolutie wordt uitgevoerd via Trotterisatie, waarbij elke stap Pauli-strings vertakt op basis van commutatierelaties met de termen van de Hamiltoniaan.

2. Clifford-geaugmenteerde Matrix Product States (CAMPS): Om de moeilijkheid aan te pakken van het evalueren van verwachtingswaarden $\langle \Psi | A(t) | \Psi \rangle$ voor verstrengelde grondtoestanden (met name in 2D), maken de auteurs gebruik van CAMPS. Deze aanpak past Clifford-transformaties toe op de Matrix Product State (MPS) om de benodigde bond-dimensie aanzienlijk te verkleinen terwijl de coëfficiëntstructuur van de Pauli-som behouden blijft. Dit maakt een efficiënte evaluatie van verwachtingswaarden mogelijk tegen complexe grondtoestanden die anders computationeel onhaalbaar zouden zijn voor standaard MPS-methoden.

3. Positief Semi-Definiete (PSD) Extensie: Erkennend dat correlatiefuncties van de vorm $C_{AA}(t) = \text{Tr}[\rho A^\dagger(t)A(0)]$ positieve functies zijn, maken de auteurs gebruik van het feit dat hun tijd-gediscretiseerde Toeplitz-matrices Positief Semi-Definiet (PSD) zijn. Gebaseerd op het werk van Carathéodory, Toeplitz en Pisarenko, wordt een rank-$r$ PSD Toeplitz-matrix uniek bepaald door $2r$ ruisvrije datapunten. De auteurs gebruiken een low-rank PSD-fit om de korte-tijd PP-data te "super-resolven", waardoor de correlatiefunctie naar latere tijden wordt uitgebreid zonder verdere operator-evolutie. Dit mitigeert effectief de "Pauli-explosie" door te vertrouwen op de observatie dat veel fysieke systemen worden gedomineerd door een klein aantal karakteristieke frequenties.

Belangrijkste Resultaten
De auteurs demonstreren de effectiviteit van deze gecombineerde aanpak op 1D en 2D Heisenberg-modellen:

• 1D Ferromagnetische Heisenberg-keten ($N=64$): De auteurs tonen aan dat zelfs met agressieve afkaping (waarbij de tijdssporen alleen kwalitatief overeenkomen met exacte resultaten), de Fourier-getransformeerde spectra de correcte fysieke kenmerken behouden (bijv. de magnon-dispersie). Het gebruik van een ad-hoc dempingsfunctie of de PSD-extensie herstelt de spectrale kenmerken ondanks het verlies van lange-tijd dynamische nauwkeurigheid in de ruwe tijdssporen.
• 1D Anti-ferromagnetische Heisenberg-keten ($N=64$): Gebruikmakend van een via DMRG voorbereide grondtoestand, vergelijken de auteurs PP met TDVP. Hoewel agressieve PP-afkap leidt tot significante frequentieverbreding en verminderde resolutie vergeleken met MPS-TDVP, legt de methode de generieke kenmerken van het spinon-continuüm succesvol vast. De studie benadrukt dat voor verstrengelde grondtoestanden de bottleneck verschuift van operator-evolutie naar de evaluatie van verwachtingswaarden.
• 2D Anti-ferromagnetische Heisenberg-rooster ($8 \times 8$): Dit vertegenwoordigt een regime waar standaard TNS-methoden te duur worden. Door PP te combineren met CAMPS (bond-dimensie $D=240$) en een rank-1 PSD-extensie, berekenen de auteurs de dynamische spin-structuurfactor $S(k, \omega)$ met slechts 7 tijdstappen. Het resulterende spectrum komt goed overeen met de lineaire spin-golf theorie en identificeert correct de Goldstone-modus bij $\Gamma$, de dispersie naar $X$, en de verzachting bij $M$.
• Schaalbaarheid: De auteurs merken op dat voor het 2D-geval de beperking niet de onbegrensde groei van operatoren op latere tijdstippen is (wat wordt vermeden door de korte-tijd extensie), maar eerder de voorbereiding van de grondtoestand en de evaluatie van verwachtingswaarden. Ze demonstreren verder in de appendix dat de methode het magnon-spectrum voor een $100 \times 100$ ferromagnetisch rooster kan herstellen.

Betekenis en Claims
Het artikel claimt dat deze gecombineerde aanpak het regime verbrekt waarin klassieke methoden betrouwbaar de real-time dynamica van interagerende quantum veel-deeltjessystemen kunnen verkennen, met name in dimensies hoger dan één.

• Overwinnen van Exponentiële Groei: De methode omzeilt de exponentiële proliferatie van Pauli-operatoren door de directe simulatie te beperken tot een kort tijdsvenster waar de ondersteuning van de operator beheersbaar is, en vervolgens de data uit te breiden met behulp van fysieke restricties (positiviteit en low-rank frequentie-dominantie).
• Efficiëntie in Hoge Dimensies: Door CAMPS te integreren, demonstreren de auteurs dat evenwichtscorrelaties voor 2D-systemen berekend kunnen worden aan de rand van wat hanteerbaar is met andere methoden, waarbij de steile groei van de bond-dimensie die gepaard gaat met standaard MPS-tijdsevolutie wordt omzeild.
• Spectrale Getrouwheid: De resultaten wijzen erop dat hoog-resolutie spectrale informatie kan worden geëxtraheerd uit extreem beperkte tijdsdata (zo weinig als 7 tijdstappen in het 2D-geval), mits de dynamica van het systeem wordt gedomineerd door een beperkt aantal frequenties.

De auteurs erkennen bescheiden de beperkingen: de low-rank aanname faalt voor systemen met een continuüm van frequenties (zoals het 1D spinon-continuüm), en het succes van de methode hangt ervan af of de afkaping voldoende nauwkeurig is om de onderliggende frequenties te identificeren. Zij positioneren dit werk niet als een vervanging voor tensornetwerken in 1D, maar als een levensvatbaar, veelbelovend pad voor het berekenen van dynamische correlatiefuncties in grote, hoog-dimensionale roosters waar de huidige state-of-the-art methoden moeite mee hebben.