Differences

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

--- wien2k:フェルミ面 [2021/04/28 21:21] – [フェルミ面の描き方] koudai
+++ wien2k:フェルミ面 [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 10: / Line 10: @@
   - XCrySDenを起動<code>
-$ xcrysden --wien_fermisurface
+$ xcrysden --wien_fermisurface /
 </code>
-  - [File]→[Open WIEN2k]→[Fermi Surface] で、計算を行っているディレクトリを選択
+    * <nowiki>--wien_fermisurface</nowiki>の後ろには作業ディレクトリのパスが入ります。上は現在いるディレクトリと同じ場合。
+    * Xcrysdenを起動して [File]→[Open WIEN2k]→[Fermi Surface] で、計算を行っているディレクトリを選択するのでも可
   - Number of kpointsを多めにとる（フェルミ面がきれいになるまで徐々に増やしていけば良い）
   - [Generate k-mesh]→[Done]
@@ Line 24: / 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>
 ===== 注意 =====