计算化学公社

标题: cp2k计算得出的.molden文件,使用Multiwfn进行CDA计算，出现CDA和ECDA符号正负不一样 [打印本页]

作者
Author: kaikai 时间: 2024-11-5 21:46
标题: cp2k计算得出的.molden文件,使用Multiwfn进行CDA计算，出现CDA和ECDA符号正负不一样
本帖最后由 kaikai 于 2024-11-5 21:46 编辑

各位老师好，最近在学习使用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

作者
Author: kaikai 时间: 2024-11-6 19:08
相同的计算参数，计算了其他的几个模型，有符号相同的，也有符号不同的

作者
Author: sobereva 时间: 2024-11-7 06:59
1 MOLOPT基组不用担心弥散的事

2 可以

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

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

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

作者
Author: kaikai 时间: 2024-11-7 08:46

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

2 可以

谢谢老师的解答

作者
Author: yubeimu 时间: 2025-11-28 15:44

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

2 可以

老师您好，查到这个帖子。发现有相似情况，想向您请教一下。
1.我的是周期性体系，H2O在板层上的吸附结果。对吸附后构型做CDA分析，发现吸附后断键形成了H和OH，这种情况下需要分成3个片段进行考虑吗？
2.同时计算了H2O\CO2竞争吸附，竞争吸附后是将CO2\H2O当作一个片段还是分开呢？
3.如果H2O仍有断键现象，形成了游离的OH-和被吸附的H+;而CO2形成了新键，这种情况应该如何考虑片段更为合适呢？

作者
Author: sobereva 时间: 2025-11-29 05:35

yubeimu 发表于 2025-11-28 15:44
老师您好，查到这个帖子。发现有相似情况，想向您请教一下。
1.我的是周期性体系，H2O在板层上的吸附结 ...

取决于你实际想分析谁与谁之间的电荷转移，每一个“谁”都要作为一个片段

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