Gaussian-to-non-Gaussian transition in supercooled fluids

Understanding Relaxation in the Kob-Andersen Liquid Based on Entropy, String, Shoving, Localization, and Parabolic Models

We assess the validity of a range of models of glass formation based on molecular dynamics simulation results of the Kob-Andersen (KA) model system under a wide range of constant volume and constant pressure conditions. These models include the Adam-Gibbs model emphasizing configurational entropy, the string model emphasizing collective particle exchange motion, the shoving model emphasizing material elasticity, the localization model emphasizing dynamical free volume, and parabolic models based on the ideas of dynamic facilitation and, alternatively, the hypothesis that glass formation involves an avoided critical point. We demonstrate that these seemingly disparate models all provide a reasonable description of structural relaxation and diffusion data for the KA model system under all simulation conditions considered. Hence, the present study points to some unity in our understanding of the relationship between leading models of glass formation, supporting inferences drawn from previous studies of polymeric glass-forming liquids.

Transport Dynamics of Water Molecules Confined between Lipid Membranes

Water molecules confined between biological membranes exhibit a distinctive non-Gaussian displacement distribution, far different from that of bulk water. Here, we introduce a new transport equation for water molecules in the intermembrane space, quantitatively explaining molecular dynamics simulation results. We find that the unique transport dynamics of water molecules stems from the lateral diffusion coefficient fluctuation caused by their longitudinal motion in the direction perpendicular to the membranes. We also identify an interfacial region where water possesses distinct physical properties, which is unaffected by changes in the intermembrane separation.

Connecting real glasses to mean-field models

We propose a novel model for a glass-forming liquid, which allows us to switch in a continuous manner from a standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducing k additional particle-particle interactions, which thus augments the effective number of neighbors of each particle. Our computer simulations of this system show that the structure of the liquid does not change with the introduction of these pseudo-neighbors and by means of analytical calculations, and we determine the structural properties related to these additional neighbors. We show that the relaxation dynamics of the system slows down very quickly with the increase in k and that the onset and the mode-coupling temperatures increase. The systems with high values of k follow the mode-coupling theory power law behavior for a larger temperature range compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamic heterogeneity decreases with the increase in k, whereas the non-Gaussian parameter is independent of it. Thus, we conclude that with the increase in the number of pseudo-neighbors, the system becomes more mean-field-like. By comparing our results with previous studies on mean-field-like systems, we come to the conclusion that the details of how the mean-field limit is approached are important since they can lead to different dynamical behavior in this limit.

Fragility and cooperative motion in a glass-forming polymer-nanoparticle composite

Polymer-nanoparticle composites play a vital role in ongoing materials development. The behavior of the glass transition of these materials is important for their processing and applications, and also represents a problem of fundamental physical interest. Changes of the polymer glass transition temperature Tg due to nanoparticles have been fairly well catalogued, but the breadth of the transition and how rapidly transport properties vary with temperature T - termed the fragility m of glass-formation - is comparatively poorly understood. In the present work, we calculate both Tg and m of a model polymer nanocomposite by molecular dynamics simulations. We systematically consider how Tg and m vary both for the material as a whole, as well as locally, for a range of nanoparticle (NP) concentrations and two polymer-NP interactions. We find large positive and negative changes in Tg and m that can be interpreted in terms of the Adam-Gibbs model of glass-formation, where the scale of the cooperative motion is identified with the scale of string-like cooperative motion. This provides a molecular perpective of fragility changes due to the addition of NPs and for glass formation more generally. We also contrast the behavior along isobaric and isochoric approaches to Tg , since these differing paths can be important to compare experiments (isobaric) and simulations (very often isochoric). Our findings have practical implications for understanding the properties of nanocomposites and fundamental significance for understanding the properties glass-forming materials more broadly.

Uncovering chromatin's contribution to the mitotic spindle: Applications of computational and polymer models

The mitotic spindle is a structure that forms during mitosis to help ensure that each daughter cell receives a full complement of genetic material. In metaphase, the spindle contains microtubules that nucleate inward from two opposing poles. Chromosomes are attached to plus-ends of these microtubules via protein structures called kinetochores. The centromere is the specific region of kinetochore attachment on the chromosome. Chromatin surrounding the centromere (pericentric chromatin) is subject to microtubule-based forces and is commonly modeled as a linear spring, where the force that it exerts is proportional to the distance that it is stretched. We have incorporated physically based models of chromatin to create more accurate and predictive models of the spindle. In addition, using fluorescence microscopy and motion analysis of fluorescently labeled chromatin spots we discovered that pericentric chromatin is restrained relative to free diffusive motion. The characterization of chromatin is crucial to understand mitotic spindle stability and to understand the cell cycle checkpoint regulating anaphase onset.

On the study of collective dynamics in supercooled liquids through the statistics of the isoconfigurational ensemble

The use of the isoconfigurational ensemble to explore structure-dynamic correlations in supercooled liquids is examined. The statistical error of the dynamic propensity and its spatial distribution are determined. The authors present the spatial distribution of the particle non-Gaussian parameter as a measure of the intermittency with which particles exhibit their propensity for motion. The ensemble average of the direction of particle motion is introduced to establish the anisotropy of the dynamic propensity.

Non-Gaussian nature of glassy dynamics by cage to cage motion

A model based on a single Brownian particle moving in a periodic effective field is used to understand the non-Gaussian dynamics in glassy systems of cage escape and subsequent recaging, often thought to be caused by a heterogeneous glass structure. The results are compared to molecular-dynamics simulations of systems with varying complexity: quasi-two-dimensional colloidlike particles, atactic polystyrene, and a dendritic glass. The model nicely describes generic features of all three topologically different systems, in particular around the maximum of the non-Gaussian parameter. This maximum is a measure for the average distance between cages.

Non-Gaussian effects, space-time decoupling, and mobility bifurcation in glassy hard-sphere fluids and suspensions

Brownian trajectory simulation methods are employed to fully establish the non-Gaussian fluctuation effects predicted by our nonlinear Langevin equation theory of single particle activated dynamics in glassy hard-sphere fluids. The consequences of stochastic mobility fluctuations associated with the space-time complexities of the transient localization and barrier hopping processes have been determined. The incoherent dynamic structure factor was computed for a range of wave vectors and becomes of an increasingly non-Gaussian form for volume fractions beyond the (naive) ideal mode coupling theory (MCT) transition. The non-Gaussian parameter (NGP) amplitude increases markedly with volume fraction and is well described by a power law in the maximum restoring force of the nonequilibrium free energy profile. The time scale associated with the NGP peak becomes much smaller than the alpha relaxation time for systems characterized by significant entropic barriers. An alternate non-Gaussian parameter that probes the long time alpha relaxation process displays a different shape, peak intensity, and time scale of its maximum. However, a strong correspondence between the classic and alternate NGP amplitudes is predicted which suggests a deep connection between the early and final stages of cage escape. Strong space-time decoupling emerges at high volume fractions as indicated by a nondiffusive wave vector dependence of the relaxation time and growth of the translation-relaxation decoupling parameter. Displacement distributions exhibit non-Gaussian behavior at intermediate times, evolving into a strongly bimodal form with slow and fast subpopulations at high volume fractions. Qualitative and semiquantitative comparisons of the theoretical results with colloid experiments, ideal MCT, and multiple simulation studies are presented.

Nature of the breakdown in the Stokes-Einstein relationship in a hard sphere fluid

Molecular dynamics simulations of high density hard sphere fluids clearly show a breakdown of the Stokes-Einstein equation (SE). This result has been conjectured to be due to the presence of mobile particles, i.e., ones which have the propensity to "hop" distances that are integer multiples of the interparticle distance. We conclusively show that the sedentary particles, i.e., ones complementary to the "hoppers," obey the SE relationship to a good approximation, even though the fluid as a whole violates the SE equation at high densities. These results support the notion that the unusual diffusive behavior of supercooled liquids is dominated by the hopping particles.

How reproducible are dynamic heterogeneities in a supercooled liquid?

The particle dynamics in a liquid exhibits a transient spatial distribution of dynamic heterogeneities. The relationship between this kinetic structure and the underlying particle configuration remains an outstanding problem. In this Letter, we present a general simulation technique for identifying the features of the dynamic heterogeneity which arise due to a specific configuration, as distinct from the random spatial variation due to the intermittent particle dynamics.

Molecular dynamics simulation of the fragile glass-former orthoterphenyl: A flexible molecule model

We present a realistic model of the fragile glass-former orthoterphenyl and the results of extensive molecular dynamics simulations in which we investigated its basic static and dynamic properties. In this model the internal molecular interactions between the three rigid phenyl rings are described by a set of force constants, including harmonic and anharmonic terms; the interactions among different molecules are described by Lennard-Jones site-site potentials. Self-diffusion properties are discussed in detail together with the temperature and momentum dependencies of the self-intermediate scattering function. The simulation data are compared with existing experimental results and with the main predictions of the mode-coupling theory.

Orientational relaxation in Brownian rotors with frustrated interactions on a square lattice

We present simulation results on the equilibrium relaxation of Brownian planar rotors based on a uniformly frustrated XY model on a square lattice. The rotational relaxation exhibits typical dynamic features of fragile supercooled liquids including the two-step relaxation. We observe a dynamic crossover from the high-temperature regime with Arrhenius behavior to the low-temperature regime with temperature-dependent activation energy. A consistent picture for the observed slow dynamics can be given in terms of the caging effect and thermal activation across potential barriers in the energy landscapes.