Dynamic phases, stratification, laning, and pattern formation for driven bidisperse disk systems in the presence of quenched disorder.

Using numerical simulations, we examine the dynamics of driven two-dimensional bidisperse disks flowing over quenched disorder. The system exhibits a series of distinct dynamical phases as a function of applied driving force and packing fraction, including a phase-separated state as well as a smectic state with liquid-like or polycrystalline features. At low driving forces, we find a clogged phase with an isotropic density distribution, while at intermediate driving forces the disks separate into bands of high and low densities with either liquid-like or polycrystalline structure in the high-density bands. In addition to the density phase separation, we find that in some cases there is a fractionation of the disk species, particularly when the disk size ratio is large. The species phase-separated regimes form a variety of patterns such as large disks separated by chains of smaller disks. Our results show that the formation of laning states can be enhanced by tuning the ratio of disk radius of the two species, due to the clumping of small disks in the interstitial regions between the large disks. This system could be experimentally realized using sterically interacting colloidal particles suspended in a viscous fluid driven over random pinning arrays or granular matter suspended in fluid moving over a random landscape.

Dynamic phases, clustering, and chain formation for driven disk systems in the presence of quenched disorder.

We numerically examine the dynamic phases and pattern formation of two-dimensional monodisperse repulsive disks driven over random quenched disorder. We show that there is a series of distinct dynamic regimes as a function of increasing drive, including a clogged or pile-up phase near depinning, a homogeneous disordered flow state, and a dynamically phase separated regime consisting of high-density crystalline regions surrounded by a low density of disordered disks. At the highest drives the disks arrange into one-dimensional moving chains. The phase separated regime has parallels with the phase separation observed in active matter systems, but arises from a distinct mechanism consisting of the combination of nonequilibrium fluctuations with density-dependent mobility. We discuss the pronounced differences between this system and previous studies of driven particles with longer-range repulsive interactions moving over random substrates, such as superconducting vortices or electron crystals, where dynamical phase separation and distinct one-dimensional moving chains are not observed. Our results should be generic to a broad class of systems in which the particle-particle interactions are short ranged, such as sterically interacting colloids or Yukawa particles with strong screening driven over random pinning arrays, superconducting vortices in the limit of small penetration depths, or quasi-two-dimensional granular matter flowing over rough landscapes.

Dynamic Control of Topological Defects in Artificial Colloidal Ice.

We demonstrate the use of an external field to stabilize and control defect lines connecting topological monopoles in spin ice. For definiteness we perform Brownian dynamics simulations with realistic units mimicking experimentally realized artificial colloidal spin ice systems, and show how defect lines can grow, shrink or move under the action of direct and alternating fields. Asymmetric alternating biasing forces can cause the defect line to ratchet in either direction, making it possible to precisely position the line at a desired location. Such manipulation could be employed to achieve mobile information storage in these metamaterials.

Clogging and jamming transitions in periodic obstacle arrays.

We numerically examine clogging transitions for bidisperse disks flowing through a two-dimensional periodic obstacle array. We show that clogging is a probabilistic event that occurs through a transition from a homogeneous flowing state to a heterogeneous or phase-separated jammed state where the disks form dense connected clusters. The probability for clogging to occur during a fixed time increases with increasing particle packing and obstacle number. For driving at different angles with respect to the symmetry direction of the obstacle array, we show that certain directions have a higher clogging susceptibility. It is also possible to have a size-specific clogging transition in which one disk size becomes completely immobile while the other disk size continues to flow.

Collective transport for active matter run-and-tumble disk systems on a traveling-wave substrate.

We examine numerically the transport of an assembly of active run-and-tumble disks interacting with a traveling-wave substrate. We show that as a function of substrate strength, wave speed, disk activity, and disk density, a variety of dynamical phases arise that are correlated with the structure and net flux of disks. We find that there is a sharp transition into a state in which the disks are only partially coupled to the substrate and form a phase-separated cluster state. This transition is associated with a drop in the net disk flux, and it can occur as a function of the substrate speed, maximum substrate force, disk run time, and disk density. Since variation of the disk activity parameters produces different disk drift rates for a fixed traveling-wave speed on the substrate, the system we consider could be used as an efficient method for active matter species separation. Within the cluster phase, we find that in some regimes the motion of the cluster center of mass is in the opposite direction to that of the traveling wave, while when the maximum substrate force is increased, the cluster drifts in the direction of the traveling wave. This suggests that swarming or clustering motion can serve as a method by which an active system can collectively move against an external drift.

Avalanches, plasticity, and ordering in colloidal crystals under compression.

Using numerical simulations we examine colloids with a long-range Coulomb interaction confined in a two-dimensional trough potential undergoing dynamical compression. As the depth of the confining well is increased, the colloids move via elastic distortions interspersed with intermittent bursts or avalanches of plastic motion. In these avalanches, the colloids rearrange to minimize their colloid-colloid repulsive interaction energy by adopting an average lattice constant that is isotropic despite the anisotropic nature of the compression. The avalanches take the form of shear banding events that decrease or increase the structural order of the system. At larger compression, the avalanches are associated with a reduction of the number of rows of colloids that fit within the confining potential, and between avalanches the colloids can exhibit partially crystalline or anisotropic ordering. The colloid velocity distributions during the avalanches have a non-Gaussian form with power-law tails and exponents that are consistent with those found for the velocity distributions of gliding dislocations. We observe similar behavior when we subsequently decompress the system, and find a partially hysteretic response reflecting the irreversibility of the plastic events.

Softening of stressed granular packings with resonant sound waves.

We perform numerical simulations of a two-dimensional bidisperse granular packing subjected to both a static confining pressure and a sinusoidal dynamic forcing applied by a wall on one edge of the packing. We measure the response experienced by a wall on the opposite edge of the packing and obtain the resonant frequency of the packing as the static or dynamic pressures are varied. Under increasing static pressure, the resonant frequency increases, indicating a velocity increase of elastic waves propagating through the packing. In contrast, when the dynamic amplitude is increased for fixed static pressure, the resonant frequency decreases, indicating a decrease in the wave velocity. This occurs both for compressional and for shear dynamic forcing and is in agreement with experimental results. We find that the average contact number Zc at the resonant frequency decreases with increasing dynamic amplitude, indicating that the elastic softening of the packing is associated with a reduced number of grain-grain contacts through which the elastic waves can travel. We image the excitations created in the packing and show that there are localized disturbances or soft spots that become more prevalent with increasing dynamic amplitude. Our results are in agreement with experiments on glass bead packings and earth materials such as sandstone and granite and may be relevant to the decrease in elastic wave velocities that has been observed to occur near fault zones after strong earthquakes, in surficial sediments during strong ground motion, and in structures during earthquake excitation.

Collective transport properties of driven Skyrmions with random disorder.

We use particle-based simulations to examine the static and driven collective phases of Skyrmions interacting with random quenched disorder. We show that nondissipative effects due to the Magnus term reduce the depinning threshold and strongly affect the Skyrmion motion and the nature of the dynamic phases. The quenched disorder causes the Hall angle to become drive dependent in the moving Skyrmion phase, while different flow regimes produce distinct signatures in the transport curves. For weak disorder, the Skyrmions form a pinned crystal and depin elastically, while for strong disorder the system forms a pinned amorphous state that depins plastically. At high drives the Skyrmions can dynamically reorder into a moving crystal, with the onset of reordering determined by the strength of the Magnus term.

Active microrheology in active matter systems: Mobility, intermittency, and avalanches.

We examine the mobility and velocity fluctuations of a driven particle moving through an active matter bath of self-mobile disks for varied density or area coverage and varied activity. We show that the driven particle mobility can exhibit nonmonotonic behavior that is correlated with distinct changes in the spatiotemporal structures that arise in the active media. We demonstrate that the probe particle velocity distributions exhibit specific features in the different dynamic regimes and identify an activity-induced uniform crystallization that occurs for moderate activity levels and is distinct from the previously observed higher activity cluster phase. The velocity distribution in the cluster phase has telegraph noise characteristics produced when the probe particle moves alternately through high-mobility areas that are in the gas state and low-mobility areas that are in the dense phase. For higher densities and large activities, the system enters what we characterize as an active jamming regime. Here the probe particle moves in intermittent jumps or avalanches that have power-law-distributed sizes that are similar to the avalanche distributions observed for nonactive disk systems near the jamming transition.

Casimir effect in active matter systems.

We numerically examine run-and-tumble active matter particles in Casimir geometries composed of two finite parallel walls. We find that there is an attractive force between the two walls of a magnitude that increases with increasing run length. The attraction exhibits an unusual exponential dependence on the wall separation, and it arises due to a depletion of swimmers in the region between the walls by a combination of the motion of the particles along the walls and a geometric shadowing effect. This attraction is robust as long as the wall length is comparable to or smaller than the swimmer run length, and is only slightly reduced by the inclusion of steric interactions between swimmers. We also examine other geometries and find regimes in which there is a crossover from attraction to repulsion between the walls as a function of wall separation and wall length.

Aspects of jamming in two-dimensional athermal frictionless systems.

In this work we provide an overview of jamming transitions in two dimensional systems focusing on the limit of frictionless particle interactions in the absence of thermal fluctuations. We first discuss jamming in systems with short range repulsive interactions, where the onset of jamming occurs at a critical packing density and where certain quantities show a divergence indicative of critical behavior. We describe how aspects of the dynamics change as the jamming density is approached and how these dynamics can be explored using externally driven probes. Different particle shapes can produce jamming densities much lower than those observed for disk-shaped particles, and we show how jamming exhibits fragility for some shapes while for other shapes this is absent. Next we describe the effects of long range interactions and jamming behavior in systems such as charged colloids, vortices in type-II superconductors, and dislocations. We consider the effect of adding obstacles to frictionless jamming systems and discuss connections between this type of jamming and systems that exhibit depinning transitions. Finally, we discuss open questions such as whether the jamming transition in all these different systems can be described by the same or a small subset of universal behaviors, as well as future directions for studies of jamming transitions in two dimensional systems, such as jamming in self-driven or active matter systems.

Dynamics and separation of circularly moving particles in asymmetrically patterned arrays.

There are many examples of driven and active matter systems containing particles that exhibit circular motion with different chiralities, such as swimming bacteria near surfaces or certain types of self-driven colloidal particles. Circular motion of passive particles can also be induced with an external rotating drive. Here we examine particles that move in circles and interact with a periodic array of asymmetric L-shaped obstacles. We find a series of dynamical phases as a function of swimming radius, including regimes where the particle motion is rectified, producing a net dc motion. The direction of the rectification varies with the swimming radius, permitting the separation of particles with different swimming radii. Particles with the same swimming radius but different chirality can also move in different directions over the substrate and be separated. The rectification occurs for specific windows of swimming radii corresponding to periodic orbits in which the particles interact one or more times with the barriers per rotation cycle. The rectification effects are robust against the addition of thermal or diffusive effects, and are in some cases even enhanced by these effects.

Colloidal lattice shearing and rupturing with a driven line of particles.

We examine the dynamics of two-dimensional colloidal systems using numerical simulations of a system with a drive applied to a thin region in the middle of the sample to produce a local shear. For a monodisperse colloidal assembly, we find a well-defined decoupling transition separating a regime of elastic motion from a plastic phase where the driven particles break away or decouple from the bulk particles and produce a shear band. For a bidisperse assembly, the onset of a bulk disordering transition coincides with the broadening of the shear band. We identify several distinct dynamical regimes that are correlated with features in the velocity-force curves. As a function of bidispersity, the decoupling force shows a nonmonotonic behavior associated with features in the noise fluctuations, power spectra, and bulk velocity profiles. When pinning is added in the bulk, we find that the shear band regions can become more localized, causing a decoupling of the driven particles from the bulk particles. For a system with thermal noise and no pinning, the shear band region becomes more extended and the average velocity of the driven particles drops at the thermal disordering transition of the bulk system.

Dynamic regimes for driven colloidal particles on a periodic substrate at commensurate and incommensurate fillings.

We numerically examine colloidal particles driven over a muffin tin substrate. Previous studies of this model identified a variety of commensurate and incommensurate static phases in which topological defects can form domain walls, ordered stripes, superlattices, or disordered patchy regimes as a function of the filling fraction. Here, we show that the addition of an external drive to these static phases can produce distinct dynamical responses. At incommensurate fillings the flow occurs in the form of localized pulses or solitons correlated with topological defect structures. Transitions between different modes of motion can occur as a function of increasing drive. We measure the average particle velocity for specific ranges of external drive and show that changes in the velocity response correlate with changes in the topological defect arrangements. We also demonstrate that in the different dynamic phases, the particles have distinct trajectories and velocity distributions. Dynamic transitions between ordered and disordered flows exhibit hysteresis, while in strongly disordered regimes there is no hysteresis and the velocity-force curves are smooth. When stripe patterns are present, transport can occur at an angle to the driving direction.

Structural transitions and dynamical regimes for directional locking of vortices and colloids driven over periodic substrates.

We examine collective dynamical locking effects for superconducting vortices and colloids interacting with square and triangular substrate arrays under a slowly rotated dc drive. A rich variety of lattice configurations associated with a series of steps in the velocity-force curves occur during the locking transitions. These include triangular, square, smectic, and disordered particle arrangements that can be identified using the structure factor. We show that the step widths vary with the ratio of the number of particles to the number of pinning sites. Unlike a static system, where matching effects occur at simple integer commensuration ratios, we find dynamical commensuration effects that arise when an integer number of particle chains flow between rows of pins. We identify two distinct types of locking as a function of substrate strength, distinguished by whether the particles flow along or between the pinning rows.

Positive and negative drag, dynamic phases, and commensurability in coupled one-dimensional channels of particles with Yukawa interactions.

We introduce a simple model consisting of two or three coupled one-dimensional channels of particles with Yukawa interactions. For the two-channel system, when an external drive is applied only to the top or primary channel, we find a transition from locked flow, where particles in both channels move together, to decoupled flow, where the particles in the secondary or undriven channel move at a slower velocity than the particles in the primary or driven channel. Pronounced commensurability effects in the decoupling transition occur when the ratio of the number of particles in the top and bottom channels is varied, and the coupling of the two channels is enhanced when this ratio is an integer or a rational fraction. Near the commensurate fillings, we find additional features in the velocity-force curves caused by the slipping of individual vacancies or incommensurations in the secondary channels. For three coupled channels, when only the top channel is driven we find a remarkably rich variety of distinct dynamic phases, including multiple decoupling and recoupling transitions. These transitions produce pronounced signatures in the velocity response of each channel. We also find regimes where a negative drag effect can be induced in one of the nondriven channels. The particles in this channel move in the opposite direction from the particles in the driven channel due to the mixing of the two different periodic frequencies produced by the discrete motion of the particles in the two other channels. In the two-channel system, we also demonstrate a ratchet effect for the particles in the secondary channel when an asymmetric drive is applied to the primary channel. This ratchet effect is similar to that observed in superconducting vortex systems when there is a coupling between two different species of vortices.

Random organization and plastic depinning.

We provide evidence that the general phenomenon of plastic depinning can be described as an absorbing phase transition, and shows the same features as the random organization which was recently studied in periodically driven particle systems [L. Corte, Nature Phys. 4, 420 (2008)]. In the plastic flow system, the pinned regime corresponds to the absorbing state and the moving state corresponds to the fluctuating state. When an external force is suddenly applied, the system eventually organizes into one of these two states with a time scale that diverges as a power law at a nonequilibrium transition. We propose a simple experiment to test for this transition in systems with random disorder.

Creating artificial ice states using vortices in nanostructured superconductors.

We demonstrate that it is possible to realize vortex ice states that are analogous to square and kagome ice. With numerical simulations, we show that the system can be brought into a state that obeys either global or local ice rules by applying an external current according to an annealing protocol. We explore the breakdown of the ice rules due to disorder in the nanostructure array and show that in square ice, topological defects appear along grain boundaries, while in kagome ice, individual defects appear. We argue that the vortex system offers significant advantages over other artificial ice systems.

Enhancing mixing and diffusion with plastic flow.

We use numerical simulations to examine two-dimensional particle mixtures that strongly phase separate in equilibrium. When the system is externally driven in the presence of quenched disorder, plastic flow occurs in the form of meandering and strongly mixing channels. In some cases, this can produce a fast and complete mixing of previously segregated particle species, as well as an enhancement of transverse diffusion even in the absence of thermal fluctuations. We map the mixing phase diagram as a function of external driving and quenched disorder parameters.