Decoupling of exchange and persistence times in atomistic models of glass formers

The role of functionalization in the translocation of peptides through multilayer graphene nanopores

The rapid translocation speed of peptides through graphene nanopores poses a challenge, hindering the accurate sensing of the biomarkers. Employing the functionalized graphene nanopores is at the forefront of reducing the translocation speed. The current work details the translocation of a negatively charged peptide endothelin-1 through a bare multilayer graphene nanopore, a hydrogen-functionalized graphene nanopore, and a hydroxyl-functionalized graphene nanopore by applying electric fields. The hydroxyl-functionalized graphene nanopore significantly reduces the peptide's translocation speed. The time required for the peptide to translocate through the hydroxyl-functionalized graphene nanopore is 2.25 times longer than in the non-functionalized graphene nanopore and 1.25 times longer than in the hydrogen-functionalized graphene nanopore. We critically analyze the factors influencing the reduced translocation speed, including the interactions between the pore and the peptide, the conformational changes of the peptide within the pore, the solvent velocity inside the pore, and the solvent's viscosity near the peptide. The altered solvent velocities within functionalized pores have a minimal role in the speed reduction of peptides. When a constant force is applied to the peptide without any electric field, the hydroxyl-functionalized graphene nanopore delivers the lowest diffusion rate. The persistence time, which serves as a measure of the solvent viscosity near the peptide, is the highest within the hydroxyl-functionalized pore. Finally, we conclude that the Coulombic interactions between the peptide and the pore play a major role in its speed reduction inside the hydroxyl-functionalized graphene nanopore.

Glassy phases of the Gaussian core model

We present results from molecular dynamics simulations exploring the supercooled dynamics of the Gaussian Core Model in the low- and intermediate-density regimes. In particular, we analyse the transition from the low-density hard-sphere-like glassy dynamics to the high-density one. The dynamics at low densities is well described by the caging mechanism, giving rise to intermittent dynamics. At high densities, the particles undergo a more continuous motion in which the concept of cage loses its meaning. We elaborate on the idea that these different supercooled dynamics are in fact the precursors of two different glass states.

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts - often in collaboration with experimental groups - to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

Excitations follow (or lead?) density scaling in propylene carbonate

Structural excitations that enable interbasin (IB) barrier crossings on a potential energy landscape are thought to play a facilitating role in the relaxation of liquids. Here, we show that the population of these excitations exhibits the same density scaling observed for α relaxation in propylene carbonate, even though they are heavily influenced by intramolecular modes. We also find that IB crossing modes exhibit a Gr[Formula: see text]neisen parameter ( γ G) that is approximately equivalent to the density scaling parameter γ TS. These observations suggest that the well-documented relationship between γ G and γ TS may be a direct result of the pressure dependence of the frequency of unstable (relaxation) modes associated with IB motion.

Kinetic frustration by limited bond availability controls the LAT protein condensation phase transition on membranes

LAT is a membrane-linked scaffold protein that undergoes a phase transition to form a two-dimensional protein condensate on the membrane during T cell activation. Governed by tyrosine phosphorylation, LAT recruits various proteins that ultimately enable condensation through a percolation network of discrete and selective protein-protein interactions. Here, we describe detailed kinetic measurements of the phase transition, along with coarse-grained model simulations, that reveal that LAT condensation is kinetically frustrated by the availability of bonds to form the network. Unlike typical miscibility transitions in which compact domains may coexist at equilibrium, the LAT condensates are dynamically arrested in extended states, kinetically trapped out of equilibrium. Modeling identifies the structural basis for this kinetic arrest as the formation of spindle arrangements, favored by limited multivalent binding interactions along the flexible, intrinsically disordered LAT protein. These results reveal how local factors controlling the kinetics of LAT condensation enable formation of different, stable condensates, which may ultimately coexist within the cell.

Single molecule demonstration of Debye-Stokes-Einstein breakdown in polystyrene near the glass transition temperature

Rotational-translational decoupling, in which translational motion is apparently enhanced over rotational motion in violation of Stokes-Einstein (SE) and Debye-Stokes-Einstein (DSE) predictions, has been observed in materials near their glass transition temperatures (T_g). This has been posited to result from ensemble averaging in the context of dynamic heterogeneity. In this work, ensemble and single molecule experiments are performed in parallel on a fluorescent probe in high molecular weight polystyrene near its T_g. Ensemble results show decoupling onset at approximately 1.15T_g, increasing to over three orders of magnitude at T_g. Single molecule measurements also show a high degree of decoupling, with typical molecules at T_g showing translational diffusion coefficients nearly 400 times higher than expected from SE/DSE predictions. At the single molecule level, higher degree of breakdown is associated with particularly mobile molecules and anisotropic trajectories, providing support for anomalous diffusion as a critical driver of rotational-translational decoupling and SE/DSE breakdown.

Experiments with high-molecular-weight polystyrene provide insights into the mechanisms behind rotational-translational decoupling in glassy systems. Specifically, particularly mobile molecules exhibiting anisotropic trajectories are found to play a key role in Debye-Stokes-Einstein breakdown.

Comparing Microscopic and Macroscopic Dynamics in a Paradigmatic Model of Glass-Forming Molecular Liquid

Glass transition is a most intriguing and long-standing open issue in the field of molecular liquids. From a macroscopic perspective, glass-forming systems display a dramatic slowing-down of the dynamics, with the inverse diffusion coefficient and the structural relaxation times increasing by orders of magnitude upon even modest supercooling. At the microscopic level, single-molecule motion becomes strongly intermittent, and can be conveniently described in terms of “cage-jump” events. In this work, we investigate a paradigmatic glass-forming liquid, the Kob–Andersen Lennard–Jones model, by means of Molecular Dynamics simulations, and compare the macroscopic and microscopic descriptions of its dynamics on approaching the glass-transition. We find that clear changes in the relations between macroscopic timescales and cage-jump quantities occur at the crossover temperature where Mode Coupling-like description starts failing. In fact, Continuous Time Random Walk and lattice model predictions based on cage-jump statistics are also violated below the crossover temperature, suggesting the onset of a qualitative change in cage-jump motion. Interestingly, we show that a fully microscopic relation linking cage-jump time- and length-scales instead holds throughout the investigated temperature range.

Breakdown of the Stokes-Einstein relation in supercooled liquids: A cage-jump perspective

The breakdown of the Stokes-Einstein relation in supercooled liquids, which is the increase in the ratio τ_ατ_D between the two macroscopic times for structural relaxation and diffusion on decreasing the temperature, is commonly ascribed to dynamic heterogeneities, but a clear-cut microscopic interpretation is still lacking. Here, we tackle this issue exploiting the single-particle cage-jump framework to analyze molecular dynamics simulations of soft disk assemblies and supercooled water. We find that τ_ατ_D∝⟨t_p⟩⟨t_c⟩, where ⟨t_p⟩ and ⟨t_c⟩ are the cage-jump times characterizing slow and fast particles, respectively. We further clarify that this scaling does not arise from a simple term-by-term proportionality; rather, the relations τ_α∝⟨t_p⟩⟨Δr_J ²⟩ and τ_D∝⟨t_c⟩⟨Δr_J ²⟩ effectively connect the macroscopic and microscopic timescales, with the mean square jump length ⟨Δr_J ²⟩ shrinking on cooling. Our work provides a microscopic perspective on the Stokes-Einstein breakdown and generalizes previous results on lattice models to the case of more realistic glass-formers.

Glass-like characteristics of intracellular motion in human cells

The motion in the cytosol of microorganisms such as bacteria and yeast has been observed to undergo a dramatic slowing down upon cell energy depletion. These observations have been interpreted as the motion being "glassy," but whether this notion is useful also for active, motor-protein-driven transport in eukaryotic cells is less clear. Here, we use fluorescence microscopy of beads in human (HeLa) cells to probe the motion of membrane-surrounded structures that are carried along the cytoskeleton by motor proteins. Evaluating several hallmarks of glassy dynamics, we show that at short length scales, the motion is heterogeneous, is nonergodic, is well described by a model for the displacement distribution in glassy systems, and exhibits a decoupling of the exchange and persistence times. Overall, these results suggest that the short length scale behavior of objects that can be transported actively by motor proteins in human cells shares features with the motion in glassy systems.

Rotational dynamics of proteins in nanochannels: role of solvent's local viscosity

Viscosity variation of solvent in local regions near a solid surface, be it a biological surface of a protein or an engineered surface of a nanoconfinement, is a direct consequence of intermolecular interactions between the solid body and the solvent. The current coarse-grained molecular dynamics study takes advantage of this phenomenon to investigate the anomaly in a solvated protein's rotational dynamics confined using a representative solid matrix. The concept of persistence time, the characteristic time of structural reordering in liquids, is used to compute the solvent's local viscosity. With an increase in the degree of confinement, the confining matrix significantly influences the solvent molecule's local viscosity present in the protein hydration layer through intermolecular interactions. This effect contributes to the enhanced drag force on protein motion, causing a reduction in the rotational diffusion coefficient. Simulation results suggest that the direct matrix-protein non-bonded interaction is responsible for the occasional jump and discontinuity in orientational motion when the protein is in very tight confinement.

A comparative study on the dynamic heterogeneity of supercooled polymers under nanoconfinement

Dynamic heterogeneity (DH) is a universal property of glass transition phenomena. In this work, we perform a comparative analysis of DH for pure polymer and polymer/nanoparticle composite systems in both film and bulk states via molecular dynamics simulations. We find that the dynamic gradient and the faster average dynamics due to the presence of a free surface are two leading factors, resulting from a nanoconfinement effect, which influence different parts of DH in a film system. The dynamic gradient results from differences in dynamics at different distances from the mobile surface, which induces a large deviation from the Gaussian distribution for the displacement distribution in the film. At the same time, the maximum string size which describes the region size for cooperative motion (dynamic correlation) can also be influenced by the dynamic gradient, although this influence is much weaker than that on the displacement distribution. On the other hand, reflecting temporal fluctuations of dynamics or temporal parts of DH, characteristic peak times of the non-Gaussian parameter and string size, and the ratio between persistent times and exchange times which describe the dynamic exchange properties, are mainly influenced by the faster dynamics on average. Our results demonstrate that measuring different properties (dynamic distribution, dynamic correlation or dynamic exchange) place an emphasis on distinct temporal and spatial parts of DH. It is necessary to use combinational measurements of these properties to give a complete picture of DH in nanoconfinement environments.

Continuous-time random-walk approach to supercooled liquids: Self-part of the van Hove function and related quantities

From equilibrium molecular dynamics (MD) simulations of a bead-spring model for short-chain glass-forming polymer melts we calculate several quantities characterizing the single-monomer dynamics near the (extrapolated) critical temperature [Formula: see text] of mode-coupling theory: the mean-square displacement g₀(t), the non-Gaussian parameter [Formula: see text] and the self-part of the van Hove function [Formula: see text] which measures the distribution of monomer displacements r in time t. We also determine these quantities from a continuous-time random walk (CTRW) approach. The CTRW is defined in terms of various probability distributions which we know from previous analysis. Utilizing these distributions the CTRW can be solved numerically and compared to the MD data with no adjustable parameter. The MD results reveal the heterogeneous and non-Gaussian single-particle dynamics of the supercooled melt near [Formula: see text]. In the time window of the early [Formula: see text] relaxation [Formula: see text] is large and [Formula: see text] is broad, reflecting the coexistence of monomer displacements that are much smaller ("slow particles") and much larger ("fast particles") than the average at time t, i.e. than [Formula: see text]. For large r the tail of [Formula: see text] is compatible with an exponential decay, as found for many glassy systems. The CTRW can reproduce the spatiotemporal dependence of [Formula: see text] at a qualitative to semiquantitative level. However, it is not quantitatively accurate in the studied temperature regime, although the agreement with the MD data improves upon cooling. In the early [Formula: see text] regime we also analyze the MD results for [Formula: see text] via the space-time factorization theorem predicted by ideal mode-coupling theory. While we find the factorization to be well satisfied for small r, both above and below [Formula: see text] , deviations occur for larger r comprising the tail of [Formula: see text]. The CTRW analysis suggests that single-particle "hops" are a contributing factor for these deviations.

Interfacial Properties and Hopping Diffusion of Small Nanoparticle in Polymer/Nanoparticle Composite with Attractive Interaction on Side Group

The diffusion dynamics of fullerene (C 60 ) in unentangled linear atactic polystyrene (PS) and polypropylene (PP) melts and the structure and dynamic properties of polymers in interface area are investigated by performing all-atom molecular dynamics simulations. The comparison of the results in two systems emphasises the influence of local interactions exerted by polymer side group on the diffusion dynamics of the nanoparticle. In the normal diffusive regime at long time scales, the displacement distribution function (DDF) follows a Gaussian distribution in PP system, indicating a normal diffusion of C 60 . However, we observe multiple peaks in the DDF curve for C 60 diffusing in PS melt, which indicates a diffusion mechanism of hopping of C 60 . The attractive interaction between C 60 and phenyl ring side groups are found to be responsible for the observed hopping diffusion. In addition, we find that the C 60 is dynamically coupled with a subsection of a tetramer on PS chain, which has a similar size with C 60 . The phenyl ring on PS chain backbone tends to have a parallel configuration in the vicinity of C 60 surface, therefore neighbouring phenyl rings can form chelation effect on the C 60 surface. Consequently, the rotational dynamics of phenyl ring and the translational diffusion of styrene monomers are found to be slowed down in this interface area. We hope our results can be helpful for understanding of the influence of the local interactions on the nanoparticle diffusion dynamics and interfacial properties in polymer/nanoparticle composites.

Static and dynamic properties of two-dimensional Coulomb clusters

We study the temperature dependence of static and dynamic responses of Coulomb interacting particles in two-dimensional confinements across the crossover from solid- to liquid-like behaviors. While static correlations that investigate the translational and bond orientational order in the confinements show the footprints of hexatic-like phase at low temperatures, dynamics of the particles slow down considerably in this phase, reminiscent of a supercooled liquid. Using density correlations, we probe long-lived heterogeneities arising from the interplay of the irregularity in the confinement and long-range Coulomb interactions. The relaxation at multiple time scales show stretched-exponential decay of spatial correlations in irregular traps. Temperature dependence of characteristic time scales, depicting the structural relaxation of the system, show striking similarities with those observed for the glassy systems, indicating that some of the key signatures of supercooled liquids emerge in confinements with lower spatial symmetries.

Dynamic heterogeneity in aqueous ionic solutions

It is well known that supercooled liquids have heterogeneous dynamics, but it is still unclear whether dynamic heterogeneity also exists in aqueous ionic solutions at room or even higher temperatures.

Quantifying Allosteric Communication via Both Concerted Structural Changes and Conformational Disorder with CARDS

Allosteric (i.e., long-range) communication within proteins is crucial for many biological processes, such as the activation of signaling cascades in response to specific stimuli. However, the physical basis for this communication remains unclear. Existing computational methods for identifying allostery focus on the role of concerted structural changes, but recent experimental work demonstrates that disorder is also an important factor. Here, we introduce the Correlation of All Rotameric and Dynamical States (CARDS) framework for quantifying correlations between both the structure and disorder of different regions of a protein. To quantify disorder, we draw inspiration from methods for quantifying "dynamic heterogeneity" from chemical physics to classify segments of a dihedral's time evolution as being in either ordered or disordered regimes. To demonstrate the utility of this approach, we apply CARDS to the Catabolite Activator Protein (CAP), a transcriptional activator that is regulated by Cyclic Adenosine MonoPhosphate (cAMP) binding. We find that CARDS captures allosteric communication between the two cAMP-Binding Domains (CBDs). Importantly, CARDS reveals that this coupling is dominated by disorder-mediated correlations, consistent with NMR experiments that establish allosteric coupling between the CBDs occurs without a concerted structural change. CARDS also recapitulates an enhanced role for disorder in the communication between the DNA-Binding Domains (DBDs) and CBDs in the S62F variant of CAP. Finally, we demonstrate that using CARDS to find communication hotspots identifies regions of CAP that are in allosteric communication without foreknowledge of their identities. Therefore, we expect CARDS to be of great utility for both understanding and predicting allostery.

Active cage model of glassy dynamics

We build up a phenomenological picture in terms of the effective dynamics of a tracer confined in a cage experiencing random hops to capture some characteristics of glassy systems. This minimal description exhibits scale invariance properties for the small-displacement distribution that echo experimental observations. We predict the existence of exponential tails as a crossover between two Gaussian regimes. Moreover, we demonstrate that the onset of glassy behavior is controlled only by two dimensionless numbers: the number of hops occurring during the relaxation of the particle within a local cage and the ratio of the hopping length to the cage size.

Timing Correlations in Proteins Predict Functional Modules and Dynamic Allostery

How protein structure encodes functionality is not fully understood. For example, long-range intraprotein communication can occur without measurable conformational change and is often not captured by existing structural correlation functions. It is shown here that important functional information is encoded in the timing of protein motions, rather than motion itself. I introduce the conditional activity function to quantify such timing correlations among the degrees of freedom within proteins. For three proteins, the conditional activities between side-chain dihedral angles were computed using the output of microseconds-long atomistic simulations. The new approach demonstrates that a sparse fraction of side-chain pairs are dynamically correlated over long distances (spanning protein lengths up to 7 nm), in sharp contrast to structural correlations, which are short-ranged (<1 nm). Regions of high self- and inter-side-chain dynamical correlations are found, corresponding to experimentally determined functional modules and allosteric connections, respectively.

Limiting Energy Dissipation Induces Glassy Kinetics in Single-Cell High-Precision Responses

Single cells often generate precise responses by involving dissipative out-of-thermodynamic-equilibrium processes in signaling networks. The available free energy to fuel these processes could become limited depending on the metabolic state of an individual cell. How does limiting dissipation affect the kinetics of high-precision responses in single cells? I address this question in the context of a kinetic proofreading scheme used in a simple model of early-time T cell signaling. Using exact analytical calculations and numerical simulations, I show that limiting dissipation qualitatively changes the kinetics in single cells marked by emergence of slow kinetics, large cell-to-cell variations of copy numbers, temporally correlated stochastic events (dynamic facilitation), and ergodicity breaking. Thus, constraints in energy dissipation, in addition to negatively affecting ligand discrimination in T cells, can create a fundamental difficulty in determining single-cell kinetics from cell-population results.

Particle jumps in structural glasses

Particles in structural glasses rattle around temporary equilibrium positions, that seldom change through a process which is much faster than the relaxation time, known as particle jump. Since the relaxation of the system is due to the accumulation of many such jumps, it could be possible to connect the single particle short time motion to the macroscopic relaxation by understanding the features of the jump dynamics. Here we review recent results in this research direction, clarifying the features of particle jumps that have been understood and those that are still under investigation, and examining the role of particle jumps in different theories of the glass transition.

Heterogeneous dynamics and its length scale in simple ionic liquid models: a computational study

We numerically investigate the dynamic heterogeneity and its length scale found in coarse-grained ionic liquid model systems.

Dynamic heterogeneity controls diffusion and viscosity near biological interfaces

At a nanometer scale, the behavior of biological fluids is largely governed by interfacial physical chemistry. This may manifest as slowed or anomalous diffusion. Here we describe how measures developed for studying glassy systems allow quantitative measurement of interfacial effects on water dynamics, showing that correlated motions of particles near a surface result in a viscosity greater than anticipated from individual particle motions. This effect arises as a fundamental consequence of spatial heterogeneity on nanometer length scales and applies to any fluid near any surface. Increased interfacial viscosity also causes the classic finding that large solutes such as proteins diffuse much more slowly than predicted in bulk water. This has previously been treated via an empirical correction to the solute size: the hydrodynamic radius. Using measurements of quantities from theories of glass dynamics, we can now calculate diffusion constants from molecular details alone, eliminating the empirical correction factor.

Using the s ensemble to probe glasses formed by cooling and aging

From length scale distributions characterizing frozen amorphous domains, we relate the s ensemble method with standard cooling and aging protocols for forming glass. We show that in a class of models where space-time scaling is in harmony with that of experiment, the spatial distributions of excitations obtained with the s ensemble are identical to those obtained through cooling or aging, but the computational effort for applying the s ensemble is generally many orders of magnitude smaller than that of straightforward numerical simulation of cooling or aging. We find that in contrast to the equilibrium ergodic state, a nonequilibrium length scale characterizes the anticorrelation between excitations and encodes the preparation history of glass states.

Dynamic phase coexistence in glass-forming liquids

One of the most controversial hypotheses for explaining the heterogeneous dynamics of glasses postulates the temporary coexistence of two phases characterized by a high and by a low diffusivity. In this scenario, two phases with different diffusivities coexist for a time of the order of the relaxation time and mix afterwards. Unfortunately, it is difficult to measure the single-particle diffusivities to test this hypothesis. Indeed, although the non-Gaussian shape of the van-Hove distribution suggests the transient existence of a diffusivity distribution, it is not possible to infer from this quantity whether two or more dynamical phases coexist. Here we provide the first direct observation of the dynamical coexistence of two phases with different diffusivities, by showing that in the deeply supercooled regime the distribution of the single-particle diffusivities acquires a transient bimodal shape. We relate this distribution to the heterogeneity of the dynamics and to the breakdown of the Stokes-Einstein relation, and we show that the coexistence of two dynamical phases occurs up to a timescale growing faster than the relaxation time on cooling, for some of the considered models. Our work offers a basis for rationalizing the dynamics of supercooled liquids and for relating their structural and dynamical properties.

Dynamic heterogeneity in crossover spin facilitated model of supercooled liquid and fractional Stokes-Einstein relation

Kinetically constrained models have gained much interest as models that assign the origins of interesting dynamic properties of supercooled liquids to dynamical facilitation mechanisms that have been revealed in many experiments and numerical simulations. In this work, we investigate the dynamic heterogeneity in the fragile-to-strong liquid via Monte Carlo method using the model that linearly interpolates between the strong liquid-like behavior and the fragile liquid-like behavior by an asymmetry parameter b. When the asymmetry parameter is sufficiently small, smooth fragile-to-strong transition is observed both in the relaxation time and the diffusion constant. Using these physical quantities, we investigate fractional Stokes-Einstein relations observed in this model. When b is fixed, the system shows constant power law exponent under the temperature change, and the exponent has the value between that of the Frederickson-Andersen model and the East model. Furthermore, we investigate the dynamic length scale of our systems and also find the crossover relation between the relaxation time. We ascribe the competition between energetically favored symmetric relaxation mechanism and entropically favored asymmetric relaxation mechanism to the fragile-to-strong crossover behavior.

Slowing down of ring polymer diffusion caused by inter-ring threading

Diffusion of long ring polymers in a melt is much slower than the reorganization of their internal structures. While direct evidence for entanglements has not been observed in the long ring polymers unlike linear polymer melts, threading between the rings is suspected to be the main reason for slowing down of ring polymer diffusion. It is, however, difficult to define the threading configuration between two rings because the rings have no chain end. In this work, evidence for threading dynamics of ring polymers is presented by using molecular dynamics simulation and applying a novel analysis method. The simulation results are analyzed in terms of the statistics of persistence and exchange times that have proved useful in studying heterogeneous dynamics of glassy systems. It is found that the threading time of ring polymer melts increases more rapidly with the degree of polymerization than that of linear polymer melts. This indicates that threaded ring polymers cannot diffuse until an unthreading event occurs, which results in the slowing down of ring polymer diffusion.

Cooperativity leads to temporally-correlated fluctuations in the bacteriophage lambda genetic switch

Cooperative interactions are widespread in biochemical networks, providing the nonlinear response that underlies behavior such as ultrasensitivity and robust switching. We introduce a temporal correlation function-the conditional activity-to study the behavior of these phenomena. Applying it to the bistable genetic switch in bacteriophage lambda, we find that cooperative binding between binding sites on the prophage DNA lead to non-Markovian behavior, as quantified by the conditional activity. Previously, the conditional activity has been used to predict allosteric pathways in proteins; here, we show that it identifies the rare unbinding events which underlie induction from lysogeny to lysis.

Coupled diffusion in lipid bilayers upon close approach

Biomembrane interfaces create regions of slowed water dynamics in their vicinity. When two lipid bilayers come together, this effect is further accentuated, and the associated slowdown can affect the dynamics of larger-scale processes such as membrane fusion. We have used molecular dynamics simulations to examine how lipid and water dynamics are affected as two lipid bilayers approach each other. These two interacting fluid systems, lipid and water, both slow and become coupled when the lipid membranes are separated by a thin water layer. We show in particular that the water dynamics become glassy, and diffusion of lipids in the apposed leaflets becomes coupled across the water layer, while the "outer" leaflets remain unaffected. This dynamic coupling between bilayers appears mediated by lipid-water-lipid hydrogen bonding, as it occurs at bilayer separations where water-lipid hydrogen bonds become more common than water-water hydrogen bonds. We further show that such coupling occurs in simulations of vesicle-vesicle fusion prior to the fusion event itself. Such altered dynamics at membrane-membrane interfaces may both stabilize the interfacial contact and slow fusion stalk formation within the interface region.

Connecting short and long time dynamics in hard-sphere-like colloidal glasses

Glass-forming materials are characterized by an intermittent motion at the microscopic scale. Particles spend most of their time rattling within the cages formed by their neighbors, and seldom jump to a different cage. In molecular glass formers the temperature dependence of the jump features, such as the average caging time and jump length, characterizes the relaxation processes and allows for a short-time prediction of the diffusivity. Here we experimentally investigate the cage-jump motion of a two-dimensional hard-sphere-like colloidal suspension, where the volume fraction is the relevant parameter controlling the slowing down of the dynamics. We characterize the volume fraction dependence of the cage-jump features and show that, as in molecular systems, they allow for a short time prediction of the diffusivity.

Time scale of dynamic heterogeneity in model ionic liquids and its relation to static length scale and charge distribution

We find a general power-law behavior: , where ζ_dh ≈ 1.2 for all the ionic liquid models, regardless of charges and the length scale of structural relaxation.

Hopping and the Stokes-Einstein relation breakdown in simple glass formers

One of the most actively debated issues in the study of the glass transition is whether a mean-field description is a reasonable starting point for understanding experimental glass formers. Although the mean-field theory of the glass transition--like that of other statistical systems--is exact when the spatial dimension d → ∞, the evolution of systems properties with d may not be smooth. Finite-dimensional effects could dramatically change what happens in physical dimensions,d = 2, 3. For standard phase transitions finite-dimensional effects are typically captured by renormalization group methods, but for glasses the corrections are much more subtle and only partially understood. Here, we investigate hopping between localized cages formed by neighboring particles in a model that allows to cleanly isolate that effect. By bringing together results from replica theory, cavity reconstruction, void percolation, and molecular dynamics, we obtain insights into how hopping induces a breakdown of the Stokes-Einstein relation and modifies the mean-field scenario in experimental systems. Although hopping is found to supersede the dynamical glass transition, it nonetheless leaves a sizable part of the critical regime untouched. By providing a constructive framework for identifying and quantifying the role of hopping, we thus take an important step toward describing dynamic facilitation in the framework of the mean-field theory of glasses.

Renewal events in glass-forming liquids

On cooling toward the glass transition temperature, glass-forming liquids display long periods of localized motion interrupted by fast "jumps" in the single-particle trajectories. Several theoretical models based on these single-particle jumps have been proposed, most prominently the continuous-time random walk (CTRW). The central assumption of the CTRW is that jumps are renewal events, i.e. that the internal clock of a particle can be reset upon a jump. In this paper, I present an easy-to-implement method to test whether jumps detected in a supercooled liquid or glass are renewal events or not. The test was applied to molecular dynamics simulations of a short-chain polymer melt, demonstrating that the jumps can in fact be treated as renewal events. The test further revealed that additional relaxation processes are present which are not accounted for in the CTRW picture, highlighting the limitations of this approach. The notion of renewal events in glass-forming systems could be a very important building block for the interpretation of aging and the glass transition. Furthermore, it could have practical implications for the study of non-equilibrium dynamics in glasses as well as mechanical rejuvenation.

Even hard-sphere colloidal suspensions display Fickian yet non-Gaussian diffusion

We scrutinize three decades of probability density displacement distribution in a simple colloidal suspension with hard-sphere interactions. In this index-matched and density-matched solvent, fluorescent tracer nanoparticles diffuse among matrix particles that are eight times larger, at concentrations from dilute to concentrated, over times up to when the tracer diffuses a few times its size. Displacement distributions of tracers, Gaussian in pure solvent, broaden systematically with increasing obstacle density. The onset of non-Gaussian dynamics is seen in even modestly dilute suspensions, which traditionally would be assumed to follow classic Gaussian expectation. The findings underscore, in agreement with recent studies of more esoteric soft matter systems, the prevalence of non-Gaussian yet Fickian diffusion.

Continuous-time random-walk approach to supercooled liquids. II. Mean-square displacements in polymer melts

The continuous-time random walk (CTRW) describes the single-particle dynamics as a series of jumps separated by random waiting times. This description is applied to analyze trajectories from molecular dynamics (MD) simulations of a supercooled polymer melt. Based on the algorithm presented by Helfferich et al. [Phys. Rev. E 89, 042603 (2014)], we detect jump events of the monomers. As a function of temperature and chain length, we examine key distributions of the CTRW: the jump-length distribution (JLD), the waiting-time distribution (WTD), and the persistence-time distribution (PTD), i.e., the distribution of waiting times for the first jump. For the equilibrium (polymer) liquid under consideration, we verify that the PTD is determined by the WTD. For the mean-square displacement (MSD) of a monomer, the results for the CTRW model are compared with the underlying MD data. The MD data exhibit two regimes of subdiffusive behavior, one for the early α process and another at later times due to chain connectivity. By contrast, the analytical solution of the CTRW yields diffusive behavior for the MSD at all times. Empirically, we can account for the effect of chain connectivity in Monte Carlo simulations of the CTRW. The results of these simulations are then in good agreement with the MD data in the connectivity-dominated regime, but not in the early α regime where they systematically underestimate the MSD from the MD.

Continuous-time random-walk approach to supercooled liquids. I. Different definitions of particle jumps and their consequences

Single-particle trajectories in supercooled liquids display long periods of localization interrupted by "fast moves." This observation suggests a modeling by a continuous-time random walk (CTRW). We perform molecular dynamics simulations of equilibrated short-chain polymer melts near the critical temperature of mode-coupling theory Tc and extract "moves" from the monomer trajectories. We show that not all moves comply with the conditions of a CTRW. Strong forward-backward correlations are found in the supercooled state. A refinement procedure is suggested to exclude these moves from the analysis. We discuss the repercussions of the refinement on the jump-length and waiting-time distributions as well as on characteristic time scales, such as the average waiting time ("exchange time") and the average time for the first move ("persistence time"). The refinement modifies the temperature (T) dependence of these time scales. For instance, the average waiting time changes from an Arrhenius-type to a Vogel-Fulcher-type T dependence. We discuss this observation in the context of the bifurcation of the α process and (Johari) β process found in many glass-forming materials to occur near Tc. Our analysis lays the foundation for a study of the jump-length and waiting-time distributions, their temperature and chain-length dependencies, and the modeling of the monomer dynamics by a CTRW approach in the companion paper [J. Helfferich et al., Phys. Rev. E 89, 042604 (2014)].

Dynamic propensity as an indicator of heterogeneity in room-temperature ionic liquids

Dynamic propensity of an RTIL system exhibits broad and asymmetric distributions, and spatial patterns of the dynamic propensity and mobility distribution are shown.

Distribution of directional change as a signature of complex dynamics

Analyses of random walks traditionally use the mean square displacement (MSD) as an order parameter characterizing dynamics. We show that the distribution of relative angles of motion between successive time intervals of random walks in two or more dimensions provides information about stochastic processes beyond the MSD. We illustrate the behavior of this measure for common models and apply it to experimental particle tracking data. For a colloidal system, the distribution of relative angles reports sensitively on caging as the density varies. For transport mediated by molecular motors on filament networks in vitro and in vivo, we discover self-similar properties that cannot be described by existing models and discuss possible scenarios that can lead to the elucidated statistical features.

Perspective: The glass transition

We provide here a brief perspective on the glass transition field. It is an assessment, written from the point of view of theory, of where the field is and where it seems to be heading. We first give an overview of the main phenomenological characteristics, or "stylised facts," of the glass transition problem, i.e., the central observations that a theory of the physics of glass formation should aim to explain in a unified manner. We describe recent developments, with a particular focus on real space properties, including dynamical heterogeneity and facilitation, the search for underlying spatial or structural correlations, and the relation between the thermal glass transition and athermal jamming. We then discuss briefly how competing theories of the glass transition have adapted and evolved to account for such real space issues. We consider in detail two conceptual and methodological approaches put forward recently, that aim to access the fundamental critical phenomenon underlying the glass transition, be it thermodynamic or dynamic in origin, by means of biasing of ensembles, of configurations in the thermodynamic case, or of trajectories in the dynamic case. We end with a short outlook.

Intracellular transport of insulin granules is a subordinated random walk

We quantitatively analyzed particle tracking data on insulin granules expressing fluorescent fusion proteins in MIN6 cells to better understand the motions contributing to intracellular transport and, more generally, the means for characterizing systems far from equilibrium. Care was taken to ensure that the statistics reflected intrinsic features of the individual granules rather than details of the measurement and overall cell state. We find anomalous diffusion. Interpreting such data conventionally requires assuming that a process is either ergodic with particles working against fluctuating obstacles (fractional brownian motion) or nonergodic with a broad distribution of dwell times for traps (continuous-time random walk). However, we find that statistical tests based on these two models give conflicting results. We resolve this issue by introducing a subordinated scheme in which particles in cages with random dwell times undergo correlated motions owing to interactions with a fluctuating environment. We relate this picture to the underlying microtubule structure by imaging in the presence of vinblastine. Our results provide a simple physical picture for how diverse pools of insulin granules and, in turn, biphasic secretion could arise.

Space-time correlated two-particle hopping in glassy fluids: structural relaxation, irreversibility, decoupling, and facilitation

The microscopic nonlinear Langevin equation (NLE) theory of correlated two-particle dynamics in dense fluids of spherical particles is extended to construct a predictive model of multiple correlated hopping and recaging events of a pair of tagged particles as a function of their initial separation. Modest coarse graining over the liquid structural disorder allows contact to be made with various definitions of irreversible particle motion within the context of a multistate Markov model. The correlated space-time hopping process that underlies structural relaxation can also be analyzed in the context of kinetically constrained models. The dependence of microscopically defined mean persistence and exchange times, their distributions, and relaxation-diffusion decoupling on hard-sphere fluid volume fraction is derived from a model in which irreversible jumps serve as the nucleating persistence event. For a subset of questions, the predictions of the two-particle theory are compared with results from the earlier single-particle NLE approach.