Hamiltonian spectra in quantum computers through the generalized eigenvalue method

Technische Samenvatting: Hamiltoniaanse Spectra in Kwantumcomputers via de Gegeneraliseerde Eigenwaarde Methode

Probleemstelling
Het bepalen van het energiediscretum van een Hamiltoniaan is een fundamenteel doel in de kwantumfysica, met name voor sterk gekoppelde veldentheorieën waar perturbatieve methoden falen. Ho'e de klassieke roosterveldtheorie maakt gebruik van Monte Carlo-evaluaties van Euclidische correlatoren om spectra te extraheren, maar deze aanpak staat voor een "mixed blessing": hoewel de grondtoestand domineert bij grote imaginaire tijd, worden geëxciteerde toestanden exponentieel onderdrukt, wat de extractie ervan bemoeilijkt door de degradatie van de signaal-ruisverhouding.

Kwantumcomputers bieden een duidelijk voordeel door het natuurlijk genereren van real-time correlatoren, waarbij alle toestanden bijdragen met ongedempte, oscillerende gewichten. Het extraheren van energie-eigenwaarden uit deze oscillerende signalen berust echter doorgaans op Fourier-transformaties. De resolutie van een Fourier-transformatie wordt beperkt door de inverse van de totale evolutietijd ($1/T$). Huidige kwantumbo hardware, die beperkt wordt door eindige coherentietijden en accumulerende gate- en Trotter-fouten, kan niet de lange evolutietijden handhaven die nodig zijn om bijna-gedegenereerde niveaus te resolven of om echte pieken te onderscheiden van ringing-artefacten via Fourier-analyse.

Methodologie
De auteurs stellen voor om de GEVP (Generalized Eigenvalue Problem) — een standaardtechniek in de Euclidische roosterveldtheorie — aan te passen aan real-time correlatoren gegenereerd op kwantumcomputers. De kernmethodologie omvat:

1. Operatorbasis en Correlatierestrictie: In plaats van één enkele correlator te analyseren, construeert de methode een matrix van correlatoren $C_{ij}(t) = \langle \Omega | O_i(t) O_j^\dagger(0) | \Omega \rangle$, waarbij $O_i$ een set imperfecte interpolerende operatoren is en $|\Omega\rangle$ een benaderde grondtoestand.
2. Theoretisch Kader: De auteurs demonstreren dat de eigenwaarden $\lambda_k(t)$ van de matrixproduct $C^{-1}(0)C(t)$ zich gedragen als $e^{-iE_k t}$, waarbij $E_k$ de energie-eigenwaarden zijn. De effectieve energie wordt geëxtraheerd via $E_k^{\text{eff}}(t) = \frac{i}{t} \log(\lambda_k(t))$. Theoretische analyse toont aan dat correcties aan de exacte energieën vervallen als $1/t$ (of sneller door destructieve interferentie), mits de operatorbasis de relevante subruimte van eigentoestanden voldoende overspant.
3. Kwantumimplementatie: Een specifiek algoritme wordt gedetailleerd voor het meten van off-diagonaal correlatoren van Pauli-operatoren op een kwantumcomputer met behulp van een enkele ancilla-qubit. De procedure omvat:
   • Het prepareren van de benaderde grondtoestand $|\Omega\rangle$ en een ancilla.
   • Het verstrengelen van het systeem met de ancilla via een gecontroleerde operatie die de operator $O_j$ bevat.
   • Het tijdsevolueren van het systeem.
   • Het meten van specifieke observabelen ($X \otimes O_i$ en $Y \otimes O_i$) om de reële en imaginaire delen van de correlator te reconstrueren.
4. Vergelijkingsbenchmark: De methode wordt gebenchmarkt tegen de discrete Fourier-transformatie (DFT) van de trace van de correlatiematrix, waarbij een dempingsfactor wordt gebruikt om Gibbs-fenomenen te mitigeren.

Belangrijkste Bijdragen en Resultaten
Het artikel valideert deze aanpak met behulp van de fuzzy $\sigma$-model (equivalent aan de Heisenberg-comb), een theorie die kenmerken deelt met fysiek relevante veldentheorieën zoals asymptotische vrijheid. De studie maakt gebruik van drie verschillende testregimes:

• Exacte Simulaties ($L=8$): Met behulp van exacte tijdsevolutie (geen ruis of Trotter-fouten), slaagde de methode erin om 8, 36 en 276 energieniveaus te resolven met drie steeds complexere operator-sets. De resultaten bereikten een precisie van minder dan één procent voor de laagste toestanden binnen een maximale evolutietijd van $T_{\text{max}} = 5.6$. De methode reproduceerde correct de degeneratiepatronen die door de $O(3)$-symmetrie van het model worden bepaald.
• Gate-Based Simulaties ($L=4$): De auteurs analyseerden de impact van Trotter-fouten. Door resultaten te extrapoleren van verschillende tijdstappen ($\Delta t$) naar $\Delta t \to 0$, bevestigden zij dat de methode robuust is, waarbij afwijkingen kwadratisch schalen met de stapgrootte.
• Kwantumhardware ($L=4$ op IonQ Forte-Enterprise): De methode werd uitgevoerd op werkelijke kwantumbo hardware. Met behulp van een enkele ancilla-qubit en een minimale operator-set, losten de auteurs vier laagliggende energieniveaus op. De geëxtraheerde energieën kwamen binnen enkele procenten overeen met de exacte diagonalisatie. De resterende discrepanties werden toegeschreven aan een combinatie van Trotter-fouten (door de grove $\Delta t = 0.4$) en hardware-ruis/decoherentie, in plaats van een falen van de GEVP-methode zelf.

Betekenis en Claims
Het artikel claimt dat de GEVP-aanpak "substantieel efficiënter" is dan Fourier-transformatie methoden voor het extraheren van spectra uit real-time kwantumsimulaties. Specifiek maakt het de resolutie van laagliggende energieniveaus mogelijk bij tijdsbereiken waar de Fourier-transformatie faalt om pieken van artefacten te onderscheiden.

De auteurs benadrukken dat de methode geen perfecte operatoren of een exacte grondtoestand vereist; het vertrouwt enkel op een operatorbasis die de toestanden van belang benaderend overspant, en een benaderde grondtoestand met een goede overlap. Hoewel dit hier is gedemonstreerd met sterke koppelingslimieten om de noodzakelijke toestanden te construeren, merken de auteurs op dat de methode algemeen is en toepasbaar op elke variationele of adiabatische preparatie van toestanden.

De primaire betekenis die wordt benadrukt, is het potentieel van deze techniek om spectroscopie mogelijk te maken in theorieën waar het "tekenprobleem" (sign problem) Euclidische benaderingen blokkeert, waarbij specif