User Tools

Site Tools


quantumespresso:respack

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
quantumespresso:respack [2020/08/01 02:44]
koudai [QuantumESPRESSOによるSCF計算]
quantumespresso:respack [2020/08/01 03:13]
koudai [分極関数の計算]
Line 176: Line 176:
   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます
   * この計算に一番時間がかかります   * この計算に一番時間がかかります
-  * ユニットセルの体積が大きいほどメモリを多く消費します。+  * バンドの数、カットオフ、k点数、ユニットセルの体積に比例して非常に多くのメモリを消費します。普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます(それ以上になる場合はスパコンの利用を検討してください)。 
 +  * 観測される物理量としての光学応答を調べたければ通常のRPAを、相互作用パラメータを求めたければ制限RPAを使用します。それぞれ収束に必要なパラメータは異なるので注意してください
  
 ===== 相互作用の計算 ===== ===== 相互作用の計算 =====
Line 202: Line 203:
   * Ecut_for_eps ... ecutwfcの1/10(デフォルト)   * 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点数を増やし、分極関数を収束させる   - SCF計算でk点数を増やし、分極関数を収束させる
  
quantumespresso/respack.txt · Last modified: 2024/02/03 18:42 by koudai