A general problem in your calculation is that the Si surface is not passivated, so there will be a lot of dangling bonds that give an incorrect electronic structure.
I would also use a Neumann condition on the vacuum side.
Can't say if this helps your specific problem, but they should be addressed, else results might not be correct.