Chaos and threshold for irreversibility in sheared suspensions

Non-equilibrium dynamic hyperuniform states

Disordered hyperuniform structures are an exotic state of matter having suppressed density fluctuations at large length-scale similar to perfect crystals and quasicrystals but without any long range orientational order. In the past decade, an increasing number of non-equilibrium systems were found to have dynamic hyperuniform states, which have emerged as a new research direction coupling both non-equilibrium physics and hyperuniformity. Here we review the recent progress in understanding dynamic hyperuniform states found in various non-equilibrium systems, including the critical hyperuniformity in absorbing phase transitions, non-equilibrium hyperuniform fluids and the hyperuniform structures in phase separating systems via spinodal decomposition.

Training physical matter to matter

Biological systems offer a great many examples of how sophisticated, highly adapted behavior can emerge from training. Here we discuss how training might be used to impart similarly adaptive properties in physical matter. As a special form of materials processing, training differs in important ways from standard approaches of obtaining sought after material properties. In particular, rather than designing or programming the local configurations and interactions of constituents, training uses externally applied stimuli to evolve material properties. This makes it possible to obtain different functionalities from the same starting material (pluripotency). Furthermore, training evolves a material in situ or under conditions similar to those during the intended use; thus, material performance can improve rather than degrade over time. We discuss requirements for trainability, outline recently developed training strategies for creating soft materials with multiple, targeted and adaptable functionalities, and provide examples where the concept of training has been applied to materials on length scales from the molecular to the macroscopic.

Jamming is a first-order transition with quenched disorder in amorphous materials sheared by cyclic quasistatic deformations

Jamming is an athermal transition between flowing and rigid states in amorphous systems such as granular matter, colloidal suspensions, complex fluids and cells. The jamming transition seems to display mixed aspects of a first-order transition, evidenced by a discontinuity in the coordination number, and a second-order transition, indicated by power-law scalings and diverging lengths. Here we demonstrate that jamming is a first-order transition with quenched disorder in cyclically sheared systems with quasistatic deformations, in two and three dimensions. Based on scaling analyses, we show that fluctuations of the jamming density in finite-sized systems have important consequences on the finite-size effects of various quantities, resulting in a square relationship between disconnected and connected susceptibilities, a key signature of the first-order transition with quenched disorder. This study puts the jamming transition into the category of a broad class of transitions in disordered systems where sample-to-sample fluctuations dominate over thermal fluctuations, suggesting that the nature and behavior of the jamming transition might be better understood within the developed theoretical framework of the athermally driven random-field Ising model.

Enhancing (quasi-)long-range order in a two-dimensional driven crystal

It has been recently shown that 2D systems can exhibit crystalline phases with long-range translational order showcasing a striking violation of the Hohenberg-Mermin-Wagner (HMW) theorem, which is valid at equilibrium. This is made possible by athermal driving mechanisms that inject energy into the system without exciting long wavelength modes of the density field, thereby inducing hyperuniformity. However, as thermal fluctuations are superimposed on the non-equilibrium driving, long-range translational order is inevitably lost. Here, we discuss the possibility of exploiting non-equilibrium effects to suppress arbitrarily large density fluctuations even when a global thermal bath is coupled to the system. We introduce a model of a harmonic crystal driven both by a global thermal bath and by a momentum conserving noise, where the typical observables related to density fluctuations and long-range translational order can be analytically derived and put in relation. This model allows us to rationalize the violation of the HMW theorem observed in previous studies through the prediction of large-wavelength phonons, which thermalize at a vanishing effective temperature when the global bath is switched off. The conceptual framework introduced through this theory is then applied to numerical simulations of a hard-disk solid in contact with a thermal bath and driven out-of-equilibrium by active collisions. Our numerical analysis demonstrates how varying driving and dissipative parameters can lead to an arbitrary enhancement of the quasi-long-range order in the system regardless of the applied global noise amplitude. Finally, we outline a possible experimental procedure to apply our results to a realistic granular system.

Anomalous relaxation and hyperuniform fluctuations in center-of-mass conserving systems with broken time-reversal symmetry

We study the Oslo model, a paradigm for absorbing-phase transition, on a one-dimensional ring of L sites with a fixed global density ρ[over ¯]; we consider the system strictly above critical density ρ_{c}. Notably, microscopic dynamics conserve both mass and center of mass (CoM), but lack time-reversal symmetry. We show that, despite having highly constrained dynamics due to CoM conservation, the system exhibits diffusive relaxation away from criticality and superdiffusive relaxation near criticality. Furthermore, the CoM conservation severely restricts particle movement, causing the mobility-a transport coefficient analogous to the conductivity for charged particles-to vanish exactly. Indeed, the steady-state temporal growth of current fluctuation is qualitatively different from that observed in diffusive systems with a single conservation law. Remarkably, far from criticality where the relative density Δ=ρ[over ¯]-ρ_{c}≫ρ_{c}, the second cumulant, or the variance, 〈Q_{i}^{2}(T,Δ)〉_{c}, of current Q_{i} across the ith bond up to time T in the steady-state saturates as 〈Q_{i}^{2}〉_{c}≃Σ_{Q}^{2}(Δ)-constT^{-1/2}; near criticality, it grows subdiffusively as 〈Q_{i}^{2}〉_{c}∼T^{α}, with 0<α<1/2, and eventually saturates to Σ_{Q}^{2}(Δ). Interestingly, the asymptotic current fluctuation Σ_{Q}^{2}(Δ) is a nonmonotonic function of Δ: It diverges as Σ_{Q}^{2}(Δ)∼Δ^{2} for Δ≫ρ_{c} and Σ_{Q}^{2}(Δ)∼Δ^{-δ}, with δ>0, for Δ→0^{+}. Using a mass-conservation principle, we exactly determine the exponents δ=2(1-1/ν_{⊥})/ν_{⊥} and α=δ/zν_{⊥} via the correlation-length and dynamic exponents, ν_{⊥} and z, respectively. Finally, we show that in the steady state the self-diffusion coefficient D_{s}(ρ[over ¯]) of tagged particles is connected to activity through the relation D_{s}(ρ[over ¯])=a(ρ[over ¯])/ρ[over ¯].

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.

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.

Dynamical Approach to the Jamming Problem

A simple dynamical model, biased random organization (BRO), appears to produce configurations known as random close packing (RCP) as BRO's densest critical point in dimension d=3. We conjecture that BRO likewise produces RCP in any dimension; if so, then RCP does not exist in d=1-2 (where BRO dynamics lead to crystalline order). In d=3-5, BRO produces isostatic configurations and previously estimated RCP volume fractions 0.64, 0.46, and 0.30, respectively. For all investigated dimensions (d=2-5), we find that BRO belongs to the Manna universality class of dynamical phase transitions by measuring critical exponents associated with the steady-state activity and the long-range density fluctuations. Additionally, BRO's distribution of near contacts (gaps) displays behavior consistent with the infinite-dimensional theoretical treatment of RCP when d≥4. The association of BRO's densest critical configurations with random close packing implies that RCP's upper-critical dimension is consistent with the Manna class d_{uc}=4.

Continuum contact process and influence of impurity on the critical behavior in absorbing-state phase transitions in two dimensions

We study via Monte Carlo simulations the influence of quenched and mobile impurities in the contact process (CP) on two-dimensional lattice and continuum systems. In the lattice system, the effect of mobile impurity was studied for the density n_{i}=0.2 and two selected values of hopping probability for impurity particles, w=0.5 and 1. In the continuum system, the CP was defined by distributing spherical impurity particles of diameter σ_{i} and number density n_{i}=0.2 and active particles of diameter unity and number density 1-n_{i} on a square substrate with periodic boundaries. In each dynamic process, a particle is selected at random; the active particle either creates with a rate λ an offspring at a distance r (1≤r≤1.5) from the active particle or annihilates with a unit rate, and the impurity particle hops a distance r (0≤r≤1), both along randomly selected directions. We found that the lattice CP shows power-law behaviors with varying critical exponents depending on the values of w. For the continuum CP with quenched impurity, the critical behavior followed the activated scaling scenario, whereas with mobile impurity usual power-law behaviors were observed but the critical exponents varied depending on the values of σ_{i}.

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.

Two-Dimensional Crystals far from Equilibrium

When driven by nonequilibrium fluctuations, particle systems may display phase transitions and physical behavior with no equilibrium counterpart. We study a two-dimensional particle model initially proposed to describe driven non-Brownian suspensions undergoing nonequilibrium absorbing phase transitions. We show that when the transition occurs at large density, the dynamics produces long-range crystalline order. In the ordered phase, long-range translational order is observed because equipartition of energy is lacking, phonons are suppressed, and density fluctuations are hyperuniform. Our study offers an explicit microscopic model where nonequilibrium violations of the Mermin-Wagner theorem stabilize crystalline order in two dimensions.

Equilibrium states corresponding to targeted hyperuniform nonequilibrium pair statistics

The Zhang-Torquato conjecture [G. Zhang and S. Torquato, Phys. Rev. E, 2020, 101, 032124.] states that any realizable pair correlation function g₂(r) or structure factor S(k) of a translationally invariant nonequilibrium system can be attained by an equilibrium ensemble involving only (up to) effective two-body interactions. To further test and study this conjecture, we consider two singular nonequilibrium models of recent interest that also have the exotic hyperuniformity property: a 2D "perfect glass" and a 3D critical absorbing-state model. We find that each nonequilibrium target can be achieved accurately by equilibrium states with effective one- and two-body potentials, lending further support to the conjecture. To characterize the structural degeneracy of such a nonequilibrium-equilibrium correspondence, we compute higher-order statistics for both models, as well as those for a hyperuniform 3D uniformly randomized lattice (URL), whose higher-order statistics can be very precisely ascertained. Interestingly, we find that the differences in the higher-order statistics between nonequilibrium and equilibrium systems with matching pair statistics, as measured by the "hole" probability distribution, provide measures of the degree to which a system is out of equilibrium. We show that all three systems studied possess the bounded-hole property and that holes near the maximum hole size in the URL are much rarer than those in the underlying simple cubic lattice. Remarkably, upon quenching, the effective potentials for all three systems possess local energy minima (i.e., inherent structures) with stronger forms of hyperuniformity compared to their target counterparts. Our methods are expected to facilitate the self-assembly of tunable hyperuniform soft-matter systems.

Rheology of periodically sheared suspensions undergoing reversible-irreversible transition

The rheology of noncolloidal suspensions under cyclic shear is studied numerically. The main findings are a strain amplitude (γ_{0}) dependent response in the shear stress and second normal stress difference (N_{2}). Specifically, we find a reduced viscosity, an enhanced intracycle shear thinning, the onset of a finite N_{2}, and its frequency doubling, all near a critical strain amplitude γ_{c} that scales with the volume fraction ϕ as γ_{c}∼ϕ^{-2}. These rheological changes also signify a reversible-irreversible transition (RIT), dividing stroboscopic particle dynamics into a reversible absorbing phase (for γ_{0}<γ_{c}) and a persistently diffusing phase (for γ_{0}>γ_{c}). We explain the results based on two flow-induced mechanisms and elucidate their connection in the context of RIT through the underlying microstructure, which tends toward hyperuniformity near γ_{0}=γ_{c}. Overall, we expect this correspondence between rheology and emergent dynamics to hold in a wide range of settings where structural organizations are dominated by volume exclusions.

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.

Microstructural changes of concentrated Newtonian suspensions in the first oscillation cycles probed with linear and non-linear rheology

Carotenuto et al. (Rheol Acta, 2021, 60, 309) recently showed that the complex viscosity of a Newtonian non-Brownian suspension is smaller than the steady shear one, whatever the imposed strain amplitude. Oscillatory shear can alter the microstructure through a shear induced particle diffusion mechanism. This mechanism needs time to show its effect and cannot be invoked to explain the observed mismatch between the steady shear and the complex viscosity. Moreover, in the limit of vanishing strain amplitudes and of very large ones, where the oscillatory shear is equivalent to consecutive steady flow reversals, the oscillatory shear should not alter the microstructure and the Cox-Merz rule should hold. With a combination of approaches exploiting the Lissajous-Bowditch plots, the Fourier transform rheology and the Sequence of Physical Processes, we investigate the microstructure changes induced in the first oscillatory cycles. The results from the different analyses agree with the microstructure rearranging mechanisms proposed by Carotenuto et al.: at small amplitudes, the oscillatory shear rotates couples of touching particles towards the flow direction, at medium amplitudes it breaks particle clusters and at very large amplitudes it reshuffles and reorients all the particles. We show that the vast majority of the microstructure rearrangement occurs soon after the flow inversion of the first cycle, while before it the microstructure is not altered. This allows us to suggest a procedure to "recover" the Cox-Merz rule: a single cycle of oscillation must be imposed and the stress response of the sole first quarter of oscillation must be analysed.

Evidence for Chaotic Behavior during the Yielding of a Soft Particle Glass

Materials comprised of deformable particles such as microgels and concentrated emulsions and foams display complex rheological behavior that includes a yielding transition from an elastic solid to viscous fluid. Most studies of this class of soft matter involve shear flows, and only a handful report both shear and normal stresses. We present measurements of the shear stress and two normal stress differences for a microgel subjected to constant shear rate flows. The shear stress evolves through the yield point in a manner indicative of simple yield stress fluid behavior. Prior to yielding, the normal stress differences are immeasurable; beyond the yield point, they evolve in a reproducibly chaotic manner.

Chaos and irreversibility of a flexible filament in periodically driven Stokes flow

The flow of Newtonian fluid at low Reynolds number is, in general, regular and time-reversible due to absence of nonlinear effects. For example, if the fluid is sheared by its boundary motion that is subsequently reversed, then all the fluid elements return to their initial positions. Consequently, mixing in microchannels happens solely due to molecular diffusion and is very slow. Here, we show, numerically, that the introduction of a single, freely floating, flexible filament in a time-periodic linear shear flow can break reversibility and give rise to chaos due to elastic nonlinearities, if the bending rigidity of the filament is within a carefully chosen range. Within this range, not only the shape of the filament is spatiotemporally chaotic, but also the flow is an efficient mixer. Overall, we find five dynamical phases: the shape of a stiff filament is time-invariant-either straight or buckled; it undergoes a period-two bifurcation as the filament is made softer; becomes spatiotemporally chaotic for even softer filaments but, surprisingly, the chaos is suppressed if bending rigidity is decreased further.

Model architecture can transform catastrophic forgetting into positive transfer

The work of McCloskey and Cohen popularized the concept of catastrophic interference. They used a neural network that tried to learn addition using two groups of examples as two different tasks. In their case, learning the second task rapidly deteriorated the acquired knowledge about the previous one. We hypothesize that this could be a symptom of a fundamental problem: addition is an algorithmic task that should not be learned through pattern recognition. Therefore, other model architectures better suited for this task would avoid catastrophic forgetting. We use a neural network with a different architecture that can be trained to recover the correct algorithm for the addition of binary numbers. This neural network includes conditional clauses that are naturally treated within the back-propagation algorithm. We test it in the setting proposed by McCloskey and Cohen and training on random additions one by one. The neural network not only does not suffer from catastrophic forgetting but it improves its predictive power on unseen pairs of numbers as training progresses. We also show that this is a robust effect, also present when averaging many simulations. This work emphasizes the importance that neural network architecture has for the emergence of catastrophic forgetting and introduces a neural network that is able to learn an algorithm.

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.

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.

Absorbing phase transitions in systems with mediated interactions

Experiments of periodically sheared colloidal suspensions or soft amorphous solids display a transition from reversible to irreversible particle motion that, when analyzed stroboscopically in time, is interpreted as an absorbing phase transition with infinitely many absorbing states. In these systems, interactions mediated by hydrodynamics or elasticity are present, causing passive regions to be affected by nearby active ones. We show that mediated interactions induce a universality class of absorbing phase transitions distinct from conserved directed percolation, and we obtain the corresponding critical exponents. We do so with large-scale numerical simulations of a minimal model for the stroboscopic dynamics of sheared soft materials and we derive the minimal field theoretical description.

Critical behavior of nonequilibrium depinning transitions for vortices driven by current and vortex density

We study the critical dynamics of vortices associated with dynamic disordering near the depinning transitions driven by dc force (dc current I) and vortex density (magnetic field B). Independent of the driving parameters, I and B, we observe the critical behavior of the depinning transitions, not only on the moving side, but also on the pinned side of the transition, which is the first convincing verification of the theoretical prediction. Relaxation times, [Formula: see text] and [Formula: see text], to reach either the moving or pinned state, plotted against I and B, respectively, exhibit a power-law divergence at the depinning thresholds. The critical exponents of both transitions are, within errors, identical to each other, which are in agreement with the values expected for an absorbing phase transition in the two-dimensional directed-percolation universality class. With an increase in B under constant I, the depinning transition at low B is replaced by the repinning transition at high B in the peak-effect regime. We find a trend that the critical exponents in the peak-effect regime are slightly smaller than those in the low-B regime and the theoretical one, which is attributed to the slight difference in the depinning mechanism in the peak-effect regime.

Shear-induced migration of confined flexible fibers

We report an experimental study of the shear-induced migration of flexible fibers in suspensions confined between two parallel plates. Non-Brownian fiber suspensions are imaged in a rheo-microscopy setup, where the top and the bottom plates counter-rotate and create a Couette flow. Initially, the fibers are near the bottom plate due to sedimentation. Under shear, the fibers move with the flow and migrate towards the center plane between the two walls. Statistical properties of the fibers, such as the mean values of the positions, orientations, and end-to-end lengths of the fibers, are used to characterize the behaviors of the fibers. A dimensionless parameter Λ_eff, which compares the hydrodynamic shear stress and the fiber stiffness, is used to analyze the effective flexibility of the fibers. The observations show that the fibers that are more likely to bend exhibit faster migration. As Λ_eff increases (softer fibers and stronger shear stresses), the fibers tend to align in the flow direction and the motions of the fibers transition from tumbling and rolling to bending. The bending fibers drift away from the walls to the center plane. Further increasing Λ_eff leads to more coiled fiber shapes, and the bending is more frequent and with larger magnitudes, which leads to more rapid migration towards the center. Different behaviors of the fibers are quantified with Λ_eff, and the structures and the dynamics of the fibers are correlated with the migration.

Directional shear jamming of frictionless ellipses

In this work we study shear reversals of dense non-Brownian suspensions composed of cohesionless elliptical particles. By numerical simulations, we show that a new fragility appears for frictionless ellipses in the flowing states, where particles can flow indefinitely in one direction at applied shear stresses but shear jam in the other direction upon shear stress reversal. This new fragility, absent in the isotropic particle case, is linked to the directional order of the elongated particles at steady shear and its reorientation at shear stress reversal, which forces the suspensions to pass through a more disordered state with an increased number of contacts in which it might get arrested.

Dynamical behaviors of self-propulsion intruder buried in granular materials

Self-propulsion intruder motion in particles is common in the field of biomimetic and exploration instrument development. In this paper, numerical simulations and laboratory experiments are conducted for the spiral upward phenomenal motion of self-propulsion spherical intruder in granular media. Dynamic particle buoyancy and particle Saffman lift are proposed to establish a dynamic model of the intruder under horizontal simple harmonic excitations. The dependencies of the net particle lift force on the horizontal displacement and the local fluidization parameters of granular are discussed. The results show that horizontal displacement of the intruder and the coordination number of particles are jointly determined by the excitation amplitude and frequency, and the intruder starts to rise when they simultaneously reach the critical value. The dynamic particle buoyancy and particle Saffman lift have clarified the ascent mechanics. Meanwhile, the motion trajectory of intruder in space is inverted conical spiral, and the vibration causes the gap filling effect after the local particle fluidization is the mechanism of the intruder floating up motion.

Critical behavior of density-driven and shear-driven reversible-irreversible transitions in cyclically sheared vortices

Random assemblies of particles subjected to cyclic shear undergo a reversible–irreversible transition (RIT) with increasing a shear amplitude d or particle density n, while the latter type of RIT has not been verified experimentally. Here, we measure the time-dependent velocity of cyclically sheared vortices and observe the critical behavior of RIT driven by vortex density B as well as d. At the critical point of each RIT, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{\mathrm {c}}$$\end{document}Bc and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{\mathrm {c}}$$\end{document}dc, 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}$$\tau $$\end{document}τ to reach the steady state shows a power-law divergence. The critical exponent for B-driven RIT is in agreement with that for d-driven RIT and both types of RIT fall into the same universality class as the absorbing transition in the two-dimensional directed-percolation universality class. As d is decreased to the average 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}$$\tau (d)$$\end{document}τ(d) shows a significant drop, indicating a transition or crossover from a loop-reversible state with vortex-vortex collisions to a collisionless point-reversible state. In either regime, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau (d)$$\end{document}τ(d) exhibits a power-law divergence at the same \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d_{\mathrm {c}}$$\end{document}dc with nearly the same exponent.

Gel breakdown in a formulated product via accumulated strain

Almost any formulated product is sufficiently complex that definitive elucidation of all interactions and microstructural evolutions is difficult at best and more likely intractable. Drilling fluids are no exception. Nevertheless, detailed experiment and comparison with simpler systems studied in the literature enable rational pictures to be deduced. We study the breakdown of a gelled formulated product, a drilling fluid, under the action of repeated deformation, i.e. weakly nonlinear oscillation. Our data may be rationalised by postulating that the fluid behaves as an arrested phase separating material whose natural slow structural evolution, aging and coarsening, is accelerated by the imposed sinusoidal strain consistent with previous work on well characterised systems. During the observed evolution the elastic modulus exhibits a maximum which appears correlated with a maximal connected heterogeneity of structure.

Lévy Walks and Path Chaos in the Dispersal of Elongated Structures Moving across Cellular Vortical Flows

In cellular vortical flows, namely arrays of counterrotating vortices, short but flexible filaments can show simple random walks through their stretch-coil interactions with flow stagnation points. Here, we study the dynamics of semirigid filaments long enough to broadly sample the vortical field. Using simulation, we find a surprising variety of long-time transport behavior-random walks, ballistic transport, and trapping-depending upon the filament's relative length and effective flexibility. Moreover, we find that filaments execute Lévy walks whose diffusion exponents generally decrease with increasing filament length, until transitioning to Brownian walks. Lyapunov exponents likewise increase with length. Even completely rigid filaments, whose dynamics is finite dimensional, show a surprising variety of transport states and chaos. Fast filament dispersal is related to an underlying geometry of "conveyor belts." Evidence for these various transport states is found in experiments using arrays of counterrotating rollers, immersed in a fluid and transporting a flexible ribbon.

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.

Random Close Packing as a Dynamical Phase Transition

Sphere packing is an ancient problem. The densest packing is known to be a face-centered cubic (FCC) crystal, with space-filling fraction ϕ_{FCC}=π/sqrt[18]≈0.74. The densest "random packing," random close packing (RCP), is yet ill defined, although many experiments and simulations agree on a value ϕ_{RCP}≈0.64. We introduce a simple absorbing-state model, biased random organization (BRO), which exhibits a Manna class dynamical phase transition between absorbing and active states that has as its densest critical point ϕ_{c_{max}}≈0.64≈ϕ_{RCP} and, like other Manna class models, is hyperuniform at criticality. The configurations we obtain from BRO appear to be structurally identical to RCP configurations from other protocols. This leads us to conjecture that the highest-density absorbing state for an isotropic biased random organization model produces an ensemble of configurations that characterizes the state conventionally known as RCP.

Memory in three-dimensional cyclically driven granular material

We perform experimental and numerical studies of a granular system under cyclic compression to investigate reversibility and memory effects. We focus on the quasistatic forcing of dense systems, which is most relevant to a wide range of geophysical, industrial, and astrophysical problems. We find that soft-sphere simulations with proper stiffness and friction quantitatively reproduce both the translational and rotational displacements of the grains. We then utilize these simulations to demonstrate that such systems are capable of storing the history of previous compressions. While both mean translational and rotational displacements encode such memory, the response is fundamentally different for translations compared to rotations. For translational displacements, this memory of prior forcing depends on the coefficient of static interparticle friction, but rotational memory is not altered by the level of friction.

Barrier-controlled nonequilibrium criticality in reactive particle systems

Nonequilibrium critical phenomena generally exist in many dynamic systems, like chemical reactions and some driven-dissipative reactive particle systems. Here, by using computer simulation and theoretical analysis, we demonstrate the crucial role of the activation barrier on the criticality of dynamic phase transitions in a minimal reactive hard-sphere model. We find that at zero thermal noise, with increasing the activation barrier, the type of transition changes from a continuous conserved directed percolation into a discontinuous dynamic transition by crossing a tricritical point. A mean-field theory combined with field simulation is proposed to explain this phenomenon. The possibility of Ising-type criticality in the nonequilibrium system at finite thermal noise is also discussed.

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.

Symmetry properties of nonlinear hydrodynamic interactions between responsive particles

Two identical particles driven by the same steady force through a viscous fluid may move relative to one another due to hydrodynamic interactions. The presence or absence of this relative translation has a profound effect on the dynamics of a driven suspension consisting of many particles. We consider a pair of particles which, to linear order in the force, do not interact hydrodynamically. If the system possesses an intrinsic property (such as the shape of the particles, their position with respect to a boundary, or the shape of the boundary) which is affected by the external forcing, hydrodynamic interactions that depend nonlinearly on the force may emerge. We study the general properties of such nonlinear response. Analysis of the symmetries under particle exchange and under force reversal leads to general conclusions concerning the appearance of relative translation and the motion's time reversibility. We demonstrate the applicability of the conclusions in three specific examples: (a) two spheres driven parallel to a wall; (b) two deformable objects driven parallel to their connecting line; and (c) two spheres driven along a curved path. The breaking of time reversibility suggests a possible use of nonlinear hydrodynamic interactions to disperse or assemble particles by an alternating force.

Hidden Order Beyond Hyperuniformity in Critical Absorbing States

Disordered hyperuniformity is a description of hidden correlations in point distributions revealed by an anomalous suppression in fluctuations of local density at various coarse-graining length scales. In the absorbing phase of models exhibiting an active-absorbing state transition, this suppression extends up to a hyperuniform length scale that diverges at the critical point. Here, we demonstrate the existence of additional many-body correlations beyond hyperuniformity. These correlations are hidden in the higher moments of the probability distribution of the local density and extend up to a longer length scale with a faster divergence than the hyperuniform length on approaching the critical point. Our results suggest that a hidden order beyond hyperuniformity may generically be present in complex disordered systems.

Transition from steady shear to oscillatory shear rheology of dense suspensions

Recent studies have highlighted that oscillatory and time-dependent shear flows might help increase the flowability of dense suspensions. While most focus has been on cross-flows we here study a simple two-dimensional suspensions where we apply simultaneously oscillatory and stationary shear along the same direction. We first show that the dissipative viscosities in this set-up significantly decrease with an increasing shear-rate magnitude of the oscillations and given that the oscillatory strain is small, in a similar fashion as found previously for cross-flow oscillations. As for cross-flow oscillations, the decrease can be attributed to the large decrease in the number of contacts and an altered microstructure as one transitions from a steady shear to an oscillatory shear dominated rheology. As subresults we find both an extension to the μ(J) rheology, a constitutive relationship between the shear stresses and the shear rate, valid for oscillatory shear flows and that shear-jamming of frictional particles at oscillatory shear dominated flows occurs at higher packing fractions compared to steady shear dominated flows.

Probing intermittency and reversibility in a dense granular suspension under shear using multiply scattered ultrasound

We study the rheology of a dense granular suspension under shear strain with the simultaneous detection of multiply scattered ultrasound through the shear band. At a low shear rate, the dissipation is rate-independent and determined by the frictional contacts between grains. Under quasistatic shear, the stress-strain curve contains elastic loading parts interrupted by stress drops. Such an intermittency is concomitant with some large decorrelation events as measured by the ultrasound probe, sensitive to the position of the grains. Under cyclic shear, the correlations between the scattered ultrasonic waves show that at low shear strain, the grains exhibit reversible motion. Beyond this linear regime, some irreversible motion of the grains is detected. Moreover, the correlation between successive ultrasound signals suggests that some specific rearrangements, which add to the homogeneous flow, take place near the maximum strain.

Hyperuniform Structures Formed by Shearing Colloidal Suspensions

In periodically sheared suspensions there is a dynamical phase transition, characterized by a critical strain amplitude γ_{c}, between an absorbing state where particle trajectories are reversible and an active state where trajectories are chaotic and diffusive. Repulsive nonhydrodynamic interactions between "colliding" particles' surfaces have been proposed as a source of this broken time reversal symmetry. A simple toy model called random organization qualitatively reproduces the dynamical features of this transition. Random organization and other absorbing state models exhibit hyperuniformity, a strong suppression of density fluctuations on long length scales quantified by a structure factor S(q→0)∼q^{α} with α>0, at criticality. Here we show experimentally that the particles in periodically sheared suspensions organize into structures with anisotropic short-range order but isotropic, long-range hyperuniform order when oscillatory shear amplitudes approach γ_{c}.

Hyperuniformity and density fluctuations at a rigidity transition in a model of biological tissues

The suppression of density fluctuations at different length scales is the hallmark of hyperuniformity. Here, we explore the presence of this hidden order in a manybody interacting model of biological tissue, known to exhibit a transition, or sharp crossover, from a solid to a fluid like phase. We show that the density fluctuations in the rigid phase are only suppressed up to a finite lengthscale. This length scale monotonically increases and grows rapidly as we approach the fluid phase reminiscent to divergent behavior at a critical point, such that the system is effectively hyperuniform in the fluid phase. Furthermore, complementary behavior of the structure factor across the critical point also indicates that hyperuniformity found in the fluid phase is stealthy.

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.

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.

Embedding orthogonal memories in a colloidal gel through oscillatory shear

It has recently been shown that in a broad class of disordered systems oscillatory shear training can embed memories of specific shear protocols in relevant physical parameters such as the yield strain. These shear protocols can be used to change the physical properties of the system and memories of the protocol can later be "read" out. Here we investigate shear training memories in colloidal gels, which include an attractive interaction and network structure, and discover that such systems can support memories both along and orthogonal to the training flow direction. We use oscillatory shear protocols to set and read out the yield strain memories and confocal microscopy to analyze the rearranging gel structure throughout the shear training. We find that the gel bonds remain largely isotropic in the shear-vorticity plane throughout the training process suggesting that structures formed to support shear along the training shear plane are also able to support shear along the orthogonal plane. Orthogonal memory extends the usefulness of shear memories to more applications and should apply to many other disordered systems as well.

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.

Directed aging, memory, and nature's greed

Disordered materials are often out of equilibrium and evolve very slowly in a rugged and tortuous energy landscape. This slow evolution, referred to as aging, is deemed undesirable as it often leads to material degradation. However, we show that aging also encodes a memory of the stresses imposed during preparation. Because of inhomogeneous local stresses, the material itself decides how to evolve by modifying stressed regions differently from those under less stress. Because material evolution occurs in response to stresses, aging can be "directed" to produce sought-after responses and unusual functionalities that do not inherently exist. Aging obeys a natural "greedy algorithm" as, at each instant, the material simply follows the path of most rapid and accessible relaxation. Our experiments and simulations illustrate directed aging in examples in which the material's elasticity transforms as desired because of an imposed deformation.

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.

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.