WIEN2k (Calculation of spin-orbit coupling)

Описание: case.in2
-----------
TEMP: EF is calculated where each eigenvalue is temperature broadened using a Fermi function with a broadening parameter of eval Ry. The total energy is corrected corresponding to T=0K. (e.g. eval=0.002 Ry gives good total energy convergence, but has no “physical“ justification)

　最初にSOを入れて計算したらエラーが出ましたので、SO無しで計算し、その結果を使ってSOを計算することにしました。
　k点は後から増やせるので、SCFの時にはかなり減らして計算しました。全エネルギーや電荷の値がフラフラと増減していたので、TEMPにマニュアルの推奨値である0.002 Ryを入力しました（確かに収束までの回数が短くなります。ブラハ教授に、計算がしにくい系はTEMPを0.002 Ryとすると計算の仕方を見せながら教えて頂いております）。
　k点はSCFを再度行なうときに増やす方法もありますが、ここではDOSのときにk点を増やすようにしました。参考にして頂けましたら幸いです。
　RKmax（WIEN2でのCutoff は 13.6 * (RKmax / Rmt ) ^2 [eV]）も後から増やせると思いますけど、どうしても計算が大変な場合は（経験的にも）6ぐらいで良いと思います。マニュアルには「Usually this value should be between 5 and 9 (APW+lo) or 6 - 10. (LAPW basis) (K2 max would be the plane wave cut-off parameter in Ry used in pseudopotential calculations.)と書いてあります。ディフォルトはLAPW basisなので、RKmax=6がマニュアルでの推奨値の下限です。sp成分が多くて系が自由電子的なら5.8や5.9という値で論文を出したこともあります。
　もう少し大きなk点やRKmaxでSCF計算をしても結果が変わらないかチェックして見た目で顕著な違い(つまり、条件を変えてギザギザやスパイク状の電子構造の変化を気にするよりは、計算した物性の解釈が変わらないということ)がなければ学会発表や論文に記載して良いと思います。
　どうしてもある構造（ピーク）が出て欲しいんだというのであれば、精度を高くする必要がありますが、そうでないなら計算条件を明記しておけば良いと思います（実験とあわなければ所詮は理論と思っておくことも大切です）。
　どこまでも精度を突詰める人は環境が整っている人で劣悪な環境の現場を経験していない人でしょう。そういう恵まれた環境にある人が相手の環境を考えて建設的な意見を出せなかった場合は、相手を尊敬しつつ、その程度の方なのだと心の中で思っておきましょう。そういう方に限って”本当の”プロフェッショナルではないです。プロなら相手の環境や低コストでより良い結果を出す方法を提言できます。
　書き忘れていましたが、spinの方向はh k l で 0 0 1の方向（つまりz軸方向に揃えて）にしてスピン軌道相互作用（SO)を計算しています（z軸方向にspinを揃えのがどの第一原理計算コードでも一般的です）。