计算化学公社

标题: 请问非简谐频率计算问题 [打印本页]

作者
Author: 最爱喵星人 时间: 2022-4-2 15:39
标题: 请问非简谐频率计算问题

附件是我计算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

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

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

谢谢大家。

作者
Author: wzkchem5 时间: 2022-4-2 16:34
没有简单关系式。非简谐频率要比简谐频率多做很多计算才能算出来，不是一个公式就能写清楚的。
如果一定要知道关系式，看doi: 10.1063/1.1824881

作者
Author: 最爱喵星人 时间: 2022-4-2 16:53

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)