This is an old revision of the document!
$ init_lapw -b -numk 2000 -rkmax 8.5 $ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100
$ init_so_lapw
$ rm *.broyd* $ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100 -so
case.klist_bandを用意してください
$ x lapw1 -band -so -p $ x lapwso $ x lapw2 -qtl -band -so -p $ x spaghetti -so -p # case.inspを編集 $ x spaghetti -so -p
うまく計算できているかどうかは、たとえばPt単体などで計算して https://arxiv.org/abs/1309.1356 の図2が再現できるかやってみると良いでしょう。 CIFファイルは https://staff.aist.go.jp/nomura-k/japanese/itscgallary.htm で手に入ります。