计算化学公社

标题: 求助CP2K对钙钛矿与水的固液界面模型做AIMD的设置 [打印本页]

作者
Author:
Weldingspock    时间: 2023-6-23 10:57
标题: 求助CP2K对钙钛矿与水的固液界面模型做AIMD的设置
各位老师好,本人想用CP2K进行AIMD的模拟,体系如下,钙钛矿与水的固液界面模型。slab下面三层原子固定,上面的弛豫。现在有几个问题想问一下各位:
1. 之前用VASP做AIMD,因为VASP默认是XYZ都加了周期性,那么我现在用CP2K做AIMD使用可以只用XY方向上的周期性吗,来防止上方的水分子在Z方向周期性存在的情况下,对slab下半部分固定原子造成影响?
2. 只在XY方向上加周期性的话,我这个体系还需要加偶极矫正吗。


希望大家解答,非常感谢!

作者
Author:
sobereva    时间: 2023-6-23 15:22
1 可以XY方向周期性结合PSOLVER MT,这样可以避免上层的水与周期镜像的slab底部的相互作用。但此时需要Z方向的盒子尺寸是现在的两倍以上,耗时也增加很多。

2 不
作者
Author:
Weldingspock    时间: 2023-6-23 15:45
sobereva 发表于 2023-6-23 15:22
1 可以XY方向周期性结合PSOLVER MT,这样可以避免上层的水与周期镜像的slab底部的相互作用。但此时需要Z方 ...

谢谢老师,我看了一下MT下不加偶极矫正可能会极大增加我这个体系的模拟时间,而且增加真空层的厚度做AIMD模拟的时候,水分子应该会全部跑到上方的真空层中。所以想问一下还有什么方法能够避免上方的水分子影响slab底部呢:
(1)现有的体系固定底部的slab原子,考虑XYZ周期性,但是不太清楚不固定的水分子会对固定的原子造成什么影响
(2) 把slab放到盒子中间,上下都加水
哪种更好一些呢
作者
Author:
sobereva    时间: 2023-6-23 15:54
Weldingspock 发表于 2023-6-23 15:45
谢谢老师,我看了一下MT下不加偶极矫正可能会极大增加我这个体系的模拟时间,而且增加真空层的厚度做AIMD ...

不用管什么影响slab底部,slab底部又不是关心的
如果说担心slab底部影响上方的水,把水层弄厚些

把slab挪到中间,只冻结slab内部部分的原子也行,此时水层需要更厚。这样相当于研究两个相同的界面,可以减小属性的统计误差,也可以视为在同样模拟时间内能增加一倍采样

作者
Author:
Weldingspock    时间: 2023-6-23 16:15
sobereva 发表于 2023-6-23 15:54
不用管什么影响slab底部,slab底部又不是关心的
如果说担心slab底部影响上方的水,把水层弄厚些

谢谢sob老师的解答




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3