is the identity tensor. The sum in (1) runs over all particles and their periodic images and is evaluated using an efficient Ewald summation method. The linear system (1) is solved iteratively using the generalized minimal residual method (GMRES) to a relative error tolerance of 10 -3 . The ensemble-averaged extinction cross-section per nanoparticle volume is