计算化学公社

标题: PSI4能量分解任务出错以及与GAUSSIAN结果相差很大的 [打印本页]

作者
Author: LIEL 时间: 2022-5-28 12:27
标题: PSI4能量分解任务出错以及与GAUSSIAN结果相差很大的
各位老师好，我是计算化学的菜鸟，最近用gaussian算结合能以及PSI4能量分解时遇到一些有问题，很困惑，想请教一下各位老师一些相关问题。最近在做聚吡咯（五个吡咯分子单元）与铬酸根的结合能计算，由于是氧化形态的聚吡咯，因此在GAUSSIAN上给聚吡咯单独加了一个正电荷。首先用B3LYP-D3(BJ)/6-311g(d,p)分别对铬酸根，聚吡咯以及两者的结合做了几何优化，再用m062x/6-311++g(d,p)进行单点能计算，得出的结合能为-11 kcal/mol。但是该结果与文献中相差很远，并且这种静电结合的能量不应该这么低（可能m062x不适合Cr，但是做其他阴离子时仍然是这个问题）。此外，进一步按sob老师的方法做PSI4进行能量分解http://sobereva.com/526，由于是开壳层计算，首先使用的是sSAPT0/jun-cc-pVDZ，但是显示对Cr没有定义，因此改用SAPT0/def2-QZVPD，但总是收敛不了，加关键词如：soscf true和soscf_max_iter 30也不行。随即采用小机组def2-SVPD，算出来的结合能（-86.13kcal/mol）虽然与gaussian算出来的结合能相差很大，但是较为符合规律。此外，想进一步优化机组以得到更准确的结果，将机组提升为def2-TZVPPD，但是出现了如下图及附件所呈现的错误。因此，想请教各位老师3个问题，希望各位老师指正：
1. 为什么gaussian与PSI4算出的结果相差这么多；
2. 用def2-TZVPPD机组分解能量出现的错误该如何解决；
3. 新手提问“PSI4能考虑溶剂化效应吗”。
PS：附件为PSI4能量分解的输入和输出文件（所用机组为def2-TZVPPD）

作者
Author: sobereva 时间: 2022-5-28 21:39
1 不上传Gaussian输入输出文件没法说
但凡说文献，必须说清楚什么文献

2 aug-cc-pVDZ对Cr有定义，没必要考虑其它的

3 PSI4能不能考虑，和能量分解能不能考虑完全是两码事。前者能，后者不能直接考虑。

那叫基组不叫机组

作者
Author: LIEL 时间: 2022-5-29 10:25

sobereva 发表于 2022-5-28 21:39
1 不上传Gaussian输入输出文件没法说
但凡说文献，必须说清楚什么文献

哇！谢谢sob老师的回答，第一次发帖，有些地方的确没注意到，以后一定注意。
已上传聚吡咯和铬酸氢根的GAUSSIAN输入文件和输出文件，并且按照您推荐的aug-cc-pVDZ基组重新能量分解了一下，仍然是87kcal/mol左右，与Gaussian算的相差很大，附件也上传了sapt输入和输出文件。
PS：文献也上传了，文献1是用Q-chem能量分解的，总结合能很高（与我的体系相比吡咯环数量和阴离子种类不同）；文献2是用VASP算的，结合能也很高，但是我按照文献2体系用gaussian粗略的算过也只有-13kcal/mol左右。

作者
Author: sobereva 时间: 2022-5-30 04:16

LIEL 发表于 2022-5-29 10:25
哇！谢谢sob老师的回答，第一次发帖，有些地方的确没注意到，以后一定注意。
已上传聚吡咯和铬酸氢根的G ...

Gaussian计算带了隐式溶剂模型，SAPT计算不带，当然相差巨大，尤其是体系带净电荷的情况

作者
Author: LIEL 时间: 2022-5-30 09:01

sobereva 发表于 2022-5-30 04:16
Gaussian计算带了隐式溶剂模型，SAPT计算不带，当然相差巨大，尤其是体系带净电荷的情况

好的，谢谢sob老师！我也觉得可能是由于溶剂的影响，就是没想到会差这么多，我再用Gaussian试试真空条件的能量。

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