Onset of irreversibility and chaos in amorphous solids under periodic shear

Interplay between an Absorbing Phase Transition and Synchronization in a Driven Granular System

Absorbing phase transitions (APTs) are widespread in nonequilibrium systems, spanning condensed matter, epidemics, earthquakes, ecology, and chemical reactions. APTs feature an absorbing state in which the system becomes entrapped, along with a transition, either continuous or discontinuous, to an active state. Understanding which physical mechanisms determine the order of these transitions represents a challenging open problem in nonequilibrium statistical mechanics. Here, by numerical simulations and mean-field analysis, we show that a quasi-2D vibrofluidized granular system exhibits a novel form of APT. The absorbing phase is observed in the horizontal dynamics below a critical packing fraction, and can be continuous or discontinuous based on the emergent degree of synchronization in the vertical motion. Our results provide a direct representation of a feasible experimental scenario, showcasing a surprising interplay between dynamic phase transition and synchronization.

Slow Fatigue and Highly Delayed Yielding via Shear Banding in Oscillatory Shear

We study theoretically the dynamical process of yielding in cyclically sheared amorphous materials, within a thermal elastoplastic model and the soft glassy rheology model. Within both models we find an initially slow accumulation, over many cycles after the inception of shear, of low levels of damage in the form strain heterogeneity across the sample. This slow fatigue then suddenly gives way to catastrophic yielding and material failure. Strong strain localization in the form of shear banding is key to the failure mechanism. We characterize in detail the dependence of the number of cycles N^{*} before failure on the amplitude of imposed strain, the working temperature, and the degree to which the sample is annealed prior to shear. We discuss our finding with reference to existing experiments and particle simulations, and suggest new ones to test our predictions.

Loss of memory of an elastic line on its way to limit cycles

Oscillatory-driven amorphous materials forget their initial configuration and converge to limit cycles. Here we investigate this memory loss under a nonquasistatic drive in a minimal model system, with quenched disorder and memory encoded in a spatial pattern, where oscillating protocols are formally replaced by a positive-velocity drive. We consider an elastic line driven athermally in a quenched disorder with biperiodic boundary conditions and tunable system size, thus controlling the area swept by the line per cycle as would the oscillation amplitude. The convergence to disorder-dependent limit cycle is strongly coupled to the nature of its velocity dynamics depending on system size. Based on the corresponding phase diagram, we propose a generic scenario for memory formation in disordered systems under finite driving rate.

Reversible, irreversible, and mixed regimes for periodically driven disks in random obstacle arrays

We examine an assembly of repulsive disks interacting with a random obstacle array under a periodic drive and find a transition from reversible to irreversible dynamics as a function of drive amplitude or disk density. At low densities and drives, the system rapidly forms a reversible state where the disks return to their exact positions at the end of each cycle. In contrast, at high amplitudes or high densities, the system enters an irreversible state where the disks exhibit normal diffusion. Between these two regimes, there can be an intermediate irreversible state where most of the system is reversible, but localized irreversible regions are present that are prevented from spreading through the system due to a screening effect from the obstacles. We also find states that we term "combinatorial reversible states" in which the disks return to their original positions after multiple driving cycles. In these states, individual disks exchange positions but form the same configurations during the subcycles of the larger reversible cycle.

Ripples in the bottom of the potential energy landscape of metallic glass

In the absence of periodicity, the structure of glass is ill-defined, and a large number of structural states are found at similar energy levels. However, little is known about how these states are connected to each other in the potential energy landscape. We simulate mechanical relaxation by molecular dynamics for a prototypical [Formula: see text] metallic glass and follow the mechanical energy loss of each atom to track the change in the state. We find that the energy barriers separating these states are remarkably low, only of the order of 1 meV, implying that even quantum fluctuations can overcome these potential energy barriers. Our observation of numerous small ripples in the bottom of the potential energy landscape puts many assumptions regarding the thermodynamic states of metallic glasses into question and suggests that metallic glasses are not totally frozen at the local atomic level.

Reversible-to-irreversible transition of colloidal polycrystals under cyclic athermal quasistatic deformation

Cyclic loading on granular packings and amorphous media exhibits a transition from reversible elastic behavior to irreversible plasticity. The present study compares the irreversibility transition and microscopic details of colloidal polycrystals under oscillatory tensile-compressive and shear strain. Under both modes, the systems exhibit a reversible to irreversible transition. However, the strain amplitude at which the transition is observed is larger in the shear strain than in the tensile-compressive mode. The threshold strain amplitude is confirmed by analyzing the dynamical properties, such as mobility and atomic strain (von Mises shear strain and the volumetric strain). The structural changes are quantified using a hexatic order parameter. Under both modes of deformation, dislocations and grain boundaries in polycrystals disappear, and monocrystals are formed. We also recognize the dislocation motion through grains. The key difference is that strain accumulates diagonally in oscillatory tensile-compressive deformation, whereas in shear deformation, strain accumulation is along the x or y axis.

Granular convergence as an iterated local map

Granular convergence is a property of a granular pack as it is repeatedly sheared in a cyclic, quasistatic fashion, as the packing configuration changes via discrete events. Under suitable conditions the set of microscopic configurations encountered converges to a periodic sequence after sufficient shear cycles. Prior work modeled this evolution as the iteration of a pre-determined, random map from a set of discrete configurations into itself. Iterating such a map from a random starting point leads to similar periodic repetition. This work explores the effect of restricting the randomness of such maps in order to account for the local nature of the discrete events. The number of cycles needed for convergence shows similar statistical behavior to that of numerical granular experiments. The number of cycles in a repeating period behaves only qualitatively like these granular studies.

Counting and Sequential Information Processing in Mechanical Metamaterials

Materials with an irreversible response to cyclic driving exhibit an evolving internal state which, in principle, encodes information on the driving history. Here we realize irreversible metamaterials that count mechanical driving cycles and store the result into easily interpretable internal states. We extend these designs to aperiodic metamaterials that are sensitive to the order of different driving magnitudes, and realize "lock and key" metamaterials that only reach a specific state for a given target driving sequence. Our metamaterials are robust, scalable, and extendable, give insight into the transient memories of complex media, and open new routes towards smart sensing, soft robotics, and mechanical information processing.

Dissipation Indicates Memory Formation in Driven Disordered Systems

Disordered and amorphous materials often retain memories of perturbations they have experienced since preparation. Studying such memories is a gateway to understanding this challenging class of systems. However, it often requires the ability to measure local structural changes in response to external drives. Here, we show that dissipation is a generic macroscopic indicator of the memory of the largest perturbation. Through experiments in crumpled sheets under cyclic drive, we show that dissipation transiently increases when first surpassing the largest perturbation due to irreversible structural changes with unique statistics. This finding is used to devise novel memory readout protocols based on global observables only. The general applicability of this approach is demonstrated by revealing a similar memory effect in a three-dimensional amorphous solid.

Coupled Dynamical Phase Transitions in Driven Disk Packings

Under the influence of oscillatory shear, a monolayer of frictional granular disks exhibits two dynamical phase transitions: a transition from an initially disordered state to an ordered crystalline state and a dynamic active-absorbing phase transition. Although there is no reason a priori for these to be at the same critical point, they are. The transitions may also be characterized by the disk trajectories, which are nontrivial loops breaking time-reversal invariance.

Mechanical annealing and memories in a disordered solid

Shearing a disordered or amorphous solid for many cycles with a constant strain amplitude can anneal it, relaxing a sample to a steady state that encodes a memory of that amplitude. This steady state also features a remarkable stability to amplitude variations that allows one to read the memory. Here, we shed light on both annealing and memory by considering how to mechanically anneal a sample to have as little memory content as possible. In experiments, we show that a "ring-down" protocol reaches a comparable steady state but with no discernible memories and minimal structural anisotropy. We introduce a method to characterize the population of rearrangements within a sample and show how it connects with the response to amplitude variation and the size of annealing steps. These techniques can be generalized to other forms of glassy matter and a wide array of disordered solids, especially those that yield by flowing homogeneously.

Anisotropic correlations of plasticity on the yielding of metallic glasses

We report computer simulations on the shear deformation of CuZr metallic glasses at zero and room temperatures. Shear bands emerge in athermal alloys at strain γ_{c}, with a finite-size effect found. The correlation of nonaffine displacement exhibits an exponential decay even after yielding in thermal alloys, but transits to a power law at γ>γ_{c} in athermal ones. The algebraic exponent is around -1 for the decay inside shear bands, consistent with the theoretical prediction in random elastic media. We quantify the anisotropic correlation with harmonic projection, finding the spectrum is weak in the exponential-decay regime, while it displays a strong polar and quadrupolar symmetry in the power-law regime. The nonvanishing quadrupolar symmetry at long distance signifies the nonlocality of plastic correlation in the athermal alloys. In contrast, the plastic correlation was found to be isotropic and localized at the yielding in the thermal alloys without shear bands.

Inter-particle adhesion induced strong mechanical memory in a dense granular suspension

Repeated/cyclic shearing can drive amorphous solids to a steady state encoding a memory of the applied strain amplitude. However, recent experiments find that the effect of such memory formation on the mechanical properties of the bulk material is rather weak. Here, we study the memory effect in a yield stress solid formed by a dense suspension of cornstarch particles in paraffin oil. Under cyclic shear, the system evolves toward a steady state showing training-induced strain stiffening and plasticity. A readout reveals that the system encodes a strong memory of the training amplitude ( γ T) as indicated by a large change in the differential shear modulus. We observe that memory can be encoded for a wide range of γ T values both above and below the yielding albeit the strength of the memory decreases with increasing γ T. In situ boundary imaging shows strain localization close to the shearing boundaries, while the bulk of the sample moves like a solid plug. In the steady state, the average particle velocity [Formula: see text] inside the solid-like region slows down with respect to the moving plate as γ approaches γ T; however, as the readout strain crosses γ T, [Formula: see text] suddenly increases. We demonstrate that inter-particle adhesive interaction is crucial for such a strong memory effect. Interestingly, our system can also remember more than one input only if the training strain with smaller amplitude is applied last.

Softening and Residual Loss Modulus of Jammed Grains under Oscillatory Shear in an Absorbing State

From a theoretical study of the mechanical response of jammed materials comprising frictionless and overdamped particles under oscillatory shear, we find that the material becomes soft, and the loss modulus remains nonzero even in an absorbing state where any irreversible plastic deformation does not exist. The trajectories of the particles in this region exhibit hysteresis loops. We succeed in clarifying the origin of the softening of the material and the residual loss modulus with the aid of Fourier analysis. We also clarify the roles of the yielding point in the softening to distinguish the plastic deformation from reversible deformation in the absorbing state.

Mean-Field Theory of Yielding under Oscillatory Shear

We study a mean field elastoplastic model, embedded within a disordered landscape of local yield barriers, to shed light on the behavior of athermal amorphous solids subject to oscillatory shear. We show that the model presents a genuine dynamical transition between an elastic and a yielded state, and qualitatively reproduces the dependence on the initial degree of annealing found in particle simulations. For initial conditions prepared below the analytically derived threshold energy, we observe a nontrivial, nonmonotonic approach to the yielded state. The timescale diverges as one approaches the yielding point from above, which we identify with the fatigue limit. We finally discuss the connections to brittle yielding under uniform shear.

Cooperative effects driving the multi-periodic dynamics of cyclically sheared amorphous solids

When subject to cyclic forcing, amorphous solids can reach periodic, repetitive states, where the system behaves plastically, but the particles return to their initial positions after one or more forcing cycles, where the latter response is called multi-periodic. It is known that plasticity in amorphous materials is mediated by local rearrangements called ``soft spots' or ``shear transformation zones'.Experiments and simulations indicate that soft spots can be modeled as hysteretic two-state entities interacting via quadrupolar displacement fields generated when they switch states and that these interactions can give rise to multi-periodic behavior. However, how interactions facilitate multi-periodicity is unknown. Here we show, using a model of random interacting two-state systems and molecular dynamics simulations, that multi-periodicity arises from oscillations in the magnitudes of the switching field of soft spots which cause soft spots to be active during some forcing cycles and idle during others. We demonstrate that these oscillations result from cooperative effects facilitated by the frustrated interactions between the soft spots. The presence of such mechanisms has implications for manipulating memory in frustrated hysteretic systems.

Reversible to Irreversible Transitions for Cyclically Driven Particles on Periodic Obstacle Arrays

We examine the collective dynamics of disks moving through a square array of obstacles under cyclic square wave driving. Below a critical density we find that system organizes into a reversible state in which the disks return to the same positions at the end of every drive cycle. Above this density, the dynamics are irreversible and the disks do not return to the same positions after each cycle. The critical density depends strongly on the angle θ between the driving direction and a symmetry axis of the obstacle array, with the highest critical densities appearing at commensurate angles such as θ=0{degree sign} and θ=45{degree sign} and the lowest critical densities falling at θ=arctan(0.618), the inverse of the golden ratio, where the flow is the most frustrated. As the density increases, the number of cycles required to reach a reversible state grows as a power law with an exponent near ν=1.36, similar to what is found in periodically driven colloidal and superconducting vortex systems.

Turbulent-like velocity fluctuations in two-dimensional granular materials subject to cyclic shear

We perform a systematic experimental study to investigate the velocity fluctuations in the two-dimensional granular matter of low and high friction coefficients subjected to cyclic shear of a range of shear amplitudes, whose velocity fields are strikingly turbulent-like with vortices of different scales. The scaling behaviors of both the transverse velocity power spectra E_T(k) ∝ k^-α_T and, more severely, the longitudinal velocity power spectra E_L(k) ∝ k^-α_L are affected by the prominent peak centered around k ≈ 2π of the inter-particle distance due to the static structure factor of the hard-particle nature in contrast to the real turbulence. To reduce the strong peak effect to the actual values of α_ν (the subscript 'ν' refers to either T or L), we subsequently analyze the second-order velocity structure functions of S(2)ν(r) in real space, which show the power-law scalings of S(2)ν(r) ∝ r^β_ν for both modes. From the values of β_ν, we deduce the corresponding α_ν from the scaling relations of α_ν = β_ν + 2. The deduced values of α_ν increase continuously with the shear amplitude γ_m, showing no signature of yielding transition, and are slightly larger than α_ν = 2.0 at the limit of γ_m → 0, which corresponds to the elastic limit of the system, for all γ_m. The inter-particle friction coefficients show no significant effect on the turbulent-like velocity fluctuations. Our findings suggest that the turbulent-like collective particle motions are governed by both the elasticity and plasticity in cyclically sheared granular materials.

Scaling theory for Wöhler plots in amorphous solids under cyclic forcing

In mechanical engineering Wöhler plots serve to measure the average number of load cycles before materials break, as a function of the maximal stress in each cycle. Although such plots have been prevalent in engineering for more than 150 years, their theoretical understanding is lacking. Recently a scaling theory of Wöhler plots in the context of cyclic bending was offered [Bhowmik, arXiv:2103.03040 (2021)]. Here we elaborate further on cyclic bending and extend the considerations to cyclic tensile loads on an amorphous strip of material; the scaling theory applies to both types of cyclic loading equally well. On the basis of atomistic simulations we conclude that the crucial quantities to focus on are the accumulated damage and the average damage per cycle. The dependence of these quantities on the loading determines the statistics of the number of cycles to failure. Finally, we consider the probability distribution functions of the number of cycles to failure and demonstrate that the scaling theory allows prediction of these distributions at one value of the forcing amplitude from measurements and another value.

Quantifying local rearrangements in three-dimensional granular materials: Rearrangement measures, correlations, and relationship to stresses

Quantifying the ways in which local particle rearrangements contribute to macroscopic plasticity is one of the fundamental pursuits of granular mechanics and soft matter physics. Here we examine local rearrangements that occur naturally during the deformation of three samples of 3D granular materials subjected to distinct boundary conditions by employing in situ x-ray measurements of particle-resolved structure and stress. We focus on five distinct rearrangement measures, their statistics, interrelationships, contributions to macroscopic deformation, repeatability, and dependence on local structure and stress. Our most significant findings are that local rearrangements (1) are correlated on a scale of three to four particle diameters, (2) exhibit volumetric strain-shear strain and nonaffine displacement-rotation coupling, (3) exhibit correlations that suggest either rearrangement repeatability or that rearrangements span multiple steps of incremental sample strain, and (4) show little dependence on local stress but correlate with quantities describing local structure, such as porosity. Our results are presented in the context of relevant plasticity theories and are consistent with recent findings suggesting that local structure may play at least as important of a role as local stress in determining the nature of local rearrangements.

Metastability as a Mechanism for Yielding in Amorphous Solids under Cyclic Shear

We consider the yielding behavior of amorphous solids under cyclic shear deformation and show that it can be mapped into a random walk in a confining potential with an absorbing boundary. The resulting dynamics is governed by the first passage time into the absorbing state and suffices to capture the essential qualitative features recently observed in atomistic simulations of amorphous solids. Our results provide insight into the mechanism underlying yielding and its robustness. When the possibility of activated escape from absorbing states is added, it leads to a unique determination of a threshold energy and yield strain, suggesting thereby an appealing approach to understanding fatigue failure.

Multiple memory formation in glassy landscapes

Cyclically sheared jammed packings form memories of the shear amplitude at which they were trained by falling into periodic orbits where each particle returns to the identical position in subsequent cycles. While simple models that treat clusters of rearranging particles as isolated two-state systems offer insight into this memory formation, they fail to account for the long training times and multiperiod orbits observed in simulated sheared packings. We show that adding interactions between rearranging clusters overcomes these deficiencies. In addition, interactions allow simultaneous encoding of multiple memories, which would not have been possible otherwise. These memories are different in an essential way from those found in other systems, such as multiple transient memories observed in sheared suspensions, and contain information about the strength of the interactions.

Multiperiodic orbits from interacting soft spots in cyclically sheared amorphous solids

When an amorphous solid is deformed cyclically, it may reach a steady state in which the paths of constituent particles trace out closed loops that repeat in each driving cycle. A remarkable variant has been noticed in simulations where the period of particle motions is a multiple of the period of driving, but the reasons for this behavior have remained unclear. Motivated by mesoscopic features of displacement fields in experiments on jammed solids, we propose and analyze a simple model of interacting soft spots-locations where particles rearrange under stress and that resemble two-level systems with hysteresis. We show that multiperiodic behavior can arise among just three or more soft spots that interact with each other, but in all cases it requires frustrated interactions, illuminating this otherwise elusive type of interaction. We suggest directions for seeking this signature of frustration in experiments and for achieving it in designed systems.

Yielding in an amorphous solid subject to constant stress at finite temperatures

Understanding the nature of the yield transition is a long-standing problem in the physics of amorphous solids. Here we use molecular dynamics simulations to study the response of amorphous solids to constant stresses at finite temperatures. We compare amorphous solids that are prepared using fast and slow quenches and show that for thermal systems, the steady-state velocity exhibits a continuous transition from very slow creep to a finite strain rate as a function of the stress. This behavior is observed for both well-annealed and poorly annealed systems. However, the transient dynamics is different in the latter and involves overcoming an energy barrier. Due to the different simulation protocol, the strain rate as a function of stress and temperature follows a scaling relation that is different from the ones that are shown for systems where the strain is controlled. Collapsing the data using this scaling relation allows us to calculate critical exponents for the dynamics close to yield, including an exponent for thermal rounding. We also demonstrate that strain slips due to avalanche events above yield follow standard scaling relations and we extract critical exponents that are comparable to the ones obtained in previous studies that performed simulations of both molecular dynamics and elastoplastic models using strain-rate control.

Topology of the energy landscape of sheared amorphous solids and the irreversibility transition

Recent experiments and simulations of amorphous solids plastically deformed by an oscillatory drive have found a surprising behavior-for small strain amplitudes the dynamics can be reversible, which is contrary to the usual notion of plasticity as an irreversible form of deformation. This reversibility allows the system to reach limit cycles in which plastic events repeat indefinitely under the oscillatory drive. It was also found that reaching reversible limit cycles can take a large number of driving cycles and it was surmised that the plastic events encountered during the transient period are not encountered again and are thus irreversible. Using a graph representation of the stable configurations of the system and the plastic events connecting them, we show that the notion of reversibility in these systems is more subtle. We find that reversible plastic events are abundant and that a large portion of the plastic events encountered during the transient period are actually reversible in the sense that they can be part of a reversible deformation path. More specifically, we observe that the transition graph can be decomposed into clusters of configurations that are connected by reversible transitions. These clusters are the strongly connected components of the transition graph and their sizes turn out to be power-law distributed. The largest of these are grouped in regions of reversibility, which in turn are confined by regions of irreversibility whose number proliferates at larger strains. Our results provide an explanation for the irreversibility transition-the divergence of the transient period at a critical forcing amplitude. The long transients result from transition between clusters of reversibility in a search for a cluster large enough to contain a limit cycle of a specific amplitude. For large enough amplitudes, the search time becomes very large, since the sizes of the limit cycles become incompatible with the sizes of the regions of reversibility.

Yielding in an Integer Automaton Model for Amorphous Solids under Cyclic Shear

We present results on an automaton model of an amorphous solid under cyclic shear. After a transient, the steady state falls into one of three cases in order of increasing strain amplitude: (i) pure elastic behavior with no plastic activity, (ii) limit cycles where the state recurs after an integer period of strain cycles, and (iii) irreversible plasticity with longtime diffusion. The number of cycles N required for the system to reach a periodic orbit diverges as the amplitude approaches the yielding transition between regimes (ii) and (iii) from below, while the effective diffusivity D of the plastic strain field vanishes on approach from above. Both of these divergences can be described by a power law. We further show that the average period T of the limit cycles increases on approach to yielding.

Shear modulus and reversible particle trajectories of frictional granular materials under oscillatory shear

In this study, we numerically investigated the mechanical responses and trajectories of frictional granular particles under oscillatory shear in the reversible phase where particle trajectories form closed loops below the yielding point. When the friction coefficient is small, the storage modulus exhibits softening, and the loss modulus remains finite in the quasi-static limit. As the friction coefficient increases, the softening and residual loss modulus are suppressed. The storage and loss moduli satisfy scaling laws if they are plotted as functions of the areas of the loop trajectories divided by the strain amplitude and diameter of grains, at least for small values of the areas.

Crystalline shielding mitigates structural rearrangement and localizes memory in jammed systems under oscillatory shear

The nature of yield in amorphous materials under stress has yet to be fully elucidated. In particular, understanding how microscopic rearrangement gives rise to macroscopic structural and rheological signatures in disordered systems is vital for the prediction and characterization of yield and the study of how memory is stored in disordered materials. Here, we investigate the evolution of local structural homogeneity on an individual particle level in amorphous jammed two-dimensional (athermal) systems under oscillatory shear and relate this evolution to rearrangement, memory, and macroscale rheological measurements. We define the structural metric crystalline shielding, and show that it is predictive of rearrangement propensity and structural volatility of individual particles under shear. We use this metric to identify localized regions of the system in which the material's memory of its preparation is preserved. Our results contribute to a growing understanding of how local structure relates to dynamic response and memory in disordered systems.

Directional clogging and phase separation for disk flow through periodic and diluted obstacle arrays

We model collective disk flow though a square array of obstacles as the flow direction is changed relative to the symmetry directions of the array. At lower disk densities there is no clogging for any driving direction, but as the disk density increases, the average disk velocity decreases and develops a drive angle dependence. For certain driving angles, the flow is reduced or drops to zero when the system forms a heterogeneous clogged state consisting of high density clogged regions coexisting with empty regions. The clogged states are fragile and can be unclogged by changing the driving angle. For large obstacle sizes, we find a uniform clogged state that is distinct from the collective clogging regime. Within the clogged phases, depinning transitions can occur as a function of increasing driving force, with intermittent motion appearing just above the depinning threshold. The clogging is robust against the random removal or dilution of the obstacle sites, and the disks are able to form system-spanning clogged clusters even under increasing dilution. If the dilution becomes too large, however, the clogging behavior is lost.

Period multiplication cascade at the order-by-disorder transition in uniaxial random field XY magnets

Uniaxial random field disorder induces a spontaneous transverse magnetization in the XY model. Adding a rotating driving field, we find a critical point attached to the number of driving cycles needed to complete a limit cycle, the first discovery of this phenomenon in a magnetic system. Near the critical drive, time crystal behavior emerges, in which the period of the limit cycles becomes an integer n > 1 multiple of the driving period. The period n can be engineered via specific disorder patterns. Because n generically increases with system size, the resulting period multiplication cascade is reminiscent of that occurring in amorphous solids subject to oscillatory shear near the onset of plastic deformation, and of the period bifurcation cascade near the onset of chaos in nonlinear systems, suggesting it is part of a larger class of phenomena in transitions of dynamical systems. Applications include magnets, electron nematics, and quantum gases.

Jamming, fragility and pinning phenomena in superconducting vortex systems

We examine driven superconducting vortices interacting with quenched disorder under a sequence of perpendicular drive pulses. As a function of disorder strength, we find four types of behavior distinguished by the presence or absence of memory effects. The fragile and jammed states exhibit memory, while the elastic and pinning dominated regimes do not. In the fragile regime, the system organizes into a pinned state during the first pulse, flows during the second perpendicular pulse, and then returns to a pinned state during the third pulse which is parallel to the first pulse. This behavior is the hallmark of the fragility proposed for jamming in particulate matter. For stronger disorder, we observe a robust jamming state with memory where the system reaches a pinned or reduced flow state during the perpendicular drive pulse, similar to the shear jamming of granular systems. We show signatures of the different states in the spatial vortex configurations, and find that memory effects arise from coexisting elastic and pinned components of the vortex assembly. The sequential perpendicular driving protocol we propose for distinguishing fragile, jammed, and pinned phases should be general to the broader class of driven interacting particles in the presence of quenched disorder.

State transition graph of the Preisach model and the role of return-point memory

The Preisach model has been useful as a null model for understanding memory formation in periodically driven disordered systems. In amorphous solids, for example, the athermal response to shear is due to localized plastic events (soft spots). As shown recently by Mungan et al. [Phys. Rev. Lett. 123, 178002 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.178002], the plastic response to applied shear can be rigorously described in terms of a directed network whose transitions correspond to one or more soft spots changing states. The topology of this graph depends on the interactions between soft spots and when such interactions are negligible, the resulting description becomes that of the Preisach model. A first step in linking transition graph topology with the underlying soft-spot interactions is therefore to determine the structure of such graphs in the absence of interactions. Here we perform a detailed analysis of the transition graph of the Preisach model. We highlight the important role played by return-point memory in organizing the graph into a hierarchy of loops and subloops. Our analysis reveals that the topology of a large portion of this graph is actually not governed by the values of the switching fields that describe the hysteretic behavior of the individual elements but by a coarser parameter, a permutation ρ which prescribes the sequence in which the individual hysteretic elements change their states as the main hysteresis loop is traversed. This in turn allows us to derive combinatorial properties, such as the number of major loops in the transition graph as well as the number of states |R| constituting the main hysteresis loop and its nested subloops. We find that |R| is equal to the number of increasing subsequences contained in the permutation ρ.

Glass Stability Changes the Nature of Yielding under Oscillatory Shear

We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our results disentangle the relative role of mechanical and thermal annealing on the mechanical properties of amorphous solids, which is relevant for various experimental situations from the rheology of soft materials to fatigue failure in metallic glasses.

Quantification of plasticity via particle dynamics above and below yield in a 2D jammed suspension

The yield transition of amorphous materials is characterized by a swift increase of energy dissipation. The connection between particle dynamics, dissipation, and overall material rheology, however, has still not been elucidated. Here, we take a new approach relating trajectories to yielding, using a custom built interfacial stress rheometer, which allows for measurement of shear moduli (G',G'') of a dense athermal suspension's microstructure while simultaneously tracking particle trajectories undergoing cyclic shear. We find an increase in total area traced by particle trajectories as the system is stressed well below to well above yield. Trajectories may be placed into three categories: reversibly elastic paths; reversibly plastic paths, associated with smooth limit cycles; and irreversibly plastic paths, in which particles do not return to their original position. We find that above yield, reversibly plastic trajectories are predominantly found near to the shearing surface, whereas reversibly elastic paths are more prominent near the stationary wall. This spatial transition between particles acting as liquids to those acting as solids is characteristic of a 'melting front', which is observed to shift closer to the wall with increasing strain. We introduce a non-dimensional measure of plastic dissipation based on particle trajectories that scales linearly with strain amplitude both above and below yield, and that is unity at the rheological yield point. Surprisingly, this relation collapses for three systems of varying degrees of disorder.

Unified phase diagram of reversible-irreversible, jamming, and yielding transitions in cyclically sheared soft-sphere packings

Self-organization, and transitions from reversible to irreversible behavior, of interacting particle assemblies driven by externally imposed stresses or deformation is of interest in comprehending diverse phenomena in soft matter. They have been investigated in a wide range of systems, such as colloidal suspensions, glasses, and granular matter. In different density and driving regimes, such behavior is related to yielding of amorphous solids, jamming, memory formation, etc. How these phenomena are related to each other has not, however, been much studied. In order to obtain a unified view of the different regimes of behavior, and transitions between them, we investigate computationally the response of soft-sphere assemblies to athermal cyclic-shear deformation over a wide range of densities and amplitudes of shear deformation. Cyclic-shear deformation induces transitions from reversible to irreversible behavior in both unjammed and jammed soft-sphere packings. Well above the minimum isotropic jamming density ([Formula: see text]), this transition corresponds to yielding. In the vicinity of the jamming point, up to a higher-density limit, we designate [Formula: see text], an unjammed phase emerges between a localized, absorbing phase and a diffusive, irreversible, phase. The emergence of the unjammed phase signals the shifting of the jamming point to higher densities as a result of annealing and opens a window where shear jamming becomes possible for frictionless packings. Below [Formula: see text], two distinct localized states, termed point- and loop-reversible, are observed. We characterize in detail the different regimes and transitions between them and obtain a unified density-shear amplitude phase diagram.

Forced deterministic dynamics on a random energy landscape: Implications for the physics of amorphous solids

The dynamics of supercooled liquids and plastically deformed amorphous solids is known to be dominated by the structure of their rough energy landscapes. Recent experiments and simulations on amorphous solids subjected to oscillatory shear at athermal conditions have shown that for small strain amplitudes these systems reach limit cycles of different periodicities after a transient. However, for larger strain amplitudes the transients become longer and for strain amplitudes exceeding a critical value the system reaches a diffusive steady state. This behavior cannot be explained using the current mean-field models of amorphous plasticity. Here we show that this phenomenology can be described and explained using a simple model of forced dynamics on a multidimensional random energy landscape. In this model, the existence of limit cycles can be ascribed to confinement of the dynamics to a small part of the energy landscape which leads to self-intersection of state-space trajectories and the transition to the diffusive regime for larger forcing amplitudes occurs when the forcing overcomes this confinement.

Rheological similarities between dense self-propelled and sheared particulate systems

Different from previous modeling of self-propelled particles, we develop a method to propel particles with a constant average velocity instead of a constant force. This constant propulsion velocity (CPV) approach is validated by its agreement with the conventional constant propulsion force (CPF) approach in the flowing regime. However, the CPV approach shows its advantage of accessing quasistatic flows of yield stress fluids with a vanishing propulsion velocity, while the CPF approach is usually unable to because of finite system size. Taking this advantage, we realize cyclic self-propulsion and study the evolution of the propulsion force with the propelled particle displacement, both in the quasistatic flow regime. By mapping the shear stress and shear rate to the propulsion force and propulsion velocity, we find similar rheological behaviors of self-propelled systems to sheared systems, including the yield force gap between the CPF and CPV approaches, propulsion force overshoot, reversible-irreversible transition under cyclic propulsion, and propulsion bands in plastic flows. These similarities suggest underlying connections between self-propulsion and shear, although they act on systems in different ways.

Absorbing-State Transitions in Granular Materials Close to Jamming

We consider a model for driven particulate matter in which absorbing states can be reached both by particle isolation and by particle caging. The model predicts a nonequilibrium phase diagram in which analogs of hydrodynamic and elastic reversibility emerge at low and high volume fractions respectively, partially separated by a diffusive, nonabsorbing region. We thus find a single phase boundary that spans the onset of chaos in sheared suspensions to the onset of yielding in jammed packings. This boundary has the properties of a nonequilibrium second order phase transition, leading us to write a Manna-like mean field description that captures the model predictions. Dependent on contact details, jamming marks either a direct transition between the two absorbing states, or occurs within the diffusive region.

Shear Band Formation in Amorphous Materials under Oscillatory Shear Deformation

The effect of periodic shear on strain localization in disordered solids is investigated using molecular dynamics simulations. We consider a binary mixture of one million atoms annealed to a low temperature with different cooling rates and then subjected to oscillatory shear deformation with a strain amplitude slightly above the critical value. It is found that the yielding transition occurs during one cycle but the accumulation of irreversible displacements and initiation of the shear band proceed over larger number of cycles for more slowly annealed glasses. The spatial distribution and correlation function of nonaffine displacements reveal that their collective dynamics changes from homogeneously distributed small clusters to a system-spanning shear band. The analysis of spatially averaged profiles of nonaffine displacements indicates that the location of a shear band in periodically loaded glasses can be identified at least several cycles before yielding. These insights are important for the development of novel processing methods and prediction of the fatigue lifetime of metallic glasses.

Annealing and rejuvenation in a two-dimensional model amorphous solid under oscillatory shear

We study the annealing and rejuvenation behavior of a two-dimensional amorphous solid model under oscillatory shear. We show that, depending on the cooling protocol used to create the initial configuration, the mean potential energy can either decrease or increase under subyield oscillatory shear. For post-yield oscillatory shear, the mean potential energy increases and is independent on the initial conditions. We explain this behavior by modeling the dynamics using a simple model of forced dynamics on a random energy landscape and show that the model reproduces the qualitative behavior of the mean potential energy and mean-square displacement observed in the particle based simulations. This suggests that some important aspects of the dynamics of amorphous solids can be understood by studying the properties of random energy landscapes and without explicitly taking into account the complex real-space interactions which are involved in plastic deformation.

Experimental signatures of a nonequilibrium phase transition near the crossover point of a Langmuir monolayer

We investigate the response of the two-dimensional (2D) continuous non-particulate film of surfactant sorbitan tristearate confined at the air-water interface under oscillatory shear deformation. The time dependence of various rheological parameters show critical-like behavior at a value of strain amplitude close to the crossover point of elastic ([Formula: see text]) and viscous ([Formula: see text]) shear moduli. Imposing oscillatory shear of different strain amplitudes ([Formula: see text]) above and below the crossover strain amplitude ([Formula: see text]) over a large number of cycles, we quantify the temporal dependence of interfacial viscous modulus, phase angle ([Formula: see text]) as well as higher harmonic components of stress. The number of shear cycles ([Formula: see text]) required for these quantities to reach the steady state value diverges near [Formula: see text]. The steady state values of the third harmonic ([Formula: see text]) show order parameter like behavior indicating the importance of higher order harmonics near the nonequilibrium transition. We further show that the energy dissipation per cycle per unit volume has a marked change near [Formula: see text], consistent with continuum level nonequilibrium shear-transformation-zone model of amorphous viscoplasticity.

Critical behavior near the reversible-irreversible transition in periodically driven vortices under random local shear

When many-particle (vortex) assemblies with disordered distribution are subjected to a periodic shear with a small amplitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}$$\end{document}d, the particles gradually self-organize to avoid next collisions and transform into an organized configuration. We can detect it from the time-dependent voltage \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{V}}{\boldsymbol{(}}{\boldsymbol{t}}{\boldsymbol{)}}$$\end{document}V(t) (average velocity) that increases towards a steady-state value. For small \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}$$\end{document}d, the particles settle into a reversible state where all the particles return to their initial position after each shear cycle, while they reach an irreversible state for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}$$\end{document}d above a threshold \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{d}}}_{{\boldsymbol{c}}}$$\end{document}dc. Here, we investigate the general phenomenon of a reversible-irreversible transition (RIT) using periodically driven vortices in a strip-shaped amorphous film with random pinning that causes local shear, as a function of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}$$\end{document}d. By measuring \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{V}}{\boldsymbol{(}}{\boldsymbol{t}}{\boldsymbol{)}}$$\end{document}V(t), we observe a critical behavior of RIT, not only on the irreversible side, but also on the reversible side of the transition, which is the first under random local shear. The relaxation time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\tau }}{\boldsymbol{(}}{\boldsymbol{d}}{\boldsymbol{)}}$$\end{document}τ(d) to reach either the reversible or irreversible state shows a power-law divergence at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{d}}}_{{\boldsymbol{c}}}$$\end{document}dc. The critical exponent is determined with higher accuracy and is, within errors, in agreement with the value expected for an absorbing phase transition in the two-dimensional directed-percolation universality class. As \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{d}}$$\end{document}d is decreased down to the intervortex spacing in the reversible regime, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\tau }}{\boldsymbol{(}}{\boldsymbol{d}}{\boldsymbol{)}}$$\end{document}τ(d) deviates downward from the power-law relation, reflecting the suppression of intervortex collisions. We also suggest the possibility of a narrow smectic-flow regime, which is predicted to intervene between fully reversible and irreversible flow.

Networks and Hierarchies: How Amorphous Materials Learn to Remember

We consider the slow and athermal deformations of amorphous solids and show how the ensuing sequence of discrete plastic rearrangements can be mapped onto a directed network. The network topology reveals a set of highly connected regions joined by occasional one-way transitions. The highly connected regions include hierarchically organized hysteresis cycles and subcycles. At small to moderate strains this organization leads to near-perfect return point memory. The transitions in the network can be traced back to localized particle rearrangements (soft spots) that interact via Eshelby-type deformation fields. By linking topology to dynamics, the network representations provide new insight into the mechanisms that lead to reversible and irreversible behavior in amorphous solids.

Classification of the reversible-irreversible transitions in particle trajectories across the jamming transition point

The reversible-irreversible (RI) transition of particle trajectories in athermal colloidal suspensions under cyclic shear deformation is an archetypal nonequilibrium phase transition which has attracted much attention recently. Most studies of the RI transitions have focused on either dilute limit or very high densities well above the jamming transition point. The transition between the two limiting cases is largely unexplored. In this paper, we study the RI transition of athermal frictionless colloidal particles over a wide range of densities, with emphasis on the region below φ_J, by using oscillatory sheared molecular dynamics simulation. We reveal that the nature of the RI transitions in the intermediate densities is very rich. As demonstrated by the previous work by Schreck et al. [Phys. Rev. E: Stat., Nonlinear, Soft Matter Phys., 2013, 88, 052205], there exist the point-reversible and the loop-reversible phases depending on the density and the shear strain amplitude. We find that, between the two reversible phases, a quasi-irreversible phase where the particles' trajectories are highly non-affine and diffusive. The averaged number of contacts of particles is found to characterize the phase boundaries. We also find that the system undergoes the yielding transition below but in the vicinity of φ_J when the strain with a small but finite strain rate is applied. This yielding transition line matches with the RI transition line separating the loop-reversible from the irreversible phases. Surprisingly, the nonlinear rheological response called "softening" has been observed even below φ_J. These findings imply that geometrical properties encoded in the sheared configurations control the dynamical transitions.

Cyclic annealing as an iterated random map

Disordered magnets, martensitic mixed crystals, and glassy solids can be irreversibly deformed by subjecting them to external deformation. The deformation produces a smooth, reversible response punctuated by abrupt relaxation "glitches." Under appropriate repeated forward and reverse deformation producing multiple glitches, a strict repetition of a single sequence of microscopic configurations often emerges. We exhibit these features by describing the evolution of the system configuration from glitch to glitch as a mapping of N states into one another. A map U controls forward deformation; a second map D controls reverse deformation. Iteration of a given sequence of forward and reverse maps, e.g., DDDDUUUU necessarily produces a convergence to a fixed cyclic repetition of states covering multiple glitches. The repetition may have a period of more than one strain cycle, as recently observed in simulations. Using numerical sampling, we characterize the convergence properties of four types of random maps implementing successive physical restrictions. The most restrictive is the much-studied Preisach model. These maps show only the most qualitative resemblance to annealing simulations. However, they suggest further properties needed for a realistic mapping scheme.

Strength of Mechanical Memories is Maximal at the Yield Point of a Soft Glass

We show experimentally that both single and multiple mechanical memories can be encoded in an amorphous bubble raft, a prototypical soft glass, subject to an oscillatory strain. In line with recent numerical results, we find that multiple memories can be formed sans external noise. By systematically investigating memory formation for a range of training strain amplitudes spanning yield, we find clear signatures of memory even beyond yielding. Most strikingly, the extent to which the system recollects memory is largest for training amplitudes near the yield strain and is a direct consequence of the spatial extent over which the system reorganizes during the encoding process. Our study further suggests that the evolution of force networks on training plays a decisive role in memory formation in jammed packings.

Critical diffusivity in the reversibility-irreversibility transition of amorphous solids under oscillatory shear

Recently it was shown that under oscillatory shear at zero temperature an amorphous solid transitions from asymptotically periodic to asymptotically diffusive steady-state at a critical maximal strain amplitude. Current understanding of the physics behind this transition is lacking. Here we show, using computer simulations, evidence that the diffusivity of the vector of coordinates of the particles comprising an amorphous solid, when subject to oscillatory shear, undergoes a second order phase transition at the reversibility-irreversibility transition point. We explain how such a transition is consistent with dissipative forced dynamics on a complex energy landscape, such as is known to exist in amorphous solids. We demonstrate that as the forcing increases, more and more state-space volume becomes accessible to the system, making it less probable for the state-space trajectory of the system to self-intersect and form a limit-cycle, which explains the slowing-down observed at the transition.