分子動力学法（molecular dynamics, MD）は、原子に働く力を計算しながら、原子の位置と速度を少しずつ時間発展させていくシミュレーション手法です。 ざっくり言えば、原子の動きをコマ送りで追いかける計算です。
計算化学の入門でよく見かける言葉ですが、最初は
- 何を計算しているのか
- 構造最適化と何が違うのか
- 実際に何が分かる手法なのか
が少し分かりにくいと思います。
この記事では、分子動力学法の基本を初学者向けに整理します。 細かい実装や高度な理論に入る前に、まずはMDでは何をしているのかをつかむことが目的です。
分子動力学法とは何か
分子動力学法は、原子配置から力を計算し、その力に従って原子を動かす方法です。 流れをかなり単純化すると、MDでは次を繰り返します。
- ある時刻での原子の位置を用意する
- その原子配置から、各原子に働く力を計算する
- 力から加速度を求める
- 加速度を使って速度と位置を少しだけ更新する
- 次の時刻でも同じことを繰り返す
これを何万ステップ、何十万ステップと進めることで、原子の運動の履歴を得ます。 各時刻の1コマ1コマの構造はスナップショットと呼ばれ、その並び全体をトラジェクトリと呼びます。
つまりMDは、1つの安定構造を探す計算というより、時間とともに系がどう動くかを追う計算です。 ここでいう系とは、シミュレーションの対象になっている原子や分子の集まり全体のことです。
原子は何に従って動くのか
原子を適当に動かしているわけではありません。 MDでは、原子に働く力をポテンシャルから計算します。
ここでいうポテンシャルは、原子配置に対してエネルギーを与える関数のことです。 また、実際のMD計算では、このポテンシャルの形やパラメータ、原子タイプ、結合情報などを含めたモデル一式を力場と呼ぶことがよくあります。 そして原子に働く力は、そのポテンシャルエネルギーの座標に対する勾配から求められます。
式で書くと、力は次の形です。
$$
F = -\nabla U
$$
U: ポテンシャルエネルギーF: 原子に働く力
この式の意味を直感的に言うと、エネルギー地形の傾きに沿って原子が動くということです。
山の斜面にボールを置くと、傾きに沿って転がります。 MDでも似た考え方で、エネルギーがどう変わるかをもとに力を計算し、その力で原子を動かします。
構造最適化との違い
初学者が最初に混同しやすいのが、構造最適化との違いです。
どちらも原子座標を更新するので似て見えますが、目的はかなり違います。
| 手法 | 主な目的 | 原子の動かし方 |
|---|---|---|
| 構造最適化 | エネルギーの低い構造を探す | エネルギーが下がる方向へ構造を調整する |
| MD | 時間発展を追う | 力に基づいて位置と速度を時間発展させる |
構造最適化が安定構造を探す計算なのに対し、MDは動いている系を追う計算です。 また、実際のMDでは少なくとも次のような設定が必要になります。
- 初期構造
- 初期速度
- ポテンシャルまたは力場
- 時間刻み
- ステップ数
- 温度や圧力の制御条件
- 周期境界条件の有無
つまりMDは、「ボタン1つで原子の動きを見る手法」ではなく、 どんなモデルで、どんな条件で、どれくらいの時間だけ系を追うかを明示して進める計算です。
MDで何が得られるのか
MDの直接の出力は、各時刻における原子の位置、速度、エネルギーなどです。 それ自体が最終目的というより、そこから後で解析を行います。
たとえば、MDの結果からは次のようなものを調べられます。
- 原子や分子がどう動いたか
- 温度やエネルギーが安定しているか
- 密度がどのくらいか
- 動径分布関数（RDF）
- 平均二乗変位（MSD）
- 拡散係数の見積もり
ここでは、MDはトラジェクトリを作り、それを解析して物性や構造の情報を得る手法と捉えておけば足ります。
MDはとても便利ですが、現実をそのまま完全に再現するわけではありません。
理由はいくつかあります。
- ポテンシャルや力場に近似が入る
- 扱える時間スケールに限界がある
- 扱える系の大きさに限界がある
- 温度や圧力の制御も理想化を含む
つまりMDは、現実そのものではなく、あるモデルのもとで時間発展を近似的に再現する方法です。 大事なのは「MDだから正しい」と思い込むことではなく、 どんな前提で計算した結果なのかを理解して使うことです。
まとめ
分子動力学法（MD）は、原子に働く力を計算しながら、原子の位置と速度を少しずつ更新していくシミュレーション手法です。
ポイントをまとめると次の通りです。
- MDは原子の時間発展を追う計算
- 力はポテンシャルエネルギーから求める
- 出力はトラジェクトリであり、そこから物性や構造を解析する
- 構造最適化とは目的が違う
- 便利な手法だが、モデルや近似に依存する
MDの入口としてまず押さえたいのは、MDは原子を動かす計算であり、その動きの履歴を解析する手法だという点です。
次の記事では、MDを使うと実際に何が分かるのかをもう少し具体的に整理します。
コメント