On the TiC and BaTiO3 calculation ( in japanese)

Posted on : April 26, 2005 (Tue) 17:09:40

by Yukihiro OKuno

富士フィルム奥野です。

おせわになっております。

先日PbTiO3の計算に関して投稿したものです。
絶縁体であるものがE＝0に状態密度をもってしまうこと
に関してですが、同様の計算をBaTiO3 (cubic 相)、TiC
に関して試行してみました。で、やはりTiC BaTiO3が
E＝0に状態密度をもってしまいました。
結果をみると、Ti のd軌道のウェイトがやけにE＝0でおおきく
なってしまっており、これが諸悪の根元かとおもっています。
計算をみていると、とくにコアとして採択している電子の数や
軌道などに問題はないようなのですが、たとえば計算結果を
得る過程などに問題があるとかいうことはないでしょうか?
私は計算過程でedelta=0.005 -> 0.003 ->0.002 などと
徐徐にちいさくして(より大きなewidthの値の収束結果を利用して)
収束させていっています。

どなたか、同様の計算などをおこなって同じ問題に遭遇されかた
などいましたら返事もらえると幸です。

以下 TiC の入力ファイル


#----------------------- input data ----------------------------#
# go/ngo/dos/dsp/spc file name
dos data/TiC

#- primitive vector --------------------------------------#
# this aux option is only effective for newest versions of AkaiKKR
aux
0.000000000000 0.500000000000 0.500000000000
0.500000000000 0.000000000000 0.500000000000
0.500000000000 0.500000000000 0.000000000000
8.160000000000

#- edelt ewidth nrl/sra mjw/vbh/vwn mag/nmag/-mag/rvrs/kick ------#
0.002 1.5000 sra vwn nmag

#- init/1st/2nd quit/update 0/1/2/../t/l/m/h/u iteration pmix -#
1st update t 100 0.01000

#- number of type -----------------------------------------------#
2
# type components rtin field l_max Z concentration----------#
C 1 0. 0.000 2 6 100.00
Ti 1 0. 0.000 2 22 100.00

# number of atoms------------------------------------------------#
2
#- position ---------------------------------------------- type -#
0.000000000000a 0.000000000000b 0.000000000000c Ti
0.500000000000a 0.500000000000b 0.500000000000c C