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

Planteamiento del Problema
El cofactor FeMo (FeMo-co) de la nitrogenasa es el sitio activo responsable de la fijación biológica de nitrógeno, convirtiendo $N_2$ en $NH_3$. A pesar de la resolución de su estructura atómica, su estructura electrónica sigue siendo objeto de intenso debate debido a la presencia de ocho iones de metales de transición y numerosos electrones desapareados. Esta complejidad ha llevado a la sugerencia de que el sistema está fuera del alcance de la simulación computacional clásica, lo que lo convierte en un referente principal para los algoritmos de computación cuántica. Específicamente, estimar la energía del estado fundamental de los modelos de FeMo-co con "exactitud química" (aproximadamente 1 kcal/mol o 1.6 milliHartrees) es una métrica estándar para evaluar los requisitos de recursos de los algoritmos cuánticos. El desafío radica en la naturaleza multiconfiguracional del sistema, donde un vasto número de configuraciones electrónicas y de espín en competencia (potencialmente decenas de miles) deben ser navegadas para identificar el verdadero estado fundamental.

Metodología
Los autores emplean un conjunto de métodos computacionales clásicos avanzados para abordar el problema del FeMo-co, centrándose específicamente en el modelo de espacio activo "LLDUC" de 76 orbitales/113 electrones, ampliamente utilizado. Su enfoque se basa en la idea de que el problema del estado fundamental puede caracterizarse como un problema de clasificación entre muchos estados en competencia, pero mayoritariamente simples, en lugar de un problema multiconfiguracional intrínsecamente intratable.

Los componentes metodológicos clave incluyen:

• Ruptura de Simetría y Elección de Orbitales: Los autores utilizan referencias de Hartree-Fock no restringido (UHF) y estados de simetría rota (BS) para reducir el entrelazamiento y aumentar el carácter de campo medio de los estados. Una innovación crítica es el uso de orbitales "no restringidos" (diferentes orbitales espaciales para los electrones de espín arriba y espín abajo) dentro de los cálculos de Grupo de Renormalización de Matriz de Densidad (DMRG), denominados UDMRG. Esto captura de manera significativa la correlación espín-carga directamente en el conjunto de bases, mejorando la convergencia.
• Protocolo de Filtrado: Para gestionar la explosión combinatoria de las configuraciones electrónicas (por ejemplo, 78,750 asignaciones de cargas de Fe(II)/Fe(III)), los autores desarrollaron un embudo de filtrado. Realizan cientos de cálculos UCCSD para clasificar las aproximaciones iniciales de UHF de baja energía, para luego seleccionar estos candidatos para cálculos de mayor nivel.
• Clúster Acoplado de Alto Orden (CC): Emplean la teoría de clúster acoplado no restringido (UCC) hasta excitaciones cuádruples (UCCSDTQ) y más allá. Debido al costo computacional del UCCSDTQ de espacio activo completo, utilizan un enfoque compuesto: calculando UCCSDT en el espacio completo y añadiendo correcciones de cálculos UCCSDTQ en subespacios activos más pequeños definidos por Orbitales Naturales Congelados (FNO).
• Extrapolación: Para alcanzar la exactitud química, extrapolan las energías de UDMRG a una dimensión de enlace infinita ($D \to \infty$) utilizando una forma funcional $\Delta E \sim \exp[-\kappa(\log D)^2]$. Del mismo modo, extrapolan los resultados de UCC desde los subespacios FNO hacia el espacio activo completo.
• Estimación de Consenso: La energía del estado fundamental final se deriva comparando los resultados extrapolados de UDMRG (límites superiores variacionales) con los resultados extrapolados de UCC (estimaciones no variacionales), proporcionando un valor de consenso con una incertidumbre estimada.

Resultos Clave

• Energía del Estado Fundamental: Los autores predicen que la energía del estado fundamental del modelo LLDUC es $E_0 = -22140.4091 \pm 0.0020$ Hartrees (o $\pm 1.25$ kcal/mol). Esto logra una precisión dentro del rango de la exactitud química.
• Isómeros de Espín Degenerados: Dentro de esta precisión, el estudio identifica que los isómeros de espín BS7-235 y BS8-237 son degenerados. Se encuentra que la familia BS8, que no fue considerada previamente como un candidato principal al estado fundamental en muchos estudios de DFT, es energéticamente competitiva con la familia BS7, tradicionalmente favorecida.
• Paisaje Electrónico: El paisaje de baja energía se caracteriza por un manifold denso de estados. Las brechas de energía entre los isómeros de espín más bajos son pequeñas (a menudo $< 1$ kcal/mol), lo que sugiere que las poblaciones térmicas de múltiples isómeros son posibles a temperatura ambiente.
• Más Allá del Modelo Mínimo: Al extender el estudio a espacios activos más grandes (hasta 404 orbitales) e incluir fluctuaciones ambientales de la proteína mediante dinámica molecular QM/MM, las características cualitativas del paisaje electrónico se mantienen robustas. BS7, BS8 y BS10 siguen siendo las clases de menor energía.
• Estados de Oxidación: El estudio sugiere que el clúster FeMo-co favorece cargas deslocalizadas. Diferentes medidas de estado de oxidación (densidad electrónica nuclear frente a población de valencia) producen diferentes ordenamientos de los centros de hierro, resaltando la complejidad de asignar estados de oxidación experimentalmente.

Significancia y Reivindicaciones
El artículo afirma demostrar que los métodos computacionales clásicos, guiados por la comprensión física de la estructura de los estados, pueden lograr la exactitud química para el modelo del cofactor FeMo. Este logro sirve a dos propósitos principales:

1. Referencia para la Ventaja Cuántica: Al proporcionar una estimación clásica rigurosa de la energía del estado fundamental, este trabajo establece una línea base para evaluar la ventaja potencial de los algoritmos cuánticos. Los autores calculan que, si bien los métodos clásicos actuales requieren recursos sustanciales (millones de horas de núcleos de CPU), son factibles en supercomputadoras existentes, mientras que se proyecta que los algoritmos cuánticos requerirán un tiempo similar o potencialmente menor dependiendo de los solapamientos de estado inicial.
2. Información Química: El trabajo desafía las suposiciones previas de que la estructura electrónica del FeMo-co es demasiado compleja para el tratamiento clásico. Sugiere que la complejidad surge del correcto ordenamiento de muchos estados en competencia de un solo referente, en lugar de dificultades multiconfiguracionales intrínsecas. Esto implica que métodos de un solo referente más simples y prácticos (como UCCSD) pueden ser suficientes para comprender el paisaje electrónico en representaciones realistas del cofactor, siempre que se identifiquen correctamente los estados iniciales.

Los autores concluyen que su protocolo proporciona un camino factible para resolver las ambigüedades en los estudios actuales de la teoría del funcional de la densidad y ofrece una vía hacia el modelado computacional preciso del mecanismo de reacción de la nitrogenasa, potencialmente combinado con aprendizaje automático para futuros estudios de alta precisión.