This shows you the differences between two versions of the page.
wien2k:フェルミ面 [2021/04/28 21:59] koudai [フェルミ面の描き方] |
wien2k:フェルミ面 [2021/06/27 22:04] |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== 概要 ====== | ||
- | |||
- | * フェルミ面を描きます。 | ||
- | * XCrySDenを使います。 | ||
- | |||
- | |||
- | ====== フェルミ面の描き方 ====== | ||
- | |||
- | SCF計算を事前に済ませておいてください。 | ||
- | |||
- | - XCrySDenを起動< | ||
- | $ xcrysden --wien_fermisurface / | ||
- | </ | ||
- | * < | ||
- | * Xcrysdenを起動して [File]→[Open WIEN2k]→[Fermi Surface] で、計算を行っているディレクトリを選択するのでも可 | ||
- | - Number of kpointsを多めにとる(フェルミ面がきれいになるまで徐々に増やしていけば良い) | ||
- | - [Generate k-mesh]→[Done] | ||
- | - [Calculate Eigenvalues [x lapw1 ]] | ||
- | * 並列化計算をする場合はParallel execution of lapw1/ | ||
- | * k点が多いと時間がかかる | ||
- | - [Calculate Fermi energy [x lapw2 -fermi]] | ||
- | - [Render Fermi surface] | ||
- | * Specify the Fermi energyはそのままにして[OK] | ||
- | * BarGraphやBandwidthから、select bandsのフェルミ準位にかかっているバンドを選択して[Selected] | ||
- | * それぞれのバンドのフェルミ面を見ることができる。全部を一気に表示したかったら[Merged Bands]を選ぶ | ||
- | * できたフェルミ面を保存したい場合はフェルミ面の表示領域で右クリックして[File]-> | ||
- | * 保存したファイルはそれぞれのバンドのフェルミ面で分かれてしまうので、自分でファイルを開いて結合してください | ||
- | - BANDGRID_3D_BANDS を結合するバンドの数にする | ||
- | - BAND: x から END_BANDGRID_3D の一個上の行までを順番にはりつけていく | ||
- | * 再び見るときは次のコマンドで開く< | ||
- | $ xcrysden --bxsf case.bxsf.band-x</ | ||
- | |||
- | ===== 注意 ===== | ||
- | |||
- | * XcrysDenのソースコードにk点数の上限が定められている | ||
- | * もし文句を言われたらXCrysDenのソースコードから wn_readbakgen.f を開いて MAX_NKPT(全ブリルアンゾーンのk点数の上限)と MAX_IRNKPT(既約ブリルアンゾーンのk点数の上限)の値を増やしてコンパイルしなおせばよい。 | ||
- | |||