Competitive Adsorption of Cadmium (II) from Aqueous Solutions onto Nanoparticles of Water Treatment Residual
There is increasing interest in using water treatment residuals (WTRs) for heavy metals removal from wastewater due to their low cost, availability, and high efficiency in removing various pollutants. In this study, novel water treatment residuals nanoparticles (nWTRs) were prepared using high energy ball milling and used for efficient removal of Cd(II) in single- and multi-ion systems. The WTR nanoparticles demonstrated high removal efficiency for Cd from aqueous solution as the adsorption capacities of nWTR were 17 and 10 times higher than those of bulk WTR in single- and multielement systems, respectively. Noticeably, Cd(II) adsorption was clearly suppressed in the multi-ion system as Cu and Pb form the most stable monohydroxo complexes. Fourier transmission infrared (FTIR) analyses suggested the participation of OH−, O-Al-O, FeOH, and FeOOH entities in the adsorption process. The stability of Cd-nWTR surface complexes is evident as less than 0. 2% of adsorbed Cd(ll) was released at the highest Cd(II) concentration load after 4 consecutive desorption cycles. Moreover, the real efficiency of nWTR for Cd(II) removal from wastewater samples studied was calculated to be 98.35%. These results highlight the potential of nWTR for heavy metals removal from wastewater.
Aquatic ecosystems contamination especially with heavy metal ions is an environmental problem worldwide. Of all contaminant metals, cadmium is considered an extremely toxic heavy metal with acceptable levels one-tenth those of most of the other toxic metals [1, 2]. Cadmium pollution arises mainly from domestic sewage and the effluents of industries involved in electroplating, smelting, alloy manufacturing, pigments, plastic, cadmium-nickel batteries, fertilizers, pesticides, pigments and dyes, textile operations, and refining [3, 4]. The main form of Cd in contaminated water is Cd(II) and the remediation technologies available to reduce Cd concentrations in contaminated water systems include ion-exchange, solvent extraction, chemical precipitation, phytoextraction, ultrafiltration, reverse osmosis, electrodialysis, and adsorption [3–6].
However, most of these technologies have shown limitations in removing the toxic contaminants from contaminated water to safe levels and they are costly, laborious, and time-consuming [7, 8]. Development of cheap and green remediation technology is urgently needed. Water treatment residuals (WTRs), waste material resulting from the treatment of surface water with Al and Fe salts, have been gaining increased attention recently as an eco-friendly and low cost adsorbent. Because of their amorphous nature, WTRs have shown strong affinity for Ni, Cu, Pb, and Hg [9–12].
Moreover, conclusive evidences have demonstrated that the smaller size fraction of WTRs greatly increases their active surface area and adsorption performance . Recently, Elkhatib et al.  developed a method to produce nanoscale sorbent from water treatment residues using ball milling. They demonstrated that the removal efficiency of the produced nanoparticles for As was 16 times higher than the removal efficiency of bulk WTR .
The greater adsorption capacity of water treatment residuals nanoparticles (nWTRs) in comparison with larger particles suggests their use as a very promising candidate for contaminated water and wastewater remediation. The potential use of bulk WTRs, inexpensive and nonhazardous waste material, was already tested for the removal of many metals in literature .
However, information is not available in the literature on the use of nanoscale WTR as super adsorbent for cadmium (II) removal. Here we report a comparative study of bulk WTRs and WTR nanoparticles for treatment of Cd(II) contaminated water. The objectives of this study were therefore (i) to determine cadmium adsorption capacity of water treatment residual nanoparticles and to investigate the pertinent adsorption mechanism and (ii) to study the effects of operational parameters including adsorbent dosages, solution pH, and competing cations on the removal of Cd by nWTRs. Such a study is important to throw the lights on nanoscale WTRs as a remediation agent for wastewater and freshwater.
2.1. Synthesis and Characterization of Water Treatment Residuals Nanparticles (nWTRs)
The bulk WTR (mWTR) was obtained from the drinking water treatment plant in Kafr El-Dawar, El-Behera, Egypt. The mWTR samples were collected from the treatment plant disposal site (open dry land in the close vicinity of the water treatment plant) using a shovel. The samples were transferred into high density polyethylene containers and transported to the laboratory at Alexandria University. The selected physical and chemical characteristics of the studied drinking water treatment residuals (WTRs) are shown in Table 1.
In the laboratory, mWTRs samples were first air-dried at room temperature for a period of 4 weeks. Relatively homogeneous material was obtained by passing dried mWTR through two different sieves having 2 mm and 51 um of pore diameters.
Synthesis of water treatment residuals nanoparticles (nWTRs) was achieved by milling subsamples of uWTR (<51 um) using Fritsch planetary mono mill . The characteristics and element contents of nWTR, mWTR, and uWTR were investigated using scanning electron microscopy (SEM) with energy dispersive X-ray (INCAx-Sight model 6587, Oxford Instruments, UK) and Autosorb-iQ surface area analyzer (Quantachrome, USA) and the results have been reported elsewhere .
2.2. Sorption Isotherms and Maximum Sorption Capacity of nWTRs
Batch sorption experiments were carried out in 50 mL polypropylene centrifuge tubes using three sizes of WTR (2 mm, less than 51 microns, and less than 100 nm). The Cd concentrations used in all batch sorption experiments ranged from 0 to 160 mg Cd/L in a background electrolyte of 0.01 M K2SO4. The WTR-Cd mixtures were equilibrated on a slowly end-over-end shaker (4 rpm) for 24 h and centrifuged at 4000 rpm for 20 min and the obtained supernatant solutions were filtered through a 0.45 μm Millipore filter. All experiments were performed in triplicate, and solutions were analyzed for Cd by inductively coupled plasma spectrometry (ICPS). The sorption data were then fitted to different isotherm models and the best fit was used to calculate the maximum sorption capacity. The solid nWTR material used in these sorption experiments was recovered and examined for sorbed Cd via scanning electron microscopy equipped with an X-ray energy dispersive spectrometer (SEM-EDX).
2.3. Competitive Adsorption Behavior
In reality, the presence of multiple competing ions is more frequent than the existence of only one kind of ion. Therefore, to assess competition for adsorption sites, multi-ion solutions were prepared and tested. Solutions containing equal amounts of Pb, Cd, and Cu (multi-ion system) ranging from 5 to 160 mg/L were prepared in 0.01 M K2SO4 background electrolyte and mixed with nWTRs. The suspensions were equilibrated for 24 h, after which aliquots were collected, centrifuged, and filtered through 0.45 um.
2.4. Sorption Kinetics and Effect of pH
A sufficient mass of nanoscale water treatment residuals was used in 0.01 M K2SO4 background electrolyte. Adsorption experiments as a function of time and pH (adsorption edge) were conducted in single-ion and multi-ion systems. The equilibrium pH of the suspension used varied from 5.0 to 9.0 with a 2.0 pH unit increment using combined pH electrode with standard pH buffers. The suspensions were equilibrated for Cd adsorption for 24 h and an aliquot was collected at 5 min, 30 min, 1 h, 4 h, 8 h, and 24 h contact time. The aliquots were centrifuged and filtered through 0.45 um Millipore filter and the supernatant was analyzed for Cd using ICP. The amount of Cd adsorbed was calculated as the difference between the initial and final metal concentrations and the obtained data were fitted into different kinetic models.
2.5. Effect of Initial Metal Concentration and nWTR Dose
The effects of initial metal concentration and nWTR dose on the extent of adsorption were investigated at the initial metal concentrations of 10.0, 50.0, 200.0, and 500.0 mg/L. The experiments were performed by adding 0.02 g, 0.05 g, or 0.10 g of nWTR sample into 10.0 mL potions of Cd solutions. The mixtures were equilibrated for 24.0 hours and then removed for centrifugation. After centrifugation, the supernatant solutions were retransferred into clean falcon tubes, acidified, and kept in the refrigerator until analysis.
2.6. Desorption Studies
In order to assess the stability of adsorption, 10.0 mL portions of high purity water were added onto Cd-loaded nWTR (after sorption experiment) and were shaken for 24 hours. After separation of liquid and solid phases by centrifugation, the solutions were acidified and kept in the refrigerator until Cd analysis.
3. Results and Discussion
3.1. Characterization of Bulk WTR and nWTR
The microstructure of WTR nanoparticles as examined by SEM indicated that WTR nanoparticles are spherical and the representative single particle sizes are typically less than 100 nm in diameter (Figure 1(a)). The major elements in nWTR as revealed by SEM-EDX are Si, Al, and Fe, and they represent around 83% of the total elements (Figure 1(a)). The X-ray diffraction analysis  suggested that amorphous iron, aluminum (hydr) oxides, and silicon oxide dominated all nWTR, with no apparent crystalline iron–Al (hydr) oxides. After the reaction with Cd ions, the SEM image of nWTR (Figure 1(b)) shows formation of a coating layer on the surface of the nWTR which indicates that the reaction occurred on the surface of nWTR. In addition, the SEM-EDX analysis spectrum (Figure 1(b)) ascertained the appearance of a cadmium peak (2.20%) amongst the elements detected in Cd-saturated nWTR. The BET specific surface area (SSA) of nWTR sample revealed that SSA of nanoscale WTR (129.0 m2 g−1) is 2-3 times higher than SSA of mWTR (53.1 m2 g−1) which demonstrates the high reactivity of nWTR and makes it ideal candidate for water treatment and desalination.
3.2. Effect of Adsorbent Dose
The effect of nWTR dose on Cd removal from solutions at different Cd concentrations is shown in Figure 2. The Cd2+ adsorption efficiency increases with increasing adsorbent dose of nWTR from 0.02 to 0.10 g. The amount of Cd2+ adsorbed increased from 900 mg kg−1 at 0.02 g to 46,000 mg kg−1 nWTR at 0.10 g. The H-type adsorption isotherm of nWTR (0.10 g/L) dose suggests a very strong metal-nWTR surface interaction . The enhanced adsorption with increasing adsorbent dose is attributed to higher surface area and the more adsorption sites available at higher adsorbent dose. Therefore, the nWTR dose of 0.10 g was proposed as optimum nWTR dose required for efficient Cd removal at different initial Cd concentrations. The nWTR dose of 0.10 g was used in subsequent experiments.
3.3. Equilibrium Studies
The data of Cd adsorption onto three different particle sizes of WTR in single- and multielement systems were analyzed using seven isotherm models: Langmuir, Freundlich, Temkin, Fowler–Guggenheim, Kiselev, Elovich, and Hill–de Boer.
Equilibrium and Kinetic Models
Equilibrium Models Freundlich Langmuir Elovich isotherm model TemkinFowler–Guggenheim (FG) Kiselev Hill–de Boer
where is relative adsorption capacity of the adsorbent (mg1−1/n)L1/n g−1, is constant and measures the strength of adsorption, is amount of Cd adsorbed per gram of adsorbent, is equilibrium concentration in the bulk solution, is maximum adsorption capacity of the adsorbent, is Langmuir constant related to the free energy of adsorption, is Elovich equilibrium constant, θ is fractional coverage, is universal gas constant (kJ mol−1 K−1), is (−H) the variation of adsorption energy (kJ mol−1), is Temkin equilibrium constant (L mg−1), is Fowler–Guggenheim equilibrium constant (L mg−1), T is temperature (K), is interaction energy between adsorbed molecules (kJ mol−1), is Kiselev equilibrium constant (L mg−1), is constant of complex formation between adsorbed molecules, is Hill–de Boer constant (L mg−1), is energetic interaction between adsorbed molecules (kJ mol−1).
Kinetics Models Elovich Intraparticle diffusion Power function First-order model
where is amount of Cd adsorbed per gram of adsorbent in time , is initial adsorption rate (mg min−1), β is a constant related to surface coverage and activation energy, is apparent sorption diffusion rate coefficient (μg g−1 min−1), t is reaction time (min), is a constant (mg g−1), v is estimated apparent sorption rate (min−1), is amount of Cd adsorption at equilibrium, is apparent adsorption rate coefficient.
The experimental data were fitted to the seven isotherm models tested using linear forms . Model parameters as well as coefficient determination () and the standard error of estimate (SE) values of all models are presented in Equilibrium and Kinetic Models. of all models tested were quite high especially for nWTR. Comparing goodness-of-fit measures between the different 7 models tested shows that Langmuir model yielded the lowest SE values for all the sorbents followed by Temkin and Fowler–Guggenheim models (Equilibrium and Kinetic Models). The high SE values of Freundlich, Hill–de Boer, Elovich, and Kiselev models indicate the low predictive ability of these models to describe Cd sorption date (Table 2). Therefore, in the following, we will refer only to the results of the Langmuir isotherm model.
3.3.1. Langmuir Isotherm Model
Cadmium sorption conformed to the Langmuir model over the entire concentration range for all three particle sizes of WTR in both single- and multielement systems (Table 2 and Figure 3). The goodness-of-fit of Langmuir model suggested the homogeneous and monolayer mode of adsorption. The calculated Langmuir values for mWTR, uWTR, and nWTR in single-ion system were 2.8, 4.6, and 47 mg Cd g−1, respectively (Table 1). In multielement system, Langmuir values were 1.75, 2.50, and 16.67 mg Cd g−1 for mWTR, uWTR, and nWTR, respectively. The calculated Langmuir value of nWTR was 17 and 10 times higher than of mWTR in single- and multielement systems, respectively. Moreover, the high value of (Langmuir adsorption coefficient related to apparent energy of adsorption) showed that Cd was much better adsorbed onto the nWTR than mWTR and uWTR (Table 2). These results demonstrated that the nWTR is a superior sorbent for Cd in comparison to mWTR and uWTR. This could be explained on the basis of high and well-developed specific surface area of nWTR.
The competing effect of multi-ion system on Cd adsorption on nWTR was investigated by measuring the adsorption capacity ratio (R) of Cd using the following formula :where and are the maximum amount of Cd adsorbed by nWTR in multi-ion and single-ion systems, respectively, and is the adsorption capacity ratio. When , Cd adsorption is enhanced by the presence of other metal ions; when , no interaction occurred; when , Cd adsorption is suppressed by other metal ions. For all adsorbents studied, adsorption capacity ratios () of Cd were less than 1 indicating the suppression of Cd adsorption due to the competition between Cd and the two ions (Pb and Cu) in the multi-ion system for the same sorption sites.
3.4. Kinetic Studies
The results of kinetic studies on Cd adsorption at 3 different solution pH values are presented in Figure 4. The process of Cd removal using nanoscale WTR was rapid in the first 15 min and then slowed down considerably . At initial solution pH 5, approximately 70% of Cd was adsorbed by nWTR within the first 15 min and slowly proceeded to 100% sorption by the end of the 24 h period. As can be seen from Figure 4, Cd removal by nWTR decreased as the pH increased from 5 to 9 and a maximum value was reached at an equilibrium pH of around 5. In the aqueous solution, the metal ions may undergo hydrolysis and the metal adsorption process involved may present as follows :
pKa value for Cd(II) is 10.1 and M(OH)2 is the dominant species at pH > 6.0, while M2+ and M(OH)+ are the dominant species at pH < 6.0. On further increase of pH, adsorption decreases due to the formation and precipitation of cadmium hydroxide .
The kinetically driven Cd adsorption data at three pH solution values (5, 7, and 9) were fitted to four kinetic models  (Equilibrium and Kinetic Models) to predict Cd adsorption rate and understand the mechanism of adsorption process. The Cd adsorption data best fitted to power function and first-order kinetic models as the determination coefficients were quite high (significant at ) and SE values were the lowest (Table 3). The high SE values of Elovich and parabolic diffusion models indicated their low predictive capability to describe Cd adsorption by nWTRs. Therefore, based on SE values, power function model is considered the most suitable model to describe Cd adsorption on nWTR (Figure 5). The power function model is empirical and can successfully be used to compare experimental results. The estimated apparent sorption rates of the power function model (v) for Cd adsorption from single-ion and multi-ion solutions are listed in Table 3. The (v) values could be related to the empirical rate coefficients (k) of the overall sorption processes over the entire reaction time range . As can be seen in Table 2, the presence of the competing ions (Pb and Cu) greatly affected the estimated apparent sorption rate (v). The estimated values for Cd in single-ion solutions were consistently higher than for Cd in multi-ion solutions indicating competition between Cd, Pb, and Cu for nWTR sorption sites. The decreased Cd adsorption by nWTR in the multi-ion system can be explained by considering the properties the other elements present. Some properties of the metal ion are (i) electronegativity, (ii) charge to radius ratio, (iii) abilities to form hydroxo complexes, and (iv) preferred adsorption site on the adsorbent. In our multielement system (Cd-Pb-Cu system), Cd(II) adsorption was suppressed as Cu and Pb form the most stable monohydroxo complexes and the least soluble hydroxides which have tendencies to adsorb more preferably than Cd(II), thereby inhibiting Cd(II) adsorption [25, 26]. From these results, it is clear that Cd(II) adsorption behavior is significantly affected by the presence of competing ions (Pb and copper). Therefore, the competitive sorption should be considered to assess the potential removal of Cd from aqueous solutions.
3.5. Fourier Transmission Infrared (FTIR) Spectrum of nWTR
FTIR spectroscopic analysis was performed with FTIR spectrophotometer model 4400 (Shimadzu Corporation, Japan) to investigate the adsorption mechanism of Cd onto nWTR. Solid adsorbent was mixed with KBr in 1 : 10 ratio and pressed under 10-ton pressure into pellet and the infrared spectra were measured. The FTIR spectra of nWTR before and after Cd adsorption are demonstrated in Figure 6. Before Cd adsorption, the two bands at 4012 cm−1 and 3416 cm−1 on nWTR surface correspond to O-H bending vibrations, while the 1634 cm−1 band can be ascribed to the bending mode of H-O-H molecule [27, 28]. In addition, the FTIR spectrum showed a small band at 543 cm−1 (O-AL-O stretching vibration) and a medium band at 1090 cm−1 which can be assigned to bending vibration of FeOH modes of feroxyhyte [29, 30]. Retention of Cd on the surface of nWTR leads to significant spectral changes. The band at 4012 cm−1 (O-H bending vibration) completely disappeared which ascertained the participation of the surface hydroxyl group on the surface of nWTR (Figure 6). Increase in intensities and shift of locations of the 3416 cm−1 band to 3377 cm−1 (O-H bending vibrations), 1636 cm−1 to 1629 cm−1 (assigned to H2O bending vibrations), 1091 cm−1 to 1029 cm−1 (corresponding to bending vibration of FeOH modes of feroxyhyte), and 543 cm−1 to 529 cm−1 (O-AL-O stretching vibration ) are noticed. Meanwhile, two new bands indicative of d′-FeOOH appeared at 794 cm−1 and 683 cm−1 [31, 32]. The band position shift and intensities increase of the previously mentioned bands are clear indication of the level of specific molecular interaction. Therefore, the participation of OH, O-Al-O, FeOH, and FeOOH entities in the adsorption process is suggested.
3.6. Fractionation of Cd-Loaded nWTR and Cd Mobility
To evaluate the potential mobility of adsorbed Cd onto nWTR relative to mWTR and uWTR, the Cd distribution in the fractions of the Cd(II)-loaded adsorbents was performed using the fractionation scheme of Tessier et al. . The sequential extractions have been commonly used for Cd speciation because of the efficiency and reproducibility of the procedure . According to the scheme, the labile Cd associated with the residual (RS) fraction is less mobile than Cd associated with the nonresidual (NORS) fraction (sum of all fractions except the residual fraction). Such an approach enables us to correlate the Cd(II) data with adsorbent fractions and identify the mobility of Cd bound to different sorbent fractions. The Cd fractionation results (Figure 7) show that Cd adsorbed on mWTR was mostly associated with the more mobile NORS fraction (66%), whereas 78.4% of Cd on nWTR was associated with the RS fraction, the less mobile fraction. Such differences could be related to particle size and crystalline effects, which in turn affect the solubility of Cd phases, present . A greater percentage of Cd associated with the residual fraction of nWTR relative to mWTR and uWTR indicates the higher capability of nWTR to immobilize Cd compared to bulk WTR (mWTR). These observations support the hypothesis that the smaller size and the higher surface area of WTR nanoparticles greatly enhance its sorptive characteristics and stability.
3.7. Stability of Cd-nWTR Surface Complexes
The percentages of Cd desorbed from Cd-loaded nWTR are shown in Figure 8. The amounts of Cd released increased from 0.02 to 0.13% as the Cd concentrations increased from 5 to 500 mg/L (at pH 7.0). It is clear that the Cd adsorbed on nWTR was not altered as less than 0. 2% of adsorbed Cd(II) was released at the highest Cd concentration load after 4 consecutive desorption cycles. These results suggest the stability of Cd-nWTR surface complexes.
3.8. Efficiency of Cadmium (II) Removal from Wastewater
The efficiency of nWTR for Cd(II) removal was determined by conducting a batch experiment on real wastewater. Four samples were collected from wastewater discharged from fruit and vegetable processing plant by using polyethylene bottles. The wastewater was characterized in terms of pH, conductivity, Cd concentration, and major and minor elements. One liter of each water sample was spiked with 100 mg/L−1 of Cd+2 and 25 mL was placed in 50 mL falcon tube along with 0.1 g nWTR. The mixture was shaken for 2 h and centrifuged at 4000 rpm for 20 min and the supernatant solutions were filtered and analyzed for Cd by ICPS.
The pH and conductivity of the wastewater were 8.1 and ds m−1, respectively. The predominant anions in the water were chloride ( mg L−1), nitrate (.54 mg L−1), and phosphate ( mg L−1) and the sodium adsorption ratio (SAR) was . Amounts of Cu ( mg L−1), Cd ( mg L−1), Cr (), As ( mg L−1), and Pb ( mg L−1) were also detected. The real efficiency of nWTR for Cd removal from wastewater samples studied was calculated and found to be 98.35%. These results highlight the potential of nWTR for heavy metals removal from wastewater.
Nanoparticles of water treatment residuals demonstrated high removal efficiency of Cd(II) from aqueous solutions under various sets of experimental conditions. However, Cd(II) adsorption by nWTR was much higher in single-ion system than in multi-ion system due to the competition of Pb+2 and Cu+2 ions with Cd2+ ions. The Cd(II) adsorption mechanism seems to be related to the interaction between Cd(II) and OH, O-Al-O, FeOH, and FeOOH entities. The real efficiency of nWTR for Cd removal from wastewater samples was found to be 98.35%. Therefore, it is expected that the obtained nWTR can be used as potential sorbents for the removal of various toxic pollutants from wastewater.
The authors declare that there are no competing interests regarding the publication of this paper.
The authors gratefully acknowledge the financial support from the Egyptian Science and Technology Development Fund (STDF 4977-2013).
P. K. Rai, “Heavy metal pollution in aquatic ecosystems and its phytoremediation using wetland plants: an eco-sustainable approach,” International Journal of Phytoremediation, vol. 10, no. 2, pp. 131–158, 2008.View at: Google Scholar
G. Drasch, M. Horvat, and M. Stoeppler, “Mercury,” in Elements and Their Compounds in the Environment: Occurrence, Analysis and Biological Relevance, E. Merian, M. Anke, M. Ihnat, and M. Stoeppler, Eds., vol. 2, pp. 931–1005, Wiley-VCH, Weinheim, Germany, 2004.View at: Google Scholar
B. J. Alloway, Heavy Metals in Soils, John Wiley & Sons, New York, NY, USA, 1995.
O. B. Akpor and M. Muchie, “Remediation of heavy metals in drinking water and wastewater treatment systems: processes and applications,” International Journal of Physical Sciences, vol. 5, no. 12, pp. 1807–1817, 2010.View at: Google Scholar
C. H. Giles, T. H. MacEwan, S. N. Nakhwa, and D. Smith, “786. Studies in adsorption. Part XI. A system of classification of solution adsorption isotherms, and its use in diagnosis of adsorption mechanisms and in measurement of specific surface areas of solids,” Journal of the Chemical Society, vol. 111, pp. 3973–3993, 1960.View at: Google Scholar
G. Wulfsberg, Principles of Descriptive Chemistry, Brookes/Cole Publishing, Montery, Calif, USA, 1987.
M. Mohapatra, L. Mohapatra, P. Singh, S. Anand, and B. Mishra, “A comparative study on Pb(II), Cd(II), Cu(II), Co(II) adsorption from single and binary aqueous solutions on additive assisted nano-structured goethite,” International Journal of Engineering, Science and Technology, vol. 2, no. 8, 2011.View at: Publisher Site | Google Scholar
M. Gotić, S. Musić, S. Popović, and L. Sekovanić, “Investigation of factors influencing the precipitation of iron oxides from Fe(II) containing solutions,” Croatica Chemica Acta, vol. 81, no. 4, pp. 569–578, 2008.View at: Google Scholar
H. D. Ruan, R. L. Frost, J. T. Kloprogge, and L. Duong, “Infrared spectroscopy of goethite dehydroxylation: III. FT-IR microscopy of in situ study of the thermal transformation of goethite to hematite,” Spectrochimica Acta—Part A, vol. 58, no. 5, pp. 967–981, 2002.View at: Publisher Site | Google Scholar