| 查看: 7632 | 回复: 18 | |||||||||||
| 【奖励】 本帖被评价14次，作者cjl1761增加金币 10.6 个 | |||||||||||
[资源]
在进行高斯计算时如何消去虚频（转载）
|
|||||||||||
|
在进行Gaussian计算时，经常会出现多个虚频个，每个都只有－30~－20左右，但是如果结构改动过大，则一不小心，虚频就会更大。消除虚频的主要方法是改变构型。其次在计算上还可以尝试：nosymm；加大循环次数；提高收敛度；iop(1/8=1)等。通常对于势能面很平，较小的虚频很难消除。 为了解决以上的这种现象,按如下的步骤一步步调整： 1）. 在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。 2）. 对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。但是有时这一方法也并没有多大作用。 3）. 如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到虚频没有。比较好的做法是分子结构输入用直角坐标，将这些产生虚频的振动坐标直接与原始坐标相加。加的时候不一定要将坐标直接相加，可以选择权重如1/3相加等。 4）. 实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做一般的电子结构计算，这样，应当没有问题。 如果上面几步仍然无法解决问题，则有可能是因为计算还没有得到最稳定的结构，那么或者是分子结构有问题，或者是计算错了，再者就是游离出现在计算的范畴。 最后如果是寻找过渡态，如果只有一个虚频，而且该虚频对应某单键旋转震动，那么产生虚频的原因就是初始结构比较接近重叠构象。这时把优化过的结构绕振动方向旋转几十度再重新优化就有可能得到正确的结果。优化中等大小的有机分子（50-100原子）的时候，应该经常碰到这种情况。 当出现两个以上虚频时，可以先看一下振动模式。如果多余的虚频，和想要的振动模式发生在不同原子上，可以先固定过渡态模式所对应的原子，然后做opt到minimal,一般可以把多余的虚频去掉，再放开做全优化。或许会有所收获。 |
回复此楼» 收录本帖的淘帖专辑推荐
 过渡态计算 | 催化与分子模拟 | 量化辅导 | 量化软件学习 |
| 量化计算有意义的讨论 | 光谱 | 量化基础 | Gaussian |
| 反应过程 |
» 猜你喜欢
厦门大学招专业学位博士，申请考核制，化学、材料、物理、计算机方向
已经有0人回复
-
有偿调查问卷，一元一份，有意请私我（不要豆包）
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有178人回复
-
求助ORCA计算单点能的输入输出文件
已经有0人回复
-
中教金源的流动池
已经有0人回复
-
RSC ADV状态问题
已经有4人回复
-
湖南大学招收锂离子电池、离子液体、高分子与凝胶材料、电催化方向博士生多名
已经有7人回复
-
南京大学蔡亮课题组诚招2026年申请-考核制博士生2-3名（电解水制氢，XAFS谱学等）
已经有32人回复
-
求助各位大佬配位化学学习教材！！！
已经有0人回复
-
大连海事大学国家级人才团队2026年博士研究生招生启事
已经有9人回复
-
内蒙古大学青年长江学者王蕾教授课题组2026年博士招聘，2-3人
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
高斯计算中zORA方法是否可以使用
已经有9人回复- 高斯计算过渡态时的能垒总是负值
已经有13人回复
- 【整理】今天整理硬盘上的资料，偶然发现的好玩的小东西就直接发上来吧
已经有105人回复
- 【求助】如何用高斯计算一个分子中的某一部分HOMO和LUMO？
已经有7人回复
- 【求助】高斯计算时无缘无故中断的问题
已经有8人回复
- 【求助】高斯频率计算死掉的原因
已经有7人回复
- 【求助】高斯计算过渡态不收敛，请问如何解决？
已经有8人回复
- 【讨论】直接进行频率计算有意义吗？
已经有9人回复
5楼2012-04-11 16:58:35
8楼2012-06-25 22:29:25
10楼2012-07-09 10:07:30
12楼2012-08-29 09:03:27
15楼2013-01-18 13:55:00
3楼2012-03-17 09:38:12
13楼2012-09-19 13:53:44
14楼2013-01-16 09:57:44
16楼2013-04-19 11:32:03
17楼2013-05-02 21:57:26
18楼2017-06-09 18:24:47
简单回复
nnipp2楼
2012-02-24 17:39
回复
五星好评 谢谢分享
gonghj4楼
2012-03-30 09:07
回复
五星好评 顶
2012-05-14 09:31
回复
五星好评
2012-06-01 16:46
回复
五星好评 
2012-07-02 16:35
回复
五星好评 
chenzz71811楼
2012-07-09 11:39
回复
五星好评 感谢分享！
yex10235041919楼
2018-01-07 13:45
回复
五星好评