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:44] koudai [電場勾配] |
wien2k:電場勾配 [2020/04/02 15:03] koudai [注意] |
||
---|---|---|---|
Line 64: | Line 64: | ||
===== SCF計算 ===== | ===== SCF計算 ===== | ||
- | - 内部構造最適化を行います< | + | - (任意)内部構造最適化を行います< |
$ init_lapw -b -numk 1000 -rkmax 7.0 | $ init_lapw -b -numk 1000 -rkmax 7.0 | ||
$ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min | $ run_lapw -ec 0.0001 -cc 0.001 -fc 0.5 -i 100 -p -min | ||
Line 71: | Line 71: | ||
$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI | $ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p -NI | ||
</ | </ | ||
+ | * 構造最適化のフラグ-minをつけて厳しい収束条件にすると、いつまでたっても計算が終わらなくなります。 | ||
===== 計算結果 ===== | ===== 計算結果 ===== | ||
Line 201: | Line 201: | ||
電気四重極モーメントは原子核の種類によって決まります。 | 電気四重極モーメントは原子核の種類によって決まります。 | ||
- | * 電気四重極モーメント一覧 https://www.sciencedirect.com/science/ | + | * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001 |
- | * 例えば43Alですと $Q=14.9~\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 250: | ||
$ grep ": | $ grep ": | ||
* k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。 | ||
- | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるのも一つの手です。 | + | * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいでしょう |