This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:電場勾配 [2020/03/28 15:33] koudai [電場勾配] |
wien2k:電場勾配 [2020/03/31 20:32] koudai [例:CaGa4] |
||
---|---|---|---|
Line 10: | Line 10: | ||
* また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください | * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください | ||
* 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。 | * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。 | ||
+ | |||
+ | ====== 例:CaGa4 ====== | ||
+ | |||
+ | 以下の論文の結晶構造から作成しました。 | ||
+ | |||
+ | * https:// | ||
+ | |||
+ | <file - CaGa4.struct> | ||
+ | blebleble | ||
+ | CXZ LATTICE, | ||
+ | MODE OF CALC=RELA unit=bohr | ||
+ | | ||
+ | ATOM -1: X=0.00000000 Y=0.00000000 Z=0.00000000 | ||
+ | MULT= 1 ISPLIT=15 | ||
+ | Ca | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | ATOM -2: X=0.40503000 Y=0.76923000 Z=0.00000000 | ||
+ | MULT= 2 ISPLIT=15 | ||
+ | -2: X=0.59497000 Y=0.23077000 Z=0.00000000 | ||
+ | Ga | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | ATOM -3: X=0.00000000 Y=0.50000000 Z=0.23283000 | ||
+ | MULT= 2 ISPLIT=15 | ||
+ | -3: X=0.00000000 Y=0.50000000 Z=0.76717000 | ||
+ | Ga | ||
+ | LOCAL ROT MATRIX: | ||
+ | | ||
+ | | ||
+ | | ||
+ | </ | ||
+ | |||
+ | この物質のWIEN2kによるNQRの研究は、以下の論文の5.1節が参考になります | ||
+ | |||
+ | * https:// | ||
+ | |||
+ | |||
+ | 電荷の収束を厳し目にして、通常のSCF計算を実行します | ||
+ | |||
+ | < | ||
+ | $ init_lapw -b -numk 1000 -rkmax 7.5 | ||
+ | $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -p | ||
+ | </ | ||
====== 電場勾配 ====== | ====== 電場勾配 ====== | ||
- | case.scf0 を開いてください。 | + | SCF計算が終了したら |
以下はCaGa4の例です。 | 以下はCaGa4の例です。 | ||
Line 35: | Line 81: | ||
: | : | ||
</ | </ | ||
- | がEFG計算の結果です。 | + | が1番目の原子のEFG計算の結果です。 |
===== 電場勾配 ===== | ===== 電場勾配 ===== | ||
Line 84: | Line 130: | ||
0.0000 | 0.0000 | ||
- | 主軸回転後の向きベクトルです。 | + | 主軸回転後の軸の向きです。 |
- | 主軸回転後のZ軸の方向(主軸方向)は (0.6989, 1, 0) 、y軸の方向は(0, | + | 主軸回転後のz軸の方向(主軸方向)は (0.6989, 1, 0) 、y軸の方向は(0, |