|qk⟩ is the q th Bloch function at point k in reciprocal space, u qk (r) is the cell periodic part of that Bloch function, and C is the matrix of crystal orbital (CO) coefficients for that k point . To simplify the notation, we write the triple sum over cell indices as a single sum over the multi-index l