Curve crossing and negative refraction in simulations of near-field coupled metallic nanoparticle arrays

We extend our previous results [R. Baer, J. Chem. Phys. 126, 014705 (2007).] to develop a simple theory of localized surface plasmon-polariton (LSPP) dispersion on regular arrays of metal nanoparticles in the weak-field and weak-damping limits. This theory describes the energy-momentum as well as the polarization-momentum properties of LSPP waves, both of which are crucial to plasmonic device design. We then explicitly compute the dispersion relation for isotropic and anisotropic two-dimensional square lattices, and show curve crossings between all three levels as well as negative refraction where the phase and group velocities (refractive indices), or at least their projection along the main axis, have different signs. The curve crossing implies that scattering between the different polarizations, and therefore different velocities, is easy at the curve crossing momenta, so that a quick change in wave packet direction can be achieved. Time-resolved wave packet dynamics simulations demonstrate negative refraction and the easy scattering over nanometer length scales. This paper also gives some computational schemes for future applications, such as a way to include source terms and how to efficiently treat dissipative effects.