Differences

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

--- quantumespresso:respack [2020/07/28 18:38]
koudai [QuantumESPRESSOによるSCF計算]
+++ quantumespresso:respack [2020/08/01 02:51]
koudai [分極関数の計算]
@@ Line 8: / Line 8: @@
   * 日本語および英語のマニュアルでは使用例とその結果が充実しているので、初心者でも特に困ることはないと思います
     * https://sites.google.com/view/kazuma7k6r
-  * f電子系はできないようです？
+  * f電子系やスピン軌道相互作用がある系には非対応。
@@ Line 88: / Line 88: @@
   py  0.2d0 0.0d0 0.0d0 0.0d0             !n: py
   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
   dzx 0.2d0 0.0d0 0.0d0 0.0d0             !n: dzx
-  dx2 0.2d0 0.0d0 0.0d0 0.0d0             !n: dx2-y2
+  dyz 0.2d0 0.0d0 0.0d0 0.0d0             !n: dyz
+  dx2 0.2d0 0.0d0 0.0d0 0.0d0             !n: dx2-y2
+  dxy 0.2d0 0.0d0 0.0d0 0.0d0             !n: dx
 &param_interpolation
-  N_sym_points =5,  !n: 計算ラインを構成する対称 k 点数
+  N_sym_points = 5,  !n: 計算ラインを構成する対称k点数
-  Ndiv         =40/ !n: 対称 k 点間の分割数
+  Ndiv         = 40/ !n: 対称k点間の分割数
-.500 0.500 0.500 !n: 対称 k 点 ; SK_sym_pts(1:3,1); L（逆格子ベクトルの分率座標）
+.500 0.500 0.500 !n: 対称k点（逆格子ベクトルの分率座標）
-.000 0.000 0.000 !n: 対称 k 点 ; SK_sym_pts(1:3,2); G
+.000 0.000 0.000 !n: 対称k点
-.500 0.000 0.500 !n: 対称 k 点 ; SK_sym_pts(1:3,3); X
+.500 0.000 0.500 !n: 対称k点
-.500 0.250 0.750 !n: 対称 k 点 ; SK_sym_pts(1:3,4); W
+.500 0.250 0.750 !n: 対称k点
-.500 0.500 0.500 !n: 対称 k 点 ; SK_sym_pts(1:3,5); L
+.500 0.500 0.500 !n: 対称k点
 &param_visualization
   flg_vis_wannier = 0, !o: 実空間ワニエ関数を計算 (しない: 0, する: 1)
@@ Line 145: / Line 145: @@
 </code>
   * 既約k点の情報が必要なので、SCF計算 (ネームリスト<nowiki>&control</nowiki>内で<nowiki>calculation='scf'</nowiki>とする計算) の直後にRESPACKを使った計算をする必要があります
-  * k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、prefix.save を削除してから再度SCF計算を実行してください
+  * k点数を変更したり、バンド計算など別の計算を行った後にRESPACKを使用する場合は、ディレクトリ prefix.save を削除してから再度SCF計算を実行してください
   * k点は automatic により生成し、シフトしていないもの (0 0 0) を使用してください
@@ Line 163: / Line 163: @@
 </code>
+  * エネルギーウインドウはワニエ基底に射影するバンドの範囲です。
+  * エネルギーインナーウィンドウ内のすべてのバンドがワニエ基底によるバンドと一致するように計算されます。指定しなくても計算可能ですが、できるだけ広い範囲を指定することでより確実なワニエ基底を得ることができます。
 ===== 分極関数の計算 =====
@@ Line 174: / Line 176: @@
   * MPIのプロセス数は MPI_num_proc_per_qcomm * MPI_num_qcomm に一致させます
   * この計算に一番時間がかかります
+  * バンドの数やカットオフ、ユニットセルの体積に比例して非常に多くのメモリを消費します。研究室にあるような普通の計算機ですとユニットセルに原子が10数個くらいの物質が限界かと思われます。
 ===== 相互作用の計算 =====
@@ Line 201: / Line 204: @@
 次の手順で収束させます(flg_calc_type=1としてEELSの結果をプロットするのが便利です)
-  - 多めの数のnbndでSCF計算を実行する（フェルミ準位より50eV程度上のバンドが入るくらいくらい）
+  - 多めの数のnbndでSCF計算を実行する（フェルミ準位より50eV程度上のバンドが入るくらい）
   - N_CALC_BANDを増やし、分極関数を収束させる。足りなくなったらnbndを増やしたSCF計算を再び行う
   - Ecut_for_eps を増やしていき、計算を収束させる。

雑多な記録

User Tools

Site Tools

Differences

Page Tools