We have performed molecular dynamics simulations to investigate the adsorption of radionuclide elements species onto substituted Montmorillonite (001) surface in the presence of different counterions. The structure and the dynamics of uranyl ion as well as its aquo, chloride ion, and carbonate complexes are analyzed. In addition, we have studied the surface energy between layered Montmorillonite sheets and the work of adhesion between radionuclide and charged Montmorillonite. The clay model used here is a Wyoming-type Montmorillonite with 0.75e negative charge per unit cell resulting from substitutions in Octahedral and Tetrahedral sheets. The system model was constructed based on CLAYFF force field potential model. To evaluate the thermodynamic work of adhesion, each surface and clay layer regions are converted to a thin film model. One and two species of radionuclide elements (UO2(H2O)(5)UO2CO3(H2O)(5), and UO2Cl2(H2O)(5)) were deposited near the clay surface in a pseudo-two-dimensional periodic cell. Analysis shows that the uranyl ion structure is preserved with two axial oxygen atoms detected at 1.8 angstrom. Radial distribution functions results indicate that average U-O-w distances are 2.45-2.61 angstrom, and 2.29-2.40 angstrom for U-O-c distance. Average U-Cl distances are 2.78-3.08 angstrom, which is relatively larger than that of Uranium atom-Oxygen atom because of electrostatic factors. (C) 2013 Elsevier B.V. All rights reserved.