Ab initio calculations were carried out to study the geometry, solvation free energy, dipole moment, molecular electrostatic potential (MESP), Mulliken and Natural charge distribution, polarizability, hyperpolarizability, Natural Bond Orbital (NBO) energetic and different molecular properties like global reactivity descriptors (chemical hardness, softness, chemical potential, electronegativity, electrophilicity index) of 2-methylimidazole. B3LYP/6-31G(d,p) level of theory was used to optimize the structure both in the gas phase and in solution. The solvation free energy, dipole moment and molecular properties were calculated by applying the Solvation Model on Density (SMD) in four solvent systems, namely water, dimethylsulfoxide (DMSO), n-octanol and chloroform. The computed bond distances, bond angles and dihedral angles of 2-methylimidazole agreed reasonably well with the experimental data except for C(2)-N(1), C(4)-C(5) and N(1)-H(7) bond lengths and N(1)-C(5)-C(4) bond angle. The solvation free energy, dipole moment, polarizability, first order hyperpolarizability, chemical potential, electronegativity and electrophilicity index of 2-methylimidazole increased on going from non-polar to polar solvents. Chemical hardness also increased with increasing polarity of the solvent and the opposite relation was found in the case of softness. These results provide better understanding of the stability and reactivity of 2-methylimidazole in different solvent systems.