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

기술 요약: 일부 이미다졸 유도체의 부식 억제 잠재력에 관한 양자 역학적 연구 및 몬테카를로 시뮬레이션

문제 제기
부식은 전 세계적으로 인프라의 실패와 금속의 퇴화, 특히 산업, 해양 및 대기 환경에서의 상당한 재정적 손실을 초래하는 중요한 경제적 및 안전상의 과제이다. 기존의 완화 전략이 존재하지만, 효과적이고 지속 가능한 부식 억제제를 찾는 연구는 여전히 매우 중요하다. 실험적인 억제 메커니즘 결정은 종종 비용이 많이 들고 시간이 오래 걸린다. 결과적으로, 특히 N, S, O와 같이 금속 표면에 흡착하기 적합한 전자 공여 헤테로 원자를 보유한 이미다졸 유도체와 같은 헤테로고리 화합물에 대해, 분자 구조와 억제 효능 사이의 상관관계를 규명할 수 있는 효율적인 이론적 모델링이 필요하다.

연구 방법론
본 연구는 세 가지 이미다졸 유도체(A, B, C로 명명)의 철(Fe) 표면에 대한 부식 억제 잠재력을 조사하기 위해 이중 계산 접근 방식을 채택하였다:

1. 밀도 범함수 이론 (DFT): 계산은 Spartan 10 소프트웨어 패키지에서 B3LYP 범함수와 6-31G(d) 기저 집합을 사용하여 수행되었다. 워크플로우에는 기하학적 구조 정의, 전역 최소 에너지 구조를 식별하기 위한 형태 탐색, 그리고 기하학적 최적화가 포함되었다. 주요 양자 화학적 기술자로서 에너지 갭($E_g$), 이온화 에너지, 전자 친화도, 전역 경도, 연성, 전기 음성도, 친전자성, 친핵성 및 전자 전달 분율($\Delta N$)을 포함한 $E_{HOMO}$, $E_{LUMO}$ 등의 프런티어 분자 궤도 에너지가 계산되었다. 활성 흡착 부위를 식별하기 위해 푸쿠이 지수(Fukui indices)와 멀리컨 전하 분포(Mulliken charge distributions)도 분석되었다.
2. 몬테카를로 (MC) 시뮬레이션: 흡착 거동은 COMPASS 힘의 장(force field)을 사용하여 Material Studio 8.0 패키지에서 시뮬레이션되었다. 시뮬레이션은 Fe(110) 표면에서 298 K, 200,000 사이클(100,000회의 평형화 단계 및 100,000회의 생성 단계) 동안 NVT 앙상블을 사용하여 수행되었다. 억제제-금속 간 상호작용의 안정성을 결정하기 위해 흡착 에너지($E_{ads}$)를 계산하였다.

주요 결과

• 전자적 특성 및 반응성:

  • 에너지 갭: 계산된 에너지 갭($E_g$)은 A(4.50 eV), B(4.28 eV), C(4.23 eV)였다. C와 B의 낮은 갭은 A에 비해 높은 반응성과 더 나은 억제 잠재력을 시사한다.
  • 공여/수용 능력: 분자 C는 가장 높은 $E_{HOMO}$ (-5.45 eV)와 가장 낮은 이온화 전위(5.45 eV)를 나타내어 가장 강력한 전자 공여체임을 보여주었다. 반대로, A와 B는 C(1.22 eV)보다 높은 전자 친화도 값(각각 1.72 eV 및 1.73 eV)을 보여 금속 표면으로부터 전자를 받아들이는 능력이 더 강함을 시示했다.
  • 전역 기술자: 분자 C는 가장 높은 친핵성(0.32)과 가장 낮은 전역 경도(2.11 eV)를 나타냈으며, 이는 흡착에 유리한 "연한(soft)" 특성을 의미한다. 분자 B는 가장 높은 전역 친전자성(3.61 eV)을 보였다.
  • 푸쿠이 지수: 분석 결과, 질소(N)와 황(S)과 같은 헤테로 원자가 주요 활성 흡착 부위임을 확인하였다. 푸쿠이 함수 값은 친핵적 및 친전자적 공격을 위한 특정 부위를 지시하였으며, 이는 억제 메커니즘에서 이러한 헤테로 원자의 역할을 확인시켜 주었다.

• 흡착 및 몬테카를로 시뮬레이션:

  • 흡착 에너지: MC 시뮬레이션은 모든 화합물에 대해 강한 음의 흡착 에너지를 나타냈으며, 이는 자발적이고 안정적인 흡착을 의미한다. 값의 순서는 C (-13719.09 kcal/mol) > B (-13039.52 kcal/mol) > A (-13032.45 kcal/mol) 순이었다.
  • 메커니즘: 매우 낮은 음의 흡착 에너지는 억제제 분자와 Fe(110) 표면 사이의 공유 결합을 포함하는 화학 흡착 메커니즘을 시사한다.
  • 배열: 시뮬레이션 결과, 분자들이 표면 피복률을 극대화하기 위해 Fe(110) 표면에서 평행한 방향을 취하는 것으로 나타났다. 이러한 배치는 억제제의 빈 궤도(LUMO)로 철의 채워진 d-궤도에서 전자가 이동하는 역공여(back-donation)에 대한 DFT의 결과를 뒷받침한다.

의의 및 주장
본 논문은 이론적 모델링이 부식 억제제 스크리닝을 위한 실행 가능하고 비용 효율적인 실험적 방법의 대안을 제공한다고 주장한다. 본 연구는 조사된 세 가지 유도체 중 분자 C가 부식 억제에 가장 유리한 특성을 보유하고 있다고 결론짓는다. 이 결론은 낮은 에너지 갭, 높은 친핵성, 그리고 가장 낮은(음의 값이 큰) 흡착 에너지에 근거하며, 이는 종합적으로 우수한 흡착 강도와 연강(mild steel)에 대한 표면 피복력을 시사한다.

저자들은 억제 메커니즘이 억제제의 HOMO에서 철의 빈 d-궤도로 전자를 공여하는 과정과, 철의 d-궤도에서 억제제의 LUMO로 전자를 역공여하는 이중 과정에 관여한다고 단언한다. 본 연구는 헤테로 원자(N, S)의 존재와 이미다졸 유도체 내의 특정 전자 분포가 부식에 대한 보호막을 형성하는 능력에 결정적인 요소임을 강조한다. 결과는 이러한 이미다졸 유도체들이 효과적인 화학적 흡착제로서 작용하며, 특히 분자 C가 금속 표면의 녹 발생을 방지하는 데 가장 높은 잠재력을 보임을 확인해 준다.