• write_inwfの最初の行を #!/usr/bin/env python2.7 としておく．
• ulimit -s unlimited
  • w2wの計算の最後でスタックサイズが巨大になるので，unlimitedにしておかないとsegmentation faultすることがある．~/.bashrcに書いておくのが良いだろう．
• 並列化するなら.machinesファイルを用意しておく．以下，並列化を仮定(オプション-p)．
• $WIENROOT/SRC_templates/case.win内の変数hr_plotをwrite_hrと書き換えておく（最新版のwannier90の仕様に合わせる）．

実行

1. prepare_w2wdir WORK で作業ディレクトリWORKをつくる．

2. cd WORK

3. init_w2w -so

   1. 基底軌道は2N個指定することになる．スピン↑と↓で同じ軌道を採るのが自然であろう．1:s 2:p ... 1:s 2:p ... など．

4. case.inwfup, case.inwfdn を用意する．inwfupにおけるスピン↓の射影は0とするのが良い．単に0と書く．

5. w2wso.shを実行する．２行目のcaseは適切に変更する．

   #!/bin/sh
   case="wan"
   cp ../.machines .
   cp $case.fermi $case.fermiup
   cp $case.fermi $case.fermidn
   cp $case.vns $case.vnsup
   cp $case.vns $case.vnsdn
   cp $case.vsp $case.vspup
   cp $case.vsp $case.vspdn
   x lapw1 -up -p
   x lapw1 -dn -p
   x lapwso -up -p
   #cp $case.inwf $case.inwfup
   #cp $case.inwf $case.inwfdn
   x w2w -so -up -p &
   x w2w -so -dn -p
   x wannier90 -so -pp
   

6. mpirun wannier90.x WORK