The intra- and intermolecular bonding in the known phase of dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate, TKX-50, has been analyzed on the basis of the experimentally determined charge density distribution from high-resolution Xray diffraction data obtained at 20 K. This was compared to the charge density obtained from DFT calculations with periodic boundary conditions using both direct calculations and derived structure factors. Results of topological analysis of the electron density corroborate that TKX-50 is best described as a layered structure linked primarily by a number of hydrogen bonds as well as by a variety of other interactions. Additional bonding interactions were identified, including a pair of equivalent 1,5-type intramolecular closed-shell interactions in the dianion. Refinement of anharmonic motion was shown to be essential for obtaining an adequate model, despite the low temperature of the study. Although generally unusual, the implementation of anharmonic refinement provided a significant improvement compared to harmonic refinement of both traditional and split-core multipole models.