Phase ordering with a global conservation law: Ostwald ripening and coalescence

Coarsening and wavelength selection far from equilibrium: A unifying framework based on singular perturbation theory

Intracellular protein patterns are described by (nearly) mass-conserving reaction-diffusion systems. While these patterns initially form out of a homogeneous steady state due to the well-understood Turing instability, no general theory exists for the dynamics of fully nonlinear patterns. We develop a unifying theory for nonlinear wavelength-selection dynamics in (nearly) mass-conserving two-component reaction-diffusion systems independent of the specific mathematical model chosen. Previous work has shown that these systems support an extremely broad band of stable wavelengths, but the mechanism by which a specific wavelength is selected has remained unclear. We show that an interrupted coarsening process selects the wavelength at the threshold to stability. Based on the physical intuition that coarsening is driven by competition for mass and interrupted by weak source terms that break strict mass conservation, we develop a singular perturbation theory for the stability of stationary patterns. The resulting closed-form analytical expressions enable us to quantitatively predict the coarsening dynamics and the final pattern wavelength. We find excellent agreement with numerical results throughout the diffusion- and reaction-limited regimes of the dynamics, including the crossover region. Further, we show how, in these limits, the two-component reaction-diffusion systems map to generalized Cahn-Hilliard and conserved Allen-Cahn dynamics, therefore providing a link to these two fundamental scalar field theories. The systematic understanding of the length-scale dynamics of fully nonlinear patterns in two-component systems provided here builds the basis to reveal the mechanisms underlying wavelength selection in multicomponent systems with potentially several conservation laws.

Spatiotemporal Mapping of Hole Nucleation and Growth during Block Copolymer Terracing with High-Speed Atomic Force Microscopy

As a platform for investigating two-dimensional phase separation, we track the structural evolution of block copolymer thin films during thermal annealing with environmentally controlled atomic force microscopy (AFM). Upon thermal annealing, block copolymer films with incommensurate thickness separate into a terraced morphology decorated with holes. With in situ imaging at 200 °C, we follow the continuous progression of terrace formation in a single region of a cylinder-forming poly(styrene-block-methyl methacrylate) thin film, beginning with the disordered morphology on an unpatterned silicon substrate and continuing through nucleation and coarsening stages. Topographic AFM imaging with nanoscale resolution simultaneously captures ensemble hole growth statistics while locally tracking polymer diffusion through measurements of the film thickness. At early times, we observe homogeneous hole nucleation and isotropic growth, with kinetics following the predictions of classical nucleation theory. At later times, however, we find anomalous hole growth which arises due to the combination of Ostwald ripening and coalescence mechanisms. In each case, our real-space observations highlight the importance of hole interactions for determining coarsening kinetics, mediated either through the interconnected phase for Ostwald ripening or through binary collision events for coalescence.

The influence of non-stationarity and interphase curvature on the growth dynamics of spherical crystals in a metastable liquid

This manuscript is concerned with the theory of nucleation and evolution of a polydisperse ensemble of crystals in metastable liquids during the intermediate stage of a phase transformation process. A generalized growth rate of individual crystals is obtained with allowance for the effects of their non-stationary evolution in unsteady temperature (solute concentration) field and the phase transition temperature shift appearing due to the particle curvature (the Gibbs-Thomson effect) and atomic kinetics. A complete system of balance and kinetic equations determining the transient behaviour of the metastability degree and the particle-radius distribution function is analytically solved in a parametric form. The coefficient of mutual Brownian diffusion in the Fokker-Planck equation is considered in a generalized form defined by an Einstein relation. It is shown that the effects under consideration substantially change the desupercooling/desupersaturation dynamics and the transient behaviour of the particle-size distribution function. The asymptotic state of the distribution function (its 'tail'), which determines the relaxation dynamics of the concluding (Ostwald ripening) stage of a phase transformation process, is derived. This article is part of the theme issue 'Transport phenomena in complex systems (part 1)'.

An in situ SAXS investigation of the formation of silver nanoparticles and bimetallic silver-gold nanoparticles in controlled wet-chemical reduction synthesis

We present a study on the formation of silver (Ag) and bimetallic silver-gold (AgAu) nanoparticles monitored by in situ SAXS as well as by ex situ TEM, XRD and UV-vis analysis in a flow reactor at controlled reaction temperature. The formation mechanism of the nanoparticles is derived from the structural parameters obtained from the experimental data. The evolution of the average particle size of pure and alloyed nanoparticles shows that the particle growth occurs initially by a coalescence mechanism. The later growth of pure silver nanoparticles is well described by Ostwald ripening and for the alloyed nanoparticles by a process with a significantly slower growth rate. Additionally, the SAXS data of pure silver nanoparticles revealed two major populations of nanoparticles, the first one with a continuous crystal growth to a saturation plateau, and the second one probably with a continuous emergence of small new crystals. The particle sizes obtained by SAXS agree well with the results from transmission electron microscopy and X-ray diffraction. The present study demonstrates the capability of an in situ investigation of synthesis processes using a laboratory based SAXS instrument. Online monitoring of the synthesis permitted a detailed investigation of the structural evolution of the system.

Oscillations in an array of bistable microelectrodes coupled through a globally conserved quantity

The dynamical behavior of an array of microelectrodes is investigated under controlled current conditions during CO electrooxidation, a bistable electrochemical reaction with an S-shaped negative differential resistance (S-NDR) current-potential curve. Under these conditions, the total current constitutes a globally conserved quantity, thus coupling all microelectrodes globally. Upon increasing the total current, the microelectrodes activate one by one, with a single microelectrode being on its intermediate S-NDR current branch and the other ones being either on their passive or their active branches. When a few coupled microelectrodes are activated, the electrochemical system exhibits spontaneous potential oscillations. Mathematical analysis shows that oscillations arise already in a two group approximation of the dynamics, the two groups consisting of 1 electrode and n - 1 electrodes with n ≥ 3, respectively, with each group being described by a single evolution equation. In this minimal representation, oscillations occur when the single electrode is on the intermediate branch and the larger group is on the active branch.

Phase diagram and aggregation dynamics of a monolayer of paramagnetic colloids

We have developed a tunable colloidal system and a corresponding theoretical model for studying the phase behavior of particles assembling under the influence of long-range magnetic interactions. A monolayer of paramagnetic particles is subjected to a spatially uniform magnetic field with a static perpendicular component and a rapidly rotating in-plane component. The sign and strength of the interactions vary with the tilt angle θ of the rotating magnetic field. For a purely in-plane field, θ=90^{∘}, interactions are attractive and the experimental results agree well with both equilibrium and out-of-equilibrium predictions based on a two-body interaction model. For tilt angles 50^{∘}≲θ≲55^{∘}, the two-body interaction gives a short-range attractive and long-range repulsive interaction, which predicts the formation of equilibrium microphases. In experiments, however, a different type of assembly is observed. Inclusion of three-body (and higher-order) terms in the model does not resolve the discrepancy. We further characterize the anomalous regime by measuring the time-dependent cluster size distribution.

Dewetting-mediated pattern formation in nanoparticle assemblies

The deposition of nanoparticles from solution onto solid substrates is a diverse subfield of current nanoscience research. Complex physical and chemical processes underpin the self-assembly and self-organization of colloidal nanoparticles at two-phase (solid-liquid, liquid-air) interfaces and three-phase (solid-liquid-air) contact lines. This review discusses key recent advances made in the understanding of nonequilibrium dewetting processes of nanoparticle-containing solutions, detailing how such an apparently simple experimental system can give rise to such a strikingly varied palette of two-dimensional self-organized nanoparticle array morphologies. Patterns discussed include worm-like domains, cellular networks, microscale rings, and fractal-like fingering structures. There remain many unresolved issues regarding the role of the solvent dewetting dynamics in assembly processes of this type, with a significant focus on how dewetting can be coerced to produce nanoparticle arrays with desirable characteristics such as long-range order. In addition to these topics, methods developed to control nanofluid dewetting through routes such as confining the geometries of drying solutions, depositing onto pre-patterned heterogeneous substrates, and post-dewetting pattern evolution via local or global manipulation are covered.

Mechanism of silver particle formation during photoreduction using in situ time-resolved SAXS analysis

Formation mechanisms of silver (Ag) particles in an aqueous ethanol solution of poly(N-vinyl-2-pyrrolidone) (PVP) by the photoreduction of AgClO(4) were investigated by means of in situ small-angle X-ray scattering (SAXS) measurements. The kinetics of association process (nucleation, growth, and coalescence) of Ag(0) atoms to produce Ag particles was successfully revealed by the quantitative SAXS analysis for the number-average of radius (R(0)), number of particles (n(Ag)), reduced standard deviation (σ(R)/R(0)), and volume fraction (ϕ(Ag)) of Ag particles produced by the photoreduction. The rate of nucleation and growth process during Ag particle formation strongly depend on the initial metal concentration. The time evolution of radius and number of Ag particles indicates that a mechanism of Ag particle formation is composed of different three processes, that is, reduction-nucleation, Ostwald ripening, and particle coalescence. In a rapid reduction-nucleation process, small nuclei or particles (average radius ∼2.5 nm) are produced by an autocatalytic reduction. After the formation of small nuclei or particles proceeds, Ostwald ripening and particle coalescence, predicted by the Lifshitz-Slyozov-Wagner theory (LSW theory), subsequently occur, resulting in the particle growth (average radius ∼11.5 nm).

Subharmonic phase clusters in the complex Ginzburg-Landau equation with nonlinear global coupling

A wide variety of subharmonic n -phase cluster patterns was observed in experiments with spatially extended chemical and electrochemical oscillators. These patterns cannot be captured with a phase model. We demonstrate that the introduction of a nonlinear global coupling (NGC) in the complex Ginzburg-Landau equation has subharmonic cluster pattern solutions in wide parameter ranges. The NGC introduces a conservation law for the oscillatory state of the homogeneous mode, which describes the strong oscillations of the mean field in the experiments. We show that the NGC causes a pronounced 2:1 self-resonance on any spatial inhomogeneity, leading to two-phase subharmonic clustering, as well as additional higher resonances. Nonequilibrium Ising-Bloch transitions occur as the coupling strength is varied.

Two-species-coagulation approach to consensus by group level interactions

We explore the self-organization dynamics of a set of entities by considering the interactions that affect the different subgroups conforming the whole. To this end, we employ the widespread example of coagulation kinetics, and characterize which interaction types lead to consensus formation and which do not, as well as the corresponding different macroscopic patterns. The crucial technical point is extending the usual one species coagulation dynamics to the two species one. This is achieved by means of introducing explicitly solvable kernels which have a clear physical meaning. The corresponding solutions are calculated in the long time limit, in which consensus may or may not be reached. The lack of consensus is characterized by means of scaling limits of the solutions. The possible applications of our results to some topics in which consensus reaching is fundamental, such as collective animal motion and opinion spreading dynamics, are also outlined.

Irregular subharmonic cluster patterns in an autonomous photoelectrochemical oscillator

Unusual subharmonic cluster patterns are observed during the oscillatory electro-oxidation of n-Si(111) under illumination. 2D in situ imaging of the electrode by means of an ellipsometric setup allows local variations in the oxide layer thickness to be monitored. The local oscillators exhibit an irregular distribution of the amplitude with the extrema locked to the constant base frequency of the total current. In addition, Ising 2-phase clustering occurs at half the base frequency. This intrinsic dynamics is described by means of a modified complex Ginzburg-Landau equation.

Nonlinear theory of nonstationary low Mach number channel flows of freely cooling nearly elastic granular gases

We employ hydrodynamic equations to investigate nonstationary channel flows of freely cooling dilute gases of hard and smooth spheres with nearly elastic particle collisions. This work focuses on the regime where the sound travel time through the channel is much shorter than the characteristic cooling time of the gas. As a result, the gas pressure rapidly becomes almost homogeneous, while the typical Mach number of the flow drops well below unity. Eliminating the acoustic modes and employing Lagrangian coordinates, we reduce the hydrodynamic equations to a single nonlinear and nonlocal equation of a reaction-diffusion type. This equation describes a broad class of channel flows and, in particular, can follow the development of the clustering instability from a weakly perturbed homogeneous cooling state to strongly nonlinear states. If the heat diffusion is neglected, the reduced equation becomes exactly soluble, and the solution develops a finite-time density blowup. The blowup has the same local features at singularity as those exhibited by the recently found family of exact solutions of the full set of ideal hydrodynamic equations [I. Fouxon, Phys. Rev. E 75, 050301(R) (2007); I. Fouxon,Phys. Fluids 19, 093303 (2007)]. The heat diffusion, however, always becomes important near the attempted singularity. It arrests the density blowup and brings about previously unknown inhomogeneous cooling states (ICSs) of the gas, where the pressure continues to decay with time, while the density profile becomes time-independent. The ICSs represent exact solutions of the full set of granular hydrodynamic equations. Both the density profile of an ICS and the characteristic relaxation time toward it are determined by a single dimensionless parameter L that describes the relative role of the inelastic energy loss and heat diffusion. At L>>1 the intermediate cooling dynamics proceeds as a competition between "holes": low-density regions of the gas. This competition resembles Ostwald ripening (only one hole survives at the end), and we report a particular regime where the "hole ripening" statistics exhibits a simple dynamic scaling behavior.

Coerced mechanical coarsening of nanoparticle assemblies

Coarsening is a ubiquitous phenomenon that underpins countless processes in nature, including epitaxial growth, the phase separation of alloys, polymers and binary fluids, the growth of bubbles in foams, and pattern formation in biomembranes. Here we show, in the first real-time experimental study of the evolution of an adsorbed colloidal nanoparticle array, that tapping-mode atomic force microscopy (TM-AFM) can drive the coarsening of Au nanoparticle assemblies on silicon surfaces. Although the growth exponent has a strong dependence on the initial sample morphology, our observations are largely consistent with modified Ostwald ripening processes. To date, ripening processes have been exclusively considered to be thermally activated, but we show that nanoparticle assemblies can be mechanically coerced towards equilibrium, representing a new approach to directed coarsening. This strategy enables precise control over the evolution of micro- and nanostructures.

Far-from-equilibrium Ostwald ripening in electrostatically driven granular powders

We report an experimental study of cluster size distributions in electrostatically driven granular submonolayers. The cluster size distribution in this far-from-equilibrium process exhibits dynamic scaling behavior characteristic of the (nearly equilibrium) Ostwald ripening, controlled by the attachment and detachment of the "gas" particles. The scaled size distribution, however, is different from the classical Wagner distribution obtained in the limit of a vanishingly small area fraction of the clusters. A much better agreement is found with the theory of Phys. Rev. E 65, 046117 (2002)] which accounts for the cluster merger.

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.

Phase-field model of island growth in epitaxy

Nucleation and growth of islands in epitaxy is simulated using a continuum phase-field model. In addition to local density of adatoms, a local phase-field variable, varying in the real space, is introduced to describe the epitaxial islands. Evolution of this phase field is determined by a time-dependent Ginzburg-Landau-like equation coupled to a diffusive transport equation of adatoms. When applied to nucleation and growth of islands in the submonolayer regime, this model reproduces both the scaling laws of island density and experimental size and spatial distributions of islands. For island growth in the multilayer regime, this phase-field model reproduces mound structures consistent with experimental images concerned. Accurate coarsening and roughening exponents of the mounds are obtained in this model. Compared with atomic models and mean-field models, this model can provide a fine visualized morphology of islands at large space and time scales of practical engineering interests.