“…We use the same angular discretization for the wave vectors of the electrons and optical phonons, which is not absolutely necessary, but greatly simplifies the implementation. Following this scheme, the modulus of optical phonon wave vectors q is discretized by q In this discretized coordinate system, the electron-phonon collision operator (4) takes the form (17) with the discretized electron distribution functions f i klm = f i (x k ,k(ε l ,ϕ m ),t) and f j kl±t η n = f i (x k ,k(ε l ±hω η ,ϕ n ),t) and the scattering element…”