Havriliak–Negami relaxation

Template:Short description Template:Use American English The Havriliak–Negami relaxation is an empirical modification of the Debye relaxation model in electromagnetism. Unlike the Debye model, the Havriliak–Negami relaxation accounts for the asymmetry and broadness of the dielectric dispersion curve. The model was first used to describe the dielectric relaxation of some polymers, by adding two exponential parameters to the Debye equation:

$\hat{\epsilon }(\omega )={\epsilon }_{\mathrm{\infty }}+\frac{\mathrm{\Delta }\epsilon }{(1+(i\omega \tau {)}^{\alpha }{)}^{\beta }},$

where ${\epsilon }_{\mathrm{\infty }}$ is the permittivity at the high frequency limit, $\mathrm{\Delta }\epsilon ={\epsilon }_s-{\epsilon }_{\mathrm{\infty }}$ where ${\epsilon }_s$ is the static, low frequency permittivity, and $\tau$ is the characteristic relaxation time of the medium. The exponents $\alpha$ and $\beta$ describe the asymmetry and broadness of the corresponding spectra.

The relation is named after Stephen Havriliak Jr. and Sinichi Negami, who published in 1967.

Depending on application, the Fourier transform of the stretched exponential function can be a viable alternative that has one parameter less.

For $\beta =1$ the Havriliak–Negami equation reduces to the Cole–Cole equation, for $\alpha =1$ to the Cole–Davidson equation.

Mathematical properties

Real and imaginary parts

The storage part ${\epsilon }^{\prime }$ and the loss part ${\epsilon }^{″}$ of the permittivity (here: $\hat{\epsilon }(\omega )={\epsilon }^{\prime }(\omega )-i{\epsilon }^{″}(\omega )$ with $(\pm i{)}^2=-1$) can be calculated as

${\epsilon }^{\prime }(\omega )={\epsilon }_{\mathrm{\infty }}+\mathrm{\Delta }\epsilon {(1+2(\omega \tau {)}^{\alpha }\cos (\pi \alpha /2)+(\omega \tau {)}^{2\alpha })}^{-\beta /2}\cos (\beta \varphi )$

and

${\epsilon }^{″}(\omega )=\mathrm{\Delta }\epsilon {(1+2(\omega \tau {)}^{\alpha }\cos (\pi \alpha /2)+(\omega \tau {)}^{2\alpha })}^{-\beta /2}\sin (\beta \varphi )$

with

$\varphi =\arctan \left(\frac{(\omega \tau {)}^{\alpha }\sin (\pi \alpha /2)}{1+(\omega \tau {)}^{\alpha }\cos (\pi \alpha /2)}\right)$

Loss peak

The maximum of the loss part lies at

${\omega }_{\mathrm{m}\mathrm{a}\mathrm{x}}={\left(\frac{\sin \left(\frac{\pi \alpha }{2(\beta +1)}\right)}{\sin \left(\frac{\pi \alpha \beta }{2(\beta +1)}\right)}\right)}^{1/\alpha }{\tau }^{-1}$

Superposition of Lorentzians

The Havriliak–Negami relaxation can be expressed as a superposition of individual Debye relaxations

$\frac{\hat{\epsilon }(\omega )-{ϵ}_{\mathrm{\infty }}}{\mathrm{\Delta }\epsilon }={\int }_{-\mathrm{\infty }}^{\mathrm{\infty }}\frac{1}{1+i\omega {\tau }_D}g(\ln {\tau }_D)d\ln {\tau }_D$

with the real valued distribution function

$g(\ln {\tau }_D)=\frac{1}{\pi }\frac{({\tau }_D/\tau {)}^{\alpha \beta }\sin (\beta \theta )}{(({\tau }_D/\tau {)}^{2\alpha }+2({\tau }_D/\tau {)}^{\alpha }\cos (\pi \alpha )+1{)}^{\beta /2}}$

where

$\theta =\arctan \left(\frac{\sin (\pi \alpha )}{({\tau }_D/\tau {)}^{\alpha }+\cos (\pi \alpha )}\right)$

if the argument of the arctangent is positive, else

$\theta =\arctan \left(\frac{\sin (\pi \alpha )}{({\tau }_D/\tau {)}^{\alpha }+\cos (\pi \alpha )}\right)+\pi$

Template:Off topic Noteworthy, $g(\ln \tau )$ becomes imaginary valued for

$\frac{\hat{\epsilon }(\omega )-{ϵ}_{\mathrm{\infty }}}{\mathrm{\Delta }\epsilon }=\frac{(i\omega \tau {)}^{\alpha \beta }}{(1+(i\omega \tau {)}^{\alpha }{)}^{\beta }}$

and complex valued for

$\frac{\hat{\epsilon }(\omega )-{ϵ}_{\mathrm{\infty }}}{\mathrm{\Delta }\epsilon }=\frac{1}{(1-(\omega \tau {)}^{2\alpha }{)}^{\beta }}$

Logarithmic moments

The first logarithmic moment of this distribution, the average logarithmic relaxation time is

$⟨\ln {\tau }_D⟩=\ln \tau +\frac{\mathrm{\Psi }(\beta )+\mathrm{E}\mathrm{u}}{\alpha }$

where $\mathrm{\Psi }$ is the digamma function and $\mathrm{E}\mathrm{u}$ the Euler constant.

Inverse Fourier transform

The inverse Fourier transform of the Havriliak-Negami function (the corresponding time-domain relaxation function) can be numerically calculated. It can be shown that the series expansions involved are special cases of the Fox–Wright function. In particular, in the time-domain the corresponding of $\hat{\epsilon }(\omega )$ can be represented as

$X(t)={\epsilon }_{\mathrm{\infty }}\delta (t)+\frac{\mathrm{\Delta }\epsilon }{\tau }{\left(\frac{t}{\tau }\right)}^{\alpha \beta -1}{E}_{\alpha ,\alpha \beta }^{\beta }(-(t/\tau {)}^{\alpha }),$

where $\delta (t)$ is the Dirac delta function and

$E_{\alpha ,\beta }^{\gamma }(z)=\frac{1}{\mathrm{\Gamma }(\gamma )}\sum _{k=0}^{\mathrm{\infty }}\frac{\mathrm{\Gamma }(\gamma +k)z^k}{k!\mathrm{\Gamma }(\alpha k+\beta )}$

is a special instance of the Fox–Wright function and, precisely, it is the three parameters Mittag-Leffler function also known as the Prabhakar function. The function $E_{\alpha ,\beta }^{\gamma }(z)$ can be numerically evaluated, for instance, by means of a Matlab code .

See also

• Debye relaxation
• Cole–Cole equation
• Cole–Davidson equation
• Curie–von Schweidler law
• Dielectric spectroscopy
• Dipole