Physica D 193, 303 (2004)

Dynamics and Melting of Stripes, Crystals, and Bubbles with Quenched Disorder

C. J. Olson Reichhardt*, C. Reichhardt, I. Martin, and A.R. Bishop

Center for Nonlinear Science and Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545.


Abstract
Two-dimensional systems in which there is a competition between long-range repulsion and short range attraction exhibit a remarkable variety of patterns such as stripes, bubbles, and labyrinths. Such systems include magnetic films, Langmuir monolayers, polymers, gels, water-oil mixtures, and two-dimensional electron systems. In many of these systems quenched disorder from the underlying substrate may be present. We examine the dynamics and stripe formation in the presence of both an applied dc drive and quenched disorder. When the disorder strength exceeds a critical value, an applied dc drive can induce a dynamical stripe ordering transition to a state that is more ordered than the originating undriven, unpinned pattern. We also consider the melting transition of the different patterns.
PACS: 71.45.Lr; 64.60.Cn; 73.20.Qt; 89.75.Kd
Keywords: Two-dimensional systems; Quenched disorder; Short- and long-range interactions


The growing interest in nonextensive thermodynamics and statistical mechanics represented in these Proceedings is stimulated by observations (from real experiments or simulations of models) of phenomena such as incomplete (e.g., fractal) coverage of phase space, rare-event statistics, intermittency, stochastic resonance, and so forth. We can hope for, and certainly seek to find, underpinning mechanisms and classifications of such systems.
We suggest that among the relevant classifications are (a) systems exhibiting mesoscopic self-assembly (usually from nonlinearity or competing scales at more microscopic scales), and (b) systems with long-range interactions (and often coexisting short-range interactions and entropic forces). Long-range interactions may be in the microscopic variables of the system (e.g., Coulomb fields), or they may be induced at mesoscopic scales because of local constraints or symmetries (e.g., vortices in fluids and magnets, patterns in convection cells and vibrated granular materials, strain-strain interactions in elastic media). It is important to distinguish topological (singularity) mesoscale excitations, which have long-range mutual interactions (such as vortices, dislocations, stress filaments) from mesoscopic collective excitations with short-range interactions (such as solitons, breathers). We expect that the former will exhibit mesoscopic statistics driven by their mutual interactions, whereas the latter should show mesoscopic statistics driven by the (colored, multiplicative) noise bath of the non-collective degrees of freedom.
We have studied aspects of anomalous statistics and dynamics in a number of physically motivated models, including: soliton diffusion [1]; relaxation in two and three dimensional complex Ginzburg-Landau equations [2]; multiscale structure and dynamics in elastic media [3]; filamentary transport in superconductor flux lattices [4]; surface morphology and growth [5]; relaxation in magnets with planar spin anisotropy [6]; glassy dynamics near a spinodal in a Van der Waals fluid [7]; charge, spin ordering and dynamics in strongly correlated electronic materials [8]; and large-scale materials deformation (fracture, friction) [9].
Motivated by the similarities in the phenomena exhibited by these (and related) model systems, we have begun studies of a "minimal" model with coexisting short- and long-range interactions. Here we describe this model, which has also been discussed in Refs. [10, 11], and briefly report some of its properties.
Recent proposals for ordering in two-dimensional (2D) electron systems consider the idea of an electron liquid crystal system [12, 13, 14]. Such a system can form stripe phases with varying degrees of disorder: An anisotropic Wigner crystal is well-ordered in both directions; a smectic phase has ordering along only one direction; in a nematic, dislocations are introduced into the smectic phase but there is still an overall preferred orientation for the system; finally, in an isotropic liquid stripe phase, the stripe pattern is filled with dislocations and has no unique orientational direction.
Stripe patterns are very general and are formed in a wide range of systems [15]. For example, the Swift-Hohenberg model produces glassy stripe configurations [16]. Other stripe-forming systems include magnetic films [17], Langmuir monolayers, polymers, gels, and water-oil mixtures [18]. One system in which stripes have attracted a particularly large amount of interest is holes in layered transition metal oxides. In a quasiclassical description, the interaction between two holes consists of both a long-range repulsion generated by the Coulomb interaction, as well as a short-range attraction that arises due to the breaking of the antiferromagnetic bonds surrounding the holes. As a result, the holes can form linear arrays termed stripes. There is experimental evidence for the presence of stripes in several of the lanthanum cuprate oxides, but it is not yet clear whether the high-Tc superconducting compounds such as YBCO contain stripes, nor what the role of the stripes would be if they exist. Experimental detection of stripes is made more difficult by the possible interaction of impurities in the material with the stripes. The effect of disorder on a stripe phase is not well understood. Therefore we have constructed a simple model containing the basic physics of the system in order to gain a better understanding of the ordering and dynamics of a stripe phase.
Our simple model consists of particles with competing long-range and short-range interactions [10, 11]. The particles experience a long-range repulsion of a Coulomb form, and a short-range attraction which we assume to be a simple exponential form. When these two interactions are combined, a new length scale emerges which depends on the density of the particles and the strength of the short-range attractive term. We further add to our model quenched disorder to represent impurities in the material, as well as an external driving field, which would be an electric field in the case of holes.
The particles are assumed to obey overdamped Langevin dynamics given by
fi = − Nc

j=1 
U(ρ) + fT + fp + fd = ηvi,
(1)
where the number of particles is Nc and we take the damping coefficient η = 1. The potential between two particles a distance r apart is
U(ρ)= q

r
B exp(−r/ξ) ,
(2)
where q=1 is the charge on each particle and ξ is the magnetic screening length. The parameter B can be used to tune the relative strength of the short-range attractive term. The long range interactions are treated as in Ref. [19]. Temperature is represented by random thermal kicks which have the property, < f(t)T > =0, < f(t)iTf(t′)jT > =2ηkB T δijδ(tt′). The disorder is modeled as random, point-like pins, represented by parabolic traps. The pinning force is given by
fip= Np

k=1 
fp

rp
|rirk(p)|Θ(rp − |rirk(p)|)
^
r
 

i,k 
,
(3)
where fp=0.2 is the pinning force, rk(p) is the location of pin k, rp=0.125ξ is the pin radius, Θ is the Heaviside step function, ri,k=(rirk(p))/|rirk(p)|, and the number of pins Np is greater than the number of particles, Np > Nc. The force from an applied field (voltage) is modeled as a uniform drive, fd=fd x. We measure the velocity signature < Vx > of the system, where Vx = (1/Nc) ∑iNcvi ·x, which would correspond to a current.
Fig1.png
Figure 1: Static phases as a function of increasing strength of the short range attraction, B, at n=0.64. (a) Wigner crystal phase at B=0. (b) Stripe phase at intermediate B=0.3. (c) Clump phase at high B=0.4.
Fig2.png
Figure 2: Static phases for increasing particle density with fixed B=1. (a) Nc=18; (b) Nc=100; (c) Nc=400; (d) Nc=1050.
We initially consider a system with no pinning or external drive. We vary the strength of the short range attraction by changing the value of B, while holding the particle density fixed at n=0.64/ξ2 for Nc=132. As shown in Fig. 1(a), taken from Ref. [11], when B is small, the long range repulsion dominates and the system forms a Wigner crystal. At high values of B, Fig. 1(c), the short range attraction dominates and the particles form clumps which become larger as the value of B is further increased. The clumps organize into a superstructure of a triangular lattice of clumps. At intermediate values of B, Fig. (b), the short and long range interactions are of comparable strength, and the system forms a disordered stripe structure. We observe the same phases if we fix the value of B and vary the density of the system. This is illustrated in Fig. 2.
Fig3.png
Figure 3: Melting temperature for the three phases as a function of B at n=2.0.
Fig4.png
Figure 4: Snapshots at different temperatures illustrating the melting of the stripe phase for B=4.75 at n=2.0. Black dots are the particles; lines indicate the trajectories of the particles over a fixed interval of time which is the same for all panels.
We have measured the melting temperature of the three phases by computing the particle diffusion at each temperature. We define the melting temperature to be the temperature at which the particles first diffuse more than a lattice constant. The melting temperature Tm as a function of B is plotted in Fig. 3 for a system with particle density n=2.0. Tm decreases in the crystal phase as B increases since the short range attractive term introduces distortions in the crystal, allowing it to melt more easily. The stripe phase has the lowest melting temperature, but Tm increases dramatically upon entering the clump phase. This is because the clump phase contains a large amount of elastic energy. Within each clump there is triangular ordering of the particles, and additionally the clumps form an ordered superstructure. We find that the clump phase simultaneously melts and dissociates for the parameters considered here. In contrast, in the stripe phase melting first occurs along the length of each stripe, as illustrated in Fig. 4. Particles diffuse freely within the stripes but do not enter the regions between the stripes until higher temperatures are applied. Thus we find evidence for a two-stage melting of the stripe phase.
Fig5.png
Figure 5: Depinning force Fdp for the three phases as a function of B at n=0.64.
We next consider the effect of disorder on the three phases. We add random point-like disorder to the system, and then apply a driving force to the particles. At small drive strengths, the particles are trapped by the disorder and do not move. We measure the threshold force required to depin the particles in each phase. As shown in Fig. 5, taken from Ref. [11], the depinning force increases with increasing B in the crystal phase, is highest in the stripe state, and then drops dramatically upon entering the clump phase. This is the same behavior, but inverted, that was observed for the melting transition in Fig. 3. This similarity is due to the fact that the softness of the stripe phase, which allows it to melt at lower temperature, also allows it to be well pinned by the randomly located pinning sites when, as considered here, there are more pinning sites than particles. The stripe state is disordered and can readily readjust to take advantage of the maximum number of pinning sites. In contrast, in the highly elastic and stiff clump phase, the particles are so constrained by the ordering of the phase that it is very difficult for them to adjust in order to find the lowest energy pinning location. As a result many of the pins are ineffective in the clump phase and the depinning force is low. This similarity between the melting and depinning signatures is particularly interesting as it implies that information about a thermal transition (melting) can be obtained from a dynamical experiment at constant temperature (depinning).
Fig6.png
Figure 6: Dynamic reordering of the clump phase for B=0.325 and n=0.64. The external driving force increases from Fd=0.0 in panel (a) to Fd=0.06 in panel (f).
At high applied drives, it is possible to (dynamically) reorder the phases, which are distorted by the presence of the pinning at low drives. Dynamical reordering of a triangular lattice of repulsively interacting particles has been observed previously in several systems, including superconducting vortices [20, 21] and Wigner crystals [22]. At high drives, the force from the pinning begins to resemble a temperature which decreases as the drive increases, until at infinite drive the effective temperature is zero, leading to the concept of dynamical freezing [20]. Dynamical freezing has already been observed in the Wigner crystal phase. We find that the clump phase also undergoes a dynamical freezing transition [11], and that at high drives the clump phase in the presence of disorder has the same form as the static clump phase in the absence of disorder, as shown in Fig. 6.
Fig7.png
Figure 7: Dynamic reordering of the stripe phase for B=0.29 and n=0.64. The external driving force increases from Fd=0.0 in panel (a) to Fd=0.06 in panel (f).
When we consider the dynamical reordering of the stripe phase, illustrated in Fig. 7, we find that the final ordered stripe phase is more ordered than the corresponding stripe phase that appears in a static system without pinning [10]. The stripes are aligned due to the driving force, rather than being randomly aligned as was the case in the thermally annealed system. This indicates that, in a stripe forming system, it may be possible to use a combination of quenched disorder and an external drive to produce well-ordered stripes. The reordering of the stripe phase occurs only if the disorder is sufficiently strong to cause the stripe to undergo disordered, plastic flow at low drives. The original ordering of the stripe must be destroyed by the disorder in order to allow for a new ordering of the stripe to form aligned with the drive. As a result of the reordering of the stripe phase, only the stripe phase shows hysteresis if the driving force is increased and then decreased again. The other two phases, clump and Wigner, do not exhibit any hysteresis.
In summary, the competition between short range attraction and long range repulsion leads to the formation of Wigner crystal, stripe, and clump phases as a function of either the strength of the short range term or the particle density. The melting temperature is lowest in the stripe phase, while the depinning force is highest in the stripe phase, and lower in the Wigner and clump phases. Quenched disorder distorts the ordered phases, but the structure of all three phases can be regained through dynamical reordering at high drives. The dynamically reordered stripe phase is more ordered than the equilibrium stripe phase in a clean system, and as a result the stripe phases exhibits hysteretic signatures which could be observed experimentally.
In the context of "nonextensive thermodynamics," it is now clearly of interest to study the various mesoscopic phases of our model (and their dynamics), and to compare their statistics and responses with predictions of Tsallis or other superstatistics.

Acknowledgments

This work was supported by the US Department of Energy under contract W-7405-ENG-36.

References

[1]
See, e.g., K.Ø. Rasmussen, S. Aubry, A.R. Bishop, and T. Tsironis, Europ. J. Phys. B 15 (2000) 169.
[2]
I.S. Aranson, A.R. Bishop, and L. Kramer, Phys. Rev. E 57 (1998) 5276; P. Kevrekidis, A.R. Bishop, and K.Ø. Rasmussen, Phys. Rev. E 65 (2001) 16122.
[3]
S.R. Shenoy, T. Lookman, A. Saxena, and A.R. Bishop, Phys. Rev. B 60 (1999) R12537; K.Ø. Rasmussen, T. Lookman, A. Saxena, A.R. Bishop, R.C. Albers, and S.R. Shenoy, Phys. Rev. Lett. 87 (2001) 5704; T. Lookman, S.R. Shenoy, K.O. Rasmussen, A. Saxena, and A.R. Bishop, Phys. Rev. B, in press.
[4]
N. Grønbech-Jensen, A.R. Bishop, and D. Domínguez, Phys. Rev. Lett. 76 (1996) 2985; C.J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 80 (1998) 2197.
[5]
I.S. Aranson, A.R. Bishop, I. Daruka, and V.M. Vinokur, Phys. Rev. Lett. 80 (1998) 1770; A. Sanchez, A.R. Bishop, and E. Moro, Phys. Rev. E 62 (2000) 3219.
[6]
S. Komineas et al., Europhys. Lett., in press.
[7]
K.-K. Loh, K. Kawasaki, A.R. Bishop, T. Lookman, A. Saxena, J. Schmalian, and Z. Nussinov, cond-mat/0206494.
[8]
B.P. Stojkovic, Z.G. Yu, A.L. Chernyshev, A.R. Bishop, A.H. Castro Neto, and N. Grønbech-Jensen, Phys. Rev. B 62 (2000) 4353; S. Teber, B.P. Stojkovic, S.A. Brazovskii, and A.R. Bishop, J. Phys. Cond. Mat. 13 (2001) 4015.
[9]
J. Röder, J.E. Hammerberg, B.L. Holian, and A.R. Bishop, Phys. Rev. B 57 (1998) 2759; O. Braun, A.R. Bishop, and J. Röder, Phys. Rev. Lett. 82 (1999) 3097.
[10]
C. Reichhardt, C.J. Olson Reichhardt, I. Martin, and A.R. Bishop, Phys. Rev. Lett. 90, 026401 (2003).
[11]
C. Reichhardt, C.J. Olson, I. Martin, and A.R. Bishop, Europhys. Lett. 61, 221 (2003).
[12]
M.M. Fogler, A.A. Koulakov, and B.I. Shklovskii, Phys. Rev. B 54 (1996) 1853.
[13]
S.A. Kivelson, E. Fradkin, and V.J. Emery, Nature 393 (1998) 550.
[14]
E. Fradkin and S.A. Kivelson, Phys. Rev. B 59 (1999) 8065.
[15]
M. Seul and D. Andelman, Science 267 (1995) 476.
[16]
D. Boyer and J. Viñals, Phys. Rev. E 65, 046119 (2002).
[17]
M. Seul and R. Wolfe, Phys. Rev. A 46 (1992) 7519.
[18]
W.M. Gelbart and A. Ben Shaul, J. Phys. Chem. 100 (1996) 13169.
[19]
N. Grønbech-Jensen, Int. J. Mod. Phys. C 7 (1996) 873.
[20]
A.E. Koshelev and V.M. Vinokur, Phys. Rev. Lett. 73 (1994) 3580.
[21]
K. Moon, R.T. Scalettar, and G.T. Zimányi, Phys. Rev. Lett. 77 (1996) 2778; C.J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 81 (1998) 3757; A.B. Kolton, D. Domínguez, and N. Grønbech-Jensen, Phys. Rev. Lett. 83 (1999) 3061.
[22]
C. Reichhardt, C.J. Olson, N. Grønbech-Jensen, and F. Nori, Phys. Rev. Lett. 86 (2001) 4354.



File translated from TEX by TTHgold, version 4.00.
Back to Home