quantumespresso:respack
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
quantumespresso:respack [2022/09/08 00:50] – [概要] koudai | quantumespresso:respack [2024/02/03 18:42] (current) – [概要] koudai | ||
---|---|---|---|
Line 11: | Line 11: | ||
* arXiv: https:// | * arXiv: https:// | ||
* 出版されたもの: | * 出版されたもの: | ||
- | * 計算は非常に重いので、研究室にあるような普通の計算機ですとユニットセルに原子が5, | + | * 計算は非常に重いので、普通のパソコンですとユニットセルに原子が5, |
* f電子系やスピン軌道相互作用がある系には非対応 | * f電子系やスピン軌道相互作用がある系には非対応 | ||
- | * 擬ポテンシャルはノルム保存型を使用する必要があります。例えば以下のサイトからダウンロードできます | + | * 擬ポテンシャルはノルム保存型を使用する必要があります。例えば以下のサイトから入手できます |
* http:// | * http:// | ||
* https:// | * https:// | ||
Line 148: | Line 148: | ||
たとえばPWscfの入力ファイルが prefix.scf.in だった場合、次のようにします。 | たとえばPWscfの入力ファイルが prefix.scf.in だった場合、次のようにします。 | ||
< | < | ||
- | $ mpirun -n 2 < prefix.scf.in > prefix.scf.out | + | $ mpirun -n 2 pw.x < prefix.scf.in > prefix.scf.out |
</ | </ | ||
* 既約k点の情報が必要なので、SCF計算 (ネームリスト< | * 既約k点の情報が必要なので、SCF計算 (ネームリスト< | ||
Line 216: | Line 216: | ||
* Ecut_for_eps ... 分極関数のカットオフ (Ry) | * Ecut_for_eps ... 分極関数のカットオフ (Ry) | ||
- | 最初は次の値を使用します | ||
- | * k点数 ... ワニエでバンドが再現できる最低限の数 | + | 有効相互作用を求める場合は、さらに Num_freq_grid(周波数のグリッド数)を増やして有効相互作用の値の収束を確認します。 |
- | * Ecut_for_eps ... 3, 4, 5 Ryあたりで計算して様子を見る | + | |
- | + | ||
- | 次の手順で収束させます | + | |
- | + | ||
- | - 多めの数のnbndでSCF計算を実行する(最初はフェルミ準位より50eV程度上のバンドが入るくらいで試すと良い) | + | |
- | - 相互作用の大きさを求めたい軌道のワニエ基底を計算する。 | + | |
- | - N_CALC_BAND に関して計算を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | + | |
- | - Ecut_for_eps に関して計算を収束させる | + | |
- | - N_CALC_BANDの収束をチェックする | + | |
- | - SCF計算でk点数を増やし、計算を収束させる | + | |
- | + | ||
- | + | ||
- | * 有効相互作用を求める場合は、さらに Num_freq_grid(周波数のグリッド数)を増やして有効相互作用の値の収束を確認しましょう。 | + | |
* 計算が重いので、有効相互作用の有効数値は2桁出すのが精一杯かもしれません | * 計算が重いので、有効相互作用の有効数値は2桁出すのが精一杯かもしれません | ||
quantumespresso/respack.1662565814.txt.gz · Last modified: 2022/09/08 00:50 by koudai