Multisurface coding simulations of the restricted solid-on-solid model in four dimensions

Kardar-Parisi-Zhang growth in ɛ dimensions and beyond

We examine anew the relationship of directed polymers in random media on traditional hypercubic versus hierarchical lattices, with the goal of understanding the dimensionality dependence of the essential scaling index β at the heart of the Kardar-Parisi-Zhang universality class. A seemingly accurate, but entirely empirical, ansatz due to Perlsman and Schwartz, proposed many years ago, can be put in proper context by anchoring the connection between these distinct lattice types at vanishing dimensionality. We graft together complementary perturbative field-theoretic and nonperturbative real-space renormalization group tools to establish the necessary connection, thereby elucidating the central mystery underlying the ansatz's uncanny apparent success, but also revealing its intrinsic limitations. Furthermore, we perform an extensive Euler integration of the KPZ equation in 3+1 dimensions which, bolstered by a separate directed polymer simulation, allows us an estimate for the critical exponent β_{3+1}^{KPZ}=0.1845(4) that greatly improves upon all previous Monte Carlo calculations in this regard and rules out the Perlsman-Schwartz value, 0.1882^{+}, in that dimension. Finally, leveraging this hybrid RG partnership permits us a versatile, more potent, tool to explore the general KPZ problem across dimensions, as well as a conjecture for its key critical exponent, β=1/2-0.22967ɛ, as ɛ→0, testable in a three-loop calculation.

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.

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.

Zero-temperature directed polymer in random potential in 4+1 dimensions

Zero-temperature directed polymer in random potential in 4+1 dimensions is described. The fluctuation ΔE(t) of the lowest energy of the polymer varies as t^{β} with β=0.159±0.007 for polymer length t and ΔE follows ΔE(L)∼L^{α} at saturation with α=0.275±0.009, where L is the system size. The dynamic exponent z≈1.73 is obtained from z=α/β. The estimated values of the exponents satisfy the scaling relation α+z=2 very well. We also monitor the end to end distance of the polymer and obtain z independently. Our results show that the upper critical dimension of the Kardar-Parisi-Zhang equation is higher than d=4+1 dimensions.

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.

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.

Numerical estimate of the Kardar-Parisi-Zhang universality class in (2+1) dimensions

We study the restricted solid on solid model for surface growth in spatial dimension d=2 by means of a multisurface coding technique that allows one to produce a large number of samples in the stationary regime in a reasonable computational time. Thanks to (i) a careful finite-size scaling analysis of the critical exponents and (ii) the accurate estimate of the first three moments of the height fluctuations, we can quantify the wandering exponent with unprecedented precision: χ(d=2)=0.3869(4). This figure is incompatible with the long-standing conjecture due to Kim and Koesterlitz that hypothesized χ(d=2)=2/5.

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.

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

We show that the theoretical machinery developed for the Kardar-Parisi-Zhang (KPZ) class in low dimensions is obeyed by the restricted solid-on-solid model for substrates with dimensions up to d=6. Analyzing different restriction conditions, we show that the height distributions of the interface are universal for all investigated dimensions. It means that fluctuations are not negligible and, consequently, the system is still below the upper critical dimension at d=6. The extrapolation of the data to dimensions d≥7 predicts that the upper critical dimension of the KPZ class is infinite.

Kardar-Parisi-Zhang equation with spatially correlated noise: a unified picture from nonperturbative renormalization group

We investigate the scaling regimes of the Kardar-Parisi-Zhang (KPZ) equation in the presence of spatially correlated noise with power-law decay D(p) ∼ p(-2ρ) in Fourier space, using a nonperturbative renormalization group approach. We determine the full phase diagram of the system as a function of ρ and the dimension d. In addition to the weak-coupling part of the diagram, which agrees with the results from Europhys. Lett. 47, 14 (1999) and Eur. Phys. J. B 9, 491 (1999), we find the two fixed points describing the short-range- (SR) and long-range- (LR) dominated strong-coupling phases. In contrast with a suggestion in the references cited above, we show that, for all values of ρ, there exists a unique strong-coupling SR fixed point that can be continuously followed as a function of d. We show in particular that the existence and the behavior of the LR fixed point do not provide any hint for 4 being the upper critical dimension of the KPZ equation with SR noise.

Restricted solid-on-solid model with a proper restriction parameter N in 4+1 dimensions

A restricted solid-on-solid growth model is studied for various restriction parameters N in d=4+1 dimensions. The interface width W grows as t^{β} with β=0.158 ± 0.006 and W follows W∼L{α} at saturation with α=0.273 ± 0.009, where L is the system size. The dynamic exponent z≈1.73 is obtained from the relation z=α/β. The estimated exponents satisfy the scaling relation α+z=2 very well. Our results indicate that the upper critical dimension of the Kardar-Parisi-Zhang equation is larger than d=4+1 dimensions. With a proper choice of the restriction parameter N, we can reduce the discrete effect of the height to the width and obtain the values of the exponents accurately.

Stiff directed lines in random media

We investigate the localization of stiff directed lines with bending energy by a short-range random potential. We apply perturbative arguments, Flory scaling arguments, a variational replica calculation, and functional renormalization to show that a stiff directed line in 1+d dimensions undergoes a localization transition with increasing disorder for d>2/3. We demonstrate that this transition is accessible by numerical transfer matrix calculations in 1+1 dimensions and analyze the properties of the disorder-dominated phase in detail. On the basis of the two-replica problem, we propose a relation between the localization of stiff directed lines in 1+d dimensions and of directed lines under tension in 1+3d dimensions, which is strongly supported by identical free-energy distributions. This shows that pair interactions in the replicated Hamiltonian determine the nature of directed line localization transitions with consequences for the critical behavior of the Kardar-Parisi-Zhang equation. We support the proposed relation to directed lines via multifractal analysis, revealing an analogous Anderson transition-like scenario and a matching correlation length exponent. Furthermore, we quantify how the persistence length of the stiff directed line is reduced by disorder.