Differences

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

--- wien2k:電場勾配 [2020/04/02 19:47]
koudai [計算結果]
+++ wien2k:電場勾配 [2020/07/10 22:01]
koudai [主軸回転]
@@ 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 184: / Line 185: @@
 ==== 電場勾配 ====
-    :EFG001:                        EFG         =     0.57135   *10**21  V / m**2
+  :EFG005:                        EFG         =     2.96921   *10**21  V / m**2
 は電場勾配 $eq = V_{zz}~\mathrm{[10^{21} V/m^2]}$ の値です。
@@ 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周波数を求める方法は核スピンの大きさで変わってくるので、教科書で確認してください。
 ==== 電気四重極テンソル ====
-.13321    0.39243    0.00000       -0.21827    0.00000    0.00000
+<code>
-.39243    0.21987    0.00000        0.00000    0.57135    0.00000
+             -2.85854   -0.03740    0.27082       -2.88543    0.00000    0.00000
-.00000    0.00000   -0.35308        0.00000    0.00000   -0.35308
+             -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行列は、元の座標系における電気四重極テンソルで、
-の左側の3×3行列は電気四重極テンソルで、
 \begin{equation}
@@ Line 222: / Line 225: @@
 の順に並んでいます。
-右側が主軸変換（対角化）後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 0.57135$, $V_{yy} = -0.35308$, , $V_{xx} = -0.21827$です。
+右側が主軸変換（対角化）後の結果で、固有値の絶対値が大きい順に$|V_{zz}| \geq |V_{yy}| \geq |V_{xx}|$と定義されますので、$V_{zz} = 2.96921$, $V_{yy} = -2.88543$, , $V_{xx} = -0.08379$です。
 ==== 主軸回転 ====
-           MAIN DIRECTIONS OF THE EFG   1.0000  0.8957  0.0000
+           MAIN DIRECTIONS OF THE EFG   1.0000  0.0986 -0.0029
-                                       -0.8957  1.0000  0.0000
+.0103 -0.0762  1.0000
-.0000  0.0000  1.0000
+                                       -0.0978  1.0000  0.0765
 主軸回転後の軸の向きです。
-主軸回転後のz軸の方向（主軸方向）は (0.8957, 1, 0) 、y軸の方向は(0,0,1)、x軸の方向は (1, -0.8957, 0) です。
+列ごとにみてください。
+主軸回転後のz軸の方向（主軸方向）は (-0.0029, 1, 0.0765) 、y軸の方向は(1, 0.0103, -0.0978)、x軸の方向は (0.0986, -0.0762, 1) です。
 ==== 非対称パラメータ ====
-  :ETA001:                         ASYMM. ETA =     0.23594
+  :ETA005:                         ASYMM. ETA =     0.94356
 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。
@@ Line 247: / Line 250: @@
 $ grep ":EFG001:" case.scf</code>で確認できます。収束が十分でないようでしたら、収束条件をより厳しくして再計算してください。
   * k点数に対しても収束が遅いので、必ずk点数に関する収束も確認してください。特に金属の場合はかなり多くのk点をとらないと収束しません。
-  * 周辺の原子位置に敏感な量ですので、実験と合わないときは構造最適化を行ってから計算してみるとよいでしょう
+  * 周辺の原子位置に敏感な量ですので、基本的には内部座標の構造最適化と組み合わせて計算します。run_lapw -p -fc 0.1 -min ...

雑多な記録

User Tools

Site Tools

Differences

Page Tools