Pt-Ni電極触媒の第一原理解析：酸素還元反応（ORR）における触媒活性#

PHASE logo

燃料電池のカソード極で進行する酸素還元反応（ORR）において、現在最も優れた単金属の電極触媒は白金（Pt）ですが、純Ptは酸素種との結合が「わずかに強すぎる」ため、理論上の最高活性（Volcano Plotの頂点）には達していません。これを克服し、白金使用量の削減と高活性化を両立する電極触媒材料として、Ptと遷移金属（Ni, Coなど）の合金系が広く研究されています。特に、合金を高温で熱処理（アニーリング）した際、表面エネルギーを最小化するために最表面にPtが偏析（Segregation）して形成される「Pt-skin（白金スキン）」構造を持つPt-Ni合金は、純Ptを大幅に凌ぐORR活性を示すことが知られています。本解析では、第一原理計算ソフトウェアAdvance/PHASEを用い、Pt-skin/Pt₃Ni合金触媒表面における酸素吸着計算を実施しました。既存の純金属のVolcano Plot解析に合金系の結果を追加・拡張した結果、第一原理計算から得られる定量的な「吸着エネルギーの相対シフト」から、合金電極触媒が高い活性を示すことが明らかになりました。

Keywords: 第一原理計算 (DFT), 燃料電池, 電極触媒, 酸素還元反応 (ORR), Pt-Ni合金, Pt-skin, Volcano Plot, 歪み効果, スピン分極

計算モデルと計算条件#

Pt₃Niのような合金電極触媒において、最高レベルのORR活性を引き出す理想的な表面構造として「Pt-skin」が挙げられます。これは、合金内部と表面の原子配列が再構築され、最表面が純Pt層で覆われた理想的な構造です [1]。

本解析ではこの優れた高活性表面を模倣し、図1に示すように、面心立方（fcc）構造の(111)面をベースとした4層のスラブモデル（2x2スーパーセル）を構築しました。最表面（第1層）は100%のPt原子とし、第2〜4層をPt₃Ni組成とした「Pt-skin/Pt₃Ni(111)」モデルを採用しています。スラブ間の真空層は20 Å確保し、下部2層の原子座標を固定して構造最適化を行いました。

酸素原子の吸着サイトは、Pt-skin表面の高対称性サイト（トップ、ブリッジ、fccホロ、hcpホロ）でそれぞれ計算を行い、最安定配置を探索しました。

Pt-skin/Pt3Niモデル
図1. Pt-skin/Pt₃Ni(111) 表面スラブモデルの上面図と斜視図（最表面がPt、下層がPt₃Ni合金）。固定原子層を半透明に表示しています。

表1. 計算条件の概要

項目	設定
擬ポテンシャル	ウルトラソフト
交換相関汎関数	RPBE
スピン分極（磁性）	考慮あり (Ferromagnetic)
波動関数のカットオフエネルギー	30 Rydberg
表面モデル	4層スラブ（底面2層固定、2x2セル）

※ Niは強磁性金属であり、合金化後も磁気モーメントを持ちます。また、酸素原子などの不対電子とのスピン相互作用を正しく記述し、反応エネルギーを正確に見積もるため、本計算ではスピン分極（Spin polarization）を明示的に考慮しています。

計算結果と考察#

Pt₃Niバルクの格子定数最適化と歪み効果#

スラブモデルの構築に先立ち、RPBE汎関数を用いてPt₃Niバルクの格子定数を最適化しました。Ni原子の初期スピン分極は2.0に設定されています。図2にPt₃Niのエネルギー-体積（E-V）曲線のフィッティング結果を示します。

Pt3NiのE-V曲線
図2. Pt₃Ni の E-V曲線 (RPBE)。最適化された平衡格子定数は 3.94 Å でした（純PtのRPBE計算値は約 4.02 Å）。

RPBEで得られたPt₃Niバルクの格子定数（3.94 Å）は、純Ptの格子定数（4.02 Å）よりも約2%小さくなります。本解析ではスラブモデルの面内格子定数としてPt₃Niの最適化値（3.94 Å）を採用しました。これにより、最表面のPt層は純Pt本来の結合長よりも押し潰された状態となり、理論的に非常に重要な「圧縮歪み（Compressive strain）」がモデルに導入されます。

酸素吸着エネルギーとVolcano Plotの更新#

Pt-skin/Pt₃Ni(111)表面において酸素を吸着させた結果、純Pt(111)と同様に fcc hollow サイトが最安定であることが確認されました。この最安定構造を用いて、単金属触媒の解析事例 [2] と同様に、Nørskovらの定義 [3] に基づき、水（H₂O）と水素（H₂）を基準とした酸素吸着エネルギー（ $\Delta E_\text{O}$ ）を算出しました。

$\Delta E_\text{O} = E_{\text{O-slab}} + E_{\text{H2}} - E_\text{slab} - E_{\text{H2O}}$

表2. 純PtとPt-Ni合金の酸素吸着エネルギー（ $\Delta E_\text{O}$ ）の比較

電極触媒モデル	$\Delta E_\text{O}$ (eV)	特徴
純Pt(111)	1.54	吸着がわずかに強すぎる
Pt-skin/Pt₃Ni(111)	1.75	吸着が適度に弱まり最適化（+0.21 eVのシフト）

計算の結果、Pt-Ni合金表面における酸素吸着エネルギーは1.75 eVとなり、純Pt（1.54 eV）と比較して結合が約0.21 eV弱まることが定量的に示されました。この結果を、計算済の5種類の単金属からなるVolcano Plot [2] にプロットしたものが図3です。

Pt-Ni合金を含むORR Volcano Plot
図3. Pt-Ni合金の結果を追加したORR Volcano Plot。Pt-skin/Pt₃Ni（赤い星）の純Ptからのエネルギーシフト量（+0.21 eV）が、文献 [4] での理想的な頂点へのシフト量と非常によく一致していることが分かります。

純PtはVolcano Plotの頂点よりも左側（結合が強すぎる領域）に位置していましたが、Pt-skin構造による合金化効果によって吸着エネルギーが右側へシフトしました。今回得られた+0.21 eVというシフト量は、文献 [4] で提唱されている理想的なエネルギーシフト量 (約+0.2 eV)と非常によく一致しており、吸着エネルギーがVolcanoの頂点（理想的な触媒活性を示すピーク）付近へ最適化されたことを示しています。これは、Pt-skin/Pt₃Ni合金電極触媒が純Ptを凌駕する高いORR活性を持つことの理論的な妥当性を示すものであり、実験結果 [1] とも整合しています。

活性向上のメカニズム考察#

本計算で得られた酸素吸着エネルギーの最適化（低下）については、以下の2つの効果が同時に作用します。

歪み効果（Strain effect）： 下層のPt₃Niの小さな格子定数に合わせるため、表面のPt原子間距離が圧縮される効果。
配位子効果（Ligand effect）： 表面Pt原子が、直下に存在する異種金属（Ni）と電子的な相互作用を行う効果。

これらの効果により、酸素の2p軌道とPtの5d軌道が混成して形成される反結合性軌道（Anti-bonding states）がフェルミ準位より下に引き下げられ、電子に占有される割合が増加します。反結合性軌道への電子の充填はO-Pt間の結合を不安定化させるため、結果として酸素吸着エネルギーがマイルド（結合が弱く）になります [5]。

Volcano Plot解析の有効性について#

本解析では、真空中のスラブモデルを用いて計算を行っています。実際の電極触媒反応は水溶液（電解液）中で進行するため、溶媒和エネルギーや界面電場の影響が存在します。また、DFTの交換相関汎関数（RPBE等）自体にも絶対値としての計算誤差は含まれます。それにもかかわらず、なぜ真空モデルの第一原理計算が実験的な触媒活性のトレンド（Volcano Plot）[3] をこれほど正確に記述できるのでしょうか。

その最大の理由は「誤差相殺効果（Error Cancellation）」にあります。本解析のように、似たような遷移金属の(111)面において、同じ酸素原子（O*）の吸着を比較する場合、水分子との相互作用（溶媒和効果）や汎関数の系統的な過大・過小評価の度合いは、どの金属表面でもほぼ同程度に発生します。そのため、純PtとPt-Ni合金の相対的なエネルギー差を議論する際には、これらの系統誤差が引き算によって相殺（キャンセル）されるのです。これが、シンプルな計算モデルを用いたVolcano Plot解析が高い信頼性と強力な予測能力を持つ理由です。

まとめ#

本解析では、第一原理計算ソフトウェアAdvance/PHASEを用いて、Pt-skin/Pt₃Ni電極触媒のスラブモデルを対象として、スピン分極を考慮したDFT計算により、合金化による酸素吸着エネルギーの低下（1.54 eV → 1.75 eV）を定量的に算出しました。算出した吸着エネルギーをVolcano Plot上にマッピングすることで、Pt-Ni合金が理想的な頂点の近くに位置することを予測し、実験結果とも整合しています。この結果は、表面の圧縮歪みとNiの配位子効果に起因するものです。本手法のような吸着エネルギーの定量的評価は、Volcano Plot解析と組み合わせることで、高価な白金に代わる次世代の高性能電極触媒のスクリーニングや合理的な材料設計において強力なアプローチとなります。

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

お問い合わせ

参考文献#

V. R. Stamenkovic, B. Fowler, B. S. Mun, G. Wang, P. N. Ross, C. A. Lucas, and N. M. Marković, "Improved Oxygen Reduction Activity on Pt₃Ni(111) via Increased Surface Site Availability", Science 315, 493 (2007).
燃料電池の酸素還元反応における触媒活性の第一原理解析：Volcano Plot
J. K. Nørskov, J. Rossmeisl, A. Logadottir, L. Lindqvist, J. R. Kitchin, T. Bligaard, and H. Jónsson, "Origin of the overpotential for oxygen reduction at a fuel-cell cathode", J. Phys. Chem. B 108, 17886 (2004).
J. Greeley, I. E. L. Stephens, A. S. Bondarenko, T. P. Johansson, H. A. Hansen, T. F. Jaramillo, J. Rossmeisl, I. Chorkendorff, and J. K. Nørskov, "Alloys of platinum and early transition metals as oxygen reduction electrocatalysts", Nat. chem. 1, 552 (2009).
V. R. Stamenkovic, B. S. Mun, M. Arenz, K. J. J. Mayrhofer, C. A. Lucas, G. Wang, P. N. Ross, and N. M. Markovic, "Trends in electrocatalysis on extended and nanoscale Pt-bimetallic alloy surfaces", Nat. Mat. 6, 241 (2007).