Infrared singularities in interface growth models

Pseudospectral versus finite-difference schemes in the numerical integration of stochastic models of surface growth

We present a comparison between finite differences schemes and a pseudospectral method applied to the numerical integration of stochastic partial differential equations that model surface growth. We have studied, in 1+1 dimensions, the Kardar, Parisi, and Zhang model (KPZ) and the Lai, Das Sarma, and Villain model (LDV). The pseudospectral method appears to be the most stable for a given time step for both models. This means that the time up to which we can follow the temporal evolution of a given system is larger for the pseudospectral method. Moreover, for the KPZ model, a pseudospectral scheme gives results closer to the predictions of the continuum model than those obtained through finite difference methods. On the other hand, some numerical instabilities appearing with finite difference methods for the LDV model are absent when a pseudospectral integration is performed. These numerical instabilities give rise to an approximate multiscaling observed in earlier numerical simulations. With the pseudospectral approach no multiscaling is seen in agreement with the continuum model.

Localization in disordered media, anomalous roughening, and coarsening dynamics of faceted surfaces

We study a surface growth model related to the Kardar-Parisi-Zhang equation for nonequilibrium kinetic roughening, but where the thermal noise is replaced by a static columnar disorder eta(x) . This model is one of the many representations of the problem of particle diffusion in trapping or amplifying disordered media. We find that probability localization in the latter translates into facet formation in the equivalent surface growth problem. Coarsening of the pattern can therefore be identified with the diffusion of the localization center. The emergent faceted structure gives rise to nontrivial scaling properties, including anomalous surface roughening in excellent agreement with an existing conjecture for kinetic roughening of faceted surfaces. In a wider context, our study sheds light onto the scaling properties in other systems displaying this kind of patterned surface.

Scaling of local slopes, conservation laws, and anomalous roughening in surface growth

We argue that symmetries and conservation laws greatly restrict the form of the terms entering the long wavelength description of growth models exhibiting anomalous roughening. This is exploited to show by dynamic renormalization group arguments that intrinsic anomalous roughening cannot occur in local growth models. However, some conserved dynamics may display superroughening if a given type of term is present.

Self-consistent mode-coupling approach to the nonlocal Kardar-Parisi-Zhang equation

The dynamic scaling of the nonlocal Kardar-Parisi-Zhang equation in the strong-coupling regime is investigated by a self-consistent mode-coupling approximation. The values of the dynamic exponent depending on nonlocal parameter rho are calculated numerically for the substrate dimension d=1, d=2, and d=3, respectively.

Dynamic finite-size scaling of the normalized height distribution in kinetic surface roughening

Using well-known simple growth models, we have studied the dynamic finite-size scaling theory for the normalized height distribution of a growing surface. We find a simple functional form that explains size-dependent behavior of the skewness and kurtosis in the transient regime, and obtain the transient- and long-time values of the skewness and kurtosis for the models. Scaled distributions of the models are obtained, and the shape of each distribution is discussed in terms of the interfacial width, skewness, and kurtosis, and compared with those for other models. Exponents eta(+) and eta(-), which characterize the form of the distribution, are determined from an exponential fitting of scaling functions. Our detailed results reveal that eta(+)+eta(-) approximately 4 for a model obeying usual scaling in contrast to eta(+)+eta(-)<4 with eta(-)=1 for a model exhibiting anomalous scaling as well as multiscaling. Since we obtain eta(+)+eta(-) approximately 4 for a model exhibiting anomalous scaling but no multiscaling, we conclude that the deviation from eta(+)+eta(-) approximately 4 is due to the presence of multiscaling behavior in a model.

Interface scaling in the contact process

A correspondence between lattice models with absorbing states and models of pinned interfaces in random media can be established by defining local height variables h(x,t) as integrals of the activity at point x up to time t. Within this context we study the interface representation of a prototypical model with absorbing states, the contact process, in dimensions 1-3. Simulations confirm the scaling relation beta(W)=1-straight theta between the interface-width growth exponent beta(W) and the exponent straight theta governing the decay of the order parameter. A scaling property of the height distribution, which serves as the basis for this relation, is also verified. The height-height correlation function shows clear signs of anomalous scaling, in accord with Lopez' analysis [Phys. Rev. Lett. 83, 4594 (1999)], but no evidence of multiscaling.