Exploring the Relaxation Landscape of a 2D Quantum Magnet on a 256-Qubit Processor

Technische Samenvatting: Verkenning van het Relaxatielandschap van een 2D Kwantummagneet op een 256-Qubit Processor

Probleemstelling
De relaxatie van kwantummaterie ver van evenwicht blijft een centraal openstaand probleem in de veel-deeltjesfysica. Hoewel analoge kwantumsimulatoren er succesvol in zijn geslaagd om theoretische voorspellingen over universele dynamische fenomenen (zoals dynamische faseovergangen en emergente hydrodynamica) te reproduceren, is hun potentieel om te fungeren als een werkelijke ontdekkingsplatform voor fysica die nog niet door de theorie is voorzien, een front dat pas net wordt verkend. Specifiek zijn de relaxatiepaden van twee-dimensionale (2D) kwantumsystemen minder goed begrepen dan die van hun 1D-tegenhangers. De auteurs onderzoeken de post-quench dynamica van het 2D transversale-veld Ising-model (TFIM), een veeleisende testomgeving die zowel thermische als kwantumfaseovergangen kenmerkt, om het relaxatielandschap in kaart te brengen en regimes te identificeren waarin klassieke simulatiemethoden falen.

Methodologie
De studie maakt gebruik van een programmeerbare analoge kwantumverwerkingseenheid (QPU) gebaseerd op een 256-qubit ($16 \times 16$) array van gevangen $^{87}$Rb-atomen. Lokale qubits worden gecodeerd in de grondtoestand ($|\downarrow\rangle$) en een hoog-geëxciteerde Rydberg-toestand ($|\uparrow\rangle$).

• Hamiltoniaan-simulatie: Het systeem simuleert de 2D TFIM-Hamiltoniaan $\hat{H}_{TFIM} = -\sum_{i<j} J_{ij}\hat{\sigma}^z_i \hat{\sigma}^z_j - h_x \sum_i \hat{\sigma}^x_i$. De interactie $J_{ij}$ komt voort uit kort-bereik van der Waals-krachten ($U \propto 1/r^6$), terwijl het transversale veld $h_x$ wordt gecontroleerd via de globale Rabi-frequentie $\Omega$.
• Protocol: Het systeem wordt geïnitialiseerd in een volledig gepolariseerde producttoestand $|\downarrow \dots \downarrow\rangle$ en gekwantumt (quenched) naar een constante $h_x$ en $J$. De ratio $h_x/J$ wordt afgestemd om over het fasediagram te scannen, variërend van de ferromagnetische (FM) fase naar de paramagnetische (PM) fase.
• Meting en Mitigatie: Dynamica wordt geprobeerd via projectieve metingen van $\hat{\sigma}^z$ met behulp van 300 experimentele shots per tijdstap. Gegevens worden gecorrigeerd voor state preparation and measurement (SPAM) fouten ($\epsilon \approx 1\%$, $\epsilon' \approx 3\%$) en gemiddeld over equivalente qubits.
• Benchmarking: De QPU-resultaten worden vergeleken met de meest geavanceerde klassieke tensor-netwerkmethoden, inclusen 2D Tensor Netwerken (2DTN), Matrix Product States (MPS) en Tree Tensor Networks (TTN), evenals Quantum Monte Carlo (QMC) voor thermische evenwichtseigenschappen.

Belangrijkste Bijdragen en Resultaten
De auteurs identificeren drie onderscheidende regimes in het relaxatielandschap van de 2D TFIM, gekenmerkt door het gedrag van lokale klassieke energie en verbonden correlaties:

1. Snelle Thermalisatie (Nabij de Dynamische Faseovergang):
   Voor quenches nabij het kritieke punt ($h_x/J \approx 2.5$), vertoont het systeem snelle thermalisatie. Correlaties relaxeren naar hun thermische verwachtingswaarde op tijdschalen van $t \approx 1/J$. De bulk-magnetisatie behoudt het teken van de initiële toestand in de FM-regio, consistent met symmetriebreking, terwijl deze rond nul oscilleert in de PM-regio, consistent met symmetrieherstel.

2. Pretherm Regime (Groot Transversaal Veld):
   Bij een grote $h_x/J$ (bijv. $h_x/J = 6.0$), betreedt het systeem een langdurig pretherm plateau. De dynamica wordt beheerst door een effectieve Floquet-Hamiltoniaan afgeleid in het roterende kader, die de vorm aanneemt van een $U(1)$-symmetrisch $ZY$-model: $\hat{H}_{ZY} = \frac{1}{2}\sum_{i<j} J_{ij}(\hat{\sigma}^z_i \hat{\sigma}^z_j + \hat{\sigma}^y_i \hat{\sigma}^y_j)$.

• Observatie: In plaats van te relaxeren naar de thermische waarde van de oorspronkelijke TFIM, worden correlaties in stand gehouden op een hogere prethermische waarde die wordt voorspeld door de thermische eigenschappen van $\hat{H}_{ZY}$.
• Correlaties: Dit regime kenmerkt zich door een sterke dynamische opbouw van verbonden correlaties die zich over grote afstanden uitstrekt, consistent met een quasi-lange-afstandsgeordende toestand bij een lage effectieve temperatuur.

3. Langzame Relaxatie Crossover Regime:
   In het intermediaire regime (bijv. $h_x/J = 3.5$), tussen de thermische en prethermische limieten, ontdekren de auteurs een onverwachte vertraging in de thermalisatie.

• Observatie: Correlaties vervallen traag naar hun thermische waarden en blijven ver van evenwicht voorbij de karakteristieke relaxatietijd ($t \gg 1/J$).
• Klassiek Falen: Cruciaal is dat dit het regime is waar de meest geavanceerde klassieke tensor-netwerkmethoden (2DTN, MPS, TTN) de controle verliezen bij late tijden. De klassieke simulaties divergeren van elkaar en falen om een gecontroleerde beschrijving te bieden vanwege de snelle groei van verstrengeling in het 2D-systeem. In tegenstelling hiertoe blijft de QPU-data consistent over verschillende systeemgroottes ($L=8$ en $L=16$) en blijft de langzame verval voortbestaan.

Betekenis en Claims
Het artikel vestigt Rydberg-atoomarrays als een robuust platform voor wetenschappelijke ontdekkingen in de niet-evenwicht kwantum-veel-deeltjesdynamica. De primaire betekenis ligt in de demonstratie dat analoge kwantumsimulatoren toegang hebben tot dynamische regimes die numeriek onberekenbaar zijn voor klassieke computers.

• Ontdekking versus Bevestiging: In tegenstelling tot eerdere studies die grotendeels voorspelde theoretische effecten reproduceerden, identificeert dit werk een "vertraging"-regime dat niet volledig door de huidige theorie wordt verklaard, wat het potentieel van kwantumsimulatoren benadrukt om nieuwe fysica te onthullen.
• Beperkingen van Benchmarking: De studie demonstreert expliciet de beperkingen van de huidige tensor-netwerkbenaderingen in 2D niet-integreerbare systemen bij late tijden, waar de groei van verstrengeling betrouwbare benaderingen met een eindige bond-dimensie onbetrouwbaar maakt.
• Robuustheid: De auteurs merken op dat hoewel experimentele ruis de kwantitatieve relaxatiesnelheden beïnvloedt, het kwalitatieve dynamische gedrag (inclusclusief het bestaan van het langzame-relaxatie-regime) robuust blijft en niet op dezelfde wijze degradeert met de systeemgrootte als fouten in klassieke simulaties (bijv. MPS).

Het werk suggereert dat de trage thermalisatie het gevolg kan zijn van de wisselwerking tussen de breuk van de leidende orde van de Floquet-Magnus expansie en de kritieke vertraging (critical slowing down) nabij de dynamische faseovergang, hoewel een definitieve theoretische verklaring voor dit specifieke crossover-regime een open vraag blijft voor toekomstig onderzoek.