Differences

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

--- quantumespresso:respack [2020/08/01 01:06]
koudai [インプットファイル]
+++ quantumespresso:respack [2020/08/01 03:13]
koudai [分極関数の計算]
@@ Line 145: / Line 145: @@
 </code>
   * 既約k点の情報が必要なので、SCF計算 (ネームリスト<nowiki>&control</nowiki>内で<nowiki>calculation='scf'</nowiki>とする計算) の直後にRESPACKを使った計算をする必要があります
-  * k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、prefix.save を削除してから再度SCF計算を実行してください
+  * k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、ディレクトリ prefix.save を削除してから再度SCF計算を実行してください
   * k点は automatic により生成し、シフトしていないもの (0 0 0) を使用してください
@@ Line 163: / Line 163: @@
 </code>
-  * エネルギーインナーウィンドウはできる限り指定しておきましょう。
+  * エネルギーウインドウはワニエ基底に射影するバンドの範囲です。
-  * 電子数は指定しても良いですが、うまくワニエ化できなくなることが多くなります。そのときは electron_number_wannier_space = 0.0 にしてください。
+  * エネルギーインナーウィンドウ内のすべてのバンドがワニエ基底によるバンドと一致するように計算されます。指定しなくても計算可能ですが、できるだけ広い範囲を指定することでより確実なワニエ基底を得ることができます。
 ===== 分極関数の計算 =====
@@ Line 176: / Line 176: @@
   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます
   * この計算に一番時間がかかります
-  * ユニットセルの体積が大きいほどメモリを多く消費します。
+  * バンドの数、カットオフ、k点数、ユニットセルの体積に比例して非常に多くのメモリを消費します。普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます（それ以上になる場合はスパコンの利用を検討してください）。
+  * 観測される物理量としての光学応答を調べたければ通常のRPAを、相互作用パラメータを求めたければ制限RPAを使用します。それぞれ収束に必要なパラメータは異なるので注意してください。
 ===== 相互作用の計算 =====
@@ Line 202: / Line 203: @@
   * Ecut_for_eps ... ecutwfcの1/10（デフォルト）
-次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です)
+次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットして収束を確認するのが便利です)
-  - 多めの数のnbndでSCF計算を実行する（フェルミ準位より50eV程度上のバンドが入るくらい）
+  - 多めの数のnbndでSCF計算を実行する（最初は非占有状態のバンドが50本くらいになるか、あるいはフェルミ準位より2,30eV程度上のバンドが入るくらいで試すと良い）
-  - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
+  - （制限RPAの場合）相互作用の大きさを求めたい軌道のワニエ基底を計算する。
-  - Ecut_for_eps を増やしていき、計算を収束させる。
+  - N_CALC_BAND に関して分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
+  - Ecut_for_eps に関して計算を収束させる。
   - SCF計算でk点数を増やし、分極関数を収束させる

雑多な記録

User Tools

Site Tools

Differences

Page Tools