24小时热门版块排行榜    

查看: 4823  |  回复: 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的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 310求调剂 +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 308求调剂 +3 墨墨漠 2026-03-21 3/150 2026-03-22 16:54 by i_cooler
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 070300化学319求调剂 +7 锦鲤0909 2026-03-17 7/350 2026-03-21 03:46 by JourneyLucky
[考研] 南京大学化学376求调剂 +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
信息提示
请填处理意见