cp2k计算得出的.molden文件,使用Multiwfn进行CDA计算，出现CDA和ECDA符号正负不一样

kaikai

各位老师好，最近在学习使用Multiwfn对cp2k计算得出的.molden文件进行CDA计算，其中模型的CDA计算结果与ECDA结果出现了符号相反的情况，CDA的计算结果与Mulliken计算结果相同，在论坛中寻找有关CDA的计算注意事项，对于出现的符号相反的情况，存在如下几个疑问：
1，在几篇帖子中提到，计算CDA时不要使用带有弥散函数的基组，在本次计算文件中，使用的DZVP-MOLOPT-SR-GTH系列是广义收缩的，是否不用考量弥散的问题？
2，结构优化使用的是默认的DZVP-MOLOPT-SR-GTH系列，进行能量计算时想着使用更高层次一点的TZVP-MOLOPT-SR-GTH系列,但是C,O,H等没有找到TZVP非SR版，使用TZVP-MOLOPT-GTH是否可以？
3，这次的计算结果文件中使用了混合基组了，在使用混合基组计算CDA时，应当注意什么问题？在cp2k输入文件中补充哪些信息？
4，目前出现模型的CDA计算结果与ECDA结果出现了符号相反的情况，建议如何改进？
麻烦老师们帮忙解答，谢谢

cp2k计算结果连接：https://www.alipan.com/s/SfuDZKcJf35

kaikai

相同的计算参数，计算了其他的几个模型，有符号相同的，也有符号不同的

sobereva

1 MOLOPT基组不用担心弥散的事

2 可以

3 保证不同体系里，相同原子用相同的基组就行了

4 CDA的目的不是算准净转移量，而主要是对转移从轨道角度对d、b进行分解讨论。和ECDA给出的转移量（本质上对应Mulliken方法算的转移量）符号不符也可能出现，不符的时候通常CDA的净转移量不合理，但不代表靠CDA讨论轨道对电子转移的贡献没实际意义。

当前的片段划分并不适合用CDA讨论什么，两个分离的H都不是个整体、不该当同一个片段。当成三个片段，其中两个片段都是H，或许结果还好点

kaikai

sobereva 发表于 2024-11-7 06:59
1 MOLOPT基组不用担心弥散的事

2 可以

谢谢老师的解答