24小时热门版块排行榜    

查看: 4813  |  回复: 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 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 求调剂 +5 Zhangbod 2026-03-21 7/350 2026-03-22 13:13 by Zhangbod
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 北科281学硕材料求调剂 +5 tcxiaoxx 2026-03-20 5/250 2026-03-20 21:35 by laoshidan
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见