Significant difference in the dynamics between strong and fragile glass formers

Transverse Viscous Transport in Classical Solid States

The transverse velocity time correlation function C[over ˜]_{T}(k,ω) with k and ω being the wave number and the frequency, respectively, is a fundamental quantity in determining the transverse mechanical and transport properties of materials. In ordinary liquids, a nonzero value of C[over ˜]_{T}(k,0) is inevitably linked to viscous material flows. Even in solids where significant material flows are precluded due to almost frozen positional degrees of freedom, our molecular dynamics simulations reveal that C[over ˜]_{T}(k,0) takes a nonzero value, whereby the time integration of the velocity field shows definite diffusive behavior with diffusivity C[over ˜]_{T}(k,0)/3. This behavior is attributed to viscous transport accompanying a small random convection of the velocity field (the inertia effect), and the resultant viscosity is measurable in the Eulerian description: the constituent particles that substantially carry momenta fluctuate slightly around their reference positions. In the Eulerian description, the velocity field is explicitly associated with such fluctuating instantaneous particle positions, whereas in the Lagrangian description, this is not the case. The present study poses a fundamental problem for continuum mechanics: reconciling liquid and solid descriptions in the limit of the infinite structural relaxation time.

Stress hyperuniformity and transient oscillatory-exponential correlation decay as signatures of strength vs fragility in glasses

We examine and compare the local stress autocorrelation in the inherent states of a fragile and a strong glass: the Kob-Andersen (KA) binary mixture and the Beest-Kramer-Santen model of silica. For both systems, local (domain-averaged) stress fluctuations asymptotically reach the normal inverse-volume decay in the large domain limit; accordingly, the real-space stress autocorrelation presents long-range power law tails. However, in the case of silica, local stress fluctuations display a high degree of hyperuniformity, i.e., their asymptotic (normal) decay is disproportionately smaller than their bond level amplitude. This property causes the asymptotic power law tails of the real-space stress autocorrelation to be swamped, up to very large distances (several nanometers), by an intermediate oscillatory-exponential decay regime. Similar contributions exist in the KA stress autocorrelation, but they never can be considered as dominating the power law decay and fully disappear when stress is coarse-grained beyond one interatomic distance. Our observations document that the relevance of power-law stress correlation may constitute a key discriminating feature between strong and fragile glasses. Meanwhile, they highlight that the notion of local stress in atomistic systems involves by necessity a choice of observation (coarse-graining) scale, the relevant value of which depends, in principle, on both the model and the phenomenon studied.

Negative density-dependence of the structural relaxation time of liquid silica: insights from a comparative molecular dynamics study

In many tetrahedral network-forming liquids, structural relaxation is anomalously accelerated by compression over relatively low pressure ranges. Here, for silica, we study this problem through comparative molecular dynamics simulations using two different models. Under compression, the network structures are compacted by slight tuning of the intertetrahedral bond angles while nearly preserving the unit tetrahedral structure. The consequent structural changes are remarkable for length scales larger than the nearest neighbor ion-pair distances. Accompanying with such structural changes, the interactions of the nearest Si-O pairs remain almost unchanged, whereas those of other ion pairs are, on average, strengthened by the degree of compression. In particular, the enhancement of the net Si-O interactions at the next nearest neighbor distance, which assist an ion in escaping from the potential well, reduces the activation energy, leading to a significant acceleration of structural relaxation. The results of our comparative molecular dynamics simulations are compatible with the scenario proposed by Angell, and further indicate that the structural relaxation dynamics cannot be uniquely determined by the configurations but strongly depends on the details of the coupling between the structure and the interaction.

Role of hydrodynamics in liquid-liquid transition of a single-component substance

Liquid-liquid transition (LLT) is an unconventional transition between two liquid states in a single-component system. This phenomenon has recently attracted considerable attention not only because of its counterintuitive nature but also since it is crucial for our fundamental understanding of the liquid state. However, its physical understanding has remained elusive, particularly of the critical dynamics and phase-ordering kinetics. So far, the hydrodynamic degree of freedom, which is the most intrinsic kinetic feature of liquids, has been neglected in its theoretical description. Here we develop a Ginzburg-Landau-type kinetic theory of LLT taking it into account, based on a two-order parameter model. We examine slow critical fluctuations of the nonconserved order parameter coupled to the hydrodynamic degree of freedom in equilibrium. We also study the nonequilibrium process of LLT. We show both analytically and numerically that domain growth becomes faster (slower), depending upon the density decrease (increase) upon the transition, as a consequence of hydrodynamic flow induced by the density change. The coupling between nonconserved order parameter and hydrodynamic interaction results in anomalous domain growth in both nucleation-growth-type and spinodal-decomposition-type LLT. Our study highlights the characteristic features of hydrodynamic fluctuations and phase ordering during LLT under complex interplay among conserved and nonconserved order parameters and the hydrodynamic transport intrinsic to the liquid state.

Growing length scale accompanying vitrification: A perspective based on nonsingular density fluctuations

A model for describing growing length scale accompanying the vitrification is introduced, in which we assume that in a subsystem whose density is above a certain threshold value, ρ_{c}, due to steric constraints, particle rearrangements are highly suppressed for a sufficiently long time period (∼structural relaxation time). We regard such a subsystem as a glassy cluster. With this assumption and without introducing any complicated thermodynamic arguments, we predict that with compression (increasing average density ρ) at a fixed temperature T in supercooled states, the characteristic length of the clusters, ξ, diverges as ξ∼(ρ_{c}-ρ)^{-2/d}, where d is the spatial dimensionality. This ξ measures the average persistence length of the steric constraints in blocking the rearrangement motions and is determined by the subsystem density. Additionally, with decreasing T at a fixed ρ, the length scale diverges in the same manner as ξ∼(T-T_{c})^{-2/d}, for which ρ is identical to ρ_{c} at T=T_{c}. The exponent describing the diverging length scale is the same as the one predicted by some theoretical models and indeed has been observed in some simulations and experiments. However, the basic mechanism for this divergence is different; that is, we do not invoke thermodynamic anomalies associated with the thermodynamic phase transition as the origin of the growing length scale. We further present arguements for the cooperative properties of the structural relaxation based on the clusters.

Universal glass-forming behavior of in vitro and living cytoplasm

Physiological processes in cells are performed efficiently without getting jammed although cytoplasm is highly crowded with various macromolecules. Elucidating the physical machinery is challenging because the interior of a cell is so complex and driven far from equilibrium by metabolic activities. Here, we studied the mechanics of in vitro and living cytoplasm using the particle-tracking and manipulation technique. The molecular crowding effect on cytoplasmic mechanics was selectively studied by preparing simple in vitro models of cytoplasm from which both the metabolism and cytoskeletons were removed. We obtained direct evidence of the cytoplasmic glass transition; a dramatic increase in viscosity upon crowding quantitatively conformed to the super-Arrhenius formula, which is typical for fragile colloidal suspensions close to jamming. Furthermore, the glass-forming behaviors were found to be universally conserved in all the cytoplasm samples that originated from different species and developmental stages; they showed the same tendency for diverging at the macromolecule concentrations relevant for living cells. Notably, such fragile behavior disappeared in metabolically active living cells whose viscosity showed a genuine Arrhenius increase as in typical strong glass formers. Being actively driven by metabolism, the living cytoplasm forms glass that is fundamentally different from that of its non-living counterpart.