くだすれFORTRAN（超初心者用）

くだすれFORTRAN（超初心者用） at TECH

408:デフォルトの名無しさん
06/06/08 15:28:18
>>405
ＯＫ
エネルギーを体積で微分すれば圧力になるのは熱力学の帰結だ．
しかしこれはFORTRANというより物理の質問だなｗ

＞せめて∑G/=0 のＧベクトルが０以外の時を全部足す方法が知りたいです。
ＦＯＲＴＲＡＮスレ的には，これは３次元逆格子ベクトルに関する和なので，
もっとも素朴には３重のＤＯ　ＬＯＯＰを回せばよいとい答えになる．
DO n1 = -n, n
　DO n2 = -n, n
　　DO n2 = -n, n
IF(n1==0 .and. n2==0 .and. n3==0) CYCLE　
云々
END DO
END DO
END DO
実際は結晶の対象性があるので，このループの１／８～１／２４位回せばいいはず．

実際の計算には逆格子ベクトルを求める必要がある．これは結晶構造に依存する．
もし今アルカリ金属が対象ならbcc構造ということになるだろう．
その逆格子ベクトルはfccと同じものになる．
これは固体物理の教科書に載っているから，一般式をつくるか，
手で求めて表を配列に入れることになる．

DO　LOOPの中にn1,n2,n3の組み合わせが逆格子ベクトルにふさわしいかどうかを
判別する式を，G=0を飛ばすのと同じように，いれてやればいい．

次ページ

続きを表示

1を表示