计算化学公社

标题: cp2k计算表面吸附离子能量可以吗？ [打印本页]

作者
Author: wwp 时间: 2022-3-8 11:41
标题: cp2k计算表面吸附离子能量可以吗？
各位老师好

我想请问一下，是否可以使用CP2K来优化计算周期性体系表面吸附离子的结构能量？只需要改变input文件中charge参数吗？

谢谢各位老师指教。

作者
Author: jiangning198511 时间: 2022-3-8 11:43
可以周期条件需要设置为XY,

作者
Author: wwp 时间: 2022-3-8 12:18

jiangning198511 发表于 2022-3-8 11:43
可以周期条件需要设置为XY,

好的谢谢老师

作者
Author: ChemG 时间: 2022-3-8 12:48
CP2K不知道是否有和VASP一样对于带电离子吸附计算不太适用的问题，之前大佬们有讨论过VASP计算表面吸附带电离子不太合适的问题 http://bbs.keinsci.com/thread-9422-1-1.html

作者
Author: wwp 时间: 2022-3-8 13:06

ChemG 发表于 2022-3-8 12:48
CP2K不知道是否有和VASP一样对于带电离子吸附计算不太适用的问题，之前大佬们有讨论过VASP计算表面吸附带电 ...

也看到了这个问题，所以考虑使用cp2k来算表面电荷的问题。真的是太麻烦了

作者
Author: jiangning198511 时间: 2022-3-8 14:42

wwp 发表于 2022-3-8 13:06
也看到了这个问题，所以考虑使用cp2k来算表面电荷的问题。真的是太麻烦了

你可以按照VASP手册中建议的那样画一下不同真空层的能量，如果是随真空层变化那就不能用，如果是恒定的那就可以使用，再者多调研一下相关文献看看使用什么方法来实现的

作者
Author: wwp 时间: 2022-3-9 11:00

jiangning198511 发表于 2022-3-8 14:42
你可以按照VASP手册中建议的那样画一下不同真空层的能量，如果是随真空层变化那就不能用，如果是恒定的 ...

老师，您好

您说的画一下不同真空层的能量，是指在改变真空层的厚度然后用vasp测试一下能量的变化吗？

作者
Author: jiangning198511 时间: 2022-3-9 11:40

wwp 发表于 2022-3-9 11:00
老师，您好

您说的画一下不同真空层的能量，是指在改变真空层的厚度然后用vasp测试一下能量的 ...

是用CP2K 并且周期性选XY, 看看是否有能量变化，如果基本不变就说明这个策略可行，当然你也可以平行对比一下VASP的结果，如果审稿人质疑方法可行性这个对比就比较说明问题另外可以找找相关的文献看看是采用什么方法

作者
Author: ChemG 时间: 2022-3-9 12:24

jiangning198511 发表于 2022-3-9 11:40
是用CP2K 并且周期性选XY, 看看是否有能量变化，如果基本不变就说明这个策略可行，当然你也可以平行对 ...

在之前的讨论中 http://bbs.keinsci.com/thread-21224-1-1.html ，卡开发发前辈提到过在相同的真空层设置下相对能量还可以比较，是不是说在表面带电的体系比较相对能量还是可以勉强可以接受的

作者
Author: wwp 时间: 2022-3-9 12:26

jiangning198511 发表于 2022-3-9 11:40
是用CP2K 并且周期性选XY, 看看是否有能量变化，如果基本不变就说明这个策略可行，当然你也可以平行对 ...

好的谢谢老师

作者
Author: wwp 时间: 2022-3-9 12:29

ChemG 发表于 2022-3-9 12:24
在之前的讨论中 http://bbs.keinsci.com/thread-21224-1-1.html ，卡开发发前辈提到过在相同的真空层设置 ...

应该是这个意思吧。带电体系真的是麻烦

作者
Author: jiangning198511 时间: 2022-3-10 10:12
本帖最后由 jiangning198511 于 2022-3-10 10:16 编辑

ChemG 发表于 2022-3-8 12:48
CP2K不知道是否有和VASP一样对于带电离子吸附计算不太适用的问题，之前大佬们有讨论过VASP计算表面吸附带电 ...

VASP中是XYZ周期，而CP2K中可以自行设置周期性边界条件，所以有可能解决这一问题，当然需要进一步文献调研

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