User Tools

Site Tools


wien2k:バンド計算

Differences

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

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
wien2k:バンド計算 [2020/03/20 17:58]
koudai [固有エネルギーの計算]
wien2k:バンド計算 [2021/06/27 22:04] (current)
Line 3: Line 3:
 バンド計算を行います。 バンド計算を行います。
  
-バンドの経路を指定するためにXCrysDenが必要です。 +  * バンドの経路を指定するためにXCrysDenが必要です。 
-また、バンドのプロットにGnuplotが必要です+
  
  
Line 24: Line 24:
 ===== 固有エネルギーの計算 ===== ===== 固有エネルギーの計算 =====
  
-指定した経路上の各k点でのエネルギーを求めます。 +  - 指定した経路上の各k点でのエネルギーを求めます。<code>
- +
-<code>+
 $ x lapw1 -band -p $ x lapw1 -band -p
 </code> </code>
 +  -(各バンドにおける軌道ごとの重みを見たい場合のみ)qtlを実行します<code>
 +$ x qtl -p
 +</code>
 +  - 実際にバンドを描いてみます。<code>
 +$ x spaghetti -p
 +</code>を実行すると、バンドプロット用の設定ファイル case.insp が設置されます。ファイルを開いて、フェルミエネルギーなどを入力します。<file - case.insp>
 +### Figure configuration
 + 5.0   3.0                         # paper offset of plot
 +10.0  15.0   3.0                   # xsize,ysize [cm], linebreak-parameter
 + 1.0                             # major ticks, minor ticks
 + 1.0                             # character height, font switch
 + 1.1      4                      # line width, line switch, color switch
 +### Data configuration            
 +-14.0  8.0  2                      # energy range, energy switch (1:Ry, 2:eV)
 +1      0.xxxx                      # Fermi switch,  Fermi-level (in Ry units)
 +1   999                            # number of bands for heavier plotting   1,1
 +0      1    0.2                    # jatom, jcol, size  of heavier plotting   
 +</file>
 +    * 9行目の 0.xxxx というところをフェルミエネルギーに置き換えます。
 +      * フェルミエネルギーは<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 にバンドが出力されます。
  
 +自分でGnuplotなどでバンド図を作りたい場合は、同時に出力される case.spaghetti_ene の4列目(波数)と5列目(エネルギー)を使用してください。
 +対称点での波数の値は case.outputsp のxmaxを見ることでわかります。
  
-===== バンドのプロット ===== 
- 
-バンドをPostScriptファイルに出力します。 
- 
- 
-<code> 
wien2k/バンド計算.1584694709.txt.gz · Last modified: 2021/06/27 22:01 (external edit)