Classical computational simulation of the FeMo-cofactor model to chemical accuracy and its implications

Problematica
Il cofattore FeMo (FeMo-co) della nitrogenasi è il sito attivo responsabile della fissazione biologica dell'azoto, convertendo $N_2$ in $NH_3$. Nonostante la risoluzione della sua struttura atomica, la sua struttura elettronica rimane oggetto di intenso dibattito a causa della presenza di otto ioni di metalli di transizione e di numerosi elettroni spaiati. Questa complessità ha portato al suggerimento che il sistema si trovi oltre la portata della simulazione computazionale classica, rendendolo un benchmark primario per gli algoritmi di calcolo quantistico. Specificamente, stimare l'energia dello stato fondamentale dei modelli FeMo-co con "accuratezza chimica" (circa 1 kcal/mol o 1,6 milliHartree) è una metrica standard per valutare i requisiti di risorse degli algoritmi quantistici. La sfida risiede nella natura multiconfigurazionale del sistema, dove un vasto numero di configurazioni elettroniche e di spin competitivi (potenzialmente decine di migliaia) deve essere navigato per identificare il vero stato fondamentale.

Metodologia
Gli autori impiegano una suite di avanzati metodi computazionali classici per affrontare il problema del FeMo-co, concentrandosi specificamente sul modello di spazio attivo "LLDUC" a 76 orbitali/113 elettroni, ampiamente utilizzato. Il loro approccio si basa sull'intuizione che il problema dello stato fondamentale possa essere caratterizzato come un problema di classificazione tra molti stati competitivi, ma in gran parte semplici, piuttosto che come un problema intrinsecamente intratabile e multiconfigurazionale.

I componenti metodologici chiave includono:

• Rottura della Simmetria e Scelta degli Orbitali: Gli autori utilizzano riferimenti di Hartree-Fock non ristretti (UHF) e stati a simmetria rotta (broken-symmetry, BS) per ridurre l'entanglement e aumentare il carattere mean-field degli stati. Un'innovazione critica è l'uso di orbitali "non ristretti" (orbitali spaziali differenti per gli elettroni spin-up e spin-down) all'interno dei calcoli Density Matrix Renormalization Group (DMRG) (UDMRG). Ciò cattura significativamente la correlazione spin-carica direttamente nel set di basi, migliorando la convergenza.
• Protocollo di Filtraggio: Per gestire l'esplosione combinatoria delle configurazioni elettroniche (ad esempio, 78.750 assegnazioni di cariche Fe(II)/Fe(III)), gli autori hanno sviluppato un imbuto di filtraggio. Effettuano centinaia di calcoli UCCSD per classificare i guess iniziali UHF a bassa energia, per poi selezionare questi candidati per calcoli di livello superiore.
• Coupled Cluster ad Alto Ordine (CC): Impiegano la teoria Coupled Cluster non ristretta (UCC) fino alle eccitazioni quadruple (UCCSDTQ) e oltre. A causa del costo computazionale della full active space UCCSDTQ, utilizzano un approccio composito: calcolano l'UCCSDT nello spazio completo e aggiungono le correzioni da calcoli UCCSDTQ in sottospazi attivi più piccoli definiti da Frozen Natural Orbitals (FNO).
• Estrapolazione: Per raggiungere l'accuratezza chimica, estrapolano le energie UDMRG a dimensione di legame infinita ($D \to \infty$) utilizzando una forma funzionale $\Delta E \sim \exp[-\kappa(\log D)^2]$. Allo stesso modo, estrapolano i risultati UCC dai sottospazi FNO allo spazio attivo completo.
• Stima di Consenso: L'energia dello stato fondamentale finale è derivata confrontando i risultati estrapolati UDMRG (limiti superiori variazionali) con i risultati estrapolati UCC (stime non variazionali), fornendo un valore di consenso con incertezza stimata.

Risultati Chiave

• Energia dello Stato Fondamentale: Gli autori prevedono che l'energia dello stato fondamentale del modello LLDUC sia $E_0 = -22140.4091 \pm 0.0020$ Hartree (ovvero $\pm 1.25$ kcal/mol). Ciò raggiunge un'accuratezza entro l'intervallo dell'accuratezza chimica.
• Isomeri di Spin Degeneri: Entro questa accuratezza, lo studio identifica che gli isomeri di spin BS7-235 e BS8-237 sono degeneri. La famiglia BS8, che non era precedentemente considerata un candidato primario per lo stato fondamentale in molti studi DFT, si rivela essere energeticamente competitiva con la tradizionalmente favorita famiglia BS7.
• Paesaggio Elettronico: Il paesaggio a bassa energia è caratterizzato da un manifold denso di stati. I gap energetici tra i più bassi isomeri di spin sono piccoli (spesso $< 1$ kcal/mol), suggerendo che le popolazioni termiche di molteplici isomeri siano possibili a temperatura ambiente.
• Oltre il Modello Minimo: Estendendo lo studio a spazi attivi più grandi (fino a 404 orbitali) e includendo le fluttuazioni ambientali proteiche tramite dinamica molecolare QM/MM, le caratteristiche qualitative del paesaggio elettronico rimangono robuste. BS7, BS8 e BS10 rimangono le classi a più bassa energia.
• Stati di Ossidazione: Lo studio suggerisce che il cluster FeMo-co favorisce cariche delocalizzate. Diverse misure degli stati di ossidazione (densità elettronica nucleare rispetto alla popolazione di valenza) producono diversi ordinamenti dei centri di ferro, evidenziando la complessità nell'assegnare gli stati di ossidazione sperimentalmente.

Significatività e Rivendicazioni
L'articolo sostiene che i metodi computazionali classici, quando guidati dall'intuizione fisica sulla struttura degli stati, possono raggiungere l'accuratezza chimica per il modello del cofattore FeMo. Questo raggiungimento serve a due scopi principali:

1. Benchmark per il Vantaggio Quantistico: Fornendo una stima classica rigorosa dell'energia dello stato fondamentale, questo lavoro stabilisce una linea di base per valutare il potenziale vantaggio degli algoritmi quantistici. Gli autori calcolano che, mentre gli attuali metodi classici richiedono risorse considerevoli (milioni di ore/core CPU), essi sono fattibili su supercomputer esistenti, mentre gli algoritmi quantistici sono progettati per richiedere tempi simili o potenzialmente inferiori, a seconda delle sovrapposizioni iniziali degli stati.
2. Intuizione Chimica: Il lavoro sfida le assunzioni precedenti secondo cui la struttura elettronica del FeMo-co sarebbe troppo complessa per il trattamento classico. Suggerisce che la complessità derivi dalla corretta classificazione di molti stati simili a un singolo riferimento piuttosto che da difficoltà multiconfigurazionali intrinseche. Ciò implica che metodi più semplici e pratici a singolo riferimento (come UCCSD) possono essere sufficienti per comprendere il paesaggio elettronico in rappresentazioni realistiche del cofattore, a condizione che vengano identificati gli stati iniziali corretti.

Gli autori concludono che il loro protocollo fornisce una via percorribile per risolvere le ambiguità negli attuali studi di teoria del funzionale della densità (DFT) e offre una strada verso la modellazione computazionale precisa del meccanismo di reazione della nitrogenasi, potenzialmente combinata con l'apprendimento automatico per futuri studi ad alta accuratezza.