Quantum Imaginary Time Evolution on an Infinite 1D Chain

Technische Samenvatting: Kwantum Imaginaire Tijdsevolutie op een Oneindige 1D-keten

Probleemstelling
De komst van Noisy Intermediate-Scale Quantum (NISQ)-apparaten heeft de belangstelling voor het simuleren van kwantumveel-deeltjessystemen aangewakkerd. Hoewel algoritmen voor real-time evolutie (bijv. TEBD, TDVP) succesvol zijn aangepast aan kwantumcircuits door gebruik te maken van de unitaire aard van kwantumpoorten, blijft het uitbreiden van deze methoden naar Imaginaire Tijdsevolutie (ITE) uitdagend. De ITE-operator, $e^{-H\tau}$, is niet-unitair, wat incompatibel is met standaard unitaire kwantumpoorten. Bestaande benaderingen voor dit probleem—Quantum ITE (QITE), Variational ITE (VITE) en Probabilistic ITE (PITE)—zijn primair beperkt tot eindige systemen. Bovendien, op echte kwanthardware, heeft de statistische ruis die inherent is aan het schatten van kostenfuncties via herhaalde metingen een significante impact op de convergentie van variationele optimalisatie. Dit werk adresseert de kloof in het uitvoeren van ITE op oneindige eendimensionale rooster-systemen met behulp van kwantumcircuits en analyseert de statistische effecten van eindige-steekproefruis op convergentie.

Methodologie
De auteurs stellen een kwantumcircuit-algoritme voor dat gebaseerd is op het Time-Dependent Variational Principle (TDVP) toegepast op een uniforme Matrix Product State (uMPS) ansatz. De kernmethodologie omvat:

1. Mapping van uMPS naar Kwantumcircuits:

   • Een translationeel invariante uMPS wordt gemapt naar een geparametriseerd kwantumcircuit (PQC). Voor een bond-dimensie $D=2^N$ mapt een enkele-site unit cell naar een unitaire poort werkend op $N+1$ qubits.
   • Om de niet-unitaire omgevingstensors ($l$ en $r$) te verwerken die vereist zijn voor de TDVP-kostenfunctie, construeren de auteurs unitaire poorten door de Hilbertruimte te vergroten (het inbedden van niet-unitaire tensors in grotere unitaire matrices).
   • De kostenfunctie, die de overlap tussen de geëvolueerde toestand en de doeltoestand vertegenwoordigt, wordt geformuleerd als een kwantumcircuit waarbij de output-waarschijnlijkheid van de all-zero toestand dient als metriek voor optimalisatie.

2. Afhandeling van Niet-Unitariteit (Twee Strategieën):

   • Unitaire Approximatie (QITE-achtig): De niet-unitaire ITE-operator wordt benaderd door een unitaire operator $e^{-iA(\tau)d\tau}$, waarbij $A(\tau)$ een Hermitische operator is die wordt bepaald door de afstand tussen de genormaliseerde ITE-toestand en zijn unitaire benadering te minimaliseren. Deze benadering vereist geen extra qubits maar introduceert approximatiefouten.
   • Ancilla-ondersteunde Inbedding (PITE-achtig): De niet-unitaire operator wordt ingebed in een grotere unitaire operator met behulp van een extra ancilla-qubit via block encoding. Deze methode vermijdt approximatiefouten maar verhoogt het aantal qubits (bijv. van 6 naar 7 qubits voor $D=2$).

3. Quantum Lanczos (QLanczos) Verfijning:

   • Om fouten van unitaire approximaties en statistische ruis te mitigeren, gebruiken de auteurs het Quantum Lanczos-algoritme.
   • In plaats van een Krylov-subruimte te genereren uit $H^k|\phi\rangle$, construeert QLanczos een subruimte uit de sequentie van ITE-toestanden $|\phi_l\rangle = e^{-l d\tau H}|\phi_0\rangle$.
   • Door een gegeneraliseerd eigenwaardeprobleem op te lossen met behulp van de overlapmatrix $S$ en de Hamiltoniaanse matrix $H$ geconstrueerd uit deze toestanden, verfijnt het algoritme de grondtoestand-energie schatting.

Belangrijkste Resultaten
De auteurs hebben hun algoritme getest met het Transverse-Field Ising Model (TFIM) op zijn kritieke punt ($H = \sum_i (-\sigma^z_i \sigma^z_{i+1} + \sigma^x_i)$) met een bond-dimensie van $D=2$.

• Statevector Simulaties:
  • De ancilla-ondersteunde methode kwam nauw overeen met klassieke Tensor Network (TN) TDVP resultaten en convergeerde naar de verwachte grondtoestand-energiedichtheid ($\approx -1.2732$).
  • De unitaire approximatiemethode vertoonde een onverwachte toename in energie bij lange evolutietijden ($\tau \approx 1$), toegeschreven aan de beperkte correlatie-domeingrootte inherent aan de oneindige-systeem opzet en de unitaire approximatiefout.
• Simulaties op Echte Kwantumapparaten (IBM Quantum):
  • Simulaties op apparaten (bijv. ibm_nairobi, ibm_auckland) toonden aan dat hoewel ruwe ITE-trajecten ruizig zijn vanwege eindige-steekproef statistiek, de QLanczos post-processing de fouten aanzienlijk verminderde, waardoor de energie-schattingen dichter bij de grondtoestandwaarde kwamen.
• Statistische Analyse van Ruis:
  • De studie analyseerde de distributie van kostenwaarden en energiedichtheden afgeleid van meting-gebaseerde schattingen.
  • Hoewel individuele kosten-schattingen een normale distributie volgen (Bernoulli-proces), vormt de geaggregeerde distributie van energiedichtheden over meerdere optimalisatieruns een rechts-scheve, Maxwell-Boltzmann-achtige envelop.
  • Deze distributie ontstaat door een combinatie van run-to-run variaties in "stall points" (waar optimalisatie stopt door ruis) en Gaussische fluctuaties binnen individuele runs. De analyse geeft aan dat optimalisatietrajecten vaak stoppen voordat ze het variationele minimum bereiken wanneer de optimizer tolerantie vergelijkbaar wordt met de eindige-steekproefruis.

Betekenis en Claims
Het artikel claimt de eerste implementatie van kwantumcircuit ITE voor oneindige 1D-systemen te demonstreren. De primaire bijdragen zijn:

1. Algoritmische Uitbreiding: Succesvolle aanpassing van het TDVP-raamwerk aan oneindige roosters met behulp van geparametriseerde kwantumcircuits, waarbij een constant aantal poorten behouden blijft ongeacht de evolutietijd.
2. Ruiskarakterisering: Het bieden van een gedetailleerde statistische studie over hoe eindige-steekproefruis de variationele convergentie in ITE beïnvloedt, inclus\u00ad bij de specifieke distributionele kenmerken van de resulterende energie-fouten.
3. Foutmitigatie: Het valideren van het nut van het QLanczos-algoritme voor het verfijnen van grondtoestand-energie schattingen op zowel simulaties als echte hardware, wat effectief compenseert voor zowel unitaire approximatiefouten als meetruis.
4. Schaalbaarheid: Hoewel de huidige implementatie zich richt op $D=2$ (wat slechts 6 of 7 qubits vereist), merken de auteurs op dat de mapping tussen uMPS en kwantumpoorten kan worden gegeneraliseerd naar grotere bond-dimensies, wat een pad biedt om complexere, sterk verstrengelde systemen in de toekomst te simuleren.

Het werk concludeert dat hoewel huidige NISQ-apparaten aanzienlijke statistische uitdagingen introduceren, de combinatie van variationele ansätze, ancilla-ondersteunde inbedding en QLanczos-verfijning een levensvatbare route biedt voor het simuleren van oneindige kwantumsystemen en het schatten van grondtoestand-eigenschappen.