Tailored coupled cluster method with sample-based quantum diagonalization: Application to titanium-based metallocene catalytic reactions for 1-hexene production

Sintesi Tecnica: Metodo Coupled Cluster su Misura con Diagonalizzazione Quantistica Basata su Campionamento

Definizione del Problema
La previsione accurata della selettività nelle reazioni catalitiche dei metalloceni a base di titanio per la produzione di 1-esene richiede il calcolo delle differenze di energia libera tra stati di transizione competitivi (inserzione dell'etilene vs. trasferimento di beta-idruro) con un'accuratezza entro 1 kcal/mol. Gli attuali metodi di struttura elettronica affrontano una dicotomia: la Teoria del Funzionale della Densità (DFT) è computazionalmente efficiente ma spesso manca dell'accuratezza necessaria a causa della dipendenza dal funzionale, mentre i metodi ad alto livello di funzione d'onda (ad es. CCSD(T)) sono computazionalmente proibitivi per sistemi con metalli di transizione realistici. Inoltre, gli approcci di computazione quantistica come la Diagonalizzazione Quantistica Basata su Campionamento (SQD) sono attualmente limitati dalle restrizioni hardware a spazi attivi troppo piccoli per catturare l'intera correlazione dinamica necessaria per risultati chimicamente accurati nei sistemi catalitici. Di conseguenza, trattare simultaneamente sia la correlazione statica (all'interno dello spazio attivo) che la correlazione dinamica (dalle orbite esterne) rimane una sfida significativa.

Metodologia
Questo studio propone un framework ibrido quantistico-classico che integra la Diagonalizzazione Quantistica Basata su Campionamento (SQD) con la teoria Tailored Coupled Cluster (TCC) per affrontare i limiti della dimensione dello spazio attivo e del trattamento della correlazione.

1. Sistema Target: Un modello molecolare minimo di un catalizzatore di titanio ansa-ciclopentadienil-arene cationico è stato utilizzato per investigare le energie relative di due stati di transizione, T3 (inserzione dell'etilene) e T4 (trasferimento di beta-idruro). Per garantire un framework orbitale coerente, lo stato T3 è stato aumentato con una molecola di etilene esplicita (T3(+C2H4)), corrispondente al numero di atomi ed elettroni di T4.
2. SQD per la Correlazione Statica: La SQD è stata impiegata per trattare la struttura elettronica dello spazio attivo. Utilizzando l'ansatz LUCJ (Local-Unitary Cluster Jastrow) sul processore quantistico IBM Heron, la SQD campiona stocasticamente le bitstring per costruire un Hamiltoniano proiettato all'interno di un sottospazio. Ciò consente il trattamento di spazi attivi (fino a 30 orbitali spaziali) che sono combinatoriamente troppo grandi per la classica CASCI (Complete Active Space Configuration Interaction) classica. Il metodo utilizza il recupero delle configurazioni auto-coerente per mitigare il rumore dell'hardware e garantire corretti numeri di elettroni e spin.
3. TCC per la Correlazione Dinamica: Le funzioni d'onda derivate da SQD hanno servito come stati di riferimento per i calcoli Tailored Coupled Cluster (TCC). In questo framework, le ampiezze di cluster per lo spazio attivo ($\hat{T}_{AS}$) sono "su misura" (fissate) sulla base dei risultati SQD, catturando la correlazione statica e le eccitazioni di ordine superiore (tripli, quadrupli) all'interno dello spazio attivo. Le restanti ampiezze per lo spazio esterno ($\hat{T}_{Ext}$) sono ottimizzate tramite le equazioni standard Coupled Cluster Singles and Doubles (CCSD) per catturare la correlazione dinamica.
4. Estensioni: Lo studio ha inoltre applicato un'estensione perturbativa dei tripli, TCC(T), per migliorare l'accuratezza. È stato utilizzato uno schema di correzione dell'energia per garantire che l'energia dello spazio attivo rifletta il livello SQD piuttosto che l'approssimazione di livello inferiore CCSD.
5. Configurazione Computazionale: Tutti i calcoli single-point post-DFT sono stati eseguiti in fase gassosa per isolare gli effetti di correlazione elettronica, sebbene le ottimizzazioni della geometria includessero gli effetti del solvente (toluene) tramite il Modello del Continuo Polarizzabile (PCM).

Risultati Chiave

• Problemi di Convergenza con la sola SQD: Quando la SQD veniva applicata isolatamente, le energie relative degli stati di transizione non potevano convergere entro gli spazi attivi attualmente sostenibili. I risultati erano sensibili alla dimensione dello spazio attivo, indicando che la correlazione statica da sola è insufficiente per questo sistema.
• Successo di SQD-TCC: L'approccio ibrido SQD-TCC ha prodotto risultati ragionevolmente convergenti per le energie relative. Incorporando la correlazione dinamica tramite il framework TCC, il metodo ha bilanciato con successo il trattamento della correlazione elettronica.
• Impatto della Correlazione Statica: Una discrepanza significativa (superiore a 1 kcal/mol) è stata osservata tra i risultati standard CCSD(T) e quelli SQD-TCC(T). Ciò evidenzia che gli effetti della correlazione statica sono critici in questo sistema a base di titanio e che trascurarli (come avviene nella standard CCSD(T) a singolo riferimento) porta a errori significativi nella previsione della selettività.
• Scalabilità dello Spazio Attivo: Lo studio ha dimostrato che la SQD consente il trattamento di spazi attivi (ad es. 30 orbitali) che sono praticamente inaccessibili ai metodi CASCI classici date le risorse computazionali disponibili (fino a 192 core CPU e 4 TiB di memoria).

Significatività e Rivendicazioni
Il documento afferma che l'integrazione di funzioni d'onda multireferenziali derivate dal calcolo quantistico con la teoria classica Coupled Cluster fornisce una via percorribile verso simulazioni chimicamente accurate per complessi sistemi catalitici. Nello specifico, il lavoro dimostra che:

1. Necessità Ibrida: Sia la correlazione elettronica statica che quella dinamica sono essenziali per descrizioni affidabili della chimica sottostante nella catalisi dei metalloceni di titanio.
2. Avanzamento Metodologico: Il framework SQD-TCC colma con successo il divario tra le attuali capacità di calcolo quantistico (spazi attivi limitati) e i requisiti per simulazioni catalitiche realistiche (grandi spazi attivi + correlazione dinamica).
3. Previsione della Selettività: Il trattamento della correlazione elettronica influenza significativamente la selettività prevista della reazione, con il metodo ibrido che offre un approccio più robusto rispetto ai metodi classici ad alto livello per sistemi con forte correlazione statica.

Gli autori concludono che, sebbene l'attuale studio si concentri sulle superfici di energia potenziale in fase gassosa, la metodologia stabilisce una base per lavori futuri che incorporino contributi di energia libera ed effetti del solvente al livello della funzione d'onda correlata.