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/04/02 13:42] koudai |
wien2k:電場勾配 [2020/04/02 15:42] koudai [例:α-CaAlF5] |
||
---|---|---|---|
Line 61: | Line 61: | ||
* https:// | * https:// | ||
+ | * この論文によると、構造最適化して計算したほうが良いということですが、今回はチュートリアルですので構造最適化前の結果を示します。 | ||
===== SCF計算 ===== | ===== SCF計算 ===== | ||
- | - 内部構造最適化を行います< | ||
- | $ init_lapw -b -numk 1000 -rkmax 7.0 | ||
- | $ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min | ||
- | </ | ||
- 収束を厳し目にして、SCF計算を再実行します< | - 収束を厳し目にして、SCF計算を再実行します< | ||
- | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI | + | $ init_lapw -b -numk 2000 -rkmax 7.0 |
+ | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p | ||
</ | </ | ||
- | |||
===== 計算結果 ===== | ===== 計算結果 ===== | ||
Line 201: | Line 198: | ||
電気四重極モーメントは原子核の種類によって決まります。 | 電気四重極モーメントは原子核の種類によって決まります。 | ||
- | * 電気四重極モーメント一覧 https://www.sciencedirect.com/science/ | + | * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001 |
- | * 例えば43Caですと $Q=4.3~\mathrm{[\times 10^{-28}m^2]}$ です | + | * 例えば27Alですと $Q=0.15~\mathrm{[\times 10^{-28}m^2]}$ |
* bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です | * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です | ||
Line 250: | Line 247: | ||
$ grep ": | $ grep ": | ||
* k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | ||
- | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手です。 | + | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいでしょう |