1 Introduction
The pollution of water resources with heavy metals has been causing worldwide concern in the last few decades. Wastewater from many industries contains one or more toxic heavy metals [1]. Mercury is one of the most toxic heavy metals commonly found in the global environment and its toxic nature has been known for centuries. It is primarily used in battery, gauges, and monometer manufacturing. It has a very high tendency to bind to proteins and due to that it mainly affects the renal and nervous systems [2]. Consequently, the selective separation of Hg2+ from natural samples needs much more attention.
In recent years, many studies have been focused on the adsorption of mercury [3–6]. Walcarius and Delacote [7] and Aguado et al. [8] used ordered mesoporous silica materials for the adsorption of mercury and reported that the adsorbent had higher adsorption capacity than other silica materials. The ordered mesoporous silica materials, due to their regular pore array with uniform and tunable mesopore diameter, high surface area and pore volume, have proved to be of considerable interest compared to other silica materials [9].
On the other side, zeolites have received much attention because of their regular porous structure with excellent stability [10], but they have poor adsorption capacity because of diffusion limitation in the microporous network [11], while mesoporous silica possesses poor stability compared to zeolites [12]. The hierarchical zeolites may be considered as desirable candidates for adsorption because they show a dramatically improved performance over the mesoporous materials. The reason is that hierarchical zeolites combine the strongly hydrothermal stability characteristic of zeolites with the excellent mass transport properties of mesoporous materials [13]. Extensive researches have been performed regarding the preparation and application of these types of zeolites in various reactions [14–18]. Nevertheless, limited researches have been accomplished on the adsorption of heavy metals using hierarchical zeolites.
In this study, thiol-functionalized beta zeolite-MCM-41 nanocomposite (SH-Beta/MCM-41) was prepared by hydrothermal procedure and post-synthesis grafting of 3-mercaptopropyltrimethoxysilane (MPTMS). The synthesized adsorbent was employed for the extraction of Hg2+ from aqueous solutions. The effect of different parameters, including pH, initial concentration, contact time and temperature on Hg2+ adsorption was investigated and optimized. Adsorption isotherms, thermodynamic parameters and kinetics of the process were studied by application of different models. The selectivity and reusability of the adsorbent was also examined.
2 Experimental
2.1 Materials and methods
Tetraethylammoniumhydroxide (20%, Merck), cetyltrimethylammoniumbromide (CTAB, Merck), aerosol-200 (Fluka) as silica source, 3-mercaptopropyltrimethoxysilane (MPTMS, ≥ 95%, Merck), toluene (≥ 99%, Merck), ethanol (≥ 99.9%, Merck) and mercury nitrate [Hg(NO3)2, ≥ 98.0%, Merck] were used without further purification. All aqueous solutions were prepared with double-distilled water ≤ 1.8 μS.cm−1. FTIR spectra of the samples were obtained using a PerkinElmer spectrum 65 FTIR spectrophotometer (wavenumber: 400–4000 cm−1, USA) with the KBr pressed-disk technique. Thermogravimetric analysis was performed using a Mettler TG-50 thermal analyzer from ambient temperature to 800 °C. The atmosphere was pure nitrogen and the heating rate was 15 °C·min−1. X-ray diffraction patterns were recorded using a Bruker (Germany), D8 ADVANCE X-ray diffractometer using Cu Kα radiation (wavelength: 1.5406 Ǻ and filter: Ni) up to 2θ = 45° at ambient temperature. Nitrogen adsorption–desorption isotherms were determined at 77 K and the specific surface area was measured by applying the BET equation to the isotherm. The pore size distribution was calculated using the adsorption branch of the isotherm and the Barrett–Joyner–Halenda (BJH) formula. For measurement on mercury, a PerkinElmer A Analyst 300 atomic absorption spectrophotometer was used.
2.2 Preparation of materials
Beta/MCM-41 composite was synthesized by the procedure suggested by Ooi et al. [19]. To prepare SH-Beta/MCM-41, 5.0 g of Beta/MCM-41 composite was refluxed in 250 mL of toluene containing 5.0 g of MPTMS at 110 °C for 24 h. The solid phase was then recovered by filtration and washed with toluene followed by ethanol. The residual organosilane was removed by the Soxhlet extraction method over ethanol for 24 h [20].
2.3 Characterization methods
FTIR spectra of the samples were obtained using a PerkinElmer (USA) spectrum 65 FTIR spectrophotometer (wavenumber: 400–4000 cm−1) by using the KBr pressed-disk technique. Thermogravimetric analysis was performed using a Mettler TG-50 thermal analyzer from ambient temperature to 800 °C. The atmosphere was pure nitrogen and the heating rate was 15 °C·min−1. X-ray diffraction patterns were taken using a Bruker (Germany), D8ADVANCE X-ray diffractometer using Cu Kα radiation (wavelength: 1.5406 Å and Ni filter) up to 2θ = 45° at ambient temperature. Nitrogen adsorption–desorption isotherms were determined at 77 K and the specific surface area was measured by applying the BET equation to the isotherm. The pore size distribution was calculated using the adsorption branch of the isotherm and the Barrett–Joyner–Halenda (BJH) formula. A PerkinElmer AAnalyst 300 cold-vapor atomic absorption spectrophotometer was used to measure mercury concentration.
2.4 Adsorption experiments
The batch technique was used to study the kinetics and equilibrium of the adsorption process. In order to evaluate the kinetic data, a known amount of the adsorbent was added in polyethylene bottles containing 25 mL of Hg2+ solution. A separate bottle was prepared for each time interval. The adsorption isotherms were constructed at four different temperatures of 25, 35, 45 and 55 °C by putting 0.050 g of the adsorbent in contact with 25 mL of the mercury solution (50 to 1500 mg·L−1). The effect of the pH was studied in the range from 1 to 4 (adjusted with HNO3). After the required time, the adsorbent was separated by filtering and the concentration of mercury in the filtrate was measured by cold-vapor atomic absorption spectrophotometry. The equilibrium uptake was calculated by the following equation [9]:
(1) |
3 Results and discussion
3.1 Adsorbent characterization
The FTIR spectra of Beta/MCM-41 and SH-Beta/MCM-41 obtained by the KBr pellet method are shown on Fig. 1. The broad band that appeared around 3422 and 3452 cm−1 was mainly attributed to the Si–OH stretching vibration of water adsorbed on the surface of the silica material. The peak at about 1630 cm−1 was also attributed to adsorbed water [21]. In both samples, the Si–OH–Al band appeared at ∼ 3630 cm−1, indicating that the zeolite structure was not changed during functionalization with mercaptosilane groups [22]. The wide bands around 1080 and the smaller band at 800 cm−1 were attributed to Si–O–Si asymmetric and symmetric stretching vibrations, respectively [23]. The weak band around 2550 cm−1 in Fig. 1b, which is typically very weak due to the aggregation of mercapto groups within the monolayer and to hydrogen binding effects [24,25], and the peak at 2933 cm−1 belonged to the stretching vibrations of CH bonds, indicating the presence of organic groups adsorbed on the surface of the adsorbent [25].
The TG and DTG curves of Beta/MCM-41 and SH-Beta/MCM-41 are shown in Fig. 2. In the thermal curves of Beta/MCM-41, a characteristic dehydration peak between 25 and 200 °C was observed (Fig. 2a) [26], while in the thermal curve of SH-Beta/MCM-41, two distinct weight loss events were observed. The first weight loss that appeared around 100 °C was attributed to the dehydration of the water molecules adsorbed on the surface of the sample, whereas the second weight loss around 330 °C was assigned to the decomposition of mercaptosilane groups (Fig. 2b) [8,27].
In the XRD pattern of the calcined Beta/MCM-41 composite, the characteristic diffraction lines of beta zeolite at 7.9° and 22.5° were clearly observed (Fig. 3a). For the functionalized sample, SH-Beta/MCM-41, similar diffraction lines were observed, with a slightly lower intensity (Fig. 3b) [28]. The diffraction lines of MCM-41 were not observed in low-angle diffraction patterns (Fig. 3c). This was in agreement with the results reported by Ooi et al. [19].
In the SEM images of Beta/MCM-41 and SH-Beta/MCM-41, massive and spherical beads were observed. The beads of Beta/MCM-41 composite had a relatively smooth surface, indicating that the composite was homogeneously formed (Fig. 4a) [19,29]. The SH-Beta/MCM-41 image revealed that the particles are more agglomerated. This was attributed to the presence of mercaptosilane groups on the surface of the composite (Fig. 4b) [9].
Fig. 5 shows the adsorption–desorption isotherms and the pore size distribution. The isotherm of Beta/MCM-41 has a typical hysteresis loop in the range P/P0 = 0.75–0.95, indicating the presence of framework-confined mesopores [30]. In addition, in the isotherm of SH-Beta/MCM-41, the gap of adsorption and desorption is low because the mercaptosilane groups are located within the mesopores and the occupied pores volume [31]. The BET surface area and the average pore size of the samples are shown in Table 1. The Beta/MCM-41 composite sample synthesized via two-step crystallization had a less pronounced mesopore size than mesoporous materials [19]. By grafting mercaptosilane groups on the sample, the BET surface area of the new composite (SH-Beta/MCM-41) decreased from 604 to 82 m2·g−1. This was attributed to the occupancy of mercaptosilane groups into the mesoporous structure of the host sample by occupation of the pores [27]. Also, due to the presence of the propyl chains of mercaptosilane groups in the pores of Beta/MCM-41, the pore volumes were reduced from 0.617 to 0.295 cm3·g−1 and 0.42 to 0.005 cm3·g−1 for the mesopores and the micropores, respectively. The presence of this group also significantly expanded the pore diameters of meso- and micropores [32] (Table 1).
N2 adsorption/desorption data.
Samples | SBET (m2·g−1) | Vp(meso) (cm3·g−1) | Vp(micro) (cm3·g−1)a | dp(meso) (nm) | dp(micro) (nm)a |
Beta/MCM-41 | 604 | 0.617 | 0.42 | 4.08 | 1.20 |
SH-Beta/MCM-41 | 82 | 0.295 | 0.005 | 14.25 | 1.90 |
a Calculated by MP-Plot method.
The micropore sizes of Beta/MCM-41 and SH-Beta/MCM-41 samples were calculated using MP-Plot in Table 1. Further MP-Plot analyses also confirm the existence of microporosity within the samples (Fig. 6).
3.2 Adsorption optimizing
The effect of initial concentration on the adsorption capacity was studied in the 50–1500 mg·L−1concentration range. As shown on Fig. 7, with increasing cation concentrations, the adsorption rate increased and then leveled off. The optimized uptake was obtained at a concentration of 1000 mg·L−1. At pH = 2, an extremely high adsorption capacity of 6.6 mequiv·g−1 was obtained, which was much higher than the one previously reported. Kao et al. [4] reported that the adsorption capacities of their adsorbent for Hg2+ binding increased in proportion to the amount of thiol groups incorporated and the maximum adsorption capacity reached 415 mg·g−1. The maximum adsorption capacity of a thiol-SNHS nano-adsorbent prepared by Rafati et al. in 330 mg·L−1 of Hg2+ solution was 160 mg·g−1 [9]. The high adsorption capacity was attributed to the presence of thiol groups on the surface of the adsorbent. The thiol group, as soft bases, has highly polarizable donor centers capable to strongly interact with low-lying orbitals of soft acids [33].
Adsorption isotherm models were used to describe experimental adsorption data. The model parameters and the underlying thermodynamic assumptions of these adsorption models can provide some insight into the sorption mechanism, the surface properties, and the affinity of the adsorbent. Therefore, obtaining the ‘best-fit’ isotherm is very important. Langmuir and Freundlich isotherms were examined. In the Langmuir isotherm model, the assumption is that the adsorption process is mono-layered [34]:
(2) |
(3) |
Isotherms parameters for Hg2+ adsorption on SH-Beta/MCM-41.
Adsorbent | Langmuir isotherm | Freundlich isotherm | ||||
R2 | KL | qm | R2 | KF | n | |
SH-Beta/MCM-41 | 0.8882 | 1.316 × 10−3 | 111.11 | 0.9692 | 6.996 | 1.539 |
In the Freundlich isotherm model, the assumption is that the adsorption process is multi-layered [34]:
(4) |
Measurement of the adsorption capacity at different pH values indicated that Hg2+ binding to SH-Beta/MCM-41 strongly depended on the pH of the solution (Fig. 10). At higher pH, the formation of metal hydroxides coincides with the absorption process. At a concentration of 1000 mg·L−1 of Hg2+, the pH of precipitation was 2.4. Therefore, the adsorption capacity obtained beyond this pH was not valid [36]. At lower pH, because of the competition for adsorption between H3O+ and Hg2+ ions on the surface of the adsorbent, the adsorption capacity toward Hg2+ was lowered. Moreover, in aqueous solutions, divalent metal ions may exist in different forms—M(OH)1+, M(OH)20, and M(OH)31−—and their abundance is dependent on the pH of the solution. The most favorite species for adsorption on the adsorbent surface is the M2+ species, which would be the dominant species at a particular pH for each individual cation [7,36]. Walcarius and Delacote investigated the effect of pH on the adsorption of mercury by thiol-functionalized mesoporous silica. By plotting the species distribution curves of Hg2+, Hg(OH)+ and Hg(OH)2 against the pH of the solution, they concluded that the maximal adsorption capacity was observed at pH = 2 because Hg2+ was the most abundant species at this pH [7].
The effect of contact time on Hg2+ adsorption was studied in 1000 mg·L−1 mercury solutions with an initial pH of 2.0 (Fig. 11). Adsorption was kinetically fast and the equilibration was attained after 5 min. At the beginning of the experiment, Hg2+ concentration was high and therefore the number of effective interactions between the adsorbates and the adsorbents was fast. Then the adsorbent surface was saturated by the cations and the rate of adsorption remained constant [37].
In order to analyze the adsorption kinetic parameters, two kinetic models, namely, pseudo-first-order and pseudo-second-order models, were applied.
The pseudo-first-order kinetic model can be expressed in linear form [9]:
(5) |
Kinetic parameters for Hg2+ adsorption on SH-Beta/MCM-41.
Kinetic models | K | R2 |
Pseudo-first-order | 0.5861 | 0.9849 |
Pseudo-second-order | 8.141 × 10−4 | 0.9989 |
The pseudo-second-order kinetic model can be expressed in linear form [9]:
(6) |
- • migration of the metal ions from the solution to the surface of the adsorbent;
- • migration of the metal ions into the pores of the adsorbent;
- • adsorption of the metal ions on the interior surface of the adsorbent [9].
The rate of the reaction is mainly controlled by the first two steps.
Removal of mercury ions from the aqueous solution occurs mainly either by electrostatic interaction (ionic interaction between positively charged metal ions and negatively charged matrices) or by chelation (donation of the lone-pair electrons of the matrices to metal ions to form coordinate bonds) (Scheme 1) [38].
The adsorption capacity of the adsorbent was measured at a constant Hg2+ concentration solution and at four different temperatures (Fig. 14). Adsorption capacity decreased from 8.2 mequiv·g−1 at 25 °C to 6.1 mequiv·g−1 at 55 °C. This is in agreement with the results reported by Arencibia et al. [39]: the mercury loading changes from 1.22 mmol·Hg·g−1 at 20 °C to 1.16 mmol·Hg·g−1 at 60 °C, which was attributed to the weakening of the bond between the cation and the sorption sites. Absorption capacity decreased with increasing temperature due to the weakening of the bond between the cation and the sorption sites. This was in agreement with the results reported by Arencibia et al. [39]. The change in Gibbs free energy, enthalpy and entropy of adsorption were calculated using the following equations [40]:
(7) |
(8) |
According to Van’t Hoff equation:
(9) |
The plot of lnKd versus 1/T was linear. The negative values of ΔH0 evidenced the exothermic nature of adsorption (Fig. 15). The negative values of ΔG0 indicated that the adsorption process led to a decrease in Gibbs free energy and that the adsorption process was feasible and spontaneous for Hg2+ at each temperature. The positive values of ΔS0 represented the increase of irregularity during the absorption process; this issue was predicted to increase the mobility of the molecules due to the temperature rise. The thermodynamic parameters for the removal of Hg2+ were given in Table 4.
Thermodynamic parameters for Hg2+ adsorption on SH-Beta/MCM-41.
ΔH0 (J·mol−1) | ΔS0·(J·mol−1·K−1) | ΔG0 (J·mol−1) | |||
−8028.91 | 28.87 | 298 K | 308 K | 318 K | 328 K |
−16632.17 | −16920.87 | −17209.57 | −17498.27 |
3.3 Reusability of the adsorbent
Reusability feature of the functionalized composite was examined for five regeneration cycles by washing the used adsorbent with a thiourea solution (5% w/v) and HNO3 3 M [6]. The optimized regeneration value was obtained with the 1:1 thiourea–HNO3 solution. The FTIR spectra taken from the sample treated with 1:1 thiourea–HNO3 showed no significant changes (Fig. 1c).
Arencibia et al. [39] compared a variety of complexing agents, such as KBr, KSCN and thiourea, and reported that thiourea had the highest complexing ability with mercury and was able to remove the higher amount of mercury from the aqueous solution. This effect can be attributed to the presence of sulfur and to the multidentate character of thiourea as a ligand, which leads to a larger mercury complexing capacity. The multidental character of thiourea and the presence of sulfur created larger mercury complexing capacity. It was a stronger nucleophile than the sulfur atom of the thiol group; therefore mercury was readily washed out into the solution [39] (Scheme 2).
The adsorption capacity of the adsorbent declined with the number of reuse cycles (Fig. 16). This was attributed to the blockage of adsorptive centers or the decomposition of some active sites during the regeneration process.
According to the hard–soft acid–base theory, Hg2+ is classified as a soft ion. Soft ions form very strong bonds with groups containing nitrogen and sulfur atoms [41]. Soft cations are tending to establish a covalent interaction and have an affinity toward the soft sulfur. Therefore, the sulfur atom in the mercaptosilane group acts as a soft base that has a favorable interaction with Hg2+ as a soft acid, and forms a hard bond (covalent bond) with Hg2+.
3.4 Selectivity of the adsorbent for mercury
The tolerance limit is defined as the interfering ion concentration causing a relative error ≥ 5% on the adsorption of the measured cation by the adsorbent [6]. To measure the tolerance limit of the interfering cations, Pb2+ and Ni2+ solutions containing 1000 mg/L of Hg2+ (optimized concentration) and of interfering cations (50–3000 mg·L−1) were prepared. Twenty-five milliliters of each solution were equilibrated with 0.010 g of the adsorbent. After equilibration, the mixture was filtered out and the amount of Hg2+ was determined. The results demonstrated that the adsorbent retained its selectivity for Hg up to a Pb2+ concentration of 1500 mg/L and a Ni2+ concentration of 2500 mg/L (Table 5).
Tolerance limits of competitive cations.
Competitive cations | Tolerance limit for adsorbent (mg·L−1) |
Pb(II) | 1500 |
Ni(II) | 2500 |
The selectivity of the adsorbent was also examined by measurement of competitive adsorption distribution coefficients Kd (mL/g), the selectivity coefficient (k) of Hg2+/M (M: Pb2+, Ni2+) by the following equations:
(10) |
(11) |
4 Conclusion
The functionalized hierarchical zeolite nanocomposite prepared in the frame of this research possessed the structural advantages of a zeolite and of a mesoporous material. It showed excellent capacity for the adsorption of mercury from aqueous solutions. The adsorption process was kinetically fast and the equilibrium state was established within 5 min. The adsorption data have been fitted to the Freundlich isotherm. The selectivity of the binding process, in particular with respect to the accessibility to the active centers in the presence of other metal ions interferences was studied. SH-Beta/MCM-41 was an excellent adsorbent to bind the mercury with high selectivity because the –SH group is highly sensitive to Hg2+ ion. Furthermore, a maximum adsorption capacity of 8.2 mequiv·g−1 was obtained, which was much higher than the earlier-reported values.
Acknowledgment
The support by the Islamic Azad University, Shahreza Branch (IAUSH) Research Council and Center of Excellence in Chemistry is gratefully acknowledged.