An efficient method of treating solvent effects in excited state molecular dynamics (ESMD) is implemented and tested by exploring the solvatochromic effects in substituted p-phenylene vinylene oligomers. A continuum solvent model is used which has very little computational overhead. This allows simulations of ESMD with solvent effects on the scale of hundreds of picoseconds for systems of up to hundreds of atoms. At these time scales, solvatochromic shifts in fluoresence spectra can be described. Solvatochromic shifts in absorption and fluorescence spectra from ESMD are compared with time-dependent density functional theory calculations and experiments. (C) 2015 Elsevier B.V. All rights reserved.