Monte Carlo (MC) simulations are frequently used to describe random processes. An important application is multiple Mie scattering in perturbed media. To significantly reduce execution time and increase the statistical accuracy of the MC simulations we have implemented a concurrent algorithm running on graphics processing units. We comment on execution time and scalability.