The absolute solvation energies (free energies and enthalpies) of the proton in ammonia are used to compute the pKa of species embedded in ammonia. They are also used to compute the solvation energies of other ions in ammonia. Despite their importance, it is not possible to determine experimentally the solvation energies of the proton in a given solvent. We propose in this work a direct approach to compute the solvation energies of the proton in ammonia from large‐sized neutral and protonated ammonia clusters. To undertake this investigation, we performed a geometry optimization of neutral and protonated ammonia 30‐mer, 40‐mer, and 50 mer to locate stable structures. These structures have been fully optimized at both APFD/6‐31++g(d,p) and M06‐2X/6‐31++g(d,p) levels of theory. An infrared spectroscopic study of these structures has been provided to assess the reliability of our investigation. Using these structures, we have computed the absolute solvation free energy and the absolute solvation enthalpy of the proton in ammonia. It comes out that the absolute solvation free energy of the proton in ammonia is calculated to be −1192 kJ mol–1, whereas the absolute solvation enthalpy is evaluated to be −1214 kJ mol–1. © 2019 Wiley Periodicals, Inc.