Differences

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

--- wien2k:電場勾配 [2020/03/28 14:38]
koudai [電場勾配]
+++ wien2k:電場勾配 [2020/03/31 20:33]
koudai [例：CaGa4]
@@ Line 3: / Line 3: @@
   * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。
-====== 注意 ======
-  * 収束は非常に遅いので<code>
+====== 例：CaGa4 ======
-$ grep ":EFG001:" case.scf</code>で必ず収束を確認してください。
-  * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。
+以下の論文の結晶構造から作成しました。
-  * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください
-  * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。
+  * https://www.sciencedirect.com/science/article/pii/0022508889901719
+<file - CaGa4.struct>
+blebleble
+CXZ LATTICE,NONEQUIV.ATOMS   3   12 C2/m
+MODE OF CALC=RELA unit=bohr
+.680586 11.559833 11.584210 90.000000 90.000000118.942000
+ATOM  -1: X=0.00000000 Y=0.00000000 Z=0.00000000
+          MULT= 1          ISPLIT=15
+Ca         NPT=  781  R0=.000050000 RMT= 2.50000     Z:  20.00000
+LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
+.0000000 1.0000000 0.0000000
+.0000000 0.0000000 1.0000000
+ATOM  -2: X=0.40503000 Y=0.76923000 Z=0.00000000
+          MULT= 2          ISPLIT=15
+      -2: X=0.59497000 Y=0.23077000 Z=0.00000000
+Ga         NPT=  781  R0=.000050000 RMT= 2.33        Z:  31.00000
+LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
+.0000000 1.0000000 0.0000000
+.0000000 0.0000000 1.0000000
+ATOM  -3: X=0.00000000 Y=0.50000000 Z=0.23283000
+          MULT= 2          ISPLIT=15
+      -3: X=0.00000000 Y=0.50000000 Z=0.76717000
+Ga         NPT=  781  R0=.000050000 RMT= 2.33        Z:  31.00000
+LOCAL ROT MATRIX:    1.0000000 0.0000000 0.0000000
+.0000000 1.0000000 0.0000000
+.0000000 0.0000000 1.0000000
+      NUMBER OF SYMMETRY OPERATIONS
+</file>
+この物質のWIEN2kによるNQRの研究は、以下の論文の5.1節が参考になります
+  * https://d-nb.info/1007751142/34
+電荷の収束を厳し目にして、通常のSCF計算を実行します
+<code>
+$ init_lapw -b -numk 1000 -rkmax 7.5
+$ run_lapw -cc 0.00001 -ec 0.00001 -i 100 -p
+</code>
 ====== 電場勾配 ======
-case.scf0 を開いてください。
+SCF計算が終了したら case.scf0 を開いてください。
 以下はCaGa4の例です。
@@ Line 35: / Line 74: @@
 :ETA001:                         ASYMM. ETA =     0.09099
 </file>
-がEFG計算の結果です。
+が1番目の原子のEFG計算の結果です。
 ===== 電場勾配 =====
@@ Line 42: / Line 81: @@
 は電場勾配 $eq = V_{zz}~\mathrm{[10^{21} V/m^2]}$ の値です。
-電気四重極モーメントを$Q~\mathrm{[10^{-30} m^2]}$とすれば、四重極結合定数は
+電気四重極モーメントを$Q~\mathrm{[10^{-28} m^2]}$とすれば、四重極結合定数は
 \begin{equation}
 C_q = \frac{e^2 q Q}{h}
-= \frac{1.60217662 \times 10^{-19}~\mathrm{[C]}}{6.62607004 \times 10^{-34}~\mathrm{[J \cdot s]} } \times Q~\mathrm{[10^{-30} m^2]} \times  V_{zz}~\mathrm{[10^{21} V/m^2]}
+= \frac{1.60217662 \times 10^{-19}~\mathrm{[C]}}{6.62607004 \times 10^{-34}~\mathrm{[J \cdot s]} } \times Q~\mathrm{[10^{-28} m^2]} \times  V_{zz}~\mathrm{[10^{21} V/m^2]}
-= 0.24179893 \times Q V_{zz}~\mathrm{[MHz]}
+= 24.179893 \times Q V_{zz}~\mathrm{[MHz]}
 \end{equation}
@@ Line 54: / Line 93: @@
   * 電気四重極モーメント一覧 https://www.sciencedirect.com/science/article/pii/S0092640X05000239
-  * 例えば43Caですと $Q=4.3~\mathrm{\times 10^{-30}m^2}$ です (bはバーンという単位で、$1~\mathrm{[b]} = $10^{-28}~\mathrm{[m^2]}$と定義されます)。
+  * 例えば43Caですと $Q=4.3~\mathrm{[\times 10^{-28}m^2]}$ です
+  * bはバーンという単位で、$1~\mathrm{[b]} = 10^{-28}~\mathrm{[m^2]}$です
@@ Line 83: / Line 123: @@
 .0000  0.0000  1.0000
-主軸回転後の向きベクトルです。
+主軸回転後の軸の向きです。
-主軸回転後のZ軸の方向（主軸方向）は (0.6989, 1, 0) 、y軸の方向は(0,0,1)、x軸の方向は (1, -0.6989, 0) です。
+主軸回転後のz軸の方向（主軸方向）は (0.6989, 1, 0) 、y軸の方向は(0,0,1)、x軸の方向は (1, -0.6989, 0) です。
@@ Line 92: / Line 132: @@
 非対称パラメータ $\eta = \frac{V_{xx} - V_{yy}}{V_{zz}}$ の値です。
+====== 注意 ======
+  * 収束は非常に遅いので<code>
+$ grep ":EFG001:" case.scf</code>で必ず収束を確認してください。
+  * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。
+  * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください
+  * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。

雑多な記録

User Tools

Site Tools

Differences

Page Tools