求助 使用CP2K计算α-MnO2表面吸附甲苯的吸附能异常偏大

11WHISPER28

各位老师好，
我想做α-MnO2表面吸附甲苯的吸附能计算，切的110面，分别几何优化了表面，吸附体系与吸附质，但是计算出的吸附能有十几ev,自己刚接触计算没有多久，其中几何优化与计算能量都没有警告与报错，以下是计算的输入文件与输出文件，对于计算这个体系应该注意哪些参数考虑，因为之前在几何优化的时候有发现scf不收敛，所以尝试了DFT+U，scf确实可以收敛了但是对于初始磁矩猜测把握不是很准确，对于不太能把握的有什么比较好的办法去得到他的初始磁矩吗？希望各位老师不吝赐教

wzkchem5

检查吸附前后的磁矩是否定性一致
常见的材料查文献即可知道是铁磁耦合还是反铁磁耦合，以及每个金属离子的磁矩是多少，据此设定初始磁矩即可

11WHISPER28

好的谢谢老师，我这个计算是按照反铁磁耦合计算的但是计算的吸附体系的能量要比slab表面加上吸附质的能量要负的更多导致吸附能有-14ev，我有将几何优化后的吸附体系基底与吸附质单独计算能量但是他们加起来的能量也与吸附体系的能量绝对差值是15ev左右，三个计算能量的文件设置都是一样的，不太清楚问题出在哪里

喝酸奶的巧乐兹

110面的不同暴露面可以对比一下，可能有的暴露面O多，有的暴露面Mn多
设置磁矩初猜时，暴露面Mn的两种不同自旋方向也可以对比一下