Montmorillonite (Mnt) is the most common 2:1 clay mineral. The negative layer charge of Mnt originates from the isomorphic substitution, including octahedral substitution (Mg for Al) and tetrahedral substitution (Al for Si; Brigatti et al., Reference Brigatti, Galan, Theng, Bergaya, Theng and Lagaly2006). Exchangeable hydrated metal cations exist in the interlayer of Mnt. Organoclays can be formed through ion exchange between these metal cations and organic cations. Organoclays have been proven to be versatile materials in industry, with a wide range of applications as adsorbents, rheological control agents and components of oil well drilling fluids (Yariv & Cross, Reference Yariv and Cross2002; Lagaly et al., Reference Lagaly, Ogawa, Dékány, Bergaya and Lagaly2006; Park et al., Reference Park, Ayoko and Frost2011; Bhatt et al., Reference Bhatt, Somani, Mody and Bajaj2013; Silva et al., Reference Silva, Sousa, Menezes, Neves, Santana and Ferreira2014; Shah et al., Reference Shah, Pan, Shukla, Shah and Chiang2018).
In the synthesis of organoclays, alkylammonium cations are favoured surface modifiers due to their affordability and ease of preparation (Gieseking, Reference Gieseking1939; Nigam et al., Reference Nigam, Setua, Mathur and Kar2004). Many experiments have been conducted to investigate the properties of alkylammonium-intercalated Mnt (Barrer & MacLeod, Reference Barrer and MacLeod1995; Bala et al., Reference Bala, Samantaraya and Srivastava2000; Scholtzová et al., Reference Scholtzová, Madejová, Jankovič and Tunega2016; Veiskarami et al., Reference Veiskarami, Sarvi and Mokhtari2016; Peng et al., Reference Peng, Zhong and Min2018). Tetraalkylammonium-modified Mnt exhibits remarkable adsorption capacity for benzene, with tetramethylammonium being particularly effective (Chun et al., Reference Chun, Sheng and Boyd2003). In other studies, Mnt modified with dodecyl trimethylammonium or dioctadecyl dimethylammonium cations has been utilized for the adsorption and degradation of chlorobenzene (Witthuhn et al., Reference Witthuhn, Klauth, Pernyeszi, Vereecken and Klumpp2006).
In recent years, the dodecyl dimethyl benzyl ammonium (DDBA) ion has been analysed by increasing numbers of researchers (Khan et al., Reference Khan, Macfie and Ray2017; Fellah et al., Reference Fellah, Hezil, Guerfi, Djellabi, Montagne and Iost2021; Ilari et al., Reference Ilari, Etcheverry, Waiman and Zanini2021; Barber & Hartmann, Reference Barber and Hartmann2022; Zeng et al., Reference Zeng, Li, Jin, Su and Yao2022). As a quaternary ammonium organic cation, DDBA is composed of a positively charged N atom linked with a lauryl group (12-C alkyl chain), a benzyl group and two methyl groups (Van de Voorde et al., Reference Van de Voorde, Lorgeoux, Gromaire and Chebbo2012; Zanini et al., Reference Zanini, Ovesenr, Hansen and Strobel2013). DDBA, as one of n-alkyl dimethyl benzyl ammonium cationic surfactants, is used to prepare organoclays (n = the number of carbon atoms in the alkyl chain). Kwolek et al. (Reference Kwolek, Hodorowicz, Stadnicka and Czapkiewicz2003) experimentally determined the adsorption capacities and structures of homologues (n-alkyl dimethyl benzyl ammonium, n = 2, 4, …, 10, 12) with different alkyl chain carbon numbers, including DDBA, in Mnt interlayers. Tahani et al. (Reference Tahani, Karroua, Van Damme, Levitz and Bergaya1999) measured the different basal spacings and interlayer structures of DDBA-intercalated Mnt with various loadings. Fellah et al. (Reference Fellah, Hezil, Guerfi, Djellabi, Montagne and Iost2021) experimentally compared the adsorption capacities of DDBA and sodium dodecyl sulfate surfactants on kaolinite surfaces. Moreover, DDBA has the characteristics of both a cationic surfactant and a disinfectant (the main component of benzalkonium; Flores et al., Reference Flores, Loveira, Yarza, Candal and Torres Sánchez2017); therefore, it was used as an environmental bactericidal and antifungal agent. DDBA-Mnt can be used as a carrier to exert a sustained antibacterial or anti-mildew effect in the environment (Yue et al., Reference Yue, Zhang, Li, Su, Jin and Qin2019). Yue et al. (Reference Yue, Zhang, Li, Su, Jin and Qin2019) measured the basal spacing of DDBA-Mnt and the release rate of DDBA-Mnt in water. Zanini et al. (Reference Zanini, Ovesenr, Hansen and Strobel2013) compared the differences in adsorption between DDBA and other benzalkonium main components (tetradecyl dimethyl benzyl ammonium) in Mnt interlayers. In other applications, DDBA-modified rectorite has been successfully utilized for the adsorption of phenol from solution and used as a filler additive to modify gel polymer electrolytes (Huang et al., Reference Huang, Ma, Liang and Yan2007, Reference Huang, Ma, Liang and Yan2008, Reference Huang, Ma, Liang and Yan2009, Reference Huang, Wang, He and Yang2010). The use of DDBA-modified Mnt to adsorb the antibiotics ampicillin and amoxicillin has also been studied (Li et al., Reference Li, Shi, Xia and Wang2021). The intercalation of DDBA in Mnt was used as a structure-directing agent in the synthesis of Fe-incorporated silica-pillared clays, which can be used for catalysis, adsorption and separation (Yang et al., Reference Yang, Liang, Gu and Mao2012, Reference Yang, Liang, Gu and Mao2013). However, although DDBA-modified clays (i.e. DDBA-intercalated clays) have been widely studied and applied, the most common model of the DDBA arrangement between clay layers requires more detail. The microscopic structure and dynamics of DDBA in the interlayer of Mnt have not been clearly revealed.
Molecular simulation has proven to be a powerful tool for elucidating the microscopic structures and dynamics of organoclay systems, providing complementary results to experimental findings (Cygan et al., Reference Cygan, Romanov and Myshakin2012, Reference Cygan, Greathouse and Kalinichev2021; Teich-McGoldrick et al., Reference Teich-McGoldrick, Greathouse, Jové-Colón and Cygan2015; Scholtzová, Reference Scholtzová, Cavallaro, Fakhrullin and Pasbakhsh2020). Minisini & Tsobnang (Reference Minisini and Tsobnang2005) used the COMPASS force field (Sun, Reference Sun1998) to simulate an alkylammonium surfactant with only one alkyl chain-intercalated Mnt system. Their simulation results accurately reproduced the X-ray diffraction (XRD) measurements of basal spacings reported by Hackett et al. (Reference Hackett, Manias and Giannelis1998). Heinz et al. (Reference Heinz, Vaia and Farmer2006) employed the PCFF force field (Sun et al., Reference Sun, Mumby, Maple and Hagler1994) to determine the cleavage energy at the interface of octadecyl ammonium and Mnt. The calculated cleavage energies and their Coulomb and van der Waals contributions were found to be in good agreement with experimental results (Giese & Van Oss, Reference Giese and Van Oss2002; Yariv & Cross, Reference Yariv and Cross2002). Lanson et al. (Reference Lanson, Mignon, Velde, Bauer, Lanson, Findling and del Valle2022) simulated the behaviour of alkylammonium in the Mnt interlayer using Amber–ClayFF force fields (Cornell et al., Reference Cornell, Cieplak, Bayly, Gould, Merz and Ferguson1995; Pearlman et al., Reference Pearlman, Case, Caldwell, Ross, Cheatham and Debolt1995; Cygan et al., Reference Cygan, Liang and Kalinichev2004) and obtained stable configurations; the calculated basal spacing were close to experimental values (Lagaly & Weiss, Reference Lagaly and Weiss1971; Lagaly, Reference Lagaly and Mermut1994). Liu et al. (Reference Liu, Lu, Wang, Zhou and Xu2007) investigated the swelling behaviour of alkylammonium-intercalated smectites as a function of carbon chain length by using the combined ClayFF–CVFF force field (Dauber-Osguthorpe et al., Reference Dauber-Osguthorpe, Roberts, Osguthorpe, Wolff, Genest and Hagler1988; Cygan et al., Reference Cygan, Liang and Kalinichev2004). The simulation results of the system basal spacings agree with the experimental findings of alkylammonium-intercalated smectites as reported by Laird et al. (Reference Laird, Scott and Fenton1989). Zhou et al. (Reference Zhou, Shen, Zhu, Zhu, He, Zhou and Yuan2014) utilized the same combined force field to derive the structural and dynamic properties of water-saturated cetyltrimethylammonium-Mnt. Their simulation findings were consistent with experimental basal spacing and the conformation of alkyl chains (Kung & Hayes, Reference Kung and Hayes1993; Zhu et al., Reference Zhu, Zhu, Zhu and Chen2008). Ghavami et al. (Reference Ghavami, Zhao, Javadi, Jangam, Jasinski and Saraei2017) determined the alterations in basal spacing of hexadecyltrimethylammonium-Mnt by adsorption of organic solvents such as naphthalene. They used the ClayFF–CVFF force field in combination with XRD, transmission electron microscopy (TEM) and Fourier-transform infrared (FTIR) spectrometry experimental methods to establish that the simulation outcomes closely matched the experimental results. These studies verified that the ClayFF–CVFF force field can accurately reproduce the microscopic properties of organic-intercalated clay minerals.
In this study, we focused on the interaction between DDBA and the Mnt surface. The interlayer microscopic structures and mobilities of DDBA-Mnts as a function of DDBA loadings were systematically investigated using classical molecular dynamics simulations. We obtain various interlayer configurations and dynamic information for different number densities of the interlayer DDBA.
Methodology
Organic-Mnt model
The chemical formula of Mnt is Na1.0(Al3.0Mg1.0)Si8O20(OH)4, where Na+ can be replaced by a monovalent DDBA cation, Na1.0–n(DDBA)n(Al3.0Mg1.0)Si8O20(OH)4, which was equal to one unit cell (n was the fraction of DDBA in the unit cell). The supercell used for simulation consisted of two clay platelets each containing 24 unit cells: 6 in the x-dimension, 4 in the y-dimension and 2 in the z-dimension. The system contained two interlayer spaces (Fig. 1a). Twenty-four isomorphic substitutions of Mg(II) for Al(III) were imposed in each sheet, which resulted in a layer charge density of –1e/unit cell. Figure 1b shows the molecular structure of DDBA. To probe the interlayer structure of DDBA-Mnt as a function of DDBA loading and to ensure that different interlayer arrangement configurations could be reproduced, we have specified the numbers of DDBA between layers as 6 (0.25 cation-exchange capacity; CEC), 13 (0.54 CEC), 24 (1.00 CEC) and 34 (1.42 CEC), for comparison with the experimental findings of Li et al. (Reference Li, Shi, Xia and Wang2021; denoted as DA-6, DA-13, DA-24, DA-34). Cl– was used to compensate for the net charge of the system. The initial basal spacing was set to 50 Å, which was large enough to accommodate the intercalated organic ions. The results presented are the averages of those obtained from each interlayer.
Data analysis
The distribution of interlayer species in space is characterized by the atomic number density distribution along the z-direction (normal to the clay basal plane), which is calculated by averaging the NVT (canonical ensemble with constant particle number, volume and temperature) trajectories, taking the plane where the oxygen atoms in the silicon–oxygen ring at the bottom of the Mnt layer are as the starting plane (z = 0).
The formula for the distribution of species number density between layers is as in Equation 1:
Here, $N\left({{{z\;-\;\Delta z} \over 2}, \;{{z\; + \;\Delta z} \over 2}\;} \right)$ is the average number of atoms appearing at the height interval $\left({{{z\;-\;\Delta z} \over 2}, \;{{z \;+\; \Delta z} \over 2}\;} \right)$ (Δz = 0.2 Å in this research) and S is the basal surface area.
The angle between the benzyl group and the clay surface is obtained by the following calculation formula (Equations 2 & 3):
Here, the vectors a and b are, respectively, the normal vector of the plane where the benzyl ring is located and the normal vector of the plane where the clay layer is located. |*| is the absolute value.
The radial distribution function (RDF) of species B around species A is as per Equation 4 (Allen & Tildesley, Reference Allen and Tildesley1987):
ρB is the average number density of species B, dnAB is the average number density of species B within the range r to r + dr from species A.
The self-diffusion coefficient D of interlayer species can be obtained using the Einstein relationship (Equation 5; Allen & Tildesley, Reference Allen and Tildesley1987; Chang et al., Reference Chang, Skipper and Sposito1997):
Here, N is the number of atoms and ri(t) is the centre-of-mass position of the ith atom at time t. d is the diffusion dimension. For example, d = 3 represents the total coefficient and d = 1 represents the coefficient components in the x-, y- or z-directions, respectively. The left-hand side of Equation 5 is the mean square displacement (MSD) term. The components of MSD in each individual direction are denoted as XX, YY and ZZ. The motion parallel to the x–y plane is obtained by calculating (XX + YY)/2.
Computational details
All molecular dynamics simulations were carried out by using the ‘Lammps’ package (Plimpton, Reference Plimpton1995). ClayFF–CVFF was employed to describe the interactions in the systems. The cut-off of short-range interactions was 10.0 Å. The electrostatic interaction was calculated using the Ewald summation method (Frenkel & Smit, Reference Frenkel and Smit2002) and the number of K-space vectors was determined to reach a precision of 1.0E–4.
The NPT (ensemble with constant particle numbers, pressure and temperature) simulations were performed at 298 K (temperature) and 1.01325 × 105 Pa (pressure) for 30.0 ns to obtain the basal spacing in equilibrium. After NPT simulations, the NVT simulations at 298 K were carried out for another 15.0 ns for production (structures and dynamics). The time step was set to 0.5 fs in each simulation. Every atom was allowed to move freely.
Results and discussion
Basal spacings
The calculated basal spacings of DDBA-Mnt are given in Table 1. The basal spacing of DDBA-Mnt increased as the loading of DDBA increased. The basal spacing of the DA-6 system was calculated to be 14.1 Å. A small amount of DDBA in the DA-6 system formed a monolayer configuration with a basal spacing similar to that of Na-Mnt, which contained two layers of water in the interlayer (Fu et al., Reference Fu, Zhang and Low1990; Sato et al., Reference Sato, Watanabe and Otsuka1992). The basal spacing of Mnt with a single layer of DDBA between layers obtained by Tahani et al. (Reference Tahani, Karroua, Van Damme, Levitz and Bergaya1999) was ~14.6 Å. This was in a good agreement with the simulation results.
a From Tahani et al. (Reference Tahani, Karroua, Van Damme, Levitz and Bergaya1999).
b From Li et al. (Reference Li, Shi, Xia and Wang2021).
The basal spacings of DA-13, DA-24 and DA-34 were 17.7, 22.5 and 28.2 Å, respectively. The basal spacings of DA-13, DA-24 and DA-34 were 17.7, 22.5 and 28.2 Å, respectively. Tahani et al. (Reference Tahani, Karroua, Van Damme, Levitz and Bergaya1999) obtained a basal spacing of 18 Å for the double layer. Li et al. (Reference Li, Shi, Xia and Wang2021) obtained 17.3, 22.4, 23.1, 23.4 and 29.7 Å basal spacings at various DDBA loadings (all samples above from these works were dried before XRD measurements so that there were almost no interlayer water molecules). For the same interlayer DDBA layered arrangement, the deviations between the calculated values and the experimental values of the basal spacings were all <5%, which means that they were in good agreement.
To illustrate the interlayer structure of DDBA-intercalated Mnt, the simulation focused more on the different configurations of the DDBA arrangement between Mnt layers and the corresponding ranges (plateaus) of basal spacings. The basal spacings obtained by simulation match the basal spacing ranges related to various interlayer configurations obtained in previous experiments, which verifies that the ClayFF–CVFF force field accurately reproduces the microstructure of DDBA-intercalated Mnt. DDBA cations formed layered configurations in the interlayer region, which explains the plateaus observed on the experimental basal spacing curves as a function of DDBA loading.
Interlayer structures
DDBA could form layered configurations in the interlayer region of Mnt. To demonstrate the layering behaviour, we calculated the atomic number density along the z-axis (Fig. 2). The ammonium group, lauryl group and benzyl group were represented by N atoms, chain C atoms of the lauryl group and ring C atoms, respectively. In the DA-6 system, all three groups presented single and clear peaks (Fig. 2a), suggesting a monolayer configuration in the interlayer region (Fig. 3a). With the increase of interlayer DDBA loadings, the lauryl group showed two, three and four clear peaks in the DA-13, DA-24, and DA-34 systems, respectively (Fig. 2b–d), corresponding to the bilayer, pseudo-trilayer and pseudo-quadrilayer configurations in the interlayer region (Fig. 3b–d). Based on the experimental measurements on DDBA-Mnt, Tahani et al. (Reference Tahani, Karroua, Van Damme, Levitz and Bergaya1999) proposed configurations of the DDBA arrangement: DDBA cations changed from lying flat on the surface to perpendicular to the surface with increasing loadings. The structure perpendicular to the surface (pseudo-quadrilayer) could be regarded as a paraffin-type structure. The appearance of the pseudo-quadrilayer (paraffin-type) was attributed to excess unattached DDBA in the middle region rather than vertically and attached to the clay surface (Lagaly et al., Reference Lagaly1986).
In the DA-13, DA-24 and DA-34 systems, the lauryl chains presented incompletely separated peaks, suggesting the formation of a non-monolayer structure and the presence of intersecting chains in the middle of the interlayer space. As shown in Fig 3d, the ammonium groups were attached to the upper and lower inner surfaces of the clay mineral, and the alkyl chains were alternately oblique to the (001) surface of the clay minerals, indicating a paraffin-type arrangement (Lagaly et al., 1986). This paraffin-type arrangement was commonly observed in the cases of excessive adsorption of long-chain organic matter into the Mnt interlayer space (Lagaly et al., Reference Lagaly, Ogawa, Dékány, Bergaya and Lagaly2006; He et al., Reference He, Ma, Zhu, Yuan and Qing2010; Huskić et al., Reference Huskić, Žagar and Žigon2012; Li et al., Reference Li, Khan, Xia, Lei, Zhu and Wang2019). This indicated that the pseudo-quadrilayer arrangement included the paraffin-type configuration of DDBA adsorbed on the upper and lower surfaces and the unattached DDBA configuration in the middle region.
The ammonium group also showed clear peaks in these systems. In the DA-13 and DA-24 systems, two peaks of the N atom of the quaternary ammonium were present (Figs 2b,c & 3b,c), and these peaks were close to the clay basal surfaces, indicating the strong electrostatic interaction between the ammonium and clay surface. In the DA-34 system, the three tiny N peaks in the middle correspond to the N atoms of DDBA in the middle of the interlayer space, which are due to excessive loading in the interlayer space.
Compared to the ammonium group and lauryl group, the peaks of the benzyl group were less well-defined. The curves indicated that the benzyl group was located preferentially on the surface as it was anchored by the ammonium group. Such configurations could be utilized to effectively adsorb aromatic pollutants. Thus, the orientation of the benzyl group relative to the clay surface might play a significant role in the adsorption of aromatic organic compounds on DDBA-intercalated Mnt.
The occurrence probabilities of the angle between the benzyl group and the clay surface are shown in Fig. 4a. As the density of DDBA increased, the average angle between the benzyl group and the clay surface gradually increased. In the DA-6 system, the benzyl groups intersected with the clay surface at a small angle, ranging from 15° to 45° (peaking at 32°; Fig. 4a). With high DDBA number density, the benzyl groups were tilted or even nearly perpendicular to the clay surface, and the angle of inclination ranged from 60° to 90° (peaking at 75°; Fig. 4b).
Interaction mechanisms
The interaction mechanisms of DDBA with clay mineral surfaces can be explored using the RDF for surface oxygen atoms of the clay layer (bridging oxygen; Ob) around the N atoms of DDBA. In the RDF plot (Fig. 5), the first peaks of the RDF curves in the DA-13, DA-24 and DA-34 systems split into two peaks with different heights. This is because the N atoms of DDBA were slightly off-centre from the surface six-member ring. As shown in Fig. 6, the two peaks were at ~3.8–4.0 and ~4.7–4.9 Å, denoting the nearby and far side surface O, respectively. By contrast, the DA-6 system shows a much sharper first peak, indicating that the N atoms of DDBA were located nearly above the centre of the surface six-member ring. These findings were supported by images of the overlapped trajectories of N atoms in DDBA projected onto the x–y plane (blue dots in Fig. 5).
The position of the quaternary ammonium N atom in DDBA deviated from the centre of the six-membered ring on the clay surface. This was due to the two asymmetric groups bonded to the N atom (i.e. the benzyl and lauryl chains were significantly different).
Mobility of interlayer DDBA
Figure 7 shows the self-diffusion coefficients derived from the MSDs (in Fig. S1). In each system, the self-diffusion coefficient of the lauryl chain was the highest and that of the N atom was the lowest. The total self-diffusion coefficient of N atoms in each system was extremely low, ranging from 0.05 × 10–12 to 0.7 × 10–12 m2 s–1. The total self-diffusion coefficients of the C atom in the lauryl chain in DA-6, DA-13, DA-24 and DA-34 were ~5.4 × 10–12, 1.9 × 10–12, 0.9 × 10–12 and 2.3 × 10–12 m2 s–1, respectively. The total self-diffusion coefficients of the benzyl C were ~0.8 × 10–12, 1.5 × 10–12, 0.3 × 10–12 and 0.6 × 10–12 m2 s–1 in the DA-6, DA-13, DA-24 and DA-34 systems, respectively.
In each system, the (XX + YY)/2 component of the self-diffusion coefficient was higher than the ZZ component because of the 2D geometry of the interlayer space. The densities of the organic phases formed by the lauryl chain were 0.22, 0.27, 0.35 and 0.37 g mL–1 in interlayer spaces of the DA-6, DA-13, DA-24 and DA-34 systems, respectively (the calculation details can be found in the Supplementary Information). Therefore, in the DA-6 system, the lauryl chain presented the highest (XX + YY)/2 component of the self-diffusion coefficient due to it having the lowest density.
Conclusion
The interlayer structures and dynamics of DDBA-Mnt were derived from systematic molecular dynamics simulations. The conclusions are as follows:
(1) DDBA shows layering behaviour in the interlayer region. With increased loadings, the basal spacing of DDBA-Mnt increased and DDBA formed monolayer, bilayer, pseudo-trilayer and pseudo-quadrilayer configurations in the interlayer space sequentially. The basal spacings of DDBA-Mnts matched well with previously reported XRD experimental ranges of basal spacings. Moreover, the pseudo-quadrilayer (paraffin-type) structure was confirmed as a common configuration with overloaded long-chain alkylammonium compounds. The ammonium group and benzyl group preferred to stay above the clay mineral surfaces, whereas the lauryl group, due to its flexible nature, was distributed evenly in the interlayer space.
(2) DDBA ions were anchored to the basal surface through ammonium groups, which were located slightly off-centre from the Si–O six-member ring. The benzyl group and lauryl group were oriented in parallel to the monolayer configuration, whereas in the multilayer configurations they were tilted and pointing towards the interlayer.
(3) The interlayer mobility of DDBA was found to be very low. Specifically, within the DDBA ions, the mobility of the amino group was the lowest, followed by the benzyl group, whereas the mobility of the lauryl group was the highest.
The Mnt surface has obvious confinement effects on DDBA, leading to the generation of a hydrophobic interlayer. The interlayer space of DDBA-intercalated Mnt can enhance the adsorption of other organic molecules due to the strongly hydrophobic alkyl chains and benzyl groups. DDBA, as the main component of disinfectants, is often used to sterilize skin wounds (e.g. benzalkonium chloride, the active ingredient in Band-Aids). Understanding the structure and mobility of DDBA cations in Mnt interlayers is helpful for the design of and research into adsorbing small-molecule drugs on DDBA-modified Mnt.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1180/clm.2023.38.
Acknowledgements
We are grateful to the High-Performance Computing Center (HPCC) of Nanjing University for performing the numerical calculations for this paper on its blade cluster system. We thank the help from Professor Xiancai Lu.
Financial support
This study was supported by the project funded by the China Postdoctoral Science Foundation (No. BX20190151) and the National Natural Science Foundation of China (No. 42172050 and 42002036). We acknowledge the financial support from the State Key Laboratory for Mineral Deposits Research at Nanjing University.
Conflicts of interest
The authors declare none.