User Tools

Site Tools


wien2k:フェルミ面

Differences

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

Link to this comparison view

wien2k:フェルミ面 [2021/04/28 22:05]
koudai [フェルミ面の描き方]
wien2k:フェルミ面 [2021/06/27 22:04]
Line 1: Line 1:
-====== 概要 ====== 
- 
-  * フェルミ面を描きます。 
-  * XCrySDenを使います。 
- 
- 
-====== フェルミ面の描き方 ====== 
- 
-SCF計算を事前に済ませておいてください。 
- 
-  - XCrySDenを起動<code> 
-$ xcrysden --wien_fermisurface / 
-</code> 
-    * <nowiki>--wien_fermisurface</nowiki>の後ろには作業ディレクトリのパスが入ります。上は現在いるディレクトリと同じ場合。 
-    * Xcrysdenを起動して [File]→[Open WIEN2k]→[Fermi Surface] で、計算を行っているディレクトリを選択するのでも可 
-  - Number of kpointsを多めにとる(フェルミ面がきれいになるまで徐々に増やしていけば良い) 
-  - [Generate k-mesh]→[Done] 
-  - [Calculate Eigenvalues [x lapw1 ]] 
-    * 並列化計算をする場合はParallel execution of lapw1/lapw2にチェックを入れる 
-    * k点が多いと時間がかかる 
-  - [Calculate Fermi energy [x lapw2 -fermi]] 
-  - [Render Fermi surface] 
-    * Specify the Fermi energyはそのままにして[OK] 
-    * BarGraphやBandwidthから、select bandsのフェルミ準位にかかっているバンドを選択して[Selected] 
-    * それぞれのバンドのフェルミ面を見ることができる。全部を一気に表示したかったら[Merged Bands]を選ぶ 
-    * できたフェルミ面を保存したい場合はフェルミ面の表示領域で右クリックして[File]->[Save Fermi Surface(s) in BXSF format] 
-      * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまうので、下のような感じになるように自分でファイルを開いて結合してください><file case.bxsf> 
- BEGIN_INFO 
-   Fermi Energy:     (フェルミ準位) 
- END_INFO 
- BEGIN_BLOCK_BANDGRID_3D 
- band_energies 
- BANDGRID_3D_BANDS 
- (バンドの数) 
- (k点の数) 
- (基準点) 
- (逆格子ベクトル1) 
- (逆格子ベクトル2) 
- (逆格子ベクトル3) 
- BAND:    1 
-(なんかデータ) 
- BAND:    2 
-(なんかデータ) 
-・・・ 
- BAND:    x 
-(なんかデータ) 
- END_BANDGRID_3D 
- END_BLOCK_BANDGRID_3D 
-</code> 
-    * 再び見るときは次のコマンドで開く<code> 
-$ xcrysden --bxsf case.bxsf.band-x</code> 
- 
-===== 注意 ===== 
- 
-  * XcrysDenのソースコードにk点数の上限が定められている 
-  * もし文句を言われたらXCrysDenのソースコードから wn_readbakgen.f を開いて MAX_NKPT(全ブリルアンゾーンのk点数の上限)と MAX_IRNKPT(既約ブリルアンゾーンのk点数の上限)の値を増やしてコンパイルしなおせばよい。 
- 
  
wien2k/フェルミ面.txt · Last modified: 2021/06/27 22:04 (external edit)