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:フェルミ面 [2018/02/25 23:03] koudai [固有エネルギーの計算] |
quantumespresso:フェルミ面 [2021/04/27 20:50] koudai [固有エネルギーの計算] |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== 概要 ====== | ====== 概要 ====== | ||
- | アルミニウムのフェルミ面の計算を行う。 | + | アルミニウムを例に金属のフェルミ面の計算を行う。 |
表示にはXCrySDenを使うので、事前にインストールすること。 | 表示にはXCrySDenを使うので、事前にインストールすること。 | ||
Line 11: | Line 11: | ||
| | ||
| | ||
- | | + | |
| | ||
/ | / | ||
Line 20: | Line 20: | ||
ntyp = 1 | ntyp = 1 | ||
| | ||
- | | + | |
- | | + | |
- | | + | |
/ | / | ||
& | & | ||
Line 38: | Line 36: | ||
====== 固有エネルギーの計算 ====== | ====== 固有エネルギーの計算 ====== | ||
- | Al.scf.inをコピーしてAl.fs.inとし、、以下の部分だけ書き換える。 | + | Al.scf.inをコピーしてAl.nscf.inとし、以下のようにcalculationとK_POINTSを書き換える。 |
+ | K_POINTSはシフトしていないものを使うこと。 | ||
k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。 | k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。 | ||
できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。 | できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。 | ||
- | <file - Al.fs.in> | + | <file - Al.nscf.in> |
- | (略) | + | & |
- | calculation = ' | + | prefix = ' |
- | (略) | + | calculation = ' |
- | occupations = ' | + | outdir = ' |
- | (略) | + | pseudo_dir = ' |
+ | / | ||
+ | & | ||
+ | ibrav = 2 | ||
+ | | ||
+ | nat = 1 | ||
+ | ntyp = 1 | ||
+ | | ||
+ | occupations = ' | ||
+ | / | ||
+ | & | ||
+ | / | ||
+ | ATOMIC_SPECIES | ||
+ | | ||
+ | ATOMIC_POSITIONS crystal | ||
+ | | ||
K_POINTS automatic | K_POINTS automatic | ||
16 16 16 0 0 0 | 16 16 16 0 0 0 | ||
</ | </ | ||
- | $ pw.x < Al.fs.in > Al.fs.out | + | $ pw.x < Al.nscf.in > Al.nscf.out |
- | + | ||
- | SCF計算の段階でk点数が十分取れるなら、SCF計算もoccupations = ' | + | |
Line 64: | Line 78: | ||
<file - FS.in> | <file - FS.in> | ||
&fermi | &fermi | ||
- | outdir = ' | + | outdir = './tmp/' |
prefix = ' | prefix = ' | ||
/ | / |