We construct model wavefunctions for the collective modes of fractional quantum Hall systems. The wavefunctions are expressed in terms of symmetric polynomials characterized by a root partition and a squeezed basis, and show excellent agreement with exact diagonalization results for finite systems. In the long wavelength limit, the model wavefunctions reduce to those predicted by the single-mode approximation, and remain accurate at energies above the continuum of roton pairs.