We here reveal the specific chromophore-solvent interactions that contribute to both inhomogenous and vibronic broadening. Solvent in the axial region of the chromophore interacts with the π-system, leading to substantial inhomogeneous broadening. Our analysis of a single axial H-bonded methanol shows that it produces a strong electric-field at the amine site with the O-H bond dipole of methanol aligned with the π-system of the chromophore, and that this one solvent molecule is responsible for the large increase in the spectral width. The methanol molecules that form strong H-bonds peripheral to the chromophore decrease the excitation energy, but do not contribute to inhomogeneous broadening or any increase in the reorganization energy, presumably because they are very stable throughout the dynamics and lead to minimal energy gap fluctuation. However, these peripheral H-bonds likely play a key role in capturing the correct vibronic broadening, as we observe a substantial increase in vibronic spectral width if these methanol molecules are treated with QM during geometry optimization and normal mode computation. Both our potential energy scans along the H-bond and our geometry optimizations in frozen solvent suggest that the MM solvent interfacing with the QM chromophore holds the solvent too close to the chromophore compared to a QM solvent treatment, which likely restricts the motion of the chromophore atoms leading to a decrease in vibronic coupling. We find that only with the QM treatment of the solvent does the CAM-B3LYP/6-31G(d) spectral shape approach that of experiment.