Dynamics and correlation length scales of a glass-forming liquid in quiescent and sheared conditions

We numerically study dynamics and correlation length scales of a colloidal liquid in both quiescent and sheared conditions to further understand the origin of slow dynamics and dynamic heterogeneity in glass-forming systems. The simulation is performed in a weakly frustrated two-dimensional liquid, where locally preferred order is allowed to develop with increasing density. The four-point density correlations and bond-orientation correlations, which have been frequently used to capture dynamic and static length scales $\xi$ in a quiescent condition, can be readily extended to a system under steady shear in this case. In the absence of shear, we confirmed the previous findings that the dynamic slowing down accompanies the development of dynamic heterogeneity. The dynamic and static length scales increase with $\alpha$-relaxation time $\tau_{\alpha}$ as power-law $\xi\sim\tau_{\alpha}^{\mu}$ with $\mu>0$. In the presence of shear, both viscosity and $\tau_{\alpha}$ have power-law dependence on shear rate in the marked shear thinning regime. However, dependence of correlation lengths cannot be described by power laws in the same regime. Furthermore, the relation $\xi\sim\tau_{\alpha}^{\mu}$ between length scales and dynamics holds for not too strong shear where thermal fluctuations and external forces are both important in determining the properties of dense liquids. Thus, our results demonstrate a link between slow dynamics and structure in glass-forming liquids even under nonequilibrium conditions.Comment: 9 pages, 17 figures. Accepted by J. Phys.: Condens. Matte

Heterogeneity in structurally arrested hard spheres

When cooled or compressed sufficiently rapidly, a liquid vitrifies into a glassy amorphous state. Vitrification in a dense liquid is associated with jamming of the particles. For hard spheres, the density and degree of order in the final structure depend on the compression rate: simple intuition suggests, and previous computer simulation demonstrates, that slower compression results in states that are both denser and more ordered. In this work, we use the Lubachevsky-Stillinger algorithm to generate a sequence of structurally arrested hard-sphere states by varying the compression rate. We find that while the degree of order, as measured by both bond-orientation and translation order parameters, increases monotonically with decreasing compression rate, the density of the arrested state first increases, then decreases, then increases again, as the compression rate decreases, showing a minimum at an intermediate compression rate. Examination of the distribution of the local order parameters and the distribution of the root-mean-square fluctuation of the particle positions, as well as direct visual inspection of the arrested structures, reveal that they are structurally heterogeneous, consisting of disordered, amorphous regions and locally ordered crystal-like domains. In particular, the low-density arrested states correspond with many interconnected small crystal clusters that form a polycrystalline network interspersed in an amorphous background, suggesting that jamming by the domains may be an important mechanism for these states

Structure, compressibility factor and dynamics of highly size-asymmetric binary hard-disk liquids

By using event-driven molecular dynamics simulation, we investigate effects of varying the area fraction of the smaller component on structure, compressibility factor and dynamics of the highly size-asymmetric binary hard-disk liquids. We find that the static pair correlations of the large disks are only weakly perturbed by adding small disks. The higher-order static correlations of the large disks, by contrast, can be strongly affected. The compressibility factor of the system first decreases and then increases upon increasing the area fraction of the small disks and separating different contributions to it allows to rationalize this non-monotonic phenomenon. Furthermore, adding small disks can influence dynamics of the system in quantitative and qualitative ways. For the large disks, the structural relaxation time increases monotonically with increasing the area fraction of the small disks at low and moderate area fractions of the large disks. In particular, "reentrant" behavior appears at sufficiently high area fractions of the large disks, strongly resembling the reentrant glass transition in short-ranged attractive colloids and the inverted glass transition in binary hard spheres with large size disparity. By tuning the area fraction of the small disks, relaxation process for the small disks shows concave-to-convex crossover and logarithmic decay behavior, as found in other binary mixtures with large size disparity. Moreover, diffusion of both species is suppressed by adding small disks. Long-time diffusion for the small disks shows power-law-like behavior at sufficiently high area fractions of the small disks, which implies precursors of a glass transition for the large disks and a localization transition for the small disks. Therefore, our results demonstrate the generic dynamic features in highly size-asymmetric binary mixtures.Comment: 9 pages, 12 figure