Next Article in Journal
Modification of DSSC Based on Polymer Composite Gel Electrolyte with Copper Oxide Nanochain by Shape Effect
Next Article in Special Issue
Facile Conversion of Polystyrene Waste into an Efficient Sorbent for Water Purification
Previous Article in Journal
Research on Silicone Rubber Sheds of Decay-Like Fractured Composite Insulators Based on Hardness, Hydrophobicity, NMR, and FTIR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Hexavalent Chromium Removal with Hydrophobically Modified Cellulose Nanofibers

by
Francisco de Borja Ojembarrena
1,*,
Jose Luis Sánchez-Salvador
1,
Sergio Mateo
1,
Ana Balea
1,
Angeles Blanco
1,
Noemí Merayo
2 and
Carlos Negro
1,*
1
Department of Chemical Engineering and Materials, Complutense University of Madrid, Avda. Complutense s/n, 28040 Madrid, Spain
2
Department of Mechanical, Chemical and Industrial Design Engineering, ETSIDI, Universidad Politécnica de Madrid, Ronda de Valencia 3, 28012 Madrid, Spain
*
Authors to whom correspondence should be addressed.
Polymers 2022, 14(16), 3425; https://doi.org/10.3390/polym14163425
Submission received: 29 July 2022 / Revised: 15 August 2022 / Accepted: 16 August 2022 / Published: 22 August 2022
(This article belongs to the Special Issue Polymeric Materials for Wastewater Treatment Applications)

Abstract

:
Cellulose nanofibers (CNF) are sustainable nanomaterials, obtained by the mechanical disintegration of cellulose, whose properties make them an interesting adsorbent material due to their high specific area and active groups. CNF are easily functionalized to optimize the performance for different uses. The hypothesis of this work is that hydrophobization can be used to improve their ability as adsorbents. Therefore, hydrophobic CNF was applied to adsorb hexavalent chromium from wastewater. CNF was synthetized by TEMPO-mediated oxidation, followed by mechanical disintegration. Hydrophobization was performed using methyl trimetoxysilane (MTMS) as a hydrophobic coating agent. The adsorption treatment of hexavalent chromium with hydrophobic CNF was optimized by studying the influence of contact time, MTMS dosage (0–3 mmol·g−1 CNF), initial pH of the wastewater (3–9), initial chromium concentration (0.10–50 mg·L−1), and adsorbent dosage (250–1000 mg CNF·L−1). Furthermore, the corresponding adsorption mechanism was identified. Complete adsorption of hexavalent chromium was achieved with CNF hydrophobized with 1.5 mmol MTMS·g−1 CNF with the faster adsorption kinetic, which proved the initial hypothesis that hydrophobic CNF improves the adsorption capacity of hydrophilic CNF. The optimal adsorption conditions were pH 3 and the adsorbent dosage was over 500 mg·L−1. The maximum removal was found for the initial concentrations of hexavalent chromium below 1 mg·L−1 and a maximum adsorption capacity of 70.38 mg·g−1 was achieved. The kinetic study revealed that pseudo-second order kinetics was the best fitting model at a low concentration while the intraparticle diffusion model fit better for higher concentrations, describing a multi-step mechanism of hexavalent chromium onto the adsorbent surface. The Freundlich isotherm was the best adjustment model.

1. Introduction

Chromium is a harmful water pollutant. Several environmental effects have been associated with the presence of trivalent and hexavalent chromium in natural water bodies. Trivalent chromium shows lower toxicity due to its low cell permeability, whereas hexavalent chromium presents acute toxicity to many species, causing carcinogenicity and mutagenicity [1,2,3], aside from being neurotoxic [4]. Furthermore, Cr(VI) is highly toxic and one of the most common environmental contaminants and is non-easily biodegradable in nature, thus staying in the environment for a long time, polluting the soil and water, with the subsequent health risks to humans and wildlife [5]. Several industries generate hexavalent chromium in their effluents including tanneries, electroplating, hardware, textile, cement, or mining industries [3,6]. The limit of discharge of this contaminant varies depending on the country, but most of European countries have low limits in water, between 0.05 and 0.1 mg·L−1 [7].
The most relevant hexavalent chromium removal treatments include adsorption, ion exchange, membrane filtration, electrochemical treatments, coagulation–flocculation, chemical precipitation, biological treatments, and photocatalysis [6,8,9,10]. Among them, adsorption is the most widely used technology due to the low operational costs, ease of operation, low environmental impact associated with the process, and its effectiveness in removing trace levels of chromium [11]. Some of these adsorption treatments are based on microbial bioadsorption [12], since some microbes are able to reduce hexavalent chromium to trivalent chromium; as the latter is much less harmful, their excretion by microbes is not an issue [5]. Common adsorbents such as activated carbons or zeolites have been successfully applied, but there is a need for new materials from renewable sources, with better biodegradability and whose raw materials have less environmental impact. Adsorbents coming from renewable sources include those obtained from waste phytomass, which could be valorized to adsorb hexavalent chromium, achieving 87.2% removal, with the possibility of regeneration for at least four cycle runs [13]. Agro-based biomass has also been used as a bioadsorbent, which contributes to significant sustainable waste management, with an adsorption capacity of more than 10 mg/g [12]. Other waste materials used to produce adsorbents are fishbone waste to produce hydroxyapatite, which successfully removed Ni2+, Cu2+, and Zn2+ (more than 95% removal at the optimum conditions) [14]. Furthermore, nanomaterials can be used as adsorbents such as the polyacrylonitrile nanofiber membrane modified with bovine serum albumin used to remove Ca2+ from process streams, achieving removal efficiencies of about 62% [15]. Other membrane structure materials used as adsorbents are, for example, cellulose acetate-based membranes with glass nanoparticles for CO2 separation [16]. The surfaces of the nanomaterial adsorbents are usually functionalized or chemically modified to improve Cr(VI) adsorption such as the attachment of active binding sites [17].
As cellulose is the most abundant biopolymer on Earth, and it is an inexpensive and renewable material, it is an excellent option as a raw material to produce new adsorbents [18,19] as well as its derivates such as sulfate cellulose [20].
Nanocelluloses are cellulosic materials with at least one dimension in the nanometer scale. There are different types of nanocelluloses including bacterial cellulose (BC), cellulose nanocrystals (CNC), and cellulose nanofibers (CNF). These materials show good adsorptive properties due to their high specific area and the presence of active groups such as carboxyl groups over their surface [21]. These nanomaterials have been proven as excellent adsorbents of different heavy metals such as cadmium, nickel, copper, lead, arsenic, iron, silver, cobalt, or mercury, thanks to their ease of functionalization and application, which enhances the adaptation to each heavy metal [22,23]. Celluloses are commonly anionic charged, caused by the pulping and bleaching processes [24]. Most of the heavy metals are cationic species, facilitating the direct application of nanocelluloses without any surface modification to obtain both a high removal yield and high adsorption capacity as well as fast adsorption rates. Khoo et al. [25] confirmed that there was a large affinity of CNC for Cu2+, Cd2+, Ni2+, Pb2+, Zn2+, Fe3+, and Ag+ adsorption, reached with none or low modifications. Compared to these cationic heavy metals, the adsorption capacity of CNC for anionic As(V) was reduced by an order of magnitude. The same trend was also observed by Liu et al. [26], who applied both untreated CNF and 2,2,6,6-tetramethylpiperidin-1-yl-oxyl (TEMPO)-mediated oxidized CNF to a list of 15 heavy metals, with hexavalent chromium the only anionic species. The authors found a maximum adsorption capacity of 87.5 mg·g−1 using TEMPO-oxidized CNF-PAN membranes to adsorb Cr(VI), which was lower than the adsorption capacity reached when cationic Pb(II) was adsorbed.
The adsorption of anionic hexavalent chromium species with cellulosic materials usually requires other strategies and one of the most common options is surface modification. These surface modifications often include a chemical reaction, which forces a change in the surface groups present in the cellulose. Cationization, oxidation, esterification, alkaline treatment, or halogenation are some of the most typical treatments [25]. These modifications have demonstrated their applicability, and some of them such as dialdehyde oxidation or the cationization of celluloses have reached the large yields and adsorption capacities of hexavalent chromium [27,28]. Nevertheless, they usually imply a purification process that needs washing steps and the loss of reagents after treatment.
Another option could be the coating of the nanocelluloses. In this case, the nanocelluloses undergo a reaction where the coating reagent is completely mixed with the nanocellulose in suspension and the nanocellulose surface becomes covered by the coating agent. As this chemical is directly added to the final nanocellulose suspension, there is no need for separation after the reaction, avoiding the loss of material. Among the possible hydrophobization reactions, silanization becomes a great option because of the maintenance of large specific area of the fibrils and the enhanced hydrophobicity of cellulose fibers [29]. This fact allows for subsequent advanced treatments such as lyophilization to obtain aerogels [30]. Other silanized materials such as graphene oxide have been successfully applied to hexavalent chromium adsorption [31,32], but little information has been developed on the application of silanized cellulose-based materials to this contaminant. The closest approach is by Jamroz et al. [33], who indicated that amino-silanized celluloses could become an option for the ultra-trace determination of hexavalent chromium in water, which is not a direct wastewater treatment.
The state-of-the-art shows that there is a lack of studies about the use of nanocellulose to remove Cr(VI) from water streams, and particularly, the efficiency of hydrophobic nanocelluloses is not known. Furthermore, the kinetic and adsorptive behavior of these modified materials have not been well-studied yet. On the other hand, in this study, the applied hydrophobization method was based on a coating process to reduce product losses.
Therefore, in this study, the main objective was the synthesis and characterization of CNF hydrogels modified with methyl trimetoxysilane (MTMS) as a hydrophobic coating agent and their application to Cr(VI) adsorption from wastewater. The novelty of this approach is related to the application of hydrophobic CNF, hydrophobized by a coating method, to improve the adsorption performance of cellulosic nanofibers. Furthermore, the proposed adsorbent material application was optimized for the removal of Cr(VI) by performing a set of experiments on a batch system to evaluate the effect of system variables, MTMS dosage, pH value, initial chromium concentration, adsorbent dosage, and contact time. The kinetic and isotherm data were modeled to increase the knowledge related to the process, making easier its future application. Furthermore, different parameters have been obtained from a wide variety of models, which are the key factors in understanding the way hexavalent chromium interacts with the surface of the tested nanocelluloses.

2. Materials and Methods

2.1. Materials

Bleached pine kraft pulp was supplied by Arauco (Chile) and used as a raw material to synthetize CNF hydrogels. All of the experiments were performed using Milli-Q quality water. The chemicals used for the production and characterization of nanocelluloses were hydroxylamine hydrochloride, silver nitrate, methyl trimetoxysilane (MTMS), and 2,2,6,6-tetramethylpiperidin-1-yl-oxyl (TEMPO) supplied by Sigma Aldrich; sodium bromide, potassium chromate, sodium chloride, sodium hydroxide pellets, and sodium hypochlorite solution (10 w/v%) supplied by Panreac; and hydrochloric acid (37% v/v) and sulfuric acid (98% v/v) were supplied by Labkem. All of the chemicals used were of analytical grade.
To analyze hexavalent chromium, standard solutions were used. Hexavalent chromium standard solution (50 mg·L−1 as Cr(VI)) was used as a calibration standard for the spectrophotometric method. Analytical reagents for hexavalent chromium determination were purchased from Macherey Nagel (Dueren, Germany) following Standard Method 3500 Cr B. Poly(diallyldimethylammonium chloride) (PDADMAC) and polyethylenesulfonate (PesNA) solutions with a concentration of 0.00025 N were used as standard titration reagents during cationic and anionic demand determination.

2.2. CNF Hydrogel Synthesis

Briefly, bleached pine kraft pulp was disintegrated at 30,000 revolutions and diluted at 1% of consistency. Then, TEMPO-mediated oxidation was carried out through the addition of 5 mmol of NaClO per gram of cellulose using TEMPO and NaBr at 0.1 mmol and 1 mmol per gram of cellulose, respectively [34,35,36,37]. The oxidized cellulose suspension, after a cleaning process, was then treated under three-step mechanical homogenization at 600 bar pressure. Finally, the homogenized suspension underwent surface modification with MTMS to perform the hydrophobization, as indicated by Zhang et al. [30].

2.3. CNF Characterization

The characterization of the CNF hydrogels differs from the typical physical-chemical characterization of adsorbents as it is a hydrogel and must be characterized as a nanocellulose suspension instead of the common solid characterization. This characterization included the determination of the consistency of the suspension as well as the anionic demand and transmittance of the CNF suspension at λ = 800 nm, as described by Balea et al. [38], and the zeta potential of cellulose nanofibers and the amount of carboxyl groups in the oxidized pulp suspension determined by conductimetric titration as described by Sanchez-Salvador et al. [39].

2.4. Experimental Procedure: Batch Adsorption of Hexavalent Chromium Solution

The experimental installation consists of a batch, stirred laboratory beakers of 250 mL filled with 100 mL of sample of synthetic wastewater placed on multiple position hot plates under temperature control. Different operating conditions were studied: contact time, pH, chromium concentration, and adsorbent dosage. In the case of the CNF hydrogel, the dosage of hydrophobic reagent MTMS in the synthesis route was also optimized regarding adsorption efficiency.
During the kinetic experiments, the adsorption contact times were considered, taking samples at 1, 5, 15, 30, 60 min, 2, 3, 4, 24, and 48 h. The samples were filtered by cellulose acetate syringe filters (pore size 0.45 μm) to separate the insoluble fraction of the adsorbent from the soluble fraction as this material did not interact or adsorb chromium while passing through. Each experiment was considered finished (equilibrium time) when the pollutant level remained constant in 2–3 samples in a row. The hexavalent chromium concentration was determined by using a spectrophotometer calibrated with diluted solutions from a standard of chromium of 50 mg·L−1 between 0.02 and 0.50 mg·L−1 and the measurement was performed at the peak wavelength λ = 540 nm, as indicated in Standard Method 3500 Cr B [40]. Synthetic hexavalent chromium solution samples were diluted when the concentration exceeded the top of the calibration curve.
The pH effect was evaluated from acidic (pH 3) to neutral and alkaline (pH 9) conditions through the addition of HCl and NaOH, both at 0.1 mol·L−1. Hexavalent chromium concentration in the kinetic and isotherm studies varied from 0.1 to 50 mg·L−1. Hydrophobized CNF hydrogels were tested in a dosage range 250–1000 mg CNF·L−1. The dosage of the hydrophobic reagent MTMS during CNF hydrogel synthesis was varied from 0 to 5 mmol MTMS·g−1 CNF hydrogel, according to the indications from Zhang et al. [30]. Each test was repeated three times.

2.5. Isotherm and Kinetic Studies

The conversion of the hexavalent chromium concentration into adsorption capacity values was determined. Efficiency of the Cr(VI) removal and adsorption capacity was calculated according to Equations (1) and (2) [41].
%   Cr VI   removal = C 0 C t C o · 100  
q Cr VI mg   Cr VI adsorbed g   CNF = C 0 C t m ads · V water
where C0 and Ct represent the initial and equilibrium concentrations of pollutants in solution, respectively (mg·L−1); Vwater is the volume (L) of the solution; mads is the adsorbent mass (g CNF). The adsorbent mass (mads) can be calculated by the following Equation (3), which relates to the volume dose of the adsorbent and the consistency measured as indicated before.
m ads   g   CNF =   V ads · Consistency  
where Vads is the volume of adsorbent and the consistency is the dried mass of adsorbent at 60 °C [42]. These compiled data were subsequently analyzed through different kinetic models.
The experimental kinetic and isotherm data were fitted through different kinetic and isotherm models to identify the adsorption mechanisms of the hexavalent chromium onto the surface of the CNF. The selected kinetic and isotherm equations and their linearized forms are summarized in Table 1 and Table 2.

3. Results and Discussion

3.1. CNF Characterization

The cellulose-based physical-chemical characterization begins with the determination of the carboxyl group content on the oxidized cellulose before homogenization. Carboxylic acids created from primary hydroxyl groups on the cellulose leads to an increase in the repulsion between individual fibrils that form the cellulose fiber, which subsequently favors the homogenization process. Therefore, the number of carboxylic groups formed is indicative of the technical feasibility of the homogenization process, particularly to avoid clogging the device. The oxidized cellulose achieved 1.099 mmol COOH·g−1. This content was higher than the values reported by Patiño-Masó et al. [48] (0.75 mmol COOH·g−1), Lu et al. [49] (0.73 mmol COOH·g−1) and Balea et al. [38] (0.50 and 0.25 mmol COOH·g−1) using bleached kraft eucalyptus pulp, bleached bagasse, recycled newspaper, and corrugated container as raw materials under equivalent NaClO dosages, respectively. These differences are related to the quantity of impurities and lignin present in the source of cellulose, which has a strong influence on the amount of carboxylic groups formed for a specific amount of NaClO [38]. Therefore, the high number of carboxylic groups that formed was due to the reduced amount of impurities and lignin in the bleached pine kraft pulp used in this case, which agrees with the good performance of the oxidized cellulose during the homogenization process to obtain the CNF, and explains the low number of homogenization cycles (three).
Once the CNF hydrogel was obtained by homogenization, it was characterized by measuring the zeta potential of the CNF suspension. Zeta potential is indicative of the stability and separation of individual cellulose nanofibers, negatively charged; therefore, as its value becomes more negative, the suspension has more stability due to the repulsion between nanofibers that avoid aggregation between them. The value of the zeta potential for this suspension was −41.33 mV, which revealed the negative charge on the nanocellulose surface and the stability of the CNF suspension. In the literature, zeta potential values below −30 mV are indicators of bilateral repulsion and colloidal stability [50,51]. Other authors have obtained similar results for nanocellulose from acid hydrolysis (−38.2 mV) and TEMPO-oxidized nanocellulose (−46.5 mV) caused by the sulfonate and carboxyl groups, respectively [52].
Aside from the zeta potential, the cationic demand of the CNF hydrogel was measured, resulting in 838.5 µeq·g−1. This value was higher than the one obtained by Balea et al. [38] for recycled fibers (200–600 µeq·g−1) and lower than those obtained by Patiño-Masó et al. [48] (1000 µeq·g−1) for bleached kraft pulp from eucalyptus under the same testing conditions. Cationic demand is indicative of the degree of defibrillation achieved by homogenization. As the cellulose surface is negatively charged, a higher amount of negatively charged groups in the suspension means that the specific surface of the cellulose material/nanomaterial was higher. Therefore, the cationic demand will be higher as the specific surface increases, which implies that a higher number of individual nanofibers was achieved during the homogenization process. The degree of nanofibrillation achieved was directly related to both the cellulose source and the cycles applied during the homogenization, which explains the higher value achieved regarding recycled fibers such as the one obtained with virgin fibers, as the amount of cellulose in virgin sources is higher than the cellulose in recycled ones.
To evaluate the amount of nanofibers in the CNF suspension, the degree of nanofibrillation in the CNF, the transmittance at 800 nm was determined [53]. The value of transmittance achieved was 95.6%, indicating that the synthetized CNF were highly nanofibrillated. As a comparison, Patiño-Masó et al. [48] obtained a 99.13% of nanofibrillation yield with a suspension of 88.0% of transmittance. This fact shows that most of the CNF can be considered disaggregated and separated in small individual nanofibers through this treatment.

3.2. Kinetics of Cr(VI) Adsorption with Hydrophobic CNF

3.2.1. Effect of MTMS Dosage

The effect of MTMS dosage, added to coat the surface of the CNF, on the hexavalent chromium adsorption was evaluated. Three MTMS dosages were tested (0, 1.5, and 3 mmol MTMS·g−1 cellulose). These adsorption batch experiments were evaluated at pH 3, 1000 mg·L−1 of adsorbent dosage, and 0.1 mg·L−1 of initial chromium concentration. Hexavalent chromium concentration in the water was measured along time, using the three different hydrophobic CNF used as adsorbents (Figure 1). The trends of chromium decreased for each dose of MTMS, indicating a variation in both the adsorption rate and the maximum adsorption capacity. Whereas the CNF hydrogel without the hydrophobic coating reached a maximum chromium removal of about 80%, after 75 h of contact time, both the MTMS doped CNF hydrogels showed a higher adsorption capacity as they adsorbed hexavalent chromium up to complete abatement, after 75 h of contact time. However, the lowest MTMS dosage applied to the CNF enhanced the adsorption rate (more than 90% of chromium was removed before 6 h of treatment) in comparison with the CNF without the MTMS coating and the CNF with the highest MTMS dosage. Other modified CNF (i.e., cationized [54], carboxylated [55,56], acid treated [57], and diethylenetriamine [58] modifications) achieved similar removal yields, over 90% of hexavalent chromium removal after 120 min of contact time using doses of the adsorbent between 0.3 and 3 g adsorbent·L−1 at acidic pH values from 1 to 5.5.
Following the calculation methodology described in Section 2.5, the adsorption capacity was determined and plotted in each experiment. Then, these data were fitted according to the different kinetic equations proposed in Table 1 (Figure 2).
The best fitting model in this case was the pseudo-second order kinetic equation, which corresponded to a saturation mechanism. The fitted kinetic parameters achieved through the evaluated models are supplied in the Appendix A (Table A1). The adjustment of the intraparticle diffusion model showed a multi-step adsorption process. The first linear adjustment of all experiments showed a high correlation coefficient and the intercept was close to the origin, meaning that internal diffusion is a rate-limiting step due to both the linearity and the number of steps, and the low boundary layer effect was considered according to the intercept [59]. Pseudo-second order kinetics and multi-step kinetic models were also indicated by Xu et al. [60], who used black wattle tannin-modified dialdehyde nanocellulose to adsorb hexavalent chromium. These authors associated this pseudo-second order fitting to the possibility of diffusion as a rate-limiting step.
Comparing the nonlinear adjustment of the pseudo-second order for each MTMS dosage applied to CNF, the maximum adsorption capacity corresponded to 0.30 mg Cr(VI)·g−1 CNF when applying the 1.5 mmol MTMS·g−1 CNF hydrogel (Figure 3), confirming the experimental results plotted in Figure 2. Non-modified CNF and hydrophobized CNF with 1.5 mmol MTMS·g−1 worked similarly, both reaching high adsorption capacities and fast saturation, while hydrophobized CNF with 3 mmol MTMS·g−1 showed slower adsorption, and saturation of the adsorbent was not found at 24 h of operation. The results of the adsorption capacity in the equilibrium and contact time to equilibrium were close to those found by other silanized cellulose applied by Jamroz et al. [33] to adsorb hexavalent chromium of 0.30 mg·g−1 at 300 min, respectively. However, these authors applied hydrophobic cellulose to measure ultra-traces of Cr(VI) in water. They also found that the pseudo-second order kinetic was the best fitting model.
The slightly MTMS-coated CNF adsorbed chromium better than the naked CNF as the CNF surface groups were hidden by the coating layer, thus reducing the interaction of these groups with water. Furthermore, this technique prevented electrostatic repulsion between the negative CNF-surface charges and the hexavalent chromate ions. Nevertheless, excessive coating of MTMS over the surface of CNF supposed a mass transfer limitation, minimizing the adsorption rate. For these reasons, the selected MTMS dosage for the resting adsorption tests was 1.5 mmol MTMS·g−1. Compared to other silanized materials applied for hexavalent chromium removal, this amount of silanization agent was considerably lower. Around 40 mmol·g−1 of silanization reagents ((3-aminopropyl)trimethoxysilane, [3-(2-aminoethylamino)propyl]trimethoxysilane, and 3-[2-(2-aminoethylamino)ethylamino]propyl-trimethoxysilane) were added to coat the graphene oxide to adsorb hexavalent chromium from water [32].

3.2.2. Effect of pH

The experimental evolution of the hexavalent chromium concentration along the batch adsorption tests by modifying the pH from pH 3 to pH 9 is shown in Figure 4. These adsorption tests were carried out by adding the best tested operating condition of MTMS dosage, 0.1 mg·L−1 of initial chromium concentration, and 1000 mg·L−1 of CNF dosage. The trend of the Cr(VI) removal rate was similar under neutral (pH 7) and alkali (pH 9) conditions, while the adsorption was faster under acidic (pH 3) conditions, with a total removal of 80% after 1 h of operation (Figure 4). The high adsorption capacity found at pH 3 was associated with the equilibrium changes of hexavalent chromium under acidic conditions. While divalent chromate is the predominant specie when pH >6, the monovalent specie is mainly present between pH 2 and 4. Therefore, the amount of adsorbate, hexavalent chromium, is doubled at acid pH because only one active site is required per anion. A similar result of pH optimization can be found with independence of the kind of adsorbent, as indicated by Owlad et al. [9] and Saha and Orvig [61]. It can be concluded that the pH effect is related to the adsorbate ionic forms, being the effect on the adsorbent negligible. After these experiments, the pH 3 condition was selected for the rest of the optimization process.
The obtained adsorption capacity data were evaluated through different kinetic models. The result of plotting each fitted kinetic to the pH 3 experimental data can be observed in Figure 5.
The adjustments demonstrate that the most representative model was the pseudo-second kinetic model (Figure 5). The fitted kinetic parameters obtained in the different pH experiments are shown in Table A2. The intraparticle model analysis showed a multi-step adsorption mechanism with independence in the pH, as seen in the MTMS optimization. The representation of the pseudo-second order kinetic adjustment to the experimental data from each experiment is seen in Figure 6. The maximum adsorption capacity was achieved while operating under pH 3 conditions, when 0.30 mg·g−1 of the adsorption capacity and 80% of the maximum removal was reached in 1 h. The operation under neutral and alkaline media was similar in adsorption rate, but with a slightly lower adsorption capacity when operating at alkaline conditions.

3.2.3. Effect of Adsorbent Dosage

Following the same experimental procedure of the previous optimization processes, the adsorbent dosage effect was also evaluated from 250 up to 1000 mg CNF·L−1 (Figure 7). The kinetic curves showed a clear tendency of increasing the adsorption rate as the dosage of CNF increased. Among the studied CNF dosages, the minimum one needed to obtain a 100% of adsorption of hexavalent chromium was 500 mg·L−1 whereas the lowest adsorbent dosage studied showed a fast saturation of the CNF and a total removal below 10% (Figure 7). The highest dosage studied, 1000 mg·L−1, achieved the complete removal of hexavalent chromium after 25 h of contact time, the same as the dosage of 500 mg·L−1. However, the higher the dosage, the faster the adsorption, as expected. In terms of the adsorption capacity, the optimum value was found while applying 500 mg·L−1, reaching the largest efficiency in chromium removal per gram of CNF.
The quantity of hydrophobic CNF needed to reach the complete removal of hexavalent chromium is in most cases lower, compared to other cellulosic materials. Xu et al. [60] reached the maximum adsorption capacity of hexavalent chromium using 500 mg·L−1 with black wattle tannin-modified dialdehyde nanocellulose. Qiu et al. [62] and Huang et al. [58] achieved the complete depletion of hexavalent chromium employing 3000 mg·L−1 of polyethylenimine facilitated ethyl cellulose and diethylenetriamine-modified hydroxypropyl methylcellulose, respectively. In the study developed by Singh et al. [63], a dosage of 5000 mg·L−1 of aminated cellulose nanocrystals adsorbed 98.33% of hexavalent chromium. These results depict that hydrophobic CNF reached equivalent yields compared to other cellulosic materials, reducing the dosage by an order of magnitude.
Compared to other kinds of adsorbents, the dosage of hydrophobic CNF can still be considered reduced. In the case of waste material adsorbents such as treated sawdust [64], pristine almond green hull [65], and Hibiscus cannabinus kenaf [66] were applied to treat hexavalent chromium under optimal conditions, 1600, 2000 and 3000 mg·L−1 were necessary, respectively. Other typical materials used to adsorb hexavalent chromium are both commercial and mango kernel-synthetized activated carbons [67,68], and chitosan microparticles and nanoparticles [69], but they require larger adsorbent dosages than the hydrophobic CNF. The optimal dosage needed for the activated carbons was 2000 mg·L−1, while 800 mg·L−1 was the minimum required dose of the chitosan microparticles and nanoparticles to reach hexavalent chromium removal.
The results of adjusting the selected kinetic models to a dosage of 500 mg·L−1 are shown in Figure 8. The optimal kinetic model corresponded to a pseudo-second order equation. All of the kinetic parameters obtained through the kinetic fittings from the dosage optimization experiments are shown in detail in Table A3. The results of the intraparticle diffusion model adjustments were similar for each dosage, corresponding to multi-step adsorption.
Figure 9 represents the final pseudo-second order kinetic curves fitted to the experimental data from adsorbent dosage evaluation. The evolution observed revealed a maximum adsorption capacity of 0.58 mg·g−1 when 500 mg CNF·L−1 was dosed compared to 0.30 mg·g−1 in the case of 1000 mg CNF·L−1. However, the pseudo-second order kinetic constant was fourteen-times higher in the case of the highest dosage. The contrast between reaching fast-equilibrium with high adsorbent dosages and increasing the adsorption efficiency can be observed through the use of other nanomaterial adsorbents such as Fe2O3 nanoparticles [70]. This fact would suggest that the optimal dosage would depend on the objective between reaching the rapid total removal or maximizing the total capacity of the adsorbent.

3.2.4. Effect of Initial Chromium Concentration

Hexavalent chromium concentration was modified to study the adsorption capacity of CNF when the initial chromium concentration was between 0.1 and 5 mg·L−1 and when the initial chromium concentration drastically increased to the range from 10 to 50 mg·L−1 (Figure 10a,b, respectively). These kinetic experiments were developed by implementing the optimal tested MTMS dosage, the optimal tested pH, and the maximum adsorbent dosage (1000 mg CNF·L−1), thus prioritizing the process kinetic over the adsorption capacity. Furthermore, this CNF hydrogel dosage was selected to prevent the CNF from extremely fast saturation.
All of the studied initial chromium concentrations showed similar trends, with an initial fast adsorption of chromium, a second stationary step, and a final slow adsorption step (Figure 10). This situation corresponded to a multilayer adsorption. At the beginning, the free active sites were easy to find, so a first layer was generated rapidly. Once the first layer was formed, there was an accumulation of anions in the CNF surface, which led to the minimization of both the driving force between the liquid and solid phases and the electrostatic force interaction. Therefore, the stationary and intermediate step is caused by the slow reduction of anionic hexavalent chromium to less toxic cationic trivalent chromium. This reduction of chromium through the oxidation of C–OH surface groups to carboxylic groups was studied deeply by Wang and Lee [19]. Afterward, the attachment of the second layer of hexavalent chromium onto the CNF surface begins once the surface gradually becomes positively charged and the attraction force of anionic hexavalent chromate from the wastewater to the cationic surface of the CNF overpasses the mass transfer limitations. From the results shown in Figure 10, the resulting kinetic curves of hydrophobic CNF adsorption when modifying the initial chromium concentration showed an increase in the adsorption capacity (from 0.30 to 70.38) as the initial concentration of chromium increased, coupled with a reduction in the adsorption efficiency (from 100% to 51%) as the initial chromium concentration increased. This fact demonstrates that both parameters, adsorption capacity and adsorption efficiency, showed opposite trends, in agreement with Pourfadakari et al. [57], as for a selected dosage, increasing the initial concentration involves both a lesser percentage of chromium removal and a higher driving force between the liquid and the solid surface. This is due to the difference in the concentrations, which enhances the efficient usage per mass unit of material, reaching a higher amount of adsorbed chromium per active site. These authors [57] also showed the change in the shape of the curves while increasing the initial concentration, where the presence of different steps became clearer as the concentration rose.
The achieved experimental data were fitted through the kinetic models applied previously, with the initial hexavalent chromium concentration of 25 mg Cr(VI)·L−1 (Figure 11). At this initial concentration, as it happened at 1 mg·L−1, the only model that allowed an adequate interpretation and simulation of the experimental results was the intraparticle diffusion model. Each kinetic parameter obtained by kinetic analysis for all the kinetic models can be checked in Table A4. The intraparticle diffusion model predicts a three-step adsorption mechanism including the first high adsorption rate step, a steady-state step, and the last slow adsorption step, confirming the trend seen in the experimental data. Evaluating the intraparticle adjustment of the first step in the experiments using an initial concentration of 5 mg·L−1 and above, the intercept is extremely high, in some cases, close to the maximum adsorption capacity. This fact indicates a strong effect of the boundary layer, which means that external diffusion limitation will play a major role in the overall adsorption rate [59].
The resulting curves of each initial chromium concentration experiment adjusted to the pseudo-second order kinetics (Figure 12a,b) showed the change in the mechanism of adsorption. At the lower initial concentration of hexavalent chromium, the adsorption mechanism was mainly saturation, whereas as the initial concentration increased, the multistep adsorption mechanism was more noticed. The final values of both equilibrium concentrations of chromium and the adsorption capacities of CNF achieved in these tests were used to plot the isotherm curves. The low correlation coefficients for the kinetic models were due to the assumption that all of the selected kinetic equations made, as they suppose the gradual increase in the adsorption capacity or continuous adsorption up to saturation of the adsorbent, but none of these models assumed a two-layer adsorption. Nevertheless, as it was previously pointed out, the kinetic data showed a third step in the adsorption process, corresponding to the adsorption at long contact times, longer than 20 h of batch adsorption. These long contact times are not feasible at an industrial scale; therefore, the selection of experimental data with shorter contact times (first and second steps) would enhance the adjustment of these models.
There are a large variety of adsorbents to remove hexavalent chromium and many results related to how the contact between water and adsorbent is made, and to the different surface modifications, which affects the mechanism of adsorption. The hydrophobic CNF hydrogel showed the slowest equilibrium time regarding other adsorbent materials based on nanocellulose or on activated carbon (Table 3), but it reached the complete abatement of chromium and good adsorption capacity. Similar results of the maximum adsorption capacity in the equilibrium using polyacrylonitrile-modified CNF membranes were reported by Yang et al. [71] (87.5 mg·g−1). The equilibrium time of the hydrophobic CNF was similar to the one of another silanized cellulosic material applied by Jamroz et al. [33] (300 min) to detect ultra-trace concentrations of hexavalent chromium in water, suggesting that silanization processes lead to larger contact times than other kinds of adsorbents (Table 4) and could be associated with relevant mass transfer limitations due to the silane reaction with celluloses. The applied dosage was in the order of magnitude of other cellulosic adsorbents such as polypyrrole-bacterial CNF and polyaniline-functionalized CNF and was lower than that of activated carbons [72,73,74,75,76]. This comparison suggests that the hydrophobized CNF hydrogel is an efficient material for hexavalent chromium adsorption from wastewater compared to other adsorbents including other nanocellulosic materials.

3.3. Isotherm Analysis

The isotherm curves were graphed with the equilibrium data of the chromium concentration. The equilibrium data of the maximum adsorption capacities and final concentration at the equilibrium stage of hexavalent chromium treated with the hydrophobic CNF hydrogel were analyzed and adjusted to different isotherm models. The isotherm data showed an exponential increase trend instead of a saturation trend (Figure 13). According to the classification established by McCabe, Smith, and Harriott [78], this concave upward curve involves a large effect of mass transfer limitations and indicates unfavorable adsorption. This mass transfer limitation can be associated with the coating of MTMS on the surface of the material, which implies more tortuosity to reach the active sites. The parameters achieved from the adjustment of the different isotherm models are summarized in Table 3. The optimal correlation parameters are found when the Freundlich isotherm is applied to experimental data, which means the distribution of the energy through an exponential equation, and considering the heterogeneous dispersion of the active sites, happens over the adsorbent surface [79]. This mechanism supposes a multilayer adsorption [80], which fits the trend seen during the adsorption kinetics in Section 3.2.
Figure 13. The isotherm experimental data and isotherm model adjustment of the Langmuir, Freundlich, Temkin, Dubinin–Raduskevich, and Sips equations.
Figure 13. The isotherm experimental data and isotherm model adjustment of the Langmuir, Freundlich, Temkin, Dubinin–Raduskevich, and Sips equations.
Polymers 14 03425 g013
Table 4. The results of the isotherm model adjustment to the adsorption equilibrium data of hexavalent chromium on a CNF hydrogel.
Table 4. The results of the isotherm model adjustment to the adsorption equilibrium data of hexavalent chromium on a CNF hydrogel.
ModelParametersValues
LangmuirIsotherm parameterskL [L·mg−1] = 21.26
qe [mg·g−1] = 0.3417
RL (C0 = 0.1 mg·L−1) [-] = 0.9670
RL (C0 = 50 mg·L−1) [-] = 5.53·10−2
Correlation parametersR2 = 0.7420
RSS = 2949.55
FreundlichIsotherm parameterskF [mg(1−1/n)-L(1/n)·g−1] = 1.3914
nF [-] = 0.8404
Correlation parametersR2 = 0.9902
RSS = 108.01
Dubinin–RaduskevichIsotherm parametersBDR [mol2·J−2] = 9.93·10−8
qmax [mg·g−1] = 27.72
Thermodynamic parametersEDR [J·mol−1] = 2243.50
Correlation parametersR2 = 0.5754
RSS = 2542.39
TemkinIsotherm parametersBT [J·mol−1] = 12.83
bT [-J·mol−1] = 188.08
AT [L·g−1] = 1.3759
Correlation parametersR2 = 0.7481
RSS = 1415.93
SipsIsotherm parametersnS [-] = 1.2442
kS [L(1/nS)·mol-(1/nS)] = 6.16·10−2
Correlation parametersR2 = 0.9023
RSS = 1529.83
The value of the parameter nF <1 indicates an unfavorable process as well as the small bond adsorbate–adsorbent compared to a favorable process [81]. A similar value of nF was reported by Dawodu et al. [82] for hexavalent chromium adsorption onto the seed coat biomass, showing a cooperative adsorption between Cr(VI) and adsorbent surface.
As expected by the shape of the isotherm, the selected models, which supposed the saturation of the adsorbent such as Langmuir, Sips, or Dubinin–Raduskevich, showed low correlation coefficients and a high residual sum of squares. These could represent the initial part of the curve, but failed the clear concave shape adjustment in the end. For this reason, these models did not offer a good estimation of the maximum adsorption capacity reached by the hydrophobic CNF.
The Langmuir’s separation factor values calculated for the minimum and maximum initial concentrations were in the interval of 0 < RL < 1, but close to the upper and lower limits of the interval, respectively. The variability of the separation factor indicates that while treating lower concentrations, the high value of RL implies a reduced affinity adsorbate–adsorbent. On the other hand, the treatment of concentrated solutions showed a reduction in the RL values close to 0, strengthening the chromium attachment onto the CNF surface [83].
The mean free energy of adsorption calculated through the Dubinin–Raduskevich model was 2.24 kJ·mol−1, similar to the typical values indicated for the physisorption of chromium (<8 kJ·mol−1) [84]. The value of the Temkin bT parameter, which is related to the heat of sorption, was 0.19 kJ·mol−1. Choudhary and Paul [85] indicated that values of bT below 8 kJ·mol−1 revealed a weak interaction chromium-CNF surface. These low bT values are related to physisorption processes, where the values of adsorption enthalpy are in the order of physical processes such as intermolecular forces.

4. Conclusions

The modified CNF hydrogels successfully removed more than 97% of Cr(VI), opening promising applications for these nanomaterials as heavy metal adsorbents. The nanofibers showed a large amount of active carboxylic and anionic groups. Thus, it was necessary to develop a surface modification through hydrophobization treatment. The inclusion of 1.5 mmol MTMS·g−1 as a hydrophobizing agent allowed for the increase in the Cr(VI) kinetic constant of adsorption k2 by 84.97%. The pseudo-second order and intraparticle diffusion kinetic models were the best fitting models. These models revealed that both the sorption rate and hexavalent chromium diffusion played a major role as rate-limiting steps. The adsorption mechanism is ruled by the multi-step adsorption of hexavalent chromium on the CNF hydrogel dominated by internal diffusion at low concentrations and external diffusion with concentrations above 5 mg·L−1. The optimized conditions were found to be pH 3 and a dosage over 500 mg·L−1. More than 97% of hexavalent chromium removal was reached, treating concentrations below 1 mg·L−1 and the maximum adsorption capacity of 70.38 mg·g−1 was achieved at 50 mg·L−1. The isotherm analysis showed that Freundlich was the best fitting model, meaning that multilayer adsorption and the heterogeneous dispersion of surface energy is the main adsorption mechanism of hexavalent chromium onto the surface of hydrophobized CNF. The Freundlich unfavorable isotherm predicts a multilayer adsorption and a weak interaction between hexavalent chromium and CNF, associated with a physical sorption mechanism. The relatively low values of mean free energy of adsorption (2.24 kJ·mol−1) and heat of sorption (0.19 kJ·mol−1) calculated through the Dubinin–Raduskevich and Temkin models are also indicators of a physical sorption process. In general terms, the hydrophobized CNF hydrogel reached relatively high adsorption capacities compared to previously developed cellulose nanomaterials and activated carbons, and the complete removal of chromium could be found, even at low adsorbent dosages.

Author Contributions

Conceptualization, C.N., N.M., and F.d.B.O.; Methodology, F.d.B.O., N.M., and C.N.; Software, F.d.B.O.; Validation, F.d.B.O., A.B. (Ana Balea), and J.L.S.-S.; Formal analysis, F.d.B.O., A.B. (Angeles Blanco), A.B. (Ana Balea), C.N., and N.M.; Investigation, F.d.B.O. and S.M.; Resources, C.N. and A.B. (Angeles Blanco); Data curation, F.d.B.O., A.B. (Ana Balea), and J.L.S.-S.; Writing—original draft preparation, F.d.B.O., A.B. (Ana Balea), and N.M.; Writing—review and editing, F.d.B.O., A.B. (Ana Balea), and C.N.; Supervision, C.N., N.M., and A.B. (Angeles Blanco); Project administration, C.N. and A.B. (Angeles Blanco); Funding acquisition, C.N. and A.B. (Angeles Blanco) All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Community of Madrid, project S2018/EMT-4459 (RETO-PROSOST-CM) and the Economy and Competitiveness Ministry of Spain, project CTQ2017-85654-C2-2-R and predoctoral fellowship PRE2018-085034.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors wish to thank the Community of Madrid and the Economy and Competitiveness Ministry of Spain for their support to projects S2018/EMT-4459 (RETO-PROSOST-CM) and CTQ2017-85654-C2-2-R as well as the support of the Economy and Competitiveness Ministry of Spain for the grant provided to F.B. Ojembarrena (PRE2018-085034).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. The results of the kinetic equation adjustment to the adsorption experimental data of 0–3 mmol MTMS·g−1 CNF.
Table A1. The results of the kinetic equation adjustment to the adsorption experimental data of 0–3 mmol MTMS·g−1 CNF.
Kinetic Model 0 mmol MTMS·g−1 CNF1.5 mmol MTMS·g−1 CNF3 mmol MTMS·g−1 CNF
Pseudo-first orderKinetic parametersk1 [h−1] = 0.2378k1 [h−1] = 8.63·10−2k1 [h−1] = 9.16·10−2
Correlation parametersR2 = 0.9653
RSS = 7.11·10−3
R2 = 0.7725
RSS = 7.84·10−2
R2 = 0.9674
RSS = 9.38·10−3
Pseudo-second orderKinetic parametersk2 [mg·g−1-h−1] = 3.7399
qe [mg·g−1] = 0.2933
k2 [mg·g−1-h−1] = 6.9175
qe [mg·g−1] = 0.3058
k2 [mg·g−1-h−1] = 1.4827
qe [mg·g−1] = 0.2610
Correlation parametersR2 = 0.9440
RSS = 1.23·10−2
R2 = 0.9661
RSS = 5.96·10−3
R2 = 0.9430
RSS = 1.08·10−2
ElovichKinetic parametersα [h·mg·g−1] = 0.8350
β [g·mg−1] = 21.5517
α [h·mg·g−1] = 6.6752
β [g·mg−1] = 28.1690
α [h·mg·g−1] = 1.6142
β [g·mg−1] = 37.7358
Correlation parametersR2 = 0.8109
RSS = 1.23·10−2
R2 = 0.8005
RSS = 7.74·10−3
R2 = 0.9540
RSS = 3.90·10−3
Intraparticle diffusionKinetic parameters: Step 1ki,1 [mg·g−1·min−0.5] = 0.1193
Ci,1 [mg·g−1] = −1.62·10−2
ki,1 [mg·g−1·min−0.5] = 0.1462
Ci,1 [mg·g−1] = −2.30·10−3
ki,1 [mg·g−1·min−0.5] = 0.1462
Ci,1 [mg·g−1] = 7.00·10−4
Correlation parametersR2 = 0.9709
RSS = 1.23·10−4
R2 = 0.9996
RSS = 3.17·10−5
R2 = 0.9918
RSS = 5.06·10−5
Kinetic parameters: Step 2ki,2 [mg·g−1·min−0.5] = 1·10−17
Ci,1 [mg·g−1] = 0.299
ki,2 [mg·g−1·min−0.5] = 4.00·10−2
Ci,2 [mg·g−1] = 0.196
ki,2 [mg·g−1·min−0.5] = 0.0239
Ci,1 [mg·g−1] = 0.0658
Correlation parametersR2 = 0.5477
RSS = 1.23·10−4
R2 = 0.9996
RSS = 4.25·10−8
R2 = 0.9839
RSS = 3.67·10−4
Kinetic parameters: Step 3 ki,3 [mg·g−1·min−0.5] = 2.10·10−3
Ci,3 [mg·g−1] = 0.2878
Correlation parameters R2 = 0.8137
RSS = 4.29·10−5

Appendix B

Table A2. The results of the kinetic equation adjustment to the adsorption experimental data of pH 3–9.
Table A2. The results of the kinetic equation adjustment to the adsorption experimental data of pH 3–9.
Kinetic Model pH 3pH 7pH 9
Pseudo-first orderKinetic parametersk1 [h−1] = 5.55·10−2k1 [h−1] = 6.09·10−2k1 [h−1] = 0.2179
Correlation parametersR2 = 0.6366
RSS = 0.167
R2 = 0.9219
RSS = 7.25·10−3
R2 = 0.9960
RSS = 1.49·10−4
Pseudo-second orderKinetic parametersk2 [mg·g−1-h−1] = 11.4405
qe [mg·g−1] = 0.3050
k2 [mg·g−1-h−1] = 0.9889
qe [mg·g−1] = 0.1821
k2 [mg·g−1-h−1] = 4.4596
qe [mg·g−1] = 0.1139
Correlation parametersR2 = 0.9980
RSS = 3.29·10−4
R2 = 0.9612
RSS = 2.77·10−3
R2 = 0.9801
RSS = 8.14·10−4
ElovichKinetic parametersα [h·mg·g−1] = 43.5512
β [g·mg−1] = 33.7838
α [h·mg·g−1] = 0.3496
β [g·mg−1] = 46.9484
α [h·mg·g−1] = 0.3013
β [g·mg−1] = 66.2252
Correlation parametersR2 = 0.8683
RSS = 5.44·10−3
R2 = 0.8989
RSS = 6.16·10−3
R2 = 0.9255
RSS = 2.18·10−3
Intraparticle diffusionKinetic parameters: Step 1ki,1 [mg·g−1·min−0.5] = 0.2346
Ci,1 [mg·g−1] = 1.38·10−2
ki,1 [mg·g−1·min−0.5] = 2.41·10−2
Ci,1 [mg·g−1] = 4.75·10−4
ki,1[mg·g−1·min−0.5] = 8.34·10−3
Ci,1 [mg·g−1] = 4.63·10−3
Correlation parametersR2 = 0.9995
RSS = 2.13·10−10
R2 = 0.9841
RSS = 5.36·10−5
R2 = 0.9055
RSS = 2.73·10−5
Kinetic parameters: Step 2ki,2 [mg·g−1·min−0.5] = 3.01·10−2
Ci,2 [mg·g−1] = 0.2178
ki,2 [mg·g−1·min−0.5] = 0.2399
Ci,2 [mg·g−1] = −0.4314
ki,2[mg·g−1·min−0.5] = 4.82·10−2
Ci,2 [mg·g−1] = −3.30·10−2
Correlation parametersR2 = 0.9965
RSS = 6.23·10−6
R2 = 0.9999
RSS = 2.13·10−8
R2 = 0.9891
RSS = 5.05·10−5
Kinetic parameters: Step 3ki,3 [mg·g−1·min−0.5] = 2.10·10−3
Ci,3 [mg·g−1] = 0.288
ki,3 [mg·g−1·min−0.5]= 7.05·10−3
Ci,3 [mg·g−1] = 0.1079
ki,3[mg·g−1·min−0.5] = 1.35·10−2
Ci,3 [mg·g−1] = 4.42·10−2
Correlation parametersR2 = 0.8137
RSS = 4.29·10−5
R2 = 0.9103
RSS = 1.96·10−4
R2 = 0.9999
RSS = 6.74·10−9
Kinetic parameters: Step 4 ki,4 [mg·g−1·min−0.5] = 0
Ci,4 [mg·g−1] = 0.1114
Correlation parameters RSS = 1.98·10−6

Appendix C

Table A3. The results of the kinetic equation adjustment to the adsorption experimental data of 250 to 1000 mg CNF·L−1.
Table A3. The results of the kinetic equation adjustment to the adsorption experimental data of 250 to 1000 mg CNF·L−1.
Kinetic Model 250 mg·L−1500 mg·L−11000 mg·L−1
Pseudo-first orderKinetic parametersk1 [h−1] = 4.49·10−2k1 [h−1] = 5.56·10−2k1 [h−1] = 6.69·10−2
Correlation parametersR2 = 0.3713
RSS = 3.31·10−3
R2 = 0.9267
RSS = 0.159
R2 = 0.6605
RSS = 0.152
Pseudo-second orderKinetic parametersk2 [mg·g−1-h−1] = 3822.25
qe [mg·g−1] = 2.76·10−2
k2 [mg·g−1-h−1] = 1.1320
qe [mg·g−1] = 0.5787
k2 [mg·g−1-h−1] = 14.3041
qe [mg·g−1] = 0.3047
Correlation parametersR2 = 0.9256
RSS = 1.55·10−4
R2 = 0.9596
RSS = 3.78·10−2
R2 = 0.9957
RSS = 9.31·10−4
ElovichKinetic parametersα [h·mg·g−1] = 229.62
β [g·mg−1] = 462.96
α [h·mg·g−1] = 3.2448
β [g·mg−1] = 15.1976
α [h·mg·g−1] = 66.5435
β [g·mg−1] = 35.9712
Correlation parametersR2 = 0.6241
RSS = 2.78·10−4
R2 = 0.9246
RSS = 4.49· 10−2
R2 = 0.8683
RSS = 5.34·10−3
Intraparticle diffusionKinetic parameters: Step 1ki,1 [mg·g−1·min−0.5] = 0.1384
Ci,1 [mg·g−1] = −1.23·10−2
ki,1 [mg·g−1·min−0.5] = 0.1397
Ci,1 [mg·g−1] = 5.45·10−2
ki,1 [mg·g−1·min−0.5] = 0.2346
Ci,1 [mg·g−1] = 1.38·10−2
Correlation parametersR2 = 0.9631
RSS = 1.51·10−4
R2 = 0.9295
RSS = 1.03·10−2
R2 = 0.9995
RSS = 2.13·10−10
Kinetic parameters: Step 2ki,2 [mg·g−1·min−0.5] = 1·10·−17
Ci,2 [mg·g−1] = 2.76·10−2
ki,2 [mg·g−1·min−0.5] = 7.03·10−2
Ci,2 [mg·g−1] = 0.2210
ki,2[mg·g−1·min−0.5] = 3.01·10−2
Ci,2 [mg·g−1] = 0.2178
Correlation parametersR2 = 0.5117
RSS = 1.44·10−11
R2 = 1.0000
RSS = 4.19·10−7
R2 = 0.9965
RSS = 6.23·10−6
Kinetic parameters: Step 3 ki,3 [mg·g−1·min−0.5] = 0
Ci,3 [mg·g−1] = 0.5658
ki,3[mg·g−1·min−0.5] = 2.10·10−3
Ci,3 [mg·g−1] = 0.2880
Correlation parameters RSS = 1.73·10−9R2 = 0.8137
RSS = 4.29·10−5

Appendix D

Table A4. The results of the kinetic equation adjustment to the adsorption experimental data of the initial hexavalent chromium concentrations from 0.1 to 50 mg·L−1.
Table A4. The results of the kinetic equation adjustment to the adsorption experimental data of the initial hexavalent chromium concentrations from 0.1 to 50 mg·L−1.
Kinetic Model 0.1 mg·L−11 mg·L−15 mg·L−1
Pseudo-first orderKinetic parametersk1 [h−1] = 8.00·10−2k1 [h−1] = 4.76·10−2k1 [h−1] = 1.5675
Correlation parametersR2 = 0.8105
RSS = 6.81·10−2
R2 = 0.8880
RSS = 1.1163
R2 = 0.7040
RSS = 35.56
Pseudo-second orderKinetic parametersk2 [mg·g−1-h−1] = 4.9284
qe [mg·g−1] = 0.3059
k2 [mg·g−1-h−1] = 3.53·10−2
qe [mg·g−1] = 1.9361
k2 [mg·g−1-h−1] = 2.99·10−2
qe [mg·g−1] = 6.9541
Correlation parametersR2 = 0.9853
RSS = 3.05·10−3
R2 = 0.8975
RSS = 0.7128
R2 = 0.5884
RSS = 54.66
ElovichKinetic parametersα [h·mg·g−1] = 8.6542
β [g·mg−1] = 31.257
α [h·mg·g−1] = 1.6140
β [g·mg−1] = 4.7996
α [h·mg·g−1] = 1.20·109
β [g·mg−1] = 4.7547
Correlation parametersR2 = 0.8606
RSS = 5.62·10−3
R2 = 0.8760
RSS = 0.5550
R2 = 0.5875
RSS = 3.5383
Intraparticle diffusionKinetic parameters: Step 1ki,1 [mg·g−1·min−0.5] = 0.1099
Ci,1 [mg·g−1] = 2.88·10−2
ki,1 [mg·g−1·min−0.5] = 0.6325
Ci,1 [mg·g−1] = −0.1448
ki,1[mg·g−1·min−0.5] = 0.3298
Ci,1 [mg·g−1] = 4.4590
Correlation parametersR2 = 0.9982
RSS = 1.68·10−4
R2 = 0.9271
RSS = 3.91·10−2
R2 = 0.8235
RSS = 19.88
Kinetic parameters: Step 2ki,2 [mg·g−1·min−0.5] = 6.00·10−3
Ci,2 [mg·g−1] = 0.2615
ki,2 [mg·g−1·min−0.5] = 2.29·10−2
Ci,2 [mg·g−1] = 0.5194
ki,2[mg·g−1·min−0.5] = 0
Ci,2 [mg·g−1]= 4.6922
Correlation parametersR2 = 0.9978
RSS = 1.66·10−6
R2 = 0.9690
RSS = 2.79·10−4
RSS = 9.68·10−9
Kinetic parameters: Step 3ki,3 [mg·g−1·min−0.5] = 0
Ci,3 [mg·g−1] = 0.3036
ki,3 [mg·g−1·min−0.5]= 0.2121
Ci,3 [mg·g−1] = −0.3775
ki,3[mg·g−1·min−0.5] = 0.617
Ci,3 [mg·g−1] = 0.4177
Correlation parametersRSS = 1.68·10−4R2 = 0.9986RSS = 2.48·10−3R2 = 0.9999RSS = 7.46·10−8
Kinetic model 10 mg/L25 mg/L50 mg/L
Pseudo-first orderKinetic parametersk1 [h−1] = 3.00·10−2k1 [h−1] = 0.1664k1 [h−1] = 3.86·10−2
Correlation parametersR2 = 0.5969
RSS = 307.88
R2 = 0.5860
RSS = 1592.74
R2 = 0.8982
RSS = 2433.18
Pseudo-second orderKinetic parametersk2 [mg·g−1-h−1] = 5.40·10−2
qe [mg·g−1] = 12.115
k2 [mg·g−1-h−1] = 8.70·10−3
qe [mg·g−1] = 34.51
k2 [mg·g−1-h−1] = 6.16·10−3
qe [mg·g−1] = 70.92
Correlation parametersR2 = 0.7242
RSS = 121.03
R2 = 0.6332
RSS = 1286.27
R2 = 0.9851
RSS = 236.03
ElovichKinetic parametersα [h·mg·g−1] = 3.82·108
β [g·mg−1] = 2.2222
α [h·mg·g−1] = 7.29·108
β [g·mg−1] = 0.8143
α [h·mg·g−1] = 153.74
β [g·mg−1] = 0.1041
Correlation parametersR2 = 0.9184
RSS = 1.5841
R2 = 0.7408
RSS = 44.33
R2 = 0.9934
RSS = 166.49
Intraparticle diffusionKinetic parameters: Step 1ki,1 [mg·g−1·min−0.5] = 1.3193
Ci,1 [mg·g−1] = 8.1756
ki,1 [mg·g−1·min−0.5] = 4.9473
Ci,1 [mg·g−1] = 21.343
ki,1[mg·g−1·min−0.5] = 27.162
Ci,1 [mg·g−1] = 4.5797
Correlation parametersR2 = 0.8392
RSS = 66.84
R2 = 0.8504
RSS = 455.52
R2 = 0.9935
RSS = 20.97
Kinetic parameters: Step 2ki,2 [mg·g−1·min−0.5] = 0
Ci,2 [mg·g−1] = 9.1085
ki,2 [mg·g−1·min−0.5] = 0.1966
Ci,2 [mg·g−1] = 24.783
ki,2[mg·g−1·min−0.5] = 8.0237
Ci,2 [mg·g−1] = 24.416
Correlation parametersRSS = 2.08·10−9R2 = 0.9204
RSS = 0.2037
R2 = 0.9987
RSS = 1.4264
Kinetic parameters: Step 3ki,3 [mg·g−1·min−0.5] = 0.3399
Ci,3 [mg·g−1] = 8.5209
ki,3 [mg·g−1·min−0.5] = 2.2280
Ci,3 [mg·g−1] = 10.785
ki,3[mg·g−1·min−0.5] = 1.2205
Ci,3 [mg·g−1] = 56.515
Correlation parametersR2 = 0.9953
RSS = 4.63·10−2
R2 = 0.9999
RSS = 2.05·10−7
R2 = 0.9443
RSS = 3.4399

References

  1. Outridge, P.M.; Scheuhammer, A.M. Bioaccumulation and toxicology of chromium: Implications for wildlife. Rev. Environ. Contam. Toxicol. 1993, 130, 31–77. [Google Scholar]
  2. Kumar, V.; Dwivedi, S.K. A review on accessible techniques for removal of hexavalent Chromium and divalent Nickel from industrial wastewater: Recent research and future outlook. J. Clean. Prod. 2021, 295, 126229. [Google Scholar] [CrossRef]
  3. Mishra, S.; Bharagava, R.N. Toxic and genotoxic effects of hexavalent chromium in environment and its bioremediation strategies. J. Environ. Sci. Health C Environ. Carcinog. Ecotoxicol. Rev. 2016, 34, 1–32. [Google Scholar] [CrossRef] [PubMed]
  4. Wise, J.P.; Young, J.L.; Cai, J.; Cai, L. Current understanding of hexavalent chromium [Cr(VI)] neurotoxicity and new perspectives. Environ. Int. 2022, 158, 106877. [Google Scholar] [CrossRef] [PubMed]
  5. Sharma, P.; Singh, S.P.; Parakh, S.K.; Tong, Y.W. Health hazards of hexavalent chromium (Cr (VI)) and its microbial reduction. Bioengineered 2022, 13, 4923–4938. [Google Scholar] [CrossRef] [PubMed]
  6. Chai, W.S.; Cheun, J.Y.; Kumar, P.S.; Mubashir, M.; Majeed, Z.; Banat, F.; Ho, S.-H.; Show, P.L. A review on conventional and novel materials towards heavy metal adsorption in wastewater treatment application. J. Clean. Prod. 2021, 296, 126589. [Google Scholar] [CrossRef]
  7. Vaiopoulou, E.; Gikas, P. Regulations for chromium emissions to the aquatic environment in Europe and elsewhere. Chemosphere 2020, 254, 126876. [Google Scholar] [CrossRef]
  8. Farooq, U.; Kozinski, J.A.; Khan, M.A.; Athar, M. Biosorption of heavy metal ions using wheat based biosorbents—A review of the recent literature. Bioresour. Technol. 2010, 101, 5043–5053. [Google Scholar] [CrossRef] [PubMed]
  9. Owlad, M.; Aroua, M.K.; Daud, W.A.W.; Baroutian, S. Removal of Hexavalent Chromium-Contaminated Water and Wastewater: A Review. Water Air Soil Pollut. 2009, 200, 59–77. [Google Scholar] [CrossRef]
  10. GracePavithra, K.; Jaikumar, V.; Kumar, P.S.; SundarRajan, P. A review on cleaner strategies for chromium industrial wastewater: Present research and future perspective. J. Clean. Prod. 2019, 228, 580–593. [Google Scholar] [CrossRef]
  11. Islam, M.A.; Angove, M.J.; Morton, D.W. Recent innovative research on chromium (VI) adsorption mechanism. Environ. Nanotechnol. Monit. Manag. 2019, 12, 100267. [Google Scholar] [CrossRef]
  12. Singh, P.; Itankar, N.; Patil, Y. Biomanagement of hexavalent chromium: Current trends and promising perspectives. J. Environ. Manag. 2021, 279, 111547. [Google Scholar] [CrossRef] [PubMed]
  13. Rambabu, K.; Thanigaivelan, A.; Bharath, G.; Sivarajasekar, N.; Banat, F.; Show, P.L. Biosorption potential of Phoenix dactylifera coir wastes for toxic hexavalent chromium sequestration. Chemosphere 2021, 268, 128809. [Google Scholar] [CrossRef] [PubMed]
  14. Hernández-Cocoletzi, H.; Salinas, R.A.; Águila-Almanza, E.; Rubio-Rosas, E.; Chai, W.S.; Chew, K.W.; Mariscal-Hernández, C.; Show, P.L. Natural hydroxyapatite from fishbone waste for the rapid adsorption of heavy metals of aqueous effluent. Environ. Technol. Innov. 2020, 20, 101109. [Google Scholar] [CrossRef]
  15. Lee, P.-X.; Liu, B.-L.; Show, P.L.; Ooi, C.W.; Chai, W.S.; Munawaroh, H.S.H.; Chang, Y.-K. Removal of calcium ions from aqueous solution by bovine serum albumin (BSA)-modified nanofiber membrane: Dynamic adsorption performance and breakthrough analysis. Biochem. Eng. J. 2021, 171, 108016. [Google Scholar] [CrossRef]
  16. Mubashir, M.; Dumée, L.F.; Fong, Y.Y.; Jusoh, N.; Lukose, J.; Chai, W.S.; Show, P.L. Cellulose acetate-based membranes by interfacial engineering and integration of ZIF-62 glass nanoparticles for CO2 separation. J. Hazard. Mater. 2021, 415, 125639. [Google Scholar] [CrossRef]
  17. Aigbe, U.O.; Osibote, O.A. A review of hexavalent chromium removal from aqueous solutions by sorption technique using nanomaterials. J. Environ. Chem. Eng. 2020, 8, 104503. [Google Scholar] [CrossRef]
  18. Varghese, A.G.; Paul, S.A.; Latha, M.S. Remediation of heavy metals and dyes from wastewater using cellulose-based adsorbents. Environ. Chem. Lett. 2019, 17, 867–877. [Google Scholar] [CrossRef]
  19. Wang, S.-L.; Lee, J.-F. Reaction mechanism of hexavalent chromium with cellulose. Chem. Eng. J. 2011, 174, 289–295. [Google Scholar] [CrossRef]
  20. Kazachenko, A.; Akman, F.; Medimagh, M.; Issaoui, N.; Vasilieva, N.; Malyar, Y.N.; Sudakova, I.G.; Karacharov, A.; Miroshnikova, A.; Al-Dossary, O.M. Sulfation of Diethylaminoethyl-Cellulose: QTAIM Topological Analysis and Experimental and DFT Studies of the Properties. ACS Omega 2021, 6, 22603–22615. [Google Scholar] [CrossRef]
  21. Negro, C.; Balea Martin, A.; Sanchez-Salvador, J.; Campano, C.; Fuente, E.; Monte, M.; Blanco, A. Nanocellulose and its potential use for sustainable industrial applications. Lat. Am. Appl. Res.-Int. J. 2020, 50, 59–64. [Google Scholar] [CrossRef]
  22. Norrrahim, M.N.F.; Kasim, N.A.M.; Knight, V.F.; Misenan, M.S.M.; Janudin, N.; Shah, N.A.A.; Kasim, N.; Yusoff, W.Y.W.; Noor, S.A.M.; Jamal, S.H. Nanocellulose: A bioadsorbent for chemical contaminant remediation. RSC Adv. 2021, 11, 7347–7368. [Google Scholar] [CrossRef] [PubMed]
  23. Reshmy, R.; Philip, E.; Madhavan, A.; Pugazhendhi, A.; Sindhu, R.; Sirohi, R.; Awasthi, M.K.; Pandey, A.; Binod, P. Nanocellulose as green material for remediation of hazardous heavy metal contaminants. J. Hazard. Mater. 2022, 424, 127516. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, F.; Hubbe, M. Charge Properties of Fibers in the Paper Mill Environment. 1. Effect of Electrical Conductivity. J. Pulp Pap. Sci. 2002, 28, 347–353. [Google Scholar]
  25. Khoo, R.Z.; Chow, W.S.; Ismail, H. Sugarcane bagasse fiber and its cellulose nanocrystals for polymer reinforcement and heavy metal adsorbent: A review. Cellulose 2018, 25, 4303–4330. [Google Scholar] [CrossRef]
  26. Liu, S.; Low, Z.-X.; Xie, Z.; Wang, H. TEMPO-Oxidized Cellulose Nanofibers: A Renewable Nanomaterial for Environmental and Energy Applications. Adv. Mater. Technol. 2021, 6, 2001180. [Google Scholar] [CrossRef]
  27. Huang, X.Y.; Dognani, G.; Hadi, P.; Yang, M.Y.; Job, A.E.; Hsiao, B.S. Cationic Dialdehyde Nanocellulose from Sugarcane Bagasse for Efficient Chromium(VI) Removal. ACS Sustain. Chem. Eng. 2020, 8, 4734–4744. [Google Scholar] [CrossRef]
  28. Ossa-Paredes, R.; Bastidas, B.; Carvajal-Barriga, E.J. Remediation of Contaminated Water with Chromium VI by Sorption in Surface-activated-nanocellulose spheroids. Pollution 2022, 8, 489–500. [Google Scholar]
  29. Ben Hammouda, S.; Chen, Z.; An, C.; Lee, K. Recent advances in developing cellulosic sorbent materials for oil spill cleanup: A state-of-the-art review. J. Clean. Prod. 2021, 311, 127630. [Google Scholar] [CrossRef]
  30. Zhang, Z.; Sèbe, G.; Rentsch, D.; Zimmermann, T.; Tingaut, P. Ultralightweight and flexible silylated nanocellulose sponges for the selective removal of oil from water. Chem. Mater. 2014, 26, 2659–2668. [Google Scholar] [CrossRef]
  31. Janik, P.; Zawisza, B.; Talik, E.; Sitko, R. Selective adsorption and determination of hexavalent chromium ions using graphene oxide modified with amino silanes. Microchim. Acta 2018, 185, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lee, J.; Park, J.-A.; Kim, H.-G.; Lee, J.-H.; Cho, S.-H.; Choi, K.; Jung, K.-W.; Lee, S.Y.; Choi, J.-W. Most suitable amino silane molecules for surface functionalization of graphene oxide toward hexavalent chromium adsorption. Chemosphere 2020, 251, 126387. [Google Scholar] [CrossRef] [PubMed]
  33. Jamroz, E.; Kocot, K.; Zawisza, B.; Talik, E.; Gagor, A.; Sitko, R. A green analytical method for ultratrace determination of hexavalent chromium ions based on micro-solid phase extraction using amino-silanized cellulose membranes. Microchem. J. 2019, 149, 104060. [Google Scholar] [CrossRef]
  34. Saito, T.; Kimura, S.; Nishiyama, Y.; Isogai, A. Cellulose Nanofibers Prepared by TEMPO-Mediated Oxidation of Native Cellulose. Biomacromolecules 2007, 8, 2485–2491. [Google Scholar] [CrossRef]
  35. Fukuzumi, H.; Saito, T.; Okita, Y.; Isogai, A. Thermal stabilization of TEMPO-oxidized cellulose. Polym. Degrad. Stab. 2010, 95, 1502–1508. [Google Scholar] [CrossRef]
  36. Balea, A.; Merayo, N.; De La Fuente, E.; Negro, C.; Blanco, Á. Assessing the influence of refining, bleaching and TEMPO-mediated oxidation on the production of more sustainable cellulose nanofibers and their application as paper additives. Ind. Crops Prod. 2017, 97, 374–387. [Google Scholar] [CrossRef]
  37. Sanchez-Salvador, J.L.; Campano, C.; Negro, C.; Monte, M.C.; Blanco, A. Increasing the Possibilities of TEMPO-Mediated Oxidation in the Production of Cellulose Nanofibers by Reducing the Reaction Time and Reusing the Reaction Medium. Adv. Sustain. Syst. 2021, 5, 2000277. [Google Scholar] [CrossRef]
  38. Balea, A.; Sanchez-Salvador, J.L.; Monte, M.C.; Merayo, N.; Negro, C.; Blanco, A. In Situ Production and Application of Cellulose Nanofibers to Improve Recycled Paper Production. Molecules 2019, 24, 1800. [Google Scholar] [CrossRef] [Green Version]
  39. Sanchez-Salvador, J.L.; Balea, A.; Monte, M.C.; Negro, C.; Miller, M.; Olson, J.; Blanco, A. Comparison Of Mechanical And Chemical Nanocellulose As Additives To Reinforce Recycled Cardboard. Sci. Rep. 2020, 10, 3778. [Google Scholar] [CrossRef]
  40. APHA. Standard Methods: 3500-Cr CHROMIUM. Chromium by colororimetry. In Standard Methods for the Examination of Water and Wastewater; APHA: Washington, DC, USA, 2018. [Google Scholar]
  41. Alsaiari, N.S.; Katubi, K.M.; Alzahrani, F.M.; Amari, A.; Osman, H.; Rebah, F.B.; Tahoon, M.A. Synthesis, Characterization and Application of Polypyrrole Functionalized Nanocellulose for the Removal of Cr(VI) from Aqueous Solution. Polymers 2021, 13, 3691. [Google Scholar] [CrossRef]
  42. Park, S.Y.; Goo, S.; Shin, H.; Kim, J.; Youn, H.J. Structural properties of cellulose nanofibril foam depending on wet foaming conditions in Pickering stabilization. Cellulose 2021, 28, 10291–10304. [Google Scholar] [CrossRef]
  43. Largitte, L.; Pasquier, R. A review of the kinetics adsorption models and their application to the adsorption of lead by an activated carbon. Chem. Eng. Res. Des. 2016, 109, 495–504. [Google Scholar] [CrossRef]
  44. Ho, Y.S.; McKay, G. A Comparison of Chemisorption Kinetic Models Applied to Pollutant Removal on Various Sorbents. Process Saf. Environ. Prot. 1998, 76, 332–340. [Google Scholar] [CrossRef] [Green Version]
  45. Wu, F.-C.; Tseng, R.-L.; Juang, R.-S. Characteristics of Elovich equation used for the analysis of adsorption kinetics in dye-chitosan systems. Chem. Eng. J. 2009, 150, 366–373. [Google Scholar] [CrossRef]
  46. Andelescu, A.; Nistor, M.A.; Muntean, S.G.; Rădulescu-Grad, M.E. Adsorption studies on copper, cadmium, and zinc ion removal from aqueous solution using magnetite/carbon nanocomposites. Sep. Sci. Technol. 2018, 53, 2352–2364. [Google Scholar] [CrossRef]
  47. Lim, L.B.L.; Priyantha, N.; Tennakoon, D.T.B.; Chieng, H.I.; Dahri, M.K.; Suklueng, M. Breadnut peel as a highly effective low-cost biosorbent for methylene blue: Equilibrium, thermodynamic and kinetic studies. Arab. J. Chem. 2017, 10, S3216–S3228. [Google Scholar] [CrossRef] [Green Version]
  48. Patiño-Masó, J.; Serra-Parareda, F.; Tarrés, Q.; Mutjé, P.; Espinach, F.X.; Delgado-Aguilar, M. TEMPO-Oxidized Cellulose Nanofibers: A Potential Bio-Based Superabsorbent for Diaper Production. Nanomaterials 2019, 9, 1271. [Google Scholar] [CrossRef] [Green Version]
  49. Lu, P.; Liu, R.; Liu, X.; Wu, M. Preparation of Self-supporting Bagasse Cellulose Nanofibrils Hydrogels Induced by Zinc Ions. Nanomaterials 2018, 8, 800. [Google Scholar] [CrossRef] [Green Version]
  50. Ghiasi, F.; Eskandari, M.H.; Golmakani, M.-T.; Hosseini, S.M.H. Development of highly stable colloidal dispersions of gelled-oil nanoparticles loaded with cuminaldehyde. J. Colloid. Interf. Sci. 2019, 541, 65–74. [Google Scholar] [CrossRef]
  51. Necula, B.S.; Apachitei, I.; Fratila-Apachitei, L.E.; Teodosiu, C.; Duszczyk, J. Stability of nano-/microsized particles in deionized water and electroless nickel solutions. J. Colloid. Interfaces.Sci. 2007, 314, 514–522. [Google Scholar] [CrossRef]
  52. Khalil, H.P.S.A.; Davoudpour, Y.; Aprilia, N.A.S.; Mustapha, A.; Hossain, S.; Islam, N.; Dungani, R. Nanocellulose-Based Polymer Nanocomposite: Isolation, Characterization and Applications. In Nanocellulose Polymer Nanocomposites; Scrivener Publishing LLC: Beverly, MA, USA, 2014; pp. 273–309. [Google Scholar]
  53. Balea, A.; Fuente, E.; Tarrés, Q.; Pèlach, M.À.; Mutjé, P.; Delgado-Aguilar, M.; Blanco, A.; Negro, C. Influence of pretreatment and mechanical nanofibrillation energy on properties of nanofibers from Aspen cellulose. Cellulose 2021, 28, 9187–9206. [Google Scholar] [CrossRef]
  54. Kim, Y.; Park, J.; Bang, J.; Kim, J.; Jin, H.-J.; Kwak, H.W. Highly efficient Cr(VI) remediation by cationic functionalized nanocellulose beads. J. Hazard. Mater. 2022, 426, 128078. [Google Scholar] [CrossRef] [PubMed]
  55. Yuan, H.; Peng, J.; Ren, T.; Luo, Q.; Luo, Y.; Zhang, N.; Huang, Y.; Guo, X.; Wu, Y. Novel fluorescent lignin-based hydrogel with cellulose nanofibers and carbon dots for highly efficient adsorption and detection of Cr(VI). Sci. Total Environ. 2021, 760, 143395. [Google Scholar] [CrossRef]
  56. Kumar, N.; Kardam, A.; Rajawat, D.S.; Jain, V.K.; Suman. Carboxymethyl nanocellulose stabilized nano zero-valent iron: An effective method for reduction of hexavalent chromium in wastewater. Mater. Res. Express 2019, 6, 1150f3. [Google Scholar] [CrossRef]
  57. Pourfadakari, S.; Jorfi, S.; Ahmadi, M.; Takdastan, A. Experimental data on adsorption of Cr(VI) from aqueous solution using nanosized cellulose fibers obtained from rice husk. Data Brief 2017, 15, 887–895. [Google Scholar] [CrossRef]
  58. Huang, H.H.; He, D.L.; Tang, Y.N.; Guo, Y.N.; Li, P.; Qv, W.; Deng, F.; Lu, F.H. Adsorption of Hexavalent Chromium from an Aqueous Phase by Hydroxypropyl Methylcellulose Modified with Diethylenetriamine. J. Chem. Eng. Data 2019, 64, 98–106. [Google Scholar] [CrossRef]
  59. Wu, F.-C.; Tseng, R.-L.; Juang, R.-S. Initial behavior of intraparticle diffusion model used in the description of adsorption kinetics. Chem. Eng. J. 2009, 153, 1–8. [Google Scholar] [CrossRef]
  60. Xu, Q.H.; Wang, Y.L.; Jin, L.Q.; Wang, Y.; Qin, M.H. Adsorption of Cu (II), Pb (II) and Cr (VI) from aqueous solutions using black wattle tannin-immobilized nanocellulose. J. Hazard. Mater. 2017, 339, 91–99. [Google Scholar] [CrossRef]
  61. Saha, B.; Orvig, C. Biosorbents for hexavalent chromium elimination from industrial and municipal effluents. Coord. Chem. Rev. 2010, 254, 2959–2972. [Google Scholar] [CrossRef]
  62. Qiu, B.; Guo, J.; Zhang, X.; Sun, D.Z.; Gu, H.B.; Wang, Q.; Wang, H.W.; Wang, X.F.; Zhang, X.; Weeks, B.L.; et al. Polyethylenimine Facilitated Ethyl Cellulose for Hexavalent Chromium Removal with a Wide pH Range. Acs Appl. Mater. Interfaces 2014, 6, 19816–19824. [Google Scholar] [CrossRef]
  63. Singh, K.; Arora, J.K.; Sinha, T.J.M.; Srivastava, S. Functionalization of nanocrystalline cellulose for decontamination of Cr(III) and Cr(VI) from aqueous system: Computational modeling approach. Clean Technol. Environ. Policy 2014, 16, 1179–1191. [Google Scholar] [CrossRef]
  64. Baral, S.S.; Das, S.N.; Rath, P. Hexavalent chromium removal from aqueous solution by adsorption on treated sawdust. Biochem. Eng. J. 2006, 31, 216–222. [Google Scholar] [CrossRef]
  65. Nasseh, N.; Taghavi, L.; Barikbin, B.; Harifi-Mood, A.R. The removal of Cr(VI) from aqueous solution by almond green hull waste material: Kinetic and equilibrium studies. J. Water Reuse Desalination 2016, 7, 449–460. [Google Scholar] [CrossRef]
  66. Omidvar Borna, M.; Pirsaheb, M.; Vosoughi Niri, M.; Khosravi Mashizie, R.; Kakavandi, B.; Zare, M.R.; Asadi, A. Batch and column studies for the adsorption of chromium(VI) on low-cost Hibiscus Cannabinus kenaf, a green adsorbent. J. Taiwan Inst. Chem. Eng. 2016, 68, 80–89. [Google Scholar] [CrossRef]
  67. Rai, M.K.; Shahi, G.; Meena, V.; Meena, R.; Chakraborty, S.; Singh, R.S.; Rai, B.N. Removal of hexavalent chromium Cr (VI) using activated carbon prepared from mango kernel activated with H3PO4. Resour.-Effic. Technol. 2016, 2, S63–S70. [Google Scholar] [CrossRef]
  68. Mohan, D.; Singh, K.P.; Singh, V.K. Removal of Hexavalent Chromium from Aqueous Solution Using Low-Cost Activated Carbons Derived from Agricultural Waste Materials and Activated Carbon Fabric Cloth. Ind. Eng. Chem. Res. 2005, 44, 1027–1042. [Google Scholar] [CrossRef]
  69. Dima, J.B.; Sequeiros, C.; Zaritzky, N.E. Hexavalent chromium removal in contaminated water using reticulated chitosan micro/nanoparticles from seafood processing wastes. Chemosphere 2015, 141, 100–111. [Google Scholar] [CrossRef]
  70. Debnath, A.; Bera, A.; Chattopadhyay, K.K.; Saha, B. Facile additive-free synthesis of hematite nanoparticles for enhanced adsorption of hexavalent chromium from aqueous media: Kinetic, isotherm, and thermodynamic study. Inorg. Nano-Met. Chem. 2017, 47, 1605–1613. [Google Scholar] [CrossRef]
  71. Yang, R.; Aubrecht, K.B.; Ma, H.; Wang, R.; Grubbs, R.B.; Hsiao, B.S.; Chu, B. Thiol-modified cellulose nanofibrous composite membranes for chromium (VI) and lead (II) adsorption. Polymer 2014, 55, 1167–1176. [Google Scholar] [CrossRef]
  72. Shao, Y.; Fan, Z.; Zhong, M.; Xu, W.; He, C.; Zhang, Z. Polypyrrole/bacterial cellulose nanofiber composites for hexavalent chromium removal. Cellulose 2021, 28, 2229–2240. [Google Scholar] [CrossRef]
  73. Jain, P.; Varshney, S.; Srivastava, S. Synthetically modified nano-cellulose for the removal of chromium: A green nanotech perspective. IET Nanobiotechnol. 2017, 11, 45–51. [Google Scholar] [CrossRef]
  74. Sun, Y.; Yue, Q.; Mao, Y.; Gao, B.; Gao, Y.; Huang, L. Enhanced adsorption of chromium onto activated carbon by microwave-assisted H3PO4 mixed with Fe/Al/Mn activation. J. Hazard. Mater. 2014, 265, 191–200. [Google Scholar] [CrossRef] [PubMed]
  75. Acharya, J.; Sahu, J.N.; Sahoo, B.K.; Mohanty, C.R.; Meikap, B.C. Removal of chromium(VI) from wastewater by activated carbon developed from Tamarind wood activated with zinc chloride. Chem. Eng. J. 2009, 150, 25–39. [Google Scholar] [CrossRef]
  76. Liu, S.X.; Chen, X.; Chen, X.Y.; Liu, Z.F.; Wang, H.L. Activated carbon with excellent chromium(VI) adsorption performance prepared by acid–base surface modification. J. Hazard. Mater. 2007, 141, 315–319. [Google Scholar] [CrossRef] [PubMed]
  77. Zhang, D.; Xu, W.; Cai, J.; Cheng, S.-Y.; Ding, W.-P. Citric acid-incorporated cellulose nanofibrous mats as food materials-based biosorbent for removal of hexavalent chromium from aqueous solutions. Int. J. Biol. Macromol. 2020, 149, 459–466. [Google Scholar] [CrossRef]
  78. McCabe, W.L.; Smith, J.C.; Harriott, P. Unit Operations of Chemical Engineering; McGraw Hill: Boston, MA, USA, 2001. [Google Scholar]
  79. Ayawei, N.; Ebelegi, A.N.; Wankasi, D. Modelling and Interpretation of Adsorption Isotherms. J. Chem. 2017, 2017, 3039817. [Google Scholar] [CrossRef]
  80. Qiu, B.; Gu, H.; Yan, X.; Guo, J.; Wang, Y.; Sun, D.; Wang, Q.; Khan, M.; Zhang, X.; Weeks, B.L.; et al. Cellulose derived magnetic mesoporous carbon nanocomposites with enhanced hexavalent chromium removal. J. Mater. Chem. A 2014, 2, 17454–17462. [Google Scholar] [CrossRef]
  81. Rao, D.G.; Senthikumar, R.; Byrne, J.A.; Feroz, S. Wastewater Treatment: Advanced Processes and Technologies; CRC Press: Boca Raton, FL, USA, 2012. [Google Scholar]
  82. Dawodu, F.A.; Akpan, B.M.; Akpomie, K.G. Sequestered capture and desorption of hexavalent chromium from solution and textile wastewater onto low cost Heinsia crinita seed coat biomass. Appl. Water Sci. 2019, 10, 32. [Google Scholar] [CrossRef] [Green Version]
  83. De Andrade Neto, J.C.; Pereira, G.J.; Morandim-Giannetti, A.d.A. Chitosan and corn stover derivative bioadsorbent: Characterization and application in hexavalent chromium adsorption processes. Cellulose 2020, 27, 6317–6331. [Google Scholar] [CrossRef]
  84. Batool, F.; Akbar, J.; Iqbal, S.; Noreen, S.; Bukhari, S.N.A. Study of Isothermal, Kinetic, and Thermodynamic Parameters for Adsorption of Cadmium: An Overview of Linear and Nonlinear Approach and Error Analysis. Bioinorg. Chem. Appl. 2018, 2018, 3463724. [Google Scholar] [CrossRef]
  85. Choudhary, B.; Paul, D. Isotherms, kinetics and thermodynamics of hexavalent chromium removal using biochar. J. Environ. Chem. Eng. 2018, 6, 2335–2343. [Google Scholar] [CrossRef]
Figure 1. The evolution of the hexavalent chromium concentration [mg·L−1] at 0.1 mg·L−1 of the initial chromium concentration, pH 3, and 1000 mg CNF·L−1 of dosage during adsorption with different doses of MTMS in the CNF hydrogels.
Figure 1. The evolution of the hexavalent chromium concentration [mg·L−1] at 0.1 mg·L−1 of the initial chromium concentration, pH 3, and 1000 mg CNF·L−1 of dosage during adsorption with different doses of MTMS in the CNF hydrogels.
Polymers 14 03425 g001
Figure 2. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, pH 3, and 1000 mg CNF·L−1 hydrophobized with 1.5 mmol MTMS·g−1 cellulose and kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Figure 2. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, pH 3, and 1000 mg CNF·L−1 hydrophobized with 1.5 mmol MTMS·g−1 cellulose and kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Polymers 14 03425 g002
Figure 3. The evolution of the kinetic adsorption experiment at the previously indicated conditions under different MTMS dosages and kinetic fitting of the pseudo-second order model.
Figure 3. The evolution of the kinetic adsorption experiment at the previously indicated conditions under different MTMS dosages and kinetic fitting of the pseudo-second order model.
Polymers 14 03425 g003
Figure 4. The evolution of the hexavalent chromium concentration [mg·L−1] during adsorption with the CNF hydrogel at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization, and 1000 mg CNF·L−1 of dosage under different pH conditions.
Figure 4. The evolution of the hexavalent chromium concentration [mg·L−1] during adsorption with the CNF hydrogel at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization, and 1000 mg CNF·L−1 of dosage under different pH conditions.
Polymers 14 03425 g004
Figure 5. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization and 1000 mg CNF·L−1 of dosage under pH 3 conditions and the kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Figure 5. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization and 1000 mg CNF·L−1 of dosage under pH 3 conditions and the kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Polymers 14 03425 g005
Figure 6. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization, and 1000 mg CNF·L−1 of dosage under different pH and kinetic fitting of the pseudo-second order model.
Figure 6. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the initial chromium concentration, 1.5 mmol MTMS·g−1 CNF applied during hydrophobization, and 1000 mg CNF·L−1 of dosage under different pH and kinetic fitting of the pseudo-second order model.
Polymers 14 03425 g006
Figure 7. The evolution of the hexavalent chromium concentration [mg·L−1] during adsorption with the CNF hydrogel at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under different adsorbent dosages.
Figure 7. The evolution of the hexavalent chromium concentration [mg·L−1] during adsorption with the CNF hydrogel at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under different adsorbent dosages.
Polymers 14 03425 g007
Figure 8. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 500 mg CNF·L−1 of dosage and the kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Figure 8. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 500 mg CNF·L−1 of dosage and the kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Polymers 14 03425 g008
Figure 9. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under different adsorbent dosages and kinetic fitting of the pseudo-second order model.
Figure 9. The evolution of the kinetic adsorption experiment at 0.1 mg·L−1 of the chromium initial concentration, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under different adsorbent dosages and kinetic fitting of the pseudo-second order model.
Polymers 14 03425 g009
Figure 10. (a) The evolution of the hexavalent chromium concentration [mg·L−1] during the adsorption with CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 0.1 to 5 mg·L−1; (b) 10 to 50 mg·L−1 initial hexavalent chromium concentrations.
Figure 10. (a) The evolution of the hexavalent chromium concentration [mg·L−1] during the adsorption with CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 0.1 to 5 mg·L−1; (b) 10 to 50 mg·L−1 initial hexavalent chromium concentrations.
Polymers 14 03425 g010
Figure 11. The evolution of the kinetic adsorption experiment of the CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 25 mg·L−1 of the initial hexavalent chromium concentration and kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Figure 11. The evolution of the kinetic adsorption experiment of the CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 25 mg·L−1 of the initial hexavalent chromium concentration and kinetic fitting of the pseudo-first, pseudo-second, Elovich, and intraparticle models.
Polymers 14 03425 g011
Figure 12. (a) The evolution of the kinetic adsorption experiment of the CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 0.1 to 5 mg·L−1; (b) 10 to 50 mg·L−1 of the initial hexavalent chromium concentration and kinetic fitting of the pseudo-second order model.
Figure 12. (a) The evolution of the kinetic adsorption experiment of the CNF hydrogel at 1000 mg·L−1 of dosage, 1.5 mmol MTMS·g−1 CNF, and pH 3 conditions under 0.1 to 5 mg·L−1; (b) 10 to 50 mg·L−1 of the initial hexavalent chromium concentration and kinetic fitting of the pseudo-second order model.
Polymers 14 03425 g012
Table 1. The nonlinear and linearized equations of the analyzed kinetic models.
Table 1. The nonlinear and linearized equations of the analyzed kinetic models.
ModelNonlinearized Equations Linearization Ref.
Pseudo-first order 1 q = q e · 1 e k 1 · t (4) ln q e q ln q e = k 1 · t (8)[43]
Pseudo-second order q = q e 2 · k 2 · t 1 + q e · k 2 · t (5) t q = 1 q e 2 · k 2 + 1 q e · t (9)[43]
1 q = 1 q e 2 · k 2 · 1 t + 1 q e (10)[44]
Elovich q = 1 β · ln t + 1 β ln α · β (6)- [45]
Weber and Morris (Intraparticle) q = k i · t 0.5 + C (7)- [46]
1 The value of qe must be previously obtained by estimation or experimentally.
Table 2. The nonlinear and linearized equations of the analyzed isotherm models.
Table 2. The nonlinear and linearized equations of the analyzed isotherm models.
ModelNonlinearized Equations Linearization Ref.
Langmuir q e = K L · q max · C e 1 + K L · C e (11) Type   I :   C e q e = 1 K L · q max + 1 q max · C e (20)[47]
Type   II :   1 q e = 1 q max + 1 K L · q max · 1 C e (21)
Type   III :   q e = 1 K L · q e C e + q max (22)
Type   IV :   q e C e = K L · q e + K L · q max (23)
R L = 1 1 + K L · C 0 (12) Type   V :   1 C e = K L · q max · 1 q e K L (24)
Freundlich q e = K F · C e 1 / n F (13) ln q e = 1 n F · ln C e + ln K F (25)
Temkin q e = B T · ln A T · C e (14) q e = B T · ln A T + B T · ln C e (26)
B T = R · T b T (15)
Dubinin–Raduskevich q e = q max · exp B DR · ε 2 (16) ln q e = ln q max B DR · ε 2 (27)
ε = R · T · ln 1 + 1 C e (17)
E = 1 2 · B DR (18)
Sips 1 q e = K S · q max · C e 1 / n S 1 + K S · C e 1 / n S (19) ln q e q max q e = 1 n s · ln C e + ln K S (28)
1 The value of qmax can be first estimated from the qmax obtained from the Langmuir model as the first input for optimization using a calculation software.
Table 3. A comparison of the hexavalent chromium adsorption through different NC and activated carbon adsorbents.
Table 3. A comparison of the hexavalent chromium adsorption through different NC and activated carbon adsorbents.
AdsorbentContact Time [min]Adsorbent Dosage [mg·L−1]Initial Cr(VI) Concentration [mg·L−1]pHqmax [mg·g−1]Maximum Removal Yield [%]Ref.
CNF from rice husk100 15003063.7692.99[57]
Polypyrrole-bacterial CNF180 2503002555.697.5[72]
Thiol-modified CNF composite20 50487.596[71]
Citric acid-incorporated CNF 120 405021123[77]
Amino-silanized cellulose membranes 300500050434.7 [33]
Polyaniline-functionalized CNC 40500302.548.9297.84[73]
Microwave-assisted H3PO4/Fe-modified activated carbon200100030334.39100[74]
ZnCl2-modified tamarind wood activated carbon70300010328.0299[75]
Acid-base surface modified activated carbon180200050 13.89 [76]
Hydrophobized CNF Hydrogel (MTMS dosage = 1.5 mmol·g−1)330 50050370.38>97.14This work
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ojembarrena, F.d.B.; Sánchez-Salvador, J.L.; Mateo, S.; Balea, A.; Blanco, A.; Merayo, N.; Negro, C. Modeling of Hexavalent Chromium Removal with Hydrophobically Modified Cellulose Nanofibers. Polymers 2022, 14, 3425. https://doi.org/10.3390/polym14163425

AMA Style

Ojembarrena FdB, Sánchez-Salvador JL, Mateo S, Balea A, Blanco A, Merayo N, Negro C. Modeling of Hexavalent Chromium Removal with Hydrophobically Modified Cellulose Nanofibers. Polymers. 2022; 14(16):3425. https://doi.org/10.3390/polym14163425

Chicago/Turabian Style

Ojembarrena, Francisco de Borja, Jose Luis Sánchez-Salvador, Sergio Mateo, Ana Balea, Angeles Blanco, Noemí Merayo, and Carlos Negro. 2022. "Modeling of Hexavalent Chromium Removal with Hydrophobically Modified Cellulose Nanofibers" Polymers 14, no. 16: 3425. https://doi.org/10.3390/polym14163425

APA Style

Ojembarrena, F. d. B., Sánchez-Salvador, J. L., Mateo, S., Balea, A., Blanco, A., Merayo, N., & Negro, C. (2022). Modeling of Hexavalent Chromium Removal with Hydrophobically Modified Cellulose Nanofibers. Polymers, 14(16), 3425. https://doi.org/10.3390/polym14163425

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop