Differences

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

--- wien2k:電場勾配 [2020/04/02 20:10] – koudai
+++ wien2k:電場勾配 [2021/06/27 22:04] (current) – external edit 127.0.0.1
@@ Line 7: / Line 7: @@
 プロジェクト名を alpha-CaAlF5とします。
-以下の論文の結晶構造から作成しました。
+以下の論文の結晶構造からstructファイルを作成しました。
   * https://scripts.iucr.org/cgi-bin/paper?du0231
@@ Line 61: / Line 61: @@
   * https://doi.org/10.1021/jp0740696
+    * スライド：http://susi.theochem.tuwien.ac.at/events/ws2017/notes/PB-Hyperfine.pdf
   * この論文によると、構造最適化して計算したほうが良いということですが、今回はチュートリアルですので構造最適化前の結果を示します。
 ===== SCF計算 =====
-  - 収束を厳し目にして、SCF計算を再実行します<code>
+  - SCF計算を実行します<code>
 $ init_lapw -b -numk 100 -rkmax 7.0
-$ run_lapw -ec 0.000001 -cc 0.000001 -i 100 -p
+$ run_lapw -ec 0.00001 -cc 0.0001 -i 100 -p
 </code>
 ===== 計算結果 =====
@@ Line 199: / Line 200: @@
   * 電気四重極モーメント一覧 https://doi.org/10.1016/j.adt.2005.04.001
-  * 例えば27Alですと $Q=0.15~\mathrm{[\times 10^{-28}m^2]}$ 程度です
+  * 例えば27Alですと $Q=0.15~\mathrm{[10^{-28}m^2]}$ 程度です
-  * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です
+  * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です。散乱断面積に使われる単位ですが、NQRでも同じくらいのオーダーの量を扱うために借用されることが多いです。
+  * 四重極結合定数からNQR周波数を求める方法は核スピンの大きさで変わってくるので、教科書で確認してください。
 ==== 電気四重極テンソル ====
-               -2.85854   -0.03740    0.27082       -2.88543    0.00000    0.00000
+<code>
-               -0.03740    2.95115    0.23490        0.00000   -0.08379    0.00000
+             -2.85854   -0.03740    0.27082       -2.88543    0.00000    0.00000
-.27082    0.23490   -0.09260        0.00000    0.00000    2.96921
+             -0.03740    2.95115    0.23490        0.00000   -0.08379    0.00000
+.27082    0.23490   -0.09260        0.00000    0.00000    2.96921
+</code>
 の左側の3×3行列は、元の座標系における電気四重極テンソルで、
@@ Line 230: / Line 234: @@
 主軸回転後の軸の向きです。
+列ごとにみてください。
 主軸回転後のz軸の方向（主軸方向）は (-0.0029, 1, 0.0765) 、y軸の方向は(1, 0.0103, -0.0978)、x軸の方向は (0.0986, -0.0762, 1) です。
@@ Line 245: / Line 250: @@
 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。
   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。
-  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいでしょう
+  * 周辺の原子位置に敏感な量ですので、基本的には内部座標の構造最適化と組み合わせて計算します。run_lapw -p -fc 0.1 -min ...