Differences

This shows you the differences between two versions of the page.

--- wien2k:スピン軌道相互作用 [2020/07/10 20:06]
koudai [常磁性体]
+++ wien2k:スピン軌道相互作用 [2021/06/27 22:04]
127.0.0.1 external edit
@@ Line 9: / Line 9: @@
   - 通常通りstructファイルを用意します
-  - スピン分極を許してSCF計算を実行します。Rkmaxやk点数は物質によって調整してください<code>
+  - SCF計算を実行します。Rkmaxやk点数は物質によって調整してください<code>
-$ instgen -ask   # 非磁性なので全部nにする
+$ init_lapw -b -numk 2000 -rkmax 8.5
-$ init_lapw -b -numk 2000 -rkmax 8.5 -sp
+$ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100
-$ runsp_c_lapw -p -cc 0.0001 -ec 0.00001 -i 100
 </code>
   - スピン軌道相互作用を入れた計算の初期設定を行います<code>
@@ Line 18: / Line 17: @@
 </code>
     * 非磁性体の場合は何を聞かれてもEnterを押すだけで良いです
-    * Do you have a spinpolarized case (and want to run symmetso) ? (y/N) y
+    * 入力ファイル case.inso と case.in1 ができます
-    * Do you want to use the new structure for SO calculations ? (y/N) y
-    * 入力ファイル case.inso ができます
   - スピン軌道相互作用を含んだSCF計算を実行します<code>
 $ rm *.broyd*
-$ runsp_lapw -p -cc 0.0001 -ec 0.00001 -i 100 -so
+$ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100 -so
 </code>
-状態密度などの計算で再度 x_lapw を実行するときは、フラグ -so をつけるようにしてください。
+  * 物理量の計算では x_lapw lapw1 のあとに必ず x_lapw lapwso を実行してください
+  * その他の場合については、x_lapw を実行する際にフラグ -so をつけてください
 ==== 例：バンド構造 ====
+case.klist_bandを用意してください
 <code>
-$ x lapw1 -band -up -p
+$ x lapw1 -p -band
-$ x lapw1 -band -dn -p
+$ x lapwso -p
-$ x lapwso -up -p
+$ x qtl -p -so
-$ x spaghetti -up -p
+$ x spaghetti -p -so
 # case.inspを編集
-$ x spaghetti -up -p
+$ x spaghetti -p -so
 </code>
+  * うまく計算できているかどうかは、たとえばPt単体などで計算して https://arxiv.org/abs/1309.1356 の図2が再現できるかやってみると良いでしょう。
+  * PtのCIFファイルは https://staff.aist.go.jp/nomura-k/japanese/itscgallary.htm で手に入ります。

雑多な記録

User Tools

Site Tools

Differences

Page Tools