Research Article  Open Access
Synthesis of Iron Doped Zeolite Imidazolate Framework8 and Its Remazol Deep Black RGB Dye Adsorption Ability
Abstract
Zeolite imidazole framework8 (ZIF8) and the iron doped ZIF8 (FeZIF8) were synthesized by the hydrothermal process. The obtained materials were characteristic of Xray diffraction (XRD), Xray photoelectron spectroscopy (XPS), scanning electron microscope (SEM), nitrogen adsorption/desorption isotherms, and atomic absorption spectroscopy (AAS). The results showed that the obtained FeZIF8 possessed the ZIF8 structure with a large specific area. ZIF8 and FeZIF8 were used for the removal of Remazol Deep Black (RDB) RGB dye from aqueous solutions. The various factors affecting adsorption such as pH, initial concentration, contact time, and temperature were investigated. The results showed that the introduction of iron into ZIF8 provided a much larger adsorption capacity and faster adsorption kinetics than ZIF8 without iron. The electrostatic interaction and ππ interaction between the aromatic rings of the RDB dye and the aromatic imidazolate rings of the adsorbent were responsible for the RDB adsorption. Moreover, the coordination of the nitrogen atoms and oxygen in carboxyl group in RDB molecules with the Fe^{2+} ions in the ZIF8 framework played a vital role for the effective removal of RDB from aqueous solution.
1. Introduction
It is wellknown that textile industries, pulp mills, and dyestuff manufacturing discharge a considerable amount of colored wastewaters, which has provoked serious environmental concerns all over the world [1]. Its removal is therefore of prime importance. Owing to their complicated chemical structures, dyes are difficult to treat with municipal waste treatment operations. There are many treatment processes such as chemical precipitation, membrane filtration, and alum coagulation of which adsorption is considered the most effective method widely employed to treat wastewater containing different classes of dyes.
Nearly 60% of the dyes used in textile plants are azo dyes which are classified as mono, di, tri, and tetraazo dyes, with the azo groups mainly bound to the benzene and naphthalene rings. Remazol Deep Black (denoted as RDB) RGB is a common diazo reactive dye in aqueous solutions and widely used in textile industries [2]. It is stable and hardly biologically degradable due to the presence of aromatic rings. Thus, much research attention has been paid to how to eliminate RDB from aqueous solutions. Various processes including adsorption, electrochemistry, and biosorption for RDB treatment have been reported. Soloman et al. [2] studied the electrochemical degradation of hydrolyzed Remazol Black. Performance of the batch recirculation system was comparatively better than the other rector configurations studied, with respect to capacity utilization and energy consumption. Cardoso et al. [3] used Brazilian pinefruit shells (Araucaria angustifolia) in natural form and as adsorbents for the removal of RDB textile dye from aqueous effluents, with a maximum sorption capacity of 74.6 mg g^{−1}. Biosorption of the azo dye by growing fungi (Aspergillus flavus) was investigated in batch reactors. The nearly complete removal of dye was found at initial concentration up to 250 mg/L and at pH 4.5 which was used as working pH value for removal of dye in all the batch studies. The removal of chemical oxygen demand (COD) was found to be 90% at 100 mg/L initial concentration of dye [4].
Metalorganic frameworks (MOFs) are hybrid materials with ordered threedimensional frameworks via strong metalligand bonds between metal cations and organic linkers [5]. Since their discoveries, MOFs have received significant attention in their potential applications in gas storage [6–13], separation [14, 15], and heterogeneous catalysis [16–19]. The zeolitic imidazolate framework (ZIF) materials have zeolite like topologies and belong to an important class of MOFs materials with interesting adsorption, separation, and catalytic properties [20–23]. Among them, ZIF8 [Zn(2methylimidazole)_{2}·2H_{2}O] constructed from 2methylimidazole ligands and Zn(II) center ions exhibits higher thermal and chemical stability than other MOFs [24]. Many studies show that ZIF8 exhibited efficient removal of arsenic and organic pollutants from aqueous solutions. Jiang et al. [25] reported that ZIF8 as a novel adsorbent for fast removal of 1Hbenzotriazole and 5tolyltriazole with regard to adsorption isotherms, kinetics, thermodynamics, desorption, and adsorbent regeneration. Lin and Chang [26] reported that the adsorption capacity of ZIF8 was much higher than that of fly ash, activated carbon, zeolites, and so forth, showing its promising potential for the removal of humic acid. PArsanilic acid which is widely used as feeding additive in the poultry and pork industries to enhance the feeding efficiency was efficiently removed by using ZIF8 as an adsorbent [27]. ZIF8 exhibited the high arsenate/arsenite adsorption (up to 50 mg·g^{−1} for As(II) and 60 mg·g^{−1} for As(V)) [28]. Zheng et al. [29] reported a free solvent synthesis of coreshell Fe_{3}O_{4}@zeolitic imidazolate frameworks8 (ZIF8) via two steps. The introduction of magnetic iron oxide into ZIF8 facilitated the separation of adsorbents by the magnetic field. Fe_{3}O_{4}@ZIF8 showed good adsorption properties for methylene blue, with a maximum adsorption capacity of 20.2 mg g^{−1}. According to these studies, ZIF8 not only is highly stable in water but also exhibits promisingly high adsorption capacities. Owing to these features, ZIF8 should be promising and feasible adsorbents to organic pollutants from aqueous solutions. On the other hand, in addition to being a catalyst, iron species can act as adsorptive sites for adsorption processes. The combination of iron and ZIF8 (denoted as FeZIF8) is expected to provide a novel adsorbent due to large accessible surface area and abundant active surface sites.
In this study, the introduction of iron into ZIF8 by a onestep process was performed with the aim of improving its adsorption ability. FeZIF8 was used as an adsorbent for removing RDB dye. In addition, the dye adsorption over ZIF8 was performed for comparison. The iron oxide incorporated into the ZIF8 framework significantly enhanced a RDB adsorption capacity compared to bare ZIF8. A possible adsorption mechanism was suggested, based on adsorption at various pH values and surface charges on FeZIF8.
2. Experimental
2.1. Materials
Zinc nitrate hexahydrate (Zn(NO_{3})_{2}·6H_{2}O, Daejung, Korea, ≥99%), iron(II) sulfate heptahydrate (FeSO_{4}·7H_{2}O, Merck, Germany > 99%), methanol (CH_{3}OH, Merck, Germany), and 2methylimidazole (C_{4}H_{6}N_{2}, SigmaAldrich, USA, 99%) were utilized in this paper. Remazol black B RGB (C_{26}H_{21}N_{5}Na_{4}O_{19}S_{6}, molecular weight = 991.82) was obtained from Thuy Duong Textile Company, Vietnam. The structure of RDB is shown in Scheme 1.
2.2. Preparation of ZIF8 and Iron Doped ZIF8 (FeZIF8)
ZIF8 and FeZIF8 were synthesized as in [30, 31]. Briefly, 2.8 mmol of zinc(II) and iron(II) (molar ratio of Fe/Zn = 0/10 or 1/9) were dissolved in 1.4 mmol of methanol. A solution consisting of 64.4 mmol of 2methylimidazole and 1.4 mol of methanol was added to the ZnFe based solution and vigorously stirred for 24 hs at ambient temperature. Nitrogen was bubbled through the solution to minimize the oxidation reaction of Fe(II) to Fe(III) species. Finally, this solution was centrifuged at 300 rpm and washed thoroughly with methanol. This washing procedure was repeated 3 times. The resultant crystals were dried overnight at 120°C. The obtained samples with the molar ratio of Fe(II)/(Zn(II) + Fe(II)) being 0/10 and 1/10 were named ZIF8 and FeZIF8, respectively. In the obtained samples, ZIF8 was white and FeZIF8 was light brown.
2.3. Determination of the Point of Zero Charge
The pH at the potential of zero charge () of ZIF8 and FeZIF8 was measured by the pH drift method [32]. To a series of 100 mL flasks, 5 mL of 0.1 M NaCl solution and 40 mL of distilled water were added. The initial pH value (pHi) of the solution was adjusted from 2.3 to 12.0 by adding either 0.1 M NaOH or 0.1 M HCl. The total volume of solution in each flask was made exactly as 50 mL by adding distilled water. The 0.01 M NaCl solutions with different pH values were obtained. Nitrogen was bubbled through the solution to eliminate the dissolved CO_{2}. Then, 0.02 grams of the sample was added to each flask and mixtures were sealed and shaken for 24 hours; then the final pH () of solution was recorded. The plot of versus initial pH_{i} was conducted. The point of intersection of curve with abscissa, at which ΔpH = 0, provided .
2.4. Adsorption Kinetics Study
Experiments were conducted in a batch process. The 3 L plastic beaker was equipped with a stainless steel flat blade impeller using an electric motor to stir the dye solution. Samples (ZIF8 or FeZIF8) (0.4 g) were vigorously mixed with 1000 mL of RDB solution in the beaker at a fixed temperature. Ten millilitres of solution was drawn at preset intervals of time through tap and the solid was removed by centrifuging process. The residual dye concentrations were determined using UVVis spectrophotometry. The experiments were conducted with various RDB concentrations ranging from 30 to 50 ppm. The amount of the dye adsorbed by the adsorbent was calculated by the following equation:where is the amount of dye adsorbed per unit of adsorbent amount (mg g^{−1}) at time, the initial dye concentration (mg L^{−1}), the dye concentration (mg L^{−1}) after the batch adsorption procedure, the volume of dye solution (L), and the mass (g) of the adsorbent. It is reasonably assumed that the adsorption of dye from solution follows reversible firstorder kinetics. The heterogeneous equilibrium between dye in solution and solid adsorbent is illustrated as follows [32]:where and are the forward and backward rate constants, respectively.
The equilibrium constant, , defined as could be expressed [26]where and are dye concentrations (mg L^{−1}) at initial and equilibrium time, respectively, and the others are described above.
Pseudofirstorder kinetic of adsorption was investigated by NatarajanKhalaf equation [33]:where is the rate constant of adsorption process.
The slope of the linear plot of versus will provide the value of .
Based on the relaxation method [34, 35], could be expressed asThen, the forward and backward rate constants could be derived from (3) and (5).
The kinetics of diffusion was studied by Webber’s intraparticlediffusion model. Webber’s intraparticlediffusion model is described in the following equation [36, 37]:where is intraparticlediffusion rate constant (mg g^{−1} min^{−0.5}) and the intercept which reflects the layer boundary effect.
The analysis of the multilinearity in pore and filmdiffusion plot using Webber’s plot was conducted by using piecewise linear regression proposed by Malash and ElKhaiary [38].
In this method, the experimental data could be fixed for one, two, or three linear segments’ line by Webber’s model. One linear segment’s line: (two parameters). Two linear segments’ line: (four parameters). Three linear segments’ line: (six parameters),where the values of , , , , , and are estimated by nonlinear regression. and called breakpoints are the boundaries between the segments. The Microsoft Excel “sign” function is defined as follows:The example for the two linear segments’ equation was expressed as follows:Then the linear equation of the first segment is , where and .
Then the linear equation of the second segment is , where and .
Nonlinear regression determines the model’s parameters by the least squares method. This is calculated by minimizing the sum of squared deviations, , by numerical optimization techniques using Solver function in Microsoft Excel. The function for minimization iswhere is experimental datum and is the value estimated by model.
The determination coefficient, , is obtained by the expressionwhere is the total sum of squares equal to ( is the mean value of ).
The comparison of models was based on Akaike’s Information Criterion (AIC) [38–40]. The determines how well the data support each model. The value of AIC can be positive or negative. The model with the lowest AICs score is most likely correct. The (for a small size sample) is calculated for each model from the following equations:where is the number of experimental points and is the parameter sum of model.
2.5. Thermodynamic and Isothermal Studies
Experimental procedure was conducted as an adsorption kinetics study. However, the temperature of the process was fixed at 298, 308, and 318 K. The activation energy, , was determined by Arrhenius equation [35]:where is the rate constant equal to the rate constant in Natarajan and Khalaf equation, the frequency factor, gas constant (8.315 J mol^{−1} K^{−1}), and absolute temperature in Kelvin.
Taking the natural logarithm of both sides of (12), one obtainsBy linear plotting ln versus , the could be obtained from slope ().
Thermodynamic parameters of activation can inform whether or not the adsorption process follows an activated complex or is prior to the final adsorption. Thermodynamic parameters of activation including the enthalpy (, entropy , and free energy Gibbs of activation for RBB adsorption kinetics were obtained by applying Eyring equation [41, 42].
The Eyring equation in its thermodynamic version is as follows: where is the rate constant equal to the rate constant in NatarajanKhalaf equation, the (1.3807 × 10^{−23} J K^{−1}) the Boltzmann constant, and (6.621 × 10^{−34} J s) the Planck constant.
Taking the natural logarithm of both sides of (14), Eyring equation in linear form was obtained:By linear plotting versus , and were obtained from the slope ( and intercept .
The Gibbs free energy of activation can be obtained by
In order to assess if the adsorption process is spontaneous or not, the thermodynamic parameters of adsorption are needed. The standard Gibbs free energy of adsorption () is given by the expression [35, 43]where , , and are the standard Gibbs free energy, enthalpy, and entropy, respectively.
is given by van’t Hoff’s equation:where is the distribution coefficient of the solute ions and equals [27, 44, 45], and the others are described above.
By replacing (18) with (17), one obtainsThe value of and was determined from the slope and intercept of the linear plot of versus .
The adsorption isotherms were developed at 25°C. The adsorption time was conducted for 24 hours to confirm saturation. Thereafter supernatant liquid was collected by centrifugation and the final dye concentrations were determined using UVVis spectrophotometry.
The experimental data were analyzed according to the Freundlich and Langmuir models.
Langmuir Isotherm. The Langmuir equation is valid for monolayer sorption onto the surface. It could be expressed as follows [37, 46]:where is the maximum monolayer capacity amount (mg g^{−1}), is Langmuir equilibrium constant (L mg^{−1}), and the others are described above.
The essential characteristics of the Langmuir isotherm can be expressed in terms of a dimensionless constant separation factor, , which is performed aswhere the value of indicates the type of isotherm: unfavorable , linear , favorable , or irreversible [47].
Freundlich Isotherm. Freundlich equation is an empirical relation based on the adsorption of adsorbates onto the heterogeneous surface. It is represented as follows [48, 49]:where is the Freundlich constant, which is a measure of adsorption capacity and an empirical parameter related to the nature and strength of the adsorption process. A large value of means that the surface is heterogeneous. For values in the range , adsorption is favorable. Values of between 2 and 10 represent good adsorption processes, whereas indicates that adsorption capacity is only slightly suppressed at lower equilibrium concentrations [50].
The parameters of models ((20) and (22)) were estimated by nonlinear regression method using Solver function in Microsoft Excel. To quantitatively compare the applicability of each model, apart from the regression coefficient , the Chisquare test was calculated as follows [51]:where and are the adsorption capacity at the equilibrium experimental condition and adsorption capacity estimated by model, respectively.
A small value of indicates that the data obtained from the model is consistent with the experimental value.
2.6. Characterization of Materials
The powder Xray diffraction (XRD) patterns were recorded by a D8 Advance, Bruker (Germany) with CuK radiation ( Å). The morphology of the obtained sample was determined by scanning electron microscope (SEM) using SEM JMS5300LV (Japan). The specific surface area of the sample was determined by nitrogen adsorption/desorption isotherms using a Micromeritics 2020 volumetric adsorption analyzer system (USA). Thermal behaviors of the obtained materials were analyzed by means of thermal analysis (TGDTA) using Labsys TG Setaram (France). The element analysis was conducted by atomic absorption spectrometry (AAS) using AA6800 Shimazu (Japan). Visible spectrophotometry was measured by Lambda 25 Spectrophotometer, PerkinElmer (Singapore) at of RDB dye (600 nm).
3. Results and Discussion
3.1. Characterization of ZIF8 and FeZIF8
Figure 1 shows XRD patterns of ZIF8 and Fe ZIF8. The XRD pattern of ZIF8 in this work agreed well with patterns from [30, 31]. There was a welldefined diffraction (011) at two theta = 7.16° in the XRD pattern of ZIF8, indicating that the crystallinity of ZIF8 in this work was relatively high. The XRD patterns of FeZIF8 also exhibited characteristic peaks of ZIF8 and no characteristic peaks of iron oxides were observed. However, intensity of these diffractions decreases as a certain amount of iron was doped into the ZIF8 framework.
SEM images of ZIF8 and FeZIF8 are presented in Figure 2. Morphology of ZIF8 consisted of spherical particles with a diameter around 100 nm while the irregular shapes of FeZIF8 with sizes around 0.1–100 nm were observed. The introduction of iron into ZIF8 caused a significant change in morphology of ZIF8.
(a)
(b)
The zinc and iron compositions were analyzed by AAS. The results are presented in Table 1. Iron percentage in the final product (0.116) was higher than the original (0.100). The pH of the synthesized gel was around 4.5. Then iron(II) was likely to incorporate completely into ZIF8; however, possibly Zn(II) was partly dissolved in the solution. This is reason why there was an increase in the percentage of iron in the final product.

The XPS spectra indicated a chemical state of element, that is, iron (Fe2p) and zinc (Zn2p). The peak of Zn 2p1/2 (1044 eV) and Zn2p3/2 (1020.96 eV) observed for both samples confirmed the existence of Zn(II) (Figure 3). For ZIF8, the peak of Fe 2p3/2 was inconspicuous, indicating that iron was a very minor component (in fact, it could not be detected). Only peak Fe2p3/2 for Fe(II) at 709.98 eV was detected implying the main iron in FeZIF8 was Fe(II). The percentage of oxidation state of iron, calculated from peak areas, was listed in Table 1. It was worth noting that the initial iron source of Fe(III) was also tested to incorporate into ZIF8 but the solid product was not obtained. This means that the presence of Fe(III) in the initial synthesized mixture was not favorable for the formation of ZIF8 structure.
(a)
(b)
Based on ZIF8 with space group of [52] the cell parameter of ZIF8 and FeZIF8 was expressed bywhere is spacing distance, is cell parameter, and , , are Miller indexes of diffraction planes.
The cell parameter of ZIF8 (16.800 Å) and FeZIF8 (16.977 Å) seems to be unchangeable. Since the charge and radii of Zn^{2+} (0.75 Å) and Fe^{2+} (0.74 Å) are similar, it is likely that Fe(II) can substitute Zn(II) in ZIF8 or disperse highly in ferrous forms.
Figure 4 shows the nitrogen adsorption/desorption isotherms of ZIF8 and FeZIF8. All samples exhibited type IV with H4 which is characteristic of mesoporous materials. FeZIF8 possessed a shape which is different from ZIF8 at high relative pressure. This result suggests that the porous structure was distorted due to the incorporation of iron oxides. ZIF8 exhibited a high specific surface area of 1380 m^{2}·g^{−1} (calculated by BET model), which was similar to that found in the previous literature [30, 31, 53].
The introduction of iron oxide into ZIF8 lowered the specific surface area, the pore diameter, and the pore volume. The specific surface areas are 1380 and 1243 m^{2}·g^{1} for ZIF8 and FeZIF8, respectively (see Table 2). This also gave evidence of the encapsulation of iron oxides within the pores of framework, which brought about the lowering of accessible void space for N_{2} gas molecules.
 
: specific surface area calculated by BET model; : specific surface area calculated by Langmuir model; : pore diameter calculated by BJH model; : pore volume. 
The thermal stability of materials was tested by TGDTA. ZIF8 and FeZIF8 were found to be highly stable up to 220°C and 250°C. Beyond this temperature the framework slowly started to decompose and a flat valley was obtained till 700°C (Figure 5). The incorporation of iron in ZIF8 seems to make the materials more stable. This behavior was also observed as TiO_{2} is doped in ZIF8 [54].
The stability of ZIF8 in water with different pHs was also studied. Figure 6 presents XRD patterns of ZIF8 which were soaked in water with pH from 2 to 12. The pH of the solution was adjusted by NaOH 0.01 M or HCl 0.01 M. The crystallinity of ZIF8 nanoparticles was retained well at pH > 3–12, which proved that ZIF8 was stable in aqueous solutions with pHs in the range of 3–12.
3.2. A Study on RDB Adsorption onto ZIF8 and FeZIF8
3.2.1. Effect of Initial RDB Concentration
The experimental results for adsorption at various concentrations (30–50 mg·L^{−1}) with contact time are shown in Figure 7. As seen from Figure 7, the adsorption capacity of adsorbent increases as initial dye concentrations go up. The RDB adsorption of FeZIF8 was higher than that of ZIF8 in the same initial concentration. The adsorption capacity of RDB onto ZIF8 increased from 30.20 to 42.14 mg g^{−1} and that onto FeZIF8 increased from 50.36 to 76.79 mg g^{−1} as the initial concentrations increased from 30 to 50 mg g^{−1}. The initial dye concentration provided a critical driving force to overcome all the mass transfer resistance of RDB between the aqueous solution and the adsorbent surface [55]. As a result, high initial RDB concentration might promote the adsorption process. Figure 7 indicates that the adsorption of RDB was fast in the earlier stage (0–100 minutes) and gradually reached the equilibrium. It is found that the time required to obtain the adsorption equilibrium was around 250 minutes.
Due to the porosity of ZIF8 and FeZIF8, intraparticle diffusion was expected in the adsorption process. This was examined by plotting RDB uptake, against the square root of time, as (6). If intraparticle diffusion is the ratelimiting step this plot will give a straight line and an intercept equal to zero. However, previous studies by various researchers showed that these plots represent multilinearity [56]. This application often suffers uncertainties caused by the multilinear nature of their plot. Therefore, the confidence intervals for parameters are determined to estimate their uncertainties. If the 95% confidence interval of intercept does not contain zero or it varies from this negative/positive value to another, intercept is significantly different from zero. It is concluded that the line does not pass through the origin or vice versa (in this study the SPSSversion 21 was used to calculate the confidence interval).
Piecewise linear regression proposed by Malash and ElKhaiary [38] was applied to analyze multilinearity of Webber’s plot. Since increasing the number of linear segments leads to an increase in the number of regression parameters, a decrease in naturally follows. For this reason, solely cannot be used to assess the goodness of fit for models. The wellknown statistical method for model comparison is Akaike’s Information Criterion (AIC) [38, 57]. This method advises which model is more likely to be correct. The smaller AIC value indicates a more compatible model. For example, Figure 8 illustrates experimental data and piecewise linear regression lines with initial concentration 50 mg L^{−1} for ZIF8 and FeZIF8. The experimental points seem to be close to regression lines for two or three linear segment lines. We could not estimate visually which one is more likely correct. The values of for RDB adsorption onto ZIF8 are 89.47 and 16.05 and 22.46 for onesegment, twosegment, and threesegment models, respectively. Those for RDB adsorption onto FeZIF8 are 98.65 and 6.67 and 10.01 for onesegment, twosegment, and threesegment models, respectively. It is decided to accept the model with two linear segments as the correct model because of the lowest value of in this model.
A comparison of models based on is presented in Table 3. Table 3 shows that Webber’s model analyzed two segments linear regression provided the lowest compared to one or threesegment models. In conclusion, the experimental data best fit with the two segments’ linear regression model. Results of piecewise two linear segments’ regression for different initial concentrations are shown in Table 4. For illustration, for 50 mg L^{−1} concentration of ZIF8 in Figure 8, the intercepts of the first line in the Webber plot were −10.19 with 95% confidence interval (−12.28 to −8.10). This value of the intercept was significantly different from zero. It means the line did not pass through the origin. The similar behaviors were observed for all the other cases. These results indicate that the adsorption of RDB dye onto ZIF8 or FeZIF8 in the first stage was controlled by film diffusion (e.g., surface adsorption and liquid film diffusion) while the second was assigned to intraparticle diffusion [57–59]. The breakpoints (see Table 4) reflect the adsorption process moving from one mechanism to another. The times of phase transition calculated by the square of breakpoints varied from 36 to 100 minutes. The adsorption could be divided into two distinct phases by the time of phase transition: (i) main adsorption of RDB molecules within 0–100 minutes of the contact times: this process was rather slow compared to other adsorbents [40, 57] in which the fast and instantaneous adsorption was observed; since the pore sizes of FeZIF8 or ZIF8 are around 1 nm [45], the adsorption of RDB with four benzene rings in which each has a critical dimension of around 0.662 nm [60] on ZIF8 limited the diffusion of the RDB molecule to the inner pore structure; however, abundant active adsorption sites could be obtained at its external surface, which is helpful to the adsorption of dye molecules; (ii) a gradual attainment of the equilibrium where only about 5–10% of the adsorption was encountered. This is due to the utilization of the all active sites on the adsorbent surface. The first linear segment indicated a boundary layer effect while the second linear segment was assigned to intraparticle diffusion [58, 59].


The rate parameter, , increases with an increase in initial concentrations. However, the effect of initial RDB on rate parameters was irregular for RDB adsorption onto FeZIF8. The irregularity of could be related to that fact that iron incorporated caused the heterogeneity of ZIF8 structure.
The experimental data usually exhibit the multilinear plots, two or more stages influencing the adsorption process. In the statistical method, the piecewise linear regression is relevant to analyzing the data with multilinearity instead of the conventional graphical analysis [61] that might misidentify objective break points or numbers of segments.
As seen from Figure 9, a possible desorption might occur where the RDB uptake appears to fluctuate or even drop a little. This behavior could be attributed to either a reversible adsorption or a back diffusion controlling mechanism [62]. The pseudofirst or secondorder kinetic models of Lagergren [33] or Ho and McKay [63], respectively, have been widely used to investigate the formal kinetics of adsorption processes [12, 27, 40, 64]. However, the present experimental data could not apply to these models due to reversible adsorption. In the present study, the pseudofirstorder kinetics model of NatarajanKhalaf was used to analyze the experimental data. The rate constant, , for the adsorption of dye was determined from the slope of NatarajanKhalaf’s plots. The results are listed in Table 5. The high coefficient of determination (0.973–0.998) implies that this model was compatible with the experimental data. The rate constant of adsorption is separated into the rate of forward and reverse adsorption. The rate constants for the forward and reverse process are also presented in Table 5. It is clear that the adsorption kinetics can be significantly improved by the introduction of iron into the ZIF8. As shown in Table 5 the rate constants for adsorption could be increased approximately 5 times by FeZIF8 and the kinetics with FeZIF8 were faster than that with ZIF8.

(a)
(b)
Adsorption thermodynamics was conducted by varying the temperature from 298 K to 318 K as shown in Figure 9.
The thermodynamic parameters including activation energy, , , and are presented in Table 6. The results showed that equilibrium adsorption capacity, , of both adsorbents increased with an increase in temperature which indicated that the process was endothermic. The equilibrium adsorption capacity of FeZIF8 is higher than that of ZIF8 for each corresponding temperature. The increasing rate constant with an increase in temperature suggests that temperature favors the adsorption process for the fast rate to proceed. The equilibrium constant, , for FeZIF8 and ZIF8 changed in the range of 1.3–1.8 and 3.4–5.7, respectively, as temperature varied from 298 to 318 K. It is worth noting that in the case of FeZIF8 is higher and increases much faster than that in the case of ZIF8. The activation energy refers to the minimum amount of energy that must be overcome to proceed with the adsorption. for ZIF8 and FeZIF8 was 48.27 kJ and 12.51 kJ, respectively. Lower for FeZIF8 indicated that FeZIF8 was more favorable for RDB adsorption than ZIF8. Low activation energy (below 42 kJ mol^{−1}) implies diffusion controlled process because the temperature dependence of pore diffusivity is relatively weak and the diffusion process refers to the movement of the solute to an external surface and not diffusivity of material along microspore surfaces in a particle [65]. with ZIF8 is slightly higher than 42 kJ mol^{−1}. However, a reversible adsorption was observed as shown in Figure 9, indicating that the ratelimiting step in this case involved a physicalchemical mechanism and not purely physical or chemical one. Therefore the RDB adsorption onto ZIF8 was also considered a diffusion controlled process and so is the RDB adsorption onto FeZIF8.

The activation parameters calculated using Eyring equation are presented in Table 7. This would express whether the adsorption process follows an active complex prior to the final sorption. The coefficient of determination was high for both ZIF8 and FeZIF8 although the Eyring equation is not strictly linear with respect to [66]. The negative value of (−29.44 J K^{−1} mol^{−1} for ZIF8 and −153.40 J K^{−1}·mol^{−1} for FeZIF8) suggests a possibility of an associative chemisorption through the formation of an activated complex between RDB molecule and adsorbent [47], [63]. Also the negative value of normally reflects that no significant change occurs in the internal structure of the adsorbent during the adsorption process [41, 66]. The values for (52. 03 kJ.mol^{−1} for ZIF8 and 9.94 for FeZIF8) suggest that these reactions are endothermic. The large, positive in both RDB adsorptions onto ZIF8 and FeZIF8 imply that these reactions require energy to convert reactants to the product and as the energy requirement is satisfied the reaction proceeds. Typically, value relates to the adsorption rate. The rate increases as decreases [41, 42]. This is seen when comparing the data from Tables 6 and 7. In Table 6, the rate constant, , of FeZIF8 is higher than that of ZIF8. Table 7 describes the trend for in which the FeZIF8 has the lower value than ZIF8.

The thermodynamic parameters, , and , of system were determined using van’t Hoff equation to assess the spontaneity of adsorption process. In Table 8, the results show that the adsorption process using ZIF8 and FeZIF8 was endothermic as indicated by the positive sign of the value. The positive value of indicates the increasing randomness at the solidliquid interface during the adsorption of RDB molecules on the adsorbent [43]. The negative values of suggest the spontaneous RDB adsorption of RDB on ZIF8 or FeZIF8. The more negative value at higher temperatures implies that the spontaneity increases with a temperature increase. As the Gibbs free energy change is negative and accompanied by the positive standard entropy change, the adsorption reaction is spontaneous with high affinity. Consistent with the findings in the kinetics, more negative values of with FeZIF8 were obtained compared to that with ZIF8. Again, this confirmed the important role of iron in ZIF8 in enhancing the RDB capacity.

3.2.2. Effect of pH
One of the important parameters controlling the adsorption process is pH. Figure 10 shows the effect of pH on the removal of RBB dye from aqueous solutions. The behaviors of pH effect for RDB adsorption onto ZIF8 and FeZIF8 were similar. However, the RDB adsorption capacity of FeZIF8 was higher than that of ZIF8. In general, the RDB adsorption capacity of adsorbents was observed to increase significantly with an increase in pH from 2.2 to 6, followed by a slight increase in pH from 6 to 10 and it decreased significantly in further pH increase. The values of for ZIF8 and FeZIF8 determined by pH drift method were around 9.2 and 9.8, respectively, and are in accordance with previous results [25, 27, 67] (the inset in Figure 10). The value of indicates that the surface of the adsorbent is positively charged when pH of the solution is below , while the surface of adsorbent becomes negatively charged at pH of the solution above .
As seen in Scheme 1, molecular structure of RDB consists of two groups of – and –, a group of NH_{2}. Commonly, equilibrium of RDB in water is expressed as in Scheme 2.
Therefore, the higher pH is the more the equilibrium of (22) shifts to the right. With regard to the correlation of the electronic charges of the adsorbent/adsorbate and solution pH values, it can be presumed that there might be an electrostatic interaction between the positively charged ZIF8 surface and the negatively charged site of RDB increase as pH of solution increases up to . The adsorption capacity of ZIF8 and FeZIF8 started to decrease significantly at pH > ~10.2, which was probably due to the electrostatic repulsion of negatively charged RDB and negative ZIF8. In addition, the low stability of ZIF8 framework in pH lower than 3 also contributes to the low RDB adsorption capacity (see Figure 6). This electrostatic interaction mechanism is similar to the adsorption mechanism for phthalic acid [67] and arsanilic acid on ZIF8 [27]. In addition, the hydrophobic and ππ interaction between the aromatic rings of the RDB and the aromatic imidazole rings of the ZIF8 framework is also thought to contribute to the RDB adsorption capacity.
The increasing RDB uptake of FeZIF8 is probably due to the adsorption through other mechanisms addition to the electrostatic interaction. The enhanced adsorbed amount with FeZIF8 probably resulted from an increased number of adsorption iron sites. It is likely that the coordination of the nitrogen atoms and oxygen in RDB molecules to the Fe^{2+} ions in the ZIF8 framework is responsible for the more efficient adsorption compared with bare ZIF8. From this discussion the possible mechanisms of RDB adsorption onto ZIF8 or FeZIF8 were illustrated in Figure 11.
3.2.3. Isotherm Adsorption Studies
The study of adsorption isotherms is helpful in determining the adsorption capacities for removal at certain dyes at fixed temperature. In the present work, the equilibrium experiments were operated as if the initial concentration ( mg·L^{−1}) was kept constant and the absorbent weight varied between 0.005, 0.007, 0.01, 0.12, 0.015, 0.017, and 0.020 g. The experimental data were analyzed according to the nonlinear form of Langmuir and Freundlich model. Figure 12 presents the experimental points and nonlinear regression curves of Langmuir and Freundlich models. The parameters of models are listed in Table 9. However, the maximum adsorption capacity is not obtained in Freundlich equation. Halsey [68] supposed that the maximum adsorption capacity, , by Freundlich equation could be expressed: calculated based on Freundlich equation is also shown in Table 9.

The determination coefficient and the Chisquare test for assessing the compatibility of experimental data with the isothermal models are listed in Table 9. The high value of and low suggest that the isothermal data of ZIF8 could be well represented by the Langmuir model. This implies a monolayer adsorption for ZIF8. For FeZIF8 both models exhibited similar values of and . Moreover, favorable characteristic parameters of for Langmuir isotherm and for Freundlich isotherm were and , which indicated that both isotherms were favorable. These results confirmed that the equilibrium data of RDB adsorption onto FeZIF8 could be well fitted by the two adsorption isotherm models. This result indicated a monolayer adsorption and the existence of heterogeneous surface in FeZIF8. It is worth noting that maximum monolayer adsorption capacity, , for ZIF8 is much higher than that for ZIF8. It was concluded that the introduction of iron into ZIF8 enhanced the RDB adsorption onto FeZIF8 in terms of kinetics and isotherm adsorption.
in Langmuir model is the equilibrium constant and describes the relation between kinetics and thermodynamics. Then, in Langmuir model is thought to be equivalent to . However, the relation of and () (see Table 5) was not clear in our study. It is explained by that fact that in the solidliquid adsorption system in the Langmuir model no longer reflects the equilibrium constant as its original meaning. In our previous study [40], as isothermal adsorption experiments were operated with constant initial concentration and variable weights of adsorbent, we found that the obtained parameters of Langmuir and Freundlich models also increase as the initial concentration increases, which should be constant at certain temperature. It is thought that these parameters turn out to be empirical coefficients rather than fixed parameters at a certain temperature.
3.2.4. Reusability of ZIF8 and FeZIF8
To estimate the reusability of ZIF8 and FeZIF8 for the removal RDB, the used adsorbents were regenerated by sonication assisted washing with 10^{−3} M NaOH solution for 6 h and drying for 24 hours at 100°C and then exploited to adsorb RDB. The adsorption capacity of the regenerated ZIF8 and FeZIF8 are presented in Figure 13. Although the RDB adsorption decreased gradually with an increase in desorption cycles, the regenerated adsorbents still exhibited good performance for their adsorption. After the desorption for three cycles, the RDB adsorption capacity of by ZIF8 and FeZIF8 reached 90% and 95% of the adsorption of the initial material, respectively. The XRD patterns of adsorbents after the three cycles seem unchangeable (see Figure 14), indicating that the adsorbents were stable in this condition.
4. Conclusion
The iron doped zeolite imidazolate framework8 was synthesized by hydrothermal process. Ferrous ion as an iron source could be directly introduced into ZIF8 to form FeZIF8 (molar mole of Fe/(Fe + Zn) = 0.116). The ferrous ions could replace partially Zn(II) in ZIF8 structure or exist in amorphous species. The ZIF8 was stable in aqueous solutions with pH in ranging from 3 to 12 for 24 hours. Both ZIF8 and FeZIF8 were used to study the RDB adsorption. The results show that the introduction of iron into ZIF8 significantly enhanced the RDB adsorption capacity compared to bare ZIF8. The study based on Webber’s intraparticle diffusion shows that the adsorption process with ZIF8 and FeZIF8 tended to follow two stages in which the first stage was filmdiffusion and the second was an intraparticlediffusion process. The kinetic parameters based on Arrhenius and Eyring equation proved that the introduction of iron into ZIF8 provided a much larger adsorption capacity and faster adsorption kinetics than ZIF8. In addition to the electrostatic interaction mechanism and the hydrophobic and ππ interaction between the aromatic rings of the RBB dye and the aromatic imidazole rings of the adsorbent for ZIF8, the coordination of the nitrogen atoms and oxygen in carboxyl group in RDB molecules to the Fe^{2+} ions in the ZIF8 framework might explain why FeZIF8 has a higher RBD adsorption capacity than ZIF8. The experimental data for ZIF8 were well correlated by Langmuir model while those for FeZIF8 were well fitted to both Langmuir and Freundlich models. The maximum monolayer adsorption capacity for FeZIF8 (193.56 mg^{−1}·g^{−1}) was approximately 1.4 times higher than that for ZIF8 (133.76 mg^{−1}·g^{−1}).
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was funded by Ministry of Education and Training, Vietnam, under the Project B2016DHH20.
References
 F. P. Van Der Zee and S. Villaverde, “Combined anaerobicaerobic treatment of azo dyes—a short review of bioreactor studies,” Water Research, vol. 39, no. 8, pp. 1425–1440, 2005. View at: Publisher Site  Google Scholar
 P. A. Soloman, C. A. Basha, M. Velan, V. Ramamurthi, K. Koteeswaran, and N. Balasubramanian, “Electrochemical degradation of Remazol Black B Dye effluent,” Clean—Soil, Air, Water, vol. 37, no. 11, pp. 889–900, 2009. View at: Publisher Site  Google Scholar
 N. F. Cardoso, R. B. Pinto, E. C. Lima et al., “Removal of remazol black B textile dye from aqueous solution by adsorption,” Desalination, vol. 269, no. 1–3, pp. 92–103, 2011. View at: Publisher Site  Google Scholar
 V. R. Ranjusha, R. Pundir, K. Kumar, M. G. Dastidar, and T. R. Sreekrishnan, “Biosorption of Remazol Black B dye (Azo dye) by the growing Aspergillus flavus,” Journal of Environ Sci Health A: Toxic/Hazardous Substances and Environmental Engineering, vol. 45, no. 10, pp. 1256–1263, 2010. View at: Publisher Site  Google Scholar
 C. Janiak and J. K. Vieth, “MOFs, MILs and more: concepts, properties and applications for porous coordination networks (PCNs),” New Journal of Chemistry, vol. 34, no. 11, pp. 2366–2388, 2010. View at: Publisher Site  Google Scholar
 O. M. Yaghi, M. O'Keeffe, N. W. Ockwig, H. K. Chae, M. Eddaoudi, and J. Kim, “Reticular synthesis and the design of new materials,” Nature, vol. 423, no. 6941, pp. 705–714, 2003. View at: Publisher Site  Google Scholar
 P. Chowdhury, C. Bikkina, and S. Gumma, “Gas adsorption properties of the chromiumbased metal organic framework MIL101,” Journal of Physical Chemistry C, vol. 113, no. 16, pp. 6616–6621, 2009. View at: Publisher Site  Google Scholar
 L. Hamon, C. Serre, T. Devic et al., “Comparative study of hydrogen sulfide adsorption in the MIL53(Al, Cr, Fe), MIL47(V), MIL100(Cr), and MIL101(Cr) metalorganic frameworks at room temperature,” Journal of the American Chemical Society, vol. 131, no. 25, pp. 8775–8777, 2009. View at: Publisher Site  Google Scholar
 Y. Li and R. T. Yang, “Hydrogen storage in metalorganic and covalentorganic frameworks by spillover,” AIChE Journal, vol. 54, no. 1, pp. 269–279, 2008. View at: Publisher Site  Google Scholar
 P. L. Llewellyn, S. Bourrelly, C. Serre et al., “High uptakes of CO2 and CH4 in mesoporous metalorganic frameworks MIL100 and MIL101,” Langmuir, vol. 24, no. 14, pp. 7245–7250, 2008. View at: Publisher Site  Google Scholar
 J. Yang, Q. Zhao, J. Li, and J. Dong, “Synthesis of metalorganic framework MIL101 in TMAOHCr(NO_{3})_{3}H_{2}BDCH_{2}O and its hydrogenstorage behavior,” Microporous and Mesoporous Materials, vol. 130, no. 1–3, pp. 174–179, 2010. View at: Publisher Site  Google Scholar
 K. Yang, Q. Sun, F. Xue, and D. Lin, “Adsorption of volatile organic compounds by metalorganic frameworks MIL101: influence of molecular size and shape,” Journal of Hazardous Materials, vol. 195, pp. 124–131, 2011. View at: Publisher Site  Google Scholar
 Z. Zhang, S. Huang, S. Xian, H. Xi, and Z. Li, “Adsorption equilibrium and kinetics of CO_{2} on chromium terephthalate MIL101,” Energy and Fuels, vol. 25, no. 2, pp. 835–842, 2011. View at: Publisher Site  Google Scholar
 R. Kitaura, K. Seki, G. Akiyam, and S. Kitagawa, “Porous coordinationpolymer crystals with gated channels specific for supercritical gases,” Angewandte Chemie—International Edition, vol. 42, no. 4, pp. 428–431, 2003. View at: Publisher Site  Google Scholar
 S. Ma, D. Sun, X.S. Wang, and H.C. Zhou, “A meshadjustable molecular sieve for general use in gas separation,” Angewandte ChemieInternational Edition, vol. 46, no. 14, pp. 2458–2462, 2007. View at: Publisher Site  Google Scholar
 D.Y. Hong, Y. K. Hwang, C. Serre, G. Férey, and J.S. Chang, “Porous chromium terephthalate MIL101 with coordinatively unsaturated sites: surface functionalization, encapsulation, sorption and catalysis,” Advanced Functional Materials, vol. 19, no. 10, pp. 1537–1552, 2009. View at: Publisher Site  Google Scholar
 Y. K. Hwang, D.Y. Hong, J.S. Chang et al., “Selective sulfoxidation of aryl sulfides by coordinatively unsaturated metal centers in chromium carboxylate MIL101,” Applied Catalysis A: General, vol. 358, no. 2, pp. 249–253, 2009. View at: Publisher Site  Google Scholar
 N. V. Maksimchuk, M. N. Timofeev, M. S. Melgunov et al., “Heterogeneous selective oxidation catalysts based on coordination polymer MIL101 and transition metalsubstituted polyoxometalates,” Journal of Catalysis, vol. 257, no. 2, pp. 315–323, 2008. View at: Publisher Site  Google Scholar
 Z. Saedi, S. Tangestaninejad, M. Moghadam, V. Mirkhani, and I. MohammadpoorBaltork, “MIL101 metalorganic framework: a highly efficient heterogeneous catalyst for oxidative cleavage of alkenes with H_{2}O_{2},” Catalysis Communications, vol. 17, pp. 18–22, 2012. View at: Publisher Site  Google Scholar
 A. Phan, C. J. Doonan, F. J. UribeRomo, C. B. Knobler, M. O'Keeffe, and O. M. Yaghi, “Synthesis, structure, and carbon dioxide capture properties of zeolitic imidazolate frameworks,” Accounts of Chemical Research, vol. 43, no. 1, pp. 58–67, 2010. View at: Publisher Site  Google Scholar
 J. R. Long and O. M. Yaghi, “The pervasive chemistry of metalorganic frameworks,” Chemical Society Reviews, vol. 38, no. 5, pp. 1213–1214, 2009. View at: Publisher Site  Google Scholar
 H. Wu, W. Zhou, and T. Yildirim, “Hydrogen storage in a prototypical zeolitic imidazolate framework8,” Journal of the American Chemical Society, vol. 129, no. 17, pp. 5314–5315, 2007. View at: Publisher Site  Google Scholar
 H. Bux, A. Feldhoff, J. Cravillon, M. Wiebcke, Y.S. Li, and J. Caro, “Oriented zeolitic imidazolate framework8 membrane with sharp H_{2}/C_{3}H_{8} molecular sieve separation,” Chemistry of Materials, vol. 23, no. 8, pp. 2262–2269, 2011. View at: Publisher Site  Google Scholar
 J. Cravillon, S. Münzer, S.J. Lohmeier, A. Feldhoff, K. Huber, and M. Wiebcke, “Rapid roomtemperature synthesis and characterization of nanocrystals of a prototypical zeolitic imidazolate framework,” Chemistry of Materials, vol. 21, no. 8, pp. 1410–1412, 2009. View at: Publisher Site  Google Scholar
 J.Q. Jiang, C.X. Yang, and X.P. Yan, “Zeolitic imidazolate framework8 for fast adsorption and removal of benzotriazoles from aqueous solution,” ACS Applied Materials and Interfaces, vol. 5, no. 19, pp. 9837–9842, 2013. View at: Publisher Site  Google Scholar
 K.Y. A. Lin and H.A. Chang, “Efficient adsorptive removal of humic acid from water using zeolitic imidazole framework8 (ZIF8),” Water, Air, and Soil Pollution, vol. 226, article 10, 2015. View at: Publisher Site  Google Scholar
 B. K. Jung, J. W. Jun, Z. Hasan, and S. H. Jhung, “Adsorptive removal of parsanilic acid from water using mesoporous zeolitic imidazolate framework8,” Chemical Engineering Journal, vol. 267, pp. 9–15, 2015. View at: Publisher Site  Google Scholar
 M. Jian, B. Liu, G. Zhang, R. Liu, and X. Zhang, “Adsorptive removal of arsenic from aqueous solution by zeolitic imidazolate framework8 (ZIF8) nanoparticles,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 465, pp. 67–76, 2015. View at: Publisher Site  Google Scholar
 J. Zheng, C. Cheng, W.J. Fang et al., “Surfactantfree synthesis of a Fe_{3}O_{4}@ZIF8 coreshell heterostructure for adsorption of methylene blue,” CrystEngComm, vol. 16, no. 19, pp. 3960–3964, 2014. View at: Publisher Site  Google Scholar
 S. Eslava, L. Zhang, S. Esconjauregui et al., “Metalorganic framework ZIF8 films as lowκ dielectrics in microelectronics,” Chemistry of Materials, vol. 25, no. 1, pp. 27–33, 2013. View at: Publisher Site  Google Scholar
 M. Zhu, S. R. Venna, J. B. Jasinski, and M. A. Carreon, “Roomtemperature synthesis of ZIF8: the coexistence of ZnO nanoneedles,” Chemistry of Materials, vol. 23, no. 16, pp. 3590–3592, 2011. View at: Publisher Site  Google Scholar
 A. Kumar, B. Prasad, and I. M. Mishra, “Adsorptive removal of acrylonitrile by commercial grade activated carbon: kinetics, equilibrium and thermodynamics,” Journal of Hazardous Materials, vol. 152, no. 2, pp. 589–600, 2008. View at: Publisher Site  Google Scholar
 N. Kannan and M. Meenakshisundaram, “Adsorption of Congo Red on various activated carbons. A Comparative Study,” Water, Air, and Soil Pollution, vol. 138, no. 1–4, pp. 289–305, 2002. View at: Publisher Site  Google Scholar
 A. J. Ahamed, V. Balakrishman, and S. Arivoli, “Kinetic and equilibrium studies of Rhodamine B adsorption by low cost activated carbon,” Archives of Applied Science Research, vol. 3, pp. 154–166, 2011. View at: Google Scholar
 P. Atkins and J. D. Paula, Physical Chemistry, Oxford University Press, New York, NY, USA, 2010.
 J. Crank, The Mathematics of Diffusion, Clarendon Press, London, UK, 1975. View at: MathSciNet
 W. J. Weber and J. C. Morris, “Kinetics of adsorption on carbon from solution,” Journal of the Sanitary Engineering Division Proceed American society of civil Engineers, vol. 89, no. 2, pp. 31–60, 1963. View at: Google Scholar
 G. F. Malash and M. I. ElKhaiary, “Piecewise linear regression: a statistical method for the analysis of experimental adsorption data by the intraparticlediffusion models,” Chemical Engineering Journal, vol. 163, no. 3, pp. 256–263, 2010. View at: Publisher Site  Google Scholar
 H. Motulsky and A. Christopoulos, Fitting Models to Biological Data Using Linear and NonLinear Regression, GraphPad Software, San Diego, Calif, USA, 2003.
 B. H. Dang Son, V. Quang Mai, D. Xuan Du, N. Hai Phong, and D. Quang Khieu, “A study on astrazon black AFDL dye adsorption onto Vietnamese diatomite,” Journal of Chemistry, vol. 2016, Article ID 8685437, 11 pages, 2016. View at: Publisher Site  Google Scholar
 T. S. Anirudhan and P. G. Radhakrishnan, “Thermodynamics and kinetics of adsorption of Cu(II) from aqueous solutions onto a new cation exchanger derived from tamarind fruit shell,” The Journal of Chemical Thermodynamics, vol. 40, no. 4, pp. 702–709, 2008. View at: Publisher Site  Google Scholar
 K. G. Scheckel and D. L. Sparks, “Temperature effects on nickel sorption kinetics at the mineralwater interface,” Soil Science Society of America Journal, vol. 65, no. 3, pp. 719–728, 2001. View at: Publisher Site  Google Scholar
 E. I. Unuabonah, K. O. Adebowale, and B. I. OluOwolabi, “Kinetic and thermodynamic studies of the adsorption of lead (II) ions onto phosphatemodified kaolinite clay,” Journal of Hazardous Materials, vol. 144, no. 12, pp. 386–395, 2007. View at: Publisher Site  Google Scholar
 Y. Liu, “Is the free energy change of adsorption correctly calculated?” Journal of Chemical and Engineering Data, vol. 54, no. 7, pp. 1981–1985, 2009. View at: Publisher Site  Google Scholar
 Q. Song, S. K. Nataraj, M. V. Roussenova et al., “Zeolitic imidazolate framework (ZIF8) based polymer nanocomposite membranes for gas separation,” Energy and Environmental Science, vol. 5, no. 8, pp. 8359–8369, 2012. View at: Publisher Site  Google Scholar
 I. Langmuir, “The constitution and fundamental properties of solids and liquids. Part I. Solids,” The Journal of the American Chemical Society, vol. 38, no. 2, pp. 2221–2295, 1916. View at: Publisher Site  Google Scholar
 T. W. Weber and R. K. Chakravorti, “Pore and solid diffusion models for fixed‐bed adsorbers,” AIChE Journal, vol. 20, no. 2, pp. 228–238, 1974. View at: Publisher Site  Google Scholar
 H. M. F. Freundlich, “Over the adsorption in solution,” Journal of Physical Chemistry, vol. 57, pp. 385–471, 1906. View at: Google Scholar
 I. Tosun, “Ammonium removal from aqueous solutions by clinoptilolite: determination of isotherm and thermodynamic parameters and comparison of kinetics by the double exponential model and conventional kinetic models,” International Journal of Environmental Research and Public Health, vol. 9, no. 3, pp. 970–984, 2012. View at: Publisher Site  Google Scholar
 I. A. W. Tan, A. L. Ahmad, and B. H. Hameed, “Adsorption of basic dye on highsurfacearea activated carbon prepared from coconut husk: equilibrium, kinetic and thermodynamic studies,” Journal of Hazardous Materials, vol. 154, no. 1–3, pp. 337–346, 2008. View at: Publisher Site  Google Scholar
 A. Asfaram, M. Ghaedi, and G. R. Ghezelbash, “Biosorption of Zn^{2+}, Ni^{2+} and Co^{2+} from water samples onto Yarrowia lipolytica ISF7 using a response surface methodology, and analyzed by inductively coupled plasma optical emission spectrometry (ICPOES),” RSC Advances, vol. 6, no. 28, pp. 23599–23610, 2016. View at: Publisher Site  Google Scholar
 O. Karagiaridi, M. B. Lalonde, W. Bury, A. A. Sarjeant, O. K. Farha, and J. T. Hupp, “Opening ZIF8: a catalytically active zeolitic imidazolate framework of sodalite topology with unsubstituted linkers,” Journal of the American Chemical Society, vol. 134, no. 45, pp. 18790–18796, 2012. View at: Publisher Site  Google Scholar
 Y. Du, R. Z. Chen, J. F. Yao, and H. T. Wang, “Facile fabrication of porous ZnO by thermal treatment of zeolitic imidazolate framework8 and its photocatalytic activity,” Journal of Alloys and Compounds, vol. 551, no. 25, pp. 125–130, 2013. View at: Publisher Site  Google Scholar
 X. Zeng, L. Huang, C. Wang, J. Wang, J. Li, and X. Luo, “Sonocrystallization of ZIF8 on electrostatic spinning TiO_{2} nanofibers surface with enhanced photocatalysis property through synergistic effect,” ACS Applied Materials and Interfaces, vol. 8, no. 31, pp. 20274–20282, 2016. View at: Publisher Site  Google Scholar
 V. Vadivelan and K. Vasanth Kumar, “Equilibrium, kinetics, mechanism, and process design for the sorption of methylene blue onto rice husk,” Journal of Colloid and Interface Science, vol. 286, no. 1, pp. 90–100, 2005. View at: Publisher Site  Google Scholar
 G. A. Seber and C. J. Wild, Nonlinear Regression, Wiley Series in Probability and Mathematical Statistics: Probability and Mathematical Statistics, John Wiley & Sons, New York, NY, USA, 1989. View at: Publisher Site  MathSciNet
 M. A. AlGhouti, M. A. M. Khraisheh, M. N. M. Ahmad, and S. Allen, “Adsorption behaviour of methylene blue onto Jordanian diatomite: a kinetic study,” Journal of Hazardous Materials, vol. 165, no. 1–3, pp. 589–598, 2009. View at: Publisher Site  Google Scholar
 A. J. Ahamed and B. J. Suganthana, “Adsorption of CrVI on commercially available ash powder—a kinetic approach,” Journal of the Indian Chemical Society, vol. 83, pp. 792–795, 2006. View at: Google Scholar
 B. Y. M. Bueno, M. L. Torem, F. Molina, and L. M. S. de Mesquita, “Biosorption of lead(II), chromium(III) and copper(II) by R. opacus: equilibrium and kinetic studies,” Minerals Engineering, vol. 21, no. 1, pp. 65–75, 2008. View at: Publisher Site  Google Scholar
 C. E. Webster, R. S. Drago, and M. C. Zerner, “Molecular dimensions for adsorptives,” Journal of the American Chemical Society, vol. 120, no. 22, pp. 5509–5516, 1998. View at: Publisher Site  Google Scholar
 R. R. Sheha and A. A. ElZahhar, “Synthesis of some ferromagnetic composite resins and their metal removal characteristics in aqueous solutions,” Journal of Hazardous Materials, vol. 150, no. 3, pp. 795–803, 2008. View at: Publisher Site  Google Scholar
 N. K. Lazaridis, T. D. Karapantsios, and D. Georgantas, “Kinetic analysis for the removal of a reactive dye from aqueous solution onto hydrotalcite by adsorption,” Water Research, vol. 37, no. 12, pp. 3023–3033, 2003. View at: Publisher Site  Google Scholar
 Y. S. Ho and G. McKay, “Pseudosecond order model for sorption processes,” Process Biochemistry, vol. 34, no. 5, pp. 451–465, 1999. View at: Publisher Site  Google Scholar
 J. X. Lin, S. L. Zhan, M. H. Fang, and X. Q. Qian, “The adsorption of dyes from aqueous solution using diatomite,” Journal of Porous Materials, vol. 14, no. 4, pp. 449–455, 2007. View at: Publisher Site  Google Scholar
 K. G. Scheckel and D. L. Sparks, “Dissolution kinetics of nickel surface precipitates on clay mineral and oxide surfaces,” Soil Science Society of America Journal, vol. 66, pp. 689–694, 2001. View at: Google Scholar
 C.C. Kan, M. C. Aganon, C. M. Futalan, and M. L. P. Dalida, “Adsorption of Mn^{2+} from aqueous solution using fe and mn oxidecoated sand,” Journal of Environmental Sciences, vol. 25, no. 7, pp. 1483–1491, 2013. View at: Publisher Site  Google Scholar
 N. A. Khan, B. K. Jung, Z. Hasan, and S. H. Jhung, “Adsorption and removal of phthalic acid and diethyl phthalate from water with zeolitic imidazolate and metalorganic frameworks,” Journal of Hazardous Materials, vol. 282, pp. 194–200, 2015. View at: Publisher Site  Google Scholar
 G. D. Halsey, “The role of surface heterogeneity in adsorption,” Advances in Catalysis, vol. 4, pp. 259–269, 1952. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Mai Thi Thanh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.