Time domain characterization of the Cole-Cole dielectric model

Sverre Holm

doi:10.2478/joeb-2020-0015

Abstract

The Cole-Cole model for a dielectric is a generalization of the Debye relaxation model. The most familiar form is in the frequency domain and this manifests itself in a frequency dependent impedance. Dielectrics may also be characterized in the time domain by means of the current and charge responses to a voltage step, called response and relaxation functions respectively. For the Debye model they are both exponentials while in the Cole-Cole model they are expressed by a generalization of the exponential, the Mittag-Leffler function. Its asymptotes are just as interesting and correspond to the Curie-von Schweidler current response which is known from real-life capacitors and the Kohlrausch stretched exponential charge response.

Keywords:: Cole-Cole; Curie-von Schweidler; Kohlrausch; Constant phase element

Introduction

The Cole Cole model [1] is a generalization of the Debye dielectric relaxation model which fits measurements in many applications including the bioimpedance field, [2, Sec. 9.2.7]. One interpretation is that it represents a distribution of relaxation processes, each described by the Debye model. Since the Debye model has a simple time domain interpretation and both the current and charge responses to a voltage step are exponential, the Cole-Cole responses can therefore be expressed as sums of exponential functions. In practice, however, this result is often too complex to lend itself to interpretation.

In recent years, there has been a development in understanding of the responses of the Cole-Cole model found in a direct way. These results depend on the Mittag-Leffler function, a generalization of the exponential which is named after Gösta Mittag-Leffler (1846–1927). This function is rightly called the “queen function of fractional calculus” [3] showing the close link between non-integer derivatives and the Cole-Cole model. The asymptotes of the Mittag-Leffler function are just as important as the function itself and is what will be emphasized here.

There are two well-established results for non-ideal dielectrics. The first is that for a long time it has been known that the current response to a step voltage for a practical non-ideal capacitor often follows the Curie-von Schweidler power law:

I (t) \propto t^{- α'},

(1)

where α′ is an order which is defined after (19). In [4] such responses are measured for many practical capacitors and the law is attributed to Curie in 1889 and von Schweidler in 1907. This is especially relevant for non-ideal dielectrics.

The second is an even older result which is due to Kohlrausch who found that the discharge of a capacitor with glass as a dielectric medium in a Leiden jar follows a stretched exponential. The charge is:

Q (t) \propto exp [- {(t / τ *)}^{α}],

(2)

where τ* is a time constant and α is given after (19). In [5] this is traced back to 1854. The result was later rediscovered by Williams and Watts [6] and it is often called the Kohlrausch-Williams-Watt model.

The purpose of this paper is to increase awareness of the time domain properties of the Cole-Cole model by collecting and interpreting some results from in particular [5, 7, 8]. The paper starts with the Debye model in order to define the relevant current and charge responses, called the response function and the relaxation function respectively. It will also be shown that both the Curie-von Schweidler power law and the Kohlrausch-Williams-Watt stretched exponential response are approximations to those of the Cole-Cole model. Finally, it is also shown that just as the Debye model corresponds to an ordinary partial differential equation for the constitutive law between the displacement field and the electric field, the Cole-Cole model corresponds to a similar equation but with non-integer, i.e. fractional derivatives.

Definitions

The constitutive relation between the displacement field, D, and the electric field, E is:

D = ε_{0} ε_{r} E = ε_{0} E + χ E = ε_{0} E + P,

(3)

where P is the polarization charge density, ɛ₀ is free space permittivity, ɛ_r is relative permittivity, and ξ is susceptibility. In the time domain, D(t), represents a charge density.

Frequency-dependency can be given either for the susceptibility [5, 9] or for the permittivity [2]. The relationship between the two is:

χ = ε_{0} (ε_{r} - 1) .

(4)

There are two reasons why we consider the permittivity here. First, in the bioimpedance field ɛ_r ≫ 1 so there is little difference in practice and D ≈ P. Second and more important, it is ɛ_r which is directly reflected in the properties of the macroscopic capacitance of the medium, and since impedance or capacitance is what can be measured, it makes sense to specify properties in terms of ɛ_r. This is also how models are justified in the bioimpedance field as [2, Sec. 3.1.3] says: “Polarization P cannot itself be measured. Dielectric theory is therefore invariably linked with the concept of a capacitor formed by two plates with the dielectric in between.”

This capacitance of such a dielectric material is

C = \frac{A}{d} ε_{0} ε_{r},

(5)

where A and d are the area and the plate distance of the capacitor. The capacitance is complex for the models considered here and its impedance is Z = (jωC)⁻¹. The frequency domain response to an input voltage is:

I (ω) = \frac{U (ω)}{Z (ω)} = j ω \frac{ε_{0} A}{d} ε_{r} (ω) U (ω) .

(6)

When the input voltage is a step, U(ω) = (jω)⁻¹, this is:

I_{step} (ω) = \frac{ε_{0} A}{d} ε_{r} (ω) .

(7)

In the time domain, the current step response is found as the inverse Fourier transform of the relative permittivity.

The current charge relation is :

I = A \cdot J = A \frac{d Q}{d t},

(8)

where J(t) is the charge density. Charge is therefore found by an integration of the result for the current plus a constant. Integration is equivalent to division by jω in the frequency domain and therefore the charge response to a step in voltage is related to

Q_{step} (ω) = \frac{I_{step} (ω)}{j ω A} = \frac{ε_{0}}{d} \frac{ε_{r} (ω)}{j ω} .

(9)

Debye model

As a reference and in order to establish terminology, the Debye model will first be analyzed for its current and charge responses. Its permittivity is

ε_{r} (ω) = ε_{\infty} + \frac{ε_{s} - ε_{\infty}}{1 + j ω τ},

(10)

where ɛ_s is the static value and ɛ_∞ < ɛ_s is the value at infinity frequency, τ is a characteristic time constant for the medium. The first term, represented by the constant ɛ_∞, represents an ideal capacitor which is in parallel with a frequency-varying part.

Time-domain characterization

The current response, (7), is

\begin{matrix} I_{step} (t) & = \frac{ε_{0} A}{d} ℱ^{- 1} {ε_{r} (ω)} \\ = \frac{ε_{0} A}{d} (ε_{\infty} δ (t) + \frac{ε_{s} - ε_{\infty}}{τ} e^{- t / τ}) . \end{matrix}

(11)

The current has an initial impulse due to the charging of an ideal capacitor followed by a current that dies out with a time constant τ.

Likewise the charge response is given by (9) or by integration of the current:

Q_{step} (t) = \frac{1}{A} \int_{0}^{t} I_{step} (u) du .

(12)

This gives

\begin{matrix} Q_{step} (t) & = \frac{ε_{0}}{d} (ε_{\infty} - (ε_{s} - ε_{\infty}) e^{- t / τ} + K) \\ = \frac{ε_{0}}{d} (ε_{s} - (ε_{s} - ε_{\infty}) e^{- t / τ}), t \geq 0, \end{matrix}

(13)

where K is a constant which is such that the initial value for the charge is proportional to ɛ_∞. The charge therefore starts with this value and ends up to be proportional to ɛ_s for large time. This is in agreement with the example in [2, Sec. 3.4.2].

Characterization of general models

The models will be given in terms of a normalized permittivity which for the Debye model is:

ε_{D} (ω) = \frac{ε_{r} (ω) - ε_{\infty}}{ε_{s} - ε_{\infty}} = \frac{1}{1 + j ω τ} .

(14)

In order to characterize subsequent models, the two descriptions of [5] will be used. The first is the response function, φ(t), which characterizes the current response. It is given as the inverse Fourier transform of the normalized relative permittivity, ɛ(ω). In the Debye example this is

ϕ (t) = ℱ^{- 1} {ε (ω)} = ℱ^{- 1} {\frac{1}{1 + j ω τ},} = \frac{1}{τ} e^{- t / τ},

(15)

which can be recognized to be the main time-varying part of (11).

The second function is the relaxation function, ψ(t), which characterizes the charge. It is given by

ψ (t) = ℱ^{- 1} {\frac{1 - ε (ω)}{j ω}} = e^{- t / τ}, t \geq 0,

(16)

which is seen in (13) with a negative sign. It is only for the simple Debye model that the two functions are the same.

The definitions of φ(t) and ψ(t) are such that they both are non-negative and non-increasing functions of time, meaning that the zero order derivative is positive and the first order derivative is negative. This pattern of sign changes repeats infinitely with the second order derivative positive and so on. This is what characterizes a completely monotone function and it ensures that it can be expressed as a continuous distribution of exponential functions and that the Laplace transform is non-negative [5]. In this way the physical realizability of the considered system is guaranteed [10].

Constitutive law

The Debye model can also be expressed as a differential equation between D and E by combining (3) with a rearranged (10):

(1 + j ω τ) D (ω) = ε_{0} ε_{s} E (ω) + j ω τ ε_{0} ε_{\infty} E (ω) .

(17)

Inverse Fourier transformation then gives:

D (t) + τ \frac{d D (t)}{d t} = ε_{0} ε_{s} E (t) + τ ε_{0} ε_{\infty} \frac{d E (t)}{d t} .

(18)

This is the electrical equivalent of the standard linear solid or Zener model in linear viscoelasticity [7, 11].

Cole-Cole model

The permittivity of the Cole-Cole model follows a more general power-law than the Debye model:

ε_{r} (ω) = ε_{\infty} + \frac{ε_{s} - ε_{\infty}}{1 + {(j ω τ)}^{1 - α'}} = ε_{\infty} + \frac{ε_{s} - ε_{\infty}}{1 + {(j ω τ)}^{α}} .

(19)

In electromagnetics, the model is sometimes presented with an exponent of 1 − α′, so that α′ = 0 corresponds to the Debye model. It may also be expressed with an exponent α = 1 − α′ where 0 < α ≤ 1 and here that convention will be followed in order to conform to [5].

The permittivity of the Cole-Cole model in normalized form is:

ε_{CC} (ω) = \frac{ε_{r} (ω) - ε_{\infty}}{ε_{s} - ε_{\infty}} = \frac{1}{1 + {(j ω τ)}^{α}} .

(20)

Time domain characterization

Section 3.1 of [5] gives the functions for the Cole-Cole model. The response function that characterizes the current response is:

ϕ_{CC} (t) = ℱ^{- 1} {ε_{CC} (ω)} = \frac{1}{τ} {(t / τ)}^{α - 1} E_{α, α} (- {(t / τ)}^{α}),

(21)

and the relaxation function which describes the charge is:

ψ_{CC} (t) = ℱ^{- 1} {\frac{1}{j ω} - \frac{ε_{CC} (ω)}{j ω}} = E_{α} (- {(t / τ)}^{α}) .

(22)

The function E_α,β is the Mittag-Leffler function which is a generalization of the exponential function. The two-parameter Mittag-Leffler function is defined by

E_{α, β} (t) = \sum_{n = 0}^{\infty} \frac{t^{n}}{⌈ (α n + β)}, 0 < α \leq 1,

(23)

where ⌈(x) is the gamma function, a generalization of the factorial for non-integer arguments and where ⌈(n) = (n − 1)! for integer arguments. Setting β = 1 gives the standard Mittag-Leffler function E_α(t) = E_α,1(t) of (22). Another special case is E₁(t) which is the exponential function. In this article, it is in particular the asymptotes of the responses which are important.

Approximation of the response function

According to [5], the response function may be approximated:

ϕ_{CC} (t) \sim {\begin{matrix} \frac{1}{τ ⌈ (α)} {(t / τ)}^{α - 1}, & t ≪ τ \\ \frac{1}{τ ⌈ (- α)} {(t / τ)}^{- α - 1}, & t ≫ τ . \end{matrix}

(24)

The small time approximation corresponds to the Curie-von Schweidler law of (1) as mentioned in [12]. An example may also be found in [4] where the current in capacitors followed the Curie-von Schweidler law over days. On first sight, this seems inconsistent with the small time approximation above, but in fact it fits well. The argument is that the capacitors were modeled by a constant phase element. The Cole-Cole model approaches such an element for (ωτ)^α ≫ 1 and ɛ_∞ = 0 and then the capacitance is:

C \approx \frac{A}{d} \frac{ε_{0} ε_{s}}{{(j ω τ)}^{α}} .

(25)

An example in [4] is a polypropylene dielectric with α′ of (19) between 0.999 and 1, i.e. α between 0 and 0.001, where α = 0 is an ideal capacitor. The factor (ωα)^α requires a very large argument to be much larger than one for such a small α, e.g. ωτ needs to be 10¹⁰⁰ for α = 0.01 in order for the factor to reach a value of say ten. Therefore, even for large frequencies, τ has to be very much larger than some days. The point is that even when capacitors followed the Curie-von Schweidler law for as long as several days, this indicates that the small argument approximation of (24) was valid.

The exact expression and the approximations are plotted in Fig. 1 for α = 0.7 using numerical code from [13, 14]. The two approximations fit very well for small time and large time respectively.

Approximation of the relaxation function

The Mittag-Leffler function with a negative argument raised to a power can also be approximated [8, 5]:

ψ_{CC} (t) \sim {\begin{matrix} exp [\frac{- {(t / τ)}^{α}}{⌈ (α + 1)}], & t ≪ τ \\ \frac{{(t / τ)}^{- α}}{⌈ (1 - α)}, & t ≫ τ . \end{matrix}

(26)

The small time approximation is the stretched exponential or Kohlrausch-Williams-Watt function of (2). For large values of t the Mittag-Leffler function approaches a power law [8].The relaxation function along with both its approximations are plotted in Fig. 2. Both approximations fit very well.

Constitutive law

The Cole-Cole model can also be expressed as a differential equation between D and E. The frequency domain relation building on (19) is:

(1 + {(j ω τ)}^{α}) D (ω) = ε_{0} ε_{s} E (ω) + {(j ω τ)}^{α} ε_{0} ε_{\infty} E (ω) .

(27)

The property of the Fourier transform which was used in transforming from (17) to (18) is that (jω)ⁿ transforms to the nth order derivative for integer n. This property is generalized in fractional calculus to non-integer orders α. This demonstrates the close relationship between power-laws and fractional calculus [11]. This results in this constitutive equation for the Cole-Cole model:

D (t) + τ^{α} \frac{d^{α} D (t)}{d t^{α}} = ε_{0} ε_{s} E (t) + τ^{α} ε_{0} ε_{\infty} \frac{d^{α} E (t)}{d t^{α}} .

(28)

This is the electrical equivalent of the fractional Zener model in linear viscoelasticity [7]. Fractional calculus has become an important research area in recent decades both in mathematics and in physics, despite its roots long ago [15]. One feature is that non-integer derivatives have memory. The details are beyond the scope of this paper, but hopefully its role in describing the Cole-Cole model may be a motivation for delving into it in e.g. the two books just cited.

Beyond the Cole-Cole model

There are several alternatives to the Cole-Cole model such as the Cole-Davidson and Havriliak-Negami models. The latter is the more general one:

\tilde{ε} (ω) = ε_{\infty} + \frac{ε_{s} - ε_{\infty}}{{(1 + {(j ω τ)}^{α})}^{β}}, 0 < α \leq 1, 0 < β < 1,

(29)

where β = 1 gives the Cole-Cole model and α = 1 gives the Cole-Davidson model. All three models yield an ideal capacitor in parallel with the constant phase element of (25) in the limit of a large τ. These models are also analyzed in [5] and the main thing to note is that the Curie-von Schweidler law and the Kohlrausch-Williams-Watt function fit the asymptotes of these models just as well as they fit the Cole-Cole model [12]. The link between the constant phase element, these early empirical results, and the Cole Cole model is therefore not unique.

Conclusion

The familiar frequency domain expression for the Cole-Cole model of order α can also be expressed in the time domain. The response function, which is related to the current response to a voltage step excitation, is expressed with a two-parameters Mittag-Leffler function. Its asymptote for small time is a power-law function which corresponds to the Curie-von Schweidler law. The relaxation function, which describes the charge response, is given by a one-parameter Mittag-Leffler function where the asymptote for small time is the stretched exponential or Kohlrausch-Williams-Watt function. For the Debye model, both these responses are given by exponential functions in time. The Debye model is also equivalent to a first-order differential equation between electric field, E, and the displacement field, D. This generalizes to a fractional differential equation with non-integer derivatives of order α for the Cole-Cole model.

Acknowledgment

I want to thank Thomas Holm for having read and commented the final manuscript.

Ethical approval

The conducted research is not related to either human or animal use.

Conflict of interest

The author states no conflict of interest.

References

1↑
Cole KS and Cole RH. Dispersion and absorption in dielectrics I. Alternating current characteristics. J. Chem. Phys. 1941; 9:341–51. DOI:
- Crossref
- Export Citation
2↑
Grimnes S and Martinsen ØG. Bioimpedance and bioelectricity basics, 3rd ed. Academic Press, 2015 :1–584. DOI:
- Crossref
- Export Citation
3↑
Mainardi F. Why the Mittag-Leffler function can be considered the queen function of the fractional calculus? Entropy 2020; 22:1359. DOI:
- Crossref
- Export Citation
4↑
Westerlund S. Dead matter has memory! Phys. Scr. 1991; 43:174. DOI:
- Crossref
- Export Citation
5↑
Garrappa R, Mainardi F, and Maione G. Models of dielectric relaxation based on completely monotone functions. Fract. Calc. Appl. Anal. 2016; 19:1105–60. DOI:
- Crossref
- Export Citation
6↑
Williams G and Watts DC. Non-symmetrical dielectric relaxation behaviour arising from a simple empirical decay function. Trans. Faraday Soc. 1970; 66:80–5. DOI:
- Crossref
- Export Citation
7↑
Mainardi F. Fractional calculus and waves in linear viscoelesticity: An introduction to mathematical models. London, UK: Imperial College Press, 2010 :1–347. DOI:
- Crossref
- Export Citation
8↑
Mainardi F. On some properties of the Mittag-Leffler function Eα(−t^α), completely monotone for t > 0 with 0 < α < 1. Discrete Cont Dyn-B 2014 :2267–78. DOI:
- Crossref
- Export Citation
9↑
Dissado L. Dielectric Response. Springer handbook of electronic and photonic materials. Springer, 2017. Chap. 10:219–45. DOI:
- Crossref
- Export Citation
10↑
Holm S and Holm MB. Restrictions on wave equations for passive media. J. Acoust. Soc. Am. 2017; 142. DOI:
- Crossref
- Export Citation
11↑
Holm S.Waves with power-law attenuation. Switzerland: Springer and ASA Press, 2019 :1–312. DOI:
- Crossref
- Export Citation
12↑
De la Fuente M, Jubindo MP, Solier J, and Tello M. The memory effect in dielectric response. J Phys C Solid State Phys 1985; 18:6547. DOI:
- Crossref
- Export Citation
13↑
Garrappa R. The Mittag-Leffler function, MATLAB Central File Exchange. 2015. Available from: https://www.mathworks.com/matlabcentral/fileexchange/48154-the-mittag-leffler-function [Accessed on: 2020 Dec 2]
14↑
Garrappa R. Numerical evaluation of two and three parameter Mittag-Leffler functions. SIAM J, Num. Anal. 2015; 53:1350–69. DOI:
- Crossref
- Export Citation
15↑
Podlubny I, Magin RL, and Trymorush I. Niels Henrik Abel and the birth of fractional calculus. Fract. Calc. Appl. Anal. 2017; 20:1068–75. DOI:
- Crossref
- Export Citation

If the inline PDF is not rendering correctly, you can download the PDF file here.

ProCite

RefWorks

Reference Manager

BibTeX

Zotero

EndNote

Article information

Received:: 2020-12-07

Published Online:: 2020-12-31

Citation Information:: Journal of Electrical Bioimpedance, Volume 11, Issue 1, Pages 101–105, eISSN 1891-5469, DOI: https://doi.org/10.2478/joeb-2020-0015.

Time domain characterization of the Cole-Cole dielectric model

Abstract

Introduction

Definitions