As above, let formula_25 denote a "translation operator" that shifts every wave function by the amount , where "ni" are integers. Because the crystal has translational symmetry, this operator commutes with the Hamiltonian operator. Moreover, every such translation operator commutes with every other. Therefore, there is a simultaneous eigenbasis of the Hamiltonian operator and every possible formula_25 operator. This basis is what we are looking for. The wave functions in this basis are energy eigenstates (because they are eigenstates of the Hamiltonian), and they are also Bloch states (because they are eigenstates of the translation operators; see Lemma above).