Re^2:KKR is not suitable for non-isotoropy structure? (in Japanese)

Posted on : September 30, 2002 (Mon) 01:47:33

by Kazuhiro Okamura

小谷先生、レスポンスありがとうございます。

>まずは，KKRのMT近似が悪いのではないか？どうやってチェックすべきか？
>とのとこですが、他の論文になっているような計算結果などを探して比較できないでしょうか？
>(もちろんasaでまず試すのは賢明な方法だと思います).
>まあ，ほんとはこういう場合は「他のLAPWなどでやってみる」というのが一番の
>お勧めだと思います.　まだ，そこまで整備できていなくて申し訳ないです.

Li3Nの計算事例としては、かのP. BlahaさんがWIEN(FLAPW)で
Hyperfine Interactions, VOL.51,NO.1/4 PAGE.917-923 1989 に発表してまして
ギャップは出現しています。

私は、LmtART(FP-LMTO)で計算してみましたが、メタルになってしまいました。

ひょってして、入力ファイルの原子位置が間違っていないかチェックしてみます。

なお、Li3Nの構造は、
P6/mmm(191), a=0.3648nm, c=0.3875nm
N 1a (0, 0, 0)
Li 1b (0, 0, 1/2)
Li 2c (1/3, 2/3, 0)
で

>--------------------
>それと相談してくださる場合、inputや出力の必要部分を，ちょっと添付してくださると
>助かる場合もあるかと思います．もちろん場合によりますが...

Akai-KKRのinputは
#---Li3N ---
go data/Li3N
hcp 6.894 1.062 1.0 90.0 90.0 120.0
0.001 1.0 nrl vwn mag 2nd
update 5 50 0.020
3
N 1 0 0.0 2 7 100
Li_a 1 0 0.0 2 3 100
Li_b 1 0 0.0 2 3 100
4
0.00000000 0.00000000 0.00000000 N
0.00000000 0.00000000 0.53111294 Li_a
0.28867513 0.50000000 0.00000000 Li_b
0.57735027 0.00000000 0.00000000 Li_b
でした。


※ 本題から外れますが、xbandで得られるLiCoO2 (hex)のsysファイルは変(bug?)でした。
LiCoO2の構造はR-3m(166)でして、hex(Z=3)とrhombo(Z=1)の2つのチョイスがありますが、
xbandでhexagonalで入力作成したsysファイルは、晶系がtrigonalとなり、原子位置はZ=1分しか
生成されません。ゆえに、Akai-KKR用のinpファイルも変になってしまいます。
（菱面体構造で作成すると正常です）

system data-file created by xband on Sun Sep 22 19:57:36 2002
LiCoO2.sys
Bravais lattice
10 trigonal primitive -3m D_3d
space group number (ITXC and AP)
166 428
structure type
UNKNOWN
lattice parameter A [a.u.]
5.329027287793
ratio of lattice parameters b/a c/a
1.000000000000 4.925531914894
lattice parameters a b c [a.u.]
5.329027287793 5.329027287793 26.248293981361
lattice angles alpha beta gamma [deg]
90.000000000000 90.000000000000 120.000000000000
primitive vectors (cart. coord.) [A]
0.577350269200 0.000000000000 1.641843970000
-0.288675134600 0.500000000000 1.641843970000
-0.288675134600 -0.500000000000 1.641843970000
number of sites NQ
4
IQ ICL basis vectors (cart. coord.) [A] RWS [a.u.] NLQ NOQ ITOQ
1 1 0.000000000000 0.000000000000 0.000000000000 3.039999999997 3 1 1
2 2 0.000000000000 0.000000000000 2.462765955000 2.620999999997 3 1 2
3 3 0.000000000000 0.000000000000 1.211680849900 1.899999999999 3 1 3
4 3 0.000000000000 0.000000000000 3.713851060100 1.899999999999 3 1 3
number of sites classes NCL
3
ICL WYCK NQCL IQECL (equivalent sites)
1 a 1 1
2 b 1 2
3 c 2 3 4
number of atom types NT
3
IT ZT TXTT NAT CONC IQAT (sites occupied)
1 3 Li 1 1.000 1
2 27 Co 1 1.000 2
3 8 O 2 1.000 3 4