The dynamics of intermolecular vibrations in molecular aggregates can be directly monitored by ultrafast spectroscopic techniques. Starting with the third order response function obtained using the path-integral approach, we calculate the pump-probe spectrum of a molecular dimer, taking into account the effects of an intermolecular vibration coupled to the solvent. The various contributions to the signal which represent distinct Liouville space paths are identified and analyzed.