Differences

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

--- quantumespresso:フェルミ面 [2018/02/25 23:06] – [pw.xによるSCF計算] koudai
+++ quantumespresso:フェルミ面 [2021/07/12 23:39] (current) – koudai
@@ Line 1: / Line 1: @@
 ====== 概要 ======
-アルミニウムのフェルミ面の計算を行う。
+アルミニウムを例に金属のフェルミ面の計算を行う。
 表示にはXCrySDenを使うので、事前にインストールすること。
@@ Line 11: / Line 11: @@
    prefix = 'Al'
    calculation = 'scf'
-   outdir = './'
+   outdir = './tmp/'
    pseudo_dir = './'
 /
@@ Line 20: / Line 20: @@
    ntyp = 1
    ecutwfc = 12.0
-   occupations = 'smearing'
+   occupations = 'tetrahedra_opt'
-   smearing = 'mp'
-   degauss = 0.03
 /
 &electrons
@@ Line 38: / Line 36: @@
 ====== 固有エネルギーの計算 ======
-Al.scf.inをコピーしてAl.fs.inとし、、以下のようにcalculation, occupations, K_POINTSを書き換える。
+Al.scf.inをコピーしてAl.nscf.inとし、以下のように編集する。
-k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。
+対称点におけるエネルギーの値は重要なので、K_POINTSはシフトしていないものを使うこと。
-できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。
-<file - Al.fs.in>
+  * calculation = 'nscf'に変更
+  * K_POINTSの数を増やす
+<file - Al.nscf.in>
 &control
    prefix = 'Al'
    calculation = 'nscf'
-   outdir = './'
+   outdir = './tmp/'
    pseudo_dir = './'
 /
@@ Line 67: / Line 67: @@
 </file>
-  $ pw.x < Al.fs.in > Al.fs.out
+k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると計算時間が非常にかかる。
+できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。
-SCF計算の段階でk点数が十分取れるなら、SCF計算もoccupations = 'tetrahedra_opt'を使っても良い（k点数に対するFermi energyの収束を必ず確認）。
+  $ pw.x < Al.nscf.in > Al.nscf.out
@@ Line 75: / Line 77: @@
-FS.inというファイルを新規に作成する
+Al.fs.inというファイルを新規に作成する
-<file - FS.in>
+<file - Al.fs.in>
 &fermi
-  outdir = './'
+  outdir = './tmp/'
   prefix = 'Al'
 /
 </file>
-  $ fs.x < FS.in > FS.out
+  $ fs.x < Al.fs.in > Al.fs.out