We calculated the extinction spectra and distributions of near-field intensity for our simulated nanoparticle networks using a recently developed many-bodied, mutual polarization method (MPM). The N induced particle dipoles p i in an electric field E 0 of magnitude E 0 = |E 0 | are computed by solving the linear system of equations for each frequency ω