Organic dyes such as methylene blue (MB) are often used in the characterization of clays and related minerals, but details of the adsorption mechanisms of such dyes are only partially understood from spectroscopic data, which indicate the presence of monomers, dimers, and higher aggregates for varying mineral surfaces. A combination of quantum (density functional theory) and classical molecular simulation methods was used to provide molecular detail of such adsorption processes, specifically the adsorption of MB onto kaolinite basal surfaces. Slab models with vacuum-terminated surfaces were used to obtain detailed structural properties and binding energies at both levels of theory, while classical molecular dynamics simulations of aqueous pores were used to characterize MB adsorption at infinite dilution and at higher concentration in which MB dimers and one-dimensional chains formed. Results for the neutral MB molecules are compared with those for the corresponding cation. Simulations of the aqueous pore indicate preferred adsorption on the hydrophobic siloxane surface, while charge-balancing chloride ions adsorb at the aluminol surface. At infinite dilution and in the gas-phase models, MB adsorbs with its primary molecular plane parallel to the siloxane surface to enhance hydrophobic interactions. Sandwiched dimers and chains are oriented perpendicular to the surface to facilitate the strong hydrophobic intermolecular interactions. Compared with quantum results, the hybrid force field predicts a weaker MB adsorption energy but a stronger dimerization energy. The structure and energetics of adsorbed MB at infinite dilution are consistent with the gas-phase binding results, which indicate that monomer adsorption is driven by strong interfacial forces rather than by the hydration properties of the dye. These results inform spectroscopic studies of MB adsorption on mineral surfaces while also revealing critical areas for development of improved hybrid force fields.

Zwitterionic porphyrin (tetrakis {4-(2-carboxyethyl)pyridinio}porphyrin (TPyCP)) with four pyridinium and carboxyl groups in the molecule was synthesized. The adsorption behaviour of TPyCP on the montmorillonite clay surface was examined in an aqueous colloidal solution. While the adsorption maximum λmax of TPyCP was 424 nm in water, it shifted to longer wavelengths on complex formation with clay. The λmax on the exfoliated clay surface was 450 nm, and that in the stacked clay sheets was 471 nm, respectively. Under alkaline conditions ([NaOH] = 2 × 10–4 M), the stacking behaviour of the clay was completely suppressed. It emerged that anionic parts in the porphyrin molecule can suppress the stacking of clay sheets. Judging from the quantitative analysis, the maximum adsorption amount of TPyCP was 100% vs. cation exchange capacity (CEC) of the clay. As the adsorption density of TPyCP increased, λmax shifted slightly to longer wavelengths due to the interactions between adjacent porphyrins. When the loading level of TPyCP was 200% vs. CEC, the stacking of TPyCP was indicated. The average stacking layer number was calculated to be 1.48. On the other hand, it is known that tetra cationic porphyrins without anionic parts do not form such stacking structures. Thus, it seems that zwitterionic porphyrin has the potential to form a three dimensional structure through electrostatic interactions between porphyrins on the clay surface.