Probes of heterogeneity in rotational dynamics: application to supercooled liquid CS2

Connection of translational and rotational dynamical heterogeneities with the breakdown of the Stokes-Einstein and Stokes-Einstein-Debye relations in water

We study the Stokes-Einstein (SE) and the Stokes-Einstein-Debye (SED) relations, Dt=kBT/6pietaR and Dr=kBT/8pietaR3, where Dt and Dr are the translational and rotational diffusivity, respectively, T is the temperature, eta the viscosity, kB the Boltzmann constant, and R the "molecular" radius. Our results are based on molecular dynamics simulations of the extended simple point charge model of water. We find that both the SE and SED relations break down at low temperature. To explore the relationship between these breakdowns and dynamical heterogeneities (DHs), we also calculate the SE and SED relations for subsets of the 7% "fastest" and 7% "slowest" molecules. We find that the SE and SED relations break down in both subsets, and that the breakdowns occur on all scales of mobility. Thus these breakdowns appear to be generalized phenomena, in contrast with a view where only the most mobile molecules are the origin of the breakdown of the SE and SED relations, embedded in an inactive background where these relations hold. At low temperature, the SE and SED relations in both subsets of molecules are replaced with "fractional" SE and SED relations, Dt approximately (tau/T)-xit and Dr approximately (tau/T)-xir, where xit approximately 0.84(<1) and xir approximately 0.75(<1). We also find that there is a decoupling between rotational and translational motion, and that this decoupling occurs in both the fastest and slowest subsets of molecules. Further, we find that, the decoupling increases upon cooling, but that the probability of a molecule being classified as both translationally and rotationally fastest also increases. To study the effect of time scale for SE and SED breakdown and decoupling, we introduce a time-dependent version of the SE and SED relations, and a time-dependent function that measures the extent of decoupling. Our results suggest that both the decoupling and SE and SED breakdowns originate at a time scale corresponding to the end of the cage regime, when diffusion starts. This is also the time scale when the DHs are more relevant. Our work also demonstrates that selecting DHs on the basis of translational or rotational motion more strongly biases the calculation of diffusion constants than other dynamical properties such as relaxation times.

Computational probes of molecular motion in the Lewis-Wahnström model for ortho-terphenyl

We use molecular dynamics simulations to investigate translational and rotational diffusion in a rigid three-site model of the fragile glass former ortho-terphenyl, at 260 K< or =T< or =346 K and ambient pressure. An Einstein formulation of rotational motion is presented, which supplements the commonly used Debye model. The latter is shown to break down at supercooled temperatures as the mechanism of molecular reorientation changes from small random steps to large infrequent orientational jumps. We find that the model system exhibits non-Gaussian behavior in translational and rotational motion, which strengthens upon supercooling. Examination of particle mobility reveals spatially heterogeneous dynamics in translation and rotation, with a strong spatial correlation between translationally and rotationally mobile particles. Application of the Einstein formalism to the analysis of translation-rotation decoupling results in a trend opposite to that seen in conventional approaches based on the Debye formalism, namely, an enhancement in the effective rate of rotational motion relative to translation upon supercooling.

On the mechanism of reorientational and structural relaxation in supercooled liquids: the role of border dynamics and cooperativity

Molecular dynamics simulation and analysis based upon the many-body potential energy landscape (PEL) are employed to characterize single molecule reorientation and structural relaxation, and their interrelation, in deeply supercooled liquid CS(2). The rotational mechanism changes from small-step Debye diffusion to sudden large angle reorientation (SLAR) as the temperature falls below the mode-coupling temperature T(c). The onset of SLAR is explained in terms of the PEL; it is an essential feature of low-T rotational dynamics, along with the related phenomena of dynamic heterogeneity and the bifurcation of slow and fast relaxation processes. A long trajectory in which the system is initially trapped in a low energy local minimum, and eventually escapes, is followed in detail, both on the PEL and in real space. During the trapped period, "return" dynamics occurs, always leading back to the trap. Structural relaxation is identified with irreversible escape to a new trap. These processes lead to weak and strong SLAR, respectively; strong SLAR is a clear signal of structural relaxation. Return dynamics involves small groups of two to four molecules, while a string-like structure composed of all the active groups participates in the escape. It is proposed that, rather than simple, nearly instantaneous, one-dimensional barrier crossings, relaxation involves activation of the system to the complex, multidimensional region on the borders of the basins of attraction of the minima for an extended period.

On the Breakdown of the Stokes−Einstein Law in Supercooled Liquids

The wavevector-dependent shear viscosity, eta(k), is evaluated for a range of temperatures in a supercooled binary Lennard-Jones liquid. The mode coupling theory of Keyes and Oppenheim (Phys. Rev. A 1973, 8, 937) expresses the self-diffusion constant, D, in terms of eta(k). Replacing eta(k) with the usual viscosity, eta identical with eta(k = 0), yields the Stokes-Einstein law. It is found that the breakdown of the SE law in this system is well described by keeping the simulated k-dependence. Simply put, bath processes on all length scales (wavevectors) contribute to D, the system is much less viscous at finite k, and thus D exceeds the SE estimate based upon eta. The functional form of eta(k) allows for the estimation of a correlation length that grows with decreasing T.

Onset of slow dynamics in difluorotetrachloroethane glassy crystal

Complementary neutron spin-echo and x-ray experiments and molecular-dynamics simulations have been performed on difluorotetrachloroethane (CFCl2-CFCl2) glassy crystal. Static, single-molecule reorientational dynamics and collective dynamics properties are investigated. Our results confirm the strong analogy between molecular liquids and plastic crystals. The orientational disorder is characterized at different temperatures and a change in the nature of rotational dynamics is observed. A careful check of the rotational diffusion model is performed using self-angular correlation functions Cl with high l values and compared to results obtained on molecular liquids composed of A-B dumbbells. Below the crossover temperature at which slow dynamics emerge, we show that some scaling predictions of the mode coupling theory hold and that alpha-relaxation times and nonergodicity parameters are controlled by the nontrivial static correlations.

Superdipole liquid scenario for the dielectric primary relaxation in supercooled polar liquids

We propose a dynamic structure of coupled dynamic molecular strings for supercooled small polar molecule liquids and accordingly we obtain the Hamiltonian of the rotational degrees of freedom of the system. From the Hamiltonian, the strongly correlated supercooled polar liquid state is renormalized to a normal superdipole liquid state. This scenario describes the following main features of the primary or alpha-relaxation dynamics in supercooled polar liquids: (1) the average relaxation time evolves from a high temperature Arrhenius to a low temperature non-Arrhenius or super-Arrhenius behavior; (2) the relaxation function crosses over from the high temperature exponential to low temperature nonexponential form; and (3) the temperature dependence of the relaxation strength shows non-Curie features. According to the present model, the crossover phenomena of the first two characteristics arise from the transition between the superdipole gas and the superdipole liquid. The model predictions are quantitatively compared with the experimental results of glycerol, a typical glass former.

Mode coupling behavior in glass-forming liquid crystalline isopentylcyanobiphenyl

Linear and nonlinear dielectric measurements of liquid crystalline chiral isopentylcyanobiphenyl (5*CB) and n -pentylcyanobiphenyl (5CB), combined with viscosity eta (T) data, are presented. The 5*CB compound glassifies on cooling in the cholesteric phase whereas 5CB crystallizes in the nematic phase. In both compounds the temperature evolution of dielectric relaxation times, the dc conductivity, and the viscosity are well described by the "critical-like" description from mode coupling theory (MCT). However, for 5*CB a unique coincidence of the MCT "critical" temperature and extrapolated temperature of the hypothetical continuous isotropic-cholesteric (T*) phase transition was found. The temperature dependence of the strong electric-field-induced changes of the dielectric permittivity exhibits a strong anomaly in the direction of negative values on approaching T* , not observed up to now. The anomaly is described by the susceptibility-related critical exponent gamma=1 . The divergence of the "nonlinear" dielectric relaxation follows a power dependence described by the exponent y=1 . This paper recalls the recent discussions on the glassy dynamics of a "hard-ellipsoid" liquid and the possible relationship between the glass transition, critical phenomena, and isotropic-nematic transition.