The structures and vibrational spectra of small methanol clusters from dimer to decamer have been calculated using a newly developed intermolecular potential which is essentially based on monomer wave functions. Special care has been taken for the description of the electrostatic interaction using a distributed multipole representation and including a penetration term. In addition, the potential model consists of repulsion, dispersion, and induction terms. Based on this potential model cluster structures have been calculated. The lowest energy dimer configuration is linear, while from trimer to decamer for the most stable structures ring configurations were found. Tetramer, hexamer, and octamer have S4-, S6-, and S(8-)symmetry, respectively. Vibrational spectra of the CO stretch and the OH stretch mode have been determined in the harmonic and in the anharmonic approximation using perturbation theory and variational calculations. Up to the tetramer the experimental spectra of the CO stretch mode are well reproduced, for larger clusters an increasing blueshift with respect to the experimental evidence is found. The experimental data for the OH stretch mode of the dimer are fairly well reproduced in all approximations, however, the spectrum of the trimer can only be reproduced using the variational calculation which includes Darling-Dennison resonance terms.