On the calculation of PbTiO3

Posted on : April 25, 2005 (Mon) 16:59:05

by Yukihiro OKuno

富士フィルム 奥野といいます。

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

AKAIKKRを利用して、PbTiO3の計算をおこなっています。(誘電体には
固溶体がおおいためCPAなどの利用で上手くできればハッピーかと考え
ています。)

が、状態密度をみてみるとエネルギー＝０(Fermiエネルギー)のところに
状態をもってしまい絶縁体にならないというおかしな結果がでてしまい
ます。(ちなみにewidth=1.0 として Tiのセミコア状態をいれなくても
同じ結果になりました。)
どうも、TIのd軌道がフェルミエネルギーのところにういてしまってきている
ようにみえます。(nrl, sra どちらの場合も)

PbZrO3の場合は状態密度をFermi-Energyのところにもたずきちんと絶縁体
になってくれます。

計算じたいはedeltを0.005より少しずつ値をちいさくさせて計算を収束
させております。

これは、どういう部分に注意が足りないからこのような事態がおこるのでしょうか?



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

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

7.370000000000

#- edelt ewidth nrl/sra mjw/vbh/vwn mag/nmag/-mag/rvrs/kick ------#
0.001 2.00 nrl vwn nmag

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

#- number of type -----------------------------------------------#
3
# type components rtin field l_max Z concentration----------#
Pb 1 0. 0.000 3 82 1.00
O 1 0. 0.000 2 8 1.00
Ti 1 0. 0.000 3 22 1.00

# number of atoms------------------------------------------------#
5
#- position ---------------------------------------------- type -#
0.000000000000a 0.000000000000b 0.000000000000c Pb
0.500000000000a 0.500000000000b 0.500000000000c Ti
0.000000000000a 0.500000000000b 0.500000000000c O
0.500000000000a 0.000000000000b 0.500000000000c O
0.500000000000a 0.500000000000b 0.000000000000c O