wien2k:wien2wannier:クイックスタート
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
|
wien2k:wien2wannier:クイックスタート [2019/02/09 20:47] koudai [WIEN2WANNIER] |
wien2k:wien2wannier:クイックスタート [2020/03/31 19:31] koudai [ワニエ軌道のプロット] |
||
|---|---|---|---|
| Line 5: | Line 5: | ||
| * 結晶構造は以下のサイトのものを使いました。 | * 結晶構造は以下のサイトのものを使いました。 | ||
| * http:// | * http:// | ||
| - | * Wannier90は並列計算に対応していません。x lapw1 で並列計算をしたい場合は .machines を作成して、オプション -p を付与してください。その場合 WIEN2WANNIERの実行時に x w2w -p としてください（並列計算は行いませんが、並列計算用に用意された出力ファイルを読み込んでくれます）。 | ||
| Line 26: | Line 25: | ||
| 事前にSrVO3というディレクトリを作って、その中で作業します。 | 事前にSrVO3というディレクトリを作って、その中で作業します。 | ||
| 結晶構造のデータから SrVO3.struct を作ってください。 | 結晶構造のデータから SrVO3.struct を作ってください。 | ||
| + | また、並列化する場合は .machines を作成してください。 | ||
| + | 以下では並列化ありの例を示します。 | ||
| 結晶構造ファイルができたら、SCF計算を実行してください。 | 結晶構造ファイルができたら、SCF計算を実行してください。 | ||
| Line 57: | Line 58: | ||
| $ prepare_w2wdir subdir | $ prepare_w2wdir subdir | ||
| + | $ cp .machines subdir/ | ||
| $ cd subdir | $ cd subdir | ||
| Line 68: | Line 70: | ||
| | | ||
| NUMBER OF K-POINTS IN WHOLE CELL: (0 allows to specify 3 divisions of G) | NUMBER OF K-POINTS IN WHOLE CELL: (0 allows to specify 3 divisions of G) | ||
| - | 1000 | + | 300 |
| - | | + | |
| Shift of k-mesh allowed. Do you want to shift: (0=no, 1=shift) | Shift of k-mesh allowed. Do you want to shift: (0=no, 1=shift) | ||
| 0 | 0 | ||
| - | 1000 | + | |
| KGEN ENDS | KGEN ENDS | ||
| 0.032u 0.008s 0:05.33 0.5% 0+0k 3352+1320io 14pf+0w | 0.032u 0.008s 0:05.33 0.5% 0+0k 3352+1320io 14pf+0w | ||
| Line 140: | Line 142: | ||
| </ | </ | ||
| - | * NUMBER OF K-POINTS IN WHOLE CELL ... ワニエ化に使うブロッホ関数のk点数を指定してください。SCF計算で使ったものよりも小さいもので十分です。 | + | * NUMBER OF K-POINTS IN WHOLE CELL ... ワニエ化に使うブロッホ関数のk点数を指定してください。最初は少ないk点数から、うまくワニエ化がいかなかったら増やしていきます。 |
| * Shift of k-mesh allowed ... k点に対称点を含めるかどうかです。通常は対称点を含むようにしますので 0 を選んでください。 | * Shift of k-mesh allowed ... k点に対称点を含めるかどうかです。通常は対称点を含むようにしますので 0 を選んでください。 | ||
| * enter energy window for findbands ... さきほど作成したバンド構造を見ながら、ワニエ化するバンドが全て含まれるエネルギーの範囲を指定してください。フェルミエネルギーを基準 (0 eV) とします。 | * enter energy window for findbands ... さきほど作成したバンド構造を見ながら、ワニエ化するバンドが全て含まれるエネルギーの範囲を指定してください。フェルミエネルギーを基準 (0 eV) とします。 | ||
| * minimal and maximal band indices ... さきほど指定したエネルギーの範囲内に含まれるバンドの指標が表示されるので、その中からワニエ化に使用するバンドの指標を指定します。 | * minimal and maximal band indices ... さきほど指定したエネルギーの範囲内に含まれるバンドの指標が表示されるので、その中からワニエ化に使用するバンドの指標を指定します。 | ||
| * next proj. ... ワニエ化する軌道を選びます。V（番号2に対応）のt2g軌道を使用します。同じ原子の軌道であれば 2: | * next proj. ... ワニエ化する軌道を選びます。V（番号2に対応）のt2g軌道を使用します。同じ原子の軌道であれば 2: | ||
| - | * write_win ... 自動的にWannier90の入力ファイルが作成されます。projectionsのところにさきほど指定した軌道を入力します< | + | * write_win ... 自動的にWannier90の入力ファイルが作成されます。通常はそのままで大丈夫ですが、必要に応じて修正します |
| - | begin projections | + | |
| - | | + | |
| - | | + | |
| - | | + | |
| - | end projections | + | |
| - | </ | + | |
| ワニエ化に使うk点での固有状態を求めます。 | ワニエ化に使うk点での固有状態を求めます。 | ||
| - | $ x lapw1 | + | $ x_lapw |
| - | WIEN2WANNIERにより、WIEN2kのデータをWannier90が使えるように変換します。 | + | WIEN2WANNIERにより、WIEN2kのデータをWannier90が使えるように変換します |
| + | （並列計算は行いませんが、フラグ-pをつけることで並列計算用に用意された出力ファイルを読み込んでくれます。以下の x_lapw wplot も同様です）。 | ||
| - | $ x w2w | + | $ x_lapw |
| - | Wannier90を使って実際にワニエ化を行います。 | + | Wannier90を使って実際にワニエ化を行います。Wannier90自体は並列化に対応していません。 |
| - | $ x wannier90 | + | $ x_lapw |
| ワニエ化したことで得られたハミルトニアンによるバンドは subdir_band.dat に出力されます。 | ワニエ化したことで得られたハミルトニアンによるバンドは subdir_band.dat に出力されます。 | ||
| Line 176: | Line 173: | ||
| ワニエ軌道がどのようになっているのかプロットしてみます。 | ワニエ軌道がどのようになっているのかプロットしてみます。 | ||
| - | $ x wplot -p # 入力ファイルの subdir.inwplot を生成します。 | + | $ x_lapw |
| $ vim subdir.inwplot | $ vim subdir.inwplot | ||
| Line 184: | Line 181: | ||
| 実空間のワニエ関数を計算します。 | 実空間のワニエ関数を計算します。 | ||
| - | $ x wplot -p | + | $ x_lapw |
| wplotで生成したファイルをXCrysDenで読み込める形式に変換します。 | wplotで生成したファイルをXCrysDenで読み込める形式に変換します。 | ||
| Line 192: | Line 189: | ||
| 得られたワニエ軌道は次のようにして見ることができます（1番目のワニエ軌道を描くときの例です）。 | 得られたワニエ軌道は次のようにして見ることができます（1番目のワニエ軌道を描くときの例です）。 | ||
| - | $ xcrysden | + | $ xcrysden --xsf subdir_1.xsf |
| - [Modify] → [Atomic Radius] → SpaceFill/ | - [Modify] → [Atomic Radius] → SpaceFill/ | ||
wien2k/wien2wannier/クイックスタート.txt · Last modified: 2021/06/27 22:04 (external edit)
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
