虚频无法消除，各种方法尝试失败

灰飞的旋律

如图，我要优化一个石墨烯双空位的结构，输入文件如下

优化结束后，出现了-17的虚频。我分别尝试了opt=tight、opt=tight和int=ultrafine组合、根据振动方向修改做标、手动大幅度调整边缘局部结构。但以上方法优化出来的结构虚频仍未减小或者消失。打算尝试以优化后的有虚频的结构做初猜继续calcall，还在计算中，但可能没啥效果。
请各位帮我分析下原因，不胜感谢！！！
另外附上石墨烯的输出文件。 DVG-1-1.log (2.57 MB, 下载次数 Times of downloads: 3)

2016-12-2 19:06 上传 Uploaded

点击下载
Click to download

我本是个娃娃

两种方法：
不使用BJ阻尼，用零阻尼
使用高斯E01

灰飞的旋律

我本是个娃娃 发表于 2016-12-2 20:37
两种方法：
不使用BJ阻尼，用零阻尼
使用高斯E01

0阻尼也没有用，试过了。E01为什么可能有效果呢？？

我本是个娃娃

灰飞的旋律 发表于 2016-12-2 20:58
0阻尼也没有用，试过了。E01为什么可能有效果呢？？

修复了BJ阻尼的一些BUG，看样子你的缺陷不可能保持平面了

灰飞的旋律

我本是个娃娃 发表于 2016-12-2 21:02
修复了BJ阻尼的一些BUG，看样子你的缺陷不可能保持平面了

谢谢您。我找找E01版本的试试

还有个问题想请教下您，如下图这个单空位。我在测试空位上对C2H3Cl的吸附能力

测试结果和构型图下


最右侧一列是吸附能，单位kcal/mol。8和9都是物理吸附，其他的我就看不明白了。要是化学吸附的话，但按照吸附能计算方法算出来吸附能却是正的，显然热力学上无法发生的。那该怎么理解这个现象呢？

我本是个娃娃

灰飞的旋律 发表于 2016-12-2 21:23
谢谢您。我找找E01版本的试试

还有个问题想请教下您，如下图这个单空位。我在测试空位上对C2H3Cl的吸 ...

这种吸附方式热力学不能实现呗，这个没什么好解释的

灰飞的旋律

我本是个娃娃 发表于 2016-12-3 11:20
这种吸附方式热力学不能实现呗，这个没什么好解释的

但是我从过渡态做IRC也会走出来1.4也就是左上角这种吸附方式。反应是C2H2和HCl加成生成氯乙烯。C2H2的这种两个C原子成键氮吸附方式吸附能虽然是负的，但也很低还不到-10。和HCll共吸附后，共吸附能也就负十几。

我本是个娃娃

灰飞的旋律 发表于 2016-12-3 20:41
但是我从过渡态做IRC也会走出来1.4也就是左上角这种吸附方式。反应是C2H2和HCl加成生成氯乙烯。C2H2的这 ...

结合并不是特别的稳定……

灰飞的旋律

我本是个娃娃 发表于 2016-12-3 21:04
结合并不是特别的稳定……

那我是不是可以理解为：C2H2单吸附不稳定；和HCl共吸附，能使吸附结构更加稳定；而生成的产物C2H3Cl却极易脱附？？

我本是个娃娃

灰飞的旋律 发表于 2016-12-3 22:01
那我是不是可以理解为：C2H2单吸附不稳定；和HCl共吸附，能使吸附结构更加稳定；而生成的产物C2H3Cl却极 ...

如果吸附能为正，热力学上是不稳定的，即使能得到稳定构型，也只是势能面上的一个势阱，而非极小点

灰飞的旋律

我本是个娃娃 发表于 2016-12-3 22:09
如果吸附能为正，热力学上是不稳定的，即使能得到稳定构型，也只是势能面上的一个势阱，而非极小点

那我这个反应C2H2+HCl——C2H3Cl也可以正常进行吧？？并且能得到产物C2HCl？

我本是个娃娃

灰飞的旋律 发表于 2016-12-4 14:06
那我这个反应C2H2+HCl——C2H3Cl也可以正常进行吧？？并且能得到产物C2HCl？

反应是可以的，可是后一个产物是什么？

灰飞的旋律

我本是个娃娃 发表于 2016-12-4 14:25
反应是可以的，可是后一个产物是什么？

不好意思写错了，产物是生成C2H3Cl。那么这个反应路径生成的产物很容易脱附？

我本是个娃娃

灰飞的旋律 发表于 2016-12-4 17:59
不好意思写错了，产物是生成C2H3Cl。那么这个反应路径生成的产物很容易脱附？

如果没有稳定的结合能，当然容易脱附