Nonlinear susceptibilities of interacting polar molecules in the self-consistent field approximation

The nonlinear stationary response of an assembly of long range interacting electric dipoles is calculated via Berne's nonlinear rotational diffusion equation [J. Chem. Phys. 62 1154 (1975)]. Analytical formulas are derived, showing that the behavior of ω and 3ω components of the nonlinear external field response spectra in such polar dielectrics strongly deviates from the Coffey-Paranjape formulas [Proc. R. Ir. Acad., Sect. A 78, 17 (1978)] as the long range dipole-dipole interactions increase, while the linear response remains qualitatively unaffected. By qualitatively comparing with recent experimental measurements on glycerol, it is further demonstrated that nonlinear dielectric response experiments provide a powerful tool for the characterization of short range intermolecular interactions in polar liquids.