To ensure adequate representation of the system, the super cell is set to be larger than 10 Å in each direction. Only the gamma-point is considered for all calculations. The energy of electronic iterations is set to reach 10 -5 eV. Additionally, a energy cutoff of for the plane-wave-basis 550 eV. Spin polarization is not considered.