Coarsening dynamics of a one-dimensional driven Cahn-Hilliard system

Steady states and coarsening in one-dimensional driven Allen-Cahn system

We study the steady states and the coarsening dynamics in a one-dimensional driven nonconserved system modeled by the so-called driven Allen-Cahn equation, which is the standard Allen-Cahn equation with an additional driving force. In particular, we derive equations of motion for the phase boundaries in a phase ordering system obeying this equation using a nearest-neighbor interaction approach. Using the equations of motion we explore kink binary and ternary interactions and analyze how the average domain size scale with respect to time. Further, we employ numerical techniques to perform a bifurcation analysis of the one-period stationary solutions of the equation. We then investigate the linear stability of the two-period solutions and thereby identify and study various coarsening modes.

Suppression of coarsening and emergence of oscillatory behavior in a Cahn-Hilliard model with nonvariational coupling

We investigate a generic two-field Cahn-Hilliard model with variational and nonvariational coupling. It describes, for instance, passive and active ternary mixtures, respectively. Already a linear stability analysis of the homogeneous mixed state shows that activity not only allows for the usual large-scale stationary (Cahn-Hilliard) instability of the well-known passive case but also for small-scale stationary (Turing) and large-scale oscillatory (Hopf) instabilities. In consequence of the Turing instability, activity may completely suppress the usual coarsening dynamics. In a fully nonlinear analysis, we first briefly discuss the passive case before focusing on the active case. Bifurcation diagrams and selected direct time simulations are presented that allow us to establish that nonvariational coupling (i) can partially or completely suppress coarsening and (ii) may lead to the emergence of drifting and oscillatory states. Throughout, we emphasize the relevance of conservation laws and related symmetries for the encountered intricate bifurcation behavior.

Unstable wave-packet perturbations in an advective Cahn-Hilliard process

A study of the perturbation dynamics in a one-dimensional advective Cahn-Hilliard system, characterized by a nonvanishing driving force, is carried out to test the stability of a uniform basic solution. The linear stability of small-amplitude perturbations is analyzed both in the case of normal Fourier modes, with a given wave number, and in the case of wave packets localized in space. The dual nature of the instability, either of convective or absolute type, is studied, revealing that the driving force creates a gap between the parametric threshold to instability of normal modes and that to instability of wave packets. When the driving force is zero, also the gap between such thresholds disappears.

Thixotropy and shear thinning of lubricated contacts with confined membranes

We have modeled the nonlinear dynamics and the rheological behavior of a system under shear containing a membrane confined between two attractive walls. The presence of the membrane induces additional tangential forces on the walls that always increase the global friction. At low shear rates, the membrane exhibits chaotic dynamics with slow coarsening leading to thixotropy, i.e. to a slow decrease of the membrane-induced tangential forces on the walls. At intermediate shear rates, the membrane profile presents stationary periodic patterns. At higher shear rates, membrane dynamics are governed by a nonlinear evolution equation which is similar to the Kuramoto-Sivashinski equation, but with a sixth-order stabilizing term. The membrane experiences chaotic dynamics without coarsening. As a consequence of the nonlinear dynamics of the membrane at intermediate and large shear rates, the system exhibits shear thinning.

Mean-field theory for coarsening faceted surfaces

A mean-field theory is developed for the scale-invariant length distributions observed during the coarsening of generic one-dimensional faceted surfaces. This theory closely follows the Lifshitz-Slyozov-Wagner theory of Ostwald ripening in two-phase systems, but the mechanism of coarsening in faceted surfaces requires the addition of convolution terms recalling work on particle coalescence, and induces an unexpected coupling between the convolution and the rate of facet loss. As a generic framework, the theory concisely illustrates how the universal processes of facet removal and neighbor merger are moderated by the system-specific mean-field velocity describing average rates of length change. For a simple, example facet dynamics associated with the directional solidification of a binary alloy, agreement between the predicted scaling state and that observed after direct numerical simulation of 40,000,000 facets is reasonable given the limiting assumption of noncorrelation between neighbors; relaxing this assumption is a clear path forward toward improved quantitative agreement with data in the future.

Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth

In the present contribution we review basic mathematical results for three physical systems involving self-organizing solid or liquid films at solid surfaces. The films may undergo a structuring process by dewetting, evaporation/condensation or epitaxial growth, respectively. We highlight similarities and differences of the three systems based on the observation that in certain limits all of them may be described using models of similar form, i.e. time evolution equations for the film thickness profile. Those equations represent gradient dynamics characterized by mobility functions and an underlying energy functional. Two basic steps of mathematical analysis are used to compare the different systems. First, we discuss the linear stability of homogeneous steady states, i.e. flat films, and second the systematics of non-trivial steady states, i.e. drop/hole states for dewetting films and quantum-dot states in epitaxial growth, respectively. Our aim is to illustrate that the underlying solution structure might be very complex as in the case of epitaxial growth but can be better understood when comparing the much simpler results for the dewetting liquid film. We furthermore show that the numerical continuation techniques employed can shed some light on this structure in a more convenient way than time-stepping methods. Finally we discuss that the usage of the employed general formulation does not only relate seemingly unrelated physical systems mathematically, but does allow as well for discussing model extensions in a more unified way.

Faceting and coarsening dynamics in the complex Swift-Hohenberg equation

The complex Swift-Hohenberg equation models pattern formation arising from an oscillatory instability with a finite wave number at onset and finds applications in lasers, optical parametric oscillators, and photorefractive oscillators. We show that with real coefficients this equation exhibits two classes of localized states: localized in amplitude only or localized in both amplitude and phase. The latter are associated with phase-winding states in which the real and imaginary parts of the order parameter oscillate periodically but with a constant phase difference between them. The localized states take the form of defects connecting phase-winding states with equal and opposite phase lag, and can be stable over a wide range of parameters. The formation of these defects leads to faceting of states with initially spatially uniform phase. Depending on parameters these facets may either coarsen indefinitely, as described by a Cahn-Hilliard equation, or the coarsening ceases leading to a frozen faceted structure.

Formation of regular structures in the process of phase separation

Phase separation under directional quenching has been studied in a Cahn-Hilliard model. In distinct contrast to the disordered patterns which develop under a homogeneous quench, periodic stripe patterns are generated behind the quench front. Their wavelength is uniquely defined by the velocity of the quench interface in a wide range. Numerical simulations match perfectly analytical results obtained in the limit of small and large velocities of the quench interface. Additional periodic modulation of the quench interface may lead to cellular patterns. The quenching protocols analyzed are expected to be an effective tool in technological applications to design nanostructured materials.

Dynamic phase separation: from coarsening to turbulence via structure formation

We investigate some new two-dimensional evolution models belonging to the class of convective Cahn-Hilliard models: (i) a local model with a scalar order parameter, (ii) a nonlocal model with a scalar order parameter, and (iii) a model with a vector order parameter. These models are applicable to phase-separating system where concentration gradients cause hydrodynamic motion due to buoyancy or Marangoni effect. The numerical study of the models shows transition from coarsening, typical of Cahn-Hilliard systems, to spatiotemporally irregular behavior (turbulence), typical of the Kuramoto-Sivashinsky equation, which is obtained in the limit of very strong driving. The transition occurs not in a straightforward way, but through the formation of spatial patterns that emerge for intermediate values of the driving intensity. As in driven one-dimensional models studied before, the mere presence of the driving force, however small, breaks the symmetry between the two separating phases, as well as increases the coarsening rate. With increasing driving, coarsening stops. The dynamics is generally irregular at strong driving, but exhibits specific structural features.

Faceting of a growing crystal surface by surface diffusion

Consider faceting of a crystal surface caused by strongly anisotropic surface tension, driven by surface diffusion and accompanied by deposition (etching) due to fluxes normal to the surface. Nonlinear evolution equations describing the faceting of 1+1 and 2+1 crystal surfaces are studied analytically, by means of matched asymptotic expansions for small growth rates, and numerically otherwise. Stationary shapes and dynamics of faceted pyramidal structures are found as functions of the growth rate. In the 1+1 case it is shown that a solitary hill as well as periodic hill-and-valley solutions are unique, while solutions in the form of a solitary valley form a one-parameter family. It is found that with the increase of the growth rate, the faceting dynamics exhibits transitions from the power-law coarsening to the formation of pyramidal structures with a fixed average size and finally to spatiotemporally chaotic surfaces resembling the kinetic roughening.

Convective Cahn-Hilliard models: from coarsening to roughening

In this paper we demonstrate that convective Cahn-Hilliard models, describing phase separation of driven systems (e.g., faceting of growing thermodynamically unstable crystal surfaces), exhibit, with the increase of the driving force, a transition from the usual coarsening regime to a chaotic behavior without coarsening via a pattern-forming state characterized by the formation of various stationary and traveling periodic structures as well as structures with localized oscillations. Relation of this phenomenon to a kinetic roughening of thermodynamically unstable surfaces is discussed.

Ripple formation through an interface instability from moving growth and erosion sources

The propagation of material interfaces is investigated under the action of a localized moving source which deposits or removes material. Among others the latter process applies to beam cutting techniques. We develop a Kuramoto-Sivashinsky-type model and find a new type of ripple forming mechanism. This theory offers a new explanation for the occurrence of striation patterns which often degrade the quality of cutting edges.

Fast coarsening in unstable epitaxy with desorption

Homoepitaxial growth is unstable towards the formation of pyramidal mounds when interlayer transport is reduced due to activation barriers to hopping at step edges. Simulations of a lattice model and a continuum equation show that a small amount of desorption dramatically speeds up the coarsening of the mound array, leading to coarsening exponents between 1/3 and 1/2. The underlying mechanism is the faster growth of larger mounds due to their lower evaporation rate.