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

Posted on : September 30, 2002 (Mon) 17:21:37

by Kobayashi Kazuaki

>ぼくは，あまり結晶構造がよくわかってないのすが．．．
>(空間群記号がじつはあまりよく読めません。。。)
>
>AkaiKKRでは，現状では，hcpと入れても,hexといれても同じです。(いいかげんです．．．)
>どうprimitive vectorを取るかだけの問題です。。。
>とにかく，標準出力をむしろ読んでもらって，きちんと原子位置がはいっているかどうか？
>primitive vectorが想定したものになってるか確認してもらう必要があります。
>
（中略）
>>>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
>>>でした。
>>
>>　上記で、ちょっと疑問に思ったのですが、最初に、Li3Nの構造が、
>>P6/mmmとしてあるのに、実際のinputはhcp(P6/mmc)となっている
>>ように見えるのですが、ここでなにか間違いがないでしょうか？。
>>
>>　論文：A. C. Ho, et al., PRB59, 6083(1999)（とその孫引き、
>>参考文献１１など）
>>によると、P6/mmm構造では、uパラメーターはなく、P6/mmc構造で
>>は、uパラメーターがあります。自信はないのですが、上記のhcp(P6/mmc)
>>としての、uパラメーターの与え方もおかしいように見えるのですが、
>>いかがでしょう。
>えっと,AkaiKKRの使用のいいかげんさによる誤解かもしれません。
>hexとかhcpは，primitive vectorを与えてるだけです．
>点群などは原子位置から，内部で計算してしまいます。。。
>
>小谷

　小谷さんコメント感謝です。疑問の一つは解決しました。
　ただ、それでもまだ上記には疑問が残ります。

　それは、前の筆者の書き込みでも書いた、最初、P6/mmmとして考える
なら、前の筆者の書き込みで参照した論文によると、uパラメーターは
ないはずで、その場合、上記の原子位置の数値はもっとすっきりした値
（例えば、0.2886云々は、0.3333というようになるはず）になると思い
ます。
　上記の数値の、0.57735(Li b)の値は、先の論文のP6/mmc(hcp)構造
でのuの値（=0.5776に酷似）から推測するに、上記の原子位置の数値は、
P6/mmc(hcp)としての入力データになっているようです。これは最初の
P6/mmmと矛盾するように見えます。

　また、もし上記でuパラメーターがある構造（P6/mmc）を採用すると
しても、この場合、おそらくセル内の原子数は４つより多いと思います。
　ここら辺は、自信も確証もないのですが、P6/mmc(hcp)の構造をもつも
ので、uパラメーターの類がある構造でセル内の原子数が４つのものがな
いように思えます（これはもっと調べる必要あり）。