• WIEN2kではSCF計算の際に、NQR周波数を求めるのに使われる電場勾配 (EFG) の計算も同時に行ってくれます。

case.scf0 を開いてください。 以下はCaGa4の例です。

CaGe4.scf0

:EFG001:                        EFG         =     0.65471   *10**21  V / m**2
                               V20  TOT/SRF=    -0.30929    -0.00945
                               V22  TOT/SRF=    -0.16361    -0.00535
                               V22M TOT/SRF=     0.44715     0.01402
                               V21  TOT/SRF=     0.00000     0.00000
                               V21M TOT/SRF=     0.00000     0.00000

              0.01496    0.44715    0.00000       -0.29757    0.00000    0.00000
              0.44715    0.34218    0.00000        0.00000    0.65471    0.00000
              0.00000    0.00000   -0.35714        0.00000    0.00000   -0.35714

         MAIN DIRECTIONS OF THE EFG   1.0000  0.6989  0.0000
                                     -0.6989  1.0000  0.0000
                                      0.0000  0.0000  1.0000
:ANG001:  ANGLE WITH OLD X-AXIS =         0.0

:ETA001:                         ASYMM. ETA =     0.00000

がEFG計算の結果です。

—

:EFG001:                        EFG         =     0.65471   *10**21  V / m**2

電場勾配 $eq = V_{zz} = 0.65471 \times 10^{21}$ [V/m^2]の値です。 $Q \mathrm{m^2}$を考えている原子核のQ値とすれば、四重極結合定数は

\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]} } V_{zz} Q = 2.4179893 \times 10^{14} \times V_{zz} Q \mathrm{[Hz]} = 2.4179893 \times 10^{8} \times V_{zz} Q \mathrm{[MHz]} \end{equation}

となります。 ほとんどの原子核でおよそ数MHz程度の値になります。