计算化学公社

标题: PS4做SAPT能力分解的问题 [打印本页]

作者
Author: dgwfly 时间: 2021-11-7 17:32
标题: PS4做SAPT能力分解的问题
本帖最后由 dgwfly 于 2021-11-7 17:35 编辑

老师好，我是用SAPT计算一个氯离子与B27N27H18平板的相互作用,示意图如下。当氯离子与B27N27H18的距离为2.4
Å时，计算得到的sSAPT0结果如下：
Special recipe for scaled SAPT0 (see Manual):
Electrostatics sSAPT0       -83.55402120 [mEh]    -52.43093988 [kcal/mol] -219.37105244 [kJ/mol]
Exchange sSAPT0             157.00930734 [mEh]    98.52482783 [kcal/mol]    412.22787963 [kJ/mol]
Induction sSAPT0             1.28083274 [mEh]    0.80373468 [kcal/mol]    3.36282589 [kJ/mol]
Dispersion sSAPT0          -25.38472096 [mEh]    -15.92915289 [kcal/mol]    -66.64757569 [kJ/mol]
Total sSAPT0                   49.35139792 [mEh]    30.96846974 [kcal/mol]    129.57207740 [kJ/mol]

用scaled SAPT0方法得到的Induction energy为正值，很奇怪;用SAPT0得到的Induction energy 为-18.05994212，和scaled 方法的差别很大。当dimer之间距离比较大的时候，结果就比较正常。不知道是我输入文件的问题，还是sSAPT0方法不适用距离太近的dimer？？附件为输入和输出文件。

(, 下载次数 Times of downloads: 2)

(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 19)

(, 下载次数 Times of downloads: 22)

作者
Author: zjxitcc 时间: 2021-11-8 10:11
用更高级SAPT算一下，或者GKS-EDA算一下，以更高级计算级别的结果为参考答案再考虑。

作者
Author: snljty 时间: 2021-11-8 12:05

zjxitcc 发表于 2021-11-8 10:11
用更高级SAPT算一下，或者GKS-EDA算一下，以更高级计算级别的结果为参考答案再考虑。

请问这么大体系更高级SAPT算的动么？

作者
Author: zjxitcc 时间: 2021-11-8 13:09

snljty 发表于 2021-11-8 12:05
请问这么大体系更高级SAPT算的动么？

PSI4估计算不动

作者
Author: dgwfly 时间: 2021-12-22 13:25

zjxitcc 发表于 2021-11-8 10:11
用更高级SAPT算一下，或者GKS-EDA算一下，以更高级计算级别的结果为参考答案再考虑。

我用小体系测试了一下，就是氯离子和B3N3H6之间的相互作用，分别用sSAPT0- jun-cc-pdvz 和最贵的SAPT2+(3)δMP2/aug-cc-pVTZ，发现结果相差巨大，而且SAPT2的induction energy就很正常，不会出现正值的情况。看来sSAPT0算含离子的体系，可能是出现问题。但大体系，用SAPT2又算不动

，好难解决这个问题

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