“…For lattice constants in the near field regime dominated by dipole-dipole interaction [3], the LSPR can be well described by the quasi-static approximation for nanosized islands on a substrate, where the near field interaction is included between the islands in a square or hexagonal lattice, also including the induced dipoles in the substrate. The Bedeaux-Vlieger model [15] is such a method, but it does not take into account polarization coupling, and hence result in blockdiagonal Mueller matrices. However, the plasmon modes occurring as a result of the lattice for larger lattice constants, can be expected to both convey polarization coupling and resonances not accounted for in the Bedeaux-Vlieger model, and also reflect the two-dimensional photonic crystal character of the surface.…”