The last decade has seen incredible growth in the quality of experiments being done on single molecule junctions. Contemporary experimental measurements have expanded far beyond simple electron transport. Measurement of vibronic effects, quantum interference and coherence effects, molecular optical response (Raman spectroscopy), and molecular spintronics are just some of the continuing areas of research in single molecule junctions. Experimental advancements demand advanced theoretical treatments, which can be used accurately within appropriate physical regimes, in order to understand measured phenomena and predict interesting directions for future study. In this dissertation we will study systems with strong intra-system interactions using a many-body states based approach. We will be focused on three related processes in molecular junctions: electron transport, electronic energy transfer, and molecular excitation. Inelastic electron transport in the regime of strong and nonlinear electron-vibration coupling within and outside of the Born-Oppenheimer regime will be investigated. To understand their appropriateness, we will compare simple semi-classical approximations in molecular redox junctions and electron-counting devices to fully quantum calculations based on many-body system states. The role of coherence and quantum interference in energy and electron transfer in molecular junctions is explored. Experiments that simultaneously measure surface enhanced Raman scattering and electron conduction have revealed a strong interaction between conducting electrons and molecular excitation. We investigate the role of the molecular response to a classical surface plasmon enhanced electric field considering the back action of the oscillating molecular dipole. Raman scattering is quantum mechanical by nature and involves strong interaction between surface plasmons in the contacts and the molecular excitation. We develop a scheme for treating strong plasmon-molecular excitation interactions quantum mechanically within nonequilibrium molecular junctions. Finally we perform preliminary calculations of the Raman spectrum of a three- ring oligophenylene vinylene terminating in amine functional groups molecule in a molecular junction and compare our results to experimental measurements. This work is the first steps towards full calculations of the optical response of current-carrying molecular junction, which should combine classical calculations of the plasmon enhanced electric field with quantum calculations for the plasmon-molecular exciton interaction and nonequilibrium Raman scattering