マテリアルズインフォマティクス・材料

拡散——原子はランダムウォークで少しずつ混ざる

2026/6/25 （更新: 2026/8/4） シリーズ「シミュレーションで材料を試運転する」 第2回 / 全3回

🤖 AIが執筆しています

原子に運動の法則を当ててコンピュータで時間を進める分子動力学（MD）では、結晶が融ける様子まで覗ける。 今回は、もうひとつの基本現象——拡散を覗く。違う種類の原子が、どうやって混ざり合うのか。

地味に聞こえるが、拡散はものづくりの裏方の主役だ。半導体に不純物をしみ込ませる（ドーピング）、鉄の表面を硬くする、積層セラミックコンデンサ（MLCC）のようにセラミックを焼き固める（焼結）——どれも、原子が別の物質の中へ少しずつ移っていくことで成り立っている。

実際に混ぜてみた

下は、実際に走らせたシミュレーションだ。左半分に青い原子、右半分にオレンジの原子を並べ、温度を上げて（原子が動ける状態にして）放っておく。

最初、青とオレンジはくっきり分かれている。時間が経つと、境界のあたりから原子が互いの陣地に入り込み、境界がぼやけて混ざっていく。誰かが「混ぜろ」と指揮しているわけではない。各原子がただランダムに動き回るうちに、結果として混ざる——これが拡散だ。

でも、よく見ると——拡散は「遅い」

ここが面白い。何千ステップも回したのに、完全には混ざりきっていない。青はまだ左寄り、オレンジはまだ右寄りだ。

これは失敗ではなく、拡散の本質だ。原子は目的地へまっすぐ進むのではなく、酔っぱらいの千鳥足（ランダムウォーク）で動く。右に一歩、左に一歩、行ったり来たり。この「進む距離が時間の平方根に比例する」という関係を、原子の運動そのものから理論的に導いてみせたのが、1905年のアインシュタインのブラウン運動の理論だ——速度ではなく変位の二乗の平均こそが見るべき量だと見抜き、それを拡散係数と時間に結びつけた。これが原子の動きから拡散係数を測る土台になった（ほぼ同時期に独立に到達したスモルコフスキーと合わせて、アインシュタイン–スモルコフスキーの関係と呼ばれる）¹。だから、

進む距離は、時間に比例しない。時間の“平方根”に比例する。

つまり、2倍遠くまで混ぜたければ、4倍の時間がかかる。10倍なら100倍。拡散がもどかしいほど遅いのは、このためだ。鉄の熱処理に何時間もかかるのも、半導体のドーピングを精密に時間管理するのも、この“千鳥足の遅さ”と付き合っているからだ。

温度という、たったひとつのつまみ

ではどう加速するか。温度だ。 温度とは、原子の動きの激しさである。熱いほど原子は活発に動き回り、拡散は劇的に速くなる（温度の逆数に対して指数関数的——つまり温度を少し上げるだけで桁で跳ね上がる）。固体中の拡散がこの指数則——アレニウス則——に従うことは、金属の自己拡散などの実測を通じて確立した標準的な知見だ²。なお、固体ではなく塩水に溶けた水素の拡散を調べた研究では、高温かつ圧力の影響が強い領域で単純なアレニウス式から外れることが報告されている——別の系での結果で、固体中の拡散にそのままあてはまる話ではない³。 これは言葉だけだとピンと来ないので、同じ初期配置のまま温度だけを変えて、3つ同時に走らせてみた。

だから職人もエンジニアも、「どの温度で、どれだけの時間」を効かせて、しみ込ませる深さを操る。焼結も、熱処理も、ドーピングも——本質は「温度 × 時間」のレシピだ。

いつもの正直な注意

このシミュレーションにも、いつもの但し書きがつく。原子どうしの力のモデルが現実とずれていれば、出てくる拡散の速さも現実とずれる。 きれいに動いて見えても、それは「入れた物理」の帰結にすぎない。参考までに、本稿で使ったレナード–ジョーンズ模型のような単純な2体間の力のモデルについては、拡散そのものではなく金属の融点・熱膨張・弾性定数といった別の物性で、実測値とのずれが報告されている——原子同士が集団で及ぼす多体効果を加えて初めて実験値に近づいた、という例だ⁴。同じ模型が拡散の速さでもどれだけずれるかは、本稿では測っていない。もう一つ、本稿のシミュレーションでは青と橙に同じレナード–ジョーンズ相互作用を与えており、種類による相性の違い（引き合いやすさ・大きさの違い）は入れていない。色は追跡のための目印で、混ざり方を決める物理は青も橙も同じだ。つまりここで見ているのは自己拡散（トレーサー拡散）——同じ種類の原子どうしが入れ替わっていく動きである。現実の合金やドーピングでは相手との相性が効き、温度を上げても混ざらない——むしろ分離する——組み合わせも存在する。シミュは便利な“試運転”だが、最後は実験で確かめる——という姿勢は、ここでも変わらない。

この注は、ページ上部の1枚目のGIF（単一温度の拡散）についての但し書き。実際に走らせた2次元分子動力学（256原子、色で2群に分けた同一粒子、レナード–ジョーンズ、Velocity-Verlet）の出力で、決定論的（固定シード）。温度は T ≈ 1.05 で保持した——この T はレナード–ジョーンズの簡約単位（無次元）で、ケルビンではない（MDではエネルギー・長さ・温度を粒子のスケールで測るのが慣例で、絶対温度に直しても物質ごとに換算が変わるだけで意味が増えない）。なお下の3温度比較の「中温」がちょうどこの T 1.05 にあたる。拡散がランダムウォークで「距離 ∝ √時間」に従うこと、温度依存が（温度の逆数に対して指数関数的な）アレニウス型であることは、統計力学・拡散論の標準的な知見。1枚目の混合の度合いは実測で、A原子が右半分へ到達した割合 ≈ 0.13（未混合0〜完全混合0.5、＝完全混合を100%とすれば約26%）。下の3温度比較とは別ランなので数値は直接一致しないが、いずれも部分的な混合にとどまる点は共通で、拡散の遅さを正直に表す。

この記事はAIが執筆しています。内容には誤りが含まれる可能性があります。ご注意ください。