24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4945  |  回复: 17
【奖励】 本帖被评价13次,作者zzgyb增加金币 12.5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zzgyb

荣誉版主 (文坛精英)


[资源] Gaussian98中消去虚频的技巧

Gaussian98中消去虚频的技巧


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

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

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

  2). 对称性的影响,很多情况下的虚频是由于分子本身的对称性造成的。这样,在优化时,如果必要,要将对称性降低,还有,输入文件有时是用内坐标。建议如果有虚频的话,将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。
  关于如何将内坐标改成直角坐标,可以直接看输出文件(例子如下):
  Full point group CS NOp 2
  Largest Abelian subgroup CS NOp 2
  Largest concise Abelian subgroup CS NOp 2
   Standard orientation:
  ---------------------------------------------------------------------
  Center Atomic Atomic Coordinates (Angstroms)
  Number Number Type X Y Z
  ---------------------------------------------------------------------
   1 46 0 -0.242506 1.028743 2.204311
   2 46 0 -0.476137 2.619863 0.000000
   3 46 0 -0.242506 1.028743 -2.204311
   4 46 0 -0.242506 -1.618028 -1.382208
   5 46 0 -0.242506 -1.618028 1.382208
   6 46 0 -1.725240 0.349678 0.000000
  Coordinates这一部分就是,如果没有,看Input orientation也行。

  3). 如果上述方法还有虚频,看一下虚频,找到强度较大的,将在频率中产生的原子的振动坐标加到相应的输入文件中。这样,重新计算。直到虚频没有。比较好的做法是分子结构输入用直角坐标,将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加,可以选择权重如1/3相加等。
  4). 实际上,如果分子柔性较大,很难找到最低点,这是电子结构计算的问题,这种情况下,需要动力学的 东西,用构象搜寻的办法解决。如:模拟退火,最陡下降法,淬火法等。将得到的能量最低的构象做一般的电子结构计算,这样,应当没有问题。

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

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

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

[ Last edited by zzgyb on 2008-1-10 at 12:17 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaoyan2418

木虫 (正式写手)


引用回帖:
Originally posted by 小笨虫 at 2007-12-7 21:59:

局长,给我们开个Gaussian的讲座吧

支持了.希望有一个讲座了.
6楼2008-03-16 12:11:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

小笨虫


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

强烈要求局长大人给我们开个Gaussian讲座
2楼2007-12-07 21:56:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小笨虫


1


局长,给我们开个Gaussian的讲座吧
3楼2007-12-07 21:59:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

虚谦

至尊木虫 (正式写手)


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

好贴啊,lz辛苦了~!真的很有用啊~!
4楼2008-03-15 20:32:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 调剂 +11 月@163.com 2026-04-11 11/550 2026-04-13 22:52 by pies112
[考研] 310求调剂 +15 666真好 2026-04-11 17/850 2026-04-13 22:27 by pies112
[考研] 327求调剂 +21 Xxjc1107. 2026-04-13 23/1150 2026-04-13 21:53 by bljnqdcc
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 求调剂 +9 吃口冰激凌 2026-04-07 9/450 2026-04-13 08:46 by 紫曦紫棋
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 085404,285分求调剂 +12 薇薇考研 2026-04-07 14/700 2026-04-09 23:10 by parmtree
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
信息提示
请填处理意见