关于slab模型的k点

lawsai47

原胞几何优化使用的是552精度，slab模型的1/埃与原胞几何优化一致，但是突然想到垂直方向有那么大的真空层，这样计算真空层所在州的k点采样准确性岂不是和其他两个轴差距很大

卡开发发

按道理真空层那个方向k点需求小，甚至直接gamma就够了，也就是说当某个真空层特别大或者k点特别大能量等性质几乎能收敛到一个稳定值，此时k多一点少一点的影响几乎很小。

lawsai47

卡开发发 发表于 2021-11-24 21:06
按道理真空层那个方向k点需求小，甚至直接gamma就够了，也就是说当某个真空层特别大或者k点特别大能量等性 ...

谢谢老师，还有一个问题，就是关于表面计算使用的cutoff energy的大小，我原胞优化的时候用的是600ev，切出来的表面弱真空层选为0是原胞体积的三倍，但计算出的能量远高于原胞的3倍，是否应该加大cutoff energy呢，但是感觉计算量会特别大

卡开发发

lawsai47 发表于 2021-11-28 20:33
谢谢老师，还有一个问题，就是关于表面计算使用的cutoff energy的大小，我原胞优化的时候用的是600ev，切 ...

切出来的表面弱真空层选为0是原胞体积的三倍

这话我没看懂，要不发图片看看吧。一般体系可能用不上那么大的动能截断。

lawsai47

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

这个是原胞的优化（一共12个原子），然后在计算110表面时，切出来的slab模型是36个原子，如果不加真空层体系的面积正好是原来的三倍，但是能量远高于原胞能量的三倍（50eV左右，平均每个原子的能量高了2eV）

卡开发发

lawsai47 发表于 2021-11-29 07:02
这个是原胞的优化（一共12个原子），然后在计算110表面时，切出来的slab模型是36个原子，如果不加真空层 ...

这个不应该，也许有东西弄错了。要不把.castep文件拿上来瞧瞧。

lawsai47

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

谢谢老师，结合优化后切出表面做的能量计算，结构是一致的

卡开发发

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原子的体系用了展宽，两个计算参数不完全一致。