Efficient classical simulation of two-dimensional long-range systems: Rydberg arrays and beyond

Technische Samenvatting: Efficiënte Klassieke Simulatie van Tweedimensionale Lange-Afstandssystemen

Probleemstelling
Variational Monte Carlo (VMC) is een primaire tool voor het benchmarken van kwantumsimulatoren, met name voor tweedimensionale (2D) systemen met lange-afstandinteracties (bijv. Rydberg-atoomarrays met $1/r^3$ of $1/r^6$ koppelingen). Er bestaat echter een aanzienlijke computationele bottleneck bij het evalueren van de "lokale energie" ($H_{\text{loc}}$) voor systemen met willekeurige all-to-all twee-lichaam-interacties.

• In standaard VMC, terwijl de bemonsteringskosten $O(N)$ zijn voor Tensor Network States (TNS) en $O(N^2)$ voor Neural Quantum States (NQS), vereist de evaluatie van $H_{\text{loc}}$ voor lange-afstand off-diagonale koppelingen het sommeren van $O(N^2)$ niet-lokale termen.
• Omdat elke term een volledige netwerkevaluatie vereist die $O(N)$ tijd kost, schaalt de totale kosten per Monte Carlo (MC) sweep als $O(N^3)$. Deze kubische schaling maakt nauwkeurige klassieke simulatie van grote 2D lange-afstandssystemen (bijv. $10 \times 10$ roosters) onuitvoerbaar.
• Daarnaast ontstaat er een technische pathologie wanneer TNS wordt geëvolueerd vanuit producttoestanden (gebruikelijk in Rydberg-experimenten): bij exacte producttoestanden ontwikkelt de Quantum Geometry Tensor (QGT) nul-richtingen als de bindingsdimensie groter is dan één, waardoor de Time-Dependent Variational Principle (TDVP) vergelijkingen singulier worden.

Methodologie
De auteurs introduceren een adaptieve-basis tensor-netwerk tVMC-framework dat is ontworpen om de computationele complexiteit van de lokale energie-evaluatie te reduceren van $O(N^3)$ naar $O(N)$, terwijl het vermogen behouden blijft om algemene inhomogene koppelingen en realistische geometrieën aan te kunnen.

1. Adaptieve Basis Bemonstering:

   • In plaats van te bemonsteren in een enkele vaste basis (bijv. de $z$-basis), wordt de Hamiltoniaan ontleed in componenten die diagonaal zijn in specifieke bases.
   • Voor het dipolaire XY-model ($H = H_{XY} + H_Z$) splitsen de auteurs de evolutie op in $H_X$ (diagonaal in de $x$-basis) en $H_Y + H_Z$ (diagonaal in de $y$-basis).
   • Er worden twee onafhankelijke Markov-ketens bijgehouden: één die bemonsteren uit $|\langle x|\psi_\theta\rangle|^2$ en een andere uit $|\langle y|U|\psi_\theta\rangle|^2$, waarbij $U$ de basisveranderingsoperator is.
   • Omdat de off-diagonale lange-afstandkoppelingen binnen elke specifieke basis verdwijnen, wordt de lokale energie-evaluatie voor elke component $O(N)$.
   • De parameterupdates worden apart berekend voor elke basis ($\dot{\theta}_X$ en $\dot{\theta}_Y$) en gecombineerd via $\dot{\theta} = \dot{\theta}_X + U^\dagger \dot{\theta}_Y$.

2. Afhandeling van Producttoestand Initiële Condities:

   • Om de singulariteit in de QGT op te lossen wanneer men start vanuit een producttoestand (waar de bindingsdimensie $D > 1$ is), hanteren de auteurs een Projected tVMC (p-tVMC) benadering voor de initiële tijdstap.
   • In plaats van de singuliere TDVP-vergelijking op te lossen, wordt de eerste stap ($\delta t$) bepaald door de fidelity $|\langle \psi_\theta | V | \psi(0) \rangle|^2$ te maximaliseren, waarbij $V$ de tweede-orde expansie van de tijdsevolutie-operator is.
   • Deze optimalisatie is efficiënt omdat de initiële toestand een producttoestand is, waardoor de bemonsteringsverdeling $|\langle s|V|\psi(0)\rangle|^2$ computationeel hanteerbaar is. Na deze initiële stap wordt de standaard tVMC hervat.

3. Exploitatie van Symmetrie:

   • Het framework legt expliciet reflectiesymmetrie langs de roosterdiagonalen op in de TNS-ansatz. Dit vermindert het aantal onafhankelijke variationele parameters met een factor vier, wat de GPU-geheugeneisen aanzienlijk verlicht en de variantie van parameterupdates vermindert.

Belangrijkste Resultaten
De methode werd toegepast om de adiabatische real-time protocol te simuleren van een $10 \times 10$ dipolair XY-model, gerealiseerd in een recent Rydberg-experiment (Chen et al., Nature 2023), een systeemgrootte die voorheen buiten het bereik van klassieke simulatie lag.

• Benchmarking: Op een $4 \times 4$ systeem convergeren de TNS-tVMC resultaten (met bindingsdimensies $D=4, 5$) systematisch naar exact diagonalisatie-resultaten, waarbij gedreven relaxatiedynamica en snelle oscillaties accuraat worden gereproduceerd.
• $10 \times 10$ Simulatie:
  • De simulatie volgde de energie-evolutie en toonde aan dat het systeem relaxeert naar een bijna stationair plateau dat iets boven de grondtoestandsenergie ligt.
  • Antiferromagnetische (AFM) Correlaties: De studie adresseerde een discrepantie in het originele experiment, waarbij de dynamisch voorbereide toestand een onderdrukte lange-afstand AFM-correlatie vertoonde vergeleken met grondtoestand Matrix Product State (MPS) berekeningen.
  • De ideale dynamica-simulatie van de auteurs (startend vanaf een perfecte Néel-toestand met een grotere detuning $\delta_0 = 19.48$) toonde aan dat het voorgeschreven adiabatische ramp proces wel ruimtelijk uitgebreide AFM-correlaties produceert die consistent zijn met de variationele grondtoestand.
  • Dit suggereert dat de onderdrukking die in het experiment werd waargenomen, waarschijnlijk te wijten is aan experimentele imperfecties (zoals beperkt laservermogen dat het gebruik van de optimale $\delta_0$ op de $10 \times 10$ array verhinderde) in plaats van een intrinsiek falen van het adiabatische protocol.
• Convergentie: De resultaten toonden een robuuste convergentie met betrekking tot de steekproefgrootte ($N_s = 81.920$ vs. $163.840$), wat de betrouwbaarheid van de adaptieve-basis benadering bevestigt.

Betekenis en Claims
Het artikel claimt een schaalbare klassieke benchmarking-tool te bieden voor lange-afstand kwantumplatforms, specifiek Rydberg en ionen-val simulatoren, die momenteel een dringende behoefte hebben aan dergelijke tools.

• Technische Prestatie: Het werk lost de $O(N^3)$ bottleneck op, wat nauwkeurige simulatie van real-time dynamica en grondtoestanden voor 2D lange-afstandssystemen mogelijk maakt met een $O(N)$ kost per MC-sweep.
• Resolutie van Pathologie: Het biedt een robuuste oplossing voor het evolueren van TNS vanuit producttoestanden, een veelvoorkomend scenario in experimentele protocollen.
• Wetenschappelijke Impact: Door een directe klassieke referentie te bieden voor het $10 \times 10$ Rydberg-experiment, verheldert dit werk de aard van de discrepantie tussen de experimentele resultaten en de theoretische grondtoestand-voorspellingen. Het demonstreert dat het ideale adiabatische protocol in staat is om de verwachte lange-afstand orde te genereren, waardoor de focus van experimentele optimalisatie verschuift naar realisatie-imperfecties.
• Toekomstperspectief: De auteurs positioneren deze methode als een brug voor klassiek-kwantum samenwerking, waarbij klassieke simulaties kunnen helpen bij het identificeren van interessante fysieke opstellingen voordat de entanglement-groei te groot wordt voor tensor netwerken, waarna de dynamica kan worden overgedragen aan kwantumsimulatoren. De methode wordt ook genoteerd als een nieuwe weg voor het bestuderen van 2D lange-afstandssystemen met sign-problemen, waar andere numerieke methoden moeite mee hebben.