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.

The intercalation of DNA into layered double hydroxides (LDHs) has various applications, including drug delivery for gene therapy and origins of life studies. The nanoscale dimensions of the interlayer region make the exact conformation of the intercalated DNA difficult to elucidate experimentally. We use molecular dynamics techniques, performed on high performance supercomputing grids, to carry out large-scale simulations of double stranded, linear and plasmid DNA up to 480 base pairs in length intercalated within a magnesium-aluminum LDH. Currently only limited experimental data have been reported for these systems. Our models are found to be in agreement with experimental observations, according to which hydration is a crucial factor in determining the structural stability of DNA. Phosphate backbone groups are found to align with aluminum lattice positions. At elevated temperatures and pressures, relevant to origins of life studies which maintain that the earliest life forms originated around deep ocean hydrothermal vents, the structural stability of LDH-intercalated DNA is substantially enhanced as compared to DNA in bulk water. We also discuss how the materials properties of the LDH are modified due to DNA intercalation.

A quantitative description is presented of the dynamical process of polymer intercalation into clay tactoids and the ensuing aggregation of polymer-entangled tactoids into larger structures, obtaining various characteristics of these nanocomposites, including clay-layer spacings, out-of-plane clay-sheet bending energies, X-ray diffractograms, and materials properties. This model of clay–polymer interactions is based on a three-level approach, which uses quantum mechanical and atomistic descriptions to derive a coarse-grained yet chemically specific representation that can resolve processes on hitherto inaccessible length and time scales. The approach is applied to study collections of clay mineral tactoids interacting with two synthetic polymers, poly(ethylene glycol) and poly(vinyl alcohol). The controlled behavior of layered materials in a polymer matrix is centrally important for many engineering and manufacturing applications. This approach opens up a route to computing the properties of complex soft materials based on knowledge of their chemical composition, molecular structure, and processing conditions.

In oil and gas drilling operations, drilling fluids perform essential tasks such as lubricating the drill bit, providing hydrostatic pressure and removing drill cuttings. One important function of the drilling fluid is to stop compacted clay minerals, commonly encountered in drilling operations, from taking up water from the drilling fluids and consequently swelling. Such a scenario can have an adverse impact on drilling operations and may lead to significantly increased oil well construction costs. With increasingly stringent environmental guidelines determining which swelling inhibitors are available for use in the oilfield as drilling fluid additives, there is a need to fully understand the mechanisms of clay hydration in order to design new swelling inhibitors which conform to evolving regulations. Using a range of computational techniques and analysis, combined with known experimental results, we have devised a set of ''rule-based'' design criteria for clay-swelling inhibitors. To achieve this, we have formulated a hydration energy parameter, which assesses the changes in energy during the step-wise progression from mono-to bi-to trilayers of water in the clay sheet galleries. This parameter can be used to rationalise and predict the swelling profiles for clays containing both cationic and neutral clay swelling inhibitors. The rules we have devised are as follows: (i) Cationic inhibitors should be able to replace sodium ions in the interlayer. (ii) Cationic inhibitors should possess a water soluble, hydrophobic backbone. (iii) Cationic inhibitors should have primary di-amine or mono-quaternary amine functionality. (iv) Cationic inhibitors should have little alcohol functionality. (v) The hydrophobic backbone of the cationic inhibitor should be long enough to form a dense monolayer in the interlayer. (vi) For neutral inhibitors, the inhibitor should be a water soluble organic molecule of low molecular weight with well defined domains of relatively high hydrophobicity and small domains of hydrophobicity. Our ''rule-based'' criteria will facilitate the rational design of improved-and more environmentally acceptable-clay swelling inhibitors for oilfield drilling operations.

Graphite's lubricating properties due to the "weak" interactions between individual layers have long been known. However, these interactions are not weak enough to allow graphite to readily exfoliate into graphene on a large scale. Separating graphite layers down to a single sheet is an intense area of research as scientists attempt to utilize graphene's superlative properties. The exfoliation and processing of layered materials is governed by the friction between layers. Friction on the macroscale can be intuitively understood, but there is little understanding of the mechanisms involved in nanolayered materials. Using molecular dynamics and a new forcefield, graphene's unusual behavior in a superlubric state is examined, and the energy dissipated between two such surfaces sliding past each other is shown. The dependence of friction on temperature and surface roughness is described, and agreement with experiment is reported. The accuracy of the simulated behavior enables the processes that drive exfoliation of graphite into individual graphene sheets to be described. Taking into account the friction between layers, a peeling mechanism of exfoliation is predicted to be of lower energy cost than shearing.