Time-Fractional Phase Field Model of Electrochemical Impedance

: In this paper, electrochemical impedance responses of subdiffusive phase transition materials are calculated and analyzed for one-dimensional cell with reﬂecting and absorbing boundary conditions. The description is based on the generalization of the diffusive Warburg impedance within the fractional phase ﬁeld 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 ﬁeld allowing us to avoid additional calculation of the activity coefﬁcient. The derived impedance spectra are applied to describe the response of supercapacitors with polyaniline/carbon nanotube electrodes.


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 pseudocapacitors with polyaniline/carbon nanotube electrodes.

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.
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.

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]) 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]: 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 Thus, combining the previous relations we obtain the equations describing the diffusion of ions and electrons in electrolyte: 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: where M a = M + M − M + +M − is the ambipolar mobility, and µ a is the effective (ambipolar) chemical potential that can be expressed as Here, we introduce the effective mixing energy g a = g(c) + k B Tc 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 where Ω is interaction parameter that is ordinarily related to critical temperature Ω = 2k 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 2 0 /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 Ca ≈ 0.69T 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]: where c(x, t) is particle concentration, V(x) is an external potential, K is a generalized diffusion coefficient and kT is the thermal energy. Here, 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: The Fourier transform of Equation (6) leads to the expression The flux density for diffusion described by (6) takes the form The subdiffusive generalization of Warburg's impedance for a semi-infinite medium can be written in the form [30] 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 α, Eliminating the electric field from the considered system of equations, we arrive at an equation with a combination of fractional derivatives:

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 α = β.
This phase-field diffusion equation is related to the following continuity equation and generalized Fick law: 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 With ordinary relation for current J = Mc∇µ a , this system of equations leads to the following time-fractional generalization of the CH-equation: For this case, the number of diffusing particles is not conserved. For the one-dimensional case, the diffusion equation 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.

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) 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).
Here, we introduce dimensionless variables and coefficients: 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ωτ .

Solution of Equation (17) has the form
with 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: The impedance function is determined by the following expression [42]: 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: while for the absorbing electrode, we have 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.
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), , , 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 Here, l = L/ k 1 .
For the absorbing case, we have Setting L → ∞, we arrive at the impedance of a semi-infinite medium: 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 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, 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: 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: with Λ(ω) = (iω) α and functions F(ω) to be specified for certain type of boundary condition. For reflecting boundary, we take and for absorbing boundary, it is Nyquist plots for Z C−CH (ω) are shown in Figures 2 and 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.

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) pseudocapacitors 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 Figures 5 and 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 C−CH refl and Z RL−CH refl defined by (22). Reasonable goodness of fit has been achieved for the simple model presented in Figure 7. 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 RL−CH refl 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.

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 C−CH refl and Z RL−CH refl 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 RL−CH refl 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 RL−CH refl by Z C−CH refl 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.