Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Unified view of avalanche criticality in sheared glasses

Plastic events in sheared glasses are considered an example of so-called avalanches, whose sizes obey a power-law probability distribution with the avalanche critical exponent τ. Although the so-called mean-field depinning (MFD) theory predicts a universal value of this exponent, τ_{MFD}=1.5, such a simplification is now known to connote qualitative disagreement with realistic systems. Numerically and experimentally, different values of τ have been reported depending on the literature. Moreover, in the elastic regime, it has been noted that the critical exponent can be different from that in the steady state, and even criticality itself is a matter of debate. Because these confusingly varying results have been reported under different setups, our knowledge of avalanche criticality in sheared glasses is greatly limited. To gain a unified understanding, in this work, we conduct a comprehensive numerical investigation of avalanches in Lennard-Jones glasses under athermal quasistatic shear. In particular, by excluding the ambiguity and arbitrariness that has crept into the conventional measurement schemes, we achieve high-precision measurement and demonstrate that the exponent τ in the steady state follows the prediction of MFD theory, τ_{MFD}=1.5. Our results also suggest that there are two qualitatively different avalanche events. This binariness leads to the nonuniversal behavior of the avalanche size distribution and is likely to be the cause of the varying values of τ reported thus far. To investigate the dependence of criticality and universality on applied shear, we further study the statistics of avalanches in the elastic regime and the ensemble of the first avalanche event in different samples, which provide information about the unperturbed system. We show that while the unperturbed system is indeed off-critical, criticality gradually develops as shear is applied. The degree of criticality is encoded in the fractal dimension of the avalanches, which starts from zero in the off-critical unperturbed state and saturates in the steady state. Moreover, the critical exponent τ is consistent with the prediction of the MFD τ_{MFD} universally, regardless of the amount of applied shear, once the system becomes critical.

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.

Isomorph theory beyond thermal equilibrium

This paper generalizes isomorph theory to systems that are not in thermal equilibrium. The systems are assumed to be R-simple, i.e., to have a potential energy that as a function of all particle coordinates R obeys the hidden-scale-invariance condition U(R_a) < U(R_b) ⇒ U(λR_a) < U(λR_b). "Systemic isomorphs" are introduced as lines of constant excess entropy in the phase diagram defined by density and systemic temperature, which is the temperature of the equilibrium state point with the average potential energy equal to U(R). The dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, the systemic temperature is equal to the bath temperature and the original isomorph formalism is recovered. The new approach rationalizes within a consistent framework previously published observations of isomorph invariance in simulations involving nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs. This paper relates briefly to granular media, physical aging, and active matter. Finally, we discuss the possibility that the energy unit defining the reduced quantities should be based on the systemic rather than the bath temperature.

Isomorph invariance of dynamics of sheared glassy systems

We study hidden scale invariance in the glassy phase of the Kob-Andersen binary Lennard-Jones system. After cooling below the glass transition, we generate a so-called isomorph from the fluctuations of potential energy and virial in the NVT ensemble: a set of density, temperature pairs for which structure and dynamics are identical when expressed in appropriate reduced units. To access dynamical features, we shear the system using the SLLOD algorithm coupled with Lees-Edwards boundary conditions and study the statistics of stress fluctuations and the particle displacements transverse to the shearing direction. We find good collapse of the statistical data, showing that isomorph theory works well in this regime. The analysis of stress fluctuations, in particular the distribution of stress changes over a given strain interval, allows us to identify a clear signature of avalanche behavior in the form of an exponential tail on the negative side. This feature is also isomorph invariant. The implications of isomorphs for theories of plasticity are discussed briefly.

Perspective: Excess-entropy scaling

This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.

Isomorph theory of physical aging

This paper derives and discusses the configuration-space Langevin equation describing a physically aging R-simple system and the corresponding Smoluchowski equation. Externally controlled thermodynamic variables like temperature, density, and pressure enter the description via the single parameter T_s/T, in which T is the bath temperature and T_s is the "systemic" temperature defined at any time t as the thermodynamic equilibrium temperature of the state point with density ρ(t) and potential energy U(t). In equilibrium, T_s ≅ T with fluctuations that vanish in the thermodynamic limit. In contrast to Tool's fictive temperature and other effective temperatures in glass science, the systemic temperature is defined for any configuration with a well-defined density, even if it is not close to equilibrium. Density and systemic temperature define an aging phase diagram, in which the aging system traces out a curve. Predictions are discussed for aging following various density-temperature and pressure-temperature jumps from one equilibrium state to another, as well as for a few other scenarios. The proposed theory implies that R-simple glass-forming liquids are characterized by the dynamic Prigogine-Defay ratio being equal to unity.

Scaling of slip avalanches in sheared amorphous materials based on large-scale atomistic simulations

Atomistic simulations of binary amorphous systems with over 4 million atoms are performed. Systems of two interatomic potentials of the Lennard-Jones type, LJ12-6 and LJ9-6, are simulated. The athermal quasistatic shearing protocol is adopted, where the shear strain is applied in a stepwise fashion with each step followed by energy minimization. For each avalanche event, the shear stress drop (Δσ), the hydrostatic pressure drop (Δσ_{h}), and the potential energy drop (ΔE) are computed. It is found that, with the avalanche size increasing, the three become proportional to each other asymptotically. The probability distributions of avalanche sizes are obtained and values of scaling exponents fitted. In particular, the distributions follow a power law, P(ΔU)∼ΔU^{-τ}, where ΔU is a measure of avalanche sizes defined based on shear stress drops. The exponent τ is 1.25±0.1 for the LJ12-6 systems, and 1.15±0.1 for the LJ9-6 systems. The value of τ for the LJ12-6 systems is consistent with that from an earlier atomistic simulation study by Robbins et al. [Phys. Rev. Lett. 109, 105703 (2012)]PRLTAO0031-900710.1103/PhysRevLett.109.105703, but the fitted values of other scaling exponents differ, which may be because the shearing protocol used here differs from that in their study.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

Hidden scale invariance in condensed matter

Recent developments show that many liquids and solids have an approximate "hidden" scale invariance that implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics in properly reduced units are invariant to a good approximation. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. Liquids and solids with isomorphs include most or all van der Waals bonded systems and metals, as well as weakly ionic or dipolar systems. On the other hand, systems with directional bonding (hydrogen bonds or covalent bonds) or strong Coulomb forces generally do not exhibit hidden scale invariance. The article reviews the theory behind this picture of condensed matter and the evidence for it coming from computer simulations and experiments.