高斯优化结构有很小的虚频

huahua0621 · 发表于 Post on 2022-4-26 20:38:50

老师好，光催化算反应路径，用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,不能描述之前的状态怎么办?谢谢老师

sobereva · 发表于 Post on 2022-4-27 05:20:12

结合能根本就不是这么算的，你的做法是严重的歪门邪道
应当复合物能量减去两个单体各自单独计算的能量

huahua0621 · 发表于 Post on 2022-4-27 09:41:55

老师可能是我表述不清，就是在反应路径的过程中加入一个物质，也是应该这样复合物减去两个单体能量算能量差吗

ChemG · 发表于 Post on 2022-4-27 10:28:42

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

KangJiaxin · 发表于 Post on 2022-4-27 11:29:24

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

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

huahua0621 · 发表于 Post on 2022-4-27 11:43:31

就是因为我的体系有单电子，所以拉远和加和的区别有点大

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 高斯优化结构有很小的虚频

浏览过的版块