Universality of fluctuations in the Kardar-Parisi-Zhang class in high dimensions and its upper critical dimension

Dimensional crossover in Kardar-Parisi-Zhang growth

Two-dimensional (2D) Kardar-Parisi-Zhang (KPZ) growth is usually investigated on substrates of lateral sizes L_{x}=L_{y}, so that L_{x} and the correlation length (ξ) are the only relevant lengths determining the scaling behavior. However, in cylindrical geometry, as well as in flat rectangular substrates L_{x}≠L_{y} and, thus, the surfaces can become correlated in a single direction, when ξ∼L_{x}≪L_{y}. From extensive simulations of several KPZ models, we demonstrate that this yields a dimensional crossover in their dynamics, with the roughness scaling as W∼t^{β_{2D}} for t≪t_{c} and W∼t^{β_{1D}} for t≫t_{c}, where t_{c}∼L_{x}^{1/z_{2D}}. The height distributions (HDs) also cross over from the 2D flat (cylindrical) HD to the asymptotic Tracy-Widom Gaussian orthogonal ensemble (Gaussian unitary ensemble) distribution. Moreover, 2D to one-dimensional (1D) crossovers are found also in the asymptotic growth velocity and in the steady-state regime of flat systems, where a family of universal HDs exists, interpolating between the 2D and 1D ones as L_{y}/L_{x} increases. Importantly, the crossover scalings are fully determined and indicate a possible way to solve 2D KPZ models.

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.

Kardar-Parisi-Zhang universality class in (d+1)-dimensions

The determination of the exact exponents of the KPZ class in any substrate dimension d is one of the most important open issues in Statistical Physics. Based on the behavior of the dimensional variation of some exact exponent differences for other growth equations, I find here that the KPZ growth exponents (related to the temporal scaling of the fluctuations) are given by β_{d}=7/8d+13. These exponents present an excellent agreement with the most accurate estimates for them in the literature. Moreover, they are confirmed here through extensive Monte Carlo simulations of discrete growth models and real-space renormalization group (RG) calculations for directed polymers in random media (DPRM), up to d=15. The left-tail exponents of the probability density functions for the DPRM energy provide another striking verification of the analytical result above.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

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.

Stirred Kardar-Parisi-Zhang Equation with Quenched Random Noise: Emergence of Induced Nonlinearity

We study the stochastic Kardar-Parisi-Zhang equation for kinetic roughening where the time-independent (columnar or spatially quenched) Gaussian random noise f(t,x) is specified by the pair correlation function ⟨f(t,x)f(t′,x′)⟩∝δ(d)(x−x′), d being the dimension of space. The field-theoretic renormalization group analysis shows that the effect of turbulent motion of the environment (modelled by the coupling with the velocity field described by the Kazantsev-Kraichnan statistical ensemble for an incompressible fluid) gives rise to a new nonlinear term, quadratic in the velocity field. It turns out that this “induced” nonlinearity strongly affects the scaling behaviour in several universality classes (types of long-time, large-scale asymptotic regimes) even when the turbulent advection appears irrelevant in itself. Practical calculation of the critical exponents (that determine the universality classes) is performed to the first order of the double expansion in ε=4−d and the velocity exponent ξ (one-loop approximation). As is the case with most “descendants” of the Kardar-Parisi-Zhang model, some relevant fixed points of the renormalization group equations lie in “forbidden zones”, i.e., in those corresponding to negative kinetic coefficients or complex couplings. This persistent phenomenon in stochastic non-equilibrium models requires careful and inventive physical interpretation.

Efimov effect at the Kardar-Parisi-Zhang roughening transition

Surface growth governed by the Kardar-Parisi-Zhang (KPZ) equation in dimensions higher than two undergoes a roughening transition from smooth to rough phases with increasing the nonlinearity. It is also known that the KPZ equation can be mapped onto quantum mechanics of attractive bosons with a contact interaction, where the roughening transition corresponds to a binding transition of two bosons with increasing the attraction. Such critical bosons in three dimensions actually exhibit the Efimov effect, where a three-boson coupling turns out to be relevant under the renormalization group so as to break the scale invariance down to a discrete one. On the basis of these facts linking the two distinct subjects in physics, we predict that the KPZ roughening transition in three dimensions shows either the discrete scale invariance or no intrinsic scale invariance.

Variational approach to KPZ: Fluctuation theorems and large deviation function for entropy production

Motivated by the time behavior of the functional arising in the variational approach to the Kardar-Parisi-Zhang (KPZ) equation, and in order to study fluctuation theorems in such a system, we have adapted a path-integral scheme that adequately fits to this kind of study dealing with unstable systems. As the KPZ system has no stationary probability distribution, we show how to proceed for obtaining detailed as well as integral fluctuation theorems. This path-integral methodology, together with the variational approach, in addition to allowing analyze fluctuation theorems, can be exploited to determine a large deviation function for entropy production.

Faceted patterns and anomalous surface roughening driven by long-range temporally correlated noise

We investigate Kardar-Parisi-Zhang (KPZ) surface growth in the presence of power-law temporally correlated noise. By means of extensive numerical simulations of models in the KPZ universality class we find that, as the noise correlator index increases above some threshold value, the surface exhibits anomalous kinetic roughening of the type described by the generic scaling theory of Ramasco et al. [Phys. Rev. Lett. 84, 2199 (2000)PRLTAO0031-900710.1103/PhysRevLett.84.2199]. Remarkably, as the driving noise temporal correlations increase, the surface develops a characteristic pattern of macroscopic facets that completely dominates the dynamics in the long time limit. We argue that standard scaling fails to capture the behavior of KPZ subject to long-range temporally correlated noise. These phenomena are not not described by the existing theoretical approaches, including renormalization group and self-consistent approaches.

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.

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.

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

We report extensive numerical simulations of growth models belonging to the nonlinear molecular beam epitaxy (nMBE) class, on flat (fixed-size) and expanding substrates (ES). In both d=1+1 and 2+1, we find that growth regime height distributions (HDs), and spatial and temporal covariances are universal, but are dependent on the initial conditions, while the critical exponents are the same for flat and ES systems. Thus, the nMBE class does split into subclasses, as does the Kardar-Parisi-Zhang (KPZ) class. Applying the "KPZ ansatz" to nMBE models, we estimate the cumulants of the 1+1 HDs. Spatial covariance for the flat subclass is hallmarked by a minimum, which is not present in the ES one. Temporal correlations are shown to decay following well-known conjectures.

Universality of (2+1)-dimensional restricted solid-on-solid models

Extensive dynamical simulations of restricted solid-on-solid models in D=2+1 dimensions have been done using parallel multisurface algorithms implemented on graphics cards. Numerical evidence is presented that these models exhibit Kardar-Parisi-Zhang surface growth scaling, irrespective of the step heights N. We show that by increasing N the corrections to scaling increase, thus smaller step-sized models describe better the asymptotic, long-wave-scaling behavior.

Scaling, cumulant ratios, and height distribution of ballistic deposition in 3+1 and 4+1 dimensions

We investigate the origin of the scaling corrections in ballistic deposition models in high dimensions using the method proposed by Alves et al. [Phys. Rev. E 90, 052405 (2014)PLEEE81539-375510.1103/PhysRevE.90.052405] in d=2+1 dimensions, where the intrinsic width associated with the fluctuations of the height increments during the deposition processes is explicitly taken into account. In the present work, we show that this concept holds for d=3+1 and 4+1 dimensions. We have found that growth and roughness exponents and dimensionless cumulant ratios are in agreement with other models, presenting small finite-time corrections to the scaling, that in principle belong to the Kardar-Parisi-Zhang (KPZ) universality class in both d=3+1 and 4+1. Our results constitute further evidence that the upper critical dimension of the KPZ class, if it exists, is larger than 4.

Width and extremal height distributions of fluctuating interfaces with window boundary conditions

We present a detailed study of squared local roughness (SLRDs) and local extremal height distributions (LEHDs), calculated in windows of lateral size l, for interfaces in several universality classes, in substrate dimensions d_{s}=1 and 2. We show that their cumulants follow a Family-Vicsek-type scaling, and, at early times, when ξ≪l (ξ is the correlation length), the rescaled SLRDs are given by log-normal distributions, with their nth cumulant scaling as (ξ/l)^{(n-1)d_{s}}. This gives rise to an interesting temporal scaling for such cumulants as 〈w_{n}〉_{c}∼t^{γ_{n}}, with γ_{n}=2nβ+(n-1)d_{s}/z=[2n+(n-1)d_{s}/α]β. This scaling is analytically proved for the Edwards-Wilkinson (EW) and random deposition interfaces and numerically confirmed for other classes. In general, it is featured by small corrections, and, thus, it yields exponents γ_{n} (and, consequently, α,β and z) in good agreement with their respective universality class. Thus, it is a useful framework for numerical and experimental investigations, where it is usually hard to estimate the dynamic z and mainly the (global) roughness α exponents. The stationary (for ξ≫l) SLRDs and LEHDs of the Kardar-Parisi-Zhang (KPZ) class are also investigated, and, for some models, strong finite-size corrections are found. However, we demonstrate that good evidence of their universality can be obtained through successive extrapolations of their cumulant ratios for long times and large l. We also show that SLRDs and LEHDs are the same for flat and curved KPZ interfaces.

From single-file diffusion to two-dimensional cage diffusion and generalization of the totally asymmetric simple exclusion process to higher dimensions

A two-dimensional constrained diffusion model is presented and characterized by numerical simulations. The model generalizes the one-dimensional single-file diffusion model by considering a cage diffusion constraint induced by neighboring particles, which is a more stringent condition than volume exclusion. Using numerical simulations we characterize the diffusion process and we particularly show that asymmetric transition probabilities lead to the two-dimensional Kardar-Parisi-Zhang universality class. Therefore, this very simple model effectively generalizes the one-dimensional totally asymmetric simple exclusion process to higher dimensions.

Strong-coupling phases of the anisotropic Kardar-Parisi-Zhang equation

We study the anisotropic Kardar-Parisi-Zhang equation using nonperturbative renormalization group methods. In contrast to a previous analysis in the weak-coupling regime, we find the strong-coupling fixed point corresponding to the isotropic rough phase to be always locally stable and unaffected by the anisotropy even at noninteger dimensions. Apart from the well-known weak-coupling and the now well-established isotropic strong-coupling behavior, we find an anisotropic strong-coupling fixed point for nonlinear couplings of opposite signs at noninteger dimensions.

Origins of scaling corrections in ballistic growth models

We study the ballistic deposition and the grain deposition models on two-dimensional substrates. Using the Kardar-Parisi-Zhang (KPZ) ansatz for height fluctuations, we show that the main contribution to the intrinsic width, which causes strong corrections to the scaling, comes from the fluctuations in the height increments along deposition events. Accounting for this correction in the scaling analysis, we obtain scaling exponents in excellent agreement with the KPZ class. We also propose a method to suppress these corrections, which consists in dividing the surface in bins of size ɛ and using only the maximal height inside each bin to do the statistics. Again, scaling exponents in remarkable agreement with the KPZ class are found. The binning method allows the accurate determination of the height distributions of the ballistic models in both growth and steady-state regimes, providing the universal underlying fluctuations foreseen for KPZ class in 2 + 1 dimensions. Our results provide complete and conclusive evidences that the ballistic model belongs to the KPZ universality class in 2+1 dimensions. Potential applications of the methods developed here, in both numerics and experiments, are discussed.