Defect structures in the growth kinetics of the Swift-Hohenberg model

Field-Theoretic Study of Salt-Induced Order and Disorder in a Polarizable Diblock Copolymer

We study a salt-doped polarizable symmetric diblock copolymer using a recently developed field theory that self-consistently embeds dielectric response, ion solvation energies, and van der Waals (vdW) attractions via the incorporation of segment polarizabilities and fixed dipoles. This field theory is amenable to direct simulation via the complex Langevin sampling technique and, thus, requires no approximations beyond the phenomenology of the underlying molecular model. We measure the shift in the order-disorder transition (ODT) of a diblock copolymer with salt-loading in field-theoretic simulations and observe rich behavior in which solvation, dilution and charge screening effects compete to determine whether the ordered or disordered phase is stabilized. At low salt concentrations, the salt behaves as a selective solvent, localizing into the high-dielectric domains and stabilizing the ordered phase. At high salt concentrations, however, the salt localization vanishes due to charge screening effects, and the salt behaves as a nonselective solvent that screens vdW attractions and stabilizes the disordered phase.

Coarsening dynamics in the Swift-Hohenberg equation with an external field

We study the Swift-Hohenberg equation (SHE) in the presence of an external field. The application of the field leads to a phase diagram with three phases, i.e., stripe, hexagon, and uniform. We focus on coarsening after a quench from the uniform to stripe or hexagon regions. For stripe patterns, we find that the length scale associated with the order-parameter structure factor has the same growth exponent (≃1/4) as for the SHE with zero field. The growth process is slower in the case of hexagonal patterns, with the effective growth exponent varying between 1/6 and 1/9, depending on the quench parameters. For deep quenches in the hexagonal phase, the growth process stops at late stages when defect boundaries become pinned.

Scaling properties of ageing orientation fluctuations in stripe phases

We investigate the non-equilibrium dynamics of an ordered stripe-forming system free of topological defects. In particular, we study the ageing and the coarsening of orientation fluctuations parallel and perpendicular to the stripes via computer simulations based on a minimal phase-field model (model B with Coulomb interactions). Under the influence of noise, the stripe orientation field develops fluctuations parallel to the stripes, with the dominant modulation length λ*∥ increasing with time t as λ*∥ ∼ t1/4 and the correlation length perpendicular to the stripes ξ⊥θ increasing as ξ⊥θ ∼ t1/2. We explain these anisotropic coarsening dynamics with an analytic theory based on the linear elastic model for stripe displacements first introduced by Landau and Peierls. We thus obtain the scaling forms and the scaling exponents characterizing the correlation functions and the structure factor of the stripe orientation field. Our results reveal how the coarsening of orientation fluctuations prevents a periodically modulated phase free of topological defects from reaching equilibrium.

Order-disorder transition in thin films of horizontally-oriented cylinder-forming block copolymers: thermal fluctuations vs. preferential wetting

We present experimental and theoretical investigations of the order-disorder transition (ODT) in thin films of cylinder-forming diblock copolymers with asymmetric wetting conditions. Grazing incidence small-angle X-ray scattering (GISAXS) was implemented to determine the ODT temperatures (TODT) for poly(styrene-b-2-vinyl pyridine) (PS-P2VP) block copolymer thin films on a P2VP-preferential silicon substrate. Specifically, films consisting of multilayers of horizontally-oriented cylindrical structures (from 1- to 9-layers) were tested. We find that films with more than 2 cylindrical layers have a TODT comparable to the bulk case. However, TODT decreases as the film becomes thinner and the monolayer system has an ODT 30 °C below the bulk. Using self-consistent field theory (SCFT), we studied the ordering in corresponding thin films with asymmetric (top and bottom surface) wetting conditions. Surprisingly, SCFT is found to predict an opposite trend in TODT with film thickness than observed experimentally. Field-theoretic simulations with complex Langevin sampling were employed to resolve this discrepancy and demonstrate that thermal fluctuations in the PS-P2VP thin-film system dominate its ordering behavior in monolayer and bilayer films.

Coarsening of stripe patterns: variations with quench depth and scaling

The coarsening of stripe patterns when the system is evolved from random initial states is studied by varying the quench depth ε, which is a measure of distance from the transition point of the stripe phase. The dynamics of the growth of stripe order, which is characterized by two length scales, depends on the quench depth. The growth exponents of the two length scales vary continuously with ε. The decay exponents for free energy, stripe curvature, and densities of defects like grain boundaries and dislocations also show similar variation. This implies a breakdown of the standard picture of nonequilibrium dynamical scaling. In order to understand the variations with ε we propose an additional scaling with a length scale dependent on ε. The main contribution to this length scale comes from the "pinning potential," which is unique to systems where the order parameter is spatially periodic. The periodic order parameter gives rise to an ε-dependent potential, which can pin defects like grain boundaries, dislocations, etc. This additional scaling provides a compact description of variations of growth exponents with quench depth in terms of just one exponent for each of the length scales. The relaxation of free energy, stripe curvature, and the defect densities have also been related to these length scales. The study is done at zero temperature using Swift-Hohenberg equation in two dimensions.

Aging of orientation fluctuations in stripe phases

Stripe patterns, observed in a large variety of physical systems, often exhibit a slow nonequilibrium dynamics because ordering is impeded by the presence of topological defects. Using computer simulations based on a well-established model for stripe formation, we show that a slow dynamics and aging occur also in stripe patterns free of topological defects. For a wide range of noise strengths, the two-time orientation correlation function follows a scaling form that is typical for systems exhibiting a growing length scale. In our case, the underlying mechanism is the coarsening of orientation fluctuations, ultimately leading to power-law spatial correlations perpendicular to the stripes. Our results show that even for the smallest amount of noise, stripe phases without topological defects do not reach equilibrium. This constitutes an important aspect of the dynamics of modulated phases.

Defect structures and ordering behaviours of diblock copolymers self-assembling on spherical substrates

One of the main differences of ordered structures constrained on curved surfaces is the nature of topological defects. We here explore the defect structures and ordering behaviours of both lamellar and cylindrical phases of block copolymers confined on spherical substrates by the Landau-Brazovskii theory, which is numerically solved by a highly accurate spectral method with a spherical harmonic basis. For the cylindrical phase, isolated disclinations and scars are generated on the spherical substrates. The number of excess dislocations in a scar depends linearly on the sphere radius. The defect fraction characterizing the ordering dynamics decays exponentially. The scars are formed from the isolated disclinations via mini-scars. For the lamellar phase, three types of defect structures (hedgehog, spiral and quasi-baseball) are identified. The disclination annihilation is the primary ordering mechanism of the lamellar phase.

Stripe patterns: role of initial state and boundary conditions

This paper presents results on stripe patterns by numerical solution of the Swift-Hohenberg equation. The focus is on the role of initial state and boundary conditions. We choose initial states which generate simple defect configurations and study their evolution. Various classes of defects are identified and their motion and relaxation is studied numerically. We first study the dynamics of a straight front and present a comparison of numerical results with some analytical results. We then study the domain-wall dynamics in configurations containing two and three domains and identify some mechanisms of their relaxation. Rates of domain-wall relaxation depend on several features like incommensuration, dislocations and orientations in neighboring domains, in addition to the curvature of the walls. For a generic class of domain walls the relaxation process has an intrinsic frustration which leads to generation of dislocations. This process also generates stripe curvature thereby making relaxation nonmonotonic. We have also generated some other topological defects and studied their evolution and the effect of boundary conditions on their stability.

Kinetics of lamellar formation on sparsely stripped patterns

Chemical epitaxy based on the self-assembly of block copolymers is viewed as a promising technique to achieve ordered patterns on a large scale. Herein, we study the kinetics of lamellar formation of block copolymers under the direction of sparsely stripped patterns using cell dynamics simulations of the time-dependent Ginzburg-Landau theory. First, a scaling law is unveiled with the ordering time of lamellae, tp, with respect to the multiples between the periods of lamellae and stripe patterns, which is consistent with the power law evolution of the correlation length existing in the bulk phase of lamellae. Second, the tolerative windows of perfect order, with deviation from integer multiples, are also estimated from the aspect of kinetics. The results of the ordering time and tolerative windows are of great interest for relevant experiments or applications. Finally, a two-stage evolution is explored during the pattern formation of chemical epitaxy by probing into the evolution of defects, which is of fundamental interest for us to understand the coarsening kinetics of block copolymers under the direction of chemical patterns.

Nucleation, growth, and coarsening of crystalline domains in order-order transitions between lamellar and hexagonal phases

Using the numerical solution of the time-dependent Ginzburg-Landau equation, we study the entire process of transformation between the lamellar and the hexagonal phases from the early-stage nucleation and growth to the late-stage coarsening regime. The metastable crystalline structure that nucleates first is identified in terms of the mean-field theory under the single-wave-number approximation. This has been borne out by the numerically efficient preparation of single-crystal structure developed via the noise-induced self-organization. We also present results for the scaling of the late-time domain growth, which is quantified by two measures: the structure factor and the orientational correlation function. In particular, the growth exponent is shown to be robust and indifferent to conservation of the order parameter.

Shear alignment of a disordered lamellar mesophase

The shear alignment of an initially disordered lamellar phase is examined using lattice Boltzmann simulations of a mesoscopic model based on a free-energy functional for the concentration modulation. For a small shear cell of width 8λ, the qualitative features of the alignment process are strongly dependent on the Schmidt number Sc=ν/D (ratio of kinematic viscosity and mass diffusion coefficient). Here, λ is the wavelength of the concentration modulation. At low Schmidt number, it is found that there is a significant initial increase in the viscosity, coinciding with the alignment of layers along the extensional axis, followed by a decrease at long times due to the alignment along the flow direction. At high Schmidt number, alignment takes place due to the breakage and reformation of layers because diffusion is slow compared to shear deformation; this results in faster alignment. The system size has a strong effect on the alignment process; perfect alignment takes place for a small systems of width 8λ and 16λ, while a larger system of width 32λ does not align completely even at long times. In the larger system, there appears to be a dynamical steady state in which the layers are not perfectly aligned--where there is a balance between the annealing of defects due to shear and the creation due to an instability of the aligned lamellar phase under shear. We observe two types of defect creation mechanisms: the buckling instability under dilation, which was reported earlier, as well as a second mechanism due to layer compression.

Novel mechanical properties in lamellar phases of liquid-crystalline diblock copolymers

Structural properties of flexible nematic diblock copolymers in the lamellar phase are investigated using a mean-field model. We address two complementary questions on the mechanics of the system: 1) How does the nematic order affect the elasticity of the one-dimensional solid? 2) What effect does the block copolymer microstructure has on the orientation of the nematic director? In the limit when the microstructure does not influence the nematic director orientation we predict a soft lamellar compression mode. When the microstructure does influence the nematic director orientation, small compressions lead to conventional elasticity, until a critical strain is reached, where there is then a transition to a softer response. On the other hand, we show that an identifiable lamellar symmetry provides a direction along which the nematic director prefers to align. Our model provides avenues to explore nonlinear properties of flexible diblock copolymers in which the monomers on both sides have mesogenic side groups.

Coarsening in potential and nonpotential models of oblique stripe patterns

We study the coarsening of two-dimensional oblique stripe patterns by numerically solving potential and nonpotential anisotropic Swift-Hohenberg equations. Close to onset, all models exhibit isotropic coarsening with a single characteristic length scale growing in time as t1/2. Further from onset, the characteristic lengths along the preferred directions x and ŷ grow with different exponents, close to 1/3 and 1/2, respectively. In this regime, one-dimensional dynamical scaling relations hold. We draw an analogy between this problem and model A in a stationary, modulated external field. For deep quenches, nonpotential effects produce a complicated dislocation dynamics that can lead to either arrested or faster-than-power-law growth, depending on the model considered. In the arrested case, small isolated domains shrink down to a finite size and fail to disappear. A comparison with available experimental results for electroconvection in nematics is presented.

Orientational order in block copolymer films zone annealed below the order--disorder transition temperature

We report measurements of rapid ordering and preferential alignment in block copolymer films zone annealed below the order-disorder transition temperature. The orientational correlation lengths measured after approximately 5 h above the glass-transition temperature ( approximately 2 microm) were an order of magnitude greater than that obtained under equivalent static annealing. The ability to rapidly process polymers with inaccessible order-disorder transition temperatures suggests zone annealing as a route toward more robust nanomanufacturing methods based on block copolymer self-assembly.

Impact of noise on domain growth in electroconvection

The growth and ordering of striped domains has recently received renewed attention due in part to experimental studies in diblock copolymers and electroconvection. One surprising result has been the relatively slow dynamics associated with the growth of striped domains. One potential source of the slow dynamics is the pinning of defects in the periodic potential of the stripes. Of interest is whether or not external noise will have a significant impact on the domain ordering, perhaps by reducing the pinning and increasing the rate of ordering. In contrast, we present experiments using electroconvection in which we show that a particular type of external noise decreases the rate of domain ordering.

Morphologies and flow patterns in quenching of lamellar systems with shear

We study the behavior of a fluid quenched from the disordered into the lamellar phase under the action of a shear flow. The dynamics of the system is described by Navier-Stokes and convection-diffusion equations with the pressure tensor and the chemical potential derived by the Brazovskii free energy. Our simulations are based on a mixed numerical method with the lattice Boltzmann equation and a finite difference scheme for Navier-Stokes and order parameter equations, respectively. We focus on cases where banded flows are observed with two different slopes for the component of velocity in the direction of the applied flow. Close to the walls the system reaches a lamellar order with very few defects, and the slope of the horizontal velocity is higher than the imposed shear rate. In the middle of the system the local shear rate is lower than the imposed one, and the system looks like a mixture of tilted lamellae, droplets, and small elongated domains. We refer to this as a region with a shear-induced structures (SIS) configuration. The local behavior of the stress shows that the system with the coexisting lamellar and SIS regions is in mechanical equilibrium. This phenomenon occurs, at fixed viscosity, for shear rates under a certain threshold; when the imposed shear rate is sufficiently large, lamellar order develops in the whole system. Effects of different viscosities have been also considered. The SIS region is observed only at low enough viscosity. We compare the above scenario with the usual one of shear banding. In particular, we do not find evidence for a plateau of the stress at varying imposed shear rates in the region with banded flow. We interpret our results as due to a tendency of the lamellar system to oppose the presence of the applied flow.

Early time evolution of Freédericksz patterns generated from states of electroconvection

We report on the early time ordering in a nematic liquid crystal subjected to a sudden change in an external ac electric field. We compare time evolution for two different initial states of electroconvection. Electroconvection is a highly driven state of a nematic liquid crystal involving convective motion of the fluid and periodic variations of the molecular alignment. By suddenly changing either the voltage or the frequency of the applied ac field, the system is brought to the same thermodynamic conditions. The time ordering of the system is characterized by the evolution of features of the power spectrum, including the average wave number, total power, and shape of the power spectrum. We observe that ordering of the system occurs faster after a sudden change in frequency than it does after a sudden change in voltage.

Coarsening of topological defects in oscillating systems with quenched disorder

We use large scale simulations to study interacting particles in two dimensions in the presence of both an ac drive and quenched disorder. As a function of ac amplitude, there is a crossover from a low drive regime where the colloid positions are highly disordered to a higher ac drive regime where the system dynamically reorders. We examine the coarsening of topological defects formed when the system is quenched from a disordered low ac amplitude state to a high ac amplitude state. When the quench is performed close to the disorder-order crossover, the defect density decays with time as a power law with alpha = 1/4 to 1/3. For deep quenches, in which the ac drive is increased to high values such that the dynamical shaking temperature is strongly reduced, we observe a logarithmic decay of the defect density into a grain boundary dominated state. We find a similar logarithmic decay of defect density in systems containing no pinning. We specifically demonstrate these effects for vortices in thin film superconductors, and discuss implications for dynamical reordering transition studies in these systems.

Growth of order in an anisotropic Swift-Hohenberg model

We have studied the ordering kinetics of a two-dimensional anisotropic Swift-Hohenberg (SH) model numerically. The defect structure for this model is simpler than for the isotropic SH model. One finds only dislocations in the aligned ordering striped system. The motion of these point defects is strongly influenced by the anisotropic nature of the system. We developed accurate numerical methods for following the trajectories of dislocations. This allows us to carry out a detailed statistical analysis of the dynamics of the dislocations. The average speeds for the motion of the dislocations in the two orthogonal directions obey power laws in time with different amplitudes but the same exponents. The position and velocity distribution functions are only weakly anisotropic.

Numerical simulations of two-dimensional magnetic domain patterns

I show that a model for the interaction of magnetic domains that includes a short range ferromagnetic and a long range dipolar antiferromagnetic interaction reproduces very well many characteristic features of two-dimensional magnetic domain patterns. In particular bubble and stripe phases are obtained, along with polygonal and labyrinthine morphologies. In addition, two puzzling phenomena, namely, the so called "memory effect" and the "topological melting" observed experimentally, are also qualitatively described. Very similar phenomenology is found in the case in which the model is changed to be represented by the Swift-Hohenberg equation driven by an external orienting field.

Numerical study of domain coarsening in anisotropic stripe patterns

We study the coarsening of two-dimensional smectic polycrystals characterized by grains of oblique stripes with only two possible orientations. For this purpose, an anisotropic Swift-Hohenberg equation is solved. For quenches close enough to the onset of stripe formation, the average domain size increases with time as t(1/2). Further from onset, anisotropic pinning forces similar to Peierls stresses in solid crystals slow down defects, and growth becomes anisotropic. In a wide range of quench depths, dislocation arrays remain mobile and dislocation density roughly decays as t(-1/3), while chevron boundaries are totally pinned. We discuss some agreements and disagreements found with recent experimental results on the coarsening of anisotropic electroconvection patterns.

Dislocation dynamics in an anisotropic stripe pattern

The dynamics of dislocations confined to grain boundaries in a striped system are studied using electroconvection in the nematic liquid crystal N4. In electroconvection, a striped pattern of convection rolls forms for sufficiently high driving voltages. We consider the case of a rapid change in the voltage that takes the system from a uniform state to a state consisting of striped domains with two different wave vectors. The domains are separated by domain walls along one axis and a grain boundary of dislocations in the perpendicular direction. The pattern evolves through dislocation motion parallel to the domain walls. We report on features of the dislocation dynamics. The kinetics of the domain motion is quantified using three measures: dislocation density, average domain wall length, and total domain wall length per area. All three quantities exhibit behavior consistent with power-law evolution in time, with the defect density decaying as t(-1/3), the average domain wall length growing as t(1/3), and the total domain wall length decaying as t(-1/5). The two different exponents are indicative of the anisotropic growth of domains in the system.

Dependence of domain wall dynamics on background wave number

We report on the growth of domains of standing waves in electroconvection in a nematic liquid crystal, focusing on the evolution of domain walls. An ac voltage is applied to the system, forming an initial state that consists of traveling striped patterns with two different orientations, zig and zag rolls. The standing waves are generated by suddenly applying a periodic modulation of the amplitude of the applied voltage that is approximately resonant with the traveling frequency of the pattern. By varying the modulation frequency, we are able to vary the steady-state, average wave number. We characterize the evolution of the domain walls as a function of the average background wave number by measuring the total area and length of domain walls present in the system as a function of time. We find that as the background wave number is varied away from the "natural" wave number for the pattern, the evolution of the domain walls occurs at a faster rate.

Model for striped growth

We introduce a model for describing the defected growth of striped patterns. This model, while roughly related to the Swift-Hohenberg model, generates a quite different mixture of defects during phase ordering. We find two characteristic lengths in the system: the scaling length L(t), and the average width of the domain walls. The growth law exponent is larger than the value of 1/2 found in typical point defect systems.