Diverging length scale of the inhomogeneous mode-coupling theory: a numerical investigation

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

Connecting relaxation time to a dynamical length scale in athermal active glass formers

Supercooled liquids display dynamics that are inherently heterogeneous in space. This essentially means that at temperatures below the melting point, particle dynamics in certain regions of the liquid can be orders of magnitude faster than other regions. Often dubbed dynamical heterogeneity, this behavior has fascinated researchers involved in the study of glass transition for over two decades. A fundamentally important question in all glass transition studies is whether one can connect the growing relaxation time to a concomitantly growing length scale. In this paper, we go beyond the realm of ordinary glass forming liquids and study the origin of a growing dynamical length scale ξ in a self-propelled "active" glass former. This length scale, which is constructed using structural correlations, agrees well with the average size of the clusters of slow-moving particles that are formed as the liquid becomes spatially heterogeneous. We further report that the concomitantly growing α-relaxation time exhibits a simple scaling law, τ_{α}∼exp(μξ/T_{eff}), with μ as an effective chemical potential, T_{eff} as the effective temperature, and μξ as the growing free energy barrier for cluster rearrangements. The findings of our study are valid over four decades of persistence times, and hence they could be very useful in understanding the slow dynamics of a generic active liquid such as an active colloidal suspension, or a self-propelled granular medium.

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.

Dynamic length scales in glass-forming liquids: an inhomogeneous molecular dynamics simulation approach

In this work, we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli, G.; et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular, the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.

Quantitative field theory of the glass transition

We develop a full microscopic replica field theory of the dynamical transition in glasses. By studying the soft modes that appear at the dynamical temperature, we obtain an effective theory for the critical fluctuations. This analysis leads to several results: we give expressions for the mean field critical exponents, and we analytically study the critical behavior of a set of four-points correlation functions, from which we can extract the dynamical correlation length. Finally, we can obtain a Ginzburg criterion that states the range of validity of our analysis. We compute all these quantities within the hypernetted chain approximation for the Gibbs free energy, and we find results that are consistent with numerical simulations.

Mode-coupling theory of the glass transition for confined fluids

We present a detailed derivation of a microscopic theory for the glass transition of a liquid enclosed between two parallel walls relying on a mode-coupling approximation. This geometry lacks translational invariance perpendicular to the walls, which implies that the density profile and the density-density correlation function depends explicitly on the distances to the walls. We discuss the residual symmetry properties in slab geometry and introduce a symmetry adapted complete set of two-point correlation functions. Since the currents naturally split into components parallel and perpendicular to the walls the mathematical structure of the theory differs from the established mode-coupling equations in bulk. We prove that the equations for the nonergodicity parameters still display a covariance property similar to bulk liquids.

Mode-coupling theory for the glass transition: test of the convolution approximation for short-range interactions

We reexamine the convolution approximation commonly used in the mode-coupling theory (MCT) of nonergodic states of classical fluids. This approximation concerns the static correlation functions used as input in the MCT treatment of the dynamics. Besides the hard-sphere model, we consider interaction potentials that present a short-range tail, either attractive or repulsive, beyond the hard core. By using accurate static correlation functions obtained from the fundamental measures functional for hard spheres, we show that the role of three-body direct correlations can be more significant than what is inferred from previous simple ansatzs for pure hard spheres. This may in particular impact the location of the glass transition line and the nonergodicity parameter.

Dynamical heterogeneity in a highly supercooled liquid: consistent calculations of correlation length, intensity, and lifetime

We have investigated dynamical heterogeneity in a highly supercooled liquid using molecular-dynamics simulations in three dimensions. Dynamical heterogeneity can be characterized by three quantities: correlation length ξ(4), intensity χ(4), and lifetime τ(hetero). We evaluated all three quantities consistently from a single order parameter. In a previous study [H. Mizuno and R. Yamamoto, Phys. Rev. E 82, 030501(R) (2010)], we examined the lifetime τ(hetero)(t) in two time intervals t = τ(α) and τ(ngp), where τ(α) is the α-relaxation time and τ(ngp) is the time at which the non-Gaussian parameter of the Van Hove self-correlation function is maximized. In the present study, in addition to the lifetime τ(hetero)(t), we evaluated the correlation length ξ(4)(t) and the intensity χ({4)(t) from the same order parameter used for the lifetime τ(hetero)(t). We found that as the temperature decreases, the lifetime τ(hetero)(t) grows dramatically, whereas the correlation length ξ(4)(t) and the intensity χ(4)(t) increase slowly compared to τ(hetero)(t) or plateaus. Furthermore, we investigated the lifetime τ(hetero)(t) in more detail. We examined the time-interval dependence of the lifetime τ(hetero)(t) and found that as the time interval t increases, τ(hetero)(t) monotonically becomes longer and plateaus at the relaxation time of the two-point density correlation function. At the large time intervals for which τ(hetero)(t) plateaus, the heterogeneous dynamics migrate in space with a diffusion mechanism, such as the particle density.

Analysis of a growing dynamic length scale in a glass-forming binary hard-sphere mixture

We examine a length scale that characterizes the spatial extent of heterogeneous dynamics in a glass-forming binary hard-sphere mixture up to the mode-coupling volume fraction ϕ(c). First, we characterize the system's dynamics. Then, we utilize a method [Phys. Rev. Lett. 105, 217801 (2010)] to extract and analyze the ensemble-independent dynamic susceptibility χ(4)(t) and the dynamic correlation length ξ(t) for a range of times between the β and α relaxation times. We find that in this time range the dynamic correlation length follows a volume fraction-independent master curve ξ(t)~ln(t). For longer times, ξ(t) departs from this master curve and remains constant up to the largest time at which we can determine the length accurately. In addition to the previously established correlation τ(α)~exp[ξ(τ(α))] between the α relaxation time, τ(α), and the dynamic correlation length at this time, ξ(τ(α)), we also find a similar correlation for the diffusion coefficient D~exp[ξ(τ(α))(θ)] with θ≈0.6. We discuss the relevance of these findings for different theories of the glass transition.

Dynamic heterogeneity in a glass forming fluid: susceptibility, structure factor, and correlation length

We investigate the growth of dynamic heterogeneity in a glassy hard-sphere mixture for volume fractions up to and including the mode-coupling transition. We use an 80,000 particle system to test a new procedure to evaluate a dynamic correlation length ξ(t): we determine the ensemble independent dynamic susceptibility χ(4)(t) and use it to facilitate evaluation of ξ(t) from the small wave vector behavior of the four-point structure factor. We analyze relations between the α relaxation time τ(α), χ(4)(τ(α)), and ξ(τ(α)). We find that mode-coupling-like power laws provide a reasonable description of the data over a restricted range of volume fractions, but the power laws' exponents differ from those predicted by the inhomogeneous mode-coupling theory. We find ξ(τ(α))~ln(τ(α)) over the full range of volume fractions studied, which is consistent with an Adam-Gibbs-type relation.