鉄－水系のプルベー（Pourbaix）図の更新：磁性体のDFT+U計算値と実験値の併用#

PHASE logo

鉄鋼材料の腐食防食設計において、水溶液中での熱力学的安定領域を示す「プルベー図（電位-pH図）」は不可欠なツールです。しかし、既存のデータベースに基づく自動生成図は、鉄のように複雑な磁性や準安定相を持つ系において、そのまま利用すると実験事実（実際の錆の成分など）と異なる結果を示す場合があります。本解析では、既存のPourbaixデータをベースとしつつ、第一原理計算ソフトウェアAdvance/PHASEによる高精度なDFT+U計算値と、信頼性の高い実験値をハイブリッドに適用してデータを更新しました。これにより、実験的に観測される腐食挙動と整合するFe-O-H系プルベー図が得られました。

Keywords: 第一原理計算, DFT+U, 強相関電子系, プールベ図 (Pourbaix Diagram), 腐食防食, 鉄鋼, Fe3O4, Magnetite, Materials Project

解析のアプローチ：データベースと独自計算の統合#

Materials Project (MP) [1] では、鉄－水 (Fe-O-H)系の豊富なデータが提供されております。また、MPで採用されているPerssonらのCorrection Scheme [2] は、固体のDFT計算結果を水溶液化学に接続し、実験値と整合させるための優れた補正手法です。そこで本事例では、MPから基本データをダウンロードした後、以下の3段階の更新処理（Update Scheme）を適用することで、プルベー図の更新を図りました。

データの更新と詳細#

1. Fe₃O₄ (Magnetite) の更新：磁性を考慮したDFT+U計算#

鉄の酸化物は強い電子相関を持つため、エネルギー評価に誤差が生じる場合があります。特にマグネタイト（Fe₃O₄）については、Advance/PHASEを用いて以下の条件で形成エネルギー（ $\Delta H$ ）を再計算し、その値を適用しました。

計算手法: GGA(PBE)+ U
U値の設定: $U_\text{eff} = 4.0 \text{ eV}$ (Fe-3d軌道に対して適用 [3])
磁気構造: Ferrimagnetic（フェリ磁性）
初期磁気モーメント設定: Fe(A-site) = -4.0 $\mu_B$ , Fe(B-site) = +4.0 $\mu_B$

DFT+Uの計算より、Fe₃O₄の磁気モーメントは 4.0 $\mu_B$ /formulaという結果が得られました。算出された形成エネルギーに対しては、Wangら [3] の手法に基づく酸素分子（O₂）のエネルギー補正を行い、MPのデータセット内のFe₃O₄の値を更新しました。

Fe₃O₄の形成エネルギーのDFT+U計算値: $-1.69$ eV/atom

2. Fe(OH)₂ の更新：実験的自由エネルギーの採用#

水酸化第一鉄（Fe(OH)₂）は、水和や溶解度積の影響を強く受けるため、DFT計算よりも実験値の方が信頼性が高い場合があります。本解析では、腐食科学分野で広く参照される実験値（標準生成ギブス自由エネルギー $\Delta G^\circ_f$ , $-486.5$ ～ $-492$ kJ/mol）[4, 5]の平均値を採用しました。

採用データ: $\Delta G^\circ_f (\text{Fe(OH)}_2) \approx -489 \text{ kJ/mol}$

この値を採用することで、高pH領域におけるFe²⁺イオンからの沈殿生成ラインを実験事実に合わせ込みました。

3. FeO (Wustite) の除外#

FeOは570℃以上で安定な高温相であり、常温（25℃）の水溶液環境下では熱力学的に不安定（FeとFe₃O₄に不均化）です [4]。物理化学的な妥当性に基づき、計算対象リストからFeOを明示的に除外するフィルタリング処理を行いました。

解析結果：更新されたFe-O-H系プルベー図#

上記の手順でデータを更新・再構築した結果得られたプルベー図（図1）を示します。

Updated Fe-O-H Pourbaix Diagram Zoomed Fe-O-H Pourbaix Diagram
図1. （上）本解析によるFe-O-H系のプルベー図。 $T=298.15\text{ K}$ , $[\text{Ion}]=10^{-6}\text{ M}$ 。赤破線は水の安定領域の上限・下限（酸素・水素発生線）を示します。（下）Fe-O-H系のプルベー図の一部拡大図。

実験事実との整合性#

図1の結果は、鉄の腐食における固体相の出現順番において、標準的な実験図 [4, 5]と良好な一致を示しています。

Fe₃O₄領域の出現: 図1の（pH 12, -0.72V）を中心とした範囲に安定相として現れており、鉄の防食に寄与する黒錆（不動態皮膜）の存在を予測できています。なお、Fe₃O₄の安定領域の幅は一般的な実験図と比較して狭くなっています。これは、Fe₂O₃（赤錆）とのエネルギー差が非常に僅差であり、計算条件や実験における実際の生成物の状態（粒子サイズ効果など）に敏感に依存するためと考えられますが [6]、腐食生成物の熱力学的な出現順序は正しく再現されています。
Fe(OH)₂領域の出現: Fe₃O₄領域の下側、より低電位側において、適切な溶解度バランスで水酸化物が生成することが確認されました。

この結果は、「データベースの値をそのまま使う」だけでなく、「DFT計算や実験知見を組み合わせてデータを更新する」というハイブリッドなアプローチが、実用的な腐食予測において有効であることを示しています。

まとめ#

本事例では、Materials Projectのデータをベースにしつつ、Advance/PHASEによる磁性体DFT+U計算値と実験値を適切に組み合わせてデータを更新することで、実験的に観測される腐食挙動と整合するFe-O-H系プルベー図を構築しました。この手法を用いることで、既存データベースの不足部分を補い、実際の腐食環境に即した信頼性の高い相図を得ることが可能です。

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

お問い合わせ

参考文献#

A. Jain et al., "Commentary: The Materials Project: A materials genome approach to accelerating materials innovation", APL Materials 1, 011002 (2013).
K. A. Persson, B. Waldwick, P. Lazic, and G. Ceder, "Prediction of solid-aqueous equilibria: Scheme to combine first-principles calculations of solids with experimental aqueous states", Phys. Rev. B 85, 235438 (2012).
L. Wang, T. Maxisch, and G. Ceder, "Oxidation energies of transition metal oxides within the GGA+U framework", Phys. Rev. B 73, 195107 (2006).
B. Beverskog and I. Puigdomenech, "Revised Pourbaix diagrams for iron at 25-300 °C", Corros. Sci. 38, 2121 (1996).
M. Pourbaix, "Atlas of Electrochemical Equilibria in Aqueous Solutions", NACE, Houston (1974).
T. Misawa, "The thermodynamic consideration for Fe-H₂O system at 25°C", Corros. Sci. 13, 659 (1973).