SCF計算を事前に済ませておいてください。
$ xcrysden --wien_fermisurface /
BEGIN_INFO Fermi Energy: (フェルミ準位) END_INFO BEGIN_BLOCK_BANDGRID_3D band_energies BANDGRID_3D_BANDS (バンドの数) (k点の数) (基準点) (逆格子ベクトル1) (逆格子ベクトル2) (逆格子ベクトル3) BAND: (最初のバンド) (データ) BAND: (次のバンド) (データ) ・・・ BAND: (最後のバンド) (データ) END_BANDGRID_3D END_BLOCK_BANDGRID_3D
$ xcrysden --bxsf case.bxsf
XcrysdenからWIEN2kを呼び出すときに -orb オプションを入れられない。 なのでGGA+Uをしている場合は先にlapw1, lapw2を計算しておく
$ x kgen # k点数を指定 $ x lapw1 -up -orb $ x lapw1 -dn -orb $ x lapw2 -up -fermi -orb # アップスピンのフェルミエネルギーを計算 $ x lapw2 -dn -fermi -orb # アップとダウンでフェルミ面が異なる場合 $ xcrysden --wien_fermisurface / # Render Fermi Surfaceを選択