The intermolecular potential energy surface of the electronic ground state of the ammonium-argon ionic dimer, NH[sub 4] [sup +] -Ar, is calculated by ab initiomethods using different levels of theory (MP2, MP4) and basis sets (aug-cc-pVXZ, X = D/T/Q). The deformation of the ammonium ion in the complex is shown to be small and its geometry is therefore fixed in these calculations to the tetrahedral structure optimized for the bare ion. The global minimum of the potential corresponds to a proton-bound structure with C[sub 3] symmetry (Rapproximately 3:4 A, Dapproximately 950cm[sup -1] ) and the barrier to internal rotation between the four equivalent minima is around 200cm[sup -1] . The three-dimensional potential is expanded in tetrahedral harmonics whose radially dependent coefficients, V(R), are compared for the considered levels of theory. The rotation-intermolecular vibration Hamiltonian is solved using a two-dimensional fixed-Rrepresentation of the calculated potentials, Vequals V(R), where the effective intermolecular separation, R[sup eff] , is determined from the experimental rotational constants of the complex. The accuracy of these parametrized potential energy surfaces is judged by their ability to reproduce the hindered rotor subband structure in the experimental upsilon[sub 3] (t[sub 2] ) infrared band of the complex. The simulations using the potentials calculated at the MP2/aug-cc-pVTZ or higher levels of theory reproduce the coarse structure of the experimental spectrum well. Further improvement could be achieved by least-squares fitting the potential parameters to the observed subband positions. The fitted V[sub 3] and V[sub 4] parameters remain in close agreement with those determined from the ab initiocalculations but the anisotropy of the potential is significantly different from that in a previous least-squares fit of V[sub 3] alone. [ABSTRACT FROM AUTHOR]