A Coarse-Grained Model for the Mechanical Behavior of Na-Montmorillonite Clay

Mesoscale Simulations Reveal How Salt Influences Clay Particles Agglomeration in Aqueous Dispersions

The aggregation of clay particles is an everyday phenomenon of scientific and industrial relevance. However, it is a complex multiscale process that depends delicately on the nature of the particle-particle and particle-solvent interactions. Toward understanding how to control such phenomena, a multiscale computational approach is developed, building from molecular simulations conducted at atomic resolution to calculate the potential of mean force (PMF) profiles in both pure and saline water environments. We document how it is possible to use such a model to develop a fundamental understanding concerning the mechanism of particle aggregation. For example, using molecular dynamics simulations conducted at the mesoscale in implicit solvents, it is possible to quantify the size and shape of clay aggregates as a function of system conditions. The approach is used to emphasize the role of salt concentration, which directly affects the potentials of the mean forces between kaolinite particles. While particle agglomeration in pure water yields large aggregates, the presence of sodium chloride in the aqueous brine leads instead to a large number of small aggregates. These results are consistent with macroscopic experimental observations, suggesting that the simulation protocol developed could be relevant for preventing pore blocking in heterogeneous porous matrixes.

Exploring the biointerfaces: ab initio investigation of nano-montmorillonite clay, and its interaction with unnatural amino acids

We investigated the interaction between biomimetic Fe and Mg co-doped montmorillonite nanoclay and eleven unnatural amino acids. Employing three different functionals (PBE-GGA, PBE-GGA + U, and HSE06), we examined the clay's structural, electronic, and magnetic properties. Our results revealed the necessity of using PBE-GGA + U with U ≥ 4 eV to accurately describe key clay properties. We identified amino acids that strongly interacted with the clay surface, with steric orientation playing a crucial role in facilitating binding. Our DFT calculations highlighted significant electrostatic interactions between the amino acids and the clay slab, with the amino group's predominant role in this interaction. These findings hold promise for designing amino acids for clay-amino acid systems, leading to innovative bio-material composites for various applications. Additionally, our ab-initio molecular dynamics simulations confirmed the stability of clay-amino acid systems under ambient conditions, and the introduction of an implicit water solvent enhanced the binding energy of amino acids on the clay surface.

Molecular insights into the temperature and pressure dependence of mechanical behavior and dynamics of Na-montmorillonite clay

Sodium montmorillonite (Na-MMT) clay mineral is a common type of swelling clay that has potential applications for nuclear waste storage at high temperatures and pressures. However, there is a limited understanding of the mechanical properties, local molecular stiffness, and dynamic heterogeneity of this material at elevated temperatures and pressures. To address this, we employ all-atomistic (AA) molecular dynamics (MD) simulation to investigate the tensile behavior of Na-MMT clay over a wide temperature range (500 K to 1700 K) and pressures (200 atm to 100 000 atm). The results show that increasing the temperature significantly reduces the tensile modulus, strength, and failure strain, while pressure has a minor effect compared to temperature, as seen in the normalized pressure-temperature plot. Mean-square displacement (MSD) analysis reveals increased molecular stiffness with increasing pressure and decreasing temperature, indicating suppressed atomic mobility. Our simulations indicate temperature-dependent dynamical heterogeneity in the Na-MMT model, supported by experimental studies and quantified local molecular stiffness distribution. These findings enhance our understanding of the tensile response and dynamical heterogeneity of Na-MMT clay under extreme conditions, aiding the development of clay minerals for engineering applications such as nuclear waste storage and shale gas extraction.

Molecular Dynamics Model to Explore the Initial Stages of Anion Exchange involving Layered Double Hydroxide Particles

A classical molecular dynamics (MD) model of fully unconstrained layered double hydroxide (LDH) particles in aqueous NaCl solution was developed to explore the initial stages of the anion exchange process, a key feature of LDHs for their application in different fields. In particular, this study focuses on the active corrosion protection mechanism, where LDHs are able to entrap aggressive species from the solution while releasing fewer corrosive species or even corrosion inhibitors. With this purpose in mind, it was explored the release kinetics of the delivery of nitrate and 2-mercaptobenzothiazole (MBT, a typical corrosion inhibitor) from layered double hydroxide particles triggered by the presence of aggressive chloride anions in solution. It was shown that the delamination of the cationic layers occurs during the anion exchange process, which is especially evident in the case of MBT^-.

First-Principles Study on Optoelectronic Properties of Fe-Doped Montmorillonite Clay

A theoretical investigation is conducted to describe optoelectronic properties of Fe-doped montmorillonite nanoclay under spin states of low spin (LS), intermediate spin (IS), and high spin (HS). Ground state electronic properties are studied using spin-polarized density functional theory calculations. The nonradiative and radiative relaxation channels of charge carriers are studied by computing nonadiabatic couplings (NACs) using an "on-the-fly" approach from adiabatic molecular dynamics trajectories. The NACs are further processed using a reduced density matrix approach with the Redfield formalism. The computational results are presented for electronic density of states, absorption spectra, charge carrier dynamics, and photoluminescence (PL) by comparing various spin multiplicities. Results on spin α and spin β components are independent and quite different because of the partial occupation of Fe 3d states. Overall, HS is the most stable with the largest Fe-O distances. One finds different nonradiative relaxation pathways in space and on the time scale for electrons and holes. The Redfield PL reveals obvious Fe 3d-3d transitions for LS and IS.