Differences

This shows you the differences between two versions of the page.

--- quantumespresso:respack [2020/08/28 14:27]
koudai [計算の収束]
+++ quantumespresso:respack [2020/09/19 19:42]
koudai [計算の収束]
@@ Line 8: / Line 8: @@
     * https://www.cond-mat.de/events/correl11/
   * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います
-  * メモリを非常に多く使うため、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません
+  * 計算は非常に重いので、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません
   * f電子系やスピン軌道相互作用がある系には非対応
   * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます
@@ Line 216: / Line 216: @@
 最初は次の値を使用します
-  * k点数 ... SCF計算で全エネルギーを収束させたもの
+  * k点数 ... ワニエでバンドが再現できる最低限の数
-  * Ecut_for_eps ... ecutwfcの1/20
+  * Ecut_for_eps ... 2, 3, 4 Ryあたりで計算して様子を見る
-次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットして収束を確認するのが便利です)
+次の手順で収束させます
-  - 多めの数のnbndでSCF計算を実行する（最初は非占有状態のバンドが50本くらいになるか、あるいはフェルミ準位より2,30eV程度上のバンドが入るくらいで試すと良い）
+  - 多めの数のnbndでSCF計算を実行する（最初はフェルミ準位より40~50eV程度上のバンドが入るくらいで試すと良い）
   - （制限RPAの場合）相互作用の大きさを求めたい軌道のワニエ基底を計算する。
-  - N_CALC_BAND に関して分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
+  - N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
-  - Ecut_for_eps に関して計算を収束させる
+  - Ecut_for_eps に関して計算を収束させる（）
   - N_CALC_BANDの収束をチェックする
-  - SCF計算でk点数を増やし、分極関数を収束させる
+  - SCF計算でk点数を増やし、計算を収束させる
-有効相互作用を求める場合は、さらに Num_freq_grid（周波数のグリッド数）を増やして有効相互作用の値の収束を確認しましょう。
+  * 有効相互作用を求める場合は、さらに Num_freq_grid（周波数のグリッド数）を増やして有効相互作用の値の収束を確認しましょう。
+  * 計算が重いので、有効相互作用の有効数値は2桁出すのが精一杯かもしれません

雑多な記録

User Tools

Site Tools

Differences

Page Tools