how to get the value of band energy , exchange energy , and so on.

Posted on : October 26, 2006 (Thu) 15:01:50

by Soichiro Sasaki

全エネルギーEtot(＝ Eband　-　Epot +Estat + Exc)を

バンドエネルギーEband、交換相関エネルギーExc、ポテンシャルエネルギーEpot、
the electrostatic energy Estat等に分解して
エネルギーを見てみたいと思うのですが、
プログラム内のどのサブルーチンのどの値をみればいいか分かりません。



バンドエネルギーEbは

banden.f 内の１９，２０行目の

write(*,1000)bnde
1000 format(' band energy for the spin state=',f10.5)

で　得たエネルギーでいいのでしょうか？


また、
交換相関エネルギーExcは
totalw.f 内の22行目の

te=u+bnd2(1)+bnd2(2)

のuで与えられているエネルギーでいいのでしょうか？




すみませんが、教えてください。宜しくお願いします。

=====================================================================
大阪大学基礎工学研究科
物質創成専攻・未来物質領域
修士　2年　佐々木　総一郎
e-mail ：sasaki@aquarius.mp.es.osaka-u.ac.jp
=====================================================================