3.1. Characterization
Slightly crosslinked styrene-divinylbenzene copolymers containing 1, 2, and 3% of DVB (CP-1DVB, CP-2DVB, CP-3DVB) were prepared by free radical suspension polymerization as shown in Scheme 1. The synthesis conditions (cylindrical shape of flask, high stirring rate and the presence of 1 wt% of suspension stabilizer) made it possible to prepare particles with relatively narrow size distributions. Hypercrosslinked copolymers (HCCP-1DVB, HCCP-2DVB, and HCCP-3DVB) were obtained by post-crosslinking of slightly crosslinked copolymers using CMME as an external agent. In this case, the selected amount of CMME introduced into the reaction (1 mol of CMME per 1 mol of the copolymer) ensures that, statistically, each phenyl group of the copolymer forms methylene bonds with two neighboring phenyl groups [17].
The FT-IR spectra of the copolymer before (CP-1DVB) and after the post-crosslinking procedure (HCCP-1DVB) are shown in Fig. 2. The absorption bands at 2853 and 2918 cm–1 can be assigned to C–H symmetric and asymmetric stretching vibrations. The absorption bands at 753 and 820 cm–1 correspond to C–H bending vibrations of mono- and disubstituted phenyl rings [31]. Intensities of the absorption bands at 753 and 820 cm–1 changed after the post-crosslinking procedure. The absorption band at 1704 cm–1 is traditionally assigned to C = O vibrations. However, in [32] it has been proved that hypercrosslinked polystyrenes have no functional groups, and this band is most probably related to vibrations of polysubstituted aromatic rings. Thus, it may be assumed that hypercrosslinked structure was formed under the synthesis conditions described above.
The nitrogen adsorption-desorption isotherms of the synthesized HCCPs are presented in Fig. 3(a). According to the IUPAC classification [33], these isotherms belong to type IV, which confirms micro/mesoporous structure of the samples. The hysteresis loop between the adsorption and desorption branches for HCCP-1DCB is assigned to the H2 type and suggests that this sample has a predominantly mesoporous structure. For HCCP-2DVB, the hysteresis loop is an intermediate between the H2 and H4 types, while the hysteresis loop for HCCP-3DVB can be assigned to the H4 type. Thus, the shapes of the hysteresis loops of HCCPs changed from the H4 to the H2 type with a decrease in the DVB content from 3 to 1%. This fact indicates the formation of mostly mesoporous structure during post-crosslinking of copolymers with a lower DVB content. The absence of a sharp rise in the isotherms in the region of high relative pressures indicates the absence of macropores in the HCCP structure. The pore size distribution curves of the synthesized HCCPs are presented in Fig. 3(b). In the case of HCCP-1DVB, the contribution of small mesopores with sizes of 2–4 nm is significantly higher in comparison with HCCP-2DVB and HCCP-3DVB. At the same time, HCCP-2DVB and HCCP-3DVB have a large number of micropores with sizes less than 2 nm.
Table 1 summarizes the information about porous structure of the obtained HCCPs. The values of the surface area, total pore volume and mesopore volume increase with a decrease in the DVB content from 3 to 1%. Sample HCCP-1DVB has the highest value of mesopore volume (0.53 cm3/g). Thus, a decrease in the DVB content in the initial copolymer promotes the formation of mesoporous structure in the post-crosslinking stage. This trend was noted in [28]; it may be explained by lower rigidity of the network with lower DVB content. As a result, the distances between chains are less regular, which contributes to lower homogeneity of post-crosslinking throughout the polymer.
Table 1
Porosity characteristics of the synthesized HCCPs
Sample
|
SBETa
(m2/g)
|
Vmicrob
(cm3/g)
|
Vmesoc
(cm3/g)
|
Vtotald
(cm3/g)
|
HCCP-1DVB
|
946
|
0.12
|
0.53
|
0.65
|
HCCP-2DVB
|
897
|
0.23
|
0.35
|
0.58
|
HCCP-3DVB
|
760
|
0.22
|
0.27
|
0.49
|
a Surface area calculated from nitrogen adsorption isotherms at 77.3 K using BET equation |
b micropore volume calculated using t-plot |
c mesopore volume calculated as Vtot - Vmicro |
d total pore volume calculated from nitrogen isotherm at P/P0 = 0.989, 77.3 K |
It is seen from the SEM images presented in Fig. 4 that fracture sections of the synthesized HCCPs do not indicate visible pore structure because samples are characterized by micro/mesoporous structure with only small mesopores. However, the scanning micrographs clearly show that the texture becomes smoother in going from HCCP-1DVB (Fig. 4(a, d)) to HCCP-3DVB (Fig. 4(c, f)). This indicates that as the content of DVB in the precursor increases from 1 to 3%, the texture of hypercrosslinked polymer becomes more uniform.
Thus, sample HCCP-1DVB with a pronounced mesoporous structure is promising for efficient adsorption of large organic molecules, in particular, rifampicin.
3.2. Adsorption isotherms and thermodynamics
Figure 5 displays the equilibrium isotherms of RIF adsorption on HCCP-1DVB, HCCP-2DVB and HCCP-3DVB at 298 K. At the equilibrium concentration (Ce) equal to 100 mg/L, the equilibrium adsorption capacity (qe) decreased in the order HCCP-1DVB (122.6 mg/g) > HCCP-2DVB (42.8 mg/g) > HCCP-3DVB (22.1 mg/g). The RIF molecule has a significant size and, therefore, can be adsorbed in mesopores. For this reason, the adsorbent with the largest mesopore volume (HCCP-1DVB) has the most adsorption capacity towards RIF.
The experimental data (Fig. 5) were analyzed using the Langmuir and Freundlich models described by equations (4) and (5), respectively [34]:
$${q}_{e}=\frac{{K}_{L}{C}_{e}{q}_{m}}{1+{K}_{L}{C}_{e}}$$
4
$${q}_{e}={K}_{F}{C}_{e}^{\frac{1}{n}}$$
5
where qm (mg/g) is the maximum adsorption capacity; Ce (mg/L) is the equilibrium concentration of RIF in solution; KL (L/mg) and KF ((mg/g)·(L/mg)1/n)) are the adsorption constants of Langmuir and Freundlich equations, respectively.
Table 2 lists the parameters obtained using these models. The adsorption equilibrium data for all samples are in better agreement with the Langmuir model due to higher correlation coefficients (R2 > 0.99). This model suggests the predominant monolayer adsorption on a homogeneous surface. The highest maximum adsorption capacity (qm) was 183.27 mg/g for HCCP-1DVB. This value is comparable with other adsorbents used for RIF adsorption or even significantly higher (Table 3).
Table 2
Parameters of adsorption isotherms for RIF on HCCPs at 298 K
|
Langmuir model
|
Freundlich model
|
Sample
|
qm
(mg/g)
|
KL
(L/mg)
|
R2
|
KF
((mg/g)·(L/mg)1/n)
|
n
|
R2
|
HCCP-1DVB
|
183.27
|
2.02 × 10− 2
|
0.994
|
26.82
|
3.18
|
0.987
|
HCCP-2DVB
|
89.81
|
0.91 × 10− 2
|
0.995
|
7.81
|
2.71
|
0.983
|
HCCP-3DVB
|
66.94
|
0.49 × 10− 2
|
0.994
|
2.501
|
2.08
|
0.981
|
Table 3
Comparison of adsorption capacity (qm) of various adsorbents towards RIF at 298 K
Adsorbent
|
qm
(mg/g)
|
Reference
|
Chitosan
|
66.91
|
[11]
|
Graphene oxide/Chitosan/Fe3O4
|
101.22
|
[12]
|
Fe3O4 nanoparticles
|
107.7
|
[10]
|
Zeolite − polypyrrole∕TiO2 nanoparticles
|
112.31
|
[8]
|
Activated carbon cocoa shells
|
135.2
|
[13]
|
β-cyclodextrin/mesoporous SiO2 nanospheres
|
246.8
|
[9]
|
HCCP-1DVB
|
183.27
|
This work
|
The effect of temperature on the RIF adsorption on HCCP-1DVB was investigated at 298, 308 and 318 K (Fig. 6). The equilibrium adsorption capacity increases with increasing temperature, suggesting that adsorption is an endothermic process. The equilibrium data were fitted to the above-mentioned models, and the corresponding parameters are given in Table 4. The Langmuir model was more suitable to fit the adsorption data with higher correlation coefficients (R2 > 0.99) as compared with the Freundlich model. According to [35, 36], the KL values can be used for calculating thermodynamics parameters, including the changes in enthalpy (ΔH, kJ/mol) and entropy (ΔS, J/(mol·K)) by the van’t Hoff Eq. (6) as follows:
where KL (L/mg) is the Langmuir constant, MRIF (822.09 g/mol) is the molecular weight of RIF, T (K) is the temperature, and R (8.3145 J/(mol·K)) is the gas constant.
The Gibbs free energy (ΔG, kJ/mol) can be calculated using Eq. (7):
The values of ΔS (143.56 J/(mol·K)) and ΔH (18.72 kJ/mol) were calculated according to [36]. The ΔG values calculated according to Eq. (7) were equal to − 24.06 kJ/mol at 298 K, – 25.51 kJ/mol at 308 K, and − 26.93 kJ/mol at 318 K. The positive ΔH value confirms that the RIF adsorption on HCCP-1DVB is an endothermic process. Since the ΔG value is negative, and the ΔS value is positive, it may be concluded that the RIF adsorption is a spontaneous and disordered process.
Table 4
Parameters of adsorption isotherms for RIF on HCCP-1DVB at 298, 308 and 318 K
|
|
Langmuir model
|
Freundlich model
|
Sample
|
T (K)
|
qm
(mg/g)
|
KL
(L/mg)
|
R2
|
KF
((mg/g)·(L/mg)1/n)
|
n
|
R2
|
|
298
|
183.27
|
2.02 × 10− 2
|
0.994
|
26.82
|
3.183
|
0.987
|
HCCP-1DVB
|
308
|
189.59
|
2.53 × 10− 2
|
0.995
|
32.41
|
3.401
|
0.986
|
|
318
|
195.88
|
3.25 × 10− 2
|
0.997
|
40.13
|
3.732
|
0.988
|
3.3. Adsorption kinetics
Figure 7 presents kinetic curves for the RIF adsorption on HCCP-1DVB, HCCP-2DVB and HCCP-3DVB. The RIF adsorption rate on HCCPs is high in the initial stage. 52.7, 27.1 and 14.9% of the equilibrium adsorption capacity are attained within the first 5 min for HCCP-1DVB, HCCP-2DVB and HCCP-3DVB, respectively. Then the adsorption rate decreases. More than 90% of the equilibrium adsorption capacity is reached within 30 min for HCCP-1DVB, within 60 min for HCCP-2DVB and HCCP-3DVB. The equilibrium is achieved within 90 min for all adsorbents.
The kinetic data can be fitted to the pseudo-first-order and pseudo-second-order models described by equations (8) and (9), respectively [37, 38]:
$${q}_{t}={q}_{e}(1-{e}^{-{k}_{1}t})$$
8
$${q}_{t}=\frac{{k}_{2}{q}_{e}^{2}t}{1+{k}_{2}{q}_{e}t}$$
9
where qt (mg/g) and qe (mg/g) are the adsorption capacity values at time t (min) and at equilibrium, respectively; k1 (L/min) and k2 (g/(mg·min)) are the constants for the pseudo-first-order and pseudo-second-order models, respectively.
The values of kinetic parameters are given in Table 5. The RIF adsorption on HCCPs is well described by the pseudo-second-order kinetic model due to higher correlation coefficients (R2 ≥ 0.996). The adsorption rate constants of the pseudo-second-order model increase in the order HCCP-1DVB (1.15 × 10− 3 (g/(mg·min)) > HCCP-2DVB (0.74 × 10− 3 (g/(mg·min)) > HCCP-3DVB (0.48 × 10− 3 (g/(mg·min)), which may be related to an increase in the mesopore volume. The adsorption capacity of RIF on HCCPs calculated using the pseudo-second-order model is in good agreement with the qmax value calculated by the Langmuir model.
Table 5
Kinetic parameters of RIF adsorption on HCCPs at 298 K
| Pseudo-first-order model | Pseudo-second-order model |
---|
Sample | qe, cal (mg/g) | k1 (min− 1) | R2 | qe, cal (mg/g) | k2 (g/(mg⋅min)) | R2 |
HCCP-1DVB | 156.26 | 11.89 × 10− 2 | 0.979 | 171.71 | 1.15 × 10− 3 | 0.997 |
HCCP-2DVB | 67.89 | 5.17 × 10− 2 | 0.986 | 80.94 | 0.72 × 10− 3 | 0.996 |
HCCP-3DVB | 47.33 | 3.04 × 10− 2 | 0.989 | 61.41 | 0.48 × 10− 3 | 0.996 |
3.4. Effect of pH on adsorption
RIF is a zwitterionic molecule (pKa1 = 1.7, pKa2 = 7.9) and thus can exist in different states depending on pH; this may affect its adsorption, since HCCPs have non-polar surfaces. The influence of pH on adsorption was studied in the pH range from 2 to 11. It was found that the adsorption capacity of HCCP-1DVB towards RIF reached maximum (162 mg/g) at pH = 7.2 (Fig. 8). This pH value is close to the isoelectric point of RIF (pI = 7.3) [39]. At the same time, the minimum values of adsorption capacity were 146 and 142 mg/g at pH equal to 2.1 and 11, respectively. Thus, the adsorption of RIF on HCCP-1DVB slightly depends on pH. It may be assumed that π-π stacking and hydrophobic interactions are the major driving forces for this process.
3.5. Effect of the adsorbent dosage on removal efficiency
The effect of the adsorbent (HCCP-1DVB) dosage on the RIF removal efficiency is illustrated in Fig. 9. The obtained data indicate that an increase in the HCCP-1DVB dosage leads to noticeable increase in removal efficiency of RIF. At an HCCP-1DVB dosage of 10 mg, the RIF removal efficiency was only 35%, but it increased up to 90% when the adsorbent dosage was equal to 50 mg. Finally, the removal efficiency was 99.8% at an adsorbent dosage of 80 mg.
The efficiency of RIF removal with HCCP-1DVB was also compared with the corresponding parameters for commercial hypercrosslinked adsorbents based on styrene-divinylbenzene copolymers (MN-200, MN-202, MN-270) of different porosities. As shown in Fig. 10, the RIF removal efficiency decreases in the order HCCP-1DVB (99.8%) > MN-200 (76.3%) > MN-202 (56.2%) > MN-270 (15.3%) when the adsorbent dosage is equal to 80 mg. The results obtained indicate that HCCP-1DVB can be used as an effective adsorbent for the removal of RIF from aqueous media.
3.6. Adsorbent reusability
An important characteristic of an adsorbent is the possibility of efficient regeneration and repeated use. Thus, HCCP-1DVB was filtered after equilibrium adsorption, and ethanol was employed as a regenerating agent. As illustrated in Fig. 11, after five adsorption-desorption cycles adsorption capacity of HCCP-1DVB was slightly reduced to 151 mg/g (93% of initial adsorption capacity), indicating that HCCP-1DVB has excellent reusability.