当前位置: 首页 > 量子化学 >计算线性三原子分子的振动频率,弯曲振动频率是个虚频?

计算线性三原子分子的振动频率,弯曲振动频率是个虚频?

作者 孙1986
来源: 小木虫 1550 31 举报帖子
+关注

使用molpro计算线性三原子分子的频率时,得到弯曲振动频率为虚频,是怎么一回事?

计算线性三原子分子的振动频率,弯曲振动频率是个虚频?
TIM图片20170804113550.png 返回小木虫查看更多

今日热帖
  • 精华评论
  • 121099

    把角度的自由度打开重新优化一下结构看看

  • 孙1986

    引用回帖:
    3楼: Originally posted by 121099 at 2017-08-05 02:20:06
    把角度的自由度打开重新优化一下结构看看

    好的,我试一下

  • 孙1986

    引用回帖:
    3楼: Originally posted by 121099 at 2017-08-05 02:20:06
    把角度的自由度打开重新优化一下结构看看

    角度自由度是打开的

  • 121099

    引用回帖:
    5楼: Originally posted by 孙1986 at 2017-08-06 15:12:44
    角度自由度是打开的...

    三原子分子结构自由度数目并不多,如果几何优化给出的是全局最小的话(所有自由度都打开),那么就可以检查一下每个频率单点计算看看是否有收敛错的单点结果。

    我个人碰到这样的情况的时候,一般选择用笛卡尔坐标系c1对称性优化,来确保优化是全局最小,然后再用优化所得笛卡尔坐标进行频率计算。

    如果用的方法是DFT的话,个人经历过很多情况设置grid对几何优化及频率计算影响很大,可以选择更小的grid数值。
    如果是WF办法的话,除了以上提到的问题暂时想不出其他的可能原因,你可以把input贴一下。

  • 孙1986

    引用回帖:
    6楼: Originally posted by 121099 at 2017-08-08 03:16:17
    三原子分子结构自由度数目并不多,如果几何优化给出的是全局最小的话(所有自由度都打开),那么就可以检查一下每个频率单点计算看看是否有收敛错的单点结果。

    我个人碰到这样的情况的时候,一般选择用笛卡尔坐 ...

    这是输入文件
    计算线性三原子分子的振动频率,弯曲振动频率是个虚频?-1
    TIM图片20170808113729.png

  • 121099

    引用回帖:
    7楼: Originally posted by 孙1986 at 2017-08-08 11:34:26
    这是输入文件

    TIM图片20170808113729.png
    ...

    用DZ优化了一下结构。
    发现如果打开夹角至175以下,则可以找到能量更低的PES极小点,
    DZ优化结果为
    R1=                  1.20567667 ANGSTROM
    R2=                  1.89890745 ANGSTROM
    A1=                123.08527193 DEGREE

    分别对应以上的数值。


    可能原因是180算是一个非常不稳定的局部极小点,几何优化步长太小跳不出这个极小点(优化步长为1度,非常小),而在频率计算时因为转换成了笛卡尔坐标系进行displancement计算,displancement为0.01,改变幅度较大,在PES的角度上来看harmonic frequency model已经不适用,所以出现虚频。
    如果一定要计算角度为180的这个局部极小点,可以考虑多点计算PES,用多项式拟合振动波函数,再计算出harmonic frequency跟其他非简谐数值等

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓