24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4935  |  回复: 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的回帖

小笨虫


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

强烈要求局长大人给我们开个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的回帖

小红豆

主管区长 (知名作家)


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

楼主厉害!
5楼2008-03-15 20:47:46
已阅   回复此楼   关注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的回帖

tianguixin521

金虫 (小有名气)


不错,值得一看
7楼2008-03-18 09:04:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxguoch

银虫 (小有名气)


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

讲的很好!!!
8楼2008-04-09 09:00:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingfei101


救人一命啊
9楼2008-05-17 13:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dolphin618

金虫 (正式写手)


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

支持一下,很有收获
10楼2008-10-14 07:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youyou421

木虫 (小有名气)


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

楼主,辛苦了!很有用啊!
11楼2008-11-28 09:41:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yyx19840628

木虫 (著名写手)


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

很好,很强大
12楼2008-12-01 14:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zzgyb 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 材料与化工调剂 10+11 下一站上岸@ 2026-04-10 36/1800 2026-04-11 10:26 by 89436494
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 285求调剂 +9 AZMK 2026-04-07 11/550 2026-04-10 15:24 by AZMK
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 青岛科技大学材料学院,环境学院调剂补录4月10日以前都可以 +3 1青科大。 2026-04-09 5/250 2026-04-10 09:58 by 翩翩一书生
[考研] 347材料专硕求调剂 +19 zj8215216 2026-04-06 19/950 2026-04-10 09:36 by 690616278
[考研] 322求调剂,08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 材料307分求大佬组收留 +17 Hll胡 2026-04-07 17/850 2026-04-09 10:53 by liuhuiying09
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 277求调剂 +4 考研调剂lxh 2026-04-06 6/300 2026-04-08 10:40 by 逆水乘风
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
信息提示
请填处理意见