24小时热门版块排行榜    

CyRhmU.jpeg
查看: 7636  |  回复: 18
【奖励】 本帖被评价14次,作者cjl1761增加金币 10.6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cjl1761

金虫 (小有名气)


[资源] 在进行高斯计算时如何消去虚频(转载)

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


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

  如果上面几步仍然无法解决问题,则有可能是因为计算还没有得到最稳定的结构,那么或者是分子结构有问题,或者是计算错了,再者就是游离出现在计算的范畴。

  最后如果是寻找过渡态,如果只有一个虚频,而且该虚频对应某单键旋转震动,那么产生虚频的原因就是初始结构比较接近重叠构象。这时把优化过的结构绕振动方向旋转几十度再重新优化就有可能得到正确的结果。优化中等大小的有机分子(50-100原子)的时候,应该经常碰到这种情况。

  当出现两个以上虚频时,可以先看一下振动模式。如果多余的虚频,和想要的振动模式发生在不同原子上,可以先固定过渡态模式所对应的原子,然后做opt到minimal,一般可以把多余的虚频去掉,再放开做全优化。或许会有所收获。
回复此楼

» 收录本帖的淘帖专辑推荐

过渡态计算 催化与分子模拟 量化辅导 量化软件学习
量化计算有意义的讨论 光谱 量化基础 Gaussian
反应过程

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daniell251

木虫 (著名写手)


引用回帖:
10楼: Originally posted by woxiangfei at 2012-07-08 14:07:30
校区虚频的方法有好多的!在进行Gaussian计算时,经常会出现多个虚频个,但是如果结构改动过大,则一不小心,虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试:nosymm;加大循环次数;提高收敛 ...

相加后的构型要不要再进行opt的优化? 还是直接认为这个构型就是优化了的构型,直接拿这个构型进行frq的分析呢 谢谢了
12楼2012-08-29 09:03:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

milanforza

新虫 (初入文坛)


★★★★★ 五星级,优秀推荐

唉,高斯菜鸟,一点都看不懂啊
5楼2012-04-11 16:58:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao_niao

铁虫 (小有名气)


★★★★★ 五星级,优秀推荐

楼主,最后一个看不懂,怎么做到“可以先固定过渡态模式所对应的原子,然后做opt到minimal,”?
8楼2012-06-25 22:29:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woxiangfei

木虫 (文坛精英)


校区虚频的方法有好多的!在进行Gaussian计算时,经常会出现多个虚频个,但是如果结构改动过大,则一不小心,虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试:nosymm;加大循环次数;提高收敛度;iop(1/8=1)等。通常对于势能面很平,较小的虚频很难消除。
消虚频的最有效的方法,就是根据虚频对应的振动位移(输出文件每一个频率下面都有一个类似坐标的3列数值,给出的就是每个原子在这个振动频率时候的移动方式),把这些位移,乘以一个适当大小的因子,直接和原来优化好的平衡构型相加,就得到新的构型了。
10楼2012-07-09 10:07:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2012-05-14 09:31   回复  
五星好评  
2012-07-02 16:35   回复  
五星好评  
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复(可上传附件)
信息提示
请填处理意见