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/31 20:32] koudai [注意] |
wien2k:電場勾配 [2020/03/31 21:02] koudai [例:CaGa4] |
||
---|---|---|---|
Line 43: | Line 43: | ||
- | 電荷の収束を厳し目にして、通常のSCF計算を実行します | + | 収束を厳し目にして、通常のSCF計算を実行します |
< | < | ||
- | $ init_lapw -b -numk 1000 -rkmax 7.5 | + | $ init_lapw -b -numk 3000 -rkmax 7.5 |
- | $ run_lapw -cc 0.0001 -ec 0.00001 -i 100 -p | + | $ run_lapw -ec 0.000001 |
</ | </ | ||
Line 74: | Line 74: | ||
: | : | ||
</ | </ | ||
- | が1番目の原子のEFG計算の結果です。 | + | が1番目(今の例ではCa)の原子のEFG計算の結果です。 |
===== 電場勾配 ===== | ===== 電場勾配 ===== | ||
Line 132: | Line 132: | ||
非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。 | ||
+ | |||
+ | |||
+ | ====== 注意 ====== | ||
+ | |||
+ | * SCF計算の各ステップでの電場勾配の値は< | ||
+ | $ grep ": | ||
+ | * 収束が十分でないようでしたら、電荷の収束条件をより厳しくして再計算してください。 | ||
+ | * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください | ||
+ | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手です。 | ||
+ |