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

原子は、どうやって混ざるのか——「拡散」を実際に見る

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

🤖 AIが下書き・人間が編集

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

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

実際に混ぜてみた

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

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

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

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

これは失敗ではなく、拡散の本質だ。原子は目的地へまっすぐ進むのではなく、酔っぱらいの千鳥足（ランダムウォーク）で動く。右に一歩、左に一歩、行ったり来たり。だから、

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

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

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

ではどう加速するか。温度だ。 温度とは、原子の動きの激しさである。熱いほど原子は活発に動き回り、拡散は劇的に速くなる（温度に対して指数関数的に）。 これは言葉だけだとピンと来ないので、同じ設定を3つの温度で同時に走らせてみた。

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

いつもの正直な注意

このシミュレーションにも、いつもの但し書きがつく。原子どうしの力のモデルが現実とずれていれば、出てくる拡散の速さも現実とずれる。 きれいに動いて見えても、それは「入れた物理」の帰結にすぎない。シミュは便利な“試運転”だが、最後は実験で確かめる——という姿勢は、ここでも変わらない。

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

この記事はAIが下書きし、人間が編集・公開しています。