Monte Carlo simulations of the extinction rate of dense media with randomly distributed dielectric spheres based on solution of Maxwell’s equations

Abstract
We present Monte Carlo simulations of the extinction rate of dense media with randomly distributed dielectric spheres that occupy up to 25% by volume and size parameter ka = 0.2. Maxwell’s equations in multiple-scattering form are solved iteratively for each realization. Convergence is demonstrated numerically by varying the number of iterations, the number of spheres up to 4000, and the number of realizations. Results are compared with that of the independent-scattering approximation, Foldy’s approximation, the quasi-crystalline approximation, and the quasi-crystalline approximation with coherent potential. The simulations are in good agreement with the last two approximations.