CP2K可以在结构优化时计算REPEAT电荷吗

thor

之前看到社长博文里面有使用CP2K计算REPEAT电荷的例子，但我不知道使用Multiwfn生成MOF结构优化脚本的时候是否可以选择同时计算REPEAT和输出.cube文件用来计算DDEC电荷？还是说必须优化完成另外开一个单点能的计算再去计算电荷呢

John_Tao

我每次都是另开一个单点。我记得选上molden文件的话，每优化一步都会生成一个。不知道电荷和cube文件会不会。

sobereva

别跟结构优化搅合在一起，应当单独用个任务
只需要计算REPEAT电荷就够了，原理上比任何其它原子电荷都更适合做MOF体系的基于力场的模拟的目的

thor

sobereva 发表于 2024-5-14 17:32
别跟结构优化搅合在一起，应当单独用个任务
只需要计算REPEAT电荷就够了，原理上比任何其它原子电荷都更适 ...

好的谢谢社长。我还有两个问题想请教一下：
1. 我看到您给出的MOF-5计算REPEAT电荷是同类型原子电荷都是相等的，这是做了什么特殊处理吗，我计算自己MOF的REPEAT电荷同类原子都是有差异，我需要特殊操作或者做平均处理吗？
2. 我用CP2K优化一个立方晶胞的MOF，优化出来的三个方向的晶格常数不相等，有一些差别，我尝试过保持对称性的操作，还是会有点差异，不过一般都差的很小，这是正常的吗，还是有办法保持一致呢。我看到别人给出DFT优化结果晶格常数都是相等

sobereva

thor 发表于 2024-5-14 17:45
好的谢谢社长。我还有两个问题想请教一下：
1. 我看到您给出的MOF-5计算REPEAT电荷是同类型原子电荷都是 ...

1 空间上（而不是化学上）等价原子才可能REPEAT电荷相等。另外，即便空间上等价，若CUTOFF比较糙的话，有可能由于静电势格点精度不够，导致拟合出来的原子电荷不等价。你可以对等价原子取平均

2 如果你用KEEP_SYMMETRY T并且&FORCE_EVAL / &SUBSYS / &CELL / SYMMETRY设了CUBIC的话，不应该算完了边长不同

thor

sobereva 发表于 2024-5-14 23:13
1 空间上（而不是化学上）等价原子才可能REPEAT电荷相等。另外，即便空间上等价，若CUTOFF比较糙的话，有 ...

好的，非常感谢社长解惑