[表示 : 全て 最新50 1-99 101- 201- 301- 401- 501- 601- 701- 801- 901- 1001- 2chのread.cgiへ]
Update time : 04/11 21:43 / Filesize : 404 KB / Number-of Response : 1002
[このスレッドの書き込みを削除する]
[＋板 最近立ったスレ＆熱いスレ一覧 : ＋板 最近立ったスレ／記者別一覧] [類似スレッド一覧]

↑キャッシュ検索、類似スレ動作を修正しました、ご迷惑をお掛けしました

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

408 名前：デフォルトの名無しさん [2006/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を飛ばすのと同じように，いれてやればいい．




[ 続きを読む ] / [ 携帯版 ]
全部読む 前100 次100 最新50 ▲ [ このスレをブックマーク！ 携帯に送る ] 2chのread.cgiへ
[＋板 最近立ったスレ＆熱いスレ一覧 : ＋板 最近立ったスレ／記者別一覧]（；´∀｀）＜404KB

read.cgi ver5.27 [feat.BBS2 +1.6] / e.0.2 (02/09/03) / eucaly.net products.
担当:undef