求助：振动配分函数的非谐校正中的非谐常数项G0的计算

xxkxycb

各位老师好，我目前在做超低温下同位素分馏计算，其中涉及到对振动配分函数进行非谐校正，有一些问题想请教一下。在二阶微扰理论下考虑非谐振动有：

其中gst仅在线性分子中考虑。想请问对于一个非线性多原子分子，公式中的非谐常数G0项可以在Gaussian的输出结果中直接给出吗？对于线性分子的
Vibrational l-type doubling非谐常数项gst应该如何计算求解？
对于线性三原子分子的非谐常数G0项可以通过如下公式求解：

我计算出的HCN分子的立方力常数四次方力常数如下，请问对于存在简并振动模式的项3a 3a 3b 3b应该如何考虑进公式呢？

        I    J    K    FI(I,J,K)   k(I,J,K)    K(I,J,K)
        1    1    1    1827.43416    37.43806     1.27249
        2    1    1     459.86433     7.53638     0.25616
        2    2    1    -105.44524    -1.38236    -0.04699
        2    2    2     623.02524     6.53373     0.22208
        3a   3a   1   -1109.73448    -5.07662    -0.17255
        3b   3b   1   -1109.73448    -5.07662    -0.17255
        3a   3a   2    -208.79593    -0.76408    -0.02597
        3b   3b   2    -208.79593    -0.76408    -0.02597

       8 3rd derivatives larger than 0.371D-04 over      20

........................................................
:                                                      :
:       QUARTIC FORCE CONSTANTS IN NORMAL MODES        :
:                                                      :
: FI =  Reduced values [cm-1]  (default input)         :
: k  =  Quartic Force Const.[AttoJ*amu(-2)*Ang(-4)]    :
: K  =  Quartic Force Const.[Hartree*amu(-2)*Bohr(-4)] :
:......................................................:

   I    J    K    L   FI(I,J,K,L) k(I,J,K,L) K(I,J,K,L)

   1    1    1    1     986.56012   204.20208     3.67285
   2    1    1    1     208.19918    34.47282     0.62004
   2    2    1    1      59.62211     7.89708     0.14204
   2    2    2    1     -37.72022    -3.99664    -0.07188
   2    2    2    2     148.48137    12.58501     0.22636
   3a   3a   1    1    -933.61301   -43.15065    -0.77612
   3b   3b   1    1    -933.61301   -43.15065    -0.77612
   3a   3a   2    1    -176.26199    -6.51689    -0.11722
   3b   3b   2    1    -176.26199    -6.51689    -0.11722
   3a   3a   2    2     -56.78766    -1.67957    -0.03021
   3b   3b   2    2     -56.78766    -1.67957    -0.03021
   3a   3a   3a   3a   1147.55738    11.84346     0.21302
   3a   3a   3b   3b    382.51913     3.94782     0.07101
   3b   3b   3b   3b   1147.55738    11.84346     0.21302