Boson sampling enhanced quantum chemistry

Riepilogo Tecnico: Chimica Quantistica Potenziata dal Campionamento di Bosoni

Enunciato del Problema
Il Problema della Struttura Elettronica (ESP) consiste nel determinare l'energia dello stato fondamentale degli Hamiltoniani molecolari per prevedere le velocità delle reazioni chimiche. Sebbene il Variational Quantum Eigensolver (VQE) sia un approccio leader per risolvere l'ESP su dispositivi Quantum a Scala Intermedia Rumorosa (NISQ), affronta sfide significative sui Sistemi Ottici Quantistici Lineari (LQOS). Gli ansatz VQE convenzionali, come il Unitary Coupled Cluster (UCC), si basano sulla dinamica di fermioni interagenti e su circuiti quantistici profondi che sono difficili da implementare su LQOS passivi a causa della mancanza di interazioni naturali fotone-fotone. Inoltre, non esiste una mappatura diretta dallo spazio di Fock di Bosoni a più modi allo spazio di più qubit che preservi la struttura dei linguaggi standard di circuiti quantistici. La sfida consiste nello sviluppare un algoritmo VQE specificamente adattato per LQOS passivi, in grado di superare la trattabilità classica mantenendo al contempo l'efficienza hardware.

Metodologia
Gli autori propongono un algoritmo ibrido quantistico-classico denominato Boson Sampling-Classic (BS-C) VQE. Questo approccio utilizza esclusivamente sistemi ottici lineari passivi (interferometri lineari) combinati con metodi di chimica computazionale classica.

1. Costruzione dell'Ansatz:
   L'ansatz BS-C è definito come $|\Phi_{BS-C}\rangle = V_C(\vec{\beta}) U_B(\vec{\alpha}) |\Phi_{ref}\rangle$.

   • Parte Quantistica ($U_B(\vec{\alpha})$): Rappresenta la dinamica di Bosoni non interagenti realizzata da un Interferometro Ottico Lineare (LOI) parametrizzato. Lo stato iniziale $|\Phi_{ref}\rangle$ è preparato utilizzando sorgenti a singolo fotone per gli orbitali occupati e il vuoto per gli orbitali virtuali (codifica single-rail).
   • Parte Classica ($V_C(\vec{\beta})$): Rappresenta operazioni di chimica computazionale classica (ad esempio, Hartree-Fock o Interazione di Configurazione) che agiscono sull'Hamiltoniano. Nello specifico, l'Hamiltoniano viene trasformato classicamente come $H_T(\vec{\beta}) = V_C^\dagger(\vec{\beta}) H_F V_C(\vec{\beta})$ prima di essere mappato su un Hamiltoniano di qubit tramite la trasformazione di Jordan-Wigner (JW).
   • Distinzione delle Risorse: A differenza dei metodi classici (HF, CISD) o dell'UCC che si basano su eccitazioni elettroniche (determinanti), la parte Bosonica introduce risorse basate sui permanenti delle matrici. Poiché il calcolo dei permanenti è #P-difficile, ciò fornisce una risorsa intrattabile classicamente distinta dalle eccitazioni fermioniche standard.

2. Funzione di Costo e Misurazione Ibrida:
   Il valore di aspettazione dell'energia non è una misurazione diretta a causa della natura non unitaria di alcuni operatori classici e dell'assenza del principio di esclusione di Pauli nei fotoni (che consente più fotoni in un singolo modo). L'energia vera è definita come un rapporto:
   $$E(\vec{\alpha}, \vec{\beta}) = \frac{\langle \Phi(\vec{\alpha}) | H_T^{JW}(\vec{\beta}) | \Phi(\vec{\alpha}) \rangle}{\langle \Phi(\vec{\alpha}) | V_C^\dagger(\vec{\beta}) V_C(\vec{\beta}) | \Phi(\vec{\alpha}) \rangle}$$
   Per misurare questo in modo efficiente, gli autori propongono una strategia di misurazione ibrida:

   • Rilevamento del Numero di Fotoni: Utilizzato per i termini che coinvolgono l'Identità ($I$) e gli operatori di Pauli-$Z$.
   • Rilevamento Omodino: Utilizzato per i termini che coinvolgono gli operatori di Pauli-$\sigma^+$ e $\sigma^-$.
     Questa combinazione permette la valutazione della funzione di costo senza richiedere la tomografia completa dello stato quantistico, sfruttando il fatto che le funzioni di kernel per le misurazioni omodine possono essere calcolate in modo efficiente.

3. Mitigazione degli Errori:
   Il documento affronta la perdita di fotoni, un canale di rumore dominante negli LQOS. Poiché la perdita di fotoni fa uscire il sistema dal sottospazio legale (dove il numero di fotoni è uguale al numero di elettroni), gli autori introducono una procedura di correzione. Eseguendo misurazioni aggiuntive del numero puro di fotoni per stimare la probabilità di specifici eventi di perdita di fotoni, derivano un fattore di normalizzazione per correggere i valori di aspettazione grezzi. Questo metodo si dimostra scalabile e non richiede una post-selezione che aumenterebbe esponenzialmente i costi di campionamento.

Contributi Chiave

• Ansatz BS-C: Un nuovo ansatz VQE che ibrida la dinamica di Bosoni non interagenti con metodi di chimica classica (HF, CISD), sfruttando la difficoltà computazionale dei permanenti per migliorare l'accuratezza.
• Protocollo di Misurazione Scalabile: Uno schema ibrido di rilevamento omodino e del numero di fotoni che valuta in modo efficiente la funzione di costo per Hamiltoniani mappati tramite JW, superando l'incompatibilità delle misurazioni omodine standard con i termini di Pauli non locali.
• Mitigazione della Perdita di Fotoni: Una tecnica di mitigazione degli errori scalabile che corregge la perdita di fotoni senza sovraccarico esponenziale, basandosi su fattori di normalizzazione statistici derivati da misurazioni ripetute.
• Efficienza Hardware: La proposta dimostra che gli LQOS passivi possono risolvere l'ESP con profondità di circuito significativamente minori ($O(M)$) rispetto all'UCCSD ($O(M^3)$) e offre tassi di campionamento molto più elevati ($10^6-10^8$ Hz) rispetto alle piattaforme a trappola ionica o superconduttive.

Risultati
Gli autori hanno eseguito esperimenti numerici su diverse molecole: LiH, BeH$_2$ e H$_4$ (in geometrie lineare e quadrata) utilizzando la base STO-3G.

• Accuratezza: Il metodo BS-C (nello specifico BS-HF per LiH e BS-CISD per gli altri) ha raggiunto l'accuratezza chimica ($1.6 \times 10^{-3}$ Hartree) su varie lunghezze di legame, superando i corrispettivi classici (HF e CISD) presi singolarmente.
• Rapporto di Proiezione: Lo studio ha analizzato il "rapporto di proiezione" (la sovrapposizione tra lo stato di Bosone e il sottospazio di Fermione codificato). I risultati hanno indicato che, sebbene il rapporto vari, può essere mantenuto a un livello sufficiente per la scalabilità scegliendo basi appropriate (aumentando il numero di orbitali rispetto agli elettroni) o limitando le eccitazioni per preservare le simmetrie.
• Prestazioni: Il metodo ha generato con successo curve di energia potenziale che corrispondevano strettamente ai benchmark Full Configuration Interaction (FCI), dimostrando l'efficacia delle risorse basate sui permanenti nel catturare la correlazione elettronica.

Significato e Affermazioni
Il documento afferma che il BS-C VQE offre una via praticabile per risolvere il Problema della Struttura Elettronica utilizzando sistemi ottici lineari passivi nell'era NISQ. Il significato risiede in:

1. Sfruttamento degli LQOS Passivi: Dimostra che gli LQOS passivi, nonostante la mancanza di interazioni naturali, possono risolvere problemi pratici sfruttando la complessità computazionale unica del campionamento di Bosoni (permanenti) per potenziare i metodi classici.
2. Vantaggio Hardware: L'approccio è altamente efficiente dal punto di vista hardware, richiedendo solo interferometri ottici lineari e offrendo tassi di campionamento di ordini di grandezza superiori rispetto ad altre piattaforme quantistiche, il che è critico per gli algoritmi variazionali.
3. Scalabilità: Le strategie ibride di misurazione e mitigazione degli errori proposte affrontano i principali colli di bottiglia del rumore e del sovraccarico di misurazione, suggerendo che il metodo è scalabile per molecole più grandi, a condizione che il rapporto di proiezione sia gestito.

Gli autori concludono che, sebbene la proposta sia alla portata della tecnologia attuale (utilizzando componenti esistenti come sorgenti a singolo fotone e interferometri riconfigurabili), sono necessarie ulteriori indagini riguardo alle prestazioni su molecole più grandi e allo sviluppo di strategie di misurazione più efficienti.