如何计算ITO(110)中氧空位缺陷

BBH

各位老师好，我最近在计算ITO（110）面的氧空位，采用了GGA+U和DFT-D2，构建的模型（此时还未构建氧空位）中包含了48个O原子，26个In原子以及4个Sn原子，进行了结构优化，然后通过文献阅读，找到了In2O3里面氧空位最容易形成的位置去除掉一个氧原子然后结构优化，但是计算已经过了快四天了，还没有收敛，想问一下是有哪里出错吗，以及计算氧空位需要考虑些什么呀，非常感谢。