We propagate the quantum mechanical density matrix, including dephasing, spontaneous emission and dissipation to study systems composed of two, three or four quantum dots in proximity to a plasmonic system when exposed to an ultrafast optical pulse [1]. The plasmonic system could be a single metal nanoparticle or an array of metal nanoparticles or some other plasmonic structure and serves to enhance the local electric fields that arise due to the pulse. We have previously studied a single quantum dot interacting with a plasmonic system and demonstrated how Fano resonances in tra