Inherent stress correlations in a quiescent two-dimensional liquid: Static analysis including finite-size effects

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.

Isotropic tensor fields in amorphous solids: Correlation functions of displacement and strain tensor fields

Generalizing recent work on isotropic tensor fields in isotropic and achiral condensed matter systems from two to arbitrary dimensions we address both mathematical aspects assuming perfectly isotropic systems and applications focusing on correlation functions of displacement and strain field components in amorphous solids where isotropy may not hold. Various general points are exemplified using simulated polydisperse Lennard-Jones particles. It is shown that the strain components in reciprocal space have essentially a complex circularly symmetric Gaussian distribution albeit weak non-Gaussianity effects become visible for large wave numbers q where also anisotropy effects become relevant. The dynamical strain correlation functions are strongly nonmonotonic with respect to q with a minimum roughly at the breakdown of the continuum limit.

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.

Yielding Is an Absorbing Phase Transition with Vanishing Critical Fluctuations

The yielding transition in athermal complex fluids can be interpreted as an absorbing phase transition between an elastic, absorbing state with high mesoscopic degeneracy and a flowing, active state. We characterize quantitatively this phase transition in an elastoplastic model under fixed applied shear stress, using a finite-size scaling analysis. We find vanishing critical fluctuations of the order parameter (i.e., the shear rate), and relate this property to the convex character of the phase transition (β>1). We locate yielding within a family of models akin to fixed-energy sandpile (FES) models, only with long-range redistribution kernels with zero modes that result from mechanical equilibrium. For redistribution kernels with sufficiently fast decay, this family of models belongs to a short-range universality class distinct from the conserved directed percolation class of usual FES, which is induced by zero modes.

Long-range correlations in elastic moduli and local stresses at the unjamming transition

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

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.

Emergence of rigidity percolation in flowing granular systems

Jammed granular media and glasses exhibit spatial long-range correlations as a result of mechanical equilibrium. However, the existence of such correlations in the flowing matter, where the mechanical equilibrium is unattainable, has remained elusive. Here, we investigate this problem in the context of the percolation of interparticle forces in flowing granular media. We find that the flow rate introduces an effective long-range correlation, which plays the role of a relevant perturbation giving rise to a spectrum of varying exponents on a critical line as a function of the flow rate. Our numerical simulations along with analytical arguments predict a crossover flow rate [Formula: see text] below which the effect of induced disorder is weak and the universality of the force chain structure is shown to be given by the standard rigidity percolation. We also find a power-law behavior for the critical exponents with the flow rate [Formula: see text].

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.

Stress hyperuniformity and transient oscillatory-exponential correlation decay as signatures of strength vs fragility in glasses

We examine and compare the local stress autocorrelation in the inherent states of a fragile and a strong glass: the Kob-Andersen (KA) binary mixture and the Beest-Kramer-Santen model of silica. For both systems, local (domain-averaged) stress fluctuations asymptotically reach the normal inverse-volume decay in the large domain limit; accordingly, the real-space stress autocorrelation presents long-range power law tails. However, in the case of silica, local stress fluctuations display a high degree of hyperuniformity, i.e., their asymptotic (normal) decay is disproportionately smaller than their bond level amplitude. This property causes the asymptotic power law tails of the real-space stress autocorrelation to be swamped, up to very large distances (several nanometers), by an intermediate oscillatory-exponential decay regime. Similar contributions exist in the KA stress autocorrelation, but they never can be considered as dominating the power law decay and fully disappear when stress is coarse-grained beyond one interatomic distance. Our observations document that the relevance of power-law stress correlation may constitute a key discriminating feature between strong and fragile glasses. Meanwhile, they highlight that the notion of local stress in atomistic systems involves by necessity a choice of observation (coarse-graining) scale, the relevant value of which depends, in principle, on both the model and the phenomenon studied.

Elastoplasticity Mediates Dynamical Heterogeneity Below the Mode Coupling Temperature

As liquids approach the glass transition temperature, dynamical heterogeneity emerges as a crucial universal feature of their behavior. Dynamic facilitation, where local motion triggers further motion nearby, plays a major role in this phenomenon. Here we show that long-ranged, elastically mediated facilitation appears below the mode coupling temperature, adding to the short-range component present at all temperatures. Our results suggest deep connections between the supercooled liquid and glass states, and pave the way for a deeper understanding of dynamical heterogeneity in glassy systems.

The role of friction in statistics and scaling laws of avalanches

We investigate statistics and scaling laws of avalanches in two-dimensional frictional particles by numerical simulations. We find that the critical exponent for avalanche size distributions is governed by microscopic friction between the particles in contact, where the exponent is larger and closer to mean-field predictions if the friction coefficient is finite. We reveal that microscopic "slips" between frictional particles induce numerous small avalanches which increase the slope, as well as the power-law exponent, of avalanche size distributions. We also analyze statistics and scaling laws of the avalanche duration and maximum stress drop rates, and examine power spectra of stress drop rates. Our numerical results suggest that the microscopic friction is a key ingredient of mean-field descriptions and plays a crucial role in avalanches observed in real materials.

Frictional Granular Matter: Protocol Dependence of Mechanical Properties

Theoretical treatments of frictional granular matter often assume that it is legitimate to invoke classical elastic theory to describe its coarse-grained mechanical properties. Here, we show, based on experiments and numerical simulations, that this is generically not the case since stress autocorrelation functions decay more slowly than what is expected from elasticity theory. It was theoretically shown that standard elastic decay demands pressure and torque density fluctuations to be normal, with possibly one of them being hyperuniform. However, generic compressed frictional assemblies exhibit abnormal pressure fluctuations, failing to conform with the central limit theorem. The physics of this failure is linked to correlations built in the material during compression from a dilute configuration prior to jamming. By changing the protocol of compression, one can observe different pressure fluctuations, and stress autocorrelations decay at large scales.

Key role of retardation and non-locality in sound propagation in amorphous solids as evidenced by a projection formalism

We investigate acoustic propagation in amorphous solids by constructing a projection formalism based on separating atomic vibrations into two, "phonon" (P) and "non-phonon" (NP), subspaces corresponding to large and small wavelengths. For a pairwise interaction model, we show the existence of a "natural" separation lengthscale, determined by structural disorder, for which the isolated P subspace presents the acoustic properties of a nearly homogenous (Debye-like) elastic continuum, while the NP one encapsulates all small scale non-affinity effects. The NP eigenstates then play the role of dynamical scatterers for the phonons. However, at variance with a conjecture of defect theories, their spectra present a finite low frequency gap, which turns out to lie around the Boson peak frequency, and only a small fraction of them are highly localized. We then show that small scale disorder effects can be rigorously reduced to the existence, in the Navier-like wave equation of the continuum, of a generalized elasticity tensor, which is not only retarded, since scatterers are dynamical, but also non-local. The full neglect of both retardation and non-locality suffices to account for most of the corrections to Born macroscopic moduli. However, these two features are responsible for sound speed dispersion and have quite a significant effect on the magnitude of sound attenuation. Although it remains open how they impact the asymptotic, large wavelength scaling of sound damping, our findings rule out the possibility of representing an amorphous solid by an inhomogeneous elastic continuum with the standard (i.e., local and static) elastic moduli.

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.

Connecting shear localization with the long-range correlated polarized stress fields in granular materials

One long-lasting puzzle in amorphous solids is shear localization, where local plastic deformation involves cooperative particle rearrangements in small regions of a few inter-particle distances, self-organizing into shear bands and eventually leading to the material failure. Understanding the connection between the structure and dynamics of amorphous solids is essential in physics, material sciences, geotechnical and civil engineering, and geophysics. Here we show a deep connection between shear localization and the intrinsic structures of internal stresses in an isotropically jammed granular material subject to shear. Specifically, we find strong (anti)correlations between the micro shear bands and two polarized stress fields along two directions of maximal shear. By exploring the tensorial characteristics and the rotational symmetry of force network, we reveal that such profound connection is a result of symmetry breaking by shear. Finally, we provide the solid experimental evidence of long-range correlated inherent shear stress in an isotropically jammed granular system.

Fluctuating Elasticity Fails to Capture Anomalous Sound Scattering in Amorphous Solids

The fluctuating elasticity (FE) model, introduced phenomenologically and developed by Schirmacher [J. Non-Cryst. Solids 357, 518 (2011)JNCSBJ0022-309310.1016/j.jnoncrysol.2010.07.052], is today the only theoretical framework available to analyze low-temperature elastic acoustic scattering in glasses. Its existing formulations, which neglect the tensorial nature of elasticity and exclude long-range disorder correlations, predict that the acoustic damping coefficients obey the standard Rayleigh scaling law: Γ∼k^{d+1}, with k the acoustic wave vector, in dimension d. However, recent numerical data, supported by the analysis of existing experimental results, show that Γ does not obey this scaling law but Γ∼-k^{d+1}lnk. Here we analyze in detail how a fully tensorial FE model can be constructed as a long wavelength approximation of the elastic response of the discrete, atomistic, problem. We show that, although it incorporates all long-range correlations, it fails to capture the observed damping in two respects: (i) it misses the anomalous scaling, and predicts the standard Rayleigh law; (ii) it grossly underestimates the amplitude of scattering by about 2 orders of magnitude. This brings clear evidence that the small scale nonaffine displacement fields, although not simply reducible to local defects, play a crucial role in acoustic wave scattering and hence cannot be ignored.

Field Theory for Amorphous Solids

Glasses at low temperature fluctuate around their inherent states; glassy anomalies reflect the structure of these states. Recently, there have been numerous observations of long-range stress correlations in glassy materials, from supercooled liquids to colloids and granular materials, but without a common explanation. Herein it is shown, using a field theory of inherent states, that long-range stress correlations follow from mechanical equilibrium alone, with explicit predictions for stress correlations in two and three dimensions. "Equations of state" relating fluctuations to imposed stresses are derived, as well as field equations that fix the spatial structure of stresses in arbitrary geometries. Finally, a new holographic quantity in 3D amorphous systems is identified.

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.