Dynamic scaling and crossover analysis for the Kuramoto-Sivashinsky equation

Out-of-time-ordered correlator in the one-dimensional Kuramoto-Sivashinsky and Kardar-Parisi-Zhang equations

The out-of-time-ordered correlator (OTOC) has emerged as an interesting object in both classical and quantum systems for probing the spatial spread and temporal growth of initially local perturbations in spatially extended chaotic systems. Here, we study the (classical) OTOC and its "light cone" in the nonlinear Kuramoto-Sivashinsky (KS) equation, using extensive numerical simulations. We also show that the linearized KS equation exhibits a qualitatively similar OTOC and light cone, which can be understood via a saddle-point analysis of the linearly unstable modes. Given the deep connection between the KS (deterministic) and the Kardar-Parisi-Zhang (KPZ, which is stochastic) equations, we also explore the OTOC in the KPZ equation. While our numerical results in the KS case are expected to hold in the continuum limit, for the KPZ case it is valid in a discretized version of the KPZ equation. More broadly, our work unravels the intrinsic interplay between noise/instability, nonlinearity, and dissipation in partial differential equations (deterministic or stochastic) through the lens of OTOC.

Divergent Stiffness of One-Dimensional Growing Interfaces

When a spatially localized stress is applied to a growing one-dimensional interface, the interface deforms. This deformation is described by the effective surface tension representing the stiffness of the interface. We present that the stiffness exhibits divergent behavior in the large system size limit for a growing interface with thermal noise, which has never been observed for equilibrium interfaces. Furthermore, by connecting the effective surface tension with a space-time correlation function, we elucidate the mechanism that anomalous dynamical fluctuations lead to divergent stiffness.

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.

One-dimensional Kardar-Parisi-Zhang and Kuramoto-Sivashinsky universality class: Limit distributions

Tracy-Widom and Baik-Rains distributions appear as universal limit distributions for height fluctuations in the one-dimensional Kardar-Parisi-Zhang (KPZ) stochastic partial differential equation (PDE). We obtain the same universal distributions in the spatiotemporally chaotic, nonequilibrium, but statistically steady state of the one-dimensional Kuramoto-Sivashinsky (KS) deterministic PDE, by carrying out extensive pseudospectral direct numerical simulations to obtain the spatiotemporal evolution of the KS height profile h(x,t) for different initial conditions. We establish, therefore, that the statistical properties of the one-dimensional (1D) KS PDE in this state are in the 1D KPZ universality class.

Coherent X-ray measurement of step-flow propagation during growth on polycrystalline thin film surfaces

The properties of artificially grown thin films are strongly affected by surface processes during growth. Coherent X-rays provide an approach to better understand such processes and fluctuations far from equilibrium. Here we report results for vacuum deposition of C₆₀ on a graphene-coated surface investigated with X-ray Photon Correlation Spectroscopy in surface-sensitive conditions. Step-flow is observed through measurement of the step-edge velocity in the late stages of growth after crystalline mounds have formed. We show that the step-edge velocity is coupled to the terrace length, and that there is a variation in the velocity from larger step spacing at the center of crystalline mounds to closely-spaced, more slowly propagating steps at their edges. The results extend theories of surface growth, since the behavior is consistent with surface evolution driven by processes that include surface diffusion, the motion of step-edges, and attachment at step edges with significant step-edge barriers.

Nonlinear dynamics of a dispersive anisotropic Kuramoto-Sivashinsky equation in two space dimensions

A Kuramoto-Sivashinsky equation in two space dimensions arising in thin film flows is considered on doubly periodic domains. In the absence of dispersive effects, this anisotropic equation admits chaotic solutions for sufficiently large length scales with fully two-dimensional profiles; the one-dimensional dynamics observed for thin domains are structurally unstable as the transverse length increases. We find that, independent of the domain size, the characteristic length scale of the profiles in the streamwise direction is about 10 space units, with that in the transverse direction being approximately three times larger. Numerical computations in the chaotic regime provide an estimate for the radius of the absorbing ball in [Formula: see text] in terms of the length scales, from which we conclude that the system possesses a finite energy density. We show the property of equipartition of energy among the low Fourier modes, and report the disappearance of the inertial range when solution profiles are two-dimensional. Consideration of the high-frequency modes allows us to compute an estimate for the analytic extensibility of solutions in [Formula: see text]. We also examine the addition of a physically derived third-order dispersion to the problem; this has a destabilizing effect, in the sense of reducing analyticity and increasing amplitude of solutions. However, sufficiently large dispersion may regularize the spatio-temporal chaos to travelling waves. We focus on dispersion where chaotic dynamics persist, and study its effect on the interfacial structures, absorbing ball and properties of the power spectrum.

Simulations of Co-GISAXS during kinetic roughening of growth surfaces

The recent development of surface growth studies using X-ray photon correlation spectroscopy in a grazing-incidence small-angle X-ray scattering (Co-GISAXS) geometry enables the investigation of dynamical processes during kinetic roughening in greater detail than was previously possible. In order to investigate the Co-GISAXS behavior expected from existing growth models, calculations and (2+1)-dimension simulations of linear Kuramoto-Sivashinsky and non-linear Kardar-Parisi-Zhang surface growth equations are presented which analyze the temporal correlation functions of the height-height structure factor. Calculations of the GISAXS intensity auto-correlation functions are also performed within the Born/distorted-wave Born approximation for comparison with the scaling behavior of the height-height structure factor and its correlation functions.

Observation and modeling of interrupted pattern coarsening: surface nanostructuring by ion erosion

We report the experimental observation of interrupted coarsening for surface self-organized nanostructuring by ion erosion. Analysis of the target surface by atomic force microscopy allows us to describe quantitatively this intriguing type of pattern dynamics through a continuum equation put forward in different contexts across a wide range of length scales. The ensuing predictions can thus be consistently extended to other experimental conditions in our system. Our results illustrate the occurrence of nonequilibrium systems in which pattern formation, coarsening, and kinetic roughening appear, each of these behaviors being associated with its own spatiotemporal range.

Dual structures of chaos and turbulence, and their dynamic scaling laws

The decay form of the time correlation function Un(t) of a state variable un(t) with a small wave number kn has been shown to take the algebraic decay 1/{1+(gammanat)2} in the initial regime t<taun(gamma) and the exponential decay alphane exp(-gammanet) in the final regime t>taun(gamma), where taun(gamma) denotes the decay time of the memory function Gamman(t). This dual structure of Un(t) is generated by the deterministic short orbits in the initial regime and the stochastic long orbits in the final regime, thus giving the outstanding features of chaos and turbulence. The kn dependence of gammana, alphane, and gammane is obtained for the chaotic Kuramoto-Sivashinsky equation, and it is shown that if kn is sufficiently small, then the dual structure of Un(t) obeys a hydrodynamic scaling law in the final regime t>tau(gamma) with scaling exponent z=2 and a dynamic scaling law in the initial regime t<taun(gamma) with scaling exponent z=1. If kn is increased so that the decay time taun(u) of Un(t) becomes equal to the decay time taun(gamma), then the decay form of Un(t) becomes the power-law decay t-3/2 in the final regime.

Single-phase-field model of stepped surfaces

We formulate a phase-field description of step dynamics on vicinal surfaces that makes use of a single dynamical field, at variance with previous analogous works in which two coupled fields are employed, namely, a phase-field proper plus the physical adatom concentration. Within an asymptotic sharp interface limit, our formulation is shown to retrieve the standard Burton-Cabrera-Frank model in the general case of asymmetric attachment coefficients (Ehrlich-Schwoebel effect). We confirm our analytical results by means of numerical simulations of our phase-field model. Our present formulation seems particularly well adapted to generalization when additional physical fields are required.

Soliton turbulence in the complex Ginzburg-Landau equation

We study spatiotemporal chaos in the complex Ginzburg-Landau equation in parameter regions of weak amplification and viscosity. Turbulent states involving many solitonlike pulses appear in the parameter range, because the complex Ginzburg-Landau equation is close to the nonlinear Schrödinger equation. We find that the distributions of amplitude and wave number of pulses depend only on the ratio of the two parameters of the amplification and the viscosity. This implies that a one-parameter family of soliton turbulence states characterized by different distributions of the soliton parameters exists continuously around the completely integrable system.

Sputter roughening of inhomogeneous surfaces: impurity pinning and nanostructure shape selection

A model is proposed for sputter roughening of inhomogeneous systems with slowly sputtered impurity particles randomly distributed in the bulk. Surface inhomogeneity, which develops as a result of coupling between the time evolution of the local surface impurity concentration and the local surface shape, is tuned by changing the dependence of the sputtering probability upon impurity concentration. In 1+1 dimensions, we find long-time scaling exponents that are consistent with Kardar-Parisi-Zhang (KPZ) values. However, for a range of surface inhomogeneity, impurity pinning results in a persistent growth regime where the surface roughens rapidly. We correlate this rapid roughening to fluctuations of the impurity concentration at the surface. Roughening in this regime leads to the formation of cones whose shape is determined by material property and sputtering flux, suggesting a unique method of nanostructure fabrication. In 2+1 dimensions, a similar variation of the roughening behavior with surface inhomogeneity is observed. For small surface inhomogeneity, there is an initial exponential roughening followed by power-law roughening with an effective growth exponent much smaller than KPZ. For larger surface inhomogeneity two power-law roughening regimes are observed, with an initial rapid roughening that crosses over to slower roughening; the effective exponent in each of these regimes increases with surface inhomogeneity. The surface morphology observed in the simulations is considerably noisier than experimental data for InP and GaSb. Our model shows noisy nonlinear pattern formation in contrast to the marked long-range hexagonal ordering seen in experiments. However, the scaling behavior is robust enough that roughening kinetics similar to that observed experimentally can be obtained depending upon the values of inhomogeneity and the strength of the nonlinear term in the model.

Short-range stationary patterns and long-range disorder in an evolution equation for one-dimensional interfaces

A local evolution equation for one-dimensional interfaces is derived in the context of erosion by ion beam sputtering. We present numerical simulations of this equation which show interrupted coarsening in which an ordered cell pattern develops with constant wavelength and amplitude at intermediate distances, while the profile is disordered and rough at larger distances. Moreover, for a wide range of parameters the lateral extent of ordered domains ranges up to tens of cells. We also provide analytical estimates for the stationary pattern wavelength and mean growth velocity.

Validity of the essential assumption in a projection operator method

The projection operator method developed by Mori involves the essential assumption that chaotic motion is successfully divided into a coherent motion and a fluctuating one. We investigate the validity of the assumption using the Kuramoto-Sivashinsky equation as a model equation of chaotic systems. It has been found that the assumption is reasonable for both long wave modes and short wave modes. We have also evaluated a value of the eddy viscosity as 9.0 by extracting the nonlinear term from the coherent part. This value is consistent with the former estimates with other methods.

Renormalization-group and numerical analysis of a noisy Kuramoto-Sivashinsky equation in 1 + 1 dimensions

The long-wavelength properties of a noisy Kuramoto-Sivashinsky (KS) equation in 1 + 1 dimensions are investigated by use of the dynamic renormalization group (RG) and direct numerical simulations. It is shown that the noisy KS equation is in the same universality class as the Kardar-Parisi-Zhang (KPZ) equation in the sense that they have scale invariant solutions with the same scaling exponents in the long-wavelength limit. The RG analysis reveals that the RG flow for the parameters of the noisy KS equation rapidly approach the KPZ fixed point with increasing strength of the noise. This is supplemented by numerical simulations of the KS equation with a stochastic noise, in which scaling behavior close to the KPZ scaling can be observed even in a moderate system size and time.

Directed percolation with long-range interactions: Modeling nonequilibrium wetting

It is argued that some phase transitions observed in models of nonequilibrium wetting phenomena are related to contact processes with long-range interactions. This is investigated by introducing a model where the activation rate of a site at the edge of an inactive island of length l is 1+a l(-sigma) . Mean-field analysis and numerical simulations indicate that for sigma>1 the transition is continuous and belongs to the universality class of directed percolation, while for 0<sigma<1 , the transition becomes first order. This criterion is then applied to discuss critical properties of various models of nonequilibrium wetting.

Mean solutions for the Kuramoto-Sivashinsky equation with incoming boundary conditions

We consider herein the Kuramoto-Sivashinsky (KS) equation with incoming boundary conditions. Using a projection operator method, we have derived a set of closed equations for the mean quantities, called a model equation, from the KS equation. One of the characteristics of the model equation is that it does not include any empirical parameters. The adequacy of the model equation is verified by comparing solutions of the model equation with time-averaged solutions obtained from the numerical simulation of the KS equation.

Dynamical scaling of sputter-roughened surfaces in 2+1 dimensions

The asymptotic scaling behavior of sputter-roughened surfaces is of great current interest. In particular, the disparately wide-ranging values of the growth exponent found experimentally, and whether sputter-roughening belongs to the Kardar-Parisi-Zhang universality class in 2+1 dimensions, are two interesting issues. We address these issues using simulations of an atomistic model. The asymptotic scaling appears to be Edwards-Wilkinson. Crossover behavior in the model leads to effective growth exponents that vary widely depending upon the regime of observation.

Spontaneous formation of dynamical patterns with fractal fronts in the cyclic lattice Lotka-Volterra model

Dynamical patterns, in the form of consecutive moving stripes or rings, are shown to develop spontaneously in the cyclic lattice Lotka-Volterra model, when realized on square lattice, at the reaction limited regime. Each stripe consists of different particles (species) and the borderlines between consecutive stripes are fractal. The interface width w between the different species scales as w(L,t) approximately L(alpha)f(t/L(z)), where L is the linear size of the interface, t is the time, and alpha and z are the static and dynamical critical exponents, respectively. The critical exponents were computed as alpha=0.49+/-0.03 and z=1.53+/-0.13 and the propagating fronts show dynamical characteristics similar to those of the Eden growth models.

Temporal and spatial persistence of combustion fronts in paper

The spatial and temporal persistence, or first-return distributions are measured for slow-combustion fronts in paper. The stationary temporal and (perhaps less convincingly) spatial persistence exponents agree with the predictions based on the front dynamics, which asymptotically belongs to the Kardar-Parisi-Zhang universality class. The stationary short-range and the transient behavior of the fronts are non-Markovian, and the observed persistence properties thus do not agree with the predictions based on Markovian theory. This deviation is a consequence of additional time and length scales, related to the crossovers to the asymptotic coarse-grained behavior.

Pseudospectral method for the Kardar-Parisi-Zhang equation

We discuss a numerical scheme to solve the continuum Kardar-Parisi-Zhang equation in generic spatial dimensions. It is based on a momentum-space discretization of the continuum equation and on a pseudospectral approximation of the nonlinear term. The method is tested in (1+1) and (2+1) dimensions, where it is shown to reproduce the current most reliable estimates of the critical exponents based on restricted solid-on-solid simulations. In particular, it allows the computations of various correlation and structure functions with high degree of numerical accuracy. Some deficiencies that are common to all previously used finite-difference schemes are pointed out and the usefulness of the present approach in this respect is discussed.

Transients due to instabilities hinder Kardar-Parisi-Zhang scaling: a unified derivation for surface growth by electrochemical and chemical vapor deposition

We propose a unified moving boundary problem for surface growth by electrochemical and chemical vapor deposition, which is derived from constitutive equations into which stochastic forces are incorporated. We compute the coefficients in the interface equation of motion as functions of phenomenological parameters. The equation features the Kardar-Parisi-Zhang (KPZ) nonlinearity and instabilities which, depending on surface kinetics, can hinder the asymptotic KPZ scaling. Our results account for the universality and the experimental scarcity of KPZ scaling in the growth processes considered.

Heuristic model for the energy spectrum of phase turbulence

We present a heuristic model for the energy spectrum of the one-dimensional phase turbulence in the steady state of the Kuramoto-Sivashinsky equation. Our model contains an energy transfer mechanism from low- to high-wave-vector modes. The energy transfer is written as the sum of local and nonlocal interactions. Our analytical results show good agreement with numerical simulations, particularly for the hump in the energy spectrum, which is mainly due to the local interactions.

Kinetic roughening in slow combustion of paper

Results of experiments on the dynamics and kinetic roughening of one-dimensional slow-combustion fronts in three grades of paper are reported. Extensive averaging of the data allows a detailed analysis of the spatial and temporal development of the interface fluctuations. The asymptotic scaling properties, on long length and time scales, are well described by the Kardar-Parisi-Zhang (KPZ) equation with short-range, uncorrelated noise. To obtain a more detailed picture of the strong-coupling fixed point, characteristic of the KPZ universality class, universal amplitude ratios, and the universal coupling constant are computed from the data and found to be in good agreement with theory. Below the spatial and temporal scales at which a crossover takes place to the standard KPZ behavior, the fronts display higher apparent exponents and apparent multiscaling. In this regime the interface velocities are spatially and temporally correlated, and the distribution of the magnitudes of the effective noise has a power-law tail. The relation of the observed short-range behavior and the noise as determined from the local velocity fluctuations is discussed.

Shock structures in time-averaged patterns for the kuramoto-sivashinsky equation

The Kuramoto-Sivashinsky equation with fixed boundary conditions is numerically studied. Shocklike structures appear in the time-averaged patterns for some parameter range of the boundary values. Effective diffusion constant is estimated from the relation of the width and the height of the shock structures.

Scaling and noise in slow combustion of paper

We present results of high resolution experiments on kinetic roughening of slow combustion fronts in paper, focusing on short length and time scales. Using three different grades of paper, we find that the combustion fronts show apparent spatial and temporal multiscaling at short scales. The scaling exponents decrease as a function of the order of the corresponding correlation functions. The noise affecting the fronts reveals short range temporal and spatial correlations, and non-Gaussian noise amplitudes. Our results imply that the overall behavior of slow combustion fronts cannot be explained by standard theories of kinetic roughening.

Scale and space localization in the Kuramoto-Sivashinsky equation

We describe a wavelet-based approach to the investigation of spatiotemporally complex dynamics, and show through extensive numerical studies that the dynamics of the Kuramoto-Sivashinsky equation in the spatiotemporally chaotic regime may be understood in terms of localized dynamics in both space and scale (wave number). A projection onto a spline wavelet basis enables good separation of scales, each with characteristic dynamics. At the large scales, one observes essentially slow Gaussian dynamics; at the active scales, structured "events" reminiscent of traveling waves and heteroclinic cycles appear to dominate; while the strongly damped small scales display intermittent behavior. The separation of scales and their dynamics is invariant as the length of the system increases, providing additional support for the extensivity of the spatiotemporally complex dynamics claimed in earlier works. We show also that the dynamics are spatially localized, discuss various correlation lengths, and demonstrate the existence of a characteristic interaction length for instantaneous influences. Our results motivate and advance the search for localized, low-dimensional models that capture the full behavior of spatially extended chaotic partial differential equations. (c) 1999 American Institute of Physics.

Collective patterns arising out of spatio-temporal chaos

We present a simple mathematical model in which a time averaged pattern emerges out of spatio-temporal chaos as a result of the collective action of chaotic fluctuations. Our evolution equation possesses spatial translational symmetry under periodic boundary conditions. Thus the spatial inhomogeneity of the statistical state arises through spontaneous symmetry breaking. The transition from a state of homogeneous spatio-temporal chaos to one exhibiting spatial order is explained by introducing a collective viscosity which relates the averaged pattern with a correlation of the fluctuations. (c) 1996 American Institute of Physics.