小正準集団

系を粒子数一定、全エネルギー一定として取り扱う。真空中の粒子の運動などはこの方法によって取り扱うのが便利である場合が多い。温度は1自由度当たりの運動エネルギーの2倍の平均値として求められる。

正準集団

系を粒子数一定、温度一定として取り扱う。実際の表面系の多くは、下地が熱浴としてはたらくために温度一定の取り扱いをする方が便利な場合が多い。

小正準集団において、系の温度を一定に保つために人為的に速度をスケールする。

具体的には、系の実際の温度（運動エネルギーの平均から求まる温度）を${T}$、設定温度を$T_c$とするときに、各時間ステップごとに

${\vec{v}'}_j={\vec{v}}_j\displaystyle\sqrt{\frac{T_c}{T}}$

のように速度${\vec{v}}_j$をスケールすると系の温度を$T_c$に保つことができる。ただし、${M}$を系の自由度の数、${K}$を系の運動エネルギーとするとき、

$k_{\rm B}T=\displaystyle\frac{2}{M}\langle K \rangle$

で定義される。ここで、$k_{\rm B}$はボルツマン定数である。

平衡状態を保つようにして温度を徐々に下げていくと粒子は最も安定な（あるいは準安定な）位置に落ち着く。このようにして物質の安定な構造を得ることができる。

速度スケーリングの方法で、1次元W型ポテンシャル中を運動する1個の粒子の最安定状態を求めるプログラムを作成する。

1. 速度スケーリングの方法で一定温度のシミュレーションを行うことができることを示す。
2. 安定な吸着位置が正方格子を組んでいる表面の上に原子をランダムに配置する。以下の三つの模型について、次のような条件のプログラムをそれぞれ作成してシミュレーションを行い、原子の位置、速度、運動エネルギー、ポテンシャルエネルギー、全エネルギーの時間変化を比較する。

• 模型
  • モデルA：隣接原子と「ばね」で繋がった原子の運動
  • モデルB：正方格子状の外場の中の原子
  • モデルC：正方格子状に並んだ下地原子と相互作用する原子の運動
• プログラムの条件
  1. 一定エネルギーのシミュレーション
  2. 一定温度のシミュレーション
  3. 徐々に温度を下げていくシミュレーション

特に、一定温度のシミュレーションと徐々に温度を下げていくシミュレーションについては、原子の平均的な位置のずれを計算し、温度に対してプロットせよ。同時にスナップショットも出力せよ。