Differences

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

--- wien2k:バンド計算 [2020/06/22 19:02]
koudai [固有エネルギーの計算]
+++ wien2k:バンド計算 [2021/06/27 22:04] (current)
@@ Line 24: / Line 24: @@
 ===== 固有エネルギーの計算 =====
-指定した経路上の各k点でのエネルギーを求めます。
+  - 指定した経路上の各k点でのエネルギーを求めます。<code>
-<code>
 $ x lapw1 -band -p
 </code>
+  -（各バンドにおける軌道ごとの重みを見たい場合のみ）qtlを実行します<code>
-（各バンドにおける軌道ごとの重みを見たい場合のみ）qtlを実行します
-<code>
 $ x qtl -p
 </code>
+  - 実際にバンドを描いてみます。<code>
-実際にバンドを描いてみます。
-<code>
 $ x spaghetti -p
-</code>
+</code>を実行すると、バンドプロット用の設定ファイル case.insp が設置されます。ファイルを開いて、フェルミエネルギーなどを入力します。<file - case.insp>
-を実行すると、バンドプロット用の設定ファイル case.insp が設置されます。
-ファイルを開いて、フェルミエネルギーなどを入力します。
-<file - case.insp>
 ### Figure configuration
 .0   3.0                         # paper offset of plot
@@ Line 54: / Line 45: @@
       1    0.2                    # jatom, jcol, size  of heavier plotting
 </file>
+    * 9行目の 0.xxxx というところをフェルミエネルギーに置き換えます。
-  * 9行目の 0.xxxx というところをフェルミエネルギーに置き換えます。
+      * フェルミエネルギーは<code>
-  * その他の設定は任意です。
-  * 最後の2行はバンドの重み付けプロットをする際に使います。
-  * jatomやjcolをどうすればよいかは[[WIEN2k:状態密度]]の解説ページも見てください。
-フェルミエネルギーは
-<code>
 $ grep ":FER" case.scf | tail -1
 :FER  : F E R M I - ENERGY(TETRAH.M.)=   0.xxxx
+</code>で確認できます。
+    * その他の設定は任意です。通常の使用ですとあとは energy range を適当な値に変更するくらいです。
+    * 最後の2行はバンドの重み付けプロットをする際に使います。qtlの計算が必要です。
+    * jatomやjcolをどうすればよいかは[[WIEN2k:状態密度]]の解説ページも見てください。
+  - 設定ファイルができたら、もう一度 x spaghetti を実行します。<code>
+$ x spaghetti -p
 </code>
-で確認できます。
+  - case.spaghetti_ps にバンドが出力されます。
-設定ファイルができたら、もう一度 x spaghetti を実行するとバンドが出力されます。
-<code>
-$ x spaghetti -p
-</code>
-case.spaghetti_ps にバンドが出力されます。
 自分でGnuplotなどでバンド図を作りたい場合は、同時に出力される case.spaghetti_ene の4列目（波数）と5列目（エネルギー）を使用してください。
 対称点での波数の値は case.outputsp のxmaxを見ることでわかります。

雑多な記録

User Tools

Site Tools

Differences

Page Tools