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/26 21:30] koudai [電場勾配] |
wien2k:電場勾配 [2020/03/28 13:56] koudai |
||
---|---|---|---|
Line 3: | Line 3: | ||
* WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | * WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。 | ||
+ | ====== 注意 ====== | ||
+ | |||
+ | * 収束は非常に遅いので< | ||
+ | $ grep ": | ||
+ | * 収束してないようでしたら、run_lapwでオプション-ccを使って、電荷の収束をより厳しい条件で再計算してください。 | ||
+ | * また、k点数に対しても収束が遅いので、必ずk点数を変化させて収束を確認してください | ||
+ | * 周辺の原子位置に敏感な量ですので、構造最適化を行った後に計算することをおすすめします。 | ||
====== 電場勾配 ====== | ====== 電場勾配 ====== | ||
Line 29: | Line 36: | ||
</ | </ | ||
がEFG計算の結果です。 | がEFG計算の結果です。 | ||
+ | |||
+ | ===== 電場勾配 ===== | ||
: | : | ||
- | 電場勾配 $eq = V_{zz} | + | は電場勾配 $eq = V_{zz} \mathrm{[10^{21} V/m^2]}$ の値です。 |
- | $Q \mathrm{m^2}$を考えている原子核のQ値とすれば、四重極結合定数は | + | 電気四重極モーメントを$Q \mathrm{[10^{-30} |
\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} Q | + | = \frac{1.60217662 \times 10^{-19} \mathrm{[C]}}{6.62607004 \times 10^{-34} \mathrm{[J \cdot s]} } \times |
- | = 2.4179893 | + | = 2.4179893 \times |
- | = 2.4179893 \times 10^{8} \times V_{zz} | + | = 0.24179893 |
\end{equation} | \end{equation} | ||
となります。 | となります。 | ||
- | ほとんどの原子核でおよそ数MHz程度の値になります。 | + | ほとんどの原子核でおよそ数〜数十MHz程度の値になります。 |
+ | |||
+ | |||
+ | ===== 電気四重極テンソル ===== | ||
+ | |||
+ | 0.01496 | ||
+ | 0.44715 | ||
+ | 0.00000 | ||
+ | |||
+ | の左側の3×3行列は電気四重極テンソルで、 | ||
+ | |||
+ | \begin{equation} | ||
+ | \begin{pmatrix} | ||
+ | V_{xx} & V_{xy} & V_{xz} \\ | ||
+ | V_{yx} & V_{yy} & V_{yz} \\ | ||
+ | V_{zx} & V_{zy} & V_{zz} | ||
+ | \end{pmatrix} | ||
+ | \end{equation} | ||
+ | |||
+ | の順に並んでいます。 | ||
+ | 右側が主軸変換(対角化)後の結果です。 | ||
+ | |||
+ | ===== 主軸回転 ===== | ||
+ | |||
+ | MAIN DIRECTIONS OF THE EFG | ||
+ | | ||
+ | 0.0000 | ||
+ | |||
+ | は主軸回転後の向きベクトルです。 | ||
+ | 固有値の絶対値が大きい順に$V_{zz}, |