pythonで超簡単な分子動力学プログラムを書く - Qiita

initial. d の各行は 原子のｘ座標、ｙ座標、速度のｘ成分、ｙ成分 の4つの項目からなり、スペースで区切られている必要がある。 物体を描画できるようにキャンバスを開き、ファイル（initial.d)から原子配置（と速度）情報を読み取って原子を初期位置に配置するプログラム例を下に示す。 後々のため、ダミーのボタンも用意している（ここではボタンが押されても特に何もしない）。 ここではinitial.dが3行なので原子数は３つとなるが、行数を増やすことで原子数が増やせる。プログラム中では各行を読み込む度に配列サイズを増やしている。 座標の単位はÅ、速度の単位はÅ／sとしている。 プログラム例を示す。 import tkinter as tk def dummy(event): print('button is pressed') def drawatom(x,y): scl=10 rad