下午茶 发表于 2019-6-5 08:45
我是算N2O吸附在CeO2上的吸附能的，按照文献把CeO2切了110面后，2*3进行扩胞，The vacuum gap thickness ...

气相分子计算时，放在大盒子里，不是把吸附结构表面移除就OK的

hakuna 发表于 2019-6-4 16:08
首先，如果计算的是气相中的分子，你的slab需要改一下，a, b方向有点短，增加对应slab边长；
其次，最好把 ...

我是算N2O吸附在CeO2上的吸附能的，按照文献把CeO2切了110面后，2*3进行扩胞，The vacuum gap thickness设置的10埃，这样盒子的a,b方向大小就确定了，我是把N2O的盒子和吸附的盒子大小一样，我试过提高收敛精度EDIFF=1E-8，EDIFFG=-1E-4，优化后虚频更大了，还有请问为什么要固定一个原子呢？

首先，如果计算的是气相中的分子，你的slab需要改一下，a, b方向有点短，增加对应slab边长；
其次，最好把分子键轴沿着某一方向放置，固定其中的某一个原子；
最后，优化时力的收敛精度不够，POTIM取小点。
无关参数都去掉

你可以看一下你们方向相关的文献上面气体分子的频率是不是从热力学数据库中取的。