Study of hydrogen transfer reactions within the framework of the non-equilibrium approach: CH3OH + CH3 → CH2OH + CH4, HCHO + CH3 → HCO + CH4 and CH4 + HO2 → CH3 + H2O2

Resumen Técnico: Estudio de las Reacciones de Transferencia de Hidrógeno dentro del Marco del Enfoque de No Equilibrio

Planteamiento del Problema
Este estudio investiga la cinética de tres reacciones de transferencia de átomo de hidrógeno en fase gaseosa:

1. $\text{CH}_3\text{OH} + \text{CH}_3 \rightarrow \text{CH}_2\text{OH} + \text{CH}_4$
2. $\text{HCHO} + \text{CH}_3 \rightarrow \text{HCO} + \text{CH}_4$
3. $\text{CH}_4 + \text{HO}_2 \rightarrow \text{CH}_3 + \text{H}_2\text{O}_2$

Aunque existen datos experimentales para estas reacciones, estudios teóricos previos utilizando la Teoría del Estado de Transición (TST) convencional y sus variaciones (por ejemplo, V-TST con Tunelamiento Multidimensional) han arrojado resultados mixtos. Específicamente, enfoques de no equilibrio aplicados previamente a la reacción asimétrica de metanol con el radical metilo resultaron en constantes de velocidad subestimadas y pendientes de dependencia de la temperatura incorrectas. Los autores señalan que no se habían realizado previamente estudios teóricos para la tercera reacción ($\text{CH}_4 + \text{HO}_2$) dentro de este marco. El problema central abordado es hasta qué punto los patrones observados en reacciones simétricas (como metano + metilo) u otros sistemas asimétricos (metano + hidroxilo) se aplican a estas reacciones de transferencia de hidrógeno asimétricas específicas, particularmente en lo que respecta al papel del tunelamiento y la reorganización del sistema.

Metodología
La investigación emplea un enfoque de no equilibrio basado en el Principio de Franck–Condon Generalizado (GFCP), adaptado para la transferencia de H en fase gaseosa. La metodología difiere de la TST estándar al tratar la reacción como un proceso de múltiples etapas que involucra un complejo de colisión formado cuando la velocidad relativa de los reactivos llega a cero.

• Mecanismo de Reacción: El modelo postula un mecanismo de tres etapas:
  1. Formación de un complejo de colisión.
  2. Aproximación de átomos pesados durante el movimiento vibracional.
  3. Tunelamiento de H a una distancia fija ($Q$) entre los átomos pesados, acompañado de una reorganización estructural simultánea (específicamente la rehibridación del átomo de carbono en el fragmento $\text{CH}_3$).
• Cálculo de la Constante de Velocidad: La constante de velocidad térmica $k(T)$ se calcula como una integral sobre la distancia $Q$ entre los átomos donante y aceptor. La ecuación incorpora la frecuencia de tunelamiento ($\nu_t$), la energía libre de activación ($\Delta G^*$) y correcciones para:
  • Asimetría: Se aplica un factor de asimetría ($\eta$) derivado del modelo WKB de Brickman al tiempo de tunelamiento para considerar la posible asimetría de la barrera.
  • Efecto Promotor: Se asume que la vibración del enlace A(H)...B en el complejo de colisión promueve la reacción al reducir la distancia de tunelamiento. Esto se modela sumando las contribuciones de varios niveles vibracionales ($j$).
  • Tiempo de Vida de la Colisión: Se aplica un factor de corrección ($m$) basado en la ecuación de Bunker para tener en cuenta el tiempo de vida finito ($\tau_{lf}$) del complejo de colisión en relación con el periodo de oscilación del enlace promotor.
• Detalles Computacionales: Los cálculos se realizaron utilizando Gaussian 03 al nivel CCSD(T)/6-311++G**//B3LYP/6-31+G**. Las frecuencias de tunelamiento se computaron utilizando la aproximación WKB con el método de Brickman. Los parámetros geométricos fueron optimizados parcialmente fijando distancias de enlace y ángulos específicos correspondientes al pozo izquierdo, la cima de la barrera y el pozo derecho.

Resultados Clave
El estudio presenta las constantes de velocidad térmica calculadas para las tres reacciones a través de un rango de temperatura de 100–2000 K, comparándolas con datos experimentales y trabajos teóricos previos.

• Reacción 1 ($\text{CH}_3\text{OH} + \text{CH}_3$): Los cálculos básicos sin efectos promotores ($k_{00}$) subestimaron los valores experimentales por un factor de 6–8 en el rango de 370–550 K. La introducción del efecto promotor (niveles vibracionales $j=0$ a $5$) y la corrección del tiempo de vida de la colisión ($k_3'$) llevaron la curva teórica a un buen acuerdo con los datos experimentales. Los resultados mostraron una estrecha concordancia con los cálculos de alta precisión de V-TST de Alecu y Truhlar [3], particularmente a temperaturas superiores a 300 K.
• Reacción 2 ($\text{HCHO} + \text{CH}_3$): Similar a la Reacción 1, el modelo base subestimó las velocidades significativamente (relación $k_{exp}/k_{00} \approx 20$ a 300 K). La inclusión de vibraciones promotoras y la corrección de tiempo de vida ($k_4'$) reprodujo con éxito la dependencia de la temperatura experimental. Los resultados fueron comparables a los obtenidos mediante el método CVT/SCT [6], lo que sugiere que el modelo de no equilibrio captura efectos que la TST/SCT estándar puede tratar de manera diferente.
• Reacción 3 ($\text{CH}_4 + \text{HO}_2$): Esta reacción es fuertemente exotérmica con una barrera alta. Se encontró que el efecto promotor era más débil en comparación con las otras dos reacciones, particularmente en niveles vibracionales altos. La constante de velocidad corregida ($k_5'$) mostró un acuerdo ligeramente mejor con el experimento que los valores no corregidos, aunque la relación entre la tasa experimental y la tasa teórica base se mantuvo alta ($\approx 7.3$ a 2000 K).
• Tunelamiento vs. Vibración: El análisis de la relación entre el tiempo de tunelamiento ($\tau_t$) y el periodo de oscilación del enlace A(H)...B ($T_{col}$) reveló que para las Reacciones 1 y 2, $\tau_t > T_{col}$ (relación $b > 1$), lo que desafía la suposición de una distancia $Q$ estrictamente fija. Para la Reacción 3, $b < 1$, lo que respalda más fuertemente la suposición de distancia fija.

Significancia y Afirmaciones
El artículo afirma que el enfoque de no equilibrio, cuando se modifica para incluir la reorganización simultánea del sistema durante el tunelamiento, la promoción vibracional y los efectos del tiempo de vida del complejo de colisión, es factible para describir la cinética de las reacciones de transferencia de H-átomo asimétricas en fase gaseosa.

Los autores sostienen que su método ofrece un equilibrio de simplicidad relativa, costo computacional moderado y una precisión satisfactoria al reproducir los datos experimentales. Concluyen que, si bien la suposición de una distancia $Q$ estrictamente fija no es totalmente válida (ya que los tiempos de tunelamiento a menudo exceden los periodos vibracionales), el concepto de un aislamiento dinámico de la vibración A(H)...B permite el uso de una distancia efectiva promedio ponderada ($Q_{eff}$) como una aproximación útil. El estudio confirma que los patrones encontrados en el trabajo previo sobre reacciones de metano-hidroxilo [7] retienen en gran medida su importancia para estos otros sistemas asimétricos, siempre que se apliquen las correcciones específicas de promoción y tiempo de vida.