Surface calculations are performed using supercell geometries with periodic boundary conditions, ensuring a vacuum gap of at least 20 Å to prevent artificial interactions between neighboring images. In this study, we focus on the (110), which has been known as the most stable facet among all low-index surfaces of rutile TiO2 .