The water-hexamethylenetetramine system displays, features of significant interest in the context of phase equilibria in molecular materials. First, it is possible to crystallize two solid phases depending oil temperature, both hexahydrate and anhydrous forms. Second, saturated aqueous solutions in equilibrium with these forms exhibit a negative dependence of solubility (retrograde) oil temperature. In this contribution, neutron scattering experiments (with isotopic substitution) of concentrated aqueous hexamethylenetetramine solutions combined with empirical potential structure refinement (EPSR) were used to investigate the time-averaged atomistic details of this system. Through the derivation of radial distribution functions. quantitative details emerge of the solution coordination, its relationship to the nature of the solid phases, and of the underlying cause of the solubility behavior of this molecule.