Universality in two-dimensional Kardar-Parisi-Zhang growth

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.

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.

Random deposition with a power-law noise model: Multiaffine analysis

We study the random deposition model with power-law distributed noise and rare-event dominated fluctuation. In this model instead of particles with unit sizes, rods with variable lengths are deposited onto the substrate. The length of each rod is chosen from a power-law distribution P(l)∼l^{-(μ+1)}, and the site at which each rod is deposited is chosen randomly. The results show that for μ<μ_{c}=3 the log-log diagram of roughness, W(t), versus deposition time, t, increases as a step function, where the roughness in each interval acts as W_{loc}(t)≈t^{β_{loc}}. The local growth exponent, β_{loc}, is zero for μ=1. By increasing the μ exponent, the value of β_{loc} is increased. It tends to the growth exponent of the random distribution model with Gaussian noise, β=1/2, at μ_{c}=3. The fractal analysis of the height fluctuations for this model was performed by multifractal detrended fluctuation analysis algorithm. The results show multiaffinity behavior for the height fluctuations at μ<μ_{c} and the multiaffinity strength is greater for smaller values of the μ exponent.

Observation of Two-Dimensional Localized Jones-Roberts Solitons in Bose-Einstein Condensates

Jones-Roberts solitons are the only known class of stable dark solitonic solutions of the nonlinear Schrödinger equation in two and three dimensions. They feature a distinctive elongated elliptical shape that allows them to travel without change of form. By imprinting a triangular phase pattern, we experimentally generate two-dimensional Jones-Roberts solitons in a three-dimensional atomic Bose-Einstein condensate. We monitor their dynamics, observing that this kind of soliton is indeed not affected by dynamic (snaking) or thermodynamic instabilities, that instead make other classes of dark solitons unstable in dimensions higher than one. Our results confirm the prediction that Jones-Roberts solitons are stable solutions of the nonlinear Schrödinger equation and promote them for applications beyond matter wave physics, like energy and information transport in noisy and inhomogeneous environments.

Ordered phases in coupled nonequilibrium systems: Dynamic properties

We study the dynamical properties of the ordered phases obtained in a coupled nonequilibrium system describing advection of two species of particles by a stochastically evolving landscape. The local dynamics of the landscape also gets affected by the particles. In a companion paper we have presented static properties of different phases that arise as the two-way coupling parameters are varied. In this paper we discuss the dynamics. We show that in the ordered phases macroscopic particle clusters move over an ergodic time scale growing exponentially with system size but the ordered landscape shows dynamics over a faster time scale growing as a power of system size. We present a scaling ansatz that describes several dynamical correlation functions of the landscape measured in steady state.

Analysis of etching at a solid-solid interface

We present a method to derive an analytical expression for the roughness of an eroded surface whose dynamics are ruled by cellular automaton. Starting from the automaton, we obtain the time evolution of the height average and height variance (roughness). We apply this method to the etching model in 1+1 dimensions, and then we obtain the roughness exponent. Using this in conjunction with the Galilean invariance we obtain the other exponents, which perfectly match the numerical results obtained from simulations. These exponents are exact, and they are the same as those exhibited by the Kardar-Parisi-Zhang (KPZ) model for this dimension. Therefore, our results provide proof for the conjecture that the etching and KPZ models belong to the same universality class. Moreover, the method is general, and it can be applied to other cellular automata models.

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.

Keldysh field theory for driven open quantum systems

Recent experimental developments in diverse areas-ranging from cold atomic gases to light-driven semiconductors to microcavity arrays-move systems into the focus which are located on the interface of quantum optics, many-body physics and statistical mechanics. They share in common that coherent and driven-dissipative quantum dynamics occur on an equal footing, creating genuine non-equilibrium scenarios without immediate counterpart in equilibrium condensed matter physics. This concerns both their non-thermal stationary states and their many-body time evolution. It is a challenge to theory to identify novel instances of universal emergent macroscopic phenomena, which are tied unambiguously and in an observable way to the microscopic drive conditions. In this review, we discuss some recent results in this direction. Moreover, we provide a systematic introduction to the open system Keldysh functional integral approach, which is the proper technical tool to accomplish a merger of quantum optics and many-body physics, and leverages the power of modern quantum field theory to driven open quantum systems.

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.

Relaxation after a change in the interface growth dynamics

The global effects of sudden changes in the interface growth dynamics are studied using models of the Edwards-Wilkinson (EW) and Kardar-Parisi-Zhang (KPZ) classes during their growth regimes in dimensions d=1 and d=2. Scaling arguments and simulation results are combined to predict the relaxation of the difference in the roughness of the perturbed and the unperturbed interfaces, ΔW^{2}∼s{c}t{-γ}, where s is the time of the change and t>s is the observation time after that event. The previous analytical solution for the EW-EW changes is reviewed and numerically discussed in the context of lattice models, with possible decays with γ=3/2 and γ=1/2. Assuming the dominant contribution to ΔW{2} to be predicted from a time shift in the final growth dynamics, the scaling of KPZ-KPZ changes with γ=1-2β and c=2β is predicted, where β is the growth exponent. Good agreement with simulation results in d=1 and d=2 is observed. A relation with the relaxation of a local autoresponse function in d=1 cannot be discarded, but very different exponents are shown in d=2. We also consider changes between different dynamics, with the KPZ-EW as a special case in which a faster growth, with dynamical exponent z_{i}, changes to a slower one, with exponent z. A scaling approach predicts a crossover time t_{c}∼s{z/z_{i}}≫s and ΔW{2}∼s{c}F(t/t_{c}), with the decay exponent γ=1/2 of the EW class. This rules out the simplified time shift hypothesis in d=2 dimensions. These results help to understand the remarkable differences in EW smoothing of correlated and uncorrelated surfaces, and the approach may be extended to sudden changes between other growth dynamics.

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.

Extremal paths, the stochastic heat equation, and the three-dimensional Kardar-Parisi-Zhang universality class

Following our numerical work [Phys. Rev. Lett. 109, 170602 (2012)] focused upon the 2+1 Kardar-Parisi-Zhang (KPZ) equation with flat initial condition, we return here to study, in depth, the three-dimensional (3D) radial KPZ problem, comparing common scaling phenomena exhibited by the pt-pt directed polymer in a random medium (DPRM), the stochastic heat equation (SHE) with multiplicative noise in three dimensions, and kinetic roughening phenomena associated with 3D Eden clusters. Examining variants of the 3D DPRM, as well as numerically integrating, via the Itô prescription, the constrained SHE for different values of the KPZ coupling, we provide strong evidence for universality within this 3D KPZ class, revealing shared values for the limit distribution skewness and kurtosis, along with universal first and second moments. Our numerical analysis of the 3D SHE, well flanked by the DPRM results, appears without precedent in the literature. We consider, too, the 2+1 KPZ equation in the deeply evolved kinetically roughened stationary state, extracting the essential limit distribution characterizing fluctuations therein, revealing a higher-dimensional relative of the 1+1 KPZ Baik-Rains distribution. Complementary, corroborative findings are provided via the Gaussian DPRM, as well as the restricted-solid-on-solid model of stochastic growth, stalwart members of the 2+1 KPZ class. Next, contact is made with a recent nonperturbative, field-theoretic renormalization group calculation for the key universal amplitude ratio in this context. Finally, in the crossover from transient to stationary-state statistics, we observe a higher dimensional manifestation of the skewness minimum discovered by Takeuchi [Phys. Rev. Lett. 110, 210604 (2013)] in 1+1 KPZ class liquid-crystal experiments.

Crystal growth inside an octant

We study crystal growth inside an infinite octant on a cubic lattice. The growth proceeds through the deposition of elementary cubes into inner corners. After rescaling by the characteristic size, the interface becomes progressively more deterministic in the long-time limit. Utilizing known results for the crystal growth inside a two-dimensional corner, we propose a hyperbolic partial differential equation for the evolution of the limiting shape. This equation is interpreted as a Hamilton-Jacobi equation, which helps in finding an analytical solution. Simulations of the growth process are in excellent agreement with analytical predictions. We then study the evolution of the subleading correction to the volume of the crystal, the asymptotic growth of the variance of the volume of the crystal, and the total number of inner and outer corners. We also show how to generalize the results to arbitrary spatial dimension.

Kardar-Parisi-Zhang universality class in (2+1) dimensions: universal geometry-dependent distributions and finite-time corrections

The dynamical regimes of models belonging to the Kardar-Parisi-Zhang (KPZ) universality class are investigated in d=2+1 by extensive simulations considering flat and curved geometries. Geometry-dependent universal distributions, different from their Tracy-Widom counterpart in one dimension, were found. Distributions exhibit finite-time corrections hallmarked by a shift in the mean decaying as t(-β), where β is the growth exponent. Our results support a generalization of the ansatz h=v(∞)t+(Γt)(β)χ+η+ζt(-β) to higher dimensions, where v(∞), Γ, ζ, and η are nonuniversal quantities whereas β and χ are universal and the last one depends on the surface geometry. Generalized Gumbel distributions provide very good fits of the distributions in at least four orders of magnitude around the peak, which can be used for comparisons with experiments. Our numerical results call for analytical approaches and experimental realizations of the KPZ class in two-dimensional systems.

Nonperturbative renormalization group for the stationary Kardar-Parisi-Zhang equation: scaling functions and amplitude ratios in 1+1, 2+1, and 3+1 dimensions

We investigate the strong-coupling regime of the stationary Kardar-Parisi-Zhang equation for interfaces growing on a substrate of dimension d = 1, 2, and 3 using a nonperturbative renormalization group (NPRG) approach. We compute critical exponents, correlation and response functions, extract the related scaling functions, and calculate universal amplitude ratios. We work with a simplified implementation of the second-order (in the response field) approximation proposed in a previous work [Phys. Rev. E 84, 061150 (2011) and Phys. Rev. E 86, 019904(E) (2012)], which greatly simplifies the frequency sector of the NPRG flow equations, while keeping a nontrivial frequency dependence for the two-point functions. The one-dimensional scaling function obtained within this approach compares very accurately with the scaling function obtained from the full second-order NPRG equations and with the exact scaling function. Furthermore, the approach is easily applicable to higher dimensions and we provide scaling functions and amplitude ratios in d = 2 and d = 3. We argue that our ansatz is reliable up to d [Symbol: see text] 3.5.

Growth inside a corner: the limiting interface shape

We investigate the growth of a crystal that is built by depositing cubes inside a corner. The interface of this crystal approaches a deterministic growing limiting shape in the long-time limit. Building on known results for the corresponding two-dimensional system and accounting for basic three-dimensional symmetries, we conjecture a governing equation for the evolution of the interface profile. We solve this equation analytically and find excellent agreement with simulations of the growth process. We also present a generalization to arbitrary spatial dimension.

Extremely large-scale simulation of a Kardar-Parisi-Zhang model using graphics cards

The octahedron model introduced recently has been implemented onto graphics cards, which permits extremely large-scale simulations via binary lattice gases and bit-coded algorithms. We confirm scaling behavior belonging to the two-dimensional Kardar-Parisi-Zhang universality class and find a surface growth exponent: β = 0.2415(15) on 2(17) × 2(17) systems, ruling out β = 1/4 suggested by field theory. The maximum speedup with respect to a single CPU is 240. The steady state has been analyzed by finite-size scaling and a growth exponent α = 0.393(4) is found. Correction-to-scaling-exponent are computed and the power-spectrum density of the steady state is determined. We calculate the universal scaling functions and cumulants and show that the limit distribution can be obtained by the sizes considered. We provide numerical fitting for the small and large tail behavior of the steady-state scaling function of the interface width.

Anomalous roughening in competitive growth models with time-decreasing rates of correlated dynamics

Lattice growth models where uncorrelated random deposition competes with some aggregation dynamics that generates correlations are studied with rates of the correlated component decreasing as a power law. These models have anomalous roughening, with anomalous exponents related to the normal exponents of the correlated dynamics, to an exponent characterizing the aggregation mechanism and to that power-law exponent. This is shown by a scaling approach extending the Family-Vicsek relation previously derived for the models with time-independent rates, thus providing a connection of normal and anomalous growth models. Simulation results for several models support those conclusions. Remarkable anomalous effects are observed even for slowly decreasing rates of the correlated component, which may correspond to feasible temperature changes in systems with activated dynamics. The scaling exponents of the correlated component can be obtained only from the estimates of three anomalous exponents, without knowledge of the aggregation mechanism, and a possible application is discussed. For some models, the corresponding Edwards-Wilkinson and Kardar-Parisi-Zhang equations are also discussed.

Morphologies from slippery ballistic deposition model: a bottom-up approach for nanofabrication

We report pattern formation using a slippery ballistic deposition (SBD) model where growth germinates from a single site or from sites distributed periodically on a lattice. By changing the sticking probability p(s) and choosing systems with different lattice constants and symmetries, we demonstrate that a variety of patterns can be generated. These patterns can be further used as scaffolds for nanofabrication. We also demonstrate that by choosing a lateral sticking probability p(l) at the base that is different than p(s), one can control both the early and late time morphologies originating from a seed. Furthermore, we indicate a possible generalization of preparing patterns to higher dimensions that in principle can have potential technological applications for preparing grooves and scaffolds of specific shapes and periodicities.

Large-scale simulations of ballistic deposition: the approach to asymptotic scaling

Extensive kinetic Monte Carlo simulations are presented for ballistic deposition (BD) in (1+1) dimensions. Asymptotic scaling is found only for lattice sizes L≳2¹². Such a large system size for the onset of scaling explains the widespread discrepancies of previous reports for exponents of BD in one and likely higher dimensions. The exponents obtained from our simulations, α=0.499±0.004 and β=0.336±0.004, are in excellent agreement with the exact values α=½ and β=1/3 for the one-dimensional Kardar-Parisi-Zhang equation. Our findings make possible a more informed exploration of exponents for BD in higher dimensions, accurate estimates of which have proven to be elusive.

Classification of (2+1) -dimensional growing surfaces using Schramm-Loewner evolution

Statistical behavior and scaling properties of isoheight lines in three different saturated two-dimensional grown surfaces with controversial universality classes are investigated using ideas from Schramm-Loewner evolution (SLE_{κ}). We present some evidence that the isoheight lines in the ballistic deposition (BD), Eden and restricted solid-on-solid (RSOS) models have conformally invariant properties all in the same universality class as the self-avoiding random walk (SAW), equivalently SLE_{8/3}. This leads to the conclusion that all these discrete growth models fall into the same universality class as the Kardar-Parisi-Zhang (KPZ) equation in two dimensions.

Conformal invariance in (2+1)-dimensional stochastic systems

Stochastic partial differential equations can be used to model second-order thermodynamical phase transitions, as well as a number of critical out-of-equilibrium phenomena. In (2+1) dimensions, many of these systems are conjectured (and some are indeed proved) to be described by conformal field theories. We advance, in the framework of the Martin-Siggia-Rose field-theoretical formalism of stochastic dynamics, a general solution of the translation Ward identities, which yields a putative conformal energy-momentum tensor. Even though the computation of energy-momentum correlators is obstructed, in principle, by dimensional reduction issues, these are bypassed by the addition of replicated fields to the original (2+1)-dimensional model. The method is illustrated with an application to the Kardar-Parisi-Zhang (KPZ) model of surface growth. The consistency of the approach is checked by means of a straightforward perturbative analysis of the KPZ ultraviolet region, leading, as expected, to its c=1 conformal fixed point.

Directed d -mer diffusion describing the Kardar-Parisi-Zhang-type surface growth

We show that d+1-dimensional surface growth models can be mapped onto driven lattice gases of d-mers. The continuous surface growth corresponds to one dimensional drift of d-mers perpendicular to the (d-1-dimensional "plane" spanned by the d-mers. This facilitates efficient bit-coded algorithms with generalized Kawasaki dynamics of spins. Our simulations in d=2, 3, 4, 5 dimensions provide scaling exponent estimates on much larger system sizes and simulations times published so far, where the effective growth exponent exhibits an increase. We provide evidence for the agreement with field theoretical predictions of the Kardar-Parisi-Zhang universality class and numerical results. We show that the (2+1)-dimensional exponents conciliate with the values suggested by Lässig within error margin, for the largest system sizes studied here, but we cannot support his predictions for (3+1)d numerically.

Scaling of clusters and winding-angle statistics of isoheight lines in two-dimensional Kardar-Parisi-Zhang surfaces

We investigate the statistics of isoheight lines of (2+1) -dimensional Kardar-Parisi-Zhang model at different level sets around the mean height in the saturation regime. We find that the exponent describing the distribution of the height-cluster size behaves differently for level cuts above and below the mean height, while the fractal dimensions of the height-clusters and their perimeters remain unchanged. The statistics of the winding angle confirms the previous observation that these contour lines are in the same universality class as self-avoiding random walks.

Mapping of (2+1) -dimensional Kardar-Parisi-Zhang growth onto a driven lattice gas model of dimers

We show that a (2+1) -dimensional discrete surface growth model exhibiting Kardar-Parisi-Zhang (KPZ) class scaling can be mapped onto a two-dimensional conserved lattice gas model of directed dimers. The KPZ height anisotropy in the surface model corresponds to a driven diffusive motion of the lattice gas dimers. We confirm by numerical simulations that the scaling exponents of the dimer model are in agreement with those of the (2+1) -dimensional KPZ class. This opens up the possibility of analyzing growth models via reaction-diffusion models, which allow much more efficient computer simulations.

Maximal- and minimal-height distributions of fluctuating interfaces

Maximal- and minimal-height distributions (MAHD, MIHD) of two-dimensional interfaces grown with the nonlinear equations of Kardar-Parisi-Zhang (KPZ, second order) and of Villain-Lai-Das Sarma (VLDS, fourth order) are shown to be different. Two universal curves may be MAHD or MIHD of each class depending on the sign of the relevant nonlinear term, which is confirmed by results of several lattice models in the KPZ and VLDS classes. The difference between MAHD and MIDH is connected with the asymmetry of the local height distribution. A simple, exactly solvable deposition-erosion model is introduced to illustrate this feature. The average extremal heights scale with the same exponent of the average roughness. In contrast to other correlated systems, generalized Gumbel distributions do not fit those MAHD and MIHD, nor those of Edwards-Wilkinson growth.

Numerical study of the Kardar-Parisi-Zhang equation

We integrate numerically the Kardar-Parisi-Zhang (KPZ) equation in 1+1 and 2+1 dimensions using a Euler discretization scheme and the replacement of (nablah)(2) by exponentially decreasing functions of that quantity to suppress instabilities. When applied to the equation in 1+1 dimensions, the method of instability control provides values of scaling amplitudes consistent with exactly known results, in contrast to the deviations generated by the original scheme. In 2+1 dimensions, we spanned a range of the model parameters where transients with Edwards-Wilkinson or random growth are not observed, in box sizes 8< or =L< or =128 . We obtain a roughness exponent of 0.37< or =alpha< or =0.40 and steady state height distributions with skewness S=0.25+/-0.01 and kurtosis Q=0.15+/-0.1 . These estimates are obtained after extrapolations to the large L limit, which is necessary due to significant finite-size effects in the estimates of effective exponents and height distributions. On the other hand, the steady state roughness distributions show weak scaling corrections and evidence of stretched exponential tails. These results confirm previous estimates from lattice models, showing their reliability as representatives of the KPZ class.

Finite-size effects in roughness distribution scaling

We study numerically finite-size corrections in scaling relations for roughness distributions of various interface growth models. The most common relation, which considers the average roughness <w(2)> as a scaling factor, is not obeyed in the steady states of a group of ballisticlike models in 2+1 dimensions, even when very large system sizes are considered. On the other hand, good collapse of the same data is obtained with a scaling relation that involves the root mean square fluctuation of the roughness, which can be explained by finite-size effects on second moments of the scaling functions. We also obtain data collapse with an alternative scaling relation that accounts for the effect of the intrinsic width, which is a constant correction term previously proposed for the scaling of <w(2)> . This illustrates how finite-size corrections can be obtained from roughness distributions scaling. However, we discard the usual interpretation that the intrinsic width is a consequence of high surface steps by analyzing data of restricted solid-on-solid models with various maximal height differences between neighboring columns. We also observe that large finite-size corrections in the roughness distributions are usually accompanied by huge corrections in height distributions and average local slopes, as well as in estimates of scaling exponents. The molecular-beam epitaxy model of Das Sarma and Tamborenea in 1+1 dimensions is a case example in which none of the proposed scaling relations work properly, while the other measured quantities do not converge to the expected asymptotic values. Thus although roughness distributions are clearly better than other quantities to determine the universality class of a growing system, it is not the final solution for this task.

Conformal curves on the WO3 surface

We have studied the isoheight lines on the WO3 surface as a physical candidate for conformally invariant curves. We have shown that these lines are conformally invariant with the same statistics of domain walls in the critical Ising model. They belong to the family of conformal invariant curves called Schramm-Loewner evolution (or SLE(kappa)), with diffusivity of kappa approximately 3. This can be regarded as the first experimental observation of SLE curves. We have also argued that Ballistic Deposition (BD) can serve as a growth model giving rise to contours with similar statistics at large scales.

Scaling behavior in corrosion and growth of a passive film

We study a simple model for metal corrosion controlled by the reaction rate of the metal with an anionic species and the diffusion of that species in the growing passive film between the solution and the metal. A crossover from the reaction-controlled to the diffusion-controlled growth regime with different roughening properties is observed. Scaling arguments provide estimates of the crossover time and film thickness as functions of the reaction and diffusion rates and the concentration of anionic species in the film-solution interface, including a nontrivial square-root dependence on that concentration. At short times, the metal-film interface exhibits Kardar-Parisi-Zhang (KPZ) scaling, which crosses over to a diffusion-limited erosion (Laplacian growth) regime at long times. The roughness of the metal-film interface at long times is obtained as a function of the rates of reaction and diffusion and of the KPZ growth exponent. The predictions have been confirmed by simulations of a lattice version of the model in two dimensions. Relations with other erosion and corrosion models and possible applications are discussed.

Surface and bulk properties of deposits grown with a bidisperse ballistic deposition model

We study roughness scaling of the outer surface and the internal porous structure of deposits generated with the three-dimensional bidisperse ballistic deposition (BBD), in which particles of two sizes are randomly deposited. Systematic extrapolation of roughness and dynamical exponents and the comparison of roughness distributions indicate that the top surface has Kardar-Parisi-Zhang (KPZ) scaling for any ratio F of the flux between large and small particles. A scaling theory predicts the characteristic time of the crossover from random to correlated growth in BBD and provides relations between the amplitudes of roughness scaling and F in the KPZ regime. The porosity of the deposits monotonically increases with F and scales as F{12} for small F, which is also explained by the scaling approach and illustrates the possibility of connecting surface growth rules and bulk properties. The suppression of relaxation mechanisms in BBD enhances the connectivity of the deposits when compared to other ballisticlike models, so that they percolate down to F approximately 0.05. The fractal dimension of the internal surface of the percolating deposits is D{F} approximately 2.9, which is very close to the values in other ballistic-like models and suggests universality among these systems.

Universal and nonuniversal features in the crossover from linear to nonlinear interface growth

We study a restricted solid-on-solid model involving deposition and evaporation with probabilities p and 1 - p, respectively, in one-dimensional substrates. It presents a crossover from Edwards-Wilkinson (EW) to Kardar-Parisi-Zhang (KPZ) scaling for p approximately 0.5. The associated KPZ equation is analytically derived, exhibiting a coefficient lambda of the nonlinear term proportional to q identical with p - 1/2, which is confirmed numerically by calculation of tilt-dependent growth velocities for several values of p. This linear lambda - q relation contrasts to the apparently universal parabolic law obtained in competitive models mixing EW and KPZ components. The regions where the interface roughness shows pure EW and KPZ scaling are identified for 0.55< or =p< or =0.8, which provides numerical estimates of the crossover times tc. They scale as tc approximately lambda -phi with phi=4.1+/-0.1, which is in excellent agreement with the theoretically predicted universal value phi=4 and improves previous numerical estimates, which suggested phi approximately 3.

Scaling of ballistic deposition from a Langevin equation

An exact lattice Langevin equation is derived for the ballistic deposition model of surface growth. The continuum limit of this equation is dominated by the Kardar-Parisi-Zhang (KPZ) equation at all length and time scales. For a one-dimensional substrate the solution of the exact lattice Langevin equation yields the KPZ scaling exponents without any extrapolation. For a two-dimensional substrate the scaling exponents are different from those found from computer simulations. This discrepancy is discussed in relation to analytic approaches to the KPZ equation in higher dimensions.

Scaling in the crossover from random to correlated growth

In systems where deposition rates are high compared to diffusion, desorption, and other mechanisms that generate correlations, a crossover from random to correlated growth of surface roughness is expected at a characteristic time t0. This crossover is analyzed in lattice models via scaling arguments, with support from simulation results presented here and in other works. We argue that the amplitudes of the saturation roughness and of the saturation time t(x) scale as t0(1/2) and t0, respectively. For models with lateral aggregation, which typically are in the Kardar-Parisi-Zhang (KPZ) class, we show that t0 approximately p(-1), where p is the probability of the correlated aggregation mechanism to take place. However, t0 approximately p(-2) is obtained in solid-on-solid models with single-particle deposition attempts. This group includes models in various universality classes, with numerical examples being provided in the Edwards-Wilkinson (EW), KPZ, and Villain-Lai-Das Sarma (nonlinear molecular-beam epitaxy) classes. Most applications are for two-component models in which random deposition, with probability 1-p, competes with a correlated aggregation process with probability p. However, our approach can be extended to other systems with the same crossover, such as the generalized restricted solid-on-solid model with maximum height difference S, for large S. Moreover, the scaling approach applies to all dimensions. In the particular case of one-dimensional KPZ processes with this crossover, we show that t0 approximately nu(-1) and nu approximately lambda(2/3), where nu and lambda are the coefficients of the linear and nonlinear terms of the associated KPZ equations. The applicability of previous results to models in the EW and KPZ classes is discussed.

Numerical study of roughness distributions in nonlinear models of interface growth

We analyze the shapes of roughness distributions of discrete models in the Kardar, Parisi, and Zhang (KPZ) and in the Villain, Lai, and Das Sarma (VLDS) classes of interface growth, in one and two dimensions. Three KPZ models in d=2 confirm the expected scaling of the distribution and show a stretched exponential tail approximately as exp(-x0.8), with a significant asymmetry near the maximum. Conserved restricted solid-on-solid models belonging to the VLDS class were simulated in d=1 and d=2. The tail in d=1 has the form exp(-x2) and, in d=2, has a simple exponential decay, but is quantitatively different from the distribution of the linear fourth-order (Mullins-Herring) theory. It is not possible to fit any of the above distributions to those of 1/f(alpha) noise interfaces, in contrast with recently studied models with depinning transitions.

Roughness of two-dimensional surfaces with global constraints

We study dynamical scaling properties of the two-dimensional surface growth models with global constraints. These include the growth model from a partition function Z = sigma{(h(r)Pi(h = h)(min) (h(max)) 1/2 (1 + z (n(h))), multiparticle-correlated surface growth models and dissociative Q-mer growth models. The equilibrium surfaces of all the models except the dimer model show the same dynamical scaling behavior W2 (L,t) = (1/2piK(G)) ln [L g (t/L(z(W)))] with z(W) = 2.5 and K(G) = 0.916 , whereas the surface in the dimer model has a correction to the scaling. The growing (eroding) surfaces have two phases. The models with z > or = 0 show the normal Kardar-Parisi-Zhang scaling behavior. In contrast the models with -1 < or = z < 0 and multiparticle-correlated growth model manifest grooved surface structures with alpha = 1. The growing surfaces of Q -mer models form rather complex facets.

What is the connection between ballistic deposition and the Kardar-Parisi-Zhang equation?

Ballistic deposition (BD) is believed to belong to the Kardar-Parisi-Zhang (KPZ) universality class. In this paper we study the validity of this belief by rigorously deriving a continuum equation from the BD microscopic rules, which deviates from the KPZ equation. We show that in one dimension and in the presence of noise the deviation is not important. This is not the case in the absence of noise. In more than one dimension and in the presence of noise we obtain an equation that superficially seems to be a continuum equation but in which the symmetry under rotations around the growth direction is broken.

Numerical study of discrete models in the class of the nonlinear molecular beam epitaxy equation

We study numerically some discrete growth models belonging to the class of the nonlinear molecular beam epitaxy equation, or the Villain-Lai-Das Sarma (VLDS) equation. The conserved restricted solid-on-solid model (CRSOS) with maximum height differences Delta H(max)=1 and Delta H(max)=2 was analyzed in substrate dimensions d=1 and d=2 . The Das Sarma and Tamborenea (DT) model and a competitive model involving random deposition and CRSOS deposition were studied in d=1. For the CRSOS model with Delta H(max)=1, we obtain the more accurate estimates of scaling exponents in d=1:roughness exponent alpha=0.94+/-0.02 and dynamical exponent z=2.88+/-0.04. These estimates are significantly below the values of one-loop renormalization for the VLDS theory, which confirms Janssen's proposal of the existence of higher-order corrections. The roughness exponent in d=2 is very near the one-loop result alpha=2/3, in agreement with previous works. The moments W(n) of orders n=2 , 3, 4 of the height distribution were calculated for all models, and the skewness S triple bond W3/W(3/2)(2) and the kurtosis Q triple bond W4/W(2)2-3 were estimated. At the steady states, the CRSOS models and the competitive model have nearly the same values of S and Q in d=1, which suggests that these amplitude ratios are universal in the VLDS class. The estimates for the DT model are different, possibly due to their typically long crossover to asymptotic values. Results for the CRSOS models in d=2 also suggest that those quantities are universal.