To study the change of interlayer structure of a Wyoming-type Na-montmorillonite as a result of the replacement of interlayer Na+ ions by cetylpyridinium (CP+) ions, a series of NPT Monte Carlo simulations of the clay mineral with different contents of CP+, Na+, Cl- ions and water in its interlayer space is carried out. In agreement with conclusions from experimental studies, the simulations show that the CP+ ions form monomolecular, bimolecular, and pseudotrimolecular layers with increasing interlayer contents. Calculated potential energies reveal that clay-organic interactions are stronger than organic-organic interactions in CP+-modified montmorillonite, which is in conformity with observations of earlier thermogravimetric experiments. The simulation results indicate that the pseudotrimolecular arrangement of CP+ ions is a prerequisite for the experimentally observed interlayer sorption of inorganic anions. Furthermore, in the interlayer space with a pseudotrimolecular layer, chloride ions favor the formation of pairs with inorganic rather than organic cations. On the basis of these findings and available experimental data, we propose that the interlayer sorption of inorganic anions from the pore space of an organically modified montmorillonite may occur not only in pairs with organic cations, as suggested earlier, but also in pairs with inorganic cations, which represents a so-far unconsidered and maybe more important mechanism of anion sorption on clay minerals.