A proper modeling is essential for calculating the surface properties of complex oxides, for example, perovskite oxides. Many studies have shown that traditional selections of slab models, either symmetric or asymmetric ones, may vary the calculated results of the polar surfaces, which will finally affect the perception of surface properties. To solve this issue, we proposed a passivation method to handle the complex polar surfaces, where the pseudo-atoms with fractional charges were used to balance non-zero net charges on one side of the surfaces. Based on this method, we revisited the low-index surfaces of LaCoO3, including (11̅02), (1̅104), and (0001) surfaces, with different terminations. The results suggested that the (11̅02)-LaO, (11̅02)-CoO2, and (0001)-LaO3 terminations were stable surfaces in different O and Co chemical environments. The oxygen vacancy formation energies (EOv) of the (11̅02)-LaO and (11̅02)-CoO2 terminations were also evaluated. Different from the unpassivated cases, the EOv of the passivated LaO and CoO2 terminations converged rapidly when slab thickness increased. The upward shift trend of the EOv after passivation was more consistent with the experimental results. The electronic structure analysis indicated that the oxygen vacancy mainly changed the valence state of the surrounding Co atoms. This work provides a more accurate calculation method, which may benefit the study on catalytic properties of complex oxides.
Read full abstract