24小时热门版块排行榜    

查看: 249  |  回复: 2
当前主题已经存档。

Jldbaobei

木虫 (著名写手)

按时毕业找到好工作

[交流] 【求助】计算过渡态得出的虚频一直很小

计算过渡态时出现两个虚频
调整计算后剩下一个虚频
但是接下来不管怎么样对键长、键角、二面角做小的调整
产生的虚频总是在20左右—50之间
命令行中也加入了nosymm scf=tight 关键字
请问虚频的调整有什么技巧吗?
怎么我弄来弄去它还是在这个区间里转啊

[ Last edited by lei0736 on 2008-11-29 at 11:18 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小不点点

木虫 (著名写手)

★ ★
lei0736(金币+2,VIP+0):谢谢
不知道你有没有根据g-view给出的振动进行调整。还有一个可能就是该结构在该区的势能面很平缓,这就需要你对结构进行较大范围的调整。
别胡闹了
2楼2008-11-16 22:34:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张文艳

交流

★ ★
lei0736(金币+2,VIP+0):谢谢
我也是看的别人说的方法,你试试有没有用吧!
在进行Gaussian计算时,经常会出现多个虚频个,每个都只有-30~-20左右,但是如果结构改动过大,则一不小心,虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试:nosymm;加大循环次数;提高收敛度;iop(1/8=1)等。通常对于势能面很平,较小的虚频很难消除。

  为了解决以上的这种现象,按如下的步骤一步步调整:

   1). 在优化时采用Scf(tight)的选项,增加收敛的标准。再去计算频率。如果还有虚频,参见下一步。

   2). 对称性的影响,很多情况下的虚频是由于分子本身的对称性造成的。这样,在优化时,如果必要,要将对称性降低,还有,输入文件有时是用内坐标。建议如果有虚频的话,将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。
3). 如果上述方法还有虚频,看一下虚频,找到强度较大的,将在频率中产生的原子的振动坐标加到相应的输入文件中。这样,重新计算。直到虚频没有。比较好的做法是分子结构输入用直角坐标,将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加,可以选择权重如1/3相加等。
   4). 实际上,如果分子柔性较大,很难找到最低点,这是电子结构计算的问题,这种情况下,需要动力学的 东西,用构象搜寻的办法解决。如:模拟退火,最陡下降法,淬火法等。将得到的能量最低的构象做一般的电子结构计算,这样,应当没有问题。
当出现两个以上虚频时,可以先看一下振动模式。如果多余的虚频,和想要的振动模式发生在不同原子上,可以先固定过渡态模式所对应的原子,然后做opt到minimal,一般可以把多余的虚频去掉,再放开做全优化。
3楼2008-11-17 09:09:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Jldbaobei 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见