“…One way to sidestep this is to assume a virtual crystal (VC) model of the oxide, and to extend the TB Hamiltonian to include this region, however, to the best of our knowledge, VC models of SiO 2 are still not well established in literature. 50,51 Although a truncated TB Hamiltonian over the silicon region can be iterated with the Poisson equation over the whole domain, there could be issues relating to charge inconsistencies which would affect the convergence. We have therefore neglected an exact quantitative description of gate screening, and have used the field as the free parameter to investigate its effects approximately.…”