24小时热门版块排行榜    

CyRhmU.jpeg
查看: 974  |  回复: 9
当前主题已经存档。

gaussian08

银虫 (小有名气)

[交流] 【讨论】虚频的问题

我用MP2方法,对两个小分子进行优化,成功优化后进行频率计算(计算路径部分:opt, freq),但是经过几次计算都发现有1到2个虚频,freq小于-20波数,我不知道是否是真正的优化结构,但是我修改了几次初始结构,都得到了小的虚频。希望有做频率研究的高手给予指导。
回复此楼
有你,有我
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

3867826

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 12-1 22:03
引用回帖:
Originally posted by gaussian08 at 2009-11-28 22:21:
我用MP2方法,对两个小分子进行优化,成功优化后进行频率计算(计算路径部分:opt, freq),但是经过几次计算都发现有1到2个虚频,freq小于-20波数,我不知道是否是真正的优化结构,但是我修改了几次初始结构, ...

但是经过几次计算都发现有1到2个虚频  你固定了没》?没固定你可以内坐标
+上虚频坐标进行消虚 。。或者改变初始构型 再计算  小分子应该可以消掉。。。的却也有很小的消不掉
拒稿,终于体会到了!乱写接受,认真写拒稿,这就是2和4的区别吗?
2楼2009-11-29 15:46:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaussian08

银虫 (小有名气)

楼上的意见很中肯,但是我现在用B3LYP/6-31G进行优化后没有得到虚频,应该说是最优结构,但是在此结构的基础上用MP2计算单点能和频率时,发现有虚频出现,怎么回事?》
有你,有我
3楼2009-11-29 16:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

caohuiming

铜虫 (初入文坛)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 11-29 19:58
引用supi的原话

首先,
1在优化时采用Scf(tight)的选项,增加收敛的标准。再去计算频率。如果还有虚频,参见下一步。
2. 对称性的影响,很多情况下的虚频是由于分子本身的对称性造成的。这样,在优化时,如果必要,要将对称性降低,还有,输入文件有时是用内坐标。建议如果有虚频的话,将内坐标改成直角坐标优化。
3.如果上述方法还有虚频,看一下虚频,找到强度较大的,将在频率中产生的原子的振动坐标加到相应的输入文件中。这样,重新计算。直到虚频没有。
4 实际上,如果分子柔性较大,很难找到最低点,这是电子结构计算的问题,这种情况下,需要动力学的 东西,用构象搜寻的办法解决。如:模拟退火,最陡下降法,淬火法等。将得到的能量最低的构象做一般的电子结构计算,这样,应当没有问题。

不要讲你还没有得到最稳定的结构,那么,是你的分子有问题,要么计算错了,要么就是游离出现代计算的范畴

[ Last edited by yjcmwgk on 2009-11-29 at 19:58 ]
4楼2009-11-29 16:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gonemrc

捐助贵宾 (正式写手)

有个问题我有点奇怪,为什么用b3lyp优化后,再用MP2来算能量与频率?为什么不继续用b3lyp? 我的拙见就是这两个方法应该取的不同的近似吧?
5楼2009-11-30 11:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 12-1 22:03
1.可以先使用hf优化计算频率,然后将chk文件保存为gjf文件再改为MP2下计算
2.初始结构你没有调整到合适的位置,继续调整,然后先小基组计算,如果没虚频,再大基组,可以节省不少时间
6楼2009-11-30 15:14:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaochao85

铜虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 12-1 22:03
一般当你改变方法时候,很容易虚频,
解决方法就是上面四楼提到的那些方法,你可以单个试试,也可以组合试试,
消虚频没有最好的方法,对不同的体系得慢慢试。
MP2适合算有氢键作用的体系。
7楼2009-12-01 10:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 12-1 22:03
算法和基组不同,自然优化的结果就不同啊。这很正常啊。
你用b3lyp优化的结构趋于收敛了,但是用mp2算就应该重新优化——虽然大体上
应该不会有太大变化,但是结构无法保证完全一致。如果不重新优化直接算频率,
有虚频是正常的。
虽然通常习惯上用b3lyp优化,然后用mp2再单点算能量,那是因为b3lyp优化收敛
速度快,而mp2能量更精确,但那是因为能量变化不大,所以近似的把b3lyp优化的
结构当成最稳定状态,但是用mp2算频率就看出来不一致了。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
8楼2009-12-01 12:29:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaussian08

银虫 (小有名气)

楼上说的很对,但是我用较高的方法,计算机能力达不到。
那讨论频率只能用较低的方法计算出的,这样做对频率的讨论有影响吗?用较低方法计算的频率准确吗?
有你,有我
9楼2009-12-01 13:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT


小木虫(金币+0.5):给个红包,谢谢回帖交流
根据你的机器能力,算到什么水平就是什么样啊。本来频率计算就不是那么精确么。
不然那个所谓的校正因子是干什么的啊。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
10楼2009-12-01 13:40:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gaussian08 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见