Differences

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

--- wien2k:スピン軌道相互作用 [2020/07/10 18:13]
koudai [手順]
+++ wien2k:スピン軌道相互作用 [2021/06/27 22:04]
127.0.0.1 external edit
@@ Line 5: / Line 5: @@
 ====== 手順 ======
+===== 常磁性体 =====
   - 通常通りstructファイルを用意します
-  - 通常のSCF計算を実行します。Rkmaxやk点数は物質によって調整してください<code>
+  - SCF計算を実行します。Rkmaxやk点数は物質によって調整してください<code>
 $ init_lapw -b -numk 2000 -rkmax 8.5
-$ run -p -cc 0.0001 -ec 0.00001 -i 100
+$ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100
 </code>
   - スピン軌道相互作用を入れた計算の初期設定を行います<code>
@@ Line 15: / Line 17: @@
 </code>
     * 非磁性体の場合は何を聞かれてもEnterを押すだけで良いです
-    * 入力ファイル case.inso ができます
+    * 入力ファイル case.inso と case.in1 ができます
   - スピン軌道相互作用を含んだSCF計算を実行します<code>
-$ run -p -cc 0.0001 -ec 0.00001 -i 100 -so
+$ rm *.broyd*
+$ run_lapw -p -cc 0.0001 -ec 0.00001 -i 100 -so
+</code>
+  * 物理量の計算では 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>
-状態密度などの計算で再度 x_lapw を実行するときは、フラグ -so をつけるようにしてください。
+  * うまく計算できているかどうかは、たとえば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