3.1 Characterization of activated carbon of Jerivá coconut
Vegetable-derived adsorbents are widely used to remove dyes (Sellaoui et al. 2020) as they are materials that contain functional groups that can influence and favor the sorption of organic compounds, such as MB. Thus, the surface functional groups of the ACJC were identified by the FTIR spectrum (Fig. 1). In the region of 3800-2700 cm-1 the highest frequencies of the hydroxyl group (O-H) are found, so the stretch at 3421 cm-1 can be attributed to this group (Deng et al. 2010; de Costa et al. 2015; Gupta et al. 2018; dias Júnior et al. 2019). The O-H group, in plant materials (such as coconut fiber), is related to cellulose, hemicellulose and lignin (Sreekala et al. 1997; Setiabudi et al. 2016; Gupta et al. 2018). It is also possible to infer the presence of the amino group (NH2) with absorptions also in the range of 3200-3600 cm-1 (Yuvaraja et al. 2014), thus suggesting the overlapping of two bands in this region. It should be noted that the amino and hydroxyl groups have adsorptive characteristics (Tang et al. 2019), corroborating the result obtained by the N2 isotherms, which point to chemisorption (Ruthven 1984). The high intensity band around 1600 cm-1 can be attributed to the vibration of the C=C groups of the aromatic rings. This band can show the presence of lignin and also the formation of aromatic compounds originated by the elimination of hydrogen and oxygen from aliphatic compounds during pyrolysis, and can provide information regarding the stretching of the C=C bond constituting the aromatic benzene molecules or rings in lignin, typical of carbonaceous material (de Sales et al. 2015; dias Júnior et al. 2019). There may also be the presence of the nitro group, since it is found in the absorption peaks of 1600-1530 cm-1 and 1390-1300 cm-1. The appearance of peaks at 1178 and 1085 cm-1 may originate from the phosphorus functional groups, resulting from the use of H3PO4 as an activated agent (Pan et al. 2019). Bands with values <1000 cm-1 are related to constituents present in wood, such as cellulose and lignin, which are eliminated at the activation temperature (600°C) (Lammers et al. 2009; Özgenç et al. 2017). Chemical modification of the carbon structure by H3PO4 and calcination may have enhanced the adsorptive properties (Luo et al. 2018; Tang et al. 2019).
The SEM allowed the study of the structural and morphological characteristics of activated carbon, such as the existence of pores and chemical composition. In Fig. 2a, the micrographs at 70x magnification of the ACJC are shown, where it is observed that the surface has a very irregular, rough and heterogeneous topography when it comes to the size of the grains. At 5000x magnification (Fig. 2b) the presence of a large volume of well-defined pores can be seen. These results are similar to those found in other studies, which developed carbons activated with ZnCl2 from Jerivá coconut (Carvalho et al. 2019; Salomón et al. 2021). Pigatto et al. (2020) also synthesized an adsorbent material having the Jerivá coconut as precursor material, however, there was no chemical activation, only calcination, when analyzing the SEM, it is clear that the surface is rough and irregular, but it is not possible to observe the presence of pores so well structured. Thus, it can be inferred that calcination was responsible for evaporating/degrading the phosphoric acid that was impregnated in the Jerivá coconut, leaving the space that was previously occupied by the chemical reagent, resulting in the opening of pores (Deng et al. 2010; Enniya et al. 2018).
Through dispersive X-ray analysis (EDS) it was possible to identify the main elements present in an ACJC sample. It is observed that the impregnation of phosphoric acid in the adsorbent was efficient, with a significant P peak at 2 keV. It was also identified the presence of C and O at 0.3 and 0.5 keV, which can influence the sorption properties, increasing the surface polarity (Barrer 1966). The presence of C, P and O was also obtained in the study by Basu et al. (2018) which performed chemical activation with H3PO4. In addition to sodium, iron, potassium, calcium and tungsten in smaller amounts on the surface. SEM micrographs and EDS analysis can be viewed in Fig. SM1.
Figure 3 shows the nitrogen adsorption/desorption isotherm and the pore size distribution. According to the IUPAC, the isotherm can be classified as type I (Ruthven 1984; Thommes et al. 2015; Thi et al. 2020), since there is no hysteresis between the curves, indicating a monoatomic approximation to a limit in the adsorption capacity that corresponds to the formation of a complete layer. Usually, this type of isotherm is found in microporous adsorbates that act by chemisorption (Thi et al. 2020). The N2 isotherm of the Jerivá coconut adsorbent developed by Salomón et al. (2021) obtained a similar result to the ACJC.
Using the N2 isotherm, the surface area, pore volume, and average pore size were calculated (Table 1) and the values were compared to those reported in the literature using this biomass. The chemical activation with H3PO4 combined with high temperature was responsible for an increase of almost 700% in the surface area of the adsorbent when compared to the thermal treatment without the addition of an activating agent (Pigatto et al. 2020). The largest surface area and pore volume of the adsorbent activated with ZnCl2 and produced by Carvalho et al. (2019) can be attributed to the calcination conditions, since the firing temperature was 450°C, while in the present study 600°C was used and Salomón et al. (2021) 700°C, suggesting that lower carbonization temperatures are more indicated.
Table 1
Texture characterization of activated carbons produced from Jerivá coconut
Activating agent | SBET (m2 g−1) | TPV (cm3 g−1) | Average pore size (nm) | Reference |
H3PO4 | 750 | 0.442 | 2.36 | This work |
ZnCl2 | 1435 | 0.876 | - | Carvalho et al. (2019) |
ZnCl2 | 782 | 0.441 | 2.26 | Salomón et al. (2021) |
- | 108 | 0.090 | - | Pigatto et al. (2020) |
The relationship between the size of the adsorbate molecule with the pores can influence the diffusion mechanisms, in the case of macropores, size differences are favored, since the radius of the dye molecule is much smaller than the average pore size (Nascimento et al. 2014). In the case of micropores, they are often attributed as the main contributors to the greater sorption capacity of small molecules, such as MB, which has a molecular diameter of 0.8 nm. According to Barton (1987), a pore diameter of the order of 1.3 nm is necessary for effective sorption, so the results of the specific surface characterization allow us to infer that the adsorbent is capable of retaining the contaminant in the micropores, macropores and on the surface of the developed material.
3.2 Experimental design
Many works in the literature carry out an isolated evaluation of each variable in order to obtain the best conditions for the sorption process (Basu et al. 2018; Othman et al. 2018; Calimli et al. 2020), however, when one parameter is fixed and the others vary, it is not possible to evaluate the interaction of the variables. Thus, experimental planning allows optimizing processes by performing an interaction of all variables at the same time and obtaining results with statistical significance. The effects of the studied variables on MB removal using ACJC were evaluated using the Pareto graph (Fig. SM2). It was possible to observe that agitation was the only statistically significant variable (+4.7415), so that the increase in agitation favors dye sorption. The agitation is a factor that interferes directly in the liquid-solid sorption (Ruthven 1984), the higher the speed with which the reactor is agitated the higher the mass transfer rate and the lower the resistance that the MB molecules will find to move to the adsorbent surface, this behavior is expected since the mass transfer coefficient in the external film (kf) is the inverse of the resistance (Nascimento et al. 2014; Cremasco 2015). The other variables were not significant, being possible to reduce operating costs, using a lower adsorbent dosage, higher dye concentrations, the aqueous solution does not need pH adjustment and room temperature can be used. These conclusions are fundamental to improve the process to minimize resources and energy consumption, thus meeting the third and fourth principles of green engineering (Anastas and Zimmerman 2003).
3.3 Mechanism of MB adsorption kinetics on ACJC
The kinetic study was carried out in order to determine the mass transfer mechanisms that control MB sorption in the ACJC. Four kinetic models were evaluated, such as fractional order Avrami, Elovich, Pseudo-first-order and Pseudo-second-order, and the parameters obtained are described in Table 2. It can be seen that the Avrami-fractional-order model obtained good adjustments to the experimental data, presenting high R2 (> 0.990) and R2adj (> 0.988) and the lowest error values, SSE (< 0.385) and HYBRID (< 0.316), for the three concentrations studied, indicating that the Avrami model is able to adequately represent the experimental results. Thus, the calculated value of qe predicted by the equation is very close to the value obtained experimentally.
Avrami's kinetic model has a wide application in studies involving the adsorption of dyes on activated carbons (Alencar et al. 2012; Adebayo et al. 2014; Lima et al. 2016), precisely for being able to adapt to different systems, not fixing a single kinetic order (nAV) (Lima et al. 2015). As mentioned before, qe indicates the sorption capacity of MB in ACJC, while the parameter nAV is associated with possible potential changes in the adsorption mechanism (Liu and Liu 2008; Lima et al. 2016), allowing many kinetic orders to describe the adsorption process over time.
Table 2
Kinetic parameters for adsorption of MB on ACJC in different concentrations.
| 51.4 mg L−1 | 77.1 mg L−1 | 103.5 mg L−1 |
Experimental data |
qe (mg g−1) | 10.26 | 15.40 | 20.67 |
Avrami fractional-order |
qe (mg g−1) | 10.52 | 15.58 | 21.19 |
kAV (min−1) | 0.3316 | 0.4851 | 0.5497 |
nAV | 0.5490 | 0.5295 | 0.4231 |
R2 | 0.9898 | 0.9965 | 0.9953 |
R2adj | 0.9880 | 0.9958 | 0.9944 |
SSE | 0.3062 | 0.2583 | 0.3855 |
HYBIRD | 0.3166 | 0.0796 | 0.0797 |
Elovich |
α (mg (g−1 min−1)) | 15.81 | 121.60 | 188.60 |
β (mg g−1) | 0.6595 | 0.5482 | 0.4253 |
R2 | 0.9608 | 0.9541 | 0.9756 |
R2adj | 0.9536 | 0.9457 | 0.9712 |
SSE | 0.5989 | 0.9541 | 0.8753 |
HYBRID | 0.5211 | 0.9457 | 0.2671 |
Pseudo-first-order |
qe (mg g−1) | 9.88 | 14.86 | 19.41 |
k1 (min−1) | 0.1538 | 0.2622 | 0.2830 |
R2 | 0.9481 | 0.9514 | 0.9141 |
R2adj | 0.9387 | 0.9426 | 0.8985 |
SSE | 0.7925 | 1.0510 | 1.7776 |
HYBRID | 3.3843 | 1.8985 | 2.3723 |
Pseudo-second-order |
qe (mg g−1) | 10.61 | 15.71 | 20.54 |
k2 (g mg−1 min−1) | 0.0232 | 0.0293 | 0.0242 |
R2 | 0.9735 | 0.9850 | 0.9705 |
R2adj | 0.9686 | 0.9823 | 0.9652 |
SSE | 0.5159 | 0.5423 | 0.9867 |
HYBRID | 1.4498 | 0.5028 | 0.7728 |
Through the results obtained, an increase in the rate (kAV) of dye sorption to the surface of the adsorbent is observed as the concentration increases, ranging from 0.3316 min−1 to 0.5497 min−1. This occurred because the driving force is directly related to the increase in concentration, providing an increase in the mass transfer potential (Feng et al. 2011). In Fig. 4, the kinetic results obtained experimentally and the adjustment are presented. It is verified that the sorption was extremely fast and the amount of MB adsorbed increases in higher concentrations of the solute, however, when there is no more significant variation in the mass transfer, the equilibrium of the process is reached (Thi et al. 2020). At the equilibrium time of 80 min, residual concentrations below 0.4 mg L−1 were obtained, which represents a removal greater than 99 % for the three concentrations studied. This result is consistent with other highly efficient adsorbents, such as Eucalyptus residue activated from H3PO4, which achieved equilibrium in 60 min (Han et al. 2020), so the equilibrium time can be considered fast compared to other adsorbents which took about 120 min (Dawood et al. 2017), 6 h (Ahmad et al. 2021). The treated water, after adsorption with ACJC, has a very high-quality level, especially when compared to separation processes by membranes (Parakala et al. 2019; Hou et al. 2020), which are highly recommended technologies for water reuse.
The study of the mass transference of the dye molecules from the fluid phase to the solid phase of the adsorbent is extremely relevant and involves the diffusion of external and internal mass and chemical reactions on the surface of the adsorbent (Prajapati and Mondal 2020). External diffusion was not evaluated because, at agitation above 150 rpm, the thickness of the film boundary layer decreases so that it can be disregarded (Wang et al. 2018). Thus, the experimental data were fitted to the intraparticle diffusion model, such as the Weber and Morris model, which assumes that the removal of adsorbate varies with the square root of time, being determinant for the velocity (Weber and Morris 1963).
According to the result shown in Figure 5, the adsorption process can be divided into three steps, following the same behavior for the three concentrations. The first zone, obtained in the first 10 minutes, corresponds to external mass transfer. A possible explanation for the higher speed (observe kid for the first zone in Table SM3) and adsorption rate may be due to the presence of available sites and some surface functional groups (Feng et al. 2012), such as OH and NH2, revealed in the characterization of the adsorbent through the FTIR. Furthermore, it is known that the external mass transfer can be affected by concentration and agitation, so that increasing the dye concentration was able to accelerate the diffusion to the solid surface (Nascimento et al. 2014). Other factors that may have contributed to the increase in adsorption speed are the pH of the solution, the particle size, and the pore size distribution.
The steeper curve suggests that the adsorption is instantaneous and that there is a strong attraction between the dye and the functional groups. The second zone was established at 40 min and according to Feng et al. (2012) can be attributed to the intraparticle diffusion process. At this stage, there is still the presence of sites available for dye sorption, although many of which have already been filled in the previous zone. The third and last zone corresponds to the equilibrium stage, where the sorption rate is much lower (observe kid for the third zone in Table SM3), until it finally ceases, due to the low concentration of adsorbate in the liquid phase (Nascimento et al. 2014; Wang et al. 2018).
The parameters of the intraparticle diffusion model are shown in Table SM3, and it is possible to observe that the model had a good adjustment to the experimental data, with high R2 (> 0.937) and R2adj (> 0.905) and low errors (SSE < 0.307 and HYBRID < 1.514). The kid constant (mg g−1 min-0.5) represents the rate of diffusion in the pores, while C (mg g−1) is related to the thickness of the boundary layer or the resistance to mass transfer in the external liquid film. From the values found for the diffusion speed, it can be seen that they are high, especially in the first and second zone, when compared to other similar works (Kavitha and Namasivayam 2007; Karagöz et al. 2008).
The values of C increased with the increase in the initial concentrations of MB, which indicates the increase in the thickness of the boundary layer at the highest dye concentrations, while the kid had the same behavior. This is because the adsorption driving force is more intense at higher concentrations (Weng et al. 2009; Feng et al. 2011). It is also possible to observe that the diffusion coefficient decreases with time, which means that the speed at which the pores are filled by the adsorbate decreases until equilibrium is reached. The value of constant C for all sorption steps was different from zero, as well as the regression curves, for both initial dye concentrations (Figure SM3), even if they were extended they will not pass through the origin. This means that some other mechanisms such as external surface sorption, adsorption equilibrium and intraparticle diffusion are involved in the process (Özcan et al. 2005; Tang et al. 2019).
3.4 Thermodynamic equilibrium study
The applicability of new adsorbents, such as Jerivá coconut, in adsorption processes depends on information related to thermodynamic equilibrium (Kebede et al. 2018), since it plays a fundamental role in the mass transfer of the contaminant from the liquid to the solid phase (Cremasco 2015), in addition to showing whether the adsorption is feasible. The study was carried out with the objective of obtaining the adsorption capacity of ACJC at temperatures of 15, 25, 35 and 45 ± 1°C and through adjustments of the BET, Langmuir, Freundlich and Sips models, the adsorption thermodynamics was clarified. Table 3 shows the parameters of each model and the adjustment coefficients. It appears that the Sips model has higher R2 (> 0.9511) and R2adj (> 0.9316), however the other models also obtained good fits, thus, through the error analysis between the values predicted by the model and the experimental data, note that the Sips model has lower results for SSE and HYBID, indicating that it is the model that best represents the adsorption of MB in ACJC.
Table 3
Adsorption equilibrium parameters for MB on ACJC performed in different temperatures (adsorbent dosage = 5 g L−1, shaking = 150 rpm, pHi = 5.85, pHf = 2.44).
| 15°C | 25°C | 35°C | 45°C |
BET |
qm (mg g−1) | 60.03 | 70.18 | 85.36 | 95.57 |
KS (L mg−1) | 0.3449 | 0.8910 | 1.2862 | 1.7985 |
KL (L mg−1) x10−4 | 4.2587 | 3.1024 | 6.4955 | 8.0202 |
R2 | 0.9749 | 0.9468 | 0.9871 | 0.9592 |
R2adj | 0.9649 | 0.9255 | 0.9819 | 0.9429 |
SSE | 4.088 | 7.527 | 4.327 | 8.739 |
HYBRID | 0.566 | 3.305 | 0.964 | 2.586 |
Freundlich |
KF (mg g−1(mg L−1)−1/n) | 30.86 | 36.49 | 48.65 | 53.03 |
nF | 7.253 | 7.970 | 7.42 | 6.694 |
R2 | 0.9801 | 0.9301 | 0.9775 | 0.9824 |
R2adj | 0.9768 | 0.9185 | 0.9738 | 0.9795 |
SSE | 3.221 | 8.597 | 5.110 | 5.346 |
HYBRID | 0.310 | 1.884 | 1.335 | 0.664 |
Langmuir |
qm (mg g−1) | 66.35 | 79.79 | 94.22 | 106.20 |
KS (L mg−1) | 0.2537 | 0.7211 | 0.9998 | 1.1180 |
R2 | 0.9553 | 0.9400 | 0.9531 | 0.8918 |
R2adj | 0.9479 | 0.9300 | 0.9453 | 0.8738 |
SSE | 4.834 | 7.951 | 7.383 | 13.575 |
HYBRID | 0.820 | 3.072 | 0.746 | 4.724 |
Sips |
qm (mg g−1) | 93.69 | 115.08 | 141.50 | 254.40 |
KS (mg L−1)−1/n | 0.365 | 0.521 | 0.492 | 0.258 |
nS | 2.897 | 4.105 | 3.416 | 4.463 |
R2 | 0.9856 | 0.9511 | 0.9827 | 0.9845 |
R2adj | 0.9799 | 0.9316 | 0.9758 | 0.9784 |
SSE | 2.583 | 5.849 | 4.225 | 4.729 |
HYBRID | 0.251 | 1.165 | 0.890 | 0.464 |
In Figure 6, it can be seen that the curves were very close to the y-axis, indicating the process is extremely favorable (Ruthven 1984). The qmax increases from 93.69 mg g−1 at 15°C to 254.40 mg g−1 when the temperature is 45°C, indicating that the process is endothermic, that is, there is an increase in adsorption capacity with increasing temperature. This behavior is expected, since the increase in temperature provides higher rates of diffusion of dye molecules through the boundary layer, due to a decrease in the viscosity of the medium and an increase in the mobility of the adsorbate. Furthermore, this increase in adsorption capacity can be attributed to the effects of swelling on the ACJC itself (Hameed and Ahmad 2009). The Sips model is empirically based and consists of a combination of Langmuir and Freundlich isotherms, thus, when the exponent nS = 1 then the model becomes similar to Langmuir, when the constant KS is close to 0, the model approaches Freundlich (Lima et al. 2015; An et al. 2017). According to Table 2, it was observed that nS > 1 and KS ≈ 0, indicating that adsorption occurs in multilayers on a heterogeneous surface.
In order to evaluate the adsorption efficiency of the ACJC MB, a survey of the evaluated dye concentration range (MB) and maximum adsorption capacity (qmax) was carried out, based on different adsorbents obtained from biomass, as shown in the Table 4. Initially, it is noticed that ACJC can be considered an adsorbent with great potential for the adsorption of dyes, such as MB, with a qmax higher than most of the investigated adsorbents. It is observed that adsorbents that had some type of surface modification by a chemical agent had higher qmax.
Table 4
Maximum adsorption capacities for MB in different biomasses as adsorbents
Adsorbents | Activator | MB (mg L−1) | qmax (mg g−1) | References |
seeds of Cedrela odorata L | - | 10-750 | 158.5 | Subratti et al. (2021) |
Eucalyptus sheathiana bark | - | 10-100 | 104.2 | Dawood et al. (2016) |
seeds of Cedrela odorata | - | 20-200 | 111.8 | Babalola et al. (2016) |
Pinus radiata | - | 10-70 | 106.4 | Dawood et al. (2017) |
Corn cob | KOH | 100-500 | 333.0 | Medhat et al. (2021) |
Peanut shell | NaOH | - | 555.6 | Ahmad et al. (2021) |
Coconut shells | ZnCl2 | 25-500 | 156.2 | Yağmur and Kaya (2021) |
Eucalyptus residue | H3PO4 | - | 977.0 | Han et al. (2020) |
Syagrus romanzoffiana | H3PO4 | 175-915 | 254.4 | This study |
Based on the equilibrium concentrations of the adsorbate and the sorption capacity of the adsorbent at different temperatures, the thermodynamic parameters were obtained, as shown in the Van’t Hoff graph (Fig. SM3). Thermodynamic analysis is important to understand how a contaminant is transferred from one phase to another and how it is distributed (Mahmoodi et al. 2010). The enthalpy change (ΔH0) was 29.99 J mol−1, the positive sign of ΔH0 indicates that the process is endothermic (Hasan et al. 2019), and the isotherms (Table 3) show that the adsorption capacity is favored by increasing temperature. Regarding the entropy variation (ΔS0), 86.39 J mol−1 K−1 was obtained, indicating the affinity between ACJC and MB. The behavior regarding the energy flow between the system and the neighborhood, through the Gibbs free energy (ΔG0), also indicates spontaneity in the opposite direction, since the values (-24.86 kJ mol−1, -25.73 kJ mol−1, -26.60 kJ mol−1 and -27.46 kJ mol−1) are negative (Mahmoodi et al. 2010). Similar thermodynamic results were found for other adsorbents developed from Syagrus romanzoffiana activated with ZnCl2 (Carvalho et al. 2019; Salomón et al. 2021).
3.5 MB Adsorption Mechanisms
The pH of the solution is one of the most important factors that regulate the sorption characteristics of active carbons, as the pH can change their surface charge and solute speciation (Wang and Chen 2015; Tang et al. 2019). Thus, this parameter was monitored during the kinetic and thermodynamic equilibrium tests. It was observed that the initial pH of the MB solution was 5.85 ± 0.20 and after contact with ACJC the pH decreased to 2.44 ± 0.32, this behavior being expected since the activation of the adsorbent occurred by impregnation of phosphoric acid. The activating agent was responsible for the acidic surface of the carbon and, even after washing with water and calcium carbonate, this characteristic did not change. As the mass transfer from the liquid medium to the solid occurred under acidic conditions and the dye pKa is 3.14, it is suggested that the sorption occurred by Van der Waals interactions. Adsorbate is in protonated form when pH < pKa and deprotonated when pH > pKa (Liu et al. 2010), thus at the beginning of the sorption process, MB is presented in deprotonated form, but over time the pH decays and stabilizes at 2.44 ± 0.32 indicating that the dye is in a protonated condition in this study.
The adsorption of a certain molecule on the surface of the adsorbent can occur through different types of interactions that affect the sorption capacity (Somsesta et al. 2020). In Fig. 7, an illustration of the possible mechanisms involved has been sketched. The first mechanism that may be acting on the adsorption of MB on ACJC is the π-π stacking interaction, since the characterization of ACJC revealed the presence of functional groups characteristic of aromatic rings and the MB molecule itself is composed of these rings, suggesting a possible multilayer adsorption (Liu et al. 2011; Iwuozor et al. 2021). Adsorbents obtained from biomass have surface groups of the OH type, which can interact with the dye molecule, more specifically with nitrogen, forming a hydrogen bond (Liu et al. 2011; Somsesta et al. 2020; Iwuozor et al. 2021). Electrostatic attractions can also contribute to adsorption, since the ACJC surface has functional oxygenated groups (with more acidic and negatively charged characteristics) and MB is a cationic dye (Somsesta et al. 2020; Jawad et al. 2021).