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/13 13:13]
koudai [概要]
quantumespresso:respack [2020/08/14 15:01]
koudai [分極関数の計算]
Line 8: Line 8:
     * https://www.cond-mat.de/events/correl11/      * https://www.cond-mat.de/events/correl11/ 
   * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います   * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います
-    https://sites.google.com/view/kazuma7k6r+  メモリを非常に多く使うため、普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます。スパコンを使えば20個くらいは行けるのではないかと思います 
   * f電子系やスピン軌道相互作用がある系には非対応   * f電子系やスピン軌道相互作用がある系には非対応
   * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます   * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます
Line 45: Line 45:
 $ sudo make install $ sudo make install
 </code> </code>
-    * makeにやや時間がかかります 
   - パスをとおします<code>   - パスをとおします<code>
 $ echo 'export PATH=$PATH:/usr/local/respack/bin/' >> ~/.bashrc $ echo 'export PATH=$PATH:/usr/local/respack/bin/' >> ~/.bashrc
Line 58: Line 57:
   - QuantumESPRESSO の pw.x によるSCF計算   - QuantumESPRESSO の pw.x によるSCF計算
   - RESPACKの calc_wannier による最局在ワニエ関数の計算。事前にSCF計算が必要   - RESPACKの calc_wannier による最局在ワニエ関数の計算。事前にSCF計算が必要
-  - RESPACKの calc_chiqw による誘電関数の計算。事前にワニエ関数の計算が必要+  - RESPACKの calc_chiqw による誘電関数の計算。cRPAを使う場合は事前にワニエ関数の計算が必要
   - RESPACKの calc_w3d による直接相互作用の計算。事前にcalc_chiqwのcRPA計算が必要   - RESPACKの calc_w3d による直接相互作用の計算。事前にcalc_chiqwのcRPA計算が必要
   - RESPACKの calc_j3d による交換相互作用の計算。事前にcalc_chiqwのcRPA計算が必要   - RESPACKの calc_j3d による交換相互作用の計算。事前にcalc_chiqwのcRPA計算が必要
Line 117: Line 116:
 &param_chiqw &param_chiqw
   Ecut_for_eps           = 3.6d0,   !o: 分極関数のカットオフ (Ry) 記述がなければSCF計算のecutwfcの値の1/10になる   Ecut_for_eps           = 3.6d0,   !o: 分極関数のカットオフ (Ry) 記述がなければSCF計算のecutwfcの値の1/10になる
-  Num_freq_grid          = 70,      !o: 計算周波数の総数+  Num_freq_grid          = 70,      !o: 計算周波数の総数(メッシュは log ω に対して切られる)
   N_CALC_BAND            = 50,      !o: 分極計算で考慮されるバンドの総数。SCF計算のnbndの数よりも小さくする。記述がなければSCF計算のnbndの値になる   N_CALC_BAND            = 50,      !o: 分極計算で考慮されるバンドの総数。SCF計算のnbndの数よりも小さくする。記述がなければSCF計算のnbndの値になる
   MPI_num_proc_per_qcomm = 2,       !o: コミュニティ当たりのプロセス数。MPIを使用する際に指定   MPI_num_proc_per_qcomm = 2,       !o: コミュニティ当たりのプロセス数。MPIを使用する際に指定
Line 188: Line 187:
   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます
   * この計算に一番時間がかかります   * この計算に一番時間がかかります
-  * バンドの数、カットオフ、k点数、ユニットセルの体積に比例して非常に多くのメモリを消費します。普通の計算機ですとユニットセルに原子が10個くらいの物質が限界かと思われます(それ以上になる場合はスパコンの利用を検討してください)。+  * バンドの数、エネルギーカットオフ、k点数、ユニットセルの体積に比例して非常に多くのメモリを消費します
   * 観測される物理量としての光学応答を調べたければ通常のRPAを、相互作用パラメータを求めたければ制限RPAを使用します。それぞれ収束に必要なパラメータは異なるので注意してください。   * 観測される物理量としての光学応答を調べたければ通常のRPAを、相互作用パラメータを求めたければ制限RPAを使用します。それぞれ収束に必要なパラメータは異なるので注意してください。
 +    * 一般に、通常のRPAの方が計算コストが大きいです 
  
 ===== 相互作用の計算 ===== ===== 相互作用の計算 =====
Line 243: Line 243:
 omega_max=200 omega_max=200
  
 +set logscale x
  
  
quantumespresso/respack.txt · Last modified: 2024/02/03 18:42 by koudai