Yuanzhi Qu‡
a,
Ling Lin‡b,
Shifeng Gaoa,
Yukun Yangb,
Hongjun Huanga,
Xin Lib,
Han Rena and
Wenjia Luo
*b
aCNPC Engineering Technology R&D Company Limited, Beijing, 102206, P. R. China
bSchool of Chemistry and Chemical Engineering, Southwest Petroleum University, Chengdu, 610500, P. R. China. E-mail: luowenjia@swpu.edu.cn
First published on 11th January 2023
Adsorption of polymers on clay in aqueous solutions has wide applications in environmental, medical, and energy-related areas, but the interactions between polymers and clay under varied conditions are still not fully understood. In this study, we investigated the adsorption mechanisms of four polymers belonging to different categories, namely anionic poly(acrylic acid) (poly-AA), cationic poly(diallyldimethylammonium chloride) (poly-DADMAC), nonionic polyacrylamide (poly-AM), and the copolymer of AA and DADMAC (poly-AADADMAC). By using molecular dynamics simulations, we compared the desorption kinetics of these polymers at different temperatures and found that poly-AA and poly-AM have the weakest and strongest adsorption abilities, respectively. Polymer adsorptions are slightly more stable at higher pressures, and high salinity favors the adsorption of charged polymers. Further analysis suggests that the adsorption of anionic poly-AA is less stable than that of cationic poly-DADMAC because the latter is attracted to the negatively charged surface by direct coulombic forces, and poly-AM is stabilized by van der Waals forces and hydrogen bonds. This study provides insights on how to enhance the adsorption affinity of polymers on a clay surface and may help the design or improvement of polymer/clay nanocomposite materials.
PCNC materials are prepared by combining polymers with clay minerals. Currently, there are two major approaches to preparing PCNC materials, which are chemical grafting and physical adsorption.15 In the former approach, clay surfaces are modified chemically and the polymers are linked to the clay surface via covalent bonds.16–19 On the other hand, polymers can also be physically adsorbed on clay via van der Waal, electrostatic, and hydrogen-bonding interactions.20–23 Our aim in this study is to investigate the general principles behind polymer/clay interactions, but the strength of chemisorption depends on the number and distributions of chemical bonds, which are determined by the specific preparation method and cannot be discussed unless the details of preparation method are provided. Therefore, we focus only on physical adsorption in this paper.
Understanding the interaction between polymer and clay is the key to improving the design and application of novel PCNC materials that can have satisfactory performance under certain working conditions. For example, additives of drilling fluids must withstand high temperatures (>400 K), high pressures (>100 MPa), and high salinity conditions (up to saturation) in deep wellbores.24–26 Polymer adsorption on clay minerals can be investigated by various experimental techniques, including Fourier transform infrared spectroscopy (FTIR), X-ray diffraction (XRD), atomic force microscopy (AFM), ζ potential, particle size distribution, to measure adsorption amount, adsorption site, polymer conformation, mechanical strength, clay basal spacing, and many other quantities.27–29 Such results can be analyzed to provide insights into the adsorption mechanisms. For example, Nakatani et al. studied the effect of polymer conformation,30 Kohay et al. investigated the influence of polycation density,31 Sengwa et al. compared polymers with different structures,32 and Kádár et al. focused on clay surface characteristics.27
Although experimental techniques are versatile, quantizable, and indispensable in the research of PCNC materials, it is usually challenging to experimentally investigate the microscopic adsorption mechanisms at the molecular level. Instead, molecular dynamics (MD) and related computational methods can be used to simulate the adsorption process of polymers and predict the properties of PCNC materials. Reviews by Pereira et al.,33 Suter et al.,34 Cygan et al.,35 and Teppen et al.36 summarized early attempts to study clay materials by MD and numerical simulations. Camara et al.,37 Greathouse et al.,38 Shahriyari et al.,39 and Yang et al.40 simulated the adsorption of salts on clay, while many other researchers focused on the adsorption of organic or inorganic small molecules on clay.41–45 Wang et al.46 studied the wettability of silica, and Yu et al.47 investigated the phase equilibrium of n-pentane in quartz pores. The formation mechanisms of CO2 hydrate48,49 or CH4 hydrate50,51 in pores of various types of clay minerals were examined using molecular simulation methods, and the effect of cations was discussed.51
The adsorption of various kinds of polymers on clay have been investigated, which include but are not limited to bitumen,52 poly(butylene 2,5-furan-dicarboxylate),53 poly(lactic-co-glycolic acid),54 chitin and chitosan,55 xyloglucans,56,57 poly(ethylene glycol),58,59 poly(ε-caprolactone),60,61 and poly(ethylene oxide).62 While most researchers chose montmorillonite (MMT) to represent clay, some other studies considered other types of clay including quartz,63 kaolinite,64 and silica.65
In many studies, the authors focused on predicting the mechanical,66–68 rheological,69,70 or thermodynamical properties71,72 of the PCNC materials. On the other hand, many other researchers investigated the adsorption mechanisms of polymers and provided insights on how to enhance the adsorption. For example, Sun et al. demonstrated that divalent Ca2+ ions can enhance the adsorption of anionic polymers,73 and cationic chitosan can enhance the adsorption of neutral poly(acrylamide).74 Quezada showed that the adsorption strength of hydrolyzed polyacrylamide is proportional to salinity but inversely proportional to surface charge concentration.63 Jamil et al. suggested that ionic clouds formed on the basal plane of clay can enhance the adsorption of polymers.75 Haouti et al. compared the adsorption strengths of two cationic molecules on clay.76 Experimental studies on the effects of temperatures, pressures, and salinity on polymer adsorption are widely available. For example, Mpofu et al. showed that temperatures up to 90 °C and salt concentrations up to 0.01 M have limited effect on polyetheramine adsorption on MMT.77 Tekin et al. investigated the adsorption of cationic polyacrylamide onto sepiolite under temperatures up to 55 °C.78 However, to the best of our knowledge, few molecular dynamics studies have been reported on comparing the adsorption mechanisms of different types of polymers on clay, especially the desorption under high temperatures.
In this study, we aim to compare the adsorption mechanisms of different types of polymers on MMT. Poly(acrylic acid) (poly-AA), poly(diallyldimethylammonium chloride) (poly-DADMAC), and polyacrylamide (poly-AM). We chose these polymers both because they are widely used in drilling fluids and also because each of them represents a different category of polymers. Specifically, poly-AA, poly-DADMAC, and poly-AM represent anionic, cation, and nonionic polymers. In addition, the copolymer of acrylic acid and diallyldimethylammonium chloride (poly-AADADMAC) was selected to represent polymers that possess both positive and negative charges. In Section 2 we describe how the molecular system was constructed and the forcefield used in this study. In Section 3 we compare the desorption kinetics of these polymers under different temperatures and investigate the conditions that can affect the adsorption stability such as salinity and pressure. Finally, we provide insights on how to enhance the adsorption affinity based on the simulation results.
Our computational model is shown in Fig. 1. In this model, two identical MMT plates were placed in the xy-plane and separated by around 6.0 nm, and the lower plate was designated to be the adsorption site for polymers. This 6.0 nm layer-to-layer distance corresponds to an inter-layer spacing of around 6.6 nm. Although experimentally measured basal spacing of MMT intercalated by polymers is usually less than 5 nm,17 we found that this 5 nm spacing is too narrow to simulate the complete desorption of polymers into the solution; the polymer must be at least 3.0 nm away from both plates. On the other hand, using an even larger spacing would dramatically increase the demand on computational resources and is unnecessary.
As stated in the Introduction section, we considered four types of polymers, poly-AA, poly-DADMAC, poly-AADADMAC, and poly-AM, in this study. At the start of each simulation, nine identical polymer chains were placed above the lower plate in a 3-by-3 fashion as shown in Fig. 1. By adjusting the number of repeating units, chains of different polymers have similar lengths of around 2 nm. Parameters of polymer chains including the numbers of atoms and molecular weights are summarized in Table 1. The structures of these four polymer chains are further illustrated in Fig. S1.†
Polymer | Repeating unitsa | Number of atomsb | Molecular weight | Charge |
---|---|---|---|---|
a All chains are isotactic.b Excluding counter ions for poly-AA and poly-DADMAC.c Poly-AADADMAC is an alternating copolymer consisting of five AA and five DADMAC units. Each AA and DADMAC unit carry a negative and positive charge, respectively. | ||||
Poly-AA | 10 | 82 | 712.57 | −10.0 |
Poly-DADMAC | 6 | 152 | 759.35 | 6.0 |
Poly-AADADMACc | 10 | 167 | 988.41 | 0.0 |
Poly-AM | 10 | 102 | 712.81 | 0.0 |
Counter ions, K+ or Cl−, were added to the system where the anionic poly-AA or the cationic poly-DADMAC are present to always keep the system charge-neutral. The space between the upper and lower plates was filled with water and salt (10 wt% KCl, unless otherwise stated), where the percent by weight (wt%) means the mass of salt divided by the mass of the salt solution (excluding polymer and clay). And the 10 wt% KCl concentration is equivalent to a mole fraction of 0.026. A periodic boundary condition was imposed on the simulation box and a vacuum space of at least 5.0 nm was left between periodic images in the z-direction. In all simulations, the lower plate was fixed in space, while the upper plate is mobile and experienced a constant downward force equivalent to a pressure of 50 MPa unless otherwise stated. Therefore, all simulations in this study were effectively performed under constant-pressure conditions although the volume of the simulation box was fixed. To simplify the problem, we did not consider the influence of pH values by adding additional H+ or OH− ions. In other words, all calculations were performed at a constant pH value of 7.
![]() | (1) |
Mixing rules are used for van der Waals parameters of unlike atoms by eqn (2)
![]() | (2) |
Long-range electrostatic interactions were summed by using the particle–particle particle-mesh (PPPM) method with an accuracy of 10−4.84 As explained in Section 2.1, although all simulations were performed under constant pressure conditions because of the external forces applied on the upper plate, the volume and shape of the periodic simulation box were fixed, therefore the MD runs were under the NVT ensemble with a Nosé–Hoover thermostat.85 All simulations were performed with a time step of 1.0 fs and a total integration time of 5.0 ns. Moreover, each simulation has been repeated three times by starting from different random velocities, and the averaged results of three independent runs were reported.
Self-diffusion coefficients (Ds) were estimated from the mean-squared displacement (MSD) based on the Einstein model (eqn (3)), where N is the number of total atoms and rk(t) is the position of atom k at time t.
![]() | (3) |
Fig. 2 shows averaged results of three independent runs. We have further provided the results of each separate run in ESI Fig. S2,† which also includes the simulations under an even lower temperature (275 K).
Both the polymer/MMT distances and adsorption percentages indicate that the adsorption of all polymers becomes less stable under higher temperatures.
As shown in Fig. 2(a), the adsorption of poly-AA is the least stable among all four polymers. Even at room temperatures (300 K), poly-AA chains steadily desorb from the surface and are on average 1.5 nm above the surface at 5.0 ns. Under higher temperatures (400–500 K), poly-AA chains desorb even faster with adsorption rates dropping to nearly zero after 3.0 ns.
The adsorption of poly-AADADMAC (Fig. 2(c)) is more stable than poly-AA (Fig. 2(a)) but less stable than poly-DADMAC (Fig. 2(b)). Poly-DADMAC can retain nearly 40% adsorption ratio at 500 K after 5 ns, while for poly-AADADMAC only 15% of chains remain on the surface under the same conditions.
Surprisingly, poly-AM shows exceptional adsorption stability as shown in Fig. 2(d). Its adsorption is completely stable at 300 K and 400 K. Only in the highest temperature (500 K), poly-AM slightly desorbs from the surface but its adsorption rate never drops below 75%. While the adsorption mechanism of poly-AA, which is presumably different from those of the other three polymers, will be further investigated in the subsequent sections, at this point we can qualitatively conclude that the relative adsorption stability of these four polymers follows the order of poly-AM > poly-DADMAC > poly-AADADMAC > poly-AA.
This order of adsorption stability loosely matches the order of self-diffusion coefficients (Ds) of polymers listed in Table 2 at least for temperatures below 400 K. For example, the Ds of poly-AA are two to five times higher than that of poly-AM. A more detailed table containing Ds values of water and salt ions is provided as Table S1.†
Polymer | Ds (10−9 m2 s−1) | ||
---|---|---|---|
300 K | 400 K | 500 K | |
a All simulations in this table are at 50 MPa, in 10 wt% KCl solution. | |||
Poly-AA | 0.1014 | 0.5059 | 2.4269 |
Poly-DADMAC | 0.0513 | 0.4728 | 1.2065 |
Poly-AADADMAC | 0.0615 | 0.4678 | 0.8268 |
Poly-AM | 0.0233 | 0.1264 | 0.9311 |
Based on Fig. 3, we conclude that pressures varying between 20 and 200 MPa have a weak influence on the desorption behaviors of polymers. Except for poly-AM, the other three polymers have slightly stronger adsorption under higher pressures. A possible reason is that the elevated density of the polymer/salt solution (approximately 1.01 g cm−3 at 20 MPa and 1.08 g cm−3 at 200 MPa, see Table S2†) under high pressures reduces the mobility of polymer chains by decreasing the Ds values of polymer chains and salt ions. For example, Table 3 suggests that the Ds of poly-AA decreases from 0.63 to 0.40 from 20 MPa to 200 MPa. Table S2† further lists the density of the polymer/salt solution and the Ds values of K+, Cl−, and water. Especially for water, its Ds strictly decreases as the pressure rises.
Polymer | Ds (10−9 m2 s−1) | ||
---|---|---|---|
20 MPa | 50 MPa | 200 MPa | |
a All simulations in this table are at 400 K, in 10 wt% KCl solution. | |||
Poly-AA | 0.6325 | 0.5059 | 0.4014 |
Poly-DADMAC | 0.5737 | 0.4728 | 0.4382 |
Poly-AADADMAC | 0.5182 | 0.4678 | 0.5069 |
Poly-AM | 0.2010 | 0.1264 | 0.1563 |
For poly-AA, Fig. 4(a) suggests that the adsorption rate is inversely proportional to the salinity, at least in the first two nanoseconds. With no salt at all, the polymer quickly detaches (distance > 3.0 nm) from the surface in just half a nanosecond, while when a high concentration of KCl (20 wt%) is present, the desorption process is significantly delayed. However, after the poly-AA chains fully desorb from the surface after 2 ns, both the adsorption rates and distances are no longer dependent on salinity.
Desorption of poly-DADMAC has a more profound dependence on salinity as depicted in Fig. 4(b). A high concentration of 20% KCl can significantly enhance the adsorption of poly-DADMAC within the entire simulated time range. However, there is no obvious difference in adsorption ratios among 0%, 5%, and 10% concentrations.
For poly-AADADMAC and poly-AM shown in Fig. 4(c) and (d), salinity has a less obvious influence on the desorption rates. In the first nanosecond, desorption rates under different salinity are almost the same (difference in desorption distance <0.5 nm). At a longer timescale (2–5 ns), the desorption distances and adsorption rates show no correlation with the salinity and are more likely a consequence of random distribution due to the Brownian motion of partially desorbed polymer chains.
Therefore, we conclude that high salinity (≥20%) can enhance the adsorption of poly-AA and poly-DADMAC, which are both charged polymers. On the other hand, salinity has no apparent influence on the desorption rates of poly-AADADMAC and poly-AM. For poly-AA, a high salinity may slow down, but not prevent its desorption at high temperatures (500 K).
It is evident from Fig. 5 that poly-AA has stronger adsorption in CaCl2 than in KCl. However, the presence of CaCl2 did not prevent the polymer from desorbing from the surface after 2.0 ns. We conclude that the divalent Ca2+ cations may enhance the adsorption of poly-AA by slowing its desorption. But from the thermodynamic point of view, KCl and CaCl2 showed no difference at the high temperature (500 K) since the polymer chains eventually desorbed in both cases.
For the other three polymers shown in Fig. 5, adsorption behaviors in KCl have no apparent difference compared with those in CaCl2 solutions at 500 K.
Additional simulations at 275 K provided in Fig. S4† suggest that similar conclusions can be drawn for low temperatures. Only poly-AA is affected by the presence of CaCl2.
Unfortunately, significant desorption of polymer chains can be observed at or above 300 K. To obtain structures in which polymer chains are adsorbed on the MMT surface, we lowered the temperature to 275 K and performed a series of simulations. These results are summarized in Fig. S3,† which is previously referenced in Section 3.1.3 as a part of the investigation into salinity dependency. At 275 K and in 10% KCl solution, all polymers including poly-AA are on average within 1.5 nm from the MMT surface. Therefore, the subsequent analysis of adsorption mechanisms will be performed based on the results under the 275 K, 10% KCl condition.
We have calculated the distributions of polymer chains and salt ions along the z-direction (perpendicular to the MMT surface) averaged over each five-nanosecond trajectory and plotted them in Fig. 6. Positions of characteristic functional groups of each polymer are chosen to represent the locations of polymer chains. Specifically, the distributions of carboxyl groups (COO−) on poly-AA, tertiary ammonium nitrogen (>N<+) on poly-DADMAC, COO− and >N<+ on poly-AADADMAC, and amide (CONH2) on poly-AM are shown. Furthermore, for each polymer, Fig. 7 shows the structure of the molecular system after 3.0 ns of MD simulations. These snapshots are intended to provide an intuitive, but less rigorous, representation of the typical structures of adsorbed polymer chains. In addition, radial distribution functions (RDFs) between polymer functional groups and salt ions are shown in Fig. S5.†
![]() | ||
Fig. 7 Snapshots of (a) poly-AA, (b) poly-DADMAC, (c) poly-AADADMAC, and (d) poly-AM adsorbed on the MMT surface after 3.0 ns of MD simulations at 275 K, 50 MPa. In each panel, the figure on the right is a close-up view of the polymer-clay interfacial region highlighted by blue squares on the left. For clarity, water and ions are shown only in the close-up views. In addition to the coloring rules defined in Fig. 1, K+ and Cl− are shown as purple and light green spheres, respectively. |
Since the carboxyl groups and the MMT surface carry like charges and repel each other, Fig. 6(a) suggests that the carboxyl groups keep a minimum distance of 0.5 nm and are most likely 1.2 nm away from the surface. Meanwhile, K+ aggregate between the polymer and the MMT surface as counter ions to form an electric double-layer (EDL). The interaction between poly-AA and the MMT is dominated by coulombic forces, which rely on the existence of cations between the polymer and the clay. This mechanism can explain the observation that both the existence of divalent Ca2+ ions and higher salinity can enhance the adsorption stability because the coulombic forces are proportional to the amounts of charges. On the other hand, under high temperatures (500 K), cations with larger kinetic energies are more likely to escape the polymer/clay interfacial region, which may lead to the desorption of poly-AA.
As shown in Fig. 6(b), positively charged poly-DADMAC is attracted to the negatively charged Na-MMT surface in the adsorbed state. The density of poly-DADMAC splits into three peaks, which can be attributed to DADMAC sidechains that point towards the surface, sideways, and upwards as suggested by the close-up view in Fig. 7(b). The first peak of >N<+ appears at around 0.5 nm, much closer than the most probable 1.2 nm distance between poly-AA and MMT surface. Furthermore, the first peak of >N<+ is followed by a peak of Cl− at 0.7 nm. The >N<+ groups from poly-DADMAC and Cl− ions form an EDL similar to poly-AA, which may stabilize the adsorbed polymer chains and can explain why high salinity also favors the adsorption of poly-DADMAC.
The adsorption mechanism of poly-AADADMAC appears to be a combination of poly-DADMAC and poly-AA. As shown in Fig. 6(c), only DADMAC units are directly adsorbed on the surface (first peak at 0.5 nm), while AA units are further away from the surface (most likely at 0.7 nm). Since the polymer chain itself is charge neutral, neither cations nor anions accumulate between the polymer and the clay surface, and the peaks of Cl− almost coincide with >N<+, and K+ peaks coincide with COO− peaks. The DADMAC and AA units form a dipole with positively charged DADMAC units facing the surface. Not surprisingly, the adsorption stability of poly-AADADMAC is stronger than poly-AA but weaker than poly-DADMAC because the adsorption affinity is contributed mainly by the DADMAC units, and poly-AADADMAC is a 1:
1 copolymer of AA and DADMAC units.
The adsorption mechanism of poly-AM is completely different from the others. The sidechains of poly-AM are polar but not ionic. The amide groups carry negative partial charges and are hydrogen-bond donors. As shown in Fig. 6(d), there are two peaks of CONH2 at 0.4 nm and 0.8 nm, which may correspond to functional groups pointing downward and upward, respectively, as demonstrated in Fig. 7(d). This short amide-surface distance (0.4 nm) leads to significant contributions from van der Waals forces.
Moreover, the amide groups form hydrogen bonds with oxygen atoms on the MMT surface. If we consider the criteria to form a hydrogen bond to be H⋯O distance less than 0.25 nm and N–H⋯O angle greater than 100°, our calculations found that each poly-AM chain on average forms 1.64 hydrogen bonds with the surface at 275 K. Under higher temperatures, the average numbers of hydrogen bonds per chain are 1.52, 1.02, and 0.78 at 300 K, 400 K, and 500 K, respectively. Since each chain contains 10 amide groups (Table 1), 8% to 16% amide groups have contributions to hydrogen bonding with the surface. Under temperatures below 400 K, each poly-AM chain is anchored to the MMT surface by at least one hydrogen bond, which may be one of the major reasons for the observed strong adsorption affinity shown in Fig. 2(d).
The COO− peak which appears at 1.2 nm at 275 K (Fig. 6(a)) gradually shifts to the right in Fig. 8(a) and becomes broader at higher temperatures, meaning that poly-AA is no longer attracted to the MMT surface and finally dissolves in the solution.
In Fig. 8(b), the second and third peaks of >N<+, which represent DADMAC units pointing sideways and upwards in the original adsorbed poly-DADMAC, gradually flatten out at 400 K and 500 K. However, the first peak at 0.5 nm, which represents >N<+ groups directly attracted to the surface, does not shift or disappear under higher temperatures. Moreover, the area under the first peak is inversely proportional to the temperature (height reduced from around 1.6 mol L−1 at 275 K to 0.8 mol L−1 at 500 K). These results suggest that the DADMAC units pointing sideways and upwards in the originally adsorbed poly-DADMAC may not have significant contributions to the adsorption stability and can easily detach from the surface, possibly because long-range coulombic forces are screened in solutions with high salinity.
Similarly, Fig. 8(c) indicates that >N<+ peaks can resist the high temperature while the peak of COO− completely disappears at 500 K. This result corroborates the conclusion that adsorption stability of poly-AADADMAC is mainly contributed by cationic DADMAC units, at least under high temperatures.
Fig. 8(d) further verifies the strong affinity of poly-AM with the MMT surface. The peak at 0.4 nm, which represents amide groups directly attached to the surface via hydrogen bonds, does not shift its position under high temperatures. But the second peak gradually disappears. The height of the first peak drops from 3.75 mol L−1 at 275 K to 2.1 mol L−1 at 500 K, which matches the change in the number of average hydrogen bonds between poly-AM and the surface from 1.64 at 275 K to 0.78 at 500 K.
![]() | ||
Fig. 9 Illustrations for the adsorption mechanisms of (a) poly-AA, (b) poly-DADMAC, (c) poly-AADADMAC, and (d) poly-AM on MMT surfaces. |
Fig. 9 summarizes the adsorption mechanisms discovered in this study. Specifically, Fig. 9(a) demonstrates the EDL structure formed by K+ and COO− groups on poly-AA. Fig. 9(b) also depicts the EDL by >N<+ and Cl−. For the neutral poly-AADADMAC, Fig. 9(c) shows that poly-AADADMAC is stabilized on the surface through the contributions from both DADMAC and AA units, and Fig. 9(d) indicates that van der Waals interactions and hydrogen bonds exist between poly-AM and the surface.
Adsorption mechanisms revealed by the simulations can also provide insights into the strategies to enhance the polymer/clay interactions:
(1) Strongly polar polymers should have higher adsorption stability than anionic or cationic polymers under high temperatures (>500 K). The introduction of hydrogen bond donors can enhance the adsorption of polymer on MMT.
(2) Cationic polymers should have stronger adsorption stability than anionic polymers since they can be attracted to the negatively charged clay surface directly through coulombic forces. However, cationic polymers are less stable than neutral polar polymers, especially under high temperatures. A possible strategy to increase the adsorption stability is to introduce both cationic and polar repeating units into the polymer chains.
(3) High salinity favors the adsorption of both anionic and cationic polymers. A possible reason is that the EDL structure formed under high salinity can better stabilize the adsorbed polymer chains. Charged polymers, especially polymers with cationic units, should be recommended in high-salinity environments.
We found that temperature significantly affects the adsorption stability of polymers as all polymers are less stable under higher temperatures. Pressure is only a minor factor of adsorption stability as polymers are marginally more sable under higher pressures. High salinity favors the adsorption of both positively and negatively charged polymers but has negligible influence on neutral polymers.
Based on the simulation results, we predict that a polymer with small and strongly polar sidechains, for example with hydrogen donor sites similar to poly-AM, would adsorb firmly on clay surfaces. Positively charged polymers would also have strong adsorption stability. On the other hand, anionic polymers would have limited adsorption stability on negatively charged clay surfaces under high temperatures.
One limitation of this study is that only montmorillonite was investigated while other types of clay minerals are also widely used in drilling fluids. Furthermore, only the basal planes of clays are considered due to periodic boundary conditions, while in reality, the clay platelet edges may also be adsorption sites. Our future work will be focused on the adsorption of polymers on other types of clay under more complex situations by considering the influence of pH values and platelet edges. Nevertheless, the simulation framework and insights proposed in this study can be extended and utilized in the design and screening of polymers with high adsorption stability on clays, which will also be the aim of our future studies.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2ra07341b |
‡ Yuanzhi Qu and Ling Lin contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2023 |