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/07/28 03:01] koudai [インプットファイル] |
quantumespresso:respack [2020/08/01 01:02] koudai [分極関数の計算] |
||
---|---|---|---|
Line 8: | Line 8: | ||
* 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います | ||
* https:// | * https:// | ||
- | * f電子系はできないようです? | + | * f電子系やスピン軌道相互作用がある系には非対応。 |
Line 76: | Line 76: | ||
DAMP = 0.1d0, | DAMP = 0.1d0, | ||
MAX_STEP_LENGTH | MAX_STEP_LENGTH | ||
- | set_inner_window | + | set_inner_window |
Lower_inner_window | Lower_inner_window | ||
Upper_inner_window | Upper_inner_window | ||
Line 88: | Line 88: | ||
py 0.2d0 0.0d0 0.0d0 0.0d0 !n: py | py 0.2d0 0.0d0 0.0d0 0.0d0 !n: py | ||
pz 0.2d0 0.0d0 0.0d0 0.0d0 !n: pz | pz 0.2d0 0.0d0 0.0d0 0.0d0 !n: pz | ||
- | dxy 0.2d0 0.0d0 0.0d0 0.0d0 !n: dx | ||
- | dyz 0.2d0 0.0d0 0.0d0 0.0d0 !n: dyz | ||
dz2 0.2d0 0.0d0 0.0d0 0.0d0 !n: dz2 | dz2 0.2d0 0.0d0 0.0d0 0.0d0 !n: dz2 | ||
dzx 0.2d0 0.0d0 0.0d0 0.0d0 !n: dzx | dzx 0.2d0 0.0d0 0.0d0 0.0d0 !n: dzx | ||
- | dx2 0.2d0 0.0d0 0.0d0 0.0d0 !n: dx2-y2 | + | |
+ | | ||
+ | dxy 0.2d0 0.0d0 0.0d0 0.0d0 !n: dx | ||
& | & | ||
N_sym_points =5, !n: 計算ラインを構成する対称 k 点数 | N_sym_points =5, !n: 計算ラインを構成する対称 k 点数 | ||
Line 126: | Line 126: | ||
flg_calc_type | flg_calc_type | ||
& | & | ||
- | calc_ifreq | + | calc_ifreq |
ix_intJ_min = 0, !o: 交換積分の計算範囲 (a1 方向始点格子 )(calc_j3d のみ ) | ix_intJ_min = 0, !o: 交換積分の計算範囲 (a1 方向始点格子 )(calc_j3d のみ ) | ||
ix_intJ_max = 0, !o: 交換積分の計算範囲 (a1 方向終点格子 )(calc_j3d のみ ) | ix_intJ_max = 0, !o: 交換積分の計算範囲 (a1 方向終点格子 )(calc_j3d のみ ) | ||
Line 146: | Line 146: | ||
* 既約k点の情報が必要なので、SCF計算 (ネームリスト< | * 既約k点の情報が必要なので、SCF計算 (ネームリスト< | ||
* k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、prefix.save を削除してから再度SCF計算を実行してください | * k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、prefix.save を削除してから再度SCF計算を実行してください | ||
+ | * k点は automatic により生成し、シフトしていないもの (0 0 0) を使用してください | ||
< | < | ||
Line 162: | Line 163: | ||
</ | </ | ||
+ | * エネルギーインナーウィンドウはできる限り指定しておきましょう。 | ||
+ | * 電子数は指定しても良いですが、うまくワニエ化できなくなることが多くなります。そのときは electron_number_wannier_space = 0.0 にしてください。 | ||
===== 分極関数の計算 ===== | ===== 分極関数の計算 ===== | ||
Line 173: | Line 176: | ||
* MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます | * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます | ||
* この計算に一番時間がかかります | * この計算に一番時間がかかります | ||
+ | * ユニットセルの体積が大きいほどメモリを多く消費します。 | ||
===== 相互作用の計算 ===== | ===== 相互作用の計算 ===== | ||
Line 200: | Line 204: | ||
次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です) | 次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です) | ||
- | - 多めの数のnbndでSCF計算を実行する(フェルミ準位より50eV程度上のバンドが入るくらいくらい) | + | - 多めの数のnbndでSCF計算を実行する(フェルミ準位より50eV程度上のバンドが入るくらい) |
- N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う | ||
- Ecut_for_eps を増やしていき、計算を収束させる。 | - Ecut_for_eps を増やしていき、計算を収束させる。 |