Universality and dependence on initial conditions in the class of the nonlinear molecular beam epitaxy equation

Scaling of surface roughness in film deposition with height-dependent step edge barriers

We perform kinetic Monte Carlo simulations of film growth in simple cubic lattices with solid-on-solid conditions, Ehrlich-Schwoebel (ES) barriers at step edges, and a kinetic barrier related to the hidden off-plane diffusion at multilayer steps. Broad ranges of the diffusion-to-deposition ratio R, detachment probability per lateral neighbor, ε, and monolayer step crossing probability P=exp[-E_{ES}/(k_{B}T)] are studied. Without the ES barrier, four possible scaling regimes are shown as the coverage θ increases: nearly layer-by-layer growth with damped roughness oscillations; kinetic roughening in the Villain-Lai-Das Sarma (VLDS) universality class when the roughness is W∼1 (in lattice units); unstable roughening with mound nucleation and growth, where slopes of logW×logθ plots reach values larger than 0.5; and asymptotic statistical growth with W=θ^{1/2} solely due to the kinetic barrier at multilayer steps. If the ES barrier is present, the layer-by-layer growth crosses over directly to the unstable regime, with no transient VLDS scaling. However, in simulations up to θ=10^{4} (typical of films with a few micrometers), low temperatures (small R, ε, or P) may suppress the two or three initial regimes, while high temperatures and P∼1 produce smooth surfaces at all thicknesses. These crossovers help to explain proposals of nonuniversal exponents in previous works. We define a smooth film thickness θ_{c} where W=1 and show that VLDS scaling at that point indicates negligible ES barriers, while rapidly increasing roughness indicates a small ES barrier (E_{ES}∼k_{B}T). θ_{c} scales as ∼exp(const×P^{2/3}) if the other parameters are kept fixed, which represents a high sensitivity on the ES barrier. The analysis of recent experimental data in the light of our results distinguishes cases where E_{ES}/(k_{B}T) is negligible, ∼1, or ≪1.

Universal interface fluctuations in the contact process

We study the interface representation of the contact process at its directed-percolation critical point, where the scaling properties of the interface can be related to those of the original particle model. Interestingly, such a behavior happens to be intrinsically anomalous and more complex than that described by the standard Family-Vicsek dynamic scaling Ansatz of surface kinetic roughening. We expand on a previous numerical study by Dickman and Muñoz [Phys. Rev. E 62, 7632 (2000)10.1103/PhysRevE.62.7632] to fully characterize the kinetic roughening universality class for interface dimensions d=1,2, and 3. Beyond obtaining scaling exponent values, we characterize the interface fluctuations via their probability density function (PDF) and covariance, seen to display universal properties which are qualitatively similar to those recently assessed for the Kardar-Parisi-Zhang (KPZ) and other important universality classes of kinetic roughening. Quantitatively, while for d=1 the interface covariance seems to be well described by the KPZ, Airy_{1} covariance, no such agreement occurs in terms of the fluctuation PDF or the scaling exponents.

One-point height fluctuations and two-point correlators of (2+1) cylindrical KPZ systems

While the one-point height distributions (HDs) and two-point covariances of (2+1) Kardar-Parisi-Zhang (KPZ) systems have been investigated in several recent works for flat and spherical geometries, for the cylindrical one the HD was analyzed for few models and nothing is known about the spatial and temporal covariances. Here, we report results for these quantities, obtained from extensive numerical simulations of discrete KPZ models, for three different setups yielding cylindrical growth. Beyond demonstrating the universality of the HD and covariances, our results reveal other interesting features of this geometry. For example, the spatial covariances measured along the longitudinal and azimuthal directions are different, with the former being quite similar to the curve for flat (2+1) KPZ systems, while the latter resembles the Airy_{2} covariance of circular (1+1) KPZ interfaces. We also argue (and present numerical evidence) that, in general, the rescaled temporal covariance A(t/t_{0}) decays asymptotically as A(x)∼x^{-λ[over ¯]} with an exponent λ[over ¯]=β+d^{*}/z, where d^{*} is the number of interface sides kept fixed during the growth (being d^{*}=1 for the systems analyzed here). Overall, these results complete the picture of the main statistics for the (2+1) KPZ class.

Height distributions in interface growth: The role of the averaging process

Height distributions (HDs) are key quantities to uncover universality and geometry-dependence in evolving interfaces. To quantitatively characterize HDs, one uses adimensional ratios of their first central moments (m_{n}) or cumulants (κ_{n}), especially the skewness S and kurtosis K, whose accurate estimate demands an averaging over all L^{d} points of the height profile at a given time, in translation-invariant interfaces, and over N independent samples. One way of doing this is by calculating m_{n}(t) [or κ_{n}(t)] for each sample and then carrying out an average of them for the N interfaces, with S and K being calculated only at the end. Another approach consists in directly calculating the ratios for each interface and, then, averaging the N values. It turns out, however, that S and K for the growth regime HDs display strong finite-size and -time effects when estimated from these "interface statistics," as already observed in some previous works and clearly shown here, through extensive simulations of several discrete growth models belonging to the EW and KPZ classes on one- and two-dimensional substrates of sizes L=const. and L∼t. Importantly, I demonstrate that with "1-point statistics," i.e., by calculating m_{n}(t) [or κ_{n}(t)] once for all NL^{d} heights together, these corrections become very weak, so that S and K attain values very close to the asymptotic ones already at short times and for small L's. However, I find that this "1-point" (1-pt) approach fails in uncovering the universality of the HDs in the steady-state regime (SSR) of systems whose average height, h[over ¯], is a fluctuating variable. In fact, as demonstrated here, in this regime the 1-pt height evolves as h(t)=h[over ¯](t)+s_{λ}A^{1/2}L^{α}ζ+⋯-where P(ζ) is the underlying SSR HD-and the fluctuations in h[over ¯] yield S_{1-pt}∼t^{-1/2} and K_{1-pt}∼t^{-1}. Nonetheless, by analyzing P(h-h[over ¯]), the cumulants of P(ζ) can be accurately determined. I also show that different, but universal, asymptotic values for S and K (related, so, to different HDs) can be found from the "interface statistics" in the SSR. This reveals the importance of employing the various complementary approaches to reliably determine the universality class of a given system through its different HDs.

Spreading fronts of wetting liquid droplets: Microscopic simulations and universal fluctuations

We have used kinetic Monte Carlo (kMC) simulations of a lattice gas to study front fluctuations in the spreading of a nonvolatile liquid droplet onto a solid substrate. Our results are consistent with a diffusive growth law for the radius of the precursor layer, R∼t^{δ}, with δ≈1/2 in all the conditions considered for temperature and substrate wettability, in good agreement with previous studies. The fluctuations of the front exhibit kinetic roughening properties with exponent values which depend on temperature T, but become T independent for sufficiently high T. Moreover, strong evidence of intrinsic anomalous scaling has been found, characterized by different values of the roughness exponent at short and large length scales. Although such a behavior differs from the scaling properties of the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class, the front covariance and the probability distribution function of front fluctuations found in our kMC simulations do display KPZ behavior, agreeing with simulations of a continuum height equation proposed in this context. However, this equation does not feature intrinsic anomalous scaling, at variance with the discrete model.

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.

Geometry dependence in linear interface growth

The effect of geometry in the statistics of nonlinear universality classes for interface growth has been widely investigated in recent years, and it is well known to yield a split of them into subclasses. In this work, we investigate this for the linear classes of Edwards-Wilkinson and of Mullins-Herring in one and two dimensions. From comparison of analytical results with extensive numerical simulations of several discrete models belonging to these classes, as well as numerical integrations of the growth equations on substrates of fixed size (flat geometry) or expanding linearly in time (radial geometry), we verify that the height distributions (HDs) and the spatial and the temporal covariances are universal but geometry-dependent. In fact, the HDs are always Gaussian, and, when defined in terms of the so-called "KPZ ansatz" [h≃v_{∞}t+(Γt)^{β}χ], their probability density functions P(χ) have mean null, so that all their cumulants are null, except by their variances, which assume different values in the flat and radial cases. The shape of the (rescaled) covariance curves is analyzed in detail and compared with some existing analytical results for them. Overall, these results demonstrate that the splitting of such university classes is quite general, being not restricted to the nonlinear ones.

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.

Local roughness exponent in the nonlinear molecular-beam-epitaxy universality class in one dimension

We report local roughness exponents, α_{loc}, for three interface growth models in one dimension which are believed to belong to the nonlinear molecular-beam-epitaxy (nMBE) universality class represented by the Villain-Lais-Das Sarma (VLDS) stochastic equation. We applied an optimum detrended fluctuation analysis (ODFA) [Luis et al., Phys. Rev. E 95, 042801 (2017)2470-004510.1103/PhysRevE.95.042801] and compared the outcomes with standard detrending methods. We observe in all investigated models that ODFA outperforms the standard methods providing exponents in the narrow interval α_{loc}^{}∈[0.96,0.98] quantitatively consistent with two-loop renormalization group predictions for the VLDS equation. In particular, these exponent values are calculated for the Clarke-Vvdensky and Das Sarma-Tamborenea models characterized by very strong corrections to the scaling, for which large deviations of these values had been reported. Our results strongly support the absence of anomalous scaling in the nMBE universality class and the existence of corrections in the form α_{loc}^{}=1-ε of the one-loop renormalization group analysis of the VLDS equation.

Circular Kardar-Parisi-Zhang interfaces evolving out of the plane

Circular KPZ interfaces spreading radially in the plane have Gaussian unitary ensemble (GUE) Tracy-Widom (TW) height distribution (HD) and Airy_{2} spatial covariance, but what are their statistics if they evolve on the surface of a different background space, such as a bowl, a mountain, or any surface of revolution? To give an answer to this, we report here extensive numerical analyses of several one-dimensional KPZ models on substrates whose size enlarges as 〈L(t)〉=L_{0}+ωt^{γ}, while their mean height 〈h〉 increases as usual [〈h〉∼t]. We show that the competition between the L enlargement and the correlation length (ξ≃ct^{1/z}) plays a key role in the asymptotic statistics of the interfaces. While systems with γ>1/z have HDs given by GUE and the interface width increasing as w∼t^{β}, for γ<1/z the HDs are Gaussian, in a correlated regime where w∼t^{αγ}. For the special case γ=1/z, a continuous class of distributions exists, which interpolate between Gaussian (for small ω/c) and GUE (for ω/c≫1). Interestingly, the HD seems to agree with the Gaussian symplectic ensemble (GSE) TW distribution for ω/c≈10. Despite the GUE HDs for γ>1/z, the spatial covariances present a strong dependence on the parameters ω and γ, agreeing with Airy_{2} only for ω≫1, for a given γ, or when γ=1, for a fixed ω. These results considerably generalize our knowledge on 1D KPZ systems, unveiling the importance of the background space on their statistics.

Kardar-Parisi-Zhang growth on one-dimensional decreasing substrates

Recent experimental works on one-dimensional (1D) circular Kardar-Parisi-Zhang (KPZ) systems whose radii decrease in time have reported controversial conclusions about the statistics of their interfaces. Motivated by this, here we investigate several one-dimensional KPZ models on substrates whose size changes in time as L(t)=L_{0}+ωt, focusing on the case ω<0. From extensive numerical simulations, we show that for L_{0}≫1 there exists a transient regime in which the statistics is consistent with that of flat KPZ systems (the ω=0 case), for both ω<0 and ω>0. Actually, for a given model, L_{0} and |ω|, we observe that a difference between ingrowing (ω<0) and outgrowing (ω>0) systems arises only at long times (t∼t_{c}=L_{0}/|ω|), when the expanding surfaces cross over to the statistics of curved KPZ systems, whereas the shrinking ones become completely correlated. A generalization of the Family-Vicsek scaling for the roughness of ingrowing interfaces is presented. Our results demonstrate that a transient flat statistics is a general feature of systems starting with large initial sizes, regardless of their curvature. This is consistent with their recent observation in ingrowing turbulent liquid crystal interfaces, but it is in contrast with the apparent observation of curved statistics in colloidal deposition at the edge of evaporating drops. A possible explanation for this last result, as a consequence of the very small number of monolayers analyzed in this experiment, is given. This is illustrated in a competitive growth model presenting a few-monolayer transient and an asymptotic behavior consistent, respectively, with the curved and flat statistics.

Radial restricted solid-on-solid and etching interface-growth models

An approach to generate radial interfaces is presented. A radial network recursively obtained is used to implement discrete model rules designed originally for the investigation in flat substrates. I used the restricted solid-on-solid and etching models as to test the proposed scheme. The results indicate the Kardar, Parisi, and Zhang conjecture is completely verified leading to a good agreement between the interface radius fluctuation distribution and the Gaussian unitary ensemble. The evolution of the radius agrees well with the generalized conjecture, and the two-point correlation function exhibits also a good agreement with the covariance of the Airy_{2} process. The approach can be used to investigate radial interfaces evolution for many other classes of universality.

Kardar-Parisi-Zhang Interfaces with Inward Growth

We study the (1+1)-dimensional Kardar-Parisi-Zhang (KPZ) interfaces growing inward from ring-shaped initial conditions, experimentally and numerically, using growth of a turbulent state in liquid-crystal electroconvection and an off-lattice Eden model, respectively. To realize the ring initial condition experimentally, we introduce a holography-based technique that allows us to design the initial condition arbitrarily. Then, we find that fluctuation properties of ingrowing circular interfaces are distinct from those for the curved or circular KPZ subclass and, instead, are characterized by the flat subclass. More precisely, we find an asymptotic approach to the Tracy-Widom distribution for the Gaussian orthogonal ensemble and the Airy_{1} spatial correlation, as long as time is much shorter than the characteristic time determined by the initial curvature. Near this characteristic time, deviation from the flat KPZ subclass is found, which can be explained in terms of the correlation length and the circumference. Our results indicate that the sign of the initial curvature has a crucial role in determining the universal distribution and correlation functions of the KPZ class.

Initial pseudo-steady state & asymptotic KPZ universality in semiconductor on polymer deposition

The Kardar-Parisi-Zhang (KPZ) class is a paradigmatic example of universality in nonequilibrium phenomena, but clear experimental evidences of asymptotic 2D-KPZ statistics are still very rare, and far less understanding stems from its short-time behavior. We tackle such issues by analyzing surface fluctuations of CdTe films deposited on polymeric substrates, based on a huge spatio-temporal surface sampling acquired through atomic force microscopy. A pseudo-steady state (where average surface roughness and spatial correlations stay constant in time) is observed at initial times, persisting up to deposition of ~10⁴ monolayers. This state results from a fine balance between roughening and smoothening, as supported by a phenomenological growth model. KPZ statistics arises at long times, thoroughly verified by universal exponents, spatial covariance and several distributions. Recent theoretical generalizations of the Family-Vicsek scaling and the emergence of log-normal distributions during interface growth are experimentally confirmed. These results confirm that high vacuum vapor deposition of CdTe constitutes a genuine 2D-KPZ system, and expand our knowledge about possible substrate-induced short-time behaviors.

Memory and Universality in Interface Growth

Recently, very robust universal properties have been shown to arise in one-dimensional growth processes with local stochastic rules, leading to the Kardar-Parisi-Zhang (KPZ) universality class. Yet it has remained essentially unknown how fluctuations in these systems correlate at different times. Here, we derive quantitative predictions for the universal form of the two-time aging dynamics of growing interfaces and we show from first principles the breaking of ergodicity that the KPZ time evolution exhibits. We provide corroborating experimental observations on a turbulent liquid crystal system, as well as a numerical simulation of the Eden model, and we demonstrate the universality of our predictions. These results may give insight into memory effects in a broader class of far-from-equilibrium systems.