24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4962  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +20 想上岸呀!! 2026-04-12 23/1150 2026-04-17 10:50 by cuisz
[考研] 收到复试调剂但是去不了 +6 小蜗牛* 2026-04-16 6/300 2026-04-17 10:05 by 涵竹刘
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 求调剂 +8 小聂爱学习 2026-04-16 10/500 2026-04-16 22:06 by 1shin_ichi
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 305求调剂 +8 玛卡巴卡boom 2026-04-11 8/400 2026-04-14 09:04 by pengliang8036
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 303求调剂 +14 SereinQ 2026-04-10 15/750 2026-04-11 20:43 by 蓝云思雨
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
信息提示
请填处理意见