Analytic first order nonadiabatic coupling matrix elements of spin-adapted open-shell time-dependent density functional theory

Este artículo presenta la derivación, implementación y evaluación de referencia de elementos de matriz de acoplamiento no adiabático analíticos de primer orden para el método X-TDDFT adaptado al espín, demostrando que reduce significativamente los errores en comparación con el U-TDDFT estándar y proporciona conocimientos cualitativamente correctos sobre la fotofísica de sistemas de capa abierta como la porfirina de cobre(II).

Lo esencial

Imagina que intentas predecir cómo reaccionará una máquina compleja, como un trompo giratorio hecho de imanes, cuando lo sacudas. En el mundo de la química, esta "máquina" es una molécula con electrones desapareados (como radicales orgánicos o complejos de metales de transición), y el "sacudido" es la luz que incide sobre ella, provocando que salte a un estado excitado.

Los científicos utilizan una herramienta llamada Teoría del Funcional de la Densidad Dependiente del Tiempo (TDDFT) para simular estas reacciones. Piensa en la TDDFT como un pronóstico del tiempo sofisticado para moléculas. Predice cómo se mueve la molécula y cómo cambia su energía.

Sin embargo, hay un problema. La TDDFT estándar (llamémosla U-TDDFT) es como un pronóstico del tiempo que asume que el viento siempre sopla en línea recta. Funciona razonablemente bien para moléculas simples, pero para las complejas con electrones "desapareados" (como nuestro trompo de imanes giratorio), se confunde. Trata los dos "espines" de los electrones (llamémoslos Espín A y Espín B) como si fueran independientes, lo que conduce a errores. Es como intentar describir un baile donde dos parejas se sostienen de la mano, pero el pronóstico asume que están bailando solos.

La Nueva Solución: X-TDDFT

Los autores de este artículo desarrollaron una actualización llamada X-TDDFT. Esto es como un nuevo modelo meteorológico que entiende que las parejas se sostienen de la mano. Obliga a las matemáticas a respetar las reglas del "espín" de la mecánica cuántica. Ya habían utilizado esto para predecir mejor la energía y la forma de estas moléculas, pero les faltaba una pieza crucial: Elementos de Matriz de Acoplamiento No Adiabático (NACME).

¿Qué es un NACME?
Imagina que la molécula es un coche conduciendo por un camino lleno de baches.

• La Energía te dice a qué velocidad va el coche.
• Los Gradientes te dicen hacia dónde se inclina el camino.
• Los NACME te dicen qué tan probable es que el coche cambie de carril o choce contra un estado diferente.

En química, este "cambio de carril" se llama Conversión Interna (CI). Es el proceso por el cual una molécula absorbe energía, se excita y luego descarga rápidamente esa energía de nuevo al estado fundamental, a menudo liberando calor en lugar de luz. Si tu cálculo de NACME es incorrecto, podrías pensar que el coche se mantendrá en su carril, cuando en realidad está a punto de desviarse violentamente hacia una zanja.

¿Qué Hicieron?

El equipo derivó las fórmulas matemáticas para calcular estas probabilidades de "cambio de carril" (NACME) utilizando su nuevo método X-TDDFT, consciente del espín. Luego lo probaron de dos maneras:

1. La Prueba Pequeña (Radical de Formaldehído): Compararon su nuevo método contra un cálculo de "estándar de oro" superpreciso (como verificar un nuevo GPS contra un mapa satelital). Descubrieron que el método antiguo (U-TDDFT) a menudo estaba equivocado por un margen enorme, a veces desviándose en un tercio a dos tercios. El nuevo método (X-TDDFT) corrigió la mayoría de estos errores, haciendo que la predicción de la velocidad a la que la molécula "se enfría" (tasa de Conversión Interna) fuera mucho más precisa. En algunos casos, el nuevo método predijo que la velocidad de enfriamiento era 100 veces más lenta de lo que predijo el método antiguo.

2. La Prueba Grande (Porfirinas de Cobre): Examinaron moléculas complejas basadas en cobre (similares a la hemo en la sangre, pero con cobre).

   • La Vieja Visión (U-TDDFT): Predijo que cuando la molécula se excita, tiene la misma probabilidad de enfriarse directamente o tomar un desvío a través de estados intermedios.
   • La Nueva Visión (X-TDDFT): Predijo que la molécula casi nunca se enfría directamente. Casi siempre toma el desvío.
   • El Resultado: Esto cambió completamente la historia de cómo se comportan estas moléculas. El método antiguo no solo obtuvo los números ligeramente mal; obtuvo la historia mal. También arruinó la comparación entre diferentes versiones de la molécula (con diferentes decoraciones químicas), haciendo que pareciera que una versión era más rápida que otra cuando la realidad era lo contrario.

La Conclusión

El artículo concluye que para moléculas con electrones desapareados (como radicales o metales de transición), no se puede confiar en las antiguas matemáticas de "línea recta" (U-TDDFT) para predecir cómo cambian de estados de energía.

Al igual que no usarías un mapa plano para navegar una cordillera, no debes usar el antiguo método TDDFT para estas moléculas complejas. El nuevo método X-TDDFT actúa como un mapa topográfico en 3D, revelando que los "carriles" (rutas de energía) son muy diferentes de lo que se pensaba anteriormente. Esto es crucial para los científicos que intentan diseñar mejores celdas solares, LED o catalizadores, porque si no sabes en qué "carril" saltará la molécula, no puedes controlar su comportamiento.

En resumen: Los autores construyeron una mejor regla para medir cómo las moléculas "saltan" entre estados de energía. Demostraron que la vieja regla era tan inexacta que estaba contando historias completamente diferentes sobre cómo funcionan estas moléculas, especialmente para aquellas que involucran cobre y otros metales de transición.

Resumen técnico

Resumen Técnico: Elementos de Matriz de Acoplamiento No Adiabático de Primer Orden Analíticos de la Teoría del Funcional de la Densidad Dependiente del Tiempo Adaptada al Espín para Sistemas de Capa Abierta

Enunciado del Problema
La teoría del funcional de la densidad dependiente del tiempo adaptada al espín (X-TDDFT) se ha establecido recientemente como un método robusto para calcular energías de excitación y gradientes analíticos de moléculas de capa abierta (como radicales y complejos de metales de transición), ofreciendo mejoras significativas sobre la TDDFT sin restricciones (U-TDDFT) al manejar adecuadamente la contaminación de espín e incluir implícitamente excitaciones dobles específicas. Sin embargo, el efecto de la adaptación al espín sobre los elementos de matriz de acoplamiento no adiabático (NACMEs), críticos para simular tasas de conversión interna (CI) y dinámicas no adiabáticas, permanecía desconocido. Aunque existen teorías rigurosas de NACMEs en U-TDDFT, estas sufren grandes errores debido a la contaminación de espín, particularmente para estados dominados por excitaciones de capa cerrada a capa vacante (CV). Las teorías existentes de NACMEs adaptadas al espín están limitadas a Hamiltonianos semiempíricos, excitaciones específicas de inversión de espín, o carecen del formalismo necesario para sistemas generales de capa abierta. En consecuencia, la simulación precisa de la fotofísica y la fotoquímica en sistemas de capa abierta utilizando marcos conservadores del espín se vio obstaculizada por la falta de una teoría consistente de NACMEs en X-TDDFT.

Metodología
Los autores derivan e implementan NACMEs analíticos de primer orden para el método X-TDDFT, específicamente para transiciones desde el estado fundamental a estados excitados (ge) y entre estados excitados (ee). La derivación sigue un enfoque de ecuación de movimiento (EOM), reemplazando los operadores estándar de excitación y desexcitación de la aproximación de fase aleatoria (RPA) por sus contrapartes adaptadas al espín.

Los componentes teóricos clave incluyen:

1. Operadores Adaptados al Espín: El formalismo utiliza un estado de referencia de Kohn-Sham de capa abierta restringida (ROKS). El operador de excitación $O^\dagger_I$ se construye utilizando excitaciones CV(0) y CV(1) adaptadas al espín, donde esta última se mezcla con excitaciones dobles específicas ($\Psi^{\bar{t}a}_{\bar{i}t}$) para garantizar la simetría de espín adecuada.
2. Formulación Lagrangiana: Para manejar la dependencia implícita de los vectores de excitación y los coeficientes de los orbitales moleculares (MO) respecto a las coordenadas nucleares, se emplea un enfoque lagrangiano. Esto permite el cálculo de derivadas analíticas sin diferenciar explícitamente los coeficientes de los MO, incorporando términos de Pulay mediante multiplicadores de Lagrange (vectores $Z$ y $W$).
3. Términos de Corrección: Los NACMEs de X-TDDFT se expresan como los NACMEs estándar de U-TDDFT más términos de corrección específicos que surgen de:
   • El uso de una referencia adaptada al espín (ROKS) en lugar de una sin restricciones (UKS).
   • La modificación del operador de excitación CV(1) para incluir excitaciones dobles.
   • Los cambios resultantes en las matrices de densidad de transición (TDMs) y en la matriz de Casida (bloques $A$ y $B$).
4. Estrategia de Implementación: La implementación está diseñada para construirse sobre códigos existentes de NACMEs en U-TDDFT. Los autores demuestran que la sobrecarga computacional adicional es aceptable, comparable a la de los cálculos de gradientes analíticos, ya que implica principalmente resolver ecuaciones de vectores Z modificadas y calcular matrices de densidad intermedias.

Contribuciones Clave

• Derivación Teórica: El artículo proporciona la primera derivación rigurosa de NACMEs analíticos para X-TDDFT conservador del espín, abarcando tanto transiciones ge como ee.
• Implementación Algorítmica: El método se implementa en el paquete de software BDF. Aprovecha la infraestructura existente de U-TDDFT, requiriendo únicamente la adición de términos de corrección específicos a las matrices de densidad de transición y a los multiplicadores de Lagrange.
• Validación: Los autores validan el método frente a cálculos de referencia multireferencia de alto nivel (XMS-CASPT2) para el catión radical del formaldehído ($CH_2O^+$).
• Aplicación a Complejos de Metales de Transición: El método se aplica a porfirinas de cobre(II) (CuP, CuOEP, CuTPP) para simular vías de relajación de estados excitados y tasas de conversión interna.

Resultados

• Precisión de la Validación ($CH_2O^+$): En comparación con los datos de referencia XMS-CASPT2, los NACMEs de X-TDDFT reducen la desviación absoluta media (MAD) de las normas de los NACMEs de U-TDDFT en aproximadamente un 32% al 65% (dependiendo del funcional). Las mejoras más significativas se observan para estados con un carácter CV(1) importante, donde U-TDDFT a menudo falla debido a la contaminación de espín. Sin embargo, los autores señalan que para ciertos estados mixtos (por ejemplo, $1^2A_2$), X-TDDFT puede ocasionalmente arrojar errores mayores que U-TDDFT, aunque esto se atribuye a una cancelación fortuita de errores en los resultados de U-TDDFT y no a un fallo de la teoría adaptada al espín.
• Porfirinas de Cobre(II):
  • Magnitudes de NACME: X-TDDFT predice magnitudes de NACME significativamente menores (hasta un orden de magnitud) para la transición $2T_1 \to 2S_0$ en comparación con U-TDDFT. Esta reducción se atribuye a la cancelación casi perfecta de las contribuciones $\alpha$ y $\beta$ en estados CV(1) puros dentro del marco adaptado al espín, una cancelación que es incompleta en U-TDDFT debido a las diferencias en la forma de los orbitales.
  • Tasas de Conversión Interna: La reducción en las normas de NACME conduce a reducciones drásticas en las tasas de CI calculadas (hasta dos órdenes de magnitud).
  • Vías de Relajación: Para CuP no sustituido, U-TDDFT predice que la vía de relajación dominante implica una CI directa al estado fundamental y una CI indirecta a través de estados de campo de ligandos en proporciones aproximadamente iguales. En contraste, X-TDDFT predice que la vía indirecta a través de estados de campo de ligandos es dominante.
  • Efectos de Sustituyentes: Las tendencias relativas de las tasas de CI para porfirinas sustituidas (CuP, CuOEP, CuTPP) difieren entre los dos métodos. U-TDDFT predice una tendencia de CuOEP < CuTPP < CuP, mientras que X-TDDFT predice CuTPP < CuOEP < CuP. Esto indica que los errores de U-TDDFT no solo son grandes, sino también no sistemáticos, fallando al capturar tendencias relativas incluso para moléculas estrechamente relacionadas.

Significado y Afirmaciones
El artículo afirma que el desarrollo de NACMEs en X-TDDFT completa el conjunto de herramientas esencial para estudiar procesos fotofísicos y fotoquímicos de moléculas de capa abierta utilizando TDDFT adaptada al espín. Los autores enfatizan que la adaptación al espín no es meramente una corrección para energías y gradientes, sino que es crucial para los NACMEs y, en consecuencia, para la precisión de las dinámicas no adiabáticas simuladas y las tasas de CI.

El estudio concluye que confiar en NACMEs de U-TDDFT para sistemas de capa abierta puede llevar a predicciones cualitativamente incorrectas de los mecanismos de relajación y a errores cuantitativos en las constantes de velocidad que son tanto grandes como no sistemáticos. Por lo tanto, los autores recomiendan el uso de NACMEs de X-TDDFT para sistemas de capa abierta en general, incluso en casos donde las energías de excitación de U-TDDFT y X-TDDFT parecen similares. El trabajo sienta las bases para futuras aplicaciones en dinámicas moleculares no adiabáticas y el estudio de fotofísica ultrarrápida en radicales y complejos de metales de transición.