# Compensating for Electrode Polarization in Dielectric Spectroscopy Studies of Colloidal Suspensions: Theoretical Assessment of Existing Methods

^{1}Environmental Fluid Mechanics, Faculty of Civil Engineering and Geosciences, Delft University of Technology, Delft, Netherlands^{2}Laboratoire PHENIX, Centre National de la Recherche Scientifique, Sorbonne Universités, UPMC Université Paris 06, Paris, France^{3}Departamento de Física Aplicada, Universidad de Granada, Granada, Spain^{4}Formerly affiliated with Leiden Institute of Chemistry, Leiden University, RA Leiden, Netherlands^{5}Department of Chemical Engineering, Delft University of Technology, Delft, Netherlands

Dielectric spectroscopy can be used to determine the dipole moment of colloidal particles from which important interfacial electrokinetic properties, for instance their zeta potential, can be deduced. Unfortunately, dielectric spectroscopy measurements are hampered by electrode polarization (EP). In this article, we review several procedures to compensate for this effect. First EP in electrolyte solutions is described: the complex conductivity is derived as function of frequency, for two cell geometries (planar and cylindrical) with blocking electrodes. The corresponding equivalent circuit for the electrolyte solution is given for each geometry. This equivalent circuit model is extended to suspensions. The complex conductivity of a suspension, in the presence of EP, is then calculated from the impedance. Different methods for compensating for EP are critically assessed, with the help of the theoretical findings. Their limit of validity is given in terms of characteristic frequencies. We can identify with one of these frequencies the frequency range within which data *uncorrected* for EP may be used to assess the dipole moment of colloidal particles. In order to extract this dipole moment from the measured data, two methods are reviewed: one is based on the use of existing models for the complex conductivity of suspensions, the other is the logarithmic derivative method. An extension to multiple relaxations of the logarithmic derivative method is proposed.

## 1. Introduction

Dielectric spectroscopy is a powerful tool to determine the electrokinetic properties of suspensions of nano- or microparticles as it can probe the suspension's response as function of the applied electric field frequency. Interfacial properties such as zeta potential and Stern layer conductance can then be derived by analyzing the dielectric spectra of the suspensions (Grosse et al., 1998; Delgado, 2002; Hollingsworth and Saville, 2003; Ohshima, 2006; Chassagne et al., 2009). It was shown recently that fitting simultaneously the dielectric spectrum of a suspension and the electrophoretic velocity of the particles composing this suspension can provide a unique set of parameters when a Stern layer conductance is necessary to fit the data in addition to the zeta potential (Chassagne et al., 2009). Most experiments on colloidal suspensions are performed in the “low frequency” regime, i.e., below 1 MHz, where the typical dispersions emerge that are associated with the colloidal particle and its double layer. Studies at higher frequencies (well above 1 MHz, see for instance Kaatze and Feldman, 2006) will not be addressed.

The determination of dielectric spectra is based on measurements of the complex conductivity of the suspension ${\stackrel{~}{K}}_{s}$ (or equivalently the complex dielectric response ${\stackrel{~}{\epsilon}}_{s}$) as a function of frequency. This is done by measuring the complex impedance ${\stackrel{~}{Z}}_{s}$ of a suspension contained in a cell with (generally two) electrodes. Although it is in principle a simple measurement, the proper determination of ${\stackrel{~}{Z}}_{s}$ is difficult due to several unwanted effects. Two types of effects are distinguished. One type of effect originates from the non-ideality of the experimental set-ups. Because of the presence of electrical circuits, wires, external components, so-called “stray impedances” will arise and have to be accounted for when the impedance of the cell, in which the suspension is inbedded, is measured. This type of effect will not be discussed here. Another type of unwanted effect stems from the distribution of the charged species in the measurement cell when an electric field is applied. It is the so-called “electrode polarization” (EP) that typically occurs at low frequencies, mostly below 10–100 kHz (Barsoukov and Macdonald, 2005). EP originates from the fact that at low frequencies ions are able to build up close to the electrodes, contributing with a large additional capacitance to the impedance ${\stackrel{~}{Z}}_{c,s}$ of the suspension, see Figure 1. The impedance ${\stackrel{~}{Z}}_{c,s}$ of the suspension measured in the cell (cleaned from all stray impedances that might have to be accounted for) can be seen as the sum of two impedances, the sought one from the bulk of the suspension ${\stackrel{~}{Z}}_{s}$ (“true” impedance) and an unwanted one due to EP, ${\stackrel{~}{Z}}_{EP}$:

As Hollingsworth (2013) recently pointed out, EP has been studied extensively. The impedance caused by the electrodes was first discussed by Kohlrausch around 1874 as a disturbance in conductivity measurements of electrolyte solutions. Theoretical work was performed by Warburg, Mandel, and Buck to name a few, see for instance Barsoukov and Macdonald (2005), Buck (1969), and van der Touw and Mandel (1971). Analytical expressions for the impedance by EP can be obtained by solving the set of relevant equations (Nernst-Planck, conservation of ions and Poisson). In most works (including the present article), the electrodes are assumed to be ideally polarizable, implying that faradaic reactions are negligible and hence that there is no charge transfer at the electrodes (blocking electrodes). We refer in particular to the work done in Hollingsworth and Saville (2003), Buck (1969), Cirkel et al. (1997), Chassagne et al. (2002), Chassagne et al. (2003) that will be used in this article. Theoretical investigation of electrode polarization with non-zero zeta potential at the electrodes has been performed numerically by White et al. (DeLacey and White, 1982), perturbation expansions have been done by Gunning et al. (1995) and semi-analytical solutions provided by Scott et al. (2000a,b).

**Figure 1. Schematic representation of the distribution of particles (ions, colloidal particles) at a given time, when an AC electric field is applied to a solution (x = e) or suspension (x = s)**. Due to the blocking nature of the electrodes, ions and colloidal particles build-up close to the electrodes at low frequencies. The equivalent impedance of the cell ${\stackrel{~}{Z}}_{c,x}$ is then made up of two contributions: the wanted “true” impedance ${\stackrel{~}{Z}}_{x}$ of the bulk of the investigated solution or suspension and the contribution due to the polarization of the electrodes ${\stackrel{~}{Z}}_{EP}.$

The measurement of the complex dielectric response of colloids is hampered seriously by EP due to the fact that the characteristic relaxation frequencies associated with the properties of the colloids show up in the low-frequency regime too. It is therefore important to be able to devise a suitable method to remove or minimize the strong EP contribution to the measured signal, so as to assess ${\stackrel{~}{Z}}_{s}$. Recent reviews of the existing methods have been given by Kaatze et al. (Kaatze and Feldman, 2006; Grosse and Delgado, 2010; Ishai et al., 2013). The present article is intended to give insight in these different methods and in particular check their range of validity, from a theoretical perspective. Characteristic frequencies are derived. We show in particular that one of these frequencies enables us to identify those cases where data *uncorrected* for EP can be used to assess ${\stackrel{~}{Z}}_{s}$.

From ${\stackrel{~}{Z}}_{s},$ the complex conductivity ${\stackrel{~}{K}}_{s}$ of the suspension can be directly derived. This conductivity is related to the interfacial properties of the colloidal particles such as zeta potential and Stern layer conductance via the dipole moments of the polarized particles and their double layers. This dipole moment is represented by the dipolar coefficient $\stackrel{~}{\beta}$ that can be modeled either analytically or numerically, see e.g., Chassagne and Bedeaux (2008), Mangelsdorf and White (1990), and Minor et al. (1998). We will show, using the analytical model presented in Supplementary Material 1 as an example, how such a model can be used to find the interfacial properties of the particles by fitting data uncorrected for EP. In a next section, we will show how an extended logarithmic derivative method allows the evaluation of $\stackrel{~}{\beta}$ from uncorrected data. In the case the data should be corrected for EP, one of the different methods to account for it should be applied first.

### 1.1. Outline of the Article

In Section 2, the purpose of dielectric spectroscopy in the context of colloids will be outlined. The complex conductivity ${\stackrel{~}{K}}_{s}$ of colloidal suspensions will be given as function of the complex dipolar coefficient $\stackrel{~}{\beta}$ as basic unknown. It is this quantity that characterizes important particles properties such as zeta potential and Stern layer conductances. Before discussing methods to find ${\stackrel{~}{K}}_{s}$ (and subsequently $\stackrel{~}{\beta}$) theoretical derivations are recalled: in Section 3, a description of EP for blocking electrodes will be given. This theory is an extension of previous ones (Buck, 1969; Cirkel et al., 1997; Chassagne et al., 2002, 2003), and is valid for the whole frequency range of interest and arbitrary mobilities and valences of ions. Two types of electrodes will be considered: planar and cylindrical, which correspond to the geometries most frequently used. Important results of the derivations given in Section 3 will be needed in Section 4. The impedance of the electrolyte (or suspension) is often represented as an equivalent circuit consisting of the combination of resistances and capacitances (Buck, 1969; Barsoukov and Macdonald, 2005). This formalism will be discussed in Section 4. It will in particular enable us to model the response of a colloidal suspension in the presence of EP. In Section 5, several methods to (a) correct for EP and (b) extract $\stackrel{~}{\beta}$ from the data will be presented and discussed. The Supplementary Material sections give details about: the mathematical derivations of some of the basic formulas reviewed, relevant aspects of the equivalent circuit and of the spatial profile of the alternating electric field within an electrolytic solution.

### 1.2. Variable Definitions

${\stackrel{~}{Z}}_{c,x}$ is the impedance of the measurement cell filled with a suspension (*x* = *s*), an electrolyte solution (*x* = *e*) or a reference electrolyte solution (*x* = *r*) (all in the presence of electrode polarization).

${\stackrel{~}{Z}}_{e}$ is the “true” impedance of the electrolyte solution (without electrode polarization).

${\stackrel{~}{Z}}_{s}$ is the “true” impedance of the suspension (without electrode polarization).

${\stackrel{~}{Z}}_{r}$ is the “true” impedance of the reference electrolyte solution (without electrode polarization).

${\stackrel{~}{Z}}_{EP}$ is the impedance accounting for electrode polarization (${\stackrel{~}{Z}}_{EP}={\stackrel{~}{Z}}_{c,x}-{\stackrel{~}{Z}}_{x}$).

These subscripts apply accordingly to the related complex conductivities and complex permittivities ($\stackrel{~}{K}$, $\stackrel{~}{\epsilon}$).

## 2. Link between Measurement and Particle's Interfacial Properties

When an alternating electric field *E*(ω) is applied to a colloidal suspension, the particles and their double layers become polarized. The charged particles will oscillate in the AC field as will the ions in the double layers, but with two different characteristic times. Such a combination of events usually produces a rich dielectric loss spectrum.

This interesting though complicated relaxation behavior can best be described by the complex dipolar coefficient $\stackrel{~}{\beta}\left(\omega \right)$ which is associated with the (complex) dipole moment of a single colloidal particle surrounded by its double layer. This dipole moment is generated by the polarization of the core material of the colloidal particle, and the anisotropic charge distributions of particle and ions in space and time, due to their respective movements. For a spherical colloidal particle of radius *a* the dipole coefficient $\stackrel{~}{\beta}$ is defined by:

where $\stackrel{~}{P}$ is the particle's dipole moment, ε_{e} the relative dielectric permittivity of the solvent (generally water) and ε_{0} the permittivity of vacuum. Since this quantity incorporates all polarization mechanisms possible, we can gather with dielectric spectroscopy important information from $\stackrel{~}{\beta}$ about the interfacial properties of the particles, in particular the zeta potential and Stern layer conductances of the particle (DeLacey and White, 1981; Delgado, 2002; Ohshima, 2006; Chassagne and Bedeaux, 2008). A general expression of $\stackrel{~}{\beta}$ can be found in Supplementary Material 1, see also Chassagne and Bedeaux (2008), Mangelsdorf and White (1990), and Minor et al. (1998). A discussion about the important characteristic frequencies associated with the system can be found in Grosse and Shilov et al.

The relation between dipolar coefficient and complex conductivity ${\stackrel{~}{K}}_{s}\left(\omega \right)$ of a suspension with low volume fraction ϕ of the dispersed particles is given by Delgado (2002):

where ε_{s}(ω) is the real part of the permittivity of the suspension and *K*_{s}(ω) the real part of its conductivity. A precise derivation of this relation can be found in Grosse. The complex conductivity of the suspending electrolyte is given by

where *K*_{e} is the Ohmic conductivity of the electrolyte solution, defined explicitly later in Equation (28). *K*_{e} equals in principle the conductivity obtained in an ideal DC conductivity experiment. We assume both *K*_{e} and ε_{e} to be frequency-independent. This hypothesis is justified for the frequencies used. In Chassagne and Bedeaux (2008), the notation (*K*_{1}, ε_{1}) was used instead of (*K*_{e}, ε_{e}) to express the conductivity and relative dielectric permittivity of the electrolyte. Instead of the complex conductivity $\stackrel{~}{K}\left(\omega \right)$ for the suspension and the electrolyte we can equally well use the complex permittivity $\stackrel{~}{\epsilon}\left(\omega \right)$ defined by:

It then follows that:

Some authors prefer to overlook the interrelation in Equations (2) and (5) between the components of $\stackrel{~}{K}$ and $\stackrel{~}{\epsilon}$ and simply write:

with the real parts denoted by a single prime and the imaginary parts by a double prime. From which we can deduce that:

Dielectric spectroscopy data are usually plotted in terms of dielectric and conductivity increments, that can be obtained from Equations (2) and (5) and are defined by:

As we are interested in the dipolar coefficient, it is convenient to write Equations (2) and (5) in terms of the real and imaginary parts of $\stackrel{~}{\beta}$:

These two relations enable us to understand an additional source for the errors that can be made in the evaluation of the dipole coefficient. Indeed, at low frequencies, where the magnitude |$\stackrel{~}{\beta}$| of the dipolar coefficient is usually the largest, we obtain:

As the conductivity of the (dilute) suspension is in general very close to the conductivity of the suspending electrolyte, *K*_{s} ≃ *K*_{e}, a small inaccuracy in the measurement of *K*_{s} and/or *K*_{e} will give rise to a large change in *Re*(β) and consequently |$\stackrel{~}{\beta}$|.

In practice ${\stackrel{~}{K}}_{s}$ (or equivently ${\stackrel{~}{\epsilon}}_{s}$) are not measured. The measured impedance of the suspension contains interfering impedances that have to be corrected for. In the introduction we distinguished between two types of effects: effects that originate from the actual electrical connections, leading to stray impedances, and effects coming from the ionic charge distribution in the cell close to the electrodes, leading to Electrode Polarization (EP) impedances. In the present article, only EP impedances will be discussed. This implies that we assume that we have properly compensated for stray impedances and that we have access to the complex impedance ${\stackrel{~}{Z}}_{c,s}$. This impedance ${\stackrel{~}{Z}}_{c,s}$ is the impedance of the cell filled with the suspension. ${\stackrel{~}{Z}}_{c,s}$ contains the EP contribution ${\stackrel{~}{Z}}_{EP},$ but no stray impedances contributions, and the contribution of the bulk of the suspension, i.e., ${\stackrel{~}{Z}}_{s}$ that we are looking for:

The general relation between impedance and conductivity (for any subscript: *x* = *c, s*, *x* = *s*, *x* = *c, e* or *x* = *e*) is given by:

where *C* is a cell constant that depends on the geometry of the electrodes. This relation is only valid in the case that the electrodes are blocking. In case they are not, there would be a remaining DC conductivity at ω = 0. From Equation (12), it follows that for blocking electrodes ${\stackrel{~}{K}}_{c,x}\left(\omega =0\right)={K}_{c,x}\left(\omega =0\right)+0\times i{\epsilon}_{0}{\epsilon}_{c,x}\left(\omega =0\right)=0$ and therefore *K*_{c, x}(ω = 0) = 0.

The cell constant *C* is defined in Supplementary Material 2 and derived in the corresponding sections for planar and cylindrical electrodes. We have:

## 3. Complex Permittivity of an Electrolyte in Presence of Electrode Polarization

In this section, the complex permittivity ${\stackrel{~}{\epsilon}}_{c,e}\left(\omega \right)$ of an electrolyte solution, in the presence of electrode polarization, is analyzed for the planar electrode case. The important hypotheses used for the derivations are:

1 - In the absence of an applied voltage the electrodes are uncharged.

2 - No charge transfer takes place at the electrodes. These “ideally polarizable electrodes” can therefore be considered as capacitor plates.

3 - The electrodes are spaced sufficiently apart so that their respective EP does not influence each other, i.e., |λ_{c}*d*| ≫ 1 (these variables are defined right below), a valid assumption for the experimental conditions encountered in impedance spectroscopy.

The solution is obtained from the set of relevant equations (Nernst-Planck, conservation of ions, Poisson) and appropriate boundary conditions. The mathematical details of the derivations for both planar and cylindrical electrodes can be found in Supplementary Material 2. The electrolyte response to the applied electric field frequency is shown to display the same characteristic features, besides geometrical aspects, in both cases.

We here recall the expressions of important variables defined in Supplementary Material 2, which are needed to evaluate the complex permittivities given below:

with

where:

ν* _{i}* is the stoichiometric coefficient of ion

*i*with

*i*= +, −

*z _{i}* is the valence of ion

*i*

*D _{i}* is the diffusion constant of ion

*i*

ω is the radial frequency of the applied electric field

κ^{−1} is the Debye length

Furthermore *d* is the distance between electrodes.

The diffusion coefficients of the ions can be obtained from their limiting conductivities ${\Lambda}_{i}^{\infty}$ (which can be found in Handbooks) from the relation:

where *N*_{a} is Avogadro's number, *e* the electric elementary charge, *k* Boltzmann's constant and *T* the temperature. The Debye length can be obtained from the relation:

where the ionic density is defined by

where *C*_{s} is the salt concentration in 10^{−3} mol/L (i.e., millimolar, mM).

In Supplementary Material 2 it is shown that the measured complex permittivity of a binary electrolyte solution, in the case of planar electrodes is given by

for all frequencies, valences and ionic strengths. One can verify that for frequencies such that $\omega \ll {D}_{\pm}{\kappa}^{2}$, this expression corresponds to the expression found in Chassagne et al. (2002), which is also valid for all types of electrolytes but for frequencies restricted to $\omega \ll {D}_{\pm}{\kappa}^{2}.$ An illustration is given in Figure 2. The condition $\omega \ll {D}_{\pm}{\kappa}^{2}$ corresponds to the one encountered mostly in experiments. For these frequencies, one also finds that

The lengthscale ${\lambda}_{c}^{-1}$ is related to the creation of the double layer close to the electrodes: for low frequencies $\left(\omega \ll {D}_{\pm}{\kappa}^{2}\right)$ one gets ${\lambda}_{c}^{-1}\simeq {\kappa}^{-1}$ and for high frequencies ${\lambda}_{c}^{-1}\simeq 0$. The lengthscale ${\lambda}_{n}^{-1}$ is related to the ionic diffusion.

**Figure 2. (A)** Conductivity K(S/m) and **(B)** relative permittivity (epsilon) of a solution of divalent salt as function of the applied electric field frequency, similar to MgCl2, for which *D*_{1} = 1.4 × 10^{−9} m^{2}/s and *D*_{2} = 2.0 × 10^{−9} m^{2} /s. The salt concentration is 0.5 mM. The electrodes are planar. The important characteristic frequencies associated with the system are ${\omega}_{0}={D}_{0}{\kappa}^{2}$ and ω_{P} = 2κ*D*_{0}/*d* where *d* = 10 × 10^{−3} m is the distance between the electrodes. The full solution found in this study, i.e., Equation (23) (red full line) is compared to the solution (blue dashed line) found in Chassagne et al. (2002), which was derived for the case ω ≪ ω_{0}. As expected, the two solutions overlap for ω ≪ ω_{0}. The green dotted line in the conductivity plot represents the theoretical conductivity, the value of which is given by Equation (29). The green dotted line in the epsilon plot represents the relative permittivity of water.

In the particular case where *D*_{−} = *D*_{+} = *D* which has been treated by several authors, see Cirkel et al. (1997), Buck (1969), and Hollingsworth and Saville (2003) one finds the simplifications:

Equation (25) is for example the expression found by Cirkel et al. [their Equation (8) in Cirkel et al., 1997]. Explicit expressions for the real and imaginary part of ${\stackrel{~}{\epsilon}}_{c,e}$ are given in Supplementary Material 4. There, we also make a comparison with the work done by Kang et al. (Kang and Dhont, 2010) on the in and out phase component of the alternating electric field. Since we have made the hypothesis that λ_{c}*d* ≫ 1 we can safely assume that *tanh*(λ_{c}*d*/2) ≃ 1 and Equation (25) is also in agreement with Equation (25) of Hollingsworth et al. in (Hollingsworth and Saville, 2003), who restricted their calculations to 1:1 electrolytes (for which *z*_{+} = −*z*_{−} = 1).

Two important relaxations frequencies are found from the analysis of Equation (23):

The frequency ω_{P} corresponds to the frequency below which charges can fully build-up a double layer close to the blocking electrodes due to the application of the alternating electric field. Below ω_{P} and beyond the double layers, the bulk electric field is consequently zero. For ω ≪ ω_{P} we get:

The conductivity *K*_{c, e} is in good approximation zero for low frequencies because the electrodes are blocking. No charge transfer is possible and the ions are all accumulated close to the electrodes.

The Maxwell-Wagner frequency ω_{0} represents the frequency below which charges can be dissipated in the system. For ω_{P} ≪ ω ≪ ω_{0} we get:

As pointed out in recent discussions (Grosse and Delgado, 2013; Hollingsworth, 2013), we note that ε_{c, e} indeed scales as ω^{−2} in this frequency range. The conductivity in this frequency range is equal to the theoretical conductivity *K*_{e} of the electrolyte solution. We didn't include ionic interactions, which would have modified the value of *K*_{e} presented here. This implies in particular that we consider electrolytes at low ionic strength, for which the conductivity varies linearly with the ionic strength (${K}_{e}\propto {\kappa}^{2}\propto {C}_{s}$ where *C*_{s} is the salt concentration). Above ω_{0} no double layer can be established at the electrodes and the electric field in the cell is everywhere equal to the applied one. For ω ≫ ω_{0} we get:

The permittivity of the electrolyte solution has reduced to the one of water (we do not consider the frequency-dependence of ε_{e}) and the conductivity is equal to the usual conductivity of the electrolyte solution since we have not considered inertial effects that might slow down the ions at high frequencies. The relaxation frequency ω_{0} is defined by the frequency above which ε_{0}ε_{e}ω > *K*_{e} i.e., ${\omega}_{0}={K}_{e}/\left({\epsilon}_{0}{\epsilon}_{e}\right)={D}_{0}{\kappa}^{2}.$ It can be verified that for $\omega \gg {D}_{\pm}{\kappa}^{2}$ the permittivity (ε_{c,e} − ε_{e}) scales as ω^{−3/2} (Cirkel et al., 1997; Hollingsworth, 2013). However, this dependence can not be observed. If one takes an electrolyte such that *D*_{−} = *D*_{+} = *D* the expansion of Equation (25) reads:

Even under extremely good conditions such that, for instance ω ≃ κ^{2}*D*, for an electrode spacing of 1 mm and a very low ionic strength of 0.1 mM NaCl, one finds ${\epsilon}_{c,e}\simeq {\epsilon}_{e}\left[1+1{0}^{-5}\right]$ and the correcting term is virtually undetectable.

Indicative values of ω_{P} and ω_{0} are given in Table 1.

**Table 1. Values of ω _{P} and ω_{0} as given by Equation (26), using a diffusion coefficient of ${D}_{0}=2.1{0}^{-9}$ m^{2}/s at room temperature, for 3 salt concentrations**.

## 4. Equivalent Circuits

Traditionally Dielectric Spectroscopy measurements have been described by equivalent electrical circuits (Barsoukov and Macdonald, 2005). A comparison between the results found from equivalent circuits and those found from solving the electrokinetic set of equations has been performed in specific cases (Buck, 1969; Barsoukov and Macdonald, 2005). We are going to show that the description of the system (any electrolyte solution or colloidal suspension) using equivalent circuits is strictly equivalent to the results found in the previous section, for all frequencies and cell geometries (planar and cylindrical). We will start with the description of electrolyte solutions. The results for colloids will be given in the following section. The solution can be described by a resistance *R*_{b} in parallel with a capacitor *C*_{b}. The electrode polarization can be seen as an additional capacitor *C*_{EP} in series with the equivalent circuit of the solution, see Figure 3.

**Figure 3. Equivalent circuit representation for a cell containing a solution or suspension, with blocking electrodes**.

The total impedance of the system ${\stackrel{~}{Z}}_{c}$ (including electrode polarization) can be evaluated from the relations:

The general expression for ${\stackrel{~}{Z}}_{c}$ is given in Supplementary Material 2. From that expression, one can verify that for ω_{0} ≫ ω ≫ ω_{P} one gets

Schwan (1992) found this result experimentally in the case of a concentrated suspension of blood cells. Like many experimentalists, he expresses his complex measured impedance in terms of a resistance and capacitance in series i.e.,

and defines:

In Figure 5 of Schwan (1992), it is shown that $C\simeq 1/\left({\omega}^{2}{R}_{b}^{2}{C}_{EP}\right)$ for low frequencies and that *C* ≃ *C*_{b} at high frequencies. As EP only contributes as a capacitance, it implies that EP only affects *C* and not *R*. In other words, in the frequency range of interest (ω_{0} ≫ ω ≫ ω_{P}) only the permittivity is affected by EP, not the conductivity, which reduces to the bulk conductivity. This is consistent with the results of Section 3.

### 4.1. Electrolyte Solutions

In Supplementary Material 3 we describe how we obtain the expressions for *C*_{EP}, *C*_{b}, and *R*_{b} from the mathematical equivalence with the analytical results of Section 3. In the case of planar electrodes, we find:

These values were also found by Buck (1969) in the case of a z-z electrolyte. One can verify that the characteristic frequencies found in Supplementary Material 3 can be related to the ones found in the Planar electrodes section by:

The signification of the frequencies (specified in Supplementary Material 2) can now be linked to circuit elements: as expected ω_{P} depends on the capacitance *C*_{EP} associated to electrode polarization, whereas ω_{0} depends on *C*_{b} associated to the bulk permittivity. Note that the “2” that appears in the right hand side part of the relation for ω_{P} can be seen as resulting from the fact that there are two electrodes in the system: each electrode contributes with a capacitance *Sε*_{0}ε_{e}κ and summing two capacitances in series yields *C*_{EP} = *Sε*_{0}ε_{e}κ/2 as in Equation (35). Some authors therefore prefer to use an alternative equivalent circuit : instead of the one represented in Figure 3, where one capacitor of capacitance *C*_{EP} = *Sε*_{0}ε_{e}κ/2 was used, one can have a circuit with two capacitors, each of capacitance *C*_{EP} = *Sε*_{0}ε_{e}κ, placed on each side of the ${\stackrel{~}{Z}}_{s}$ element, as was sketched in Figure 1.

Following the same procedure for the cylindrical electrodes (not detailed here), one finds the equivalence, using the circuit represented in Figure 3:

The analytical expressions of Section 3 and the equivalent circuit expressions give the same results for all frequencies, ionic strengths, and type of salts (not shown). Buck (1969) also derived an equivalent circuit in the case of spherical concentric electrodes for z-z electrolytes. We refer to him for this geometry which will not be further discussed, as most experiments nowadays are performed with either planar or cylindrical electrodes.

### 4.2. Colloidal Suspensions

For the sake of argument, we will consider planar electrodes. Instead of a simple electrolyte solution, as studied in the previous section, we now like to find an equivalent circuit in the case a complex electrolyte solution such as a colloidal suspension. We make the hypothesis that the equivalent circuit can be represented as the one given in Figure 3: we would like that when the concentration of colloidal particles goes to zero, one would find again Equation (35). The equivalent circuit elements are taken to be:

Note that in the general case the circuit elements now become frequency-dependent: even though the *R* and *C* in Equation (37) have the dimensions of resistance and capacitance, they represent no “real” resistor and capacitors, and no “real” electronic equivalent circuit can be made of them. This was not the case for electrolyte solutions, where all *R* and *C* were frequency-independent, and which could therefore be substituted by “real” resistor and capacitors. The equivalent circuit would then reproduce the behavior of the electrolyte between blocking electrodes. The expressions for ε_{s} and *K*_{s} are given in Equation (8). An important check for the validity of the proposed equivalent circuit is that in the absence of electrode polarization, the analytical expression for ${\stackrel{~}{K}}_{s}$ and the equivalent circuit expression should give the same results: this is illustrated in Figure 4. This is not surprising, as one can see by comparing Equations (12) and (31).

**Figure 4. (A)** Conductivity K(S/m) and **(B)** relative permittivity (epsilon) as function of frequency of a suspension of 100 nm colloidal spheres (ϕ = 1%, *eζ*/*kT* = 4) in a 1 mM electrolyte solution of monovalent salt solution for which ${D}_{1}=2\times 1{0}^{-9}$ m^{2}/s and ${D}_{2}=3\times 1{0}^{-9}$ m^{2}/s. Red curve: the suspension in the absence of electrode polarization, corresponding to ε_{s} and *K*_{s} from Equation (8). Dashed blue curve: the equivalent circuit model corresponds to the theoretical prediction provided that one takes *R*_{b} = 1/*K*_{s}, *C*_{b} = ε_{0}ε_{s} and *C*_{EP} = 0.

If there are no colloidal particles, we get from Equation (5) that ${\stackrel{~}{\epsilon}}_{s}\left(\omega \right)={\stackrel{~}{\epsilon}}_{e}\left(\omega \right)$ and we indeed find again Equation (35), provided that the unknown ε* equals ε_{e}. The permittivity ε* is clearly linked to the EP phenomena by construction. One could argue that EP originates mainly from the contribution of the ions and not that of the colloidal particles. This implies that the dielectric permittivity ε* in the small slab of liquid adjacent to the electrodes should be equal to ε_{e} whether in the presence or not of colloidal particles. If the colloidal particles contribute to ε*, ε* should be close to ε_{s}(ω). The relation between ε_{s} and ε_{e} is given by Equation (8), and these permittivities do not differ much, especially not at low volume fractions. However, even a small difference in permittivities will prevent that the subtraction method is applicable when ${\epsilon}^{*}\ne {\epsilon}_{e}$ (See section 5). We note that both *C*_{b} and *R*_{b} depend on *d*, whereas *C*_{EP} does not. This fact will allow to account for EP by the variable separation method (See section 5).

In Figure 4, we have shown the case where *C*_{EP} = 0 (no electrode polarization). Two relaxations features can be observed that are related to the characteristic frequencies associated with the polarization of the particle and its double layer. The double layer around a colloidal particle and the double layer at the electrodes have the same relaxation frequency ω_{0} and the frequency associated with the relaxation of the particle is given by ${\omega}_{a}={D}_{0}/{a}^{2}$ where *a* is the radius of the colloidal particle and *D*_{0} is given in Equation (14). If we take electrode polarization into account, the dielectric permittivity at low frequency is 10^{6} times larger than that without electrode polarization. This is illustrated in Figure 5. In this example we have used ${\epsilon}^{*}={\epsilon}_{s}$.

**Figure 5. (A)** Conductivity K(S/m) and **(B)**: relative permittivity (epsilon) as function of frequency. Suspension with same properties as the one given in Figure 4. The suspending electrolyte solution (blue line) and suspension (red line) in the presence of electrode polarization. The magenta dashed line corresponds to the solution found in Figure 4 for no EP. The equivalent circuit of the suspension in the presence of electrode polarization was constructed by taking *R*_{b} = 1/*K*_{s}, *C*_{b} = ε_{0}ε_{s} and *C*_{EP} = ε_{0}ε_{s}κ*d*/2 with *d* = 10 mm.

In the frequency range of interest, i.e., for ω_{0} ≫ ω ≫ ω_{P} Equation (32) holds. A new frequency can be defined by

This frequency is discussed in Supplementary Material 4. One can verify from Equation (32) that for ω ≫ ω_{b} the influence of electrode polarization can be neglected and

Equivalently,

from which $\stackrel{~}{\beta}$ can be directly deduced. This is illustrated in Figure 6 which is an enlargement of Figure 5.

**Figure 6. Enlargement of Figure 5**. Above ω_{b} one finds that ε_{c,s} ≃ ε_{s} i.e., that the EP plays no role anymore. Note that for any frequency above ω_{P} the relation *K*_{c,s} ≃ *K*_{s} holds (see Figure 5). This implies that for frequencies above ω_{b} one has ${\stackrel{~}{K}}_{c,s}\simeq {\stackrel{~}{K}}_{s}$ (or equivalently ${\stackrel{~}{\epsilon}}_{c,s}\simeq {\stackrel{~}{\epsilon}}_{s}$). Similarly, above ω_{b} one has ${\stackrel{~}{K}}_{c,e}\simeq {\stackrel{~}{K}}_{e}$ which implies in particular for the present figure that ε_{c,s} ≃ ε_{e}.

The characteristic frequency ω_{b}, representing the frequency above which EP becomes negligible, increases with ionic strength in the following way:

An indication for the values of ω_{b} is given in Table 2.

**Table 2. Values of ω _{b} as given by Equation (41), using a typical diffusion coefficient of ${D}_{0}=2.1{0}^{-9}$ m^{2}/s at room temperature, for 3 salt concentrations**.

Interestingly, one can verify that in most cases encountered in practice, one has ω_{b} ≲ (ω_{a}, ω_{0}). We will come back to this point in the following section.

As above ω_{b} one has ${\stackrel{~}{K}}_{c,s}\simeq {\stackrel{~}{K}}_{s}$ this also implies that the value of ε* only plays a role below ω_{b}. This is illustrated in Figure 7. In order to better distinguish the features of the curves, we have chosen a large volume fraction (20%). This is permitted for the theoretical considerations presented, as using Equation (2) one can verify that $|\varphi \stackrel{~}{\beta |}<1$.

**Figure 7. (A)** Conductivity *K*_{c,s}(S/m) and **(B)** relative permittivity (epsilon) ε_{c,s} of a suspension of 100 nm colloidal spheres (ϕ = 20%, *eζ*/*kT* = 4) in a 0.23 mM electrolyte solution of monovalent salt solution for which ${D}_{1}=2\times 1{0}^{-9}$ m^{2}/s and ${D}_{2}=1.98\times 1{0}^{-9}$ m^{2}/s. The spacing between electrodes is 10 mm. Red curve: the case where ε* = ε_{e}; Dashed magenta curve: the case where ε* = ε_{s}.

## 5. Accounting for Electrode Polarization

In this section, first three methods used to compensate experimentally for electrode polarization will be analyzed. In the last section, a general method to extract $\stackrel{~}{\beta}$ directly from the data will be presented. The frequency range of validity of each method is discussed.

### 5.1. Subtraction Method

This method relies on the fact that EP can be modeled by a capacitance *C*_{EP} (see Section 3) associated with a relative permittivity ${\epsilon}^{*}={\epsilon}_{e}$. The aim is to find a reference electrolyte solution (symbolized by the subscript *r*) that will give as best as possible the same *C*_{EP} as the one of the investigated suspension. From the section about equivalent circuits and from general definitions [see Equation (12) and Supplementary Material 1], we can write that (for *x* = *r* or *x* = *s*):

For the sake of argument, we will consider planar electrodes and use *C* = *S*/*d*. We then have:

The ionic strength is determined by the amount of charge carriers in the system, which implies that conductivity and ionic strength are linked, see Equation (28) for electrolyte solutions for example. For the reference electrolyte, it is therefore found that ${K}_{r}\propto {\kappa}_{r}^{2}$. On the other hand, the conductivity for colloidal suspensions is given by Equation (8) i.e., at low frequencies:

Equation (44) gives the relation for the conductivity of the suspension after it has been mixed with the suspending electrolyte. The conductivity ${K}_{e}={\epsilon}_{0}{\epsilon}_{e}{\kappa}_{s}^{2}{D}_{0}$ therefore represents the conductivity of the suspending electrolyte *including* all other ions (counterions and impurities) originating from the original (concentrated) suspension used to be mixed with the suspending electrolyte. In some cases, one can assume that the original suspension is devoid of impurities, and that the amount of counterions is negligeable compared to the amount of ions stemming from the electrolyte suspension. This last condition is usually true when dilute suspensions of large (> 10 nm) colloidal particles are studied. One can then safely assume that κ_{s} = κ_{e} i.e., the ionic strength of the suspension (symbolized by the subscript *s*) is given by the ionic strength of the suspending electrolyte (symbolized by the subscript *e*). In the case of nanocolloids, i.e., when there is a relative large total interface between particles and water and hence a non-negligible amount of counterions, or in the case that the original suspension of colloids (before it is added to the suspending electrolyte) contains impurities, i.e., ions that are not counterions, one can have κ_{s} ≠ κ_{e}.

In practice, the value of κ_{s} can be difficult to obtain, because measuring the conductivity of a suspension does not give direct access to κ_{s}. A procedure, based on experimental results, has consequently been developed to get a reference electrolyte solution for which in close approximation *C*_{EP,s} ≃ *C*_{EP,r} without the need of κ_{s} being known. This procedure is now outlined.

In the frequency range of interest, i.e., for ω_{0} ≫ ω ≫ ω_{P} we get from Equations (32), (35), and (37) for the suspension and the reference solution:

As discussed in the previous section, the conductivity in this frequency range is not affected by EP, only the permittivity is, in the range [ω_{P}, ω_{b}]. Above ω_{b} one has ε_{c,s}(ω) ≃ ε_{s}(ω). The subtraction method is therefore only applied to the permittivities, and should compensate for EP in the range [ω_{P}, ω_{b}]. The reference electrolyte solution is made of an electrolyte of the same type as the electrolyte present in the suspension (ex: if NaCl is the suspending electrolyte, NaCl will be taken as reference electrolyte). The concentration of the reference electrolyte is chosen such that the conductivity of the reference electrolyte *K*_{r} equals:

which implies that:

In the range ω_{P} < ω < ω_{b} the conductivity *K*_{s}(ω_{P} < ω < ω_{b}) is constant. For frequencies smaller than ω_{P} the conductivity *K*_{c,s} is going to zero, as the electrodes are blocking, but above ω_{P} one has *K*_{c,s}(ω) = *K*_{s}(ω), as was shown in Section 3. One can verify from Equation (44) that below ω_{b} one has *K*_{s}(ω) = *K*_{s} (i.e., *K*_{s} does not depend on frequency) provided that $\text{Re}\left(\stackrel{~}{\beta}\right)$ remains constant.$\text{}\text{Re}\left(\stackrel{~}{\beta}\right)$ is indeed a constant at low frequencies, and will start to exhibit relaxation phenomena when the lowest characteristic frequency associated with the colloidal particle is reached. This frequency is given by

where *a* is the radius of the colloidal particle and *D*_{0} is given in Equation (14). This relaxation frequency occurs in the frequency regime [10^{4}–10^{6}] rad/s for particle sizes between [25–250] nm. In most cases encountered in experiments one has therefore ω_{a} > ω_{b} and $\text{Re}\left(\stackrel{~}{\beta}\right)$ is constant below ω_{b} implying that *K*_{s} is constant. For extreme cases, like very large particles at low ionic strength, the situation ω_{a} ≃ ω_{b} could occur in which case an option can be to increase *d* so as to lower ω_{P} and take *K*_{r} = *K*_{s}(ω ≳ ω_{P}).

As ϕ is supposed to be small, one has in the range ω_{P} < ω < ω_{b}, in good approximation, ${K}_{s}\simeq {K}_{e}={\epsilon}_{0}{\epsilon}_{e}{\kappa}_{s}^{2}{D}_{0}$ which implies that κ_{s} ≃ κ_{r}. One can now get rid of the EP contribution in the whole frequency range ω_{0} ≫ ω ≫ ω_{P}. By substracting ε_{c, r} from ε_{c,s} one obtains for ε_{s}(ω):

As stated at the beginning of this section, the subtraction method will not work when ${\epsilon}^{*}\ne {\epsilon}_{e}$ as this would imply that *C*_{EP,s} ≠ *C*_{EP,e}. If one assumes that ${\epsilon}^{*}={\epsilon}_{s}$ even though ε_{s} and ε_{e} differ by less than 10% this will give huge errors in the subtraction method, as Figure 8 illustrates. Note that around ω_{p} the values for the permittivity are of the order of 10^{7} (as can be seen in Figure 5). Substracting such huge numbers goes at the cost of accuracy: this explains the deviations observed the theoretical prediction (red line) and substraction method result (dashed blue line) at low frequencies.

**Figure 8. Relative permittivity (epsilon) of a suspension of 100 nm colloidal spheres (ϕ = 20%, eζ/kT = 4) in a 0.23 mM electrolyte solution of monovalent salt solution for which ${D}_{1}=2\times 1{0}^{-9}$ m^{2}/s and ${D}_{2}=1.98\times 1{0}^{-9}$ m^{2}/s**. The spacing between electrodes is 10 mm.

**(A)**The case where ε* = ε

_{e}.

**(B)**The case where ε* = ε

_{s}Red curves: ε

_{s}Blue curves: ε

_{s}found by using the subtraction procedure

As discussed after Equation (45), two main frequency regions can be distinguished in the spectrum. In the range [ω_{P}, ω_{b}] electrode polarization is dominant and we obtain from Equation (45):

where ${\epsilon}_{x}={\epsilon}^{*}$ for suspensions and ε_{x} = ε_{e} for electrolytes. One might wonder whether it is possible to compensate for electrode polarization by dividing the dielectric spectrum of the suspension with the one of the solvent. This gives in the range [ω_{P}, ω_{b}]:

This ratio will probably be close to one, considering the discussion given above about the values of these variables. Clearly, this way of elimination of the electrode polarization contribution does not give access to ε_{s} in the frequency range where electrode polarization is prevailing. Above ω_{b} we have shown that electrode polarization does not play a significant role, and therefore we obtain from Equation (45):

Combining real and imaginary terms we obtain from Equations (12, 45) for the complex ratio, above ω_{b}:

From Equations (2) and (5) one can see that above ω_{b} the ratios $\text{Re}\left({\stackrel{~}{K}}_{c,s}/{\stackrel{~}{K}}_{c,e}\right),$ or alternatively $\text{Re}\left({\stackrel{~}{\epsilon}}_{c,s}/{\stackrel{~}{\epsilon}}_{c,e}\right),$ provide direct information about $\text{Re}\left(\stackrel{~}{\beta}\right)$ whereas $\text{Im}\left({\stackrel{~}{K}}_{c,s}/{\stackrel{~}{K}}_{c,e}\right),$ or $\text{Im}\left({\stackrel{~}{\epsilon}}_{c,s}/{\stackrel{~}{\epsilon}}_{c,e}\right),$ provide direct information about $\text{Im}\left(\stackrel{~}{\beta}\right).$ By contrast, the ratios $\text{Re}\left({\stackrel{~}{\epsilon}}_{c,s}\right)/\text{Re}\left({\stackrel{~}{\epsilon}}_{c,e}\right)$ and $\text{Im}\left({\stackrel{~}{\epsilon}}_{c,s}\right)/\text{Im}\left({\stackrel{~}{\epsilon}}_{c,e}\right)$ are less appealing because they each depend on a mix of $\text{Re}\left(\stackrel{~}{\beta}\right)$ and $\text{Im}\left(\stackrel{~}{\beta}\right)$.

### 5.2. Variable Electrode Separation Method

This technique is, in principle, extremely well-suited for compensating for electrode polarization. We start again from the general definitions:

For the sake of argument, we again consider planar electrodes and hence *C* = *S*/*d*. Moreover:

Combining these equations leads to:

From the experimental data one can then estimate ${\stackrel{~}{K}}_{s}^{-1}$ (and hence ${\stackrel{~}{K}}_{s}$) from a linear regression of ${\stackrel{~}{K}}_{c,s}^{-1}$ as function of 1/*d*. This method is applicable for the whole range of frequency, however it requires the combination of both the real and imaginary part of ${\stackrel{~}{K}}_{c,s}$ in the analysis. In the frequency range of interest, i.e., for ω_{0} ≫ ω ≫ ω_{P}, the variable electrode separation technique can be applied to the imaginary part of ${\stackrel{~}{K}}_{c,s}$ (which is related to ε_{c,s}) only. Indeed, we have found that in the range of interest [see Equation (45)]:

It is therefore possible to “clean” the signal from electrode polarization by a linear regression of ε_{c,s} as function of 1/*d*, from which ε_{s} is then easily obtained. Although 2 electrode separations suffice in theory, in practice, 3 or 4 electrode separations are usually used to minimize the experimental error. The technique was introduced by Fricke and Curtis (1937), and used, for suspensions, in combination with the subtraction method (described in the previous section) by the Dutch groups (van der Touw and Mandel, 1971; van der Touw et al., 1975; Cirkel et al., 1997; Chassagne et al., 2002, 2003). This technique is also used by the groups in Princeton and Granada (Hollingsworth and Saville, 2003; Jimenez et al., 2007). For more recent work on the topic, and experimental limitations of the techniques, we refer to Hollingsworth and Saville (2004), Beltramo and Furst (2012), and Beltramo and Furst (2012).

### 5.3. Four Electrode Cell Method

The idea of using 4 electrodes stems from Schwan (1992). Four electrode cells are designed such that the two inner electrodes are able to measure a voltage difference that is devoid of EP. The outer electrodes are the current carrying electrodes, where EP occurs and causes an extra voltage drop close to these electrodes. Figure 9 sketches the measuring cell with 4 electrodes and equipped with special electronics. For the voltage sensing often needle-like electrodes are used so as to minimize their influence on the ionic flow. The inner electrodes should be designed such that the electric current generated by the outer electrodes is not disturbed. This requires that the voltage difference measured at the inner electrodes ${\stackrel{~}{V}}_{i}$ should be at virtually zero current (${\stackrel{~}{I}}_{i}\simeq 0$), so that no current leaks away in the probing circuit. This implies that the impedance ${\stackrel{~}{Z}}_{i}$ of the inner electrodes should be virtually infinite (since ${\stackrel{~}{V}}_{i}={\stackrel{~}{Z}}_{i}{\stackrel{~}{I}}_{i}$). The measurement of ${\stackrel{~}{V}}_{i}$ can be achieved with the use of operational amplifiers. The current $\stackrel{~}{I}$ that is flowing in the whole cell can be measured with commercial frequency response analyzers or with impedance analyzers. We have:

where ${\stackrel{~}{Z}}_{c}$ is the impedance of the cell (with electrode polarization) and ${\stackrel{~}{Z}}_{ci}$ the impedance of the fluid between the inner electrodes. Since EP should not play a role at the inner measuring electrodes, we have ${\stackrel{~}{Z}}_{ci}={\stackrel{~}{Z}}_{s}$ i.e., the measured impedance ${\stackrel{~}{Z}}_{ci}$ gives directly the sought conductivity (or permittivity) of the suspension:

Note that the cell constant *C* can be found separately from the impedance of an electrolyte with known electrical properties as the correct geometry of the probing electrodes might be unknown. We then simply have $C={\stackrel{~}{K}}_{e}\times {\stackrel{~}{Z}}_{ci}$, where ${\stackrel{~}{K}}_{e}={K}_{e}+i\omega {\epsilon}_{0}{\epsilon}_{e}$ is known, and ${\stackrel{~}{Z}}_{ci}$ is measured.

**Figure 9. Schematic representation of a 4-electrode cell; the voltage difference at the outer electrodes, where EP takes place, is ${\stackrel{~}{V}}_{o}=\Delta {\stackrel{~}{V}}_{out}={\stackrel{~}{Z}}_{o}\stackrel{~}{I}$ whereas the voltage difference at the inner (probing) electrodes is ${\stackrel{~}{V}}_{i}=\Delta {\stackrel{~}{V}}_{in}$; as EP should be minimized at the inner electrodes, the voltage difference measured at the inner electrodes ${\stackrel{~}{V}}_{i}$ should be done with virtually zero inner current (${\stackrel{~}{I}}_{i}\simeq 0$), implying that the impedance ${\stackrel{~}{Z}}_{i}$ of the inner electrodes should be virtually infinite**. One then obtains, from the measurement of ${\stackrel{~}{V}}_{i}$ and $\stackrel{~}{I}$ : ${\stackrel{~}{Z}}_{s}={\stackrel{~}{V}}_{i}/\stackrel{~}{I}.$

Clearly, both the cell design and the electronic are more sophisticated than those of the common 2 electrode cell. Further details are given in Hayakawa et al. (1975) and Pelc et al. (2011). In the early publications gain/phase detectors were not available, so they were home built. Also, high technical skills in order to properly compensate for parasitic stray capacitances and inductions were required (Hayakawa et al., 1975; van der Touw et al., 1975). A new technique van der Ploeg and Mandel (1991) was proposed in 1991 to reduce the measurement time, and hence the possible drifts in conductivity. This technique is no longer in use as current set-ups are improved in such a way that the measurement time is not anymore an issue. Nowadays setting up 4 electrode experiments has become easier (Pelc et al., 2011). Four electrode cells are now also commercially available from e.g., Novocontrol.

It was already Schwan (1992) who warned that in practice ${\stackrel{~}{Z}}_{ci}$ could still contain some contribution of EP. This arises from the parasitic capacitance between the 4 electrodes. This extra parasitic EP can be compensated for by using the subtraction method explained above. This has been done by authors like Mandel, Saville, Hayakawa, and Kijlstra (Hayakawa et al., 1975; van der Touw et al., 1975; Myers and Saville, 1989; Kijlstra and Wegh, 1994). Critical remarks on this issue have also been made recently by Grimnes et al., Mazzeo et al., and Lvovich (Mazzeo, 2009; Lvovich, 2012; Grimnes and Martinsen, 2015). Another concern is that the 4 electrode system and its additional electronics is less suitable for high frequencies.

### 5.4. Logarithmic Derivative Method

This method was introduced by Jimenez et al. (2002, 2007) in order to compensate for EP, following the work of van Turnhout et al. (Wübbenhorst and van Turnhout, 2002). The logarithmic derivative method was originally derived to study the dielectric relaxation in polymers, but Jimenez et al. showed that it can be applied to colloidal suspensions as well, see also Kaatze et al. (Kaatze and Feldman, 2006).

These authors use the following definitions for the complex permittivity/conductivity of the system, which are different from Equations (6) :

implying in particular that:

As discussed earlier in the present article (see Section 2), the conductivity ${\stackrel{~}{K}}_{c}$ is zero at ω = 0 in the case of blocking electrodes; this implies that in our case *K*_{DC} = 0. We therefore obtain the equivalence, using Equations (6), (12), and (60):

The logarithmic method is based on the use of a new variable defined by Wübbenhorst and van Turnhout (2002):

where the subscript “*D*” refers to “derivative.” This new variable ${\epsilon}_{D}^{\prime \prime}\left(\omega \right)$ enables to better distinguish the characteristics relaxation frequencies of the system.

Despite its notation, ${\epsilon}_{D}^{\prime \prime}$ differs from ε″. For a reference electrolyte solution of conductivity *K*_{r} one has, in the frequency range of interest:

This implies that:

This result is in agreement with the discussion in Wübbenhorst and van Turnhout (2002) about the expected quadratic dependence of ${\epsilon}_{D}^{\prime \prime}$ on ε″ for sharp Debye like loss peaks. This quadratic dependence implies that ∂ε′/∂ ln ω provides a higher spectral resolution than ε″. It makes this quantity better suited for the decomposition of overlapping relaxations.

Let us say a few words about the origin of the use of ∂ε′/∂ ln ω (Steeman and van Turnhout, 1994). It was ∂ε′/∂ ln ω that enabled the authors of Steeman and van Turnhout (1994) to detect in its very spectrum for solid polymers just beyond the main-chain or α-relaxation, the space charge or ρ-relaxation, which arises from the motion of ions at high temperature and low frequencies. By contrast, the ρ-relaxation could not be seen at all in the ε″-spectra. Its better revealing power is furthermore due to another special feature of ∂ε′/∂ ln ω, namely that it does not contain any contribution of ohmic conduction. The conduction loss does add up however in the ε″ data, in which it often overshadows the genuine relaxation losses at low frequencies.

If the relaxation of a system stems from a variety of processes, then we better model the total response with a logarithmic distribution of relaxation times. Both ∂ε′/∂ ln ω and ε″ can be linked to this distribution (Steeman and van Turnhout, 1994). In fact, they provide approximations to it, −∂ε′/∂ ln ω gives a first order and ε″ a zero order estimate. This explains once more that ∂ε′/∂ ln ω has a higher resolving power than ε″.

Both approximations become the better, the broader the distribution. For broad dielectric distributions we can therefore envisage −∂ε′/∂ ln ω as being a close approximation to π/2ε″. It was for this reason that the symbol ${\epsilon}_{D}^{\prime \prime}$ was used for (−π/2)∂ε′/∂ ln ω.

In Figure 10 we have considered a colloidal suspension. One can see that indeed $\left({\epsilon}_{D,c,s}^{\prime \prime},{\epsilon}_{D,s}^{\prime \prime}\right)$ display clear relaxation phenomena, whereas these dispersions cannot be seen in the variables $\left({\epsilon}_{s}^{\prime},{\epsilon}_{s}^{\prime \prime}\right)$ and $\left({\epsilon}_{c,s}^{\prime},{\epsilon}_{c,s}^{\prime \prime}\right)$.

**Figure 10. Comparison between ${\epsilon}_{D,s}^{\prime \prime}$, ${\epsilon}_{D,c,s}^{\prime \prime}$, $\left({\epsilon}_{s}^{\prime},{\epsilon}_{s}^{\prime \prime}\right)$ and $\left({\epsilon}_{c,s}^{\prime},{\epsilon}_{c,s}^{\prime \prime}\right)$ as function of frequency; suspension of 100 nm colloidal spheres (ϕ = 1%, eζ/kT = 4) in a 1 mM electrolyte solution of monovalent salt solution for which ${D}_{1}=2\times 1{0}^{-9}$ m^{2}/s and ${D}_{2}=1.98\times 1{0}^{-9}$ m^{2}/s**. The spacing between the electrodes is 10 mm. Using the logarithmic derivative method (in blue) enables to better distinguish the relaxation processes associated with the colloidal particles.

This high-resolution property of ${\epsilon}_{D}^{\prime \prime}$ therefore enables us to identify the two important relaxations frequencies occurring in a colloidal suspension, namely ${\omega}_{a}={D}_{0}/{a}^{2}$ and ω_{0} [see Equations (26) and (48)] where *a* is the radius of the colloidal particle and *D*_{0} is given in Equation (14). The relaxation frequency ω_{a} shows up in the frequency regime [10^{4}–10^{6}] rad/s for particle sizes between [25–250] nm. The second relaxation frequency, ${\omega}_{0}={D}_{0}{\kappa}^{2},$ is associated with the colloidal particle double layer. The characteristic length scale related to this frequency is the Debye length κ^{−1}. Because the electrolyte that determines the conductivity of the bulk electrolyte is the same as the one in which the colloidal particles are suspended, ω_{0} is a characteristic frequency both for the electrolyte solution (as in this case ${\omega}_{0}={K}_{e}/({\epsilon}_{0}{\epsilon}_{e})={D}_{0}{\kappa}^{2}$) and for the colloidal particles. We have seen in Table 1that in general ω_{0} ≥10^{6} rad/s. Note that for nanoparticles the two characteristic frequencies might overlap: ${\omega}_{a}\simeq {\omega}_{0}\simeq 1{0}^{7}$ rad/s for 10 nm particles in an 1 mM salt. An important relaxation feature present in the spectrum is due to EP. EP is linked to ω_{b} as defined in the previous section, see Equation (38). We have already noted that in many experimental conditions ω_{b} ≲ (ω_{a}, ω_{0}). We will show that in the cases where this inequality applies, it is possible to make use of this fact by performing the fits on the measured data *uncorrected* for EP.

In the following Sections 5.4.1 and 5.4.2, we show how the logarithmic derivative method can be applied when the relation between ${\stackrel{~}{\epsilon}}_{s}$ and $\stackrel{~}{\beta}$ is known (section The direct fitting method) and when it is not known (The general fitting method). This last section is based on the work of Jimenez et al. (2007).

The “experimental data” needed in these sections to perform the fits are generated numerically using an equivalent circuit model, with the elements given in Equation (37). The necessary parameters are given in the legends of the figures. The “experimental data” mimic the response $\left({\epsilon}_{c,s}^{\prime},{\epsilon}_{c,s}^{\prime \prime}\right)$ of a colloidal suspension (including EP). From the ${\epsilon}_{c,s}^{\prime}$ thus obtained, one can derive

#### 5.4.1. The Direct Fitting Method (with a Dipolar Coefficient Model)

Theories have been developed to express $\stackrel{~}{\beta}$ in terms of the relevant parameters of the colloidal particles, namely particle size, zeta potential or surface charge and Stern layer parameters. These theories have been derived for the case of colloidal suspensions consisting of spherical, homogeneously charged particles (DeLacey and White, 1981; Delgado, 2002; Ohshima, 2006; Chassagne and Bedeaux, 2008). Some work has also been performed on homogeneously charged spheroids, see Chassagne and Bedeaux (2008) and references within. An analytical expression for $\stackrel{~}{\beta}$ in the case of spheres can be found in Supplementary Material 1. Other (numerical) models can be found in Mangelsdorf and White (1990) and Minor et al. (1998). We are going to use the analytical expression given in Supplementary Material 1, but we emphasize that any other theory can be applied, even numerical ones, even though the fitting procedure becomes more challenging in this case. Equation (8) is used to link ${\stackrel{~}{\epsilon}}_{s}$ (${\stackrel{~}{K}}_{s}$) to $\stackrel{~}{\beta}$. The logarithmic derivative ${\epsilon}_{D,s}^{\prime \prime}$ can be calculated numerically from ${\epsilon}_{s}^{\prime}\left(\stackrel{~}{\beta}\right)$ using the relation

The ${\epsilon}_{D,s}^{\prime \prime}$ thus obtained is a function of $\stackrel{~}{\beta}$ and therefore a function of the relevant parameters of the colloidal particles, the zeta potential in particular.

In Figure 11 we show two examples of the fit of ${\epsilon}_{D,c,s}^{\prime \prime}$ (black curve) using the ${\epsilon}_{D,s}^{\prime \prime}$ defined just above (red curve). The blue curve indicates the ${\epsilon}_{D,s}^{\prime \prime}$ calculated using the equivalent circuit defined in Section 4.2 with *C*_{EP} = 0 which is equivalent to take an infinite electrode spacing.

**Figure 11. Top:** Logarithmic derivative ${\epsilon}_{D}^{\prime \prime}$ for a suspension consisting colloidal spheres (ϕ = 1%, *eζ*/*kT* = 4) in a 1 mM electrolyte solution of monovalent salt for which *D*_{1} = 2 × 10^{−9} m^{2}/s and *D*_{2} = 3 × 10^{−9} m^{2}/s. The electrode spacing is 10 mm. (left): 25 nm particles and (right): 250 nm particles. The fit with the dipolar coefficient was done between [ω_{b} = 2.5 × 10^{6}-10^{10}] rad/s using (ζ, *a*) as adjustable parameters. **Bottom**: Recalculated $\stackrel{~}{\beta}$ (blue) and original beta (red).

For the two cases, which span a large particle size distribution at moderate ionic strength (1 mM), we are in the situation where ω_{a} ≳ ω_{b} (left figure) and ω_{a} ≃ ω_{b} (right figure). For all curves, the fits were always performed using the data in the frequency range above ω_{b}. Different situations were tested: we checked that is is possible to fit the data using (a) the zeta potential ζ as single adjustable parameter, (b) ζ and the particle size *a* as adjustable parameters. Both gave the same excellent fit as in Figure 11. The fits were also possible adding small random errors to the data (not shown). We checked that it is possible to use 3 adjustable parameters: ζ, *a* and a Stern layer parameter *St*, defined in Supplementary Material 1, if a Stern layer conductance is introduced in $\stackrel{~}{\beta}$ and similarly introduced in the numerically generated “experimental data” (not shown).

Interestingly, when ω_{a} ≃ ω_{b} (see Figure 11), the direct fitting method gives very good results despite the fact that only a part of the first relaxation peak is fitted. This is because ${\epsilon}_{s}^{\prime}\left(\stackrel{~}{\beta}\right)$ contains all possible relaxations processes of the colloid, and is therefore versatile enough to give a proper fit with few parameters.

#### 5.4.2. The General Fitting Method

This procedure, introduced by Jimenez et al. (2002, 2007) is based on two steps. The first one is to remove the electrode polarization (EP) contribution from the logarithmic derivative ${\epsilon}_{D,c,s}^{\prime \prime}$ data in order to assess ${\epsilon}_{D,s}^{\prime \prime}$, and the second step is to fit ${\epsilon}_{D,s}^{\prime \prime}$ using a using a Havriliak-Negami (HN) relaxation function. This method is quite general, and can in theory be applied to a variety of complex fluids, such as polymeric or colloidal suspensions. The use of the Havriliak-Negami (HN) function originates from work done to describe empirically the dielectric relaxation of polymers (Havriliak and Negami, 1967). Jimenez et al. showed it is applicable to colloidal suspensions, and we refer to their work for more details and examples of application of the method (Jimenez et al., 2002, 2005, 2007).

In order to remove the EP contribution to ${\epsilon}_{D,c,s}^{\prime \prime}$ Jimenez et al. introduce the relation:

where *m* should be close to 2 [see Equation (64)], as *Aω*^{−m} represents the EP contribution (*A* is a parameter to be fitted to the data). The contribution *Aω*^{−m} is fitted to the data in the low-frequency part (assuming that $A{\omega}^{-m}\gg {\epsilon}_{D,s}^{\prime \prime}$ in this part). The permittivity ${\epsilon}_{D,s}^{\prime \prime}$, (named ${\epsilon}_{D,cor}^{\prime \prime}$ by Jimenez et al.) represents the ${\epsilon}_{D}^{\prime \prime}$ of the suspension in the absence of EP. It is obtained by subtraction over the whole frequency range (${\epsilon}_{D,s}^{\prime \prime}={\epsilon}_{D,c,s}^{\prime \prime}-A{\omega}^{-m}$).

It is clear that the logarithmic derivative method gives ${\epsilon}_{s}^{\prime}={\epsilon}_{s}$ without its ω-independent part. The same holds for ${\epsilon}_{HN,s}^{\prime}$ i.e., the ${\epsilon}_{s}^{\prime}$ found from the logarithmic derivative of corresponding HN function. The expression for the HN function reads:

(Note that for *b* = 1 the function is called the Cole-Cole function while for *a* = *b* = 1 the function is called the Debye function). We can therefore define:

The parameters Δε_{HN, s} = ε_{HN, s}(0) − ε_{HN, s}(∞), τ, *a* and *b* are to be fitted to the data (ε_{HN, s}(∞), being a constant, is not included in the fit). The missing constant is found from ${\epsilon}_{HN,s}\left(\infty \right)={\epsilon}_{s}^{\prime}\left(\infty \right)$ because ${\epsilon}_{s}^{\prime}$ reaches a constant at high frequencies. Indeed, we then have from Equation (2): ${\epsilon}_{s}^{\prime}\left(\omega \to \infty \right)={\epsilon}_{e}\left(1+3\varphi \text{Re}\left(\stackrel{~}{\beta}\left(\omega \to \infty \right)\right)\right)={\epsilon}_{\infty}^{\prime}$. From the parameters found from the fit of ${\epsilon}_{D,s}^{\prime \prime}$ by ${\epsilon}_{D,HN,s}^{\prime \prime}$ one then finds $\left[{\epsilon}_{HN,s}^{\prime}-{\epsilon}_{HN}\left(\infty \right)\right]=\left[{\epsilon}_{s}^{\prime}-{\epsilon}^{\prime}\left(\infty \right)\right]$ using Equation (68). From Equation (68), we also find that ${\epsilon}_{HN,s}^{\prime \prime}\left(\omega \to 0\right)=0$. This implies that we can make the equivalence :

where *K*_{r} = *K*_{s}(ω → 0). The electrolyte with a conductivity having this property has been termed reference electrolyte in the Section 5.1.

We define the increment $\Delta {\stackrel{~}{K}}_{HN,s}$ which can be found directly from the parameters of the fit using Equation (68):

where ${\stackrel{~}{K}}_{r}\left(\omega \right)={K}_{r}+i\omega {\epsilon}_{0}{\epsilon}_{e}$ and where we have used, for a colloidal sphere with relative permittivity ε_{2} [see Supplementary Material 1 for a general expression of $\stackrel{~}{\beta}\left(\omega \right)$]

From Equation (2) we deduce that (we recall that ${\stackrel{~}{K}}_{e}={\stackrel{~}{K}}_{1}$):

Combining Equations (71) and (73), it is possible to find the required dipole coefficient $\stackrel{~}{\beta}$ by using the relation:

Note that ${\stackrel{~}{K}}_{r}$ and ${\stackrel{~}{K}}_{1}$ can be obtained by measuring the conductivities of the corresponding electrolytes. In order to estimate β_{inf} only the value of ε_{2} is required (the value of the relative permittivity of water, i.e., ε_{e}, is known). It is also possible to determine β_{inf} experimentally from

As discussed at the beginning of this section, there are two characteristic frequencies associated with the colloidal particle's relaxations. It is therefore possible to fit both relaxation frequencies when they are experimentally accessible, by using a sum of two HN functions such that:

which implies the use of 8 fitting parameters: Δε_{HN, 1}, τ_{1}, *a*_{1}, *b*_{1}, Δε_{HN, 2}, τ_{2}, *a*_{2}, *b*_{2}. Equations (71)–(74) can subsequently be applied to find $\stackrel{~}{\beta}$. We have checked that this procedure gives good results, see Figure 12). In this example, we have fitted the uncorrected measured data, in the frequency range above ω_{b}. For the smallest particle used (25 nm) ω_{a} ≳ ω_{b} the fit is very good. For larger colloidal particles (250 nm particles), we find that despite the fact that ω_{a} ≃ ω_{b} and that the double HN is less versatile than the dipolar coefficient used in the direct fitting method discussed in the section above, the result of the fit is surprisingly good, especially at low frequencies. This is due to the fact that the value of $\stackrel{~}{\beta}$ at low frequencies depends nearly exclusively on the value of *K*_{r} = *K*_{s}(ω → 0) which is not affected by EP, and does not depend on $\Delta {\stackrel{~}{K}}_{HN}$:

The correct fit of the second peak (associated with the double layer relaxation and ω_{0}) is, on the other hand, very important to get the correct relaxation of $\stackrel{~}{\beta}$ in the high frequency range.

**Figure 12. Top:** Logarithmic derivative ${\epsilon}_{D}^{\prime \prime}$ for a suspension consisting of colloidal spheres(ϕ = 1%, *eζ*/*kT* = 4) in a 1 mM electrolyte solution of monovalent salt for which *D*_{1} = 2 × 10^{−9} m^{2}/s and *D*_{2} = 3 × 10^{−9} m^{2}/s. (left): 25 nm particles and (right): 250 nm particles. The electrode spacing is 10 mm, and the 2 HN fit was done between [2.5.10^{6}–10^{10}] rad/s. **Bottom**: recalculated $\stackrel{~}{\beta}$ [blue line, from Equation (74) and original beta (red)]: despite the inaccuracy in the 2 HN fit, the error in $\stackrel{~}{\beta}$ is small, see explanation in text.

## 6. Conclusions

In this article, we have discussed the dielectric spectroscopy response of electrolytes and colloidal suspensions. We concentrated ourselves on the elimination of the phenomenon called “electrode polarization” (EP) that can overshadow the measured signal. This phenomenon is caused by the build-up of charges close to the electrodes in the suspension (or electrolyte solution). We have not discussed the other source of unwanted effects, namely stray impedances that originate from the non-ideality of the experimental set-up, and that also affect the measured signal. From the data, important parameters for the characterization of colloidal particles can be assessed. These parameters are found from the evaluation of the dipolar coefficient $\stackrel{~}{\beta}$ which is obtained from the cleaned and/or fitted data.

It was found that several characteristic frequencies could be identified and linked to EP effects. An important frequency is ${\omega}_{b}=\sqrt{{\omega}_{0}{\omega}_{P}}$ with ω_{P} = 2κ*D*_{±}/*d* and ${\omega}_{0}={\kappa}^{2}{D}_{\pm}$ where κ is the inverse of the Debye length, *D*_{±} the characteristic ionic diffusion coefficient and *d* the distance between electrodes (planar electrodes) or the smallest of the two radii in the case of cylindrical electrodes. Above ω_{b} the measured data are not affected anymore by EP. This enables, for many cases encountered in practice, to assess the properties of the colloidal particles by analyzing directly the data uncorrected for EP. This is shown in the last section. In that section, two fitting procedures are described: first, we consider “the direct fitting method” that makes use of existing dipolar coefficient models for colloidal suspensions, and then we explain the “general fitting method” which is based on the logarithmic derivative method and makes use of Havriliak-Negami (HN) relaxation functions. This last method is especially suited when the fluid under investigation is complex, i.e., when no good model for its behavior as function of the frequency exists. The correction of the experimental data for EP with the logarithmic derivative thus is quite successful. It has the advantage that it can be achieved simply by invoking ∂ε′/∂ ln ω in the data analysis. No special adaptations have to be made to the measurement or the equipment.

Also some experimental techniques enable us to compensate for EP. In this article, we tested them from a theoretical perspective. As outlined in the corresponding sections, even though a model could be ideal to compensate for EP in a given range of frequency, there are often technical limitations that could limit their validity. However, from the theoretical modeling, we have found that:

1 - The “subtraction method” can compensate for electrode polarization above the critical frequency ω_{P} provided that the dielectric permittivity ε* in the small slab of liquid adjacent to the electrodes is given by ε_{e}. If this permittivity differs from ε_{e}, the subtraction method is not suited for eliminating EP. The fact that ε* may or may not be affected by the presence of colloidal particles in the slab of liquid adjacent to the electrodes is still an open question.

2 - The “variable electrode separation method” is well-suited for compensating for electrode polarization over the whole frequency range.

3 - The “4-electrode cell method” gives directly access to the wanted signal, devoid of EP over the whole frequency range.

## Author Contributions

All the authors made substantial contributions to the conception of the work; They all helped revising it critically for important intellectual content; They all approved the final version to be published; They agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

Partial financial support of the Agence Nationale de la Recherche in the frame of the project Celadyct (ANR-12-BS08-0017-01) is gratefully acknowledged. Financial support of this investigation by Junta de Andalucía, Spain (grant No. PE2012-FQM0694) and MINECO, Spain (project No. FIS2013-47666-C3-1-R) is gratefully acknowledged. The authors also would like to thank the TU Delft Library for the financial support that made it possible to publish this article in open access. One of the authors, CC would like to warmly thank Dr. S. Gourdin-Bertin for many enlightening discussions, while JT is grateful to Prof. Sudhölter and Prof. Kreutzer for allowing him a stay in the Dept. of Chemical Engineering. All the authors would like to thank the referees. Especially, the detailed comments made by referee 1 enabled them to improve the article, while the comments of referee 3 prompted them to add Supplementary Material 4.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fchem.2016.00030

## References

Barsoukov, E., and Macdonald, J. R. (2005). *Impedance Spectroscopy Theory, Experiments and Applications, 2nd Edn*. NewYork, NY: Wiley. doi: 10.1002/0471716243

Beltramo, P. J., and Furst, E. M. (2012). Transition from dilute to concentrated electrokinetic behavior in the dielectric spectra of a colloidal suspension. *Langmuir* 28, 10703–10712. doi: 10.1021/la301876w

Beltramo, P. J., and Furst, E. M. (2012). A simple, single-measurement methodology to account for electrode polarization in the dielectric spectra of colloidal dispersions. *Chem. Lett.* 41, 1116–1118. doi: 10.1246/cl.2012.1116

Buck, R. P. (1969). Diffuse layer charge relaxation at the ideally polarized electrode. *J. Electroanal. Chem.* 23, 219–240. doi: 10.1016/S0022-0728(69)80212-3

Chassagne, C., and Bedeaux, D. (2008). The dielectric response of a colloidal spheroid. *J. Coll. Interface Sci.* 326, 240–253. doi: 10.1016/j.jcis.2008.06.055

Chassagne, C., Bedeaux, D., van der Ploeg, J. P. M., and Koper, G. J. M. (2002). Theory of electrode polarization: application to parallel plate cell dielectric spectroscopy experiments. *Colloids Surf. A* 210, 137–145. doi: 10.1016/S0927-7757(02)00379-5

Chassagne, C., Bedeaux, D., van der Ploeg, J. P. M., and Koper, G. J. M. (2003). Polarization between concentric cylindrical electrodes. *Physica A* 326, 129–140. doi: 10.1016/S0378-4371(03)00286-3

Chassagne, C., Mietta, F., and Winterwerp, J. C. (2009). Electrokinetic study of kaolinite suspensions. *J. Coll. Interface Sci.* 336, 352–359. doi: 10.1016/j.jcis.2009.02.052

Cirkel, P. A., van der Ploeg, J. P. M., and Koper, G. J. M. (1997). Electrode effects in dielectric spectroscopy of colloidal suspensions. *Physica A* 235, 269–278. doi: 10.1016/S0378-4371(96)00347-0

DeLacey, E. H. B., and White, L. R. (1981). Dielectric response and conductivity of dilute suspensions of colloidal particles. *J. Chem. Soc. Faraday Trans.* 77, 2007–2039. doi: 10.1039/f29817702007

DeLacey, E. H. B., and White, L. R. (1982). The polarization impedance of an ideally polarizable plane electrode. *J. Chem. Soc. Faraday Trans.* 78, 457–479. doi: 10.1039/f29827800457

Delgado, A. V. (ed.). (2002). *Interfacial Electrokinetics and Electrophoresis.* New York, NY: Dekker.

Fricke, H., and Curtis, H. J. (1937). The dielectric properties of water-dielectric interphases. *J. Phys. Chem.* 41, 729–745. doi: 10.1021/j150383a011

Grimnes, S., and Martinsen, O. J. (2015). *Bioimpedance and Bioelectricity. 3rd Edn.* Amsterdam: Elsevier.

Grosse, C., and Delgado, A. V. (2010). Dielectric dispersion in aqueous colloidal systems. *Curr. Opin. Coll. Interface Sci.* 15, 145–159. doi: 10.1016/j.cocis.2009.11.004

Grosse, C., and Delgado, A. V. (2013). Corrigendum to dielectric dispersion in aqueous colloidal systems. *Curr. Opin. Coll. Interface Sci.* 18, 161–163. doi: 10.1016/j.cocis.2013.02.007

Grosse, C., Tirado, M., Pieper, W., and Pottel, R. (1998). Broad frequency range study of the dielectric properties of suspensions of colloidal polystyrene particles in aqueous electrolyte solutions. *J. Coll. Interface Sci.* 205, 26–41. doi: 10.1006/jcis.1998.5587

Grosse, C. (2002). “Relaxation mechanism of homogenous particles and cells suspended in aqueous electrolyte solutions,” in *Interfacial Electrokinetics and Electrophoresis*, ed A. V. Delgado (NewYork, NY: Dekker), 277–327.

Gunning, J., Chan, D. Y. C., and White, L. R. (1995). The impedance of the planar diffuse double layer: an exact low-frequency theory. *J. Coll. Interface Sci.* 170, 522–537. doi: 10.1006/jcis.1995.1131

Havriliak, S., and Negami, S. (1967). A complex plane representation of dielectric and mechanical relaxation processes in some polymers. *Polymer* 8, 161–210. doi: 10.1016/0032-3861(67)90021-3

Hayakawa, R., Kanda, H., Sakamoto, M., and Wada, Y. (1975). New apparatus for measuring the complex dielectric constant of a highly conductive material, Japan. *J. Appl. Phys.* 14, 2039–2052. doi: 10.1143/JJAP.14.2039

Hollingsworth, A. D., and Saville, D. A. (2003). Broad frequency range dielectric spectrometer for colloidal suspensions: cell design, calibration, and validation. *J. Coll. Interface Sci.* 257, 65–76. doi: 10.1016/S0021-9797(02)00029-2

Hollingsworth, A. D., and Saville, D. A. (2004). Dielectric spectroscopy and electrophoretic mobility measurements interpreted with the standard electrokinetic model. *J. Coll. Interface Sci.* 272, 235–245. doi: 10.1016/j.jcis.2003.08.032

Hollingsworth, A. D. (2013). Remarks on the determination of low-frequency measurements of the dielectric response of colloidal suspensions. *Curr. Opin. Coll. Int. Sci.* 18, 157–159. doi: 10.1016/j.cocis.2013.01.002

Ishai, P. B., Talary, M. S., Caduff, A., Levy, E., and Feldman, Y. (2013). Electrode polarization in dielectric measurements: a review. *Meas. Sci. Technol.* 24:102001. doi: 10.1088/0957-0233/24/10/102001

Jimenez, M. L., Arroyo, F. J., Carrique, F., and Delgado, A. V. (2007). Surface conductivity of colloidal particles: experimental assessment of its contributions. *J. Coll. Interface Sci.* 316, 836–843. doi: 10.1016/j.jcis.2007.07.016

Jimenez, M. L., Arroyo, F. J., Carrique, F., Kaatze, U., and Delgado, A. V. (2005). Determination of stagnant layer conductivity in polystyrene suspensions-temperature effects. *J. Coll. Interface Sci.* 281, 503–509. doi: 10.1016/j.jcis.2004.08.093

Jimenez, M. L., Arroyo, F. J., van Turnhout, J., and Delgado, A. V. (2002). Analysis of the dielectric permittivity of suspensions by means of the logarithmic derivative of its real part. *J. Coll. Interface Sci.* 249, 327–335. doi: 10.1006/jcis.2001.8141

Kaatze, U., and Feldman, Y. (2006). Broadband dielectric spectrometry of liquids and biosystems. *Meas. Sci. Technol.* 17, R17–R35. doi: 10.1088/0957-0233/17/2/R01

Kang, K., and Dhont, J. K. G. (2010). Electric-field induced transitions in suspensions of charged colloidal rods. *Soft Matter* 6, 273–286. doi: 10.1039/B910046F

Kijlstra, J., and Wegh, R. A. J. (1994). van Leeuwen HP. Impedance spectroscopy of colloids. *J. Electroanal. Chem.* 366, 37–42. doi: 10.1016/0022-0728(93)02121-W

Loeb, A. L., Overbeek, J. T. H. G., and Wiersema, P. H. (1961). *The Electrical Double Layer around a Spherical Colloidal Particle*. Cambridge, MA: MIT Press.

Lvovich, V. F. (2012). *Impedance Spectroscopy-Applications to Electrochemical and Dielectric Phenomena*. Hoboken, NJ: Wiley. doi: 10.1002/9781118164075

Mangelsdorf, C. S., and White, L. R. (1990). Effects of Stern-layer conductance on electrokinetic transport properties of colloidal particles. *J. Chem. Soc. Faraday Trans.* 86, 2859–2870. doi: 10.1039/ft9908602859

Mazzeo, B. A. (2009). Parasitic capacitance influence of potential sensing electrodes on 4 electrode liquid impedance measurements. *J. Appl. Phys.* 105:094106. doi: 10.1063/1.3124365

Minor, M., Van Leeuwen, H. P., and Lyklema, J. (1998). Low-frequency dielectric response of polystyrene latex dispersions. *J. Coll. Interface Sci.* 206, 397–406. doi: 10.1006/jcis.1998.5619

Myers, D. F., and Saville, D. A. (1989). Dielectric spectroscopy of colloidal suspensions. I Dielectric spectrometer. *J. Coll. Interface Sci.* 131, 448–460. doi: 10.1016/0021-9797(89)90187-2

Ohshima, H., Healy, T. W., and White, L. R. (1982). Accurate analytic expressions for the surface charge density/surface potential relationship and double-layer potential distribution for a spherical colloidal particle. *J. Coll. Interface Sci.* 90, 17–26. doi: 10.1016/0021-9797(82)90393-9

Pelc, D., Marion, S., and Basletic, M. (2011). Four-contact impedance spectroscopy of conductive liquid samples. *Rev. Sci. Instrum.* 82:073907. doi: 10.1063/1.3609867

Rica, R. A., Jimenez, M. L., and Delgado, A. V. (2010). Electric permittivity of concentrated suspensions of elongated suspensions goethite particles. *J. Coll. Interface Sci.* 343, 564–573. doi: 10.1016/j.jcis.2009.11.063

Schwan, H. P. (1992). Linear and nonlinear electrode polarization and biological materials. *Ann. Biomed. Eng.* 20, 269–288. doi: 10.1007/BF02368531

Scott, M., Paul, R., and Kaler, K. V. I. S. (2000a). Theory of frequency-dependent polarization of general planar electrodes with zeta potentials of arbitrary magnitude in ionic media, Pt 1. Theoretical foundations and general results. *J. Coll. Interface Sci.* 230, 377–387. doi: 10.1006/jcis.2000.7125

Scott, M., Paul, R., and Kaler, K. V. I. S. (2000b). Theory of frequency-dependent polarization of general planar electrodes with zeta potentials of arbitrary magnitude in ionic media, Pt 2. Applications and results from homogeneous and array systems of electrodes. *J. Coll. Interface Sci.* 230, 388–395. doi: 10.1006/jcis.2000.7126

Shilov, V. N., Delgado, A. V., Gonzalez-Caballerro, F., and Grosse, C. Donath, E. (2002). “Suspensions in an alternating electric field: dielectric and electrorotation spectroscopies,” in *Interfacial Electrokinetics and Electrophoresis*, ed A. V. Delgado, (New York, NY: Dekker), 329–367.

Steeman, P. A. M., and van Turnhout, J. (1994). Fine structure in the parameters of dielectric and viscoelastic relaxations. *Macromolecules* 27, 5421–5427. doi: 10.1021/ma00097a023

van der Ploeg, J. P. M., and Mandel, M. (1991). A low frequency relative permittivity meter for electrically conducting liquids. *Meas. Sci. Technol.* 2, 389–395. doi: 10.1088/0957-0233/2/4/019

van der Touw, F., and Mandel, M. (1971). Plane-parallel condenser with variable electrode spacing for determination of electric permittivity of highly conducting liquids below 1 MHz. *J. Chem. Soc. Trans. Faraday Soc.* 67, 1336–1342. doi: 10.1039/tf9716701336

van der Touw, F., Selier, G., and Mandel, M. (1975). Measuring device for the determination of the electric permittivity of conducting liquids in the frequency range 2-500 kHz.II The measuring cell. *J. Phys. E* 8, 844–846. doi: 10.1088/0022-3735/8/10/016

Keywords: colloidal suspension, complex conductivity and permittivity, electrode polarization

Citation: Chassagne C, Dubois E, Jiménez ML, van der Ploeg JPM and van Turnhout J (2016) Compensating for Electrode Polarization in Dielectric Spectroscopy Studies of Colloidal Suspensions: Theoretical Assessment of Existing Methods. *Front. Chem*. 4:30. doi: 10.3389/fchem.2016.00030

Received: 07 December 2015; Accepted: 17 June 2016;

Published: 19 July 2016.

Edited by:

Alireza Mashaghi, Harvard University, USAReviewed by:

Rafael Roa, Helmholtz-Zentrum Berlin, GermanyAndrew David Hollingsworth, New York University, USA

Jan Karel George Dhont, Forschungszentrum Jülich, Germany

Copyright © 2016 Chassagne, Dubois, Jiménez, van der Ploeg and van Turnhout. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Claire Chassagne, c.chassagne@tudelft.nl