Viscous flow and jump dynamics in molecular supercooled liquids. II. Rotations

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Random walk model of mobility in glass formers

A mechanism responsible for the generic features of the mean squared displacement and the decay of the orientational autocorrelator of a molecule in a glass forming liquid is poorly understood, where such a mechanism would be critical for creating the theory of glass transition. A discrete random walk model is proposed where, instead of a straight line, a walk is along a tortuous path consisting of blocks of switchback ramps. Subdiffusive regime, short-term dynamic heterogeneity, and existence of the α- and β-relaxation processes emerge naturally from the model. The model suggests that slowing of the rate of relaxation may be due to an increase in the number of the switchback ramps per block rather than due to growth of an energy barrier as typically assumed.

Rotational and translational dynamics in dense fluids of patchy particles

We explore the effect of directionality on rotational and translational relaxation in glassy systems of patchy particles. Using molecular dynamics simulations, we analyze the impact of two distinct patch geometries, one that enhances the local icosahedral structure and the other one that does not strongly affect the local order. We find that in nearly all investigated cases, rotational relaxation takes place on a much faster time scale than translational relaxation. By comparing to a simplified dynamical Monte Carlo model, we illustrate that rotational diffusion can be qualitatively explained as purely local motion within a fixed environment, which is not coupled strongly to the cage-breaking dynamics required for translational relaxation. Nonetheless, icosahedral patch placement has a profound effect on the local structure of the system, resulting in a dramatic slowdown at low temperatures, which is strongest at an intermediate "optimal" patch size.

Fast Vibrational Modes and Slow Heterogeneous Dynamics in Polymers and Viscous Liquids

Many systems, including polymers and molecular liquids, when adequately cooled and/or compressed, solidify into a disordered solid, i.e., a glass. The transition is not abrupt, featuring progressive decrease of the microscopic mobility and huge slowing down of the relaxation. A distinctive aspect of glass-forming materials is the microscopic dynamical heterogeneity (DH), i.e., the presence of regions with almost immobile particles coexisting with others where highly mobile ones are located. Following the first compelling evidence of a strong correlation between vibrational dynamics and ultraslow relaxation, we posed the question if the vibrational dynamics encodes predictive information on DH. Here, we review our results, drawn from molecular-dynamics numerical simulation of polymeric and molecular glass-formers, with a special focus on both the breakdown of the Stokes-Einstein relation between diffusion and viscosity, and the size of the regions with correlated displacements.

Coupling and decoupling between translational and rotational dynamics in supercooled monodisperse soft Janus particles

We perform dynamics simulations to investigate the translational and rotational glassy dynamics in a glass-forming liquid of monodisperse soft Janus particles. We find that, with decreasing temperature, the mean-square angular displacement shows no clear plateau in the caging region, in contrast with the apparent caging behavior of translational motion. By defining a reorientational mean-square angular displacement, the caging behavior of rotational motion can be recognized. On approaching the glass transition (decreasing temperature), the coupling between translational and rotational relaxation increases, while the coupling between translational and rotational diffusion decreases, whereas the coupling between translational and reorientational diffusion increases. The strong decoupling between translational and rotational diffusion is due to the suppressed translational mobility but promoted rotational mobility of soft Janus particles. We think that the low-T SE and SED decoupling is mainly attributed to hopping motion of soft Janus particles, whereas the high-T SE and SED decoupling is mainly attributed to collective cage motion of soft Janus particles. Our results demonstrate that interaction anisotropy has a critical effect on the translational and rotational dynamics of soft Janus particles.

The effect of polymorphism on the structural, dynamic and dielectric properties of plastic crystal water: A molecular dynamics simulation perspective

We have employed molecular dynamics simulations based on the TIP4P/2005 water model to investigate the local structural, dynamical, and dielectric properties of the two recently reported body-centered-cubic and face-centered-cubic plastic crystal phases of water. Our results reveal significant differences in the local orientational structure and rotational dynamics of water molecules for the two polymorphs. The probability distributions of trigonal and tetrahedral order parameters exhibit a multi-modal structure, implying the existence of significant local orientational heterogeneities, particularly in the face-centered-cubic phase. The calculated hydrogen bond statistics and dynamics provide further indications of the existence of a strongly heterogeneous and rapidly interconverting local orientational structural network in both polymorphs. We have observed a hindered molecular rotation, much more pronounced in the body-centered-cubic phase, which is reflected by the decay of the fourth-order Legendre reorientational correlation functions and angular Van Hove functions. Molecular rotation, however, is additionally hindered in the high-pressure liquid compared to the plastic crystal phase. The results obtained also reveal significant differences in the dielectric properties of the polymorphs due to the different dipolar orientational correlation characterizing each phase.

Peculiar Molecular Shape and Size Dependence of the Dynamics of Fluids Confined in a Small-Pore Metal-Organic Framework

Force-field-based molecular dynamics simulations were deployed to systematically explore the dynamics of confined molecules of different shapes and sizes, that is, linear (CO₂ and N₂) and spherical (CH₄) fluids, in a model small pore system, that is, the metal-organic framework SIFSIX-2-Cu-i. These computations unveil an unprecedented molecular symmetry dependence of the translational and rotational dynamics of fluids confined in channel-like nanoporous materials. In particular, this peculiar behavior is reflected by the extremely slow decay of the Legendre reorientational correlation functions of even-parity order for the linear fluids, which is associated with jump-like orientation flips, while the spherical fluid shows a very fast decay taking place on a subpicosecond time scale. Such a fundamental understanding is relevant to diverse disciplines such as in chemistry, physics, biology, and materials science, where diatomic or polyatomic molecules of different shapes/sizes diffuse through nanopores.

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.

Rotational dynamics of simple asymmetric molecules

Molecular dynamic simulations were carried out on rigid diatomic molecules, which exhibit both α (structural) and β (secondary) dynamics. The relaxation scenarios range from onset behavior, in which a distinct α process emerges on cooling, to merging behavior, associated with two relaxation peaks that converge at higher temperature. These properties, as well as the manifestation of the β peak as an excess wing, depend not only on thermodynamic conditions, but also on both the symmetry of the molecule and the correlation function (odd or even) used to analyze its dynamics. These observations help to reconcile divergent results obtained from different experiments. For example, the β process is more intense and the α-relaxation peak is narrower in dielectric relaxation spectra than in dynamic light scattering or NMR measurements. In the simulations herein, this follows from the weaker contribution of the secondary relaxation to even-order correlation functions, related to the magnitude of the relevant angular jumps.

Orientational glass in mixtures of elliptic and circular particles: structural heterogeneities, rotational dynamics, and rheology

Using molecular dynamics simulation with an angle-dependent Lennard-Jones potential, we study orientational glass with quadrupolar symmetry in mixtures of elliptic particles and circular impurities in two dimensions. With a mild aspect ratio (= 1.23) and a mild size ratio (= 1.2), we realize a plastic crystal at relatively high temperature T. With further lowering T, we find a structural phase transition for very small impurity concentration c and pinned disordered orientations for not small c. The ellipses are anchored by the impurities in the planar alignment. With increasing c, the orientation domains composed of isosceles triangles gradually become smaller, resulting in orientational glass with crystal order. In our simulation, the impurity distribution becomes heterogeneous during quenching from liquid, which then produces rotational dynamic heterogeneities. We also examine rheology in orientational glass to predict a shape memory effect and a superelasticity effect, where a large fraction of the strain is due to collective orientation changes.

Characteristics of the Johari-Goldstein process in rigid asymmetric molecules

Molecular dynamics simulations were carried out on a Lennard-Jones binary mixture of rigid (fixed bond length) diatomic molecules. The translational and rotational correlation functions, and the corresponding susceptibilities, exhibit two relaxation processes: the slow structural relaxation (α dynamics) and a higher frequency secondary relaxation. The latter is a Johari-Goldstein (JG) process, by its definition of involving all parts of the molecule. It shows several properties characteristic of the JG relaxation: (1) merging with the α relaxation at high temperature; (2) a change in temperature dependence of its relaxation strength on vitrification; (3) a separation in frequency from the α peak that correlates with the breadth of the α dispersion; and (4) sensitivity to volume, pressure, and physical aging. These properties can be used to determine whether a secondary relaxation in a real material is an authentic JG process, rather than more trivial motion involving intramolecular degrees of freedom. The latter has no connection to the glass transition, whereas the JG relaxation is closely related to structural relaxation, and thus can provide new insights into the phenomenon.

Formation of double glass in binary mixtures of anisotropic particles: dynamic heterogeneities in rotations and displacements

We study glass behavior in a mixture of elliptic and circular particles in two dimensions at low temperatures using an orientation-dependent Lennard-Jones potential. The ellipses have a mild aspect ratio (∼1.2) and tend to align at low temperatures, while the circular particles play the role of impurities disturbing the ellipse orientations at a concentration of 20%. These impurities have a size smaller than that of the ellipses and attract them in the homeotropic alignment. As a result, the coordination number around each impurity is mostly 5 or 4 in glassy states. We realize double glass, where both the orientations and the positions are disordered but still hold mesoscopic order. We find a strong heterogeneity in the flip motions of the ellipses, which sensitively depends on the impurity clustering. In our model, a small fraction of the ellipses still undergo flip motions relatively rapidly even at low temperatures. In contrast, the nonflip rotations (with angle changes not close to ±π) are mainly caused by the cooperative configuration changes involving many particles. Then, there arises a long-time heterogeneity in the nonflip rotations closely correlated with the dynamic heterogeneity in displacements.

Heterogeneity in Single-Molecule Observables in the Study of Supercooled Liquids

Bulk approaches to studying heterogeneous systems obscure important details, as they report average behavior rather than the distribution of behaviors in such environments. Small-molecule and polymeric supercooled liquids, which display heterogeneity in their dynamics without an underlying structural heterogeneity that sets those dynamics, are important constituents of this category of condensed matter systems. A variety of approaches have been devised to unravel ensemble averaging in supercooled liquids. This review focuses on the ultimate subensemble approach, single-molecule measurements, as they have been applied to the study of supercooled liquids. We detail how three key experimental observables (single-molecule probe rotation, translation, and fluorescence lifetime) have been employed to provide detail on dynamic heterogeneity in supercooled liquids. Special attention is given to the potential for, but also the challenges in, discriminating spatial and temporal heterogeneity and detailing the length scales and timescales of heterogeneity in these systems.

Dynamical Line-Shifts in High-Field Electron Spin Resonance: Applications to Polymer Physics

High-field high-frequency Electron Paramagnetic Resonance spectroscopy (HF 2 -EPR) is a powerful tool to investigate, with ultra-high angular resolution, the rotational dynamics of complex systems like polymers, viscous fluids and glasses. Usually, information is drawn by detailed numerical analysis of the overall lineshape. Here, we present a simplified analytical model of the line shifts due to the rotational dynamics of the paramagnetic centre. The model captures the basic features of the reorientation process (time scale and size of the angular jump). It is compared with experimental results concerning the reorientation of a paramagnetic guest molecule dissolved in polystyrene. We find that, if the rotational model to describe the reorientation of the radical is consistent, the best-fit parameters yield equally acceptable best-fits of the overall spectrum by numerical simulations and dynamical line shifts by independent analytic expressions.

Molecular probe dynamics reveals suppression of ice-like regions in strongly confined supercooled water

The structure of the hydrogen bond network is a key element for understanding water's thermodynamic and kinetic anomalies. While ambient water is strongly believed to be a uniform, continuous hydrogen-bonded liquid, there is growing consensus that supercooled water is better described in terms of distinct domains with either a low-density ice-like structure or a high-density disordered one. We evidenced two distinct rotational mobilities of probe molecules in interstitial supercooled water of polycrystalline ice [Banerjee D, et al. (2009) ESR evidence for 2 coexisting liquid phases in deeply supercooled bulk water. Proc Natl Acad Sci USA 106: 11448–11453]. Here we show that, by increasing the confinement of interstitial water, the mobility of probe molecules, surprisingly, increases. We argue that loose confinement allows the presence of ice-like regions in supercooled water, whereas a tighter confinement yields the suppression of this ordered fraction and leads to higher fluidity. Compelling evidence of the presence of ice-like regions is provided by the probe orientational entropy barrier which is set, through hydrogen bonding, by the configuration of the surrounding water molecules and yields a direct measure of the configurational entropy of the same. We find that, under loose confinement of supercooled water, the entropy barrier surmounted by the slower probe fraction exceeds that of equilibrium water by the melting entropy of ice, whereas no increase of the barrier is observed under stronger confinement. The lower limit of metastability of supercooled water is discussed.

Rotational dynamics of a dipolar supercooled liquid

We study the rotational dynamics of a supercooled molecular liquid by means of molecular dynamics simulations. The system under investigation is composed of rigid diatomic molecules with an associate dipole moment. At room temperature, orientational correlations decrease rapidly with increasing distances. Upon cooling, angles between dipole moments of molecules within the first coordination shell decrease. As for the dynamical properties, rotational diffusion coefficients decrease with temperature at a smaller rate than translational diffusion coefficients do, and the critical temperature associated with the former is lower than the one corresponding to their translational counterparts. Translation and rotation about an inertial axis are uncorrelated, whereas some coupling between translation and dipole reorientation is obtained.

Molecular dynamics simulation of the Johari-Goldstein relaxation in a molecular liquid

Molecular dynamics simulations were carried out to investigate the reorientational motion of a rigid (fixed bond length), asymmetric diatomic molecule in the liquid and glassy states. In the latter the molecule reorients via large-angle jumps, which we identify with the Johari-Goldstein (JG) dynamics. This relaxation process has a broad distribution of relaxation times, and at least deeply in the glassy state, the mobility of a given molecule remains fixed over time; that is, there is no dynamic exchange among molecules. Interestingly, the JG relaxation time for a molecule does not depend on the local density, although the nonergodicity factor is weakly correlated with the packing efficiency of neighboring molecules. In the liquid state the intensity of the JG process increases significantly, eventually subsuming the slower α relaxation. This evolution of the JG motion into structural relaxation underlies the correlation of many properties of the JG and α dynamics.

Manifestations of probe presence on probe dynamics in supercooled liquids

Experimental studies that follow behavior of single probes embedded in heterogeneous systems are increasingly common. The presence of probes may perturb the system, and such perturbations may or may not affect interpretation of host behavior from the probe observables typically measured. In this study, the manifestations of potential probe-induced changes to host dynamics in supercooled liquids are investigated via molecular dynamics simulations. It is found that probe dynamics do not necessarily mirror host dynamics as they exist either in the probe-free or probe-bearing systems. In particular, for a binary supercooled liquid, we find that smooth probes larger than the host particles induce increased translational diffusion in the host system; however, the diffusion is anisotropic and enhances caging of the probe, suppressing probe translational diffusion. This in turn may lead experiments that follow probe diffusion to suggest Stokes-Einstein behavior of the system even while both the probe-free and probe-bearing systems exhibit deviations from that behavior.

ESR evidence for 2 coexisting liquid phases in deeply supercooled bulk water

Using electron spin resonance spectroscopy (ESR), we measure the rotational mobility of probe molecules highly diluted in deeply supercooled bulk water and negligibly constrained by the possible ice fraction. The mobility increases above the putative glass transition temperature of water, T(g) = 136 K, and smoothly connects to the thermodynamically stable region by traversing the so called "no man's land" (the range 150-235 K), where it is believed that the homogeneous nucleation of ice suppresses the liquid water. Two coexisting fractions of the probe molecules are evidenced. The 2 fractions exhibit different mobility and fragility; the slower one is thermally activated (low fragility) and is larger at low temperatures below a fragile-to-strong dynamic cross-over at approximately 225 K. The reorientation of the probe molecules decouples from the viscosity below approximately 225 K. The translational diffusion of water exhibits a corresponding decoupling at the same temperature [Chen S-H, et al. (2006) The violation of the Stokes-Einstein relation in supercooled water. Proc Natl Acad Sci USA 103:12974-12978]. The present findings are consistent with key issues concerning both the statics and the dynamics of supercooled water, namely the large structural fluctuations [Poole PH, Sciortino F, Essmann U, Stanley HE (1992) Phase behavior of metastable water. Nature 360:324-328] and the fragile-to-strong dynamic cross-over at approximately 228 K [Ito K, Moynihan CT, Angell CA (1999) Thermodynamic determination of fragility in liquids and a fragile-to-strong liquid transition in water. Nature 398:492-494].

Orientational relaxation in a dispersive dynamic medium: generalization of the Kubo-Ivanov-Anderson jump-diffusion model to include fractional environmental dynamics

The Ivanov-Anderson model (and an earlier treatment by Kubo) envisages a decay of the orientational correlation by random but large amplitude molecular jumps, as opposed to infinitesimal small jumps assumed in Brownian diffusion. Recent computer simulation studies on water and viscous liquids have shown that large amplitude motions may indeed be more of a rule than exception. Existing theoretical studies on jump diffusion mostly assume an exponential (Poissonian) waiting time distribution for jumps, thereby again leading to an exponential decay. Here we extend the existing formalism of Ivanov and Anderson to include an algebraic waiting time distribution between two jumps. As a result, the first (l=1) and second (l=2) rank orientational time correlation functions show the same long time power law, but their short time decay behavior is quite different. The predicted Cole-Cole plot of dielectric relaxation reproduces various features of non-Debye behavior observed experimentally. We also developed a theory where both unrestricted small jumps and large angular jumps coexist simultaneously. The small jumps are shown to have a large effect on the long time decay, particularly in mitigating the effects of algebraic waiting time distribution, and in giving rise to an exponential-like decay, with a time constant, surprisingly, less than the time constant that arises from small amplitude decay alone.

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.

Glassy orientational dynamics of rodlike molecules near the isotropic-nematic transition

We investigate the single-particle orientational dynamics of rodlike molecules across the isotropic-nematic transition in computer simulations of a family of model systems of thermotropic liquid crystals. Several remarkable features of glassy dynamics are on display including nonexponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. In order to obtain a quantitative measure of glassy dynamics in line with the established methods in supercooled liquids, we construct a relaxation time versus scaled inverse temperature plot and demonstrate that one can indeed define a "fragility index" for thermotropic liquid crystals that depends on density and aspect ratio. The values of the fragility parameter are surprisingly in the range observed for glass-forming liquids. A plausible correlation between the energy landscape features and the observed fragility is discussed.

Orientational relaxation in a discotic liquid crystal

We investigate orientational relaxation of a model discotic liquid crystal, consisting of disclike molecules, by molecular dynamics simulations along two isobars starting from the high temperature isotropic phase. The two isobars have been so chosen that (a) the phase sequence isotropic- (I-) nematic- (N-) columnar (C) appears upon cooling along one of them and (b) the sequence isotropic- (I-) columnar- (C) along the other. While the orientational relaxation in the isotropic phase near the I-N phase transition in system (a) shows a power law decay at short to intermediate times, such power law relaxation is not observed in the isotropic phase near the I-C phase boundary in system (b). In order to understand this difference (the existence or the absence of the power law decay), we calculated the growth of the orientational pair distribution functions (OPDFs) near the I-N phase boundary and also near the I-C phase boundary. We find that the OPDF shows a marked growth in long range correlation as the I-N phase boundary is approached in the I-N-C system (a), but such a growth is absent in the I-C system, which appears to be consistent with the result that I-N phase transition in the former is weakly first order while the I-C phase transition in the latter is not weak. As the system settles into the nematic phase, the decay of the single-particle second-rank orientational time correlation function follows a pattern that is similar to what is observed with calamitic liquid crystals and supercooled molecular liquids.

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.

Fractional Stokes-Einstein and Debye-Stokes-Einstein relations in a network-forming liquid

We study the breakdown of the Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) relations for translational and rotational motion in a prototypical model of a network-forming liquid, the ST2 model of water. We find that the emergence of fractional SE and DSE relations at low temperature is ubiquitous in this system, with exponents that vary little over a range of distinct physical regimes. We also show that the same fractional SE relation is obeyed by both mobile and immobile dynamical heterogeneities of the liquid.

Power law relaxation and glassy dynamics in Lebwohl-Lasher model near the isotropic-nematic phase transition

Orientational dynamics in a liquid crystalline system near the isotropic-nematic (IN) phase transition is studied using molecular dynamics simulations of the well-known Lebwohl-Lasher model. As the IN transition temperature is approached from the isotropic side, we find that the decay of the orientational time correlation functions (OTCF) slows down noticeably, giving rise to a power law decay at intermediate time scales. The angular velocity time correlation function also exhibits a rather pronounced power law decay near the IN boundary. In the mean squared angular displacement at comparable time scales, we observe the emergence of a subdiffusive regime which is followed by a superdiffusive regime before the onset of the long-time diffusive behavior. We observe signature of dynamical heterogeneity through pronounced non-Gaussian behavior in orientational motion particularly at lower temperatures. This behavior closely resembles what is usually observed in supercooled liquids. We obtain the free energy as a function of orientational order parameter by the use of the transition matrix Monte Carlo method. The free energy surface is flat for the system considered here and the barrier between isotropic and nematic phases is vanishingly small for this weakly first-order phase transition, hence allowing large scale, collective, and correlated orientational density fluctuations. This might be responsible for the observed power law decay of the OTCFs.

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.

Complete breakdown of the Debye model of rotational relaxation near the isotropic-nematic phase boundary: effects of intermolecular correlations in orientational dynamics

The Debye-Stokes-Einstein (DSE) model of rotational diffusion predicts that the orientational correlation times tau l vary as [l(l+1)]-1, where l is the rank of the orientational time correlation function (given in terms of the Legendre polynomial of rank l). One often finds significant deviation from this prediction, in either direction. In supercooled molecular liquids where the ratio tau 1/tau 2 falls considerably below 3 (the Debye limit), one usually invokes a jump diffusion model to explain the approach of the ratio tau 1/tau 2 to unity. Here we show in a computer simulation study of a standard model system for thermotropic liquid crystals that this ratio becomes much less than unity as the isotropic-nematic phase boundary is approached from the isotropic side. Simultaneously, the ratio tau 2/eta, eta, being the shear viscosity of the liquid, becomes much larger than the hydrodynamic value near the I-N transition. We also analyze the breakdown of the Debye model of rotational diffusion in ratios of higher order orientational correlation times. We show that the breakdown of the DSE model is due to the growth of orientational pair correlation and provide a mode coupling theory analysis to explain the results.

Dynamic arrest in a liquid of symmetric dumbbells: Reorientational hopping for small molecular elongations

We present extensive equilibrium and out-of-equilibrium molecular-dynamics simulations of a liquid of symmetric dumbbell molecules, for constant packing fraction, as a function of temperature and molecular elongation. We compute diffusion constants as well as odd and even orientational correlators. The notations odd and even refer to the parity of the order l of the corresponding Legendre l polynomial, evaluated for the orientation of the molecular axis relative to its initial position. Rotational degrees of freedom of order l are arrested if, in the long-time limit, the corresponding orientational l correlator does not decay to zero. It is found that for large elongations translational and rotational degrees of freedom freeze at the same temperature. For small elongations only the even rotational degrees of freedom remain coupled to translational motions and arrest at a finite common temperature. On the contrary, the odd rotational degrees of freedom remain ergodic at all investigated temperatures. Hence, in the translationally arrested state, each molecule remains trapped in the cage formed by its neighboring molecules, but is able to perform 180 degrees rotations, which lead to relaxation only for the odd orientational correlators. The temperature dependence of the characteristic time of these residual rotations is well described by an Arrhenius law. Finally, we discuss the evidence in favor of the presence of the type-A transition for the odd rotational degrees of freedom, as predicted by mode-coupling theory for small molecular elongations. This transition is distinct from the type-B transition, associated with the arrest of the translational and even rotational degrees of freedom for small elongations, and with all degrees of freedom for large elongations. Odd orientational correlators are computed for small elongations at very low temperatures in the translationally arrested state. The obtained results suggest that hopping events restore the ergodicity of the odd rotational degrees of freedom at temperatures far below the A transition.

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.

Interconversion of mechanical and dielectrical relaxation measurements for dicyclohexylmethyl-2-methyl succinate

A comparison between results of dielectrical relaxation and dynamic mechanical spectroscopies is carried out for the alpha-relaxation of the ester dicyclohexyl methyl-2-methyl succinate (DCMMS). The results for the dielectric permittivity and the shear modulus measurements are presented according to the empirical Havriliak-Negami (HN) equation. By using the time-temperature principle a master curve in each case was obtained for several temperatures. The comparative analysis presented here is based on the assumption of a relationship between rotational and shear viscosities. The former one is associated to the dielectrical relaxation, whereas the latter is associated to mechanical relaxation. Both viscosities are not necessarily equal in general, and we assume that the difference between them is an important factor to appropriately compare the dielectrical and mechanical results.

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.

Molecular dynamics study of the thermal and the density effects on the local and the large-scale motion of polymer melts: scaling properties and dielectric relaxation

Results from a molecular dynamics simulation of a melt of unentangled polymers are presented. The translational motion, the large-scale and the local reorientation processes of the chains, as well as their relations with the so-called "normal" and "segmental" dielectric relaxation modes are thoroughly investigated in wide temperature and pressure ranges. The thermodynamic states are well fitted by the phenomenological Tait equation of state. A global time-temperature-pressure superposition principle of both the translational and the rotational dynamics is evidenced. The scaling is more robust than the usual Rouse model. The latter provides insight but accurate comparison with the simulation calls for modifications to account for both the local chain stiffness and the nonexponential relaxation. The study addresses the issue whether the temperature or the density is a dominant control parameter of the dynamics or the two quantities give rise to comparable effects. By examining the ratio /alphatau//alphaP between the isochronic and isobaric expansivities, one finds that the temperature is dominant when the dynamics is fast. If the relaxation slows down, the fluctuations of the free volume increase their role and become comparable to those of the thermal energy. Detectable cross-correlation between the "normal-mode" and the "segmental" dielectric relaxations is found and contrasted with the usual assumption of independent modes.

Anomalous glassy relaxation near the isotropic-nematic phase transition

Dynamical heterogeneity in a system of Gay-Berne ellipsoids near its isotropic-nematic (I-N) transition, and also in an equimolar mixture of Lennard-Jones spheres and Gay-Berne ellipsoids in deeply supercooled regime, is probed by the time evolution of non-Gaussian parameters (NGP). The appearance of a dominant second peak in the rotational NGP near the I-N transition signals the growth of pseudonematic domains. Surprisingly, such a second peak is instead observed in the translational NGP for the glassy binary mixture. Localization of orientational motion near the I-N transition is found to be responsible for the observed anomalous orientational relaxation.