Symmetry-Protected Basin Localization in Variational Quantum Eigensolvers

Riepilogo Tecnico: Localizzazione in Bacini Protetti da Simmetria negli Eigensolver Quantistici Variazionali

Enunciato del Problema
Gli Eigensolver Quantistici Variazionali (VQE) affrontano un ostacolo fondamentale nei regimi molecolari fortemente correlati: il paesaggio energetico non convesso è partizionato dalla simmetria molecolare in bacini di attrazione concorrenti. La modalità di fallimento primaria non si verifica durante il ciclo di ottimizzazione, ma all'inizializzazione.

1. Inizializzazione Casuale: Approssima un insieme unitario simile a quello di Haar, portando a "altipiani sterili" (barren plateaus) dove la varianza del gradiente svanisce esponenzialmente ($O(1/2^N)$) a causa della concentrazione della misura, rendendo la discesa locale non informativa.
2. Inizializzazione a Campo Medio (Hartree–Fock): Fallisce nei regimi fortemente correlati (oltre il punto di Coulson–Fischer) perché rompe spontaneamente la simmetria di spin. Ciò sposta lo stato variazionale verso varietà a simmetria rotta con sovrapposizione trascurabile con il settore dello stato fondamentale correlato, deviando l'ottimizzazione verso un minimo locale concorrente.

In entrambi i casi, la misura di inizializzazione è spostata dal bacino dello stato fondamentale prima che venga compiuto il primo passo del gradiente, causando il fallimento dell'algoritmo prima ancora che inizi l'ottimizzazione.

Metodologia
Gli autori introducono un precondizionatore condizionato alla geometria, indicato come $P_{eq}: \mathbb{R} \mapsto \theta_0$, che mappa la geometria nucleare direttamente sui parametri del circuito. Questa mappa è vincolata dalla covarianza $SE(3)$ dell'Hamiltoniano molecolare.

• Quadro Teorico: Il precondizionatore è definito come un vincolo olonomo che commuta con il gruppo di simmetria fisica dell'Hamiltoniano. Richiedendo che $| \psi(P_{eq}(g \cdot R)) \rangle \simeq \hat{U}(g) | \psi(P_{eq}(R)) \rangle$ per ogni $g \in SE(3)$, la mappa rimuove i gradi di libertà ridondanti di rotazione e traslazione. Ciò restringe la misura di inizializzazione indotta allo spazio quoziente $\mathbb{R}^{3M}/SE(3)$ e forza lo stato iniziale $\theta_0$ nella componente consistente con la simmetria del bacino dello stato fondamentale.
• Implementazione: $P_{eq}$ è realizzato utilizzando una rete di passaggio messaggi equivariante sotto $SE(3)$ (tipo MACE), costruita a partire dai principi di accoppiamento tensoriale utilizzati nei modelli atomistici. L'architettura della rete rispetta i vincoli di località del circuito quantistico (adattando il campo ricettivo al cono causale del circuito).
• Addestramento: La rete è addestrata su benchmark di vettori di stato (ottenuti tramite basin-hopping e raffinamento L-BFGS-B) per prevedere rappresentanti parametrici attrattivi per i bacini. Viene utilizzata una funzione obiettivo consapevole del gauge per gestire le ridondanze parametriche, combinando l'ancoraggio dei parametri con termini di fedeltà dello stato.
• Protocollo Ibrido: Per sistemi disordinati, gli autori propongono un metodo ibrido che combina il targeting deterministico del bacino equivariante con una fuga stocastica per risolvere trappole su scala fine all'interno del bacino localizzato.

Contributi Chiave

1. Passaggio dal Controllo della Concentrazione al Controllo della Curvatura: Il documento stabilisce un cambiamento teorico nelle statistiche dei gradienti. I VQE standard soffrono di gradienti controllati dalla concentrazione (varianza soppressa dal volume dello spazio di Hilbert). La localizzazione in bacini tramite $P_{eq}$ cambia il regime in gradienti controllati dalla curvatura, dove la varianza è determinata dalla curvatura locale dell'Hessiano ($\kappa_j^2 \sigma_j$) all'interno di un intorno fortemente convesso, indipendentemente dalla dimensione completa dello spazio di Hilbert.
2. Selezione del Bacino Pre-Shot: La procedura esegue l'identificazione del bacino prima del ciclo quantistico limitato dal numero di shot. La mappa classica seleziona il bacino compatibile con la simmetria, e il circuito quantistico è riservato esclusivamente al raffinamento della correlazione all'interno di quel bacino selezionato.
3. Derivazione Formale: Gli autori derivano la Proposizione 1 (e il Teorema S1 nel Materiale Supplementare), dimostrando formalmente che se la misura di inizializzazione è supportata all'interno di un bacino quadratico fortemente convesso, la varianza del gradiente scala con la curvatura locale piuttosto che esponenzialmente con la dimensione del sistema.

Risultati
Il metodo è stato valutato su sei benchmark molecolari allungati ($N_2$, $H_2O$, $BeH_2$, $CO$, $LiH$, $H_8$) e catene disordinate di $H_{10}$ utilizzando simulazioni di vettori di stato:

• Riduzione dell'Errore di Inizializzazione: $P_{eq}$ ha ridotto gli errori di inizializzazione Hartree–Fock di fattori compresi tra 38$\times$ e 6250$\times$.
• Accuratezza:
  • Ha raggiunto un'accuratezza di inizializzazione sotto il millihartree (sub-mHa) per $CO$, $LiH$e$H_8$.
  • Ha collocato $N_2$, $H_2O$ e $BeH_2$ all'interno del bacino correlato su scala mHa.
  • In $N_2$ allungato ($R=2.5$ Å), l'inizializzazione Hartree–Fock è caduta in un bacino a simmetria rotta con $\Delta E > 0.2$ Ha, mentre $P_{eq}$ ha collocato lo stato nel bacino dello stato fondamentale correlato.
• Robustezza al Disordine: Nelle catene disordinate di $H_{10}$, il protocollo ibrido (targeting equivariante + fuga stocastica) ha raggiunto una probabilità di successo unitaria con un budget di ottimizzazione fisso, mentre l'inizializzazione casuale ha sofferto di distribuzioni di errore a code pesanti e frequenti errori di identificazione del bacino.
• Trasferibilità: Il metodo ha dimostrato trasferibilità verso sistemi poliatomici (ad esempio, stiramento O–H in $H_2O$), mantenendo un'accuratezza su scala mHa attraverso geometrie in cui l'inizializzazione casuale falliva di ordini di grandezza.

Significato e Affermazioni
Il documento afferma che la barriera primaria per i VQE nei sistemi fortemente correlati è la selezione del bacino, non l'inefficienza dell'ottimizzatore. Imponendo l'equivarianza $SE(3)$ sulla mappa di inizializzazione, gli autori dimostrano che:

• La scoperta del bacino delocalizzato può essere sostituita dalla discesa locale controllata dalla curvatura.
• Il fenomeno dell'"altipiano sterile" può essere visto come una proprietà delle misure di inizializzazione delocalizzate; restringere la misura a un supporto localizzato nel bacino e consistente con la simmetria aggira la soppressione esponenziale dei gradienti.
• L'approccio disaccoppia il costo della ricerca geometrica (gestita classicamente e a basso costo) dal costo del raffinamento della correlazione quantistica. Ciò sposta il carico computazionale dal ciclo quantistico limitato dal numero di shot a una mappa classica protetta dalla simmetria, offrendo un aumento di velocità in termini di tempo reale (wall-clock speedup) di $10^4$–$10^6$ per la fase di inizializzazione.

Gli autori mantengono un tono modesto riguardo alla portata, notando che la costruzione è specifica per Hamiltoniani di struttura elettronica con simmetria $SE(3)$ e ansatz che preservano la località, e che l'estensione a sistemi periodici o Hamiltoniani con accoppiamento spin-orbita richiede l'adattamento delle mappe equivarianti ai gruppi di simmetria pertinenti.