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 22:48] koudai [pw.xによるSCF計算] |
quantumespresso:フェルミ面 [2021/04/27 20:50] koudai [フェルミ面の計算] |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== 概要 ====== | ====== 概要 ====== | ||
- | アルミニウムのフェルミ面の計算を行う。 | + | アルミニウムを例に金属のフェルミ面の計算を行う。 |
表示にはXCrySDenを使うので、事前にインストールすること。 | 表示にはXCrySDenを使うので、事前にインストールすること。 | ||
Line 11: | Line 11: | ||
| | ||
| | ||
- | | + | |
| | ||
/ | / | ||
Line 29: | Line 29: | ||
| | ||
K_POINTS automatic | K_POINTS automatic | ||
- | 8 8 8 0 0 0 | + | 6 6 6 1 1 1 |
</ | </ | ||
$ pw.x < Al.scf.in > Al.scf.out | $ 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はシフトしていないものを使うこと。 |
- | <file - Al.fs.in> | + | k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。 |
- | (略) | + | できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。 |
- | calculation = ' | + | |
- | (略) | + | <file - Al.nscf.in> |
+ | & | ||
+ | prefix = ' | ||
+ | calculation = ' | ||
+ | outdir = ' | ||
+ | | ||
+ | / | ||
+ | & | ||
+ | ibrav = 2 | ||
+ | | ||
+ | nat = 1 | ||
+ | ntyp = 1 | ||
+ | | ||
+ | | ||
+ | / | ||
+ | & | ||
+ | / | ||
+ | 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 |
+ | |||
Line 54: | Line 74: | ||
- | FS.inというファイルを新規に作成する | + | Al.fs.inというファイルを新規に作成する |
- | <file - FS.in> | + | <file - Al.fs.in> |
&fermi | &fermi | ||
- | outdir = ' | + | outdir = './tmp/' |
prefix = ' | prefix = ' | ||
/ | / | ||
</ | </ | ||
- | $ fs.x < FS.in > FS.out | + | $ fs.x < Al.fs.in > Al.fs.out |