Spectral Gaps via Imaginary Time

Technische Samenvatting: Spectrale Gaten via Imaginaire Tijd

Probleemstelling
De spectrale kloof, gedefinieerd als $\Delta E = E_1 - E_0$ (het energieverschil tussen de grondtoestand en de eerste aangeslagen toestand), is een kritische parameter in diverse fysieke contexten, variërend van fluorescentielengtes en excitonische energietransfer in materialen tot adiabatische looptijden en state-transfer fideliteiten in quantuminformatie. Het bepalen van deze kloof is computationeel uitdagend; in het slechtste geval is het beslissen of een lokale Hamiltoniaan een kloof heeft onbeslisbaar, en het schatten van laagliggende spectra is over het algemeen NP-hard. Bestaande numerieke methoden, zoals Fermionische Quantum Monte Carlo (die lijdt onder het tekenprobleem) en Tensor-Netwerk methoden (beperkt door de bond-dimensie), kennen significante afruilbetrachtingen met betrekking tot systeemgrootte, basisgrootte en meetoverhead. Terwijl Imaginaire-Tijd Evolutie (ITE) een robuuste klassieke tool is voor het vinden van grondtoestanden, blijft het extraheren van specifieke spectrale kloven zonder de volledige dichtheid van toestanden te resolveren of quantumfase-estimatie uit te voeren, een specifieke uitdaging.

Methodologie
De auteurs stellen een methode voor om de spectrale kloof $\Delta E$ direct te schatten vanuit de ratio van twee verwachtingswaarden van geneste commutatoren, geëvalueerd langs een imaginaire-tijd traject.

1. Theoretisch Kader:
   Voor een einddimensionele Hamiltoniaan $H$ met spectrale decompositie $H = \sum E_n \Pi_n$, en een initiële toestand $|\phi_0\rangle$ met ondersteuning op zowel de grondtoestand ($\Pi_0$) als de eerste aangeslagen subspace ($\Pi_1$), evolueert het systeem in imaginaire tijd $\tau$ als $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$.
   De methode maakt gebruik van een "coördinerende" observabele $O$ die de grondtoestand en de eerste aangeslagen subspace koppelt. De auteurs definiëren de $M$-de geneste commutator $[H, O]_M$.
   Het kernresultaat (Theorem 1) demonstreert dat naarmate $\tau \to \infty$, de ratio van de verwachtingswaarden convergeert naar het kwadraat van de energie-gap:
   $$\frac{\langle [H, O]_{M+2} \rangle_\tau}{\langle [H, O]_M \rangle_\tau} = (E_0 - E_1)^2 + \mathcal{O}(e^{-\tau(E_2 - E_1)})$$
   Deze estimator is agnostisch voor degeneraties in hogere aangeslagen niveaus omdat het steunt op projectorformuleringen in plaats van individuele eigenvectoren.

2. Logaritmische Estimators:
   Het paper leidt ook logaritmische estimators af om $E_0 + E_1$ te extraheren (uit de helling van $\ln|\langle [H, O]_M \rangle_\tau|$ versus $\tau$) en de tweede kloof $E_2 - E_1$ (uit de residu van de ratio-fit), wat alternatieven biedt wanneer numerieke deling van kleine getallen instabiel is.

3. Quantum Implementatie (ITQDE):
   Om dit op quantumhardware te implementeren, waar unitaire evolutie inherent is maar imaginaire tijd niet, maken de auteurs gebruik van Imaginary-Time Quantum Dynamical Emulation (ITQDE). Ze gebruiken een Hubbard-Stratonovich transformatie om de Euclidische kernel $e^{-\tau H^2}$ uit te drukken als een Gaussische gemiddelde van real-time unitaries. Deze kernel wordt gediscretiseerd met behulp van Gauss-Hermite kwadratuur, waarbij de evolutie wordt benaderd als een gewogen som van unitaries.
   De verwachtingswaarden worden vervolgens gereduceerd tot gewogen sommen van de reële delen van overlaps van de vorm $\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$, waarbij $P_\alpha$ Pauli-strings zijn. Deze overlaps kunnen worden gemeten met standaard Hadamard-tests. De auteurs merken op dat vanwege de symmetrie van de kernel, alleen de reële delen van deze overlaps vereist zijn.

Belangrijkste Resultaten
De auteurs hebben hun methode getest middels klassieke simulaties (via de QuSpin library) op twee standaardmodellen:

• 1D Transverse-Field Ising Model (TFIM): Voor systeemgroottes $L=4$, nam de relatieve fout in de geschatte kloof exponentieel af met $\tau$ (bereikde $\sim 10^{-8}$) voordat deze bij zeer grote $\tau$ weer toenam door de onderdrukking van de eerste aangeslagen toestand onder de numerieke ruis.
• 1D Fermi-Hubbard Model: Een vergelijkbare exponentiële afname in fout werd waargenomen voor een half-gevulde keten ($L=4$), waarbij relatieve fouten van $\sim 10^{-4}$ werden bereikt. De methode handelde de interagerende fermionische setting met niet-triviale spin- en ladingsstructuren succesvol af.
• Quantum Simulatie: Een proof-of-principle simulatie werd uitgevoerd op de TFIM met de Qiskit Aer backend met $10^6$ shots. Hoewel de methode compatibel was met het ITQDE-framework, was de minimale fout die werd bereikt $\sim 10^{-1}$. De auteurs schrijven dit toe aan de "bandwidth tradeoff": de Gaussische filter die nodig is om hogere toestanden te onderdrukken, drijft ook de relevante overlaps richting de shot-noise floor als de imaginaire tijd $\tau$ en de kwadratuurorde $K$ niet zorgvuldig worden afgestemd.

Betekenis en Claims
Het paper claimt een "bescheiden maar genuien voordeel in structuur" te bieden voor spectrale schatting.

• Eenvoud en Omvang: De methode biedt een manier om lage-energie kloven te schatten zonder individuele eigenfuncties te resolveren, quantumfase-estimatie uit te voeren, of de volledige dichtheid van toestanden te reconstrueren. Het zet spectrale vragen om in vragen over relaxatie langs een één-parameter flow.
• Robuustheid: De estimator is robuust tegen degeneraties in hogere energieniveaus en is toepasbaar op zowel klassieke ijle lineaire algebra als quantum sum-of-terms Hamiltoniaanse.
• Bescheidenheid: De auteurs stellen expliciet dat deze benadering de worst-case complexiteit van het lokale Hamiltoniaan-probleem niet omzeilt. Het is geen vervanging voor volledige spectrale reconstructie of fase-estimatie wanneer een gedetailleerde dichtheid van toestanden vereist is.
• Toekomstverwachting: Dit werk dient als een proof-of-principle. De auteurs merken op dat het concurrerend maken van de methode op quantumhardware vereist dat de volledige bandwidth-analyse uit eerder werk wordt geïmporteerd om parameters ($\tau, K$) te optimaliseren, betere coördinerende observabelen worden ontwor than, en meer meet-efficiënte overlap-estimators worden toegepast. Het paper positioneert deze kloof-georiënteerde perspectief als een tool die aansluit bij bredere programma's die effectieve Hamiltoniaanse gebruiken om differentiaalvergelijkingen op te lossen, waarbij relaxatiesnelheden en kloven vaak de primaire grootheden van belang zijn.