Abstract

We present an alternative representation of integer and fractional electrical elements in the Laplace domain for modeling electrochemical systems represented by equivalent electrical circuits. The fractional derivatives considered are of Caputo and Caputo-Fabrizio type. This representation includes distributed elements of the Cole model type. In addition to maintaining consistency in adjusted electrical parameters, a detailed methodology is proposed to build the equivalent circuits. Illustrative examples are given and the Nyquist and Bode graphs are obtained from the numerical simulation of the corresponding transfer functions using arbitrary electrical parameters in order to illustrate the methodology. The advantage of our representation appears according to the comparison between our model and models presented in the paper, which are not physically acceptable due to the dimensional incompatibility. The Markovian nature of the models is recovered when the order of the fractional derivatives is equal to 1.

1. Introduction

Electrochemical Impedance Spectroscopy (EIS) is widely used to investigate the interfacial and bulk properties of materials, interfaces of electrode-electrolyte, and the interpretation of phenomena such as electrocatalysis, corrosion, or behavior of coatings on metallic substrates. This technique relates directly measurements of impedance and phase angle as functions of frequency, voltage, or current applied. The stimulus is an alternating current signal of low amplitude intended to measure the electric field or potential difference generated between different parts of the sample. The relationship between the data of the applied stimulus and the response obtained as a function of frequency provides the impedance spectrum of samples studied [1]. Transfer function analysis is a mathematical approach to relate an input signal (or excitation) and the system’s response. The ratio formed by the pattern of the output and the input signal makes it possible to find the zeros and poles, respectively. To analyze the behavior of the transfer function in the frequency domain, several graphical methods were used, such as Bode plots that provide a graphical representation of the magnitude and phase versus frequency of the transfer function and the Nyquist diagrams that are polar plots of impedance modulus and phase lag. It is very common in the literature to analyze impedance results by a physical model; this model is expressed by a mathematical representation and usually is represented by equivalent electrical circuits composed of elements such as resistors, capacitors, inductors, the Constant Phase Element (CPE), and Warburg elements (which is a special case of the CPE). Although the resulting model is not necessarily unique, it describes the system with great precision in the range of frequencies studied [1].

Fractional calculus (FC) is the investigation and treatment of mathematical models in terms of derivatives and integrals of arbitrary order [25]. In the last years, the interest in the field has considerably increased due to many practical potential applications [613]. In the literature, a number of definitions of the fractional derivatives have been introduced, namely, the Hadamard, Erdelyi-Kober, Riemann-Liouville, Riesz, Weyl, Grünwald-Letnikov, Jumarie, and the Caputo representation [25]. For example, for the Caputo representation, the initial conditions are expressed in terms of integer-order derivatives having direct physical significance [14], this definition is mainly used to include memory effects. Recently, Caputo and Fabrizio in [15] present a new definition of fractional derivative without a singular kernel; this derivative possesses very interesting properties, for instance, the possibility to describe fluctuations and structures with different scales. Furthermore, this definition allows for the description of mechanical properties related to damage, fatigue, and material heterogeneities. Properties of this new fractional derivative are reviewed in detail in [16].

Measurements of properties of materials, interfaces of electrode-electrolyte, corrosion, tissue properties of protein fibers, semiconductors and solid-state devices, fuel cells, sensors, batteries, electrochemical capacitors, coatings, and electrochromic materials have shown that their impedance behavior can only be modeled by using Warburg elements in conjunction with resistors, dispersive inductors, or CPE elements [17]. In this context, FC allows the investigation of the nonlocal response of electrochemical systems, this being the main advantage when compared with classical calculus. Some researches concerning EIS introduce FC; for example, the authors of [18] study the electrical impedance of vegetables and fruits from a FC perspective; the experiments are developed for measuring the impedance of botanical elements and the results are analyzed using Bode and Nyquist diagrams. In [19] a technique is presented to extract the parameters that characterize a dispersion Cole-Cole impedance model. Oldham in [20] describes how fractional differential equations have influenced the electrochemistry. Other applications of fractional calculus in electrochemical impedance are given in [2126].

Unlike the work of the authors mentioned above, in which the pass from an ordinary derivative to a fractional one is direct, Gómez-Aguilar et al. in [12] analyze the ordinary derivative operator and try to bring it to the fractional form in a consistent manner. Following this idea we present an alternative representation of integer and fractional electrical elements in the Laplace domain for modeling electrochemical systems represented by equivalent electrical circuits; the order of the fractional equation is . In this representation an auxiliary parameter is introduced; this parameter characterizes the existence of the fractional temporal components and relates the time constant of the system.

The paper is organized as follows: Section 2 explains the basic concepts of the FC, Section 3 presents the examples considered and the interpretation of typical diagrams, and the conclusions are given in Section 4.

2. Introduction to Fractional Calculus

The use of Caputo Fractional Derivative (CD) in Physics is gaining importance because of the specific properties: the derivative of a constant is zero and the initial conditions for the fractional order differential equations can be given in the same manner as for the ordinary differential equations with a known physical interpretation [4].

The CD is defined as follows [4]:where is a CD with respect to , is the order of the fractional derivative, and represents the gamma function.

The Laplace transform of the CD has the following form [4]:

The Caputo-Fabrizio fractional derivative (CF) is defined as follows [15, 16]:where is a CF with respect to and is a normalization function such that ; in this definition the derivative of a constant is equal to zero, but, unlike the usual Caputo definition (1), the kernel does not have a singularity at .

If and , the CF fractional derivative, , of order , is defined by

The Laplace transform of (3) is defined as follows [15, 16]:For this representation in the time domain it is suitable to use the Laplace transform [15, 16].

From this expression we have

3. Electrochemical Impedance and the Interpretation of Typical Diagrams

Regarding the equivalent circuits, there is a diversity of models used. The most common is to adjust the system to a simple model or one that includes a bilayer structure, either with RC groups in parallel or in series model, although there are cases where it is appropriate to include the Warburg impedance element to consider possible diffusive processes on the surface. Generally, using a complex equivalent circuit is not necessary to obtain a good characterization of the real system; commonly a simple electrical circuit is the first choice and increases its complexity when knowledge of the electrochemical behavior of the system also increases.

The Cole impedance model is based on replacing the ideal capacitor in the Debye model [27, 28]. Cole model is represented by a series resistor , a capacitor , and a resistor in parallel . In general it reflects the electrical resistance of the interface sample-electrode and maintains a negligible value with respect to ; is the order of the power that best fits the model obtained, , giving an ideal capacitor when it is . Using the algebraic representation of the circuit can be said to represent the total impedance aswhere .

On electric structures RC type, bias resistor represents the charge transfer resistance and capacitance of double layer; the CPE is a component that models the behavior of a double layer capacitor in actual electrochemical cells, that is, an imperfect capacitor, and the impedance is represented asEquation (8) describes the deviation from ideal capacitors. Considering and the constant (the inverse of capacitance), this equation describes a capacitor. For a CPE, the exponent is less than one [29].

In electrochemical systems the diffusion can create an impedance called Warburg impedance [26], commonly used to describe phenomena such as diffusion, adsorption, or desorption of electroactive substances at interfaces metal/coating; this impedance depends on perturbation frequency: at high frequency a small Warburg impedance results and at low frequency a higher Warburg impedance is generated. The Warburg impedance parameter indicates the existence of diffusive processes that can be related to the release of the dissolved species. This parameter only reports the blocking ability of the passive layer, so it is not possible to know the nature of the species which are dissolved by electrochemical impedance spectroscopy technique. The equation for the infinite thickness of the Warburg impedance is given bywhere is a Warburg coefficient. On a Nyquist plot the infinite Warburg impedance appears as a diagonal line with a slope of ; on a Bode plot, the Warburg impedance exhibits a phase shift of 45° [29].

Equation (9) is valid if the diffusion layer has an infinite thickness. Quite often this is not the case. If the diffusion layer is bounded, the impedance at lower frequencies no longer obeys (9). For the Warburg impedance with finite thickness, we get the formwhere is the Nernst diffusion layer thickness; is an average value of the diffusion coefficients of the diffusing species. This equation is more general and is called finite Warburg [29].

3.1. Equivalent Representation Based on Fractal Capacitors

One of the problems of the fractional representation is the correct sizing of the physical parameters involved in the differential equation, to be consistent with dimensionality and following [12] we introduce an auxiliary parameter in the following way:orwhere is an integer and when . Expressions (11) and (12) become a classical derivative; the auxiliary parameter has the dimension of time (seconds). This nonlocal time is called the cosmic time in the literature [30]. Another physical and geometrical interpretation of the fractional operators is given in Moshrefi-Torbati and Hammond [31]. Parameter characterizes the fractional temporal structures (components that show an intermediate behavior between a conservative system and dissipative; such components change the time constant of the system) of the fractional temporal operator [32]. In the following we will apply this idea to construct the fractional equivalent circuits and examples are analyzed.

3.2. Polarizable Electrode

The model of polarizable electrode, also known as faradaic reaction, provides a simple description of the impedance of an electrochemical reaction on electrode surfaces. The equivalent circuit is represented by Figure 1.

Considering initial conditions equal to zero, the equivalent impedance is found by the following equation in the complex frequency domain:Applying Kirchhoff laws to the circuit of Figure 1, we haveBefore applying the Laplace transform of (15) the following considerations must be taken into account:Consider (11); in the Caputo sense (1), (17) becomesSubstituting (19) into (15), we obtainApplying the Laplace transform (2) to (21), we obtainFinally from (23) we have the fractional impedance of the circuitwhere .

Consider (12), in the Caputo-Fabrizio sense (3), (17) becomesSubstituting (26) into (15), we obtainApplying the Laplace transform (5) to (28) we obtainFinally from (30) we have the fractional impedance of the circuitwhere and represent the fractional components of the system. Equations (24) and (31) are the result of applying the fractional temporal operator of Caputo and Caputo-Fabrizio type in (17) for the current in the capacitor; this general representation includes an arbitrary constant, , which can be considered its own electrochemical parameter. In the particular case of , (24) is reduced to the Cole model (7). Equations (24) and (31) presented here preserve the dimensionality of the studied system for any value of the exponent of the fractional derivative. On the other hand, if in (24) and (31), we obtain an ideal RC circuit and the Markovian nature of the model is recovered.

Consider that the values of the parameters of the circuit shown in Figure 1 correspond to , , and and . Figures 2(a), 2(b), 2(c), and 2(d) show the Nyquist and Bode plots for (24).

Now consider that the values of the parameters of the circuit shown in Figure 1 correspond to , , and and . Figures 3(a), 3(b), 3(c), and 3(d) show the Nyquist and Bode plots for (31).

Several studies use equivalent circuits considering pure capacitances for adjusting the impedance spectra and thus describe phenomena as deterioration of materials due to its porosity or the effects of exposure and surface preparation of substrates on the impedances [3336]. The electrical capacitance has information on the conductive properties of materials and their chemical composition. Recently, several authors have replaced this pure capacitance by a CPE, because this has been considered an improvement in the settings of the theoretical equations regarding experimental results [3739], for example, in predicting the useful life of coatings and surface roughness heterogeneity resulting from the presence of impurities, fractality, dislocations, adsorption of inhibitors, or formation of porous layers [4047]. Replacing the capacitor by a CPE adds a pseudocapacitive constant to the circuit and causes the attenuation of the imaginary part of the impedance causing a flattening of the semicircles in the Nyquist diagrams (this is their main feature); see Figures 2(a) and 2(c) for the Caputo approach and Figures 3(a) and 3(c) for the Caputo-Fabrizio approach. These figures allow seeing the effect of varying the values of the pseudocapacitive constants: in the range we have the behavior of a CPE and when the values correspond to pure capacitances. In this range, the order of the derivative transforms the capacitor into a CPE. When increasing the value of the pseudocapacitive constants, this also increases the imaginary part of the impedance and this changes the time constants of the system [32], which correspond to creating fractal structures for particular values of . In the Caputo-Fabrizio approach it is noted that the value of the resistance series is modified, showing the existence of heterogeneities in this component. In Figures 2(b) and 2(d) for the Caputo approach and Figures 3(b) and 3(d) for the Caputo-Fabrizio approach, the Bode plots exhibit changes in the cutoff frequency, hence the phase shift and the decrease of the magnitude. For the range of frequency 10 to 10000 rad/s a shift in the magnitude and phase is shown, which implies that the proposed circuit better defines these frequencies. We see that, by increasing the frequency from the cutoff frequency, the magnitude decreases, while for frequencies below the cutoff magnitude it is almost constant. This means that the current through the circuit is very high with decreasing frequency and is very low when the frequency increases. The current flowing through these resistors at low frequencies presents dissipative effects that correspond to the nonlinear situation of the physical process (realistic behavior that is nonlocal in time), for example, the ohmic friction, which raises the temperature and therefore the kinetic energy of the molecules of the system [48]. Concerning the phase, we have that by increasing frequency the displacement current and polarization are very much increased. In this context, the need for CPE has been attributed by some authors to the presence of roughness, corrosion products, changes in the morphology of the material, or heterogeneous surfaces.

3.3. Heterogeneous Reaction

This model describes a heterogeneous reaction which occurs in two stages with absorption of intermediates products and absence of diffusion limitations. This circuit has been commonly used to model a porous electrode coating or a defective electrolyte interface and has recently been applied in the assessment of passive metal-electrolyte or metal-electrolyte interfaces with hard coating [49, 50]. The equivalent circuit is represented by Figure 4.

This model has also been used in the study of coated metals. In this case, represents the capacitance of the coating; this value is smaller than a capacitance represented by a CPE; represents the pore resistance; it is the resistance of ion conducting paths developed in the coating. These paths may be physical pores filled with electrolyte; this electrolyte solution can be very different than the bulk solution outside of the coating; (double layer capacitance) and (charge transfer reaction) represent the interface between this pocket of solution and the bare metal.

Following the same methodology from previous example, the fractional impedance in the Caputo sense of this circuit is

In the Caputo-Fabrizio sense the fractional impedance is given bywhere , in (32) and (33), and represents the fractional components of the system. This general representation includes an arbitrary constant, , which can be considered its own electrochemical parameter. In the case of the physical parameters involved are and for they are , the time constant of the system. Equations (32) and (33) presented here preserve the dimensionality of the studied system for any value of the exponent of the fractional derivative; when , we recover the Markovian nature of the model.

Consider that the values of the parameters of the circuit shown in Figure 4 correspond to , , , , , and . Figures 5(a), 5(b), 5(c), and 5(d) show the Nyquist and Bode plots for (32).

Now consider that the values of the parameters of the circuit shown in Figure 4 correspond to , , , , , and . Figures 6(a), 6(b), 6(c), and 6(d) show the Nyquist and Bode plots for (33).

Some authors make changes to the ideal capacitors including elements of Warburg [49, 50]; these elements are represented in our circuit maintaining the order of the derivatives fixed at a value of . For Caputo and Caputo-Fabrizio approach, Figures 7(a) and 7(c), respectively, show the resulting Nyquist diagrams considering the order of the derivative for both capacitors at . From these figures it can be seen that for high frequency the impedance is small because the reagents must move away from the surface. The low-frequency disturbances allow the movement of distant molecules.

The circuit shown in Figure 4 has also been used with success in the description of titanium alloys, considering carbon coatings, silicon oxide, and titanium oxide; in such applications it is generally considered that the coating presents a porous outer sublayer and a denser internal sublayer [49, 51]. In Figure 7(c) it is observed that the resistance value decreases, which can provide valuable information on the characteristics of a coating film.

3.4. Limited Diffusion, Warburg Element

This model describes the polarization of an electrode considering limiting the diffusion. This circuit models a cell where polarization is due to a combination of kinetic and diffusion processes, is the resistance of the electrolyte solution, is the imperfect capacitor, is the electron-transfer resistance, and is the Warburg element due to diffusion of the redox couple to the interface from the bulk of the electrolyte [52]. The equivalent circuit is represented by Figure 8.

Following the same methodology from previous example, the fractional impedance in the Caputo sense of this circuit is

In the Caputo-Fabrizio sense the fractional impedance is given bywhere , in (34) and (35), and and represent the fractional components of the system. This general representation includes an arbitrary constant, , which can be considered its own electrochemical parameter. In the case of the physical parameters involved are for and for are , the time constant of the system. For (34) and (35) when , we recover the Markovian nature of the model.

The capacitor and Warburg element shown in Figure 8 are replaced by two fractional capacitors and , respectively. The capacitor is assigned an exponent , causing replacement of pure capacitance and causing thus an impedance according to an infinite diffusion layer or Warburg impedance. This element arises from one-dimensional diffusion of an ionic species to the electrode [18]. The capacitance of fractional order can be modeled in the range , causing that this pure capacitance represents a CPE or imperfect capacitance. Equations (34) and (35) presented here preserve the dimensionality of the studied system for any value of the exponent of the fractional derivative. The equivalent circuit shown in Figure 8 has been used to model phenomena of corrosion stability in the physiological environment such as the simulated body fluid; represents the electrolyte resistance, represents the coating pore resistance, and the capacitances represent the frequency-dependent electrochemical phenomena, such as the coating impedance and passive oxide film capacitance. CPE is used in these models to compensate the nonhomogeneity in the system [53].

Consider that the values of the parameters of the circuit shown in Figure 8 correspond to , , , , and . Figures 9(a), 9(b), 9(c), and 9(d) show the Nyquist and Bode plots for (34).

Now consider that the values of the parameters of the circuit shown in Figure 8 correspond to , , , , and . Figures 10(a), 10(b), 10(c), and 10(d) show the Nyquist and Bode plots for (35).

Replacing capacitor by a CPE introduces a pseudocapacitive constant at circuit and causes the attenuation of the imaginary part of the impedance causing a flattening of the semicircles in the Nyquist diagrams; see Figure 9(a) for the Caputo approach and Figure 10(a) for the Caputo-Fabrizio approach. Increasing the value of the pseudocapacitive constant implies an increase in the contribution of the imaginary part of the impedance and this changes the system’s time constants [32], creating fractal structures which correspond to particular values of . In Figures 9(b), 9(d), 10(b), and 10(d) for the Caputo and Caputo-Fabrizio approach, respectively, the Bode plots show variations in the cutoff frequency, hence the phase shift and the decrease of the magnitude.

If , then we have a model that geometrically describes the activation surface inhomogeneities or deviations from linear diffusion processes. In our example . This occurs naturally when the diffusion occurs in a dilute solution or in the case that the spread does not obey the laws of Fick; see Figures 9(c) and 10(c).

4. Conclusions

An impedance spectrum which is obtained in response to the small amplitude signal excitation is often interpreted in terms of an equivalent electrical circuit. This one is based on a physical model that may represent and characterize elements whose electrochemical properties and structural features or the physicochemical processes are taking place in the studied being. The resulting spectrum described in a wide frequency range might present one, two, or more time constants, depending on the monolayer or multilayer structure, porosity, or diffusive limitations caused by the charge transfer process. The existence of a time constant indicates a homogeneous layer structure; if two different constants are presented, this may indicate the existence of two sublayers.

In the field of biomaterials, the electrochemical impedance spectroscopy technique proves to be useful in characterizing roughness or heterogeneous surfaces, in coatings analysis, as cell suspensions, in studying the adsorption of protein, and in characterizing the performance of materials, mainly with regard to electrocatalysis and corrosion. In this context, FC allows the investigation of the nonlocal response of electrochemical systems. FC has been used successfully to modify many existing models of physical processes; the representation of equivalent models in integer-order derivatives provided a good approximation of the electrochemical response of the model.

On the basis of Cole’s proposal to add an extra degree of freedom in order to solve the RC circuits for characterization purposes and improve the correlation in the adjustment to experimental data, we have developed analytical arguments to derive this result based on integration in weighted individual relaxation processes. However, the distributions of relaxation times involve complex functions and are difficult to measure. This study has shown a pattern in which the Cole type behavior appears as a result of competition between a capacitive and resistive behavior within the sample, characterized by the fractional order derivative of the applied voltage. This combination of stored and dissipated energy is conveniently based on the representation of linear viscoelastic behavior; this dissipation is known as internal friction. In the literature it is common to characterize based on least-squares fit of equivalent electrical circuit models on experimental data, including Cole models. From the description of the fractional differential equation models it can be noted that the representation of Cole models is derived as a particular solution to the RC circuit under FC.

Since some authors replace the derivative of a fractional entire order in a purely mathematical context, the physical parameters involved in the differential equation do not have the dimensionality obtained in the laboratory and from the physical point of view of engineering this is not entirely correct. In this representation an auxiliary parameter is introduced; this parameter characterizes the existence of the fractional temporal components and related the time constant of the system; the solutions presented here preserve the dimensionality of the studied system for any value of the exponent of the fractional derivative. The advantage of this alternative representation in when comparison with the models presented in the literature is the physical compatibility of the solutions.

The Caputo representation has the disadvantage that their kernel had singularity; this kernel includes memory effects and therefore this definition cannot accurately describe the full effect of the memory. Due to this inconvenience, Caputo and Fabrizio in [15] present a new definition of fractional derivative without singular kernel, the Caputo-Fabrizio fractional derivative. The two definitions of fractional derivatives must apply conveniently depending on the nature of the system and the choice of the fractional derivative depends upon the problem studied and on the phenomenological behavior of the system.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to thank Mayra Martínez for the interesting discussions. J. F. Gómez-Aguilar acknowledges the support provided by CONACYT: catedras CONACYT para jovenes investigadores 2014.