Next Article in Journal
Editorial to the Special Issue “Propagation of Coronal Mass Ejections”
Previous Article in Journal
Quiescent and Active Galactic Nuclei as Factories of Merging Compact Objects in the Era of Gravitational Wave Astronomy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Random Walk on a Rough Surface: Renormalization Group Analysis of a Simple Model

by
Nikolay V. Antonov
1,2,*,
Nikolay M. Gulitskiy
1,
Polina I. Kakin
1,* and
Dmitriy A. Kerbitskiy
1
1
Department of Physics, Saint Petersburg State University, Universitetskaya nab. 7/9, St. Petersburg 199034, Russia
2
N.N. Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, Dubna 141980, Moscow Region, Russia
*
Authors to whom correspondence should be addressed.
Universe 2023, 9(3), 139; https://doi.org/10.3390/universe9030139
Submission received: 26 January 2023 / Revised: 27 February 2023 / Accepted: 2 March 2023 / Published: 7 March 2023
(This article belongs to the Section Field Theory)

Abstract

:
The field-theoretic renormalization group is applied to a simple model of a random walk on a rough fluctuating surface. We consider the Fokker–Planck equation for a particle in a uniform gravitational field. The surface is modeled by the generalized Edwards–Wilkinson linear stochastic equation for the height field. The full stochastic model is reformulated as a multiplicatively renormalizable field theory, which allows for the application of the standard renormalization theory. The renormalization group equations have several fixed points that correspond to possible scaling regimes in the infrared range (long times and large distances); all the critical dimensions are found exactly. As an example, the spreading law for the particle’s cloud is derived. It has the form R 2 ( t ) t 2 / Δ ω with the exactly known critical dimension of frequency Δ ω and, in general, differs from the standard expression R 2 ( t ) t for an ordinary random walk.

1. Introduction

Over decades, stochastic growth processes, kinetic roughening phenomena, and fluctuating surfaces or interfaces have been attracting constant attention. The most prominent examples include the deposition of a substance on a surface and the growth of the corresponding phase boundary; propagation of flame, smoke, and solidification fronts; growth of vicinal surfaces and bacterial colonies; erosion of landscapes and seabed profiles; molecular beam epitaxy; and many others, see [1,2,3,4,5,6,7,8,9,10,11,12,13] and references therein.
Another vast area of research is that of diffusion and random walks in a random environment, such as disordered, inhomogeneous, porous, or turbulent media, see, e.g., [14,15,16,17,18,19,20,21,22,23,24].
In this paper, we study a simple model of a random walk on a rough fluctuating surface. We consider the Fokker–Planck equation for a particle in a uniform gravitational field. The surface is modeled by the generalized Edwards–Wilkinson linear stochastic equation for the height field [1]. The generalized model involves two arbitrary exponents, ε and η , related to the spectrum and the dispersion law of the height field, respectively. A detailed description of the model and its relation to various special cases is given in Section 2.
Using the general Martin–Siggia–Rose–de Dominicis–Janssen theorem, the original stochastic problem is reformulated as a certain field-theoretic model. This allows one to apply the well-developed formalism of Feynman diagrammatic techniques, renormalization theory, and a renormalization group (RG). The model is shown to be multiplicatively renormalizable, so that the RG equation can be derived in a standard way. The corresponding renormalization constants and the RG functions (anomalous dimensions and β functions) are explicitly calculated in the leading one-loop order of the RG perturbation theory. These issues are discussed in Section 3 and Section 4.
The RG equations have two Gaussian (free) fixed points and two nontrivial ones. Those points are infrared (IR) attractive depending on the values of the parameters ε and η , which implies the existence of scaling (self-similar) asymptotic regimes in the IR range (long times and large distances) for the various response and correlation functions of the model (Section 4). The critical dimensions for those regimes are found exactly as functions of ε and η . As an indicative application, the time dependence of the mean-square radius of a cloud of randomly walking particles is obtained (Section 5). It is described by a power law with the exponent that depends on the fixed point, is known exactly as a function of ε and η , and for nontrivial points, it differs from the ordinary random walk: R 2 ( t ) t .
Some implications and possible generalizations are discussed in Section 6.

2. Description of the Model

We consider the random walk of a point particle on a two-dimensional rough surface embedded into the ( d + 1 ) -dimensional space. The particle is located on the surface on the height h ( t , x ) , where x ( t ) = { x i ( t ) } is the particle’s coordinate projection on the d-dimensional substrate, i = 1 d . Thus, d is an arbitrary (for generality) dimension of the substrate x space.
While the coordinates x i with i = 1 , , d determine the particle’s location, the coordinate x i with i = ( d + 1 ) is not treated as an independent Cartesian coordinate but is restricted to the surface, h ( t , x ) x ( d + 1 ) . This setup excludes the possibility for a particle to “jump off” or “escape” the surface; such interesting phenomena are not included in our simple model.
The basic stochastic equation of motion of a particle located at the point x ( t ) = { x i ( t ) } in an external drift field F has the form [14,15,16,17,18,19,20,21,22,23,24]:
t x i = F i ( t , x ) + ζ i , ζ i ( t ) ζ j ( t ) ζ = 2 ν 0 δ ( t t ) .
Here, ζ i = ζ i ( t ) is a Gaussian noise with zero mean and a given pair correlation function, ν 0 > 0 will play a role of the diffusion coefficient, and F is an external drift field (a force or an advecting velocity, depending on the specific context).1
The probability distribution function P ( t , x ) satisfies the Fokker–Planck equation
t + i ( F i ν 0 i ) P ( t , x ) = 0 .
Here and below, summation over repeated indices is implied.
From physics reasons, the drift F in a gravitational field should obey the O ( d ) symmetry and the invariance with the respect to the shift h h + const, so that it should be built of gradients of the field h. Thus, in the simplest linear approximation, it is taken in the form
F i = λ 0 i h ,
with the parameter λ 0 > 0 proportional to the particle’s mass m and the gravitational acceleration g. Possible higher-order corrections to the linear approximation should also be constructed from gradients of the field h and obey the O ( d ) symmetry, for example, i ( l h l h ) n . They have higher canonical dimensions in comparison with (3) (see Section 3), are IR irrelevant (in the sense of Wilson), and should be dropped in the analysis of IR scaling.
Careful interpretation of the gradient i h for a rough surface and the very existence (in a rigorous mathematical sense) of corresponding continuous equations is a serious problem. Recently, important progress was achieved for the Kardar–Parisi–Zhang (KPZ) model, see [25,26,27,28,29] and references therein.
From a more practical physical point of view, there is a very small microscopical scale a below which the field h becomes smooth and differentiable. In our and the KPZ’s treatment, this scale is tacitly set to zero, so that the field becomes rough. There is an apparent analogy with the well-known dissipative anomaly in turbulence, see, e.g., [30]. Practically, in this study, we use the formal perturbation theory, where this problem does not arise, and the expression (3) is applied without further comments.
The simplest model of surface roughening, proposed within the context of landscape erosion, is the one due to Edwards and Wilkinson [1]. In the continuous formulation, it is described by the diffusion-type stochastic equation for the height field h = h ( t , x ) :
{ t κ 0 2 } h ( t , x ) = f ( t , x ) ,
where κ 0 > 0 is (a kind of) surface tension coefficient, 2 = i i is the Laplace operator, and f is a Gaussian random noise with zero mean and a given pair correlation function. The most popular choices are the white noise
f ( t , x ) f ( t , x ) f = D 0 δ ( t t ) δ ( x x )
with the positive amplitude D 0 > 0 and the quenched noise; the simplified version of the latter is
f ( t , x ) f ( t , x ) f = D 0 δ ( x x ) .
In this paper, we consider a generalized equation
{ t + κ 0 k 2 η } h ( t , x ) = f ( t , x ) ,
written here in the symbolic notation with k being the wave number,2 while the correlation function is taken in a power-like form:
f ( t , x ) f ( t , x ) f = D 0 δ ( t t ) d k ( 2 π ) d k 2 d y exp { i k ( x x ) } .
Here, η and y are arbitrary exponents and d is the dimension of space. Clearly, the choice η = 0 , 2 d y = 0 corresponds to the model (4) and (5); as we will see, the model (4) and (6) can also be obtained from (7) and (8).
For a linear stochastic equation with a Gaussian additive random noise, the field h is also a Gaussian field defined by its pair correlation function. For the model (7) and (8), the latter has the following form in the Fourier ( ω k ) representation
D h ( ω , k ) = D 0 k 2 d y ω 2 + [ κ 0 k 2 η ] 2 = g 0 u 0 ν 0 3 k 2 d η ε ω 2 + [ u 0 ν 0 k 2 η ] 2 .
In the second relation, we introduced the new variables: the exponent ε and the amplitudes g 0 , u 0 , defined by the relations
ε = y η , κ 0 = u 0 ν 0 , D 0 = g 0 u 0 ν 0 3 .
They are convenient, in particular, because the equal-time correlation function
D h ( k ) = d ω 2 π D ( ω , k ) g 0 ν 0 2 k d ε
involves the parameters g 0 , ε , while the dispersion law
ω ( k ) u 0 ν 0 k 2 η
is expressed only via u 0 , η .
The choice η 0 can be justified by the ideas of self-organized criticality (SOC), according to which the evolution of a sandpile surface is not an ordinary diffusion-type process but involves several discrete steps: expectation period, reaching a threshold, and avalanche, see, e.g., [31,32,33].
Thus, according to [32], self-organized critical dynamical systems give rise to the so-called 1 / f α noise because the characteristic size of an avalanche is related to its lifetime via a power law s t 1 + γ , where the exponent γ is the rate at which the event propagates across the system, see also, e.g., Section 1.3.2 in the book [31] and papers [33,34]. In the ω k representation, this corresponds to the dispersion law (12) with the exponent η = ( 1 + 2 γ ) / ( 1 + γ ) . It is also worth noting that the 1 / f α noise appears also in models of random walks in random media, see, e.g., [14,15].
The model (9) includes two special cases interesting on their own. In the limit u 0 and g 0 = g 0 / u 0 fixed, the function D ( ω , k ) becomes independent of the frequency ω , and the field h ( t , x ) becomes white in time. Indeed, one obtains in the (t k ) representation
D ( t t , k ) = δ ( t t ) g 0 ν 0 2 k 2 d ε + η .
Here, the exponent 0 < ( ε η ) < 2 plays a role of a Hölder’s exponent that measures “roughness” of the field h (“Batchelor limit” ( ε η ) 2 corresponds to a smooth field).
In the limit u 0 0 and g 0 fixed, the function D h ( k ) in (11) remains finite, so that (9) tends to
D h ( ω , k ) = π δ ( ω ) g 0 ν 0 2 k d ε ,
which corresponds to the time-independent (quenched or frozen) field h. Surprisingly enough, for ε = 4 d , this reproduces the model (4) and (6) where one has D h δ ( ω ) / k 4 .
Substituting the gravitational force (3) with the random height field from (7) and (8) into the Fokker–Planck equation (2) turns the latter into a stochastic equation in its own right. It has the form
t θ = ν 0 2 θ + λ 0 i ( θ i h ) + f ,
where the random field θ ( t , x ) can be interpreted as the density of walking particles, while the role of the (deterministic) probability distribution function P ( t , x ) is now conveyed to the linear response function, see Equation (44) in Section 5.
This completes formulation of the problem.

3. Field-Theoretic Formulation and Renormalization of the Model

According to the general theorem (see, e.g., Section 5.3 in the monograph [35]), the full stochastic problem (8), (15) is equivalent to the field-theoretic model for the doubled set of fields Φ = { θ , h , θ , h } with the de Dominicis–Janssen action functional:
S ( Φ ) = θ t θ + ν 0 2 θ + λ 0 i ( θ i h ) + S h ( h , h ) ,
S h ( h , h ) = 1 2 h D f h + h t + κ 0 k 2 η h .
Here, D f is the correlator (8), θ is the density field, h is the height field, and θ , h are the corresponding Martin–Siggia–Rose response fields; all the needed integrations over their arguments x = { t , x } and summations over repeated indices are implied. The field-theoretic formulation means that various correlation and response functions of the original stochastic problem are represented by functional averages with the weight exp S ( Φ ) . The field h can easily be removed by Gaussian integration, then S h ( h , h ) would be replaced with S h ( h ) = h D h 1 h / 2 with D h from (9), but the expanded representation (17) is more convenient for the renormalization purposes. The constant λ 0 can be removed by rescaling of the fields h , h and other parameters. Thus, in the following, with no loss of generality, we set λ 0 = 1 .
The model (16) and (17) corresponds to Feynman diagrammatic technique with bare propagators θ θ 0 , h h 0 , h h 0 (the latter does not enter into relevant diagrams) and the only vertex θ i ( θ i h ) .
It is well known that an analysis of ultraviolet (UV) divergences is based on an analysis of canonical dimensions, see, e.g., [35] (Sections 1.15 and 1.16). In contrast to conventional static models, dynamic ones have two independent scales: a time scale [ T ] and a spatial scale [ L ] , see [35] (Sections 1.17 and 5.14). Thus, the canonical dimension of any quantity F (a field or a parameter) is determined by two numbers: the frequency dimension d F ω and the momentum dimension d F k :
F T d F ω L d F k .
The dimensions are found from obvious normalization conditions
d k k = d x k = 1 , d k ω = d x ω = 0 , d ω k = d t k = 0 , d ω ω = d t ω = 1
and from the requirement that all terms in the action functional be dimensionless with respect to both the canonical dimensions separately. The total canonical dimension is defined as d F = d F k + 2 d F ω (the coefficient 2 follows from the relation t 2 in the free theory). In the renormalization procedure, d F plays the same role as the conventional (momentum) dimension does in static models, see Section 5.14 in [35].
Canonical dimensions of all the fields and parameters of our model are given in Table 1. It also involves renormalized parameters (without subscript “0”) and the reference mass μ , an additional parameter of the renormalized theory; they all will appear later on.
Note that for the fields θ , θ all these dimensions can be unambiguously defined only for the product θ θ . Formally, this follows from the invariance of the action functional (16) under the dilatation θ λ θ , θ λ 1 θ .
As can be seen from Table 1, the model becomes logarithmic (both coupling constants g 0 , u 0 become dimensionless) for η = y = 0 (or equivalently for ε = y = 0 ) and arbitrary d.3 According to the general strategy of renormalization, the exponents η , y, or ε that “measure” the deviation from the logarithmicity should be treated as formal small parameters of the same order. The UV divergences manifest themselves as singularities at y 0 , etc., in the correlation functions; in the one-loop approximation, they have the form of simple poles.
The total canonical dimension of a certain 1-irreducible Green’s function is given by
d Γ = ( d + 2 ) Φ d Φ N Φ ,
where N Φ are the numbers of the fields Φ = { θ , h , θ , h } entering the Green’s function and d Φ are their total canonical dimensions.
The formal index of divergence δ Γ is the total dimension of the Green’s function in the logarithmic theory ( y = η = 0 ), that is, δ Γ = d Γ | y = η = 0 . Superficial UV divergences, whose removal requires introducing counterterms, can be present in the Green’s function Γ if δ Γ is a non-negative integer.
When analyzing the divergences in the model (16) and (17), the following additional considerations should be taken into account, see, e.g., [35] (Section 5.15) and [36] (Section 1.4).
(i) For any dynamic model of this type, all the 1-irreducible functions without the response fields contain closed circuits of retarded propagators θ θ 0 and vanish. Thus, it is sufficient to consider the functions with N θ + N h 1 .
(ii) For all non-vanishing functions, N θ = N θ (otherwise no diagrams can be constructed). Formally, this is a consequence of the invariance of the action functional (16) with respect to dilatation θ λ θ , θ λ 1 θ .
(iii) Using integration by parts, one derivative in the vertex can be moved onto the field θ , i.e., θ i ( θ i h ) ( i θ ) ( i h ) θ . Thus, in any 1-irreducible diagram, each external field θ or h “releases” the external momentum, and the real index of divergence decreases by the corresponding number of units, i.e., δ = δ N θ N h . Furthermore, these fields enter the counterterms only in the form of spatial gradients. This observation excludes the counterterms θ t θ and ( θ θ ) 2 , the latter allowed by the formal index for d 2 .
(iv) It is clear that the fields θ , θ do not affect the statistics of the field h. In the field-theoretic terms, this “passivity” means that any 1-irreducible Green’s function with N θ = 0 , N θ > 0 , and N h + N h > 0 vanishes: no corresponding diagrams can be constructed.
Taking into account these considerations, one obtains:
δ = ( d + 2 ) d ( N θ + N h ) , δ = ( d + 2 ) ( d + 1 ) N θ N h d N h
(we recall that N θ = N θ , so that only N θ is indicated).
Then, the straightforward analysis shows that the superficial divergences in our model are present only in the 1-irreducible functions θ θ and θ θ h , and the corresponding counterterms necessarily contract to the forms θ 2 θ ( δ = 2 , δ = 1 ) and ( i θ ) ( i h ) θ ( δ = 2 , δ = 0 ). Such terms are already present in the action (16), which means that our model (16) and (17) is multiplicatively renormalizable with only two independent renormalization constants Z 1 and Z 2 .
The renormalized action has the form
S R ( Φ ) = θ t θ + Z 1 ν 2 θ + Z 2 i ( θ i h ) + S h R ( h , h ) ,
which is naturally reproduced as renormalization of the field h and the coefficient ν 0 ; no renormalization of the product θ θ is needed:
ν 0 = ν Z ν , Z ν = Z 1 , Z h = Z 2 , Z θ θ = 1 .
The functional (17) is not renormalized, S h R ( h , h ) = S h ( h , h ) , but it should be expressed in renormalized variables, taking into account Equations (8) and (9):
g 0 = g μ y Z g , u 0 = u μ η Z u , κ 0 = κ Z κ ,
where the renormalization mass μ is introduced so that renormalized couplings g and u are completely dimensionless. Then, it follows from the absence of renormalization of S h that
Z h Z h = 1 , Z h 2 Z g Z u Z ν 3 = 1 , Z u Z ν = Z κ = 1 .
Along with (21), this finally gives the following relations:
Z g = Z 2 2 Z 1 1 , Z u = Z 1 1 , Z ν = Z 1 .
We calculated the renormalization constants Z 1 and Z 2 in the leading one-loop approximation (the first order of the perturbative expansion in g). It is sufficient to find them for η = 0 , because the anomalous dimensions in the minimal subtraction (MS) renormalization scheme are independent of the parameters such as η and y, while the exponent y alone provides UV regularization. Then, one obtains:
Z 1 = 1 g y C d 2 d ( u 1 ) ( u + 1 ) 2 , Z 2 = 1 + g y C d 2 d 1 ( u + 1 ) 2 ,
with the higher-order corrections in g. Here, C d = S d / ( 2 π ) d , S d = 2 π d / 2 / Γ ( d / 2 ) is the surface area of the unit sphere in d-dimensional space. It is convenient to absorb overall factors into the coupling constant g, which gives
Z 1 = 1 g y ( u 1 ) ( u + 1 ) 2 , Z 2 = 1 + g y 1 ( u + 1 ) 2 .
For η 0 , the expressions (25) and (26) would be infinite sums, see, e.g., [37].

4. RG Equations, RG Functions, and Fixed Points

Because our model is multiplicatively renormalizable, the corresponding RG equations are derived in a standard fashion. In particular, for a certain renormalized (full or connected) Green’s function W R , the RG equation reads
D μ + β g g + β u u γ ν D ν Φ N Φ γ Φ W R ( g , u , ν , μ ; ) = 0 .
Here, the ellipsis stands for other variables (times and coordinates or frequencies and momenta), x = / x , D x = x x for any variable x, and the sum runs over all fields Φ = { θ , h , θ , h } .
The coefficients in the RG differential operator (27)—the anomalous dimensions γ and the β functions—are defined as
γ α = D ˜ μ ln Z α for any α , β g = D ˜ μ g , β u = D ˜ μ u ,
where D ˜ μ is the differential operation D μ at fixed bare (unrenormalized) parameters, see, e.g., Sections 1.24 and 1.25 in the monograph [35].
From (21)–(24) and definition (28), it follows that
γ θ θ = 0 , γ h = γ h = γ 2 , γ g = 2 γ 2 γ 1 , γ u = γ ν = γ 1 ,
β g = g [ ε γ g ] , β u = u [ η γ u ] .
From (30) and the one-loop result (26), one obtains
γ 1 = g u 1 ( u + 1 ) 2 , γ 2 = g 1 ( u + 1 ) 2 ,
β g = g ε + g 2 u ( u + 1 ) 2 , β u = u η + g u 1 ( u + 1 ) 2 ,
with the higher-order corrections in g.
The IR asymptotic behavior of the Green’s functions is determined by IR attractive fixed points of the corresponding RG equations. The coordinates of fixed points g * , u * are found from the requirement that all the β functions vanish simultaneously:
β g ( g * , u * ) = β u ( g * , u * ) = 0 .
The type of a fixed point is determined by the matrix of derivatives Ω i j = i β j ( g * ) at the given point g i = { g , u } : for an IR attractive point, all the eigenvalues should have positive real parts.
An analysis of the expressions (32) reveals four fixed points:
(i) Gaussian (free) fixed point:
g * = 0 , u * = 0 ;
(ii) Nontrivial fixed point:
g * = 2 ( ε η ) 2 ε 2 η , u * = ε ε 2 η .
The point (i) is IR attractive for ε < 0 , η < 0 , while the point (ii) is IR attractive for ε > 0 , η < ε / 2 .
Two more points are found in the following way. In order to explore the limiting case u with g / u fixed, we have to pass to new variables: g g / u and w 1 / u . For this case, we obtain
β g = g η ε + g w + 1 , β w = w η + g w 1 ( w + 1 ) 2 .
Finding the zeros of the β functions, we find two additional fixed points:
(iii) Gaussian (free) fixed point:
g * = 0 , w * = 0 ;
(iv) Nontrivial fixed point:
g * = ε η , w * = 0 .
The point (iii) is IR attractive if ε > 0 , ε / 2 < η < ε , and the point (iv) is IR attractive if ε < 0 , η > 0 or ε > 0 , η > ε .
The general stability pattern of the fixed points in the ε η plane is shown in Figure 1.
In the one-loop approximation, the regions of IR stability for all the points are given by sectors that cover the full plane without gaps or overlaps between them.
Some remarks are in order. Clearly, the Gaussian points correspond to cases in which the dynamics of the field θ are not affected by the statistics of the height field h (only in the leading order of the IR asymptotic behavior!). In these cases, we deal with an ordinary random walk.
The point (iv) corresponds to the limiting case (13) when the field h, in comparison with θ , behaves as if it was δ correlated in time.
However, we did not find a nontrivial point that would correspond to the frozen limit (14). This follows from the fact that the function β g in (32) becomes trivial for u 0 : β g = ε g . A similar triviality was observed earlier in models of diffusion in time-independent potential vector fields where it was shown to be exact in all orders of perturbation theory [18,19]. Because those models have a close formal resemblance with the limit (14) of our model and its special case (4) and (6), we believe that in the latter cases β g is also trivial exactly.

5. Critical Dimensions and Scaling Behavior

The existence of IR attractive fixed points of the RG equations implies the existence of the scaling behavior of the correlation functions in the IR range.
In dynamical models, the critical dimension of any quantity F (a field or a parameter) is given by the expression (see, e.g., Sections 5.16 and 6.7 in [35] and Section 2.1 in [36])
Δ F = d F k + Δ ω d F ω + γ F * , Δ ω = 2 γ ν *
(with the standard normalization convention that Δ k = Δ x = 1 ). Here and below, γ * denotes the value of the anomalous dimension γ at a fixed point.
For the Gaussian points (i) and (iii), one has
Δ θ θ = d , Δ ω = 2 .
For the fixed point (ii), one obtains the exact results from the relation (29) and definition (30):
Δ θ θ = d , Δ ω = 2 η .
As already mentioned, the point (iv) corresponds to the limit (13), where the propagator h h 0 becomes δ -correlated in time. As a result, closed circuits of retarded propagators θ θ 0 appear in almost all diagrams relevant for a renormalization procedure and they therefore vanish. The only exception is the one-loop diagram contributing to Z 1 . Thus, one has Z 2 = 1 identically, while Z 1 is given exactly by the one-loop expression, cf. the discussion of Kraichnan’s rapid-change model of passive scalar advection [38]. Then, one readily derives the exact expressions for the critical dimensions:
Δ θ θ = d , Δ ω = 2 ε + η .
As an illustrative application, consider the mean-square distance of a random walker on a rough surface. For such a particle that started moving at t = 0 from the origin x = 0 , it is given by
R 2 ( t ) = d x x 2 θ ( t , x ) θ ( 0 , 0 ) ,
where t > 0 is a later time and x is the corresponding current position. Substituting the scaling representation for the linear response function
P ( t , x ) = θ ( t , x ) θ ( 0 , 0 ) r Δ θ θ F ( t r Δ ω )
gives
R 2 ( t ) t ( d + 2 Δ θ θ ) / Δ ω .
Taking into account the exact relation Δ θ θ = d , valid for all fixed points (i)-(iv), one arrives at the spreading law
R 2 ( t ) t 2 / Δ ω ,
with the exact expressions Δ ω = 2 for the points (i), (iii), Δ ω = 2 η for (ii), and Δ ω = 2 ε + η for (iv).

6. Conclusions

We studied a model of a random walk of a particle on a rough fluctuating surface described by the Fokker–Planck equation for a particle in a constant gravitational field while the surface was modeled by the (generalized) Edwards–Wilkinson model. The full stochastic problem, (2), (3), (7) and (8), is mapped onto a multiplicatively renormalizable field-theoretic model (16) and (17).
The corresponding RG equations reveal two Gaussian (free) and two nontrivial fixed points, which means that the system exhibits various types of IR scaling behavior (long times and large distances). Although the practical calculation is confined within the leading one-loop approximation, the main critical dimensions are found exactly.
As an illustrative example, we considered the mean-square displacement of a walking particle (in another interpretation, the radius of particles’ cloud). It shows that the particle is not trapped in a finite area but travels all across the system with a spreading law similar to the ordinary random walk but, in general, with different exponents, see (46) and the text below.
As one can see, even a comparatively simple model demonstrates interesting types of IR behavior. Thus, it is interesting to study more involved situations. There are several directions for possible generalizations.
Linear stochastic equations such as (4) and (7) (corresponding to Gaussian statistics for the height field) can be replaced by nonlinear models, such as the KPZ [2] or Pavlik’s [5,8] ones.
Although our expressions (41) and (42) for the critical dimensions are exact, they are derived within perturbation theory based on the assumption that the expansion parameters ε and η are small. Then, it is supposed that the one-loop pattern of fixed points is qualitatively correct. However, in some cases, a crossover in the scaling behavior occurs for finite values of parameters analogous to ε and η [37,39]. In the field-theoretic approach, this effect can be related to the appearance of composite operators with negative dimensions [37]. This issue requires a special investigation.
On some occasions, the motion of a particle is not an ordinary random walk (1) but is described, e.g., by Lévy flights, see, e.g., [21]. This possibility is supported by the ideas of self-organized criticality that the underlying surface evolves via avalanches [31,32,33,34], while the particle can slide upon the surface. If so, it is natural to replace the Laplace operator in the Fokker–Planck equation (2) with a fractional derivative, 2 k 2 k 2 η , with a certain new exponent η .
It is especially interesting to include anisotropy (as a consequence of an overall tilt of the surface). This can be done by describing the field h by the Pastor-Satorras–Rothman model for an eroding landscape [9,10] or the Hwa–Kardar model of a running sandpile [40,41].
This work remains for the future and is partly in progress.

Author Contributions

All authors contributed equally to this paper. All authors have read and agreed to the published version of the manuscript.

Funding

The work of P.I.K. was supported by the Foundation for the Advancement of Theoretical Physics and Mathematics “BASIS”, project 22-1-3-33-1, and by the Ministry of Science and Higher Education of the Russian Federation, agreement no. 075-15-2022-287.

Acknowledgments

The authors are indebted to M.A. Reiter for the discussion.

Conflicts of Interest

The authors declare no conflict of interest.

Notes

1
Here and below, the subscript 0 refers to bare parameters which will be renormalized in the following.
2
Detailed discussion of fractional derivatives can be found in [21].
3
Although u 0 is not an expansion parameter in perturbation theory, its renormalized counterpart is dimensionless, enters into renormalization constants and RG functions, and should be treated on equal footing with g 0 . We also recall that λ 0 = 1 .

References

  1. Edwards, S.F.; Wilkinson, D.R. The Surface Statistics of a Granular Aggregate. Proc. R. Soc. 1982, 381, 17–31. [Google Scholar]
  2. Kardar, M.; Parisi, G.; Zhang, Y.-C. Dynamic scaling of growing interfaces. Phys. Rev. Lett. 1986, 56, 889–892. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yan, H.; Kessler, D.A.; Sander, L.M. Roughening phase transition in surface growth. Phys. Rev. Lett. 1990, 64, 926–929. [Google Scholar] [CrossRef] [PubMed]
  4. Yan, H.; Kessler, D.A.; Sander, L.M. Kinetic Roughening in Surface Growth. MRS Online Proc. Libr. 1992, 278, 237–247. [Google Scholar] [CrossRef]
  5. Pavlik, S.I. Scaling for a growing phase boundary with nonlinear diffusion. JETP 1994, 79, 303, [Translated from the Russian: ZhETF 1994, 106, 553.]. [Google Scholar]
  6. Halpin-Healy, T.; Zhang, Y.-C. Kinetic roughening phenomena, stochastic growth, directed polymers and all that. Aspects of multidisciplinary statistical mechanics. Phys. Rep. 1995, 254, 215–414. [Google Scholar] [CrossRef]
  7. Barabási, A.-L.; Stanley, H.E. Fractal Concepts in Surface Growth; Cambridge University Press: Cambridge, MA, USA, 1995. [Google Scholar]
  8. Antonov, N.V.; Vasil’ev, A.N. The quantum-field renormalization group in the problem of a growing phase boundary. JETP 1995, 81, 485, [Translated from the Russian: ZhETF 1995, 108, 885.]. [Google Scholar]
  9. Pastor-Satorras, R.; Rothman, D.H. Stochastic equation for the erosion of inclined topography. Phys. Rev. Lett. 1998, 80, 4349–4352. [Google Scholar] [CrossRef] [Green Version]
  10. Pastor-Satorras, R.; Rothman, D.H. Scaling of a slope: The erosion of tilted landscapes. J. Stat. Phys. 1998, 93, 477–500. [Google Scholar] [CrossRef]
  11. Antonov, N.V.; Kakin, P.I. Scaling in erosion of landscapes: Renormalization group analysis of a model with infinitely many couplings. Theor. Math. Phys. 2017, 190, 193–203. [Google Scholar] [CrossRef] [Green Version]
  12. Duclut, C.; Delamotte, B. Nonuniversality in the erosion of tilted landscapes. Phys. Rev. E 2017, 96, 012149. [Google Scholar] [CrossRef] [Green Version]
  13. Song, T.; Xia, H. Kinetic roughening and nontrivial scaling in the Kardar–Parisi–Zhang growth with long-range temporal correlations. J. Stat. Mech. 2021, 2021, 073203. [Google Scholar] [CrossRef]
  14. Marinari, E.; Parisi, G.; Ruelle, D.; Windey, P. Random Walk in a Random Environment and 1/f Noise. Phys. Rev. Lett. 1983, 50, 1223. [Google Scholar] [CrossRef]
  15. Marinari, E.; Parisi, G.; Ruelle, D.; Windey, P. On the interpretation of 1/f noise. Commun. Math. Phys. 1983, 89, 1–12. [Google Scholar] [CrossRef]
  16. Fisher, D.S. Random walks in random environments. Phys. Rev. A 1984, 30, 960. [Google Scholar] [CrossRef]
  17. Fisher, D.S.; Friedan, D.; Qiu, Z.; Shenker, S.J.; Shenker, S.H. Random walks in two-dimensional random environments with constrained drift forces. Phys. Rev. A 1985, 31, 3841–3845. [Google Scholar] [CrossRef]
  18. Kravtsov, V.E.; Lerner, I.V.; Yudson, V.I. The Einstein relation and exact Gell-Mann-Low function for random walks in media with random drifts. Phys. Lett. A 1986, 119, 203–206. [Google Scholar] [CrossRef]
  19. Honkonen, J.; Pis’mak, Y.M.; Vasil’ev, A.N. Zero beta function for a model of diffusion in potential random field. J. Phys. A Math. Gen. 1988, 21, L835–L841. [Google Scholar] [CrossRef]
  20. Bouchaud, J.-P.; Georges, A. Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications. Phys. Rep. 1990, 195, 127–293. [Google Scholar] [CrossRef]
  21. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  22. Zeitouni, O. Random walks in random environment. In Computational Complexity; Meyers, R., Ed.; Springer: New York, NY, USA, 2012. [Google Scholar]
  23. Révész, P. Random Walk in Random and Non-Random Environments, 3rd ed.; World Scientific Book: Singapore, 2013. [Google Scholar]
  24. Haldar, A.; Basu, A. Marching on a rugged landscape: Universality in disordered asymmetric exclusion processes. Phys. Rev. Res. 2020, 2, 043073. [Google Scholar] [CrossRef]
  25. Hairer, M. Solving the KPZ equation. Ann. Math. 2013, 178, 559–664. [Google Scholar] [CrossRef] [Green Version]
  26. Hairer, M.; Shen, H. A central limit theorem for the KPZ equation. arXiv 2016, arXiv:1507.01237. [Google Scholar] [CrossRef] [Green Version]
  27. Hairer, M. Exactly solving the KPZ equation. In Random Growth Models. Proceedings of Symposia in Applied Mathematics; Damron, M., Rassoul-Agha, F., Seppäläinen, T., Damron, M., Rassoul-Agha, F., Seppäläinen, T., Eds.; American Mathematical Society: Providence, RI, USA, 2018; Volume 75, p. 75. [Google Scholar]
  28. Corwin, I.; Shen, H. Some recent progress in singular stochastic partial differential equations. Bull. Am. Math. Soc. 2020, 57, 409–454. [Google Scholar] [CrossRef] [Green Version]
  29. Barraquand, G.; Corwin, I. Stationary measures for the log-gamma polymer and KPZ equation in half-space. arXiv 2022, arXiv:2203.11037. [Google Scholar]
  30. Falkovich, G.; Gawȩdzki, K.; Vergassola, M. Particles and fields in fluid turbulence. Rev. Mod. Phys. 2001, 73, 913–975. [Google Scholar] [CrossRef]
  31. Pruessner, G. Self-Organized Criticality: Theory, Models and Characterisation; Cambridge University Press: Cambridge, MA, USA, 2012. [Google Scholar]
  32. Bak, P.; Tang, C.; Wiesenfeld, K. Self-organized criticality: An explanation of the 1/f noise. Phys. Rev. Lett. 1987, 59, 381–384. [Google Scholar] [CrossRef]
  33. Bak, P.; Tang, C.; Wiesenfeld, K. Self-organized criticality. Phys. Rev. A 1988, 38, 364–374. [Google Scholar] [CrossRef]
  34. Maslov, S.; Tang, C.; Zhang, Y.-C. 1/f Noise in Bak-Tang-Wiesenfeld Models on Narrow Stripes. Phys. Rev. Lett. 1999, 83, 2449–2452. [Google Scholar] [CrossRef] [Green Version]
  35. Vasiliev, A.N. The Field Theoretic Renormalization Group in Critical Behaviour Theory and Stochastic Dynamics; Chapman & Hall/CRC: Boca Raton, FL, USA, 2004; [Translated from the Russian: St Petersburg Institute of Nuclear Physics: Gatchina, Russia, 1999; ISBN 5-86763-122-2. [Google Scholar]
  36. Adzhemyan, L.T.; Antonov, N.V.; Vasil’ev, A.N. The Field Theoretic Renormalization Group in Fully Developed Turbulence; Gordon & Breach: London, UK, 1999. [Google Scholar]
  37. Antonov, N.V. Anomalous scaling regimes of a passive scalar advected by the synthetic velocity field. Phys. Rev. E 1999, 60, 6691–6707. [Google Scholar] [CrossRef] [Green Version]
  38. Adzhemyan, L.T.; Antonov, N.V.; Vasil’ev, A.N. Renormalization group, operator product expansion, and anomalous scaling in a model of advected passive scalar. Phys. Rev. E 1998, 58, 1823–1835. [Google Scholar] [CrossRef] [Green Version]
  39. Avellaneda, M.; Majda, A. Mathematical models with exact renormalization for turbulent transport II: Non-Gaussian statistics, fractal interfaces, and the sweeping effect. Commun. Math. Phys. 1992, 146, 139–204. [Google Scholar] [CrossRef]
  40. Hwa, T.; Kardar, M. Dissipative transport in open systems: An investigation of self-organized criticality. Phys. Rev. Lett. 1989, 62, 1813–1816. [Google Scholar] [CrossRef]
  41. Hwa, T.; Kardar, M. Avalanches, hydrodynamics and great events in models of sandpiles. Phys. Rev. A 1992, 45, 7002–7023. [Google Scholar] [CrossRef]
Figure 1. Regions of stability of the fixed points (i)–(iv).
Figure 1. Regions of stability of the fixed points (i)–(iv).
Universe 09 00139 g001
Table 1. Canonical dimensions for the action functional (16) and (17).
Table 1. Canonical dimensions for the action functional (16) and (17).
F θ θ h h ν 0 , ν g 0 u 0 g, u μ , m
d F k d d + 2 2 2 ε η 01
d F ω 0 1 110000
d F dd00 ε η 01
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Antonov, N.V.; Gulitskiy, N.M.; Kakin, P.I.; Kerbitskiy, D.A. Random Walk on a Rough Surface: Renormalization Group Analysis of a Simple Model. Universe 2023, 9, 139. https://doi.org/10.3390/universe9030139

AMA Style

Antonov NV, Gulitskiy NM, Kakin PI, Kerbitskiy DA. Random Walk on a Rough Surface: Renormalization Group Analysis of a Simple Model. Universe. 2023; 9(3):139. https://doi.org/10.3390/universe9030139

Chicago/Turabian Style

Antonov, Nikolay V., Nikolay M. Gulitskiy, Polina I. Kakin, and Dmitriy A. Kerbitskiy. 2023. "Random Walk on a Rough Surface: Renormalization Group Analysis of a Simple Model" Universe 9, no. 3: 139. https://doi.org/10.3390/universe9030139

APA Style

Antonov, N. V., Gulitskiy, N. M., Kakin, P. I., & Kerbitskiy, D. A. (2023). Random Walk on a Rough Surface: Renormalization Group Analysis of a Simple Model. Universe, 9(3), 139. https://doi.org/10.3390/universe9030139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop