Convergence of Monte Carlo algorithm for solving integral equations in light scattering simulations

The light scattering process can be modeled mathematically using the Fredholm integral equation. This equation is usually solved after its discretization and transformation into the system of algebraic equations. Volume integral equations can be also solved without discretization using the Monte Carlo algorithm, but its application to the light scattering simulations has not been sufficiently studied. Here we present the implementation of this algorithm for one and three-dimensional light scattering computations and discuss its applicability in this field. We show that the Monte Carlo algorithm can provide valid and accurate results but, due to its convergence properties, it might be difficult to apply for problems with large volumes or refractive indices of scattering objects.