Derivation of the implicit current componentsLet us consider a simple one-electron reaction at a stationary, planar, macroscopic electrode of a dissolved redox couple, free of any adsorption phenomena and coupled chemical reactions:$${\text{Red}}\left( {{\text{sol}}} \right) \rightleftarrows {\text{Ox}}\left( {{\text{sol}}} \right) + {\text{e}}$$
(1)
The solutions for the surface concentrations of redox species are well known17:$$c_{{{\text{Red}}}} \left( {0,t} \right) = c^{*} – I_{con} \left( t \right)$$
(2)
and$$c_{{{\text{Ox}}}} \left( {0,t} \right) = I_{con} \left( t \right)$$
(3)
here \(I_{con} \left( t \right) = \mathop \smallint \limits_{0}^{t} \frac{I\left( \tau \right)}{{FA\sqrt D }}\frac{d\tau }{{\sqrt {\pi \left( {t – \tau } \right)} }}\) is the convolution, or semi-integral, c* is the bulk concentration of Red species initially present in the solution, D is the common diffusion coefficient, and other symbols have their common meaning. Equations (2) and (3) follow from the solution of the diffusion equation only, thus reflecting solely the process of diffusion. In another words, convolution integral \(I_{con} \left( t \right)\) is mathematical representation the diffusion effect in the course of the voltammetric experiment. It is however important to recall that the diffusion is driven by the concentration gradients prompted by the electrode reaction.If the electrode reaction is kinetically controlled, then the Butler-Volmer equation holds at the electrode surface, relating the current with the surface concentrations \(c_{{\text{i}}} \left( {0,t} \right)\), electrode potential, and kinetic parameters:$$\frac{I\left( t \right)}{{FA}} = k_{s} \exp (\beta \phi )\left[ {c_{{{\text{Red}}}} \left( {0,t} \right) – c_{{{\text{Ox}}}} \left( {0,t} \right)\exp ( – \phi \left( t \right))} \right]$$
(4)
Specifically, ks is the standard rate constant (cm s-1), β is the anodic electron transfer coefficient, and \(\phi \left( t \right) = \frac{F}{RT}\left( {E\left( t \right) – E^{{ \circ^{\prime}}} } \right)\) is dimensionless potential defined versus the formal potential (\(E^{{ \circ^{\prime}}}\)) of the electrode reaction. Combining kinetic Eq. (4), which explains the phenomena at the electrode surface, with Eq. (2) and (3), which describe the mass transfer phenomena prompted by the electrode reaction, provide the theoretical basis for understanding electrochemical experiments under chronoamperometric and voltammetric conditions.Equation (4) is derived assuming, that implicit anodic Ia(t) and cathodic Ic(t) current components determine the net I(t), experimentally measurable, current:$$I\left( t \right) = I_{a} \left( t \right) + I_{c} \left( t \right)$$
(5)
where$$I_{a} \left( t \right) = FAk_{s} \exp (\beta \phi \left( t \right))\left[ {c_{{{\text{Red}}}} \left( {0,t} \right))} \right]$$
(6)
$$I_{c} \left( t \right) = – FAk_{s} {\text{exp}}\left[ { – \left( {1 – \beta } \right)\phi \left( t \right)} \right]c_{{{\text{Ox}}}} \left( {0,t} \right)$$
(7)
Note the negative sign in front of the cathodic current in Eq. (7), thus the net current I(t) in Eq. (5) is represented as a sum of the anodic and cathodic current components. Clearly, substituting Eq. (2) into (6), and (3) into (7), the anodic and cathodic current components can be related with the total, net current through the convolution integral \(I_{con} \left( t \right)\) as follows:$$I_{a} \left( t \right) = FAk_{s} \exp (\beta \phi \left( t \right))\left( {c^{*} – I_{con} \left( t \right)} \right)$$
(8)
$$I_{c} \left( t \right) = – FAk_{s} {\text{exp}}\left[ { – \left( {1 – \beta } \right)\phi \left( t \right)} \right]I_{con} \left( t \right)$$
(9)
For the sake of clarity let us emphasise that though convolutive integral is required to calculate implicit current components \(I_{a} \left( t \right)\) and \(I_{c} \left( t \right)\), they are just common, not convolutive currents, making clear difference of the present method with the convolutive voltammetry, where the convolutive current \(I_{con} \left( t \right)\) is analysed as a function of the potential13,14,15,16.Considering the last two Eqs. (8) and (9), it follows that the ratio of Ia(t)/Ic(t) results in a relation in which kinetic parameters (ks and β) cancel out:$$\frac{{I_{a} \left( t \right)}}{{I_{c\left( t \right)} }} = – \frac{{c^{*} – I_{con} \left( t \right)}}{{\exp \left( { – \phi \left( t \right)} \right)I_{con} \left( t \right)}}$$
(10)
Equation (10) is the central equation of the current study, revealing that the ratio of current components is clearly contingent upon the potential alone. Solving the system of Eqs. (5) and (10) straightforwardly yields the anodic and cathodic current components at any moment during the voltammetric experiment. This necessitates the initial calculation of the convolution integral and prior knowledge of the formal potential of the redox couple, without any prior knowledge of the electrode kinetic parameters (ks and β). Once being revealed, the implicit current components lead straightforwardly to the intrinsic kinetic parameters, as clearly shown by Eqs. (8) and (9). Such a theoretical consideration is generally valid for any voltammetric and chronoamperometic experiment.Let us now consider theoretically a simple experiment under conditions of linear sweep cyclic voltammetry. The response can be calculated in a dimensionless current domain, according to the numerical step-function method18:$$\Psi_{m} = \frac{{\kappa \exp \left( {\beta \phi_{m} } \right)\left[ {1 – \frac{{2\sqrt {dEF} }}{{\sqrt {\pi RT} }}\left( {1 + \exp \left( { – \phi_{m} } \right)} \right)\mathop \sum \nolimits_{j = 1}^{m – 1} \Psi_{j} S_{m – j + 1} } \right]}}{{1 + \frac{{2\kappa \sqrt {dEF} }}{{\sqrt {\pi RT} }}\left( {1 + \exp \left( { – \phi_{m} } \right)} \right)}}$$
(11)
where the dimensionless current is defined as \(\psi \left( t \right) = \frac{{I\left( t \right)\sqrt {RT} }}{{FAc^{*} \sqrt {DvF} }}\), \(\kappa = k_{s} \sqrt {\frac{RT}{{DvF}}}\) is the dimensionless electrode kinetic parameter, dE is the potential increment, and \(S_{m} = \sqrt m – \sqrt {m – 1}\) is the numerical integration parameter with a serial number m. Once the net current ψ is known (either being calculated with the aid of Eq. (11), or in the real experiment, measuring the current at each time t of the experiment), the system of Eqs. (5) and (10) can be used to obtain the cathodic and anodic implicit current components in a dimensionless form:$$\Psi_{c,m} = \frac{{\Psi_{m} b_{m} }}{{b_{m} – a_{m} }}$$
(12)
$$\Psi_{a,m} = \Psi_{m} – \Psi_{c,m}$$
(13)
where \(a_{m} = 1 – \frac{{2\sqrt {dEF} }}{{\sqrt {RT\pi } }}\Psi con_{m} ,\) \(b_{m} = \frac{{2\sqrt {dEF} \exp \left( { – \varphi_{m} } \right)}}{{\sqrt {RT\pi } }}\Psi con_{m}\), and \(\Psi con_{m} = \mathop \sum \limits_{j = 1}^{m} \Psi_{j} S_{m – j + 1}\) is the numerical form of the convolution integral.The initial outcome of such calculations is illustrated in Fig. 1, where the blue curves represent typical dimensionless conventional cyclic voltammograms for different degree of electrochemical reversibility, encompassing sluggish (electrochemically irreversible), quasireversible, and fast, apparently electrochemically reversible electrode reaction, calculated by means of Eq. (11). Once the net currents are known, and using only the formal potential values, one can calculate the implicit anodic and cathodic current components in the course of the voltammetric experiment by means of Eqs. (12) and (13). For each half-cycles of the voltammetric experiment corresponding implicit anodic and cathodic current components are calculated. Figure 1 illustrates a revised version of cyclic voltammograms, showcasing the sole anodic (red curves) and cathodic (green curves) current components associated with the forward (anodic) and reverse (cathodic) half-cycles of the virtual voltammetric experiment.Figure 1Typical dimensionless cyclic voltammograms for different degree of electrochemical reversibility (blue curves), including the sole anodic (Ψa; red curves) and cathodic (Ψc; green curves) implicit current components corresponding to the anodic and cathodic half-potential circles, respectively. The electrode kinetic parameter \(\kappa\) is 0.1 (A); 0.5 (B), and 5.5 (C). Panel (D) shows current components corresponding to the forward scan of the cyclic voltammetric experiment only, for electrode kinetic parameter values corresponding to panels (A)–(C). The electron transfer coefficient is β = 0.5, potential increment is dE = 0.1 mV, and temperature is T = 298.15 K.In instances where the electrode reaction is sluggish, characterized by the electrode kinetic parameter \(\kappa\) = 0.1 (log(κ) ≤ − 1), the sole current components align seamlessly with the net current, as evidenced in Fig. 1A. This alignment could serve as a real criterion for a precise defining of the electrochemical irreversibility in terms of the electrode kinetics. Under these conditions, the implicit anodic component ψa is identical with the net current in the course of the forward, anodic potential scan, implying that the contribution of the implicit cathodic current component in the course of the forward scan is insignificant. An analogous discussion can be given for the reverse, cathodic potential scan.The relationship between the current components and the net current undergoes significant changes for typically quasireversible and fast (apparently electrochemically reversible) electrode reactions, as depicted in Fig. 1B,C, respectively. Specifically, for a quasireversible scenario (− 1 ≤ log(κ) ≤ 1), the current components already surpass the net current, as illustrated in Fig. 1B (\(\kappa\) = 0.5; log(\(\kappa\)) = − 0.3). During the forward, anodic potential scan, the implicit anodic current component (Ψa) is larger than the net current. This implies that the implicit cathodic current component is also significant, resulting in a lower overall net current due to their difference. It is important to note that the cathodic current component (Ψc) corresponding to the forward potential scan is not shown on the plot for the sake of presentation simplicity. Similarly, for analogous reasons, the implicit cathodic current component corresponding to the reverse potential scan of the experiment is larger than the net cathodic current (green curve in Fig. 1B).The discrepancy between the net and current components becomes significant when the electrode reaction is very fast (log(κ) > 0.5), as depicted in Fig. 1C. For κ = 5.5, the peak current of the implicit anodic current component (Ψa) is almost an order of magnitude larger than the net anodic peak current of the conventional cyclic voltammogram (compare red and blue curves in Fig. 1C). This indicates that the electrode reaction Red ⇢ Ox + e proceeds rapidly, generating a highly intensive anodic current component. However, simultaneously, the reverse reaction Ox + e ⇢ Red occurs at a comparable rate, resulting in an intensive implicit cathodic current component (not shown). Thus, the overall effect is represented by a relatively low net current in the course of the forward potential scan. A similar discussion is relevant to the reverse, cathodic potential scan of the experiment.For the sake of simplification, the analysis of the implicit anodic and cathodic current components could be restricted only to the initial, forward scan of the cyclic voltammetric experiment, as shown in Fig. 1D. Such analysis reveals the degree of electrochemical reversibility in both qualitative and quantitative terms.Analysis of the implicit current components in the apparently reversible kinetic regionTo unlock the genuine potential of the proposed methodology, one must scrutinize electrode reactions that are very fast and apparently electrochemically reversible. The conventional cyclic voltammograms depicted in Fig. 2A illustrate this progression, ranging from a very fast quasireversible electrode process (κ = 2.53) to apparently reversible electrode processes (κ > 5). Within this kinetic range, the conventional cyclic voltammograms appear to be identical, seem to become independent of the electrode kinetics and predominantly affected by the mass transfer only.Figure 2(A) Conventional cyclic voltammograms of very fast quasireversible electrode reaction characterized with the electrode kinetic parameter \(\kappa\) of 2.53; 5.1, and 50.1 (curves are overlapped). Panels (B–D) showcase differential current component (Ψdif) associated with the anodic half-cycle of the cyclic voltammogram simulated for an increasing degree of the electrochemical reversibility. The electrode kinetic parameter \(\kappa\) is 2.53 (B), 5.1 (C), and 50.1 (D). The other conditions are the same as for Fig. 1.For an in-depth analysis of the behavior of implicit current components in the apparently reversible kinetic region, we propose an approach which resembles the reasoning of differential voltammetric techniques19,20. For this purpose, we use the implicit current components associated with the forward, anodic potential scan of the cyclic voltammetry only, and by analogy with differential voltammetric techniques, we define a differential current as \(I_{dif} \left( t \right) = I_{a} \left( t \right) – I_{c} \left( t \right)\). Recalling the definitions of the current components through Eqs. (8) and (9), the expression for the differential current is given by:$$I_{dif} \left( t \right) = FAk_{s} e^{\beta \phi \left( t \right)} \left[ {c^{*} – I_{con} \left( t \right)\left( {1 – {\text{e}}^{ – \phi \left( t \right)} } \right)} \right]$$
(14)
The reason for introducing the differential current is that it has unique features, which are not attributed to the conventional, net current. Specifically, at the formal potential (\(E = E^{{ \circ^{\prime}}} ; \phi = { }0\)), the differential current equals \(I_{dif} = FAc^{*} k_{s}\). Obviously, under these conditions the convolution integral is cancelled out in Eq. (14), removing the effect of the mass transport in this artificial current component. In addition, the differential current is becoming independent on the electron transfer coefficient. It follows that \(I_{dif} \left( {E = E^{{ \circ^{\prime}}} } \right)\) represents the exchange current for equilibrium conditions of the studied redox couple. Expressing this in a dimensionless form is \(\Psi_{dif} \left( {E = E^{{ \circ^{\prime}}} } \right) = \kappa\).For better clarity of these profound features of the differential current, which could be calculated only by previous knowledge of the implicit anodic and cathodic current components, let us write an equation for the common, net current, analogous to Eq. (14):$$I\left( t \right) = FAk_{s} e^{\beta \phi \left( t \right)} \left[ {c^{*} – I_{con} \left( t \right)\left( {1 + {\text{e}}^{ – \phi \left( t \right)} } \right)} \right]$$
(15)
The difference between Eqs. (14) and (15) lies solely in the sign within the small bracket. It is easy to see that at \(E = E^{{ \circ^{\prime}}} ( \phi = { }0\)) the convolution integral (i.e., the effect of diffusion) cannot be cancelled out in the common net current I(t) (Eq. (15)). Thus, it remains affected by both diffusion and electrode kinetics, and their intricate complex interplay is the sole origin of complexity in the interpretation of the voltammetric data in general.The analysis of an electrode reaction by means of differential current is illustrated in Fig. 2B–D, depicting the evolution of the differential current during the forward scan of the cyclic voltammetry. Notably, when the electrode reaction is rapid and apparently reversible, the differential current component calculated aligns at the formal potential (Fig. 2C,D). Its peak current provides a direct indication of the electrode kinetic parameter \(\kappa\) (note the difference in the scale of the ordinate in Fig. 2B–D). Therefore, the data shown in Fig. 2 clearly suggest that the proposed methodology enables the estimation of rapid electrode kinetics using standard macroscopic electrodes, which seems to be a unique virtue of the present methodology. On the other hand, the peak potential of the differential current depends on the electron transfer coefficient, enabling straightforward estimation of this parameter. For instance, for a very fast reaction associated with \(\kappa\) = 50, the dependence of the peak potential of the differential current and the anodic electron transfer coefficient β is linear, associated with the following linear regression line: Ep = 0.115 β -0.056 (R 2 = 0.999), over the range of β values from 0.2 to 0.8. Therefore, the analysis of the differential current component provides a means for simple, fast, and unambiguous estimation of both standard rate constant and the electron transfer coefficient.Limitations of the methodThe distinctive feature of the proposed method for elucidating the kinetic parameters of exceedingly rapid, ostensibly reversible electrode reactions must be scrutinized critically in relation to a genuinely reversible electrode reaction. Let us revisit the remarkably fast electrode reaction associated with the electrode kinetic parameter κ = 50. This high value of the electrode kinetic parameter corresponds to an electrode reaction with, for instance, a standard rate constant of ks = 0.1 cm s−1, diffusion coefficient of D = 5 × 10−6 cm2 s−1 under experiment scan rate of v = 20 mV s−1. Consequently, the electrode reaction falls comfortably within the apparent reversible kinetic region, where the peak currents of the conventional cyclic voltammogram are virtually independent of the electrode kinetics. In addition, let us now presume that the electron transfer behaves as a truly thermodynamically reversible process, described by the Nernst equation. The cyclic voltammograms corresponding to these two cases are compared in Fig. 3A and appear to be identical over the entire potential window. Applying Eqs. (12) and (13) to the voltammogram associated with \(\kappa\) = 50 yields well-defined anodic and cathodic current components, while the peak current of the differential current component is equal to \(\kappa\), as discussed previously (cf. Fig. 2D). However, attempting to apply the same formulas (12) and (13) to the voltammogram of a truly reversible electrode reaction obeying the Nernst equation at any potential of the voltammetric experiment (e.g., ferrocene/ferrocenium redox couple) does not yield any reasonable data. Consequently, the net current of a truly reversible electrode reaction cannot be broken down into separate anodic and cathodic current components.Figure 3(A) Conventional cyclic voltammograms of very fast, apparently reversible electrode reaction characterized with the electrode kinetic parameter \(\kappa\) = 50.1 (red circles) and true reversible electrode reaction (blue solid line). (B) Relative variance between the currents of the voltammograms from panel (A) calculated for the first anodic half-cycle, calculated as \(\frac{{\left| {{\varPsi}\left( {\kappa = 50} \right) -{\varPsi}\left( {{rev}} \right)} \right|}}{{\left| {{\varPsi}\left( {\kappa = 50} \right)} \right|}}100{\text{\% }}\).In a meticulous consideration of these theoretical data, a critical question arises: How identical are the net current data of a truly reversible electrode reaction and an apparently reversible electrode reaction, as implied by Fig. 3A? The relative variance in the currents of these two voltammograms for the first half of the cyclic voltammetry is plotted in Fig. 3B. It is evident that they slightly differ (up to 1%) at the very beginning of the rising part of the voltammogram. Specifically, at low potentials, when the driving force of the electrode reaction is low enough, the electrode kinetics of the quasireversible electrode reaction still plays a role (red circles in Fig. 3A), being a slower step relative to diffusional mass transfer. Such an effect is completely absent in a truly reversible electrode reaction (black line in Fig. 3A). Therefore, these small differences between the very fast, apparently reversible, and the truly reversible reaction are of critical importance.Considering that even tiny differences in the net currents play a decisive role in the successful application of the method, it is reasonable to anticipate that the noise from the current, an inevitable artefact in a real experiment, will hold significant importance. In order to investigate this effect through simulations, uniformly distributed noise was introduced to the recurrent formula (11) for calculating the net current, and the results are presented in Fig. 4. The noisy current was simulated by the extending the recurrent formula (11) as follows:$$\Psi \left( {noise} \right)_{m} = \frac{{\kappa \exp \left( {\beta \phi_{m} } \right)\left[ {1 – \frac{{2\sqrt {dEF} }}{{\sqrt {\pi RT} }}\left( {1 + \exp \left( { – \phi_{m} } \right)} \right)\mathop \sum \nolimits_{j = 1}^{m – 1} \Psi(noise)_{j} S_{m – j + 1} } \right]}}{{1 + \frac{{2\kappa \sqrt {dEF} }}{{\sqrt {\pi RT} }}\left( {1 + \exp \left( { – \phi_{m} } \right)} \right)}} + Noise_{m}$$
(22)
where Noisem is a function generating random numbers with a uniform distribution within the interval defined by the noise percentage, relative to the current measured at the lower foot of the idealized voltammograms (i.e., percentage of the current calculated by the recurrent formula (11) at potentials at least 200 mV less than the formal potential). The added noise at the beginning of the voltammogram propagates throughout the simulation due to the recurrent nature of Eq. (22). The analysis presented in Fig. 4 corresponds to an experimental scenario examining a very fast electrode reaction (κ = 50) under varying noise levels ranging from 1 to 5%. Panels A–C illustrate the anodic and cathodic current components derived from the raw, noisy data, referring to the forward scan of the cyclic voltammetric experiment. Evidently, processing the raw data accentuates the noise in both the anodic and cathodic current components, as a consequence of the convolution (summation) procedure. However, the overall morphology, peak currents, and peak potentials of both implicit current components remain virtually unaffected. To substantiate this claim, Fig. 4D presents the implicit current components estimated from the idealized, noise-free net currents. The peak currents and potentials, crucial for estimating the kinetic parameters, are virtually identical between the noisy and idealized data, demonstrating the robustness of the method.Figure 4Anodic (Ψa, red curves) and cathodic (Ψc, green curves) implicit current components calculated from a noisy net current of a very fast, apparently reversible electrode reaction with the electrode kinetic parameter of κ = 50. For the sake of simplicity only the first half-cycle of the conventional cyclic voltammograms was considered. The level of noise, with a uniform distribution, was 1% (A), 3% (B) and 5% (C) relative to the net current measured at the lower foot of the cyclic voltammogram. Panel (D) displays the anodic and cathodic current components calculated from idealized, noise-free net currents. Panels (E) compares implicit current components estimated from net current with 1% noise, after smoothing the row data and smoothing the implicit current components (doted lines), with data calculated from idealized, noise-free net currents (full lines). The blue lines correspond to differential (ψdif) component. The smoothing of the noise net current was performed by applying “medsmooth” function of Mathcad23, with a weight of 5 in averaging the raw data.Yet, the noise effect should be mitigated by smoothing the net currents before applying Eqs. (12) and (13), as well as by filtering and smoothing the obtained data for both anodic and cathodic implicit current components, while a plethora of tools are available for such a purpose21,22. An example of such procedure is given in Fig. 4E, including the initial smoothing of the raw data with the noise of 1%, and filtering and smoothing of the obtained anodic and cathodic implicit current components by means of the function “medsmooth” of the Mathcad software23. Let us note that the “medsmooth” function performs median filtering. It is the most robust among smoothing functions because it is less susceptible to influence from outliers. However, it tends to smooth out sharp features in the data. One can contrast this approach with Gaussian kernel smoothing or localized least-squares smoothing methods. In addition, Fig. 4E compares the processed data (dotted lines) with the idealized, noise free data, illustrating the effectiveness of the smoothing procedure. Importantly, Fig. 4E also displays the differential current component, suggesting that in the presence of noise, the proposed method leads to reliable estimation of the standard rate constant (from the peak current of the differential component) and the electron transfer coefficient (from the peak potential of the differential component).Experimental verificationThe characterization of electrode kinetics, besides employing scanning electrochemical microscopy24, commonly utilizes cyclic voltammetry and electrochemical impedance spectroscopy25. Various methodologies are employed when applying linear sweep potential in studies, as detailed by Nicholson26, Klingler-Kochi27, or Gileadi28. Additionally, the development of (ultra)microelectrodes was prompted for the purpose of evaluating heterogeneous electron transfer kinetics of very fast processes by changing the regime of diffusion to be faster relative to the electron transfer29.Electrode reaction of hexaammineruthenium(III) complex is one of the fastest electrode reactions of a dissolved redox couple24, being a perfect example of apparently electrochemically reversible process under common (slow scan rates) voltammetric conditions. Typical cyclic voltammogram for the reduction of ruthenium system is shown in Fig. 5A, exhibiting all characteristics of an electrochemically reversible electrode reaction of a dissolved redox couple. Over the scan rate interval from 20 to 500 mV/s, the peak potential difference remains constant at 58.5 ± 3 mV, with a mid-potential of − 0.129 ± 0.003 V, and the peak current ratio of 1.025 ± 0.033. Both cathodic and anodic peak currents depend linearly on the square-root of the scan rate with a linear regression coefficient of R2 = 0.999 and 0.998, respectively. Invoking the Randles–Sevcik Eq. 30, the slope of these lines enabled to estimate the real electrode surface area to be 7.1 × 10–6 m2, assuming a common diffusion coefficient of 9 × 10–10 m2 s-131,32,33.Figure 5Typical cyclic (A) and square-wave voltammogram (B) of 1 mmol/L [Ru(NH3)6]Cl3 at GCE in 1 mol/L KNO3 supporting electrolyte at scan rate 20 mV/s and potential increment of 2.44 mV for (A) and frequency of 25 Hz, SW amplitude of 25 mV and step potential 5 mV, showing forward (cathodic; red), backward (anodic; blue) and net SW component (black). All curves have been smoothed and in the case of CV background corrected. The left inset of (B) shows the dependence of the dimensionless net peak current in SWV on the of the inverse square-root of the frequency. Red circles correspond to the normalized experimental data, while the black line represents the theoretical data simulated for the standard rate constant of 0.05 cm s-1 and electron transfer coefficient β = 0.5. The other conditions for normalisation of the experimental currents are: diffusion coefficient D = 9 × 10–6 cm2 s−1, stoichiometric number of electrons n = 1, and temperature 293.15 K. The right inset of (B) shows the dependency of the amplitude-normalized net peak currents on the SW amplitude. Red circles show experimental data while the black line corresponds to the simulations performed for identical kinetic parameters as in the left inset. The experiments and simulations are conducted at the frequency of 100 Hz.An attempt to shift the electrode reaction to the quasireversible kinetic region was achieved by applying square-wave voltammetry (SWV)17. Initially, the electrode reaction was studied at moderate frequency of 25 Hz and a low amplitude of 25 mV, when the electrode kinetic is expected to be within the electrochemically reversible kinetic region (Fig. 5B). The peak potential of the net component is -0.128 V, which is virtually identical with the mid-potential measured with CV, suggesting the reversible character of the electrode reaction. The other voltammetric features of the response shown in Fig. 5B, such as peak potential difference of the cathodic and anodic components, as well as the half-peak width of the net peak align well with the theoretical expectation for a very fast, apparently reversible electrode reaction17. Increasing the frequency up to 400 Hz does not affect the net peak potential, which remains at the value of -0.128 V; yet, the peak current does not increase strictly linearly with the square-root of the frequency, implying that the effect of the quasireversible kinetic behaviour commences to be effective, in particular, at frequencies larger than 100 Hz. Let us note that the frequency of 100 Hz corresponds to a typical time window of the voltammetric experiment of 5 ms in SWV.17 Assuming the same time interval in CV over potential window of, e.g., 0.4 V, would corresponds to the scan rate of 160 V/s.Thus, at the SW frequencies f > 100 Hz, the electrode reaction is gradually shifted from the reversible to the quasireversible kinetic region. The theory of SWV suggest that the dimensionless net peak current is becoming sensitive to the frequency, whereas the net peak potential remains unaffected up to the irreversible kinetic region, which is in excellent accord with the experimental data. The left inset of Fig. 5B shows the dependency of the normalized (dimensionless net peak current) on the parameter \(\log \left( {\frac{1}{\sqrt f }} \right)\). Note that the latter parameter \(\frac{1}{\sqrt f }\) corresponds to the dimensionless electrode kinetic parameter in SWV \(\left( {\kappa = \frac{{k_{s} }}{{\sqrt {Df} }}} \right)\). The red circles in the left inset of Fig. 5B represent the experimental data, which aligns with the theoretically predicted dependence of ψp vs. log(\(\frac{1}{\sqrt f }\)) (black line)17. The linearly increasing part of the function ψp vs. log(\(\frac{1}{\sqrt f }\)) is typical for the quasireversible kinetic region. An attempt to fit the experiment with the simulations, results in the best fit with a standard rate constant of ks = (0.05 ± 0.01) cm s-1. In addition, kinetic measurements in SWV can be done by means of the amplitude variation at the constant frequency as elaborated in34. Applying the method of the amplitude-based quasireversible maximum35, the fitting of the experimental with the theoretical data is presented in the right inset of Fig. 5B, achieved again with the standard rate constant being within the interval of 0.05 ± 0.02 cm s−1.This general kinetic characterisation of the ruthenium redox system confirms the well-known fact that it is one of the fastest electrode reactions of a dissolved redox couple studied36. Under moderate experimental conditions, as those in Fig. 5A, the electrode reaction is virtually reversible process. The electrode kinetics could be estimated in CV by applying conventional analysis by varying the scan rate over 160 V/s; the artefacts arising from such extreme scan rates are quite predictable, which in reality practically preclude kinetic measurements at a macroscopic electrode. Yet the methodology proposed in this study promises kinetic estimation at mild experimental conditions, i.e., macroscopic electrode and low scan rate, i.e., in conditions when the artefacts arising from uncompensated resistance and charging current can be reasonably neglected. Specifically, for the present experimental conditions and the scan rate of 0.02 V/s, the charging current is estimated to be within the interval from 0.02 to 0.08 µA, whereas the uncompensated resistance in 1 mol/L KNO3 used as a supporting electrolyte and the distance between the working and reference electrode of 1 cm is estimated to be in the range of 9 mΩ, which is insignificant for distortion of the applied potential sweep.Let us first study the evolution of the classical convolution current as a function of the potential, corresponding to the initial, cathodic potential scan of the cyclic voltammetric experiment of ruthenium system, which is shown in Fig. 6A. The sigmoid curve is constructed from experimental currents after normalisation, smoothing and background correction for the experiment conducted at low scan rate of 20 mV/s. The convolution current is virtually free of noise, with a half-wave potential of − 0.129 V, which is in excellent agreement with the formal potential estimated from the conventional CV and SWV. Having in hand the convolution currents and applying Eqs. (12) and (13), the genuine anodic and cathodic implicit current components corresponding to the forward scan of the cyclic voltammetry, have been calculated and depicted in Fig. 6B. The shape of the curves is in excellent agreement with the theoretical predictions for a very fast electrode reaction (compare Fig. 6B with 1D). The intensity of both cathodic and anodic current components is an order of magnitude higher than the conventional voltammogram, shown with a blue line in Fig. 6B. The peak potential of both cathodic and anodic implicit current components is identical, located at the formal potential of the electrode reaction, while the peak current ratio is virtually one. Such evolution of current components is typical for very fast, apparently reversible electrode processes, with an electron transfer coefficient equal to 0.5, as illustrated in Fig. 1D. Following the proposed methodology, a differential current component is constructed as a \(\Psi_{dif} = \Psi_{c} – \Psi_{a}\), which is depicted as an inset of Fig. 6B. The peak current of the differential component is 7.42. Thus, the dimensionless electrode kinetic parameter for the ruthenium system is \(\kappa\) = 7.42 leading to the standard rate constant \(k_{s} = \kappa \sqrt {\frac{DvF}{{RT}}}\) = 0.02 cm s-1. To the best of our knowledge, this is the only method enabling to estimate such a high standard rate constant at a macroscopic electrode and a single scan rate of 20 mV/s.Figure 6(A) Convolution dimensionless current as a function of the electrode potential and (B) dimensionless cathodic (green circles) and anodic (red circles) implicit current components as function of the relative potential E vs. E0′ of the ruthenium system, corresponding to the forward half-cycle of the cyclic voltammogram shown in Fig. 6B (solid blue line). The inset of panel (B) shows the differential current component calculated as a difference between the cathodic and anodic implicit current components. All experimental conditions are identical as in Fig. 5A.To simplify, the previous estimation assumed identical diffusion coefficients for both redox species in the ruthenium system. In a more rigorous approach, different diffusion coefficients should be employed, as discussed in our previous study.12 In such cases, the implicit anodic and cathodic current components additionally depend on the ratio ρ = DR/DO, where DR and DO are the diffusion coefficients of the reduced and oxidised form, respectively. The work of Compton et al.31 reports values of 5.3 × 10–6 cm2/s and 7.3 × 10–6 cm2/s for Ru(III) and Ru(II), respectively in 0.5 mol/L KNO3 as a supporting electrolyte. Assuming these values applied under the current experimental conditions, the diffusion coefficient ratio is ρ = 1.38, which insignificantly affects the magnitude of the convolution currents Ψcon calculated for each potential in the cyclic voltammograms of the ruthenium system.More critically, when diffusion coefficients are unequal, the mid-potential Em estimated from the cyclic voltammograms is \(E_{m} = E^{0^{\prime}} + \frac{RT}{F}ln\sqrt \rho\), differing from the formal potential. Conversely, precise knowledge of the formal potential is crucial for estimating implicit current components (cf. Eqs. (12) and (13) and the definition of parameters a and b therein) through the values of the dimensionless potential \(\phi\). Adjusting for this, the formal potential of the ruthenium system is estimated at -0.133 V, which is 4 mV more negative than the previously used value of -0.129 V. Revising the standard rate constant by incorporating the ρ parameter12 and using the corrected formal potential yields again values in the range of 0.02 ± 0.005 cm s−1.In Table 1 we are giving a brief summary of literature values for standard rate constants of the ruthenium system. As often happens, the kinetic data differ significantly depending on the method used. Here, our data are supported by an independent estimation using SWV.Table 1 A selection of standard rate constants estimated for ruthenium redox system.
Genuine anodic and cathodic current components in cyclic voltammetry
