quantumespresso:respack
Differences
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:53] koudai [分極関数の計算] |
quantumespresso:respack [2020/08/01 03:04] 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計算を実行する（最初は非占有状態のバンドが50本くらいになるか、あるいはフェルミ準位より数10eV程度上のバンドが入るくらいで試すと良い） |
| + | - （制限RPAの場合）相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | ||
| - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | ||
| - Ecut_for_eps を増やしていき、計算を収束させる。 | - Ecut_for_eps を増やしていき、計算を収束させる。 | ||
quantumespresso/respack.txt · Last modified: 2024/02/03 18:42 by koudai
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
