# Better Resolved Low Frequency Dispersions by the Apt Use of Kramers-Kronig Relations, Differential Operators, and All-In-1 Modeling

- Department of Chemical Engineering, Sect. Organic Materials and Interfaces, Delft University of Technology, Delft, Netherlands

The dielectric spectra of colloidal systems often contain a typical low frequency dispersion, which usually remains unnoticed, because of the presence of strong conduction losses. The KK relations offer a means for converting ε′ into ε″ data. This allows us to calculate conduction free ε″ spectra in which the l.f. dispersion will show up undisturbed. This interconversion can be done on line with a moving frame of logarithmically spaced ε′ data. The coefficients of the conversion frames were obtained by kernel matching and by using symbolic differential operators. Logarithmic derivatives and differences of ε′ and ε″ provide another option for conduction free data analysis. These difference-based functions actually derived from approximations to the distribution function, have the additional advantage of improving the resolution power of dielectric studies. A high resolution is important because of the rich relaxation structure of colloidal suspensions. The development of all-in-1 modeling facilitates the conduction free and high resolution data analysis. This mathematical tool allows the apart-together fitting of multiple data and multiple model functions. It proved also useful to go around the KK conversion altogether. This was achieved by the combined approximating ε′ and ε″ data with a complex rational fractional power function. The all-in-1 minimization turned out to be also highly useful for the dielectric modeling of a suspension with the complex dipolar coefficient. It guarantees a secure correction for the electrode polarization, so that the modeling with the help of the differences ε′ and ε″ can zoom in on the genuine colloidal relaxations.

## Introduction

Dielectric spectroscopy is a powerful technique to study electrokinetic phenomena, it determines the frequency dependence of the real and imaginary part of the permittivity, ε′ and ε″. The measurements can span a very broad frequency range. In electrokinetic spectroscopic studies the low frequency part is quite relevant. At low frequencies the dissipative loss in ε″ caused by ohmic conduction becomes prominent, often to such an extent that it overshadows the genuine dispersion or pure relaxation losses of the colloidal particles. The unwanted ohmic loss might thus hamper the resolution of nearby relaxation processes seriously. So, in order to reach a high resolution the contribution of the ohmic conduction to ε″ should be eliminated. Conduction may also reduce the diagnostic selectivity of electrochemical sensors, so cleaning ε″ from conduction might be beneficial here as well.

Figure 1 illustrates the problem caused by ohmic conduction. In order to reveal the hidden peak due to the true relaxation losses of the colloidal particles the ohmic conduction should be removed from the measured ε″-data.

**Figure 1. Simulated observed ε″-spectrum and that of ε″ calculated by converting ε′-data into conduction free ε″-ones**. Only the latter shows us the authentic loss peak we look for. The ε′-spectrum remains unaffected by conduction. We assumed that the electrode polarization, which may also turn up at low frequencies, can be neglected, cf. Figure 3.

Clearly, the disturbing ohmic energy loss will not show up in the real part of the permittivity ε′, because ε′ is a measure of the number of dipoles and ions involved, but not of the energy required for their motion. Now one of the two Kramers-Kronig dispersion relations offers a means to calculate from the ε′ data, the genuine relaxation loss in the imaginary part ε″. Unfortunately such a KK interconversion is made difficult by the fact that the KK relations are singular integrals. We will describe two methods to approximate the unwieldy integrals accurately in an easy and fast way. This has the advantage that the ε′ to ε″ conversion can be coded as one-liners and can be used on line. The ε″ data obtained from ε′ will be denoted ε″_{kk} or ε″_{cf} (with cf being short for conduction free). The fast interconversion basically consists of a *frame* of coefficients by which ε″_{kk} can be calculated from a limited set of logarithmically spaced ε′ data clustered around the conversion frequency. Our compact KK relations in the form of moving conversion frames can in reverse be used to find the strength of the conduction loss from the difference of the observed and the converted ε″ data. This information is for example crucial for assessing the onset of percolation when the conducting phase in a colloidal mixture becomes co-continuous. A versatile way to avoid the problematic KK integrals is kernel matching. This mathematical tool relies on the fact that the permittivity can be considered to originate from a continuous distribution of elementary Debye relaxation processes. The distributed ε′ and ε″ do obey the KK relations but these improper integrals are no longer needed, because the conversion can now be accomplished by approximating the Debye kernel of the ε″ distribution-integral with a sum of logarithmic spaced Debye kernels of the ε′ distribution-integral. This singularity free approach thus provides the desired conversion panel with which ε″_{kk} can be uncovered readily by moving the panel along the observed logarithmically spaced ε′ data. The other method we have explored for the fast evaluation of the KK integrals is the approach of “integration by differentiation.” In this route the KK integrals are replaced by symbolic differential operators. It turned out that the logarithmic differential operators cannot be used in a broad sense. However, one operator, viz. the cot-operator scheme could be made useful for calculating *dε*′∕*dlnω* from a narrow window of logarithmically spaced ε″data. This logarithmic derivative can of course also be calculated from ε′ proper, in that case it will be automatically conduction free. The ‘loss’ peaks appearing in the *dε*′∕*dlnω* spectra have the advantage of being sharper than their corresponding ε″ counterparts, which implies that the resolution increases^{1}. The various conversions we will discuss are summarized in Figure 2. In addition to ${\epsilon}_{kk}^{\u2033}$ and *dε*′∕*dlnω*, we will also pay attention to the special features of the logarithmic differences of ε′ and ε″.

**Figure 2. Options explored for converting ε′ and ε″ data**. We will show that the differences of ε′ and ε″ offer new possibilities for improving the data analysis. All conversions are done with panels and executed as one-liners.

Another option to increase the resolution of lf-dispersions or actually of all dispersions is to compute the underlying distribution function by inverting the ε′ and/or ε″ data. For the latter we clearly should take its pure relaxation part, i.e., ε″_{kk} or ε″_{cf}. Recall that a single Debye relaxation is narrowed down upon inversion to an extremely sharp delta distribution. The ensuing distribution spectra will therefore show the highest resolution possible. We will discuss a simple way to accomplish the complex Stieltjes inversion of the ε′ and ε″ data by making use of rational polynomials in fractional powers. A joint data-analysis by fitting ε′ and ε″ “apart-together” leads to an improvement in resolution power as well. Such a paired simultaneous modeling, dubbed *all-in-one* modeling, wherein the relaxation parameters are kept the same in the ε′ and ε″ fit functions, but whereby the conduction is contained only in the ε″ fit function, requires a two-way switch in the co-fit procedure which links the proper data to the proper fit formulae and thus assures that the two non-linear least squares minimizations are always done in parallel. Although we will focus on the data handling of the complex permittivity, the fast conversion methods developed can also be applied directly to other complex electrokinetic quantities like the dielectric modulus and the impedance or to the magnitude and phase.

Apart from the ohmic conduction, another disturbance is prominently active in the dielectric spectroscopy of aqueous colloidal systems. The point being that the ions of the electrolyte tend to accumulate near the electrodes of the measuring cell. This causes electrode polarization (EP), which may enhance ε′ and ε″ strongly. This phenomenon too may therefore overshadow the true l.f. relaxation of the colloidal particles. We will briefly describe in Section Improving the Resolution by All-In-1 Modeling of the Real and Imaginary Data, how the effect of this nuisance can be eliminated as well as. Several methods to accomplish this are discussed recently in more detail in an upcoming article of this journal (Chassagne et al., submitted) and in van Turnhout et al. (2016).

Figure 3 illustrates the impact of the EP. The local ion motions in tune with the a.c. voltage create a special relaxation peak, which we have called ρ or space charge peak. Clearly, the calculated ${\epsilon}_{kk}^{\u2033}$ spectrum will also reveal the presence of this specific space charge relaxation peak.

**Figure 3. The l.f. ρ or space charge relaxation, which is due to the rather slow up and down motion of the ion clouds near the electrodes, remains hidden in the ε″ curve**. It too can be disentangled from the conduction contribution by calculating ε″_{kk}. The conduction, which stems from the gross ion motion, appears alongside the EP, if the electrodes are not fully blocking.

**Figure 4. Kramers (left), in particular Kronig (right) was relatively young when he proposed one of the KK integrals**. At that time he added de Laer to his last name. His name should not be spelled Krönig.

## KK Relations-the Basics of the Interconversion of Real and Imaginary Data

The KK dispersion relations were proposed by Kramers and Kronig about a century ago. A candid review of their history has been given by Bohren (2010). The relations were given as integrals that relate the real and imaginary part of all dispersion phenomena. Kramers proposed both integrals for the first time in Copenhagen (Kramers, 1927). Kronig proposed one of them a year later (Kronig, 1926). Both did not refer in their accounts to the Hilbert transform and it inverse, which are basically the same.

Electrical engineers speak of “real part sufficiency,” which implies that the imaginary part can be calculated from the real part of the response of an electric network. They present the integrals in their textbooks without reference to KK (see e.g., Turtle, 1958). It took some time before the KK relations were actually used in the research of dielectric phenomena. Early pioneers were for instance: Bayard (1935), Gorter and Kronig (1936), Gross (1943), Brachman and Macdonald (1954).

The KK relations can be derived by making use of Cauchy's integral theorem (Kremer and Schönhals, 2002). They have the following form for dielectric dispersions or relaxations

where ε′ and ε″ are the real and imaginary part of the complex permittivity ε*(ω) = ε′(ω) − iε″(ω). Similar integrals hold for the other quantities used to describe dielectric relaxations such as the complex dielectric modulus, m*(ω) = 1∕ε*(ω), the magnitude, |ε*(ω)|, and phase θ(ω) = atn[ε″(ω)/ε′(ω)]. We shall restrict ourselves mainly to the interconversion of ε′ to ε″ and vice versa. In this process we will derive so-called conversion frames that can equally well be used to convert m′ into m″, etc.

Considering the integrals given, it is not surprising that in practice the KK relations are still not broadly used. This is due to the fact that their kernels become singular for ω → ω_{o}. However, this singularity can be removed by rewriting them to

Nonetheless the conversion of ε′ into ε″-values with an integral like Equation (2.2) is a tedious job. The more so because we want to know the converted values across the whole frequency scan. We have therefore discarded numerical integration and have followed two different approaches, see Figure 5.

**Figure 5. On the left, numerical integration is not the favorite tool for KK conversion**. It can be better realized with kernel matching and with symbolic differential operators. The conversion produces useful spin offs, they are depicted on the right. The higher resolution results from the elimination of the conduction from the ε″-data.

Before outlining this in Sections Easy to Use Conversion Frames Obtained by Matching Debye Kernels and Computing Conversion Frames with Symbolic Differential Operators, we should point out an important salient property of the KK relations, which is often overlooked. By the very fact that ε″ can be calculated from ε′, we miss out any possible contribution by ohmic conduction. In other words by invoking the KK relations we get a special set of ε″-data that are conduction free. We should therefore label the converted data with cf or kk and denote them as ${\epsilon}_{kk}^{\u2033}$ or ${\epsilon}_{cf}^{\u2033}$. Since the ohmic conduction causes a loss component of ${\epsilon}_{c}^{\u2033}(\omega )=\gamma \u2215({\epsilon}_{o}\omega )$ we have

where ε″(ω) represents the measured ε″-data, ${\epsilon}_{c}^{\u2033}(\omega )$ the conductive part in ε″, γ the ohmic conductivity and ε_{o} the permittivity of vacuum.

The KK relations can be of help in the removal of the dissipative ohmic loss only from ε″. They cannot be recruited for its removal from the other dielectric quantities mentioned. Let us illustrate this for the magnitude-phase relation. Writing ε*(ω) in its polar Euler form

and taking the natural logarithm we get:

with ε_{a}(ω) equaling the absolute value |ε*(ω)| and θ(ω) = atn[ε″(ω)/ ε′(ω)].

The real and imaginary parts lnε_{a} and θ obey KK relations similar to those of Equation (2.1). This means that it is possible to compute θ from lnε_{a} values, but this does not imply that in doing so the ${\epsilon}_{c}^{\u2033}$ contribution to θ is nullified. The reason being that the real part lnε_{a} also contains a contribution from ${\epsilon}_{c}^{\u2033}$. Recall that ${\epsilon}_{a}(\omega )=\sqrt{{\epsilon}^{\prime}{(\omega )}^{2}+{\epsilon}^{\u2033}{(\omega )}^{2}}$, so ε′ and ε″ are now mixed up in the real part.

Zahner (Germany) has implemented in their software this “logarithmic” version of the KK relations (Schiller et al., 2001; Lasia, 2014). It embodies an algorithm to calculate lnε_{a} from measured θ values. We will return to this in Sect. 4. The logarithmic version of the KK relations is often used by electrical engineers (Turtle, 1958).

Figure 6 shows the various conversions we will discuss below. For each we will give the coefficients of the appropriate conversion frames. The main focus will be on getting one liners for the conversion of ε′ to ε″_{kk} and ε″ to *dε*′∕*dlnω*. Although, as just said, we will also shortly touch upon the conversion of θ to lnε_{a} or for that matter of ε″ to ε′.

**Figure 6. Illustration of the different conversions covered**. The resulting ε″_{kk} and *dε*′∕*dlnω* will improve the l.f. spectral resolution, both with respect to ohmic conduction and to EP.

## Easy to Use Conversion Frames Obtained by Matching Debye Kernels

The direct use of the unwieldy KK integrals can be avoided by making use of the integral equations based on the description of relaxation phenomena with a distribution in relaxation times.

The most elementary dielectric dispersion phenomenon is the Debye relaxation, which Debye derived for dipole relaxations in gases and liquids. We can view the intricate relaxation of colloidal systems as being buildup from a sum of Debye processes

where ε_{∞} is ε′(∞), τ_{k} the relaxation time of process k. The relaxation strength Δε_{k} can be expressed as Δε_{k} = ε_{sk} − ε_{∞} with $\sum _{k=1}^{n}{\epsilon}_{sk}={\epsilon}^{\prime}(0)$. Often the individual Debye processes may come very close, so we better go a step further and think of the total dielectric response arising from a continuous distribution of Debye processes f (τ)

in which we can split the complex Debye kernel into a real and imaginary part

We can prove that these integrals obey the KK relations by inserting them into Equation (2.1).

We now have 2 additional relations between ε′ and ε″ that are linked up by the distribution function. Expressions like Equation (3.3) are common in the theory of viscoelastic phenomena. It were Ninomiya and Ferry (1959) who were the first to suggest a powerful trick to manipulate the various viscoelastic interrelations. By the way they did not deal with the interconversion of Equation (3.3) as such, but addressed e.g., the conversion of time to frequency responses and vice versa.

The clue, as sketched in Figure 7, is to approximate the respective Debye kernels. We have

where x = ωτ. Hence for the conversion of ε′ to ε″ we should approximate ε″ with a sum of κ′(x). A good option is

**Figure 7. Kernel matching is a powerful method for computing the coefficients of conversion frames, in this case for ε′ → ε″ _{kk}**. With the frames we can obtain ε″

_{kk}with one liners, cf. Equations (3.7) and (3.8).

The subtraction of ½ makes the κ′ kernel antisymmetric around *x* = 1 and assures that the a_{k}'s become symmetric. By prescribing the logarithmic spacing of the κ′ terms at *h* = 2, we have reduced the least squares approximation to a linear one, which is easy to do. We found:

These coefficients differ somewhat from the ones, we have given earlier (Steeman and van Turnhout, 1997; Wübbenhorst and van Turnhout, 2002; van Turnhout, 2004). We call this set of coefficients a *conversion* frame. With such a frame or panel we can convert a set of log-spaced ε′ data easily to a ${\epsilon}_{kk}^{\u2033}$value. In fact

with a_{4k} = {0.40378, 0.300158, −0.17918, 0.16192}. The data acquisition of modern test equipment (e.g., from Novocontrol, HP-Agilent, now Keysight Technologies, etc) allow us to measure at log-spaced intervals, whereby the choice of a spacing *h* = 2 poses no problem.

Sampling the list of measured ε′ by sliding the frame from the beginning to the end step by step along all data creates a new set of ε″ data that are conduction free. We notice that for each ε″ value we need 4 ε′ data upfront and 4 at the rear. This means that we will lose 4 ε″ at the start and 4 at the end.

For that reason we have also calculated asymmetric conversion frames by starting the kernel-approximation of Equation (3.5) not from *k* = −4, but from −3, −2, −1. Using these start up frames we only miss the very first ${\epsilon}_{kk}^{\u2033}$ value. By starting the sum in Equation (3.5) from *k* = −7, −6, −5, we get special end frames by which we lose out only one ${\epsilon}_{kk}^{\u2033}$-value at the end of the frequency scan.

Most dielectric measurements cover a very wide ω range. Then the loss of a few ε″_{kk} data at the start and the end of the range is hardly serious. In that case we can restrict the conversion to that with the central frame and simply use Equation (3.7).

Table 1 lists the coefficients of the ε′ to ε″ frames. We have only included the central frame and those for the start, because those at the end are in fact the same as those at the start albeit in reversed order and with a minus-sign.

Implemented in Mathematica, the frames with the coefficients b_{1k} to b_{4k} from top to bottom can be used as follows:

At the start from ω_{1} we stick to the first 9 ε′ values to get ε″(2ω_{1}), ε″(4ω_{1}), ε″(16ω_{1}). Next we can march on and drop the first ε′(ω_{1}) and add a new ε′-value at the end of the central frame to get ε″(32ω_{1}), ε″(64ω_{1}), etc. This can be continued till we reach the end of the total of n ε′-data at ω_{n}. We then stick to the last 9 ε′ data to get ε″[ω_{n}/16], ε″[ω_{n}/8], ε″[ω_{n}/2]. The present frames give a more accurate conversion than the ones give a more accurate conversion than the ones given earlier (Steeman and van Turnhout, 1997).

Figure 8 shows the accuracy achieved with the central panel for the ε″ data of a Debye relaxation (the sharpest relaxation possible). The recalculated values deviate a bit at the wings. This deviation becomes much less if the begin and end panels are used as well, as suggested in Equation (3.8).

**Figure 8. KK conversion tested for a Debye relaxation with Equation (3.7)**. The dashed line pertains to ε″_{kk} calculated.

Brather has also discussed KK conversion based on kern matching (Brather, 1979). He proposed the use of a series of ε′ terms. These large sums can be applied less easily used than our conversion frames, which allow conversion on line. He has also not discussed the special sums needed to begin and finish the conversion.

We should realize that the ohmic conduction will not contribute at all at the high end of the frequency range, because its loss ${\epsilon}_{c}^{\u2033}$ drops off with 1/ω. At the end we may therefore also fill in the missing ${\epsilon}_{kk}^{\u2033}$ with the measured ε″ values and thus skip the use of the special end frames.

In Steeman and van Turnhout (1997) we have also carried out the KK conversion of ε′ to ε″ by solving iteratively a triangular set of ε″ terms. It will be clear that the one-liners of Equations (3.7) and (3.8) are much easier to implement in the data analysis.

Clearly, the calculation of conversion frames by kernel matching can easily be extended to other conversions. Such as the conversion of a few ε″ values to *dε*′∕*dlnω*. a mixed conversion of ε″ and ε′ to *dε*′∕*dlnω* is possible, possible, see Figure 9.

**Figure 9. A few options to obtain dε′∕dlnω via KK conversion**. Any contribution of possible conduction loss should of course be avoided if we recruit ε″ data. This can be effected by implying ε″

_{cf}values. This special option is obtained by invoking a constraint. ε′ midway refers to sampling of the ε′ data at half the interval, cf. Equation (7.1). It produces a higher accuracy than the central discretization of ε′ with Equation (3.11).

The logarithmic derivative of ε′ is a compelling quantity, because it has a higher resolution of nearby peaks than ε″. However, we should be careful not to infect the resulting *dε*′∕*dlnω* by ohmic conduction from ε″. This can easily be prevented by imposing a constraint on the linear l.s.q. minimization of the *dε*′∕*dlnω* kernel. We therefore matched

with as constraint $\sum _{k=-1}^{1}{c}_{k}({2}^{-k}\u2215x+{2}^{k}x)=0$. This leads to the following mixed conversion

This mixed conversion illustrates the versatility of kernel matching. Yet it is still compatible with the KK relations. Apparently, we only need 2 ε′ and 3 ε″ values to get an accurate estimate of *dε*′∕*dlnω*.

Admittedly, we can use as alternative to Equation (3.9) the 5 term central logarithmic difference of ε′. This gives:

with d_{k} = {−1, 8, 0, −8, 1}/(12lnh). This numerical derivative, which derives from a 3rd degree logarithmic polynomial, is for *h* = 2 less accurate than (3.10) for sharp Debye like peaks.

The accuracy obtained with 3 options for calculating *dε*′∕*dlnω* is demonstrated in Figure 10. The central difference frame performs less than the mixed frame of Equation (3.9) and also less than the “halfway” frame of Equation (7.1). An interval of *h* = 2 roughly equals 10^{0.3}, taking *h* = 10^{0.2} improves the accuracy of the central difference markedly. It then scores virtually as good as the other two.

**Figure 10. Accuracy of 3 options for computing dε′∕dlnω for a Debye relaxation**. The commonly used central frame of Equation (3.11) is evidently not the best near the maximum for

*h*= 2. The offset at the maximum disappears if we reduce h to 10

^{0.2}.

## Computing Conversion Frames with Symbolic Differential Operators

There is a saying that one can do “integration by differentiation,” this also applies to the KK relations. The KK integrals are in fact logarithmic convolution integrals. This can be shown be rewriting Equation (2.1) to

in which ε′(*uω*_{o}) now depends on the product of u and ω_{o}. Integrals with such a function are called convolution integrals. They can be evaluated in a special way (Hirschman and Widder, 1955).

For pursuing this, we introduce the logarithmic derivative D_{l}. This allows us to move the x-position of a function f to hx using the symbolic exponential operation

This is the logarithmic variant of a better known operation with the normal derivative D, which produces a linear shift

Equation (4.2) can be proven be developing both sides in a Taylor expansion around *h* = 1.

The exponential logarithmic derivative operation of Equation (4.2) can be replaced by a logarithmic shift E_{l}

which produces the equality

The action of both logarithmic operators for transforming f(x) into f(hx) is portrayed in Figure 11 on the left.

**Figure 11. On the left, two symbolic tools to transform f(x) to f(hx)**. On the right, trick to link ε′(uω_{o}) symbolically to ε′(ω_{o}) and thus push the latter out of the convolution integral of KK in Equation (4.8).

It may be illuminating to use the symbolic equality of Equation (4.5) to find the panel of coefficients for the 1st derivative of ε′: D_{l}ε′ = *dε*′∕*dlnω*. Expanding lnE_{l} in a series of up to the 4th order around E_{l} = 1 results in

with the coefficients d_{k} equaling those given in Equation (3.11). Note that we pre-multiplied with ${E}_{l}^{2}$ in order to get the coefficients of the central difference. Series expansions like that of Equation (4.6) are easy to do in symbolic programs like Maple and Mathematica.

With the help of e Equation (4.2), as sketched on the right in Figure 11, we may replace ε′(*uω*_{o}) by

inserting this in Equation (4.1) gives

in other words we can push ε′(ω_{o})−ε_{∞} out of the KK-integral. The resulting integral with u as variable is known analytically. We thus get

We speak of a *symbolic* differential operator. In a similar way we find as symbolic differential operator for the ε″ to ε′ conversion

We have thus indeed succeeded in replacing the KK integrals by KK differential operators.

Admittedly, the KK relations can be applied for the analysis of dispersions in many areas. One of these areas are ultrasonic relaxation studies, where their differential forms did not remain unnoticed (Waters et al., 2003).

The use of these unusual symbolic operators looks at first sight not easy. The job was facilitated a lot by inserting for D_{l} = lnE_{l}/lnh in the tan and cot operator. Let us illustrate this for Equation (4.10), see also Figure 12. If we take the logarithmic derivative of both sides, we get

**Figure 12. Inserting D _{l} = lnE_{l}/ln2 in a symbolic operator provides a versatile tool to get an expansion in E_{l}, which in turn provides an easy-to-use conversion frame, in this case for ε″ to dε′∕dlnω**.

In order to arrive at a conversion frame for 5 logarithmically spaced ε″-data with *h* = 2, we expand Equation (4.11) as follows around E_{l} = 1

which yields for which yields for *dε*′∕*dlnω*,

We did a similar expansion for only 3 ε″ terms. The *dε*′∕*dlnω* calculated with neither the 3 nor the 5 terms was satisfactorily accurate for a Debye relaxation (one produced an undershoot, the other an overshoot). We therefore took the average of the coefficients. This average results in

with these coefficients the *dε*′∕*dlnω* calculated with Equation (4.13) turns out to be optimal.

We should stress that Equation (4.13) should *not* be used if the ε″-data contain ohmic conduction. But like kernel matching, the route via differential operators can be corrected for that. We will return to this in Section Simple Routes to a Higher Spectral Resolution: Logarithmic Derivatives or Differences of ε′ and ε″. We strongly recommend to rather use the frame from that section to convert measured ε″ data to *dε*′∕*dlnω*.

Being less ambitious, we could limit ourselves to using only just the very 1st term of the expansion of *D*_{l}cot(π*D*_{l}∕2)in D_{l} near 0. This equals 2/π, which gives as crude approximation

This approximation merely gives a reasonable estimate of *dε*′∕*dlnω* for broad relaxations. This modest zero-order estimate for *dε*′∕*dlnω* is furthermore not conduction free.

We mentioned in Section Introduction that Zahner (Germany) has incorporated a KK conversion in their software, viz. for calculating lnε_{a}(ω) from an integral of the phase θ(ω) (Schiller et al., 2001; Lasia, 2014). This conversion follows directly from Equation (4.10) by expanding the cot operator in D_{l} up to the 1st term

Like the inverse of D, 1/D, the inverse operator 1/Dl stands for integration, in this case logarithmic integration. Realizing that lnε_{a} corresponds to ε′ and θ to ε″ we get

The full integration of the θ values down to ω_{o} can be avoided if we calculate the logarithmic difference of lnε_{a}. A close approximation for this difference can be derived by series matching not the cot-operator itself, but the cot operator minus the integration operator

The pre-multiplication with 1/E_{l}-E_{l} is required because

We have done the series matching in a special way by matching to a sum in E${}_{l}^{2}$ rather than E_{l}. The reason being that this gives a more accurate approximation. This leads to

which gives *a*_{k} = {−1, 2, −1}π/(12lnh) and results for Δlnε_{a} with *h* = 2 in

By combining this with the logarithmic version of the Simpson rule we get

Checking Equation (4.21) for a Debye relaxation, we found that we could improve the results by multiplying the a_{k} coefficients with 0.94. This produces for *h* = 2 the following conversion

with *a*_{sk} = {−0.355035, 0.14709, 1.29843, 0.14709, −0.355035}. If we wish we can use Equation (4.22) recursively and obtain by starting from a given lnε_{a} value at high ω, values for lnε_{a} at successively lower ω's.

This new approximation for lnε_{a} can clearly also be used to convert ε″-data to Δε′(ω) ones

An interesting spin off of this conversion of ε″ to Δε′(ω) is that it allows us to calculate the ε′ response of models of which only ε″ is analytically known. Two typical examples are the Fuoss-Kirkwood and Jonscher models (see e.g., Tschoegl, 1989; Kremer and Schönhals, 2002).

There is hardly any need to use Equation (4.23) for converting measured ε″ data to Δε′. However, if we do, then we should realize that casual ohmic conduction may contaminate the Δε′ values, because the conversion frame is not conduction free. We mentioned above that by contrast the presence of ohmic conduction does not affect the proper conversion of θ to Δlnε_{a}.

By employing kernel matching with constrains we have found a conversion frame in which any conduction loss is canceled. This reads

with *a*_{fk} = {−0.485235, 0.465031, 0.899672, 0.465031, −0.485235}. Imposing the conduction free constraint causes the conversion to Δε′ to become a bit less accurate.

Surely, we could have found the frame for converting ε″ to Δε′values also directly from the cot operator itself, by simply matching

This gives a 3 term frame, we combined it with a 5 term frame to get optimal results for a Debye relaxation. This yields a combined frame with the following coefficients

*a*_{ck} = {−0.421343, 0.321676, 1.08188, 0.321676, −0.421343}. It turned out that the frame based on *a*_{sk} performs somewhat better to find Δε′ from ε″ than *a*_{ck}. Overall the panel of *a*_{sk} shows the highest accuracy for the KK-conversion of ε″ to Δε′.

Shtrauss has also discussed the use of conversion frames (Shtrauss, 2005, 2006). He derives the frames via the Mellin transform. This approach has some resemblance with the symbolic route. Shtrauss calculated his frames (which he calls functional filters) with l.s.q. If he would have used series expansion then he would have found the same coefficients as we have presented in our symbolic panels.

## Uncovering the l.f. Dispersion by Calculating the Conduction Free ε″ Losses

The removal of the contribution by ohmic loss to the l.f. dispersion has become an easy task with the availability of a fast KK conversion of ε′ to ε″. We just have to include the conversion frames of Table 1 in the data handling. They require just a few lines of code as Equation (3.8) shows. In particular the central part of the ${\epsilon}_{kk}^{\u2033}$ data is easily computed with Equation (3.7).

The ohmic conduction originates from the continuous flow of ions toward the electrodes. A flow that is driven by the applied a.c. field. In some colloidal systems percolation might happen at a certain critical concentration, when the conducting phase becomes co-continuous. This will show up in a sharp rise in the ohmic conduction. This can be monitored by calculating the conduction from the difference between the observed ε″ and the computed ${\epsilon}_{kk}^{\u2033}$, by using Equation (2.3) in reverse

Figure 13 displays the result of the computed conduction for a Debye relaxation with conduction loss present. The conduction is recovered closely up to a quite high frequency.

**Figure 13. One of the spin offs of the KK conversion, the retrieval of the conduction from the measured and the converted ε″ data**.

## Conversion without KK by Approximating ε′ and ε″ Data with a Complex Rational Function

In commercial software the KK conversion is avoided by assuming that the data can be modeled with a complex empirical function, ε′ and ε″ of which can each be specified explicitly. Usually the HN function is preferred for this purpose (Kremer and Schönhals, 2002). It reads

in which the last term accounts for ohmic conduction. The HN function has two shape or peak broadening parameters a and b. It can thus model the real and imaginary part of the experimental data of a relaxation phenomenon often successfully. Usually, a and b are assumed to remain below 1, but the proper constraints are 0 < a < 1 and 0 < ab < 1. The latter implies that b may exceed 1.

Imposing a specific model function is unsatisfactory, because it does not derive from the underlying process. This pertains the more so, if the system under study embodies a variety of relaxation phenomena. This is in fact often the case for colloidal systems (Grosse, 2002; Chassagne and Bedeaux, 2008; Delgado et al., 2014).

We have therefore followed a different model free approach and have approximated ε′ and ε″ with the real and imaginary part of a complex ratio in fractional power sums. We have also included a conduction term. The leads to the following expression for ε*

where c is a fractional power and ε${}_{r}^{*}$ represents the dispersion or relaxation part. $-\Im \left[{\epsilon}_{r}^{*}(\omega )\right]$ now plays the role of ε″_{kk}. An upper limit of *n* = 2 is often sufficient. The possibilities offered by approximating ε′ and ε″ data with a ratio of complex fractional polynomials in ω are indicated in Figure 14^{2}.

**Figure 14. A ratio of complex fractional sums provide a convenient means to model relaxation spectra**. Multiple peaks pose no problem and we get easy access to conduction free ε″ data and *dε*′∕*dlnω* without any involvement of the KK relations.

We could model ε′ and ε″ spectra with several relaxation peaks quite accurately with Equation (6.2). We usually started the l.s.q. minimization with *c* = 0.5 as initial value. The capability of Equation (6.2) to model multiple peaks in one go, is due to the fact that we do not need to include a relaxation time like we have to do for the common empirical functions like that of HN. They are meant to fit one peak only, whereby the relaxation time τ dictates the peak position on the scale (ω_{m} about 1/τ). We should take care to put in Equation (6.2) the degree in the degree in the denominator one higher than that of the numerator.

In order to recover both the relaxation part ε${}_{r}^{*}$ and the conduction term we have combined a two component input (the ε′ and ε″ data) with a two way fit [viz. to the real and imaginary part of Equation (6.2)] when the observed relaxation spectra show multiple peaks. We further advise to fit instead to the featureless ε′ or in addition to it the unstructured ε′, logarithmic difference Δε′(ω), defined by

Further details of this new joint multifunctional fit will be given in Section Improving the resolution by all-in-1 modeling of the real and imaginary data.

After this what we call all-in-1 modeling, ε″_{kk} simply follows from

The ε″_{kk} obtained from the fractional power approximation covers directly the whole ω range, no special calculations are needed at the start and the end of the range.

Having an analytical expression for ε*, we can also use Equation (6.2) to calculate its logarithmic derivative. This can be given in closed form, the real part of which leads to

where m is the upper limit in the numerator sum of ε*, and n that in the denominator sum. We may usually take n = m + 1, with m say *m* = 2. The use of the all-in-1 l.s.q. fit as the basis to obtain ε″_{kk} and *dε*′∕*dlnω*, has the advantage that these quantities become less prone to experimental errors.

The accuracy reached with Equation (6.5) and *m* = 2 and *n* = 3 for the merger of a double HN relaxation with ε_{∞} = 0.6, Δε_{1} = 0.75, a_{1} = 0.6, b_{1} = 0.7, τ_{1} = 1, Δε_{2} = 0.5, a_{2} = 0.9, b_{2} = 0.7, and τ_{2} = 7 is shown in Figure 15. We see a nice lining up with the exact *dε*′∕*dlnω* curve. This holds by the way also for the fits to the ε′ and ε″ data. It is interesting to note that the *dε*′∕*dlnω* curve with its better resolution indeed hints, albeit it vaguely, to the presence of two relaxations.

**Figure 15. After an all-in-1 rational fit dε′∕dlnω is directly analytically available from Equation (6.5)**. The accuracy of such a l.s.q. based

*dε*′∕

*dlnω*is clearly high. It coincides with the exact curve across a broad ω range.

We successfully tried some other model free approximations for ε*. These too allowed us to go around the KK conversion; they will be discussed elsewhere. The fractional power approximation for ε* of Equation (6.2) will be used in Section Improving the Resolution by All-In-1 Modeling of the Real and Imaginary Data to obtain a close analytical approximation to the distribution function of the underlying relaxation processes.

## Simple Routes to a Higher Spectral Resolution: Logarithmic Derivatives or Differences of ε′ and ε″

The removal of the ohmic conduction from ε″ is essential for getting insight in the low frequency part of the dielectric spectrum. But also the rest of the spectrum contains a lot of information about the diversity of processes possible. They occur alongside each other, and they often overlap.

We will obtain the best fingerprint if the various phenomena are resolved as good as possible. The most common way is to consider the ε″ spectra. For sure this demonstrates with its peaks much more structure than the monotonous ε′-spectra. But better options are at hand. The derivative *dε*′∕*dlnω* e.g., has a better resolution than ε″. However, an accurate numerical calculation of *dε*′∕*dlnω* from *h* = 2 spaced ε′-data is less easy than it seems for sharp peaks. Since most colloidal systems are liquids most of their relaxations are rather sharp (i.e., Debye like).

In addition to the options to compute *dε*′∕*dlnω* presented in Sections Computing Conversion Frames with Symbolic Differential Operators and Uncovering the l.f. Dispersion by Calculating the Conduction Free ε″ Losses, we like to draw attention to the use the midway central difference. This can be obtained by expanding

around E_{l} = 1, which gives d_{hk} = {−1, 27, −27, 1}/(24lnh). The higher accuracy of the midway logarithmic central difference for sharp peaks is due to the fact that it is based on ε′'s closer to the peak than the normal central difference Equation (4.6). That the halfway difference is more accurate has also been pointed out by Shtrauss (2006). The high accuracy of the half-spaced difference for *dε*′∕*dlnω* with *h* = 2 has already been shown in Figure 10.

The possible uncertainties of the numerical calculation of *dε*′∕*dlnω* from the observed ε′ data for sharp peaks brought us to use a well-defined alternative, the simple symmetric logarithmic difference Δε′ of ε′

If the input data are not available at a spacing of h = 2 then we can take more generally

As depicted in Figure 16, we can also opt for logarithmic differences of ε″. Those differences should show like Δε′ a peak for each relaxation time present (when ω_{m} = 1/τ). Another prerequisite is that they should be conduction free, so that they cancel any possible contribution from ohmic conduction in the experimental ε″ data. This gives us two options an asymmetric one with two terms and a symmetric one with three terms

**Figure 16. We will show that additional functions like the logarithmic difference of Δε′ and the conduction free ε″ difference Δε″ _{cf} will increase the resolution**. Interestingly, the handy Δε′ performs almost as good as

*dε*′∕

*dlnω*.

Here too we can use a spacing different from *h* = 2

It goes without saying that in the data analysis we should take h the same in all parts of the l.s.q. minimization. Hence, we should choose the same h in the sampling of the experimental data and in the theoretical model representation of Δε′ and Δε″. In order to simplify the minimization we best model Δε′ and Δε″ by taking:

The reason being that the analytical expression for ε* is often much simpler than that for ε′ and ε″. By using Equations (7.2.2–7.4.2) we let the l.s.q. routine compute the real and imaginary part of the differences, which is the most efficient way.

The increase in resolution by using the ε′ and ε″ differences is shown in Figure 17. Their peaks are like that of *dε*′∕*dlnω* much narrower than that of ε″ or ε″_{kk} for a Debye relaxation. However, they can be calculated much easier, while like *dε*′∕*dlnω* they have the added advantage that they too remove the conduction. It is further gratifying to notice that Δε′ enhances the resolution virtually as good as *dε*′∕*dlnω*.

**Figure 17. On the left peak shape of dε′∕dlnω compared to that of ε″**. On the right the peaks of the logarithmic differences Δε′ and Δε″

_{cf}, in particular the latter is much spikier than ε″

_{kk}. The curves are given for a Debye relaxation and

*h*= 2. The rise in ε″ in the low ω region is caused by ohmic conduction, it clearly lessens the resolution.

An enhanced resolution becomes of course more acute if the relaxation processes tend to overlap as often occurs in colloidal systems. Such a merger is also the case for the 2 HN relaxations considered in Figure 15. In Figure 18 we compare the ensuing spectra of Δε′ and Δε″_{cf} of these associated HN relaxations with the traditionally used spectrum of ε″_{h}, which at low frequencies also contains the additional conduction loss. The 2 underlying HN relaxations can in particularly be conceived in the Δε″_{cf} spectrum. The distinction in the Δε′ spectrum is less, but better than in the ε″-curve, which merely shows one united peak.

**Figure 18. Another illustration that Δε′ and Δε″ _{cf} increase the resolution of dielectric spectroscopy**. We simulated the different spectra of 2 nearby HN relaxations, which fuse together into a single ε″

_{h}peak. By contrast, the new differences do hint, even for these very close by relaxations, to the presence of 2 relaxations. This is in part due to the absence of the conduction loss in Δε′ and Δε″

_{cf}.

We can next make use of Δε″_{cf} to get an approximation for *dε*′∕*dlnω* from a frame of 5 ε″ values via the cot operator that is conduction free. We warned that this was not the case with the frame of Equation (4.13).

We recall from Section Computing Conversion Frames with Symbolic Differential Operators.

We further have for Δε″_{cf}

By replacing ε″(ω) by Δε″_{cf} in Equation (7.4) we get accordingly

By series matching the r.h.s of Equation (7.7) for one term and three terms in E_{l} we get as coefficients: {0.63662} and {0.183437, 0.269745, 0.183437}. We combined these coefficients into one set so that the optimal approximation for *dε*′∕*dlnω* of a Debye relaxation was achieved. This gives for the best coefficients of the conversion of Δε″_{cf} from Equation (7.4) at respectively ω/2,ω and 2 ω to *dε*′∕*dlnω*: {0.0984351, 0.43975, 0.0984351}.

We can now turn back directly to the measured ε″-data. This results in the following 5 term frame for the conduction free conversion of the *observed* ε″-data into the logarithmic derivative of ε′

with e_{fk} = {−0.19687, −0.38732, 1.80501, −0.38732, −0.19687}.

The accuracy of the special conduction free conversion of ε″ to *dε*′∕*dlnω* for a Debye relaxation can be seen in Figure 19. It behaves about equally well as the unconstrained conversion of Equation (4.13). However, it performs slightly less at the wings than the mixed conversion of ε′ and ε″ of Equation (3.10), which is also conduction free.

**Figure 19. The KK conversion allows the calculation of dε′∕dlnω from ε″ data**. The figure compares the plain, direct ε″ conversion of Equation (4.13) and the special one of Equation (7.8), by taking ε″ data of a Debye relaxation polluted by conduction.

## Enhancing the Resolution by Calculating Simple Approximations to the Distribution Function

If several concurrent relaxation processes are active, as is often the case in colloidal systems, then we can sort out the entire relaxation behavior comprehensively with a distribution function in relaxation times.

A Debye relaxation appears in such a distribution as a line spectrum. This is clear if we substitute a δ function in Equation (2.1) we then get

This shows that the distribution function offers us the highest resolution possible. A powerful method to obtain the distribution from a complex function like ε* was proposed as long ago as 1894 by Stieltjes in his correspondence with Hermite about his pioneering work on continued fractions. His inversion relies on inserting for w → i/τ. This transforms the kernel in Equation (3.2) into a singularity and leads owing to Cauchy's integral theorem to

where g(lnτ) is the so-called logarithmic distribution function τf(τ).

In view of Equation (8.2) we can now take up the complex rational fractional power approximation of Equation (6.2) and simply find

We can further make use of Equation (4.2) and the fact that lni = iπ/2 to derive in a simple way the symbolic differential operators for the inversion of ε′ and ε″. By expressing ε*(iω) in ε′ and ε″ we get:

This gives in view of Equation (8.2) for the symbolic inversion via ε′ and ε″

It may be worthwhile to point out that the sin and cos inversion operator and the tan and cot KK operator are compatible. In fact by dividing the sin and cos operations in Equation (8.5) for g(lnτ) we just get our KK tan operator of Equation (4.9).

By using the product or Mittag-Leffler approximations in D_{l} for the sin and cos operator, we can show that crude approximations to g(lnτ) can be obtained from:

with ω → 1/τ (see e.g., Tschoegl, 1989). If we discretize the derivatives in Equation (8.6), then we obtain as differences

We have manipulated the discretization of the derivatives of ε″ somewhat so that with just 2 and 3 terms full nullification for ohmic conduction was achieved.

We used these logarithmic differences already in Section Simple Routes to a Higher Spectral Resolution: Logarithmic Derivatives or Differences of ε′ and ε″. It is now no longer surprising that we found them to promote the spectral resolution, as we have already demonstrated in Figures 17, 18. The mix of ε″+ *dε*″/dlnω was proposed by Kaatze (2003), as a means to remove the ohmic conduction. In fact ε″(ω)−*d*^{2}ε″∕*d*ln^{2}ω is a better choice, not only because it is symmetric, but also because it has a higher resolution. Obviously, our differences of ε″ are much easier calculated than the combination of ε″ with its derivatives, while their performance almost matches that of the differential expressions.

The performance of the difference approximations for decomposing the overlapping HN relaxations used earlier in Figure 15 compared to that of the distribution obtained with Equation (8.3) is shown in Figure 20. Since Stieltjes inversion is based on ω → i/τ, we have rather plotted g(lnα), where α = 1/τ. In this way we can use the same scale for all spectra.

**Figure 20. Evidence of the strong resolution power of the distribution function**. Using Equation (8.3) we computed the distribution from a rational all-in-1 modeling to ε′and ε″of the 2 overlapping HN relaxations employed before in Figures 15, 18. The underpinning relaxations can also be distinguished in the spectra of the differences Δε′and Δε″_{cf}, this is in particular true for the latter.

Clearly the exact complex Stieltjes inversion performs best. The curve of the distribution function reveals convincingly the 2 HN relaxations present. The next best resolution is provided by Δε″_{cf}. The simple difference approximations on the other hand have the advantage that they can be applied directly to the measured data, and thus readily provide straight insight via plots in the relaxation behavior under study. The location of the peaks in the plots supply good estimates for the various relaxation times that are active. These estimates, which often lie decades apart, can be used as initial values in the final mathematical modeling.

## Improving the Resolution by All-in-1 Modeling of the Real and Imaginary Data

It is often illuminating if overlapping relaxations can be separated visually in plots. We showed in Section Enhancing the Resolution by Calculating Simple Approximations to the Distribution Function that this can be achieved by plotting approximations to the distribution function like Δε′ and Δε″_{cf} or even more so by plotting the distribution function itself by invoking the complex Stieltjes inversion.

However, we should realize that appropriate *fitting* offers a much more powerful tool to separate nearby relaxations. The common approach is to model a multiple relaxation with a sum of HN functions. By modifying Equation (6.1) to a sum, we have

Even a limited sum of 2 HN functions and one conduction term leads to no less than 10 unknowns. This hampers a proper mathematical separation of the 2 processes, the more so because the l.s.q. minimization is a nonlinear one.

The usual approach is to minimize the deviations between the measured and the HN-model values of ε′ and ε″ in one sum as follows

However, a better choice is to split the sum in two and minimize

This apart-together or all-in-1 fitting can be achieved by making use of a mathematical switch, which combines the proper data, say ε′, with the proper model i.e., ${\epsilon}_{hn}^{\prime}=\Re \left[{\epsilon}_{hn}^{*}\right]$, etc. Such a coupled minimization assures that the HN parameters in ε′_{hn} and ε″_{hn} are not allowed to differ, which clearly is a necessity. Figure 21 outlines the basics of the all-in-1modeling.

**Figure 21. Scheme of the all-in-1 or apart-together l.s.q. modeling**. The strategy is to link the mix of input data to the mix of models via an extra dummy variable, in fact an index, that acts as a switch and thus takes care of the correct link.

Most commercial software, such as Maple, Mathematica, Matlab, have a nonlinear l.s.q. procedure. We should preferably use the ones with the Levenberg-Marquardt routine built-in. They should also allow for the minimization of more than one variable. A typical example is the routine FindFit of Mathematica, in which the switching can be implemented with several commands. One of the switches we use is based on the conditional if statement, which can be applied in other software as well. Another one is e.g., Kronecker's delta.

FindFit has the following structure FindFit [input data, model, unknowns, variables]. It gives the values of the unknowns as output. Assuming that we do the all-in-1 modeling for ε′ and ε″, with two HN's and one conduction term, we have

We should then activate FindFit as follows,

input data: Join[Table[{ω_{k}, 1, ε′(ω_{k})}, {k, n}], Table[{ω_{k}, 2, ε″(ω_{k})}, {k, n}].

model: If[i = 1, ε′_{hn}(ω), 0] + If[i = 2, ε″_{hn}(ω), 0].

unknowns: ε_{∞}, Δε_{1}, a_{1}, b_{1}, τ_{1}, Δε_{2}, a_{2}, b_{2}, τ_{2} and γ.

variables: ω, i.

We have thus added to ω an extra variable i. This is an index that controls the if switch and makes the proper choice between the input data and the ε′ and ε″ models, for i = 1 it selects the ε′ data and links it to the ε′-model, while for i = 2 it selects and combines the ε″-data and the ε″-model. This is the reason that the input data comprise along with ω_{k} and the function values for ε′(ω_{k}) and ε″(ω_{k}), the index i as 2nd variable. In doing so we have transformed the common univariate minimization in ω, into a quasi-bivariate (or two-dimensional) minimization in ω and i. Although ε′_{hn} and ε″_{hn} can be written out analytically, we better model them as suggested earlier with

in order to keep the minimization simple.

The all-in-1 modeling can of course be extended to more input data and more model functions. In particular data and models for Δε′ and Δε″ are attractive, because they have intrinsically a better resolution power. They are therefore more powerful in separating well the genuine colloid relaxations from that of the electrode polarization. Both have the additional advantage over ε″ that they are not affected by ohmic conduction. Furthermore, the models for both are simple for e.g., Δε′ we frankly have

Figure 22 illustrates what can be achieved with an all-in-1 modeling of 2 overlapping HN relaxations, using Equation (9.4). Visually it seems that we are only dealing with 1 relaxation. However, by using not only ε′ and ε″, but Δε′ and Δε″_{cf} as input as well, we could unravel the 2 underlying HN's nicely, despite the fact that we have contaminated the data in this simulation with random relative errors of ±1%. We can see on the r.h.s. of Figure 22 that the recalculated ε″_{h1} and ε″_{h2} relaxation peaks closely resemble the original exact curves. The data were generated for: ε_{∞} = 0.6, Δε_{1} = 0.75, a_{1} = 0.6, b_{1} = 0.7, τ_{1} = 1, Δε_{2} = 0.5, a_{2} = 0.9, b_{2} = 0.7, τ_{2} = 7, and γ = 0.003, with the all-in-1 modeling we recovered these values as: ε_{∞} = 0.60, Δε_{1} = 0.75, a_{1} = 0.58, b_{1} = 0.74,τ_{1} = 1, Δε_{2} = 0.50, a_{2} = 0.91, b_{2} = 0.66, τ_{2} = 7.22, and γ = 0.003.

**Figure 22. Decomposing 2 overlapping HN relaxations with all-in-1 modeling**. Despite the fact that the 2 HN relaxations almost blend into one we could recover, as the results on the right show, each of the two genuine HN relaxations comfortably. Even though the simulated data were disturbed with random errors of ±1%.

The HN function is just an empirical model function. A much better way to characterize the various dielectric properties of colloidal systems is by making use of the complex dipolar coefficient β*.

This function is defined as follows, if we have a mixture of two components, say of spherical particles dispersed in a medium, then we can model according to Maxwell and Wagner the complex permittivity of this mixture ${\epsilon}_{m}^{*}$ by Grosse (2002) and Steeman and van Turnhout (2002)

where ${\epsilon}_{b}^{*}$ is the complex permittivity of the bulk, ${\epsilon}_{p}^{*}$ that of the particles, ϕ their volume fraction and β*(ω) the complex dipolar coefficient. We can rewrite Equation (9.6) to an explicit expression for ${\epsilon}_{m}^{*}$

In colloidal suspensions the bulk consists of the suspending electrolyte so ε${}_{b}^{*}=$ ε${}_{e}^{*}$. We see from Equation (9.6) that β*(ω) is determined by ${\epsilon}_{p}^{*}$ and ${\epsilon}_{e}^{*}$. The electrolyte is full of ions which are pushed by the alternating electric field to and fro the electrodes and the (insulating) particles. The ionic motion is however in part also controlled by diffusion. Near the electrodes the ions form electric double layers, which give rise to a strong EP. The motion of the ionic clouds near the particles give rise to special relaxations. The permittivity of the particles ε${}_{p}^{*}$ therefore is not simply due to the dipole relaxations and possible ohmic conduction. We should further realize that the particles may carry a charge (this depends on their ζ potential).

All in all it is quite a challenge to account for the many processes possible. Expressions for β*(ω) have e.g., been given by Grosse (2002) and Chassagne and Bedeaux (2008). The β*-model from Chassagne and Bedeaux (2008) was not fully explicit. Fortunately, it has been possible to achieve this lately, further details are given in Chassagne et al. (submitted).

Replacing ${\epsilon}_{m}^{*}$ by ${\epsilon}_{s}^{*}$, the permittivity of the suspension, and ${\epsilon}_{b}^{*}$ by ${\epsilon}_{e}^{*}$, the permittivity of the electrolyte, we get

where ${\epsilon}_{e}^{*}(\omega )={\epsilon}_{w}-i{\epsilon}_{w}\u2215(\omega {\tau}_{e})$, with τ_{e} = ε_{o}ε_{w}∕γ_{e}. However, if we measure the capacity of a cell with a colloidal suspension, then it includes a large contribution caused by the EP. This huge contribution manifests itself at low frequencies, because the ions move rather slowly.

The usual way to account for EP is to use a series circuit model as depicted in Figure 23. This consists of a large capacitor filled with electrolyte devoid of its ions. This represents the part formed by the two electric double layers near the electrodes. The DL capacitor is tied up in series with the complex bulk capacitor of the suspension, which is a mixture of the electrolyte and the colloidal particles. This leads to

where ${\epsilon}_{cs}^{*}$ denotes the complex permittivity of the suspension as measured with the cell, ε_{w} the (real) permittivity of water and ${\epsilon}_{s}^{*}$ the true complex permittivity of the suspension. We may neglect 2d_{l} against d in Equation (9.9), because the Debye length d_{l} is extremely small for aqueous systems. Note that d_{l} is known analytically, and can be calculated using handbook values e.g., from Weast (1987).

**Figure 23. Series model for the total complex capacitance of a colloidal suspension in a measuring cell**. The model on the right shows the true complex capacitance of the suspension ${C}_{s}^{*}$ in series with the DL capacitance, which is responsible for the EP. The circuit in the middle shows the series model for the electrolyte without the colloidal particles. The model on the left depicts the ion layers near the electrodes, when the top electrode is negative and the bottom one positive. We have also pointed to a continuous flow of ions, This ongoing flow is only possible if the electrodes are partially blocking. In our theoretical analysis, we assume that the electrodes are fully blocking, the EP is then at its strongest.

To simplify matters we suppose that the electrodes in the series model are totally blocking any charge exchange. We can generalize the model to partial blocking, the EP then becomes less dominant, but this will not be pursued here (see e.g., Coelho, 1983). Although it should be stressed that it is only for partial blocking that we will observe ohmic conduction losses.

We rather like to illustrate the use of Equations (9.9) and (9.10) in all-in-1 modeling. Obviously, we are mainly, not to say only, interested in the true response of the colloidal particles free of EP. This can again be achieved by using Δε′ and $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$ as input. They should be calculated from the ${\epsilon}_{cs}^{\prime}$ and ${\epsilon}_{cs}^{\u2033}$ values measured. However, we can model Δε′ and $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$ by taking straightaway ${\epsilon}_{s}^{*}$ from Equation (9.9). We need not to use the more complicated theoretical for Δε′ and $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$ from ${\epsilon}_{cs}^{*}$ measured with the EP included. This simplification is allowed in the ω region where the EP has died out. This ω region becomes readily visible from the Δε′ and Δε${}_{cf}^{{}^{\u2033}}$ plots. So in the apart-together model with built-in switch we take:

where ${\epsilon}_{s}^{*}$ follows from Equation (9.9) wherein ε${}_{e}^{*}$ can be specified for most electrolytes using the conductivities given in a handbook like (Weast, 1987). Instead of $\Delta {\epsilon}_{fc}^{\u2033}$ we can also use the symmetric difference Δε″_{cf}. The basic parameters in β* are the zeta potential and the particle radius a_{o} (see Chassagne et al., submitted). It are these important unknowns that will come out of the all-in-1 fit.

An unique aspect of the β* function over the HN fit is that it allows the fitting of multiple relaxation processes with just *one* function, provided that it is explicit and that in the β* model all underlying relaxations are properly accounted for Chassagne et al. (submitted).

The real and imaginary part of the dipolar coefficient of 50 nm particles with a zeta-potential of 4 dispersed in an 1:1 electrolyte is shown in Figure 24. We generated these curves by taking the same parameter values for the electrolyte as mentioned in Chassagne et al. (submitted) Interestingly, we could retrieve the given particle size and zeta-potential accurately by performing the all-in-1 fit as outlined above. The recalculated β′ and β″ curves coincide nicely with the exact ones across the whole ω range. We found that this remains so, if we include random errors in the input data.

**Figure 24. The β′and β″ spectra for a colloidal mixture of 50 nm sized spherical particles with ε _{p} = 2 carrying a charge equivalent to a ζ potential of 4 suspended in an 1:1 electrolyte**.

The dipolar coefficient represented in Figure 24 actually gives rise to two distinct relaxation phenomena at high frequencies. These relaxations do not show up in the ε″ curve, because the ε″ loss is dominated and controlled by that of the electrolyte and so is very high. However, they do emerge brightly in the logarithmic differences of ε′ and ε″.

This is demonstrated in Figure 25, which shows on the left Δε′ and on the right $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$ for the same colloidal suspension with ζ = 4 and particle radius a_{o} = 25 nm. They unravel that such a suspension shows right after the EP relaxation, two clear relaxation peaks due to the particles. Both could be nicely fitted with one and the same β* model using Δε′ and $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$. Not surprisingly, we could therefore ascertain the parameter values used in the simulated data quite accurately. It may be worth pointing out that the steep rise at low frequencies in Δε′_{cs} and Δε${}_{cf,cs}^{{}^{\u2033}}$, which quantities both derive from the measured data, is brought about by the EP in the measuring cell. Comparing the Δε′_{cs} and Δε${}_{cf,cs}^{{}^{\u2033}}$ curves we see that this rise starts earlier, i.e., at higher ω's, for $\Delta {\epsilon}_{cs}^{\prime}$. This means that even in the differences of the measured cell data the genuine relaxation peaks come out quite reliably. This is specially apparent in the Δε${}_{cf,cs}^{{}^{\u2033}}$ spectrum. It is by virtue of this close overlap between Δε′_{s} vs. Δε′_{cs} and Δε${}_{cf,s}^{{}^{\u2033}}$ vs. Δε${}_{cf,cs}^{{}^{\u2033}}$ that the all-in-1 modeling is so successful.

**Figure 25. Plots of Δε′, Δε${}_{\mathrm{\text{cf}}}^{{}^{\u2033}}$, and ε″ for the measured data (labeled with subscript cs) and the true values of the suspension itself (labeled with s)**. The ε″ loss is very high, so we have rather plotted 10^{−6}ε″. The data were generated with the β* model of Figure 24. The ε″ curves do not reveal the factual colloidal relaxations present. By contrast the difference curves do so indisputably, in particular the Δε${}_{\text{cf}}^{{}^{\u2033}}$ curves exhibit two prominent peaks.

Clearly, we would have found a very similar close agreement between the measured and true curves of the logarithmic derivative *dε*′∕*dlnω* beyond the EP region. The possibility to remove with a more advanced data analysis the EP-contribution was earlier established and explored by Jimenez et al. (2002). Instead of *dε*′∕*dlnω* we now rather advocate and prefer the use of Δε′ and $\Delta {\epsilon}_{cf}^{{}^{\u2033}}$, because these well-defined quantities are much easier to calculate from the data and to model. It should be stressed that actually an analytical model for *dε*′/dlnω of Equation (9.9) cannot be given, because β*(ω) can only only be represented by a very complicated expression, cf. Chassagne et al. (submitted)^{3}.

## Discussion

The use of the singular KK integrals for the interconversion between the real and imaginary part is still limited. They are mainly used in impedance spectroscopy to discover spurious effects in the data analysis caused by measuring errors (Barsoukov and Macdonald, 2005; Orazem and Tribollet, 2008).

By taking up the real part sufficiency based on the KK relations we have managed to convert ε′ into ε″ data. In this way conduction free ε″ losses can be calculated that allow us to uncover the genuine l.f. colloidal dispersions that otherwise remain hidden below strong conduction losses. This dissipative energy loss originates from the continuous flow of ions in the electrolyte toward the electrodes. It has nothing to do with the true relaxations of the colloidal particles.

We have strived for a fast and easy KK conversion. This was achieved by replacing the KK integrals by conversion frames, which consist of a window of coefficients by which a set of logarithmically spaced ε′ data are multiplied. By moving such frames step by step along all ε′ data we get a new set of ε″ fully cleaned from ohmic conduction. The cleaned ε″ data can only come from the genuine relaxation losses. The best method to derive the easy to use conversion frames is kernel matching. This versatile method can also provide panels for other conversions such as from ε″ to *dε*′∕*dlnω*. We took care to do the latter in such a way that this conversion remains unaffected by possible conduction in the observed ε″-data. This requires a special frame of logarithmically spaced ε″-data.

As other source to generate conversion frames we explored symbolic differential operators. They were found by replacing the conventional KK integration by differentiation. This option is possible because the KK integrals can be seen as logarithmic convolution integrals. We invoked the logarithmic shift operator E_{l} to expand the cot operator in an apt way in conversion frames for *dε*′∕*dlnω*. We further showed that this operator plays a key role in Zahner's software for the logarithmic form of the KK conversion (Schiller et al., 2001; Lasia, 2014).

Since the dielectric response of colloidal systems often involves several processes (Grosse, 2002; Chassagne and Bedeaux, 2008; Delgado et al., 2014; Chassagne et al., submitted), we paid also attention to an improvement of the resolution power for the entire dielectric spectrum. Major gains in this respect could be made by using derivatives or differences in ε′ and ε″. The logarithmic differences made up of a few terms are the easiest to use, while they perform almost as good as the derivatives. Our approach has been to use differences based on approximations to the distribution function.

The ensuing distribution in fact outperforms all other methods and produces the highest resolution possible. We calculated the distribution to a high accuracy by making use of Stieltjes's complex inversion. For this purpose we approximated ε′ and ε″ by a complex rational fractional power function. This model free function has the additional advantage that it also allows a direct interconversion between ε′ and ε″. A KK conversion is no longer needed.

A better resolution shows up visually in plots vs. frequency. An even higher resolution of nearby relaxations can be achieved by mathematical modeling. A common approach is to describe each relaxation process present in an empirical way by a Cole-Cole function or Havriliak-Negami function. The latter is the most general. A sum of a few HN functions contains a lot of adjustable parameters. This makes the nonlinear l.s.q. minimization a formidable task.

The accuracy of the parameter estimation and thus of the reliability of the resolution could be improved a lot by developing a joint “apart-together” minimization of the real and imaginary parts. We have called this all-in-1 modeling and incorporated it in Mathematica's l.s.q. one liner routine FindFit. The trick is that we have built in a mathematical switch that links the proper input data to the proper part of the model function. The switch is controlled by an index i that is included as dummy variable in FindFit, which transforms it from an univariate ω-minimization into a pseudo bivariate ω,i-fit.

We successfully retrieved with the all-in-1 modeling those parameters that are the most crucial in the dielectric response of colloidal suspensions, viz. the charge and size of the particles. We could recover these vital unknowns by using the explicit β-model discussed in Chassagne et al. (submitted) for the dipolar coefficient of a suspension.

The ohmic conduction is not the only disturbance that poses a problem to the characterization of colloidal systems with DS. The electrode polarization may be an even larger thread to the unraveling of the true nature of the l.f. dispersion. The elimination of this hindrance by mere data fitting is touched upon in Section Improving the resolution by all-in-1 modeling of the real and imaginary data. Its removal is recently dealt with to a greater extent in Chassagne et al. (submitted) and van Turnhout et al. (2016), see also Ishai et al. (2013).

The KK relations can be applied to many areas (Peiponen et al., 1999; King, 2009; van Dalen et al., 2013). They have basically the same form in all areas. This means that the fast conversion tools based on conversion frames we have discussed may equally well be applied in the analysis of many other spectroscopic methods.

## Conclusions

By a ready-made panel-based conversion of ε′ to ε″ data the l.f. resolution of dielectric spectra of colloidal systems can be improved. The crux bringing about this improvement is that ε′ data will never contain any contribution of the ohmic conduction current, which is in phase with the applied a.c. voltage, whereas the current from the real or ε′ part is out phase.

A high resolution is also an issue when several relaxation processes occur in near unison. Simple avenues are proposed to enhance the resolution of nearby colloidal peaks. In addition to the use of 2 or 3 term differences of ε′ and ε″, a new way of modeling is advocated. The methods designed furthermore facilitate, by virtue of their resolving power, the proper correction of the data of colloids for electrode polarization.

## Author Contributions

The author devised and tested the new methods described in the work. He contrived its intellectual content. He approves the final version. He is accountable for all aspects of it in ensuring that questions related to the accuracy or integrity are appropriately investigated and resolved.

## Conflict of Interest Statement

The author declares 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

The author is very grateful to the former and present head of the Chemical Engineering Dept. Profs. E. J. R. Sudhölter and M. T. Kreutzer for the friendly hospitality provided to him. He is also obliged to Prof. E. J. R. Sudhölter (who heads the Sect. Organic Materials and Interfaces) for his useful comments. He further likes to thank the TU Delft Library for financial support to publish this article in open access. The remarks made by the referees are appreciated too.

## Abbreviations

cos(π*D*_{l}/2), cos-operator; cot(π*D*_{l}/2), cot-operator; sin(π*D*_{l}/2), sin-operator; tan(π*D*_{l}/2), tan-operator; d, electrode distance; d_{l}, Debye length; D, normal derivative; D_{l}, logarithmic derivative; E, normal shift; E_{l}, logarithmic shift; f(τ), distribution function; g(lnτ) = τf(τ), logarithmic distribution function; h, logarithmic spacing; i, imaginary unit; m*(ω), complex dielectric modulus; z*(ω), complex impedance; β*(ω), complex dipolar coefficient; ε*(ω), complex permittivity; ε_{a}(ω), magnitude or absolute value; ${\epsilon}_{cs}^{*}$ (ω), measured complex cell permittivity of suspension; ε_{dl}, real permittivity of double layer; ${\epsilon}_{e}^{*}$(ω), complex permittivity electrolyte; ${\epsilon}_{ep}^{*}$ (ω), complex permittivity electrode polarization; ε${}_{s}^{*}$(ω), true permittivity of suspension; ε_{w}, real permittivity of water; ε′, real part of permittivity; ε″, imaginary part of permittivity which may include conduction losses; ${\epsilon}_{cf-rf}^{\u2033}$, ε″ free of conduction via ratio of fractional sums; ε″_{kk} = ε″_{cf}, ε″ free of conduction via KK; γ, ohmic conductivity; γ_{e}, ohmic conductivity of electrolyte; κ′(x), Debye's real kernel; κ″(x), Debye's imaginary kernel; τ, relaxation time; τ_{e}, (ohmic) relaxation time electrolyte; *dε*′∕*dlnω*, logarithmic derivative of ε′; *dε*″∕*dlnω*, logarithmic derivative of ε″; Δε′(ω), logarithmic difference of ε′; $\Delta {\epsilon}_{fc}^{\u2033}$(ω), asymmetric conduction free ε″ difference; Δε″_{cf}(ω), symmetric conduction free ε″ difference.

## Footnotes

1. ^Although *dε*′∕*dlnω* is a negative quantity, because ε′ drops with ω, we found it more convenient to plot and treat it below as a positive quantity. So in most cases we will actually deal with the absolute value of d ε′∕*dlnω*.

2. ^We may envisage the use of fractional polynomials in iω as a modeling of the dielectric response with a parallel circuit of several Warburg impedances. By combining such polynomials in a ratio, we only need sums with a limited number of terms.

3. ^This means that under these circumstances the theoretical model function for *dε*′/*dlnω* has to be calculated via a logarithmic difference with an interval h → 1, say of h = 10^{0.05} or less.

## References

Barsoukov, E., and Macdonald, J. R. (eds.). (2005). *Impedance Spectroscopy. Theory, Experiments and Applications, 2nd Edn.* New York, NY: Wiley.

Bayard, M. (1935). Relations entre les parties réelles et imaginaires des impedances et determination des impedances en function de l'une des parties. *Rev. Gen. Electricité* 37, 659–664.

Bohren, C. F. (2010). What did Kramers and Kronig do and how did they do it? *Eur. J. Phys*. 31, 573–577. doi: 10.1088/0143-0807/31/3/014

Brachman, M. K., and Macdonald, J. R. (1954). Relaxation-time distribution functions and the Kramers-Kronig relations. *Physica* 20, 266–270.

Brather, A. (1979). Numerisch einfache Beziehungen zwischen Verlust- und Speicherkomponente der dynamischen dielektrischen Faktors. *Coll. Polym. Sci*. 257, 467–477.

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

Delgado, A. V., Rica, R. A., Arroyo, F. J., Ahualli, S., and Jiménez, M. L. (2014). “Ac electrokinetics in concentrated suspensions,” in *Colloid and Interface Chemistry for Nanotechnology*, eds P. Kralchevsky, R. Miller et al. (Boca Raton, FL: CRC Press), 313–350.

Gorter, C. J., and Kronig, R. de L. (1936). On the theory of absorption and dispersion in paramagnetic and dielectric media. *Physica* 3, 1009–1020.

Gross, B. (1943). On an integral transformation of general circuit theory. *Amer. Math. Monthly* 50, 90–93.

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

Hirschman, I. I., and Widder, D. V. (1955). *The Convolution Transform.* Princeton: Princeton Univ. Press.

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, 21. doi: 10.1088/0957-0233/24/10/102001

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. Interf. Sci.* 249, 327–335. doi: 10.1006/jcis.2001.8141

Kaatze, U. (2003). Logarithmic derivative complex permittivity spectrometry. *Meas. Sci. Technol.* 14, N55–N58. doi: 10.1088/0957-0233/14/9/401

Kramers, H. A. (1927). “La diffusion de la lumiere par les atomes,” in *Atti Congr. Intern. Fisici, Trans. Volta Centenary Congr., Como*, 545–557 (earlier unpublished communication to Royal Danish Academy in 1925).

Kronig, R. de L. (1926). On the theory of dispersions of x-rays. *J. Opt. Soc. Am. Rev. Sci. Instrum*. 12, 547–557 (in later publications he dropped de L. which stands for his mother's name de Laer, sometimes his name is misspelled as Krönig).

Lasia, A. (2014). *Electrochemical Impedance Spectroscopy and Its Applications*. New York, NY: Springer.

Ninomiya, K., and Ferry, J. D. (1959). Some approximate equations useful in the phenomenological treatment of linear viscoelastic data. *J. Colloid Sci*. 14, 36–48.

Peiponen, K.-E., Vartainen, E. M., and Asakura, T. (1999). *Dispersion, Complex Analysis and Optical Spectroscopy*. Berlin: Springer.

Schiller, C. A., Richter, F., Gülzow, E., and Wagner, N. (2001). Validation and evaluation of electrochemical impedance spectra of systems with states that change with time. *J. Phys. Chem. Phys*. 3, 374–378. doi: 10.1039/b007678n

Shtrauss, V. (2005). Digital signal processing for relaxation data conversion. *J. Non-Crystal. Solids* 351, 2911–2916. doi: 10.1016/j.jnoncrysol.2005.04.087

Shtrauss, V. (2006). FIR Kramers-Kronig transformers for relaxation data conversion. *Signal Processing* 86, 2887–2900. doi: 10.1016/j.sigpro.2005.12.005

Steeman, P. A. M., and van Turnhout, J. (1997). A numerical Kramers-Kronig transform for the calculation of dielectric losses free from Ohmic conduction. *Coll. Polym. Sci.* 275, 106–115.

Steeman, P. A. M., and van Turnhout, J. (2002). “Dielectric properties of inhomogeneous media,” in *Broadband Dielectric Spectroscopy*, eds F. Kremer and A. Schönhals (Berlin: Springer), 495–522.

Tschoegl, N. W. (1989). *The Phenomenological Theory of Linear Viscoelastic Behaviour*. Berlin: Springer.

van Dalen, K. N., Slob, E., and Schoemaker, C. (2013). Generalized minimum-phase relations for memory functions associated with wave phenomena. *Geophys. J. Int.* 195, 1620–1629. doi: 10.1093/gji/ggt297

van Turnhout, J. (2004). “The Kramers-Kronig relations: fast methods for their evaluation,” in *Lecture at 3rd Int. Conf. Broadband Dielectric Spectroscopy and its Applications.* Delft: Delft University of Technology.

van Turnhout, J., Sudhölter, E. J. R., Jiménez, M. L., and Delgado, A. V. (2016). Removal of electrode polarization of colloidal suspensions with differences and all-in-a modeling.

Waters, K. R., Hughes, M. S., Mobley, J., and Miller, J. G. (2003). Differential forms of the Kramers-Kronig dispersion relations. *IEEE Trans. Ultrason. Ferroelectr. Freq. Contr*. 50, 68–76. doi: 10.1109/TUFFC.2003.1176526

Keywords: all-in-1modeling, electrode polarization, KK conversion frames, logarithmic derivatives and differences, matching Debye kernels, multivariate apart-together fitting, spectral resolution, symbolic differential operators

Citation: van Turnhout J (2016) Better Resolved Low Frequency Dispersions by the Apt Use of Kramers-Kronig Relations, Differential Operators, and All-In-1 Modeling. *Front. Chem*. 4:22. doi: 10.3389/fchem.2016.00022

Received: 19 March 2016; Accepted: 25 April 2016;

Published: 11 May 2016.

Edited by:

Raúl A. Rica, ICFO - The Institute of Photonic Sciences, SpainReviewed by:

Hamit Yurtseven, Middle East Technical University, TurkeyFrancisco José Arroyo, University of Jaen, Spain

Copyright © 2016 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: J. van Turnhout, j.vanturnhout@tudelft.nl