计算化学公社

标题: 溶剂中BSSE校正出现问题 [打印本页]

作者
Author: captain 时间: 2016-4-14 20:58
标题: 溶剂中BSSE校正出现问题

两分子相互作用能的计算

气相中直接用 counterpoise
不经BSSE校正这俩分子相互作用能 -147.46 kcal/mol
经过BSSE校正这俩分子相互作用能 -136.15 kcal/mol

实际上(经BSSE)相互作用能没有算出来的(不经BSSE)那么负很正常

液相水中不能直接用counterpoise 所以设定Ghost原子手动计算

不经BSSE校正这俩分子相互作用能 -67.39 kcal/mol
经过BSSE校正这俩分子相互作用能 -82.71 kcal/mol

怎么经过BSSE校正相互作用能更负了呢?跟溶剂有关系吗?

请大神们指点迷津！

作者
Author: sobereva 时间: 2016-4-14 21:38
溶剂模型下计算相互作用能没有严格的办法考虑BSSE问题，相对可行的做法是使用气相下得到的BSSE校正能，加到溶剂下相互作用能上去。
溶剂模型下用Bq不合适，因为Bq原子在溶剂模型计算时也会产生相应的孔洞，而不仅仅是表现相邻原子的基组了，这会导致BSSE校正能成为无意义的负值。

作者
Author: captain 时间: 2016-4-14 21:56

sobereva 发表于 2016-4-14 21:38
溶剂模型下计算相互作用能没有严格的办法考虑BSSE问题，相对可行的做法是使用气相下得到的BSSE校正能，加到 ...

原来如此明白了老师

作者
Author: jlm 时间: 2023-7-4 21:13

sobereva 发表于 2016-4-14 21:38
溶剂模型下计算相互作用能没有严格的办法考虑BSSE问题，相对可行的做法是使用气相下得到的BSSE校正能，加到 ...

老师，我想一下，溶剂下的相互作用能怎么去算呢？

作者
Author: wzkchem5 时间: 2023-7-4 21:55

jlm 发表于 2023-7-4 14:13
老师，我想一下，溶剂下的相互作用能怎么去算呢？

sob老师已经说了一种方法了。其他几个可能的方法：
（1）把基组选得足够大，使得不用做Counterpoise结果也可以用
（2）做两个不同基组的计算，外推到基组极限
（3）加F12校正
（4）用gCP代替Counterpoise来做BSSE校正，gCP可以和溶剂一起用，缺点是gCP的参数只对某些特定基组有定义，如果你用的基组gCP不支持的话就算不了，或者能计算但结果不可靠

作者
Author: jlm 时间: 2023-7-5 16:13

wzkchem5 发表于 2023-7-4 21:55
sob老师已经说了一种方法了。其他几个可能的方法：
（1）把基组选得足够大，使得不用做Counterpoise结果 ...

好的，谢谢老师

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