24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

孙1986

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 1046.9
散金: 1087
红花: 19
帖子: 455
在线: 243.1小时
虫号: 1405621
注册: 2011-09-17
专业: 原子和分子物理

[交流] 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？已有9人参与

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

TIM图片20170804113550.png

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

临渊羡鱼，不如退而结网

1楼2017-08-04 11:30:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

121099

银虫 (小有名气)

应助: 6 (幼儿园)
金币: 327.9
散金: 3
红花: 2
帖子: 76
在线: 76.4小时
虫号: 1299508
注册: 2011-05-18
专业: 理论和计算化学

★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+4 2017-08-11 16:16:48

引用回帖:

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跟其他非简谐数值等。