A New Accurate Equation for Estimating the Baseline for the Reversal Peak of a Cyclic Voltammogram
A New Accurate Equation for Estimating the Baseline for the Reversal Peak of a Cyclic Voltammogram
Journal of Electrochemical Science and Technology. 2016. Dec, 7(4): 293-297
Copyright © 2016, The Korean Electrochemical Society
  • Received : October 12, 2016
  • Accepted : November 10, 2016
  • Published : December 31, 2016
Export by style
Cited by
About the Authors
Sung-Hoon Oh
Byoung-Yong Chang

Here we propose a new equation by which we can estimate the baseline for measuring the peak current of the reverse curve in a cyclic voltammogram. A similar equation already exists, but it is a linear algebraic equation that over-simplifies the voltammetric curve and may cause unpredictable errors when calculating the baseline. In our study, we find a quadratic algebraic equation that acceptably reflects the complexity included in a voltammetric curve. The equation is obtained from a laborious numerical analysis of cyclic voltammetry simulations using the finite element method, and not from the closed form of the mathematical equation. This equation is utilized to provide a virtual baseline current for the reverse peak current. We compare the results obtained using the old linear and new quadratic equations with the theoretical values in terms of errors to ascertain the degree to which accuracy is improved by the new equation. Finally, the equations are applied to practical cyclic voltammograms of ferricyanide in order to confirm the improved accuracy.
1. Introduction
Cyclic voltammetry is the most useful technique for electrochemical experiments because it can be carried out very easily and the cyclic voltammogram (CV) obtained therefrom includes rich information about the electrochemical reactions that are characterized by the relationship between the potential and the current. Even though the i-E (current-potential) curve appears monotonic, several important parameters are involved, as shown in the following equation of a reversible faradaic reaction,
PPT Slide
Lager Image
where i p is the peak current; F , the Faraday constant; R , the gas constant; T , the temperature; n , the number of electrons transferred; A , the electrode area; D ox ,the diffusion coefficient of an oxidized species; C O * , the concentration of the oxidized species; and ν , the scanrate as conventionally mentioned in a textbook [1] .
However, the drawback of such an information-rich and easy to perform technique is the laborious analysis required to resolve the plots that appear simple but are complex in reality. One of the examples is performing the same experiments whilst varying only the scan rate. Fig. 1 shows a good example of finding the diffusion coefficient of an oxidized species ( D ox ) from the slope of the cathodic peak currents in Fig. 1(b) [2] . In the same way, α , n , and bulk concentrations can be evaluated as demonstrated in the literature [3 - 4] .
PPT Slide
Lager Image
(a) Simulated cyclic voltammograms with increasing potential scan rates (in the direction of the arrow). Simulation conditions: k0=1.0 cm/s, α = 0.5, n=1, COx= 1 mM, CRed= 0 mM, DOx= 5×10−6cm2/s, scan rate = 10 mV/s − 3.2 V/s. (b) Plots of cathodic and anodic peak currents of CVs when the switching potentials (Eλ) are −0.2, −0.4, and −1.0 V vs. E0.
It is interesting that while the CV comprises forward and reversed potential scans, there are only a few studies on the reverse currents, except for those examining the reversibility of the electrochemical reaction of interest [5] . One of the possible reasons for this is the effect of the switching potential ( Eλ ) on the reverse current. Fig. 1(b) compares the anodic peak currents obtained with Eλ = −0.2 V, −0.4 V, and −1.0 V vs . E 0 , respectively. Even though the forward scans give the same peak current values, the reversal scans do not. Thus, we could speculate that all the plots have strong linear relationships with ν 1/2 , as seen in eq. (1); yet the exact values of parameters cannot be determined without information that would have to be obtained by conducting other experiments.
Nevertheless, an equation found in the literature [1 , 6] describes the ratios of the forward and reverse current peaks of a CV. When only the oxidized species (Ox) is present and its cathodic peak current ( i pc ) is generated during the forward linear potential sweep, the reversal anodic peak current ( i pa ) should be measured from the baseline of the decaying curve (refer to i pa in Fig. 5 ). A Nernstian wave with stable products will give i pa / i pc = 1, regardless of Eλ . However, obtaining a baseline for measuring i pa is not always possible. In such cases, the baseline can be estimated by the following equation:
PPT Slide
Lager Image
Even though eq. (2) provides a good approximation for measuring the reversal peak current, the linear algebraic equation over-simplifies the voltammetric curve, even if the i-E curve does not appear linear. This means that the tolerance of deviation should be taken into account.
Here, we propose a new equation that provides more accuracy by describing the current relationships through a quadratic algebraic equation. In order to acquire the equation, CV simulations were carried out by varying the parameters in eqs. (1) and (2), and by numerical analysis with varying currents.
2. Experimental Section
- 2.1 Simulation
The finite element method (FEM) was employed to conduct the digital simulation of CVs as has been described elsewhere [7 - 10] . To summarize, the molecular species Ox and Red and their concentrations C ox and C red were digitized as finite elements of distance ( d ) and time ( t ) to form C O ( d , t ) and C R ( d , t ). For a convenient calculation in MATLAB, they were framed in a 2D-matrix, and redefined as C O ( x , y ) and C R ( x , y ) [10] . Here, x and y are dimensionless parameters defined by x = d d , and y = t t , respectively, as is frequently done in electrochemistry simulations [1 , 7 - 10] .
On the electrode surface, the balance between the rate of charge transfer and the flux of Ox and Red is calculated by the following equation with k f and k b as the rate constants of the forward and backward reactions, and D M , a dimensionless diffusion coefficient calculated by DΔt / Δx 2
PPT Slide
Lager Image
In the bulk solution, Ox and Red molecules diffuse in and out according to the concentration gradient. It is formulated by Fick’s laws and calculated in the matrix as shown below:
PPT Slide
Lager Image
Finally, the current was calculated by re-dimensionalizing the flux, J ( y ).
- 2.2 Fitting procedure
The quadratic equation was obtained by fitting the calculated currents to a quadratic equation. The fitting was done by a pre-defined function in Matlab. The code is as follows:
P = polyfit( i sp,0 / i pc , i ap,0 / i pc , 2).
- 2.3 Cyclic voltammetry experiments
A conventional 3-electrode electrochemical cell was used. The working, counter, and reference electrodes were Au disk, Pt mesh, and Ag|AgCl, respectively. The electrolyte solution consisted of 5 mM ferricyanide (Sigma Aldrich), 100 mM KNO 3 (Sigma Aldrich) as the supporting electrolyte, and doubly distilled water as the solvent. Linear sweep and cyclic voltammetry experiments were carried out using SP200 (Bio-Logic). The voltammograms were recorded vs . a home-made Ag|AgCl reference electrode at a scan rate of 20 mV/sec.
3. Results and Discussion
Fig. 2(a) shows the change in the reverse current with switching potential ( Eλ ). As Eλ increases (negatively), the currents at Eλ ( i sp,0 ) decay downwards along the curve. Meanwhile, the uncorrected anodic peak current ( i pa,0 ) increases in the anodic direction. Both i sp,0 and i pa,0 are normalized by the cathodic peak current ( i pc ) and are plotted in Fig. 2(b) . Thus, fitting to a quadratic algebraic equation provides a good result.
PPT Slide
Lager Image
(a) Simulated cyclic voltammograms with increasing Eλ from −0.1 to −0.9 V vs. E0. The definitions of isp,0 and ipa,0 are graphically presented for easy understanding. (b) Plots of isp,0 and ipa,0 (normalized by ipc) are fitted to a quadratic equation (red solid line) by eq. (5).
PPT Slide
Lager Image
Additional simulations with different diffusion coefficients give the same results. If the baseline at the potential of the anodic peak for measuring i ap is i ap,1 , i ap = i ap,0 + i ap,1 , then i ap,0 is replaced with i ap i ap,1 in eq. (5). If products are stable and no additional reaction is involved, i pa / i pc = 1, regardless of E λ . Finally, we reach an equation that can calculate the virtual baseline current ( i ap,1 ) from the current at the switching potential ( i sp,0 ), irrespective of when the potential scan is reversed, as long as the current decays sufficiently, which is assumed to be around 70 % of the peak current of the forward sweep.
The new equation is compared with the existing one, eq. (2). Fig. 3 shows the baseline currents calculated according to i sp,0 / i pc changing from 0.70 to 0.00. In a real system, the potential sweep can be applied within a range of −2 V to +2 V vs. E 0 in order to avoid interfering with faradaic reactions such as electrolysis of the solvent. Thus, we can practically consider the variance of i sp,0 / i pc from 0.14 to 0.70. In the figure, the calculated values by both equations are almost in agreement with each other for i sp,0 / i pc from 0.45 to 0.70 which corresponds to 210 mV to 100 mV away from E 0 . However, when the potential shifts further, the old formula deviates from the theoretical values while the new one does not. Therefore, the new equation is confirmed to provide more accurate virtual baseline currents for measuring the anodic peak current.
PPT Slide
Lager Image
The red and blue lines are calculated by eq. (2) and (5), respectively, with isp,0 changing from 0.0 to 1.0. The filled circles are obtained from the CV simulations by increasing Eλ from −0.1 to −2.0 V vs. E0.
The anodic peak current estimated using the equations are compared with those measured using theoretically simulated cyclic voltammograms. Fig. 4 plots the values of i pa,eq(2) / i pa and i pa,eq(5) / i pa . Here, i pa,eq(2) and i pa,eq(5) are the estimated values using eqs. (2) and (5), respectively, and i pa is the theoretical value. In other words, the estimated currents are normalized by the theoretical currents. Both equations provide good approximations, but when the over-potential increases, eq. (2) shows an increasing deviation from the theoretical values. Meanwhile, eq. (5) keeps i pa,eq(5) / i pa = 1 even though the potential increases within the tolerance at less than ± 0.5%. This confirms that the new equation provides a more accurate baseline for measuring the peak current of the reverse scanning curve.
PPT Slide
Lager Image
ipa,eq(2) and ipa,eq(5) are ipa measured from the baseline currents calculated by eqs. (2) and (5) are normalized by ipa,theo . They are plotted along the values of Eλ. The shaded zone covers ± 0.5% of relative error from the true value (1.00).
A practical application of this method is demonstrated with an aqueous solution of ferricyanide and the experimental results are shown in Fig. 5 . Curve (a) is made by a linear sweep of potential from 0.30 to −0.35 V vs . Ag|AgCl and curve (b) is recorded with a cyclic potential scan switched at 0.00 V vs . Ag|AgCl. Curve (c) is the baseline current extracted from the current of curve (a) decaying continuously from the switching potential of curve (b). It is observed that until the potential reaches the switching potential (0.00 V), curves (a) and (c) are identical, which confirms the validity of curve (c) as a baseline for curve (b), though it is estimated from curve (a). Eventually, the value of i pa / i pc is measured to be 1.008 which is slightly higher than the ideal value of 1.000 for a reversible faradaic reaction. At the same time, the value of i pa / i pc is measured with the virtual baseline current to be 1.007 which is almost the same as the measured one. Both values have < 1% error and the error may be caused by the non-faradaic current such as the charging current for the electric double layer. On the other hand, the old equation calculates 1.014 which is very close to 1.007 but not as accurate as 1.008.
PPT Slide
Lager Image
(a) A linear sweeping voltammogram and (b) a cyclic voltammogram an aqueous solution of 5.0 mM ferricyanide at a scan rate of 20 mV/s. (c) a baseline curve made by the mirror image of curve (a) around Eλ of curve (b). ipa is measured from the baseline curve to the uncorrected peak current as shown by the arrow.
4. Conclusion
We have explained how a new equation was developed from CV simulations and confirmed its validity with experimental results. The new equation describes that the currents at the scan-switching potential and the reversal peak potential are related to each other through a quadratic algebraic equation, providing a virtual baseline current for correcting the measured peak current to a true value. Both simulation and experimental results confirm that the correction is quite accurate; the error is less than ± 0.5%.
The proposed method will be successfully applied to cases where making a baseline is impossible because a high over-potential may activate interfering reactions such as the electrolysis of solvents or supporting electrolytes. In those cases, the potential sweep should be switched back to an appropriate potential; even then acquisition of the baseline current curve by the further potential may become experimentally impossible. At such moments, our equation can provide a virtual baseline to approximate the true reversal peak current. Our proposed method is expected to be valuably utilized in such restricted experimental conditions.
This research was supported by a Research Grant of Pukyong National University (2015).
Bard A.J. , Faulkner L.R. (2001) Electrochemical Methods: Fundamentals and Applications Wiley New York
Kamel M. K. 2014 J. Electrochem. Sci. Technol. 5 (1) 23 - 31    DOI : 10.5229/JECST.2014.5.1.23
Rashid M. , Jun T.-S , Kim Y.S. 2014 J. Korean Electrochem. Soc. 17 (1) 18 - 25    DOI : 10.5229/JKES.2014.17.1.18
Siadat S. O. R. 2015 J. Electrochem. Sci. Technol. 6 (4) 111 - 115    DOI : 10.5229/JECST.2015.6.4.111
Whitson P.E. , VandenBorn H. W. , Evans D. H. 1973 Anal. Chem. 45 (8) 1298 - 1306    DOI : 10.1021/ac60330a016
Nicholson R. S. 1965 Anal. Chem. 37 (11) 1351 - 1355    DOI : 10.1021/ac60230a016
Rudolph M. , Reddy D.P. , Feldberg S.W. 1994 Anal. Chem. 66 (10) 589A - 600A    DOI : 10.1021/ac00082a725
Gosser D. K. (1993) Wiley-VCH New York
Streeter I. , Compton R. G. 2008 J. Phys. Chem. C 112 (35) 13716 - 13728    DOI : 10.1021/jp804442m
Chang B.-Y. 2015 J. Electrochem. Sci. Technol. 6 (4) 146 - 151    DOI : 10.5229/JECST.2015.6.4.146