Integrating factor

Havriliak–Negami relaxation is an empirical modification of the Debye relaxation model, accounting 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.

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.

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

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

and

${\epsilon }^{″}(\omega )={(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)$

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}$

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

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

with the 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$

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.

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. Unfortunately, there are no numerical algorithms available for the computation of such functions with the required parameters.

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.

See also

• Cole–Cole equation
• Dielectric spectroscopy
• Dipole