Fast and scalable algorithm for the simulation of multiple Mie scattering in optical systems.

The Monte Carlo simulation of light propagation in optical systems requires the processing of a large number of photons to achieve a satisfactory statistical accuracy. Based on classical Mie scattering, we experimentally show that the independence of photons propagating through a turbid medium imposes a postulate for a concurrent and scalable programming paradigm of general purpose graphics processing units. This ensures that, without rewriting code, increasingly complex optical systems can be simulated if more processors are available in the future.