Differences

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

--- quantumespresso:フェルミ面 [2018/02/25 22:44]
koudai 作成
+++ quantumespresso:フェルミ面 [2021/04/27 20:50]
koudai [フェルミ面の計算]
@@ Line 1: / Line 1: @@
 ====== 概要 ======
-アルミニウムのフェルミ面の計算を行う。
+アルミニウムを例に金属のフェルミ面の計算を行う。
 表示にはXCrySDenを使うので、事前にインストールすること。
@@ Line 11: / Line 11: @@
    prefix = 'Al'
    calculation = 'scf'
-   outdir = './'
+   outdir = './tmp/'
    pseudo_dir = './'
 /
@@ Line 33: / Line 33: @@
   $ pw.x < Al.scf.in > Al.scf.out
 ====== 固有エネルギーの計算 ======
-Al.scf.inをコピーしてAl.fs.inとし、、以下の部分だけ書き換える。
+Al.scf.inをコピーしてAl.nscf.inとし、以下のようにcalculationとK_POINTSを書き換える。
-k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間がかなりかかるので、最初は4*4*4くらいから始めて、できあがったフェルミ面を見ながら徐々に増やしていくとよい。
+K_POINTSはシフトしていないものを使うこと。
-経験的には16*16*16くらい必要だが、物質によってはかなり時間がかかることも。
-<file - Al.fs.in>
+k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。
+できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。
-（略）
+<file - Al.nscf.in>
-calculation = 'nscf'
+&control
-（略）
+   prefix = 'Al'
+   calculation = 'nscf'
+   outdir = './tmp/'
+   pseudo_dir = './'
+/
+&system
+   ibrav = 2
+   celldm(1) = 7.50
+   nat = 1
+   ntyp = 1
+   ecutwfc = 12.0
+   occupations = 'tetrahedra_opt'
+/
+&electrons
+/
+ATOMIC_SPECIES
+ Al  26.98  Al.pz-vbc.UPF
+ATOMIC_POSITIONS crystal
+ Al  0.0  0.0  0.0
 K_POINTS automatic
 16 16 0 0 0
 </file>
-  $ pw.x < Al.fs.in > Al.fs.out
+  $ pw.x < Al.nscf.in > Al.nscf.out
@@ Line 56: / Line 74: @@
-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

雑多な記録

User Tools

Site Tools

Differences

Page Tools