Variational Monte Carlo, diffusion Monte Carlo, and stereographic projection path integral simulations are performed on eight selected species from the (NH(3))(n), (ND(3))(n), (NH(2)D)(n), and (NH(3))(n-1)(ND(3)) clusters. Each monomer is treated as a rigid body with the rotation spaces mapped by the stereographic projection coordinates. We compare the energy obtained from path integral simulations at several low temperatures with those obtained by diffusion Monte Carlo, for two dimers, and we find that at 4 K, the fully deuterated dimer energy is in excellent agreement with the ground state energy of the same. The ground state wavefunction for the (NH(3))(2-5) clusters is predominantly localized in the global minimum of the potential energy. In all simulations of mixed isotopic substitutions, we find that the heavier isotope is almost exclusively the participant in the hydrogen bond.