计算化学公社

标题: 高斯优化结构有很小的虚频 [打印本页]

作者
Author: huahua0621 时间: 2022-4-26 20:38
标题: 高斯优化结构有很小的虚频
老师好，光催化算反应路径，用B3LYP-D3/6-31G**&SDD(Ir)优化结构，在处理反应过程中有新的物质比如A加上到B的时候，都是先优化A和B在一起的复合物，然后把A-B拉开10-20埃的距离，优化，用两种结构的能量差作为结合能，画在反应路径里。问题是每次拉开10-20埃的距离优化的结构总有虚频，虚频不大，只有-0.1左右最大-1.2，试着在坐标上加上虚频下面xyz的方向再进行优化，还是没有消除，尝试了加0.1倍的也没有用。想问问：1.这样算结合能在反应路径上表示对不对，因为如果拉远10埃两个分子会重新回到复合物的距离；2.这个虚频会不会影响结果，是不是一定要消除，怎么消除；3.反应过程两个物质靠近的过程中可能有电荷转移，比如本身是 A*和B，但是优化复合物后就成了A和B*,拉开后还是A和B*,不是最初的A*和B,不能描述之前的状态怎么办?谢谢老师

作者
Author: sobereva 时间: 2022-4-27 05:20
结合能根本就不是这么算的，你的做法是严重的歪门邪道
应当复合物能量减去两个单体各自单独计算的能量

作者
Author: huahua0621 时间: 2022-4-27 09:41
老师可能是我表述不清，就是在反应路径的过程中加入一个物质，也是应该这样复合物减去两个单体能量算能量差吗

作者
Author: ChemG 时间: 2022-4-27 10:28
这种做法确实比较离谱。你可以每步分开考虑呀，比如第一步反应的产物C+D，C继续反应而D脱离出体系。那么你第二步反应参考的初始相对能量还取第一步产物的相对能量，就跟D无关了，再考虑第二部反应的过渡态及产物能量和第一步产物能量（也是第二部反应复合物的能量零点）的相对值。也就是每一步反应是单独考虑的，要是每步都积累一点分子，反应十来步怕是成了团簇了。

作者
Author: KangJiaxin 时间: 2022-4-27 11:29

huahua0621 发表于 2022-4-27 09:41
老师可能是我表述不清，就是在反应路径的过程中加入一个物质，也是应该这样复合物减去两个单体能量算能量 ...

是，这和你过程中开始一开始没有关系。算两个物质结合能应当是复合物减去两个单体的和，你拉的再远计算的时候还是按照一个整体去算的。

作者
Author: huahua0621 时间: 2022-4-27 11:43
就是因为我的体系有单电子，所以拉远和加和的区别有点大

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