Molecular Simulation of Cesium Adsorption at the Basal Surface of Phyllosilicate Minerals

A better understanding of the thermodynamics of radioactive cesium uptake at the surfaces of phyllosilicate minerals is needed to understand the mechanisms of selective adsorption and help guide the development of practical and inexpensive decontamination techniques. In this work, molecular dynamics simulations were carried out to determine the thermodynamics of Cs+ adsorption at the basal surface of six 2:1 phyllosilicate minerals, namely pyrophyllite, illite, muscovite, phlogopite, celadonite, and margarite. These minerals were selected to isolate the effects of the magnitude of the permanent layer charge (⩽2), its location (tetrahedral vs. octahedral sheet), and the octahedral sheet structure (dioctahedral vs. trioctahedral). Good agreement was obtained with the experiments in terms of the hydration free energy of Cs+ and the structure and thermodynamics of Cs+ adsorption at the muscovite basal surface, for which published data were available for comparison. With the exception of pyrophyllite, which did not exhibit an inner-sphere free energy minimum, all phyllosilicate minerals showed similar behavior with respect to Cs+ adsorption; notably, Cs+ adsorption was predominantly inner-sphere, whereas outer-sphere adsorption was very weak with the simulations predicting the formation of an extended outer-sphere complex. For a given location of the layer charge, the free energy of adsorption as an inner-sphere complex varied linearly with the magnitude of the layer charge. For a given layer charge location and magnitude, adsorption at phlogopite (trioctahedral sheet structure) was much less favorable than at muscovite (dioctahedral sheet structure) due to electrostatic repulsion between adsorbed Cs+ and the H atom of the OH- ion directly below the six-membered siloxane ring cavity. For a given layer charge magnitude and octahedral sheet structure, adsorption to celadonite (octahedral sheet layer charge) was favored over adsorption to muscovite (tetrahedral sheet layer charge) due to the increased distance to the surface K+ ions and the decreased distance to the O atom of the OH- ion directly below the surface cavity.