| 查看: 4808 | 回复: 17 | |||
| 【奖励】 本帖被评价13次，作者zzgyb增加金币 12.5 个 | |||
| 当前主题已经存档。 | |||
[资源]
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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 计算化学 |
» 猜你喜欢
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
物理化学论文润色/翻译怎么收费?
已经有101人回复
-
请问四氢呋喃溶解的聚合物用甲醇沉淀时，如何使沉淀过程加速？
已经有2人回复
-
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材
已经有0人回复
-
分子动力学模拟合作需求
已经有2人回复
-
诚邀有志之士加入江苏大学环境学院合成生物学课题组！(长期接收学硕专硕、硕博连读)
已经有16人回复
-
一志愿南大理论与计算化学方向（070300）342分，求调剂
已经有1人回复
-
欢迎化学、化工、物理、材料、或计算机背景的考研咨询
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2007-12-07 21:56:29
3楼2007-12-07 21:59:25
4楼2008-03-15 20:32:30
5楼2008-03-15 20:47:46
6楼2008-03-16 12:11:01
7楼2008-03-18 09:04:04
8楼2008-04-09 09:00:57
9楼2008-05-17 13:05:40
10楼2008-10-14 07:50:30
11楼2008-11-28 09:41:34
12楼2008-12-01 14:04:33
