where J is the electrical current operator corresponding to Eq. (1), obtained through replacing c αk → c αk+eA and taking J = ∂H/∂A. To lowest order in λ or in 1/N ,. Within the memory matrix method, the dynamics is projected onto a set of "slow operators" that are nearly conserved by the Hamiltonian. We briefly review the method in Appendix A 1. In our model, the operators n αk = c † αk c αk are nearly-conserved in the limit of either small λ or large N [58]. The optical conductivity can then be cast in the form:where χ Jx,αk ≡´β 0 dτ J x (τ )n αk (0) and χ αk,βk ≡