计算化学公社

标题: 关于slab模型的k点 [打印本页]
作者
Author: lawsai47    时间: 2021-11-24 15:43
标题: 关于slab模型的k点
原胞几何优化使用的是552精度，slab模型的1/埃与原胞几何优化一致，但是突然想到垂直方向有那么大的真空层，这样计算真空层所在州的k点采样准确性岂不是和其他两个轴差距很大
作者
Author: 卡开发发    时间: 2021-11-24 21:06
按道理真空层那个方向k点需求小，甚至直接gamma就够了，也就是说当某个真空层特别大或者k点特别大能量等性质几乎能收敛到一个稳定值，此时k多一点少一点的影响几乎很小。
作者
Author: lawsai47    时间: 2021-11-28 20:33
卡开发发 发表于 2021-11-24 21:06
按道理真空层那个方向k点需求小，甚至直接gamma就够了，也就是说当某个真空层特别大或者k点特别大能量等性 ...

谢谢老师，还有一个问题，就是关于表面计算使用的cutoff energy的大小，我原胞优化的时候用的是600ev，切出来的表面弱真空层选为0是原胞体积的三倍，但计算出的能量远高于原胞的3倍，是否应该加大cutoff energy呢，但是感觉计算量会特别大
作者
Author: 卡开发发    时间: 2021-11-28 23:36
lawsai47 发表于 2021-11-28 20:33
谢谢老师，还有一个问题，就是关于表面计算使用的cutoff energy的大小，我原胞优化的时候用的是600ev，切 ...
切出来的表面弱真空层选为0是原胞体积的三倍

这话我没看懂，要不发图片看看吧。一般体系可能用不上那么大的动能截断。
作者
Author: lawsai47    时间: 2021-11-29 07:02
卡开发发 发表于 2021-11-28 23:36
这话我没看懂，要不发图片看看吧。一般体系可能用不上那么大的动能截断。

这个是原胞的优化（一共12个原子），然后在计算110表面时，切出来的slab模型是36个原子，如果不加真空层体系的面积正好是原来的三倍，但是能量远高于原胞能量的三倍（50eV左右，平均每个原子的能量高了2eV）
作者
Author: 卡开发发    时间: 2021-11-29 10:13
lawsai47 发表于 2021-11-29 07:02
这个是原胞的优化（一共12个原子），然后在计算110表面时，切出来的slab模型是36个原子，如果不加真空层 ...

这个不应该，也许有东西弄错了。要不把.castep文件拿上来瞧瞧。
作者
Author: lawsai47    时间: 2021-11-29 10:37
卡开发发 发表于 2021-11-29 10:13
这个不应该，也许有东西弄错了。要不把.castep文件拿上来瞧瞧。

谢谢老师，结合优化后切出表面做的能量计算，结构是一致的
作者
Author: 卡开发发    时间: 2021-11-29 15:34
本帖最后由 卡开发发 于 2021-11-29 15:36 编辑
lawsai47 发表于 2021-11-29 10:37
谢谢老师，结合优化后切出表面做的能量计算，结构是一致的
标记为原胞：
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   12
                      Total number of species in cell =    2
                        Max number of any one species =    8
标记为110面：
                           -------------------------------
                                     Cell Contents
                           -------------------------------
                         Total number of ions in cell =   72
                      Total number of species in cell =    2
                        Max number of any one species =   48

不计入色散修正的情况：标记为原胞体系12个原子，能量是-14972.25478644eV，标记为110面的体系是72原子，能量是-89741.87201524eV，两个能量比值大概在5.994左右。
计入色散修正的情况：12原子体系能量为-14974.46001410eV，72原子体系为-89755.24858565eV，大概也在5.994被左右。
这个结果显然正常，但问题是72原子可能用的是fixed occ，所以没有展宽修正，而12原子的体系用了展宽，两个计算参数不完全一致。




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