Long-range spatial correlations of particle displacements and the emergence of elasticity

Topological linking determines elasticity in limited valence networks

Understanding the relationship between the microscopic structure and topology of a material and its macroscopic properties is a fundamental challenge across a wide range of systems. Here we investigate the viscoelasticity of DNA nanostar hydrogels-a model system for physical networks with limited valence-by coupling rheology measurements, confocal imaging and molecular dynamics simulations. We discover that these networks display a large degree of interpenetration and that loops within the network are topologically linked, forming a percolating network-within-network structure. Below the overlapping concentration, the fraction of branching points and the pore size determine the high-frequency elasticity of these physical gels. At higher concentrations, we discover that this elastic response is dictated by the abundance of topological links between looped motifs in the gel. Our findings highlight the emergence of 'topological elasticity' as a previously overlooked mechanism in generic network-forming liquids and gels and inform the design of topologically controllable material behaviours.

Displacement correlations in a two-dimensional colloidal liquid and their relationship with shear strain correlations

Correlations of the displacement field in a two-dimensional model colloidal liquid are studied numerically and analytically. By calculating the displacement correlations and the shear strain correlations from the numerical data of particle simulations, the displacement field is shown to have nontrivial correlations, even in liquids that are only slightly glassy with the area fraction as low as 0.5. It is suggested analytically and demonstrated numerically that the displacement correlations are more informative than the shear correlations: the former behaves logarithmically with regard to the spatial distance at shorter scales, while the corresponding information is missing from the shear correlations. The logarithmic behavior of the displacement correlations is interpreted as manifesting a long-lived aspect of the cage effect.

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.

Testing Theories of the Glass Transition with the Same Liquid but Many Kinetic Rules

We study the glass transition by exploring a broad class of kinetic rules that can significantly modify the normal dynamics of supercooled liquids while maintaining thermal equilibrium. Beyond the usual dynamics of liquids, this class includes dynamics in which a fraction (1-f_{R}) of the particles can perform pairwise exchange or "swap moves," while a fraction f_{P} of the particles can move only along restricted directions. We find that (i) the location of the glass transition varies greatly but smoothly as f_{P} and f_{R} change and (ii) it is governed by a linear combination of f_{P} and f_{R}. (iii) Dynamical heterogeneities (DHs) are not governed by the static structure of the material; their magnitude correlates instead with the relaxation time. (iv) We show that a recent theory for temporal growth of DHs based on thermal avalanches holds quantitatively throughout the (f_{R},f_{P}) diagram. These observations are negative items for some existing theories of the glass transition, particularly those reliant on growing thermodynamic order or locally favored structure, and open new avenues to test other approaches, as we illustrate.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

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.

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.

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.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

Microscopic ergodicity breaking governs the emergence and evolution of elasticity in glass-forming nanoclay suspensions

We report a study combining x-ray photon correlation spectroscopy (XPCS) with in situ rheology to investigate the microscopic dynamics and mechanical properties of aqueous suspensions of the synthetic hectorite clay Laponite, which is composed of charged, nanometer-scale, disk-shaped particles. The suspensions, with particle concentrations ranging from 3.25 to 3.75 wt %, evolve over time from a fluid to a soft glass that displays aging behavior. The XPCS measurements characterize the localization of the particles during the formation and aging of the soft-glass state. The fraction of localized particles, f_{0}, increases rapidly during the early formation stage and grows more slowly during subsequent aging, while the characteristic localization length r_{loc} steadily decreases. Despite the strongly varying rates of aging at different concentrations, both f_{0} and r_{loc} scale with the elastic shear modulus G^{'} in a manner independent of concentration. During the later aging stage, the scaling between r_{loc} and G^{'} agrees quantitatively with a prediction of naive mode coupling theory. Breakdown of agreement with the theory during the early formation stage indicates the prevalence of dynamic heterogeneity, suggesting the soft solid forms through precursors of dynamically localized clusters.

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

A simple and rigorous approach to obtain stress correlations in viscoelastic liquids (including supercooled liquid and equilibrium amorphous systems) is proposed. The long-range dynamical correlations of local shear stress are calculated and analyzed in 2-dimensional space. It is established how the long-range character of the stress correlations gradually emerges as the relevant dynamical correlation length l grows in time. The correlation range l is defined by momentum propagation due to acoustic waves and vorticity diffusion which are the basic mechanisms for transmission of shear stress perturbations. We obtain the general expression defining the time- and distance-dependent stress correlation tensor in terms of material functions (generalized relaxation moduli). The effect of liquid compressibility is quantitatively analyzed; it is shown to be important at large distances and/or short times. The revealed long-range stress correlation effect is shown to be dynamical in nature and unconnected with static structural correlations in liquids (correlation length ξ_s). Our approach is based on the assumption that ξ_s is small enough as reflected in weak wave-number dependencies of the generalized relaxation moduli. We provide a simple physical picture connecting the elucidated long-range fluctuation effect with anisotropic correlations of the (transient) inherent stress field, and discuss its implications.

Stress auto-correlation tensor in glass-forming isothermal fluids: From viscous to elastic response

We develop a generalized hydrodynamic theory, which can account for the build-up of long-ranged and long-lived shear stress correlations in supercooled liquids as the glass transition is approached. Our theory is based on the decomposition of tensorial stress relaxation into fast microscopic processes and slow dynamics due to conservation laws. In the fluid, anisotropic shear stress correlations arise from the tensorial nature of stress. By approximating the fast microscopic processes by a single relaxation time in the spirit of Maxwell, we find viscoelastic precursors of the Eshelby-type correlations familiar in an elastic medium. The spatial extent of shear stress fluctuations is characterized by a correlation length ξ which grows like the viscosity η or time scale τ ∼ η, whose divergence signals the glass transition. In the solid, the correlation length is infinite and stress correlations decay algebraically as r-d in d dimensions.

Collective Motion of Repulsive Brownian Particles in Single-File Diffusion with and without Overtaking

Subdiffusion is commonly observed in liquids with high density or in restricted geometries, as the particles are constantly pushed back by their neighbors. Since this “cage effect” emerges from many-body dynamics involving spatiotemporally correlated motions, the slow diffusion should be understood not simply as a one-body problem but as a part of collective dynamics, described in terms of space–time correlations. Such collective dynamics are illustrated here by calculations of the two-particle displacement correlation in a system of repulsive Brownian particles confined in a (quasi-)one-dimensional channel, whose subdiffusive behavior is known as the single-file diffusion (SFD). The analytical calculation is formulated in terms of the Lagrangian correlation of density fluctuations. In addition, numerical solutions to the Langevin equation with large but finite interaction potential are studied to clarify the effect of overtaking. In the limiting case of the ideal SFD without overtaking, correlated motion with a diffusively growing length scale is observed. By allowing the particles to overtake each other, the short-range correlation is destroyed, but the long-range weak correlation remains almost intact. These results describe nested space–time structure of cages, whereby smaller cages are enclosed in larger cages with longer lifetimes.

Elastic moduli of a Brownian colloidal glass former

The static, dynamic and flow-dependent shear moduli of a binary mixture of Brownian hard disks are studied by an event-driven molecular dynamics simulation. Thereby, the emergence of rigidity close to the glass transition encoded in the static shear modulus [Formula: see text] is accessed by three methods. Results from shear stress auto-correlation functions, elastic dispersion relations, and the elastic response to strain deformations upon the start-up of shear flow are compared. This enables one to sample the time-dependent shear modulus [Formula: see text] consistently over several decades in time. By that a very precise specification of the glass transition point and of [Formula: see text] is feasible. Predictions by mode coupling theory of a finite shear modulus at the glass transition, of α-scaling in fluid states close to the transition, and of shear induced decay in yielding glass states are tested and broadly verified.

Emergence of Long-Ranged Stress Correlations at the Liquid to Glass Transition

A theory for the nonlocal shear stress correlations in supercooled liquids is derived from first principles. It captures the crossover from viscous to elastic dynamics at an idealized liquid to glass transition and explains the emergence of long-ranged stress correlations in glass, as expected from classical continuum elasticity. The long-ranged stress correlations can be traced to the coupling of shear stress to transverse momentum, which is ignored in the classic Maxwell model. To rescue this widely used model, we suggest a generalization in terms of a single relaxation time τ for the fast degrees of freedom only. This generalized Maxwell model implies a divergent correlation length ξ∝τ as well as dynamic critical scaling and correctly accounts for the far-field stress correlations. It can be rephrased in terms of generalized hydrodynamic equations, which naturally couple stress and momentum and furthermore allow us to connect to fluidity and elastoplastic models.

Strain Pattern in Supercooled Liquids

Investigations of strain correlations at the glass transition reveal unexpected phenomena. The shear strain fluctuations show an Eshelby-strain pattern [∼cos(4θ)/r^{2}], characteristic of elastic response, even in liquids, at long times. We address this using a mode-coupling theory for the strain fluctuations in supercooled liquids and data from both video microscopy of a two-dimensional colloidal glass former and simulations of Brownian hard disks. We show that the long-ranged and long-lived strain signatures follow a scaling law valid close to the glass transition. For large enough viscosities, the Eshelby-strain pattern is visible even on time scales longer than the structural relaxation time τ and after the shear modulus has relaxed to zero.

Calculation of displacement correlation tensor indicating vortical cooperative motion in two-dimensional colloidal liquids

As an indicator of cooperative motion in a system of Brownian particles that models two-dimensional colloidal liquids, a displacement correlation tensor is calculated analytically and compared with numerical results. The key idea for the analytical calculation is to relate the displacement correlation tensor, which is a kind of four-point space-time correlation, to the Lagrangian two-time correlation of the deformation gradient tensor. Tensorial treatment of the statistical quantities, including the displacement correlation itself, allows capturing the vortical structure of the cooperative motion. The calculated displacement correlation also implies a negative long-time tail in the velocity autocorrelation, which is a manifestation of the cage effect. Both the longitudinal and transverse components of the displacement correlation are found to be expressible in terms of a similarity variable, suggesting that the cages are nested to form a self-similar structure in the space-time.

Cooling-rate dependence of kinetic and mechanical stabilities of simulated glasses

Glasses created through vapor deposition on a substrate maintained at a proper temperature possess higher kinetic and mechanical stabilities than glasses created by cooling at a constant rate. Molecular dynamics simulations are being increasingly used to understand why vapor deposition improves glasses' stability. There are, however, few detailed molecular dynamics studies of the dependence of the properties of glasses cooled at a constant rate on the rate of cooling. Thus, there is no clear benchmark for comparing ultrastable simulated glasses to simulated glasses prepared through cooling at a constant rate. Here, we examine the dependence of the properties of simulated glasses on the cooling rate used in their preparation. We examine the kinetic stability by measuring the time it takes for a glass to transform back to a liquid upon heating and heterogeneous dynamics during heating. We also examine properties of the energy landscape, and we evaluate mechanical stability by calculating the shear modulus of the glass. The methods outlined here can be used to assess kinetic and mechanical stabilities of simulated glasses generated using specialized algorithms and provide a benchmark for those algorithms.