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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Molpro/Molcas » 计算线性三原子分子的振动频率，弯曲振动频率是个虚频？
321/4返回列表1234下一页
查看: 3975  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 291求调剂 +9 关忆北. 2026-04-14 9/450 2026-04-16 22:49 by cfdbai
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考研] 26药学专硕105500求调剂 +6 喽哈加油 2026-04-13 7/350 2026-04-16 14:31 by zhouxiaoyu
[考研] 一志愿A区211，22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 调剂求收留 +34 果然有我 2026-04-10 35/1750 2026-04-15 13:05 by 西北望—风沙
[考研] 材料工程281还有调剂机会吗 +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by 西北望—风沙
[考研] 化工学硕294分，求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 药学305求调剂 +10 玛卡巴卡boom 2026-04-10 10/500 2026-04-14 15:55 by zs92450
[考研] 求调剂 +12 何气正 2026-04-13 13/650 2026-04-14 14:47 by zs92450
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 296求调剂 +14 汪！？！ 2026-04-10 16/800 2026-04-12 10:48 by zhouyuwinner
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
信息提示
关闭
请填处理意见
关闭