Large-Scale Molecular Dynamics Study of Montmorillonite Clay:  Emergence of Undulatory Fluctuations and Determination of Material Properties

Abstract

Large-scale atomistic simulations, which we define as containing more than 100 000 atoms, are becoming more commonplace as computational resources increase and efficient classical molecular dynamics algorithms are developed. With the advent of grid-computing methods, it is now possible to simulate even larger systems efficiently. Using this new technology, we have simulated montmorillonite clay systems containing up to approximately ten million atoms whose dimensions approach those of a realistic clay platelet. This considerably extends the spatial dimensions of microscopic simulation into a domain normally encountered in mesoscopic simulation. The simulations exhibit emergent behavior with increasing size, manifesting collective thermal motion of clay sheet atoms over lengths greater than 150 Å. This motion produces low-amplitude, long-wavelength undulations of the clay sheets, implicitly inhibited by the small system sizes normally encountered in atomistic simulation. The thermal bending fluctuations allow us to calculate material properties, which are hard to obtain experimentally due to the small size of clay platelets. Montmorillonite is commonly used as a filler in clay−polymer nanocomposites, and estimation of the elastic properties of the composite requires accurate knowledge of the elastic moduli of the components. We estimate the bending modulus to be 1.6 × 10-17 J, corresponding to an in-plane Young's modulus of 230 GPa. We encounter a clay sheet persistence length of approximately 1400 Å, which dampens the undulations at long wavelengths for the largest system in our study.

Linking preferred orientations to elastic anisotropy in Muderong Shale, Australia

“…For other phyllosilicates, e.g., kaolinite and illite-smectite, there are only model elastic properties available. It is likely that elastic properties vary with details of chemical composition, ion exchange capacity, polytype, and water content (e.g., Suter et al, 2007;Mazo et al, 2008). The volume, shape, and distribution of organic content (kerogen) are other factors that can affect the elastic properties (e.g., Vernik and Nur, 1992).…”

Use of UKLight as a fast network for data transport from Grid infrastructures

“…Iperf is a network tool which we used to test UDP and TCP bandwidth between networked computers 5 . Iperf UDP tests provided the maximum bandwidth of the connection before packet loss is seen.…”

“…Operating over at least tens of thousands of atoms, emergent mesoscopic properties are observed in full atomistic detail. Large-scale simulations of clay nanocomposites reveal long wavelength, low amplitude thermal undulations [5]. Small simulation sizes implicitly inhibit long wavelength clay sheet flexibility due to the periodic boundaries used in condensed matter molecular dynamics, which effectively pin the clay sheet at the edges of the simulation cell.…”

Characterizing and Modeling Mechanical Properties of Nanocomposites-Review and Evaluation

“…In 1997, Cornwell et al used molecular dynamics to predict the elastic properties of single-walled carbon nanotubes [101]. In recent years, molecular dynamics simulation has been extensively used in predicting mechanical properties of carbon nanotubes and nanotubes reinforced composites [102][103][104][105][106][107][108][109], graphite/epoxy nanocomposites [110][111][112], and other nanocomposites [113][114][115][116][117][118][119].…”