请问 非简谐频率计算问题

最爱喵星人

附件是我计算CO2非简谐频率的输出文件。有个问题请教一下大家。

以下是立方力常数和四次方力常数，3a和3b对应着简并的弯曲振动。


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

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

     1      1      1            264.87753      1.36572      0.04642
     2      2      1            506.82404      4.58759      0.15593
     3a     3a     1           -158.41000     -0.40183     -0.01366
     3b     3b     1           -158.41000     -0.40183     -0.01366

Num. of 3rd derivatives larger than  0.371D-04: 4 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      40.90120      1.34523      0.02420
     2      2      1      1      80.40071      4.64230      0.08350
     2      2      2      2     161.10262     16.33005      0.29372
     3a     3a     1      1     -47.63350     -0.77076     -0.01386
     3b     3b     1      1     -47.63350     -0.77076     -0.01386
     3a     3a     2      2    -113.42440     -3.22200     -0.05795
     3b     3b     2      2    -113.42440     -3.22200     -0.05795
     3a     3a     3a     3a     74.08378      0.58976      0.01061
     3a     3a     3b     3b     24.69459      0.19659      0.00354
     3b     3b     3b     3b     74.08378      0.58976      0.01061

在最后输出fundamental bands时
Fundamental Bands
-----------------
    Mode(n,l)    Status       E(harm)  E(anharm)     Ba(x)      Ca(y)
H      1(1,+0)   active      1371.890  1406.119   0.389995   0.389995
          2(1,+0)   active      2408.411  2362.154   0.388094   0.388094
          3(1,-1)   active        674.939   671.428   0.391912   0.391912
          3(1,+1)   active       674.939   671.428   0.391912   0.391912

弯曲振动的简谐频率和非简谐频率差很小。

请问非简谐频率与简谐频率，非简谐常数之间的关系式是什么？

谢谢大家。

wzkchem5

没有简单关系式。非简谐频率要比简谐频率多做很多计算才能算出来，不是一个公式就能写清楚的。
如果一定要知道关系式，看doi: 10.1063/1.1824881

最爱喵星人

wzkchem5 发表于 2022-4-2 16:34
没有简单关系式。非简谐频率要比简谐频率多做很多计算才能算出来，不是一个公式就能写清楚的。
如果一定要 ...

非常感谢，文章非常有用！