Simple extended constant capacitance surface complexation models have been developed to represent the adsorption of polyaromatic dyes (9-aminoacridine, 3,6-diaminoacridine, azure A and safranin O) to kaolinite, and the competitive adsorption of the dyes with Cd. The formulation of the models was based on data from recent publications, including quantitative adsorption measurements over a range of conditions (varying pH and concentration), acid-base titrations and attenuated total reflectance-Fourier transform infrared spectroscopic data. In the models the dye molecules adsorb as aggregates of three or four, forming outer-sphere complexes with sites on the silica face of kaolinite. Both electrostatic and hydrophobic interactions are implicated in the adsorption processes. Despite their simplicity, the models fit a wide range of experimental data, thereby supporting the underlying hypothesis that the flat, hydrophobic, but slightly charged silica faces of kaolinite facilitate the aggregation and adsorption of the flat, aromatic, cationic dye molecules.