Sie sind hier: ICP » R. Hilfer » Publikationen

11 Experimental Evidence

[218.2.1] If fractional time evolutions from eq. (27) with \alpha<1 must be expected on general grounds, then they should be observable in experiment. [218.2.2] Numerous experimental examples of anomalous dynamics or strange kinetics have been identified (see [17, 18, 19, 20] and the present volume for reviews). [218.2.3] Here the example of dielectric \alpha-relaxation in glasses is briefly discussed[53, 54], because it provides experimental data over up to 19 decades in time [55], and because the explanation of its excess wing has been a matter of debate.

Figure 1: Five different fits to the real part \varepsilon^{\prime}(f) of the complex dielectric function of propylene carbonate at T=193K as a function of frequency f. Experimental data represented by crosses are from Ref. [56]. The fitting functions corresponds to an exponential (Debye), stretched exponential (KWW), Cole-Davidson [57, 58] Havriliak-Negami [59, 60] and the fractional dynamics (FD) relaxation from (40) The range over which the data were fitted is indicated by dashed vertical lines in the figure. For clarity the data were displaced vertically by half a decade each. The original location of the data corresponds to the curve labelled FD.
Figure 2: Five different fits as in Figure 1 for the imaginary part \varepsilon^{{\prime\prime}}(f) of the complex dielectric function of propylene carbonate at T=193K as a function of frequency f.

[218.3.1] For every induced time evolution on {G} with time scale \tau>0 and fractional order \beta(\tau)

{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau)}}^{{\tau t}}={\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau)}}^{{\tau _{1}t}}{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau)}}^{{(\tau-\tau _{1})t}}={\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau)}}^{{\tau _{1}t}}{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau _{1}+\tau _{2})}}^{{\tau _{2}t}} (36)

holds generally with \tau=\tau _{1}+\tau _{2}. [218.3.2] A physical system typically shows different physical phenomena on different time scales \tau _{i}. [218.3.3] In [53, 54] it was assumed that the second factor in eq. (36) becomes approximately fractional in the sense that

{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau _{1}+\tau _{2})}}^{{\tau _{2}t}}\approx{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\alpha(\tau _{2})}}^{{\tau _{2}t}} (37)

holds in the weak* or strong topology with

\lim _{{\tau _{2}\to 0}}\alpha(\tau _{2})=\beta(\tau). (38)

[218.3.4] The resulting composite time evolution {\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\beta(\tau)}}^{{\tau _{1}t}}{\rule[-2.0pt]{0.0pt}{10.0pt}_{{{G}}}\! T}_{{\alpha(\tau _{2})}}^{{\tau _{2}t}} was studied in [53, 54] for the case \beta(\tau)=1. [218.3.5] Rescaling this composite operator as in the case of Debye relaxation and computing the infinitesimal generator yields the fractional differential equation [53, 54]

\tau _{1}\frac{\mathrm{d}f}{\mathrm{d}t}+\tau _{2}^{\alpha}\mathrm{A}_{\alpha}f=-f (39)

[page 219, §0]    with \mathrm{A}_{\alpha} from eq. (33) and inital value f(0)=1. [219.0.1] Its solution is

f(t)=\mathrm{E}_{{(1,1-\alpha),1}}\left(-\frac{t}{\tau _{1}},-\frac{\tau _{2}^{\alpha}}{\tau _{1}}t^{{1-\alpha}}\right) (40)

where

\mathrm{E}_{{(a_{1},a_{2}),b}}(z_{1},z_{2})=\sum _{{k=0}}^{\infty}\underset{\ell _{1}+\ell _{2}=k}{\sum _{{\ell _{1}\geq 0}}\sum _{{\ell _{2}\geq 0}}}\frac{k!}{\ell _{1}!\ell _{2}!}\frac{z_{1}^{{\ell _{1}}}z_{2}^{{\ell _{2}}}}{\Gamma(b+a_{1}\ell _{1}+a_{2}\ell _{2})} (41)

with a_{1},a_{2}>0 and b,z_{1},z_{2}\in\mathbb{C} is the binomial Mittag-Leffler function [61].

[page 220, §1]    [220.1.1] The complex frequency dependent susceptibility is obtained from the the normalized relaxation function as \chi(u)=1-u\widehat{f}(u) where \widehat{f}(u) is the Laplace transform of f(t) and u=\text{i}\omega is the imaginary circular frequency [54, p. 402, eq.(18)]. [220.1.2] The real part of the complex dielectric susceptibility for propylene carbonate at temperature T=193K is plotted in Figure 1, its imaginary part in Figure 2. [220.1.3] These figures are taken from [53]. [220.1.4] Crosses represent experimental data. [220.1.5] Different fit functions are shifted by half a decade for better visibility. [220.1.6] The range over which the data were fitted is indicated by two dashed vertical lines. [page 221, §0]    [221.0.1] The curve labelled FD (fractional dynamics) is the susceptibility corresponding to the relaxation function in eq. (40) It reproduces the high frequency wing even outside the range of its fit. [221.0.2] This is not the case for the other four curves curves, labelled Debye, KWW, CD and HN. They correspond to four popular fit functions for dielectric relaxation [55, 62]. [221.0.3] The curve Debye corresponds to a simple exponential function, KWW (Kohlrausch-Williams-Watts) is a stretched exponential relaxation function. [221.0.4] The relaxation functions for the two remaining cases, CD (Cole-Davison) and HN (Havriliak-Negami) were given for the first time in [58, 60].

[221.1.1] Figure 3 from [54] shows the real and imaginary part of the dielectric susceptibility for glycerol as its temperature varies over the glass transition range from T=323K to T=184K. [221.1.2] The fits are based on a trinomial fractional relaxation function as detailed in [54, 61].

Figure 3: Separate fits for real part (upper figure) and imaginary part (lower figure) of the complex dielectric susceptibility \chi(\nu)=\chi^{\prime}(\nu)+\text{i}\chi^{{\prime\prime}}(\nu) of glycerol for temperatures T=323, 303, 295, 289, 273, 263, 253, 243, 234, 223, 213, 204, 195, 184 K (from right to left) as function of frequency \nu (from [54]). The experimental data are taken from Ref. [56], the fit uses a generalized composite fractional relaxation model. For details see [54].