コンテンツにスキップ

ステンレス鋼の耐食性（錆びにくさ）：第一原理計算による不動態化メカニズムの解明#

Advance/PHASE Logo

ステンレス鋼は、鉄（Fe）を主成分としながらもクロム（Cr）を添加することで、過酷な環境下でも錆びない（腐食しない）優れた特性を持っています。この耐食性の起源は、表面に形成されるナノメートルオーダーの緻密な保護膜「不動態皮膜」にあります。ステンレス鋼は長年の実験的知見に基づき広く実用化されていますが、その優れた耐食性の背後にある「なぜCrが主役となるのか」という根本的なメカニズムを原子・電子の振る舞いから理論的に裏付けることは、より高性能な次世代の耐食合金を論理的に設計するための強力な足がかりとなります。本事例では、第一原理計算ソフトウェア Advance/PHASE を用い、bcc鉄(001)表面におけるCr原子の挙動と酸素との相互作用を解析しました。計算コストを抑えつつ、エネルギー論、吸着状態、表面磁性の観点から、ステンレス鋼の不動態化の根源を理論的に裏付けました。

Keywords: 第一原理計算, DFTシミュレーション, ステンレス鋼, 耐食性, 不動態化, 表面偏析, 吸着エネルギー, スピン分極, 磁性状態

背景・計算方法：単純化モデルによる現象の本質抽出#

鉄の錆（腐食）は、水分と酸素が存在する環境下でFe原子が電子を失い、酸化物へと変化するプロセスです。ステンレス鋼では、添加されたCrがFeよりも優先的に酸素と結合し、緻密で剥がれにくいCr2O3を主体とする不動態皮膜を形成することで、内部のFeを保護します [1]。この「なぜCrが？」という根本的な疑問に対し、原子レベルの相互作用からアプローチします。

計算にはbcc鉄（フェライト系）の(001)表面を模したスラブモデルを採用しました。具体的には、Fe(001) 2x2スーパーセル（表面の局所的なCr濃度25 at%に相当）を用い、定性的なエネルギー差の比較と計算コストのバランスから5原子層としています。ウルトラソフト型擬ポテンシャルを採用し、GGA-PBE汎関数を用いました。格子定数は最適化された2.8424 Åとし、波動関数のカットオフエネルギーは25 Rydberg、k点サンプリングは6x6x1としています。FeとCrの磁性を考慮するためスピン分極計算を行い、底面2層をバルク格子に固定することで内部環境を模擬し、上面3層と吸着原子を構造最適化（緩和）させました。このモデルを用いることで、計算コストを抑えつつ、複雑な腐食現象の本質が抽出可能となります。

1. エネルギー論：Crの表面偏析と保護膜形成の駆動力#

まず、ステンレス鋼内部のCrがどのように表面へ移動し（表面偏析）、なぜFeよりも強固な皮膜を作る起点となるのかをエネルギー論の観点から評価しました。

Cr in bulk or surface

図1. 左：Crがバルク内部にある状態、右：Crが表面に移動し酸素と結合した状態。固定層（底面2層）の原子が半透明で表示されています。

【結果】Crの表面偏析：強力な熱力学的駆動力#

図1に示すように、酸素が吸着した環境において、Cr原子が鉄の内部（第3層）にある状態と、表面層に移動した状態の全エネルギーを比較しました。
その結果、Crが表面層にあるモデルの方がエネルギーが約0.66 eV低く、安定であることが示されました。これは、「酸素が存在する環境下では、Crは鉄の内部にとどまるよりも、自ら表面に移動して酸素と結びつく方が圧倒的に安定である」という熱力学的な駆動力を明らかにしています。第一原理計算を用いた先行研究では、酸素が存在しない真空中においてはCrの表面偏析は起こりにくいことが報告されていますが [2]、本解析により、「酸素の存在」こそがCrを表面に引き出し、自己修復的な不動態皮膜を形成する真のトリガーであることが理論的に裏付けられました。

【結果】酸素吸着エネルギー：Feより強固なCr-O結合#

次に、図2に示す純鉄モデルと表面Crモデルに対して、酸素原子の吸着エネルギー（）を算出しました。基準として、O2分子ではなく、計算誤差の少ない水分子（H2O）と水素分子（H2）のエネルギー差を用いています（仮想的な反応 ）。

  • 純鉄モデル: -0.81 eV
  • 表面Crモデル: -0.96 eV

oxygen adsorption on Fe(001) surface

図2. 酸素吸着構造。左：純鉄表面モデル、右：表面Cr置換モデル。Cr-O結合の方がエネルギー的に有利です。

計算結果は、Crが表面にある場合の方が約0.15 eVも強く酸素を吸着（結合）していることを示しています。「CrはFeよりも酸素との親和性が高く、より強固な結合を作る」という化学的性質が、原子レベルでも如実に再現されました。この強固なCr-O結合が、緻密で剥がれにくい不動態皮膜形成の初期過程（起点）を駆動します。

吸着エネルギーの基準の選び方

一般的に、酸素の吸着エネルギーを考える際、気相の酸素分子（O2）を基準としがちです。しかし、DFT計算（特にGGA汎関数）において、O2分子の結合エネルギーを過大評価してしまう傾向（過剰結合誤差）があります。そのため、現実の腐食環境（水溶液環境）に近く、かつPBE汎関数でも高精度に計算できる水分子（H2O）と水素分子（H2）を基準とすることで、この誤差を回避し、より信頼性の高い吸着エネルギーの比較が可能となります。

2. 表面磁性：反平行スピンと電子的な相互作用#

エネルギー論的な強固さに加え、FeとCrが混在する表面でのユニークな電子状態（磁性）についても解析しました。FeとCrはいずれも磁性原子ですが、Fe母相中にCrが存在する場合、Crの磁気モーメントはFeに対して「反平行（逆向き）」に向く性質があります。Advance/PHASEでの計算設定においても、この反平行スピンを初期状態とすることで、安定な基底状態へと収束させました。

Surface Magnetism State with Cr Segregation

図3. Cr偏析による表面磁性状態（スピン密度の断面図）。黄色：Fe原子領域（上向きスピン）、青色：Cr原子領域（下向きスピン）。視認性のため、Fe(001)4x4スーパーセルで表示しています。

【結果と考察】反平行スピンが反映する電子的な再配置#

図3は、表面付近の磁気モーメント（スピン）の分布を可視化したものです（色の違いがスピンの向きと強さに対応）。Fe原子領域（黄）とCr原子領域（青）が、互いに逆向きの磁気モーメントを持っている様子が鮮明に見て取れます。
このCrの反平行スピン配置は、Crのd軌道電子がFeのd軌道電子と複雑に相互作用（混成）している結果であり、単に原子が並んでいるだけではない、電子レベルでの強固な結合状態を示唆しています。この磁性状態を正確に捉えることは、Fe-Cr合金系の電子構造設計において不可欠な視点です。

まとめ#

第一原理計算ソフトウェアAdvance/PHASEを用いた本解析により、bcc鉄(001)表面におけるCr原子の挙動と酸素吸着メカニズムが、原子・電子レベルで明らかになりました。エネルギー論的な解析の結果、酸素存在下においてCrが表面へ自発的に移動（表面偏析）する熱力学的駆動力（約0.66 eV）が生じること、ならびにFeよりも強固なCr-O結合（Feより約0.15 eV 安定）を形成することが定量的に評価されました。これらは、ステンレス鋼が緻密で安定な不動態皮膜を自律的に形成し、錆を防ぐメカニズムの根幹をなすと考えられます。さらに表面磁性の解析から、Cr原子がFe原子に対してスピンを反平行に向けている状態が確認され、複雑なd軌道電子間の混成が示唆されました。

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

お問い合わせ

参考文献#

  1. C.-O.A. Olsson, and D. Landolt, "Passive films on stainless steels—chemistry, structure and growth", Electrochim. Acta 48, 1093 (2003).
  2. W. T. Geng, "Cr segregation at the Fe-Cr surface: A first-principles GGA investigation", Phys. Rev. B 68, 233402 (2003).

関連ページ#