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

Boltzmann's H-function for molecules with orientational degrees of freedom: Emergence of unique features

Boltzmann's H-function H(t), often regarded as an analog of time-dependent entropy, holds a venerable place in the history of science. However, accurate numerical evaluation of H(t) for particles other than atoms is rare. To remove this lacuna, we generalize Boltzmann's H-function to a gas of molecules with orientational degrees of freedom and evaluate H(t) from the time-dependent joint probability distribution function f(p, L, t) for linear (p) and angular (L) momenta, evolving from an initial nonequilibrium state, by molecular dynamics simulations. We consider both prolate- and oblate-shaped particles, interacting via the well-known Gay-Berne potential and obtain the relaxation of the generalized molecular H(t) from initial (t = 0) nonequilibrium states. In the long-time limit, the H function saturates to its exact equilibrium value, which is the sum of translational and rotational contributions to the respective entropies. Both the translational and rotational components of H(t) decay nearly exponentially with time; the rotational component is more sensitive to the molecular shape that enters through the aspect ratio. A remarkable rapid decrease in the rotational relaxation time is observed as the spherical limit is approached, in a way tantalizingly reminiscent of Hu-Zwanzig hydrodynamic prediction with the slip boundary condition. In addition, we obtain H(t) analytically by solving the appropriate translational and rotational Fokker-Planck equation and obtain a modest agreement with simulations. We observe a remarkable signature of translation-rotation coupling as a function of molecular shape, captured through a physically meaningful differential term that quantifies the magnitude of translation-rotation coupling.

What is special about how roaming chemical reactions traverse their potential surfaces? Differences in geodesic paths between roaming and non-roaming events

With the notable exception of some illustrative two-degree-of-freedom models whose surprising classical dynamics has been worked out in detail, theories of roaming have largely bypassed the issue of when and why the counterintuitive phenomenon of roaming occurs. We propose that a useful way to begin to address these issues is to look for the geodesic (most efficient) pathways through the potential surfaces of candidate systems. Although roaming manifests itself in an unusual behavior at asymptotic geometries, we found in the case of formaldehyde dissociation that it was the pathways traversing the parts of the potential surface corresponding to highly vibrationally excited reactants that were the most revealing. An examination of the geodesics for roaming pathways in this region finds that they are much less tightly defined than the geodesics in that same region that lead directly to dissociation (whether into closed-shell products or into radical products). Thus, the broader set of options available to the roaming channel gives it an entropic advantage over more conventional reaction channels. These observations suggest that what leads to roaming in other systems may be less the presence of a localized "roaming transition state," than the existence of an entire region of the potential surface conducive to multiple equivalent pathways.

The inherent dynamics of isotropic- and nematic-phase liquid crystals

The geodesic (shortest) pathways through the potential energy landscape of a liquid can be thought of as defining what its dynamics would be if thermal noise were removed, revealing what we have called the "inherent dynamics" of the liquid. We show how these inherent paths can be located for a model liquid crystal former, showing, in the process, how the molecular mechanisms of translation and reorientation compare in the isotropic and nematic phases of these systems. These mechanisms turn out to favor the preservation of local orientational order even under macroscopically isotropic conditions (a finding consistent with the experimental observation of pseudonematic domains in these cases), but disfavor the maintenance of macroscopic orientational order, even in the nematic phase. While the most efficient nematic pathways that maintain nematic order are indeed shorter than those that do not, it is apparently difficult for the system to locate these paths, suggesting that molecular motion in liquid-crystal formers is dynamically frustrated, and reinforcing the sense that there are strong analogies between liquid crystals and supercooled liquids.

The inherent dynamics of a molecular liquid: geodesic pathways through the potential energy landscape of a liquid of linear molecules

Because the geodesic pathways that a liquid follows through its potential energy landscape govern its slow, diffusive motion, we suggest that these pathways are logical candidates for the title of a liquid's "inherent dynamics." Like their namesake "inherent structures," these objects are simply features of the system's potential energy surface and thus provide views of the system's structural evolution unobstructed by thermal kinetic energy. This paper shows how these geodesic pathways can be computed for a liquid of linear molecules, allowing us to see precisely how such molecular liquids mix rotational and translational degrees of freedom into their dynamics. The ratio of translational to rotational components of the geodesic path lengths, for example, is significantly larger than would be expected on equipartition grounds, with a value that scales with the molecular aspect ratio. These and other features of the geodesics are consistent with a picture in which molecular reorientation adiabatically follows translation-molecules largely thread their way through narrow channels available in the potential energy landscape.

Potential energy landscape and inherent dynamics of a hard-sphere fluid

Hard-sphere models exhibit many of the same kinds of supercooled-liquid behavior as more realistic models of liquids, but the highly nonanalytic character of their potentials makes it a challenge to think of that behavior in potential energy landscape terms. We show here that it is possible to calculate an important topological property of hard-sphere landscapes, the geodesic pathways through those landscapes, and to do so without artificially coarse-graining or softening the potential. We show, moreover, that the rapid growth of the lengths of those pathways with increasing packing fraction quantitatively predicts the precipitous decline in diffusion constants in a glass-forming hard-sphere mixture model. The geodesic paths themselves can be considered as defining the intrinsic dynamics of hard spheres, so it is also revealing to find that they (and therefore the features of the underlying potential energy landscape) correctly predict the occurrence of dynamic heterogeneity and nonzero values of the non-Gaussian parameter. The success of these landscape predictions for the dynamics of such a singular model emphasizes that there is more to potential energy landscapes than is revealed by looking at the minima and saddle points.

Single molecule probing of dynamics in supercooled polymers

Fluorescence experiments with single BODIPY molecules embedded in a poly(methyl acrylate) matrix have been performed at various temperatures in the supercooled regime. By using pulsed excitation, fluorescence lifetime and linear dichroism time trajectories were accessible at the same time. Both observables have been analyzed without data binning. While the linear dichroism solely reflects single particle dynamics, the fluorescence lifetime observable depends on the molecular environment, so that the dynamics from the polymer host surrounding a chromophore contributes to this quantity. We observe that the lifetime correlation decays slightly faster than polarization correlation, indicating the occurrence of large angular reorientations. Additionally, dichroism time trajectories have been adducted to reveal directly the geometry of rotational dynamics. We identify small but also significantly larger rotational jumps being responsible for the overall molecular reorientation.

Dynamics in the metabasin space of a Lennard-Jones glass former: connectivity and transition rates

Using simulations, we construct the effective dynamics in metabasin space for a Lennard-Jones glass former. Metabasins are identified via a scheme that measures transition rates between inherent structures and generates clusters of inherent structures by drawing in branches that have the largest transition rates. This construction is fundamentally different from the stochastic approach based on molecular-dynamics trajectories. The effective dynamics in this metabasin space is shown to be Markovian but to differ significantly from the simplest trap models. We specifically show that retaining information about the connectivity in this metabasin space is crucial for reproducing the slow dynamics observed in this system.

Exploring the potential energy landscape of glass-forming systems: from inherent structures via metabasins to macroscopic transport

In this review a systematic analysis of the potential energy landscape (PEL) of glass-forming systems is presented. Starting from the thermodynamics, the route towards the dynamics is elucidated. A key step in this endeavor is the concept of metabasins. The relevant energy scales of the PEL can be characterized. Based on the simulation results for some glass-forming systems one can formulate a relevant model system (ideal Gaussian glass-former) which can be treated analytically. The macroscopic transport can be related to the microscopic hopping processes, using either the strong relation between energy (thermodynamics) and waiting times (dynamics) or, alternatively, the concepts of the continuous-time random walk. The relation to the geometric properties of the PEL is stressed. The emergence of length scales within the PEL approach as well as the nature of finite-size effects is discussed. Furthermore, the PEL view is compared to other approaches describing the glass transition.

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.

Thermodynamic and dynamic anomalies for dumbbell molecules interacting with a repulsive ramplike potential

Using collision driven discrete molecular dynamics, we investigate the thermodynamics and dynamics of systems of 500 dumbbell molecules interacting by a purely repulsive ramplike discretized potential, consisting of n steps of equal size. We compare the behavior of the two systems, with n=18 and n=144 steps. Each system exhibits both thermodynamic and dynamic anomalies, a density maximum and the translational and rotational mobilities show anomalous behavior. Starting with very dense systems and decreasing the density, both mobilities first increase, reach a maximum, then decrease, reach a minimum, and finally increase; this behavior is similar to the behavior of SPC/E (Simple Point Charge-Extended) water. The regions in the pressure-temperature plane of translational and rotational mobility anomalies depend strongly on n. The product of the translational diffusion coefficient and the orientational correlation time increases with temperature, in contrast with the behavior of most liquids.

Decoupling phenomena in supercooled liquids: signatures in the energy landscape

A significant deviation from the Debye model of rotational diffusion in the dynamics of orientational degrees of freedom in an equimolar mixture of ellipsoids of revolution and spheres is found to begin at a temperature at which the average inherent structure energy of the system starts falling with drop in temperature. We argue that this onset temperature corresponds to the emergence of the process as a distinct mode of orientational relaxation. Further, we find that the coupling between rotational and translational diffusion breaks down at a still lower temperature where a change occurs in the temperature dependence of the average inherent structure energy.

Relation between rotational and translational dynamic heterogeneities in water

We use molecular dynamics simulations to probe the rotational dynamics of the extended simple point charge model of water for a range of temperatures down to 200 K, 6 K above the mode coupling temperature. We find that rotational dynamics is spatially heterogeneous; i.e., there are clusters of molecules that rotate significantly more than the average for a given time interval, and we study the size and the temporal behavior of these clusters. We find that the position of a rotational heterogeneity is strongly correlated with the position of a translational heterogeneity, and that the fraction of molecules belonging to both kinds of heterogeneities increases with decreasing temperature. We further find that although the two types of heterogeneities are not identical, they are related to the same physical picture.

Signatures of the fast dynamics in glassy polystyrene: First evidence by high-field Electron Paramagnetic Resonance of molecular guests

The reorientation of one small paramagnetic molecule (spin probe) in glassy polystyrene (PS) is studied by high-field electron paramagnetic resonance spectroscopy at two different Larmor frequencies (190 and 285 GHz). Two different regimes separated by a crossover region are evidenced. Below 180 K the rotational times are nearly temperature independent with no apparent distribution. In the temperature range of 180-220 K a large increase of the rotational mobility is observed with the widening of the distribution of correlation times which exhibits two components: (i) a deltalike, temperature-independent component representing the fraction of spin probes w which persist in the low-temperature dynamics; (ii) a strongly temperature-dependent component, to be described by a power distribution, representing the fraction of spin probes 1-w undergoing activated motion over an exponential distribution of barrier heights g(E). Above 180 K a steep decrease of w is evidenced. The shape and the width of g(E) do not differ from the reported ones for PS within the errors. For the first time the large increase of the rotational mobility of the spin probe at 180 K is ascribed to the onset of the fast dynamics detected by neutron scattering at T(f)=175+/-25 K.

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.