クラスター展開法を用いたα黄銅の規則・不規則相転移シミュレーション#

Advance/PHASE Logo

Cu-Zn合金（黄銅）は、優れた加工性と電気・熱伝導性を持ち、広く利用されている実用合金です。特に電子濃度によって相が決定されるヒューム・ロザリー合金の代表例として知られています。 $\beta$ 黄銅の規則-不規則相転移は高温で生じるため有名ですが、Cuリッチな $\alpha$ 黄銅（Zn 25%付近、FCC構造）においても、低温で $L1_2$ 構造などの長距離規則相（ $\alpha^\prime$ 相）へ転移することが第一原理計算から理論的に予測されています。しかし、転移温度が300 K以下と低いため原子拡散が遅く、特殊な触媒還元プロセスによるナノ粒子としての報告を除き、バルク状態での実験的な観測は困難とされています。本事例では、第一原理計算ソフトウェアAdvance/PHASEと合金熱力学ツールATATを連携させ、絶対零度での電子状態計算から、有限温度でのCuリッチ側のマクロな相転移挙動（等温線）およびミクロな原子配列の可視化までを一貫して解析したマルチスケール・シミュレーションを行います。

Keywords: 第一原理計算 (DFT), クラスター展開法 (CE), モンテカルロ法 (MC), ATAT, Cu-Zn合金, α黄銅, 相転移

理論的背景：マルチスケール合金熱力学#

クラスター展開法 (Cluster Expansion Method)#

有限温度における合金の相状態を予測するには、膨大な原子配置パターンに対するエネルギーを評価する必要があります。すべての配置に対して第一原理計算を行うことは計算コストの観点から不可能です。そこで、少数の代表的な構造に対する第一原理計算の結果から、任意の原子配置のエネルギーを高精度に予測する手法がクラスター展開法 [1] です。

クラスター展開によるエネルギーの定式化:
合金の任意の配置 $\sigma$ のエネルギー $E(\sigma)$ は、有効相互作用パラメータ（ECI）$V_c$ とクラスター相関関数 $\Phi_c(\sigma)$ の和として以下のように展開されます。 $$E(\sigma) = \sum_{c} m_c V_c \Phi_c(\sigma)$$

ここで、$c$ はペア（対）やトリプレット（3体）などのクラスターの種類、$m_c$ はその多重度を表します。DFTで計算された数十〜数百の構造エネルギーから、ATAT [2] を用いて最適化された $V_c$ を決定します。

モンテカルロ法による有限温度シミュレーション#

構築されたクラスター展開モデルを用いることで、巨大なスーパーセル（数千〜数万原子）のエネルギーを高速に計算可能になります。これにメトロポリス法に基づくモンテカルロ（MC）法を適用し、温度（ $T$ ）や化学ポテンシャル（ $\mu$ ）を変化させることで、相転移温度の特定やヒステリシス挙動の観測といったマクロな熱力学特性を導き出します。

計算手順#

本事例では、ATATによる候補構造の生成から、Advance/PHASEのインプットファイル作成・ジョブ投入・結果解析、そしてクラスター展開モデル（ECI）の構築と精度評価に至る一連の計算サイクルを自動化しており、複雑な合金系の相空間探索を効率的に実行しています。

第一原理計算: 自動化ワークフローに沿って、純Cuから純Znまでの多様なFCC構造（約100種類以上）を生成し、PBEsol (GGA)を用いて構造最適化および全エネルギー計算を実施します。Cu-Zn系はCuリッチ側でFCC構造（ $\alpha$ 相）、Cu:Zn=1:1付近でBCC構造（ $\beta$ 相または $\beta^\prime$ 相）、Znリッチ側でHCP構造（ $\eta$ 相）をとることが知られていますが [3]、本事例ではCuリッチ側に焦点を当てているため、FCCベースでの探索を行っています。
モデル構築: ATATを用いて自動抽出された計算結果を随時学習し、交差検証（Cross-Validation）スコアが最小となる高精度なクラスター展開モデル（ECI）を構築します。
等温線（Isotherm）の計算: グランドカノニカルアンサンブルにて化学ポテンシャル $\mu$ をスイープし、300 Kおよび100 KにおけるZn濃度 $x_{Zn}$ の変化（ヒステリシスとプラトー）を解析します。
ミクロ構造の可視化: カノニカルアンサンブル（濃度固定）のMCシミュレーションで得られた不規則相と、DFT計算で決定された基底状態である規則相の原子配列を比較します。

計算結果と考察#

1. 高精度なクラスター展開モデルの構築#

100種類以上の構造データを用いてフィッティングを行った結果、CVスコアは 0.017 eV/atom (17 meV/atom) となりました。この予測誤差は100 Kの熱エネルギー（約 8.6 meV）を上回っているため、相転移温度の定量的予測には不確実性が残りますが、室温以下の低温域において、規則・不規則相転移が駆動される物理的メカニズムや、その定性的な傾向を議論するための基盤として十分に機能します。図1(a)の生成熱（Convex Hull）を見ると、 $x_{Zn} \approx 0.6$ 付近（ $\gamma$ 相領域）に系全体の基底状態（最安定状態）が存在する一方で、本解析で着目する $\alpha$ 相領域の $x_{Zn} = 0.25$ も凸包（Convex Hull）上の安定点として同定されており、Cu₃Zn 規則相の安定性が明確に表現されています。また、図1(b)の残差プロットでは大半の構造が $\pm 0.02$ eV以内の誤差に収まり、図1(c)のECIの距離依存性も物理的に妥当な減衰を示しています。

Convex Hull of Cu-Zn system
図1(a) 生成熱と凸包 (Convex Hull)

Residual error of the fit
図1(b) DFTエネルギーとモデル予測値の残差

ECI vs. cluster radius
図1(c) 有効相互作用パラメータ (ECI) のクラスター半径依存性

2. マクロな熱力学挙動：等温線（Isotherm）の温度依存性#

得られたモデルを用いて、化学ポテンシャルを変化させた際のZn濃度の応答（等温線）を計算しました。図2(a)の300 Kでは、化学ポテンシャルを増加させた場合（Zn挿入）と減少させた場合（Zn除去）の経路が完全に一致し、ヒステリシスを持たない滑らかな曲線を描いています。これは系が不規則固溶体（ランダムな $\alpha$ 相）であることを示しています。

一方、図2(b)の100 Kでは、300 Kの滑らかな曲線とは異なり、ステップ状の細かい変化とわずかなヒステリシスが現れました。本計算では、100 Kという極低温における原子拡散の遅れ（速度論的なトラップ）を考慮して、モンテカルロ計算の緩和ステップ数を30,000 stepsに設定してシミュレーションを行いました。その結果、理想的な $x_{Zn} = 0.25$ （Cu₃Zn規則相）での完全な平坦部（プラトー）を形成するには至っていませんが、300 Kと比較して明らかな等温線の形状変化（ステップの出現）が観測されたことは、室温以下の低温域において、系が不規則相から規則相へと向かう強い熱力学的な駆動力を持っていることを定性的に裏付けています。

Isotherm at 300K
図2(a) 300 K における等温線（不規則相）

Isotherm at 100K
図2(b) 100 K における等温線（規則相の形成）

ミクロな原子配列の変化として、DFTで決定された基底状態であるCu₃Znの規則相と、300 K（濃度固定）のモンテカルロ計算から得られた不規則相（ランダム固溶体）の原子配列スナップショットを図3に示します。本シミュレーションにより、Cu₃Znの規則-不規則転移温度は100 K〜300 Kの間にあることが示唆されました。これはTurchiらが予測した「 $\alpha^\prime$ 相の安定領域は300 K以下」という結果と定性的に一致します [4]。特殊な触媒還元プロセスによるナノ粒子としての報告 [5] を除き、実験的にこの規則相がバルクで観測されないのは、室温以下の低温では原子の熱振動エネルギーが活性化障壁を下回り、拡散（並び替え）が事実上停止してしまうためです。本解析は、実験では到達困難な「熱力学的な真の平衡状態」を計算機上で明らかにする第一原理合金熱力学の強力さを示しています。

Disordered and Ordered phases of Cu3Zn
図3. Cu₃Znの原子配列比較。左：規則相、右：不規則相（ランダム固溶体）。

まとめ#

本事例では、第一原理計算ソフトウェアAdvance/PHASEによる高精度なDFT計算データをATATに連携させることで、Cu-Zn合金のクラスター展開モデルを構築しました。モンテカルロシミュレーションを通じて等温線におけるヒステリシスやステップ状の変化（規則相形成の兆候）を確認し、規則相・不規則相の相転移温度が300 K以下に存在することを特定しました。これは既存の知見と整合しています。この手法は、未知の金属間化合物の探索、状態図の作成、および新規合金の材料設計に有効なアプローチとなります。

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

お問い合わせ

参考文献#

A. van de Walle and G. Ceder, "Automating first-principles phase diagram calculations", J. Phase Equilib. 23, 348 (2002).
A. van de Walle, M. Asta, and G. Ceder, "The alloy theoretic automated toolkit: A user guide," CALPHAD 26, 539 (2002).
C. M. Eastman, Q. Zhang, and J.-C. Zhao, "Diffusion Coefficients and Phase Equilibria of the Cu-Zn Binary System Studied Using Diffusion Couples", J. Phase Equilib. Diffus. 41, 642 (2020).
P. E. A. Turchi, M. Sluiter, F. J. Pinski, D. D. Johnson, D. M. Nicholson, G. M. Stocks, and J. B. Staunton, "First-Principles Study of Phase Stability in Cu-Zn Substitutional Alloys", Phys. Rev. Lett. 67, 1779 (1991).
S. Derrouiche, H. Lauron-Pernot, and C. Louis, "Synthesis and Treatment Parameters for Controlling Metal Particle Size and Composition in Cu/ZnO Materials—First Evidence of Cu₃Zn Alloy Formation", Chem. Mater. 24, 2282 (2012).