The equations that describe the time-evolution of transverse and longitudinal 15N magnetisations in tetrahedral ammonium ions, 15NH 4 +, are derived from the Bloch-Wangsness-Redfield density operator relaxation theory. It is assumed that the relaxation of the spin-states is dominated by (1) the intra-molecular 15N– 1H and 1H– 1H dipole–dipole interactions and (2) interactions of the ammonium protons with remote spins, which also include the contribution to the relaxations that arise from the exchange of the ammonium protons with the bulk solvent. The dipole–dipole cross-correlated relaxation mechanisms between each of the 15N– 1H and 1H– 1H interactions are explicitly taken into account in the derivations. An application to 15N-ammonium bound to a 41 kDa domain of the protein DnaK is presented, where a comparison between experiments and simulations show that the ammonium ion rotates rapidly within its binding site with a local correlation time shorter than approximately 1 ns. The theoretical framework provided here forms the basis for further investigations of dynamics of AX 4 spin systems, with ammonium ions in solution and bound to proteins of particular interest.