24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yjcmwgk

禁虫 (文坛精英)

密度泛函·小卒

QC强帖: 38
应助: 165 (高中生)
贵宾: 16.999
金币: 553.8
散金: 39880
红花: 395
沙发: 2
帖子: 10406
在线: 1528.1小时
虫号: 448368
注册: 2007-11-01
专业: 理论和计算化学
管辖: 科研工具资源

[交流] 【调查】想做个Gaussian的FAQ了。大家一起参与意见吧（回帖者均可获得金币）

已经进入征集答案阶段。活动地址见：
http://muchong.com/bbs/viewthread.php?tid=1561636&fpage=1

[ Last edited by yjcmwgk on 2009-9-23 at 17:27 ]

回复此楼

» 猜你喜欢

26年秋季博士申请已经有0人回复
推荐给英语教学者的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有59人回复
物理化学论文润色/翻译怎么收费? 已经有62人回复
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》已经有32人回复
核磁分析软件MestReNova打开文件时报错已经有0人回复
在职博后不能申请博后基金了，那么在职博后意义何在？已经有2人回复
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生已经有0人回复
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）已经有0人回复
求助Cu2+1O的CIF文件(PDF: 05-0667) 已经有1人回复

1楼 2009-09-22 21:52:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

虚谦

至尊木虫 (正式写手)

应助: 1 (幼儿园)
金币: 10056.1
红花: 1
帖子: 841
在线: 64.1小时
虫号: 512837
注册: 2008-02-27
性别: GG
专业: 理论和计算化学

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+4,VIP+0): 10-3 15:08

引用回帖:

Originally posted by 3867826 at 2009-10-2 18:42:
1 消虚的方法和经验？
2 固定的方法？算频率是否固定
3 IOP命令的介绍下一些常用的？

1.在进行Gaussian计算时，经常会出现多个虚频个，每个都只有－30~－20左右，但是如果结构改动过大，则一不小心，虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试：nosymm；加大循环次数；提高收敛度；iop(1/8=1)等。通常对于势能面很平，较小的虚频很难消除。
  为了解决以上的这种现象,按如下的步骤一步步调整：  1）. 增加收敛标准
在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。  2）. 降低对称性
对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。  3）. 坐标中和法
如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中产生的原子的振动坐标加到相应的输入文件中。这样，重新计算，直到虚频没有。比较好的做法是分子结构输入用直角坐标，将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加，可以选择权重如1/3相加等。  4）. 其他方法
实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做一般的电子结构计算，这样，应当没有问题。  如果上面几步仍然无法解决问题，则有可能是因为计算还没有得到最稳定的结构，那么或者是分子结构有问题，或者是计算错了，再者就是游离出现在计算的范畴。  最后如果是寻找过渡态，如果只有一个虚频，而且该虚频对应某单键旋转震动，那么产生虚频的原因就是初始结构比较接近重叠构象。这时把优化过的结构绕振动方向旋转几十度再重新优化就有可能得到正确的结果。优化中等大小的有机分子（50-100原子）的时候，应该经常碰到这种情况。当出现两个以上虚频时，可以先看一下振动模式。如果多余的虚频，和想要的振动模式发生在不同原子上，可以先固定过渡态模式所对应的原子，然后做opt到minimal,一般可以把多余的虚频去掉，再放开做全优化。或许会有所收获。

2.关于怎么固定前面的回帖已经讲过了，你自己去看看！
3.iop命令你可以去下一个gaussian的资料，里面是关于iop命令的！google就可以找到！