Ding Hongbo1,2 Pan Zhongxiao1 and Renato Seeber2
1 Department of Applied Chemistry. University of Science and
Technology of China,Hefei,230026,China P.R. Email:
2 Department of Chemistry, University of Modena and Reggio Emilia,Modena,41100,Italy
The bistable model for the process of active-passive transition of iron in sulfuric acid was given. With this model, stochastic resonance phenomenon was simulated.
Keywords active-passive transition, stochastic resonance, bistability
Passivation of iron in sulfuric acid was first observed by Flade . Since the seventies, it has been shown that polarization curves for iron rotating disk electrodes, obtained with a potentiostat, in 1M sulfuric acid have displayed a hysteresis loop [2~4]. From the nonlinear dynamics point of view, this can be treated as possessing bistable characteristics.
The concept of stochastic resonance (SR) was first put forward in the seminal paper by Benzi and collaborators wherein they address the problem of the periodically recurrent ice ages. The basic ideas under the concept [5,6] are: for a given nonlinear system which possess the characteristic of bistability (or more generally, a form of threshold), when theres some coherence among the nonlinear condition of the system, the input signal and the input noise, an extra dose of noise can in fact help rather that hinder the performance of the system, a kind of phenomenon of the coherent effect between the stochastic force and the signal. Theres now even a term called signal-to-noise ratio (SNR) to quantify the effect. The concept of stochastic resonance has not only theoretical importance but potential application.
In the research area of electrochemical reactions, much attention has been attracted to nonlinear dynamics, such as oscillation and surface pattern. Among them, there are also some reports on electrochemical bistable systems [4,8,9]. In this paper, those experimental results from Epelboin et al was analyzed, and the bistable model for the process of active-passive transition for iron in sulfuric acid was described. With this model, stochastic resonance phenomenon was simulated. Through this way, the new nonlinear dynamics concept was introduced to electrochemists and this work might help researchers to find SR in this system experimentally.
Figure 1, which was taken from Epelboin et al, shows a hysteresis loop in the current-voltage curve. For a positive voltage sweep, the curve a-b-d-e-h will be obtained, while for a negative sweep, the curve h-e-f-b-a is obtained.
Fig.1 Current-voltage curve for a 5mm diam. disk in 1M sulfuric acid, rotating at 750rpm, taken from Epelboin et al
In the figure, the vertical line "bf" corresponds to Flade potential. From the theoratical point of view, both the positive and negative polarization curve should be the same, that is: a-b-f-h or h-f-b-a. Due to the effect of ohmic potential drop, the hysteresis loop was introduced. Therefore, for a certain cross-sectional rotating disk electrode, when the controlling parameter of V(polarization potential) are in the range of "f" and "e" , the polarization curve will be in two possible stable steady states "bd" and "fe"(This corresponds to the active and the passive state of the electrode system respectively). Therefore, from the nonlinear dynamics point of view, it can be viewed as a typical bistable state[10,11]. In certain time, the system will be in one certain state which is determined by the initial condition. With the experimental data, the differential equation describing the behavior of the system can be modeled as follows:
This differential equation possess two stable steady solution [10,11]: I=2 and I=0 (corresponds to active state and passive state respectively). In addition, it has also one unstable steady solution: I=1. However, due to the effect of the inherent stochastic process of electrode reaction [11,12], its impossible for the system to be in this state. Therefore, the system will only be in either the active state or passive state.
The mechanism of SR is simple to explain. Consider a heavily damped particle mass m and viscous friction , moving in a symmetric double-well potential V(x). The particle is subject to fluctuational forces that are, for example, induced by coupling to a heat bath. Such a model is archetypal for investigations in reaction-rate theory.The fluctuational forces cause transitions between the neighboring potential wells with a rate given by the famous Kramers rate  i.e.,
with being the squared angular frequency of the potential minima at , and the squared angular frequency at the top of the barrier, located at xb, V is the height of the potential barrier separating the two minima. The noise strength D=kB/T is related to the temperature T.
If we apply a weak periodic forcing to the particle, the double-well potential is tilted asymmetrically up and down, periodically raising and lowering the potential barrier. Although the periodic forcing is too weak to let the particle roll periodically from one potential well into the other one, noise induced hopping between the potential wells can become synchronized with the weak periodic forcing. This statistical synchronization takes place when the average waiting time TK(D)=1/rK between two noise-induced inter-well transitions is comparable with half the period TO of the periodic forcing. This yields the time-scale matching condition for stochastic resonance i.e.,
2 TK(D)= TO(3)
In short, stochastic resonance in a symmetric double-well potential manifests itself by a synchronization of activated hopping events between the potential minima with the weak periodic forcing. For a given period of the forcing TO, the time scale matching can be fulfilled by tuning the noise level Dmax to the value determined by Eq.(3).
In summary, the effect requires three basic ingredients, (i) an energetic activation barrier or, more generally, a form of threshold; (ii) a weak coherent input; (iii) a source of noise that is inherent in the system, or that ads to the input. Given these these features, the response of the system undergoes resonance-like behavior as a function of the noise level; hence the name stochastic resonance. The underlying mechanism is fairly simple and robust. As a consequence, SR has been observed in a large variety of systems, including chemical reactions .
SR can be envisioned as a particular problem of signal extraction from background noise. Its quite natural that a number of authors tried to characterize SR within the formalism of data analysis, most notably by introducing the notion of signal-to noise ratio (SNR). When noise amplitude fulfilled the coherent condition of equation (3), SNR will achieve its maximum value.
According to the theoretical discussion above, the adopted stochastic differential equation with this simulation work is:
Here, A refers to the amplitude of the input sinusoid current, w refers to angular frequency. H(t) is input noise current. In this simulation, after some modification of the parameters, the amplitude of sinusoid current and the angular frequency were chosen as A=0.38 and w=0.002, then theres only one variable: the amplitude of the noise. Here, the angular frequency of w was chosen with a very low value based on the fact that the relaxation time for the electrode system are relatively very long.
Fig.2 The optimum output of the system after modulating the noise amplitude to an optimum value
According to figure 2, when the noise amplitude was modified to an optimum value, the state of the system undergone periodic hopping with respect to every periodic signal. In this circumstance, the signal gives the system an optimum modulation.
In order to further discuss the relationship between the output of the system and the noise amplitude, the notion of signal-to-noise ratio (SNR) was employed. With the help of FFT technique, the time series output signal was translated into the frequency domain. For simplicity, here, the definition of SNR=S/N was taken. Figure 3 is the result of SNR as a function of the noise variance (analogous to noise amplitude).
Fig.3 SNR-H curve for the output of the system
According to figure 3, it can be seen that theres a peak at the position of the optimum value of noise. On the left hand of the peak, with the increment of noise amplitude, SNR increases; while on the right hand of the peak, with the further increment of noise amplitude, theres over output of noise, SNR decreases.
The appearances of hysteresis loop on the current-voltage curves for the process of active-passive transition of iron rotating disk electrode in 1M sulfuric acid can be described as a symmetric bistable system. Under this bistable model, the SR phenomenon can be simulated. This work might be help for further experimental verifications.
1. F.Flade, Z.Phys.Chem.,76,513(1911)
2. I. Epelboin, C. Gabrielli, M. Keddam, J.C. Lestrande and H. Takenouti., J. Electrochem. Soc., 119, 1632 (1972)
3. C.G. Law and J. Newman, J. Electrochem. Soc., 126, 2150 (1979)
4. P.P. Russell and J. Newman, J. Electrochem. Soc., 130, 547 (1983)
5. R. Benzi, G. Parisi, A. Sutera and A. Vulpiani., J. Phys. A, 14, 453 (1981)
6. L. Gammaitoni, P. Hanggi, P. Jung and F. Marchesoni, R. Mod. Phys., 70, 223 (1998)
7. J.L. Hudson and T.T. Tsotsis., Chem. Eng. Sci., 49, 1493 (1994)
8. D.Gilroy and B.E. Conway, Can. J. Chem., 46, 875 (1968)
9. Y.H. Xu and M. Schell, J. Phys. Chem., 94, 7137 (1990)
10. G. Nicolis and I. Prigogine, "Exploring Complexity", 1988 (Chinese version)
11. G. Hu, Stochastic Forces and Nonlinear Systems, Shanghai Scientific and Technological Education Publishing House, 1994
12. U. Bertocci and F. Huet, Corrosion, 51, 131 (1995)
13. P. Hanggi, P. Talkner and M. Borkorec, Rev. Mod. Phys, 62, 251 (1990)
14. H. Kramers, Physica (Utrecht), 7, 284 (1940)
15. L.F. Wang, C.S. Yin, Q. Ren and Z.X. Pan, Chinese J. Anal. Chem., 27 (1999, in press)
16. C.N.Cao, Principles for Metal Corrosion Electrochemistry, Publishing House of Chemical Industry, 1985 (in Chinese)