Universal fluctuations of global geometrical measurements in planar clusters.

We characterize universal features of the sample-to-sample fluctuations of global geometrical observables, such as the area, width, length, and center-of-mass position, in random growing planar clusters. Our examples are taken from simulations of both continuous and discrete models of kinetically rough interfaces, including several universality classes, such as Kardar-Parisi-Zhang. We mostly focus on the scaling behavior with time of the sample-to-sample deviation for those global magnitudes, but we have also characterized their histograms and correlations.

Shape effects in the fluctuations of random isochrones on a square lattice.

We consider the isochrone curves in first-passage percolation on a 2D square lattice, i.e., the boundary of the set of points which can be reached in less than a given time from a certain origin. The occurrence of an instantaneous average shape is described in terms of its Fourier components, highlighting a crossover between a diamond and a circular geometry as the noise level is increased. Generally, these isochrones can be understood as fluctuating interfaces with an inhomogeneous local width which reveals the underlying lattice structure. We show that once these inhomogeneities have been taken into account, the fluctuations fall into the Kardar-Parisi-Zhang universality class with very good accuracy, where they reproduce the Family-Vicsek Ansatz with the expected exponents and the Tracy-Widom histogram for the local radial fluctuations.

Anomalous ballistic scaling in the tensionless or inviscid Kardar-Parisi-Zhang equation.

The one-dimensional Kardar-Parisi-Zhang (KPZ) equation is becoming an overarching paradigm for the scaling of nonequilibrium, spatially extended, classical and quantum systems with strong correlations. Recent analytical solutions have uncovered a rich structure regarding its scaling exponents and fluctuation statistics. However, the zero surface tension or zero viscosity case eludes such analytical solutions and has remained ill-understood. Using numerical simulations, we elucidate a well-defined universality class for this case that differs from that of the viscous case, featuring intrinsically anomalous kinetic roughening (despite previous expectations for systems with local interactions and time-dependent noise) and ballistic dynamics. The latter may be relevant to recent quantum spin chain experiments which measure KPZ and ballistic relaxation under different conditions. We identify the ensuing set of scaling exponents in previous discrete interface growth models related with isotropic percolation, and show it to describe the fluctuations of additional continuum systems related with the noisy Korteweg-de Vries equation. Along this process, we additionally elucidate the universality class of the related inviscid stochastic Burgers equation.

Non-KPZ fluctuations in the derivative of the Kardar-Parisi-Zhang equation or noisy Burgers equation.

The Kardar-Parisi-Zhang (KPZ) equation is a paradigmatic model of nonequilibrium low-dimensional systems with spatiotemporal scale invariance, recently highlighting universal behavior in fluctuation statistics. Its space derivative, namely the noisy Burgers equation, has played a very important role in its study, predating the formulation of the KPZ equation proper, and being frequently held as an equivalent system. We show that, while differences in the scaling exponents for the two equations are indeed due to a mere space derivative, the field statistics behave in a remarkably different way: while the KPZ equation follows the Tracy-Widom distribution, its derivative displays Gaussian behavior, hence being in a different universality class. We reach this conclusion via direct numerical simulations of the equations, supported by a dynamic renormalization group study of field statistics.

Gaussian statistics as an emergent symmetry of the stochastic scalar Burgers equation.

Symmetries play a conspicuous role in the large-scale behavior of critical systems. In equilibrium they allow us to classify asymptotics into different universality classes, and out of equilibrium, they sometimes emerge as collective properties which are not explicit in the "bare" interactions. Here we elucidate the emergence of an up-down symmetry in the asymptotic behavior of the stochastic scalar Burgers equation in one and two dimensions, manifested by the occurrence of Gaussian fluctuations even within the time regime controlled by nonlinearities. This robustness of Gaussian behavior contradicts naive expectations due to the detailed relation-including the lack of up-down symmetry-between the Burgers equation and the Kardar-Parisi-Zhang equation, which paradigmatically displays non-Gaussian fluctuations described by Tracy-Widom distributions. We reach our conclusions via a dynamic renormalization group study of the field statistics, confirmed by direct evaluation of the field probability distribution function from numerical simulations of the dynamical equation.

Nonuniversality of front fluctuations for compact colonies of nonmotile bacteria.

The front of a compact bacterial colony growing on a Petri dish is a paradigmatic instance of non-equilibrium fluctuations in the celebrated Eden, or Kardar-Parisi-Zhang (KPZ), universality class. While in many experiments the scaling exponents crucially differ from the expected KPZ values, the source of this disagreement has remained poorly understood. We have performed growth experiments with B. subtilis 168 and E. coli ATCC 25922 under conditions leading to compact colonies in the classically alleged Eden regime, where individual motility is suppressed. Non-KPZ scaling is indeed observed for all accessible times, KPZ asymptotics being ruled out for our experiments due to the monotonic increase of front branching with time. Simulations of an effective model suggest the occurrence of transient nonuniversal scaling due to diffusive morphological instabilities, agreeing with expectations from detailed models of the relevant biological reaction-diffusion processes.

Strong anisotropy in two-dimensional surfaces with generic scale invariance: nonlinear effects.

We expand a previous study [Phys. Rev. E 86, 051611 (2012)] on the conditions for occurrence of strong anisotropy in the scaling properties of two-dimensional surfaces displaying generic scale invariance. In that study, a natural scaling ansatz was proposed for strongly anisotropic systems, which arises naturally when analyzing data from, e.g., thin-film production experiments. The ansatz was tested in Gaussian (linear) models of surface dynamics and in nonlinear models, like the Hwa-Kardar (HK) equation [Phys. Rev. Lett. 62, 1813 (1989)], which are susceptible of accurate approximations through the former. In contrast, here we analyze nonlinear equations for which such approximations fail. Working within generically scale-invariant situations, and as representative case studies, we formulate and study a generalization of the HK equation for conserved dynamics and reconsider well-known systems, such as the conserved and the nonconserved anisotropic Kardar-Parisi-Zhang equations. Through the combined use of dynamic renormalization group analysis and direct numerical simulations, we conclude that the occurrence of strong anisotropy in two-dimensional surfaces requires dynamics to be conserved. We find that, moreover, strong anisotropy is not generic in parameter space but requires, rather, specific forms of the terms appearing in the equation of motion, whose justification needs detailed information on the dynamical process that is being modeled in each particular case.

Circular Kardar-Parisi-Zhang equation as an inflating, self-avoiding ring polymer.

We consider the Kardar-Parisi-Zhang equation for a circular interface in two dimensions, unconstrained by the standard small-slope and no-overhang approximations. Numerical simulations using an adaptive scheme allow us to elucidate the complete time evolution as a crossover between a short-time regime with the interface fluctuations of a self-avoiding ring or two-dimensional vesicle, and a long-time regime governed by the Tracy-Widom distribution expected for this geometry. For small-noise amplitudes, scaling behavior is only of the latter type. Large noise is also seen to renormalize the bare physical parameters of the ring, akin to analogous parameter renormalization for equilibrium three-dimensional membranes. Our results bear particular importance on the relation between relevant universality classes of scale-invariant systems in two dimensions.

Strong anisotropy in two-dimensional surfaces with generic scale invariance: Gaussian and related models.

Among systems that display generic scale invariance, those whose asymptotic properties are anisotropic in space (strong anisotropy, SA) have received relatively less attention, especially in the context of kinetic roughening for two-dimensional surfaces. This is in contrast with their experimental ubiquity, e.g., in the context of thin-film production by diverse techniques. Based on exact results for integrable (linear) cases, here we formulate a SA ansatz that, albeit equivalent to existing ones borrowed from equilibrium critical phenomena, is more naturally adapted to the type of observables that are measured in experiments on the dynamics of thin films, such as one- and two-dimensional height structure factors. We test our ansatz on a paradigmatic nonlinear stochastic equation displaying strong anisotropy like the Hwa-Kardar equation [Phys. Rev. Lett. 62, 1813 (1989)], which was initially proposed to describe the interface dynamics of running sand piles. A very important role to elucidate its SA properties is played by an accurate (Gaussian) approximation through a nonlocal linear equation that shares the same asymptotic properties.

Dynamic effects induced by renormalization in anisotropic pattern forming systems.

The dynamics of patterns in large two-dimensional domains remains a challenge in nonequilibrium phenomena. Often it is addressed through mild extensions of one-dimensional equations. We show that full two-dimensional generalizations of the latter can lead to unexpected dynamic behavior. As an example we consider the anisotropic Kuramoto-Sivashinsky equation, which is a generic model of anisotropic pattern forming systems and has been derived in different instances of thin film dynamics. A rotation of a ripple pattern by 90° occurs in the system evolution when nonlinearities are strongly suppressed along one direction. This effect originates in nonlinear parameter renormalization at different rates in the two system dimensions, showing a dynamic interplay between scale invariance and wavelength selection. Potential experimental realizations of this phenomenon are identified.

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.

Kardar-Parisi-Zhang asymptotics for the two-dimensional noisy Kuramoto-Sivashinsky equation.

We study numerically the Kuramoto-Sivashinsky equation forced by external white noise in two space dimensions, that is a generic model for, e.g., surface kinetic roughening in the presence of morphological instabilities. Large scale simulations using a pseudospectral numerical scheme allow us to retrieve Kardar-Parisi-Zhang (KPZ) scaling as the asymptotic state of the system, as in the one-dimensional (1D) case. However, this is only the case for sufficiently large values of the coupling and/or system size, so that previous conclusions on non-KPZ asymptotics are demonstrated as finite size effects. Crossover effects are comparatively stronger for the two-dimensional case than for the 1D system.

Unstable nonlocal interface dynamics.

Nonlocal effects occur in many nonequilibrium interfaces, due to diverse physical mechanisms like diffusive, ballistic, or anomalous transport, with examples from flame fronts to thin films. While dimensional analysis describes stable nonlocal interfaces, we show the morphologically unstable condition to be nontrivial. This is the case for a family of stochastic equations of experimental relevance, paradigmatically including the Michelson-Sivashinsky system. For a whole parameter range, the asymptotic dynamics is scale invariant with dimension-independent exponents reflecting a hidden Galilean symmetry. The usual Kardar-Parisi-Zhang nonlinearity, albeit irrelevant in that parameter range, plays a key role in this behavior.

Coupling of morphology to surface transport in ion-beam-irradiated surfaces: normal incidence and rotating targets.

Continuum models have proved their applicability to describe nanopatterns produced by ion-beam sputtering of amorphous or amorphizable targets at low and medium energies. Here we pursue the recently introduced 'hydrodynamic approach' in the cases of bombardment at normal incidence, or of oblique incidence onto rotating targets, known to lead to self-organized arrangements of nanodots. Our approach stresses the dynamical roles of material (defect) transport at the target surface and of local redeposition. By applying results previously derived for arbitrary angles of incidence, we derive effective evolution equations for these geometries of incidence, which are then numerically studied. Moreover, we show that within our model these equations are identical (albeit with different coefficients) in both cases, provided surface tension is isotropic in the target. We thus account for the common dynamics for both types of incidence conditions, namely formation of dots with short-range order and long-wavelength disorder, and an intermediate coarsening of dot features that improves the local order of the patterns. We provide for the first time approximate analytical predictions for the dependence of stationary dot features (amplitude and wavelength) on phenomenological parameters, that improve upon previous linear estimates. Finally, our theoretical results are discussed in terms of experimental data.

Unified moving-boundary model with fluctuations for unstable diffusive growth.

We study a moving-boundary model of nonconserved interface growth that implements the interplay between diffusive matter transport and aggregation kinetics at the interface. Conspicuous examples are found in thin-film production by chemical vapor deposition and electrochemical deposition. The model also incorporates noise terms that account for fluctuations in the diffusive and attachment processes. A small-slope approximation allows us to derive effective interface evolution equations (IEEs) in which parameters are related to those of the full moving-boundary problem. In particular, the form of the linear dispersion relation of the IEE changes drastically for slow or for instantaneous attachment kinetics. In the former case the IEE takes the form of the well-known (noisy) Kuramoto-Sivashinsky equation, showing a morphological instability at short times that evolves into kinetic roughening of the Kardar-Parisi-Zhang (KPZ) class. In the instantaneous kinetics limit, the IEE combines the Mullins-Sekerka linear dispersion relation with a KPZ nonlinearity, and we provide a numerical study of the ensuing dynamics. In all cases, the long preasymptotic transients can account for the experimental difficulties in observing KPZ scaling. We also compare our results with relevant data from experiments and discrete models.

Nonlinear ripple dynamics on amorphous surfaces patterned by ion beam sputtering.

Erosion by ion-beam sputtering (IBS) of amorphous targets at off-normal incidence frequently produces a (nanometric) rippled surface pattern, strongly resembling macroscopic ripples on aeolian sand dunes. A suitable generalization of continuum descriptions of the latter allows us to describe theoretically for the first time the main nonlinear features of ripple dynamics by IBS, namely, wavelength coarsening and nonuniform translation velocity, that agree with similar results in experiments and discrete models. These properties are seen to be the anisotropic counterparts of in-plane ordering and (interrupted) pattern coarsening in IBS experiments on rotating substrates and at normal incidence.

Single tensionless transition in the Laplacian roughening model.

We report large scale Monte Carlo simulations of the equilibrium discrete Laplacian roughening (dLr) model, originally introduced as the simplest one accommodating the hexatic phase in two-dimensional melting. The dLr model is also relevant to surface roughening in molecular beam epitaxy (MBE). Our data suggest a single phase transition, possibly of the Kosterlitz-Thouless type, between a flat low-temperature phase and a rough, tensionless, high-temperature phase. Thus, earlier conclusions on the order of the phase transition and on the existence of a hexatic phase are seen as due to finite size effects, the phase diagram of the dLr model being similar to that of a continuum analog previously formulated in the context of surface growth by MBE.

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.