Quantum Mechanical Studies And Monte Carlo Simulations On The Corrosion Inhibitive Potentials Of Some Imidazole Derivatives

Resumen Técnico: Estudios de Mecánica Cuántica y Simulaciones de Monte Carlo sobre los Potenciales de Inhibición de la Corrosión de Algunos Derivados de Imidazol

Planteamiento del Problema
La corrosión representa un desafío económico y de seguridad significativo a nivel mundial, causando fallas en la infraestructura y pérdidas financieras sustanciales debido a la degradación de metales, particularmente en entornos industriales, marinos y atmosféricos. Si bien existen estrategias de mitigación tradicionales, la búsqueda de inhibidores de corrosión efectivos y sostenibles sigue siendo crítica. La determinación experimental de los mecanismos de inhibición suele ser costosa y lenta. En consecuencia, existe la necesidad de un modelado teórico eficiente para correlacionar las estructuras moleculares con la eficacia de inhibición, específicamente para compuestos heterocíclicos como los derivados de imidazol, que poseen heteroátomos donadores de electrones (N, S, O) adecuados para la adsorción en superficies metálicas.

Metodología
Este estudio empleó un enfoque computacional dual para investigar los potenciales de inhibición de la corrosión de tres derivados de imidazol (designados como A, B y C) sobre una superficie de hierro (Fe):

1. Teoría del Funcional de la Densidad (DFT): Los cálculos se realizaron utilizando el funcional B3LYP con el conjunto de bases 6-31G(d) en el paquete de software Spartan 10. El flujo de trabajo incluyó la definición de la geometría, búsquedas conformacionales para identificar las estructuras de energía mínima global y optimización geométrica. Se calcularon descriptores de la química cuántica clave, incluyendo energías de orbitales moleculares frontera ($E_{HOMO}$, $E_{LUMO}$), brechas de energía ($E_g$), potencial de ionización, afinidad electrónica, dureza global, suavidad, electronegatividad, electrofilicidad, nucleofilicidad y la fracción de transferencia de electrones ($\Delta N$). También se analizaron los índices de Fukui y las distribuciones de carga de Mulliken para identificar los sitios activos de adsorción.
2. Simulaciones de Monte Carlo (MC): El comportamiento de adsorción se simuló utilizando el paquete Material Studio 8.0 con el campo de fuerza COMPASS. Las simulaciones se llevaron a cabo en una superficie de Fe(110) utilizando un ensamble NVT a 298 K durante 200,000 ciclos (100,000 pasos de equilibración y 100,000 pasos de producción). Las energías de adsorción ($E_{ads}$) se calcularon para determinar la estabilidad de la interacción inhibidor-metal.

Resultados Clave

• Propiedades Electrónicas y Reactividad:

  • Brechas de Energía: Las brechas de energía ($E_g$) calculadas fueron 4.50 eV (A), 4.28 eV (B) y 4.23 eV (C). La brecha menor en C y B sugiere una mayor reactividad y un mejor potencial de inhibición en comparación con A.
  • Capacidades Donadoras/Aceptoras: La Molécula C exhibió el $E_{HOMO}$ más alto (-5.45 eV) y el potencial de ionización más bajo (5.45 eV), lo que indica que es el donador de electrones más fuerte. Por el contrario, las Moléculas A y B mostraron valores de afinidad electrónica más altos (1.72 eV y 1.73 eV, respectivamente) en comparación con C (1.22 eV), lo que sugiere que A y B tienen una mayor capacidad para aceptar electrones de la superficie metálica.
  • Descriptores Globales: La Molécula C mostró la nucleofilicidad más alta (0.32) y la dureza global más baja (2.11 eV), lo que implica un carácter más "suave" favorable para la adsorción. La Molécula B exhibió la electrofilicidad global más alta (3.61 eV).
  • Índices de Fukui: El análisis identificó heteroátomos, específicamente Nitrógeno y Azufre, como los principales sitios activos para la adsorción. Los valores de la función de Fukui indicaron sitios específicos para ataques nucleofílicos y electrofílicos, confirmando el papel de estos heteroátomos en el mecanismo de inhibición.

• Adsorción y Simulaciones de Monte Carlo:

  • Energías de Adsorción: Las simulaciones de MC revelaron energías de adsorción fuertemente negativas para todos los compuestos, lo que indica una adsorción espontánea y estable. Los valores siguieron el orden: C (-13719.09 kcal/mol) > B (-13039.52 kcal/mol) > A (-13032.45 kcal/mol).
  • Mecanismo: Las energías de adsorción altamente negativas sugieren un mecanismo de quimisorción que involucra enlaces covalentes entre las moléculas del inhibidor y la superficie de Fe(110).
  • Configuración: Las simulaciones mostraron que las moléculas adoptan una orientación paralela sobre la superficie de Fe(110), maximizando la cobertura superficial. Este arreglo respalda los hallazgos de DFT de retrodonación, donde los electrones de los orbitales d llenos del hierro se transfieren a los orbitales vacíos (LUMO) de los inhibidores.

Significancia y Afirmaciones
El artículo afirma que el modelado teórico proporciona una alternativa viable y rentable a los métodos experimentales para el cribado de inhibidores de corrosión. El estudio concluye que, entre los tres derivados investigados, la Molécula C posee las características más favorables para la inhibición de la corrosión. Esta conclusión se basa en su menor brecha de energía, su mayor nucleofilicidad y su energía de adsorción más negativa, lo que colectivamente sugiere una fuerza de adsorción y una cobertura superficial superiores.

Los autores sostienen que el mecanismo de inhibición involucra un proceso dual: la donación de electrones desde el HOMO del inhibidor hacia los orbitales d vacantes del hierro, y la retrodonación de electrones desde los orbitales d del hierro hacia el LUMO del inhibidor. El estudio destaca que la presencia de heteroátomos (N, S) y la distribución electrónica específica dentro de los derivados de imidazol son factores críticos en su capacidad para formar películas protectoras contra la corrosión. Los resultados confirman que estos derivados de imidazol actúan como adsorbentes químicos efectivos, demostrando la Molécula C el mayor potencial para prevenir la oxidación en superficies metálicas.