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

Sintesi Tecnica: Studi di Meccanica Quantistica e Simulazioni Monte Carlo sul Potenziale Inibitore della Corrosione di Alcuni Derivati dell'Imidazolo

Problematica
La corrosione rappresenta una sfida economica e di sicurezza significativa a livello globale, causando il cedimento delle infrastrutture e sostanziali perdite finanziarie attraverso la degradazione dei metalli, particolarmente in ambienti industriali, marini e atmosferici. Sebbene esistano strategie di mitigazione tradizionali, la ricerca di inibitori della corrosione efficaci e sostenibili rimane critica. La determinazione sperimentale dei meccanismi di inibizione è spesso costosa e richiede molto tempo. Di conseguenza, esiste la necessità di una modellazione teorica efficiente per correlare le strutture molecolari con l'efficacia inibitoria, specificamente per i composti eterociclici come i derivati dell'imidazolo, che possiedono eteroatomi donatori di elettroni (N, S, O) adatti all'adsorbimento sulle superfici metalliche.

Metodologia
Questo studio ha impiegato un doppio approccio computazionale per investigare i potenziali di inibizione della corrosione di tre derivati dell'imidazolo (designati come A, B e C) su una superficie di ferro (Fe):

1. Teoria del Funzionale della Densità (DFT): I calcoli sono stati eseguiti utilizzando il funzionale B3LYP con il set di basi 6-31G(d) nel pacchetto software Spartan 10. Il flusso di lavoro ha incluso la definizione della geometria, la ricerca conformazionale per identificare le strutture a energia minima globale e l'ottimizzazione geometrica. Sono stati calcolati i principali descrittori chimico-quantistici, tra cui le energie degli orbitali molecolari di frontiera ($E_{HOMO}$, $E_{LUMO}$), i gap energetici ($E_g$), il potenziale di ionizzazione, l'affinità elettronica, la durezza globale, la morbidezza, l'elettronegatività, l'elettrofilia, la nucleofilia e la frazione di trasferimento elettronico ($\Delta N$). Sono stati analizzati anche gli indici di Fukui e le distribuzioni di carica di Mulliken per identificare i siti attivi di adsorbimento.
2. Simulazioni Monte Carlo (MC): Il comportamento di adsorbimento è stato simulato con il pacchetto Material Studio 8.0 utilizzando il campo di forza COMPASS. Le simulazioni sono state condotte su una superficie Fe(110) in un ensemble NVT a 298 K per 200.000 cicli (100.000 passi di equilibratura e 100.000 passi di produzione). Le energie di adsorbimento ($E_{ads}$) sono state calcolate per determinare la stabilità dell'interazione inibitore-metallo.

Risultati Chiave

• Proprietà Elettroniche e Reattività:

  • Gap Energetici: I gap energetici ($E_g$) calcolati sono stati di 4,50 eV (A), 4,28 eV (B) e 4,23 eV (C). Il gap inferiore in C e B suggerisce una maggiore reattività e un migliore potenziale inibitorio rispetto ad A.
  • Capacità di Donatore/Accettore: La Molecola C ha esibito l'$E_{HOMO}$ più elevata (-5,45 eV) e il potenziale di ionizzazione più basso (5,45 eV), indicando che è il più forte donatore di elettroni. Al contrario, le Molecole A e B hanno mostrato valori di affinità elettronica più elevati (1,72 eV e 1,73 eV rispettivamente) rispetto a C (1,22 eV), suggerendo che A e B hanno una maggiore capacità di accettare elettroni dalla superficie metallica.
  • Descrittori Globali: La Molecola C ha mostrato la nucleofilia più alta (0,32) e la durezza globale più bassa (2,11 eV), implicando un carattere più "morbido" favorevole all'adsorbimento. La Molecola B ha esibito l'elettrofilia globale più elevata (3,61 eV).
  • Indici di Fukui: L'analisi ha identificato gli eteroatomi, specificamente l'Azoto e lo Zolfo, come i principali siti attivi per l'adsorbimento. I valori della funzione di Fukui hanno indicato siti specifici per attacchi nucleofili ed elettrofili, confermando il ruolo di questi eteroatomi nel meccanismo di inibizione.

• Adsorbimento e Simulazioni Monte Carlo:

  • Energie di Adsorbimento: Le simulazioni MC hanno rivelato energie di adsorbimento fortemente negative per tutti i composti, indicando un adsorbimento spontaneo e stabile. I valori hanno seguito l'ordine: C (-13719,09 kcal/mol) > B (-13039,52 kcal/mol) > A (-13032,45 kcal/mol).
  • Meccanismo: Le energie di adsorbimento altamente negative suggeriscono un meccanismo di chemisorbimento che coinvolge legami covalenti tra le molecole inibitrici e la superficie Fe(110).
  • Configurazione: Le simulazioni hanno mostrato che le molecole adottano un orientamento parallelo sulla superficie Fe(110), massimizzando la copertura superficiale. Questa disposizione supporta i risultati DFT di retrodonazione, dove gli elettroni dagli orbitali d riempiti del ferro si trasferiscono agli orbitali vuoti (LUMO) degli inibitori.

Significato e Rivendicazioni
L'articolo sostiene che la modellazione teorica fornisce un'alternativa valida e conveniente ai metodi sperimentali per lo screening degli inibitori della corrosione. Lo studio conclude che, tra i tre derivati investigati, la Molecola C possiede le caratteristiche più favorevoli per l'inibizione della corrosione. Questa conclusione si basa sul suo gap energetico inferiore, la sua nucleofilia più elevata e la sua energia di adsorbimento più negativa, che collettivamente suggeriscono una forza di adsorbimento e una copertura superficiale superiori sull'acciaio dolce.

Gli autori affermano che il meccanismo di inibizione prevede un processo duplice: la donazione di elettroni dall'HOMO dell'inibitore agli orbitali d vuoti del ferro, e la retrodonazione di elettroni dagli orbitali d del ferro al LUMO dell'inibitore. Lo studio evidenzia che la presenza di eteroatomi (N, S) e la specifica distribuzione elettronica all'interno dei derivati dell'imidazolo sono fattori critici nella loro capacità di formare film protettivi contro la corrosione. I risultati confermano che questi derivati dell'imidazolo agiscono come efficaci adsorbenti chimici, con la Molecola C che dimostra il più alto potenziale per prevenire la ruggine sulle superfici metalliche.