Physics:Cole–Cole equation

The Cole–Cole equation is a relaxation model that is often used to describe dielectric relaxation in polymers.

It is given by the equation

[math]\displaystyle{ \varepsilon^*(\omega) = \varepsilon_\infty + \frac{\varepsilon_s - \varepsilon_\infty}{1+(i\omega\tau)^{1 - \alpha}} }[/math]

where [math]\displaystyle{ \varepsilon^* }[/math] is the complex dielectric constant, [math]\displaystyle{ \varepsilon_s }[/math] and [math]\displaystyle{ \varepsilon_\infty }[/math] are the "static" and "infinite frequency" dielectric constants, [math]\displaystyle{ \omega }[/math] is the angular frequency and [math]\displaystyle{ \tau }[/math] is a dielectric relaxation time constant.

The exponent parameter [math]\displaystyle{ \alpha }[/math], which takes a value between 0 and 1, allows the description of different spectral shapes. When [math]\displaystyle{ \alpha=0 }[/math], the Cole-Cole model reduces to the Debye model. When [math]\displaystyle{ \alpha\gt 0 }[/math], the relaxation is stretched. That is, it extends over a wider range on a logarithmic [math]\displaystyle{ \omega }[/math] scale than Debye relaxation.

The separation of the complex dielectric constant [math]\displaystyle{ \varepsilon(\omega) }[/math] was reported in the original paper by Kenneth Stewart Cole and Robert Hugh Cole as follows:

[math]\displaystyle{ \varepsilon' = \varepsilon_\infty+ (\varepsilon_s-\varepsilon_\infty)\frac{1+(\omega\tau)^{1-\alpha} \sin\alpha\pi/2}{1+2(\omega\tau)^{1-\alpha}\sin\alpha\pi/2+(\omega\tau)^{2(1-\alpha)}} }[/math]

[math]\displaystyle{ \varepsilon''=\frac{(\varepsilon_s-\varepsilon_\infty)(\omega\tau)^{1-\alpha}\cos\alpha\pi/2}{1+2(\omega\tau)^{1-\alpha}\sin\alpha\pi/2+(\omega\tau)^{2(1-\alpha)}} }[/math]

Upon introduction of hyperbolic functions, the above expressions reduce to:

[math]\displaystyle{ \varepsilon'=\varepsilon_\infty+\frac{1}{2}(\varepsilon_0-\varepsilon_\infty)\left[1-\frac{\sinh((1-\alpha)x)}{\cosh((1-\alpha)x)+\sin(\alpha\pi/2)}\right] }[/math]

[math]\displaystyle{ \varepsilon''=\frac{1}{2}(\varepsilon_0-\varepsilon_\infty)\frac{\cos(\alpha\pi/2)}{\cosh((1-\alpha)x)+\sin(\alpha\pi/2)} }[/math]

Here [math]\displaystyle{ x = \ln(\omega\tau) }[/math].

These equations reduce to the Debye expression when [math]\displaystyle{ \alpha =0 }[/math].

The Cole-Cole equation's time domain current response corresponds to the Curie–von Schweidler law and the charge response corresponds to the stretched exponential function or the Kohlrausch–Williams–Watts (KWW) function, for small time arguments.

Cole–Cole relaxation constitutes a special case of Havriliak–Negami relaxation when the symmetry parameter [math]\displaystyle{ \beta=1 }[/math], that is, when the relaxation peaks are symmetric. Another special case of Havriliak–Negami relaxation where [math]\displaystyle{ \beta\lt 1 }[/math] and [math]\displaystyle{ \alpha=1 }[/math] is known as Cole–Davidson relaxation. For an abridged and updated review of anomalous dielectric relaxation in disordered systems, see Kalmykov.

See also

• Debye relaxation
• Cole–Davidson relaxation
• Havriliak–Negami relaxation
• Curie–von Schweidler law

Further reading

• "Dispersion and Absorption in Dielectrics II. Direct Current Characteristics". Journal of Chemical Physics 10 (2): 98–105. 1942. doi:10.1063/1.1723677. Bibcode: 1942JChPh..10...98C. 
• "Microscopic Models for Dielectric Relaxation in Disordered Systems". Physical Review E 70 (4): 041103. 2004. doi:10.1103/PhysRevE.70.041103. PMID 15600393. Bibcode: 2004PhRvE..70d1103K.