3.1. A Fractal Perspective of the EIS of GHDPE/ KCl Aqueous Solution Composite
The interfacial region is formed when the electrical circuit in the cell is closed and then, the electrochemical potential gradient associated to the charge migration remains commonly normal to the electrode surface. Therefore, the fractal dimension () of the electrode surface must be of the same order as those calculated through this Equation Error! Reference source not found. for diffusion-controlled processes. In the limiting case, for a well-polished inert metallic electrode the fractal dimension of the electrode surface would approximate the Euclidean dimension .
Our hypothesis leads to
of the electrode/solution surface is associated to the intersection of two transport cluster (
Figure 1). One corresponds to the flow of electrons through the internal cluster of the condensed composite material fractal dimension (
dC) and, the other, to the ion transport of formation of the interfacial region (double layer). The electrochemical surface is the intersection of the composite percolative cluster with the interfacial cluster formed by the vector gradients of the electrochemical potentials at the interfacial region.
Hence, the fractal dimension d
F of the surface electrode is defined as the intersection of a 2D plane embedded into three-dimensional object
[
59] as is the studied electrode:
The fractal dimension of the polymer+graphite composite material
depends increasingly on the probability of occupancy of the particles within the material p, i.e., on their volume content
(
Figure 2). Therefore, consistent with Equation (3) , the fractal dimension of the fractal dimension of the electrode surface
should increase with graphite content to a value close to 2.
A modulated ac perturbation on the composite-solution system allows to obtain valuable information on the electrical percolation through the system at each stabilization potential and applied frequency . As is commonly accepted, alternating current flows through the composite material of the working electrode, by three different ways: polarization, tunnel effect and directly by electron-hopping through the clusters formed by the conductive particles in contact.
By perturbing the dissolution-electrode system with an applied potential at a given frequency (
:
the modulated linear electrochemical response is out of phase in time with a phase angle (
) according to the sinusoidal intensity function:
Immediately, one can conclude that the measure of phase angle points of the nature of the percolation of ac perturbation through the material:
decreases from a limit value of
rad in the absence of graphite to values typical of a resistive behavior when the content of conducting particles is greater than the first percolation threshold [
4]. This implies that the impedance transfer function
provides information on the conductivity of the composite material and the permittivity
of the interfacial region. Then, both electrical magnitudes are dependent on the volumetric content of graphite
. Therefore, the dependences of real and imaginary components of impedance shows an interesting direct information about the system. However, the electrochemical processes that occur in the presence of electroactive substances, including the presence of oxygen in the solution, complicate the interpretation of the electrochemical impedance of the system (
Figure S3, results shown as Supplementary Material).
of graphite: 0.185 (black), 0.229 (red), 0.267 (blue), 0.308 (pink), 0.356 (purple). E0 = 0.4V; ΔE = 10 mV; 1M KCl. T = 298K.
The experimental values are fitted to a simple equivalent circuit R-CPE:
where
, and
if the value of the ohmic drop through the cell is negligible compared to that of the composite material while
and
are fitting coefficients of the real
Z’ and imaginary
Z” of the impedance:
The fractal dimension is related to this exponent but it depends on the geometry of the cell and its elements [
55] and the electrode model considered as well as the particular electrode process. In any case,
must be related to the CPE exponent
which, is due to the existence of a time constants distribution on the electroactive surface, as is generally accepted [
60,
61,
62]. Although this topic is extraordinarily complex, a simple relationship could be postulated as a first approximation to the problem for a planar polished electrode and if it is sufficiently far from the another parallel plane counter electrode:
Then, the postulated Equation (9) is consistent with the insertion of the electrochemical surface of fractal dimension with a virtual plane of the interfacial region or double layer embedded in a 3D space between the surface and a virtual second plane of the interfacial capacitor. In this case the electrode surface was placed horizontally to avoid gravitational convection. A platinum, foil placed parallel to the working electrode was also used as an auxiliary electrode. In this case the fractal dimension of the cluster formed by the lines of force in the interfacial region, being a monodirectional transport by migration, .
From other perspective, the interfacial region can be modelized of a framework of graphite particles in contact with the solution in a surface 2D of a geometrical fractal dimension. Also, this interfacial region can be modelized by means an electrical framework of micro-capacitors of an overall ≃ because above the percolation threshold = 0.14 the values of n are close to unity.
Assuming that the electrode surface consists of squares, whose minimum size is
and maximum
, and assuming that the double layer capacities
associated with the electrode/dissolution interface are proportional to the number of graphite particles on the surface and their occupied area
, the fractality of the electrode surface can be expressed as follows according to the scaling unit λ (an arbitrary value of area):
It can be inferred that
, which is the same as
As the capacitance is an extensive magnitude proportional to the surface area it is easy to induce the Equation (11) since in the limit case of an ideal polish electrode surface = 2 and n = 1, corresponding a pure capacitor of de double layer.
Despite GHDPE not shows hygroscopic chemical interaction nor porous, this hypothesis is quite controversial and represents a drastic approach proposed for this particular model under study. In fact, the correlation between the fractal dimension of fractal electrodes with the exponent of the CPE has been rigorously considered for different electrode models in the literature (
Table 1).
The Equation (3) supposes that the interfacial region has an apparent percolative structure formed by the resultant vector of the electrochemical potential gradient, which is considered normal to the plane containing the electroactive graphite particles of the composite surface in contact with the solution. Therefore, the values of are related to the electroactive area occupied , which is theoretically lesser than the Euclidean geometric one for a well-polished electrode , because the heterogenous electrode surface forme by conducting graphite particles dispersed between non-conductive areas of plastic
Following the same ideal way, it is possible to propose a second hypothesis: If we consider that the associated dimension fractal of the interfacial region depends on the capacitance distribution, analogously in Equation (9) the exponent of the CPE can be considered a fractal dimension of the intersection cluster of the electroactive surface with the cluster of electric charges of the plane-parallel equivalent capacitor, both embedded in a raised 3D space. This assumption considers that capacitance associated to the interfacial region
is the magnitude of the experimental observation, however more general is considerer into acount the dimensions of the magnitudes. Physical objects [
59] considered fractal are defined as fractal geometric objects in a range [
,
]. The number of elements
of a range is proportional to
:
but the concept can be extended to capacitance quantities. If it is accepted that there are a number of squares on the electrode surface
contributing to the experimental capacitance in each range of graphite content in the composite analogously to the equation will be fulfilled:
It follows to the Equation(9).
From this perspective, the fractal dimension of the interfacial region is that of the cluster associated with the electrochemical potential gradients normal to the electrode surface:
This means that by adjusting the impedance measurements to an equivalent circuit consisting of a resistance in series with the CPA, it is possible to estimate
values (
Table 2). In the extreme ideal case of a well-polished inert metal electrode, the
exponent would approach unity in accordance with double layer model of an ideally polarizable electrode, so the fractal dimension of the electrode surface it would coincides with the Euclidean fractal dimension of a surface
≈ 2. While in the case of the real composite, although the electrode surface is polished, the constant
of the CPA will depend on the area of the electroactive fraction as well as the on the contribution of the excess concentration of the electrolyte in the interfacial region, which depends also on the ionic strength of the electrolyte in consistence with the proposed model.
In the electrochemical set-up of the experiments carried out, in the absence of convection, and the counter electrode being a platinum foil parallel to the working electrode, the fractal dimension of the percolation cluster through the interfacial region is the intersection of the fractal of the surface with this migration transport cluster immersed in 3D of the solution.
3.2. Relation between the Experimental Imaginary Part of the Impedance and the Fractal Dimensions
Laser diffractometry has been used to measure the size dispersion of graphite particles. (
Figure 4).
= 16.2 ± 0.2 μm.
The first percolation threshold
≈ 0.14 have been calculated from the dependence of the resistance
on
measured by dc and EIS (See
Figure S2 in supplementary materials). The second percolation threshold corresponds approximately to the volume fraction of graphite that the material loses its consistency
Then, the theoretical occupational factor of GHDPE is
f ≈ (0.14+0.41) = 0.55, and the probability of percolation. This value agrees the theoretical 0.247±0.005 for a percolation of links for a simple cubic structure [
59,
81,
82]. This value is approximately the one obtained by a simulation. (
Figure 2,
Figure 5) Therefore, the coordination index,
z, of a percolation tree is related to the critical threshold by
z = 1/p
c. Substituting the experimental value
pc ≈ 0.254, we find for GHDPE that
z = 3.94 ≈ 4. This value is consistent with the calculated occupancy factor in accordance with the anisotropic conductance of the graphite particles.
= 16.2 μm , respectively. Only a limited number of surface graphite particles pertain a conductive cluster (red squares in b). The overall number of squares of the electrode surface is 92226. The overall number of cubes of the electrodes samples are approximately 56444026. Simulated = 0.154 (3D) and = 0.44 (2D) where calculated.
EIS provides valuable information in the case of the GHDPE system immersed in an electrolyte dissolution. From the real component of the impedance, extrapolating to infinite frequency, the values of could be calculated. From the module of the imaginary component, the capacitive character of the electrode is extracted.
From the limit at high frequencies of the real impedance component, the resistance
is calculated, but in addition, also the CPE parameters, which provide information on the distribution of potentials and time constants on the electrode surface, which is directly related to the fractal dimension
of the electrode surface.
increases as the proportion of surface occupied by graphite is greater, in agreement with the predictions of the Percolation Theory and the experimental data obtained by EIS (
Figure 6a). The parabolic curve starts from the point corresponding to the percolation threshold
≈ 0.14 to the theorical maximum capacitance associated with the double layer of the electrode-solution interface
on the volume percentage of graphite. b) Dependence on the ratio of the overall electrical active surface of the graphite particles
SCG and their separation
. At the same conditions of
Table 2.
From the values of the CPE coefficients for the samples with different contents of conductive filler, two qualitatively different ranges are distinguished over the percolation threshold: in a range which the values of
and
grow on the graphite content, and another, in which both parameters vary little (
Table 2 and
Figure 6b). Despite the probability of current passage across the bulk over
= 0.14 is
= 1 in both ranges, the potential dispersion at the multi-micro electrodes of the graphite particles on the surface of the composite is lesser over
= 0.25, showing that the surface reaches a single behavior in the interfacial region because the calculated
remains practically constant (
Table 2). These results are consistent with the initial hypothesis that the electrode surface is the result of the intersection of the cluster of the composite with the cluster associated at the perpendicular plane interfacial region.
It is necessary to consider the CPE to explain the capacitive behavior of composites due to the heterogeneity of the surface. In this way the magnitude
explains the percolation through the bulk of the composite, the CPE gives basic information on the heterogeneity of the surface, and hypothetically on the 2D distribution of time constants of the processes that take place in the interfacial region. The Equation
Error! Reference source not found. and Equation
Error! Reference source not found. are based on the fractal modeling of the three-dimensional percolative cluster of the solid composite
(see last column of
Table 2 ) and are close to the theoretically [
59] calculated for an infinite cluster
= 2.523.
If the coefficient has been associated with the fractal dimension of the surface of the electrode, the physical meaning of the constant is clear when the parameter tends to unity. For graphite contents, ≥ 0.25 the capacity of the double layer can be considered as corresponding to the equivalent associated
capacitors of
, formed by the individual graphite particles on the surface in contact with the aqueous KCl solution. Therefore, it can be simplified, assuming that the interfacial region shows a phenomenological capacitance that depends on an apparent dielectric permittivity ε of the dissolution and an apparent thickness
of the double layer, named frequently reciprocal Debye length:
being
, the area of a selected electroactive experimental surface of the electrode. Therefore, the value of this parameter can be considered proportional to the fraction of area occupied by the graphite particles located on the surface of the electrode. At low concentration of KCl, his Debye length varies with the square root of the inverse ionic strength if the interfacial region is considered as an equivalent parallel plate capacitor [
83,
84,
85].
In accordance with the assumption of Equation (9), the exponent of the CPA can be considered directly as a fractal dimension of the intersection cluster of the electrochemical surface cluster of
of with the electrical charges of the other capacitor plate of
, both embedded in a 3D space. That justifies the hypothesis embodied in Equation (7) in the simple model proposed, where the CPE is directly related to the electrical properties of the interfacial region only. However, if the CPE included in the equivalent circuit is directly associated to the solid electrode as well as the case of deposited films [
86] or porous and rough electrodes (
Table 2), the interpretation of the constant-phase-element is in general more complicated [
41,
87,
88]. In this composite, the real and imaginary components of the impedance depend on both CPE parameters following Equations (9) and (10). This suggest that the proposed model based on the self-definition of the R-CPE circuit is apparently useful for describing the electrochemical behavior of the composite electrode/solution for graphite contents
. But, this hypothesis, based in in the intuitive idea that the imaginary component of the impedance depends only of the interphase electrode/solution requires a further discussion, because the surface morphology of the electrode depends directly on the percolation inside the material and, therefore, in this system, lacking porosity and significant faradaic processes, the CPE associated to the imaginary part must depend on the distribution of graphite particles on the surface and also inside the material. That is, on the particle/solution interfaces at the surface, as well as to a large extent on the graphite particle/polyethylene interface inside the material.