应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？
321/4返回列表1234下一页
查看: 3927  |  回复: 31
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

孙1986

金虫 (正式写手)

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

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

计算线性三原子分子的振动频率，弯曲振动频率是个虚频？
TIM图片20170804113550.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
临渊羡鱼，不如退而结网
1楼 2017-08-04 11:30:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wangf44

禁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
13楼2017-08-12 22:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

gengxin60
2楼2017-08-04 21:30:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+3 2017-08-06 15:09:08
把角度的自由度打开重新优化一下结构看看
赞一下 回复此楼
3楼2017-08-05 02:20:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙1986

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

好的，我试一下
回复此楼
临渊羡鱼，不如退而结网
4楼2017-08-06 15:09:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙1986

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

角度自由度是打开的
赞一下 回复此楼
临渊羡鱼，不如退而结网
5楼2017-08-06 15:12:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
孙1986: 金币+3 2017-08-08 11:29:08
引用回帖:
5楼: Originally posted by 孙1986 at 2017-08-06 15:12:44
角度自由度是打开的...

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

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

如果用的方法是DFT的话，个人经历过很多情况设置grid对几何优化及频率计算影响很大，可以选择更小的grid数值。
如果是WF办法的话，除了以上提到的问题暂时想不出其他的可能原因，你可以把input贴一下。
赞一下 回复此楼
6楼2017-08-08 03:16:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙1986

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

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

这是输入文件
计算线性三原子分子的振动频率，弯曲振动频率是个虚频？-1
TIM图片20170808113729.png
赞一下 回复此楼
临渊羡鱼，不如退而结网
7楼2017-08-08 11:34:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

121099

银虫 (小有名气)
★ ★ ★ ★ ★
小木虫: 金币+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跟其他非简谐数值等。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

孙1986
8楼2017-08-11 04:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙1986

金虫 (正式写手)
送红花一朵
引用回帖:
8楼: Originally posted by 121099 at 2017-08-11 04:43:28
用DZ优化了一下结构。
发现如果打开夹角至175以下，则可以找到能量更低的PES极小点，
DZ优化结果为
R1=                  1.20567667 ANGSTROM
R2=                  1.89890745 ANGSTROM
A1=                ...

这是一个排斥激发态的的结构，可能确实不好处理。但是按照楼主的讲解不太会处理，不过楼主的讲解倒是给了我一点启发，我先把角度稍微改一下，再试试。
赞一下 回复此楼
临渊羡鱼，不如退而结网
9楼2017-08-11 16:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangf44

禁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽

» 本帖已获得的红花（最新10朵）

孙1986
10楼2017-08-11 16:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孙1986 的主题更新
321/4返回列表1234下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工考研调剂 +9 孅華 2026-03-22 9/450 2026-03-25 13:09 by cmz0325
[考研] 材料调剂 +3 iwinso 2026-03-23 3/150 2026-03-25 11:29 by greychen00
[考研] 289求调剂 +9 怀瑾握瑜l 2026-03-20 9/450 2026-03-25 11:02 by userper
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +6 niko- 2026-03-22 6/300 2026-03-25 09:14 by aa331100
[考研] B区考研调剂 +4 yqdszhdap－ 2026-03-22 5/250 2026-03-25 08:51 by baoball
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 0854AI CV方向招收调剂 +3 章小鱼567 2026-03-23 3/150 2026-03-24 20:25 by 汪！？！
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考博] 26申博自荐 +3 whh869393 2026-03-24 3/150 2026-03-24 09:55 by 21018060
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 一志愿武理材料工程348求调剂 +6 ￣^￣゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 考研调剂 +3 呼呼？~+123456 2026-03-21 3/150 2026-03-21 20:04 by 无际的草原
[考研] 一志愿南大，0703化学，分数336，求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 279求调剂 +5 红衣隐官 2026-03-21 5/250 2026-03-21 14:59 by lature00
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 288求调剂 +16 于海海海海 2026-03-19 16/800 2026-03-20 22:28 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
信息提示
关闭
请填处理意见
关闭