Numerical quadrature in the Brillouin zone for periodic Schrödinger operators