Simultaneous Determination of Zn ( II ) , Cd ( II ) , Pb ( II ) , and Cu ( II ) Using Differential Pulse Anodic Stripping Voltammetry at a Bismuth Film-Modified Electrode

e simultaneous electrochemical determination of Zn(II), Cd(II), Pb(II), and Cu(II) in the aqueous solution has been developed on the basis of the bismuth film glassy carbon electrode (GCE) using differential pulse anodic stripping voltammetry (DP-ASV). e bismuth film electrode (BiFE) was prepared by adding 500 ppb bismuth(III) directly to the sample solution and simultaneously depositing bismuth and the metal analytes on GCE. e optimal operational parameters, namely, accumulation potential (–1.6 V), accumulation time (110 s), pulsed amplitude (0.07 V), and scan rate (0.021 V·s), were found using a Box–Behnken design. Under the optimum conditions, a linear relationship exists between the current and the concentration of Zn(II), Cd(II), Pb(II), and Cu(II) in the range between 5.0 ppb and 110.0 ppb with the detection limits of 1.07 for Zn(II), 0.93 ppb for Cd(II), 0.65 ppb for Pb(II), and 0.94 ppb for Cu(II) calculated on the basis of a signalto-noise ratio equal to 3 (S/N � 3). e interference experiments show that Co(II), Ni(II), and Fe(III) have a little influence on the DP-ASV signals of Zn(II), Cd(II), Pb(II), and Cu(II). In addition, a high reproducibility was indicated from small relative standard deviations (1.03%, 1.74%, 1.32%, and 4.74%) for 25 repeated measurements of 15 ppb copper, lead, cadmium, and zinc solutions. BiFE was successfully applied to determine Zn(II), Cd(II), Pb(II), and Cu(II) in river samples, and the results are in a good agreement with those determined with graphite furnace atomic absorption spectrometry (GFAAS).


Introduction
Heavy metal pollution (Hg(II), Pb(II), Cd(II), and Ni(II)) has become one of the most critical environmental problems today.Unlike organic pollutants, heavy metals are nonbiodegradable.ey are found in industrial wastewater, and their discharge to the environment is alarming owing to their acute toxicity to aquatic and terrestrial life, including humans [1].Of these, excessive intake of copper (over 1.0 mg•L −1 in drinking water) accumulates in the livers of humans and animals, which causes hemochromatosis and gastrointestinal catarrh diseases [1].Lead is a highly toxic and cumulative poison and accumulates mainly in the bones, brain, kidneys, and muscles.Lead poisoning in humans causes severe damage to the kidneys, nervous and reproductive systems, liver, and brain [2,3].Cadmium (Cd) is one of the most toxic pollutants in aquatic systems.It cannot be removed easily from water systems by self-purification and could pose a serious threat to the environment and human health [4].e longtime accumulation (30 years) of Cd through the food chain can lead to serious osteoporosislike bone diseases [5].Zinc, an important element for living organisms, when exceeding the specific limits, accumulates in nature as it cannot be biodegraded.Zinc is associated with short-term "metal-fume fever," nausea, diarrhea, depression, lethargy, and neurological signs, such as seizures and ataxia [6,7].erefore, the determination of heavy metals in the human environment is very necessary and important.ere are several methods to determine heavy metals in water.e most widely used methods for heavy metal determination include atomic absorption spectrometry (AAS) [8], flame atomic absorption spectrometry (FAAS) [9], inductivelycoupled plasma mass spectroscopy (ICP-MS) [10], graphite furnace atomic absorption spectroscopy (GF-AAS), and X-ray fluorescence [11], and they are employed currently for trace analysis of heavy metals.Among these methods, anodic stripping voltammetry approaches are an emerging and attractive method because the main advantages of this technique can be ordered as high selectivity, sensitivity, in situ operation, and low cost.
In 2000, a new type of electrode, the bismuth filmmodified electrode (denoted as BiFE), was found as an alternative to mercury film electrodes (MFEs) [12], which are highly toxic.e significant advantage of BiFEs is that they are environmentally friendly since the toxicity of bismuth and its salts is negligible.However, the advantageous analytical properties of BiFEs in the voltammetric analysis, roughly comparable to those of MFEs, are assigned to the property of bismuth to form "fused alloys" with heavy metals, which is analogous to the amalgams that mercury form [13]. Wang et al. [12] reported that bismuth-coated carbon electrodes display an attractive stripping voltammetric performance that is similar to that of common mercury film electrodes and could be utilized to simultaneously determine Cu(II), Cd(II), and Pb(II) in aqueous solutions.Wu et al. [14] reported the simultaneous determination of Pb, Cd, and Zn using DP-ASV at a bismuth/ poly(p-aminobenzene sulfonic acid) film electrode.Lead and cadmium in rice samples are determined by means of ASV using a bismuth film/crown ether/Nafion-modified screenprinted carbon electrode [15].Li et al. [16] reported the simultaneous detection of ultratrace lead and cadmium by square wave stripping voltammetry with in situ depositing bismuth at the Nafion-medical stone-doped disposable electrode, and the assay results of heavy metals in wastewater with the proposed method were in a good agreement with the atomic absorption spectroscopy method.To the best of our knowledge, a few articles have attempted to study the simultaneous determination of Zn(II), Cd(II), Pb(II), and Cu(II) by BiFE.
is article demonstrates the attractive stripping voltammetric behavior of BiFE.BiFE was prepared by simultaneously accumulating bismuth and the target heavy metals (Cu, Pb, Cd, and Zn). e Box-Behnken design was conducted to optimise the operational parameters.BiFE displays well-defined, sharp, and highly reproducible stripping peaks for low (ppb) concentrations of copper, lead, cadmium, and zinc over a low background current.
e simultaneous determination of copper, lead, cadmium, and zinc in real samples was also addressed.

Experimental
2.1.Materials.All solutions were prepared using doubledistilled water.e zinc, cadmium, lead, copper, and bismuth standard stock solutions (1000 ppm Zn, Cd, Pb, and Cu or 5 (Bi) wt.% nitric acid) were obtained from Aldrich and diluted as required.An acetate buffer solution (0.10 M, pH 4.5) served as the supporting electrolyte.

Apparatus.
Stripping voltammetry was performed with a CPA-HH5 Computerized Polarography Analyzer (Vietnam).A glassy carbon electrode (GCE) (3 mm diameter, BAS) served as a working electrode, with Ag/AgCl (3 M KCl) and platinum wire acting as the reference and counter electrodes, respectively.Graphite furnace atomic absorption spectroscopy (GF-AAS) was performed in ZEEnit 700P (Analytik Jena, Germany) as a reference method.

Procedure.
e glassy carbon electrode was polished using 1.0 μm Al 2 O 3 slurry, rinsed with ethanol/water and dried.e analytes (Zn, Cd, Pb, and Cu) were accumulated on the surface of GCE at a potential of -1.4 V (E acc ) for 120 s of accumulation (t acc ).During this step, the electrode was rotated at a constant rate of 1000 rpm.After that, the electrode rotation was ceased; then, cyclic voltammograms (CVs) were recorded from -1.5 V to +0.1 V (forward potential scan) and then from +0.1 V to -1.5 V (reverse potential scan) at a scan rate of 0.2 V•s −1 .
Stripping voltammetric measurements were conducted with an in situ accumulation of the bismuth film and target metals (Zn, Cd, Pb, and Cu).
e three electrodes were immersed in a 10 mL electrochemical cell containing 0.1 M acetate buffer (pH 4.5) and an appropriate amount of bismuth.e accumulation potential at -1.4 V was applied to the carbon working electrode under stirring.After the accumulation step (120 s), stirring was stopped, and after 10 s, the voltammogram was recorded by applying a positivegoing differential pulse voltammetric potential scan (with an amplitude of 50 mV and a potential step of 6 mV).e scan was terminated at 0.1 V.All experiments were carried out at room temperature.

Operational Parameters of Box-Behnken Central Design (BBD).
e operational parameters were optimised by using BBD.e effect of accumulation potential (x 1 ), accumulation time (x 2 ), pulse amplitude (x 3 ), and scan rate (x 4 ) on the peak current of BiEF was studied.e number of experiments (N) required for the performance of this design is defined as N � 2k • (k − 1) + C 0 , where k is the number of factors (k � 4) and C 0 is the number of replicates at the center point [17].us, a total of 27 experiments were carried out for optimising these four variables.Each independent variable is considered at three levels: low, medium, and high, coded as -1, 0, and +1, respectively (Table 1).e center points were used to determine the error.On the basis of the experimental data, a second-order polynomial model is 2 Advances in Materials Science and Engineering obtained.
e relationship between the response and the variables could be presented as follows: where y is the predicted response value (anodic peak current, I pa ); x 1 , x 2 , x 3 , and x 4 are the independent variables; b 0 is the intercept term (the mean value of response at the center points);  1.

Sample Preparation.
One milliliter of the sample solution was adjusted to pH 5 with NaOH or HCl and diluted to 5 mL with a 0.1 M acetate buffer (pH 5), and then it was placed in the cell without deaeration.Stripping voltammetric measurements were performed with an in situ accumulation of bismuth and the target metals (Zn, Cd, Pb, and Cu), and the accumulation was carried out for 110 s at -1.6 V, while the solution was stirred.After the solution was kept standstill for 10 s, DP-ASV was recorded.e metal analytes in the river water were determined by means of GF-AAS for the sake of comparison.

Effect of Bismuth Concentration.
e influence of the bismuth ion concentration on the anodic current peak (I pa ) of the metals was studied in the range 0-1000 ppb for a solution containing 15 ppm Zn(II), Cd(II), Pb(II), and Cu(II) in a 0.1 M acetate buffer (pH 4.5) at bare GCE.Very small and distorted stripping peaks of the metal analytes are observed at bare GCE, and the signal of Zn is not obtained (Figure 1(a)).When bismuth (500 ppb) was added to the sample, distinct and undistorted stripping peaks for the analytes are achieved, indicating that bismuth is responsible for the attractive electrochemistry behavior.It is supposed that the amount of the Bi(III) ion decides the thickness of the Bi film, whereas the film thickness does not influence the peak position of any metals.Figure 1(b) presents the effect of Bi concentration on the anodic peak current.e anodic peak of cadmium increases significantly with the bismuth concentration from 0 to 400 ppb and then levels off.A similar trend is observed for Pb below 600 ppb, and the anodic peak of lead fluctuates slightly after this concentration.e copper and zinc peaks increase significantly with the bismuth concentration up to 500 ppb, and the rise becomes slight afterwards.
e results suggest that the appropriate combination of sensitivity, peak sharpness, and background current (especially close to the Zn peak) was performed for Bi(III) concentration of 500 ppb.

Effect of pH.
e effect of pH on the stripping peak current of Zn(II), Cd(II), Pb(II), and Cu(II) at BiFE was studied (Figure 2(a)).e optimal pH ranges are 5.0-6.0 for Zn(II) and Cd(II) and 4.0-6.0 for Pb(II) and Cu(II) (Figure 2(b)).If the acidity is higher or lower than the optimum ranges, the responses of Zn(II) and Cd(II) decrease.For Pb(II) and Cu(II) peak currents, I pa decreases slightly in pH from 4.0 to 6.0 and declines significantly with the further increasing pH.e best compromise of I pa could be obtained at pH 5.0-6.0.In the present study, we choose pH of 5.0 for the next experiments.
Figure 2(c) provides the linear plots of pH dependence of peak potential (E p ). e linear equations are expressed as follows: e linear relationship of E p vs. pH exhibits a high statistical significance with high determination coefficients (0.941-0.990, p < 0.05).e slopes of the lines are very close to the theoretical value of 0.059 for the Nerstian equation, indicating that the redoxidation involves the equal number of protons and electrons.

Effect of CV Scan Rate.
e relationship between the peak current and the scan rate could provide important information about the electrochemical mechanism.erefore, the effect of scan rate on the peak potential (E p ) and peak current (I p ) was investigated using CV, as shown in Figure 3(a).If the electrooxidation reaction is reversible, E p is independent on v, and vice versa.As can be seen from Advances in Materials Science and Engineering 3 4 Advances in Materials Science and Engineering Figure 3(a), the anodic peak potential (E p ) shifts to a higher potential with the increase in the scan rate, and then it is suggested that the electron transfer in analyte electrooxidation is irreversible.e anodic peak current (I pa ) increases with the increase in the scan rate from 0.2 to 0.6 V•s −1 , as shown in Figure 3(b), suggesting that the electron transfer reaction involves a surface-confined process.
e relationship between the peak potential (E pc ) and the natural logarithm of the scan rates is described by Laviron's equation [18]: where α is the charge transfer coefficient, k s is the heterogeneous electron transfer rate constant of the surface-confined redox couple, n is the number of electrons transferred, v is the scan rate (V•s −1 ), and E 0 is the formal redox potential, T � 298 K, R � 8.314 J•mol•K −1 , and e plots of E pa of the analytes versus ln v are shown in Figure 3(c).e linear regression equations are as follows: Advances in Materials Science and Engineering E pa,Cu � (-0.0210 ± 0.0001) +(0.024 ± 0.001) • ln ], R 2 � 0.997.
e value of α•n in equation ( 3) for Zn, Cd, Pb, and Cu can be derived from the slope of regression equations ( 4)-( 7), and they are 1.11, 1.02, 1.11, and 1.05, respectively.e value of α is assumed to be equal to 0.5 for an irreversible system [19].erefore, the number of electrons transferred (n) in the electrooxidation of Zn, Cd, Pb, and Cu is 2.23, 2.05, 2.23, and 2.10, respectively.Consequently, with n � 2, the oxidation mechanism for Zn, Cd, Pb, and Cu could involve two electrons and two protons at pH 5. e reactions at the electrode are as follows: In order to decide whether the electrooxidation reaction is controlled by adsorption or diffusion, the plots of peak current (I p ) against v 1/2 are drawn (Figure 3(d)).If the linear plot of I p vs. v 1/2 crosses the origin, the process is controlled by diffusion [20].In the range 0.2 to 0.6 V•s −1 , the linear regression equations of I pa for Zn, Cd, Pb, and Cu oxidation versus v 1/2 are as follows: I pa,Pb � (2.26 ± 0.02) +(0.19 ± 0.04 Although the plots of I p versus v 1/2 are linear with high determination coefficients (0.958-0.981, p ≥ 0.001), the intercepts do not cross the origin because the 95% confidence intervals for the intercepts do not contain zero (see the intercepts of equations ( 9)-( 12)). is means that the electrode process of analyte electrooxidation is not controlled by diffusion.

Optimising the Operational Parameters Using BBD.
A total of 27 experiments consisting of 24 standard cube points and 3 center points in the cube were carried out.e design of experiments and response values (I p : anodic stripping current) of Cu, Pb, Cd, and Zn corresponding to each run is listed in Table 2.
e analysis of the data using Minitab 16 gives the following regression equations for I p : e values of the coefficient of determination are 0.999, 0.998, 0.999, and 0.998 for Zn, Cd, Pb, and Cu, respectively, indicating an excellent agreement between experimental and predicted values.e positive and negative signs in each equation imply the synergistic and antagonistic effect of the variables.Although the magnitude of the coefficients in the equations varies, their sign is the same, indicating that the 6 Advances in Materials Science and Engineering impact of each variable on the intensity of the peak current has the same tendency.e tests for the significance of the regression model for the response of analytes were adjusted, and the ANOVA tests are presented in Table 3. e nonsignificant lack-of-fit (0.654 for Cu, 0.541 for Pb, 0.086 for Cd, and 0.132 for Zn) confirms the validity of these quadratic models in the present study.From the statistical point of view, the coefficient with the p value > 0.05 is insignificant. is means that their contribution to the peak currents is negligible, and these effects should be discarded from the model.us, the predictive models become as follows: For the sake of simplicity, the I p signal of Zn was used to optimise the experimental conditions because it is the lowest of those of the four metals at the same concentration.e profile for predicted values in Minitab 16 was performed for the optimisation process.
e optimisation design matrix (Figure 4) represents the maximal I pa,Zn (1.59 μA) at the optimal conditions: accumulation potential -1.6 V, accumulation time 110 s, pulse amplitude 0.07 V, and scan rate 0.021 V•s −1 .To evaluate the repeatability of the model, the experiment under optimal conditions was repeated three times, and the obtained I pa,Zn is 1.55, 1.57, and 1.60 μA. e one-sample t-test proves that there is no statistical difference among these values (t(2) � 1.147, p � 0.37 > 0.05).ese optimal operational parameters were utilized for further experiments.

Interference Study.
e interference study was performed by adding various foreign substances to a standard solution containing 15 ppb of Zn(II), Cd(II), Pb(II), and Cu(II) and 500 ppb Bi(III) in the 0.1 M acetate buffer (pH 5) at the accumulation potential of −1.6 V and accumulation time of 110 s. e ratios of interference to Zn(II), Cd(II), Pb(II), and Cu(II) for the ±5.0%signal change are 15 for Co(II), 15 for Ni(II), and 45 for Fe(III) (Tables S1-S3).

Precision and Limit of Detection (LOD). An appropriate precision is another attractive feature of bismuth electrodes.
A series of 25 repetitive measurements of a solution containing 15 ppb zinc, cadmium, lead, and copper results in highly reproducible stripping peaks with a relative standard deviation of 1.74%, 1.29%, 1.07%, and 0.76%, respectively, less than 1/2 RDS H (Table S4).Such a good precision is attributed to the reproducible film renewal accrued from the in situ bismuth film.Compared with the "stabilization" period required for in situ plating of mercury-film electrodes [21], BiFE displays a highly stable response starting with the first run.
To verify the practicality of BiFE for the simultaneous determination of Zn(II), Cd(II), Pb(II), and Cu(II), the stripping process of these four metal ions in the mixture was also investigated (Figure 5).In each measurement, the concentration of one species changes while keeping that of the others constant.It can be seen from Figure 5(a) that the peak current of Cu(II) is positively proportional to its concentration from 5 to 110 ppb when keeping the concentration of Zn(II), Cd(II), and Pb(II) constant; meanwhile, the response of Zn(II), Cd(II), and Pb(II) is practically unaltered with the increasing Cu(II) content (a slight increase in the peak heights of lead and cadmium with increasing zinc concentration is merely due to the additive signal contribution of the slightly overlapping copper peak).A similar situation is observed with Zn(II), Cd(II), and Pb(II) (Figures 5(b)-5(d)).No interference can be detected for the determination of Zn(II), Cd(II), Pb(II), or Cd(II) by the coexistence of the other three species.erefore, it can be concluded that, in the quaternary mixtures containing Zn(II), Cd(II), Pb(II), and Cu(II), the stripping peaks of the four analytes are clearly separated from each other.e simultaneous determination of Zn, Cd, Pb, and Cu was also performed using the DP-ASV method.e electrochemical signals of Zn, Cd, Pb, and Cu were recorded simultaneously when increasing the concentration of the analytes in the 0.1 acetate buffer solution with pH � 5.    Advances in Materials Science and Engineering very close to each other.is suggests that the simultaneous voltammetric determination of Cu, Pb, Cd, and Zn is possible at BiFE in the mixed samples without any cross interference.is linearity is also relevant for the determination of Zn, Cd, Pb, and Cu in water.
A comparison of the developed method with other voltammetric methods for Zn, Cd, Pb, and Cu determination is listed in Table 4.It could be noticed that LOD of Zn, Cd, Pb, and Cu from this method is lower or comparable with that obtained on the basis of modified electrodes in the literature.Overall, BiFE proves to be an effective electrode for the simultaneous determination of Zn, Cd, Pb, and Cu.

Practical Application.
Water samples taken randomly from five rivers in Quang Binh province, Vietnam, namely, Cau Rao, Kien Giang, Nhat Le, and Son and Gianh river were used to determine the concentration of copper, lead, cadmium, and zinc using the proposed method (DP-ASV) and GF-AAS for the sake of comparison.Table 5 lists the content of Zn, Cd, Pb, and Cu in the samples determined by using these two methods.e paired t-test was used to compare the data.With the significant level α � 0.05, the data show that the amounts of Zn, Cd, Pb, and Cu determined with the proposed method are not statistically different from those determined with GF-AAS (Cu: t(4) � 0.291; p � 0.785; Pb: t(4) � 0.347; p � 0.746; Cd: t(4) � 0.975; p � 0.385; Zn: t(4) � 1.545; p � 0.197).
e recovery measurements were also performed to evaluate the accuracy of the method.e data provide good average recoveries, ranging from 92% to 105% for Cu, 91% to 102% for Pb, 96% to 108% for Cd, and 92% to 106% for Zn, suggesting that the developed method does not suffer from any significant effects of matrix interference.

Conclusion
In this article, the differential pulse anodic stripping voltammetry was used to determine Zn, Cd, Pb, and Cu in the rivers at a bismuth film-modified electrode.Simultaneous determination of Zn, Cd, Pb, and Cu using DP-ASV on BiFE displayed low detection limit, excellent long-time stability, Advances in Materials Science and Engineering and high reproducibility.e limits of detection were found to be 1.07, 0.93, 0.65, and 0.94 ppb for Zn, Cd, Pb, and Cu, respectively.
e proposed method was successfully employed for the determination of Zn, Cd, Pb, and Cu in river water in an agreement with GF-AAS analysis.
b 1 , b 2 , b 3 , and b 4 are the linear coefficients; b 12 , b 13 , and b 14 and b 23 , b 24 , and b 34 are the cross-product coefficients; and b 11 , b 22 , b 33 , and b 44 are the quadratic-term coefficients.e design and analysis of the experiment were carried out using Minitab version 16. e definition of the factors and their level values are shown in Table

Figure 4 :
Figure 4: (a) Profiles for predicated values; (b) reliability function for peak current of Zn.

Table 1 :
Factors in BBD and their levels.

Table 2 :
Design matrix and responses for full factorial design.