Universal and nonuniversal features in the crossover from linear to nonlinear interface growth

Growth instabilities shape morphology and genetic diversity of microbial colonies

Cellular populations assume an incredible variety of shapes ranging from circular molds to irregular tumors. While we understand many of the mechanisms responsible for these spatial patterns, little is known about how the shape of a population influences its ecology and evolution. Here, we investigate this relationship in the context of microbial colonies grown on hard agar plates. This a well-studied system that exhibits a transition from smooth circular disks to more irregular and rugged shapes as either the nutrient concentration or cellular motility is decreased. Starting from a mechanistic model of colony growth, we identify two dimensionless quantities that determine how morphology and genetic diversity of the population depend on the model parameters. Our simulations further reveal that population dynamics cannot be accurately described by the commonly-used surface growth models. Instead, one has to explicitly account for the emergent growth instabilities and demographic fluctuations. Overall, our work links together environmental conditions, colony morphology, and evolution. This link is essential for a rational design of concrete, biophysical perturbations to steer evolution in the desired direction.

Universality and crossover behavior of single-step growth models in 1+1 and 2+1 dimensions

We study the kinetic roughening of the single-step (SS) growth model with a tunable parameter p in 1+1 and 2+1 dimensions by performing extensive numerical simulations. We show that there exists a very slow crossover from an intermediate regime dominated by the Edwards-Wilkinson class to an asymptotic regime dominated by the Kardar-Parisi-Zhang (KPZ) class for any p<1/2. We also identify the crossover time, the nonlinear coupling constant, and some nonuniversal parameters in the KPZ equation as a function p. The effective nonuniversal parameters are continuously decreasing with p but not in a linear fashion. Our results provide complete and conclusive evidence that the SS model for p≠1/2 belongs to the KPZ universality class in 2+1 dimensions.

Nonuniversal effects in mixing correlated-growth processes with randomness: interplay between bulk morphology and surface roughening

To construct continuum stochastic growth equations for competitive nonequilibrium surface-growth processes of the type RD+X that mixes random deposition (RD) with a correlated-growth process X, we use a simplex decomposition of the height field. A distinction between growth processes X that do and do not create voids in the bulk leads to the definition of the effective probability p(eff) of the process X that is a measurable property of the bulk morphology and depends on the activation probability p of X in the competitive process RD+X. The bulk morphology is reflected in the surface roughening via nonuniversal prefactors in the universal scaling of the surface width that scales in p(eff). The equation and the resulting scaling are derived for X in either a Kardar-Parisi-Zhang or Edwards-Wilkinson universality class in (1+1) dimensions and are illustrated by an example of X being a ballistic deposition. We obtain full data collapse on its corresponding universal scaling function for all p∈(0;1]. We outline the generalizations to (1+n) dimensions and to many-component competitive growth processes.

Height distribution of equipotential lines in a region confined by a rough conducting boundary

This work considers the behavior of the height distributions of the equipotential lines in a region confined by two interfaces: a cathode with an irregular interface and a distant flat anode. Both boundaries, which are maintained at distinct and constant potential values, are assumed to be conductors. The morphology of the cathode interface results from the deposit of 2 × 10(4) monolayers that are produced using a single competitive growth model based on the rules of the Restricted Solid on Solid and Ballistic Deposition models, both of which belong to the Kadar-Parisi-Zhang (KPZ) universality class. At each time step, these rules are selected with probability p and q = 1 - p. For several irregular profiles that depend on p, a family of equipotential lines is evaluated. The lines are characterized by the skewness and kurtosis of the height distribution. The results indicate that the skewness of the equipotential line increases when they approach the flat anode and this increase has a non-trivial convergence to a delta distribution that characterizes the equipotential line in a uniform electric field. The morphology of the equipotential lines is discussed; the discussion emphasizes their features for different ranges of p that correspond to positive, null and negative values of the coefficient of the non-linear term in the KPZ equation.

Relaxation after a change in the interface growth dynamics

The global effects of sudden changes in the interface growth dynamics are studied using models of the Edwards-Wilkinson (EW) and Kardar-Parisi-Zhang (KPZ) classes during their growth regimes in dimensions d=1 and d=2. Scaling arguments and simulation results are combined to predict the relaxation of the difference in the roughness of the perturbed and the unperturbed interfaces, ΔW^{2}∼s{c}t{-γ}, where s is the time of the change and t>s is the observation time after that event. The previous analytical solution for the EW-EW changes is reviewed and numerically discussed in the context of lattice models, with possible decays with γ=3/2 and γ=1/2. Assuming the dominant contribution to ΔW{2} to be predicted from a time shift in the final growth dynamics, the scaling of KPZ-KPZ changes with γ=1-2β and c=2β is predicted, where β is the growth exponent. Good agreement with simulation results in d=1 and d=2 is observed. A relation with the relaxation of a local autoresponse function in d=1 cannot be discarded, but very different exponents are shown in d=2. We also consider changes between different dynamics, with the KPZ-EW as a special case in which a faster growth, with dynamical exponent z_{i}, changes to a slower one, with exponent z. A scaling approach predicts a crossover time t_{c}∼s{z/z_{i}}≫s and ΔW{2}∼s{c}F(t/t_{c}), with the decay exponent γ=1/2 of the EW class. This rules out the simplified time shift hypothesis in d=2 dimensions. These results help to understand the remarkable differences in EW smoothing of correlated and uncorrelated surfaces, and the approach may be extended to sudden changes between other growth dynamics.

Langevin equations for competitive growth models

Langevin equations for several competitive growth models in one dimension are derived. For models with crossover from random deposition (RD) to some correlated deposition (CD) dynamics, with small probability p of CD, the surface tension ν and the nonlinear coefficient λ of the associated equations have linear dependence on p due solely to this random choice. However, they also depend on the regularized step functions present in the analytical representations of the CD, whose expansion coefficients scale with p according to the divergence of local height differences when p→0. The superposition of those scaling factors gives ν~p(2) for random deposition with surface relaxation (RDSR) as the CD, and ν~p, λ~p(3/2) for ballistic deposition (BD) as the CD, in agreement with simulation and other scaling approaches. For bidisperse ballistic deposition (BBD), the same scaling of RD-BD model is found. The Langevin equation for the model with competing RDSR and BD, with probability p for the latter, is also constructed. It shows linear p dependence of λ, while the quadratic dependence observed in previous simulations is explained by an additional crossover before the asymptotic regime. The results highlight the relevance of scaling of the coefficients of step function expansions in systems with steep surfaces, which is responsible for noninteger exponents in some p-dependent stochastic equations, and the importance of the physical correspondence of aggregation rules and equation coefficients.

Conformal invariance in (2+1)-dimensional stochastic systems

Stochastic partial differential equations can be used to model second-order thermodynamical phase transitions, as well as a number of critical out-of-equilibrium phenomena. In (2+1) dimensions, many of these systems are conjectured (and some are indeed proved) to be described by conformal field theories. We advance, in the framework of the Martin-Siggia-Rose field-theoretical formalism of stochastic dynamics, a general solution of the translation Ward identities, which yields a putative conformal energy-momentum tensor. Even though the computation of energy-momentum correlators is obstructed, in principle, by dimensional reduction issues, these are bypassed by the addition of replicated fields to the original (2+1)-dimensional model. The method is illustrated with an application to the Kardar-Parisi-Zhang (KPZ) model of surface growth. The consistency of the approach is checked by means of a straightforward perturbative analysis of the KPZ ultraviolet region, leading, as expected, to its c=1 conformal fixed point.

Parameter-free scaling relation for nonequilibrium growth processes

We discuss a parameter-free scaling relation that yields a complete data collapse for large classes of nonequilibrium growth processes. We illustrate the power of this scaling relation through various growth models, such as the competitive growth model with random deposition and random deposition with surface diffusion or the restricted solid-on-solid model with different nearest-neighbor height differences, as well as through a deposition model with temperature-dependent diffusion. The scaling relation is compared to the familiar Family-Vicsek relation, and the limitations of the latter are highlighted.

Numerical study of the Kardar-Parisi-Zhang equation

We integrate numerically the Kardar-Parisi-Zhang (KPZ) equation in 1+1 and 2+1 dimensions using a Euler discretization scheme and the replacement of (nablah)(2) by exponentially decreasing functions of that quantity to suppress instabilities. When applied to the equation in 1+1 dimensions, the method of instability control provides values of scaling amplitudes consistent with exactly known results, in contrast to the deviations generated by the original scheme. In 2+1 dimensions, we spanned a range of the model parameters where transients with Edwards-Wilkinson or random growth are not observed, in box sizes 8< or =L< or =128 . We obtain a roughness exponent of 0.37< or =alpha< or =0.40 and steady state height distributions with skewness S=0.25+/-0.01 and kurtosis Q=0.15+/-0.1 . These estimates are obtained after extrapolations to the large L limit, which is necessary due to significant finite-size effects in the estimates of effective exponents and height distributions. On the other hand, the steady state roughness distributions show weak scaling corrections and evidence of stretched exponential tails. These results confirm previous estimates from lattice models, showing their reliability as representatives of the KPZ class.