1 Department of Chemistry and Biochemistry and Centre for Research in Molecular Modeling (CERMM), Concordia University, 7141 Sherbrooke Street West, Montréal, Québec H4B 1R6, Canada.

A polarizable model for ammonia is optimized based on the ab initio properties of the NH3 molecule and the NH3-NH3 and NH3-H2O dimers calculated at the MP2 level. For larger (NH3)m, NH3(H2O)n, and H2O(NH3)n clusters (m = 2-7 and n = 1-4), the model yields structural and binding energies in good agreement with ab initio calculations without further adjustments. It also reproduces the structure, density, heat of vaporization, self-diffusion coefficient, heat capacity, and isothermal...