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/01 02:51] koudai [分極関数の計算] |
quantumespresso:respack [2020/08/01 03:00] koudai [計算の収束] |
||
---|---|---|---|
Line 176: | Line 176: | ||
* MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます | * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます | ||
* この計算に一番時間がかかります | * この計算に一番時間がかかります | ||
- | * バンドの数やカットオフ、ユニットセルの体積に比例して非常に多くのメモリを消費します。研究室にあるような普通の計算機ですとユニットセルに原子が10数個くらいの物質が限界かと思われます。 | + | * バンドの数、カットオフ、k点数、ユニットセルの体積に比例して非常に多くのメモリを消費します。普通の計算機ですとユニットセルに原子が10数個くらいの物質が限界かと思われます。 |
+ | * 観測される物理量としての光学応答を調べたければ通常のRPAを、相互作用パラメータを求めたければ制限RPAを使用します。それぞれ収束に必要なパラメータは異なるので注意してください。 | ||
===== 相互作用の計算 ===== | ===== 相互作用の計算 ===== | ||
Line 204: | Line 205: | ||
次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です) | 次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です) | ||
- | - 多めの数のnbndでSCF計算を実行する(フェルミ準位より50eV程度上のバンドが入るくらい) | + | - 多めの数のnbndでSCF計算と、制限RPAの場合はワニエ基底を実行する(フェルミ準位より50eV程度上のバンドが入るくらい) |
- N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | ||
- Ecut_for_eps を増やしていき、計算を収束させる。 | - Ecut_for_eps を増やしていき、計算を収束させる。 |