アルミニウムのフェルミ面の計算を行う。 表示にはXCrySDenを使うので、事前にインストールすること。

Al.scf.in

&control
   prefix = 'Al'
   calculation = 'scf'
   outdir = './'
   pseudo_dir = './'
/
&system
   ibrav = 2
   celldm(1) = 7.50
   nat = 1
   ntyp = 1
   ecutwfc = 12.0
   occupations = 'smearing'
   smearing = 'mv'
   degauss = 0.03
/
&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とし、、以下のようにcalculation, occupations, K_POINTSを書き換える。 k点の数はなるべく多い方が正確なフェルミ面が描けるが、増やしすぎると時間が非常にかかる。 できあがったフェルミ面を見ながらk点数を徐々に増やしていくとよい。

Al.fs.in

&control
   prefix = 'Al'
   calculation = 'nscf'
   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
 16 16 16 0 0 0

$ pw.x < Al.fs.in > Al.fs.out

SCF計算の段階でk点数が十分取れるなら、SCF計算もoccupations = 'tetrahedra_opt'を使っても良い（k点数に対するFermi energyの収束を必ず確認）。

FS.inというファイルを新規に作成する

FS.in

&fermi
  outdir = './'
  prefix = 'Al'
/

$ fs.x < FS.in > FS.out

$ xcrysden --bxsf Al_fs.bxsf