jpqiu 发表于 2024-5-14 17:25
老师你好,如果是加抗衡离子的话,我有三个疑问,
1.我做过相应的吸附能测试, 带电分子的吸附能随真空层的 ...

没具体结构不好说
先参考此文的讨论Phys.Chem.Chem.Phys., 20, 8456 (2018)

sobereva 发表于 2024-4-8 02:30
用三维周期边界条件算这个
如果不用抗衡离子，那就做个吸附能随真空层厚度的收敛性测试，确保真空层厚度 ...

老师你好,如果是加抗衡离子的话,我有三个疑问,
1.我做过相应的吸附能测试, 带电分子的吸附能随真空层的厚度是线性关系的,并不能收敛,不知道是哪里出了问题.
2.抗衡离子的位置其实是会影响吸附能的大小的,对比不同体系的吸附能时候,我们怎么去确保这个抗衡离子带来的影响是固定的呢,貌似解决不了带电分子和中性分子吸附能的对比?
3.抗衡离子在计算的时候它的电荷并不是Local的,所以分子和金属层的电荷分布是会受到它的影响的吧?

晓声寒月 发表于 2024-4-7 08:19
老师，是的，切完两个表面之后，目前严格来说就是 一维周期性了；那么如果只考虑垂直屏幕的周期性，确实 ...

用三维周期边界条件算这个
如果不用抗衡离子，那就做个吸附能随真空层厚度的收敛性测试，确保真空层厚度足够

sobereva 发表于 2024-4-7 06:08
看不出来你算的是什么东西，似乎只是在垂直于屏幕方向是周期性的。如果就是这样，这连二维周期性也不是，只 ...

老师，是的，切完两个表面之后，目前严格来说就是 一维周期性了；那么如果只考虑垂直屏幕的周期性，确实后面的问题就没必要了；
但就是想问问如果硬是不当成一维（抱歉不是想成心钻牛角尖，一方面是我们小方向下很多人就把它当三维周期性，另一方面是一维的Psolver似乎不是很好？），那么原则上，是不是应该对001方向（上下）的偶极和100方向（侧面）的负电荷加以处理？
最后就是想问问，如果硬是不处理侧面的负电荷，是不是对计算吸附能也无伤大雅（误差抵消？）

看不出来你算的是什么东西，似乎只是在垂直于屏幕方向是周期性的。如果就是这样，这连二维周期性也不是，只是一维周期性而已。没法回答你PSOLVER和偶极校正的问题。

CP2K没有SCPC校正，Makov-Payne校正自己手动计算就完了。第三届北京科音CP2K第一性原理计算培训班里讲“第一性原理计算的理论基础”的部分明确给了其具体公式和手动计算的例子，结果和文献里一致。本来CP2K对三维周期性带电体系自动就用了均匀补偿背景电荷的处理，Makov-Payne校正只不过是能让结果随盒子尺寸收敛更快而已。