Hydrodynamic theory of the glass transition

Q-dependent collective relaxation dynamics of glass-forming liquid Ca0.4K0.6(NO3)1.4 investigated by wide-angle neutron spin-echo

The relaxation behavior of glass formers exhibits spatial heterogeneity and dramatically changes upon cooling towards the glass transition. However, the underlying mechanisms of the dynamics at different microscopic length scales are not fully understood. Employing the recently developed wide-angle neutron spin-echo spectroscopy technique, we measured the Q-dependent coherent intermediate scattering function of a prototypical ionic glass former Ca_0.4K_0.6(NO₃)_1.4, in the highly viscous liquid state. In contrast to the structure modulated dynamics for Q < 2.4 Å⁻¹, i.e., at and below the structure factor main peak, for Q > 2.4 Å⁻¹, beyond the first minimum above the structure factor main peak, the stretching exponent exhibits no temperature dependence and concomitantly the relaxation time shows smaller deviations from Arrhenius behavior. This finding indicates a change in the dominant relaxation mechanisms around a characteristic length of 2π/(2.4 Å⁻¹) ≈ 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

Length scale dependence is important for understanding the collective relaxation dynamics in glass-forming liquids. Here, the authors find in liquid Ca_0.4K_0.6(NO₃)_1.4 a change in the dominant relaxation mechanisms around 2.6 Å, below which the relaxation process exhibits a temperature independent distribution and more Arrhenius-like behavior.

Qualitatively different collective and single-particle dynamics in a supercooled liquid

The equations of fluctuating nonlinear hydrodynamics for a two component mixture are obtained with a proper choice of slow variables which correspond to the conservation laws in the system. Using these nonlinear equations we construct the basic equations of the mode coupling theory (MCT) and consequent ergodic-nonergodic (ENE) transition in a binary mixture. The model is also analyzed in the one component limit of the mixture to study the dynamics of a tagged particle in the sea of identical particles. According to the existing MCT, dynamics of the single-particle correlation is slaved to that of the collective density fluctuations and, hence, both correlations freeze simultaneously at the ENE transition. We show here from a nonperturbative approach that at the ENE transition, characterized by the freezing of the long time limit of the dynamic correlation of collective density fluctuations to a nonzero value, the tagged-particle correlation still decays to zero. Our result implies that the point at which simulation or experimental data of the self-diffusion constant extrapolate to zero would not correspond to the ENE transition of simple MCT.

Numerical study of long-time dynamics and ergodic-nonergodic transitions in dense simple fluids

Since the mid-1980s, mode-coupling theory (MCT) has been the de facto theoretic description of dense fluids and the transition from the fluid state to the glassy state. MCT, however, is limited by the approximations used in its construction and lacks an unambiguous mechanism to institute corrections. We use recent results from a new theoretical framework--developed from first principles via a self-consistent perturbation expansion in terms of an effective two-body potential--to numerically explore the kinetics of systems of classical particles, specifically hard spheres governed by Smoluchowski dynamics. We present here a full solution for such a system to the kinetic equation governing the density-density time correlation function and show that the function exhibits the characteristic two-step decay of supercooled fluids and an ergodic-nonergodic transition to a dynamically arrested state. Unlike many previous numerical studies--and in stark contrast to experiment--we have access to the full time and wave-number range of the correlation function with great precision and are able to track the solution unprecedentedly close to the transition, covering nearly 15 decades in scaled time. Using asymptotic approximation techniques analogous to those developed for MCT, we fit the solution to predicted forms and extract critical parameters. We find complete qualitative agreement with known glassy behavior (e.g. power-law divergence of the α-relaxation time scale in the ergodic phase and square-root growth of the glass form factors in the nonergodic phase), as well as some limited quantitative agreement [e.g. the transition at packing fraction η*=0.60149761(10)], consistent with previous static solutions under this theory and with comparable colloidal suspension experiments. However, most importantly, we establish that this new theory is able to reproduce the salient features seen in other theories, experiments, and simulations but has the advantages of being derived from first principles and possessing a clear mechanism for making systematic corrections.

Renormalized dynamics of the Dean-Kawasaki model

We study the model of a supercooled liquid for which the equation of motion for the coarse-grained density ρ(x,t) is the nonlinear diffusion equation originally proposed by Dean and Kawasaki, respectively, for Brownian and Newtonian dynamics of fluid particles. Using a Martin-Siggia-Rose (MSR) field theory we study the renormalization of the dynamics in a self-consistent form in terms of the so-called self-energy matrix Σ. The appropriate model for the renormalized dynamics involves an extended set of field variables {ρ,θ}, linked through a nonlinear constraint. The latter incorporates, in a nonperturbative manner, the effects of an infinite number of density nonlinearities in the dynamics. We show that the contributing element of Σ which renormalizes the bare diffusion constant D(0) to D(R) is same as that proposed by Kawasaki and Miyazima [Z. Phys. B Condens. Matter 103, 423 (1997)]. D(R) sharply decreases with increasing density. We consider the likelihood of a ergodic-nonergodic (ENE) transition in the model beyond a critical point. The transition is characterized by the long-time limit of the density correlation freezing at a nonzero value. From our analysis we identify an element of Σ which arises from the above-mentioned nonlinear constraint and is key to the viability of the ENE transition. If this self-energy would be zero, then the model supports a sharp ENE transition with D(R)=0 as predicted by Kawasaki and Miyazima. With the full model having nonzero value for this self-energy, the density autocorrelation function decays to zero in the long-time limit. Hence the ENE transition is not supported in the model.

Brownian dynamics: from glassy to trivial

We endow a system of interacting particles with two distinct, local, Markovian, and reversible microscopic dynamics that both converge to the Boltzmann-Gibbs equilibrium of standard liquids. While the first, standard, one leads to glassy dynamics, we use field-theoretical techniques to show that the latter displays no sign of glassiness. The approximations we use, akin to the mode-coupling approximation, are famous for magnifying glassy aspects of the dynamics, supposedly through the neglect of activated events. Despite this, the modified dynamics seem to stick to standard liquid relaxation. This finding singles out as applying to a realistic system of interacting particles in low dimensions and questions the role of the dynamical rules used to explore a given static free-energy landscape. Moreover, our peculiar choice of dynamical rules offers the possibility of a direct connection with replica theory, and our findings therefore call for a clarification of the interplay between replica theory and the underlying dynamics of the system.

Critical dynamics of an isothermal compressible nonideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible nonideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.

Critical test of the mode-coupling theory of the glass transition

In its common implementation, the mode-coupling theory of the glass transition predicts the time evolution of the intermediate scattering functions in viscous liquids on the sole basis of the structural information encoded in two-point density correlations. We provide a critical test of this property and show that the theory fails to describe the strong differences of dynamical behavior seen in two model liquids characterized by very similar pair-correlation functions. Because we use "exact" static information provided by numerical simulations, our results are a direct indication that some important information about the dynamics of viscous liquids is not captured by pair correlations and is thus not described by the mode-coupling theory, even in the temperature regime where the theory is usually applied.

Density nonlinearities in field theories for a toy model of fluctuating nonlinear hydrodynamics of supercooled liquids

We study a zero-dimensional version of the fluctuating nonlinear hydrodynamics (FNH) of supercooled liquids originally investigated by Das and Mazenko (DM) [Shankar P. Das and Gene F. Mazenko Phys. Rev. A 34, 2265 (1986)]. The time-dependent density-like and momentum-like variables are introduced with no spatial degrees of freedom in this toy model. The structure of nonlinearities takes the similar form to the original FNH, which allows one to study in a simpler setting the issues raised recently regarding the field theoretical approaches to glass forming liquids. We study the effects of density nonlinearities on the time evolution of correlation and response functions by developing field theoretic formulations in two different ways: first by following the original prescription of DM and then by constructing a dynamical action which possesses a linear time-reversal symmetry as proposed recently. We show explicitly that, at the one-loop order of the perturbation theory, the DM-type field theory does not support a sharp ergodic-nonergodic transition, while the other admits one. The simple nature of the toy model in the DM formulation allows us to develop numerical solutions to a complete set of coupled dynamical equations for the correlation and response functions at the one-loop order.

Fluctuating nonlinear hydrodynamics does not support an ergodic-nonergodic transition

Despite its appeal, real and simulated glass forming systems do not undergo an ergodic-nonergodic (ENE) transition. We reconsider whether the fluctuating nonlinear hydrodynamics (FNH) model for this system, introduced by us in 1986, supports an ENE transition. Using nonperturbative arguments, with no reference to the hydrodynamic regime, we show that the FNH model does not support an ENE transition. Our results support the findings in the original paper. Assertions in the literature questioning the validity of the original work are shown to be in error.

Random diffusion model

We study the random diffusion model. This is a continuum model for a conserved scalar density field varphi driven by diffusive dynamics. The interesting feature of the dynamics is that the bare diffusion coefficient D is density dependent. In the simplest case, D=D[over ]+D_{1}deltavarphi , where D[over ] is the constant average diffusion constant. In the case where the driving effective Hamiltonian is quadratic, the model can be treated using perturbation theory in terms of the single nonlinear coupling D1 . We develop perturbation theory to fourth order in D1 . The are two ways of analyzing this perturbation theory. In one approach, developed by Kawasaki, at one-loop order one finds mode-coupling theory with an ergodic-nonergodic transition. An alternative more direct interpretation at one-loop order leads to a slowing down as the nonlinear coupling increases. Eventually one hits a critical coupling where the time decay becomes algebraic. Near this critical coupling a weak peak develops at a wave number well above the peak at q=0 associated with the conservation law. The width of this peak in Fourier space decreases with time and can be identified with a characteristic kinetic length which grows with a power law in time. For stronger coupling the system becomes metastable and then unstable. At two-loop order it is shown that the ergodic-nonergodic transition is not supported. It is demonstrated that the critical properties of the direct approach survive, going to higher order in perturbation theory.

Effective separation of forces in a mode coupling theory of self-diffusion

A mode coupling theory (MCT) expression for the self-diffusion coefficient follows simply when the soft fluctuating intermolecular forces are projected along a collective densitylike variable. The projected forces separate into two parts: from the gradient of the direct correlation function (dcf), and from the short range forces. The time correlation function of the dcf-derived forces is related to the excess entropy, as shown by Ali [J. Chem. Phys. 124, 144504 (2006)], and this relationship is evaluated for two variations of MCT. As for hard spheres, the derivation of an analogous MCT is beset by a number of singularities that kinetic theory could not remove. A justifiable MCT for hard sphere fluids may not exist.

Do current-density nonlinearities cut off the glass transition?

Extended mode-coupling theories for dense fluids predict that nonlinear current-density couplings cut off the singular "ideal glass transition" present in the standard mode-coupling theory where such couplings are ignored. We suggest here that, rather than allowing for activated processes as sometimes supposed, contributions from current-density couplings are always negligible close to a glass transition. We discuss in schematic terms how activated processes can nonetheless cut off the transition by causing the memory function to become linear in correlators at late times.

Supercooled liquids under shear: theory and simulation

We analyze the behavior of supercooled fluids under shear both theoretically and numerically. Theoretically, we generalize the mode-coupling theory of supercooled fluids to systems under stationary shear flow. Our starting point is the set of generalized fluctuating hydrodynamic equations with a convection term. A nonlinear integrodifferential equation for the intermediate scattering function is constructed. This theory is applied to a two-dimensional colloidal suspension. The shear rate dependence of the intermediate scattering function and the shear viscosity is analyzed. We have also performed extensive numerical simulations of a two-dimensional binary liquid with soft-core interactions near, but above, the glass transition temperature. Both theoretical and numerical results show the following. (i) A drastic reduction of the structural relaxation time and the shear viscosity due to shear. Both the structural relaxation time and the viscosity decrease as gamma(-nu) with an exponent nu< or =1, where gamma; is the shear rate. (ii) Almost isotropic dynamics regardless of the strength of the anisotropic shear flow.

Self-diffusion in a system of interacting Langevin particles

The behavior of the self-diffusion constant of Langevin particles interacting via a pairwise interaction is considered. The diffusion constant is calculated approximately within a perturbation theory in the potential strength about the bare diffusion constant. It is shown how this expansion leads to a systematic double expansion in the inverse temperature beta and the particle density rho. The one-loop diagrams in this expansion can be summed exactly and we show that this result is exact in the limit of small beta and rhobeta constants. The one-loop result can also be resummed using a semiphenomenological renormalization group method which has proved useful in the study of diffusion in random media. In certain cases the renormalization group calculation predicts the existence of a diverging relaxation time signaled by the vanishing of the diffusion constant, possible forms of divergence coming from this approximation are discussed. Finally, at a more quantitative level, the results are compared with numerical simulations, in two dimensions, of particles interacting via a soft potential recently used to model the interaction between coiled polymers.