Depinning with dynamic stress overshoots: a hybrid of critical and pseudohysteretic behavior

Scaling theory for the statistics of slip at frictional interfaces

Slip at a frictional interface occurs via intermittent events. Understanding how these events are nucleated, can propagate, or stop spontaneously remains a challenge, central to earthquake science and tribology. In the absence of disorder, rate-and-state approaches predict a diverging nucleation length at some stress σ^{*}, beyond which cracks can propagate. Here we argue for a flat interface that disorder is a relevant perturbation to this description. We justify why the distribution of slip contains two parts: a power law corresponding to "avalanches" and a "narrow" distribution of system-spanning "fracture" events. We derive novel scaling relations for avalanches, including a relation between the stress drop and the spatial extension of a slip event. We compute the cut-off length beyond which avalanches cannot be stopped by disorder, leading to a system-spanning fracture, and successfully test these predictions in a minimal model of frictional interfaces.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate

Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P_{6} and the collective drift velocity V_{x}, three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P_{6} and V_{x}, as well as the observation of hysteresis of V_{x} while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

How collective asperity detachments nucleate slip at frictional interfaces

Sliding at a quasi-statically loaded frictional interface can occur via macroscopic slip events, which nucleate locally before propagating as rupture fronts very similar to fracture. We introduce a microscopic model of a frictional interface that includes asperity-level disorder, elastic interaction between local slip events, and inertia. For a perfectly flat and homogeneously loaded interface, we find that slip is nucleated by avalanches of asperity detachments of extension larger than a critical radius [Formula: see text] governed by a Griffith criterion. We find that after slip, the density of asperities at a local distance to yielding [Formula: see text] presents a pseudogap [Formula: see text], where θ is a nonuniversal exponent that depends on the statistics of the disorder. This result makes a link between friction and the plasticity of amorphous materials where a pseudogap is also present. For friction, we find that a consequence is that stick-slip is an extremely slowly decaying finite-size effect, while the slip nucleation radius [Formula: see text] diverges as a θ-dependent power law of the system size. We discuss how these predictions can be tested experimentally.

Depinning and nonequilibrium dynamic phases of particle assemblies driven over random and ordered substrates: a review

We review the depinning and nonequilibrium phases of collectively interacting particle systems driven over random or periodic substrates. This type of system is relevant to vortices in type-II superconductors, sliding charge density waves, electron crystals, colloids, stripe and pattern forming systems, and skyrmions, and could also have connections to jamming, glassy behaviors, and active matter. These systems are also ideal for exploring the broader issues of characterizing transient and steady state nonequilibrium flow phases as well as nonequilibrium phase transitions between distinct dynamical phases, analogous to phase transitions between different equilibrium states. We discuss the differences between elastic and plastic depinning on random substrates and the different types of nonequilibrium phases which are associated with specific features in the velocity-force curves, fluctuation spectra, scaling relations, and local or global particle ordering. We describe how these quantities can change depending on the dimension, anisotropy, disorder strength, and the presence of hysteresis. Within the moving phase we discuss how there can be a transition from a liquid-like state to dynamically ordered moving crystal, smectic, or nematic states. Systems with periodic or quasiperiodic substrates can have multiple nonequilibrium second or first order transitions in the moving state between chaotic and coherent phases, and can exhibit hysteresis. We also discuss systems with competing repulsive and attractive interactions, which undergo dynamical transitions into stripes and other complex morphologies when driven over random substrates. Throughout this work we highlight open issues and future directions such as absorbing phase transitions, nonequilibrium work relations, inertia, the role of non-dissipative dynamics such as Magnus effects, and how these results could be extended to the broader issues of plasticity in crystals, amorphous solids, and jamming phenomena.

Frictional dynamics of viscoelastic solids driven on a rough surface

We study the effect of viscoelastic dynamics on the frictional properties of a (mean-field) spring-block system pulled on a rough surface by an external drive. When the drive moves at constant velocity V, two dynamical regimes are observed: at fast driving, above a critical threshold V(c), the system slides at the drive velocity and displays a friction force with velocity weakening. Below V(c) the steady sliding becomes unstable and a stick-slip regime sets in. In the slide-hold-slide driving protocol, a peak of the friction force appears after the hold time and its amplitude increases with the hold duration. These observations are consistent with the frictional force encoded phenomenologically in the rate-and-state equations. Our model gives a microscopical basis for such macroscopic description.

Bulk metallic glasses deform via slip avalanches

For the first time in metallic glasses, we extract both the exponents and scaling functions that describe the nature, statistics, and dynamics of slip events during slow deformation, according to a simple mean field model. We model the slips as avalanches of rearrangements of atoms in coupled shear transformation zones (STZs). Using high temporal resolution measurements, we find the predicted, different statistics and dynamics for small and large slips thereby excluding self-organized criticality. The agreement between model and data across numerous independent measures provides evidence for slip avalanches of STZs as the elementary mechanism of inhomogeneous deformation in metallic glasses.

Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with 10(3)-10(6) particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic deformation is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

Distribution of velocities and acceleration for a particle in Brownian correlated disorder: inertial case

We study the motion of an elastic object driven in a disordered environment in presence of both dissipation and inertia. We consider random forces with the statistics of random walks and reduce the problem to a single degree of freedom. It is the extension of the mean-field Alessandro-Beatrice- Bertotti-Montorsi (ABBM) model in presence of an inertial mass m. While the ABBM model can be solved exactly, its extension to inertia exhibits complicated history dependence due to oscillations and backward motion. The characteristic scales for avalanche motion are studied from numerics and qualitative arguments. To make analytical progress, we consider two variants which coincide with the original model whenever the particle moves only forward. Using a combination of analytical and numerical methods together with simulations, we characterize the distributions of instantaneous acceleration and velocity, and compare them in these three models. We show that for large driving velocity, all three models share the same large-deviation function for positive velocities, which is obtained analytically for small and large m, as well as for m=6/25. The effect of small additional thermal and quantum fluctuations can be treated within an approximate method.

Continuous depinning transition with an unusual hysteresis effect

We identify a strange hysteresis which occurs in models of the depinning transition with both elastic and transient overshoot stresses. This hysteresis occurs generically despite the fact that the phase transition is still second order. We calculate the size of the hysteresis gap exactly in a large class of models. Because it is caused by irrelevant perturbations, the hysteresis does not show up in the field theory of the continuum limit. But since it is a real phenomenon it may be a cause for the hysteresis observed in natural depinning. We therefore discuss its experimental and numerical signatures which include microscopic nucleation.

Higher correlations, universal distributions, and finite size scaling in the field theory of depinning

Recently we constructed a renormalizable field theory up to two loops for the quasistatic depinning of elastic manifolds in a disordered environment. Here we explore further properties of the theory. We show how higher correlation functions of the displacement field can be computed. Drastic simplifications occur, unveiling much simpler diagrammatic rules than anticipated. This is applied to the universal scaled width distribution. The expansion in d=4-epsilon predicts that the scaled distribution coincides to the lowest orders with the one for a Gaussian theory with propagator G(q)=1/q(d+2 zeta), zeta being the roughness exponent. The deviations from this Gaussian result are small and involve higher correlation functions, which are computed here for different boundary conditions. Other universal quantities are defined and evaluated: We perform a general analysis of the stability of the fixed point. We find that the correction-to-scaling exponent is omega=-epsilon and not -epsilon/3 as used in the analysis of some simulations. A more detailed study of the upper critical dimension is given, where the roughness of interfaces grows as a power of a logarithm instead of a pure power.

Transition from stable to unstable growth by an inertial force

We introduce a simple growth model where the growth of the interface is affected by an inertial force and a white noise. The magnitude of the inertial force is controlled by a constant p between 0 and 1. An inertial force increases continuously from 0, as p does from 0 to 1. In our model, the interface starts growing from a flat state. When p<p(c), the interface width in our model increases continuously from 0 as time elapses, but it saturates to a constant value in the long time limit. The saturated values of the interface width are the same for different values of p if p<p(c). When p>p(c), however, the interface width increases continuously without saturation as time elapses. We explain via simple calculation how this interesting phenomenon occurs in our model. We find p(c)=0.5 from the calculation. This critical value is in excellent agreement with the critical value p(c)=0.50(1) found from the simulations of our model.