3.1. Characterization
3.1.1 Fourier Transform Infrared Spectra Analysis
The Fourier transform infrared spectrum of the chitosan and the modified chitosan before adsorption is given in Fig. 1a and 1b respectively. From Fig. 1a, it can be observed from that the functional groups responsible for the adsorption on the surface of the adsorbent are the hydroxyl group and amino group at 3427cm− 1. The peaks observed at 3072, 2911, and 2734 cm− 1 were assigned to the stretching vibrations of C–H bond in methyl group while 1383 cm− 1 could be associated to the bending vibration of C–H bond in the methyl group and the band at 878 cm-1 can represents the symmetrical vibration in a chain of P-O-P and to P-C phosphorus-containing compound as seen in the figure. These functional groups are negatively charged and can attract positively charged lead metal ions. While Fig. 1b, it can be observed that the functional groups responsible for adsorption on the surface of the adsorbent are, a hydroxyl group and amino group at 3427 cm− 1region.The peaks observed at 3072, 2911 and 2734 cm− 1 were assigned to the stretching vibrations of C–H bond in methyl group while 1383 cm− 1 was assigned to the bending vibration of C–H bond in methyl group[35].
From Fig. 2a and 2b, the decrease in the number of peaks after adsorption, formation of new absorption bands, the change in absorption intensity, and the shift in wavenumber of functional groups could be attributed to complexation between metal ions and binding sites of adsorbents. The binding mechanism involved sharing of electron pair between electron donor atoms (O and N) and metal ion. FTIR suggests that both hydroxyl and amine groups are indeed the main adsorption sites in the chitosan and the modified chitosan(Kyzas and Bikiaris 2015).
Table 3
Kinetics and Isotherm Models Equations Explored
Kinetic Models
|
Equations
|
Isotherm Models
|
Equations
|
Pseudo- First Order
|
\(\frac{{d}_{qt}}{dt}= {K}_{1}\left({q}_{e}-{q}_{t}\right)\)
|
Langmuir
|
\(\frac{{C}_{e}}{{q}_{e}}=\frac{1}{{K}_{L}{q}_{m}}+\frac{{C}_{e}}{{q}_{e}}\)
|
Pseudo- Second Order
|
\(\frac{dqt}{dt}= {K}_{2}{\left({q}_{e}-{q}_{t}\right)}^{2}\)
|
Freundlich
|
\(\text{log}{q}_{e}=\text{log}K+\frac{1}{n} \times \text{log}{C}_{e}\)
|
Elovich
|
\(\left(\frac{dqt}{dt}\right)=a\text{exp}\left(-bqt\right)\)
|
Temkin
|
\({q}_{e}\) = Bln\({A}_{T}\) + Bln\({C}_{e}\)
|
Bhattacharya and Venkobachar
|
\(\text{ln}\left(1-Ut\right)=-\beta t\)
|
Dubinin–Radushkevich
|
ln\({q}_{e}\)= ln\({q}_{o}\)- Kdε2
|
Weber and Morris
|
\({q}_{t}\) = Kid t1/2 + C
|
Lineweaver-Burk
|
\(\frac{1}{{q}_{e}}= \frac{1}{{q}_{m}}+ \frac{1}{{q}_{m}{K}_{L}{C}_{e}}\)
|
Boyd Model
|
Log(\({q}_{m}\)– \({q}_{t}\)) = log(\({q}_{m}\)) –\(\left(\frac{{D}_{i}}{2.303}\right)\)
|
|
|
3.1.2. Scanning Electron Microscope (SEM)
Figures 3a, 3b, 4a and 4b show that the investigated adsorbent contains many pores in which metal ions can be trapped and adsorbed. Because of its uniform surface area, the material's shape can enhance metal adsorption[20]. Furthermore, SEM analysis revealed the significant changes on the surface morphology after modification with phosphoric acid as shown by Fig. 4a. The surface of the phosphoric acid modified chitosan is rougher than that of the chitosan. So based on the morphology, it can be concluded that this material presents an adequate morphological profile to retain metal ions.
3.1.3. X-ray Diffraction (XRD) Analysis
XRD patterns for the chitosan and modified chitosan are shown in Figs. 5a and 5b. The obtained characteristic diffraction peaks for these biosorbents are consistent with those in the literature (Kyzas and Deliyanni, 2013; Balau et al., 2004). In the X–ray diffraction spectrogram of chitosan, a broad peak at 2θ = 20.2° present due to the amorphous state of chitosan with other three characteristic peaks at 20.6°, 20.9° and 30.4° attributable to the indices were observed. In the spectrogram for the phosphoric acid modified chitosan, only a broad peak at 2θ = 20° was observed. The average crystallite size, D (nm), ofthe modified chitosan was calculated about 10.5 nm, using the Debye–Scherrer equation [37]:
$$D=\frac{{K}_{s}.{\lambda }}{B.cos{\theta }}$$
6
Where Ks is a constant (Ks = 0.9 for CuKa), λ (nm) is the wavelength (0.15405 nm for CuKa), B is the peak width of half–maximum (rad) and θ is the diffraction angle.
Phosphoric acid reduces the crystallinity of the raw chitosan and improves it via cross linking, thereby making modified chitosan less crystalline as confirmed by XRD analysis[39]. Crystallinity lowers adsorption capacity because there is less surface area and few open active sites for metal binding compared to amorphous adsorbents with large surface area and more open active sites. This reduction in crystallinity enhances the probability for an adsorbate to interact with the adsorbent due to loosen inter-chain structure and bonding, thus expose more active sites for adsorption which increased the adsorption capacity of the modified chitosan [40].
3.1.4. Differential Scanning CalorimetryThermogram
The thermogram for the chitosan and the modified chitosan are given in Fig. 6a and 6b respectively. From Fig. 6a, a broad endotherm over the temperature range of 28 and 98°C (with peak at 83.42°C) was seen, followed by a sharp endotherm which peaked at 204.71°C. The diffuse endotherm of chitosan which peaked at 83.42°C can be ascribed to loss of the absorbed water (Dhawade and Jagtap 2012; Jana et al., 2015)and 0.009824J of heat per milligram sample of the adsorbent was absorbed in the process. The second endotherm can be ascribed to polymer melting. The area under the curve (0.0021143 J/mg) represents the latent heat of melting while the peak (204.71°C) represents the melting point of the polymer [41]. From Fig. 6b, the first endotherm of the modified chitosan corresponds to relaxation transition, which peaked at 193.81oC. Endothermic relaxation is a second order reaction just like glass transition [43]. Meanwhile, second endotherm which can be attributed to polymer (phosphoric acid modified chitosan) melting picked at 240oC. It is obvious that the modified chitosan demonstrated significant stability over a wide temperature range (80-170oC).
3.2. Development of Model Equation
The design matrix containing the studied factors, their ranges and the response which is the removal efficiency (%) of Pb(II), (Ypmc) and (\({\mathbf{Y}}_{\mathbf{c}\mathbf{s}}\)) is presented in Table 2. The results obtained were then evaluated with CCD in RSM for the development of the model’s regression equation. A quadratic second order polynomial expression as suggested by the software fitted the data suitably. The model equation was selected in conformity with the sequential model sum of the square that is based on the highest order polynomial where the additional terms of the model are significant, and the model is not aliased. Both the regression coefficients (R2) and adjusted R2wereused to validate the quality of the proposed models; values closer to 1.0 confirm a good agreement between the predicted and experimental data. Thus, the correlation between experimental and predicted response is evident as indicated by the model’s R2 and adjusted R2 values of 0.94 and 0.90, which were within the desirability range. The final model’s equation for the removal efficiency of lead by modified chitosan, (\({\mathbf{Y}}_{\mathbf{p}\mathbf{m}\mathbf{c}}\)) and raw chitosan, (\({\mathbf{Y}}_{\mathbf{c}\mathbf{s}}\)) are given by Eqs. (7) and (8) respectively:
\({\mathbf{Y}}_{\mathbf{p}\mathbf{m}\mathbf{c}}\) = 74.57829+6.1607035*A + 8.2275766*B + 2.3218341*C-0.104926*D + 0.2963844*E-3.320276*A2-5.156986*B2-2.437222*C2-2.73161*D2-0.116199*E2 + 2.2462946*AB-1.761432*AC + 1.3925446*AD-2.74808*AE + 2.3489324*BC + 2.7949554*BD-0.14192*BE + 2.1239324*CD-0.725443*CE + 1.1368304*DE (7)
\({\mathbf{Y}}_{\mathbf{c}\mathbf{s}}\) = 63.72168 + 6.603598*A + 7.481349*B -0.31528*C-1.61259D-0.36278*E-0.92096*A2-4.42379*B2-0.5086*C2 + 1.594576*D2-0.8856*E2 + 2.963686*AB-0.95303*AC + 1.135561*AD-2.39256*AE + 1.960528*BC + 1.990689*BD-1.91994*BE + 1.752403*CD-0.00822*CE + 1.208189*DE (8)
The cooperative and counter effects of the studied factors were elucidated by the negative and positive coefficients before the terms in the model equation. A negative coefficient value implies that the term negatively affects Pb2+adsorption (i.e. the removal efficiency decreases), whereas a positive coefficient values mean that the term increase Pb2+adsorption in the tested range (Muluh, 2017).
Table 4
Statistical parameters obtained from the analysis of variances (ANOVA) for the
models for Pb (II) Ion % removal from waste water by the raw and modified chitosan.
Statistical Parameters
|
chitosan
|
Modified chitosan
|
Std. Dev.
|
4.95
|
4.88
|
Mean
|
59.32
|
62.72
|
Coefficient of variation, CV
|
8.34
|
7.78
|
PRESS
|
2684.98
|
2762.28
|
R- squared (R2)
|
0.91
|
0.92
|
R2 adjusted
|
0.84
|
0.87
|
Predicted R2
|
0.77
|
0.69
|
Adequate Precision
|
19.39
|
19.90
|
Table 5a: Analysis of variances (ANOVA) and lack-of- fit test for response surface quadratic model for removal of Pb (II) ions from aqueous solution by the modified chitosan
Source
|
Sum of squares
|
DF
|
Mean square
|
F-Value
|
Prob.>F
|
Comment
|
Model
|
8209.80
|
20
|
410.49
|
17.25
|
< 0.0001
|
Significant
|
A
|
1641.90
|
1
|
1641.90
|
69.00
|
< 0.0001
|
|
B
|
2928.41
|
1
|
2928.40
|
123.07
|
< 0.0001
|
|
C
|
228.50
|
1
|
228.49
|
9.60
|
0.0043
|
|
D
|
0.47
|
1
|
0.48
|
0.02
|
0.8885
|
|
E
|
3.80
|
1
|
3.80
|
0.16
|
0.6923
|
|
A2
|
612.18
|
1
|
612.18
|
25.73
|
< 0.0001
|
|
B2
|
1476.80
|
1
|
1476.80
|
62.07
|
< 0.0001
|
|
C2
|
330.25
|
1
|
330.24
|
13.88
|
0.0008
|
|
D2
|
414.35
|
1
|
414.35
|
17.41
|
0.0002
|
|
E2
|
0.75
|
1
|
0.75
|
0.03
|
0.8603
|
|
AB
|
161.19
|
1
|
161.19
|
6.77
|
0.0144
|
|
AC
|
96.43
|
1
|
96.43
|
4.05
|
0.0535
|
|
AD
|
61.95
|
1
|
61.95
|
2.60
|
0.1174
|
|
AE
|
241.26
|
1
|
241.26
|
10.13
|
0.0035
|
|
BC
|
171.48
|
1
|
171.48
|
7.20
|
0.0119
|
|
BD
|
249.56
|
1
|
249.56
|
10.49
|
0.0030
|
|
BE
|
0.64
|
1
|
0.64
|
0.03
|
0.8705
|
|
CD
|
140.20
|
1
|
140.19
|
5.89
|
0.0216
|
|
CE
|
16.36
|
1
|
16.36
|
0.68
|
0.4138
|
|
DE
|
41.29
|
1
|
41.29
|
1.74
|
0.1981
|
|
Residual
|
690.00
|
29
|
23.79
|
|
|
|
Lack of Fit
|
689.62
|
22
|
31.35
|
0.83
|
< 0.0001
|
Not significant
|
Pure Error
|
0.38
|
7
|
0.05
|
|
|
|
Cor Total
|
8899.80
|
49
|
|
|
|
|
Table 5b: Analysis of variances (ANOVA) and lack-of- fit test for response surface quadratic model for removal of Pb (II) ion from aqueous solution by the chitosan
Source
|
Sum of squares
|
DF
|
Mean square
|
F-Value
|
Prob.>F
|
Comment
|
Model
|
6852.60
|
20
|
342.63
|
14.01
|
< 0.0001
|
Significant
|
A
|
1886.46
|
1
|
1886.46
|
77.15
|
< 0.0001
|
|
B
|
2421.29
|
1
|
2421.29
|
99.02
|
< 0.0001
|
|
C
|
4.21
|
1
|
4.21
|
0.17
|
0.6811
|
|
D
|
112.49
|
1
|
112.49
|
4.60
|
0.0405
|
|
E
|
5.69
|
1
|
5.69
|
0.23
|
0.6331
|
|
A2
|
47.09
|
1
|
47.09
|
1.93
|
0.1758
|
|
B2
|
1086.72
|
1
|
1086.72
|
44.44
|
< 0.0001
|
|
C2
|
14.38
|
1
|
14.38
|
0.59
|
0.4493
|
|
D2
|
141.19
|
1
|
141.19
|
5.77
|
0.0229
|
|
E2
|
43.55
|
1
|
43.55
|
1.78
|
0.1924
|
|
AB
|
280.59
|
1
|
280.59
|
11.48
|
0.0020
|
|
AC
|
28.23
|
1
|
28.23
|
1.15
|
0.2915
|
|
AD
|
41.19
|
1
|
41.19
|
1.68
|
0.2045
|
|
AE
|
182.87
|
1
|
182.87
|
7.48
|
0.0105
|
|
BC
|
119.46
|
1
|
119.45
|
4.88
|
0.0351
|
|
BD
|
126.59
|
1
|
126.59
|
5.18
|
0.0305
|
|
BE
|
117.75
|
1
|
117.75
|
4.82
|
0.0364
|
|
CD
|
95.44
|
1
|
95.44
|
3.90
|
0.0578
|
|
CE
|
0.02
|
1
|
0.01
|
8.59
|
0.9927
|
|
DE
|
46.63
|
1
|
46.63
|
1.90
|
0.1778
|
|
Residual
|
709.13
|
29
|
24.45
|
|
|
|
Lack of Fit
|
690.19
|
22
|
31.37
|
0.60
|
0.0014
|
Not significant
|
Pure Error
|
18.93
|
7
|
2.70
|
|
|
|
Cor Total
|
7561.73
|
49
|
|
|
|
|
3.3. Statistical Analysis
The fitness of the model was investigated using the analysis of variance (ANOVA) at 95% confidence interval. The ANOVA of Pb2+ removal efficiency by modified chitosan and chitosan are presented in Table 5a and 5b respectively. The significance and insignificance of each term in the model were determined by the Fisher’s F-test and P-value. The model F-value of 17.25 and 14.01 respectively for modified chitosan and chitosan, which implies that the model is significant and the P-value is less than 0.05[23], [44]. A high value of the adjusted determination coefficient (R2-Adj = 0.84 and 0.87 for chitosan and modified chitosan respectively) was estimated. This result means that 84% and 87% of the total variation on Pb(II) adsorption data can be described by the selected model. Based on the P-values presented in Table 5a, the significant model terms were A, B, C, D, E,B2, D2, E2, AB, BC, BD, BE and CD with AC, AD, CE, DE, A2, C2 and E2 insignificant to the response. In order to simplify the model Chitosan, the insignificant terms (AC, AD, CE, DE, A2, C2 and E2) can be eliminated. And based on the P-values presented in Table 5b, the significant model terms were A, B, C, D, E, A2, B2, C2, D2, E2, AB, BC, BD, BE and CD with AC, AD, CE and DE insignificant to the response. In order to simplify the model for modified chitosan, the insignificant terms (AC, AD, CE and DE) can be eliminated. Adequate precision (AP) is basically a measure of signal to noise ratio. Ratios greater than 4 indicate that the model is adequate and can be used to navigate the design space[23], [44]. In this study, the AP ratio of 19.39 and 19.90 respectively of chitosan and modified chitosan, indicate an adequate signal and thus the model can be used to predict the responses. The coefficient of variation (CV) is the ratio of the standard deviation of the mean expressed as a percentage. For a model to be considered reliable and reproducible, it must have a CV less than 10% [23]. In this investigation, the CV value of 8.34% and 7.78% respectively of chitosan and modified chitosan reflect a good precision and reliability of the experiments.
3.4. Effects of individual variables and their interactions
It can be inferred from Tables 2, 5a and 5b that the individual effects of factors (contact time, pH, adsorbent dose and initial concentration) on lead uptake by chitosan and modified chitosan were more dominant than the effect of temperature, based on their high F-values of 0.23 and 0.16 (initial concentration), 0.17 and 9.60 (adsorbent dose), 99.09 and 123.08 (pH), 77.15 and 69.01 (contact time) and 4.60 and 0.02 (temperature) respectively of the prepared adsorbents. The two factor interaction effects of initial concentration and pH (F-value of 3.36 and 12.28), initial concentration and contact time (40.70 and 5.64) respectively of the prepared adsorbent are more pronounced in the interaction. From these results, it can be concluded that the adsorption of lead from aqueous solution by chitosan and phosphoric acid modified chitosan is highly dependent on initial concentration of adsorbate and the combined effects of the contact time or pH and initial concentration of adsorbate, due to their very high F-values. Figure 7 to 9 depicts the three-dimensional response surfaces of the interaction effects of the process variables on the adsorption of lead and can be noticed that an increase in the removal efficiency of Pb (II) by prepared chitosan and the modified chitosan can be observed with decrease in initial concentration, increase in the mass of adsorbent, in the contact time and increase in pH to a maximum of about 10. The interaction effect of initial concentration and pH can be observed to show a lopsided effect on Lead removal efficiency.
3.5. Adsorption Isotherm Studies
Adsorption isotherms are useful in determining the nature of the interaction between the adsorbate and the adsorbent. To determine the most befitting isotherm model for the sorption of lead, the equilibrium data obtained using optimized values for the studied parameters were fitted to the Langmuir, Freundlich, Temkins, Dubinin-Radushkevich and Lineweaver-Burk isotherm models. The results of the isotherm parameters from the plots are presented in Table 6, and Freundlich isotherm model (Fig. 10) with coefficient of determination (R2 = 0.9286 and 0.9323 for chitosan and modified chitosan respectively) and Lineweaver-Burk isotherm model (Fig. 11) with coefficient of determination (R2 = 0.9518 and 0.9664 respectively for chitosan and modified chitosan) appeared to be much more relevant in fitting the experimental data than the other models studied.
The fitness of the Freundlich model implies that the adsorption of Pb2+ ions from bulk solution assumes heterogeneous surface energies. It can be observed from Table 6 that the n values gotten from the Freundlich plot were greater than unity for the studied concentration range; this suggests that the adsorption conditions were favourable. All the RL (dimensionless separation factor) values from the Langmuir plot lie between 0 and 1 indicating that the adsorption process is favourable for the under studied conditions [45]
Table 6: Adsorption isotherm parameters for different models studied
Isotherm Models
|
Modified chitosan
|
Chitosan
|
Langmuir model
|
|
|
qm (mg/g)
|
2.600
|
0.2413
|
KL
|
0.0374
|
0.0248
|
R2
|
0.7801
|
0.6660
|
RL
|
0.3484
|
0.4462
|
Freundlich model
|
|
|
N
|
1.718
|
1.701
|
Kf
|
0.166
|
0.154
|
R2
|
0.9323
|
0.9286
|
Temkin model
|
|
|
(Lmg-1)
|
0.423
|
0.370
|
B
|
0.496
|
0.407
|
R2
|
0.8021
|
0.7969
|
Dubinin-Radushkevich model
|
|
|
(mg/g)
|
1.224
|
1.153
|
E (KJmol-1)
|
176776.7
|
176776.7
|
(mol2KJ-2)10-6
|
4 10-6
|
4 10-6
|
R2
|
0.7909
|
0.7791
|
Lineweaver-Burk
|
|
|
R2
|
0.9664
|
0.9518
|
qm (mg/g)
|
2.600
|
0.2413
|
3.6. Adsorption Kinetic Studies
Four linearized forms of kinetic models viz; pseudo-first order, pseudo-second order, Elovich kinetic model, Bhattacharya and Venkobachar kinetic model and Weber and Morris intra-particle diffusion kinetic models have been used to analyze the collected experimental data and the linear plots obtained for the different kinetic models are presented in Fig. 15to20 with the kinetic parameters and coefficients of determination values obtained summarized in Table 7. The pseudo-second-order kinetics model fitted well to the adsorption data as observed from the high R2 value (R2 > 0.99). More so, it best predicted the value for qe than those of first-order model given in Table 7.Therefore, it can be stated that the rate-controlling step is chemisorption, involving valence forces through sharing or exchange of electrons between the adsorbent surface and adsorbate ions with no involvement of a mass transfer in solution[46]. A deviation of plot from origin was seen in the intra-particle diffusion model indicating that it is not the sole rate-limiting step of the reaction. This deviation from the origin could be because of difference in the rate of mass transfer in the initial and final step of the adsorption process [47]. The C (mg g− 1) values of the intra-particle diffusion model indicate the boundary layer thickness of Pb2+ on the adsorbent surface. The C values from the parameters were observed to be smaller than the experimental qe values (Table 7) and this means that Pb2+ uptake on these adsorbents involved both surface adsorption and partitioning inside adsorbent components such organic matter and aliphatic amino groups, among other things.
Table 7: Summary of parameters for kinetic models studied

3.5. Process Optimization/Validation
The function of desirability with preparation conditions was applied. Based on economical perspective of production process, the variables were set within the studied range and goals were targeted to achieve maximum possible removal percentage by the software. The optimum condition chosen by the software was: Contact time, pH, Adsorbent dose, Temperature and Adsorbate initial concentration of 150min, 10, 2.0g, 55°C, and 50mg/L respectively to obtain highest removal efficiency. At this optimum condition, removal percentage obtained by chitosan and the modified chitosan biosorbents were 84 and 87% respectively. The experimental values obtained were in good agreement with the predicted results by the software with relatively small errors of 0.38 and 2.83% respectively. A plot of studentized residuals versus run order was tested and the residuals were scattered randomly around ± 3sigma (Fig. 21). This was an indication of better fit for both the models with the experimental data and shows that the model fits well to optimize the independent variables for the removal of Pb(II)[44].