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

Resumo Técnico: Estudos de Mecânica Quântica e Simulações de Monte Carlo sobre os Potenciais de Inibição de Corrosão de Alguns Derivados de Imidazol

Definição do Problema
A corrosão representa um desafio econômico e de segurança significativo globalmente, causando falhas em infraestruturas e perdas financeiras substanciais através da degradação de metais, particularmente em ambientes industriais, marinhos e atmosféricos. Embora existam estratégias de mitigação tradicionais, a busca por inibidores de corrosão eficazes e sustentáveis permanece crítica. A determinação experimental dos mecanismos de inibição é frequentemente cara e demorada. Consequentemente, há uma necessidade de modelagem teórica eficiente para correlacionar estruturas moleculares com a eficácia de inibição, especificamente para compostos heterocíclicos como derivados de imidazol, que possuem heteroátomos doadores de elétrons (N, S, O) adequados para adsorção em superfícies metálicas.

Metodologia
Este estudo empregou uma abordagem computacional dupla para investigar os potenciais de inibição de corrosão de três derivados de imidazol (designados A, B e C) em uma superfície de ferro (Fe):

1. Teoria do Funcional da Densidade (DFT): Os cálculos foram realizados utilizando o funcional B3LYP com o conjunto de bases 6-31G(d) no pacote de software Spartan 10. O fluxo de trabalho incluiu definição de geometria, buscas conformacionais para identificar estruturas de energia mínima global e otimização geométrica. Descritores de química quântica fundamentais foram calculados, incluindo energias de orbitais moleculares de fronteira ($E_{HOMO}$, $E_{LUMO}$), lacunas de energia ($E_g$), potencial de ionização, afinidade eletrônica, dureza global, suavidade, eletronegatividade, eletrofilicidade, nucleofilicidade e a fração de transferência de elétrons ($\Delta N$). Índices de Fukui e distribuições de carga de Mulliken também foram analisados para identificar sítios ativos de adsorção.
2. Simulações de Monte Carlo (MC): O comportamento de adsorção foi simulado usando o pacote Material Studio 8.0 com o campo de força COMPASS. As simulações foram conduzidas em uma superfície Fe(110) em um ensemble NVT a 298 K por 200.000 ciclos (100.000 etapas de equilíbrio e 100.000 etapas de produção). As energias de adsorção ($E_{ads}$) foram calculadas para determinar a estabilidade da interação inibidor-metal.

Principais Resultados

• Propriedades Eletrônicas e Reatividade:

  • Lacunas de Energia: As lacunas de energia ($E_g$) calculadas foram 4,50 eV (A), 4,28 eV (B) e 4,23 eV (C). A lacuna menor em C e B sugere maior reatividade e melhor potencial de inibição em comparação a A.
  • Capacidades de Doador/Aceitador: A molécula C exibiu o $E_{HOMO}$ mais alto (-5,45 eV) e o menor potencial de ionização (5,45 eV), indicando que é o doador de elétrons mais forte. Inversamente, as moléculas A e B mostraram valores de afinidade eletrônica mais altos (1,72 eV e 1,73 eV, respectivamente) comparados a C (1,22 eV), sugerindo que A e B têm uma maior capacidade de aceitar elétrons da superfície metálica.
  • Descritores Globais: A molécula C exibiu a maior nucleofilicidade (0,32) e a menor dureza global (2,11 eV), implicando um caráter mais "suave" favorável à adsorção. A molécula B exibiu a maior eletrofilicidade global (3,61 eV).
  • Índices de Fukui: A análise identificou heteroátomos, especificamente Nitrogênio e Enxofre, como os principais sítios ativos para adsorção. Os valores da função de Fukui indicaram sítios específicos para ataques nucleofílicos e eletrofílicos, confirmando o papel desses heteroátomos no mecanismo de inibição.

• Adsorção e Simulações de Monte Carlo:

  • Energias de Adsorção: As simulações de MC revelaram energias de adsorção fortemente negativas para todos os compostos, indicando adsorção espontânea e estável. Os valores seguiram a ordem: C (-13719,09 kcal/mol) > B (-13039,52 kcal/mol) > A (-13032,45 kcal/mol).
  • Mecanismo: As energias de adsorção altamente negativas sugerem um mecanismo de quimiossorção envolvendo ligação covalente entre as moléculas inibidoras e a superfície Fe(110).
  • Configuração: As simulações mostraram que as moléculas adotam uma orientação paralela na superfície Fe(110), maximizando a cobertura superficial. Este arranjo apoia os achados de DFT de retrodoação, onde elétrons dos orbitais d preenchidos do ferro transferem-se para os orbitais vazios (LUMO) dos inibidores.

Significância e Alegações
O artigo alega que a modelagem teórica fornece uma alternativa viável e de baixo custo aos métodos experimentais para a triagem de inibidores de corrosão. O estudo conclui que, entre os três derivados investigados, a Molécula C possui as características mais favoráveis para a inibição da corrosão. Esta conclusão baseia-se em sua menor lacuna de energia, maior nucleofilicidade e na energia de adsorção mais negativa, que coletivamente sugerem uma força de adsorção e cobertura superficial superiores em aço doce.

Os autores afirmam que o mecanismo de inibição envolve um processo duplo: a doação de elétrons do HOMO do inibidor para os orbitais d vazios do ferro, e a retrodoação de elétrons dos orbitais d do ferro para o LUMO do inibidor. O estudo destaca que a presença de heteroátomos (N, S) e a distribuição eletrônica específica dentro dos derivados de imidazol são fatores críticos em sua capacidade de formar filmes protetores contra a corrosão. Os resultados confirmam que esses derivados de imidazol atuam como adsorventes químicos eficazes, com a Molécula C demonstrando o maior potencial para prevenir a oxidação em superfícies metálicas.