arxiv:cond-mat/ v1 [cond-mat.soft] 11 Dec 2006

Transcription

1 Grain boundary dynamic in tripe phae of non potential ytem Zhi-Feng Huang Department of Phyic and tronomy, Wayne State Univerity, Detroit, MI and McGill Intitute for dvanced Material, and Department of Phyic, McGill Univerity, Montreal, QC H3 2T8, Canada arxiv:cond-mat/ v1 [cond-mat.oft] 11 Dec 2006 Jorge Viñal McGill Intitute for dvanced Material, and Department of Phyic, McGill Univerity, Montreal, QC H3 2T8, Canada (Dated: July 25, 2018) We decribe numerical olution of two non potential model of pattern formation in nonequilibrium ytem to addre the motion and decay of grain boundarie eparating domain of tripe configuration of different orientation. We firt addre wavenumber election becaue of the boundary, and poible decay mode when the periodicity of the tripe phae i different from the elected wavenumber for a tationary boundary. We dicu everal decay mode including long wavelength undulation of the moving boundary a well a the formation of localized defect and their ubequent motion. We find three different regime a a function of the ditance to the tripe phae threhold and initial wavenumber, and then correlate thee finding with domain morphology during domain coarening in a large apect ratio configuration. I. INTRODUCTION tripe configuration i a particular example of a modulated configuration that exhibit tructural periodicity along only one patial direction. In thermodynamic equilibrium, thee phae have been widely oberved in a variety of ytem [1], including, for example, microphae eparated block copolymer[2]. Outide of equilibrium, tripe pattern are alo widepread, including Rayleigh-Bénard convection in fluid [3], or electroconvection in nematic liquid crytal [4]. When the lateral ize of the ytem of interet i large compared with the wavelength of the modulation (large apect ratio limit), a tripe configuration in either cla of ytem i uually characterized by a tranient multidomain and defected configuration with only local order. Domain coarening i believed to be controlled by the motion of exiting topological defect in the configuration uch a grain boundarie, dilocation, and diclination. In thi paper we focu on non potential (or no variational) model ytem and invetigate the procee of defect formation and motion near grain boundarie, and how that there exit everal regime in parameter pace that are characterized by qualitatively different coarening behavior. major topic of interet in domain coarening of tripe pattern i the growth law for the characteritic linear cale of the configuration. Thi ha been addreed both in theoretical tudie [5, 6, 7, 8, 9] and in experiment in nematic [4, 10] and block copolymer [11]. Thi characteritic length cale R(t) i expected to obey a power law growth in time: R(t) t x, with a value of the coarening exponent x which range from 1/5 to 1/2 in the tudie to date. The actual value appear to depend on ǫ (the dimenionle ditance from the onet of the tripe phae) [8, 12], thermal noie [5], and the type of length cale under tudy [6, 9]. lo, domain morphologie a well a the value of x have been found to be different in potential ytem (in which the evolution i olely driven by minimization of a potential) than in the non potential cae [6]. Numerical tudie of the Swift-Hohenberg model of Rayleigh-Bénard convection (a potential ytem) ugget that x = 1/3 in the limit ǫ 1 a a reult of grain boundary motion[7]. ǫ increae, the value of x i een to decreae from 1/3, fact that wa argued to be a manifetation of defect pinning [12]. On the other hand, the generalized Swift-Hohenberg model ha been tudied a an example of a non potential model [6], and it ha been hown that the evolution of tripe pattern i aymptotically dominated by the motion of iolated dilocation. Furthermore, recent reult from a direct olution of the Bouineq equation in a Rayleigh-Bénard configuration ugget that multiple length cale hould be introduced in the cae of non potential ytem [9]. Depite a ignificant body of reearch on domain coarening of modulated pattern, the mechanim reponible for the oberved propertie and the value of the coarening exponent, epecially for non potential ytem, are not yet undertood. We report in thi paper our analyi of the relation between the elected wavenumber at grain boundarie and the critical wavenumber, and conequence on both early time defect formation and aymptotic coarening. Our tarting point i a planar 90 grain boundary eparating two tripe domain of mutually perpendicular orientation, a hown in Fig. 1. We determine the elected wavevector for a tationary boundary configuration q in two different non potential model, and how that for other wavenumber the mode of decay of the configuration can be claified according to it wavenumber and the dimenionle ditance from threhold ǫ. For mall ǫ, an initially planar grain boundary eparating two domain, one of which at leat ha q q, i een to propagate, either a a planar front, or by developing long wavelength undulation, depending on the value of q. For larger value of ǫ we find that the

2 2 q y B q x FIG. 1: 90 tilt grain boundary (GB) configuration of a tripe phae, with wave number q x and for two domain and B of mutually perpendicular tripe orientation. grain boundary configuration decay through the formation of diclination around the boundary and their ubequent motion. We then correlate our finding with the oberved domain coarening morphology of configuration obtained from initially uniform but untable tate. The rationale for the comparion i the obervation that in an extended ytem an early time tranient configuration following the decay of an untable tate will comprie a ditribution of locally ordered domain with wavenumber approximately equal to the critical wavenumber for intability q 0, but with different orientation. Depending on whether q 0 i cloe to q or not, the reulting boundarie will be locally cloe to being tationary, or quickly decay leaving behind different type of defect. The enuing domain coarening i hown to be different depending on the (near) tability of grain boundarie that pontaneouly form in the early tage following the decay of the initially untable configuration. Our tudy i baed on a numerical olution of two different non potential extenion of the Swift-Hohenberg model of Rayleigh-Bénard convection [3], a decribed in Sec. II. One of them i made non potential on account of the pecific choice of nonlinearity, wherea the other allow for mean flow repreenting the effect of vertical vorticity encountered in Rayleigh-Bénard convection in mall Prandtl number fluid [13]. They are decribed in Section III and IV repectively. lthough wave number election for a tationary grain boundary i different in both model, imilar propertie concerning grain boundary motion and decay have been found. They correlate with two type of domain coarening behavior, either grain boundary dominated or diclination/dilocation dominated, a deduced from the reult on grain boundary decay that will be dicued in Sec. V. II. NON POTENTIL MODEL EQUTIONS The model ytem we have invetigated are two dimenional generalized Swift-Hohenberg equation which have been ued a a phenomenological decription of tripe pattern [3, 13]: t ψ +U ψ = ǫψ ( 2 +q 2 0) 2 ψ gψ 3 +3(1 g) ψ 2 2 ψ, (1) where ψ(x,y,t) i a dimenionle real order parameter field, and the control parameter ǫ i the ditance from the threhold for intability of the olution ψ = 0 (at ǫ = 0). The wavenumber q 0 (q 0 = 1 in the dimenionle unit that we are uing) i the critical wave number for intability. n advection term with a drift velocity U(x,y,t) = (U x,u y ) ha been added to Eq. (1) to introduce mean flow effect reulting from the coupling to vertical vorticity in the fluid [13, 14]. The velocity U(x,y,t) i given in term of the vorticity potential ζ(x,y,t): which atifie, U = (ζẑ) = ( y ζ, x ζ), (2) [ t σ( 2 c 2 ) ] 2 ζ = g m [ ( 2 ψ) ψ ] ẑ, (3) with σ beingthe Prandtlnumberofthe fluid, c 2 introducedto phenomenologicallymodel the effect ofnolip boundary condition at the top and bottom plate of the convection cell [14], and the coupling parameter g m choen to be a decreaing function of the Prandtl number. Following Ref. 13 a Gauian filtering operator F γ i introduced on the r.h.. of Eq. (3) to reduce high wavevector contribution to the vertical vorticity. In Fourier pace, F γ act by multiplying the Fourier tranform of the filtered function by a factor ˆF γq = e γ2 q 2,

3 where q 2 = q 2 x +q2 y, with wavevector component q x and, and γ repreent the filtering radiu. Different type of model equation for decribing tripe pattern have been included in Eq. (1) (3), depending on the value of g and g m. For g = 1 and g m = 0 (o that U = 0), the nonlinearity in the model i given by the cubic term ψ 3 only, leading to the original Swift-Hohenberg equation [15], which i a potential model. non potential model without mean flow follow from g = g m = 0 (referred to a model 1 in thi paper), which include a nonlinear term of the form ψ 2 2 ψ. For g = 1 and g m 0 (and hence U 0), we obtain the Swift-Hohenberg model upplemented by mean flow, alo a non potential ytem (model 2). The initial configuration conidered in the firt part of our tudy i a 90 tilt boundary eparating two tripe domain, a hown in Fig. 1. tationary tate i known to exit [16] if unique value of the wavenumber q x (for the tripe parallel to the grain boundary in domain ) and (for the tripe perpendicular to the boundary in domain B) are elected in both domain. For the original Swift-Hohenberg model, thee wavenumber are q x = = q 0. In our numerical tudy, we conider intead a ymmetric pair of grain boundarie along the x direction o that periodic boundary condition can be adopted. The peudopectral algorithm decribed in Ref. 17 ha been ued to numerically olve the model equation, with a time tep t = 0.2 and a grid pacing x = y = λ 0 /8 (i.e., 8 grid point per tripe wavelength λ 0 = 2π/q 0 ). Mot of our reult hown below correpond to a ytem ize grid node, with pot check with larger ize ( ) to verify our reult. 3 III. WVEVECTOR SELECTION ND GRIN BOUNDRY DECY - MODEL 1. Wavenumber election Model 1 include the nonlinear term ψ 2 2 ψ but doe not have mean flow. The elected wave number (qx and qy ) for a tationary grain boundary configuration have been already calculated by Teauro and Cro [16]. Their analyi of the appropriate amplitude equation for thi cae yield q x = q 0 ǫ 2q 3 0 up to order ǫ for the elected wavenumber in domain, and qy = q 0 up to order ǫ 3/4 for domain B. They alo conducted a limited numerical tudy of thi model to ugget that qy q ǫ (although thi reult i baed on calculation involving only two different value of ǫ). We have numerically calculated the elected wave number for the ame configuration but with much larger ytem ize and longer integration time. Our reult for qx are in agreement with Eq. (4), but the calculated value of qy are different from thoe of previou tudie, a preented in Fig. 2 and 3. The deviation of from the critical wavenumber q 0 i een to obey two different power law in ǫ. For mall ǫ < ǫ 0, with ǫ 0 around 0.2, we find (Fig. 3) while for larger ǫ we obtain, q y q ǫ, (5) q y q ǫ 3/2. (6) Note that although both reult differ from the numerical calculation of Ref. 16, they are not inconitent with an analyi baed on amplitude equation. (4) B. Decay of the grain boundary configuration Given the reult of Eq. (4) (6), and alo our numerical reult in Fig. 2, it i clear that the wave number elected by a grain boundary i different from q 0, the mot untable mode near onet. Thu the locally dominant wavenumber during the early tage of coarening could be ignificantly different from (qx, ), and therefore it appear to be of value to tudy the decay of grain boundarie under thoe condition. Of coure, grain boundarie would not be fully formed at thoe early tage, but we hope that tudying untable trajectorie away from well defined grain boundary configuration will illutrate difference in oberved morphologie between potential and non potential model. We tart by noting that in the configuration which we conider there i no phae conervation in the x direction. Hence if q x qx i choen initially, the bulk wave number of domain can readjut to qx through lamella formation or detruction near the boundary [16]. We will not conider thi cae. More complicated phenomena can occur for qy, and thi i the focu of our tudy: In the calculation that follow we alway keep q x = qx in domain, while we explore a range of value of in domain B.

4 4 ε I Zigzag (bulk) =q 0 (linear intability) III 0.1 Eckhau (bulk) II ε T FIG. 2: Stability diagram for the non potential model 1 a a function of ǫ and (the wave number of the tripe perpendicular to the grain boundary). The plu ymbol correpond to the unique wavevector (q x,) = (q x = q 0 ǫ/(2q 3 0),q y) for a tationary grain boundary configuration. lo hown (a dot-dahed line) are tability boundarie for bulk Eckhau and zig-zag intabilitie. Three regime for grain boundary decay are indicated, and alo the critical wave number q 0 for linear intability. 0.1 ε 3/2 q 0 - S 0.01 ε 0.23 ε ε 3/ ε FIG. 3: Wave number q y for a tationary grain boundary (model 1) a a function of ǫ. The figure how two different type of dependency on ǫ: q 0 q y ǫ (for mall ǫ) and ǫ 3/2 (for larger ǫ). < >, ε < ε T >, ε > ε T B q x = q x q x = q x (a) Regime I (b) Regime II (c) Regime III FIG. 4: Grain boundary decay in the three regime hown in Fig. 2 for model 1, with q x = q x. The direction of boundary motion i indicated by arrow. (a) Regime I: < q y for all ǫ; the example hown here correpond to = 0.85, ǫ = 0.3 and at t = (b) Regime II: > q y and ǫ < ǫ T; we have choen = q 0 = 1, ǫ = 0.04, and t = (c) Regime III: > q y and ǫ > ǫ T; the configuration hown here correpond to = q 0 = 1, ǫ = 0.5, and t = 100.

5 5 B = q x q x = (a) t=0 (b) t=500 FIG. 5: Grain boundary decay (for model 1) due to a low tranvere modulation of domain in an initial tationary configuration (with q x = q x and = q y), with ǫ = 0.5, the modulation wavenumber Q = q x/16, and the magnitude of the phae modulation δ = λ 0. indicated in Fig. 2, we find that our reult can be qualitatively claified according to the value of and ǫ. In regime I defined by < qy, the two grain boundarie move from domain toward domain B until they merge. Both boundarie remain approximately planar during the proce. Such a motion i found for all the value of ǫ invetigated, with a particular example of ǫ = 0.3 and = 0.85 hown in Fig. 4a. On the other hand, a longitudinal ditortion of the tripe in domain B around the grain boundary region i oberved if > qy. In addition, the evolution following the appearance of the initial ditortion i qualitatively different depending on the value of ǫ. In regime II of Fig. 2 with ǫ < ǫ T (where the value of ǫ T i around 0.1 to 0.2 a etimated from our numerical olution), long range undulation of diffue grain boundarie are oberved, a hown in Fig. 4b. The tripe of domain B are curved around the boundarie, the latter moving toward each other at the expene of domain B until they finally merge. Depite the undulation, the grain boundary configuration in regime II i preerved during the proce. In regime III on the other hand, the boundary i harp, and the longitudinal ditortion lead to the formation of convex diclination in the boundary region (Fig. 4c for ǫ > ǫ T ). The boundary itelf alo move toward the bulk of domain B a indicated by the arrow in the figure, reulting again in the eventual diappearance of the grain boundary configuration. The longitudinal ditortion that we oberve in both regime, including the formation of an array of diclination in regime III, are due to wavenumber mimatch between the initial value and the tationary value qy <. In regime II (Fig. 4b), a periodic longitudinal modulation of tripe B appear, reminicent of an Eckhau ditortion that promote wavenumber reduction. Since the velocity of the boundary depend on the local wavenumber of tripe B ahead of it, a modulation of the boundary poition alo reult from thi wavenumber variation. In regime III, on the other hand, the boundary i harp, and cro roll, dilocation and diclination form at the edge of domain B leading alo to wavenumber reduction. The bulk of domain and B, however, remain largely unperturbed. In order to verify that thi pontaneou ditortion originate in the boundary region, we have tudied the decay of the configurationhown in Fig. 5a in which low tranveremodulation of wavenumber Q q x and phae magnitude δ havebeen added to the initial condition ofq x = qx in domain. In thi cae, = qy in the bulk ofdomain B but not at the grain boundarie. For large ǫ (regime III), the grain boundarie decay through the formation of diclination, a hown in Fig. 5b. The ame decay wa hown in Fig. 4c correponding to the cae qy in the whole domain B. imilar analyi can be carried out in regime II (mall ǫ), with the ame concluion. We note, however, an intereting apect of the decay a compared to potential model. The grain boundary move from the modulated domain toward the unditorted domain B, exactly the oppoite of what i expected for a potential model given our choice of wavenumber. With the conideration that wavenumber q x and in the bulk of both domain are equal to the tationary value qx and of the correponding potential model, the extra energy that i tored in the modulated tripe in domain would lead to boundary motion in the oppoite direction [7]. IV. WVEVECTOR SELECTION ND GRIN BOUNDRY DECY - MODEL 2 The econd non potential model tudied, a decribed in Sec. II, include the o called mean flow. They are known to have a ignificant effect on the tability of a tripe configuration [13], on defect dynamic [18, 19], and allow the appearance of piral defect chao [20]. We conider here the decay of the grain boundary configuration tudied above. We have choen model parameter value that are typical for fluid ytem in Eq. (1) (3): σ = 1, c 2 = 2, g m = 10,

6 6 ε < ε T ε > ε T B > q x q x = (a) Regime II (b) Regime III FIG. 6: Grain boundary decay for model 2 a a reult of mean flow, with g m = 10, σ = 1, c 2 = 2, and q x = q x = q 0. (a) Regime II: > q y and ǫ < ǫ T; the parameter ued in thi cae are = 1.02, ǫ = 0.04, and t = (b) Regime III: > q y and ǫ > ǫ T; the configuration hown correpond to = 1.02, ǫ = 0.3, and t = 600. and γ = λ 0 /2 for filtering. We have firt verified numerically that the unique wavenumber for a tationary grain boundary are qx = = q 0, the ame a in the potential Swift-Hohenberg model. Thi reult follow from the obervation that the lowet order contribution of mean flow to the amplitude equation of a grain boundary i O(ǫ 5/4 ) [18]. The known tability diagram of a tripe configuration i imilar to that hown in Fig. 2 for model 1, except that qy = q 0 = 1, and the type and detailed location of the tability boundarie are different (e.g., the Eckhau boundary in Fig. 2 hould be replaced by the kewed-varicoe boundary [13]). In term of our tudy of the decay of a grain boundary, we alo find in thi model three different regime which are qualitatively imilar to thoe of model 1. Rigid motion of the grain boundary i oberved when < qy for all ǫ, while a longitudinal ditortion of the tripe in domain B appear when > qy. We can alo qualitatively ditinguih between regime II (ǫ < ǫ T ) and III (ǫ > ǫ T ), with two typical example given in Fig. 6. wa the cae in model 1, we oberve long wavelength ditortion of the tripe in domain B in regime II (Fig. 6a), wherea regime III lead to the formation and ubequent motion of diclination, and with them, the decay of the grain boundary configuration (Fig. 6b). Diclination alo appear a a reult of local deviation of wavenumber from qy in the grain boundary region. Thi ha been verified by adding low tranvere modulation on tripe while keeping q = q in the bulk of both domain, with reult that are qualitatively imilar to thoe of Fig. 5. The driving force behind the appearance of longitudinal ditortion, and with them either long wavelength undulation or the formation of diclination i again local wavenumber reduction around the boundary. Thi i analogou to what i oberved in model 1, and ha been already decribed in detail in Sec. IIIB. However, once the undulation have been formed or the defect nucleated, mean flow are effective in their ubequent motion. Thi can be een from the patial ditribution of the vorticity potential ζ, or the correponding drift velocity U obtained in our numerical olution. In particular, we note that mean flow play a key role in the motion of diclination in regime III, a regime in which diclination become pinned in their abence (g m 0). V. DISCUSSION ND CONCLUSIONS The reult decribed on the decay of a grain boundary configuration in non potential model how the exitence of a range of parameter in which the decay lead to the formation of defect (motly dilocation and diclination) in the boundary region. Thi occur for a configuration in which the tripe parallel to the boundary plane have the elected wavenumber q x = qx, wherea the tripe perpendicular to the boundary are under light compreion > qy. lo, defect are only een above a certain threhold ǫ > ǫ T. Below the threhold, undulation and motion of the boundary follow, but no additional defect are produced. Thee reult can help interpret qualitative feature of domain coarening behavior, both in potential and non potential model. In the cae of non potential model 1 (with a ψ 2 2 ψ nonlinearity), the tationary wave number qx and qy are different from q 0, the mot untable wave number (ee Eq. (4) (6) and Fig. 2). Following the decay of an initially diordered configuration (ψ = 0), and the appearance of locally ordered domain with wavenumber approximately equal to q 0, grain boundarie would be expected to decay according to the mechanim decribed for regime II and III depending on the value of ǫ. Whether coarening take place in the regime in which the decay of boundarie lead to the formation of dilocation and diclination (regime III), or to their motion (regime II)

7 7 (a) ε < ε T (b) ε > ε T FIG. 7: Multidomain configuration for model 1 howing different qualitative coarening behavior for (a) ǫ < ǫ T (with ǫ = 0.04 and at t = 1000) and (b) ǫ > ǫ T (with ǫ = 0.5 and at t = 500). (a) ε < ε T (b) ε > ε T FIG. 8: Multidomain configuration of tripe phae in non potential model 2 (with mean flow), for (a) ǫ < ǫ T (ǫ = 0.04 and t = 1000) and (b) ǫ > ǫ T (ǫ = 0.5 and t = 1000). Other parameter are the ame a thoe in Fig. 6: g m = 10, σ = 1, and c 2 = 2. one might anticipate qualitatively different late time domain coarening. Figure 7 how the reult of a numerical integration of model 1 from a initial random configuration of ψ for two particular value of ǫ, both below and above ǫ T. For ǫ < ǫ T (Fig. 7a), grain boundarie are viible in the late time configuration a expected from the reult hown in Fig. 4b. Qualitatively different behavior i oberved for ǫ > ǫ T due to the early tage formation of dilocation or diclination around the emerging grain boundarie. Our reult in thi cae are hown in Fig. 7b, with the reulting diclination/dilocation dominated morphology. Thee reult are conitent with earlier obervation in Ref. 6, and more recent reearch on the Bouineq equation in Ref. 9. The average wave number of the evolving configuration change from q 0 at early time following the intability, toward a long time value that i cloer to the tationary wave number of iolated dilocation. One anticipate that domain coarening will be qualitatively different from the cae of mall ǫ, a uggeted by the reult of Ref. 6 howing a complex caling behavior at large ǫ. In the caeof model 2 which include the coupling to mean flow, the elected wavenumberq x and approximately coincide with the critical wavenumber for intability q 0, and therefore the effect decribed above will be weaker, but not completely abent a local deviation between the tripe wavenumber and the elected wavenumber till occur in region in which the tripe are curved. Near threhold (ǫ < ǫ T ), the morphology of a tranient tripe configuration i hown in Fig. 8a obtained by numerical olution of the governing equation from a random initial condition. tripe pattern qualitatively imilar to both model 1 and the potential Swift-Hohenberg model i found, characterized by differently oriented domain eparated by 90 grain boundarie. On the other hand, when ǫ exceed a certain threhold ǫ T, grain boundarie are expected to be abent a decribed in Sec. IV. The ubequent evolution of the ytem in thi cae of i much more complicated reulting in peritent dynamic (piral defect chao [20]), a hown in Fig. 8b.

8 In both non potential cae dicued in thi paper we would expect qualitatively different domain coarening behavior depending on the interplay between the critical wavenumber for intability and the elected wavenumber at a grain boundary. Thi i in contrat with the potential Swift-Hohenberg model in which grain boundary dynamic i known to dominate for mall ǫ until pinning effect at finite ǫ effectively arret motion [8, 12]. In thi model, the elected wavenumber for a tationary grain boundary configuration are q x = q y = q 0 [16], with q = q 0 alo marking the tability boundary of the bulk zig-zag intability. Therefore a diagram analogou to that of Fig. 2 would only have regime II and III characterizing grain boundary motion. In both > q y, longitudinal ditortion around the boundary can alo be identified a a reult of local wavenumber reduction toward q y = q 0. In analogy with the non potential cae, in regime II (with ǫ < ǫ T ) diffue boundarie with long range undulation of B tripe are oberved. For ǫ > ǫ T, the formation of convex diclination around the boundary can alo be een in numerical calculation. However, depite their emergence thee diclination become pinned along with the grain boundarie, mainly due to the nonadiabatic effect ariing from the coupling between fat cale of bae tripe pattern and low cale of the aociated envelope [12]. Thi reult in a glay configuration with very low dynamic. Interetingly, thi pinning phenomenon i only oberved in the potential Swift-Hohenberg model, but not in the two non potential model tudied here. In thi latter cae, grain boundarie quickly decay at large ǫ leaving behind a domain morphology dominated by diclination and dilocation (Fig. 7b), or piral defect (Fig. 8b). In ummary, we have tudied wavenumber election, motion, and decay of grain boundarie in a tripe configuration baed on a pair of two-dimenional generalized Swift-Hohenberg model equation that are non potential in nature. Three different regime of grain boundary dynamic have been identified, with longitudinal ditortion emerging around the grain boundary and manifeting themelve either a long range undulation of the moving boundary or in the formation and ubequent motion of defect uch a convex diclination. When an extended ytem i taken from an initial tate below threhold to ǫ > 0, the initial decay of the untable tate will lead to the formation of locally ordered domain with a ditribution of orientation. Domain interface will form in accordance with the analyi given above, and therefore qualitative difference in domain coarening are expected depending on the value of ǫ. 8 cknowledgment Thi reearch ha been upported by the National Science Foundation under grant DMR , and by NSERC Canada. [1] M. Seul and D. ndelman, Science 267, 476 (1995). [2] F. S. Bate and G. H. Fredrickon, Phy. Today 52, 32 (1999). [3] M. C. Cro and P. C. Hohenberg, Rev. Mod. Phy. 65, 851 (1993). [4] L. Purvi and M. Dennin, Phy. Rev. Lett. 86, 5898 (2001). [5] K. Elder, J. Viñal, and M. Grant, Phy. Rev. Lett. 68, 3024 (1992). [6] M. C. Cro and D. I. Meiron, Phy. Rev. Lett. 75, 2152 (1995). [7] D. Boyer and J. Viñal, Phy. Rev. E 63, (2001). [8] D. Boyer and J. Viñal, Phy. Rev. E 64, (R) (2001). [9] M. R. Paul, K.-H. Chiam, M. C. Cro, and P. F. Ficher, Phy. Rev. Lett. 93, (2004). [10] C. Kamaga, F. Ibrahim, and M. Dennin, Phy. Rev. E 69, (2004). [11] C. Harrion, D.. damon, Z. Cheng, J. M. Sebatian, S. Sethuraman, D.. Hue, R.. Regiter, and P. M. Chaikin, Science 290, 1558 (2000). [12] D. Boyer and J. Viñal, Phy. Rev. E 65, (2002). [13] H. S. Greenide and M. C. Cro, Phy. Rev. 31, 2492 (1985). [14] P. Manneville, J. Phyique 44, 759 (1983). [15] J. Swift and P. C. Hohenberg, Phy. Rev. 15, 319 (1977). [16] G. Teauro and M. C. Cro, Phil. Mag. 56, 703 (1987). [17] M. C. Cro, D. I. Meiron, and Y. Tu, Chao 4, 607 (1994). [18] G. Teauro and M. C. Cro, Phy. Rev. 34, 1363 (1986). [19] T. Walter, W. Pech, and E. Bodenchatz, Chao 14, 933 (2004). [20] H.-W. Xi, J. D. Gunton, and J. Viñal, Phy. Rev. Lett. 71, 2030 (1993).