量子化学計算（NWChem）との連携による表面防錆剤分子の性能評価#

Advance/PHASE Logo

防錆剤（腐食抑制剤）の開発において、分子が金属表面にどのように吸着し、保護膜を形成するかを原子・電子レベルで理解することは重要です。マテリアルズ・インフォマティクス（MI）の進展により、膨大な分子構造から有望な候補を効率的に見つけ出す「ハイスループット・スクリーニング」が求められています。本解析では、第一原理計算ソフトウェアAdvance/PHASEの表面計算機能と、オープンソースの量子化学計算ソフトNWChemを連携させ、低コストかつ物理的整合性の高い防錆剤性能評価フローを構築しました。

Keywords: 第一原理計算 (DFT), NWChem, ハイスループットスクリーニング, 仕事関数, ESM法, 電子移動モデル ( $\Delta N$ ), アゾール化合物

理論背景と計算手順#

電子移動モデル（ΔN）#

防錆剤分子が金属表面に吸着する際、分子から金属へと電子が移動して化学結合が形成されます。この際の電子の受け渡し能力を定量化する指標として、PearsonのHSAB則に基づき、ParrおよびPearsonによって1983年に提唱された「移行電子数（ $\Delta N$ ）」が用いられます。この指標は、Lukovitsら [1] やKokalj [2] らの研究によって腐食抑制剤の性能評価に導入され、その物理的妥当性と有用性が広く検証されてきました。

移行電子数の算出式
$\Delta N = \frac{\Phi_{\text{metal}} - \chi_{\text{inh}}}{2(\eta_{\text{metal}} + \eta_{\text{inh}})} \approx \frac{\Phi_{\text{metal}} - \chi_{\text{inh}}}{2\eta_{\text{inh}}}$

ここで、 $\Phi_{\text{metal}}$ は金属の仕事関数、 $\chi_{\text{inh}}$ は分子の電気陰性度、 $\eta_{\text{inh}}$ は分子の化学的硬度です。バルク金属では $\eta_{\text{metal}} \approx 0$ と近似できます。なお、計算上 $\chi_{\text{inh}}$ と $\eta_{\text{inh}}$ は、分子の最高被占軌道（HOMO）および最低空軌道（LUMO）のエネルギーから以下のように近似的に算出されます。

$\chi_{\text{inh}} \approx -(E_{\text{HOMO}} + E_{\text{LUMO}})/2$

$\eta_{\text{inh}} \approx (E_{\text{LUMO}} - E_{\text{HOMO}})/2$
Electron Transfer Model
図1. 電子移動モデルの模式図。金属のフェルミ準位（ $-\Phi_{\text{metal}}$ ）よりも、分子の電子を保持するポテンシャル（ $-\chi_{\text{inh}}$ ）が高い場合、電子は分子から金属へと流れ込みます。

$\Delta N$ から何が分かるか？
$\Delta N$ が正の値（ $\Delta N > 0$ ）である場合、分子から金属の空軌道へ電子が供与され、配位結合的な吸着が進行しやすいことを示します。一般に、同族の分子間や類似構造を持つ異性体間では、 $\Delta N$ が大きいほど金属表面との相互作用（電子供与による結合）が強く、防錆性能が高まる傾向があると解釈されます [1]。

本手法の鍵は、金属側の物性である仕事関数をAdvance/PHASEで、分子側の物性をNWChem [3] でそれぞれ個別に高精度に算出することにあります。これにより、吸着系の高コストな計算を回避しつつ、合理的なスクリーニングが可能となります。

ESM法による高精度な仕事関数の算出#

分子の電子授受能力を評価する電子移動モデルにおいて、基準となる金属表面の「仕事関数（ $\Phi$ ）」の精度は結果に直結します。通常の周期境界条件（PBC）を用いたスラブモデルの計算では、スラブ間に生じる人工的な双極子相互作用により、真空準位の決定が不正確になることがあります。

Advance/PHASEに実装されているESM（Effective Screening Medium）法を利用すると、真空領域のポテンシャルを完全に平坦化することが可能です。本解析では、Cu(100)表面の仕事関数として、通常のPBC計算による4.10 eVではなく、実験値（4.65 eV）に非常に近い、ESM法による算出値（4.60 eV）を採用しました (既存事例 [4] 参照)。

NWChemを用いたハイスループット量子化学計算#

30種類の有機分子（アゾール系防錆剤、各種ヘテロ環、および比較用の代表的な単環式芳香族化合物）を対象に、Advance/PHASEに同梱されているPythonライブラリpymatgen [5] に含まれるpymatgen.io.nwchem モジュールを用いた自動実行パイプラインを構築し、分子の電子状態パラメータ（ $\chi$ , $\eta$ ）を算出しました。ハイブリッド汎関数（B3LYP/6-311++G**）を用いることで、絶対エネルギーの正確な評価を行っています。

振動計算による「真の安定構造」の検証#

自動スクリーニングを構築する上で極めて重要なのが、真の安定構造の確認です。構造最適化計算は、エネルギー勾配がゼロになった点（停留点）で終了しますが、それが「真の安定状態（極小点）」なのか「自然界に存在し得ない鞍点（遷移状態）」なのかは判別できません。本パイプラインでは、構造最適化の後に振動計算を自動実行し、「虚数振動」が存在しない（すべて正の振動数である）ことを確認しています。これにより、抽出されたHOMO/LUMOエネルギーが「物理的に意味のある真の安定構造」のものであることを保証し、MIにおけるデータセットの信頼性を担保しています。

計算結果と考察#

1. 検証済みの28分子におけるスクリーニング結果#

SMILES記号から生成した30種類の分子を計算した結果、2種類の分子において虚数振動が確認されたため、本スクリーニングフローでは除外としました。なお、実際のより高度な探索ワークフローにおいては、この虚数振動モードに沿って構造を微小変位させ、真の安定構造を自動探索（再最適化）するフェーズをパイプラインに組み込むことも可能です。図2は振動計算をパスした28分子の物性値の分布を示します。

Chi vs Eta Scatter Plot
図2. 28種の有機分子における絶対電気陰性度（ $\chi$ ）と化学的硬度（ $\eta$ ）の分布。赤の破線はESM法によるCu(100)の仕事関数（ $\Phi = 4.60$ eV）を示します。この線より左側（オレンジ背景）の領域にある分子が $\Delta N > 0$ となり、防錆効果が期待できる候補となります。

2. 移行電子数（ΔN）による有望候補の序列化#

Advance/PHASEのESM法で算出した Cu(100) の仕事関数（ $\Phi = 4.60$ eV）を基準とした $\Delta N$ の算出結果を図3に示します。不正確なPBC値（4.10 eV）を用いた場合、実績のある防錆剤でさえ $\Delta N$ が負となりモデルが破綻してしまいますが、ESM法による高精度な基準を設けることで、物理的に妥当なスクリーニングが可能となっています。

Delta N Bar Plot
図3. 28分子の移行電子数（ $\Delta N$ ）ランキング。緑色は防錆剤として代表的、あるいは特筆すべき高い電子供与性を示した分子を表しています。

【考察】
本スクリーニングの結果、ATA（3-amino-1,2,4-triazole）やチオール基を持つMBTが高い電子供与ポテンシャル（ $\Delta N$ ）を示し、これらが銅に対して優れた防錆性能を持つという実験事実と概ね整合する結果が得られました。また、アニリン（aniline）が全分子中で最高の値を示しましたが、これは芳香族アミンとしての高い電子供与性を反映したものです。ただし、アニリンは強い毒性を持ち、分子サイズが小さいため強固な疎水性バリアを形成しにくいといった実用上の課題があり、実際の材料選定ではこうした物理的制約も併せて考慮する必要があります。

一方で、代表的な防錆剤であるBTAH（1h_benzotriazole）の $\Delta N$ は中程度に留まりました。本モデルの予測（ATA > BTAH）と実際の実験事実（BTAH ≥ ATA）[6] の間に見られるこのような逆転現象は、本手法の評価指標が「単体分子」の電子状態のみに依存していることに起因します。実際の金属表面では、BTAHのように平面構造を活かした $\pi$ 相互作用や、金属イオンとの錯体形成が吸着安定性に大きく寄与する場合があり、これらは分子単体のパラメータだけでは捉えきれません。

発展的な考察：単一評価の限界と吸着シミュレーション#

今回の $\Delta N$ モデルは、低コストに有望候補を絞り込む「ハイスループット・スクリーニング」として有効です。しかし、先行研究 [6, 7] でも指摘されている通り、実環境における性能をより正確に予測するためには、以下の要因を詳細に評価する段階的なアプローチが不可欠です。

幾何学的な吸着構造の影響： 分子が表面に平行に吸着して $\pi$ 電子系が関与するか、垂直に吸着して特定の官能基が結合するかといった、立体的な吸着形態の違い
溶媒および競合効果： 実際の腐食環境における水分子の排除能力や、塩素イオン（Cl^-）等の腐食因子との吸着サイトの奪い合い

したがって、本スクリーニングで上位に抽出された ATA や MBT、および BTAH などの有望分子に対しては、次ステップとしてDFTやMDを用いた『吸着シミュレーション』を実行することが推奨されます。界面での電荷密度差分や安定サイトを詳細に解析することで、電子論的メカニズムに基づいた確実な材料設計が可能となります。

まとめ#

本事例では、量子化学計算パッケージNWChemによる自動ハイスループット計算と、平面波・擬ポテンシャルに基づいた第一原理計算ソフトウェアAdvance/PHASEのESM法による高精度な表面電位評価を組み合わせ、防錆剤スクリーニングフローを実証しました。振動計算による安定性確認とハイブリッド汎関数による高精度解析を組み合わせることで、電子移動モデル（ $\Delta N$ ）の評価精度を向上させ、防錆剤のハイスループット・スクリーニングにおける本手法の有用性を示しました。それぞれのソフトウェアの強みを活かした本手法の連携は、実験コストを削減しつつ、メカニズムに基づいた新材料設計を加速させることが期待されます。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献#

I. Lukovits, E. Kálmán, and F. Zucchi, "Corrosion Inhibitors—Correlation between Electronic Structure and Efficiency", Corrosion 57, 3 (2001).
A. Kokalj, "On the HSAB based estimate of charge transfer between adsorbates and metal surfaces", Chem. Phys. 393, 1 (2012).
E. Aprà et al., "NWChem: Past, present, and future", J. Chem. Phy. 152, 184102 (2020).
ESM法による仕事関数の高精度計算
S. P. Ong, W. D. Richards, A. Jain, G. Hautier, M. Kocher, S. Cholia, D. Gunter, V. Chevrier, K. A. Persson, and G. Ceder, "Python Materials Genomics (pymatgen) : A Robust, Open-Source Python Library for Materials Analysis", Comp. Mater. Sci. 68, 314 (2013).
A. Kokalj, "Is the analysis of molecular electronic structure of corrosion inhibitors sufficient to predict the trend of their inhibition performance?", Electrochimica Acta 56, 745 (2010).
A. Kokalj, "Molecular modeling of organic corrosion inhibitors: Calculations, pitfalls, and conceptualization of molecule–surface bonding", Corrosion Science 193, 109650 (2021).