第一原理計算

質問

以下に示すのは、計算シミュレーションの記述なのですが
これは一体何の計算を行うものでしょうか？
ご面倒かと思いますが、回答お願いします

# FCC Fe
# Test the parallelism over the spins, k points and bands
# If ndtset=2, perform the computation of the Gamma dynamical matrix.
ndtset 1 # Uncomment this line
#ndtset 2 # or this line
#Dataset 1 : SCF
kptopt1 1
tolvrs1 1.0d-13

#Dataset 2 : RF at q=0 0 0 (this is always needed for IFCs)
getwfk2 1
kptopt2 2
nqpt2 1
qpt2 0.0 0.0 0.0
rfatpol2 1 1
rfdir2 1 1 1
rfphon2 1
tolvrs2 1.0d-5
#Common data
ngkpt 12 12 12 # 12 12 12 gives nkpt=182 for dataset 1
# Other possible values : 8 8 8 (giving nkpt=60 for dataset 1)
# 16 16 16 (nkpt will be on the order of 500)
ecut 35
acell 3*7.00
ixc 1
natom 1
nband 12
nline 5
nsppol 2
spinat 0.0 0.0 3.0
nstep 30
ntypat 1
occopt 7
rprim 0 .5 .5 .5 0 .5 .5 .5 0
nshiftk 4
shiftk 0.5 0.5 0.5
0.5 0.0 0.0
0.0 0.5 0.0
0.0 0.0 0.5
timopt 2
typat 1
xred 0.0 0.0 0.0
znucl 26.0

回答No.4

takanohana10

固体物理2003にも実験ノート 第一原理バンド計算プログラム「ABINIT」 (三上昌義)という特集が組まれていますね。参考まで。

回答No.3

takanohana10

たしか応用物理, 第75巻 (2006)にABINITの紹介記事が載っていました。第一原理計算について一通り書いてあるので、読んでみてください。アドバイスできるのはこのくらいかと。

回答No.2

takanohana10

どの程度ABINITのこと知っていますか?バンド計算とか、密度汎関数法とか用語だけでも知っていますか?研究室に所属するなどしていて、ABINITを知っている先輩などいますか?理論計算をするけど、ABINITをやる人間はあなたが始めてですか?

補足

今回授業の実験で使い、先生も初めてのことであるので。
このように導入として書かれていました。

To test the speed-up on a number of processor larger than 10, one
should use the files t_kpt+spin.in and t_kpt+spin.files .
The parallelisation over k-points and spins can be tested.

One is advised first to execute the above-mentioned test I ,
and see whether the output file is correct, see the fldiff file,
where the automatic analysis is performed. The test "kpt+spin"
contains parameters that make the run much longer, and much
more suitable for parallelisation, still being quite realistic.

回答No.1

takanohana10

ううむ・・・・パッケージ名も示さず、ソースを解読させようとは・・・・感心できませんね。

とりあえずFCC構造のFeについてのバンド計算ということ。”parallelism over the spins”というからスピンを考慮した計算ということですね。SCF計算ですから、繰り返し計算してポテンシャルが収束するように計算するコードではないでしょうか?

もしやABINITですかね?・・・・なんでこんな類推してあげなきゃいけないんだ。

補足

すいません、ABINITです。
k_kpt+spinというファイルでした。