Orbital-optimized spin-adapted multistate contracted VQE for excited states and properties on quantum hardware

Questo articolo introduce il metodo orbital-optimized multistate contracted VQE (oo-MC-VQE), che utilizza operatori adattati allo spin per calcolare efficientemente gli stati fondamentale e eccitati insieme alle loro proprietà su hardware quantistico, bilanciando accuratezza e complessità del circuito attraverso una scalabilità lineare dei parametri con il numero di stati.

L'essenziale

Immagina di cercare di accordare una massiccia e complessa orchestra per suonare una sinfonia perfetta. Nel mondo della chimica, l'"orchestra" è una molecola, e la "musica" è il modo in cui i suoi elettroni si muovono e interagiscono. Per capire come una molecola assorbe la luce (il che ci fornisce i colori e guida processi come la fotosintesi), gli scienziati devono calcolare le note esatte che questi elettroni suonano.

Per molto tempo, calcolare questo per molecole con molti elettroni è stato come cercare di risolvere un puzzle che diventa esponenzialmente più difficile man mano che si aggiungono pezzi. I computer classici (quelli che usiamo oggi) alla fine incontrano un muro e non riescono a risolvere questi puzzle per molecole complesse.

Questo articolo presenta un nuovo modo per risolvere questi puzzle utilizzando i computer quantistici, che sono macchine speciali progettate per gestire naturalmente questo tipo di complessità. Ecco una semplice suddivisione di ciò che gli autori hanno fatto e scoperto:

1. Il Problema: Accordare Molte Note Contemporaneamente

Di solito, gli scienziati cercano di accordare l'orchestra per suonare perfettamente anche solo una nota (lo stato fondamentale). Ma per capire come una molecola reagisce alla luce, abbiamo bisogno di conoscere molte note diverse (stati eccitati) contemporaneamente.

• La Sfida: Se provi ad accordare l'orchestra per 10 canzoni diverse simultaneamente, le istruzioni (il circuito del computer) diventano incredibilmente lunghe e complicate. Se le istruzioni sono troppo lunghe, il computer quantistico si confonde a causa del "rumore" (disturbi o errori), e la musica va in pezzi.
• Il Compromesso: Hai bisogno di un circuito complesso per ottenere una risposta accurata, ma un circuito complesso ha maggiori probabilità di fallire su macchine rumorose come quelle attuali.

2. La Soluzione: Un Direttore d'Orchestra Intelligente e Simmetrico

Gli autori hanno sviluppato un nuovo metodo chiamato oo-MC-VQE. Immagina questo come un "direttore d'orchestra intelligente" per l'orchestra quantistica.

• Adattato allo Spin (Spin-Adapted): Nella chimica quantistica, gli elettroni hanno una proprietà chiamata "spin" (come trottole che ruotano). Gli autori hanno costruito il loro metodo in modo che il direttore mantenga sempre le trottole in rotazione nel modo corretto e simmetrico. Questo evita che la musica si "scordi" a causa di errori di simmetria.
• Ottimizzato per gli Orbitali (Orbital Optimized): Hanno anche permesso al direttore di riorganizzare lo schema dei posti a sedere dei musicisti (gli orbitali) per far suonare meglio la musica prima ancora di iniziare il complesso lavoro di accordatura.
• Contratto Multistato (Multistate Contracted): Invece di cercare di accordare 10 canzoni con 10 manuali di istruzioni separati e massicci, hanno trovato un modo per utilizzare un unico set di istruzioni condivuso ed efficiente che funzioni per tutte le canzoni contemporaneamente.

3. La Scoperta: Crescita Lineare

Una delle grandi domande era: Se voglio calcolare 10 stati invece di 1, ho bisogno di 10 volte più potenza di calcolo?

• Il Risultato: Gli autori hanno scoperto che la risposta è sorprendentemente semplice. La quantità di "sforzo" del computer (parametri del circuito) necessaria cresce in modo lineare. Se raddoppi il numero di stati che vuoi calcolare, aumenti approssimativamente solo il doppio della lunghezza del manuale di istruzioni. È un'ottima notizia perché significa che il metodo è scalabile.

4. Il Test nel Mondo Reale: Suonare su un Palco Rumoroso

Gli autori non hanno solo simulato questo su un computer perfetto; hanno effettivamente eseguito il loro metodo su hardware quantistico reale (computer quantistici IBM).

• La Configurazione: Hanno testato due piccole molecole: la Formaldeide (un prodotto chimico comune) e il catione triidrogeno ($H_3^+$).
• Il Problema del Rumore: I veri computer quantistici sono come un palco con una folla rumorosa e luci che sfarfallano. I risultati erano disordinati senza aiuto.
• La Soluzione: Hanno utilizzato tecniche di "mitigazione dell'errore". Immagina questo come un ingegnere del suono che usa un software per filtrare il rumore della folla e lo sfarfallio delle luci dopo l'esecuzione.
• Il Risultato:
  • Per la Formaldeide, il metodo ha funzionato abbastanza bene. Anche con il rumore, potevano vedere chiaramente i "picchi" nello spettro di assorbimento (i colori che la molecola assorbe).
  • Per l' $H_3^+$, il rumore è stato un problema maggiore, spostando significativamente i risultati. Gli autori hanno notato che questo accade perché la matematica di questa specifica molecola è più sensibile al rumore (come uno strumento delicato che si scorda facilmente).
  • Conclusione Chiave: Sebbene i numeri non fossero perfetti sulle macchine reali, la forma dei risultati era corretta. Potevano comunque vedere le caratteristiche principali del comportamento della molecola.

Riassunto

L'articolo dimostra che, utilizzando un approccio intelligente e simmetrico, gli scienziati possono calcolare il comportamento degli elettroni eccitati nelle molecole utilizzando gli attuali computer quantistici, che sono imperfetti. Hanno dimosto che calcolare più stati non richiede una quantità impossibile di risorse e che, con alcuni trucchi di "cancellazione del rumore", è possibile ottenere informazioni chimiche utili dai dispositivi quantistici reali già oggi.

Cosa NON hanno affermato:
L'articolo non sostiene che questo metodo possa immediatamente progettare nuovi pannelli solari, curare malattie o creare nuovi materiali. Si concentra strettamente sul dimostrare che il metodo funziona per calcolare gli spettri su hardware quantistico. Qualsiasi applicazione futura è implicata dal campo generale, ma non è un'affermazione specifica di questo studio.

Sommario tecnico

Sintesi Tecnica: VQE multistato contratto con ottimizzazione orbitale e adattamento di spin

Problematica
La computazione robusta e bilanciata degli stati elettronici eccitati e delle loro proprietà di transizione rimane una sfida significativa, particolarmente in regimi che comportano una forte riorganizzazione elettronica, degenerazioni prossime o carattere multiconfigurazionale. Sebbene i metodi a risposta lineare (LR) siano efficaci in molti contesti, essi trattano gli stati eccitati come perturbazioni di un riferimento dello stato fondamentale, portando a limitazioni concettuali in prossimità delle intersezioni coniche e in sistemi fortemente correlati. Inoltre, gli attuali approcci di Variational Quantum Eigensolver (VQE) per gli stati eccitati, come il subspace-search VQE (SS-VQE) e il multistate contracted VQE (MC-VQE), affrontano spesso un compromesso tra l'accuratezza nella descrizione di un numero crescente di stati elettronici e l'espressività (complessità del circuito) dell'ansatz sottostante. Inoltre, l'attuale hardware quantistico è limitato dal rumore e dalla decoerenza, il che restringe la profondità del circuito e rende necessari lo sviluppo di metodi che siano sia efficienti che compatibili con le tecniche di mitigazione dell'errore.

Metodologia
Gli autori introducono il metodo orbital-optimized multistate contracted variational quantum eigensolver (oo-MC-VQE), che incorpora operatori adattati allo spin per calcolare gli stati fondamentale ed eccitati, nonché le proprietà specifiche dello stato e di transizione.

• Adattamento di Spin: Per garantire che la simmetria di spin degli stati di riferimento sia conservata durante tutta l'ottimizzazione VQE, il metodo utilizza un pool di operatori adattati allo spin. Nello specifico, gli autori impiegano operatori di singola eccitazione e operatori di doppia eccitazione di coppia (pair-double excitation) adattati allo spin singoletto. Ciò evita la necessità di termini di penalità di spin e garantisce la corretta simmetria di spin di tutti gli stati inclusi senza introdurre parametri ridondanti.
• Framework Multistato Contrattato: Il metodo ottimizza simultaneamente più stati elettronici utilizzando un funzionale di energia mediata sugli stati (state-averaged). Ogni stato è costruito utilizzando un circuito di preparazione dello stato fisso e privo di parametri (inizializzando basi ortogonali) seguito da un circuito correlatore parametrizzato (utilizzando un ansatz Unitary Product State, o UPS) che introduce la correlazione elettronica. L'ottimizzazione mira all'energia mediata sugli stati, seguita da una diagonalizzazione classica della matrice hamiltoniana multistato per risolvere le energie e le funzioni d'onda dei singoli stati.
• Ottimizzazione Orbitale: Il framework viene esteso per includere l'ottimizzazione orbitale (oo-MC-VQE). Le rotazioni orbitali sono trattate come una funzione di costo classica, permettendo l'ottimizzazione simultanea dei parametri del circuito ($\theta$) e dei parametri di rotazione orbitale ($\kappa$). Gli autori dimostrano che, per l'ansatz specifico utilizzato, i parametri di circuito delle singole eccitazioni adiacenti alla parametrizzazione orbitale sono ridondanti e possono essere rimossi, riducendo la complessità del circuito.
• Valutazione degli Elementi di Matrice: Per valutare gli elementi di matrice tra diversi stati (necessari per la matrice hamiltoniana multistato e le proprietà di transizione), gli autori impiegano una tecnica di decomposizione di superposizione. Questo approccio evita circuiti di preparazione dello stato profondi per gli stati multi-determinanti, valutando superposizioni più semplici di due determinanti, riducendo così il numero di gate CNOT a costo di un aumento delle misurazioni del valore di aspettazione.
• Mitigazione dell'Errore: L'implementazione integra Mansatz0, una tecnica di mitigazione dell'errore di lettura e di gate basata sull'ansatz; per i calcoli multistato che coinvolgono superposizioni di stati di riferimento, gli autori propongono un'estensione (M0+) che caratterizza gli errori per tutti gli unici stati di riferimento multi-determinanti, pur notando la scalabilità esponenziale di questo approccio.

Contributi Chiave

1. oo-MC-VQE adattato allo spin: Lo sviluppo di una formulazione adattata allo spin per il VQE a media sugli stati con ottimizzazione orbitale, garantendo la conservazione della corretta simmetria di spin attraverso più stati.
2. Analisi di Scalabilità: Un'indagine esplicita sul compromesso tra accuratezza e complessità del circuito. Gli autori riscontrano che il numero di parametri del circuito richiesti per ottenere risultati accurati scala approssimativamente in modo lineare con il numero di stati inclusi nella media sugli stati.
3. Implementazione su Hardware Quantistico: Un'esplicita implementazione in circuito quantistico del metodo, che include la riduzione dei gate CNOT tramite superposizione di determinanti e l'integrazione di tecniche di mitigazione dell'errore.
4. Dimostrazione su Dispositivi Reali: Esecuzione del metodo su hardware quantistico reale (IBM ibm marrakesh e ibm aachen) per calcolare gli spettri di assorbimento per sistemi molecolari benchmark (formaldeide e $H_3^+$).

Risultati

• Complessità del Circuito: Le simulazioni su LiH, $H_2O$ e $NH_3$ utilizzando gli ansatz ADAPT e oo-ADAPT mostrano che il numero di parametri richiesti aumenta approssimativamente in modo lineare con il numero di stati. Ad esempio, ottenere un errore inferiore a $10^{-4}$ Hartree per $H_2O$ ha richiesto 27, 32, 24 e 30 parametri per stato per 1, 2, 3 e 4 stati, rispettivamente.
• Impatto dell'Ottimizzazione Orbitale: Sebbene l'ottimizzazione orbitale migliori la convergenza nelle prime fasi della procedura ADAPT, ha un effetto trascurabile sul numero finale di parametri del circuito per una convergenza più stretta. In alcuni casi (ad esempio, $NH_3$ con 3 o 4 stati), l'inclusione delle singole eccitazioni nell'ottimizzazione orbitale ha avuto un effetto deleterio sul conteggio dei parametri.
• Performance dell'Hardware Quantistico:
  • Rumore Simulato: Su hardware emulato utilizzando modelli di rumore IBM, le tecniche di mitigazione dell'errore (M0 e M0+) hanno prodotto risultati quasi identici alle simulazioni senza rumore. Il sistema $H_3^+$ è stato più avversamente influenzato dal rumore rispetto alla formaldeide, a causa di un numero di condizionamento più elevato della matrice hamiltoniana multistato.
  • Hardware Reale: I calcoli su dispositivi reali hanno mostrato che, sebbene i risultati mitigati dall'errore preservassero i trend qualitativi e le caratteristiche primarie degli spettri di assorbimento, l'accuratezza quantitativa era significativamente inferiore rispetto ai dispositivi simulati. Questa discrepanza evidenzia i limiti degli attuali modelli di rumore riguardo al rumore tempo-dipendente e al crosstalk. Tuttavia, il metodo è riuscito a estrarre intuizioni chimiche significative nonostante il rumore dell'hardware.

Significatività e Rivendicazioni
Il lavoro sostiene che il metodo oo-MC-VQE fornisca una descrizione elettronica bilanciata, adatta a scenari impegnativi come le intersezioni coniche, affrontando le limitazioni intrinseche degli approcci a risposta lineare. Gli autori sottolineano che la scalabilità lineare dei parametri del circuito con il numero di stati suggerisce che il metodo sia praticabile per mirare a molteplici stati eccitati su hardware quantistico near-term. Pur riconoscendo che l'accuratezza quantitativa rimane una sfida sugli attuali dispositivi, il lavoro dimostra che è possibile estrarre intuizioni chimiche significative, come gli spettri di assorbimento, combinando il metodo con appropriate tecniche di mitigazione dell'errore. Lo studio funge da prova di concetto per l'esecuzione di calcoli multistato, adattati allo spin e con ottimizzazione orbitale su veri dispositivi quantistici.