This is an old revision of the document!
アルミニウムのフェルミ面の計算を行う。 表示にはXCrySDenを使うので、事前にインストールすること。
&control prefix = 'Al' calculation = 'scf' outdir = './' 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 6 6 6 1 1 1
$ pw.x < Al.scf.in > Al.scf.out
Al.scf.inをコピーしてAl.fs.inとし、、以下の部分だけ書き換える。 k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間がかなりかかるので、最初は4*4*4くらいから始めて、できあがったフェルミ面を見ながら徐々に増やしていくとよい。 経験的には16*16*16くらい必要だが、物質によってはかなり時間がかかることも。
(略) calculation = 'nscf' (略) K_POINTS automatic 16 16 16 0 0 0
$ pw.x < Al.fs.in > Al.fs.out
$ xcrysden --bxsf Al_fs.bxsf