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.

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.

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.

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.

Disordering transitions and peak effect in polydisperse particle systems.

We show numerically that in a binary system of Yukawa particles, a dispersity-driven disordering transition occurs. In the presence of quenched disorder this disordering transition coincides with a marked increase in the depinning threshold, known as a peak effect. We find that the addition of poorly pinned particles can increase the overall pinning in the sample by increasing the amount of topological disorder present. If the quenched disorder is strong enough to create a significant amount of topological disorder in the monodisperse system, addition of a poorly pinned species generates further disorder but does not produce a peak in the depinning force. Our results indicate that for binary mixtures, optimal pinning occurs for topological defect fraction densities from 0.2 to 0.25. For defect densities below this range, the system retains orientational order. We determine the effect of the pinning density, strength, and radius on the depinning peak and find that the peak effect is more pronounced in weakly pinning 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.

Stripes, clusters, and nonequilibrium ordering for bidisperse colloids with repulsive interactions.

We show that two-dimensional bidisperse assemblies of colloids with strictly repulsive interactions exhibit stripe, cluster, and partially crystallized states when driven over a quenched random substrate. The nonequilibrium states on a substrate are significantly more ordered than equilibrium states both with and without substrates. A minimum substrate strength is necessary to induce the nonequilibrium pattern formation. Our results suggest that a combination of driving and quenched disorder offers a new approach to controlling pattern formation in colloid mixtures.

Point-defect dynamics in two-dimensional colloidal crystals.

We study the topological configurations and dynamics of individual point defect vacancies and interstitials in a two-dimensional crystal of colloids interacting via a repulsive Yukawa potential. Our Brownian dynamics simulations show that the diffusion mechanism for vacancy defects occurs in two phases. The defect can glide along the crystal lattice directions, and it can rotate during an excited topological transition configuration to assume a different direction for the next period of gliding. The results for the vacancy defects are in good agreement with recent experiments. For interstitial point defects, which were not studied in the experiments, we find several of the same modes of motion as in the vacancy defect case along with two additional diffusion pathways. The interstitial defects are more mobile than the vacancy defects due to the more two-dimensional nature of the diffusion of the interstitial defects.

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.

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.

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.

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.

Electrophoresis of DNA on a disordered two-dimensional substrate.

We propose a method for electrophoretic separation of DNA in which adsorbed polymers are driven over a disordered two-dimensional substrate which contains attractive sites for the polymers. Using simulations of a model for long polymer chains, we show that the mobility increases with polymer length, in contrast to gel electrophoresis techniques, and that separation can be achieved for a range of length scales. We demonstrate that the separation mechanism relies on steric interactions between polymer segments, which prevent substrate disorder sites from trapping more than one DNA segment each. Since thermal activation does not play a significant role in determining the polymer mobility, band broadening due to diffusion can be avoided in our separation method.

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.

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.

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.

Structure and melting of two-species charged clusters in a parabolic trap.

We consider a system of charged particles interacting with an unscreened Coulomb repulsion in a two-dimensional parabolic confining trap. The static charge on a portion of the particles is twice as large as the charge on the remaining particles. The particles separate into a shell structure with those of greater charge situated farther from the center of the trap. As we vary the ratio of the number of particles of the two species, we find that for certain configurations, the symmetry of the arrangement of the inner cluster of singly charged particles matches the symmetry of the outer ring of doubly charged particles. These matching configurations have a higher melting temperature and a higher thermal threshold for intershell rotation between the species than the nonmatching configurations.

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.

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.

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.