This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
wien2k:電場勾配 [2020/03/27 11:14] koudai [電場勾配] |
wien2k:電場勾配 [2020/03/27 11:26] koudai [概要] |
||
---|---|---|---|
Line 2: | Line 2: | ||
* WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | ||
+ | * 収束は非常に遅いので< | ||
+ | $ grep ": | ||
+ | で必ず収束を確認してください。 | ||
+ | * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。 | ||
====== 電場勾配 ====== | ====== 電場勾配 ====== | ||
Line 39: | Line 42: | ||
\begin{equation} | \begin{equation} | ||
C_q = \frac{e^2 q Q}{h} | 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]} } V_{zz} | + | = \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 |
- | = 2.4179893 \times V_{zz} | + | = 2.4179893 \times |
- | = 0.24179893 \times V_{zz} | + | = 0.24179893 \times |
\end{equation} | \end{equation} | ||
Line 54: | Line 57: | ||
0.00000 | 0.00000 | ||
- | は電気四重極テンソルで、 | + | の左側の3×3行列は電気四重極テンソルで、 |
\begin{equation} | \begin{equation} | ||
Line 65: | Line 68: | ||
の順に並んでいます。 | の順に並んでいます。 | ||
- | 右側が何もしていない結果、左側が主軸変換(対角化)後の結果です。 | + | 右側が主軸変換(対角化)後の結果です。 |
---- | ---- |