# Moving Vortex Phases, Dynamical Symmetry Breaking, and Jamming for Vortices in Honeycomb Pinning Arrays

## Abstract

We show using numerical simulations that vortices in honeycomb pinning arrays can exhibit a remarkable variety of dynamical phases that are distinct from those found for triangular and square pinning arrays. In the honeycomb arrays, it is possible for the interstitial vortices to form dimer or higher -mer states which have an additional orientational degree of freedom that can lead to the formation of vortex molecular crystals. For filling fractions where dimer states appear, a novel dynamical symmetry breaking can occur when the dimers flow in one of two possible alignment directions. This leads to transport in the direction transverse to the applied drive. We show that dimerization produces distinct types of moving phases which depend on the direction of the driving force with respect to the pinning lattice symmetry. When the dimers are driven along certain directions, a reorientation of the dimers can produce a jamming phenomenon which results in a strong enhancement in the critical depinning force. The jamming can also cause unusual effects such as an increase in the critical depinning force when the size of the pinning sites is reduced.

###### pacs:

74.25.Qt## I Introduction

Vortex matter in type-II superconductors has been extensively studied as a unique system of many interacting particles in which nonequilibrium phase transitions can be accessed readily (1); (2); (3); (4); (5); (6); (7); (8). In the absence of driving or quenched disorder, the vortex-vortex interactions favor a triangular crystalline ordering. If the sample contains sufficiently strong quenched disorder in the form of randomly placed pinning sites, the vortex lattice ordering can be lost as the vortices adjust their positions to accommodate to the pinning landscape (1); (2); (3). Under an applied drive such as the Lorentz force from a current, the vortices remain immobile or pinned for low driving forces; however, there is a threshold applied force above which the vortices begin to move over the quenched disorder. For strong disorder, the initial moving state is highly inhomogeneous with the vortices flowing in meandering and fluctuating channels, and there is a coexistence between pinned vortices and flowing vortices (1); (2). At higher drives the vortices move more rapidly, the effectiveness of the quenched disorder is reduced, and the fluctuations experienced by the vortices become anisotropic due to the directionality of the external drive (4). The vortex-vortex interactions become more important at the higher drives when the quenched disorder becomes ineffective, and a dynamical transition can occur into a moving smectic state where the vortices regain partial order in one direction (5); (6); (7). Here, the system has crystalline order in the direction transverse to the vortex motion and liquid-like order in the direction of vortex motion. Depending on the dimensionality and the strength of the pinning, it is also possible for the vortices to reorganize in both directions at high drives to form a moving anisotropic crystal (4); (5); (6); (7); (8); (9). The existence of these different phases and transitions between the phases can be inferred from signatures in transport (3) and noise fluctuations (10); (11), and the moving phases have also been imaged directly using various techniques (7); (8).

In addition to the naturally occurring randomly placed pinning sites, it is also possible for artificial pinning sites to be created in a periodic structure (12). Recent advances in nanostructuring permit the creation of a wide variety of periodic pinning landscapes where the periodicity, shape, size and density of the pinning sites can be well controlled. Distinct types of pinning arrays such as square (13); (14); (15); (16); (17); (18); (19); (20) triangular (21); (22), rectangular (23); (24), honeycomb (25); (26), kagomé (25), quasicrystalline (27), and partially ordered (28) structures have been created. In these arrays the type of vortex structure that forms is determined by whether the vortex lattice is commensurate with the underlying pinning array. Commensurate arrangements appear at integer multiples of the matching field , which is the magnetic field at which the vortex density matches the pinning density, and in general, ordered vortex states occur at matching or rational fractional values of (13); (14); (15); (29); (30); (31); (32). In samples where only one vortex can be captured by each pinning site, the vortices that appear above the first matching field sit in the interstitial regions between the pinning sites, and these interstitial vortices can adopt a variety of crystalline configurations (14); (15); (16); (17); (19); (24); (29); (30); (32).

Since a number of distinct ordered and partially ordered vortex states can be created in periodic pinning arrays, a much richer variety of dynamical vortex behaviors occur for periodic pinning than for random pinning arrays (15); (16); (17); (33); (34); (35); (36); (37); (38); (39); (40); (41); (42); (43); (44); (45); (46); (47); (48); (49). Several of the dynamical phases occur due to the existence of highly mobile interstitial vortices which channel between the pinned vortices (15); (17); (35); (36); (37); (39); (40); (42); (43); (45); (49). As a function of applied drive, various types of moving phases occur, including interstitial vortices moving coherently between the pinning sites in one-dimensional paths (15); (17); (33); (34); (36); (35); (39); (43) or periodically modulated winding paths (33); (35); (44); (45), disordered regimes where the vortex motion is liquidlike (33); (35); (39); (43), and regimes where vortices flow along the pinning rows (33); (38); (46); (47); (48). Other dynamical effects, such as rectification of mixtures of pinned and interstitial vortices, can be realized when the periodic pinning arrays are asymmetric (49).

Most of the studies of vortex ordering and dynamics in periodic pinning arrays have been performed for square and triangular arrays. Experiments with honeycomb and kagomé pinning arrays revealed interesting anomalies in the critical current at nonmatching fields which are as pronounced as the anomalies observed at matching fields in triangular pinning arrays (25); (26). A honeycomb pinning array is constructed by removing every third pinning site from a triangular pinning array, producing a periodic arrangement of triangular interstitial sites. In Figs. 1(a,b) we illustrate a triangular pinning array along with the honeycomb pinning array that results after the removal of one third of the pinning sites. The matching anomalies in the experiments coincide with fields , with an integer. At these fields, the vortex density would match with the regular triangular pinning array. At the matching anomalies for , a portion of the vortices are located in the large interstitial regions of the honeycomb lattice (26), as illustrated in Fig. 1(c) for . The overall vortex lattice structure is triangular and a strong peak in the depinning force occurs at this field (50).

Recently, we used numerical simulations to demonstrate that vortices in honeycomb pinning arrays have a rich equilibrium phase diagram as a function of vortex density (50), with matching anomalies at integer and half integer matching fields that are in agreement with experiments. The large interstitial sites created by the missing pinning sites can capture multiple interstitial vortices which form cluster states of vortices. For , dimer states with form, while for higher fields trimer and higher order -mer states form. At the integer and half-integer matching fields, the -mer states can assume a global orientational ordering which may be of ferromagnetic or antiferromagnetic type; herringbone structures can also form, similar to those observed for colloidal particles on periodic substrates (51); (52); (53); (54); (55) and molecules on atomic substrates (56). These orientationally ordered states have been termed vortex molecular crystals. Certain vortex molecular crystals have ground states that are doubly or triply degenerate, such as the dimer state illustrated in Fig. 1(d) at where the dimers align in one of three equivalent directions (50). As the temperature is increased, the -mers undergo a transition from an ordered state to an orientationally disordered state in which the -mers rotate randomly but remain confined to the interstitial pinning sites. The rotating states have been termed vortex plastic crystals. At matching fields where vortex plastic crystals form, the anomalies in the critical current disappear (50). The predictions from the simulations are in general agreement with the experimental observation of the loss of certain higher order matching anomalies at higher temperatures (26). The formation of -mers that can be aligned along degenerate symmetry directions has also been predicted for kagomé pinning arrays where every other pinning site is removed from every other row of a triangular lattice (50); (57).

The formation of dimer states in the honeycomb pinning array leads to a variety of novel dynamical phases, including a spontaneous dynamical symmetry breaking effect in which the moving vortices organize into one of two equivalent states which have a component of translation perpendicular to the applied drive in either the positive or negative direction (58). The transverse response appears when the external driving force is applied halfway between the two symmetric directions of aligned dimer motion. The dynamical symmetry breaking occurs when the equilibrium ground states have no global symmetry breaking. At , the ground state is symmetry broken and the dynamical moving state has the same broken symmetry as the ground state. For incommensurate fillings, when the dimer alignment is disrupted, there is no global symmetry breaking in the ground state, and instead a dynamical symmetry breaking occurs due to the applied drive.

In this work we map the dynamical phase diagram for vortices in honeycomb arrays. We focus on the states to understand where dynamical symmetry breaking occurs and to examine what other types of moving phases are possible. We study how the dynamical phases change for driving along different axes of the pinning lattice. We find that very different kinds of dynamics occur when the driving direction is varied, and that the value of the depinning threshold is strongly directionally dependent. We also find that a novel jamming phenomenon can occur due to the formation of the dimer states. For certain directions of drive, the dimers are anti-aligned with the drive, causing the dimers to become blocked in the interstitial regions.

Although our results are specifically for vortices in type-II superconductors, the general features of this work should also be relevant for other interacting particle systems where a periodic substrate is present. Examples of such systems include vortices on periodic substrates in Bose-Einstein condensates (BEC), where different kinds of crystalline phases can occur which depend on the strength of the substrate (59); (60). It should be possible to observe different types of vortex flow states in BEC systems (61). Our results are also relevant for colloids on periodic substrates, where an orientational ordering of colloidal molecular crystals occurs which is very similar to that of the vortex molecular crystal states (51); (52); (53); (54); (55); (62). Other related systems include charged balls on periodic substrates (63) and models of sliding friction (64).

## Ii Simulation

We use the same simulation employed in the previous study of vortex equilibrium states in honeycomb pinning arrays (50). We consider a 2D system of dimensions and with periodic boundary conditions in the and directions. The sample contains vortices, giving a vortex density of which is proportional to the external magnetic field. In addition, there are pinning sites placed in a honeycomb arrangement with a pinning density of . The field at which the number of vortices equals the number of pinning sites is defined to be the matching field .

The dynamics of vortex located at position is governed by the following overdamped equation of motion:

(1) |

Here the damping constant is , where is the thickness of the superconducting sample, is the superconducting coherence length, is the normal state resistivity of the material, and is the elementary flux quantum. The vortex-vortex interaction force is

(2) |

where is the modified Bessel function, is the London penetration depth, , is the distance between vortex and vortex , and the unit vector . In this work all length scales are measured in units of and forces in units of . The vortex vortex interaction decreases sufficiently rapidly at large distances that a long range cutoff is placed on the interaction force at to permit more efficient computation times. We have found that the cutoff does not affect the results for the fields and forces we consider here.

The pinning force originates from individual nonoverlapping attractive parabolic traps of radius which have a maximum strength of . In this work we consider the limit where only one vortex can be captured per pinning site, with the majority of the results obtained for . The exact form of the pinning force is:

(3) |

Here, , is the location of pinning site , the unit vector , and is the Heaviside step function.

The external drive represents the Lorentz force from an applied current which is perpendicular to the driving force and is applied uniformly to all the vortices. We apply the drive at various angles to the symmetry axes of the honeycomb pinning array. The thermal force originates from Langevin kicks with the properties and . Unless otherwise noted, the thermal force is set to zero. The initial vortex configurations are obtained by simulated annealing, and the external force is then applied gradually in increments of every 1000 simulation time steps. For the range of pinning forces used in this work, we find that this force ramp rate is sufficiently slow that transients in the vortex dynamics do not affect the overall velocity-force curves. We obtain the velocity-force curves by summing the velocities in the (longitudinal) direction, , and the (transverse) direction, , where . In Fig. 1(c,d) we illustrate the pinning sites and vortex configurations after simulated annealing for [Fig. 1(c)] and [Fig. 1(d)]. Here and . In our previous work, Ref. (58), the drive was applied along the -direction for the geometry in Fig. 1.

## Iii Dynamics and Transverse Response For Driving in the Longitudinal Direction

We first consider the case for driving in the or longitudinal direction, , for the system shown in Fig. 1(d) with , , and . In Figs. 2(a,b) we plot and versus . At this filling there are four distinct dynamical phases, with the pinned (P) phase occurring at low . The depinning threshold occurs near when the interstitial vortices become depinned. For a system with random pinning and , there would be no transverse velocity response; the system would have and only would be finite. In contrast, for the honeycomb pinning array there is a finite velocity both in the positive direction and in either the or direction. In Fig. 2(b) the transverse response is negative, indicating that the vortices are moving at a negative angle to the axis for . Figure 3(a) illustrates the vortex motion at , where the vortices flow in one-dimensional paths oriented at to the axis. In Fig. 3(b) a snapshot of the vortex positions shows that the vortex lattice remains ordered in the moving phase, indicating that the vortices are flowing in a coherent manner. We term the phase shown in Fig. 3(a) the symmetry broken (SB) phase, since the flow can be tilted in either the positive or negative -direction.

At and , the interstitial vortices form an aligned dimer configuration with a three-fold degenerate ground state in which the dimers can be oriented along the -direction, as in Fig. 1(d), or along or to the -direction, as shown in previous work (58). When a driving force is applied to the or ground states, the vortices depin and flow along or , respectively. In these cases, the symmetry breaking in the moving state is not dynamical in nature but reflects the symmetry breaking within the ground state. If the dimers are initially aligned along the -direction in the ground state, an applied drive induces an instability in the pinned phase and causes the dimers to rotate into the or directions, as we discuss in further detail below. In this case the symmetry breaking is dynamical in origin.

In Fig. 2(a,b) we find pronounced oscillations in both and just above the depinning threshold . These oscillations are not intrinsic features but are due to the fact that at the interstitial vortex lattice is perfectly ordered, so the interstitial vortices move in a coherent fashion as shown in Fig. 3(a). At depinning, the interstitial vortices are slowly moving through a periodic potential created by vortices that remain trapped at the pinning sites. This periodic potential causes the moving interstitial vortices to develop an oscillating velocity. In Fig. 4(a), the instantaneous time traces of the vortex velocity and at constant show strong velocity oscillations. The oscillations are also visible in Fig. 2 at low drives due to our choice of averaging time spent at each value of the driving current. If the averaging time is increased, the oscillations in Fig. 2 disappear. We note that the locations of the boundaries between the different phases are not affected by the value of the velocity averaging time. At incommensurate fields, there is enough dispersion in the velocity of the moving interstitial vortices that the coherent velocity oscillations are no longer distinguishable.

As increases, the net vortex velocity in the SB phase increases linearly until , where there is an abrupt increase in . Fig. 2(a,b) shows that this increase coincides with a jump in to a zero average, indicating that the vortices are moving only in the -direction on average. In Fig. 3(c) we illustrate the disordered vortex trajectories that occur in this phase at . The vortices are continually depinning and being repinned, and the order in the vortex lattice is lost, as shown in Fig. 3(d). We term this the random (R) phase. It resembles random dynamical phases that have previously been observed for vortices in square pinning arrays when the interstitial vortices begin to depin vortices from the pinning sites (29). Figure 2 shows that there are pronounced random fluctuations in and in phase R, and also that does not increase linearly with but has a curvature consistent with , where is the threshold value for the SB-R transition. In the SB phase, the number of moving vortices is constant and is equal to the number of interstitial vortices, while in the R phase the number of moving vortices increases with .

At , the system organizes into a one-dimensional flowing state where the vortex motion is locked along the pinning rows, as shown in Fig. 3(e,f) for . The onset of this phase also coincides with the decrease of fluctuations in and the loss of fluctuations in , as shown in Fig. 4(c). For , all of the vortices are mobile and Fig. 2(a) illustrates that the versus curve becomes linear again. We term this the moving locked (ML) phase, since the vortex motion is effectively locked along the pinning sites. When the vortices are rapidly moving, the pinning sites have the same effect as a flashing one-dimensional trough that channels the vortices (33); (38). The vortices assume a smectic structure in the ML phase, since different rows have different numbers of vortices, resulting in the formation of aligned dislocations. The ML phase is essentially the same state found in square pinning arrays at high drives when (33).

In previous studies of square pinning arrays with strong pinning, the initial motion of the vortices for occurred in the form of one-dimensional channels between the vortices trapped at the pinning sites (33). In the honeycomb pinning array, similar flow occurs in the SB phase as shown in Fig. 3(a). For in the honeycomb array, the initial interstitial flow for depinning in the direction occurs via the flow of individual vortices in a zig-zag pattern around the pinned vortices. Since there is no dimer ordering for these fillings, no transverse response occurs for . For , the interstitial vortices begin to form dimer states when two interstitial vortices are captured in a single large interstitial site. The dimers can lower their orientational energy by aligning with each other in both the ground state and the moving states. Dimers can only remain aligned in the moving state if they are channeling along one of the symmetry axes of the pinning lattice. If the dimers were to move strictly in the -direction, they would be forced directly into the pinned vortex in the pinning site to the right of each large interstitial site. This would destabilize the rodlike dimers. Instead, the dimers maintain their integrity by moving along to the -axis. Within the moving state, if one of the dimers were to move along while the remaining dimers were moving along , the two interstitial vortices comprising the first dimer would be forced close together, destabilizing the dimer state due to the repulsive vortex-vortex interactions. Instead, all of the dimers move in the same direction.

The SB-R transition occurs when the combined forces on the pinned vortices from the external drive and the moving dimers are strong enough to depin the pinned vortices. At the closest approach in the direction between a dimer and a pinned vortex, the frontmost dimer vortex is a distance from the pinned vortex and the rear dimer vortex is a distance from the pinned vortex, where is the lattice constant of the undiluted triangular pinning lattice. In addition to the force from the dimerized vortices, the pinned vortex experiences an opposing force from the neighboring pinned vortex a distance away. In a simple approximation, the driving force needed to depin a vortex at a pinning site is thus . Setting gives , close to the value of for the SB-R transition in Fig. 2. Once the pinned vortices depin, the system enters the random (R) phase, and since is still considerably less than , it is possible for vortices to be pinned temporarily in phase R.

Studies of square pinning arrays have shown that after the onset of a random dynamical phase, the vortices can organize into a more ordered phase of solitonlike pulse motion along the pinning rows, followed by a phase in which all of the vortices channel along the pinning rows (33). At the transition to the one-dimensional pulse like motion, a larger fraction of the vortices are pinned compared to the random phase, so a drop in with increasing occurs, giving a negative differential conductivity. In the honeycomb pinning arrays for the parameters we have chosen here, we do not observe one-dimensional pulse motion or negative differential conductivity for driving along the -direction. For the one-dimensional pulse motion or the ML phase motion seen in Fig. 3(e,f) to occur, the vortices must be moving at a sufficiently high velocity for the pinning sites to act like a flashing trough. When the vortices move along the pinning rows, the vortex lattice structure adopts a highly anisotropic configuration which would be unstable at . During the period of time when a vortex passes through a pinning site, the vortex is pulled toward the center of the pinning row, which stabilizes the one-dimensional motion. When the vortex is moving between the pinning sites, it can drift away from the one-dimensional path until it encounters another pinning site. In Ref. (33), it was shown that for square pinning arrays, increasing the pinning radius stabilized the one-dimensional flow down to lower values of . In the honeycomb pinning array, the one-dimensional flow is less stable due to the fact that the vortices must move over the much wider large interstitial site, giving the vortices more time to drift away from the pinning row. Since this means that a larger value of is required to stabilize the one-dimensional motion, it should be more difficult in general to observe the onset of one-dimensional solitonlike motion or negative differential conductivity in the honeycomb pinning arrays than in the square pinning arrays.

### iii.1 Fluctuations and Noise Characteristics

In order to characterize the moving phases more quantitatively, in Fig. 4 we show time traces and of the vortex velocities at fixed for the different phases for the system in Fig. 2. In the symmetry broken (SB) phase at , shown in Fig. 4(a), is greater than by or about . Here both components of the velocity show a pronounced oscillation which arises when the interstitial vortices move in a coherent fashion over the periodic potential substrate created by the immobile vortices in the pinning sites. In Fig. 5(a), we plot the corresponding power spectrum of obtained from

(4) |

There is a pronounced peak in at the frequency of the velocity oscillation in the SB phase, indicating that mode locking effects could appear at when the symmetry breaking flow occurs. In square pinning arrays, experiments (17) and simulations (36) revealed Shapiro step-like mode locking of interstitial moving vortices at . In the honeycomb lattice, since there is also a strong oscillation in in the SB phase, we expect that transverse mode locking could occur if an additional ac drive is applied in the -direction. Such mode locking would appear as steps in both and versus in the SB phase. Transverse phase locking, which produces steps that are distinct from Shapiro steps, has been observed for the motion of vortices in square arrays (37). In general, if the vortices already have an intrinsic velocity oscillation in the transverse direction, then pronounced transverse phase locking is possible.

In Fig. 4(b) we plot the time trace of and for the random (R) phase at . In this case , and although both and show fluctuations, no oscillations or washboard frequencies appear. In Fig. 5(b) we show the corresponding for , where we find a broad band noise feature consistent with disordered plastic flow (3); (6); (11). Since there are no coherent velocity oscillations, mode locking should be absent in the random (R) phase.

In Fig. 4(c) we plot and in the ML phase at , where has been shifted down by a factor of 3 for clarity. There are no visible fluctuations in due to the one-dimensional nature of the flow, but there are small periodic oscillations in generated by the motion of the vortices over the periodic substrate. Due to the fact that different one-dimensional rows contain different numbers of vortices, producing dispersion in the vortex velocities, the oscillation in is not as pronounced as in the SB phase. The corresponding power spectrum in Fig. 5(c) contains a rich variety of peaks due to the wide range of frequencies present in this phase. The main peak is smaller in magnitude than that found for the SB phase. As increases, the frequency at which the first peak occurs also increases. It should be possible to generate phase locking in the ML phase; however, it would likely not be as pronounced as in the SB phase. These results suggest that noise fluctuations can be a useful technique for exploring the presence of different dynamical phases in periodic pinning arrays.

### iii.2 Dynamical Symmetry Breaking in the Pinned Phase

As previously noted, the ground state at is three-fold degenerate. When the dimers are aligned at either or to the -axis in the ground state, the subsequent SB flow is aligned in the same direction as the ground state. It is also possible for the dimers to align in the -direction, as shown in Fig. 6(a). At , the dimers and the vortices in the pinning sites are shifted slightly to the right due to the applied drive. As is further increased, a symmetry breaking transition occurs within the pinned phase. For the dimers remain aligned the -direction; however, at , the rotational instability illustrated in Fig. 6(b) occurs. The dimers rotate in such a way that they end up aligned in the direction. The interstitial vortex at the bottom of the dimer moves in the direction and by a smaller amount in the direction, while the vortex at the top of the dimer moves in the direction and by a smaller amount in the direction. There is also a slight shift of the vortices in the pinning sites that are closest to the bottom of each dimer. In Fig. 6(c) the rotation process is completed and the dimers are aligned in a new direction, . The vortices remain pinned until , at which point the system enters the SB phase. At finite temperatures, the dimer realignment occurs at even lower values of . The rearrangement can also be observed as a jump in and as shown in Fig. 7, where there is a positive spike in and a negative spike in near , in agreement with the motion shown in Fig. 6(b). We term this a dimer polarization effect since the driving force induces an alignment of the dimers. In runs with slightly different initial conditions, the dimers may align along the direction with 50% probability.

### iii.3 Dynamics for

We next consider the effect of changing the vortex density for fillings where interstitial dimers are present and the SB phase occurs. In Fig. 8(a) we illustrate the vortex positions for , where a mixture of monomers and dimers appear in the large interstitial sites. At this filling, the overall orientational ordering of the dimers is lost in the ground state, and the dimers are oriented only in local patches. For , a mixture of interstitial dimers and trimers is present, as shown in Fig. 8(b) for , and the orientational ordering is again lost. In Ref. (58) it was shown that the SB state still occurs at incommensurate fields as long as some dimers are present. If is suddenly increased from zero to a finite value at which only the interstitial vortices depin, the moving state for the incommensurate fields organizes into a dynamically symmetry broken state where all of the dimers flow along or . At the incommensurate fields, only the dimers undergo dynamical symmetry breaking; the monomers and trimers continue to move in the direction of the drive, with some fluctuations in the transverse direction.

In Fig. 9 we plot for the system in Fig. 2 at , , and . In Figs. 9(a,b), the same four phases described above are labeled. The SB phase has opposite sign in Fig. 9(a) and Fig. 9(b); the dynamical symmetry breaking can occur in either direction since there is no symmetry breaking in the ground state. If slightly different initial conditions are used, such as by changing the initial annealing procedure, the dynamical symmetry breaking has equal probability to occur in the positive or negative direction, as shown previously (58). In Figs. 9(a) and (b) the initial portion of the SB phase has fluctuations in due to the fact that we are increasing at a finite rate and there is a transient time for the moving state to fully organize into the SB state, as studied previously (58). The transient time increases as increases. If we decrease , the fluctuations at depinning are reduced; however, the boundary between the phases does not shift.

In Fig. 9(c) at , since there are only trimer states present. The large oscillations in occur when the system forms a completely ordered trimer ground state (50) and the vortex motion is highly coherent, similar to the effect seen in Fig. 2. For the rate at which we sample and average versus , the periodic fluctuating vortex velocity is visible. For the system enters a partially moving locked phase where a portion of the vortices move along the pinning rows. There are, however, too many vortices to form straight one-dimensional chains of the type shown in Fig. 3(e) for . A buckling instability of the chains occurs since the amount of anisotropy that would occur if one-dimensional chains formed is too large for the vortex lattice to sustain. Instead, a partially moving locked (PML) phase forms with a disordered moving vortex lattice. This result is interesting since it indicates that moving vortex phases do not always organize into ordered states. A time trace of a PML state at fixed shows much weaker velocity oscillations than those shown in Fig. 4(c) for the ML state. This suggests that phase locking with PML states will be very weak or absent. In previous work on phase locking for square arrays, it was shown that the phase locking is most pronounced at commensurate fields where the moving vortex structures are more ordered (17); (36).

By performing a series of simulations for varied , measuring the features in the velocity force curves and observing the vortex structures, we construct the dynamical phase diagram of vs shown in Fig. 10. The depinning force marking the end of the pinned (P) phase shows peaks at , 2.0, and , corresponding to the commensurate and ordered ground states reported previously (50). The SB-R transition line is fairly flat as a function of with an enhancement to higher values of occurring near , while at the incommensurate fields, monomers or trimers create fluctuations that cause the vortices at the pinning sites to depin at slightly lower values of . For , upon increasing the random state organizes into a ML state where all the vortices move along the pinning rows as shown in Fig. 3(e), while for , the random state organizes into the PML state. The width of the random phase, as determined by the fluctuations in the velocity, increases and persists to higher values of for increasing at . For and , where dimers are no longer present, the SB phase is lost and a new set of dynamical phases arises which we discuss in more detail below.

### iii.4 Effect of Changing the Pinning Strength

We next consider the effect of changing the pinning strength when . The four phases in Fig. 2 occurred in a sample with . As is varied, we find several different kinds of ordering within the pinned phase that affect the dynamics which occur at finite . In Fig. 11 we plot the threshold depinning force as a function of . For the pinning is weak enough that the vortex-vortex interactions dominate over the pinning energy and a nearly triangular vortex lattice forms, as shown in Fig. 12(a). In this arrangement, half of the pinning sites are still occupied, so the vortex lattice is partially pinned (PP) and there is a finite depinning threshold. This type of partially pinned vortex lattice was observed in previous simulations on honeycomb pinning lattices (50), and similar partially pinned vortex lattice states have been predicted for square pinning arrays (32) and observed for metallic particles on periodic structures (63). The depinning transition from the PP state is elastic, and all the vortices depin simultaneously to form the moving triangular crystal (MC) shown in Fig. 13(a). In the MC, half of the vortices move in one-dimensional paths along the pinning rows while the remaining vortices move through the interstitial regions with a small transverse oscillation.

Figure 11 shows that the pinned ordered dimer state (P) forms for . Over the range , the depinning from state P does not occur by the initial flow of the interstitial vortices into the SB phase, unlike the case shown earlier for . Instead, for , both the interstitial vortices and the vortices at the pinning sites depin simultaneously and rearrange into the moving crystal (MC) state shown in Fig. 13(a). We also find a peak in at . This peak occurs due to both the change in the pinning configuration and a change in the depinning process. For , only half of the pinning sites are occupied and the vortex lattice depins elastically. At , all of the pinning sites are now occupied in the P state, but the vortex lattice still depins elastically. The pinning energy that must be overcome to depin the lattice is increased compared to the PP state, leading to an increase in . For , the depinning process is plastic and only the interstitial vortices flow at depinning to form the SB state. Since the plastic depinning process does not require pinned vortices to depin, the threshold force drops, producing the peak in at .

For in Fig. 11, the depinning threshold slowly increases with increasing , and a transition in the pinned vortex structure occurs at . For , the dimers in the pinned state are no longer aligned but form a pinned herringbone (PHB) type structure such as that shown in Fig. 12(b), where the dimers are tilted in the same direction in one row and tilted in the opposite direction in the adjacent rows. Herringbone ordering of dimers has previously been observed for colloidal dimers on triangular lattices (33) and for vortices in kagomé arrays at filling (50). At depinning, the PHB state does not form a SB phase but instead forms the winding interstitial phase shown in Fig. 13(b). The dimers break apart into two monomers, with one monomer passing around the pinned vortices in the positive direction and the other monomer passing the pinned vortices in the negative direction. We term this state the moving interstitial (MI) phase.

In contrast to the herringbone state, the aligned dimer or ferromagnetic ordering of the P state occurs when the pinned vortices adjacent to the dimers are able to undergo a periodic distortion within the pinning sites, reducing the interaction energy between the pinned and interstitial dimer vortices and permitting the dimer alignment. If the dimers are aligned along , as in Fig. 5(a) of Ref. (50), the two pinned vortices closest to each interstitial vortex in the dimer both splay outwards away from the direction. When is increased, the pinned vortices are pulled toward the center of each pinning site and are no longer able to distort in order to accommodate the aligned dimer state, so the pinned herringbone (PHB) state forms instead. This result suggests that there may be other types of ground state ordering for vortices in honeycomb and kagomé arrays in addition to those that have been reported previously. It may be possible to use the size and shape of the pinning sites as a means of controlling the type of crystalline structure that forms.

In Fig. 14, we plot versus for , 1.125, and 2.125. Figure 14(a) shows the elastic depinning process for the PP state at , which moves directly into the MC phase after depinning. At in Fig. 14(b), increases linearly with through the SB phase. The slope of increases in the R phase, and the velocity-force relationship becomes linear in the ML phase. For , Fig. 14(c) indicates that the depinning occurs in two steps. The first depinning transition of the interstitial vortices only takes the system from the P phase into the moving interstitial (MI) phase, while at the second depinning transition, the pinned vortices depin and the sample enters the ML phase. Unlike the behavior at in Fig. 14(b), at the intermediate random (R) phase is lost and is replaced by a sharp jump into the ML phase.

By conducting a series of simulations we construct the dynamical phase diagram as a function of and , as shown in Fig. 15. At high , the MC phase forms for , while for the ML phase appears instead. The SB phase exists for , and the SB-R boundary shifts to higher with increasing until it terminates at . For , the PHB state occurs at low drive, and the system depins into the MI phase. The MI phase also extends as far down as , where the system passes from the SB phase into a narrow window of the random (R) phase with increasing before the vortices organize into the MI phase. As continues to increase, the vortices at the pinning sites depin and the system passes through a second narrow window of the R phase until the vortices organize into the ML phase. At high , the R phase becomes vanishingly small and the system passes directly from the MI to the ML phase. The transition into the ML state increases linearly with increasing , while the depinning force saturates with increasing .

### iii.5 Changing and

We next examine the effects of changing the pinning radius in a system with fixed and . In Fig. 16 we show the dynamic phase diagram for versus obtained from a series of simulations. For , we find the same three moving phases, SB, R, and ML, as in Fig. 2. For , the pins are large enough to permit double vortex occupancy at the individual pinning sites; in this case, a new set of phases appears which we do not consider in this work. For , the R-ML transition decreases in with increasing since the larger pinning sites make it easier for the vortices to localize along the pinning rows and flow in the one-dimensional motion of the ML phase. For at low , we find the partially pinned PP state illustrated in Fig. 12(a). The onset of the PP phase coincides with a drop in at . Unlike the PP phase that occurs at low in Fig. 15, which depins elastically, the PP phase at small depins into a moving interstitial phase that is distinct from the moving interstitial phase shown in Fig. 13(b). In Fig. 17(a) we illustrate the vortex trajectories in the phase which we term the moving interstitial phase 2 (MI2). The vortices flow in winding interstitial channels; however, unlike the MI phase, in the MI2 phase only half of the pinning sites are occupied. In Fig. 17(b) we plot vs for a system with . A clear two-step depinning transition occurs, with increasing linearly with increasing in the MI2 phase. At high and , the pinning sites are too small for the ML phase to occur, and instead the vortices flow in the MC phase illustrated in Fig. 13(a).

We next consider samples with fixed , , and , but vary the value of by changing the pinning density . This alters the average spacing between neighboring vortices. Up to this point we have used . In Fig. 18 we illustrate the dynamic phase diagram for versus . As increases, increases since the depinning of the interstitial vortices is determined by the potential created by the vortices located at the pinning sites, and as the vortex density increases, the depth of the interstitial pinning potential also increases. The R-ML transition shifts to higher with decreasing . Since the distance between the pinning sites increases with decreasing , the moving vortices spend less time in the pinning sites. This destabilizes the ML phase and the vortices must move at higher velocities for the effective trough potential to be able to stabilize the ML phase. For , the SB phase is lost, since at this pinning density the interactions between the interstitial vortices and the pinned vortices become sufficiently strong that the depinning of the interstitial vortices also causes the pinned vortices to depin. As a result, the system passes directly from the P phase to the R phase with increasing . For , the vortex-vortex interaction becomes weak enough that the system depins into the MI phase illustrated in Fig. 13(b). This also coincides with an increase in the value of at which the R phase appears, since the moving interstitial vortices in the MI phase do not approach the pinned vortices as closely as they do in the SB phase.

### iii.6 Effect of Finite Temperature

We next consider the effect of finite temperature on the system in Fig. 2 with , , , and . In Ref. (50), we showed that a transition can occur at finite temperature in which the vortex -mer states lose their orientational ordering and begin to rotate while remaining confined within the large interstitial sites. This state was termed a vortex plastic crystal. In Ref. (58) we demonstrated that the SB phase disappears in the vortex plastic crystal state. In Fig. 19 we plot the dynamical phase diagram of vs for the same system in Fig. 2. In our units, the dimers melt at . Above the melting temperature, there is appreciable creep of the interstitial vortices as they hop from one large interstitial site to another, as illustrated in Fig. 20(a) for and . As is further increased the system enters the random (R) phase shown in Fig. 20(b) for . At higher drives, the vortices begin to localize along the pinning rows in one-dimensional channels; however, there is still appreciable hopping from one row to another as shown in Fig. 20(c) for . At even higher drives, the ML is recovered as illustrated in Fig. 20(d) for . For , the ML phase is lost and the high flow is in the random (R) phase shown in Fig. 20(c). In the inset to Fig. 20 we demonstrate that at , the sharp features in the velocity force curve seen at in Fig. 2 disappear.

## Iv Dynamics for Driving in the Transverse Direction

We now consider the case where is applied along the -direction, , for the same system as in Fig. 1(d) with , , and . A different set of dynamic phases appear that are distinct from those found for driving in the -direction. In particular, the SB phase is lost and the dimers align in the -direction with the initial depinning occurring in one-dimensional interstitial flow paths. In Fig. 21 we plot the velocity force curves for 2.15, and 2.32.

Figure 21(a) shows the three phases, pinned (P), one-dimensional moving interstitial (1DMI), and moving locked (ML) that occur at . In the P phase, if the ground state contains dimers which are aligned at either or to the -axis, a polarization effect is induced by the applied drive similar to the effect discussed earlier. In this case, however, the dimers shift such that they are aligned in the -direction. The one-dimensional moving interstitial (1DMI) state which appears above depinning is illustrated in Fig. 22(a) at , where the interstitial vortices move between the vortices in the pinning sites. Near there is a sharp depinning transition for the vortices in the pinning sites. After this depinning transition occurs, the vortices very rapidly rearrange into a moving locked (ML) phase where the vortices move along the pinning sites, as shown in Fig. 22(b) for . In the ML phases for driving along the and the directions, the vortices travel in one-dimensional channels along a row or a column of pinning sites, respectively. The pinning rows followed by the vortices for -direction driving are evenly spaced in the -direction, so the vortices flow through the centers of the pinning sites. In contrast, the pinning columns followed by the vortices for -direction driving are unevenly spaced in the direction due to the symmetry of the honeycomb lattice. As a result, for the -direction driving the vortices do not flow through the centers of the pinning sites, but are instead shifted to the right and left of the pinning sites in alternate columns, as seen in Fig. 22(b). This produces a more even spacing between the columns of moving vortices. Since different columns contain different numbers of vortices, the ML phase for -direction driving has smectic type characteristics. The velocity-force curves for have the same general form as the curve in Fig. 21(a) and show the same three phases.

For the appearance of trimer states disrupts the 1DMI flow since the trimers cannot align completely in the -direction. This produces random (R) vortex flow at depinning, as shown in Fig. 21(b), with diffusive vortex motion occurring along the -direction. There are more pronounced fluctuations in in the R phase, and the velocity-force curve is nonlinear and lower than the extrapolated linear behavior in the 1DMI phase that begins near . The trimers can block the one-dimensional channels of flow shown in Fig. 22(a), lowering the number of mobile vortices. At higher drives, the trimers depin, straighten into a linear configuration, and flow in the 1DMI phase. At , the vortices in the pinning sites depin, resulting in a transition from the 1DMI phase to the R phase. For sufficiently high drives, the ML phase forms. For driving along the -axis at , we showed in Fig. 9(c) that the ML phase is lost due to a buckling transition of the one-dimensional chains of vortex motions, and that a partially moving locked (PML) phase forms instead when a portion of the vortices move through the interstitial regions. For driving along the -axis at , the ML state remains stable for much higher values of than for the -axis case. A comparison of Fig. 22(b) and Fig. 3(e) shows that the interstitial region crossed by the vortices in the moving channels is not as wide for -axis driving as for the -axis driving, resulting in more stable -axis ML flow. As is further increased, the random regime grows until the 1DMI phase is completely lost, as shown in Fig. 21(c) for . In Fig. 23 we plot the dynamical phase diagram for versus , highlighting the onset of the different phases. The transition to the ML phase shifts to higher values of with increasing since the ML vortex channels become increasingly anisotropic as the number of vortices in the sample increases.

### iv.1 Dynamics as a function of and Dimer Jamming

We now consider the vortex dynamics in a system with fixed , , and for varying with . As noted above, for a partially pinned (PP) vortex lattice forms. In Fig. 24(a), the velocity-force curve for shows that the depinning of the PP phase is elastic and occurs in a single step transition to a moving crystal (MC) phase, as seen earlier in Fig. 14(a) for driving along the -axis. In the MC phase, half of the vortices move along the pinning sites. For , at low drive the system is in the pinned (P) phase of orientationally ordered dimers, and as the drive increases, Fig. 24(b) shows that the same 1DM1 and ML phases illustrated in Fig. 21(a) appear. The rapid rearrangement of the vortices from the 1DMI phase to the ML phase results in a small jump near which marks the 1DMI-ML transition. For strong pinning , Fig. 24(c) indicates that the ground state forms the pinned herringbone (PHB) phase illustrated in Fig. 12(b). The same pinned state appears for -direction driving at strong pinning, as shown in Fig. 14(c). In Fig. 24(c), the velocity-force curve at shows the abrupt nature of the depinning transition from the PHB phase to the 1DMI phase, which differs from the smoother depinning transition that occurs from the P phase to the 1DMI phase in Fig. 24(b). The depinning threshold increases markedly with increasing once the system enters the PHB state. In Fig. 25 we plot the dynamic phase diagram for versus . Near the transition from the PP to the P phase, there is a peak in similar to the peak observed at the PP-P transition for driving in the -direction in Fig. 11. For the strong enhancement of the depinning threshold in the PHB state can be seen clearly.

In Fig. 26(a) we illustrate the vortex positions just before depinning for . Even though the drive is applied in the -direction, the dimers have aligned with the -direction. When the dimers are oriented along the -axis, they cannot fit through the easy-flow one-dimensional channel between the pinning sites, but instead are essentially jammed by the two pinned vortices at the top edge of the large interstitial site. In the ordered dimer pinned (P) phase, the dimers all reorient in the same direction under an applied drive. In contrast, in the pinned herringbone (PHB) phase, the dimers rotate in opposite directions under an applied drive, so when the drive is applied along the -direction the dimers end up aligning in the -direction. In Fig. 26(b) only the vortex positions from Fig. 26(a) are shown to indicate more clearly the shift of the dimers in the positive -direction. This vortex configuration, which we term the jammed (J) state, has a structure that is distinct from that of the pinned herringbone phase shown in Fig. 12(b). The jammed state configuration exists only in the presence of the applied drive. For the dimers return to the pinned herringbone (PHB) state. In the jammed (J) state, the critical current is up to three times larger than in the state where the dimers are aligned in the -direction. We also note that at incommensurate fields for , the net vortex flow is reduced since some of the dimers align in the -direction and effectively block the motion of other vortices along the -direction.

In order to better characterize the enhancement of in the jammed state, in Fig. 27 we plot the critical depinning force in the -direction, , and in the -direction, , versus . In the inset of Fig. 27 we show the ratio versus . In the partially pinned (PP) phase, , while in the pinned aligned dimer (P) phase, is slightly higher than since the vortices can depin more readily into the 1DMI phase in the -direction. In the jammed state that forms from the PHB phase, is 3.1 times higher than for the same value of .

### iv.2 Effects of Changing and

In Fig. 28 we plot the dynamical phase diagram versus for driving in the -direction with , , and . For , the system depins into the 1DMI phase and makes a transition to the ML phase at higher drives. For , the system forms the partially pinned (PP) phase where only half of the pinning sites are occupied. The PP phase depins into a moving interstitial phase (MI2Y) that resembles the MI2 state observed for driving in the -direction in Fig. 17(b), where half the vortices depin while the other half remain pinned. The MI2Y phase is oriented from the MI2 phase. At for , the vortices at the pinning sites begin to depin and repin, giving a regime of the random (R) phase until becomes large enough for all the vortices to depin into the ML phase. For , the jammed (J) state discussed in Fig. 26 occurs due to the formation of dimers aligned in the -direction, which is associated with a marked increase in . As discussed earlier, the pinned herringbone (PHB) phase and jammed (J) state occur when becomes high enough that the vortices in the pinning sites cannot shift to allow for dimer ordering to occur. Similarly, as is reduced, the vortices in the pinning sites have less room to adjust for dimer ordering, so the PHB state forms. The jamming also produces the counterintuitive effect that as increases above , the depinning threshold decreases.

In Fig. 29 we show the dynamical phase diagram for versus , which determines the value of , for , , and . As increases, the critical depinning force into the 1DMI phase increases since the repulsion from the pinned vortices experienced by the interstitial vortices increases as the average vortex-vortex spacing decreases. The transition from the 1DMI phase to the ML phase shifts to higher values of as decreases since the distance between the pinning sites which stabilize the ML flow increases.

## V Discussion

Our results are for honeycomb pinning arrays where it was shown in previous work that -merization of the interstitial vortices into vortex molecular crystal states occurs for . Many of the dynamical effects presented in this work are due to the -merization effect. In kagomé pinning arrays, similar types of vortex molecular crystal states appear, so we expect that many of the same types of dynamic phases described here will also occur for kagomé pinning arrays, although we do expect that there will be certain differences as well. In the kagomé pinning array, the vortex dimer state appears at and has a herringbone ordering even for large, weak pins. There are no easy flow channels along to the -axis, so the symmetry breaking flows should be absent. Additionally, since there is no easy flow channel in the -direction, the anisotropic depinning dynamics may be different as well.

We have only considered in this work. At higher fields, a wide array of vortex molecular crystal states occur that should also have interesting dynamical phases. Since the low matching fields are more robust, observing the dynamics near these low fields experimentally is more feasible. Although our results are specifically for pinning sites with single vortex occupation, similar dynamics should occur if the first few matching fields have multiple vortices at the pinning sites. In this case, the effective dimerization of the interstitial vortices would be shifted to higher magnetic fields.

Although true phase transitions are associated only with equilibrium phenomena, the nonequilibrium phases considered here have many analogies to equilibrium phases. For example, several of the transitions between the nonequilibrium phases have a continuous type behavior, while in other cases the transitions are sharp, indicative of a first order nature. Future studies could explore the possible emergence of a growing correlation length near the transitions to see whether they exhibit the true power law behavior associated with continuous phase transitions or whether they show crossover behavior. For transitions that exhibit first order characteristics, it would be interesting to prepare a small patch of pinning sites with different characteristics that could act as a nucleation site for one of the phases in order to understand whether there is a length scale analogous to a critical nucleus size.

We also note that the dynamics we observe should be general to systems with similar geometries and repulsively interacting particles. For example, in colloidal systems, square pinning arrays with flat regions between the pinning sites (muffin-tin potentials) have been fabricated, and in these systems the interstitial colloids are much more mobile than in washboard-type pinning potentials. Honeycomb pinning arrays could be created using similar techniques for this type of system.

## Vi Summary

We have shown that vortices in honeycomb pinning arrays exhibit a rich variety of dynamical phases that are distinct from those found in triangular and square pinning arrays. The honeycomb pinning arrays allow for the appearance of -mer type states that have orientational degrees of freedom. We specifically focused on the case where dimer states appear. At , the dimers can have a ferromagnetic type of ordering which is three-fold degenerate. At depinning, the dimers can flow in the direction in which they are aligned. For the case of driving along the -axis, the dimers flow at to the applied drive, giving a transverse velocity response. At incommensurate fields where dimers are present, even though the orientational order is lost, the moving states can dynamically order into a broken symmetry state where the vortices flow with equal probability at either or to the -axis. As the driving in the -direction increases, there is a depinning transition for the vortices in the pinning sites, and the transverse response is lost when the vortices either flow in a random phase or channel along the pinning sites. As a function of pinning force, we find other types of vortex lattice ordering at zero driving, including a partially pinned lattice and a herringbone ordering of the dimers. These other orderings lead to new types of dynamical phases, including an elastic depinning for weak pinning where all the vortices depin simultaneously into a moving crystal phase, and an ordered interstitial flow in which the moving dimers break apart. The transitions between these flow phases appear as clear steps in the velocity force curves, and we have mapped the dynamical phase diagrams for various system parameters. We also showed that the different phases have distinct fluctuations and noise characteristics. When the temperature is high enough, the dimer states lose their orientational ordering and begin to rotate within the interstitial sites. This destroys the symmetry breaking flow; however, the moving locked phase can still occur at high drives.

The transition in the vortex ground state ordering as a function of pinning force causes the critical depinning force for driving in the and -directions to differ. When driving along the -direction, the initial depinning occurs in the form of one-dimensional interstitial channels, and at high drives the vortices can form an anisotropic moving locked phase. We find a large enhancement of the depinning force in the -direction associated with the pinned herringbone phase when the dimers align in the -direction and creates a jamming effect. The jammed state can enhance the critical depinning force by a factor of three, and can also arise for decreasing pinning size. We expect that many of the general features we observe will carry over to the higher matching fields in the honeycomb pinning arrays and in kagomé arrays since ordered -mers states occur for the kagomé lattice as well.

## Vii Acknowledgments

This work was carried out under the auspices of the National Nuclear Security Administration of the U.S. Department of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396.

### References

- H.J. Jensen, A. Brass, and A.J. Berlinsky, Phys. Rev. Lett. 60, 1676 (1988); M.C. Faleski, M.C. Marchetti, and A.A. Middleton, Phys. Rev. B 54, 12427 (1996); C.J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 80, 2197 (1998); K.E. Bassler, M. Paczuski, and E. Altshuler, Phys. Rev. B 64, 224517 (2001); E. Olive and J.C. Soret, ibid. 77, 144514 (2008).
- A.E. Koshelev and V.M. Vinokur, Phys. Rev. Lett. 73, 3580 (1994).
- S. Bhattacharya and M.J. Higgins, Phys. Rev. Lett. 70, 2617 (1993).
- P. Le Doussal and T. Giamarchi, Phys. Rev. B 57, 11356 (1998).
- L. Balents, M.C. Marchetti, and L. Radzihovsky, Phys. Rev. Lett. 78, 751 (1997); Phys. Rev. B 57, 7705 (1998).
- K. Moon, R.T. Scalettar, and G.T. Zimányi, Phys. Rev. Lett. 77, 2778 (1996); S. Spencer and H.J. Jensen, Phys. Rev. B 55, 8473 (1997); C.J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 81, 3757 (1998); A.B. Kolton, D. Domínguez, and N. Grønbech-Jensen, ibid. 83, 3061 (1999); H. Fangohr, S.J. Cox, and P.A.J. de Groot, Phys. Rev. B 64, 064505 (2001).
- F. Pardo, F. de la Cruz, P.L. Gammel, E. Bucher, and D.J. Bishop, Nature (London) 396, 348 (1998).
- A.M. Troyanovski, J. Aarts, and P.H. Kes, Nature (London) 399, 665 (1999).
- Q.-H. Chen and X. Hu, Phys. Rev. Lett. 90, 117005 (2003); A.K. Ghosh, P. Olsson, and S. Teitel, ibid. 97, 267002 (2006).
- A.C. Marley, M.J. Higgins, and S. Bhattacharya, Phys. Rev. Lett. 74, 3029 (1995).
- Y. Togawa, R. Abiru, K. Iwaya, H. Kitano and A. Maeda, Phys. Rev. Lett. 85, 3716 (2000); S. Okuma, J. Inoue, and N. Kokubo, Phys. Rev. B 76, 172503 (2007).
- A.T. Fiory, A.F. Hebard, and S. Somekh, Appl. Phys. Lett. 32, 73 (1978).
- V.V. Metlushko, M. Baert, R. Jonckheere, V.V. Moshchalkov, and Y. Bruynseraede, Sol. St. Commun. 91, 331 (1994).
- M. Baert, V.V. Metlushko, R. Jonckheere, V.V. Moshchalkov, and Y. Bruynseraede, Phys. Rev. Lett. 74, 3269 (1995).
- K. Harada, O. Kamimura, H. Kasai, T. Matsuda, A. Tonomura, and V.V. Moshchalkov, Science 274, 1167 (1996).
- E. Rosseel, M. Van Bael, M. Baert, R. Jonckheere, V.V. Moshchalkov, and Y. Bruynseraede, Phys. Rev. B 53, R2983 (1996).
- L. Van Look, E. Rosseel, M.J. Van Bael, K. Temst, V.V. Moshchalkov, and Y. Bruynseraede, Phys. Rev. B 60, R6998 (1999).
- V.V. Metlushko, L.E. DeLong, M. Baert, E. Rosseel, M.J. Van Bael, K. Temst, V.V. Moshchalkov, and Y. Bruynseraede, Europhys. Lett. 41, 333 (1998).
- V. Metlushko, U. Welp, G.W. Crabtree, R. Osgood, S.D. Bader, L.E. DeLong, Z. Zhang, S.R.J. Brueck, B. Ilic, K. Chung, and P.J. Hesketh, Phys. Rev. B 60, R12585 (1999).
- S.B. Field, S.S. James, J. Barentine, V. Metlushko, G. Crabtree, H. Shtrikman, B. Ilic, and S.R.J. Brueck, Phys. Rev. Lett 88, 067003 (2002); A.N. Grigorenko, S.J. Bending, M.J. Van Bael, M. Lange, V.V. Moshchalkov, H. Fangohr, and P.A.J. de Groot, ibid. 90, 237001 (2003).
- J.I. Martín, M. Vélez, J. Nogués, and I.K. Schuller, Phys. Rev. Lett. 79, 1929 (1997).
- U. Welp, Z.L. Xiao, J.S. Jiang, V.K. Vlasko-Vlasov, S.D. Bader, G.W. Crabtree, J. Liang, H. Chik, and J.M. Xu, Phys. Rev. B 66, 212507 (2002).
- J.I. Martín, M. Vélez, A. Hoffmann, I.K. Schuller, and J.L. Vicent, Phys. Rev. Lett. 83, 1022 (1999).
- G. Karapetrov, J. Fedor, M. Iavarone, D. Rosenmann, and W.K. Kwok, Phys. Rev. Lett. 95, 167002 (2005).
- D.J. Morgan and J.B. Ketterson, Phys. Rev. Lett. 80, 3614 (1998).
- T.C. Wu, J.C. Wang, L. Horng, J.C. Wu, and T.J. Yang, J. Appl. Phys. 97, 10B102 (2005).
- V.R. Misko, S. Savel’ev, and F. Nori, Phys. Rev. B 74, 024522 (2006); M. Kemmler, C. Gürlich, A. Sterck, H. Pöhler, M. Neuhaus, M. Siegel, R. Kleiner, and D. Koelle, Phys. Rev. Lett. 97, 147003 (2006); A.V. Silhanek, W. Gillijns, V.V. Moshchalkov, B.Y. Zhu, J. Moonens, and L.H.A. Leunissen, Appl. Phys. Lett. 89, 152507 (2006).
- J.E. Villegas, M.I. Montero, C.-P. Li, and I.K. Schuller, Phys. Rev. Lett. 97, 027002 (2006); J. Eisenmenger, M. Oettinger, C. Pfahler, A. Plettl, P. Walther, and P. Ziemann, Phys. Rev. B 75, 144514 (2007); C. Reichhardt and C.J. Olson Reichhardt, ibid. 76, 094512 (2007).
- C. Reichhardt, C.J. Olson, and F. Nori, Phys. Rev. B 57, 7937 (1998).
- C.J. Olson Reichhardt, A. Libál, and C. Reichhardt, Phys. Rev. B 73, 184519 (2006).
- G.R. Berdiyorov, M.V. Milosevic, and F.M. Peeters, Phys. Rev. B 74, 174512 (2006).
- W.V. Pogosov, A.L. Rakhmanov, and V.V. Moshchalkov, Phys. Rev. B 67, 014532 (2003); G.R. Berdiyorov, M.V. Milosevic, and F.M. Peeters, ibid. 76, 134508 (2007).
- C. Reichhardt, C.J. Olson, and F. Nori, Phys. Rev. Lett. 78, 2648 (1997); Phys. Rev. B 58, 6534 (1998).
- C. Reichhardt and N. Grønbech-Jensen, Phys. Rev. B 63, 054510 (2001).
- C. Reichhardt, G.T. Zimányi, and N. Grønbech-Jensen, Phys. Rev. B 64, 014501 (2001).
- C. Reichhardt, R.T. Scalettar, G.T. Zimányi, and N. Grønbech-Jensen, Phys. Rev. B 61, R11914 (2000).
- C. Reichhardt, A.B. Kolton, D. Domínguez, and N. Grønbech-Jensen, Phys. Rev. B 64, 134508 (2001); V.I. Marconi, A.B. Kolton, D. Domínguez, and N. Grønbech-Jensen, ibid. 68, 104521 (2003).
- G. Carneiro, Phys. Rev. B 66, 054523 (2002).
- B.Y. Zhu, L. Van Look, V.V. Moshchalkov, B.R. Zhao, and Z.X Zhao, Phys. Rev. B 64, 012504 (2001).
- L. Van Look, B.Y.Zhu, R. Jonckheere, B.R. Zhao, Z.X. Zhao, and V.V. Moshchalkov, Phys. Rev. B 66, 214511 (2002).
- M. Benkraouda, I.M. Obaidat, U. Al Khawaja, and N.M.J. Mulaa, Supercond. Sci. Technol. 19, 368 (2006).
- Q.H. Chen, G. Teniers, B.B. Jin, and V.V. Moshchalkov, Phys. Rev. B 73, 014506 (2006).
- V.R. Misko, S. Savel’ev, A.L. Rakhmanov, and F. Nori, Phys. Rev. B 75, 024509 (2007).
- T.C. Wu, P.C. Kang, L. Horng, J.C. Wu, and T.J. Yang, J. Appl. Phys. 95, 6696 (2004).
- R. Cao, T.C. Wu, P.C. Kang, J.C. Yang, T.J. Yang, and L. Horng, Sol. St. Commun. 143, 171 (2007).
- K.D. Fisher, D. Stroud, and L. Janin, Phys. Rev. B 60, 15371 (1999); V.I. Marconi and D. Domínguez, Phys. Rev. Lett. 82, 4922 (1999).
- C. Reichhardt and G.T. Zimányi, Phys. Rev. B 61, 14354 (2000); V. Gotcheva, A.T.J. Wang, and S. Teitel, Phys. Rev. Lett. 92, 247005 (2004).
- R. Surdeanu, R.J. Wijngaarden, R. Griessen, J. Einfeld, and R. Wördenweber, Europhys. Lett. 54, 682 (2001).
- C.C. de Souza Silva, J. Van de Vondel, B.Y. Zhu, M. Morelle, and V.V. Moshchalkov, Phys. Rev. B 73, 014507 (2006); Q. Lu, C.J. Olson Reichhardt, and C. Reichhardt, ibid. 75, 054502 (2007); L. Dinis, E.M. Gonzalez, J.V. Anguita, J.M.R. Parrondo, and J.L. Vicent, New. J. Phys. 9, 366 (2007).
- C. Reichhardt and C.J. Olson Reichhardt, Phys. Rev. B 76, 064523 (2007).
- C. Reichhardt and C.J. Olson, Phys. Rev. Lett. 88, 248301 (2002).
- M. Brunner and C. Bechinger, Phys. Rev. Lett. 88, 248302 (2002).
- R. Agra, F. van Wijland, and E. Trizac, Phys. Rev. Lett. 93, 018304 (2004).
- A. Sarlah, T. Franosch, and E. Frey, Phys. Rev. Lett. 95, 088302 (2005).
- S. El Shawish, J. Dobnikar, and E. Trizac, Soft Matter 4, 1491 (2008).
- A.B. Harris and A.J. Berlinsky, Can. J. Phys. 57, 1852 (1979); E.J. Nicol, C. Kallin, and A.J. Berlinsky, Phys. Rev. B 38, 556 (1988); P. Zeppenfeld, J. Goerge, V. Diercks, R. Halmer, R. David, G. Comsa, A. Marmier, C. Ramseyer, and C. Girardet, Phys. Rev. Lett. 78, 1504 (1997).
- M.F. Laguna, C.A. Balseiro, D. Domínguez, and F. Nori, Phys. Rev. B 64, 104505 (2001).
- C. Reichhardt and C.J. Olson Reichhardt, Phys. Rev. Lett. 100, 167002 (2008).
- H. Pu, L.O. Baksmaty, S. Yi, and N.P. Bigelow, Phys. Rev. Lett. 94, 190401 (2005); J.W. Reijnders and R.A. Duine, Phys. Rev. A 71, 063607 (2005).
- S. Tung, V. Schweikhard, and E.A. Cornell, Phys. Rev. Lett. 97, 240402 (2006).
- K. Kasamatsu and M. Tsubota, Phys. Rev. Lett. 97, 240404 (2006).
- P.T. Korda, G. C. Spalding, and D.G. Grier, Phys. Rev. B 66, 024504 (2002); K. Mangold, P. Leiderer, and C. Bechinger, Phys. Rev. Lett. 90, 158302 (2003).
- G. Coupier, M. Saint Jean, and C. Guthmann, Phys. Rev. B 75, 224103 (2007).
- O.M. Braun, T. Dauxois, M.V. Paliy, and M. Peyrard, Phys. Rev. Lett. 78, 1295 (1997); E. Granato and S.C. Ying, Phys. Rev. B 69, 125403 (2004); J. Tekic, O.M. Braun, and B. Hu, Phys. Rev. E 71, 026104 (2005).