Improved scheme for modeling a spaser made of identical gain elements

Simulation of the most general configuration of a spaser where each gain element is assumed to have distinct behavior is almost intractable with current computing technologies. Hence, here we consider a representative, yet simple enough system with all possible interactions allowed, in which all gain elements have identical features including two energy levels, to gain deeper understanding of the operational characteristics. We reason in an intuitively sound, yet mathematically rigorous, way how to exploit the symmetry of the most efficient method currently available to reduce the size of the problem by order OM , where M is the size of the truncated plasmon number state. Finally, we explore some physical results as predicted by our method.