The solvation free energy of proton in methanol was calculated by B3LYP flavor of density functional calculations in combination with the Poisson-Boltzmann continuum solvation model. In order to check the adequacy of the computation level, the free energies of clustering in the gas phase were compared with the experimental data. The solvents were taken into account in a hybrid manner, i.e. one to five molecules of methanol were explicitly considered while other solvent molecules were represented with an implicit solvation model.