Classical molecular dynamics simulations were performed for the smectite clay hectorite with charge-balancing Cs+ cations using a newly developed structural model with a disordered distribution of Li/Mg substitutions in the octahedral sheet and the fully flexible CLAYFF force field. Calculations for systems with interlayer galleries containing 0–19 H2O/Cs+ suggest that the monolayer hydrate is the only stable state at all relative humidities at ambient pressure and temperature, in agreement with experimental results and previous molecular calculations. The basal spacing of this structure is also in good agreement with experimental values. In contrast to previous molecular modeling results, however, the new simulations show that interlayer Cs+ occurs on 2 different inner sphere adsorption sites: above the center of ditrigonal cavities and above Si tetrahedra. Unlike previous simulations, which employed a rigid clay model and fixed orientations of the structural −OH groups, the present results are obtained ...