• フェルミ面を描きます。
• XCrySDenを使います。

SCF計算を事前に済ませておいてください。

1. XCrySDenを起動

   $ xcrysden --wien_fermisurface /

   • --wien_fermisurfaceの後ろには作業ディレクトリのパスが入ります。上は現在いるディレクトリと同じ場合。
   • Xcrysdenを起動して [File]→[Open WIEN2k]→[Fermi Surface] で、計算を行っているディレクトリを選択するのでも可
2. Number of kpointsを多めにとる（フェルミ面がきれいになるまで徐々に増やしていけば良い）
3. [Generate k-mesh]→[Done]
4. [Calculate Eigenvalues [x lapw1 ]]
   • 並列化計算をする場合はParallel execution of lapw1/lapw2にチェックを入れる
   • k点が多いと時間がかかる
5. [Calculate Fermi energy [x lapw2 -fermi]]
6. [Render Fermi surface]
   • Specify the Fermi energyはそのままにして[OK]
   • BarGraphやBandwidthから、select bandsのフェルミ準位にかかっているバンドを選択して[Selected]
   • それぞれのバンドのフェルミ面を見ることができる。全部を一気に表示したかったら[Merged Bands]を選ぶ
   • できたフェルミ面を保存したい場合はフェルミ面の表示領域で右クリックして[File]→[Save Fermi Surface(s) in BXSF format]
   • 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまう？
   • 再び見るときは次のコマンドで開く

     $ xcrysden --bxsf case.bxsf.band-x

• XcrysDenのソースコードにk点数の上限が定められている
• もし文句を言われたらXCrysDenのソースコードから wn_readbakgen.f を開いて MAX_NKPT（全ブリルアンゾーンのk点数の上限）と MAX_IRNKPT（既約ブリルアンゾーンのk点数の上限）の値を増やしてコンパイルしなおせばよい。