计算化学公社

标题: 虚频无法消除，各种方法尝试失败 [打印本页]

作者
Author: 灰飞的旋律 时间: 2016-12-2 19:04
标题: 虚频无法消除，各种方法尝试失败
本帖最后由灰飞的旋律于 2016-12-2 19:06 编辑

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

作者
Author: 我本是个娃娃 时间: 2016-12-2 20:37
两种方法：
不使用BJ阻尼，用零阻尼
使用高斯E01

作者
Author: 灰飞的旋律 时间: 2016-12-2 20:58

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

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

作者
Author: 我本是个娃娃 时间: 2016-12-2 21:02

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

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

作者
Author: 灰飞的旋律 时间: 2016-12-2 21:23
本帖最后由灰飞的旋律于 2016-12-2 21:41 编辑

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

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

还有个问题想请教下您，如下图这个单空位。我在测试空位上对C2H3Cl的吸附能力
(, 下载次数 Times of downloads: 50)
测试结果和构型图下
(, 下载次数 Times of downloads: 47)
(, 下载次数 Times of downloads: 52)
最右侧一列是吸附能，单位kcal/mol。8和9都是物理吸附，其他的我就看不明白了。要是化学吸附的话，但按照吸附能计算方法算出来吸附能却是正的，显然热力学上无法发生的。那该怎么理解这个现象呢？

作者
Author: 我本是个娃娃 时间: 2016-12-3 11:20

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

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

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

作者
Author: 灰飞的旋律 时间: 2016-12-3 20:41

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

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

作者
Author: 我本是个娃娃 时间: 2016-12-3 21:04

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

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

作者
Author: 灰飞的旋律 时间: 2016-12-3 22:01

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

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

作者
Author: 我本是个娃娃 时间: 2016-12-3 22:09

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

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

作者
Author: 灰飞的旋律 时间: 2016-12-4 14:06

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

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

作者
Author: 我本是个娃娃 时间: 2016-12-4 14:25

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

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

作者
Author: 灰飞的旋律 时间: 2016-12-4 17:59

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

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

作者
Author: 我本是个娃娃 时间: 2016-12-4 18:14

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)