Differences

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

--- wien2k:スピン軌道相互作用 [2020/07/10 19:49] – [手順] koudai
+++ wien2k:スピン軌道相互作用 [2021/07/21 02:11] (current) – koudai
@@ Line 5: / Line 5: @@
 ====== 手順 ======
+===== 常磁性体 =====
   - 通常通り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 16: / 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_c_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 -p -band
+$ x lapwso -p
+$ x qtl -p -so
+$ x spaghetti -p -so
+# case.inspを編集
+$ x spaghetti -p -so
+</code>
+  * 下の図は非磁性体のPtのバンドです
+    * 参考：https://arxiv.org/abs/1309.1356
+  * PtのCIFファイルは https://staff.aist.go.jp/nomura-k/japanese/itscgallary.htm で手に入ります。
+{{ :wien2k:pt-band-so.png?400 |}}