Structural and thermodynamics properties of organo-modified montmorillonite clay

•We used molecular dynamics simulation based on two force fields (CLAYFF and CHARMM).•We have evaluated the thermodynamics and structural properties of organoclay.•The basal spacing and density are in good agreement with available experimental data.•The effect of the cation exchange capacity (CEC) on the basal spacing and the thermodynamics properties was considered.•The calculated isothermal bulk modulus agrees well with the density value of organoclays with two different CEC.

Polymer clay nanocomposites (PCNs) have been seen as the most novel materials in engineering applications since they exhibit significant improvement in mechanical and physical properties. Indeed, with few amount of organoclay, PCNs exhibit enhanced mechanical, optical, thermal and liquid or gas barrier properties compared to pure polymers and to their counterpart microcomposites. Thus, organoclays are extensively used as precursors in the preparation of PCNs. They are the best candidate in reinforcing PCNs because of the lightweight and the high availability of clay minerals in the nature. However, structure and physical phenomena arising at molecular level in organoclays, and subsequently in PCNs, are not completely or difficultly accessible with existing experimental techniques. In this work, molecular dynamics (MD) simulation was conducted using the combination of two force fields (CLAYFF and CHARMM) to evaluate the thermodynamics and structural properties of organoclay such as heat capacities, isothermal bulk modulus, density, basal spacing and chains arrangement in the interlayer spacing. Our results regarding the basal spacing and density are in fairly good agreement with available experimental data. This allows us to validate the use of the two force fields to represent interactions in organoclays. The effect of the cation exchange capacity (CEC) on the basal spacing and the thermodynamics properties is assessed. We found, through our MD simulation, that the calculated isothermal bulk modulus is in good agreement with the density value of organoclays with two different CEC.