Urotropin (U) and azelaic acid (AA) form 1:1 co-crystals (UA) that give rise to a rather complex diffraction pattern, the main features of which are diffuse rods and bands in addition to the Bragg re¯ections. UA is characterized by solvent inclusions, parasite phases, and high vacancy and dislocation densities. These defects compounded with the pronounced tendency of U to escape from the crystal edi®ce lead to at least seven exotic phase transitions (many of which barely manifest themselves in a differential scanning calorimetry trace). These involve different incommensurate phases and a peritectoid reaction in the recrystallization regime (T h b 0X6). The system may be understood as an OD (order±disorder) structure based on a layer with layer group Pcc2 and cell a o 9 4.7, b 9 26.1 and c 9 14.4 A Ê . At 338 K the layer stacking is random, but with decreasing temperature the build-up of an orthorhombic MDO (maximal degree of order) structure with cell a 1 = 2a o , b 1 = b, c 1 = c and space group Pcc2 is begun (at $301 K). The superposition structure of the OD system at T = 286 (1) K with space group Bmmb and cell a = 2a o , b = b and c = ca2 owes its cohesion to van der Waals interactions between the AA chains and to three types of hydrogen bonds of varied strength between UÐU and UÐAA. Before reaching completion, this MDO structure is transformed, at 282 K, into a monoclinic one with cell a m = Àa o c/4, b m = b, c m = À2a o ca2, space group P2 1 ac, spontaneous deformation $2, and ferroelastic domains. This transformation is achieved in two steps: ®rst a furtive triggering transition, which is not yet fully understood, and second an improper ferroelastic transition. At $233 K, the system reaches its ground state (cell a M = a m , b M = b, c M = c m and space group P2 1 ac) via an irreversible transition. The phase transitions below 338 K are described by a model based on the interaction of two thermally activated slip systems. The OD structure is described in terms of a three-dimensional Monte Carlo model that involves ®rst-and second-neighbour interactions along the a axis and ®rst-neighbour interactions along the b and c axes. This model includes random shifts of the chains along their axes and satisfactorily accounts for most features that are seen in the observed diffraction pattern.