Kinetic Monte Carlo (KMC) is a modern technique for modelling charge transport processes in materials. Researchers are highly interested in investigating the properties of materials for solar cells, batteries and LEDs. Furthermore, researchers wish to investigate the behaviour of these materials over relatively large timescales. KMC is a stochastic approach that models charge movement as a series of discrete jumps that occur at a particular rate. To accurately compute the probability that a jump occurs the electrostatic system energy must be computed. However, the computation of the electrostatic energy has been a prohibitively expensive operation for all but small system sizes. In this talk we describe the Fast Multipole Method (FMM) which is an existing method to compute electrostatic interactions with optimal computational complexity. Furthermore we describe our extension to the FMM that allows KMC iterations to be performed with optimal computational complexity with controllable error. We demonstrate that with our approach device sized systems may be modelled with accurate electrostatic interactions at a reasonable computational cost.