We compute electroluminescent signal in a current carrying single molecule junction using a superoperator formalism. Liouville space loop diagrams are used to identify all density matrix pathways that emit photons via the electroluminescence process. A frequency resolved spectrum is expressed in terms of the various Fock space states of the isolated molecule that participate in the creation and subsequent recombination of exciton. Application is made to a multilevel Coulomb blockade model system and to a gold鈥揵enzene-1,4-dithiol鈥揼old molecular junction.