Anomalous dielectric relaxation in strong ac external fields

Dielectric relaxation of complex polar fluids is considered in the context of the anomalous diffusion characterized by a fractional parameter alpha < or = 1 (subdiffusion). An infinite hierarchy of three-term differential-recurrence equations governing the time evolution of the electric polarization is established by following a purely phenomenological procedure. The matrix-continued fraction method is used to derive the exact numerical solution of the stationary regime for an assembly of nonelectrically interacting, polar symmetric-top molecules in presence of a strong ac electric field. The results so obtained are valid to any order in the field strength parameter gamma1, thus extending previous theories applicable to fields of very small amplitudes only. This is illustrated by Cole-Cole diagrams and three-dimensional relaxation spectra for the first- and third-harmonic components of the electric polarization as a function of alpha, gamma1, and the angular frequency.