The effect of solvent in the first two (pi) *$IMP(pi) excited states of indole and 3-methylindole, (1La and 1Lb), has been studied using a hybrid theoretical method that couples molecular dynamics and a semiempirical molecular orbital procedure. This method yields information about the mechanism and the time scales involved in bulk solvent reorganization after excitation of the indole chromophore to the 1La, (or 1Lb), state, as well as the solvent induced inhomogeneous broadening. The fluorescence red shifts predicted in several solvents, (water, methanol, butanol, and dimethyl ether), agree reasonably well with experimental values. These time resolved calculations also predict inversion of the excited states with respect to vacuum and indicate that the solvent relaxation has two components: the first one is inertial in character, with a Gaussian shape and a half-width at half-maximum of approximately 15 fs for water, and of 100 to 300 fs for the other solvents. The second component shows an exponential decay behavior and seems related to the longitudinal relaxation time of the solvent. The correlation times for this component are approximately 170 fs for water and a few picoseconds for the other solvents. Calculations performed in butanol at 0 K indicate that we can expect a fast, (1 to 5 ps), red shift in the fluorescence, as large as 1000 to 2000 cm-1, even when the solvent is rigid, because the inertial response is still possible in this environment. Absorption shifts are not calculated well using this method because it does not take into account the electronic polarizability of the solvent, which would be responsible for half of the absorption shift.