Damped harmonic oscillators are utilized to calculate Franck-Condon factors within displaced harmonic oscillator approximation. This is practically done by scaling unperturbed Hessian matrix that represents local modes of force constants for molecule in gaseous phase, and then by diagonalizing perturbed Hessian matrix it results in direct modification of Huang-Rhys factors which represent normal modes of solute molecule perturbed by solvent environment. Scaling parameters are empirically introduced for simulating absorption and fluorescence spectra of an isolated solute molecule in solution. The present method is especially useful for simulating vibronic spectra of polycyclic aromatic hydrocarbon molecules in which hydrogen atom vibrations in solution can be scaled equally, namely the same scaling factor being applied to all hydrogen atoms in polycyclic aromatic hydrocarbons. The present method is demonstrated in simulating solvent enhanced X (1)Ag ↔ A(1)B1u absorption and fluorescence spectra of perylene (medium-sized polycyclic aromatic hydrocarbon) in benzene solution. It is found that one of six active normal modes v10 is actually responsible to the solvent enhancement of spectra observed in experiment. Simulations from all functionals (TD) B3LYP, (TD) B3LYP35, (TD) B3LYP50, and (TD) B3LYP100 draw the same conclusion. Hence, the present method is able to adequately reproduce experimental absorption and fluorescence spectra in both gas and solution phases.