3.1. Characteristics of Adsorbents
The morphologies of the as-prepared samples were studied based on SEM images.
Figure 1A shows that the PBC had a hierarchical porous structure similar to that of carbon foam. Large amounts of chelating agents were attached to the surface of E-PBC, which resulted in significant aggregation (
Figure 1B). Rectangular crystals spread all over the surface of PBC in O-PBC (
Figure 1C). The hierarchical porous structure of PBC was retained in H-PBC and projections were observed on its surface (
Figure 1D). The existence of the porous structure agreed well with the result obtained by BET analysis.
Figure 2A shows the N
2 adsorption–desorption isotherms determined for the as-prepared samples. The adsorption–desorption curves for PBC and H-PBC were type IV adsorption/desorption isotherms with H4 type hysteresis loops, and they suggested the existence of a mesoporous structure according to a previous study [
17]. By contrast, the curves for E-PBC and O-PBC were type I adsorption/desorption isotherms due to their microporous filling. As shown in
Figure 2B, most of the pores in PBC and H-PBC were mesopores. The BET surface areas determined for PBC, H-PBC, O-PBC, and E-PBC were 360.41, 344.17, 3.66, and 1.64 m
2/g, respectively (
Table 1). Compared with PBC, the surface area was clearly much lower for E-PBC, possibly because the micropores were covered by EDTA. The surface area of O-PBC was also significantly reduced. The reason is similar to EPBC. Interestingly, H-PBC maintained a high surface area in a similar manner to PBC, thereby indicating that modification with H
3PO
4 retained the porous structure of the raw material, as well as maintaining the high surface area of the biochar through the formation of mesopores.
The element analysis of the samples is shown in
Table 1. It can be found that E-PBC and O-PBC have less C content than PBC, while H-PBC has more C content than PBC. The O content and H content showed a similar order. For H-PBC, the presence of P content proved that H
3PO
4 successfully binds to biochar. The H/C atomic ratio shows the aromaticity of a biochar due to hydrogenation and dehydrogenation process [
18]. For H-PBC, E-PBC, and O-PBC, the H/C atom ratio decreased indicated that the samples were more aromatic and stable. The (N+O)/C ratio is related to the hydrophobicity and polarity of biochar [
19]. For E-PBC, and O-PBC, the (N+O)/C atom ratio decreased revealed that increase in hydrophobicity of biochar after EDTA and NaOH modification, while H-PBC showed an opposite order.
The XRD patterns obtained for the PBC samples are illustrated in
Figure 2C. For PBC, the characteristic peaks at 29.81°, and 43° corresponded to the (112), and (101) planes, respectively [
17], which confirmed that PBC sample exhibited certain degree of graphitization. For H-PBC, the peak observed at 23° corresponded to the (002) planes of the graphite lattice [
20] and the sharp peaks at 27.2°, 33.7°, and 35.2° were other diffraction peaks on H-PBC. According to the literature [
4,
9,
14], amorphous materials do have high heavy metals adsorption capacity and as a result of H-PBC has excellent potential ability to remove heavy metals from aqueous solution. For O-PBC, a weak peak around 48° was related to (100) diffraction and another weak peak at around 43° was related to (101) diffraction by the activated biochar [
21]. For E-PBC, the presence of new peaks at 24.72° confirming the reaction between PBC and EDTA to form amorphous polymer composite [
22].
FTIR analysis was conducted to identify the functional groups on the surfaces of the prepared PBC samples and the results are shown in
Figure 2D. For PBC, the adsorption peak around 3430 cm
–1 could be attributed to the stretching of -NH groups in acidic and aliphatic compounds. The peak at 1049 cm
–1 was due to the stretching vibration of C-O [
17]. For O-PBC, the signals detected at 3430 cm
–1 and 1442 cm
–1 corresponded to the stretching vibrations of -OH and chromene, thereby, suggesting that O-PBC was successfully synthesized and the abundant functional groups may be beneficial to the adsorption of heavy metals. For H-PBC, the peaks at 1223 cm
–1, 1185 cm
–1, and 500 cm
–1 could be assigned to the stretching vibrations of C-O-C, P-O-C [
13], and C-P bonds [
23], respectively. Those functional groups are attributed to a long-range chain order owing to the strong hydrogen bonding and polar intra- and inter-molecular interactions [
14], which might strengthen binding to heavy metals. In addition, the bands at 1610 cm
–1 and 833 cm
–1 could be assigned to the stretching vibrations of aromatic C=C and =C-H, respectively [
14], thereby indicating the high aromaticity of H-PBC. For E-PBC, the absorption peaks near 2925 cm
–1, 1627 cm
–1, and 1390 cm
–1 could be attributed to the asymmetric stretching of -CH
2- [
24], C=O [
22], and -NH
2, thereby indicating that EDTA was successfully grafted onto PBC.
These results indicate that the E-PBC, O-PBC, and H-PBC were successfully synthesized according to previous studies [
14,
22,
25].
3.3. Thermodynamic Analysis
Temperature is an important parameter in the adsorption process. The effects of temperature on the adsorption of Ni(II) and Pb(II) are shown in
Figure S1D. The heavy metal removal efficiencies gradually improved as the temperature increased. Moreover, H-PBC had the highest capacity for adsorbing Ni(II) and Pb(II). When the temperature reached 45 °C, the Ni(II) and Pb(II) removal efficiencies with H-PBC were 10.6% and 39.4%, respectively. The thermodynamic behavior was investigated further using Equations (7)–(9) [
28]:
where
is the Gibbs free energy (kJ/mol),
is the entropy change (J/mol·K),
is the enthalpy change (kJ/mol),
R is the general gas constant (J/mol·K),
T is the solution temperature (K),
Ce is the concentration of the compound at equilibrium (mg/L), and
CAe is the amount adsorbed after reaching equilibrium at a specific temperature (mg/g).
The values of
and
were calculated from the plot of ln
Kc versus 1/
T (
Figure S3). The thermodynamic parameters are presented in
Table 2 and
Table 3. The negative values of
for the adsorption of Ni(II) and Pb(II) indicated the spontaneous nature of the adsorption process. The adsorption of heavy metals was more efficient at high temperatures because
decreased as the temperature increased.
was positive in the cases of Ni(II) and Pb(II) adsorption, thereby indicating the endothermic nature of their adsorption. The positive values of
reflected increased randomness on the adsorbents/solution surface and the affinity of Ni(II) and Pb(II) for the adsorbents [
23].
3.4. Adsorption Kinetics
Pseudo-first order and pseudo-second order models were employed to describe the kinetic adsorption processes. As shown in
Figure 3A,B, the pseudo-first order model did not obtain a good linear relationship, but the experimental data were in agreement with the pseudo-second order model. For Ni(II), the correlation coefficients (R
2 values) obtained for the pseudo-first order models for PBC, H-PBC, E-PBC, and O-PBC were 0.967, 0.960, 0.871, and 0.812 (
Table 4), respectively, which were lower than those for the pseudo-second order models (0.994, 0.995, 0.993, and 0.983). The experimental results obtained for Pb(II) were consistent with those for Ni(II). For Pb(II), the correlation coefficients (
R2 values) obtained for the pseudo-first order models for PBC, H-PBC, E-PBC, and O-PBC were 0.883, 0.822, 0.880, and 0.867 (
Table 5), respectively, which were lower than those for the pseudo-second order models (0.989, 0.987, 0.993, and 0.993). Furthermore, the values calculated for the pseudo second order model of
qe were consistent with the experimental q
e,exp values based on Ni(II) and Pb(II). The pseudo-second order kinetic constants (
K2 values) were small for Ni(II) and Pb(II), thereby indicating that the adsorption rate decreased as the contact time increased and it was proportional to the amount of adsorption sites [
4]. These results suggest that the adsorption process was controlled by chemical adsorption rather than physical adsorption.
3.5. Adsorption Isotherms
The Langmuir and Freundlich models were used to simulate the adsorption equilibria for heavy metals by the adsorbents. The fitted models and experimental data are shown in
Figure 4, and the calculated parameters are presented in
Table 6 and
Table 7. Langmuir model obtained a high coefficient of determination (
R2 ≥ 0.99) and it was more suitable than the Freundlich model, where it suggested the homogeneous adsorption of Ni(II) and Pb(II). The maximum monolayer adsorption capacities for H-PBC were determined as 64.94 mg/g for Ni(II) and 243.90 mg/g for Pb(II). According to the Freundlich isotherm, the value of 1/
n < 1 indicated the favorable adsorption of heavy metals. Moreover, the dimensionless separation factor R
L was employed to check the favorability of the adsorption process. Based on
RL = 1/(1 + K
LC
0) [
29], the Ni(II) and Pb(II) adsorption processes were thermodynamically favorable because the values of R
L ranged from 0 to 1.
The comparison of the maximum adsorption capacities of Ni(II) and Pb(II) using previously reported adsorbents are summarized in
Table 8 and
Table 9, showing that the adsorption capacity of H-PBC was slightly higher than other adsorbents, such as sawdust [
30,
31], and straw [
32,
33]. Different adsorbents exhibit different adsorption properties for heavy metals, which may be caused by the variances in raw material composition and preparation methods, resulting in the composition and pore structure of biochar. Consequently, H-PBC is assumed to be a promising adsorption material to remove heavy metals from water.
3.6. Adsorption Mechanism
The adsorption mechanisms for toxic metals on modified-biochar involve the comprehensive influence several types of reactions, such as ions exchange, precipitation, complexation, and physical adsorption [
40,
41]. For H-PBC, the larger S
BET value and the experimental data (
R2 > 0.87) fitted well by Freundlich model, indicated that physical adsorption was involved in the process. Moreover, FITR analysis demonstrated the presence of carboxyl and hydroxyl groups, which might facilitate the adsorption process via ion exchange and hydrogen bonding. Jiang et al., [
8] synthesized H
3PO
4-HC via a hydrothermal carbonization process and found that the oxygen-containing functional groups on H
3PO
4-HC, such as C=O groups, played the main role in Pb(II) adsorption. Additionally, the electrostatic interaction may also affect adsorption process, which could prove in adsorption process under different pH conditions. For E-PBC, the main adsorption mechanism may be complexation. Li et al., [
15] found that the formation of -NH
2Pb when a new peak was detected in N spectra after using EDTA-modified biochar to adsorb Pb(II) and attributed the complexation between -NH
2 and Pb(II). FITR analysis demonstrated the presence of -NH
2 in E-PBC. For O-PBC, FITR analysis proved the existed of -OH, C-O, which could enhance the adsorption process via precipitation and complexation.