MDで観るトライボロジー〜原子が擦れるその場をシミュレーション〜

[スポンサーリンク]

化学者のつぶやき

MDで観るトライボロジー〜原子が擦れるその場をシミュレーション〜

[スポンサーリンク]

軽くて強いのに摩耗に弱いチタン合金は、航空機や人工関節に広く使われています。原子が擦れ合うその場を分子動力学（MD）で覗き、摩耗の弱さの一端と、新たに原子膜上で滑り出す条件を見出した研究を紹介します。

“Atomic-Scale Investigation of the Tribological Behaviors of Titanium Alloy Interfaces with 2D Nanomaterials”,
K. Sun, S. Yu, L. Yao, Y. Xu, Langmuir 2025, 41, 16938–16951.
DOI: 10.1021/acs.langmuir.5c01352

強いのに摩耗に弱いチタン合金

チタン合金（Ti-6Al-4V、通称TC4）は軽いわりに機械的強度が高く、生体にもなじむため、航空機や人工関節、船舶部品に広く使われています。ところが熱が逃げにくく、弾性率や表面硬さも低いため、擦れ合うしゅう動部では傷や摩耗が起きやすいという泣きどころを抱えています。

この弱点を補うため、「表面を硬くして守る」という工夫が重ねられてきました。Ni/Al₂O₃＋TiO₂ナノ粒子を含む硬いNiコーティングを電気アークで表面に埋め込み、表層に硬いマルテンサイト組織をつくる手法¹や、樹脂（ポリウレタン）の層で覆って耐摩耗性を高める手法²が報告されています。こうした対策の効きめは、主に巨視的な摩擦試験で確かめられてきました。

ただし、そこから得られるのは、摩擦係数や摩耗量、摩耗痕の形といった平均化された結果が中心です。擦れている瞬間に界面で原子がどう動くかまでは、直接追いにくい。この”見えない原子/分子の世界”を計算機で動的にとらえたのが、2025年にLangmuir誌へ報告されたMDシミュレーション研究です³。

コンピュータの中の摩擦試験機

著者らは、チタン合金TC4と相手材の超硬合金YG8（WC-8%Co）を向かい合わせた仮想の摩擦試験機をMDで組みました³。材料ごとに原子間相互作用の表し方は違うので、原子間ポテンシャルをTC4はEAM、YG8はTersoff型、グラフェンはAIREBO、黒リンはSWと使い分け、異種原子間はLennard-Jonesでつないでいます。そのため、このモデルは化学反応そのものよりも、押し込みやすべりに伴う力学的接触をみるためのものです。

もちろん、ここで扱うのは長さがオングストローム（Å、原子間距離のスケール）、力がナノニュートン（nN）という極微の世界です。計算は緩和・圧縮・すべり摩擦の3段階で進め、速度・温度・押し込み深さの効き方を、少ない試行数で効率よく調べるBox–Behnken計画（応答曲面法の一種）で系統的に振りました（Fig. 1）。

Fig. 1. MDによる摩擦モデル．(a)速度・温度・押し込み深さ，(b)表面粗さ，(c，d)グラフェン(Gr)・黒リン(BP)の影響を調べる系．TC4基板にYG8圧子を押し込む構成（出典：文献[3] Fig. 1）．

応答曲面法で整理すると、この計算条件で摩擦力と摩耗原子数を最も大きく左右したのは、速度でも温度でもなく「どれだけ深く押し込むか」でした。報告された摩擦力の回帰式では、押し込み深さの項の係数が8.10と、速度（2.90）や温度（1.84）を大きく上回ります。計算では、押し込みが深くなるほど界面の応力集中が強まり、転位の発生と界面温度の上昇が同時に進むことが示されました。

二次元ナノ材料はどう滑るのか

次に著者らは、界面にグラフェン（Gr）と黒リン（black phosphorus, BP）を1枚ずつ挟んで潤滑効果を調べています。乾燥摩擦に比べ、GrでもBPでも摩擦係数は低下し、なかでもBPは「貼りついては滑って」を周期的に繰り返すスティックスリップ（stick–slip）を示しました。平均摩擦係数と時間変動はBPよりGrのほうが小さく、Grのほうがより低く安定した摩擦を示しました（ただしGr表面にかかる応力はBPより高いと報告されています）（Fig. 2）。

Fig. 2. 同一荷重での摩擦係数(COF)の時間変化．(a)乾燥摩擦・Gr・BPの比較，(b)GrとBPの重ね描き．BPは周期的なスティックスリップ（stick–slip）を示す（出典：文献[3] Fig. 10）．

もう一つの鍵は、モデルで比較された二次元材料の上下2層の重ね角（0°/90°）にありました。2枚のカードをぴたりと角をそろえて重ねる（回転角0°）と、原子の周期がそろって滑りにくくなります。一方、90°配置では周期がずれて噛み合いが弱まり、層と層がすべる「層間すべり」が起きて摩擦係数も下がります。細かな凹凸をもつ板を2枚重ねたとき、山と谷が正面から合うと引っかかりやすく、位置をずらすと引っかかりが減ってすべりやすくなる、というイメージです。

この差は、層と層をつなぐvan der Waals力が層内の結合よりずっと弱い、という二次元材料に共通する性質で説明できます。本計算条件では、黒リンは90°配置の荷重20 nNで層間すべりを見せますが、30 nNでは構造そのものが壊れてしまいます。一方、グラフェンはより高い荷重まで形を保ち、50 nN程度で層間すべりが顕著になるなど、材料ごとに「どう滑り」、「どこまで耐えるか」が違うことも見せてくれます（Fig. 3）。

Fig. 3. 二次元材料の変位．(a)BP，(b)Gr．整合接触(AA, 0°)では全体すべりのみ，非整合接触(AZ, 90°)で層間すべりが生じる。ここでいう全体すべりは上下2層がほぼ一体として動くこと，層間すべりは上下層どうしが相対的にずれることを指す（出典：文献[3] Fig. 14）．

「ずらすと滑る」はどこへつながるか

本論文が示した「重ね角しだいで滑りが変わる」性質は、今回のような計算モデルの中だけの特殊な話ではありません。結晶格子の周期が合わない面同士では、各原子が受ける横向きの力が互いに打ち消し合い、摩擦が非常に小さくなることがあります。この考え方は構造超潤滑（structural superlubricity）と呼ばれ、2004年には、グラファイト同士の接触で、格子がそろうと摩擦が大きく、向きをずらして非整合にすると摩擦が大きく下がることが実験的に示されました⁴。

この「ずれを使う」発想は、その後少しずつ大きなスケールや実材料寄りの設計へ広がってきました。2012年には、マイクロメートルサイズのグラファイト構造で超低摩擦に由来する自己復帰運動が観察され、非整合接触による超潤滑が微小機械材料の設計に関わりうることが示されました⁵。さらに、グラフェン、ナノダイヤモンド、ダイヤモンドライクカーボンを組み合わせ、摩擦中に生じるグラフェンナノスクロールが低摩擦に寄与することを実験とシミュレーションで示した研究⁶や、グラフェン/MoS₂ヘテロ構造コーティングで、構造の不均一性と荷重下での端の自己不動態化を潤滑向上に結びつける研究⁷も報告されています。つまり、2D材料は「薄くて滑る膜」ではなく、格子のずれ、層間すべり、端の化学状態まで含めて設計する材料になりつつあります。

原子の並びをずらして、摩擦を設計する

もちろん、今回のMD結果をそのまま実機の性能と読み替えることはできません。
Methodsを見る限り、本研究は長さ（Å）・時間（ps）・力（nN）を理想化したモデルであり、実環境での酸化やトライボフィルム（摩擦で生じる移着膜）の形成、欠陥や不純物の影響は含まれていません。応答曲面の当てはまり（決定係数0.9944）が高いのも、統計モデルが今回の計算データをよく再現したという意味であり、それだけで現実の摩擦の物理法則が証明されたわけではありません。つまり、この計算が示しているのは「この条件なら実機でも必ず低摩擦になる」という性能予測ではなく、押し込み、応力集中、原子層の噛み合い、層間すべりが摩擦をどう変えるかという原子スケールの筋道です。構造超潤滑も、大きな面積や実環境では、欠陥、汚染、端の結合、弾性緩和、局所的な整合接触によって失われやすい点に注意が必要です。

その意味で、押し込みで応力が集中し、原子層が噛み合い、ある角度でふっと滑り出す。その見えない瞬間を分解して原子/分子レベルでも見られることは、MDならではの強みです。

この研究の面白さは、摩擦に対して、冒頭にて示したような「表面を硬くして守る」という発想だけでなく、どこで滑らせるかを原子配列から設計する問題として見せた点にあると考えます。原子膜をどの角度で重ね、端の化学状態をどう整え、どの荷重まで構造を保たせるか。擦れる界面を原子スケールで眺めることで、固体潤滑は経験則だけでなく、配置と構造の設計問題として扱えるようになるでしょう。