分子動力学シミュレーター 戻る
分子シミュレーション

分子動力学シミュレーター

レナード=ジョーンズポテンシャルとVerlet積分で粒子間の引力・斥力を計算。温度スライダーを動かして気相・液相・固相の転移をリアルタイムで観察できます。

パラメータ

操作

気相
計算結果
運動温度 T*
運動エネルギー
ポテンシャル
全エネルギー
粒子シミュレーション(青=低温 → 赤=高温)
エネルギー履歴
理論・主要公式
$$V(r)=4\varepsilon\left[\left(\frac{\sigma}{r}\right)^{12}-\left(\frac{\sigma}{r}\right)^{6}\right]$$

引力最小点 rmin = 21/6σ ≈ 1.122σ

分子動力学シミュレーターとは

🙋
分子動力学シミュレーションって、教科書で見る気体・液体・固体の状態変化を、実際にパソコンで再現できるんですか?
🎓
その通り!このシミュレーターは、たくさんの粒子(分子)の動きをコンピュータで計算して、温度や密度を変えるとどうなるかをリアルタイムに見せてくれるんだ。例えば、上の「温度」スライダーをぐっと下げてみて。粒子がバラバラに動いていたのが、だんだん集まって規則正しく並び始める、つまり「気体→固体」への転移が目の前で起こるよ。
🙋
え、そうなんですか!でも、粒子同士が引き合ったり反発したりする力は、どうやって決めてるんですか?
🎓
そこがこのシミュレーションの肝だね。粒子間の力は「レナード=ジョーンズポテンシャル」という式で決めている。パラメータの「ε(イプシロン)」は引力の強さ、「σ(シグマ)」は粒子の大きさを表すんだ。εを大きくすると粒子がくっつきやすくなり、σを大きくすると大きい粒子をシミュレートできる。実際に動かして、液体ができる条件を探してみよう。
🙋
粒子の数「Number of Particles」も変えられるけど、数を増やすと何が変わるんですか?計算が重くなるだけ?
🎓
いいところに気がついたね。数を増やすと、少数の粒子では見えなかった「巨視的」な現象が現れやすくなるんだ。例えば、粒子数が少ないとただの塊にしか見えないけど、数を数百に増やして温度を下げると、きれいな六方最密充填(ハニカム構造)のような結晶が自然に形成される様子が観察できるよ。相転移の様子がよりはっきりするんだ。

よくある質問

温度スライダーで系全体の運動エネルギー(粒子速度)を調整できます。低温では粒子が規則正しく並んで固相(結晶)に、中温では流動的な液相に、高温では自由に飛び回る気相へと変化する様子をリアルタイムで観察できます。
現バージョンでは粒子数は固定です。ただし、粒子の初期配置や温度を変えることで、同じ粒子数でも気体・液体・固体の異なる状態を再現できます。粒子数を増やすと計算負荷が高くなるため、パフォーマンスを考慮して設定されています。
粒子間に働く距離に応じた力を表すモデルです。近距離では強い斥力(衝突防止)、中距離では弱い引力(凝集)、遠距離ではほぼゼロになります。このシンプルな式で、分子間のファンデルワールス力を模擬しています。
温度が高すぎると粒子の速度が大きくなり、計算の時間刻み(タイムステップ)に対して移動量が過大になるためです。その場合、温度スライダーを少し下げるか、ブラウザをリロードして初期状態に戻してください。

実世界での応用

ナノ材料設計:カーボンナノチューブやグラフェンシートの機械的・熱的性質を、原子レベルでのシミュレーションで予測します。LJポテンシャルは炭素原子間のファンデルワールス相互作用をモデル化するのに使われます。

創薬・タンパク質折りたたみ:薬剤分子が標的タンパク質のどの部位に結合するか、またタンパク質がどのように立体構造を形成するかを、水分子も含めた大規模な分子動力学シミュレーションで解析します。生体分子の相互作用の基礎モデルとしてLJポテンシャルが活用されます。

潤滑剤・界面科学:金属表面と潤滑油分子の間の摩擦や潤滑メカニズムを解明します。分子がどのように配列し、せん断力に抵抗するかをシミュレーションすることで、低摩擦・高耐久な潤滑剤の開発に貢献します。

ガス貯蔵・分離材料:多孔性材料(MOFなど)の細孔内でのガス分子(水素、二酸化炭素など)の吸着・拡散挙動をシミュレーションします。材料の構造を変えて、どれだけ効率的にガスを貯められるか、事前に評価できます。

よくある誤解と注意点

まず、「シミュレーション結果がそのまま実験データ」と思わないことが大事だ。このシミュレーターは、あくまで「レナード=ジョーンズ(LJ)ポテンシャルという単純化されたモデル」の世界を再現している。例えば、水の相転移を厳密に再現したいなら、水分子の極性(プラスとマイナスの偏り)を考慮したより複雑な力のモデルが必要になる。LJモデルはアルゴンなどの希ガス原子の挙動を理解するための「入り口」だと心得よう。

次に、パラメータ設定で陥りがちなのが「極端な値」を設定してしまうこと。例えば、温度を一気に0K(絶対零度)近くに下げると、粒子がほとんど動かなくなり、シミュレーションが「凍りついて」見えてしまう。逆に、ε(引力の強さ)を極端に大きくすると、粒子が強固に固まって動きがなくなり、逆に相転移の「過程」が見えにくくなる。まずはデフォルト値から少しずつ(例えば温度を10%ずつ、εを1.2倍ずつなど)変化させて、その影響を観察するのがコツだ。

最後に、「粒子数が多ければ多いほど良い」という思い込み。確かに、数百〜数千粒子にすると液体や固体の「塊」としての振る舞いが明確になるが、その分計算負荷は増大する。実務では、求められる精度と計算リソースのバランスが重要。例えば、ナノスケールの微粒子の凝集を調べるのであれば、数十粒子の挙動を詳細に追跡する方が本質的なケースもある。このツールで粒子数を変えながら、現象の本質を捉えるのに必要な最小限の粒子数はどれくらいか、感覚を掴んでおくといい。

使い方ガイド

  1. 粒子数スライダー(nMDNum)を2〜1000の範囲で設定。アルゴン原子シミュレーションの場合は256個から開始
  2. 温度スライダー(tempMDNum)を50〜5000Kで調整。低温で固体状態、中温で液体、高温で気体相への遷移を観察
  3. レナード=ジョーンズパラメータのεスライダー(epsMDNum)をkJ/molで設定し、分子間相互作用エネルギーを定義
  4. σスライダー(sigmaMDNum)をÅ単位で0.3〜0.5に設定し、粒子径を決定。Verlet積分により時間ステップ0.005psで計算実行

具体的な計算例

アルゴン(Ar)の気液共存状態を再現する場合、粒子数256個、温度94.4K、ε=119.8kJ/mol、σ=3.405Åを設定。初期密度0.8g/cm³で10000ステップシミュレーション実行すると、ポテンシャルエネルギーは約−600kJ/molに収束。104K未満では液体相(配位数10〜11)、120K以上では気体相(配位数2〜3)へ転移し、相図の三重点近傍の物理化学的挙動を定量確認できます

実務での注意点