3.2. IR Spectroscopy
According to the IR spectroscopy data (
Figure 2), several characteristic absorption bands were observed in the spectrum of the initial bentonite. These included bands at 978 and 1197 cm
−1, which corresponded to free and associated forms of the Si–OH groups. The absorption band at 1638 cm
−1 was related to the stretching and deformation vibrations of the OH group. The bands in the region of 798–838 cm
−1 were caused by the vibrations of the Si-O bond in tetrahedral molecules, while the bands in the range 426–439 cm
−1 corresponded to the vibrations of this bond within the tetrahedron [
31]. The bands between 3000 and 3700 cm
−1 represented the stretching and deformation of OH groups in free or bound water molecules, while the absorption at 2512 cm
−1 corresponded to the formation of chelate H-bridges with OH groups. Additionally, characteristic absorption bands for Al-O groups in Al
2O
3 (1111 cm
−1) and Si-O-Al stretches (645–669 cm
−1) were also observed as well as those related to O-C-O vibrations in CO
32− ions (1197 cm
−1) [
32].
The IR spectra of bentonite samples treated with amphoteric surfactants showed absorption bands that were characteristic of the original mineral form. However, there was a decrease in the number of free -OH groups, which was likely due to the distribution of surfactant molecules within the interlayer spaces of bentonite and the displacement of water that was previously located there. The absorption bands at 1465, 2853–2927, and 1378 cm
−1 corresponded to vibrations of the CH
2 bond and CH
3, respectively, indicating the presence of organic compounds within the structure. Additionally, the appearance of absorption bands at 1542, 1320, and 1583 cm
−1 was related to amino groups and carboxylate acid residues (
Figure 2).
For organoclays synthesized with the participation of nonionic surfactants, a similar situation was observed as with the modification by amphoteric surfactants (
Figure 3).
In addition, there was a shift to the long-wave region of hydroxyls (1608 cm−1), carbonyls (-C=O 1718 cm−1), and absorption bands 2855–2929 cm−1, belonging to asymmetric and symmetric stretching vibrations of the CH2− group of the adsorbed surfactant.
3.3. The Acid–Base Properties of the Original and Modified Bentonite
When characterizing the cation exchange capacity of bentonite and the organoclays, the acid–base properties of the samples were determined. The most significant factor for the cation exchange process is the exchange acidity, which is determined in a salt extract. H
+ and Al
3+ ions from this type of acidity are the first to enter into exchange reactions and compete with exchange cations. The results of the acid–base property determination are presented in
Table 2.
The pH values of the aqueous and salt suspensions of the original and modified forms of bentonite were alkaline, ranging from 8.19 ± 0.02 to 10.92 ± 0.04. The original form of bentonite had a slightly alkaline environment (aqueous pH of 8.38 ± 0.04). The salt extract had a pH of 9.78 ± 0.01, indicating the presence of some exchange acidity due to the presence of mobile exchange ions in the material.
The modification of bentonite with different types of surfactants (amphoteric and nonionic) affected the dehydration of the surface functional groups of the bentonite. The pH of the solutions of the aqueous and salt extracts from the modified samples shifted toward the alkaline region to varying degrees
3.4. Isotherm Studies
An accurate representation of the dynamic adsorption process of a solute on an adsorbent depends on the description of its equilibrium distribution between the two phases. Equilibrium is reached when the amount of solute adsorbed on the adsorbent is equal to the amount disturbed in the solution. The equilibrium concentration of the solution remains constant. The graph displaying the dependence of the concentration of the adsorbed solute in the solid phase and the concentration in the liquid phase shows the equilibrium adsorption isotherm [
33,
34,
35,
36].
Sorption equilibrium isotherms are an important tool for analyzing and designing adsorption systems and provide a way to assess and interpret thermodynamic parameters. While isotherms may accurately match the experimental data under certain conditions, they can completely disagree under other conditions. This paper presents the results of fitting adsorption isotherms using the Langmuir, Freundlich, and BET models as they allowed us to find the best possible fit between the experimental and theoretical data.
3.4.1. The Langmuir Isotherm
The Langmuir adsorption isotherm theory assumes that a monolayer of the adsorbate covers a uniform surface of the adsorbent. This is represented graphically by a plateau in the isotherm. The model states that at equilibrium, the adsorption reaches a saturation point where no further absorption can occur. Equation (2) was used to analyze the experimental data on adsorption processes within the context of the Langmuir model.
where Q is the amount of absorbed cations, mg/g; Q
∞ is the maximum adsorption value of an element, mmol/g; K
L is the Langmuir constant, L/mmol; C
eq is the concentration of an element in an equilibrium solution, mmol/L.
The obtained sorption isotherms are shown in
Figure 4.
The Langmuir adsorption isotherms for Pb
2+ ions on bentonite and its modified forms with the participation of amphoteric and nonionic surfactants followed the L-type according to the Giles classification [
37]. The shape of the isotherm at low concentrations indicates that the binding strength to the sorbent is relatively strong. At the beginning of the process, there is an excess of available functional groups for binding and a high concentration of metal in the solution. As equilibrium is approached, functional groups become occupied by metal ions, and filling the adsorption surface becomes more difficult due to repulsive forces between the bound metal ions and those still present in solution [
38]. The parameters of the Langmuir isotherm can be found in
Table 3.
It was found that high correlation coefficients (above 0.97) indicated the applicability of the Langmuir model for describing the adsorption of lead ions on the surface of the initial bentonite and organoclays. Based on the value of maximum adsorption (A∞), the sorbents studied with the use of amphoteric surfactants can be arranged in the following order: bentonite, organoclay containing sodium cocoamphodiacetate, and organoclay containing disodium cocoamphodiacetate. The parameter KL, which describes the strength of the bond between metal ions and functional centers on the sorptive surface, increased in the series: bentonite < organoclay with disodium cocoamphodiacetate < organoclay with sodium cocoiminodipropionate.
In the case of modification with nonionic surfactants, the sorbents analyzed can be arranged in order of increasing A
∞ value as follows: organoclay with cocamide diethanolamine < bentonite < organoclay with lauramine oxide < organoclay with alkyl polyglucoside. Although the maximum calculated adsorption capacity of the organomineral complex involving diethanolamine cocamide may have been lower than that of bentonite, it was still absorbed quite effectively by this organoclay at relatively low concentrations of metal (
Figure 4b). Based on the K
L value, the order is: bentonite < organoclay with alkyl polyglycoside < organoclay with lauramine oxide < organoclay with cocamide diethanolamine.
In terms of maximum adsorption, organoclays synthesized using both amphoteric and nonionic surfactants formed a series in ascending order: organoclay with cocamide diethanolamine (nonionic surfactant) < bentonite < organoclay with lauramine oxide (nonionic surfactant) < organoclay with sodium cocoiminodipropionate (amphoteric surfactant) < organoclay with disodium cocoamphodiacetate (amphoteric surfactant) < organoclay with alkyl polyglucosides (nonionic surfactant).
The Langmuir constant varied in the following order for both types of surfactants: bentonite < organoclay with alkyl polyglucoside (nonionic surfactant) < organoclay with lauramine oxide (nonionic surfactant) < organoclay with disodium cocoamphodiacetate (amphoteric surfactant) < organoclay with sodium cocoiminodipropionate (amphoteric surfactant) < organoclay with diethanolamine cocamide (nonionic surfactant). Therefore, based on the KL value, both amphoteric and nonionic surfactants improved the bond strength between bluestone and the organometallic complex when compared to the original bentonite.
The data obtained on the adsorption of Pb
2+ ions by the initial form of bentonite and the synthesized organoclays from aqueous solutions can be compared with the literature data for existing analogs (
Table 4):
Comparison of the obtained data with analogous materials (within the studied range of concentrations and pH) led to the conclusion that the synthesized organoclay materials could be used as efficient sorbents for Pb2+ ions.
3.4.2. The Freundlich Isotherm
The Freundlich model is used to describe adsorption behavior in heterogeneous systems. The surface of natural sorbents is believed to be non-uniform, with several different sorption sites characterized by different energy levels [
42]. The main limitation of the Freundlich model lies in the fact that the adsorption capacity is not fully saturated, leading to a continuous increase in adsorption. Consequently, the parameters of the Freundlich equation can only provide information about the degree of surface non-uniformity, and indirectly, the number of adsorption sites [
43].
The two-parameter Freundlich sorption model is described by the exponential Equation (3):
where Q is the amount of absorbed cations, mmol/g; K
F is the Freundlich constant, l/mmol; C
eq is the concentration of an element in an equilibrium solution, mmol/L; 1/n is the empirical exponent.
This expression is characterized by a heterogeneity factor of 1/n, which allows the Freundlich isotherm to be used to describe heterogeneous systems. Theoretically, an infinite amount of adsorption can occur using this expression [
34].
The Freundlich isotherms and parameters for this model are presented in
Figure 5 and
Table 5.
Based on the correlation coefficients, the Freundlich model is less suitable for describing lead ion adsorption processes on organoclay surfaces than the Langmuir model. The coefficient KF, which measures absorption capacity, increased in the following order when the organoclays were modified with amphoteric surfactants: bentonite < organoclay with sodium cocoiminodipropionate < organoclay with disodium cocoamphodiacetate. In the case of nonionic surfactant modification, the sorbents analyzed can be arranged in order of increasing KF value as follows: bentonite < organoclay with cocamide diethanolamine < organoclay with lauramine oxide < organoclay with alkyl polyglucoside. The dimensionless parameter 1/n can be used to identify the energy heterogeneity of reaction centers on sorption surfaces, and its value can vary between 0 and 1. When 1/n is close to 0, there is high heterogeneity among the sorption sites, while when it is closer to 1, the sites become more homogeneous. The specified parameter for all analyzed sorbents varied within the range from 0.295 ± 0.003 to 0.413 ± 0.005, indicating an increase in the heterogeneity of the sorption layer.
The high correlation coefficients and 1/n values indicate that the Pb2+ ions were more effectively adsorbed onto the organoclays compared to the original bentonite form. This was likely due to a redistribution of electron density within the adsorbent atoms, leading to the formation of active adsorption sites on the surface.
3.4.3. The BET Isotherm
Both the Langmuir and Freundlich models have drawbacks: the data on sorption equilibrium at different concentrations cannot be accurately described by a single set of parameters. The Langmuir model allows us to calculate the limiting adsorption capacity, which is useful in systems where all adsorption sites have the same energy level and ability to bind heavy metals. However, the Freundlich model is better suited for describing processes occurring on heterogeneous surfaces with multilayer adsorption. It cannot calculate the maximum adsorption capacity, but it can provide information about the heterogeneity of the surface. The BET model addresses the limitations of both models by considering energy heterogeneity on porous surfaces. Adsorbed molecules often interact with each other, leading to a violation of stoichiometry and limiting adsorption to the formation of multiple layers. The molecules of the first layer adhere to the surface of the adsorbent due to the intermolecular interaction between the adsorbent and adsorbate. Each adsorbed molecule in the first adsorption layer can act as a center for adsorption of molecules in the second layer, and so on. This process leads to the formation of subsequent sorption layers, for which physical sorption is possible. Therefore, the value of A∞ becomes more accurate compared to the Langmuir model.
To describe polymolecular adsorption, the Brunauer–Emmett–Teller (BET) theory was used. The isotherm equation can be expressed as follows:
where Q is the adsorption value, mmol/g; Q
∞ is the maximum adsorption capacity of monolayer, mmol/g; C is a constant for a given adsorption system, directly related to the heat and entropy of adsorption; C, C
0 is the equilibrium and initial concentration of lead (Pb
2+) ions, mmol/L.
The BET adsorption equation in linear form is shown in Equation (5):
The results of the experiment on the adsorption of lead ions onto clay particles were analyzed using the BET equation, which was applied in a linear form (
Figure 6).
From the linear dependencies shown in
Figure 6, the values of the maximum adsorption (maximum adsorption capacity of the monolayer) (A
∞) and the constant (c) were determined. This was conducted by using the tangent of the angle of inclination of the straight lines obtained, and the size of the segments cut off on the y-axis (
Table 6).
It has been established that the correlation coefficients for describing the sorption process using the BET equation are lower than those obtained using the Freundlich and Langmuir models. This is due to the fact that the BET model is only applicable in the region of high concentrations, and therefore is not suitable for describing sorption processes at low concentrations. Based on the value of the maximum adsorption (A∞), the order of the studied sorbent materials can be arranged as follows: bentonite < organoclay with cocamide diethanolamine < organoclay with lauramine oxide < organoclay with sodium cocoiminodipropionate < organoclay with disodium cocoamphodiacetate < organoclay with alkyl polyglucoside. The constant of the BET equation—c, which characterizes the ratio of the energy of intermolecular interactions between the sorbent (initial and modified clays) and Pb2+ ions, increased in the following order: organoclay with lauramine oxide < organoclay with cocamide diethanolamine < organoclay with sodium cocoiminodipropionate < organoclay with disodium cocoamphodiacetate < organoclay with alkyl polyglucoside < bentonite. Additionally, for organoclays derived from lauramine oxide and cocamide diethanolamine, the values of c were less than 1, at −0.146 ± 0.004 and −0.071 ± 0.001, respectively. When c is less than 1, even before the formation of a monolayer of adsorbed molecules is complete, the formation of a polymolecular layer begins.
Since the sorption process is best described by the Langmuir equation (with the highest coefficients of determination), the adsorption equilibrium constants obtained were used to calculate the standard Gibbs free energy (Formula (6),
Table 7).
where ΔG is the Gibbs energy, kJ/mol; R is the universal gas constant, J/(mol K); T is the temperature, 298 K, kJ/mol, K is the Langmuir constant.
The results obtained indicate that the adsorption of Pb2+ ions on both the original bentonite form and on the synthesized organoclays occurred spontaneously in the forward direction.
In our case, we observed the highest correlation coefficients when using the Langmuir model to describe the adsorption processes. This model assumes monomolecular adsorption, heterogeneity of the adsorbent surface, and predominant interaction due to chemical forces. A number of studies have also shown that the adsorption of Cu(II) onto organoclays with zwitterionic surfactants is better described by the Langmuir model compared to the Freundlich model for adsorption isotherms [
44]. At the same time, [
45] suggested that the adsorption of Pb
2+ onto montmorillonite modified with different chain length betaines could be well-described by the Freundlich adsorption isotherm. The Langmuir and Freundlich R
2 calculated from Cu
2+, Pb
2+, Ni
2+, Cd
2+, Fe
2+, and Zn
2+ isotherms at the metal adsorption on natural zeolite (clinoptilolite) modified by nonionic surfactant Triton X-100 showed that for Pb, Cd, and Fe, the adsorption data were fitted equally well by both the Langmuir equation and the Freundlich equation. Adsorption data for Cu, Ni, and Zn were fitted better by the Freundlich adsorption isotherm [
46].
3.5. Determination of the Zeta Potential of Bentonite and Organoclays
The electrokinetic or zeta potential (ζ potential) is the voltage on the slip plane, which is the difference in charge between the bulk fluid and the fluid layer containing ions that are attached to the surface of a particle. For molecules and nanoparticles, a higher ζ potential means greater stability, as the solution or suspension will be less likely to aggregate. Conversely, when the ζ potential is lower, the attraction between particles exceeds the repulsive force, disrupting the stability of the suspension. Therefore, colloids with high ζ potentials are electrically stable, while those with low ζ potentials tend to coagulate or form flocs [
47]. For stability, it does not matter whether the variance has a positive or negative sign; however, the sign of the variance can have a huge impact on the application of the variance.
It is now recognized that a zeta potential of at least 30 mV (optimum > 60 mV) is required for complete electrostatic stabilization, while potentials between 5 and 15 mV are in the region of limited flocculation, and those between 5 and 3 mV correspond to maximum flocculation [
48]. It is important to note that the magnitude of the charge on the surface of the nanoparticles depends on the pH of the solution.
Table 8 presents the change in the zeta potential of the original bentonite and its modified versions, which were obtained using amphoteric and nonionic surfactants.
The measured zeta potential of the initial bentonite sample was −37 ± 1 mV, indicating resistance to aggregation of its negatively charged particles. When the bentonite was modified with amphoteric and nonionic surfactants, the negative values of the zeta potential increased, ranging from −49 ± 2 to −73 ± 1 mV. This indicates an increase in the stability of the bentonite suspension and a greater possibility of electrostatic interaction between the negatively charged bentonite particles and positively charged lead ions.
3.6. Possible Mechanisms of Lead Cation Sorption on Synthesized Organoclays
Modern concepts of the structure of the interface between solid mineral particles and the solution as well as the mechanisms of interaction between the solid and liquid phases are based on two main methodological approaches [
49,
50]. The first of these uses the principles and concepts of colloid chemistry, which underlie the theory of the electrical double layer [
51]. The second approach is based on the theoretical principles of coordination compound chemistry. When using this approach, the reactions of interaction between solid soil particles and solution components are considered as processes of the formation of surface complexes [
52].
According to the first approach, the more efficient absorption of lead cations by the colloidal particles of organoclays compared to the colloidal particles of the original clay mineral can be explained by the increase in the negative charge of the colloidal particles. The zeta potential, which is not equal to the adsorption potential or surface potential in the double electric layer, is often used to assess the properties of this layer. An increase in the zeta potential indicates an increase in colloid charge, which is associated with specific adsorption. In this type of adsorption, ions are adsorbed on the surface of the solid phase without releasing an equivalent number of ions into the solution, leading to the acquisition of an electric charge by the solid phase. This leads to the fact that, near the surface, due to the action of electrostatic forces of attraction, an equal number of ions with opposite charges are grouped, forming a double electric layer.
According to the Paneth–Fajans–Hahn adsorption rule, ions that are capable of completing their crystal lattice or forming a poorly soluble compound with one of the ions present in the crystal are specifically adsorbed onto the surface of a crystalline solid from an electrolyte solution. Lead is unlikely to be incorporated into the crystal structure of modified bentonite, but the formation of precipitates of lead hydroxides or salts is possible (see below). The observed acidification of the equilibrium solution during lead adsorption, compared to the organoclay solution, may indicate a predominant exchange nature of metal adsorption by organoclays. In general, the addition of salts containing higher-valence cations (less hydrated) can displace hydrogen and other metal ions from the surface of the micelles [
53].
In bentonite, the permanent negative charge is quite high, accounting for 90–95% of the total charge. This high negative charge density is due to a significant degree of isomorphic substitution, where Fe
2+ or Mg
2+ replace the Al
3+ element in the octahedral layer, and some Al
3+ ions replace Si
4+ ions in the tetrahedral layer. These permanent charges are mainly located on the basal surface and account for most of the surface charge density of bentonite. Positive charges are mainly located at the edges of the clay structures and only account for 5–10% of the total charge [
54]. According to the second approach, positively charged fragments of amphoteric surfactants interact with negatively charged groups on the surface of bentonite, while their own anionic fragments remain unoccupied and can interact with lead cations. In this scenario, each molecule of the amphoteric surfactant contains two carboxyl groups.
During the synthesis of organoclays, surfactant molecules are added and adsorbed on the clay surface as well as penetrate (at least partially) into the interlayer space. The authors in [
54] investigated the effect of the length of the carbon chain in the amphoteric surfactant, betaine, on the interlayer spacing of clays. They found that the increase in basal spacing was most significant when betaine molecules were first introduced into the galleries, but then changed only slightly as the carbon chain length increased. Therefore, an increase in betaine molecule size forced them to reorganize within the existing space rather than expanding the galleries further. The adsorption of amphoteric surfactants into interlayer spaces occurs in a monolayer, and it should be noted that we refer to small carbon chains within the molecule in this case.
The authors in [
55], based on the results of infrared (IR) and X-ray photoelectron spectroscopy (XPS) studies and sorption experiments, proposed several mechanisms of interaction between Cd
2+ and bentonite or organoclay including electrostatic interactions, ion exchange, surface complexation, and chelation. Chelation is also responsible for the adsorption of Cd
2+ onto organoclays that contain amphoteric surfactants.
Theoretically, electrostatic interactions between molecules of nonionic surfactants and clay minerals should be absent. However, the synthesis of organoclays using nonionic surfactants is a widely used process [
20,
56,
57]. Compared to ionic surfactants, clays modified with nonionic compounds have a hydrophobic surface without changing the degree of surface charge [
56]. Organoclays synthesized with nonionic surfactants are a type of organomineral complex whose formation is not based on ion exchange. Instead, the authors suggest that the nonionic surfactant molecules interact with the silicate through their functional groups and with water molecules associated with exchange cations in clay minerals through ion-dipole and hydrogen bonds [
58]. It was also noted that nonionic surfactants adsorb onto clay surfaces due to the combination of polar attractions, Van der Waals forces, and weak C-H…O (clay) bonds [
59]. This can lead to the displacement of water molecules and the formation of organic layers within the interlayer space. X-ray diffraction analysis of systems containing nonionic surfactants and bentonite has shown that these surfactants can cause a slight increase in the spacing between the layers of smectite [
58].
In our opinion, the main mechanism of interaction between organoclays based on nonionic surfactants and lead cations involves the formation of complex compounds including chelate complexes between the structures of surfactant molecules adsorbed on the clay and the cations. Hypothetical modeling of Pb
2+ adsorption onto synthesized organoclays can be represented as follows (
Figure 7).
During the sorption experiment, a decrease in the pH of the equilibrium solutions was observed as a result of the interaction between the initial working solutions of lead nitrate and the mineral or organic clays based on both amphoteric and nonionic surfactants. The pH of the solutions decreased from 8.51 ± 0.05 to 5.11 ± 0.08, depending on the type of surfactant used for modification. This was compared to the pH values of the salt extracts, which are shown in
Table 2. It is known that at pH > 6, it is possible for lead cations to form the Pb(OH)
+ complex or the Pb(OH)
2 precipitate, as shown in
Figure 7(2) and
Figure 7(3), respectively [
60,
61]. The XRD patterns also indicated the formation of lead carbonate hydroxide (hydrocerussite) as a result of Pb adsorption on Na-bentonite under similar experimental conditions [
62].