Verlet法による調和振動子のMD計算

3.プログラムの実行

  • プログラムをビルド、実行する。
  • プロジェクト・ディレクトリの中の、プロジェクト名と同じディレクトリの中に “harmonic.csv”というファイルが生成されている。ここに、運動エネルギー、ポテンシャルエネルギー、ハミルトニアンの時間変化が記録されている。
  • Microsoft Office Excelを起動し、“harmonic.csv”を開く。
  • 「挿入」→「散布図」を選択すると、運動エネルギー、ポテンシャルエネルギー、ハミルトニアンの時間変化のグラフが表示される(下図参照)。
  • ハミルトニアンがほぼ一定に保たれていることを確認しよう。
  • 周期が理論値と一致していることを確認しよう。

excel1