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 | ||
quantumespresso:respack [2020/09/19 19:02] koudai [概要] |
quantumespresso:respack [2020/09/19 19:37] koudai [計算の収束] |
||
---|---|---|---|
Line 216: | Line 216: | ||
最初は次の値を使用します | 最初は次の値を使用します | ||
- | * k点数 ... SCF計算で全エネルギーを収束させたもの | + | * k点数 ... ワニエでバンドが再現できる最低限の数 |
- | * Ecut_for_eps ... ecutwfcの1/ | + | * Ecut_for_eps ... 2, 3, 4 Ryあたりで計算して様子を見る |
- | 次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットして収束を確認するのが便利です) | + | 次の手順で収束させます |
- | - 多めの数のnbndでSCF計算を実行する(最初は非占有状態のバンドが50本くらいになるか、あるいはフェルミ準位より2,30eV程度上のバンドが入るくらいで試すと良い) | + | - 多めの数のnbndでSCF計算を実行する(最初はフェルミ準位より30~40eV程度上のバンドが入るくらいで試すと良い) |
- (制限RPAの場合)相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | - (制限RPAの場合)相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | ||
- | - N_CALC_BAND に関して分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | + | - N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う |
- Ecut_for_eps に関して計算を収束させる | - Ecut_for_eps に関して計算を収束させる | ||
- N_CALC_BANDの収束をチェックする | - N_CALC_BANDの収束をチェックする | ||
- | - SCF計算でk点数を増やし、分極関数を収束させる | + | - SCF計算でk点数を増やし、計算を収束させる |
- | 有効相互作用を求める場合は、さらに Num_freq_grid(周波数のグリッド数)を増やして有効相互作用の値の収束を確認しましょう。 | + | |
+ | * 有効相互作用を求める場合は、さらに Num_freq_grid(周波数のグリッド数)を増やして有効相互作用の値の収束を確認しましょう。 | ||
+ | * 計算が重いので、有効相互作用の有効数値は2桁出すのが精一杯かもしれません | ||