The covalent and ionic clusters of ammonium nitrate and hydroxyl ammonium nitrate are characterized using density functional theory and second-order vibrational perturbation theory. The most stable structures are covalent acid-base pairs for the monomers and ionic acid-base pairs for the dimers. The hydrogen-bonding distances are greater in the ionic dimers than in the covalent monomers, and the stretching frequencies are significantly different in the covalent and ionic clusters. The anharmonicity of the potential energy surfaces is found to influence the geometries, frequencies, and nuclear magnetic shielding constants for these systems. The inclusion of anharmonic effects significantly decreases many of the calculated vibrational frequencies in these clusters and improves the agreement of the calculated frequencies with the experimental data available for the isolated neutral species. The calculations of nuclear magnetic shielding constants for all nuclei in these clusters illustrate that quantitatively accurate predictions of nuclear magnetic shieldings for comparison to experimental data require the inclusion of anharmonic effects. These calculations of geometries, frequencies, and shielding constants provide insight into the significance of anharmonic effects in ionic materials and provide data that will be useful for the parametrization of molecular mechanical force fields for ionic liquids. Anharmonic effects will be particularly important for the study of proton transfer reactions in ionic materials.