Spectral Gaps via Imaginary Time

Sintesi Tecnica: Gap Spettrali tramite Tempo Immaginario

Problematica
Il gap spettrale, definito come $\Delta E = E_1 - E_0$ (la differenza di energia tra lo stato fondamentale e il primo stato eccitato), è un parametro critico in diversi contesti fisici, che spaziano dalle lunghezze d'onda di fluorescenza e il trasferimento di energia eccitonica nei materiali, fino ai tempi di esecuzione adiabatici e alla fedeltà del trasferimento di stato nella informazione quantistica. Determinare questo gap è computazionalmente impegnativo; nel caso peggiore, decidere se un Hamiltoniano locale sia gappato è indecidibile, e stimare gli spettri a bassa energia è generalmente NP-hard. I metodi numerici esistenti, come il Monte Carlo Quantistico Fermionico (che soffre del problema del segno) e i metodi basati su Reti Tensoriali (limitati dalla dimensione del legame), affrontano significativi compromessi riguardanti la dimensione del sistema, la dimensione della base e l'overhead di misurazione. Sebbene l'Evoluzione in Tempo Immaginario (ITE) sia uno strumento classico robusto per trovare stati fondamentali, estrarre gap spettrali specifici senza risolvere l'intera densità degli stati o eseguire la stima della fase quantistica rimane una sfida specifica.

Metodologia
Gli autori propongono un metodo per stimare il gap spettrale $\Delta E$ direttamente dal rapporto tra due valori di aspettazione di commutatori annidati valutati lungo una traiettoria in tempo immaginario.

1. Quadro Teorico:
   Per un Hamiltoniano $H$ a dimensione finita con decomposizione spettrale $H = \sum E_n \Pi_n$, e uno stato iniziale $|\phi_0\rangle$ con supporto sia sul sottospazio fondamentale ($\Pi_0$) che su quello del primo stato eccitato ($\Pi_1$), il sistema evolve in tempo immaginario $\tau$ come $|\psi(\tau)\rangle = e^{-\tau H}|\phi_0\rangle$.
   Il metodo utilizza un osservabile "coordinante" $O$ che accoppia i sottospazi fondamentale e del primo stato eccitato. Gli autori definiscono il commutatore annidato $M$-esimo $[H, O]_M$.
   Il risultato centrale (Teorema 1) dimostra che per $\tau \to \infty$, il rapporto tra i valori di aspettazione converge al quadrato del gap di energia:
   $$\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)})$$
   Questo stimatore è agnostico rispetto alle degenerazioni nei livelli eccitati superiori perché si basa su formulazioni di proiettore piuttosto che su singoli autovettori.

2. Stimatori Logaritmici:
   Il documento deriva anche stimatori logaritmici per estrarre $E_0 + E_1$ (dalla pendenza di $\ln|\langle [H, O]_M \rangle_\tau|$ rispetto a $\tau$) e il secondo gap $E_2 - E_1$ (dal residuo del fit del rapporto), fornendo alternative quando la divisione numerica di numeri piccoli risulta instabile.

3. Implementazione Quantistica (ITQDE):
   Per implementare il metodo su hardware quantistico, dove l'evoluzione unitaria è nativa ma il tempo immaginario no, gli autori utilizzano l'Emulazione Dinamica Quantistica in Tempo Immaginario (ITQDE). Utilizzano una trasformazione di Hubbard-Stratonovich per esprimere il kernel euclideo $e^{-\tau H^2}$ come una media gaussiana di unitari in tempo reale. Questo kernel viene discretizzato usando la quadratura di Gauss-Hermite, approssimando l'evoluzione come una somma pesata di unitari.
   I valori di aspettazione vengono quindi ridotti a somme pesate di parti reali di sovrapposizioni della forma $\langle \phi_0 | U_j^\dagger P_\alpha U_k | \phi_0 \rangle$, dove $P_\alpha$ sono stringhe di Pauli. Queste sovrapposizioni possono essere misurate utilizzando standard test di Hadamard. Gli autori notano che, a causa della simmetria del kernel, sono richieste solo le parti reali di queste sovrapposizioni.

Risultati Chiave
Gli autori hanno testato il metodo tramite simulazioni classiche (tramite la libreria QuSpin) su due modelli standard:

• Modello di Ising a Campo Trasverso (TFIM) 1D: Per dimensioni di sistema $L=4$, l'errore relativo del gap stimato è diminuito esponenzialmente con $\tau$ (raggiungendo $\sim 10^{-8}$) prima di aumentare nuovamente a $\tau$ molto grandi a causa della soppressione del primo stato eccitato sotto il rumore numerico.
• Modello di Fermi-Hubbard 1D: Un decadimento esponenziale dell'errore è stato osservato similmente per una catena a semipieno ($L=4$), con errori relativi che raggiungono $\sim 10^{-4}$. Il metodo ha gestito con successo l'ambiente fermionico interagente con strutture spin e carica non banali.
• Simulazione Quantistica: Una simulazione di prova di concetto è stata condotta su TFIM utilizzando il backend Aer di Qiskit con $10^6$ shot. Sebbene il metodo fosse compatibile con il framework ITQDE, l'errore minimo raggiunto è stato di $\sim 10^{-1}$. Gli autori attribuiscono ciò al "compromesso della larghezza di banda": il filtro gaussiano necessario per sopprimere gli stati superiori spinge anche le sovrapposizioni rilevanti verso il livello del rumore di shot, se $\tau$ e l'ordine della quadratura $K$ non sono calibrati con cura.

Significato e Rivendicazioni
Il documento sostiene di offrire un "vantaggio strutturale modesto ma genuino" per la stima spettrale.

• Semplicità e Ambito: Il metodo offre un modo per stimare i gap a bassa energia senza risolvere i singoli autostati, eseguire la stima della fase quantistica o ricostruire l'intera densità degli stati. Converte le questioni spettrali in questioni di rilassamento lungo un flusso a un parametro.
• Robustezza: Lo stimatore è robusto contro le degenerazioni nei livelli di energia superiori e si applica sia alla teoria lineare sparsa classica che agli hamiltoniani quantistici a somma di termini.
• Modestia: Gli autori dichiarano esplicitamente che questo approccio non elude la complessità del caso peggiore del problema dell'hamiltoniano locale. Non è un sostituto della ricostruzione spettrale completa o della stima della fase quando è richiesta una densità degli stati dettagliata.
• Prospettive Future: Il lavoro serve come prova di concetto. Gli autori notano che rendere il metodo competitivo sull'hardware quantistico richiede l'importazione dell'intera analisi della larghezza di banda da lavori precedenti per ottimizzare i parametri ($\tau, K$), progettare osservabili coordinanti migliori ed impiegare stimatori di sovrapposizione più efficienti in termini di misurazione. Il documento posiziona questa prospettiva focalizzata sui gap come uno strumento che si integra con programmi più ampi che utilizzano hamiltoniani effettivi per risolvere equazioni differenziali, dove i tassi di rilassamento e i gap sono spesso le quantità primarie di interesse.