Sintesi Tecnica: Un Algoritmo Quantistico Senza Penalità per Trovare Autostati Energetici

Enunciato del Problema

Trovare autovalori e autostati di Hamiltoniani a molti corpi è una sfida fondamentale in fisica e scienza computazionale. I metodi classici affrontano una scalabilità esponenziale (memoria $O(2^{2n})$ e tempo $O(2^{3n})$) all'aumentare del numero di qubit $n$. Sebbene gli algoritmi quantistici offrano una via d'uscita, gli approcci esistenti presentano limitazioni significative:

• Variational Quantum Eigensolver (VQE): Un approccio ibrido quantistico-classico che spesso soffre di "barren plateaus" (gradienti che svaniscono) durante l'ottimizzazione. Trovare stati eccitati richiede tipicamente termini di penalità per sopprimere i contributi dello stato fondamentale, aumentando la complessità e il potenziale di errore.
• Quantum Annealing: Richiede tempi di evoluzione lunghi per garantire l'adiabaticità e fatica con piccoli gap energetici tra lo stato fondamentale e gli stati eccitati.
• Evoluzione nel Tempo Immaginario (ITE) & QITE: Sebbene promettenti, l'ITE standard affronta problemi di scalabilità nel realizzare operatori non unitari. Il Quantum ITE (QITE) spesso si basa su ottimizzazione variazionale o termini di penalità per gli stati eccitati, reintroducendo le sfide dei barren plateaus e richiedendo ansatz complessi.
• Altri Metodi: Tecniche come la Diagonalizzazione con Filtro Quantistico (QFD) e la Stima della Fase Quantistica Filtrata (FQPE) dipendono fortemente da stati di riferimento specifici o funzioni di filtro ottimizzate.

Gli autori identificano la necessità di un algoritmo interamente quantistico in grado di trovare sia stati fondamentali che eccitati senza funzioni di penalità, passaggi variazionali o loop ibridi classico-quantistici, evitando così i barren plateaus.

Metodologia

L'algoritmo proposto è un metodo interamente quantistico e senza penalità basato sul campionamento stocastico di una matrice di densità di stato misto combinato con una simulazione basata sugli stati.

1. Decomposizione dell'Hamiltoniano e Mappatura della Matrice di Densità

L'algoritmo assume che l'Hamiltoniano target $H$ possa essere decomposto in due componenti risolvibili, $H = h_1 + h_2$, dove gli autovalori ($\lambda^j_i$) e gli autostati ($|\psi^j_i\rangle$) di $h_1$ e $h_2$ sono noti e facili da preparare.

• Gli autori spostano e ridimensionano gli autovalori di $h_1$ e $h_2$ per creare una distribuzione di probabilità valida.
• Costruiscono una matrice di densità di stato misto $\rho$ equivalente all'Hamiltoniano spostato:
  $$\rho = \frac{1}{C} (H + (c_1 + c_2)I) = \sum_{i,j} p^j_i |\psi^j_i\rangle\langle\psi^j_i|$$
  dove $p^j_i$ sono probabilità derivate dagli autovalori spostati. Gli autostati di $\rho$ sono identici a quelli di $H$.

2. Evoluzione Stocastica nel Tempo Immaginario (ITE) per lo Stato Fondamentale

Per trovare lo stato fondamentale, l'algoritmo esegue ITE su $\rho$ utilizzando il campionamento stocastico:

• Campionamento: Ad ogni passo temporale $\eta$, un proiettore $P^j_i = |\psi^j_i\rangle\langle\psi^j_i|$ viene campionato dalla distribuzione $\{p^j_i\}$.
• Operatore di Evoluzione: Invece di implementare il pieno $e^{-\eta \rho}$ non unitario, l'algoritmo implementa l'operatore $I - \eta P^j_i$ per il proiettore campionato.
• Implementazione del Circuito: Ciò è ottenuto tramite un circuito unitario con post-selezione che coinvolge un qubit ancilla. L'unitaria $U$ è una porta controllata-SU(2) su $n$ qubit. Se il sistema è nello stato corrispondente al proiettore campionato, l'ancilla subisce una rotazione; altrimenti, rimane invariata. Post-selezionando l'ancilla nello stato $|0\rangle$ si ottiene l'operazione desiderata $I - \eta P^j_i$.
• Convergenza: Ripetendo questo processo, il sistema converge verso lo stato fondamentale di $\rho$ (e quindi di $H$) senza ottimizzazione variazionale.

3. Simulazione Basata sugli Stati per gli Stati Eccitati

Per trovare stati eccitati, l'algoritmo deve sopprimere i contributi degli stati a energia inferiore già trovati (ad esempio, lo stato fondamentale $P_g$).

• Sfida: Implementare direttamente proiettori per stati fondamentali altamente entangled è difficile.
• Soluzione: Gli autori utilizzano un protocollo di simulazione basata sugli stati. Questo coinvolge un sistema duale di $n$ qubit (ancilla) preparato nello stato fondamentale noto $P_g$, un qubit di controllo e il sistema principale.
• Meccanismo: Viene applicata una porta controllata-SWAP tra il sistema principale e l'ancilla, condizionata dal qubit di controllo. Misurare il qubit di controllo nella base $|+\rangle$ implementa efficacemente l'evoluzione $e^{-\eta P_g}$ (approssimata come $I - \eta P_g$), proiettando fuori la componente dello stato fondamentale.
• Procedura di Sollevamento: Per trovare il primo stato eccitato, la probabilità di campionare il proiettore dello stato fondamentale viene artificialmente "sollevata" (aumentata) nella distribuzione della matrice di densità. L'algoritmo esegue quindi ITE su questa distribuzione modificata. Quando viene campionato il proiettore dello stato fondamentale, viene utilizzata la simulazione basata sugli stati per filtrarlo. Questo processo può essere iterato per trovare stati eccitati superiori.

4. Variante per Hamiltoniani Generali

Per Hamiltoniani non facilmente divisibili in due parti risolvibili, gli autori propongono una variante in cui $H$ è trattata come una combinazione lineare di operatori di Pauli ($H = \sum \alpha_i O_i$). Il campionamento stocastico viene applicato direttamente a questi operatori di Pauli, con implementazioni circuitali simili per i passaggi di evoluzione.

Risultati Chiave

Gli autori hanno validato l'algoritmo utilizzando simulazioni classiche del Modello di Ising a Campo Trasverso (una catena 1D con condizioni al contorno periodiche), che funge da benchmark perché può essere diviso in componenti risolvibili nella base $Z$ e nella base $X$.

• Stato Fondamentale (sistema a 4 qubit): Ha raggiunto una fedeltà del 99,9% con lo stato fondamentale esatto.
• Stati Eccitati (sistema a 4 qubit):
  • 1° stato eccitato: fedeltà del 98,6%.
  • 2° stato eccitato: fedeltà del 96,8%.
  • 3° stato eccitato: fedeltà del 96,6%.
• Test di Scalabilità (sistema a 10 qubit): L'algoritmo ha trovato con successo lo stato fondamentale (fedeltà 98,6%) e il primo stato eccitato (fedeltà 97,8%). Ha anche risolto il terzo autostato più basso (fedeltà 90,2%) nonostante la quasi-degenerazione tra i primi due stati.
• Variante di Pauli: Anche la variante di campionamento degli operatori di Pauli ha dimostrato alte fedeltà (97,1%–99,5%) per i primi quattro autostati nel caso a 4 qubit.
• Scalabilità delle Risorse: L'analisi teorica indica che l'algoritmo scala con $O(n)$ porte a due qubit per passo, un miglioramento significativo rispetto ai metodi ITE basati su Trotter che scalano come $O(n^2)$.

Significato e Affermazioni

Il documento afferma di presentare un algoritmo interamente quantistico che affronta il problema degli autostati senza gli svantaggi dei metodi ibridi o variazionali attuali:

1. Senza Penalità: Trova stati eccitati senza aggiungere termini di penalità a una funzione di perdita, evitando la necessità di ansatz complessi e il rischio di barren plateaus.
2. Nessun Passo Ibrido: Elimina la necessità di loop di ottimizzazione classica, rendendolo una procedura puramente quantistica.
3. Ampia Applicabilità: L'assunzione che $H$ possa essere divisa in parti risolvibili è argomentata come non restrittiva, coprendo molti modelli fisici e persino problemi QMA-completi.
4. Robustezza: Il metodo dimostra la capacità di trovare sequenzialmente molteplici autostati a bassa energia, anche in sistemi con livelli energetici quasi degeneri.

Gli autori riconoscono che il tasso di successo complessivo dipende dalla post-selezione, che potrebbe essere una sfida sui dispositivi quantistici intermedi su scala rumorosa (NISQ) attuali a causa di errori di porta e decoerenza. Tuttavia, sostengono che il metodo offra un percorso concreto per utilizzare circuiti quantistici per affrontare problemi di autostati intrattabili classicamente, fungendo da aggiunta fondamentale alla cassetta degli attrezzi del calcolo quantistico.