A four-site intermolecular potential function for NH3 has been developed and tested in Monte Carlo statistical mechanics simulations of the liquid at its boiling point (-33.35°C) and 1 atm. The potential yields good thermodynamic results for liquid ammonia, while the structures of the liquid are characterized through radial distribution functions and hydrogen-bonding analyses. The results indicate that each ammonia forms on average three hydrogen bonds and the liquid contains winding chains of hydrogen bonded monomers. Roughly linear hydrogen bonds predominates in the liquid. The results were also compared with those obtained using the five-site model developed by Ferrario et al. Although the computed heat of vaporization, liquid density, and the radial distribution functions are in good accord with our results and experimental data, to our surprise, a significant difference in the dimer and liquid structures exists between the two theoretical models. The structural difference is due to the occurrence of unrealistic dimer structures with donation of hydrogen bonds to the opposite side of the ammonia lone pair in the liquid from the model of Ferrario et al. In contrast, both ab initio calculations and the present four-site model predict no stable complexes for these structures. These findings emphasize the importance of specific consideration of bimolecular interactions in developing potential functions for fluid simulations and suggest that erroneous results might be obtained if the potential functions are fitted to reproduce the condensed-phase properties alone.