3.1 Hydrogels characterization
The preparation process of three series of hydrogels was illustrated in Figure. 1. APEG is an important raw material of polycarboxylic acid-type water-reducing agent for cement paste (Li and Wang 2021; Sun et al. 2011). When it is introduced into P(AA-co-AM) hydrogels, its ether and hydroxyl groups are prone to form hydrogen bonds with carboxylic acid groups and carbonyl groups. PVP is a water-soluble linear polymer, whose carbonyl groups can also form hydrogen bonds with carboxylic acid groups of PAA (Han et al. 2021; Wang et al. 2013). According to Re values in Table 1, HGA2, HGB2, and HGC3 hydrogels were selected for further study.
Figure 2a showed the FTIR spectra of HGA2, HGB2, and HGC3 hydrogels. The wide absorption band detected at 3446 cm-1 was due to the stretching vibration of O-H and N-H (Tang et al. 2021; Wang et al. 2020). The bands between 2856 and 2971 cm-1 were assigned to the C-H symmetric and asymmetric stretching vibrations (Sarmah and Karak 2020; Zhang et al. 2020). The absorption band at1722 cm-1 was attributed to the C = O asymmetric in the carboxylate anion (Wang et al. 2020). The absorption band at1663 cm-1 was responsible for the stretching vibration of C = O in amide groups or pyrrolidone (Wang et al. 2021a; Zhang et al. 2020). The absorption bands at 1454, and 1404 cm-1 were related to -COO- stretching vibrations of carboxyl groups (Thakur et al. 2016; Wang et al. 2013). The absorption band at1283 cm-1 corresponded to the C-N vibration of pyrrolidone (Wang et al. 2021a). The presence of bands at 1163 cm-1 was associated with the C-N stretching vibration of MBA (Melo et al. 2018; Sarmah and Karak 2020), supporting the presence of the cross-linker.
The cross-section morphology images of HGA2, HGB2, and HGC3 hydrogels before adsorption were presented in Fig. 2b-d. It could be clearly observed that three samples had well-connected 3D porous structures with a pore diameter in the range of 10–30 µm, which provided a large surface area and enough channels for the adsorption and diffusion of water and MB molecules. Moreover, the HGB2 and HGC3 hydrogels indicated thicker pore-wall than the HGA2 hydrogel, which might be caused by the additional non-covalent crosslinks.
3.2 Swelling and adsorption behaviors
The moderate swelling capacity of hydrogels is necessary and beneficial to the transportation and adsorption of dye molecules. Nevertheless, excessive amounts of water in hydrogels are also disadvantaged for desorption and re-usage. Figure 3a showed the swelling ratio of HGA2, HGB2, and HGC3 hydrogels. Their SR values of them were 582.6, 213.2, and 234.5 g g-1, respectively. It was reasonable to deduce that the additional non-covalent crosslinks in HGB2 and HGC3 hydrogels increased their crosslinking densities and consequently reduced their swelling capacities.
The pH of MB solution had an important influence on the adsorption of hydrogels. As shown in Fig. 3b, the removal efficiency of three samples exhibited the highest values in the pH environment of 7–9, which was consistent with other reports (Sarmah and Karak 2020; Wang et al. 2020). Close to neutrality or under weak alkaline conditions, Re values could reach higher than 93%, which meant the blue MB solution could be changed to a transparent solution with slight cyan after adsorption (Fig. 3c). In the acidic environment, anionic adsorption sites in hydrogels were mostly occupied by protons (H+), which suppressed the adsorption of MB molecules (Melo et al. 2018). While in a weak alkaline environment, the most of anionic adsorption sites in hydrogels were exposed to MB molecules due to deprotonation of the carboxylic groups, which facilitated the electrostatic attractions between MB and carboxylate (-COO-) (Chen et al. 2021). When the pH further exceeded 9, the charge screening effect of the excess Na+ led to the reduction of their Re values (Sarmah and Karak 2020).
The effect of initial MB concentration on the adsorption of HGA2, HGB2, and HGC3 hydrogels was also explored. As illustrated in Fig. 3d-f, the Re values of them decreased with the increasing MB concentration from 10 to 1000 mg L-1. Such a finding might be attributed to the availability of active adsorption sites. At a lower initial MB concentration, there were enough active adsorption sites for MB molecules, which resulted in higher removal efficiency of MB. In contrast, at higher initial MB concentration, the available adsorption sites for MB molecules became insufficient due to the competition between more and more MB molecules, which led to lower removal efficiency of MB. On the contrary, the adsorption capacities of three hydrogels increased with the increasing MB concentration from 10 to 800 mg L-1. Then the adsorption capacities of HGA2 and HGB2 hydrogels decreased with the further increase of MB concentration. It's noteworthy to mention that the adsorption capacity of HGC3 hydrogel continued to ascend even in the initial MB concentration of 2000 mg L-1, indicating an outstanding adsorption ability for MB. The highest adsorption amounts of three hydrogels were 1260, 913, and 1996.8 mg g-1, respectively.
3.3 Influence of contact time and kinetics parameters
The influence of contact time on the MB adsorption capacity of three hydrogels was shown in Fig. 4a. It could be seen that the adsorption rates were very fast in the initial 100 min and then slowed down in the following adsorption stage. To better understand the adsorption kinetics, the pseudo-first-order model (Zhou et al. 2014), pseudo-second-order model (Omer et al. 2021; Ren et al. 2022), Elovich model (Melo et al. 2018), liquid film diffusion model (He et al. 2018), and intra-particle diffusion model (Tang et al. 2021) were adopted to analyze these adsorption data (Fig. 4b-f). The linear forms of these kinetic models were provided in equations (4)-(8) respectively.
The pseudo-first-order model:
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, qt (mg g-1) is the adsorption capacity of hydrogels at time t (min), and K1 (min-1) is the rate constant of the model.
The pseudo-second-order model:
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, qt (mg g-1) is the adsorption capacity of hydrogels at time t (min), and K2 (g mg-1 min-1) is the rate constant of the model.
Elovich model:
$${q_t}=\frac{1}{\beta }\ln (\alpha \beta )+\frac{1}{\beta }\ln t$$
6
where α (mg g-1 min-1) is the initial adsorption rate, qt (mg g-1) is the adsorption capacity of hydrogels at time t (min), and β (g mg-1) is a constant related to the desorption process (Liu et al. 2022).
Intra-particle diffusion model:
$${q_t}={K_i} \times {t^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0pt} 2}}}+C$$
7
where qt (mg g-1) is the adsorption capacity of hydrogels at time t (min), Ki (mg g-1 min-1/2) is the intra-particle diffusion rate constant, and C (mg g-1) is a constant related to the thickness of the boundary layer.
Liquid film diffusion model:
$$\ln (1 - \frac{{{q_t}}}{{{q_e}}})= - {K_d} \times t$$
8
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, qt (mg g-1) is the adsorption capacity of hydrogels at time t (min), and Kd (min-1) is the liquid film diffusion rate constant.
All the above-involved kinetic parameters and the correlation coefficient (R2) values were listed in Table 2.
Table 2
Kinetic parameters for the MB adsorption of HGA2, HGB2, and HGC3 hydrogels.
Models
|
Parameters
|
HGA2
|
HGB2
|
HGC3
|
Pseudo-first-order
|
R2
|
0.869
|
0.972
|
0.837
|
K1 (min-1)
|
0.00196
|
0.00243
|
0.00136
|
qe (mg g-1)
|
4.457
|
8.607
|
3.867
|
Pseudo-second-order
|
R2
|
0.999
|
0.998
|
0.999
|
K2 (g mg-1 min-1)
|
0.00284
|
0.00100
|
0.00390
|
qe (mg g-1)
|
18.741
|
19.316
|
18.592
|
h0 (mg g-1 min-1)
|
0.997
|
0.373
|
1.348
|
Elovich
|
R2
|
0.868
|
0.986
|
0.996
|
α (mg g− 1 min− 1)
|
80.786
|
4.546
|
272.336
|
β (g mg− 1)
|
0.584
|
0.414
|
1.055
|
Intra-particle diffusion
|
|
|
|
|
Step I
|
R2
|
0.957
|
0.986
|
0.969
|
Ki (mg g-1 min-1/2)
|
0.696
|
0.545
|
0.296
|
C (mg g-1)
|
7.446
|
7.020
|
13.287
|
Step II
|
R2
|
0.776
|
0.992
|
0.953
|
Ki (mg g-1 min-1/2)
|
0.0471
|
0.247
|
0.0673
|
C (mg g-1)
|
16.394
|
11.237
|
16.254
|
Liquid film diffusion
|
R2
|
0.776
|
0.948
|
0.954
|
Kd (min-1)
|
0.00252
|
0.00305
|
0.0022
|
As revealed in the fitting diagrams of the pseudo-first-order and second-order kinetics (Fig. 4b-c) and the R2 values of them in Table 2, compared to the pseudo-first-order model, the pseudo-second- order kinetic model demonstrated R2 values close to 1 and a remarkable consistency between the theoretical qe values in Table 2 and the experimental equilibrium amounts in Fig. 4a. Therefore, MB adsorption of the three hydrogels could be well described by the pseudo-second-order kinetic model, which indicated that the adsorption process of MB in these hydrogels was dominantly controlled by the chemisorption via the exchange or sharing of electrons between MB molecules and polymer chains of hydrogels, such as electrostatic attraction, hydrogen-bonding interaction, and n-π interaction (Chen et al. 2021; Sarmah and Karak 2020; Thakur et al. 2016; Ulu et al. 2022; Zhou et al. 2014). Based on the kinetic parameters of the pseudo-second- order kinetic model in Table 2, the initial adsorption rate (h0) of the three samples could be obtained from the following Eq. (9) (Viana et al. 2020):
$${h_0}={K_2} \times {q_e}^{2}$$
9
It could be found that HGC3 had the highest h0 value and HGB2 had the lowest h0 value among the three hydrogels. The lowest h0 value of HGB2 was probably related to its lowest SR value. In addition, according to the Elovich model, which describes the chemical adsorption mechanism for heterogeneous adsorption processes, the initial adsorption rates (α) of the three hydrogels in Table 2 also had a similar order. It could be inferred that the introduction of PVP might be more feasible to enhance the initial adsorption rate of HGA hydrogels than the introduction of APEG.
To further understand the adsorption mechanism of MB in these hydrogels, the intra-particle diffusion model and the liquid film diffusion model were also used to evaluate these kinetic results from the perspective of MB molecule diffusion. It could be viewed that the fitting curves of the intra-particle diffusion model in Fig. 4e exhibited two linear regions with different slopes for three samples. Furthermore, these fitting plots did not pass through zero. These results indicated the overall adsorption process might be controlled by multiple mechanisms (Sarmah and Karak 2020). The step I straight lines with large slopes originated from the rapid transfer of MB molecules on the external surface of hydrogels through boundary layer diffusion, while the step II straight lines with small slopes associated with the intra-particle or pore diffusion of MB molecules in hydrogels (Ren et al. 2022; Zhou et al. 2014). On the contrary, the fitting plots of the liquid film diffusion model displayed poor R2 values lower than those of the intra-particle diffusion model (Table 2) and did not pass through the origin (Fig. 4f), which meant that the liquid film diffusion model was not applicable for the MB adsorption in the three samples.
3.4 Adsorption isotherms
Figure 5a showed the adsorption isotherms of HGA2, HGB2, and HGC3 hydrogels for MB. It could be found that the adsorption capacities of hydrogels increased with the increase of the MB equilibrium concentration. To better understand the interaction between these hydrogels and MB molecules, the equilibrium data were evaluated by Langmuir (Zhou et al. 2014), Freundlich (Maijan et al. 2020; Omer et al. 2021), and Dubinin-Radushkevich (D-R) (Ren et al. 2022) isotherm models. The plotting curves were displayed in Fig. 5b-d. The linear forms of these isotherm models were expressed in equations (10)-(12), respectively. All the calculated isotherm model parameters and their R2 values of HGA2, HGB2, and HGC3 hydrogels for the MB adsorption were recorded in Table 3.
Table 3
Isotherm model parameters for the MB adsorption of HGA2, HGB2, and HGC3 hydrogels.
Models
|
Parameters
|
HGA2
|
HGB2
|
HGC3
|
Langmuir
|
R2
|
0.994
|
0.999
|
0.997
|
KL (L mg-1)
|
0.0133
|
0.0092
|
0.0053
|
qm (mg g-1)
|
1838.2
|
1207.3
|
2483.2
|
RL
|
0.07 ~ 0.883
|
0.098 ~ 0.916
|
0.159 ~ 0.95
|
Freundlich
|
R2
|
0.97772
|
0.98286
|
0.99447
|
Kf [(mg g-1) (L mg-1)1/n]
|
28.03629
|
22.2564
|
19.60334
|
1/n
|
0.80359
|
0.6702
|
0.79926
|
Dubinin-Radushkevich
|
R2
|
0.676
|
0.479
|
0.488
|
KDR (mol2 kJ-2)
|
0.766
|
0.522
|
0.916
|
qm (mg g-1)
|
430.8
|
253.9
|
376.2
|
E (kJ mol-1)
|
0.808
|
0.979
|
0.739
|
Langmuir model:
$$\frac{{{C_e}}}{{{q_e}}}=\frac{1}{{{q_m} \times {K_L}}}+\frac{{{C_e}}}{{{q_m}}}$$
10
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, qm (mg g-1) is the maximum adsorption capacity of hydrogels, Ce (mg L-1) is the MB equilibrium concentration, and KL (L mg-1) is the Langmuir constant, which associates with the free energy and affinity of adsorption.
Freundlich model:
$$\ln {q_e}=\ln {K_f}+\frac{1}{n}\ln {C_e}$$
11
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, Ce (mg L-1) is the MB equilibrium concentration, and Kf ((mg g-1) (L mg-1)1/n) is the Freundlich constant, and 1/n is the heterogeneity factor.
D-R model:
$$\ln {q_e}=\ln {q_m} - {K_{DR}} \times {\varepsilon ^2}$$
12
$$\varepsilon =RT\ln (1+\frac{1}{{{C_e}}})$$
13
where qe (mg g-1) is the equilibrium adsorption capacity of hydrogels, qm (mg g-1) is the maximum theoretical adsorption capacity of hydrogels, Ce (mg L-1) is the MB equilibrium concentration, KDR (kJ2 mol-2) is the adsorption constant related to the adsorption energy, R (8.314 J K-1 mol-1) is the universal gas constant, and T (K) is the testing temperature in Kelvin. Additionally, the mean adsorption free energy (E, kJ mol-1) of MB can be obtained from Eq. (14) (Ren et al. 2022).
$$E=\frac{1}{{\sqrt {2{K_{DR}}} }}$$
14
Based on the E value, the nature of the isotherm can be judged as physisorption (E < 8 kJ mol-1) and chemisorption (E = 8 ~ 16 kJ mol-1).
The applicability of these isotherm models to the adsorption of HGA2, HGB2, and HGC3 hydrogels for MB could be determined through R2 values. As illustrated in Table 3, the Langmuir model demonstrated R2 values between 0.994 and 0.999, which were larger than those of the Freundlich and D-R models. This revealed that the Langmuir model was more appropriate to describe the MB adsorption process on the three hydrogels, which suggested that the adsorption of MB molecules on these hydrogels occurred in a monolayer coverage (Malatji et al. 2020). The separation constant (RL) of the Langmuir model could be calculated via Eq. (15) (Melo et al. 2018).
$${R_L}=\frac{1}{{1+{K_L}{C_0}}}$$
15
As noted in Table 3, all the RL values were in the range of 0 ~ 1, which meant that the MB adsorption on the three hydrogels was a favorable process. In the case of the Langmuir model, the maximum theoretical adsorption capacities of them respectively were 1838.2, 1207.3, and 2483.2 mg g-1, which were higher than or comparable to that of other hydrogel-based adsorbents in literature (Table 4). This comparison demonstrated that the HGC3 hydrogel might be an effective adsorbent for MB removal.
Table 4
Comparative assessment of adsorption capacities with other hydrogel-based adsorbents.
Adsorbents
|
qm (mg g-1)
|
Reference
|
HGA2
|
1838.2
|
This work
|
HGB2
|
1207.3
|
This work
|
HGC3
|
2483.2
|
This work
|
TiO2/sodium alginate crosslinked polyacrylic acid composite hydrogel
|
2257.36
|
(Thakur et al. 2016)
|
Chitosan/κ-carrageenan/acid-activated bentonite composite
|
18.8
|
(Ulu et al. 2022)
|
Partially hydrolyzed polyacrylamide/cellulose nanocrystal nanocomposite hydrogel
|
358.42
|
(Zhou et al. 2014)
|
Hydrophobic starch based amphoteric hydrogel
|
133.65
|
(Sarmah and Karak 2020)
|
Cellulose nanowhiskers/ chitosan-g-poly(acrylic acid) hydrogel
|
2074
|
(Melo et al. 2018)
|
Polyacrylamide/sodium alginate hydrogel
|
90.9
|
(İsmail and Gökçe Kocabay 2021)
|
Functionalized cotton charcoal/chitosan biomass-based hydrogel
|
590.72
|
(Fan et al. 2022)
|
Villi-like poly(acrylic acid) based hydrogel
|
2249
|
(Tang et al. 2021)
|
Polyacrylamide/chitosan/Fe3O4 composite hydrogel
|
1603
|
(Zhang et al. 2020)
|
Semi‑interpenetrated PVA/PAMPS hydrogel
|
14.7
|
(Omer et al. 2021)
|
CMC-cl-pAA hydrogel
|
1109.55
|
(Malatji et al. 2020)
|
P(NIPAM-co-AAc)/MoS2 hydrogel
|
1258
|
(Yang et al. 2021)
|
CMC/ polyacrylamide/MMT hydrogel
|
410
|
(Ma et al. 2021)
|
3.5 Regeneration
From a practical point of view, the regeneration and reusability of adsorbents are very desirable issues. Regeneration experiments of three hydrogels were performed and the results were exhibited in Fig. 6. The Re values of them remained higher than 80% even after 5 consecutive adsorption/desorption cycles. As a consequence, the HGC3 hydrogel could be employed as a promising recyclable adsorbent for MB removal from wastewater due to its attractive properties, such as low SR value, fast initial adsorption rate, large adsorption capacity, and excellent reusability.