微扰方式考虑自旋－轨道耦合

注意：

• 必须使用Unrestricted
• 必须设置对称性为NoSym(如果使用图形界面设置参数，在这种情况下，软件会自动设置为NoSym，因此推荐尽量使用图形界面设置参数)
• 默认计算所有的g、A因子
• 必须使用全电子基组（Main > Frozen Core: None）
• 对于重原子核，使用Gaussian型有限核（Details > Relativity > Nuclear model: Gaussian），也即是原子核不再作为点电荷处理，而是有一定的大小
• 必须使用特别大的基组，推荐QZ4P，尤其推荐J > QZ4P-J
• 杂化泛函中，只有PBE0可用
• 不能使用metaGGA和metahybrids泛函

选中要计算的原子Ti，点击右边窗口的加号：

SCM > Output，在下方搜索框搜索：TOTAL EPR Delta g-matrix，可以找到。

===  TOTAL EPR Delta g-matrix (parts per thousand)
 
                        ***********************************
                        CARTESIAN AXIS REPRESENTATION
 
                        ==== total Delta g-matrix
                        -----------------------------------
                             118.025    -0.819    13.574
                              -0.819   125.910    -2.383
                              13.574    -2.383   103.853
                        -----------------------------------
                                  isotropic =    115.929
 
                        ***********************************
                        -----------------------------------
                        PRINCIPAL AXIS REPRESENTATION
 
                        ==== Principal components:
 
                              95.540   124.187   128.061
 
                        ==== Principal Axis System:
 
                              -0.515     0.606     0.607
                               0.053     0.729    -0.682
                               0.856     0.319     0.407
                        -----------------------------------

该值是相对于参考值 2.002319（通用与所有体系）的变化量。因此 g 因子=2.002319+115.929‰=2.118248

搜索：Total calculated hyperfine coupling，可以得到Cartesian Axis Representation与Principal Axis Representation下，A的值（au、MHz两种单位均列出），例如前者：

 Total calculated hyperfine coupling:    
  isotropic    A(au)       2.219                           A(MHz)     -66.837
       1.458       0.000       0.000            -43.924       0.000      -0.001
       0.000       2.059       0.000              0.000     -62.032       0.000
       0.000       0.000       3.139              0.001       0.000     -94.554
 
        span        skew       asymm               span        skew       asymm
       1.681       0.285       0.653            -50.630      -8.577     -19.682

其中各向同性A=-66.837MHz。