By using large-scale molecular dynamics simulations, the dynamics of two-dimensional (2D) supercooled liquids turns out to be dependent on the system size, while the size dependence is not pronounced in three-dimensional (3D) systems. It is demonstrated that the strong system-size effect in 2D amorphous systems originates from the enhanced fluctuations at long wavelengths which are similar to those of 2D crystal phonons. This observation is further supported by the frequency dependence of the vibrational density of states, consisting of the Debye approximation in the low-wave-number limit. However, the system-size effect in the intermediate scattering function becomes negligible when the length scale is larger than the vibrational amplitude. This suggests that the finite-size effect in a 2D system is transient and also that the structural relaxation itself is not fundamentally different from that in a 3D system. In fact, the dynamic correlation lengths estimated from the bond-breakage function, which do not suffer from those enhanced fluctuations, are not size dependent in either 2D or 3D systems.

The violation of the Stokes-Einstein (SE) relation D ~ (η/T)⁻¹ between the shear viscosity η and the translational diffusion constant D at temperature T is of great importance for characterizing anomalous dynamics of supercooled water. Determining which time scales play key roles in the SE violation remains elusive without the measurement of η. We provide comprehensive simulation results of the dynamic properties involving η and D in the TIP4P/2005 supercooled water. This enabled the thorough identification of the appropriate time scales for the SE relation Dη/T. In particular, it is demonstrated that the temperature dependence of various time scales associated with structural relaxation, hydrogen bond breakage, stress relaxation, and dynamic heterogeneities can be definitely classified into only two classes. That is, we propose the generalized SE relations that exhibit “violation” or “preservation.” The classification depends on the examined time scales that are coupled or decoupled with the diffusion. On the basis of the classification, we explain the physical origins of the violation in terms of the increase in the plateau modulus and the nonexponentiality of stress relaxation. This implies that the mechanism of SE violation is attributed to the attained solidity upon supercooling, which is in accord with the growth of non-Gaussianity and spatially heterogeneous dynamics.

We report an extensive and systematic investigation of the multi-point and multi-time correlation functions to reveal the spatio-temporal structures of dynamic heterogeneities in glass-forming liquids. Molecular dynamics simulations are carried out for the supercooled states of various prototype models of glass-forming liquids such as binary Kob–Andersen, Wahnström, soft-sphere, and network-forming liquids. While the first three models act as fragile liquids exhibiting super-Arrhenius temperature dependence in their relaxation times, the last is a strong glass-former exhibiting Arrhenius behavior. First, we quantify the length scale of the dynamic heterogeneities utilizing the four-point correlation function. The growth of the dynamic length scale with decreasing temperature is characterized by various scaling relations that are analogous to the critical phenomena. We also examine how the growth of the length scale depends upon the model employed. Second, the four-point correlation function is extended to a three-time correlation function to characterize the temporal structures of the dynamic heterogeneities based on our previous studies [K. Kim and S. Saito, Phys. Rev. E 79, 060501–R (2009)10.1103/PhysRevE.79.060501; K. Kim and S. Saito, J. Chem. Phys. 133, 044511 (2010)10.1063/1.3464331]. We provide comprehensive numerical results obtained from the three-time correlation function for the above models. From these calculations, we examine the time scale of the dynamic heterogeneities and determine the associated lifetime in a consistent and systematic way. Our results indicate that the lifetime of the dynamical heterogeneities becomes much longer than the α-relaxation time determined from a two-point correlation function in fragile liquids. The decoupling between the two time scales is remarkable, particularly in supercooled states, and the time scales differ by more than an order of magnitude in a more fragile liquid. In contrast, the lifetime is shorter than the α-relaxation time in tetrahedral network-forming strong liquid, even at lower temperatures.

A supercooled liquid generally exhibits marked shear-thinning behavior, but its detailed mechanism remains elusive. Here we study the dynamics of structural rearrangements in supercooled liquids under shear, using two-dimensional (2D) molecular dynamics simulation. To elucidate the relationship between heterogeneous dynamics and the rheological behavior, we extend the four-point correlation function, which has been used for analyzing “dynamic heterogenity” in a quiescent condition, to a system under steady shear. In the Newtonian regime, the rearrangement dynamics is strongly heterogeneous in space, but remains isotropic. Contrary to this, in the non-Newtonian regime, where marked shear-thinning behavior appears, we find a novel dynamic effect: The mobile region tends to form anisotropic “fluidized bands.” This finding suggests a link between nonlinear rheology and inhomogeneization of flow.

We propose a numerical method to simulate electrohydrodynamic phenomena in charged colloidal dispersions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier-Stokes equations so that the electrohydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters; however, our simulation predicts notable deviations in the case of dense dispersions.