分子動力学によるセリウムの超精密ダイヤモンド切削の解明

概要

構造相転移と転位ダイナミクスの結合が、セリウムのナノスケール変形の理解を複雑にしている。本研究では、分子動力学（MD）シミュレーションを用いて、超精密ダイヤモンド切削時のセリウムの変形メカニズムを解明した。経験的に導かれたCe–CeおよびCe–Cポテンシャルを組み込んだMDモデルを構築し、2つのfccセリウム相の弾性応答を評価した。シミュレーションの結果、塑性変形は転位すべりが支配的であり、径方向分布関数（RDF）解析では、加工表面と切りくずにおいてγ→δ相転移はわずかしか観察されなかった。パラメータスタディにより、高品質な加工面を得るための最適な切削条件が特定された。

背景

セリウム（Ce、原子番号58）は、独自の機械的・物理的・化学的特性を持つランタノイド金属として重宝されている。表面形態は、セリウム部品の性能、耐食性、水素貯蔵挙動に決定的な影響を与える[1–3]。超精密ダイヤモンド切削は、亜表面損傷を最小限に抑えた超平滑面を実現できる[6,7]が、セリウム切削に関する実験的・理論的研究はほとんど存在しない。超精密切削では工具刃先半径が切込み深さと同程度となるため、材料特性がプロセスを支配し、セリウム切削メカニズムは特に困難である。

Ceは4f電子の非局在化によって駆動される豊富な圧力-温度相図を示す。常圧かつ110 K未満ではα-Ce（fcc）が安定し、45–275 Kでβ-Ce（dhcp）、270–999 Kでγ-Ce（fcc）、999–1071 Kでδ-Ce（bcc）へと転移する[8–11]。特筆すべきは、295 K、8 kbarでのγ→α転移で体積が約20%収縮すること[12–14]、またγ→δ転移は高温高圧で誘発されることである。これらの相転移は結合と弾性特性を変化させ、変形挙動に影響を与える。金属Ceは転位活動によって支配される延性材料でもある[15]。しかし、ダイヤモンド切削時の相転移と転位の相互作用はほぼ未解明のままである。

MDシミュレーションは、原子スケールでの基礎的な加工メカニズムを探る強力な手法である。先行研究では、Cu、Al、Fe、高エントロピー合金のナノスクラッチや直交切削が調べられ、転位核生成、表面パイルアップ、摩擦挙動への洞察が得られている[16–27]。これらの研究は、結晶配向、工具すくい角、相互作用ポテンシャルの重要性を浮き彫りにした。しかし、セリウムの機械加工をMDで調べた研究はなく、本研究の動機となっている。

手法

ダイヤモンド切削のMDモデル

図1にシミュレーション設定を示す：単結晶Ceワーク（41 × 25 × 31 nm、約1 × 10^6個のγ-Ce原子）と、逃げ角9°のダイヤモンド切削工具（0.1 × 10^6個のC原子）。ワークは固定された2 nmの底層と、1 fsタイムステップでニュートン力学に従う移動層に分割。周期境界条件は縦方向のみ適用。3つの結晶配向—(010)、(110)、(111)—を検討。7つのすくい角（−30°、−20°、−10°、0°、10°、20°、30°）をテスト。ダイヤモンド工具は剛体として扱うため、C–C相互作用は省略。

原子間相互作用は以下のようにモデル化：Ce–Ce力にはShengら[29]の埋め込み原子法（EAM）パラメータ化を使用し、fcc Ceの弾性定数を正確に再現。Ce–C相互作用にはモースポテンシャルを採用：

$$E_{\text{tot}}=\sum_{ij}D_0\left[e^{-2\alpha(r-r_0)}-2e^{-\alpha(r-r_0)}\right]$$

パラメータ：D₀=0.087 eV、α=5.14、r₀=2.93 Å、カットオフ1.0 nm[30]。

系はまず30 K、0 barでNPTアンサンブルで平衡化し、次に一定速度100 m/s、切込み深さ4 nmでNVTアンサンブル下でダイヤモンド切削を実行。切削力は切削方向に沿って分解（図中の赤矢印が切削方向）。シミュレーション寸法と切削速度は実験的な超精密条件より桁違いに小さいが、本MDモデルはプロセスの本質的な物理を捉えている。共通近隣解析（CNA）でfcc、hcp、bcc、欠陥原子を識別[31]。LAMMPS（1 fsステップ）でシミュレーション実行、OVITOで可視化[32,33]。

セリウム相の特性評価

γ、α、β、ε、δの5つのセリウム相を文献値[8–11]の格子パラメータでモデル化。バルク構造を単軸引張、せん断、圧縮し、弾性定数と機械的特性を抽出（表2）。計算されたγ-Ceのヤング率（24.17 GPa）は実験的ナノインデンテーションデータ（36.7 GPa）[10]と整合し、弾性定数は非弾性中性子散乱結果[34]と一致。より高密度なα-Ceはγ-Ceより著しく高い剛性を示す。

径方向分布関数（RDF）解析で相間の構造差を特性化。図2aに各相の明確なRDFピークを示す（最近接距離に対応：γ-Ce 3.64、5.13、6.30 Å；α-Ce 3.41、4.85、5.92 Å）。γ-Ceを均一圧縮して20%体積収縮させると、α-Ceと一致するRDFが得られる（図2b）、これによりγ→α転移を確認[12–14]。

結果と考察

セリウムの加工メカニズム

0°すくい角でのCe(010)ダイヤモンド切削MDシミュレーションでは、3つの力のレジームが明らかになった（図3）。初期には接着力により負の切削力と法線力が発生；弾性変形が続き、力が急上昇。切削長さ約2.3 nmで急激な低下が見られ、塑性の開始を示す。これは自由表面近くで核生成した1/6<112>ショッキ部分転位が{111}面で滑ることで支配される。切削力と法線力は10–35 nmで安定し、工具前方と下方に転位が蓄積。35 nm以降、工具がワークから離れるため力が低下。全域で法線力は切削力より小さい。

図4は主要な切削長さにおける欠陥構造と表面形態の進化を示す。ショッキ部分転位は転位コアで境界づけられたスタッキングフォルトを生成。切削が進むと、前面から部分転位の密雲が発生し、すくい面に沿って切りくずを形成（図4f）。工具後方の転位は上面へ上昇し、表面回復を促進。工具が左境界に近づくと転位伝播が阻害され、切りくず体積が増大（図4g）。切りくず分離後、自由表面での消滅により転位密度が減少（図4h）。

欠陥ゾーンのRDF解析では、限定的なγ→δ相転移が確認された。表面下の欠陥領域のRDFはδ-Ceピークと一致し、局所加熱によるわずかなγ→δ転移を示唆（図5a）。切りくずでも同様のδ-Ceシグネチャが現れる（図5b）が、全体量は微々たるもので、相転移が切削中に支配的でないことを示す。

すくい角の影響

図6に示すように、切削力と法線力はすくい角の増大とともに減少し、マーチャントのせん断面エネルギー論[37]と一致。大きなすくい角は力の差を減らし、加工安定性を向上させる。

欠陥解析（図7a,b）では、負のすくい角（−30°）が正の角度（30°）より高い転位密度を示す。転位セグメント数（図7c）はすくい角とともに単調減少し、塑性変形の複雑さ低下と加工硬化の低減を反映。

表面形態（図8）もこれを裏付ける：負のすくい角では非対称で顕著な表面パイルアップが生じる一方、30°すくい角では最小のパイルアップ体積と対称的なプロファイルを実現。したがって、30°すくい角が低切削力、最小転位密度、最高の表面仕上げを得る最適条件である。

結晶配向の影響

最適な30°すくい角下で、Ce(010)、Ce(110)、Ce(111)配向を比較（図9）。Ce(010)が最も低い切削力と転位数を示す。Ce(111)は転位数が最少だが、{111}面が表面と平行なため表面パイルアップを促進し、切削抵抗を増大させる。加工面形態（図10）では、Ce(010)が最小パイルアップ体積と最も対称的なプロファイル、Ce(110)が最大パイルアップを示す。

結論

EAMとモースポテンシャルを組み込んだMDフレームワークを開発し、セリウムの超精密ダイヤモンド切削をモデル化した。Ce相の弾性定数、機械的特性、相転移挙動を検証し、モデルの予測能力を確認。シミュレーションにより、転位核生成と滑りが塑性を支配し、γ→δ相転移はわずかしか観察されないことを明らかにした。最適加工条件—30°すくい角と(010)結晶配向—が、最低切削力、最小転位密度、最高表面品質をもたらす。