计算化学公社

标题: 计算V6O7(OCH3)12外重组能出现巨大的偏差 [打印本页]

作者
Author: 果果123 时间: 2023-10-11 15:59
标题: 计算V6O7(OCH3)12外重组能出现巨大的偏差
按照帖子http://sobereva.com/333  计算了V6O7(OCH3)12 netural电荷和自旋多重度设置为-1  6 charged的电荷和自旋多重度设置为0  5
输入文件和结果文件的链接如下    链接：https://pan.baidu.com/s/1HNZ61g74ZXINVjCwR6oAhQ?pwd=o8qr  提取码：o8qr
第一步计算之后能量是 E= -2334.9649505 Hartree 第二步计算的能量是 E = -2492.2838537 Hartree
也分别计算了第三步和第四步能量均在 -2334 Hartree 附近  结构和chk文件均使用的上一步计算过后的想请各位老师帮忙看一下哪里有问题需要改正
（因为对于多金属氧簇 CPCM溶剂化模型描述的最好所以不得不使用CPCM溶剂化模型）

作者
Author: sobereva 时间: 2023-10-11 16:10
我从没听说过“多金属氧簇 CPCM溶剂化模型描述的最好”
至少从原理上我不理解怎么会比默认的更严格的IEFPCM更好

先甭管非平衡溶剂问题，就直接算个溶剂下的电离能，看结果是多少。确保计算级别一致、收敛到了稳定波函数

作者
Author: wzkchem5 时间: 2023-10-11 16:32
对于结构优化、频率计算，IEFPCM、CPCM往往比SMD好，但即使这种情况下IEFPCM也比CPCM好；对于单点计算，SMD比IEFPCM、CPCM好。所以仅当程序不支持IEFPCM且你做的是结构优化、频率计算时，才可以说CPCM最好。建议再三确认你这个结论是从哪里看来的，是不是你断章取义了或者原作者过度引申了
不过即便CPCM不是最好也不足以解释这么大的能量差别。至于这个能量差别问题，建议检查一下（1）是否有基组线性相关问题（检查是否加了弥散函数，以及检查重叠矩阵的最小本征值是否太小）；（2）是否收敛到了正确的波函数（检查自旋布居和自旋密度是否符合化学常识）；（3）溶剂模型是否有数值不稳定性（尝试增大溶剂cavity的格点数、改变cavity的半径等等）

作者
Author: 风起~ 时间: 2023-10-12 09:28
应该是CPCM更加好收敛，
对于多核簇，IEFPCM乱七八糟的事情太多了

作者
Author: 果果123 时间: 2023-10-15 11:45

sobereva 发表于 2023-10-11 16:10
我从没听说过“多金属氧簇 CPCM溶剂化模型描述的最好”
至少从原理上我不理解怎么会比默认的更严格的IEFP ...

好的感谢sob老师回复我没有具体去考察只是看到了计算同类型体系均采用了这个CPCM的溶剂化模型后续再考究一下这个结论之前有计算过这个过程的内重组能在正常值得范围内后续采用了PCM溶剂化模型去计算发现也出现了巨大的数值误差

作者
Author: sobereva 时间: 2023-10-15 16:02

风起~ 发表于 2023-10-12 09:28
应该是CPCM更加好收敛，
对于多核簇，IEFPCM乱七八糟的事情太多了

迄今为止我没发现比CPCM更严格的IEFPCM有什么明显问题，或者是遇到什么问题时光靠改用CPCM就能解决的。如果有确实能证明CPCM比IEFPCM算此类问题更好的例子体系，欢迎给出，供大家交流讨论

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