Electron-molecule scattering via R-matrix variational algorithms on a quantum computer

Sintesi Tecnica: Scattering Elettrone-Molecola tramite Algoritmi Variazionali R-matrix su un Computer Quantistico

Enunciato del Problema
Gli urti tra elettroni e molecole sono fondamentali per i processi naturali e le tecnologie come la lavorazione al plasma. Sebbene il metodo R-matrix costituisca una strategia computazionale standard per questi problemi, affronta sfide significative di scalabilità quando applicato a sistemi complessi. Il problema della regione interna, che comporta la risoluzione dell'Hamiltoniana a molti elettroni all'interno di un volume confinato per catturare la correlazione e lo scambio elettronico, scala in modo sfavorevole con il numero di elettroni e la dimensione del set di basi. Gli approcci convenzionali spesso fanno affidamento su metodi approssimati o diventano computazionalmente proibitivi per soluzioni ad alta accuratezza. Inoltre, gli algoritmi variazionali standard di chimica quantistica mirano tipicamente allo stato fondamentale, mentre lo scattering R-matrix richiede la determinazione simultanea di molteplici autostati (inclusi pseudostati eccitati e stati di tipo continuo) all'interno di un settore di simmetria specifico per definire le condizioni al contorno dello scattering.

Metodologia
Gli autori propongono un quadro computazionale quantistico che mappa il problema della regione interna R-matrix su un Hamiltoniana di qubit utilizzando il Variational Quantum Eigensolver (VQE) e le sue varianti. L'approccio coinvolge diversi componenti chiave:

1. Mappatura dell'Hamiltoniana e Vincoli: L'Hamiltoniana per $(N+1)$ elettroni per la regione interna è mappata sui qubit tramite la trasformazione di Jordan–Wigner. Viene imposto un vincolo fisico critico: al massimo un elettrone può occupare qualsiasi orbitale continuo. Ciò è realizzato strutturalmente all'interno dell'ansatz variazionale utilizzando una formulazione di operatore di proiezione del numero, garantendo che lo stato preparato $|\psi(\theta)\rangle$ soddisfi $P|\psi(\theta)\rangle = |\psi(\theta)\rangle$, dove $P$ proietta sulle configurazioni di occupazione valide.
2. Architettura del Circuito: L'algoritmo utilizza un sistema a tre registri: un registro target ($t$), un registro continuo ($c$) e un registro selettore ($a$). Il circuito prepara stati di prova che coprono una base troncata di funzioni di canale antisimmetrizzate e configurazioni legate ($L^2$).
   • Simmetria di Spin: L'ansatz impone la simmetria totale dello spin ($S, M$) utilizzando rotture di Clebsch–Gordan (CG) fisse (implementate come rotazioni Givens a due qubit) per accoppiare lo spin del target con lo spin dell'elettrone in arrivo.
   • Ottimizzazione Sussessiva del Sottospazio (SSO): Per recuperare lo spettro completo degli autostati senza l'overhead della diagonalizzazione del sottospazio o dei termini di penalità, gli autori impiegano un protocollo SSO. Ciò comporta la decomposizione di una rotazione $SO(k)$ in una cascata di rotazioni Givens a due qubit. L'ottimizzazione procede in sequenza: il più basso autostato viene ottimizzato per primo e i suoi parametri vengono fissati. I round successivi ottimizzano il successivo autostato all'interno del complemento ortogonale, diagonalizzando efficacemente a blocchi l'Hamiltoniana uno stato alla volta. Ciò permette l'uso di un semplice valore di aspettazione dell'Hamiltoniana $\langle H \rangle$ come funzione di costo, evitando la necessità di termini dell'Hamiltoniana al quadrato ($H^2$) o misurazioni di sovrapposizione.
3. Lettura e Ampiezze al Confine: Viene introdotto un protocollo di somma coerente per stimare i valori di aspettazione sul registro di sistema selezionando postivamente il registro selettore, riducendo il numero di valori di aspettazione distinti richiesti. Crucialmente, i parametri ottimali del circuito codificano direttamente i coefficienti di espansione ($a_{ijk}$) della funzione d'onda. Questi coefficienti sono utilizzati per calcolare le ampiezze al confine R-matrix ($w_{ik}(a)$), che servono come condizioni di matching per il calcolo dello scattering nella regione esterna.

Contributi Chiave

• Prima Formulazione Quantistica della Regione Interna R-matrix: A conoscenza degli autori, questa è la prima applicazione di algoritmi quantistici allo scattering elettrone-molecola e la prima formulazione del problema della regione interna R-matrix su un computer quantistico.
• Circuito che Impone la Simmetria: Il lavoro introduce il primo utilizzo di un circuito che impone la simmetria di Clebsch–Gordan all'interno di un ansatz variazionale per lo scattering, utilizzando blocchi Givens fissi per mantenere la simmetria di spin.
• Ottimizzazione Efficiente Multi-Stato: Il protocollo SSO proposto differisce dai metodi esistenti come SSVQE e MCVQE imponendo l'ortogonalità attraverso la struttura del circuito piuttosto che tramite termini di penalità o misurazioni di sovrapposizione. Ciò elimina la necessità di misurare elementi fuori diagonale dell'Hamiltoniana o di eseguire la diagonalizzazione classica del sottospazio, riducendo la funzione di costo a semplici valutazioni di $\langle H \rangle$.
• Estrazione Diretta dei Dati di Scattering: Il metodo dimostra che i parametri variazionali ottimali producono direttamente le ampiezze al confine R-matrix, bypassando la necessità di un'ulteriore post-elaborazione per estrarre osservabili di scattering.

Risultati
L'algoritmo è stato dimostrato su un problema modello di scattering di elettroni dalla molecola di idrogeno ($H_2$) utilizzando un simulatore classico privo di rumore.

• Configurazione del Sistema: La simulazione ha utilizzato una base minima (STO-3G) con $N=2$ elettroni target, risultando in uno spazio di prova di dimensione $k=5$ (quattro rami di canale aperto e una configurazione legata) all'interno del settore di simmetria $S=1/2, M=-1/2, B_{1u}$.
• Prestazioni: L'approccio SSO ha recuperato con successo lo spettro completo di cinque autostati con errori energetici inferiori a $10^{-7}$ Hartree.
• Efficienza delle Risorse: Il metodo SSO ha richiesto solo 573 valutazioni della funzione di costo coinvolgenti 92 termini di Pauli. Al contrario, un'ottimizzazione simultanea che minimizza la somma delle varianze ha richiesto 13.345 valutazioni di 1.080 termini di Pauli, e un approccio a singolo stato con Hamiltoniana ripiegata ha richiesto 1.397 valutazioni.
• Profondità del Circuito: Il circuito non ottimizzato per questo esempio minimale ha utilizzato 7 qubit, 217 porte CNOT e una profondità di 314.

Significato e Affermazioni
Il documento afferma di fornire una via concreta per applicare l'hardware quantistico allo scattering elettrone-molecola, affrontando specificamente il collo di bottiglia della regione interna del metodo R-matrix. Gli autori sottolineano che il loro approccio:

1. Evita l'archiviazione di grandi matrici Hamiltoniane rappresentando gli stati di prova sui qubit.
2. Riduce l'overhead di misurazione evitando osservabili dell'Hamiltoniana al quadrato e utilizzando l'ottimizzazione sequenziale.
3. Fornisce un output "nativo del problema" dove i parametri del circuito mappano direttamente su quantità fisiche di scattering (ampiezze al confine).

Gli autori notano con modestia che, sebbene non affermino un vantaggio quantistico asintotico provato rispetto ai migliori metodi iterativi classici al momento, il metodo offre un paradigma ibrido quantistico-classico vitale. Sottolineano che l'approccio è particolarmente promettente all'aumentare della dimensione dello spazio di prova, dove la diagonalizzazione classica scala come $O(k^3)$, mentre l'approccio quantistico scala il conteggio dei parametri e la profondità come $O(k^2)$. Il lavoro getta le basi per future estensioni a target più grandi e rappresentazioni continue più sofisticate, integrando potenzialmente processori quantistici in codici di scattering consolidati come UKRmol+.