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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？
321/4返回列表1234下一页
查看: 3925  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 303求调剂 +4 元夕元 2026-03-20 5/250 2026-03-25 00:20 by 元夕元
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 材料专硕找调剂 +5 哈哈哈吼吼吼哈 2026-03-23 5/250 2026-03-24 19:07 by 了了了了。。
[考研] 085602 289分求调剂 +5 WWW西西弗斯 2026-03-24 5/250 2026-03-24 18:51 by jhhcooi
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 求调剂 +6 研研，接电话 2026-03-24 7/350 2026-03-24 17:01 by barlinike
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 一志愿北京化工大学材料与化工 264分各科过A区国家线 +3 哈哈157349 2026-03-21 3/150 2026-03-24 14:11 by zhyzzh
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 求材料，环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123！
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 一志愿东华大学控制学硕320求调剂 +3 Grand777 2026-03-21 3/150 2026-03-21 19:23 by 简之-
[考研] 材料与化工（0856）304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 一志愿华南师大 070300（化学）304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
信息提示
关闭
请填处理意见
关闭