Long-range stress correlations in viscoelastic and glass-forming fluids

Long ranged stress correlations in the hard sphere liquid

The smooth emergence of shear elasticity is a hallmark of the liquid to glass transition. In a liquid, viscous stresses arise from local structural rearrangements. In the solid, Eshelby has shown that stresses around an inclusion decay as a power law r-D, where D is the dimension of the system. We study glass-forming hard sphere fluids by simulation and observe the emergence of the unscreened power-law Eshelby pattern in the stress correlations of the isotropic liquid state. By a detailed tensorial analysis, we show that the fluctuating force field, viz., the divergence of the stress field, relaxes to zero with time in all states, while the shear stress correlations develop spatial power-law structures inside regions that grow with longitudinal and transverse sound propagation. We observe the predicted exponents r-D and r-D-2. In Brownian systems, shear stresses relax diffusively within these regions, with the diffusion coefficient determined by the shear modulus and the friction coefficient.

General Relations between Stress Fluctuations and Viscoelasticity in Amorphous Polymer and Glass-Forming Systems

Mechanical stress governs the dynamics of viscoelastic polymer systems and supercooled glass-forming fluids. It was recently established that liquids with long terminal relaxation times are characterized by transiently frozen stress fields, which, moreover, exhibit long-range correlations contributing to the dynamically heterogeneous nature of such systems. Recent studies show that stress correlations and relaxation elastic moduli are intimately related in isotropic viscoelastic systems. However, the origin of these relations (involving spatially resolved material relaxation functions) is non-trivial: some relations are based on the fluctuation-dissipation theorem (FDT), while others involve approximations. Generalizing our recent results on 2D systems, we here rigorously derive three exact FDT relations (already established in our recent investigations and, partially, in classical studies) between spatio-temporal stress correlations and generalized relaxation moduli, and a couple of new exact relations. We also derive several new approximate relations valid in the hydrodynamic regime, taking into account the effects of thermal conductivity and composition fluctuations for arbitrary space dimension. One approximate relation was heuristically obtained in our previous studies and verified using our extended simulation data on two-dimensional (2D) glass-forming systems. As a result, we provide the means to obtain, in any spatial dimension, all stress-correlation functions in terms of relaxation moduli and vice versa. The new approximate relations are tested using simulation data on 2D systems of polydisperse Lennard-Jones particles.

Universal stress correlations in crystalline and amorphous packings

We present a universal characterization of stress correlations in athermal systems, across crystalline to amorphous packings. Via numerical analysis of static configurations of particles interacting through harmonic as well as Lennard-Jones potentials, for a variety of preparation protocols and ranges of microscopic disorder, we show that the properties of the stress correlations at large lengthscales are surprisingly universal across all situations, independent of structural correlations, or the correlations in orientational order. In the near-crystalline limit, we present exact results for the stress correlations for both models, which work surprisingly well at large lengthscales, even in the amorphous phase. Finally, we study the differences in stress fluctuations across the amorphization transition, where stress correlations reveal the loss of periodicity in the structure at short lengthscales with increasing disorder.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

Strain correlation functions in isotropic elastic bodies: large wavelength limit for two-dimensional systems

Strain correlation functions in two-dimensional isotropic elastic bodies are shown both theoretically (using the general structure of isotropic tensor fields) and numerically (using a glass-forming model system) to depend on the coordinates of the field variable (position vector r in real space or wavevector q in reciprocal space) and thus on the direction of the field vector and the orientation of the coordinate system. Since the fluctuations of the longitudinal and transverse components of the strain field in reciprocal space are known in the long-wavelength limit from the equipartition theorem, all components of the correlation function tensor field are imposed and no additional physical assumptions are needed. An observed dependence on the field vector direction thus cannot be used as an indication for anisotropy or for a plastic rearrangement. This dependence is different for the associated strain response field containing also information on the localized stress perturbation.

Correlations of tensor field components in isotropic systems with an application to stress correlations in elastic bodies

Correlation functions of components of second-order tensor fields in isotropic systems can be reduced to an isotropic fourth-order tensor field characterized by a few invariant correlation functions (ICFs). It is emphasized that components of this field depend in general on the coordinates of the field vector variable and thus on the orientation of the coordinate system. These angular dependencies are distinct from those of ordinary anisotropic systems. As a simple example of the procedure to obtain the ICFs we discuss correlations of time-averaged stresses in isotropic glasses where only one ICF in reciprocal space becomes a finite constant e for large sampling times and small wave vectors. It is shown that e is set by the typical size of the frozen-in stress components normal to the wave vectors, i.e., it is caused by the symmetry breaking of the stress for each independent configuration. Using the presented general mathematical formalism for isotropic tensor fields this finding explains in turn the observed long-range stress correlations in real space. Under additional but rather general assumptions e is shown to be given by a thermodynamic quantity, the equilibrium Young modulus E. We thus relate for certain isotropic amorphous bodies the existence of finite Young or shear moduli to the symmetry breaking of a stress component in reciprocal space.

Different types of spatial correlation functions for non-ergodic stochastic processes of macroscopic systems

Focusing on non-ergodic macroscopic systems, we reconsider the variances [Formula: see text] of time averages [Formula: see text] of time-series [Formula: see text]. The total variance [Formula: see text] (direct average over all time series) is known to be the sum of an internal variance [Formula: see text] (fluctuations within the meta-basins) and an external variance [Formula: see text] (fluctuations between meta-basins). It is shown that whenever [Formula: see text] can be expressed as a volume average of a local field [Formula: see text] the three variances can be written as volume averages of correlation functions [Formula: see text], [Formula: see text] and [Formula: see text] with [Formula: see text]. The dependences of the [Formula: see text] on the sampling time [Formula: see text] and the system volume V can thus be traced back to [Formula: see text] and [Formula: see text]. Various relations are illustrated using lattice spring models with spatially correlated spring constants. .

Molecular simulations and hydrodynamic theory of nonlocal shear stresscorrelations in supercooled fluids

A supercooled fluid close to the glass transition develops nonlocal shear stress correlations that anticipate the emergence of elasticity. We performed molecular dynamics simulations of a binary Lennard-Jones mixture at different temperatures and investigated the spatiotemporal autocorrelation function of the shear stressfor different wavevectors, q, from a locally measured and Fourier-transformed stress tensor. Anisotropic correlations are observed at non-zero wavevectors, exhibiting strongly damped oscillations with a characteristic frequency ω(q). A comparison with a recently developed hydrodynamic theory [Maier et al., Phys. Rev. Lett. 119, 265701 (2017)] shows a remarkably good quantitative agreement between the particle-based simulations and the theoretical predictions.

Theory of length-scale dependent relaxation moduli and stress fluctuations in glass-forming and viscoelastic liquids

The spatiotemporal correlations of the local stress tensor in supercooled liquids are studied both theoretically and by molecular dynamics simulations of a two-dimensional (2D) polydisperse Lennard-Jones system. Asymptotically exact theoretical equations defining the dynamical structure factor and all components of the stress correlation tensor for low wave-vector q are presented in terms of the generalized (q-dependent) shear and longitudinal relaxation moduli, G(q, t) and K(q, t). We developed a rigorous approach (valid for low q) to calculate K(q, t) in terms of certain bulk correlation functions (for q = 0), the static structure factor S(q), and thermal conductivity κ. The proposed approach takes into account both the thermostatting effect and the effect of polydispersity. The theoretical results for the (q, t)-dependent stress correlation functions are compared with our simulation data, and an excellent agreement is found for qb̄≲0.5 (with b̄ being the mean particle diameter) both above and below the glass transition without any fitting parameters. Our data are consistent with recently predicted (both theoretically and by simulations) long-range correlations of the shear stress quenched in heterogeneous glassy structures.

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.

Simple models for strictly non-ergodic stochastic processes of macroscopic systems

We investigate simple models for strictly non-ergodic stochastic processes [Formula: see text] (t being the discrete time step) focusing on the expectation value v and the standard deviation [Formula: see text] of the empirical variance [Formula: see text] of finite time series [Formula: see text]. [Formula: see text] is averaged over a fluctuating field [Formula: see text] ([Formula: see text] being the microcell position) characterized by a quenched spatially correlated Gaussian field [Formula: see text]. Due to the quenched [Formula: see text]-field [Formula: see text] becomes a finite constant, [Formula: see text], for large sampling times [Formula: see text]. The volume dependence of the non-ergodicity parameter [Formula: see text] is investigated for different spatial correlations. Models with marginally long-ranged [Formula: see text]-correlations are successfully mapped on shear stress data from simulated amorphous glasses of polydisperse beads.

Relaxation moduli of glass-forming systems: temperature effects and fluctuations

Equilibrium and dynamical properties of a two-dimensional polydisperse colloidal model system are characterized by means of molecular dynamics (MD) and Monte Carlo (MC) simulations. We employed several methods to prepare quasi-equilibrated systems: in particular, by slow cooling and tempering with MD (method SC-MD), and by tempering with MC dynamics involving swaps of particle diameters (methods Sw-MD, Sw-MC). It is revealed that the Sw-methods are much more efficient for equilibration below the glass transition temperature T_g leading to denser and more rigid systems which show much slower self-diffusion and shear-stress relaxation than their counterparts prepared with the SC-MD method. The shear-stress relaxation modulus G(t) is obtained based on the classical stress-fluctuation relation. We demonstrate that the α-relaxation time τ_α obtained using a time-temperature superposition of G(t) shows a super-Arrhenius behavior with the VFT temperature T₀ well below T_g. We also derive novel rigorous fluctuation relations providing isothermic and adiabatic compression relaxation moduli in the whole time range (including the short-time inertial regime) based on correlation data for thermostatted systems. It is also shown that: (i) the assumption of Gaussian statistics for stress fluctuations leads to accurate predictions of the variances of the fluctuation moduli for both shear (μ_F) and compression (η_F) at T⪆T_g. (ii) The long-time (quasi-static) isothermic and adiabatic moduli increase on cooling faster than the affine compression modulus η_A, and this leads to a monotonic temperature dependence of η_F which is qualitatively different from μ_F(T) showing a maximum near T_g.

Rheology based estimates of self- and collective diffusivities in viscous liquids

The self-diffusion coefficient of viscous liquids is estimated on the basis of a simple analysis of their rheological shear spectra. To this end, the Almond-West approach, previously employed to access single-particle diffusivities in ionic conductors, is generalized for application to molecular dynamics in supercooled liquids. Rheology based estimates, presented for indomethacin, ortho-terphenyl, and trinaphthylbenzene, reveal relatively small, yet systematic differences when compared with diffusivity data directly measured for these highly viscous liquids. These deviations are discussed in terms of mechanical Haven ratios, introduced to quantify the magnitude of collective translational effects that have an impact on the viscous flow.

Fluctuations of non-ergodic stochastic processes

We investigate the standard deviation [Formula: see text] of the variance [Formula: see text] of time series [Formula: see text] measured over a finite sampling time [Formula: see text] focusing on non-ergodic systems where independent "configurations" c get trapped in meta-basins of a generalized phase space. It is thus relevant in which order averages over the configurations c and over time series k of a configuration c are performed. Three variances of [Formula: see text] must be distinguished: the total variance [Formula: see text] and its contributions [Formula: see text], the typical internal variance within the meta-basins, and [Formula: see text], characterizing the dispersion between the different basins. We discuss simplifications for physical systems where the stochastic variable x(t) is due to a density field averaged over a large system volume V. The relations are illustrated for the shear-stress fluctuations in quenched elastic networks and low-temperature glasses formed by polydisperse particles and free-standing polymer films. The different statistics of [Formula: see text] and [Formula: see text] are manifested by their different system-size dependences.

Ensemble fluctuations matter for variances of macroscopic variables

Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average [Formula: see text] and the standard deviation [Formula: see text] of the variance [Formula: see text] of time series [Formula: see text] of a stochastic process x(t) measured over a finite sampling time [Formula: see text]. Assuming a stationary, Gaussian and ergodic process, [Formula: see text] is given by a functional [Formula: see text] of the autocorrelation function h(t). [Formula: see text] is shown to become large and similar to [Formula: see text] if [Formula: see text] corresponds to a fast relaxation process. Albeit [Formula: see text] does not hold in general for non-ergodic systems, the deviations for common systems with many microstates are merely finite-size corrections. Various issues are illustrated for shear-stress fluctuations in simple coarse-grained model systems.

Stress correlation function and linear response of Brownian particles

We determine the non-local stress autocorrelation tensor in an homogeneous and isotropic system of interacting Brownian particles starting from the Smoluchowski equation of the configurational probability density. In order to relate stresses to particle displacements as appropriate in viscoelastic states, we go beyond the usual hydrodynamic description obtained in the Zwanzig-Mori projection-operator formalism by introducing the proper irreducible dynamics following Cichocki and Hess, and Kawasaki. Differently from these authors, we include transverse contributions as well. This recovers the expression for the stress autocorrelation including the elastic terms in solid states as found for Newtonian and Langevin systems, in case that those are evaluated in the overdamped limit. Finally, we argue that the found memory function reduces to the shear and bulk viscosity in the hydrodynamic limit of smooth and slow fluctuations and derive the corresponding hydrodynamic equations.

Emergent solidity of amorphous materials as a consequence of mechanical self-organisation

Amorphous solids have peculiar properties distinct from crystals. One of the most fundamental mysteries is the emergence of solidity in such nonequilibrium, disordered state without the protection by long-range translational order. A jammed system at zero temperature, although marginally stable, has solidity stemming from the space-spanning force network, which gives rise to the long-range stress correlation. Here, we show that such nonlocal correlation already appears at the nonequilibrium glass transition upon cooling. This is surprising since we also find that the system suffers from giant anharmonic fluctuations originated from the fractal-like potential energy landscape. We reveal that it is the percolation of the force-bearing network that allows long-range stress transmission even under such circumstance. Thus, the emergent solidity of amorphous materials is a consequence of nontrivial self-organisation of the disordered mechanical architecture. Our findings point to the significance of understanding amorphous solids and nonequilibrium glass transition from a mechanical perspective.

Theory of applying shear strains from boundary walls: Linear response in glasses

We construct a linear response theory of applying shear deformations from boundary walls in the film geometry in Kubo's theoretical scheme. Our method is applicable to any solids and fluids. For glasses, we assume quasi-equilibrium around a fixed inherent state. Then, we obtain linear-response expressions for any variables including the stress and the particle displacements, even though the glass interior is elastically inhomogeneous. In particular, the shear modulus can be expressed in terms of the correlations between the interior stress and the forces from the walls. It can also be expressed in terms of the inter-particle correlations, as has been shown in the previous literature. Our stress relaxation function includes the effect of the boundary walls and can be used for inhomogeneous flow response. We show the presence of long-ranged, long-lived correlations among the fluctuations of the forces from the walls and the displacements of all the particles in the cell. We confirm these theoretical results numerically in a two-dimensional model glass. As an application, we describe emission and propagation of transverse sounds after boundary wall motions using these time-correlation functions. We also find resonant sound amplification when the frequency of an oscillatory shear approaches that of the first transverse sound mode.