The structures and infrared spectra of protonated ammonia clusters NH(4+)(NH3)n, for n < or = 8, are investigated using density functional-theory (DFT) calculations and semiempirical DFT/molecular dynamics simulations. For n < 5 the clusters are found to be mostly stable up to 100 K, while the larger clusters (n > or = 5) isomerize. Temperature effects are taken into account by performing ab initio molecular dynamics simulations with the computationally tractable self-consistent charges density functional tight-binding method. The infrared spectra at 10 K for the most stable isomers for n = 3-8 compare qualitatively with predissociation experiments, and using a common scaling factor almost quantitative agreement is found. For n > or = 6 the notion of multiple isomers present under the experimental conditions is supported. Of the 13 stable structures for n = 8 only three are found to survive at 100 K. All other clusters isomerize. Cluster structures are inferred from the analysis of the cumulative radial distribution function of the ammonia molecules surrounding the NH(4+) core. The infrared spectra are found to be typical for the structure of the clusters, which should help to relate the experimentally measured infrared spectra to the number and identity of the contributing isomers. For clusters that reorganize to a more stable isomer during the dynamics, the infrared spectrum is generally similar to that of the stable isomer itself. The clusters are found to preferably form globular structures, although chain-like arrangements are also among the low-energy configurations.