第一原理分子動力学プログラム STATE Senri Wiki

STMの計算方法 †

準備 †

通常通りSCF計算と構造最適化を行い電荷密度を求める。

入力ファイル †

入力ファイルでパラメータを例えば

のように変更して計算するとnfchgt_r.dataが出力される。

STM計算 †

STM.fというスクリプトで

を指定してコンパイル。

$ ifort -o STM.fep STM.f

STM.fepというバイナリにnfchgt_r.dataを読みこませると

$ ./STM.fep < nfchgt_r.data

占有状態および非占有状態の等電荷面がそれぞれocc.dataとunoc.dataに出力される。

STM像の描画 †

例えばgnuplotを使うと

$ gnuplot
gnuplot> set pm3d map
gnuplot> set size ratio -1
gnuplot> splot "occ.data"