Differences

This shows you the differences between two versions of the page.

--- wien2k:フェルミ面 [2021/04/28 22:05] – [フェルミ面の描き方] koudai
+++ wien2k:フェルミ面 [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 25: / Line 25: @@
     * それぞれのバンドのフェルミ面を見ることができる。全部を一気に表示したかったら[Merged Bands]を選ぶ
     * できたフェルミ面を保存したい場合はフェルミ面の表示領域で右クリックして[File]->[Save Fermi Surface(s) in BXSF format]
-      * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまうので、下のような感じになるように自分でファイルを開いて結合してください＞<file case.bxsf>
+      * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまうので（ファイル名はcase.bxsf.band-xとなる。xは対応するバンドのインデックス）、下のような感じになるように自分でファイルを開いて結合してください<file - case.bxsf>
  BEGIN_INFO
    Fermi Energy:     （フェルミ準位）
@@ Line 38: / Line 38: @@
  （逆格子ベクトル2）
  （逆格子ベクトル3）
- BAND:    1
+ BAND:    （最初のバンド）
-（なんかデータ）
+（データ）
- BAND:    2
+ BAND:    （次のバンド）
-（なんかデータ）
+（データ）
 ・・・
- BAND:    x
+ BAND:    （最後のバンド）
-（なんかデータ）
+（データ）
  END_BANDGRID_3D
  END_BLOCK_BANDGRID_3D
-</code>
+</file>
     * 再び見るときは次のコマンドで開く<code>
-$ xcrysden --bxsf case.bxsf.band-x</code>
+$ xcrysden --bxsf case.bxsf</code>
+===== 軌道ポテンシャルを入れている場合 =====
+XcrysdenからWIEN2kを呼び出すときに -orb オプションを入れられない。
+なのでGGA+Uをしている場合は先にlapw1, lapw2を計算しておく
+<code>
+$ 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を選択
+</code>
 ===== 注意 =====