计算化学公社

标题: 求助：二维Cu-MOF单点能计算SCF不收敛 [打印本页]

作者
Author: John_Tao 时间: 2024-3-22 23:38
标题: 求助：二维Cu-MOF单点能计算SCF不收敛
本帖最后由 John_Tao 于 2024-3-22 23:38 编辑

各位老师好，

我在用CP2K计算CuHHTP这个二维MOF时，使用PBE DZVP-MOLOPT-SR-GTH DFT-D3(BJ)进行单点能计算，但是SCF不收敛，请问老师们如何解决？

结构是文献中报道的结构（J. Am. Chem. Soc. 2020, 142, 12367-12373）。

附上了输入输出文件。

谢谢！

作者
Author: sobereva 时间: 2024-3-23 10:43
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

算REPEAT电荷的目的留出的真空区必须更大

作者
Author: John_Tao 时间: 2024-3-23 14:20
好的，谢谢sob老师回复！

作者
Author: John_Tao 时间: 2024-3-27 00:26

sobereva 发表于 2024-3-23 10:43
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.ht ...

sob老师好，我想追问一下：1.我想计算这个二维MOF的DOS，需要Z轴扩出真空层吗？还是应该直接用文献中实验测到的晶胞参数（层间距为3.3埃）进行计算呢； 2. 我是否需要XY方向扩建超胞呢~

作者
Author: sobereva 时间: 2024-3-27 00:40

John_Tao 发表于 2024-3-27 00:26
sob老师好，我想追问一下：1.我想计算这个二维MOF的DOS，需要Z轴扩出真空层吗？还是应该直接用文献中实验 ...

如果你想要计算层状晶体状态的MOF，而不是固-气界面的状态，就没有扩真空区的事
如果你要基于gamma点计算得到DOS，应当XY方向扩胞。CP2K结合Multiwfn产生DOS图时可以考虑k点，此时用原胞结合k点就完了，计算效率明显更高也更推荐。

作者
Author: John_Tao 时间: 2024-3-27 00:44

sobereva 发表于 2024-3-27 00:40
如果你想要计算层状晶体状态的MOF，而不是固-气界面的状态，就没有扩真空区的事

好的，感谢老师。那计算电荷也不需要扩真空层嘛~

作者
Author: sobereva 时间: 2024-3-27 01:27

John_Tao 发表于 2024-3-27 00:44
好的，感谢老师。那计算电荷也不需要扩真空层嘛~

取决于你算什么电荷
如果不是拟合静电势电荷，并且你要得到体相状态的原子电荷，没真空层的事
如果你要得到REPEAT等拟合静电势电荷，用于体相状态诸如吸附小分子的模拟，也没真空层的事
如果你要得到REPEAT等拟合静电势电荷，用于COF表面吸附的模拟，需要真空层

作者
Author: John_Tao 时间: 2024-3-27 01:51

sobereva 发表于 2024-3-27 01:27
取决于你算什么电荷
如果不是拟合静电势电荷，并且你要得到体相状态的原子电荷，没真空层的事
如果你要 ...

get！谢谢老师

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