求助：二维Cu-MOF单点能计算SCF不收敛

John_Tao · 发表于 Post on 2024-3-22 23:38:02

本帖最后由 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）。

附上了输入输出文件。

谢谢！

sobereva · 发表于 Post on 2024-3-23 10:43:53

CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

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

John_Tao · 发表于 Post on 2024-3-23 14:20:54

好的，谢谢sob老师回复！

John_Tao · 发表于 Post on 2024-3-27 00:26:43

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方向扩建超胞呢~

sobereva · 发表于 Post on 2024-3-27 00:40:48

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

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

John_Tao · 发表于 Post on 2024-3-27 00:44:28

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

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

sobereva · 发表于 Post on 2024-3-27 01:27:50

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

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

John_Tao · 发表于 Post on 2024-3-27 01:51:15

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

get！谢谢老师

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助：二维Cu-MOF单点能计算SCF不收敛