Differences

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

--- wien2k:フェルミ面 [2021/04/28 21:28]
koudai [フェルミ面の描き方]
+++ wien2k:フェルミ面 [2021/04/29 21:00]
koudai [軌道ポテンシャルを入れている場合]
@@ Line 25: / Line 25: @@
     * それぞれのバンドのフェルミ面を見ることができる。全部を一気に表示したかったら[Merged Bands]を選ぶ
     * できたフェルミ面を保存したい場合はフェルミ面の表示領域で右クリックして[File]->[Save Fermi Surface(s) in BXSF format]
-    * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまう？
+      * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまうので（ファイル名はcase.bxsf.band-xとなる。xは対応するバンドのインデックス）、下のような感じになるように自分でファイルを開いて結合してください<file - case.bxsf>
+ 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
+</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>
 ===== 注意 =====

雑多な記録

User Tools

Site Tools

Differences

Page Tools