Fokker-Planck equation dielectric relaxation

Chapter 8 by W. T. Coffey, Y. P. Kalmykov, and S. V. Titov, entitled Fractional Rotational Diffusion and Anomalous Dielectric Relaxation in Dipole Systems, provides an introduction to the theory of fractional rotational Brownian motion and microscopic models for dielectric relaxation in disordered systems. The authors indicate how anomalous relaxation has its origins in anomalous diffusion and that a physical explanation of anomalous diffusion may be given via the continuous time random walk model. It is demonstrated how this model may be used to justify the fractional diffusion equation. In particular, the Debye theory of dielectric relaxation of an assembly of polar molecules is reformulated using a fractional noninertial Fokker-Planck equation for the purpose of extending that theory to explain anomalous dielectric relaxation. Thus, the authors show how the Debye rotational diffusion model of dielectric relaxation of polar molecules (which may be described in microscopic fashion as the diffusion limit of a discrete time random walk on the surface of the unit sphere) may be extended via the continuous-time random walk to yield the empirical Cole-Cole, Cole-Davidson, and Havriliak-Negami equations of anomalous dielectric relaxation from a microscopic model based on a... 

We shall now almost exclusively concentrate on the fractal time random walk excluding inertial effects and the discrete orientation model of dielectric relaxation. We shall demonstrate how in the diffusion limit this walk will yield a fractional generalization of the Debye-Frohlich model. Just as in the conventional Debye relaxation, a fractional generalization of the Debye-Frohlich model may be derived from a number of very different models of the relaxation process (compare the approach of Refs. 22, 23, 28 and 34—36). The advantage of using an approach based on a kinetic equation such as the fractional Fokker-Planck equation (FFPE) however is that such a method may easily be extended to include the effects of the inertia of the dipoles, external potentials, and so on. Moreover, the FFPE (by use of a theorem of operational calculus generalized to fractional exponents and continued fraction methods) clearly indicates how many existing results of the classical theory of the Brownian motion may be extended to include fractional dynamics. 

Thus we have demonstrated how the empirical Havriliak-Negami equation [Eq. (11)] can be obtained from a microscopic model, namely, the fractional Fokker-Planck equation [Eq. (101)] applied to noninteracting rotators. This model can explain the anomalous relaxation of complex dipolar systems, where the anomalous exponents ct and v differ from unity (corresponding to the classical Debye theory of dielectric relaxation) that is, the relaxation process is characterized by a broad distribution of relaxation times. Hence, the empirical Havriliak-Negami equation of anomalous dielectric relaxation which has been... 

In the present section, it is demonstrated how the linear response of an assembly of noninteracting polar Brownian particles to a small external field F applied parallel and perpendicular to the bias field Fo may be calculated in the context of the fractional noninertial rotational diffusion in the same manner as normal rotational diffusion [8]. In order to carry out the calculation, it is assumed that the rotational Brownian motion of a particle may be described by a fractional noninertial Fokker-Planck (Smoluchowski) equation, in which the inertial effects are neglected. Both exact and approximate solutions of this equation are presented. We shall demonstrate that the characteristic times of the normal diffusion process, namely, the integral and effective relaxation times obtained in Refs. 8, 65, and 67, allow one to evaluate the dielectric response for anomalous diffusion. Moreover, these characteristic times yield a simple analytical equation for the complex dielectric susceptibility tensor describing the anomalous relaxation of the system. The exact solution of the problem reduces to the solution of the infinite hierarchies of differential-recurrence equations for the corresponding relaxation functions. The longitudinal and transverse components of the susceptibility tensor may be calculated exactly from the Laplace transform of these relaxation functions using linear response theory [72]. 

---