Recent years have seen significant interest in the mechanical properties of clay–polymer hybrids due to their suitability for possible application as sustainable materials in green chemistry. The objective of the present study was to investigate the mechanical properties of clay–polymer hybrids and their corresponding pristine smectite clay minerals. The density functional theory (DFT) method, employing the D3 scheme for corrections of dispersion interactions, was used to calculate elastic constants (Cij) of models of pristine smectites, particularly montmorillonite, beidellite, saponite, and hectorite, and their hybrids built on the polymer poly(2-methyl-2-oxazoline), PMeOx. Following that, the elastic moduli, encompassing the bulk modulus (KVRH), shear modulus (GVRH), Young’s modulus (EVRH), and Poisson’s ratio (ν), were calculated. The results revealed a reduction in elastic constants and elastic moduli following the intercalation of smectite clay minerals with the PMeOx polymer. The findings highlighted a distinctive ranking of mechanical properties among pristine smectite clay minerals and clay–polymer hybrids, with hectorite and its hybrid (Htr-PMeOx) demonstrating better performance compared with saponite, montmorillonite, and beidellite and their respective hybrids.

A detailed structural characterization of organo-clays is a key in understanding their properties. In this work, mono-, di-, tri-, and tetra-butylammonium (nBA; n = 1–4) cations intercalated in the layered clay mineral montmorillonite (Mnt) have been studied for the first time by combining a theoretical approach based on density functional theory (DFT) and infrared spectroscopy. The DFT calculations revealed the detailed structure and position of nBA cations in the interlayer space. A relation between the basal spacing (d001 parameter) and the cation size and structure was found, and explained with respect to the structure, composition, and size of the organic cations. Hydrogen bonds between -NHx/-CH3/-CH2 groups of the nBA cations and oxygen atoms of the basal planes of the Mnt layers were found to be an important factor for the arrangement and energetic stabilization of cations in the interlayer space. The N-H-O hydrogen bonds are stronger than C-H-O hydrogen bonds and the stabilization decreases with decreased number of bands. Analysis of DFT-calculated vibrational modes helped in understanding a problematic region of the experimental infrared spectra (4000–3000 cm-1), in which assignment of all vibrational modes unambiguously was not possible because of a significant overlap of broad bands.

A hybrid of montmorillonite (Mnt) and rifampicin (RIF) was synthesized and the structure and stability of the drug carrier system clarified. Density functional theory calculations involving dispersion corrections (DFT-D3) were performed to characterize interactions acting in the interlayer space of montmorillonite intercalated with rifampicin. The structure and stability of the RIF-Mnt intercalated complex were determined. Calculations revealed the deformation of the molecular structure of rifampicin after intercalation into the Mnt interlayer space due to the clay environment. The ansa chain of RIF was bent in the interlayer space compared with the structure of the RIF molecule in the monocrystal. RIF was keyed into the Mnt surface by means of numerous hydrogen bonds of weak to moderate strength. The calculated vibrational spectrum from ab initio molecular dynamics (AIMD) was in good agreement with the FTIR measured spectra and helped to analyze the overlapped vibrational bands. Based on analysis of structural stability, theoretical calculations revealed that Mnt is a suitable drug carrier for delayed release of the RIF drug. Batch adsorption experiments showed the large adsorption capacity of montmorillonite for RIF.

An ab initio theoretical approach has been used to calculate optimized geometries and the relative energies of various compositional arrangements in structures of dioctahedral smectites based on models consisting of two unit-cells. These calculations indicate that the energy differences between structures having vacancies in sites with cis- or trans-OH coordination are small and that their relative energies vary with the chemical nature of the substitutions. For example, a cis-OH coordination for the vacancy was the most stable when the interlayer charge originated from substitution of Al for Si in the tetrahedral sheet, whereas the trans-coordination was the more stable for most cases of substitution in the octahedral sheet, an exception being Fe(II) for Al where the cis-OH coordination was favored. It seems likely, therefore, that long-range structural disorder will be a common phenomenon in natural phyllosilicate specimens.

Theoretical models of the mechanical properties of hydrated smectites, saturated with a variety of cations, are of much value in determining the potential for their use in various applications, including clay-polymer nanocomposites, but the development of such models is still in its infancy. The purpose of this study was to calculate the effects of divalent cations on the structural and mechanical elasticity of montmorillonite under different degrees of hydration. A theoretical study of the swelling and hydration behavior of montmorillonite was, therefore, undertaken using density functional theory (DFT) to investigate the basal spacing behavior of the homoionic montmorillonite with varying amounts of water in the interlayer space. The effect of the degree of the hydration of divalent interlayer cations (Mg2+/Ca2+/ Sr2+/Ba2+) on the structure expansion of the interlayer space was analyzed. In addition, the results obtained were compared to calculations performed on the montmorillonite model with a monovalent cation (Na+). The basal spacing (d001) is governed by the size and the degree of hydration of the countercations. The structures containing divalent cations are more compact than structures with monovalent cations. Ba-exchanged montmorillonite was found to have the largest d001 value for any degree of hydration (‘dry,’ one water layer, or two layers). The basal spacings of ‘dry’ montmorillonite exchanged with small cations, Mg2+ and Ca2+, are very similar. In hydrated models, the d001 expansion correlates with the ionic radius of the interlayer cation. The dependence of the total electronic energy on the volume expansion was calculated. From the energetic curves, bulk modulus (B0) was obtained by fitting in order to show how the compliance of the material depends on the type of interlayer cation and on the degree of hydration. With increasing water content in the interlayer space, the bulk modulus decreased, suggesting that the c-axial compression becomes easier with increasing hydration of the clay mineral. The values of the bulk modulus in hydrated systems are less sensitive to the type of the interlayer cation.

Vibrational spectra of two kaolinite-dimethylsulfoxide intercalates, obtained using inelastic neutron scattering (INS), were analyzed with a view to understanding the dynamics of the hydrogen atoms in the structure. The main focus was on the spectral region 0–1700 cm−1, which is difficult to analyze using optical spectroscopy. The experimental vibrational spectra of kaolinite:dimethylsulfoxide and kaolinite:d6-dimethylsulfoxide collected using two different spectrometers were interpreted by means of the solid-state DFT calculations. Calculated spectra were obtained by both normal-mode analysis and molecular dynamics going beyond the harmonic approximation. The Al-O-H bending modes were found to be spread over the large interval 100−1200 cm−1, with the dominant contributions located between 800 and 1200 cm−1. The shape of the individual hydrogen spectrum depends on whether or not the respective hydrogen atom is involved in an O-H⋯O hydrogen bond and on its strength. The modes corresponding to the in-plane movements of the inner-surface hydrogen atoms are well defined and always appear at the top of the intervals of energy transfer. In contrast, the modes generated by the out-of-plane movements of the hydrogen atoms are spread over large energy intervals extending down to the region of external (lattice) modes. The C-H modes are concentrated mainly in the three regions 1200–1450 cm−1, 800–1100 cm−1, and 0–400 cm−1. While the first two regions are typical of the various deformational modes of methyl groups, the low-energy region is populated by the modes corresponding to the movements of the whole dimethylsulfoxide molecule.

Kaolinite is often used as a base for the synthesis of new organo-mineral nanomaterials designed for applications in industry and in environmental protection. To make the mineral structure more likely to interact with organic molecules, a kaolinite-methanol complex (KM) can be used. In the present study, different experimental procedures were tested to investigate the formation of the KM. The kaolinitedimethyl sulfoxide intercalation compound (KDS), either wet or dried, was used as a pre-intercalate. The samples obtained were characterized using X-ray diffraction, Fourier-transform infrared spectroscopy, CHNS elemental analysis, 13C CP-magic angle spinning nuclear magnetic resonance (MAS NMR), and 27Al and 29Si MAS NMR techniques. The method of density functional theory with dispersion corrections (DFT-D2) was used to explain the structure and to simulate the vibrational spectra of KM. Theoretical results were compared with experimental data. The most effective formation of the KM (d001 = 11.1 Å — wet; d001 = 8.7 Å — dried) was observed when the dried KDS precursor was used. In such conditions the degree of intercalation reached ~98% after 24 h of reaction time. As indicated by the CHNS elemental analysis, ~1/6 of the inner-surface OH groups were grafted by OCH3 groups. The esterification reaction was less efficient at higher temperatures or when wet KDS was used. In the latter case, the excess of very polar dimethyl sulfoxide molecules prevented intercalation of methanol and further grafting. Detailed analysis of the results of theoretical simulations revealed that the reaction of the KDS with methanol led to the formation of kaolinite with both grafted methoxy groups and intercalated methanol, and water molecules in the interlayer space. The spectra calculated revealed the contribution of individual vibrational modes into the complex bands, i.e. the energy of C-H vibrations was in the order: νasCHmet > νasCHmtx > νsCHmet > νsCHmtx.

The stability of organoclays prepared from smectites and organic cations depends on the type of used cation, among other factors. This study provides a prediction of the structure, stability and dynamic properties of organoclays based on montmorillonite (Mt) intercalated with two types of organic cations – tetrabutylammonium (TBA) and tetrabutylphosphonium (TBP) – using first-principle density functional theory. The results obtained from simulations were also used in the interpretation of the experimental infrared spectrum of the TBP-Mt organoclay. Analysis of interatomic distances showed that weak C–O···H hydrogen bonds were important in the stabilization of both TBA- and TBP-Mt models, with slightly stronger hydrogen bonds for the TBP cation. Calculated intercalation and adsorption reaction energies (ΔEint//ΔEads*/ΔEads**) confirmed that TBP-Mt structures (–72.4//–32.8/–53.8 kJ/mol) were considerably more stable than TBA-Mt structures (–56.7//–22.6/–37.4 kJ/mol). The stronger interactions of the alkyl chains of the TBP cation with Mt basal surfaces in comparison to those of the TBA cation were also correlated with the positions of the calculated bands of the C–H stretching vibrations.