Preferential interlayer adsorption sites in phyllosilicate clay edge models by molecular dynamics simulation

Phyllosilicate clay minerals have been proposed as a possible buffer material to be used in deep geological repositories containing high-level waste and used nuclear fuel. This work precisely characterizes ion interactions with two types of adsorption sites present in these clays: Mg _Al’ substitutions and undercoordinated edge surface atoms. A number of unique structural models were considered to represent the diverse local environments that ions in these systems are likely to encounter. Using molecular dynamics simulation with the CLAYFF potential, the spatial distribution, interlayer composition, and residence times of Na ⁺ and Cl ⁻ ions as radionuclide analogs in pyrophyllite and montmorillonite clay models were investigated to identify the most favorable conditions for sequestration. The most significant factor impacting ion adsorption was found to be the localization of charge density at substitution sites. In a montmorillonite system in which substitution sites were distributed evenly to produce a low charge density, sequestration performance was found to be comparable to pyrophyllite.