Next Article in Journal
A Special Study of the Mixed Weighted Fractional Brownian Motion
Next Article in Special Issue
A Data-Driven Memory-Dependent Modeling Framework for Anomalous Rheology: Application to Urinary Bladder Tissue
Previous Article in Journal
Population Forecast of China’s Rural Community Based on CFANGBM and Improved Aquila Optimizer Algorithm
Previous Article in Special Issue
Incorporating Rheological Nonlinearity into Fractional Calculus Descriptions of Fractal Matter and Multi-Scale Complex Fluids
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Time-Fractional Phase Field Model of Electrochemical Impedance

by
Pavel E. L’vov
1,*,
Renat T. Sibatov
2,3,*,
Igor O. Yavtushenko
1 and
Evgeny P. Kitsyuk
3
1
Laboratory of Diffusion Processes, Ulyanovsk State University, 432017 Ulyanovsk, Russia
2
Department of Theoretical Physics, Moscow Institute of Physics and Technology (MIPT), 141700 Dolgoprudny, Russia
3
Scientific–Manufacturing Complex “Technological Centre”, 124498 Moscow, Russia
*
Authors to whom correspondence should be addressed.
Fractal Fract. 2021, 5(4), 191; https://doi.org/10.3390/fractalfract5040191
Submission received: 13 September 2021 / Revised: 25 October 2021 / Accepted: 26 October 2021 / Published: 30 October 2021
(This article belongs to the Special Issue Fractal Media and Fractional Viscoelasticity)

Abstract

:
In this paper, electrochemical impedance responses of subdiffusive phase transition materials are calculated and analyzed for one-dimensional cell with reflecting and absorbing boundary conditions. The description is based on the generalization of the diffusive Warburg impedance within the fractional phase field approach utilizing the time-fractional Cahn–Hilliard equation. The driving force in the model is the chemical potential of ions, that is described in terms of the phase field allowing us to avoid additional calculation of the activity coefficient. The derived impedance spectra are applied to describe the response of supercapacitors with polyaniline/carbon nanotube electrodes.

1. Introduction

Electrochemical impedance spectroscopy (EIS) is one of the most informative methods for studying charge carrier transport characteristics in different types of materials [1,2]. The basic approach of the method implies the measurement of impedance dependence on applied voltage frequency with low value of amplitude. EIS is widely applied to study electrochemical systems and solid-state devices [2]. In case of normal diffusion of charge carriers in semi-infinite planar material, the Warburg impedance is usually observed [3,4]: Z ( ω ) ( i ω ) 1 / 2 .
In general case, the model of Warburg impedance cannot be applied to materials that exceed the basic assumptions of the approach. Presence of lattice defects in material, nonlinear effects in diffusion process, various boundary conditions, etc. can cause the significant deviation of Z ( ω ) for real systems from Warburg impedance. Therefore, some extensions of the Wargurg model have been developed (see, e.g., in [3,4]).
Charge transport in electrolytic systems is described by the system of electrochemical transport equations. The chemical and electric potential gradients play the role of driving forces in the system. If there are two types of ions only, the coupled equations for the concentrations of ions can be simplified by the ambipolar diffusion approximation, and the problem can be reduced to ambipolar diffusion equation with nonlinear diffusion coefficient that is expressed via the activity coefficients of ions. Application of the approach demands additional calculation of the activity coefficients.
Most microscopic models of electrochemical power sources are based on a system of diffusion equations that do not take into account collective transport and the anomalous type of ion diffusion, and are unable to explain the nonlinear nature of the responses of the devices under consideration. The mathematical basis of anomalous self-similar diffusion is formed by equations with fractional-order derivatives [5,6]. The phase field method is an efficient and informative method to model and describe transport phenomena and phase evolution of the microstructure in electrochemical power sources [7,8,9].
Application of the phase-field theory to the analysis of intercalation and transport of ions in electrochemical systems was launched by Guyer et al. [7,8] who developed the model of electrodeposition of ions at the sharp electrode/electrolyte interface. The dynamics of intercalant concentration was described in terms of the Cahn–Hilliard equation and regular solution approximation for interaction of the solute and solvent.
Fleck et al. [10] developed the model of diffusion limited solid–solid phase transformation during lithium intercalation in LiFePO 4 nanoparticles. The model describes propagation of the lithium concentration front in a single particle by the Cahn–Hilliard equation with anisotropic mobility, whereas transformation of the particle shape is described in terms of the Allen–Cahn equation taking into account elastic energy. In [11], real-time tracking of lithium transport and conversion in FeF 2 nanoparticles by nanoscale imaging reveals a surprisingly fast conversion process in individual particles, with a morphological evolution resembling spinodal decomposition. One of the most popular methods to describe the spinodal decomposition is based on the Cahn–Hilliard equation [12,13].
In present work, the frequency dependences of the generalized Warburg impedance are derived for the case of ion transport controlled by the subdiffusive Cahn–Hilliard equation with fractional derivatives. The paper is organized as follows. Next section provides brief survey on fractional Cahn–Hilliard equations. Then, we discuss the Cahn–Hilliard equation for ionic transport and justify its fractional generalization within bipolar phase-field diffusion approximation. In Section 4, monofractional phase-field subdiffusion is discussed. In Section 5, we apply a linearized version of the fractional Cahn–Hilliard equation to derive the response to a constant input signal with small harmonic perturbation component and derive expressions for impedance under different conditions. Section 6 presents the application of the model to the description of impedance spectra of pseudo-capacitors with polyaniline/carbon nanotube electrodes.

2. Brief Survey on Fractional Cahn–Hilliard Equations

The classical models of the phase field are formulated in terms of partial differential equations of integer order, implying the locality of spatial interactions and the absence of memory effects. However, in the original Cahn–Hilliard formulation [12], non-local interactions were included in the physical model of the phase field, and later the model became local due to the assumption of slow spatial changes. There are different ways to include memory effects into phase-field model. One of the approaches is presented in [14], where the authors use the Jäckle and Frisch model of phase transition with memory and the Cahn–Hilliard model. As a result, they derive an analogue of the telegraph equation in the phase-field theory.
Other approaches utilize the fractional calculus [6,15,16]. A method based on kinetic equations with fractional derivatives [17,18,19,20,21] is often used to describe nonlocal transport. There are several works, where fractional phase field models are developed [22,23,24] to include nonlocal interactions and/or memory effects.
In a number of works, equations with derivatives of fractional order are derived from assumptions about the microscopic mechanism of the physical process. Let us list the main methods for deriving the fractional differential subdiffusion equation. Equations with fractional time derivative were obtained by Nigmatullin from the consideration of diffusion on fractal structures simulating porous and disordered media [25,26]. Another method is based on an asymptotic transition from the integral equations of random walks of the Scher and Montroll model [27] with a power-law distribution of waiting times to fractional-differential ones (for more details, see in [28]). The Scher and Montroll model is known to describe well the main features of dispersive transport. Based on the Langevin equation, Metzler and Klafter [5] developed a three-level description of subdiffusion with a transition to macroscopic fractional differential dynamics. In [5], a non-Markov generalization of the Chapman–Kolmogorov equation for random processes with continuous time, controlled by the distribution of waiting times, was investigated. The asymptotically power-law distribution (with “heavy tails”) of waiting times leads to the fractional Klein–Kramers equation, from which fractional generalizations of the Rayleigh and Fokker–Planck equations are derived. Fractional equations can be obtained from the balance equations of capture-emission of the multiple capture model with an exponential density of localized states (for more details, see in [29,30]). In [31], anomalous transport of charge carriers in two-dimensional systems with coexisting energetic and structural disorder is described by fractional diffusion model of multiple trapping in a comb-like geometry. All of the mentioned methods are based on assumptions about the transfer mechanism. These assumptions have to be substantiated for an individual case, for a specific material and for specific experimental conditions. In recent years, phase field models based on equations with fractional derivatives have been actively developed. When considering long-range interactions between particles, phase field models are introduced based on the Allen–Cahn equations or the Cahn–Hilliard equation with fractional spatial operators [22,32,33]. However, it should be noted that a direct sequential derivation of these fractional equations is not provided in the mentioned works. A serious problem is also associated with physically substantiated adequate boundary conditions for a real nonlocal problem. The form of these conditions fundamentally influences the character of the generalized diffusion Warburg impedance studied in this work.
In highly inhomogeneous porous media with a percolation structure, particles can be localized in restricted regions, and the residence times in these regions are often characterized by wide distributions with power-law ‘tails’. As is known, such type of distribution can lead to anomalous diffusion [27] described by an equation with a time fractional derivative [17]. A phase field model with a fractional time derivative of variable order has been successfully used to describe the evolution of structural damage and fatigue [34]. In [35], the scale behavior of the time-fractional Kardar–Parisi–Zhang equation in the interface growth model was investigated to study the kinetic roughening of surfaces. In [23], generalizations of the Allen–Cahn and Cahn–Hillard equations with a fractional time derivative are investigated to describe phase field propagation with account for the subdiffusion behavior of transport in heterogeneous porous materials.
In this work, we study the charge transport in the electrode–electrolyte system, where the conduction is provided by monopolar or bipolar transport of charges. The driving force is the chemical potential of ions that is described in terms of phase-field that allows us to avoid the additional calculation of the activity coefficient. In this approach, we should obtain the bipolar modification of well known Cahn–Hilliard equation that is widely applied in simulation of first order phase transitions [13,36]. Furthermore, Warburg impedance theory should be modified by the substitution of normal diffusion equation by the fractional Cahn–Hilliard equation for subdiffusive ionic systems in bipolar phase-field diffusion approximation.

3. Cahn–Hilliard Equation for Ionic Transport and Its Fractional Generalization

First, we consider a model of bipolar ion transport based on the classical Cahn–Hillard equation, and then introduce a fractional-differential generalization.
Let us consider the diffusion of charge carriers that have positive and negative signs. The evolvement of the carriers’ composition field can be described in terms of continuity equations (see, e.g., in [36])
c + t = j + j + e
c t = j j e
where c + c + ( r , t ) and c c ( r , t ) are the concentration fields of cations and anions, respectively; j + and j are the flux densities of the ions related to the process of diffusion; j + e = q + M + c + E and j e = q M c E are the flux densities defined by the applied electric field with strength E = φ coupled with electric potential φ , q + and q absolute values of ion charges; and M + and M are their mobilities. Furthermore, we suggest that ions do not interact with each other, so we exclude process of their recombination and possible diffusion of neutral particles.
We assume that diffusion of cations is described in terms of by the chemical diffusion j + = M + c + μ + [4], where μ + is the chemical potential of cations that is described on the basis of phase-field theory [37]:
μ + = g c + κ 2 c + .
Here, g ( c + ) is mixing energy depending on concentration of positive ions and κ is the gradient energy coefficient. In this approach, the chemical potential of positive ions depends on the second derivative of their concentration. The term introduces the principal difference from classical electrolyte theory based on the normal diffusion equation.
If the negative carriers are electrons we can define their motion as the process of normal diffusion that allows us to express the corresponding current density as j = D c , where D = M k B T , M is the mobility of electrons.
Thus, combining the previous relations we obtain the equations describing the diffusion of ions and electrons in electrolyte:
c + t = ( M + c + μ + q + M + c + E )
c t = ( M k B T c + q M c E )
Let us consider the charge carriers transport in the system in approximation of ambipolar diffusion. Following the approximation, we suggest the equality of absolute values of carrier charges q + = q = e and apply the condition of the local electric neutrality c = c + = c . A combination of equations with respect to the above suggestions leads us to equation describing charge transport in the system:
c t = ( M a c μ a )
where M a = M + M M + + M is the ambipolar mobility, and μ a is the effective (ambipolar) chemical potential that can be expressed as
μ a = g a c k 2 c .
Here, we introduce the effective mixing energy g a = g ( c ) + k B T c ln c . Ordinarily, ions have the mobility much less then mobility of electrons ( M + M ) , therefore the ambipolar mobility is approximately equal to mobility of ions M a M + . Usually, the mixing energy g ( c ) in binary alloys is described in terms of regular solutions approximation [12] that allows writing the function as
g a c = Ω c 1 c + k B T 2 c ln c + 1 c ln 1 c ,
where Ω is interaction parameter that is ordinarily related to critical temperature Ω = 2 k B T c . Equation (2) is similar to the well-known Cahn–Hilliard (CH) equation that is widely applied in simulation of first order phase transitions. In contrast to ordinarily application, the CH equation is employed for description of ions’ concentration field with respect to electronic subsystem taking into account ambipolar diffusion approximation. Furthermore, note that gradient energy coefficient is related with interaction parameter by the expression κ = Ω r 0 2 / 3 , where r 0 is the intermolecular distance [13].
The form of effective energy (3) can be interpreted as mixing energy with additional entropic term that can cause the effective shifting of the corresponding critical temperature, binodal and spinodal lines described by the function g a ( c ) . Figure 1 shows the effective transformation of the phase diagram with respect to influence of electronic subsystem. The critical temperature of the regular solution with respect to ionization and influence of electronic subsystem equals to T C a 0.69 T C that corresponds to composition c = 59.6 at %.
Thus, we conclude that Cahn–Hilliard equation with variable mobility (3) can be applied for systems consisting of cations and anions in ambipolar diffusion approximation.
Different specific mechanisms can lead to anomalous diffusion processes [5,6]. The features of anomalous diffusion are often described within the continuous time random walk (CTRW) model. The decoupled version of CTRW model implies that different jump lengths and waiting (trapping) times between two successive jumps are independent random variables. Fractional-order diffusion equation can be obtained as the asymptotic equation of CTRW with power law distributions of random trapping times τ : P { τ > t } t ν ( 0 < ν < 1 ). Considering CTRW with space dependent jump probabilities, Barkai et al. [38] derived a time-fractional Fokker–Planck equation (FFPE) for the case, when the mean waiting time diverges. FFPE describes anomalous diffusion in an external force field, it is derived in the following form [39,40]:
c ( x , t ) t = 0 D t 1 ν K x c ( x , t ) x + c ( x , t ) k T V x ,
where c ( x , t ) is particle concentration, V ( x ) is an external potential, K is a generalized diffusion coefficient and k T is the thermal energy.
Here,
0 D t 1 ν c ( x , t ) = 1 Γ ( ν ) t 0 t c ( x , τ ) ( t τ ) 1 ν d τ , 0 < ν 1 ,
is the Riemann–Liouville derivative of fractional order 1 ν [15]. Fundamental solutions to Equation (5) can be found in [39,40].
The fractional diffusion equation for the simplest case V ( x ) = 0 has the following form:
c ( x , t ) t = K 0 D t 1 ν 2 c ( x , t ) x 2 .
The Fourier transform of Equation (6) leads to the expression
( i ω ) ν c ˜ ( x , ω ) = K Δ c ˜ ( x , ω ) .
The flux density for diffusion described by (6) takes the form
J ( x , t ) = K 0 D t 1 ν c ( x , t ) x .
The subdiffusive generalization of Warburg’s impedance for a semi-infinite medium can be written in the form [30]
Z ( j ω ) = B ( i ω ) ( 1 ν / 2 ) .
Here, B is a frequency-independent parameter.
Based on fractional Equations (5) and (6), different generalizations of diffusion impedances can be obtained for different geometries and boundary conditions (see in [30,41,42,43,44] and references therein). In recent paper [45], effects produced by ion subdiffusion in lithium-ion cell are described within the fractional-order generalization of the single-particle model.
The similar derivation can be applied to subdiffusive Cahn–Hilliard equation with time-fractional derivatives. The subdiffusion of cations is characterized by fractional order α ,
0 C D t α c + = ( M + c + μ + q + M + c + E )
and transport of anions is characterized by fractional order β
0 C D t β c = ( M c μ + q M c E ) .
Here,
0 C D t α c ( x , t ) = 1 Γ ( 1 α ) 0 t c ( x , τ ) / τ ( t τ ) α d τ , 0 < α < 1
is the fractional Caputo derivative.
Eliminating the electric field from the considered system of equations, we arrive at an equation with a combination of fractional derivatives:
1 M + 0 C D t α + 1 M 0 C D t β c = c μ a
Denote the time-operator in the left-hand side by D t α , β , and rewrite the equation in the compact form:
D t α , β c = c μ a .

4. Mono-Fractional Phase-Field Subdiffusion

In the present work, we consider the role of anomalous phase field diffusion on impedance of phase transition material. In application, we limit ourselves to the case α = β . With M = M + M / ( M + + M ) the latter equation can be rewritten in the form
0 C D t α c = M c μ a .
This phase-field diffusion equation is related to the following continuity equation and generalized Fick law:
c t = div J , J = 0 R L D t 1 α c μ ¯ a
Monofractional Equation (11) can also be used to describe transport of intercalant atoms in solid particles. For example, when lithium is intercalated into some materials (LiFePO 4 or LiMn 2 O 4 ), the separation of the lithiated and non-lithiated phases can take place. In [11], real-time tracking of lithium transport and conversion in FeF2 nanoparticles by nanoscale imaging reveals a surprisingly fast conversion process in individual particles, with a morphological evolution resembling spinodal decomposition. The description of spinodal decomposition is usually based on the Cahn–Hilliard equation [13].
Following Bisquert and Compte [46], we also consider another type of the continuity equation
0 R L D t α c = div J .
With ordinary relation for current J = M c μ a , this system of equations leads to the following time-fractional generalization of the CH-equation:
0 R L D t α c = M c μ a .
For this case, the number of diffusing particles is not conserved. For the one-dimensional case, the diffusion equation
0 R L D t α c ( x , t ) = D 2 c ( x , t ) x 2
has clear interpretation [29,30]. In the multiple trapping model, this equation describes evolution of delocalized carriers concentration. In [46], the case (14) is classified as ‘anomalous diffusion Ia’ (AD-Ia), while case (6) is denoted as ‘anomalous diffusion Ib’.
In the following section, we apply linearization procedure to derive the response to a constant input signal with small harmonic perturbation component.

5. Impedance of Fractional Phase-Field Diffusion

Let us consider definition of impedance in electrode where conduction is related to mobility of ions in ambipolar diffusion approximation. The problem is analogous to calculation of diffusion impedance performed by Warburg. In this study, we employ CH Equation (3) with effective mixing energy (4) that describes charge transport in ambipolar diffusion approximation. The CH equation is widely used for the analysis of first-order phase transitions in solid solutions [47]. The CH Equation (3) can be easily linearized, if we apply the linear expansion of nonlinear term
g c g c c M + 2 g c 2 c M δ c
where c M is the average concentration of ions, δ c c c M is the small perturbation of ions’ concentration provided by the applied periodic voltage U = φ ( 0 , t ) φ ( L , t ) . Leaving terms of the first order, we derive the fractional Cahn–Hilliard equation in the linearized form
0 D τ α δ c = * 2 k 1 δ c k 2 * 2 δ c .
Here, we introduce dimensionless variables and coefficients:
r * = r a , * = 1 a , τ = t c M M a k B T a 2 1 / α ,
k 1 = 1 k B T 2 g a c 2 | c M , k 2 = κ k B T a 2 .
The CH equation can describe formation of phases by mechanisms of spinodal decomposition or nucleation [37]. Therefore, we assume that c M has the value not far from the equilibrium concentration in area of stable and metastable states. Furthermore, the restriction defines the sign of k 1 that has positive values ( k 1 > 0 ). Then, we apply to electrode edge ( x = 0 ) periodic potential φ ( 0 , t ) e i ω τ . The potential at another edge of electrode ( x = L , L ) is supposed to be zero φ ( L , τ ) = 0 . We can expect that chemical potential and concentration of cations changes with the same frequency: μ ( x , t ) e i ω τ , δ c e i ω τ .
The Fourier transformation of Equation (16) gives
Λ ( ω ) δ c ˜ ( x , ω ) = * 2 k 1 δ c ˜ ( x , ω ) k 2 * 2 δ c ˜ ( x , ω )
with
Λ ( ω ) = i ω α .
Solution of Equation (17) has the form
δ c ˜ ( x , ω ) = C 1 cosh λ + x + C 2 sinh λ + x + B 1 cosh λ x + B 2 sinh λ x
with
λ + = 1 + 1 4 χ Λ ( ω ) 2 χ k 1 , λ = 1 1 4 χ Λ ( ω ) 2 χ k 1 .
Constants C 1 , C 2 , B 1 and B 2 are determined from the boundary conditions for a particular problem.
The applied potential difference determines the carrier concentration at one edge of the cell in which phase diffusion occurs. On the other side, carriers can be locked or ejected. These boundary conditions are commonly referred to as reflective and absorbing, respectively [30]. Following the work in [48], we consider the impedance of 1D system, with reflective boundary condition, and absorbing condition, but for subdiffusive Cahn–Hilliard equation with time-fractional derivatives. One of the boundary conditions usually includes the current. This current is related to the ion concentrations through non-Fickian law in the case of anomalous phase-field diffusion. As it follows from Equation (17), the Fourier image of current is determined by the following relation:
J ˜ ( ω ) = i ω Λ ( ω ) k 1 x δ c ˜ ( x , ω ) k 2 3 x 3 δ c ˜ ( x , ω ) .
The impedance function is determined by the following expression [42]:
Z ˜ ( ω ) = κ L δ μ ˜ ( ω ) J ˜ ( ω ) .
Impedance response is calculated here under the following initial conditions
δ c ( t 0 , x ) = δ c 0 ( t 0 , x ) + δ c ω ( t 0 , x ) = 0 ,
and the following boundary conditions. For the reflecting electrode, the boundary conditions are as follows:
x = 0 : δ c ω x = 0 , 3 δ c ω x 3 = 0 ,
x = L : δ c ω x = 0 ; k 2 i ω Λ ( ω ) 3 δ c ω x 3 = J ,
while for the absorbing electrode, we have
x = 0 : δ c ω = 0 , 2 δ c ω x 2 = 0 ,
x = L : δ c ω x = 0 ; k 2 i ω Λ ( ω ) 3 δ c ω x 3 = J .
As mentioned above, these boundary conditions are typical when considering diffusion impedance. In the case of phase diffusion, they are more variable for interpretation. In particular, consider the reflecting boundary condition. One of interpretations implies charge injection from the left flat electrode at point L, phase diffusion of ions in the region [ 0 , L ] with a blocking electrode at point 0. Another interpretation is based on the method of images, which is applicable since the considered equation is linearized, spatially symmetric and local. Diffusion area can be extended to [ L , L ] and the boundary condition at zero can be replaced by the following condition at point L .
x = L : δ c ω x = 0 ; k 2 i ω Λ ( ω ) 3 δ c ω x 3 = J .
These boundary conditions correspond to the problem of intercalation/deintercalation of intercalant particles (e.g., lithium and sodium) into electrode particle specified by interval [ L , L ] . The Cahn–Hilliard equation in this case describes phase-field diffusion of intercalated phase.
The generalized impedances take the form similar to those presented in [48], but the presence of Λ ( ω ) accounts for bipolar phase field subdiffusion. For the reflecting boundary condition, we obtain the following expressions for constants in solution (18),
C 2 = B 2 = 0 ,
C 1 = Λ ( ω ) i ω k 2 J L λ + 1 λ + 2 λ 2 1 sinh λ + L = Λ ( ω ) i ω k 2 J L λ + χ k 1 1 4 χ Λ ( ω ) 1 sinh λ + L ,
B 1 = Λ ( ω ) i ω k 2 J L λ 1 λ + 2 λ 2 1 sinh λ L = Λ ( ω ) i ω k 2 J L λ χ k 1 1 4 χ Λ ( ω ) 1 sinh λ L ,
λ + 2 λ 2 = 1 4 χ Λ ( ω ) χ k 1 .
Details of derivation for the integer-order case are provided in Appendix of work [48]. In the time-fractional case, the derivation looks very similar. For the reflecting case, we obtain
Z refl Λ ( ω ) i ω 1 4 χ Λ ( ω ) 1 + 1 4 χ Λ ( ω ) 3 / 2 coth l 2 χ 1 1 4 χ Λ ( ω ) 1 1 4 χ Λ ( ω ) 3 / 2 coth l 2 χ 1 + 1 4 χ Λ ( ω ) .
Here,
l = L / k 1 .
For the absorbing case, we have
Z abs Λ ( ω ) i ω 1 4 χ Λ ( ω ) 1 + 1 4 χ Λ ( ω ) 3 / 2 tanh l 2 χ 1 1 4 χ Λ ( ω ) 1 1 4 χ Λ ( ω ) 3 / 2 tanh l 2 χ 1 + 1 4 χ Λ ( ω ) .
Setting L , we arrive at the impedance of a semi-infinite medium:
Z inf Λ ( ω ) i ω 1 4 χ Λ ( ω ) 1 + 1 4 χ Λ ( ω ) 3 / 2 1 1 4 χ Λ ( ω ) 3 / 2 .
Figure 2 shows Nyquist plots corresponding to the fractional ( 0 < α = β 1 ) phase-field subdiffusion for different values of α (a) and χ (b). The slope of the lines is determined not only by fractional exponent α , but also by parameter χ . Warburg impedance is reproduced for χ = 0 and α = 1 .
Figure 3 demonstrates the Nyquist plots for different boundary conditions. The case α = β and γ + = γ = 1 is considered. There are different cases (a–d) for various values of χ and α with reflecting opposite electrode. The role of phase field behavior is enhanced for larger values of χ . When χ 0 , we arrive at ordinary diffusive ( α = 1 ) or subdiffusive ( α < 1 ) impedance responses. Subdiffusion parameter influences the shapes of arcs and slope of low-voltage tail. The phase-field character of diffusion (with increasing χ ) makes the transitions between the regions more abrupt.
The classical finite-length Warburg impedances take place at χ = 0 and α = 1 . For normal diffusion, Λ ( ω ) = i ω . If χ = 0 , for the reflecting and absorbing cases, we have
Z refl Λ ( ω ) i ω coth l Λ ( ω ) , Z abs Λ ( ω ) i ω tanh l Λ ( ω )
with
Λ ( ω ) = γ + i ω α + γ i ω β .
In other words, we arrive at the finite-length impedance for bipolar subdiffusion.
In case of monopolar subdiffusion α = β , we arrive at the impedance obtained by Bisquert and Compte [46] for subdiffusion,
Z refl ( i ω ) α / 2 1 coth ( i ω / ω d ) α / 2 , Z abs ( i ω ) α / 2 1 tanh ( i ω / ω d ) α / 2 ,
where ω d is a constant dependent on l and subdiffusion coefficient.
Thus, phase-field analogs of anomalous diffusion models AD-Ib and AD-Ia considered by Bisquert and Compte [46] can be described by the following equations:
0 C D t α c = M c μ a , 0 R L D t α c = M c μ a .
We denote these fractional generalizations of the Cahn–Hilliard (CH) equation as C-CH and RL-CH, respectively, where the first letter denotes the type of fractional time derivative used in the model (Caputo or Riemann–Liouville). The corresponding impedances can be written in the following form:
Z C CH Λ ( ω ) i ω F ( ω ) , Z RL CH 1 Λ ( ω ) F ( ω )
with
Λ ( ω ) = ( i ω ) α
and functions F ( ω ) to be specified for certain type of boundary condition. For reflecting boundary, we take
F refl ( ω ) = 1 + 1 4 χ Λ ( ω ) 3 / 2 coth l 2 χ 1 1 4 χ Λ ( ω ) 1 1 4 χ Λ ( ω ) 3 / 2 coth l 2 χ 1 + 1 4 χ Λ ( ω ) 1 4 χ Λ ( ω ) .
and for absorbing boundary, it is
F abs ( ω ) = 1 + 1 4 χ Λ ( ω ) 3 / 2 tanh l 2 χ 1 1 4 χ Λ ( ω ) 1 1 4 χ Λ ( ω ) 3 / 2 tanh l 2 χ 1 + 1 4 χ Λ ( ω ) 1 4 χ Λ ( ω ) .
Nyquist plots for Z C CH ( ω ) are shown in Figure 2 and Figure 3 and described above. Nyquist plots corresponding to the RL-CH model with reflecting boundary condition for different values of α and χ are shown in Figure 4. The slope of the low-frequency part of the Nyquist diagram is determined not only by the subdiffusion exponent, but also by the χ parameter.

6. Application to Pseudo-Capacitors with Polyaniline/Carbon Nanotube Composite Electrodes

In this section, we consider an example of application of the developed model to description of impedance spectra for electrochemical supercapacitors (ESCs) with polyaniline/carbon nanotube composite electrodes. ESCs are high-capacity energy storage devices for applications requiring high power density. They are distinguished by charge storage mechanism: (1) electric double-layer (EDL) supercapacitors, where the capacitance arises due to the accumulation of charges on the electrode-electrolyte interface; (2) pseudo-capacitors with fast reversible Faraday redox reactions at the electrodes; and (3) hybrid supercapacitors that use EDL and Faraday capacitive mechanisms. To increase the active area in supercapacitors, electrodes with a highly developed surface are used. These electrodes can include porous materials based on activated carbon, carbon nanotubes, graphene oxide, foamed metals or conducting polymers and other materials.
In [44], we studied the behavior of planar microsupercapacitors (PMSC) with electrodes of MWCNT array under sinusoidal excitation, step voltage input and linear voltage input. The fractional differential model for planar supercapacitor proposed in [44] successfully describes the measured impedance spectra. Charge storage mechanism in PMSC implies the EDL formation, and the model used in [44] is based on the linear fractional diffusion equation. Here, we apply fractional phase-field diffusion model to pseudo-capacitors with PANI/MWCNT electrodes.
Polyaniline (PANI) was one of the first electroactive materials that opened the way to the fundamental and practical development of pseudo-capacitors [49]. In recent times, PANI has attracted additional attention due to its invaluable advantages in composition with other electroactive nanomaterials [50]. Short description of the system considered in our work is given in [51].
Electrodes of the pseudo-capacitor studied in this work contain vertically aligned MWCNT array (see Figure 5 and Figure 6a). The organic compound was grown by the pyrolysis method. The area of the plates is of 0.5 cm 2 . The MWCNT layer was covered with a thin layer of polyaniline (emeraldine form), obtained by the chemical method of oxidation of an aniline solution. The thickness of PANI layer is ~250 nm. A solution of polyvinyl alcohol (PVA) and phosphoric acid was used as a layer between the plates. A schematic representation of the structure is presented in Figure 6a. The parameters of the studied system were estimated by cyclic voltammetry and impedance spectroscopy. The measurements were made with the potentiostat-galvanostat P-40X (Electrochemical Instruments). Voltammograms are measured in the range from 0.5 V to 0.5 V at a potential sweep speed of 20 and 100 mV/s. Impedance spectra were taken in the frequency range from 0.1 Hz to 50 kHz with an amplitude of 50 mV.
Typical voltammogram for the studied pseudo-capacitor is demonstrated in Figure 6b. The capacity determined from voltammogram is of 6.73 × 10 2 F.
We fit measured frequency dependencies of the PANI-MWCNT pseudo-capacitor impedance by a simple equivalent circuit containing elements Z refl C CH and Z refl RL CH defined by (22). Reasonable goodness of fit has been achieved for the simple model presented in Figure 7. Fitted parameters are following: α = β = 0.8 , χ = 7 , l = 1 , C 0 = 10 4 F, R = 3 Ohm and W = 7 Ohm·Hz 1 / 2 , L = 2.5   μ m. Changing Z refl RL CH by Z refl C CH makes matching difficult and requires significant complication of the equivalent circuit and supplementing the model by additional elements. Comparison of measured and fitted spectra is demonstrated in Figure 8.
The used equivalent circuit can be reasonably interpreted. In the configuration of the supercapacitor shown in Figure 6a, the transport can be modeled by a one-dimensional transfer model. The capacity of EDL formed around the MWCNTs is modeled by ideal capacitor. The open Warburg impedance describes diffusion in the interelectrode space. The series resistor takes into account the resistance in the tubes, polymer and electrolyte. Element Z refl RL CH describes phase-field ion diffusion in PANI filling the MWCNT array with reflecting boundary condition at the base of nanotube array. The RL-CH model implies non-conserving ion density that can be related to the EDL formation by fraction of ions during phase diffusion inside the MWCNT array.
The spectra modeled by circuit presented in Figure 7 for listed parameters are sensitive to impact of phase-field diffusion, that is confirmed by variation of Bode plots for different values of χ (Figure 9). For chosen parameters, for small χ values, the peak in frequency dependence of the phase shift disappears, and the height of this peak can indicate the impact of phase-field transport. However, the latter statement requires verification and separate study under different conditions and for a larger number of systems.

7. Conclusions

In the present work, within the framework of the phase field approach based on the fractional Cahn–Hilliard equation, the model of diffusion formation of the electric double layer is revised and the frequency dependences of the generalized Warburg impedance are derived for the case of ion transport controlled by the subdiffusion Cahn–Hilliard equation with fractional derivatives. We studied charge transfer in an electrode, where conductivity is provided by monopolar or bipolar ionic charge transfer. The driving force is the chemical potential of ions, which is described in terms of the phase field, that allows us to avoid additional calculation of the activity coefficient. We apply linearization procedure to derive the response to a constant input signal with small harmonic perturbation component, and nonlinear fractional Cahn–Hilliard equation becomes linear. Then, we solve it by means of the integral transformation method. Impedance responses for a semi-infinite medium, and for finite-length cell with reflecting and absorbing boundary conditions for time-fractional Cahn–Hilliard equation are calculated and described.
The model is applied to description of impedance spectra of pseudo-capacitors with polyaniline/carbon nanotube composite electrodes. We fitted measured frequency dependencies of the PANI-MWCNT pseudo-capacitor impedance by a simple equivalent circuit containing elements Z refl C CH and Z refl RL CH defined by (22). Reasonable goodness of fit has been achieved for the simple model presented in Figure 7, that can be provided by clear interpretation. The configuration of our pseudo-capacitor implies one-dimensional transfer model. The capacity of EDL formed around MWCNTs is modeled by ideal capacitor. The open Warburg impedance describes diffusion in interelectrode space. The series resistor takes into account the resistance in the tubes, polymer and electrolyte. Element Z refl RL CH describes phase-field ion diffusion in PANI filling the MWCNT array with reflecting boundary condition at the base of nanotube array. The RL-CH model implies non-conserving ion density that can be related to the EDL formation by fraction of ions during phase diffusion inside the MWCNT array. Changing of Z refl RL CH by Z refl C CH makes matching difficult and requires significant complication of the equivalent circuit and supplementing the model by additional elements. We show that the assumption of phase-field diffusion plays an important role in this approximation. The peak in frequency dependence of phase is sensitive to the χ -value and can indicate the impact of phase-field diffusion.

Author Contributions

Conceptualization, R.T.S. and P.E.L.; methodology, R.T.S. and P.E.L.; validation, I.O.Y. and E.P.K.; investigation, R.T.S. and P.E.L.; resources, I.O.Y. and E.P.K.; data curation, E.P.K.; writing—original draft preparation, R.T.S. and P.E.L.; writing—review and editing, R.T.S. and P.E.L.; visualization, I.O.Y. and E.P.K.; project administration, R.T.S. and E.P.K.; funding acquisition, R.T.S. and E.P.K. All authors have read and agreed to the published version of the manuscript.

Funding

P.E.L., R.T.S. and I.O.Y. acknowledge the support from the Russian Science Foundation (project No. 19-71-10063), and E.P.K. acknowledges the support from the Ministry of Science and Higher Education of the Russian Federation (project No. 0591-2021-0002) and the Russian Federation President program (grant MK-1823.2021.4).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barsoukov, E.; Macdonald, J.R. Impedance Spectroscopy Theory, Experiment, and Applications, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
  2. Orazem, M.E.; Tribollet, B. Electrochemical Impedance Spectroscopy; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2008; pp. 383–389. [Google Scholar]
  3. Taylor, S.R.; Gileadi, E. Physical interpretation of the Warburg impedance. Corrosion 1995, 51, 664–671. [Google Scholar] [CrossRef]
  4. Huang, J. Diffusion impedance of electroactive materials, electrolytic solutions and porous electrodes: Warburg impedance and beyond. Electrochim. Acta 2018, 281, 170–188. [Google Scholar] [CrossRef]
  5. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  6. Uchaikin, V.V.; Sibatov, R. Fractional Kinetics in Solids: Anomalous Charge Transport in Semiconductors, Dielectrics, and Nanosystems; World Scientific: Singapore, 2013. [Google Scholar]
  7. Guyer, J.E.; Boettinger, W.J.; Warren, J.A.; McFadden, G.B. Phase field modeling of electrochemistry. I. Equilibrium. Phys. Rev. E 2004, 69, 021603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Guyer, J.E.; Boettinger, W.J.; Warren, J.A.; McFadden, G.B. Phase field modeling of electrochemistry. II. Kinetics. Phys. Rev. E 2004, 69, 021604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ferguson, T.R.; Bazant, M.Z. Nonequilibrium thermodynamics of porous electrodes. J. Electrochem. Soc. 2012, 159, A1967. [Google Scholar] [CrossRef]
  10. Fleck, M.; Federmann, H.; Pogorelov, E. Phase-field modeling of Li-insertion kinetics in single LiFePO4-nano-particles for rechargeable Li-ion battery application. Comput. Mater. Sci. 2018, 153, 288–296. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, F.; Yu, H.C.; Chen, M.H.; Wu, L.; Pereira, N.; Thornton, K.; Van der Ven, A.; Zhu, Y.; Amatucci, G.G.; Graetz, J. Tracking lithium transport and electrochemical reactions in nanoparticles. Nat. Commun. 2012, 3, 1201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Cahn, J.W.; Hilliard, J.E. Free energy of a nonuniform system. I. Interfacial free energy. J. Chem. Phys. 1958, 28, 258–267. [Google Scholar] [CrossRef]
  13. Cahn, J.W. On spinodal decomposition. Acta Metall. 1961, 9, 795–801. [Google Scholar] [CrossRef]
  14. Chen, C.; Ciucci, F. Generalized transport model for phase transition with memory. Phys. Lett. A 2013, 377, 2668–2672. [Google Scholar] [CrossRef]
  15. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Integrals and Derivatives of Fractional Order and Some of Their Applications; Nauka i Tekhnika: Minsk, Belarus, 1987. (In Russian) [Google Scholar]
  16. Zaslavsky, G.M. Chaos, fractional kinetics, and anomalous transport. Phys. Rep. 2002, 371, 461–580. [Google Scholar] [CrossRef]
  17. Schneider, W.R.; Wyss, W. Fractional diffusion and wave equations. J. Math. Phys. 1989, 30, 134–144. [Google Scholar] [CrossRef]
  18. Mainardi, F. Fractional relaxation-oscillation and fractional diffusion-wave phenomena. Chaos Solitons Fractals 1996, 7, 1461–1477. [Google Scholar] [CrossRef]
  19. Hilfer, R. Fractional diffusion based on Riemann-Liouville fractional derivatives. J. Phys. Chem. B 2000, 104, 3914–3917. [Google Scholar] [CrossRef] [Green Version]
  20. Uchaikin, V.V.; Sibatov, R.T. Fractional differential kinetics of dispersive transport as the consequence of its self-similarity. JETP Lett. 2007, 86, 512–516. [Google Scholar] [CrossRef]
  21. Uchaikin, V.V.; Ambrozevich, A.S.; Sibatov, R.T.; Ambrozevich, S.A.; Morozova, E.V. Memory and nonlinear transport effects in charging-discharging of a supercapacitor. Tech. Phys. 2016, 61, 250–259. [Google Scholar] [CrossRef]
  22. Akagi, G.; Schimperna, G.; Segatti, A. Fractional Cahn–Hilliard, Allen–Cahn and porous medium equations. J. Differ. Equ. 2016, 261, 2935–2985. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, H.; Cheng, A.; Wang, H.; Zhao, J. Time-fractional Allen–Cahn and Cahn–Hilliard phase-field models and their numerical investigation. Comput. Math. Appl. 2018, 76, 1876–1892. [Google Scholar] [CrossRef]
  24. Tripathi, N.K.; Das, S.; Ong, S.H.; Jafari, H.; Al Qurashi, M.M. Solution of time-fractional Cahn–Hilliard equation with reaction term using homotopy analysis method. Adv. Mech. Eng. 2017, 9, 1687814017740773. [Google Scholar] [CrossRef]
  25. Nigmatullin, R.R. The realization of the generalized transfer equation in a medium with fractal geometry. Phys. Status Solidi (b) 1986, 133, 425–430. [Google Scholar] [CrossRef]
  26. Nigmatullin, R.R. Fractional integral and its physical interpretation. Theor. Math. Phys. 1992, 90, 242–251. [Google Scholar] [CrossRef]
  27. Scher, H.; Montroll, E.W. Anomalous transit-time dispersion in amorphous solids. Phys. Rev. B 1975, 12, 2455. [Google Scholar] [CrossRef]
  28. Saichev, A.I.; Zaslavsky, G.M. Fractional kinetic equations: Solutions and applications. Chaos Interdiscip. J. Nonlinear Sci. 1997, 7, 753–764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Sibatov, R.T.; Uchaikin, V.V. Fractional differential kinetics of charge transport in unordered semiconductors. Semiconductors 2007, 41, 335–340. [Google Scholar] [CrossRef]
  30. Bisquert, J. Fractional diffusion in the multiple-trapping regime and revision of the equivalence with the continuous-time random walk. Phys. Rev. Lett. 2003, 91, 010602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Sibatov, R.; Shulezhko, V.; Svetukhin, V. Fractional derivative phenomenology of percolative phonon-assisted hopping in two-dimensional disordered systems. Entropy 2017, 19, 463. [Google Scholar] [CrossRef]
  32. Ainsworth, M.; Mao, Z. Analysis and approximation of a fractional Cahn–Hilliard equation. SIAM J. Numer. Anal. 2017, 55, 1689–1718. [Google Scholar] [CrossRef] [Green Version]
  33. Weng, Z.; Zhai, S.; Feng, X. A Fourier spectral method for fractional-in-space Cahn–Hilliard equation. Appl. Math. Model. 2017, 42, 462–477. [Google Scholar] [CrossRef]
  34. Caputo, M.; Fabrizio, M. Damage and fatigue described by a fractional derivative model. J. Comput. Phys. 2015, 293, 400–408. [Google Scholar] [CrossRef]
  35. Xia, H.; Tang, G.; Ma, J.; Hao, D.; Xun, Z. Scaling behaviour of the time-fractional Kardar-Parisi-Zhang equation. J. Phys. Math. Theor. 2011, 44, 275003. [Google Scholar] [CrossRef]
  36. L’vov, P.E.; Svetukhin, V.V. Simulation of the first order phase transitions in binary alloys with variable mobility. Model. Simul. Mater. Sci. Eng. 2017, 25, 075006. [Google Scholar] [CrossRef]
  37. Provatas, N.; Elder, K. Phase-Field Methods in Materials Science and Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2011. [Google Scholar]
  38. Barkai, E.; Metzler, R.; Klafter, J. From continuous time random walks to the fractional Fokker–Planck equation. Phys. Rev. E 2000, 61, 132. [Google Scholar] [CrossRef] [PubMed]
  39. Barkai, E. Fractional Fokker-Planck equation, solution, and application. Phys. Rev. E 2001, 63, 046118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Metzler, R.; Barkai, E.; Klafter, J. Anomalous diffusion and relaxation close to thermal equilibrium: A fractional Fokker-Planck equation approach. Phys. Rev. Lett. 1999, 82, 3563. [Google Scholar] [CrossRef] [Green Version]
  41. Macdonald, J.R.; Evangelista, L.R.; Lenzi, E.K.; Barbero, G. Comparison of impedance spectroscopy expressions and responses of alternate anomalous Poisson? Nernst? Planck diffusion equations for finite-length situations. J. Phys. Chem. C 2011, 115, 7648–7655. [Google Scholar] [CrossRef]
  42. Ciuchi, F.; Mazzulla, A.; Scaramuzza, N.; Lenzi, E.K.; Evangelista, L.R. Fractional diffusion equation and the electrical impedance: Experimental evidence in liquid-crystalline cells. J. Phys. Chem. C 2012, 116, 8773–8777. [Google Scholar] [CrossRef]
  43. Sibatov, R.T.; Uchaikin, V.V. Fractional kinetics of charge carriers in supercapacitors. In Volume 8 Applications in Engineering, Life and Social Sciences, Part B; De Gruyter: Berlin, Germany, 2019; pp. 87–118. [Google Scholar]
  44. Kitsyuk, E.P.; Sibatov, R.T.; Svetukhin, V.V. Memory effect and fractional differential dynamics in planar microsupercapacitors based on multiwalled carbon nanotube arrays. Energies 2020, 13, 213. [Google Scholar] [CrossRef] [Green Version]
  45. Sibatov, R.T.; Svetukhin, V.V.; Kitsyuk, E.P.; Pavlov, A.A. Fractional differential generalization of the single particle model of a lithium-ion cell. Electronics 2019, 8, 650. [Google Scholar] [CrossRef] [Green Version]
  46. Bisquert, J.; Compte, A. Theory of the electrochemical impedance of anomalous diffusion. J. Electroanal. Chem. 2001, 499, 112–120. [Google Scholar] [CrossRef]
  47. Pego, R.L. Front migration in the nonlinear Cahn-Hilliard equation. Proc. R. Soc. Lond. A Math. Phys. Sci. 1989, 422, 261–278. [Google Scholar]
  48. Ciucci, F.; Lai, W. Electrochemical impedance spectroscopy of phase transition materials. Electrochim. Acta 2012, 81, 205–216. [Google Scholar] [CrossRef]
  49. Eftekhari, A.; Li, L.; Yang, Y. Polyaniline supercapacitors. J. Power Sources 2017, 347, 86–107. [Google Scholar] [CrossRef]
  50. Gupta, V.; Miura, N. Polyaniline/single-wall carbon nanotube (PANI/SWCNT) composites for high performance supercapacitors. Electrochim. Acta 2006, 52, 1721–1726. [Google Scholar] [CrossRef]
  51. Yavtushenko, I.O.; Sibatov, R.T.; Somov, A.I.; Makhmud-Akhunov, M.Y. Fractional circuit model for supercapacitors with polyaniline/carbon nanotube composite-based electrodes. J. Phys. Conf. Ser. 2020, 1695, 012039. [Google Scholar] [CrossRef]
Figure 1. Phase diagram for binary solution (1, 2) and binary solution with respect to electronic subsystem in approximation of ambipolar diffusion (3, 4). Curves 1 and 3 show the temperature dependence of equilibrium composition. The dashed lines 2 and 4 show the metastability limit.
Figure 1. Phase diagram for binary solution (1, 2) and binary solution with respect to electronic subsystem in approximation of ambipolar diffusion (3, 4). Curves 1 and 3 show the temperature dependence of equilibrium composition. The dashed lines 2 and 4 show the metastability limit.
Fractalfract 05 00191 g001
Figure 2. Nyquist plots corresponding to the fractional ( 0 < α = β 1 ) phase-field subdiffusion for different values of α (a) and χ (b).
Figure 2. Nyquist plots corresponding to the fractional ( 0 < α = β 1 ) phase-field subdiffusion for different values of α (a) and χ (b).
Fractalfract 05 00191 g002
Figure 3. Nyquist plots corresponding to the phase-field subdiffusion ( α = β , χ = 1 ) for different values of α with reflecting boundary condition, l = 1 : (a) α = 1.0 , (b) α = 0.9 , (c) α = 0.8 .
Figure 3. Nyquist plots corresponding to the phase-field subdiffusion ( α = β , χ = 1 ) for different values of α with reflecting boundary condition, l = 1 : (a) α = 1.0 , (b) α = 0.9 , (c) α = 0.8 .
Fractalfract 05 00191 g003
Figure 4. Nyquist plots corresponding to the RL-CH model for different values of α (a) and χ (b) with reflecting boundary condition ( l = 1 ).
Figure 4. Nyquist plots corresponding to the RL-CH model for different values of α (a) and χ (b) with reflecting boundary condition ( l = 1 ).
Fractalfract 05 00191 g004
Figure 5. SEM image of a grown MWCNT array on a titanium plate used by us for preparation of PANI-MWCNT pseudo-capacitor.
Figure 5. SEM image of a grown MWCNT array on a titanium plate used by us for preparation of PANI-MWCNT pseudo-capacitor.
Fractalfract 05 00191 g005
Figure 6. Schematic representation of the cell with polyaniline/vertically aligned carbon nanotube array electrodes (a,b). Voltammogram for pseudocapacitors for potential sweep speed of 100 mV/s.
Figure 6. Schematic representation of the cell with polyaniline/vertically aligned carbon nanotube array electrodes (a,b). Voltammogram for pseudocapacitors for potential sweep speed of 100 mV/s.
Fractalfract 05 00191 g006
Figure 7. Equivalent circuit model of PANI-MWCNT pseudo-capacitor. Element Z refl RL is described by impedance from (22) for the case of the RL-CH model with reflecting boundary condition.
Figure 7. Equivalent circuit model of PANI-MWCNT pseudo-capacitor. Element Z refl RL is described by impedance from (22) for the case of the RL-CH model with reflecting boundary condition.
Fractalfract 05 00191 g007
Figure 8. Nyquist (a) and Bode (b) plots of measured and modeled impedance spectra for PANI-MWCNT pseudo-capacitor at room temperature.
Figure 8. Nyquist (a) and Bode (b) plots of measured and modeled impedance spectra for PANI-MWCNT pseudo-capacitor at room temperature.
Fractalfract 05 00191 g008
Figure 9. Bode plots for the equivalent circuit with Z refl RL CH for different values of χ .
Figure 9. Bode plots for the equivalent circuit with Z refl RL CH for different values of χ .
Fractalfract 05 00191 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

L’vov, P.E.; Sibatov, R.T.; Yavtushenko, I.O.; Kitsyuk, E.P. Time-Fractional Phase Field Model of Electrochemical Impedance. Fractal Fract. 2021, 5, 191. https://doi.org/10.3390/fractalfract5040191

AMA Style

L’vov PE, Sibatov RT, Yavtushenko IO, Kitsyuk EP. Time-Fractional Phase Field Model of Electrochemical Impedance. Fractal and Fractional. 2021; 5(4):191. https://doi.org/10.3390/fractalfract5040191

Chicago/Turabian Style

L’vov, Pavel E., Renat T. Sibatov, Igor O. Yavtushenko, and Evgeny P. Kitsyuk. 2021. "Time-Fractional Phase Field Model of Electrochemical Impedance" Fractal and Fractional 5, no. 4: 191. https://doi.org/10.3390/fractalfract5040191

APA Style

L’vov, P. E., Sibatov, R. T., Yavtushenko, I. O., & Kitsyuk, E. P. (2021). Time-Fractional Phase Field Model of Electrochemical Impedance. Fractal and Fractional, 5(4), 191. https://doi.org/10.3390/fractalfract5040191

Article Metrics

Back to TopTop