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/09/19 19:40] koudai [計算の収束] |
quantumespresso:respack [2022/09/08 00:50] koudai [概要] |
||
---|---|---|---|
Line 5: | Line 5: | ||
* 誘電関数から求まる動的物理量(電子エネルギー欠損スペクトル、光学伝導度、反射率) | * 誘電関数から求まる動的物理量(電子エネルギー欠損スペクトル、光学伝導度、反射率) | ||
* 制限乱雑位相近似 (constraint random-phase approximation, | * 制限乱雑位相近似 (constraint random-phase approximation, | ||
- | * cRPAの理論はこちらのテキストの7章を参照してください | + | * cRPAの理論はこちらのテキストの7章がわかりやすいです |
* https:// | * https:// | ||
* 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | ||
- | | + | * 最新版 (20200113) の使用方法は論文を参照してください(いろいろ変更点があるので、昔のマニュアルを使うと混乱が生じるかもしれません) |
+ | * arXiv: https:// | ||
+ | * 出版されたもの: | ||
+ | | ||
* f電子系やスピン軌道相互作用がある系には非対応 | * f電子系やスピン軌道相互作用がある系には非対応 | ||
- | * 擬ポテンシャルはノルム保存型のものにのみ対応しています。例えば以下のサイトからダウンロードできます | + | * 擬ポテンシャルはノルム保存型を使用する必要があります。例えば以下のサイトからダウンロードできます |
* http:// | * http:// | ||
* https:// | * https:// | ||
Line 115: | Line 118: | ||
1 2 3 4 5 6 7 8 9 !o: N_write_wannierを指定したときの、計算したいワニエ軌道の番号 | 1 2 3 4 5 6 7 8 9 !o: N_write_wannierを指定したときの、計算したいワニエ軌道の番号 | ||
& | & | ||
- | Ecut_for_eps | + | Ecut_for_eps |
Num_freq_grid | Num_freq_grid | ||
N_CALC_BAND | N_CALC_BAND | ||
Line 126: | Line 129: | ||
ttrhdrn_dmna | ttrhdrn_dmna | ||
ttrhdrn_dmnr | ttrhdrn_dmnr | ||
- | flg_cRPA | + | flg_cRPA |
flg_calc_type | flg_calc_type | ||
& | & | ||
Line 145: | Line 148: | ||
たとえばPWscfの入力ファイルが prefix.scf.in だった場合、次のようにします。 | たとえばPWscfの入力ファイルが prefix.scf.in だった場合、次のようにします。 | ||
< | < | ||
- | $ mpirun -n 24 < prefix.scf.in > prefix.scf.out | + | $ mpirun -n 2 < prefix.scf.in > prefix.scf.out |
</ | </ | ||
* 既約k点の情報が必要なので、SCF計算 (ネームリスト< | * 既約k点の情報が必要なので、SCF計算 (ネームリスト< | ||
Line 169: | Line 172: | ||
* OpenMPによる並列化のみ対応しています | * OpenMPによる並列化のみ対応しています | ||
- | * エネルギーウインドウはワニエ基底に射影するバンドの範囲です。 | + | * エネルギーウインドウはワニエ基底に射影するバンドの範囲です。射影しない軌道のバンドが少し混じっていても構いません。 |
- | * エネルギーインナーウィンドウ内のすべてのバンドがワニエ基底によるバンドと一致するように計算されます。指定しなくても計算可能ですが、できるだけ広い範囲を指定することでより確実なワニエ基底を得ることができます。 | + | * 特にエネルギーインナーウィンドウの範囲を指定すると、その中にあるすべてのバンドがワニエ基底によるバンドと一致するように計算されます。 |
===== 分極関数の計算 ===== | ===== 分極関数の計算 ===== | ||
Line 210: | Line 212: | ||
分極関数の計算で収束させなければならないパラメータは次の3つです | 分極関数の計算で収束させなければならないパラメータは次の3つです | ||
- | * k点数 | + | * k点数(SCF計算) |
* N_CALC_BAND ... バンドの数 | * N_CALC_BAND ... バンドの数 | ||
* Ecut_for_eps ... 分極関数のカットオフ (Ry) | * Ecut_for_eps ... 分極関数のカットオフ (Ry) | ||
Line 217: | Line 219: | ||
* k点数 ... ワニエでバンドが再現できる最低限の数 | * k点数 ... ワニエでバンドが再現できる最低限の数 | ||
- | * Ecut_for_eps ... 2, 3, 4 Ryあたりで計算して様子を見る | + | * Ecut_for_eps ... 3, 4, 5 Ryあたりで計算して様子を見る |
次の手順で収束させます | 次の手順で収束させます | ||
- | - 多めの数のnbndでSCF計算を実行する(最初はフェルミ準位より40~50eV程度上のバンドが入るくらいで試すと良い) | + | - 多めの数のnbndでSCF計算を実行する(最初はフェルミ準位より50eV程度上のバンドが入るくらいで試すと良い) |
- | - (制限RPAの場合)相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | + | - 相互作用の大きさを求めたい軌道のワニエ基底を計算する。 |
- N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | - N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | ||
- Ecut_for_eps に関して計算を収束させる | - Ecut_for_eps に関して計算を収束させる |