Ballistic deposition on surfaces

Spectral boundary of the asymmetric simple exclusion process: Free fermions, Bethe ansatz, and random matrix theory

In nonequilibrium statistical mechanics, the asymmetric simple exclusion process (ASEP) serves as a paradigmatic example. We investigate the spectral characteristics of the ASEP, focusing on the spectral boundary of its generator matrix. We examine finite ASEP chains of length L, under periodic boundary conditions (PBCs) and open boundary conditions (OBCs). Notably, the spectral boundary exhibits L spikes for PBCs and L+1 spikes for OBCs. Treating the ASEP generator as an interacting non-Hermitian fermionic model, we extend the model to have tunable interaction. In the noninteracting case, the analytically computed many-body spectrum shows a spectral boundary with prominent spikes. For PBCs, we use the coordinate Bethe ansatz to interpolate between the noninteracting case to the ASEP limit and show that these spikes stem from clustering of Bethe roots. The robustness of the spikes in the spectral boundary is demonstrated by linking the ASEP generator to random matrices with trace correlations or, equivalently, random graphs with distinct cycle structures, both displaying similar spiked spectral boundaries.

Finite-size and finite-time scaling for kinetic rough interfaces

We consider discrete models of kinetic rough interfaces that exhibit space-time scale invariance in height-height correlation. We use the generic scaling theory of Ramasco et al. [Phys. Rev. Lett. 84, 2199 (2000)0031-900710.1103/PhysRevLett.84.2199] to confirm that the dynamical structure factor of the height profile can uniquely characterize the underlying dynamics. We apply both finite-size and finite-time scaling methods that systematically allow an estimation of the critical exponents and the scaling functions, eventually establishing the universality class accurately. The finite-size scaling analysis offers an alternative way to characterize the anomalous rough interfaces. As an illustration, we investigate a class of self-organized interface models in random media with extremal dynamics. The isotropic version shows a faceted pattern and belongs to the same universality class (as shown numerically) as the Sneppen model (version A). We also examine an anisotropic version of the Sneppen model and suggest that the model belongs to the universality class of the tensionless Kardar-Parisi-Zhang (tKPZ) equation in one 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.

Logarithmic or algebraic: Roughening of an active Kardar-Parisi-Zhang surface

The Kardar-Parisi-Zhang (KPZ) equation sets the universality class for growing and roughening of nonequilibrium surfaces without any conservation law and nonlocal effects. We argue here that the KPZ equation can be generalized by including a symmetry-permitted nonlocal nonlinear term of active origin that is of the same order as the one included in the KPZ equation. Including this term, the 2D active KPZ equation is stable in some parameter regimes, in which the interface conformation fluctuations exhibit sublogarithmic or superlogarithmic roughness, with nonuniversal exponents, giving positional generalized quasi-long-ranged order. For other parameter choices, the model is unstable, suggesting a perturbatively inaccessible algebraically rough interface or positional short-ranged order. Our model should serve as a paradigmatic nonlocal growth equation.

A triangular model of fractal growth with application to adsorptive spin-coating of polymers

Over the last 40 years, applied mathematicians and physicists have proposed a number of mathematical models that produce structures exhibiting a fractal dimension. This work has coincided with the discovery that objects with fractal dimension are relatively common in the natural and human-produced worlds. One particularly successful model of fractal growth is the diffusion limited aggregation (DLA) model, a model as notable for its simplicity as for its complex and varied behavior. It has been modified and used to simulate fractal growth processes in numerous experimental and empirical contexts. In this work, we present an alternative fractal growth model that is based on a growing mass that bonds to particles in a surrounding medium and then exerts a force on them in an iterative process of growth and contraction. The resulting structure is a spreading triangular network rather than an aggregate of spheres, and the model is conceptually straightforward. To the best of our knowledge, this model is unique and differs in its dynamics and behavior from the DLA model and related particle aggregation models. We explore the behavior of the model, demonstrate the range of model output, and show that model output can have a variable fractal dimension between 1.5 and 1.83 that depends on model parameters. We also apply the model to simulating the development of polymer thin films prepared using spin-coating which also exhibit variable fractal dimensions. We demonstrate how the model can be adjusted to different dewetting conditions as well as how it can be used to simulate the modification of the polymer morphology under solvent annealing.

Current fluctuations in a partially asymmetric simple exclusion process with a defect particle

We study an exclusion process on a ring comprising a free defect particle in a bath of normal particles. The model is one of the few integrable cases in which the bath particles are partially asymmetric. The presence of the free defect creates localized or shock phases according to parameter values. We use a functional approach to Bethe equations resulting from a nested Bethe ansatz to calculate exactly the mean currents and diffusion constants. The results agree very well with Monte Carlo simulations and reveal the main modes of fluctuation in the different phases of the steady state.

Marangoni spreading on liquid substrates in new media art

With the advent of new media art, artists have harnessed fluid dynamics to create captivating visual narratives. A striking technique known as dendritic painting employs mixtures of ink and isopropanol atop paint, yielding intricate tree-like patterns. To unravel the intricacies of that technique, we examine the spread of ink/alcohol droplets over liquid substrates with diverse rheological properties. On Newtonian substrates, the droplet size evolution exhibits two power laws, suggesting an underlying interplay between viscous and Marangoni forces. The leading edge of the droplet spreads as a precursor film with an exponent of 3/8, while its main body spreads with an exponent of 1/4. For a weakly shear-thinning acrylic resin substrate, the same power laws persist, but dendritic structures emerge, and the texture of the precursor film roughens. The observed roughness and growth exponents (3/4 and 3/5) suggest a connection to the quenched Kardar-Parisi-Zhang universality class, hinting at the existence of quenched disorder in the liquid substrate. Mixing the resin with acrylic paint renders it more viscous and shear-thinning, refining the dendrite edges and further roughening the precursor film. At larger paint concentrations, the substrate becomes a power-law fluid. The roughness and growth exponents then approach 1/2 and 3/4, respectively, deviating from known universality classes. The ensuing structures have a fractal dimension of 1.68, characteristic of diffusion-limited aggregation. These findings underscore how the nonlinear rheological properties of the liquid substrate, coupled with the Laplacian nature of Marangoni spreading, can overshadow the local kinetic roughening of the droplet interface.

Physical properties of a generalized model of multilayer adsorption of dimers

We investigate the transport properties of a complex porous structure with branched fractal architectures formed due to the gradual deposition of dimers in a model of multilayer adsorption. We thoroughly study the interplay between the orientational anisotropy parameter p_{0} of deposited dimers and the formation of porous structures, as well as its impact on the conductivity of the system, through extensive numerical simulations. By systematically varying the value of p_{0}, several critical and off-critical scaling relations characterizing the behavior of the system are examined. The results demonstrate that the degree of orientational anisotropy of dimers plays a significant role in determining the structural and physical characteristics of the system. We find that the Einstein relation relating to the size scaling of the electrical conductance holds true only in the limiting case of p_{0}→1. Monitoring the fractal dimension of the interface of the multilayer formation for various p_{0} values, we reveal that in a wide range of p_{0}>0.2 interface shows the characteristic of a self-avoiding random walk, compared to the limiting case of p_{0}→0 where it is characterized by the fractal dimension of the backbone of ordinary percolation cluster at criticality. Our results thus can provide useful information about the fundamental mechanisms underlying the formation and behavior of wide varieties of amorphous and disordered systems that are of paramount importance both in science and technology as well as in environmental studies.

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.

Accessibility of the surface fractal dimension during film growth

Fractal properties on self-affine surfaces of films growing under nonequilibrium conditions are important in understanding the corresponding universality class. However, measurement of the surface fractal dimension has been intensively investigated and is still very problematic. In this work, we report the behavior of the effective fractal dimension in the context of film growth involving lattice models believed to belong to the Kardar-Parisi-Zhang (KPZ) universality class. Our results, which are presented for growth in a d-dimensional substrate (d=1,2) and use the three-point sinuosity (TPS) method, show universal scaling of the measure M, which is defined in terms of discretization of the Laplacian operator applied to the height of the film surface, M=t^{δ}g[Θ], where t is the time, g[Θ] is a scale function, δ=2β, Θ≡τt^{-1/z}, β, and z are the KPZ growth and dynamical exponents, respectively, and τ is a spatial scale length used to compute M. Importantly, we show that the effective fractal dimensions are consistent with the expected KPZ dimensions for d=1,2, if Θ≲0.3, which include a thin film regime for the extraction of the fractal dimension. This establishes the scale limits in which the TPS method can be used to accurately extract effective fractal dimensions that are consistent with those expected for the corresponding universality class. As a consequence, for the steady state, which is inaccessible to experimentalists studying film growth, the TPS method provided effective fractal dimension consistent with the KPZ ones for almost all possible τ, i.e., 1≲τ<L/2, where L is the lateral size of the substrate on which the deposit is grown. In the growth of thin films, the true fractal dimension can be observed in a narrow range of τ, the upper limit of which is of the same order of magnitude as the correlation length of the surface, indicating the limits of self-affinity of a surface in an experimentally accessible regime. This upper limit was comparatively lower for the Higuchi method or the height-difference correlation function. Scaling corrections for the measure M and the height-difference correlation function are studied analytically and compared for the Edwards-Wilkinson class at d=1, yielding similar accuracy for both methods. Importantly, we extend our discussion to a model representing diffusion-dominated growth of films and find that the TPS method achieves the corresponding fractal dimension only at steady state and in a narrow range of the scale length, compared to that found for the KPZ class.

Kardar-Parisi-Zhang Physics and Phase Transition in a Classical Single Random Walker under Continuous Measurement

We introduce and study a new model consisting of a single classical random walker undergoing continuous monitoring at rate γ on a discrete lattice. Although such a continuous measurement cannot affect physical observables, it has a nontrivial effect on the probability distribution of the random walker. At small γ, we show analytically that the time evolution of the latter can be mapped to the stochastic heat equation. In this limit, the width of the log-probability thus follows a Family-Vicsek scaling law, N^{α}f(t/N^{α/β}), with roughness and growth exponents corresponding to the Kardar-Parisi-Zhang (KPZ) universality class, i.e., α_{KPZ}^{1D}=1/2 and β_{KPZ}^{1D}=1/3, respectively. When γ is increased outside this regime, we find numerically in 1D a crossover from the KPZ class to a new universality class characterized by exponents α_{M}^{1D}≈1 and β_{M}^{1D}≈1.4. In 3D, varying γ beyond a critical value γ_{M}^{c} leads to a phase transition from a smooth phase that we identify as the Edwards-Wilkinson class to a new universality class with α_{M}^{3D}≈1.

Hierarchical deposition and scale-free networks: A visibility algorithm approach

The growth of an interface formed by the hierarchical deposition of particles of unequal size is studied in the framework of a dynamical network generated by a horizontal visibility algorithm. For a deterministic model of the deposition process, the resulting network is scale free with dominant degree exponent γ_{e}=ln3/ln2 and transient exponent γ_{o}=1. An exact calculation of the network diameter and clustering coefficient reveals that the network is scale invariant and inherits the modular hierarchical nature of the deposition process. For the random process, the network remains scale free, where the degree exponent asymptotically converges to γ=3, independent of the system parameters. This result shows that the model is in the class of fractional Gaussian noise through the relation between the degree exponent and the series' Hurst exponent H. Finally, we show through the degree-dependent clustering coefficient C(k) that the modularity remains present in the system.

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.

Universality and crossover behavior of single-step growth models in 1+1 and 2+1 dimensions

We study the kinetic roughening of the single-step (SS) growth model with a tunable parameter p in 1+1 and 2+1 dimensions by performing extensive numerical simulations. We show that there exists a very slow crossover from an intermediate regime dominated by the Edwards-Wilkinson class to an asymptotic regime dominated by the Kardar-Parisi-Zhang (KPZ) class for any p<1/2. We also identify the crossover time, the nonlinear coupling constant, and some nonuniversal parameters in the KPZ equation as a function p. The effective nonuniversal parameters are continuously decreasing with p but not in a linear fashion. Our results provide complete and conclusive evidence that the SS model for p≠1/2 belongs to the KPZ universality class in 2+1 dimensions.

Transverse temperature interfaces in the Katz-Lebowitz-Spohn driven lattice gas

We explore the intriguing spatial patterns that emerge in a two-dimensional spatially inhomogeneous Katz-Lebowitz-Spohn (KLS) driven lattice gas with attractive nearest-neighbor interactions. The domain is split into two regions with hopping rates governed by different temperatures T>T_{c} and T_{c}, respectively, where T_{c} indicates the critical temperature for phase ordering, and with the temperature boundaries oriented perpendicular to the drive. In the hotter region, the system behaves like the (totally) asymmetric exclusion processes (TASEP), and experiences particle blockage in front of the interface to the critical region. To explain this particle density accumulation near the interface, we have measured the steady-state current in the KLS model at T>T_{c} and found it to decay as 1/T. In analogy with TASEP systems containing "slow" bonds, we argue that transport in the high-temperature subsystem is impeded by the lower current in the cooler region, which tends to set the global stationary particle current value. This blockage is induced by the extended particle clusters, growing logarithmically with system size, in the critical region. We observe the density profiles in both high- and low-temperature subsystems to be similar to the well-characterized coexistence and maximal-current phases in (T)ASEP models with open boundary conditions, which are respectively governed by hyperbolic and trigonometric tangent functions. Yet if the lower temperature is set to T_{c}, we detect marked fluctuation corrections to the mean-field density profiles, e.g., the corresponding critical KLS power-law density decay near the interfaces into the cooler region.

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.

Statistical mechanics of a single active slider on a fluctuating interface

We study the statistical mechanics of a single active slider on a fluctuating interface, by means of numerical simulations and theoretical arguments. The slider, which moves by definition towards the interface minima, is active as it also stimulates growth of the interface. Even though such a particle has no counterpart in thermodynamic systems, active sliders may provide a simple model for ATP-dependent membrane proteins that activate cytoskeletal growth. We find a wide range of dynamical regimes according to the ratio between the timescales associated with the slider motion and the interface relaxation. If the interface dynamics is slow, the slider behaves like a random walker in a random environment, which, furthermore, is able to escape environmental troughs by making them grow. This results in different dynamic exponents to the interface and the particle: the former behaves as an Edward-Wilkinson surface with dynamic exponent 2, whereas the latter has dynamic exponent 3/2. When the interface is fast, we get sustained ballistic motion with the particle surfing a membrane wave created by itself. However, if the interface relaxes immediately (i.e., it is infinitely fast), particle motion becomes symmetric and goes back to diffusive. Due to such a rich phenomenology, we propose the active slider as a toy model of fundamental interest in the field of active membranes and, generally, whenever the system constituent can alter the environment by spending energy.

Bacterial range expansions on a growing front: Roughness, fixation, and directed percolation

Directed percolation (DP) is a classic model for nonequilibrium phase transitions into a single absorbing state (fixation). It has been extensively studied by analytical and numerical techniques in diverse contexts. Recently, DP has appeared as a generic model for the evolutionary and ecological dynamics of competing bacterial populations. Range expansion-the stochastic reproduction of bacteria competing for space to be occupied by their progeny-leads to a fluctuating and rough growth front, which is known from experiment and simulation to affect the underlying critical behavior of the DP transition. In this work, we employ symmetry arguments to construct a pair of nonlinear stochastic partial differential equations describing the coevolution of surface roughness with the composition field of DP. Macroscopic manifestations (phenomenology) of these equations on growth patterns and genealogical tracks of range expansion are discussed; followed by a renormalization group analysis of possible scaling behaviors at the DP transition.

Anomalous one-dimensional fluctuations of a simple two-dimensional random walk in a large-deviation regime

The following question is the subject of our work: could a two-dimensional (2D) random path pushed by some constraints to an improbable "large-deviation regime" possess extreme statistics with one-dimensional (1D) Kardar-Parisi-Zhang (KPZ) fluctuations? The answer is positive, though nonuniversal, since the fluctuations depend on the underlying geometry. We consider in detail two examples of 2D systems for which imposed external constraints force the underlying stationary stochastic process to stay in an atypical regime with anomalous statistics. The first example deals with the fluctuations of a stretched 2D random walk above a semicircle or a triangle. In the second example we consider a 2D biased random walk along a channel with forbidden voids of circular and triangular shapes. In both cases we are interested in the dependence of a typical span 〈d(t)〉∼t^{γ} of the trajectory of t steps above the top of the semicircle or the triangle. We show that γ=1/3, i.e., 〈d(t)〉 shares the KPZ statistics for the semicircle, while γ=0 for the triangle. We propose heuristic derivations of scaling exponents γ for different geometries, justify them by explicit analytic computations, and compare with numeric simulations. For practical purposes, our results demonstrate that the geometry of voids in a channel might have a crucial impact on the width of the boundary layer and, thus, on the heat transfer in the channel.

Competing Universalities in Kardar-Parisi-Zhang Growth Models

We report on the universality of height fluctuations at the crossing point of two interacting (1+1)-dimensional Kardar-Parisi-Zhang interfaces with curved and flat initial conditions. We introduce a control parameter p as the probability for the initially flat geometry to be chosen and compute the phase diagram as a function of p. We find that the distribution of the fluctuations converges to the Gaussian orthogonal ensemble Tracy-Widom distribution for p<0.5, and to the Gaussian unitary ensemble Tracy-Widom distribution for p>0.5. For p=0.5 where the two geometries are equally weighted, the behavior is governed by an emergent Gaussian statistics in the universality class of Brownian motion. We propose a phenomenological theory to explain our findings and discuss possible applications in nonequilibrium transport and traffic flow.

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.

Width Scaling of an Interface Constrained by a Membrane

We investigate the shape of a growing interface in the presence of an impenetrable moving membrane. The two distinct geometrical arrangements of the interface and membrane, obtained by placing the membrane behind or ahead of the interface, are not symmetrically related. On the basis of numerical results and an exact calculation, we argue that these two arrangements represent two distinct universality classes for interfacial growth: while the well-established Kardar-Parisi-Zhang (KPZ) growth is obtained in the "ahead" arrangement, we find an arrested KPZ growth with a smaller roughness exponent in the "behind" arrangement. This suggests that the surface properties of growing cell membranes and expanding bacterial colonies, for example, are fundamentally distinct.

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 modes in d-dimensional directed polymers

We define a stochastic lattice model for a fluctuating directed polymer in d≥2 dimensions. This model can be alternatively interpreted as a fluctuating random path in two dimensions, or a one-dimensional asymmetric simple exclusion process with d-1 conserved species of particles. The deterministic large dynamics of the directed polymer are shown to be given by a system of coupled Kardar-Parisi-Zhang (KPZ) equations and diffusion equations. Using nonlinear fluctuating hydrodynamics and mode coupling theory we argue that stationary fluctuations in any dimension d can only be of KPZ type or diffusive. The modes are pure in the sense that there are only subleading couplings to other modes, thus excluding the occurrence of modified KPZ-fluctuations or Lévy-type fluctuations, which are common for more than one conservation law. The mode-coupling matrices are shown to satisfy the so-called trilinear condition.

Size Sorting on the Rubble-Pile Asteroid Itokawa

Photographs of the asteroid Itokawa reveal unexpectedly strong size segregation between lowlands populated almost entirely by small pebbles and highlands consisting of larger boulders. We propose that this segregation may be caused by a simple and unexplored effect: pebbles accreting onto the asteroid rebound from boulders, but sink into pebbly regions. By number, overwhelmingly more particles on Itokawa are pebbles, and collisions involving these pebbles must unavoidably cause pebbly regions to grow. We carry out experiments and simulations that demonstrate that this mechanism of size sorting based on simple counting of grains produces strong lateral segregation that reliably obeys an analytic formula.

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.

Random-field Ising model on isometric lattices: Ground states and non-Porod scattering

We use a computationally efficient graph cut method to obtain ground state morphologies of the random-field Ising model (RFIM) on (i) simple cubic (SC), (ii) body-centered cubic (BCC), and (iii) face-centered cubic (FCC) lattices. We determine the critical disorder strength Δ_{c} at zero temperature with high accuracy. For the SC lattice, our estimate (Δ_{c}=2.278±0.002) is consistent with earlier reports. For the BCC and FCC lattices, Δ_{c}=3.316±0.002 and 5.160±0.002, respectively, which are the most accurate estimates in the literature to date. The small-r behavior of the correlation function exhibits a cusp regime characterized by a cusp exponent α signifying fractal interfaces. In the paramagnetic phase, α=0.5±0.01 for all three lattices. In the ferromagnetic phase, the cusp exponent shows small variations due to the lattice structure. Consequently, the interfacial energy E_{i}(L) for an interface of size L is significantly different for the three lattices. This has important implications for nonequilibrium properties.

Surface growth on percolation networks by a conserved-noise restricted solid-on-solid growth model

Surface growth by the conserved-noise restricted solid-on-solid model is investigated on diluted lattices, i.e., on percolation networks that are embedded in two spatial dimensions. The growth exponent β and the roughness exponent α are defined, respectively, by the mean-square surface width via W(2)(t)∼t(2β) and the mean-square saturated width via W(sat)(2)(L)∼L(2α), where L is the system size. These are measured on both an infinite network and a backbone network and the results are compared with power-counting predictions obtained using the fractional Langevin equation. While the Monte Carlo results on deterministic fractal substrates show excellent agreement with the predictions [D. H. Kim and J. M. Kim, Phys. Rev. E 84, 011105 (2011)], the results on critical percolation networks deviate by 8%-12% from these predictions.

Understanding the mechanisms of solid-water reactions through analysis of surface topography

The topography of a reactive surface contains information about the reactions that form or modify the surface and, therefore, it should be possible to characterize reactivity using topography parameters such as surface area, roughness, or fractal dimension. As a test of this idea, we consider a two-dimensional (2D) lattice model for crystal dissolution and examine a suite of topography parameters to determine which may be useful for predicting rates and mechanisms of dissolution. The model is based on the assumption that the reactivity of a surface site decreases with the number of nearest neighbors. We show that the steady-state surface topography in our model system is a function of, at most, two variables: the ratio of the rate of loss of sites with two neighbors versus three neighbors (d(2)/d(3)) and the ratio of the rate of loss of sites with one neighbor versus three neighbors (d(1)/d(3)). This means that relative rates can be determined from two parameters characterizing the topography of a surface provided that the two parameters are independent of one another. It also means that absolute rates cannot be determined from measurements of surface topography alone. To identify independent sets of topography parameters, we simulated surfaces from a broad range of d(1)/d(3) and d(2)/d(3) and computed a suite of common topography parameters for each surface. Our results indicate that the fractal dimension D and the average spacing between steps, E[s], can serve to uniquely determine d(1)/d(3) and d(2)/d(3) provided that sufficiently strong correlations exist between the steps. Sufficiently strong correlations exist in our model system when D>1.5 (which corresponds to D>2.5 for real 3D reactive surfaces). When steps are uncorrelated, surface topography becomes independent of step retreat rate and D is equal to 1.5. Under these conditions, measures of surface topography are not independent and any single topography parameter contains all of the available mechanistic information about the surface. Our results also indicate that root-mean-square roughness cannot be used to reliably characterize the surface topography of fractal surfaces because it is an inherently noisy parameter for such surfaces with the scale of the noise being independent of length scale.

Generalized model for solid-on-solid interface growth

We present a probabilistic cellular automaton (PCA) model to study solid-on-solid interface growth in which the transition rules depend on the local morphology of the profile obtained from the interface representation of the PCA. We show that the model is able to reproduce a wide range of patterns whose critical roughening exponents are associated to different universality classes, including random deposition, Edwards-Wilkinson, and Kardar-Parisi-Zhang. By means of the growth exponent method, we consider a particular set of the model parameters to build the two-dimensional phase diagram corresponding to a planar cut of the higher dimensional parameter space. A strong indication of phase transition between different universality classes can be observed, evincing different regimes of deposition, from layer-by-layer to Volmer-Weber and Stransk-Krastanov-like modes. We expect that this model can be useful to predict the morphological properties of interfaces obtained at different surface deposition problems, since it allows us to simulate several experimental situations by setting the values of the specific transition probabilities in a very simple and direct way.

Skewness in (1+1)-dimensional Kardar-Parisi-Zhang-type growth

We use the (1+1)-dimensional Kardar-Parisi-Zhang equation driven by a Gaussian white noise and employ the dynamic renormalization-group of Yakhot and Orszag without rescaling [J. Sci. Comput. 1, 3 (1986)]. Hence we calculate the second- and third-order moments of height distribution using the diagrammatic method in the large-scale and long-time limits. The moments so calculated lead to the value S=0.3237 for the skewness. This value is comparable with numerical and experimental estimates.

Aging of the (2+1)-dimensional Kardar-Parisi-Zhang model

Extended dynamical simulations have been performed on a (2+1)-dimensional driven dimer lattice-gas model to estimate aging properties. The autocorrelation and the autoresponse functions are determined and the corresponding scaling exponents are tabulated. Since this model can be mapped onto the (2+1)-dimensional Kardar-Parisi-Zhang surface growth model, our results contribute to the understanding of the universality class of that basic system.

(2+1)-Dimensional directed polymer in a random medium: scaling phenomena and universal distributions

We examine numerically the zero-temperature (2+1)-dimensional directed polymer in a random medium, along with several of its brethren via the Kardar-Parisi-Zhang (KPZ) equation. Using finite-size and KPZ scaling Ansätze, we extract the universal distributions controlling fluctuation phenomena in this canonical model of nonequilibrium statistical mechanics. Specifically, we study point-point, point-line, and point-plane directed polymer geometries, scenarios which yield higher-dimensional analogs of the Tracy-Widom distributions of random matrix theory. Our analysis represents a robust, multifaceted numerical characterization of the 2+1 KPZ universality class and its limit distributions.

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.

High-order harmonic generation employing field enhancement by metallic fractal rough surfaces

We numerically investigate high-order harmonic generation (HHG) in noble gases in the vicinity of fractal structures of metallic rough surfaces described by the restricted solid-on-solid model. The calculated intensity enhancement factors in the range of 10³ enables HHG up to the 50-th order with low pump intensity down to tens of GW/cm². The increased interaction volume of "hot spots" in the case of grazing incidence of s-polarized pump pulses leads to an efficiency of harmonics in the plateau region of about 10⁻⁷.

Smoothly varying hopping rates in driven flow with exclusion

We consider the one-dimensional totally asymmetric simple exclusion process (TASEP) with position-dependent hopping rates. The problem is solved, in a mean-field adiabatic approximation, for a general (smooth) form of spatial rate variation. Numerical simulations of systems with hopping rates varying linearly against position (constant rate gradient), for both periodic and open-boundary conditions, provide detailed confirmation of theoretical predictions, concerning steady-state average density profiles and currents, as well as open-system phase boundaries, to excellent numerical accuracy.

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.

Slow relaxation and aging kinetics for the driven lattice gas

We numerically investigate the long-time behavior of the density-density autocorrelation function in driven lattice gases with particle exclusion and periodic boundary conditions in one, two, and three dimensions using precise Monte Carlo simulations. In the one-dimensional asymmetric exclusion process on a ring with half the lattice sites occupied, we find that correlations induce extremely slow relaxation to the asymptotic power law decay. We compare the crossover functions obtained from our simulations with various analytic results in the literature and analyze the characteristic oscillations that occur in finite systems away from half-filling. As expected, in three dimensions correlations are weak and consequently the mean-field description is adequate. We also investigate the relaxation toward the nonequilibrium steady state in the two-time density-density autocorrelations, starting from strongly correlated initial conditions. We obtain simple aging scaling behavior in one, two, and three dimensions, with the expected power laws.

Density profiles of loose and collapsed cohesive granular structures generated by ballistic deposition

Loose granular structures stabilized against gravity by an effective cohesive force are investigated on a microscopic basis using contact dynamics. We study the influence of the granular Bond number on the density profiles and the generation process of packings, generated by ballistic deposition under gravity. The internal compaction occurs discontinuously in small avalanches and we study their size distribution. We also develop a model explaining the final density profiles based on insight about the collapse of a packing under changes of the Bond number.

Growth and Structure of Pyrogenic Silica

Fractal analysis is used to interpret small-angle X-ray and neutron scattering data on commercial fumed silica powders. Rough primary particles are found in high surface area powders. A model based on ballistic polymerization, sintering, and diffusion-limited aggregation is proposed to explain the observed powder structure.

Ballistic deposition patterns beneath a growing Kardar-Parisi-Zhang interface

We consider a (1+1)-dimensional ballistic deposition process with next-nearest-neighbor interactions, which belongs to the Kardar-Parisi-Zhang (KPZ) universality class. The focus of our analysis is on the properties of structures appearing in the bulk of a growing aggregate: a forest of independent clusters separated by "crevices." Competition for growth (mutual screening) between different clusters results in "thinning" of this forest, i.e., the number density c(h) of clusters decreases with the height h of the pattern. For the discrete stochastic equation describing the process we introduce a variational formulation similar to that used for the randomly forced continuous Burgers equation. This allows us to identify the "clusters" and crevices with minimizers and shocks in the Burgers turbulence. Capitalizing on the ideas developed for the latter process, we find that c(h) ∼ h(-α) with α=2/3. We compute also scaling laws that characterize the ballistic deposition patterns in the bulk: the law of transversal fluctuations of cluster boundaries and the size distribution of clusters. It turns out that the intercluster interface is superdiffusive: the corresponding exponent is twice as large as the KPZ exponent for the surface of the aggregate. Finally we introduce a probabilistic concept of ballistic growth, dubbed the "hairy" Airy process in view of its distinctive geometric features. Its statistical properties are analyzed numerically.

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.

Surface pattern formation and scaling described by conserved lattice gases

We extend our 2+1 -dimensional discrete growth model [Odor, Phys. Rev. E 79, 021125 (2009)] with conserved, local exchange dynamics of octahedra, describing surface diffusion. A roughening process was realized by uphill diffusion and curvature dependence. By mapping the slopes onto particles, two-dimensional nonequilibrium binary lattice model emerges, in which the (smoothing or roughening) surface diffusion can be described by attracting or repelling motion of oriented dimers. The binary representation allows simulations on very large size and time scales. We provide numerical evidence for Mullins-Herring or molecular-beam epitaxy class scaling of the surface width. The competition of inverse Mullins-Herring diffusion with a smoothing deposition, which corresponds to a Kardar-Parisi-Zhang (KPZ) process, generates different patterns: dots or ripples. We analyze numerically the scaling and wavelength growth behavior in these models. In particular, we confirm by large size simulations that the KPZ type of scaling is stable against the addition of this surface diffusion, hence this is the asymptotic behavior of the Kuramoto-Sivashinsky equation as conjectured by field theory in two dimensions, but has been debated numerically. If very strong, normal surface diffusion is added to a KPZ process, we observe smooth surfaces with logarithmic growth, which can describe the mean-field behavior of the strong-coupling KPZ class. We show that ripple coarsening occurs if parallel surface currents are present, otherwise logarithmic behavior emerges.

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.

Surface growth on diluted lattices by a restricted solid-on-solid model

An influence of diluted sites on surface growth has been investigated, using the restricted solid-on-solid model. It was found that, with respect to equilibrium growth, the surface width and the saturated width exhibited universal power-law behaviors, i.e., W approximately t(beta) and W(sat) approximately L(zeta), regarding all cases with respect to the concentration of diluted sites x=1-p , with p being the occupation probability on each lattice site. For x < x(c) (=1-p(c), p(c) being the percolation threshold), the growth appeared to be similar to that of a regular lattice, both in two and three dimensions. For x=x(c), the growth yielded nontrivial exponents which were different from those on a regular lattice. In nonequilibrium growth, a considerable amount of diluted sites (x < or = x(c)) appeared to yield nonuniversal growth, unlike the case of a regular lattice. The cause of nonuniversal growth dynamics has been investigated, considering the growth on a backbone cluster and on lattices constructed with periodically and randomly diluted subcells.