User Tools

Site Tools


wien2k:wien2wannier:クイックスタート

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:wien2wannier:クイックスタート [2019/02/09 20:49]
koudai [概要]
wien2k:wien2wannier:クイックスタート [2021/06/27 22:04] (current)
Line 5: Line 5:
   * 結晶構造は以下のサイトのものを使いました。   * 結晶構造は以下のサイトのものを使いました。
     * http://www.crystallography.net/cod/1512117.html     * http://www.crystallography.net/cod/1512117.html
-  * Wannier90は並列計算に対応していません。lapw1 で並列計算をしたい場合は .machines を作成して、オプション -p を付けて実行してください。その場合は 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/.machines
   $ cd subdir   $ cd subdir
  
Line 68: Line 70:
             symmetry operations without inversion             symmetry operations without inversion
   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 
- length of reciprocal lattice vectors:   0.866   0.866   0.866  10.000  10.000  10.000+ length of reciprocal lattice vectors:   0.866   0.866   0.866   6.694   6.694   6.694
   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  k-points generated, ndiv=          10          10          10+         216  k-points generated, ndiv=           6                     6
 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:
 </code> </code>
  
-  *  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:dxz,dyz,dxy と一気に指定することも可能です。   * next proj. ... ワニエ化する軌道を選びます。V(番号2に対応)のt2g軌道を使用します。同じ原子の軌道であれば 2:dxz,dyz,dxy と一気に指定することも可能です。
-  * write_win ... 自動的にWannier90の入力ファイルが作成されます。projectionsところにさきほど指定した軌道を入力します<file> +  * write_win ... 自動的にWannier90の入力ファイルが作成されます。通常はそのままで大丈夫でが、必要に応じて修正します
-begin projections +
-   2:dxz +
-   2:dyz +
-   2:dxy +
-end projections +
-</file>その他は必要に応じて修正します+
  
 ワニエ化に使うk点での固有状態を求めます。 ワニエ化に使うk点での固有状態を求めます。
  
-  $ lapw1+  $ x_lapw lapw1 -p
  
-WIEN2WANNIERにより、WIEN2kのデータをWannier90が使えるように変換します。+WIEN2WANNIERにより、WIEN2kのデータをWannier90が使えるように変換します 
 +(並列計算は行いませんが、フラグ-pをつけることで並列計算用に用意された出力ファイルを読み込んでくれます。以下の x_lapw wplot も同様です)
  
-  $ w2w+  $ x_lapw w2w -p
  
-Wannier90を使って実際にワニエ化を行います。+Wannier90を使って実際にワニエ化を行います。Wannier90自体は並列化に対応していません
  
-  $ wannier90+  $ x_lapw wannier90
  
 ワニエ化したことで得られたハミルトニアンによるバンドは subdir_band.dat に出力されます。 ワニエ化したことで得られたハミルトニアンによるバンドは subdir_band.dat に出力されます。
Line 176: Line 173:
 ワニエ軌道がどのようになっているのかプロットしてみます。 ワニエ軌道がどのようになっているのかプロットしてみます。
  
-  $ wplot -p           # 入力ファイルの subdir.inwplot を生成します。+  $ x_lapw wplot -p      # 入力ファイルの subdir.inwplot を生成します。
   $ vim subdir.inwplot   # 必要に応じて修正します。   $ vim subdir.inwplot   # 必要に応じて修正します。
  
Line 184: Line 181:
 実空間のワニエ関数を計算します。 実空間のワニエ関数を計算します。
  
-  $ wplot -p+  $ x_lapw wplot -p
  
 wplotで生成したファイルをXCrysDenで読み込める形式に変換します。 wplotで生成したファイルをXCrysDenで読み込める形式に変換します。
Line 192: Line 189:
 得られたワニエ軌道は次のようにして見ることができます(1番目のワニエ軌道を描くときの例です)。 得られたワニエ軌道は次のようにして見ることができます(1番目のワニエ軌道を描くときの例です)。
  
-  $ xcrysden xcrysden --xsf subdir_1.xsf+  $ xcrysden --xsf subdir_1.xsf
  
   - [Modify] → [Atomic Radius] → SpaceFill/Ball factor を0.5に設定 → [OK]   - [Modify] → [Atomic Radius] → SpaceFill/Ball factor を0.5に設定 → [OK]
wien2k/wien2wannier/クイックスタート.1549712977.txt.gz · Last modified: 2021/06/27 22:01 (external edit)