Slab models of rutile TiO2 (110) surface: DFT and DFT+U calculations

The F2B slab model weakens the oscillations and calculations on F2B model quickly converge. However, the F2B model leads to artificial narrowness of band gap. Besides, when the number of layers increases, surface energy obtained from all three slab models approaches similar values. In particular, values of surface energy from DFT calculations converge to the experimental range for all three slab models.