Differences

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

--- quantumespresso:respack [2020/09/19 19:23]
koudai [計算の収束]
+++ quantumespresso:respack [2020/09/20 19:25]
koudai [計算の収束]
@@ Line 8: / Line 8: @@
     * https://www.cond-mat.de/events/correl11/
   * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います
-  * 計算は非常に重いので、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません
+  * 計算は非常に重いので、研究室にあるような普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらい行けるかもしれません
   * f電子系やスピン軌道相互作用がある系には非対応
-  * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます
+  * 擬ポテンシャルはノルム保存型を使用する必要があります。例えば以下のサイトからダウンロードできます
     * http://www.pseudo-dojo.org/ ... optimized norm-conserving Vanderbilt PPs
     * https://www.quantum-espresso.org/pseudopotentials/hartwigesen-goedecker-hutter-pp ... Hartwigesen-Goedecker-Hutter PPs
@@ Line 115: / Line 115: @@
 2 3 4 5 6 7 8 9    !o: N_write_wannierを指定したときの、計算したいワニエ軌道の番号
 &param_chiqw
-  Ecut_for_eps           = 3.6d0,   !o: 分極関数のカットオフ (Ry) 記述がなければSCF計算のecutwfcの値の1/10になる
+  Ecut_for_eps           = 3.0d0,   !o: 分極関数のカットオフ (Ry) 記述がなければSCF計算のecutwfcの値の1/10になる
   Num_freq_grid          = 70,      !o: 計算周波数の総数（メッシュは log ω に対して切られる）
   N_CALC_BAND            = 50,      !o: 分極計算で考慮されるバンドの総数。SCF計算のnbndの数よりも小さくする。記述がなければSCF計算のnbndの値になる
@@ Line 126: / Line 126: @@
   ttrhdrn_dmna           = 1.0d-3,  !o: テトラへドロン計算での縮退判定パラメータ 1 (eV)
   ttrhdrn_dmnr           = 1.0d-3,  !o: テトラへドロン計算での縮退判定パラメータ 2 (eV)
-  flg_cRPA               = 0,       !o: 制限 RPA オプション (0: 通常 RPA, 1: 制限 RPA)
+  flg_cRPA               = 0,       !o: 制限 RPA オプション (0: 通常RPA, 1: 制限RPA)
   flg_calc_type          = 0/       !o: 0: すべてのq, 1: q=0のみ, 2: 任意のq（マニュアル14章）
 &param_calc_int
@@ Line 210: / Line 210: @@
 分極関数の計算で収束させなければならないパラメータは次の3つです
-  * k点数
+  * k点数（SCF計算）
   * N_CALC_BAND ... バンドの数
   * Ecut_for_eps ... 分極関数のカットオフ (Ry)
@@ Line 216: / Line 216: @@
 最初は次の値を使用します
-  * k点数 ... SCF計算で全エネルギーを収束させたもの
+  * k点数 ... ワニエでバンドが再現できる最低限の数
-  * Ecut_for_eps ... 2, 3, 4 Ryあたりで計算して様子を見る
+  * Ecut_for_eps ... 3, 4, 5 Ryあたりで計算して様子を見る
 次の手順で収束させます
-  - 多めの数のnbndでSCF計算を実行する（最初はフェルミ準位より2,30eV程度上のバンドが入るくらいで試すと良い）
+  - 多めの数のnbndでSCF計算を実行する（最初はフェルミ準位より50eV程度上のバンドが入るくらいで試すと良い）
-  - （制限RPAの場合）相互作用の大きさを求めたい軌道のワニエ基底を計算する。
+  - 相互作用の大きさを求めたい軌道のワニエ基底を計算する。
   - N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
   - Ecut_for_eps に関して計算を収束させる
@@ Line 228: / Line 228: @@
   - SCF計算でk点数を増やし、計算を収束させる
-有効相互作用を求める場合は、さらに Num_freq_grid（周波数のグリッド数）を増やして有効相互作用の値の収束を確認しましょう。
+  * 有効相互作用を求める場合は、さらに Num_freq_grid（周波数のグリッド数）を増やして有効相互作用の値の収束を確認しましょう。
+  * 計算が重いので、有効相互作用の有効数値は2桁出すのが精一杯かもしれません

雑多な記録

User Tools

Site Tools

Differences

Page Tools