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

Resumen Técnico: Método de Coupled Cluster Adaptado con Diagonalización Cuántica Basada en Muestreo

Planteamiento del Problema
La predicción precisa de la selectividad en reacciones catalíticas de metalocenos basados en titanio para la producción de 1-hexeno requiere computar las diferencias de energía libre entre estados de transición competitivos (inserción de etileno frente a transferencia de beta-hidruro) con una precisión dentro de 1 kcal/mol. Los métodos actuales de estructura electrónica enfrentan una dicotomía: la Teoría del Funcional de la Densidad (DFT) es computacionalmente eficiente pero a menudo carece de la precisión necesaria debido a la dependencia del funcional, mientras que los métodos de alta resolución de la función de onda (p. ej., CCSD(T)) son computacionalmente prohibitivos para sistemas realistas de metales de transición. Además, los enfoques de computación cuántica como la Diagonalización Cuántica Basada en Muestreo (SQD) están actualmente limitados por restricciones de hardware a espacios activos demasiado pequeños para capturar la correlación dinámica completa requerida para resultados de precisión química en sistemas catalíticos. En consecuencia, tratar tanto la correlación estática (dentro del espacio activo) como la correlación dinámica (de orbitales externos) simultáneamente sigue siendo un desafío significativo.

Metodología
Este estudio propone un marco híbrido cuántico-clásico que integra la Diagonalización Cuántica Basada en Muestreo (SQD) con la teoría de Coupled Cluster Adaptado (TCC) para abordar las limitaciones del tamaño del espacio activo y el tratamiento de la correlación.

1. Sistema Objetivo: Se utilizó un modelo molecular mínimo de un catalizador de titanio ansa-ciclopentadienil-areno catiónico para investigar las energías relativas de dos estados de transición, T3 (inserción de etileno) y T4 (transferencia de beta-hidruro). Para asegurar un marco orbital consistente, el estado T3 fue aumentado con una molécula de etileno explícita (T3(+C2H4)), igualando el recuento de átomos y electrones de T4.
2. SQD para Correlación Estática: Se empleó SQD para tratar la estructura electrónica del espacio activo. Utilizando el ansatz LUCJ (Local-Unitary Cluster Jastrow) en el procesador cuántico IBM Heron, la SQD muestrea estocásticamente bitstrings para construir un Hamiltoniano proyectado dentro de un subespacio. Esto permite el tratamiento de espacios activos (hasta 30 orbitales espaciales) que son combinatoriamente demasiado grandes para la Interacción de Configuración de Espacio Activo Completo (CASCI) clásica. El método utiliza la recuperación de configuración autoconsistente para mitigar el ruido del hardware y asegurar los números correctos de electrones y espín.
3. TCC para Correlación Dinámica: Las funciones de onda derivadas de SQD sirvieron como estados de referencia para los cálculos de Tailored Coupled Cluster (TCC). En este marco, las amplitudes de cluster para el espacio activo ($\hat{T}_{AS}$) están "adaptadas" (fijas) basándose en los resultados de SQD, capturando la correlación estática y excitaciones de orden superior (triples, cuadruples) dentro del espacio activo. Las amplitudes restantes para el espacio externo ($\hat{T}_{Ext}$) se optimizan mediante las ecuaciones estándar de Coupled Cluster Singles and Doubles (CCSD) para capturar la correlación dinámica.
4. Extensiones: El estudio aplicó además una extensión de triples perturbativos, TCC(T), para mejorar la precisión. Se utilizó un esquema de corrección de energía para asegurar que la energía del espacio activo refleje el nivel de SQD en lugar de la aproximación de menor nivel CCSD.
5. Configuración Computacional: Todos los cálculos de punto único post-DFT se realizaron en fase gaseosa para aislar los efectos de la correlación electrónica, aunque las optimizaciones de geometría incluyeron efectos de solvente (tolueno) mediante el Modelo de Continuo Polarizable (PCM).

Resultados Clave

• Problemas de Convergencia con SQD por sí sola: Cuando se aplicó SQD de forma aislada, las energías relativas de los estados de transición no pudieron converger dentro de los espacios activos actualmente factibles. Los resultados fueron sensibles al tamaño del espacio activo, indicando que la correlación estática por sí sola es insuficiente para este sistema.
• Éxito de SQD-TCC: El enfoque híbrido SQD-TCC produjo resultados razonablemente convergidos para las energías relativas. Al incorporar la correlación dinámica a través del marco TCC, el método logró equilibrar con éxito el tratamiento de la correlación electrónica.
• Impacto de la Correlación Estática: Se observó una discrepancia significativa (superior a 1 kcal/mol) entre los resultados de CCSD(T) estándar y los de SQD-TCC(T). Esto resalta que los efectos de la correlación estática son críticos en este sistema basado en titanio y que omitirlos (como en el CCSD(T) de referencia única estándar) conduce a errores significativos en la predicción de la selectividad.
• Escalabilidad del Espacio Activo: El estudio demostró que la SQD permite el tratamiento de espacios activos (p. ej., 30 orbitales) que son prácticamente inaccesibles para los métodos CASCI clásicos dados los recursos computacionales disponibles (hasta 192 núcleos de CPU y 4 TiB de memoria).

Significancia y Reivindicaciones
El artículo afirma que la integración de funciones de onda multireferencia derivadas de la computación cuántica con la teoría de coupled cluster clásica proporciona una vía viable hacia simulaciones de precisión química para sistemas catalíticos complejos. Específicamente, el trabajo demuestra que:

1. Necesidad Híbrida: Tanto las correlaciones electrónicas estáticas como las dinámicas son esenciales para descripciones fiables de la química subyacente en la catálisis de metalocenos de titanio.
2. Avance Metodológico: El marco SQD-TCC cierra con éxito la brecha entre las capacidades actuales de la computación cuántica (espacios activos limitados) y los requisitos para simulaciones catalíticas realistas (espacios activos grandes + correlación dinámica).
3. Predicción de Selectividad: El tratamiento de la correlación electrónica influye significativamente en la selectividad predicha de la reacción, ofreciendo el método híbrido un enfoque más robusto que los métodos clásicos de alto nivel para sistemas con fuerte correlación estática.

Los autores concluyen que, si bien el presente estudio se centra en superficies de energía potencial en fase gaseosa, la metodología establece una base para trabajos futuros que incorporen contribuciones de energía libre y efectos de solvente al nivel de la función de onda correlacionada.