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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？
321/4返回列表1234下一页
查看: 3993  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然面上和省基金B类撒花 +17 花田半亩～白 2026-04-21 17/850 2026-04-23 07:30 by john198197
[考研] 有没有学校收留 +3 蒋昌鹏qtj 2026-04-20 3/150 2026-04-22 20:25 by 学员JpLReM
[考研] 312求调剂 +3 山河似你温柔 2026-04-22 3/150 2026-04-22 20:17 by 学员JpLReM
[考博] 华师大读博 +3 xq83 2026-04-22 5/250 2026-04-22 10:42 by xq83
[论文投稿] 急需审稿人！！！ +3 陆小果画大饼 2026-04-21 3/150 2026-04-21 23:54 by jzy_123456
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 收到复试调剂但是去不了 +8 小蜗牛* 2026-04-16 8/400 2026-04-18 11:15 by zixin2025
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
信息提示
关闭
请填处理意见
关闭