Source author record

C. Reichhardt

C. Reichhardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

55 published item(s)

preprint2026arXiv

Reversible to Irreversible Transitions in Pattern-Forming Systems with Cyclic Interactions

Transitions from reversible to irreversible or fluctuating states above a critical density and shear amplitude have been extensively studied in non-thermal cyclically sheared suspensions and amorphous solids. Here, we propose that the same type of reversible to irreversible transition occurs for a system of particles with competing short-range attraction and long-range repulsion, which can form crystals, stripes, and bubbles as the ratio of attraction to repulsion varies. By oscillating the strength of the attractive part of the potential, we find that the system can organize into either time-periodic states consisting of nondiffusive complex closed orbits, or into a diffusive fluctuating state. A critical point separates these states as a function of the maximum strength of the attraction, oscillation frequency, and particle density. We also find a re-entrant behavior of the reversible state as a function of the strength of the attraction and the oscillation frequency.

preprint2022arXiv

Soliton Motion in Skyrmion Chains

Using a particle-based model we examine the depinning motion of solitons in skyrmion chains in quasi-one dimensional (1D) and two-dimensional (2D) systems containing embedded 1D interfaces. The solitons take the form of a particle or hole in a commensurate chain of skyrmions. Under an applied drive, just above a critical depinning threshold the soliton moves with a skyrmion Hall angle of zero. For higher drives, the entire chain depins, and in a 2D system we observe that both the solitons and chain move at zero skyrmion Hall angle and then transition to a finite skyrmion Hall angle as the drive increases. In a 2D system with a 1D interface that is at an angle to the driving direction, there can be a reversal of the sign of the skyrmion Hall angle from positive to negative. Our results suggest that solitons in skyrmion systems could be used as information carriers in racetrack geometries that would avoid the drawbacks of finite skyrmion Hall angles. The soliton states become mobile at significantly lower drives than the depinning transition of the skyrmion chains themselves.

preprint2019arXiv

Noise spectra in the reversible-irreversible transition in amorphous solids under oscillatory driving

We study the stress fluctuations in simulations of a two-dimensional amorphous solid under a cyclic drive. It is known that this system organizes into a reversible state for small driving amplitudes and remains in an irreversible state for high driving amplitudes, and that a critical driving amplitude separates the two regimes. Here we study the time series of the stress fluctuations below and above the reversible-irreversible transition. In the irreversible regime above the transition, the power spectrum of the stress fluctuations is broad and has a $1/f^α$ shape with $1< α<2$. We find that the low frequency noise power peaks near the stress at which dc yielding occurs, which is consistent with the behavior expected in systems undergoing a non-equilibrium phase transition.

preprint2019arXiv

Quenched Dynamics of Artificial Spin Ice: Coarsening versus Kibble-Zurek

Artificial spin ices are ideal frustrated model systems in which to explore or design emergent phenomena with unprecedented characterization of the constituent degrees of freedom. In square spin ice, violations of the ice rule are topological excitations essential to the kinetics of the system, providing an ideal testbed for studying the dynamics of such defects under varied quench rates. In this work we describe the first test of the Kibble-Zurek mechanism and critical coarsening in colloidal square and colloidal hexagonal ice under quenches from a weakly interacting liquid state into a strongly interacting regime. As expected, for infinitely slow quenches, the system is defect free, while for increasing quench rate, an increasing number of defects remain in the sample. For square ice, we find regimes in which the defect population decreases as a power law with decreasing quench rate. A detailed scaling analysis shows that for a wide range of parameters, including quench rates that are accessible by experiments, the behavior is described by critical coarsening rather than by the Kibble-Zurek mechanism, since the defect-defect interactions are long ranged. For quenches closer to the critical point, however, there can be a competition between the two mechanisms.

preprint2016arXiv

Emergent Geometric Frustration of Artificial Magnetic Skyrmion Crystals

Magnetic skyrmions have been receiving growing attention as potential information storage and magnetic logic devices since an increasing number of materials have been identified that support skyrmion phases. Explorations of artificial frustrated systems have led to new insights into controlling and engineering new emergent frustration phenomena in frustrated and disordered systems. Here, we propose a skyrmion spin ice, giving a unifying framework for the study of geometric frustration of skyrmion crystals in a non-frustrated artificial geometrical lattice as a consequence of the structural confinement of skyrmions in magnetic potential wells. The emergent ice rules from the geometrically frustrated skyrmion crystals highlight a novel phenomenon in this skyrmion system: emergent geometrical frustration. We demonstrate how skyrmion crystal topology transitions between a non-frustrated periodic configuration and a frustrated ice-like ordering can also be realized reversibly. The proposed artificial frustrated skyrmion systems can be annealed into different ice phases with an applied current induced spin-transfer torque, including a long range ordered ice rule obeying ground state, as-relaxed random state, biased state and monopole state. The spin-torque reconfigurability of the artificial skyrmion ice states, difficult to achieve in other artificial spin ice systems, is compatible with standard spintronic device fabrication technology, which makes the semiconductor industrial integration straightforward.

preprint2015arXiv

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 non-monotonic behavior that is correlated with distinct changes in the spatial-temporal 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 that 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 which how power-law distributed sizes that are similar to the avalanche distributions observed for non-active disk systems near the jamming transition.

preprint2015arXiv

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 compressions, 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 even smectic 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.

preprint2015arXiv

Doped Colloidal Artificial Ice

We examine square and kagome artificial spin ice for colloids confined in arrays of double-well traps. Unlike magnetic artificial spin ices, colloidal and vortex artificial spin ice realizations allow creation of doping sites through double occupation of individual traps. We find that doping square and kagome ice geometries produces opposite effects. For square ice, doping creates local excitations in the ground state configuration that produce a local melting effect as the temperature is raised. In contrast, the kagome ice ground state can absorb the doping charge without generating non-ground-state excitations, while at elevated temperatures the hopping of individual colloids is suppressed near the doping sites. These results indicate that in the square ice, doping adds degeneracy to the ordered ground state and creates local weak spots, while in the kagome ice, which has a highly degenerate ground state, doping locally decreases the degeneracy and creates local hard regions.

preprint2015arXiv

Quantized Transport for a Skyrmion Moving on a Two-Dimensional Periodic Substrate

We study the dynamics and velocity-force curves of a skyrmion moving over a two-dimensional periodic substrate using simulations of a particle-based skyrmion model, focusing on the role of the non-dissipative Magnus term. In the overdamped limit, the skyrmion depins into a sliding state with a Hall angle of zero. When a Magnus term is included, the Hall angle is nonzero in the absence of a substrate. On a periodic substrate, the Hall angle varies with the drive amplitude. Due to the substrate symmetry the Hall angle does not change continuously with drive, but forms a series of discrete steps at rational ratios of the skyrmion velocity components perpendicular and parallel to the drive direction, when the skyrmion motion locks to symmetry directions of the substrate for fixed intervals of the drive amplitude. On each step, the Hall angle is constant and the skyrmion motion is orderly. Transitions between locked phases generate pronounced cusps in the velocity-force curves, as well as regions of negative differential mobility. The number of observable locking steps increases when the relative strength of the Magnus term increases. We find an overshoot phenomena where the Hall angle exceeds the clean limit value, as well as an acceleration effect where the skyrmion moves faster over a substrate than through a clean sample. These effects are robust for different types of periodic substrates. With a simple model for a skyrmion interacting with a single pin, we can capture the behavior of the Hall angle. The Magnus term induces a curvature in the skrymion orbit as it moves through the pin, resulting in a side step phenomenon that decreases with increasing drive. When the Magnus term is large, the range of impact parameters that permit the skyrmion to be trapped by the pin is small, which is a reasons that strong Magnus force effects reduce the pinning in skyrmion systems.

preprint2015arXiv

Reversible Ratchet Effects for Vortices in Conformal Pinning Arrays

A conformal transformation of a uniform triangular pinning array produces a structure called a conformal crystal which preserves the six-fold ordering of the original lattice but contains a gradient in the pinning density. Here we use numerical simulations to show that vortices in type-II superconductors driven with an ac drive over gradient pinning arrays produce the most pronounced ratchet effect over a wide range of parameters for a conformal array, while square gradient or random gradient arrays with equivalent pinning densities give reduced ratchet effects. In the conformal array, the larger spacing of the pinning sites in the direction transverse to the ac drive permits easy funneling of interstitial vortices for one driving direction, producing the enhanced ratchet effect. In the square array, the transverse spacing between pinning sites is uniform, giving no asymmetry in the funneling of the vortices as the driving direction switches, while in the random array, there are numerous easy-flow channels present for either direction of drive. We find multiple ratchet reversals in the conformal arrays as a function of vortex density and ac amplitude, and correlate the features with a reversal in the vortex ordering, which is greater for motion in the ratchet direction. The enhanced conformal pinning ratchet effect can also be realized for colloidal particles moving over a conformal array, indicating the general usefulness of conformal structures for controlling the motion of particles.

preprint2015arXiv

Transverse ac Driven and Geometric Ratchet Effects for Vortices in Conformal Crystal Pinning Arrays

With a series of conformal arrays it is possible to create asymmetric substrates, and it was previously shown that when an ac drive is applied parallel to the asymmetry direction, a pronounced ratchet effect occurs with a net dc flow of vortices in the same direction as the ac drive. Here we show that when the ac drive is applied perpendicular to the substrate asymmetry direction, it is possible to realize a transverse ratchet effect where a net dc flow of vortices is generated perpendicular to the ac drive. The conformal transverse ratchet effect is distinct from previous versions of transverse ratchets in that it occurs due to the generation of non-Gaussian transverse vortex velocity fluctuations by the plastic motion of vortices, so that the system behaves as a noise correlation ratchet. The transverse ratchet effect is much more pronounced in the conformal arrays than in random gradient arrays, and is absent in square gradient arrays due the different nature of the vortex flow in each geometry. We show that a series of reversals can occur in the transverse ratchet effect due to changes in the vortex flow across the pinning gradient as a function of vortex filling, pinning strength, and ac amplitude. We also consider the case where a dc drive applied perpendicular to the substrate asymmetry direction generates a net flow of vortices perpendicular to the dc drive, producing a geometric or drift ratchet that arises due to non-Gaussian dynamically generated fluctuations. The drift ratchet is more efficient than the ac driven ratchet and also exhibits a series of reversals for varied parameters. Our results should be general to a wide class of systems undergoing nonequilibrium dynamics on conformal substrates, such as colloidal particles on optical traps.

preprint2014arXiv

Avalanches and disorder-induced criticality in artificial spin ices

We show that both square and kagome artificial spin ice systems exhibit disorder-induced nonequilibrium phase transitions, with power law avalanche distributions at the critical disorder level. The different nature of geometrical frustration in the two lattices produces distinct types of critical avalanche behavior. For the square ice, the avalanches involve the propagation of locally stable domain walls separating the two polarized ground states, and the scaling collapse agrees with an interface depinning mechanism. In contrast, avalanches in the fully frustrated kagome ice exhibit pronounced branching behaviors that resemble those found in directed percolation. The kagome ice also shows an interesting crossover in the power-law scaling of the avalanches at low disorder. Our results show that artificial spin ices are ideal systems in which to study nonequilibrium critical point phenomena.

preprint2014arXiv

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 separation 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 run length.

preprint2014arXiv

Pinning, Ordering, and Dynamics of Vortices in Conformal Crystal and Gradient Pinning Arrays

We numerically investigate magnetization, pinning, ordering, and dynamics of vortices interacting with pinning arrangements which have a density gradient. We focus on conformal crystal structures obtained by conformally transforming a spatially uniform periodic array, as well as non-conformal gradient structures and structures with quasiperiodic order. The conformal structures feature a density gradient and local ordering. Using magnetization simulations we find that conformal pinning arrays exhibit enhanced pinning compared to non-conformal gradient arrays as well as compared to random, periodic, and quasiperiodic arrays, for a broad range of fields. The effectiveness of conformal arrays arises from the continuum of length scales introduced into the arrays by the conformal transformation, allowing for a broad range of local commensuration effects. At higher vortex fillings above the range of conformal effectiveness, we show that a non-conformal rectangular gradient array exhibits strong pinning due to a novel commensuration effect and vortex ordering. Using transport simulations where vortices are driven along the gradient and at an angle to the gradient, we confirm the effectiveness of conformal pinning at increasing the critical current. For a rotated drive, the gradient arrays produce a strong vortex guidance effect in the direction perpendicular to the gradient.

preprint2013arXiv

Comment on "Giant Plasticity of a Quantum Crystal"

In their Letter, Haziot et al. [Phys. Rev. Lett. 110 (2013) 035301] report a novel phenomenon of giant plasticity for hcp Helium-4 quantum crystals. They assert that Helium-4 exhibits mechanical properties not found in classical plasticity theory. Specifically, they examine high-quality crystals as a function of temperature and applied strain, where the shear modulus reaches a plateau and dissipation becomes close to zero; both quantities are reported to be independent of stress and strain, implying a reversible dissipation process and quantum tunneling. In this Comment, we show that these signatures can be explained with a classical model of thermally activated dislocation glide without the need to invoke quantum tunneling or dissipationless motion. Recently, we proposed a dislocation glide model in solid Helium-4 containing the dissipation contribution in the presence of other dislocations with qualitatively similar behavior [Zhou et al., Philos. Mag. Lett. 92 (2012) 608].

preprint2013arXiv

Comparing the Dynamics of Skyrmions and Superconducting Vortices

Vortices in type-II superconductors have attracted enormous attention as ideal systems in which to study nonequilibrium collective phenomena, since the self-ordering of the vortices competes with quenched disorder and thermal effects. Dynamic effects found in vortex systems include depinning, nonequilibrium phase transitions, creep, structural order-disorder transitions, and melting. Understanding vortex dynamics is also important for applications of superconductors which require the vortices either to remain pinned or to move in a controlled fashion. Recently, topological defects called skyrmions have been realized experimentally in chiral magnets. Here we highlight similarities and differences between skyrmion dynamics and vortex dynamics. Many of the previous ideas and experimental setups that have been applied to superconducting vortices can also be used to study skyrmions. We also discuss some of the differences between the two systems, such as the potentially large contribution of the Magnus force in the skyrmion system that can dramatically alter the dynamics and transport properties.

preprint2013arXiv

Dynamic Regimes for Driven Colloidal Particles on a Periodic Substrate at Commensurate and Incommensurate Fillings

We examine colloidal particles driven over a periodic muffin tin substrate using numerical simulations. In the absence of a driving force this system exhibits a rich variety of commensurate and incommensurate static phases in which topological defects can form domain walls, ordered stripes, superlattices, and disordered patchy regimes as a function of the filling fraction. When an external drive is applied, these different static phases generate distinct dynamical responses. At incommensurate fillings the flow generally occurs in the form of localized pulses or solitons correlated with the motion of the topological defect structures. We also find dynamic transitions between different types of moving states that are associated with changes in the velocity force curves, structural transitions in the topological defect arrangements, and modifications of the velocity distributions and particle trajectories. We find that the 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 along the stripe direction rather than along the driving direction. Structural dynamic transitions can also occur within the pinned regimes when the applied drive causes distortions of the interstitially pinned particles while the particles trapped at pinning sites remain immobile.

preprint2013arXiv

Dynamics and Separation of Circularly Moving Particles in Asymmetric 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.

preprint2013arXiv

Frustrated colloidal ordering and fully packed loops in arrays of optical traps

We propose that a system of colloidal particles interacting with a honeycomb array of optical traps that each contain three wells can be used to realize a fully packed loop model. One of the phases in this system can be mapped to Baxter's three-coloring problem, offering an easily accessible physical realization of this problem. As a function of temperature and interaction strength, we find a series of phases, including long range ordered loop or stripe states, stripes with sliding symmetries, random packed loop states, and disordered states in which the loops break apart. Our geometry could be constructed using ion trap arrays, BEC vortices in optical traps, or magnetic vortices in nanostructured superconductors.

preprint2013arXiv

Vortex States in Archimedean Tiling Pinning Arrays

We numerically study vortex ordering and pinning in Archimedean tiling substrates composed of square and triangular plaquettes. The two different plaquettes become occupied at different vortex densities, producing commensurate peaks in the magnetization at non-integer matching fields. We find that as the field increases, in some cases the fraction of occupied pins can decrease due to the competition between fillings of the different plaquette types. We also identify a number of different types of vortex orderings as function of field at integer and non-integer commensurate fillings.

preprint2013arXiv

Vortex Transport and Pinning Effectiveness in Conformal Pinning Arrays

We examine the current driven dynamics for vortices interacting with conformal crystal pinning arrays and compare to the dynamics of vortices driven over random pinning arrays. We find that the pinning is enhanced in the conformal arrays over a wide range of fields, consistent with previous results from flux gradient-driven simulations. At fields above this range, the effectiveness of the pinning in the moving vortex state can be enhanced in the random arrays compared to the conformal arrays, leading to crossing of the velocity-force curves.

preprint2012arXiv

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 particles in the driven region break away or decouple from the particles in the bulk, producing a shear band. For a bidisperse assembly, we find that 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.

preprint2012arXiv

Domain and Stripe Formation Between Hexagonal and Square Ordered Fillings of Colloidal Particles on Periodic Pinning Substrates

Using large scale numerical simulations, we examine the ordering of colloidal particles on square periodic two-dimensional muffin-tin substrates consisting of a flat surface with localized pinning sites. We show that when there are four particles per pinning site, the particles adopt a hexagonal ordering, while for five particles per pinning site, a square ordering appears. For fillings between four and five particles per pinning site, we identify a rich variety of distinct ordering regimes, including disordered grain boundaries, crystalline stripe structures, superlattice orderings, and disordered patchy arrangements. We characterize the different regimes using Voronoi analysis, energy dispersion, and ordering of the domains. We show that many of the boundary formation features we observe occur for a wide range of other fillings. Our results demonstrate that grain boundary tailoring can be achieved with muffin-tin periodic pinning substrates.

preprint2012arXiv

Jamming in Systems With Quenched Disorder

We numerically study the effect of adding quenched disorder in the form of randomly placed pinning sites on jamming transitions in systems that jam at a well defined point J in the clean limit. Quenched disorder decreases the jamming density and introduces a depinning threshold. The onset of a finite threshold coincides with point J at the lowest pinning densities, but for higher pinning densities there is always a finite threshold even well below jamming. We find that proximity to point J strongly affects the transport curves and noise fluctuations, and observe a change from plastic behavior below jamming, where the system is highly heterogeneous, to elastic depinning above jamming. Many of the general features we find are related to other systems containing quenched disorder, including the peak effect observed in vortex systems.

preprint2012arXiv

Static and Dynamic Phases for Vortex Matter with Attractive Interactions

Exotic vortex states with long range attraction and short range repulsion have recently been proposed to arise in superconducting hybrid structures and multi-band superconductors. Using large scale simulations we examine the static and dynamic properties of such vortex states interacting with random and periodic pinning. In the absence of pinning this system does not form patterns but instead completely phase separates. When pinning is present there is a transition from inhomogeneous to homogeneous vortex configurations similar to a wetting phenomenon. Under an applied drive, a dynamical dewetting process can occur from a strongly pinned homogeneous state into pattern forming states. We show that a signature of the exotic vortex interactions under transport measurements is a robust double peak feature in the differential conductivity curves.

preprint2012arXiv

Statics and Dynamics of Yukawa Cluster Crystals on Ordered Substrates

We examine the statics and dynamics of particles with repulsive Yukawa interactions in the presence of a two-dimensional triangular substrate for fillings of up to twelve particles per potential minimum. We term the ordered states Yukawa cluster crystals and show that they are distinct from the colloidal molecular crystal states found at low fillings. As a function of substrate and interaction strength at fixed particle density we find a series of novel crystalline states that we characterize using the structure factor. For fillings greater than four, shell and ring structures form at each potential minimum and can exhibit sample-wide orientational order. A disordered state can appear between ordered states as the substrate strength varies. Under an external drive, the onsets of different orderings produce clear changes in the critical depinning force, including a peak effect phenomenon that has generally only previously been observed in systems with random substrates. We also find a rich variety of dynamic ordering transitions that can be observed via changes in the structure factor and features in the velocity-force curves. The dynamical states encompass a variety of moving structures including one-dimensional stripes, smectic ordering, polycrystalline states, triangular lattices, and symmetry locking states. Despite the complexity of the system, we identify several generic features of the dynamical phase transitions which we map out in a series of phase diagrams. Our results have implications for the structure and depinning of colloids on periodic substrates, vortices in superconductors and Bose-Einstein condensates, Wigner crystals, and dusty plasmas.

preprint2012arXiv

Strongly Enhanced Vortex Pinning by Conformal Crystal Arrays

Conformal crystals are non-uniform structures created by a conformal transformation of regular two-dimensional lattices. We show that gradient-driven vortices interacting with a conformal pinning array exhibit substantially stronger pinning effects over a much larger range of field than found for random or periodic pinning arrangements. The pinning enhancement is partially due to matching of the critical flux gradient with the pinning gradient, but the preservation of the sixfold ordering in the conformally transformed hexagonal lattice plays a crucial role. Our results can be generalized to a wide class of gradient-driven interacting particle systems such as colloids on optical trap arrays.

preprint2011arXiv

Active Matter on Asymmetric Substrates

For collections of particles in a thermal bath interacting with an asymmetric substrate, it is possible for a ratchet effect to occur where the particles undergo a net dc motion in response to an ac forcing. Ratchet effects have been demonstrated in a variety of systems including colloids as well as magnetic vortices in type-II superconductors. Here we examine the case of active matter or self-driven particles interacting with asymmetric substrates. Active matter systems include self-motile colloidal particles undergoing catalysis, swimming bacteria, artificial swimmers, crawling cells, and motor proteins. We show that a ratchet effect can arise in this type of system even in the absence of ac forcing. The directed motion occurs for certain particle-substrate interaction rules and its magnitude depends on the amount of time the particles spend swimming in one direction before turning and swimming in a new direction. For strictly Brownian particles there is no ratchet effect. If the particles reflect off the barriers or scatter from the barriers according to Snell's law there is no ratchet effect; however, if the particles can align with the barriers or move along the barriers, directed motion arises. We also find that under certain motion rules, particles accumulate along the walls of the container in agreement with experiment. We also examine pattern formation for synchronized particle motion. We discuss possible applications of this system for self-assembly, extracting work, and sorting as well as future directions such as considering collective interactions and flocking models.

preprint2011arXiv

Bidirectional Sorting of Flocking Particles in the Presence of Asymmetric Barriers

We numerically demonstrate bidirectional sorting of flocking particles interacting with an array of asymmetric barriers. Each particle aligns with the average swimming direction of its neighbors according to the Vicsek model and experiences additional steric interactions as well as repulsion from the fixed barriers. We show that particles preferentially localize to one side of the barrier array over time, and that the direction of this rectification can be reversed by adjusting the particle-particle exclusion radius or the noise term in the equations of motion. These results provide a conceptual basis for isolation and sorting of single- and multi-cellular organisms which move collectively according to flocking-type interaction rules.

preprint2011arXiv

Dynamically Induced Locking and Unlocking Transitions in Driven Layered Systems with Quenched Disorder

Using numerical simulations, we examine a simple model of two or more coupled one-dimensional channels of driven particles with repulsive interactions in the presence of quenched disorder. We find that this model exhibits a remarkably rich variety of dynamical behavior as a function of the strength of the quenched disorder, coupling between channels, and external drive. For weaker disorder, the channels depin in a single step. For two channels we find dynamically induced decoupling transitions that result in coexisting pinned and moving phases as well as moving decoupled phases where particles in both channels move at different average velocities and slide past one another. Decoupling can also be induced by changing the relative strength of the disorder in neighboring channels. At higher drives, we observe a dynamical recoupling or locking transition into a state with no relative motion between the channels.This recoupling produces unusual velocity-force signatures, including negative differential conductivity. The depinning threshold shows distinct changes near the decoupling and coupling transitions and exhibits a peak effect phenomenon of the type that has been associated with transitions from elastic to plastic flow in other systems. We map several dynamic phase diagrams showing the coupling-decoupling transitions and the regions in which hysteresis occurs. We also examine the coexistence regime for channels with unequal amounts of quenched disorder. For multiple channels, multiple coupling and decoupling transitions can occur; however, many of the general features found for the two channel system are still present. Our results should be relevant to depinning in layered geometries in systems such as vortices in layered or nanostructured superconductors and Wigner or colloidal particles confined in nano-channels; they are also relevant to the general understating of plastic flow.

preprint2011arXiv

Hysteresis and Return Point Memory in Artificial Spin Ice Systems

We investigate hysteresis loops and return point memory for artificial square and kagome spin ice systems by cycling an applied bias force and comparing microscopic effective spin configurations throughout the hysteresis cycle. Return point memory loss is caused by motion of individual defects in kagome ice or of grain boundaries in square ice. In successive cycles, return point memory is recovered rapidly in kagome ice. Memory is recovered more gradually in square ice due to the extended nature of the grain boundaries. Increasing the amount of quenched disorder increases the defect density but also enhances the return point memory since the defects become trapped more easily.

preprint2011arXiv

Jamming in Granular Polymers

We examine the jamming transition in a two-dimensional granular polymer system using compressional simulations. The jamming density ϕ_c decreases with increasing length of the granular chain due to the formation of loop structures, in excellent agreement with recent experiments. The jamming density can be further reduced in mixtures of granular chains and granular rings, also as observed in experiment. We show that the nature of the jamming in granular polymer systems has pronounced differences from the jamming behavior observed for polydisperse two-dimensional disk systems at Point J. This result provides further evidence that there is more than one type of jamming transition.

preprint2011arXiv

Multi-Step Ordering in Kagome and Square Artificial Spin Ice

We show that in colloidal models of artificial kagome and modified square ice systems, a variety of ordering and disordering regimes occur as a function of biasing field, temperature, and colloid-colloid interaction strength, including ordered monopole crystals, biased ice rule states, thermally induced ice rule ground states, biased triple states, and disordered states. We describe the lattice geometries and biasing field protocols that create the different states and explain the formation of the states in terms of sublattice switching thresholds. For a system prepared in a monopole lattice state, we show that a sequence of different orderings occurs for increasing temperature. Our results also explain several features observed in nanomagnetic artificial ice systems under an applied field.

preprint2011arXiv

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 non-driven 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.

preprint2011arXiv

Structural Transitions and Dynamical Regimes for Directional Locking of Particles Driven over Periodic Substrates

We numerically investigate collective ordering and disordering effects for vortices in type-II superconductors interacting with square and triangular substrate arrays under a dc drive that is slowly rotated with respect to the fixed substrate. A series of directional locking transitions occur when the particle motion locks to symmetry directions of the substrate, producing a series of steps in the velocity-force curves. The locking transitions coincide with structural transitions between triangular, square, smectic, or disordered particle arrangements, which can be identified using the structure factor. We show that the widths of the locking steps pass through local minima and maxima as a function of the ratio of the number of particles to the number of substrate minima. Unlike a static system, where matching effects occur for simple integer commensuration ratios, our system exhibits dynamical commensuration effects where an integer number of particle chains flow between one-dimensional lines of substrate minima. As the system enters and exits the locking steps, order-disorder transitions in the structure of the moving particle assembly occur. We identify two distinct symmetry locking regimes as a function of substrate strength which produce different locking step characteristics. For weak substrates, all the particles are in motion and a portion of the particles flow through the substrate minima, leading to structural transitions at certain driving angles. For strong substrates, some particles are permanently pinned while the remaining particles flow around them. At the crossover between these two regimes of substrate strength, the locking steps are destroyed due to the onset of chaotic plastic flow which produces pronounced changes in the transport characteristics. Similar effects occur for colloidal particles driven over square and triangular substrate arrays.

preprint2010arXiv

Anisotropic Sliding Dynamics, Peak Effect, and Metastability in Stripe Systems

A variety of soft and hard condensed matter systems are known to form stripe patterns. Here we use numerical simulations to analyze how such stripe states depin and slide when interacting with a random substrate and with driving in different directions with respect to the orientation of the stripes. Depending on the strength and density of the substrate disorder, we find that there can be pronounced anisotropy in the transport produced by different dynamical flow phases. We also find a disorder-induced "peak effect" similar to that observed for superconducting vortex systems, which is marked by a transition from elastic depinning to a state where the stripe structure fragments or partially disorders at depinning. Under the sudden application of a driving force, we observe pronounced metastability effects similar to those found near the order-disorder transition associated with the peak effect regime for three-dimensional superconducting vortices. The characteristic transient time required for the system to reach a steady state diverges in the region where the flow changes from elastic to disordered. We also find that anisotropy of the flow persists in the presence of thermal disorder when thermally-induced particle hopping along the stripes dominates. The thermal effects can wash out the effects of the quenched disorder, leading to a thermally-induced stripe state. We map out the dynamical phase diagram for this system, and discuss how our results could be explored in electron liquid crystal systems, type-1.5 superconductors, and pattern-forming colloidal assemblies.

preprint2010arXiv

Commensurability, Jamming, and Dynamics for Vortices in Funnel Geometries

We use numerical simulations to examine vortex states and dynamics in periodic funnel geometries where a drive is applied in the easy flow direction. We show that this system exhibits a number of different commensurability effects when the vortex configurations match to both the periodicity of the array and the geometry of the funnels. The vortex configurations in this system are generally different from those observed for single isolated triangular superconducting samples due to the coupling of vortices in adjacent funnels. At certain matching fields, peaks in the critical current are absent due to the particular vortex configurations that occur at these fields. We find that the overall depinning force increases with increasing vortex density as a result of the enhanced vortex-vortex interactions caused by a crowding effect at the funnel tips. When a system becomes less mobile as a result of increased particle interactions, it is said to exhibit a jamming behavior. Under an applied drive we observe a series of elastic and plastic vortex flow phases which produce pronounced features such as jumps or dips in the transport curves. In all of the flow phases, only one vortex can pass through the funnel tip at a time due to the vortex-vortex repulsion forces. As a consequence of this constraint, we observe the remarkable result that the sum of the vortex velocities at a fixed drive remains nearly constant with increasing magnetic field B rather than increasing linearly. This result is similar to the behavior of sand in an hourglass. We also show how noise fluctuations can be used to distinguish the different flow phases. Our results should be readily generalizable to other systems of particles flowing in periodic funnel geometries, such as colloids or Wigner crystals.

preprint2010arXiv

Dynamical Ordering and Directional Locking For Particles Moving Over Quasicrystalline Substrates

We use molecular dynamics simulations to study the driven phases of particles such as vortices or colloids moving over a decagonal quasiperiodic substrate. In the regime where the pinned states have quasicrystalline ordering, the driven phases can order into moving square or smectic states, or into states with aligned rows of both square and triangular tiling which we term dynamically induced Archimedean-like tiling. We show that when the angle of the drive is varied with respect to the substrate, directional locking effects occur where the particle motion locks to certain angles. It is at these locking angles that the dynamically induced Archimedean tiling appears. We also demonstrate that the different dynamical orderings and locking phases show pronounced changes as a function of filling fraction.

preprint2010arXiv

Structural Transitions, Melting, and Intermediate Phases for Stripe and Clump Forming Systems

We numerically examine the properties of a two-dimensional system of particles which have competing long range repulsive and short range attractive interactions as a function of density and temperature. For increasing density, there are well defined transitions between a low density clump phase, an intermediate stripe phase, an anticlump phase, and a high density uniform phase. To characterize the transitions between these phases we propose several measures which take into account the different length scales in the system. For increasing temperature, we find an intermediate phase that is liquid-like on the short length scale of interparticle spacing but solid-like on the larger length scale of the clump, stripe, or anticlump pattern. This intermediate phase persists over the widest temperature range in the stripe state when the local particle lattice within an individual stripe melts well below the temperature at which the entire stripe structure breaks down, and is characterized by intra-stripe diffusion of particles without inter-stripe diffusion. This is followed at higher temperatures by the onset of inter-stripe diffusion in an anisotropic diffusion phase, and then by breakup of the stripe structure. We identify the transitions between these regimes through diffusion, specific heat, and energy fluctuation measurements, and find that within the intra-stripe liquid regime, the excess entropy goes into disordering the particle arrangements within the stripe rather than affecting the stripe structure itself. The clump and anticlump phases also show multiple temperature-induced diffusive regimes which are not as pronounced as those of the stripe phase.

preprint2010arXiv

The Effect of Pinning on Drag in Coupled One-Dimensional Channels of Particles

We consider a simple model for examining the effects of quenched disorder on drag consisting of particles interacting via a Yukawa potential that are placed in two coupled one-dimensional channels. The particles in one channel are driven and experience a drag from the undriven particles in the second channel. In the absence of pinning, for a finite driving force there is no pinned phase; instead, there are two dynamical regimes of completely coupled or locked flow and partially coupled flow. When pinning is added to one or both channels, we find that a remarkably rich variety of dynamical phases and drag effects arise that can be clearly identified by features in the velocity force curves. The presence of quenched disorder in only the undriven channel can induce a pinned phase in both channels. Above the depinning transition, the drag on the driven particles decreases with increasing pinning strength, and for high enough pinning strength, the particles in the undriven channel reach a reentrant pinned phase which produces a complete decoupling of the channels. We map out the dynamic phase diagrams as a function of pinning strength and the density of pinning in each channel. Our results may be relevant for understanding drag coupling in 1D Wigner crystal phases, and the effects we observe could also be explored using colloids in coupled channels produced with optical arrays, vortices in nanostructured superconductors, or other layered systems where drag effects arise.

preprint2009arXiv

Coherent and Incoherent Vortex Flow States in Crossed Channels

We examine vortex flow states in periodic square pinning arrays with one row and one column of pinning sites removed to create an easy flow crossed channel geometry. When a drive is simultaneously applied along both major symmetry axes of the pinning array such that vortices move in both channels, a series of coherent flow states develop in the channel intersection at rational ratios of the drive components in each symmetry direction when the vortices can cross the intersection without local collisions. The coherent flow states are correlated with a series of anomalies in the velocity force curves, and in some cases can produce negative differential conductivity. The same general behavior could also be realized in other systems including colloids, particle traffic in microfluidic devices, or Wigner crystals in crossed one-dimensional channels.

preprint2009arXiv

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.

preprint2009arXiv

Fluctuation Theorem in Driven Nonthermal Systems with Quenched Disorder

We demonstrate that the fluctuation theorem of Gallavotti and Cohen can be used to characterize the class of dynamics that arises in nonthermal systems of collectively interacting particles driven over random quenched disorder. By observing the frequency of entropy-destroying trajectories, we show that there are specific dynamical regimes near depinning in which this theorem holds. Hence the fluctuation theorem can be used to characterize a significantly wider class of nonequilibrium systems than previously considered. We discuss how the fluctuation theorem could be tested in specific systems where noisy dynamics appear at the transition from a pinned to a moving phase such as in vortices in type-II superconductors, magnetic domain walls, and dislocation dynamics.

preprint2009arXiv

Fluctuations, Jamming, and Yielding for a Driven Probe Particle in Disordered Disk Assemblies

Using numerical simulations we examine the velocity fluctuations of a probe particle driven with constant force through a two-dimensional disordered assembly of disks which has a well-defined jamming point J at a density of ϕ_J=0.843. As ϕincreases toward ϕ_J, the average velocity of the probe particle decreases and the velocity fluctuations show an increasingly intermittent or avalanchelike behavior. When the system is within a few percent of the jamming density, the velocity distributions are exponential, while when the system is less than a percent away from jamming, the velocity distributions have a non-exponential or power law character. The velocity power spectra exhibit a crossover from a Lorentzian form to a 1/f shape near jamming. We extract a correlation exponent νwhich is in good agreement with recent shear simulations. For ϕ> ϕ_J, there is a critical threshold force F_c that must be applied for the probe particle to move through the sample which increases with increasing ϕ. The onset of the probe motion above ϕ_J occurs via a local yielding of the particles around the probe particle which we term a local shear banding effect.

preprint2009arXiv

Pattern Switching and Polarizability for Colloids in Optical Trap Arrays

We show that colloidal molecular crystal states interacting with a periodic substrate, such as an optical trap array, and a rotating external field can undergo a rapid pattern switching in which the orientation of the crystal changes. In some cases, a martensitic-like symmetry switching occurs. It is also possible to create a polarized state where the colloids in each substrate minima develop a director field which smoothly rotates with the external drive, similar to liquid crystal behavior. These results open the possibility for creating novel types of devices using photonic band gap materials, and should be generalizable to a variety of other condensed matter systems with multiple particle trapping.

preprint2009arXiv

Random Organization and Plastic Depinning

We provide evidence that plastic depinning falls into the same class of phenomena as the random organization which was recently studied in periodically driven particle systems [L. Corte et al., Nature Phys. 4, 420 (2008)]. In the plastic flow system, the pinned regime corresponds to the quiescent 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 colloidal systems and superconducting vortex systems with random disorder.

preprint2009arXiv

Shear Banding and Spatiotemporal Oscillations in Vortex Matter in Nanostructured Superconductors

We propose a simple nanostructured pinning array geometry where a rich variety of complex vortex shear banding phenomena can be realized. A single row of pinning sites is removed from a square pinning array. Shear banding effects arise when vortex motion in the pin-free channel nucleates motion of vortices in the surrounding pinned regions, creating discrete steps in the vortex velocity profile away from the channel. Near the global depinning transition, the width of the band of moving vortices undergoes oscillations or fluctuations that can span the entire system. We use simulations to show that these effects should be observable in the transport properties of the system. Similar large oscillations and shear banding effects are known to occur for sheared complex fluids in which different dynamical phases coexist.

preprint2009arXiv

Switching and Jamming Transistor Effect for Vortex Matter in Honeycomb Pinning Arrays with ac Drives

We show that a remarkable variety of dynamical phenomena, including switching, polarization, symmetry locking, and dynamically induced pinning, can occur for vortices in type-II superconductors in the presence of a honeycomb pinning array and an ac or combined ac and dc drive. These effects occur at the second matching field where there are two vortices per pinning site, and arise due to the formation of vortex dimer states in the interstitial regions of the honeycomb array. The orientation of the pinned and moving vortex dimers can be controlled externally by the application of a drive. We term this a polarization effect and demonstrate that it can lock or unlock the vortex motion into different symmetry directions of the underlying pinning lattice. If the moving vortices are locked into one direction, the motion can be switched into a different direction by applying an additional bias drive, producing sharp jumps in the transverse and longitudinal velocities. Further, the dc vortex motion in one direction can be controlled directly by application of a force in the perpendicular direction. When the moving dimers reorient, we find a remarkable dynamical pinning effect in which the dimers jam when they become perpendicular to the easy flow direction of the pinning lattice. Since application of an external field can be used to switch off the vortex flow, we term this a jamming transistor effect. These effects do not occur in triangular pinning arrays due to the lack of the n-merization of the vortices in this case. The switching and dynamical pinning effects demonstrated here may be useful for the creation of new types of fluxtronic devices.

preprint2008arXiv

Collective Sliding States for Colloidal Molecular Crystals

We study the driving of colloidal molecular crystals over periodic substrates such as those created with optical traps. The n-merization that occurs in the colloidal molecular crystal states produces a remarkably rich variety of distinct dynamical behaviors, including polarization effects within the pinned phase and the formation of both ordered and disordered sliding phases. Using computer simulations, we map the dynamic phase diagrams as a function of substrate strength for dimers and trimers on a triangular substrate, and correlate features on the phase diagram with transport signatures.

preprint2003arXiv

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.

preprint2003arXiv

Temperature and ac Effects on Charge Transport in Metallic Arrays of Dots

We investigate the effects of finite temperature, dc pulse, and ac drives on the charge transport in metallic arrays using numerical simulations. For finite temperatures there is a finite conduction threshold which decreases linearly with temperature. Additionally we find a quadratic scaling of the current-voltage curves which is independent of temperature for finite thresholds. These results are in excellent agreement with recent experiments on 2D metallic dot arrays. We have also investigated the effects of an ac drive as well as a suddenly applied dc drive. With an ac drive the conduction threshold decreases for fixed frequency and increasing amplitude and saturates for fixed amplitude and increasing frequency. For sudden applied dc drives below threshold we observe a long time power law conduction decay.

preprint2002arXiv

Novel Colloidal Crystalline States on Two Dimensional Periodic Substrates

We show using numerical simulations that a rich variety of novel colloidal crystalline states are realized on square and triangular two dimensional periodic substrates which can be experimentally created using crossed laser arrays. When there are more colloids than potential substrate minima, multiple colloids are trapped at each substrate minima and act as a single particle with a rotational degree of freedom, giving rise to a new type of orientational order. We call these states colloidal molecular crystals. A two-step melting can also occur in which individual colloidal molecules initially rotate, destroying the overall orientational order, followed by the onset of inter-well colloidal hopping.

preprint1999arXiv

Dynamic Vortex Phases and Pinning in Superconductors with Twin Boundaries

We investigate the pinning and driven dynamics of vortices interacting with twin boundaries using large scale molecular dynamics simulations on samples with near one million pinning sites. For low applied driving forces, the vortex lattice orients itself parallel to the twin boundary and we observe the creation of a flux gradient and vortex free region near the edges of the twin boundary. For increasing drive, we find evidence for several distinct dynamical flow phases which we characterize by the density of defects in the vortex lattice, the microscopic vortex flow patterns, and orientation of the vortex lattice. We show that these different dynamical phases can be directly related to microscopically measurable voltage - current V(I) curves and voltage noise. By conducting a series of simulations for various twin boundary parameters we derive several vortex dynamic phase diagrams.

preprint1996arXiv

Plastic Flow, Voltage Bursts, and Vortex Avalanches in Superconductors

We use large-scale parallel simulations to compute the motion of superconducting magnetic vortices during avalanches triggered by small field increases. We find that experimentally observable voltage bursts correspond to pulsing vortex movement along branched channels or winding chains, and relate vortex flow images to features of statistical distributions. As pin density is increased, a crossover occurs from interstitial motion in narrow easy-flow winding channels with typical avalanche sizes, to pin-to-pin motion in broad channels, characterized by a very broad distribution of sizes. Our results are consistent with recent experiments.

preprint1996arXiv

Vortex Plastic Flow, $B(x,y,H(t)), M(H(t)), J_c(B(t))$, Deep in the Bose Glass and Mott-Insulator Regimes

We present simulations of flux-gradient-driven superconducting vortices interacting with strong columnar pinning defects as an external field $H(t)$ is quasi-statically swept from zero through a matching field $B_ϕ$. We analyze several measurable quantities, including the local flux density $ B(x,y,H(t))$, magnetization $M(H(t))$, critical current $J_{c}(B(t))$, and the individual vortex flow paths. We find a significant change in the behavior of these quantities as the local flux density crosses $B_ϕ$, and quantify it for many microscopic pinning parameters. Further, we find that for a given pin density $J_c(B)$ can be enhanced by maximizing the distance between the pins for $ B < B_ϕ $.