Molecular dynamics simulations of uranyl and uranyl carbonate adsorption at aluminosilicate surfaces

Adsorption at mineral surfaces is a critical factor controlling the mobility of uranium(VI) in aqueous environments. Therefore, molecular dynamics (MD) simulations were performed to investigate uranyl(VI) adsorption onto two neutral aluminosilicate surfaces, namely, the orthoclase (001) surface and the octahedral aluminum sheet of the kaolinite (001) surface. Although uranyl preferentially adsorbs as a bidentate inner-sphere complex on both surfaces, the free energy of adsorption on the orthoclase surface (-15 kcal mol(-1)) is significantly more favorable than that at the kaolinite surface (-3 kcal mol(-1)), which is attributed to differences in surface functional groups and the ability of the orthoclase surface to release a surface potassium ion upon uranyl adsorption. The structures of the adsorbed complexes compare favorably with X-ray absorption spectroscopy results. Simulations of the adsorption of uranyl complexes with up to three carbonate ligands revealed that uranyl complexes coordinated to up to two carbonate ions are stable on the orthoclase surface whereas uranyl carbonate surface complexes are unfavored at the kaolinite surface. Combining the MD-derived equilibrium adsorption constants for orthoclase with aqueous equilibrium constants for uranyl carbonate species indicates the presence of adsorbed uranium complexes with one or two carbonates under alkaline conditions, in support of current uranium(VI) surface complexation models.