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/08/17 18:57] koudai [分極関数の計算] |
quantumespresso:respack [2020/09/19 19:20] koudai [計算の収束] |
||
---|---|---|---|
Line 8: | Line 8: | ||
* https:// | * https:// | ||
* 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | ||
- | * メモリを非常に多く使うため、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません | + | * 計算は非常に重いので、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません |
* f電子系やスピン軌道相互作用がある系には非対応 | * f電子系やスピン軌道相互作用がある系には非対応 | ||
* 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます | * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます | ||
Line 194: | Line 194: | ||
& | & | ||
- | 事前に制限RPA法(calc_chiqw で flg_cRPA=1)を使った計算が必要です。 | + | 事前に制限RPA法(calc_chiqw で flg_cRPA=1)を使ったすべてのq (flg_calc_type = 0) での計算が必要です。 |
* 直接相互作用< | * 直接相互作用< | ||
Line 208: | Line 208: | ||
====== 計算の収束 ======= | ====== 計算の収束 ======= | ||
- | 分極関数の計算で収束させなければならないパラメータは次の2つです | + | 分極関数の計算で収束させなければならないパラメータは次の3つです |
* k点数 | * k点数 | ||
Line 217: | Line 217: | ||
* k点数 ... SCF計算で全エネルギーを収束させたもの | * k点数 ... SCF計算で全エネルギーを収束させたもの | ||
- | * Ecut_for_eps ... ecutwfcの1/ | + | * Ecut_for_eps ... 2, 3, 4 Ryあたりで計算して様子を見る |
- | 次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットして収束を確認するのが便利です) | + | 次の手順で収束させます |
- | - 多めの数のnbndでSCF計算を実行する(最初は非占有状態のバンドが50本くらいになるか、あるいはフェルミ準位より2, | + | - 多めの数のnbndでSCF計算を実行する(最初はフェルミ準位より2, |
- (制限RPAの場合)相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | - (制限RPAの場合)相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | ||
- N_CALC_BAND に関して分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | - N_CALC_BAND に関して分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う |