Simulation of oxygen vacancies at the Si–SiO2 interface