固体相の第一原理計算データを用いたプルベー(Pourbaix)図の更新: Cu-Au-O-H系#

PHASE logo

材料の腐食防食設計において、水溶液中での金属の安定領域を示す「プルベー図（電位-pH図）」は欠かせないツールです。しかし、新規合金や複雑な金属間化合物の場合、実験データが存在しないことが多く、理論的な予測が強く求められています。本解析では、第一原理計算ソフトウェアAdvance/PHASEの高精度な固体相データと、Materials Project の広範なイオン・酸化物データを統合し、Cu-Au合金系における更新版プルベー図を作成しました。この手法により、任意の合金組成における腐食リスクの定量的評価が可能となります。

Keywords: 第一原理計算, DFTシミュレーション, プールベ図, 電位-pH図, 腐食予測, マテリアルズ・インフォマティクス, Materials Project, Cu-Au

プルベー図の原理と構築のアプローチ#

プルベー図は、電位（ $E$ ）とpHを軸として、水溶液中で最も安定な化学種（金属、イオン、酸化物など）を図示したものです。各境界線は、反応のギブス自由エネルギー変化（ $\Delta G$ ）がゼロになる条件から導かれるネルンストの式に基づいています。

プルベー図作成のベースデータ#

プルベー図を構築するためには、系に含まれる可能性のある全ての化学種のエネルギーデータが必要です。具体的には、以下の2つのカテゴリのデータが不可欠となります。

水溶液中のイオンのエネルギー： 金属が水中に溶出した状態（例：Cu²⁺, Au⁺ など）の自由エネルギー。水分子による複雑な水和構造（溶媒和）を考慮する必要があるため、通常の固体物理用バンド計算コードでの算出は高コストであり、一般的には信頼性の高い実験値や専用のデータベースを利用します。
固体相のエネルギー： 純金属、酸化物、水酸化物、そして今回対象とする合金（例：CuAu）のエネルギー。これらは第一原理計算によって高精度に求めることが可能な領域です。

プルベー図は、これらの「固体のエネルギー」と「イオンのエネルギー」を比較し、あるpHと電位環境下でどちらの状態がエネルギー的に低い（安定である）かを判定することで描かれます。

データ統合の課題と解決策#

しかし、第一原理計算を用いて独自のプルベー図を作成する場合、単純に異なるソースのデータを組み合わせるだけでは物理的な整合性が取れません。主な課題は以下の通りです。

計算と実験の混在： イオンの実験データ（またはMaterials Project [1] のデータベース値）と、手元の固体計算データを混在させる必要があります。
エネルギー基準の不整合： 異なる計算コード（例: Advance/PHASEとVASP）や設定間で、エネルギーの基準点（Reference）が異なるため、値を単純に数値比較することができません。

本解析では、これらの課題を解決するために「形成エネルギーブリッジ法」を採用しました。これは、Advance/PHASEで計算した「固体の形成エネルギー（ $\Delta H_\text{form}$ ）」を、Materials Project（MP）のエネルギー体系に適合させる手法です。

エネルギー統合のロジック:

MP互換の合金エネルギー $E_\text{compatible}$ は以下の式で算出します。

$E_\text{compatible} = \Delta H_\text{form}\text{(PHASE)} + \Sigma_i\left[ n_i E_\text{ref}\text{(MP)}\right]$

ここで、 $\Delta H_\text{form}\text{(PHASE)}$ はAdvance/PHASEで計算された形成エネルギー、 $E_\text{ref}\text{(MP)}$ はMaterials Projectのデータベースから抽出した単体（Cu, Au）の基準エネルギー、 $n_i$ は化学式中の各元素の原子数です。この手法により、第一原理計算データの精度を保ったまま、MPが持つ膨大な水溶液データと整合させることができます。

解析手順と計算データ#

1. 固体相の計算（Advance/PHASE）#

Cu-Au系の主要な金属間化合物（Cu₃Au, CuAu, CuAu₃）について、構造最適化およびエネルギー計算を行いました [2] 。

※PBEsol汎関数の採用について：

本解析では、交換相関汎関数としてPBEsolを採用しました。一般にMaterials Project等はPBE汎関数を用いていますが、Cu-Au系においては、PBEsolの方が実験的に観測される準安定相（特にCuAu₃など）の格子定数や相対安定性をより正確に再現できることが確認されています [2]。本手法では、固体の記述に適したPBEsol計算と、水溶液データの豊富なMPデータを適切に接続することで、実用的な予測精度を確保しています。

2. データの統合と更新（Python + pymatgen）#

オープンソースライブラリpymatgen [3] を使用し、以下の手順でデータを統合しました。

基準値の抽出： MPのPourbaixデータセットから、補正済みの単体エネルギー（Cu, Au）を取得し、基準点としました。
データの差し替え： MPデータベース内の既存のCu-Au合金データを、Advance/PHASEで計算した高精度データに置換しました。
環境データの結合： 水溶液中のイオン（Cu²⁺, Au⁺等）や酸化物（CuO, Au₂O₃等）は、MPのデータをそのまま利用しました。

3. プルベー図の作成#

Pythonライブラリ pymatgen [3] を使用し、以下の環境条件でプルベー図を作成しました。

温度： 298.15 K (25℃)
イオン濃度： $10^{-6}$ mol/L (一般的な腐食判定基準)

解析結果：Cu-Au-O-H系の更新版プルベー図#

電位-pH図 (Pourbaix Diagram)#

作成されたCu-Au系（Cu:Au=1:1組成）のプルベー図を図1に示します。

※図はクリックで拡大表示可能です。

図1.　Advance/PHASEの計算結果を用いて作成したCu-Au系 (Cu:Au=1:1) のプルベー図。計算条件は $T=298.15\text{ K}, [\text{Ion}]=10^{-6}\text{ M}$ 。中央下部の領域は固体合金（CuAu）が安定な領域（不感域）を示します。破線 (water lines)は水の安定領域を示します。

結果の解釈： 図中の「CuAu(s)」領域は、熱力学的に合金が安定な状態です。この領域外では、Cu成分の溶出（ $\text{Cu}^{2+}$ 等）や酸化物の形成が予測されます。

安定性ヒートマップ (Stability Heatmap)#

合金の組成によって耐食性は大きく変化します。ここでは、Auリッチな準安定相 CuAu₃ (Au: 75%) と、Cuリッチな安定相 Cu₃Au (Cu: 75%) について、それぞれの分解エネルギー（腐食への駆動力）を比較しました。

※用語の区別について：

「安定相・準安定相」という言葉は固体としての熱力学的な安定性（形成エネルギーの低さ）を指しますが、以下のヒートマップの色は水溶液中での腐食のしやすさ（分解エネルギーの高さ）を表しています。両者は必ずしも一致しない点にご注目ください。

準安定相 CuAu₃ (Au-rich)#

Stability Heatmap of Metastable CuAu3
図2.　準安定相 CuAu₃ の分解エネルギーヒートマップ。青色は分解エネルギーが低く、準安定状態を維持しやすい（耐食性が高い）領域を示します。

安定相 Cu₃Au (Cu-rich)#

Stability Heatmap of Cu3Au
図3.　安定相 Cu₃Auの分解エネルギーヒートマップ。図2と比較して、電位上昇に伴いより急速に赤色（高分解エネルギー）へ変化しており、腐食リスクが高いことを示す。

比較考察：

CuAu₃（図2）： 固体としては準安定相ですが、ヒートマップ上では広いpH・電位範囲で青色の領域（Low Driving Force）が確認されます。Au含有率が高いため、水溶液中での分解への駆動力が小さく、相対的に高い耐食性を示すことがわかります。
Cu₃Au（図3）： 固体としては熱力学的に安定相ですが、図2と比較すると、より低い電位（0V付近）から色が黄色〜赤色に変化しています。これはCu成分が多いため、わずかな電位上昇で $\text{Cu}^{2+}$ として溶出しようとする駆動力が強く働くことを示しています。すなわち、「固体として安定であること」と「水溶液中で腐食しないこと」は同義ではないことが、定量的に可視化されています。

補足：気相分子を含む系の注意点#

本事例のCu-Au系は金属間化合物ですが、酸化物などを計算対象とする場合は注意が必要です。一般にDFT計算では、O₂などの気相分子のエネルギー計算に誤差が含まれることが知られています。そのため、酸化物の形成エネルギーを算出する際は、Materials Project等の推奨する補正スキーム（O correctionなど）を適用するか、あるいは固体の酸化物基準で化学ポテンシャルを決定するなどの工夫が必要です [4]。本手法は、それらの補正を行った後の形成エネルギー値さえあれば、同様に適用可能です。

まとめ#

本解析では、第一原理計算ソフトウェアAdvance/PHASEによるDFT計算データと、Materials Projectのデータを連携させ、Cu-Au系の腐食特性を評価しました。単なる安定領域の図示にとどまらず、分解エネルギーのヒートマップ解析を行うことで、準安定相であっても腐食環境下での安定性が期待できる領域や、組成による耐食性の違い（CuAu₃ vs Cu₃Au）を明確に示すことができました。この手法は、実験データが乏しい新規合金のスクリーニングにおいて極めて有効です。

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

お問い合わせ

参考文献#

A. Jain et al., "Commentary: The Materials Project: A materials genome approach to accelerating materials innovation", APL Materials 1, 011002 (2013).
Cu-Au二元系金属間化合物のConvex Hullと熱力学的安定性の第一原理計算
S. P. Ong et al., "Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis", Comput. Mater. Sci. 68, 314 (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).