You are an expert at summarizing long articles. Proceed to summarize the following text: 
the creation of new types of crystalline or partially ordered states and the dynamics of assemblies of interacting particles have attracted much attention both in terms of the basic science of self - assembly and dynamic pattern formation as well as for applications utilizing self - assembly processes . one of the most extensively studied systems exhibiting this behavior is assemblies of colloidal particles , where the equilibrium structures can be tuned by changing the directionality of the colloid - colloid interactions @xcite . since it can be difficult to control and tune the exact form of the interaction , another approach is to use colloids with well defined interactions that are placed on some type of ordered substrate . optical trapping techniques are one of the most common methods of creating periodic substrates for colloids @xcite . studies of colloidal ordering and melting for one dimensional ( 1d ) periodic substrate arrays have revealed ordered colloidal crystalline structures as well as smectic type structures where the colloids are crystalline in one direction and liquidlike in the other @xcite . these experiments show that the substrate strength strongly influences the type of colloidal structure that forms and that as the substrate strength increases , the resulting enhancement of fluctuations can induce a transition from crystalline to smectic order @xcite . numerous different colloidal crystalline structures can also appear on 1d substrate arrays of fixed periodicity when the colloid density is varied @xcite . more recent studies addressed colloidal ordering on two - dimensional ( 2d ) periodic substrates @xcite . in these systems , the filling factor @xmath0 is defined as the number of colloids per substrate minimum . for integer fillings @xmath1 , the colloids in each minimum can form an effective rigid @xmath2-mer , such as a dimer or trimer @xcite . the @xmath2-mers have an orientational degree of freedom , and depending on the effective interaction between neighboring @xmath2-mers , all the @xmath2-mers may align into a ferromagnetically ordered state , sit perpendicularly to their neighbors in an antiferromagnetically ordered state , or form other orientationally ordered states . for square and triangular substrate arrays , @xmath2-mer states have been studied up to @xmath3 @xcite ; however , it is not known what structures would form at higher fillings when the simple picture of rigid @xmath2-mers no longer applies . studies of the ordering of bidisperse colloidal assemblies with two different charges on 2d periodic substrates produced novel ordered phases , while a pattern switching could be induced by application of an external field @xcite . similar pattern switching also occurs for colloids with monodisperse charges under external driving @xcite . the colloidal @xmath2-mer states have been termed `` colloidal molecular crystals , '' and for conditions under which they loose their orientational ordering , they are referred to as `` colloidal plastic crystals . '' colloidal molecular crystals appear for integer fillings @xmath1 . at fractional fillings such as @xmath4 , @xmath5 , or @xmath6 , it is possible for ordered composite states to form containing two coexisting species of @xmath2-mers ; however , for other fractional fillings , the system is frustrated and the @xmath2-mer states are disordered @xcite . other studies have shown that novel orderings appear when the 2d substrate array has quasicrystalline order @xcite .    once colloidal crystals have formed on a substrate , the driven dynamics can be explored by applying an external driving field to the sample . a variety of dynamical locking phases can occur in which the colloids preferentially flow along symmetry directions of the substrate @xcite . as the filling fraction @xmath0 is varied , a series of peaks in the critical force needed to depin the colloids occurs at integer values of @xmath0 indicating the existence of commensurability effects @xcite . recent experiments with strongly interacting colloids on two periodic arrays show that kink - type dynamics can occur near @xmath7 @xcite . it would be interesting to explore higher fillings where new types of dynamics could emerge @xcite .    many of the same types of phenomena found in colloidal molecular crystals can also be realized for other systems that can be modeled as interacting particles in the presence of a 2d periodic substrate . for example , the antiferromagnetic ordering of dimer colloidal molecular crystals on a square substrate was reproduced using vortices in bose - einstein condensates confined by optical traps with two vortices per trap @xcite as well as with vortices trapped by large pinning sites in type - ii superconductors @xcite . experimental and numerical studies of molecular ordering on periodic substrates show similar orderings @xcite . other systems where similar states could be realized include classical electrons or dusty plasmas with some form of substrate as well as crystalline cold atoms on optical lattices .    to our knowledge , previous studies of yukawa interacting colloidal molecular crystals have focused only on systems with up to four colloids per trap . for the case of three colloids per trap , only a limited number of studies have considered the dynamics , and even in this limit there are several new features that we describe for the first time in this work . we show that at high fillings , the rigid @xmath2-mer picture breaks down and new cluster and ring states form . several general features of the statics and dynamics emerge at these larger fillings which are independent of the specific filling . one of the key findings in our work is the development of orientationally ordered shell structures at higher fillings for the 2d arrays . the development of particle shell structures was studied previously for repulsive particles in isolated individual traps , in systems that include classical charges @xcite , wigner islands @xcite , dusty plasmas @xcite , colloids @xcite , charged balls in traps @xcite , and vortices in confined geometries @xcite . the shell ordering in these systems can be altered by the shape or type of trap as well as by competing attractive and repulsive interactions between the particles @xcite . in our system , the particles within each shell can exhibit an additional ordering due to interactions with particles in the neighboring traps . we show that as a function of substrate strength , a rich variety of colloidal cluster crystals can be created depending on the filling , and that certain shell structures are more stable than others . the number of particles in the shells and the number of shells depend on the substrate strength , and in certain cases the particles form ring structures instead of shells . we also find reentrant ordered phases as a function of substrate strength . for weak and strong substrates , the system is ordered ; however , for substrates of intermediate strength , the system is generally disordered . we find that the depinning threshold can show distinct changes at the boundaries between these different phases as a function of substrate strength , and remarkably , we find that in some cases it is possible for the depinning force to decrease with increasing substrate strength . the velocity - force curves contain clear signatures of the different phases corresponding to different modes of depinning . we consider a 2d system with periodic boundary conditions in the @xmath8 and @xmath9 directions . the sample contains @xmath10 particles interacting with a triangular sinusoidal substrate containing @xmath11 potential minima . we focus on the case where the number of colloids @xmath10 is an integer multiple of @xmath12 such that @xmath13 , where @xmath0 is the filling factor and @xmath2 is an integer . a single particle @xmath14 responds to forces according to the following overdamped equation of motion : @xmath15 the particle - particle interaction force @xmath16 and the particle interactions are of a yukawa form , @xmath17 , where @xmath18 , @xmath19 is the solvent dielectric constant , @xmath20 is the effective charge , @xmath21 is the screening length , @xmath22 is the unit of length which is of order a micron , and @xmath23 . we model the substrate as a triangular lattice , @xmath24\ ] ] where @xmath25 , @xmath26 , @xmath27 , and @xmath28 . @xmath29 is the amplitude of the substrate potential . the initial colloid positions are obtained by simulated annealing , where the colloids start in a high temperature molten state and are gradually cooled to a finite temperature or to @xmath30 . we have also considered other initialization procedures , including placing a commensurate number of colloids in each substrate minima and letting the system relax . this procedure gives very similar initial states as the simulated annealing , especially for strong substrates ; however , for weaker substrates the configuration obtained through simulated annealing is generally more disordered , as we discuss in more detail later . after annealing , we investigate the transport properties by applying an external drive @xmath31 . we increment @xmath32 from @xmath33 to some final value in steps of @xmath34 , and we wait @xmath35 simulation time steps between increments to avoid transient effects . the depinning threshold is obtained by measuring the average colloid velocity @xmath36 . the value of @xmath34 is modified according to the strength of the substrate . for weaker substrates we use smaller values of @xmath37 in order to obtain an accurate depinning threshold . we first consider the cases @xmath39 , 4 , and @xmath40 . earlier numerical studies addressed the statics and dynamics for @xmath41 and @xmath42 and showed that for strong triangular substrates , @xmath41 produced a herringbone state that transitioned under an applied drive into an aligned ferromagnetic state where the particles move in 1d channels @xcite . for weaker substrates a triangular lattice formed and the depinning was elastic , while for intermediate substrate strengths the depinning was plastic and the herringbone state depinned into a fluctuating state that reordered at higher drives . there was a strong increase in the critical depinning force between the elastic and plastic depinning regimes similar to the peak effect phenomenon found at the transition between elastic and plastic depinning for vortices in type - ii superconductors . for @xmath43 the pinned state was always triangular and for weak substrates the particles depinned into either a moving crystal or a moving smectic state . at these low fillings it was also shown that for weak substrates the critical depinning force @xmath44 , as expected for elastic depinning , while for strong substrates , @xmath45 , as expected for single particle depinning .     for different substrate strengths of ( a ) @xmath46 and ( b ) @xmath47 . in ( b ) the trimer states that form at each trap are aligned . ( c ) the structure factor @xmath48 corresponding to panel ( a ) at @xmath49 . ( d ) @xmath48 corresponding to panel ( b ) at @xmath47 shows long range order and some features at larger @xmath50 values corresponding to the smaller real space length scale of the trimers . ] we now consider several features at @xmath43 that were not reported in the earlier studies , including the velocity - force curves and an additional dynamical phase . figures  [ fig1](a ) and [ fig1](b ) show the particle configurations at @xmath43 for @xmath46 and @xmath47 , respectively , while figs . [ fig1](c ) and [ fig1](d ) show the corresponding structure factor @xmath48 . for @xmath46 in fig . [ fig1](a , c ) , a triangular lattice forms with one dominant length scale that is the distance between the particles . for @xmath51 in fig . [ fig1](b , d ) , each substrate minimum clearly captures three particles and the structure factor has the signature of two length scales . there is a well defined hexagonal pattern of maxima in @xmath48 surrounding the origin in fig . [ fig1](d ) that was not present in fig . [ fig1](c ) which is associated with the longer length scale of the substrate minima locations . the hexagonal void arrangement at larger @xmath50 values in fig .  [ fig1](d ) is associated with the smaller length scale of the trimer colloid arrangement within each substrate minimum . in general , as @xmath29 increases further , the structure of the lattice remains the same as shown in fig . [ fig1](c , d ) except that the trimers gradually reduce in size . vs external drive @xmath52 for the @xmath43 system in fig . [ fig1 ] for @xmath53 , 2 , 3 , 4 , 5 , and 6 , from left to right . for @xmath54 the three different phases are labeled for the pinned ( p ) , random flow ( r ) , and moving smectic ( ms ) phases . for the @xmath55 curve there is no random phase but two different elastic flow phases @xmath56 and @xmath57 occur which are labeled with lines . ]     and @xmath54 for the @xmath58 sample from fig . the trajectories are rapidly changing with time and mixing of the particles occurs . ( c ) the corresponding @xmath48 shows sixfold peaks at small @xmath50 due to the overall triangular ordering of the particles on the substrate as well as a ring structure at larger @xmath50 due to the liquidlike nature of the particle positions . ( b ) particle positions ( dots ) and trajectories ( lines ) over a fixed time interval in the moving smectic phase at @xmath59 for the same sample , where the particles follow well - defined winding quasi-1d paths . ( d ) the corresponding @xmath48 shows more ordering , including more pronounced peaks along @xmath60 that are consistent with a smectic phase . ]    in fig . [ fig2 ] we plot the average particle velocity @xmath61 versus the external drive @xmath52 for @xmath53 , 2 , 3 , 4 , 5 , and @xmath62 . for the @xmath54 curve we have labeled the different dynamical phases . the pinned phase ( p ) transitions into a rapidly fluctuating phase we term the random phase ( r ) . the random phase is characterized by rapid mixing of the particles , as illustrated in fig . [ fig3](a ) for @xmath54 at @xmath63 , where the particle trajectories show significant excursions in the direction transverse to the drive and are continuously changing with time . in fig . [ fig3](b ) , the corresponding @xmath48 taken from a snapshot of the particle positions in the random phase shows features of an anisotropic liquid . there is triangular ordering in @xmath48 at smaller @xmath50 values due to the underlying triangular substrate ; however , the higher order peaks of fig . [ fig1](d ) that indicated long range order in the pinned state are replaced in fig . [ fig3](b ) by a ring structure characteristic of a liquid . we also find a larger number of peaks along @xmath60 compared to @xmath64 due to the anisotropy induced by the @xmath8-direction drive . for higher drives the system transitions into a state with no transverse diffusion , as illustrated in fig . [ fig3](b , d ) for a sample with @xmath65 at @xmath66 . in this state , which we term the moving smectic ( ms ) phase , the motion is confined to winding quasi-1d flows as shown in fig . [ fig3](b ) . the plot of @xmath48 in fig . [ fig3](d ) indicates that the ms state has stronger partial ordering than the random phase but that there is still a tendency for the system to be more ordered along the @xmath60 direction . the r to ms transition is correlated with a decrease in @xmath61 with increasing @xmath52 , producing a negative @xmath67 signature known as negative differential conductivity ( ndc ) . simulations and experiments on superconducting vortices moving over periodic pinning sites have revealed ndc at transitions from random or turbulent flows to ordered 1d flows @xcite . the ndc in the vortex system is much more pronounced than the ndc we observe in fig . [ fig2 ] , and this is likely due to the difference in the type of substrate used in the two systems . for the vortex system , artificially fabricated pinning sites produce a short range muffin - tin potential that permits some vortices to sit completely outside of the pinning sites in the flat interstitial regions . in the random phase , a large number of vortices move through the interstitial areas and do not interact directly with the pins , but at the transition to 1d flow , these vortices suddenly fall into flowing channels that pass through the pinning sites , abruptly increasing the effective amount of drag exerted by the pinning sites on the vortices . in the sinusoidal triangular substrate we consider here , there are no interstitial regions and the moving particles are always experiencing some drag from the pinning . additionally , in the vortex system ndc was only associated with fillings just above @xmath7 , the first matching field @xcite , while in fig . [ fig2 ] we find no ndc until @xmath68 . we note that the ndc for the triangular substrate is much more robust than ndc in the vortex system since it appears for a much wider range of fillings , including incommensurate fillings . in a recent vortex experiment @xcite ndc did not occur until @xmath69 ; however , this effect could not be reproduced in vortex simulations using only muffin - tin type pinning potentials . this could mean that for the particular system considered in ref . @xcite , despite the fact that the pinning sites are localized , longer range interactions between the pins or distortions of the substrate may have caused the substrate to behave more like the sinusoidal type of substrate we consider here rather than like a muffin - tin potential . sample in fig . [ fig2 ] at @xmath55 . ( a ) in phase ef1 at @xmath70 , the particles form a triangular lattice and flow through winding channels but there is no particle mixing . ( b ) in phase ef2 at @xmath71 , the trajectories become straight . the transition between these two phases appears as a kink in @xmath61 near @xmath72 in fig . [ fig2 ] . ]    for @xmath73 we find the same three phases described for @xmath65 , with the transitions between phases shifting to higher @xmath32 with increasing @xmath29 . in fig . [ fig2 ] we identify two distinct phases observed for a weaker substrate with @xmath55 . here the random phase is absent and the system depins elastically , with no sharp jump in @xmath61 at depinning as found for @xmath74 when the flow above depinning is plastic . even though the depinning transition for @xmath75 is elastic , we still find transitions between distinct dynamical phases that can be identified through signatures in the velocity - force curve . in fig . [ fig2 ] we mark the two elastic flow phases ef1 and ef2 that appear on either side of a discontinuity in @xmath76 centered at @xmath72 . in both phases the particle lattice is triangular and has an @xmath48 similar to that shown in fig .  [ fig1](b ) . the ef1 and ef2 phases also appear for @xmath47 with the kink in @xmath61 shifted up to higher @xmath32 . in fig . [ fig4](a ) we plot the particle trajectories in phase ef1 for the @xmath75 system . the triangular particle lattice moves in a zig - zag pattern in order to avoid passing over the potential maxima of the triangular substrate . at higher drives the particles no longer avoid the potential maxima and enter phase ef2 where they move in the direction of drive , as illustrated in fig . [ fig4](b ) .     for samples with @xmath77 . ( a , b ) at @xmath46 the particles form a triangular lattice . ( c , d ) at @xmath78 a quadrimer state forms . the quadrimers have a ferromagnetic alignment . ( e , f ) at @xmath79 the quadrimer state has no orientational ordering , producing ring structures in @xmath48 . ]    for @xmath77 , when the substrate is weak a triangular lattice of particles forms , as illustrated in fig . [ fig5](a , b ) for a sample with @xmath49 . weak substrates still affect the lattice by breaking rotational symmetry and causing the particle lattice to preferentially orient in a direction determined by the filling factor . at @xmath80 , one of the symmetry axes of the particle lattice is aligned with the @xmath8 axis as in fig . [ fig5](a ) , while in fig . [ fig1](a ) at @xmath58 it was aligned with the @xmath9 axis . as the substrate strength increases for @xmath80 , quadrimer arrangements of the particles form in each substrate minimum , and the orientation of neighboring quadrimers varies with @xmath29 . for @xmath78 , fig . [ fig5](c ) shows that all the quadrimers are aligned in the same direction . the corresponding @xmath48 in fig . [ fig5](d ) shows sixfold peaks at small @xmath50 from the triangular substrate , while at larger @xmath50 a fourfold void structure appears due to the square ordering of each quadrimer . at @xmath81 in fig . [ fig5](e ) , the individual quadrimers remain intact but their long range ferromagnetic orientational ordering is lost , producing a ring structure in @xmath48 as shown in fig . [ fig5](f ) . we have tried several different methods of preparing the initial configuration for @xmath80 and @xmath82 but have not found a ferromagnetically ordered state . the orientational ordering of dimer and trimer states was previously shown to arise due to quadrupole or higher pole moment interactions between neighboring @xmath2-mers . for quadrimers , it is likely that the higher pole moments play a more important role . if the pole moment favors ferromagnetic alignment , than as @xmath29 increases the size of the quadrimer decreases , reducing the pole moment responsible for the orientational ordering . previous studies found that as the substrate strength increases , the temperature at which dimers and trimers lose their orientational ordering drops due to the decreasing effective pole moment . we note that the configurations in fig . [ fig5 ] are all obtained at @xmath83 . when we anneal the system from a high temperature and slowly decrease the temperature to zero , we still obtain the same states illustrated in fig . [ fig5 ] , suggesting that the energy of the ferromagnetic state must be only slightly smaller than that of the rotationally disordered state . as @xmath29 increases , the picture of rigid quadrimers begins to break down , as shown in fig . [ fig5](e ) where the quadrimers become slightly distorted . this could further decrease the contributions of different pole moments to the orientational ordering or may even lead to competing interactions between neighboring quadrimers , resulting in a ground state with ordering produced by very long range interactions , as found in certain spin ice systems . for @xmath84 we find states similar to those shown in fig . [ fig5](e , f ) with increasing distortion in the square ordering , until for high enough @xmath29 the distortion becomes strong enough to permit some particles to sit at the center of a substrate minimum . the dynamic phases for @xmath80 are similar to those illustrated in fig . [ fig2 ] for the @xmath58 system , with elastic depinning for low @xmath29 crossing over to plastic depinning at higher @xmath29 , followed by a dynamical reordering transition at higher @xmath32 .     for samples with @xmath85 . ( a , b ) at @xmath86 the particles have triangular ordering that is not commensurate with the substrate , as indicated by the smearing of the peaks in @xmath48 . ( c , d ) at @xmath87 the system is disordered . ( e , f ) at @xmath88 the system forms what we term a jack state with one particle located at the center of each potential minimum surrounded by four outer particles in a square configuration . in large regions the squares form a herringbone structure , as indicated by the dashed lines ; however , some localized disorder is present in the sample . ( g , h ) at @xmath82 , an aligned pentagon state with long range orientational order forms . for @xmath84 the same structure persists but the pentagons shrink in size . ]    for @xmath89 at @xmath33 , as a function of @xmath29 we find disordered phases interspersed among ordered phases . this is also the first filling at which the @xmath2-mer assumption breaks down and a shell structure of the particles in each substrate minimum begins to form , in contrast to the case of isolated trapped clusters of repulsively interacting particles , which first develop equilibrium shell structure at @xmath90 @xcite but have a metastable state with one particle in the center for @xmath91 @xcite . in fig . [ fig6](a ) for @xmath92 at @xmath91 , a triangular lattice containing a small tilt distortion forms . the smearing of the corresponding @xmath48 in fig . [ fig6](b ) indicates that for weak substrates at this filling , the particle lattice is not commensurate with the substrate . at @xmath87 , shown in fig . [ fig1](c , d ) , the lattice is disordered as indicated by the ring pattern in @xmath48 . when we vary the initialization protocol , we always observe a disordered lattice for @xmath93 . for @xmath94 the system orders again as shown in fig . [ fig6](e , f ) for @xmath95 . each substrate minimum contains what we term a jack pattern consisting of one particle located at the center of the minimum surrounded by four particles in a square arrangement . as indicated by the dashed lines in fig . [ fig6](e ) , the jacks are tilted at @xmath96 or @xmath97 to the @xmath8 axis in every other row , producing a herringbone structure . the ordering is not complete as there are local regions where the jack structure breaks down , producing some smearing in @xmath48 as plotted in fig . [ fig6](f ) . the structure can be viewed as the beginning of a shell ordering where , in each substrate minimum , the first shell contains a single particle and the next shell contains four particles . for @xmath98 the interactions between particles in neighboring substrate minima are reduced and the pentagon ground state structure expected for particles in an isolated trap emerges . the pentagons exhibit ferromagnetic order and are all aligned in the same direction , as illustrated in fig . [ fig6](g , h ) for @xmath82 . for increasing @xmath29 the pentagons become smaller ; however , we observe no further change in the structure .     vs force @xmath32 curves for samples with @xmath89 . ( a ) at @xmath99 we observe a pinned phase ( p ) , a random phase ( r ) , and a moving crystal phase ( mc ) . ( b ) at @xmath88 the system depins into a moving channel phase ( c ) illustrated in fig .  [ fig8](a ) . this is followed by the r phase and the moving smectic phase ( ms ) . ( c ) at @xmath100 the c phase is absent and the system transitions directly from the r phase to the ms phase . ] ( a ) the moving channel ( c ) phase at @xmath88 and @xmath101 . ( b ) the moving smectic ( ms ) phase at @xmath100 and @xmath102 . ] we can relate the different pinned structures at @xmath91 to features in the velocity force curves and to dynamic phases . for @xmath103 there is an elastic depinning into a moving triangular lattice . for @xmath104 the disordered phase depins into a fluctuating random phase similar to that observed for @xmath43 , and at higher drives the system can reorder into a moving crystal ( mc ) phase . in fig . [ fig7](a ) we plot @xmath61 versus @xmath52 for a sample with @xmath99 . the random flow phase is characterized by large fluctuations , and is followed by a sharp transition near @xmath105 into a mc phase with small fluctuations . the transition between these two phases can also be detected by analyzing the noise fluctuations of the velocity response . in the random phase the fluctuations are characterized by a broad band noise signal with @xmath106 characteristics where @xmath107 to @xmath108 , while in the mc phase there is a narrow band signal with a characteristic frequency that increases with increasing @xmath52 , similar to what has been observed for dynamically reordered phases of superconducting vortices moving over periodic substrates @xcite . for @xmath109 we find no dynamical reordering and the particle arrangement remains disordered up to the highest drives we considered ; however , it is possible that a transition into a moving ordered state could occur at very high drives . for @xmath110 , where the pinned system is in the ordered jack state , the depinning is elastic and the particles depin into the ordered winding channel flow phase ( c ) illustrated in fig .  [ fig8](a ) . as @xmath32 increases , there is a crossover to a random phase followed by another transition at high drives into a moving smectic phase ( ms ) of the type shown in fig . [ fig8](b ) . the velocity signatures associated with these transitions appear in fig . [ fig7](b ) for a sample with @xmath111 . the c phase is associated with small velocity fluctuations and is followed first by the strongly fluctuating r phase and then by the ms phase at higher @xmath32 . for @xmath112 the c phase is lost and the sample depins from the ordered pentagon state into a random flow state , as illustrated in fig .  [ fig7](c ) for a sample with @xmath113 . this is followed by the formation of a moving smectic state at higher @xmath32 . as @xmath29 increases , the extent of both the pinned region and the random fluctuating state increases .     vs @xmath29 for a sample with @xmath89 . regions where the p , c , mc , r , and ms phases occur are marked . ( b ) the depinning curve @xmath114 from panel ( a ) plotted on a log - log scale . dashed lines mark the transitions between the pinned triangular crystal ( pc ) , pinned disordered or pinned glass ( pg ) , pinned jack ( pj ) , and pinned pentagon ( pp ) states . ( c ) the critical depinning force @xmath115 vs @xmath29 for @xmath116 , 3 , 4 , and @xmath40 , from top to bottom . upper dashed line : a linear fit @xmath117 . lower dashed line : a fit to @xmath44 . these curves indicate that , in general , the system exhibits collective depinning behavior . ]    by conducting a series of simulations , we map the dynamic phase diagram shown in fig . [ fig9](a ) . the disordered r phase separates the moving crystal and moving smectic phases . there is also a domelike region where the winding channel c phase occurs . in fig . [ fig9](b ) we plot the depinning line on a log - log scale with dashed lines marking the regions where the pinned phases illustrated in fig .  [ fig6 ] occur . for @xmath118 , the depinning force @xmath115 monotonically increases with increasing @xmath29 . there is a small decrease in @xmath114 at the onset of the pinned disordered or pinned glass ( pg ) phase , followed by a plateau in @xmath114 . near @xmath51 , @xmath115 begins to increase rapidly with increasing @xmath29 in the pinned jack state , while in the pinned pentagon state @xmath114 still increases with increasing @xmath29 but with a reduced slope . the behavior of @xmath114 near the onset of the pinned disordered phase is somewhat unusual since it indicates that even though the substrate strength is increasing , the depinning force does not increase . when the system is disordered , there are numerous dislocations present in the particle structure that produce weak spots which flow first at the depinning transition . in contrast , in the depinning of a crystalline state , there are no weak spots since the lattice structure is the same everywhere . this behavior is the opposite of the so - called peak effect phenomenon found for vortices in type - ii superconductors , where a sudden increase in @xmath115 occurs when the vortex lattice becomes disordered . the peak effect is generally believed to arise due to the softening of the vortex lattice when dislocations are present . in the softer lattice , the vortices can more freely move to adjust to a random pinning potential and maximize the pinning force . in our system , the pinning potential is not random , so the crystalline states are more strongly pinned than the disordered states . this has been demonstrated clearly in studies of commensurate - incommensurate transitions for both vortices and colloids in periodic substrates , where crystalline states form at integer values of @xmath0 . near but not at commensuration , numerous defects appear in the commensurate lattice , and these defects cause a reduction in @xmath114 . as a result , @xmath114 passes through a series of peaks at integer values of @xmath0 as the filling fraction is varied . in the system we are considering here , the @xmath89 state is unusual since even through the system is at an integer filling , we find a regime where the pinned state is noncrystalline . this is in contrast to the lower integer fillings @xmath7 , 2 , 3 , and 4 , where all the pinned states are ordered . this disordering effect and the plateau or drop in @xmath115 at intermediate substrate strengths in the disordered pinned regions is a general feature of @xmath119 systems at integer fillings , and both features becomes more prominent for higher fillings . it may appear from fig . [ fig9](b ) that the plateau in @xmath114 occurs in the disordered region and not in the pinned jack state ; however , the plateau persists into a window of the pinned jack state since the pj state contains some defects .    in fig . [ fig9](c ) we plot @xmath115 vs @xmath29 for samples with @xmath116 , @xmath42 , @xmath120 , and @xmath40 . the dynamics and configurations of the @xmath41 system were described in detail in previous work @xcite . for @xmath116 we find that @xmath45 , as indicated by the upper dashed line . this is expected in the limit of single particle depinning and indicates that particle - particle interactions are unimportant at depinning for this filling . for the higher fillings @xmath121 , we instead find @xmath122 as indicated by the lower dashed line . this is characteristic of collective depinning transitions , where the interactions among the particles within each minimum play an important role in the depinning process . in the limit of very strong substrates , the depinning generally occurs when one particle is pushed near the saddle point of the substrate by the other particles in the potential minimum . this single particle suddenly escapes and triggers a cascade of escape events in the rest of the system . at intermediate substrate strengths , the higher order fillings all generally show a change in the slope of @xmath114 versus @xmath29 associated with a change in the @xmath2-mer structure . versus @xmath29 . ( a ) at @xmath124 , a sharp increase in @xmath115 occurs at the transition from a pinned triangular lattice to a disordered state . dashed lines indicate regions where different pinned phases appear . triangular : triangular lattice ; disordered : disordered state ; pentagon / triangular : a monomer at the center of each well surrounded by five particles arranged in a pentagon ; hexagon : a hexagonal ring in each substrate minimum . ( b ) at @xmath125 , the triangular to disordered transition is associated with a small dip and a plateau in @xmath115 with increasing @xmath29 . dashed lines indicate the regions where different pinned states appear . triangular : triangular lattice ; disordered : disordered state ; pentagon / dimer : a composite pentagon - dimer arrangement in each minimum ; hexagon / monomer : a composite hexagon - monomer arrangement in each minimum ; ring : a single ring of particles in each minimum found at large @xmath29 . ]    . in panels ( c , d , e ) , dashed lines indicate the orientation of the pentagons in selected traps . ( a ) at @xmath46 and @xmath33 a triangular state forms . ( b ) at @xmath126 and @xmath33 a disordered state occurs . ( c ) at @xmath127 and @xmath33 there is a pentagon / monomer state with no orientational order . ( d ) at @xmath128 and @xmath33 we find an aligned pentagon / monomer state . ( e ) at @xmath128 and @xmath129 , just below the depinning transition , the pentagon / monomer state is polarized by the applied drive . ( f ) at @xmath130 and @xmath33 , a hexagonal ring state with orientational order forms . dashed lines indicate the two hexagon alignments in representative traps . ( g ) @xmath48 for the disordered phase in panel ( b ) . ( h ) @xmath48 for the hexagonal ring state in panel ( f ) . ]    in fig . [ fig10](a ) we plot @xmath115 versus @xmath29 for a sample with @xmath124 . a large jump in @xmath115 occurs at @xmath131 at the transition from a pinned triangular lattice , illustrated in fig . [ fig11](a ) , to a disordered pinned state , illustrated in fig . [ fig11](b ) . for @xmath132 , the system becomes disordered , as shown in fig . [ fig11](b ) for @xmath133 , and has a ringlike @xmath48 signature , shown in fig . [ fig11](g ) . in the disordered region the system depins into a plastic random flow state , while for @xmath134 the system depins elastically into a mc state . at this filling , @xmath114 is depressed in the triangular state since it is not possible for all the particles to simultaneously be in a triangular lattice and sit at substrate minima locations , so that some particles are instead located on substrate maxima . when the substrate strength increases enough , the particles are forced off of the substrate maxima and the system disorders , allowing more particles to be located closer to substrate minima . for @xmath135 , a shell structure emerges . the system forms a pentagonal / triangular lattice composite in which one particle sits at each substrate minimum and is surrounded by five particles in a pentagon arrangement , as illustrated in fig . [ fig11](c ) for @xmath127 . in this case the ordering is not ferromagnetic and the tilt of the pentagons varies from site to site ; however , for @xmath136 the pentagon - monomer states become aligned , as shown in fig . [ fig11](d ) for @xmath137 . an applied drive can induce a polarization of the pentagons in the driving direction , as illustrated in fig . [ fig11](e ) for @xmath128 under an @xmath8 direction drive just before depinning . for @xmath138 a ring state forms in which the center of the substrate minimum is empty and the particles form a hexagon . the rings are not all aligned ; instead , two orientations coexist in the sample as shown in fig . [ fig11](f ) . this produces twelve - fold modulations in the structure factor as illustrated in fig . [ fig11](h ) . as @xmath29 increases further , the ring structure persists but decreases in diameter . in general we find ring structures rather than shell structures for @xmath139 at large @xmath29 . the different phases are highlighted in fig . [ fig10](a ) .    . ( a ) the disordered state at @xmath54 . ( b)at @xmath128 a pentagon / dimer composite appears in which the dimers are partially aligned . dashed lines indicate two different pentagon / dimer orientations . ( c ) at @xmath140 a hexagon / monomer composite forms . dashed lines indicate two different hexagon orientations . ( d ) at @xmath141 a ring structure appears with seven particles per ring . ] figure  [ fig10](b ) shows that a similar set of pinned phases forms in a sample with @xmath125 . there is still a triangular lattice for small @xmath29 ; however , at @xmath142 when the system enters the disordered phase , @xmath114 shows a small decrease followed by a plateau up to the end of the disordered phase at @xmath143 . for @xmath144 the triangular lattice is more commensurate with the substrate and thus better pinned than the triangular lattice at low @xmath29 in the @xmath90 system . in the disordered phase , illustrated in fig . [ fig12](a ) for @xmath65 , @xmath48 contains a ringlike structure similar to that in fig . [ fig11](g ) . for @xmath145 the system forms a shell structure with an outer pentagon and two dimerized inner particles , as shown in fig . [ fig12](b ) for @xmath128 . the pentagon / dimer structures are only locally aligned , leading to smearing in the corresponding @xmath48 . figure  [ fig12](c ) illustrates the configuration at @xmath140 where an array of hexagons appears with a monomer at the center of each hexagon . the hexagons have one of two possible orientations , similar to the @xmath90 hexagon state shown in fig . [ fig11](f ) . for @xmath146 shown in fig . [ fig12](d ) , the shell ordering is replaced by a ring structure where each ring has seven particles and the ring radius decreases with increasing @xmath29 . as also found for the @xmath124 case , it is possible for an external drive to orient the structures close to depinning . the dynamical phases for @xmath147 and @xmath144 are very similar to those observed at @xmath89 . the system depins elastically at small @xmath29 in the triangular ordering regime and depins plastically in the disordered pinned regime . a disordered flow regime extends to divergingly large drives for regions of @xmath29 where a disordered pinned phase is present . for values of @xmath29 above the disordered pinned phase , the strongly driven system dynamically orders into a moving smectic phase , and the driving force at which this transition occurs increases with increasing @xmath29 . for the higher fillings of @xmath90 and @xmath144 , there are some additional features in the @xmath148 curves at high @xmath29 that did not appear at @xmath91 . for example , in the fluctuating phase , the system shows a coexistence regime with some particles that remain pinned while other particles move along 1d channels or troughs without rearranging the pinned particles . this feature becomes more noticeable at higher @xmath0 .    . in ( a , b ) , dashed circles indicate the centers of the substrate minima . ( a ) at @xmath55 a triangular lattice forms that is rotated by @xmath149 with respect to the underlying substrate . the dashed line indicates the direction along which 8 particles fit between consecutive substrate minima . ( b ) at @xmath95 , dimerization of a subset of the particles occurs and the dimers are aligned approximately @xmath150 from the @xmath8 axis . the dashed line indicates the lattice alignment direction , and a representative hexagon surrounding one dimer is shown . ( c ) at @xmath140 , seven particles surround a single monomer in each minimum . ( d ) at @xmath151 a ring state appears with eight particles per ring . ]    for @xmath152 there are several new features in the dynamic phases due to the fact that there is no intermediate disordered pinned phase . this results from a matching of the registry of the @xmath153 lattice with the underlying substrate . at @xmath75 in fig . [ fig13](a ) , the particles form a triangular lattice that is slightly anisotropic and rotated by @xmath149 with respect to the underlying substrate . this orientation permits as many particles as possible to sit close to substrate minima , and as indicated by the dashed line in fig . [ fig13](a ) , along one particle lattice symmetry direction , there are eight particles filling the space between consecutive substrate minima . this arrangement differs from the disordered lattice found at the eighth matching field in a sample with a muffin - tin potential @xcite . as @xmath29 increases , each substrate minimum traps two particles that dimerize and are surrounded by six other particles , as illustrated in fig . [ fig13](b ) . the dimers are rotationally ordered , and the overall particle lattice is rotated by @xmath154 with respect to the substrate , as indicated by the dashed line in fig .  [ fig13](b ) . orientationally ordered dimers did not form for @xmath91 , 6 , or 7 . for @xmath155 the dimer state is lost and replaced by a monomer at each substrate minimum surrounded by an outer shell of seven particles , as shown in fig . [ fig13](c ) for @xmath156 , and for @xmath157 a ring state with eight particles per ring occurs , as illustrated in fig . [ fig13](d ) for @xmath158 . vs @xmath29 for a sample with @xmath153 . the pinned ( p ) phases include : a triangular lattice ( triangular ) , a triangular lattice containing dimers ( triangular / dimer ) , heptagons surrounding monomers ( heptagon / monomer ) , and a ring state ( ring ) . there is a large increase in @xmath115 at the onset of the dimerized state . the transition flow region is labeled tr . the mc , ms , and r flow phases are also marked . ] ( upper curve ) and @xmath159 ( lower curve ) vs @xmath52 for a sample with @xmath152 at @xmath95 . p : pinned phase ; tr : transverse flow phase ; r : random flow phase with motion oriented along the @xmath8 direction only . there is a sharp transition between the tr and r flow phases . inset : particle positions ( dots ) and trajectories ( lines ) over a fixed time interval for the same sample in the transverse flow phase just above depinning at @xmath160 . dashed circles indicate the centers of the substrate minima . the particles do not flow with the applied drive along the @xmath8-direction but instead travel at @xmath161 to the @xmath8 axis along one of the principal axes of the particle lattice . ] the lack of disorder for @xmath153 is clearly shown in the dynamical phase diagram in fig . [ fig14 ] . there are no disordered pinned phases , and the random flow phase ( r ) does not extend out to large @xmath32 but is instead bounded on the high @xmath32 side by the ms flow state . for @xmath162 , the triangular lattice depins elastically in the direction of drive into an mc state . there is a sharp increase in @xmath114 at the onset of the dimerized state , and the mc flow is replaced by ms flow at higher drives . for @xmath163 , the system depins transverse to the driving direction along one of the symmetry directions of the particle lattice at @xmath161 from the @xmath8 axis , as illustrated in the inset of fig . [ fig15 ] . as @xmath52 increases , there is a sharp transition out of the transverse flow regime into a random flow phase with an average velocity oriented in the driving direction , as shown in the main panel of fig . [ fig15 ] where we plot @xmath164 and @xmath165 versus @xmath32 . in the transverse flow regime for @xmath166 , @xmath165 is negative and @xmath164 is positive , and both @xmath167 and @xmath165 increase linearly with increasing drive until , at @xmath168 , @xmath165 drops to zero and @xmath164 jumps up at the transition to the random flow regime . note that the orientation of the pinned lattice is degenerate , so for samples prepared with different random seeds , the particles would flow along @xmath161 in half the cases and along @xmath169 in the other half of the cases . the transverse flow regime forms an intermediate state between the pinned and random flow regimes . flows where the particles do not move in the direction of drive but at an angle to the drive have previously been reported for colloids on triangular substrates at a filling of @xmath170 @xcite . in that study , the drive was rotated @xmath171 relative to the substrate symmetry direction from the case considered here . dimerization of all the particles in the system created an orientational degree of freedom that determined the flow direction of the colloids just above depinning . numerical studies of vortices in type - ii superconductors with kagome and honeycomb pinning arrays also produced similar flows at angles to the driving direction in cases where two vortices were located in the center of the kagome or honeycomb plaquettes , forming dimer states @xcite .     at @xmath172 , a triangular lattice forms . ( b ) for @xmath173 at @xmath174 , each substrate minimum contains a linear trimer state aligned along the @xmath8 direction surrounded by dimer states , as indicated by the dashed lines . ( c ) for @xmath173 at @xmath175 , triangular trimers sit at each substrate minimum surrounded by a hexagon of six particles . ( d ) for @xmath176 at @xmath46 , a partially ordered triangular state occurs . ( e ) for @xmath176 at @xmath128 , aligned quadrimers form at each substrate minimum surrounded by six particles . ( f ) for @xmath176 at @xmath177 , aligned triangular trimers appear that are surrounded by seven particles . ( g ) @xmath115 vs @xmath29 for @xmath173 ( upper curve ) and @xmath176 ( lower curve ) . ]    in samples with @xmath178 , triangular configurations appear at low @xmath29 , as shown in fig . [ fig16](a ) for @xmath75 . this is the same configuration found at the ninth matching field for vortices in triangular pinning arrays in ref . @xcite . as @xmath29 increases , the particles align with the @xmath8 axis in an unusual new structure illustrated in fig . [ fig16](b ) for @xmath143 . the system breaks into dimer and linear trimer states . each substrate minimum contains a linear trimer flanked by two dimers , all oriented in the @xmath8 direction , while between adjacent minima there are elongated dimers oriented in the @xmath9 direction . for @xmath177 , fig . [ fig16](c ) shows that a superlattice of aligned triangular trimers appears at each substrate minimum , surrounded by an outer shell of six particles forming a hexagon . as @xmath29 increases further , we find a transition to a state with a monomer at each substrate minimum surrounded by an eight - particle shell , followed by a ring state at high @xmath29 ( not shown ) . for the parameters we considered , we never observed a state with two dimerized particles at the centers of the substrate minima surrounded by seven outer shell particles . either this state appears only for an extremely narrow range of @xmath29 , or it is simply too high in energy to form at all . in general we find that at the higher fillings , certain shell combinations are not observed . these are often , but not always , associated with incommensurate structures such as odd - even shells like the dimer - heptamer structure for @xmath178 .    for @xmath179 at @xmath46 , fig . [ fig16](d ) indicates that a partially ordered triangular lattice appears . this is a floating disordered solid that is only very weakly pinned by the substrate . at @xmath180 , the system transitions to a pinned disordered solid or pinned glass state ( not shown ) , where the last vestiges of triangular ordering are lost . for @xmath181 , the particles begin to develop an incipient quadrimer structure around each substrate minimum . the quadrimers are still very extended and the lattice is only partially localized by the substrate minima . at @xmath182 , the quadrimers become well localized within the minima and the state illustrated in fig .  [ fig16](e ) for @xmath183 emerges , with aligned quadrimers surrounded by hexagons . at @xmath184 , aligned triangular trimers sit at each substrate minimum as illustrated in fig . [ fig16](f ) . this is similar to the state shown in fig . [ fig16](c ) except each trimer is surrounded by seven particles in the outer ring instead of six . a similar set of orderings as those found for @xmath178 occurs for @xmath179 as @xmath29 increases until a ring state forms at high @xmath29 . [ fig16](g ) we plot @xmath114 versus @xmath29 for @xmath178 and @xmath179 . for @xmath173 , the triangular lattice in fig . [ fig16](a ) transitions into the aligned linear trimer state in fig . [ fig16](b ) at @xmath143 , and there is a kink in @xmath114 associated with this transition . for @xmath176 , the low-@xmath29 , partially ordered lattice illustrated in fig .  [ fig16](d ) has a very low @xmath115 since it is still floating above the substrate . a rapid increase in @xmath114 at @xmath180 is correlated with the transition to a pinned disordered solid . @xmath114 gradually increases as the quadrimer structure gradually organizes , until at @xmath185 , when the quadrimers become well localized within the substrate maxima into the state shown in fig . [ fig16](e ) , @xmath114 kinks upward again . near @xmath186 there is another upward kink in @xmath114 at the formation of the ring state , which depin as if they have no orientational degree of freedom . for higher order fillings @xmath187 , we generally observe the same features described for the @xmath188 fillings , so we believe that we have identified the generic features that arise in this type of system . for higher fillings it is likely that instead of only two shells of particles , three or more particle shell structures will form . the orientational ordering of particles trapped in adjacent substrate minima is also likely to become more fragile as the number of shells increases and the outer shell becomes more and more circular . thermal effects on multiple shell structures are beyond the scope of this work ; however , we expect that the shell structures would have multiple disordering transitions as a function of temperature that are distinct from the transitions observed at lower fillings in colloidal molecular crystals @xcite . another factor that could influence the ordering is the shape of the substrate minima . here we considered a sinusoidal substrate ; however , with optical trapping , other forms of substrates could be created which may favor or disfavor certain shell structures from forming .     at @xmath189 , each substrate minimum contains an inner shell of five particles and an outer shell of six particles . ( b ) for @xmath190 at @xmath46 , a triangular lattice appears . ( c ) for @xmath190 at @xmath191 , there are four particles in the inner shells surrounded by eight particles . ( d ) for @xmath190 at @xmath192 , there are five particles in the inner shells surrounded by seven particles . ]    in fig . [ fig17 ] we show some representative features of higher filling states . at @xmath193 and @xmath194 in fig . [ fig17](a ) , an aligned pentagonal inner ring forms at each substrate minimum with an outer hexagonal ring . for @xmath195 , a triangular lattice appears for @xmath49 as illustrated in fig . [ fig17](b ) , an oriented structure with four particles in the inner shell and nine in the outer shell forms at @xmath196 as shown in fig . [ fig17](c ) , and at @xmath140 there are five particles in the inner shell and seven in the outer shell as seen in fig . [ fig17](d ) . we have examined the static configurations and driven dynamics of particles interacting with a repulsive yukawa potential , such as colloids , in the presence of a triangular substrate for varied integer fillings of up to twelve particles per substrate minimum . under static conditions , we observe a rich variety of crystalline structures for fillings that have not been previously reported for this model . we show that for fillings @xmath197 , shell structures can form in the substrate minima with well - defined numbers of particles in the inner and outer shells . these shell clusters can exhibit orientational ordering across the sample . as we vary the substrate strength for a fixed filling level , we find that a series of different particle lattices occur with different numbers of particles in the shells and different amounts of orientational ordering . in the limit of strong substrates we observe a ring lattice . several fillings exhibit disordered particle configurations ; however , these fillings can show a reentrant ordering behavior as the substrate strength is increased . under an applied drive , the different particle orderings alter the critical depinning force , producing either a strong increase in the depinning force or a plateau and even decrease in the depinning force at structural transitions induced by increasing the substrate strength . we also find a remarkably rich variety of dynamical phases that can be correlated with the structure of the static pinned configuration . as the drive increases , we find moving crystalline states and moving smectic states as well as disordered flow states . our results are relevant to colloids on periodic substrates , vortices in nanostructured superconductors , and other commensurate - incommensurate systems , including aspects of friction and nonequilibrium transport . this work was carried out under the auspices of the nnsa of the u.s . doe at lanl under contract no . de - ac52 - 06na25396 . tkachenko , phys . lett .  * 89 * , 148303 ( 2002 ) ; p.l . biancaniello , a.j . kim , and j.c . crocker , _ ibid _ * 94 * , 058302 ( 2005 ) ; m.e . leunissen , r. dreyfus , f.c . cheong , d.g . grier , r. sha , n.c . seeman , and p.m. chaikin , nature mater . * 8 * , 590 ( 2009 ) . g.  malescio and g.  pellicane , nature mater . * 2 * , 97 ( 2003 ) ; m.  rechtsman , f.  stillinger , and s.  torquato , phys . e * 73 * , 011406 ( 2006 ) ; m.a . glaser , g.m . grason , r.d . kamien , a.  kosmrlj , c.d . santangelo , and p.  ziherl , epl * 78 * , 46004 ( 2007 ) . korda , m.b . taylor , and d.g . grier , phys . lett . * 89 * , 128301 ( 2002 ) ; m.p . macdonald , g.c . spalding , and k. dholakia , nature ( london ) * 426 * , 421 ( 2003 ) ; a.m. lacasta , j.m . sancho , a.h . romero , and k. lindenberg , phys . * 94 * , 160601 ( 2005 ) ; m. balvin , e. sohn , t. iracki , g. drazer , and j. frechette , phys .  rev . lett . * 103 * , 078301 ( 2009 ) .            c.  reichhardt and n.  grnbech - jensen , phys . lett .  * 85 * , 2372 ( 2000 ) ; c. reichhardt and c.j . olson reichhardt , phys . b * 76 * , 064523 ( 2007 ) ; g.r . berdiyorov , m.v . milosevic , and f.m . peeters , new j.  phys .  * 11 * , 013025 ( 2009 ) . hermann , c. rohr , m. balb ' as gambra , a. malecki , m.s . malarek , e. frey , and t. franosch , phys . b * 82 * , 165451 ( 2010 ) ; c. rohr , m.b . gambra , k. gruber , c. hohl , m.s . malarek , l.j . scherer , e.c . constable , t. franosch , and b.a . hermann , chem . commun . * 47 * , 1800 ( 2011 ) . grigorieva , w.  escoffier , j.  richardson , l.y . vinnikov , s.  dubonos , and v.  oboznov , phys . lett .  * 96 * , 077005 ( 2006 ) ; v.r . misko , b.  xu , and f.m . peeters , phys . b * 76 * , 024516 ( 2007 ) ; n. kokubo , s. okayasu , a. kanda , and b. shinozaki , phys . rev . b * 82 * , 014501 ( 2010 ) . zhao , v.r . misko , f.m . peeters , v.  oboznov , s.v . dubonos , and i.v . grigorieva , phys . b * 78 * , 104517 ( 2008 ) ; m.m . doria , a.r . de c. romaguera , and f.m . peeters , phys . b * 81 * , 104529 ( 2010 ) ; b. xu , m.v . milosevic , s - h . lin , f.m . peeters , and b. jank ' o , phys . lett * 107 * , 057002 ( 2011 ) .
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 .