Dmol3计算CoN4负载于石墨烯的形成能时候的问题

丁真xxx

在计算，CoN4负载在石墨烯的形成能的时候，套用文献给的公式，跟自己计算的模型图片如下面两图：




在计算中，使用DMol3模块，分别计算1.总体能量    
2.C32石墨烯的能量（超胞为4*4）
3.N2的化学式 （直接在doml3导入两个NN，几何优化后的能量除以2就是每个N的能量，我这引入了四个N原子，所以能量*4）
4.引入Co块体，几何优化后的总能除以Co原子的个数，求出Co的化学式
5.引起单层石墨烯，几何优化除以C的个数，得出C的化学式（本模型剔除了6个，C的化学式*6）
这样求得模型的形成能为0.034Ha左右，跟实际形成能差了大约0.2Ha（最稳定的结构的创新能应该的-0.2Ha）请问这样可能哪里出了问题？所有计算都用的FIne，K点取得441.
困扰一天了！！！求各位大神解答

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “ 第一性原理，Dmol3，形成能计算” 改了，以后务必注意，下次将删帖+扣分处理。