Front. Phys.Frontiers in PhysicsFront. Phys.2296-424XFrontiers Media S.A.10.3389/fphy.2016.00052PhysicsOriginal Researchd-Dimensional KPZ Equation as a Stochastic Gradient Flow in an Evolving Landscape: Interpretation and Time Evolution of Its Generating FunctionalWioHoracio S.1*RodríguezMiguel A.1GallegoRafael2RevelliJorge A.3AlésAlejandro4DezaRoberto R.41Instituto de Física de Cantabria, Universidad de Cantabria and CSICSantander, Spain2Mathematics Department, Gijón Campus, Universidad de OviedoGijón, Spain3Instituto de Física Enrique Gaviola, Universidad Nacional de Córdoba and CONICET, FaMAF-UNCCórdoba, Argentina4Instituto de Investigaciones Físicas de Mar del Plata, Universidad Nacional de Mar del Plata and CONICET, FCEyN-UNMdPMar del Plata, Argentina
Edited by: Manuel Asorey, University of Zaragoza, Spain
Reviewed by: Rodolfo Cuerno, Universidad Carlos III de Madrid, Spain; Fernando Albuquerque Oliveira, University of Brasília, Brazil
*Correspondence: Horacio S. Wio wio@ifca.unican.es
This article was submitted to Mathematical Physics, a section of the journal Frontiers in Physics
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
The deterministic KPZ equation has been recently formulated as a gradient flow. Its non-equilibrium analog of a free energy—the “non-equilibrium potential” Φ[h], providing at each time the landscape where the stochastic dynamics of h(x,t) takes place—is however unbounded, and its exact evaluation involves all the detailed histories leading to h(x,t) from some initial configuration h0(x,0). After pinpointing some implications of these facts, we study the time behavior of 〈Φ[h]〉t (the average of Φ[h] over noise realizations at time t) and show the interesting consequences of its structure when an external driving force F is included. The asymptotic form of the time derivative Φ˙[h] is shown to be valid for any substrate dimensionality d, thus providing a valuable tool for studies in d > 1.
non-equilibrium growthscaling lawsvariational principlesstochastic methodsFIS2014-59462-PMinisterio de Economía y Competitividad10.13039/5011000033291. Introduction
The KPZ equation for kinetic interface roughening (KIR) [1–4],
∂th(x,t)=ν∇2h(x,t)+λ2[∇h(x,t)]2+ξ(x,t),
where h(x, t) is the interface height and ξ(x, t) a Gaussian noise with
〈ξ(x,t)〉=0,〈ξ(x,t)ξ(x′,t′)〉=2Dδd(x−x′)δ(t−t′),
is nowadays a paradigm of systems exhibiting non-equilibrium critical scaling [5]. In fact—besides standing out as a representative of a large and robust class of microscopic KIR models1, from which the phenomenological parameters ν, λ, D can be computed—it is intimately related to two highly non-trivial problems:
Through the Hopf–Cole transformation
ϕ(x,t)=exp[λ2νh(x,t)],
(1) is isomorphic [8] to the diffusion equation with multiplicative noise obeyed by the restricted partition function of directed polymers in random media (DPRM). Thus, the DPRM problem belongs to the KPZ universality class, and its progress reinforces that of KIR, as much as vice versa.
Via v = −∇h, (1) can be mapped (for λ = 1) into the Burgers equation for a randomly stirred vorticity-free fluid [9, 10]. As a consequence (being the non-linear term in the latter, part of the substantial derivative), λ has to be invariant under scale changes. The invariance of λ under scale changes leads to the remarkable relation
α+z=2
(a signature of the KPZ universality class) in any substrate dimensionality d.
Equation (2) is often attributed to the Burgers equation's Galilean invariance (GI), which in turn translates into KPZ equation's invariance under changes in tilt. This opinion has been repeatedly challenged and in fact, the numerically computed exponents obey (2) in a discrete version of (1) where both Galilean invariance and the (1d–peculiar) fluctuation–dissipation theorem are explicitly broken [11, 12]2.
Since the first and third terms in the r.h.s. of (1) were already present in the model by Edwards and Wilkinson (EW), the innovation came from the tilt-dependent local growth velocity in the second term (or rather, by its interplay with the third one). Under additive uncorrelated local Gaussian white noises, plane (eventually moving) interfaces—which are stable for ξ(x, t) = 0 [14, 15]—develop (in both EW and KPZ models and for large enough systems) into statistically self-affine fractals, whose typical roughness width scales as a power β (called the growth exponent) of the elapsed time. β is in turn the ratio between the interface's Hölder (or roughness) exponent α and the dynamic exponent z, governing the growth in time of the correlation length. In the well understood d = 1 case, α turns out to be that of simple random walk (namely α=12) in both models. However, (2) imposes that z (and thus β) be appreciably different (at least after some crossover time, needed for the second term to take over the first). The initial success of (1) was then due to the consistency of βKPZ with numerical results on microscopic KIR models [2–4].
By about half its lifetime so far, the field was mature enough to approach not only the second moment of the fluctuations in h, but its full statistics [16–20]. In the 1d case, by proposing for h(x, t) the asymptotic behavior
h(x,t)~v∞t+(Γt)1/3χ(x′),
with A:=ν2D, Γ:=A2λ2, and x′:=Ax2(Γt)-2/3, one can see that the stochastic variable χ obeys the Tracy–Widom (TW) statistics of the largest eigenvalue of a random-matrix ensemble, here determined by the dynamics of the substrate3. In the stationary state however, temporal correlations are governed by the Baik–Rains (BR) F0 limit distribution [23], not related to random-matrix theory (RMT).
In the last few years, a handful of exact solutions to (1) have arisen: whereas most of them were inspired by RMT [18, 19, 24, 25], one was born right within the field of stochastic differential equations [26]. On the other hand, the field is by now so mature that experiments can decide on the statistics [20, 27–29]. And pretty much the same occurs with the numerics: in a recent review article [21], the statistics of KPZ itself has been found to agree with the results transposed from DPRM [16–20, 30, 31].
It is the purpose of this article to contribute to lighten up the abovementioned developments from the perspective of a recent variational formulation of the KPZ equation [11, 12, 32–34]. In the following section we derive general consequences from some properties of the free-energy-like functional Φ[h]. Next, we undertake a numerical study of the time dependence of Φ[h] (featuring an estimator of the EW-to-KPZ crossover time and its dependence on parameter λ) in d = 1, 2, and 3, and derive some consequences from the forced case. Then we relate the time behavior of Φ˙[h] with the asymptotic form in (3) and finally, we collect our conclusions.
2. The functional
The intuitive power of disposing of a potential landscape in the study of dynamical systems—as provided by energy in conservative systems, or free energy in equilibrium thermodynamics—and even in other fields (given its closeness to concepts such as “information” or “fitness function”), is beyond doubt. In gradient dissipative systems, it is immediate to find a Lyapunov function to play that role. In general relaxational ones, meeting the integrability conditions may be hard and moreover, the obtained potential landscape partially loses its intuitive power. A breakthrough proposal [35] to circumvent the first drawback—relying on Langevin (thus Fokker–Planck) dynamics having stationary states—was made three decades ago. Unfortunately, KIR is but an example of those cases where this proposal cannot be applied.
In the literature, it is commonly assumed that the KPZ equation cannot be directly obtained from a Hamiltonian. However, as shown in Wio [32], there is no trouble in expressing (1) as a stochastic gradient flow,
∂th(x,t)=-δΦ[h]δh(x,t)+ξ(x,t).
The functional Φ is defined as
Φ[h(x,t)]=∫dx[ν2(∇h)2−λ2∫h0(x,0)h(x,t)dψ(∇ψ)2],
where h0(x, 0) is an arbitrary initial pattern (usually assumed to be constant, in particular h0 = 0). The first term is clearly the Landau–Ginzburg free-energy functional
ΦEW[h(x,t)]=ν2∫dx[∇h(x,t)]2.
associated to the (equilibrium) EW process. Unfortunately, the second one has not an explicit density. Even though at time t, Φ depends only on the field h(x, t), its evaluation requires knowing the detailed history that led from h0 to h(x, t). In other words, retrieving information on Φ (such as its landscape at certain time or the time dependence of its mean value) requires averaging not simply over field configurations h(x) at time t, but over (statistically weighted) trajectories of the field configuration.
Being the KPZ equation a stochastic gradient flow, the functional Φ governing its deterministic component provides the landscape where the stochastic dynamics of h(x) takes place at time t, and in the absence of noise fulfills explicitly the Lyapunov property Φ˙[h]=-[δΦ[h]δh(x,t)]2≤0. Unfortunately, this does not make Φ into a Lyapunov functional, since (as shown in the next paragraph) it is unbounded from below. As it is well-known, Equation (1) has the joint symmetry λ ↔ −λ, h ↔ −h. In the case of Equation (5), the change λ ↔ −λ—together with the (∇h)2 dependence—implies a specular reflection h ↔ −h with respect to the h axis.
In fact—given a reference configuration h0(x, t0)—an expansion of Φ[h] in powers of δh(x, t) = h(x, t) − h0(x, t0) exactly terminates at the third order:
and the definitions (∇h0)i : = (h0)i+1 − (h0)i−1, (∇2h0)i:=(h0)i+1+(h0)i-1-2(h0)i.
Note that δhi > ν/λ would produce a negative parabolic potential, which implies local instabilities. Such a diffusive instability in [ΔΦ2]i tells us that the above procedure makes sense only for short times i.e., it provides an “instantaneous landscape” of the potential. However, this is all we need to extract important information on the qualitative aspects of the dynamics. On one hand—by associating the height δhi at site i with the coordinate of a particle submitted to the potential [ΔΦ]i—we retrieve a mechanical picture that helps understanding the local dynamics. We notice that [ΔΦ]i is at most linear in δhi. But whereas the coefficient of the first term is a constant, the other two terms are some kind of “mean-field potentials.”
Figures 1, 2 depict different contributions to [ΔΦ]i, as functions of δhi−1 and δhi+1. Figure 1 shows that whereas for EW, [ΔΦ2]i has a concave parabolic shape with constant curvature, for KPZ this curvature depends on the sign of δhi through the contribution from δ3Φ[h(x, t)] in (6). Figure 2 illustrates the breakdown of the “horizontal” symmetry δhi+1 ↔ δhi−1 by showing 3d plots of [ΔΦ3]i when (∇h0)i is either positive or negative.
Left: 3d plot of [ΔΦ]i for the EW case, with ν = 1. Right: projection on the plane δhi+1 = −δhi−1. Solid black line: EW case, with ν = 1. Remaining lines: KPZ case (λ = 1) illustrating the breakdown of the “vertical” symmetry δhi ↔ − δhi (red dash-dotted line: δhi = 1; blue dashed line: δhi = −1).
3d plots of [ΔΦ3]i—the third term of the local potential variation for particle i—in the KPZ case, adopting ν = λ = 1. Left: [∇h0]i = 1; Right: [∇h0]i = −1. Both cases exhibit breakdown of the “horizontal” symmetry δhi+1 ↔ δhi−1, but with different signs.
Rather than the quantitative information these figures may provide, some essential aspects—as the local symmetries involved in each case and the dynamical consequences of their breakdown—are worth discussing. We focus on three facts from (7):
Since the coefficient of the first-order term is symmetric (similarly to what occurs in the EW case, where this coefficient is proportional to the curvature), the fluctuating local linear field it produces has no consequence whatsoever on the global effect of the potential. A breakdown of this symmetry produces a final global effect in the potential that, if it is uniform, gives rise to a linear potential.
When the symmetry under “vertical” flip (namely, changing the sign of δhi) is broken—as it occurs in the second term of (7), one would expect some asymmetry in the probability density of δhi.
When the symmetry under “horizontal” flip (namely, changing δhi+1 ↔ δhi−1) holds, diffusive propagation of fluctuations (and thus z = 2) is expected, as occurs in the EW case. If it is broken by including local ballistic effects—as in the third term of (7)—superdiffusive propagation of fluctuations becomes possible, offering an explanation to the origin of the value z = 3/2 observed in KPZ dynamics.
Note that in the EW case these symmetries hold, while in the KPZ all of them are broken. These facts lend themselves to insightful interpretations that we study in what follows.
Before delving however into the thorough analysis of Φ[h] to which this article is devoted, we find it convenient to introduce a genuine Lyapunov functional for KPZ, to see its pros and cons. The functional
F[h]=ν2∫exp[λνh(x,t)][∇h(x,t)]2dx,
of Wio [32], allows the KPZ equation to be written as
∂∂th(x,t)=-Γ[h]δF[h]δh(x,t)+ξ(x,t),
with Γ[h] given by
Γ[h]=exp[-λνh(x,t)].
It is straightforward to verify that (i) the Lyapunov property F˙[h]≤0 holds, (ii) the EW free-energy functional is recovered for λ = 0, and (iii) 〈F[h]〉t behaves asymptotically as f(t)exp[λνv∞t], with f(t) far milder than an exponential [2, 3]. Moreover, numerical simulations show that v∞ ∝ λ, so the asymptotic behavior results ~exp[λ2t].
Since (9) is not of gradient form, F[h] does not provide an intuitive landscape for the stochastic dynamics of h(x, t). So hereafter, we focus on the analysis of Φ[h]. As stated before, its form provokes several thoughts:
2.1. On normalization
The first-order contribution to the potential in (7) is composed of two terms: one is proportional to the curvature of the initial shape [∇2h0]i and the other (which appears only in the KPZ case) is always positive and given by [(∇h0)2]i. This term breaks the symmetry observed in the EW case and induces a potential which is linear in δhi. The dynamics is then analog to that of a particle in a constant gravitational field [36]. In this case, the mean position of the particle moves with constant velocity and its fluctuations obey a diffusive behavior. The probability density of the particle keeps evolving (becoming flatter around its mean position) and a stationary solution (and even normalization of this probability) is only possible if there exists an external boundary.
The equilibrium particle distribution in a constant gravitational field would in principle be the stationary limit of the (exact) solution of the Fokker–Planck equation (FPE) associated to the Langevin problem
x˙=F+ξ(t),〈ξ(t)〉=0,〈ξ(t)ξ(t′)〉=Dδ(t-t′),
with the initial condition P(x, t → t0|x0, t0) = δ(x − x0), if this stationary distribution could be normalized! As known, the physics of this simple problem dictates that there is always a boundary. In fact, the correct interval has a finite limit in the sense of the decreasing potential, thus allowing for normalization [37, 38].
In this case, it is not a stationary but an asymptotic solution to (10) what we look for. Forcing the condition ∂tP[h] = 0 leads to Pas[h] ∝ exp(−Φ[h]/D).
As in the simple analogy before, the pdf in the potential indicated by (7) cannot be normalized except by considering an interval (−∞, hm].
2.2. On the statistics
The symmetry under “vertical” flip is broken by the second term of (7) in the KPZ case. When δhi is positive (see Figure 1), the curvature 2(ν-λ3δhi) of the parabolic potential decreases with respect to the EW value, which means that larger fluctuations are possible of the distance between neighbors. Conversely, when δhi is negative, the opposite effect occurs. The final global effect will be some asymmetry in the probability distribution of hi around its mean value, which is in agreement with recent exact results in 1d, where it has been shown that the pdf is not a Gaussian but a Tracy–Widom distribution [16–19]. Recall that since δhi > ν/λ would produce a negative parabolic potential—meaning local instabilities—the expansion is valid only for short times, such that the potential remains positive.
It is thus worthwhile stressing that the well-known pdf [2, 3]
Pstat[h]∝exp{-ν2D∫0L(∂xh)2dx},
makes sense only for a finite 1d system with periodic boundary conditions and for times larger than t ~ Lz, the saturation one.
2.3. Potential landscape's long memory
The second term in (5), namely
−λ2∫dx∫h0(x,0)h(x,t)dψ(∇ψ)2,
is not a functional integral of the form ∫∏x dψ(x) [39], but a kind of “grand mean” of (∇h)2. For a suitably small time τ, the integral over ψ at fixed x can be evaluated by resort to the Mean Value Theorem, yielding
∫h0(x,0)h(x,τ)dψ(∇ψ)2≈[h(x,τ)−h(x,0)][∇h˜(x)]2,
where h~(x) is some intermediate value between h(x, τ) and h(x, 0). For finite time t, a discretization with step τ and index μ yields
with h~μ(x) intermediate between hμ(x) and hμ+1(x), and h˙(x, s): = limτ→0[hμ+1(x) − hμ(x)]/τ. This allows us to write the potential in the form
Φ[h]=∫dx{ν2(∇h)2−λ2∫0tdsh˙(x,s)[∇h(x,s)]2},
which highlights the fact that even though the field h(x, t) obeys a Fokker–Planck equation (thus it is Markovian, albeit without stationary regime), it evolves at every time in an adaptive potential landscape, which itself evolves according to the whole history of h(x, t) and has thus a very long (“infinite”) memory. The already mentioned dependence of the asymptotic statistics on the dynamics of the substrate might probably be traced back to this fact [16–20].
3. Numerical results for Φ[h]
Spectral methods pay back the overhead of Fourier transforming h(x) forth and back every time t, with a discrete approximation of ∇h of the order of the system's size. They have proved to be more stable and reliable than pure real-space finite-differences schemes in the integration of some non-linear growth equations [40–42]. Using a spectral method to integrate (1), we have analyzed the time behavior of 〈Φ[h]〉t (the average of Φ over noise realizations, at time t) as given by (14).
Figure 3 displays the time dependence of the NEP's average over 100 samples, for systems in 1d (size 1,024), 2d (size 1282), and 3d (size 643), and several values of λ. For any λ > 0 there is a maximum, where the non-linear (KPZ) term overcomes the linear (EW) one. Only for λ = 0 (corresponding to EW) is there true saturation. The λ = 0.01 case shown here (visually resembling a saturation behavior) attains its maximum outside the plotted range.
Time behavior of Φ[h], averaged over 100 samples, in (A) 1d (size 1024), (B) 2d (size 1282), and (C) 3d (size 643). The symbols (denoting a subset of the simulation points) indicate the values adopted for λ in each curve. ■ : λ = 0.01, ▴ : λ = 0.10, ◇ : λ = 0.20, ⋆ : λ = 0.30, • : λ = 0.50, ▾ : λ = 1.00. Inset: Detail showing the existence of maxima in 2d for any λ > 0 (the same occurs in 1d and 3d).
Figure 4—as well as (16) and (17) below—shows that past this maximum, 〈Φ[h(t)]〉t ~ A − Bt with B ~ λ2. This result shows that—due to the correlations and in an effective way—the NEP behaves as having only a linear dependence on h (just as in the toy model discussed before [37, 38]).
Mean λ-dependence over 100 samples of the NEP's asymptotic slope. Here, the symbols denote the substrate's dimensionality: ◇ 1d (size 1024), ■ 2d (size 1282), and • 3d (size 643). The dashed lines (besides guiding the eye) are best fits with a λb (b = 2.01 in d = 1 and 2, and b = 2.03 in 3d).
Figure 5 shows the dependence on λ of the time txovr at which the maximum occurs, i.e., at which the non-linear (KPZ) contribution overcomes the linear (EW) one, making it a proxy of the EW-to-KPZ crossover time [43–45]. Although roughly compatible with a txovr~λ-α law (α = 3 in Guo et al. [43] and α = 4 in Forrest and Toral, [44] for d = 1) and in agreement with the results of previous studies exploiting the time behavior of the stochastic action [46, 47], this dependence becomes milder as d increases.
Mean λ-dependence over 100 samples of the time of occurrence of the NEP's maximum, in •: 1d (size 1024), ■: 2d (size 1282), and ◇: 3d (size 643). Dashed lines: best fits with a λ−b yield b = 1.14 in 1d, b = 0.98 in 2d, and b = 0.79 in 3d.
3.1. External field
By including an external driving force F, it is possible to capture some new aspects of the dynamics. Equation (5) now reads
Figure 6 shows (on a much longer timescale than in Figure 3, so that the already discussed initial rise due to the EW term cannot be appreciated) the NEP's time behavior for both positive and negative values of F. The F > 0 case brings no novelty; but when negative enough forces are applied, a reversion is observed of the NEP's time behavior, corresponding to a reversion of the front motion 〈h〉 and supporting the analogy to an activation-like behavior (which could be guessed from the cubic-like form of the NEP). As already stated however, this behavior is not of the exponential (Kramers-like) form. Moreover the abovementioned diffusive instability, now written as
Φ[h]=∫dx[(ν2−λδh6)(∇δh)2−Fδh],
limits its validity to a range of F. If F < 0, the evolution toward negative values of δh(x, t) proceeds mainly along the almost zero-measure subset of directions where ∇δh(x, t) = 0. However, as soon as—due to the presence of noise—δh(x, t) slightly depart from those directions, the effect of the NEP is such that it will drive the system in the opposite direction! At any rate, this will only happen for F larger than some threshold value. Thus, for F = 0, the time behaviors of 〈Φ[h]〉 and 〈h〉 are correlated, even though the dispersion of the NEP [48] and that of the front differ (the former continuously increases with time, while the latter saturates for t ≥ Lz). These aspects will be fully clarified in the next section. For the time being let us remind that in the EW case and without external force, 〈h〉 = 0. If a force F is applied, the EW NEP acquires the form
ΦEW[h(x,t)]=ν2∫dx{[∇h(x,t)]2−Fh(x,t)}.
NEP's long-time behavior in 1d (λ = 1), from 100 samples of size 256. Upper frame: For F < −0.3, 〈Φ[h]〉 increases with time in the observation interval. Lower frame: detail showing that for F just slightly lower than −0.3, 〈Φ[h]〉 first increases with time (on a much longer timescale than the initial rise due to the EW term) and then begins to decrease, supporting the hypothesis of an activation-like behavior.
When the linear term largely exceeds the quadratic one, the EW equation becomes h˙ = νF, so 〈h〉 → νFt. This effect seems to indicate that including an external force reinforces the statistical weight of those directions such that ∇δh(x, t) = 0, which for F = 0 have (almost) zero measure. As a matter of fact, neither the inclusion of the external force F nor the effects shown in Figure 6 change the scaling behavior.
4. Some asymptotic results
As known in the EW case [2–4] and shown in Figure 3, the average 〈 ∫ dx(∇h)2〉 over noise realizations saturates to a constant C0 at long times. Using in (15) the discrete expression for Φ[h]—namely (12), with the approximation in (13b)—we can write
In the discrete representation, (3) implies δh ~ v∞τ at long times, where τ corresponds to δt and v∞ is the asymptotic front velocity. Inserting this expression into the previous one and replacing 〈 ∫ dx(∇h)2〉 by its bound C0, we get
〈Φ[h]〉=ν2C0-(λ2C0+F)v∞t.
This result, valid in any dimensionality, completely agrees with the numerical ones presented above. In particular, 〈Φ[h]〉 will decrease at long times only if λC0 + F ≥ 0. If λC0 + F < 0, 〈Φ[h]〉 will increase asymptotically, as shown in Figure 6. In addition, it is well-known that v∞ ∝ λ [3]. Hence (particularly in the F = 0 case) we find the λ2 dependence for the slope of 〈Φ[h]〉 vs. t, as obtained in simulations.
We close this section by looking at the asymptotic behavior from another point of view. Here we restrict to d = 1, the continuum representation, and the F = 0 case. As already indicated, the statistics of the stochastic variable χ in (3), is of paramount importance [16, 17]. In addition to the forms of Φ[h] given by (5) and (14), we can also obtain an expression for Φ˙[h] as
as the 2nd and 3rd terms decay very fast. On the other hand,
(∂xh)2~A24(Γt)−2/3(χ′)2,∂x2h~A24(Γt)−1χ″.
Exploiting these relations we obtain
Φ˙[h]~-∫dx′A2[νχ″+λ2(χ′)2]v∞.
If we look at long times, where the non-linear contribution dominates, we have
Φ˙[h]~-∫dx′Aλ4(χ′)2v∞.
Since v∞ ∝ λ, we get right away the expected result: the constant goes as λ2.
We now show why 〈Φ[h]〉 and 〈h〉 are correlated. At long times, when 〈∫dx(∇δh)2〉→C0 and 〈h〉 → v∞t, (13b)—now in terms of δh—“factorizes” out and thus δ3Φ[h(x,t)]→-λC0〈h〉. Hence the second term in (5)—which already dominates over the first one—behaves in the same way as the term in EW corresponding to an external force:
For any substrate dimensionality d, the KPZ equation can be expressed as a gradient flow by means of the functional Φ[h(x, t)] in (5). This is a state functional, i.e., it depends explicitly only on h(x, t). However, evaluating Φ[h] for a given field configuration h(x, t) requires integrating over its whole story, with initial condition h0(x). So Φ[h] carries information on the whole growth story. If (5) is discretized according to (∇h)2→12[(hi+1-hi)2+(hi-1-hi)2] [11, 12], a cubic polynomial in the {hi}—containing hi3 terms—is obtained.
Expansion (6) provides an “instantaneous landscape” of Φ[h(x, t)]. Even using (∇h)2→(hi+1-hi-1)2, its discretization also yields a cubic polynomial (in the {δhi}) but now the absence of δhi3 terms enables a clearer understanding of “horizontal” and “vertical” symmetry breakdowns when λ ≠ 0—which may in turn be related to the asymmetry of the Tracy–Widom distribution—and brings an explanation to the shift in the dynamical exponent z (a change from diffusive to superdiffusive behavior).
Even though Φ[h] fulfills the Lyapunov property, its functional form is for short times that of a cubic polynomial, and is thus unbounded below. The KPZ equation turns out to be a high-dimensional (non Kramers) escape-like problem, and its asymptotic pdf cannot be normalized on the whole h–space (except for fixed and normalized forms of h0). We have moreover provided a bridge between the observed asymptotic Tracy–Widom statistics, and the behavior of Φ (which as argued, provides an instantaneous landscape for the stochastic dynamics).
Given an ensemble of stories h(x, t), one can attempt to calculate the time dependence of the lower moments of Φ[h]. Our focus in this work is on 〈Φ[h]〉t, the average of Φ[h] over noise realizations at time t. We have shown both numerically (Figures 3, 4) and analytically in (17) that in the absence of external forcing (F = 0), 〈Φ[h]〉t → A − Bt as t → ∞ (in fact, surprisingly soon!), with B ~ λ2. We have moreover traced back the origin of this asymptotic linear dependence on t to two known effects: (a) 〈h〉 → v∞t, (b) 〈∫dx(∇h)2〉→C0, which lead to factorization of the “KPZ term” in 〈Φ[h]〉t.
The time tmax at which 〈Φ[h]〉t goes through a maximum, can be regarded as an estimation of the crossover time txovr from the EW regime to the KPZ one. We have found that tmax~λ-1, in agreement with previous results obtained by exploiting the time behavior of the stochastic action [46, 47].
Other issues worth remarking are that—because of the abovementioned linear dependence of 〈Φ[h]〉t on t—the behaviors of 〈Φ[h]〉t and 〈h〉 are correlated, and that (both with and without external forcing) the NEP's time behavior indicates a (non Kramers) activation-like phenomenon. In a forthcoming work [48] we undertake a thorough analysis of the time behavior of the first few higher moments of Φ[h], namely the NEP's dispersion σ(t), skewness γ1(t) and kurtosis K(t). An interesting point (left to further work) would be to “interpret” the aging phenomenon [49, 50] from the NEP perspective.
We have discussed some strengths of a novel tool for describing the KPZ dynamics, the NEP approach. This innovative framework is expected to contribute answering questions that remain open by today, and handling the KPZ problem from still another perspective [13, 36].
Author contributions
All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of interest statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Support by MINECO (Spain), under project No. FIS2014-59462-P, is acknowledged by HW, MR, and RG; support by CONICET, UNMdP and UNC (Argentina), by RD and JR. The authors thank J.M. López, R. Toral, C. Escudero, P. Colet and E. Hernández-García for fruitful discussions.
ReferencesKardarMParisiGZhangYC. Dynamic scaling of growing interfaces. (1986) 56:889–92. 10.1103/PhysRevLett.56.88910033312Halpin-HealyTZhangYC. Kinetic roughening phenomena, stochastic growth, directed polymers and all that: aspects of multidisciplinary statistical mechanics. (1995) 254:215–414. 10.1016/0370-1573(94)00087-JBarabásiALStanleyHE. . Cambridge, UK: Cambridge University Press (1995).KrugJ. Origins of scale invariance in growth processes. (1997) 46:139–282. 10.1080/00018739700101498LesneALaguësM. . Berlin: Springer (2012).RamascoJJLópezJMRodríguezMA. Generic dynamic scaling in kinetic roughening. (2000) 84:2199–202. 10.1103/PhysRevLett.84.219911017243WioHSDezaRRLópezJM. . Singapore: World Scientific (2012).HoldenHLindstrømTØksendalBUbøeJZhangTS. The Burgers equation with noisy force and the stochastic heat equation. (1994) 19:119–41. 10.1080/03605309408821011FogedbyH. Morphology and scaling in the noisy Burgers equation: soliton approach to the strong coupling fixed point. (1998) 80:1126–9. 10.1103/PhysRevLett.80.1126FogedbyH. Canonical phase-space approach to the noisy Burgers equation. (1999) 60:4950–3. 10.1103/PhysRevE.60.495011970362WioHSRevelliJADezaRREscuderoCdeLa Lama MS. KPZ equation: galilean-invariance violation, consistency, and fluctuation–dissipation issues in real-space discretization. (2010) 89:40008. 10.1209/0295-5075/89/4000820866543WioHSRevelliJADezaRREscuderoCdeLa Lama MS. Discretization-related issues in the Kardar–Parisi–Zhang equation: consistency, Galilean-invariance violation, and fluctuation–dissipation relation. (2010) 81:066706. 10.1103/PhysRevE.81.06670620866543RodriguesEAMelloBAOliveiraFA. Growth exponents of the etching model in high dimensions. (2015) 48:035001. 10.1088/1751-8113/48/3/035001KapralRLiviROppoGLPolitiA. Dynamics of complex interfaces. (1994) 49:2009–22. 10.1103/PhysRevE.49.200925518729XuZHanHWuX. Numerical method for the deterministic Kardar–Parisi–Zhang equation in unbounded domains. (2006) 1:479–93.PrähoferMSpohnH. Statistical self-similarity of one-dimensional growth processes. (2000) 279:342–52. 10.1016/S0378-4371(99)00517-8PrähoferMSpohnH. Universal distributions for growth processes in 1+1 dimensions and random matrices. (2000) 84:4882–5. 10.1103/PhysRevLett.84.488210990822SasamotoTSpohnH. One-dimensional Kardar–Parisi–Zhang equation: an exact solution and its universality. (2010) 104:230602. 10.1103/PhysRevLett.104.23060220867222CalabresePLeDoussal P. Exact solution for the Kardar–Parisi–Zhang equation with flat initial conditions. (2011) 106:250603. 10.1103/PhysRevLett.106.25060321770622TakeuchiKA. Crossover from growing to stationary interfaces in the Kardar–Parisi–Zhang class. (2013) 110:210604. 10.1103/PhysRevLett.110.21060423745853Halpin-HealyTTakeuchiKA. A KPZ cocktail—shaken, not stirred…Toasting 30 years of kinetically roughened surfaces. (2015) 160:794–814. 10.1007/s10955-015-1282-123473312CarrascoISSTakeuchiKAFerreiraSCOliveiraTJ. Interface fluctuations for deposition on enlarging flat substrates. (2014) 16:123057. 10.1088/1367-2630/16/12/123057BaikJRainsEM. Limiting distributions for a polynuclear growth model with external sources. (2000) 100:523–41. 10.1023/A:1018615306992CorwinI. The Kardar–Parisi–Zhang equation and universality class. (2012) 1:1130001. 10.1142/S201032631130001424229127ImamuraTSasamotoT. Exact solution for the stationary Kardar–Parisi–Zhang equation. (2012) 108:190603. 10.1103/PhysRevLett.108.19060323003020HairerM. Solving the KPZ equation. (2013) 178:559–664. 10.4007/annals.2013.178.2.4TakeuchiKASanoM. Universal fluctuations of growing interfaces: evidence in turbulent liquid crystals. (2010) 104:230601. 10.1103/PhysRevLett.104.23060120867221TakeuchiKASanoMSasamotoTSpohnT. Growing interfaces uncover universal fluctuations behind scale invariance. (2011) 1:34. 10.1038/srep0003422355553TakeuchiKASanoM. Evidence for geometry-dependent universal fluctuations of the Kardar–Parisi–Zhang interfaces in liquid-crystal turbulence. (2012) 147:853–90. 10.1007/s10955-012-0503-0Halpin-HealyT. (2+1)-dimensional directed polymer in a random medium: scaling phenomena and universal distributions. (2012) 109:170602. 10.1103/PhysRevLett.109.17060223215169Halpin-HealyT. Extremal paths, the stochastic heat equation, and the three-dimensional Kardar–Parisi–Zhang universality class. (2013) 88:042118. 10.1103/PhysRevE.88.04211824229127WioHS. Variational formulation for the KPZ and related kinetic equations. (2009) 19:2813–21. 10.1142/S0218127409024505WioHSEscuderoCRevelliJADezaRRdeLa Lama MS. Recent developments on the Kardar–Parisi–Zhang surface-growth equation. (2011) 369:396–411. 10.1098/rsta.2010.025921149379WioHSDezaRREscuderoCRevelliJA. Invited review: KPZ. Recent developments via a variational formulation. (2013) 5:050010. 10.1098/rsta.2010.0259GrahamR. Weak noise limit and nonequilibrium potentials of dissipative dynamical systems, In: TirapeguiEVillarroelD editors. . Dordrecht: D. Reidel (1987). p. 271–90.SakaHBengtsonRDReichlLE. Relaxation of Brownian particles in a gravitational field. (2009) 77:240–3. 10.1119/1.303902922740596vanKampen NG. . Amsterdam: North-Holland (1990).GardinerCW. . Berlin: Springer (2004).LangoucheFRoekaertsDTirapeguiE. . Dordrecht: Reidel (1982).GallegoRCastroMLópezJM. Pseudospectral versus finite-difference schemes in the numerical integration of stochastic models of surface growth. (2007) 76:051121. 10.1103/PhysRevE.76.05112118233637GiadaLGiacomettiARossiM. Pseudospectral method for the Kardar–Parisi–Zhang equation. (2002) 65:036134. 10.1103/PhysRevE.65.03613411909192GallegoR. Predictor–corrector pseudospectral methods for stochastic partial differential equations with additive white noise. (2011) 218:3905–17. 10.1016/j.amc.2011.09.038GuoHGrossmannBGrantM. Crossover scaling in the dynamics of driven systems. (1990) 41:7082–5. 10.1103/PhysRevA.41.70829903136ForrestBMToralR. Crossover and finite-size effects in the (1+1)-dimensional Kardar–Parisi–Zhang equation. (1993) 70:703–20. 10.1007/BF01053591OliveiraTJDechoumKRedinzJAAarãoReis FDA. Universal and nonuniversal features in the crossover from linear to nonlinear interface growth. (2006) 74:011604. 10.1103/PhysRevE.74.01160416907104WioHSDezaRRRevelliJAEscuderoC. A novel approach to the KPZ dynamics. . (2013) 44:889–98. 10.5506/APhysPolB.44.88927786254WioHSRevelliJAGallegoRRodríguezMADezaRR. Estimation of the crossover time in arbitrary dimensions for the KPZ equation, through the analysis of its stochastic action. (2016).WioHSGallegoRRodríguezMAAlésADezaRRRevelliJA. Generating functional for the KPZ equation in arbitrary dimensions: time evolution of its higher cumulants. (2016).BustingorryS. Aging dynamics of non-linear elastic interfaces: the Kardar–Parisi–Zhang equation. (2007) 10:P10002. 10.1088/1742-5468/2007/10/P10002HenkelMNohJDPleimlingM. Phenomenology of aging in the Kardar–Parisi–Zhang equation. . (2012) 85:030102(R). 10.1103/PhysRevE.85.03010222587028
1We deliberately leave aside KIR models exhibiting diverse kinds of anomalous scaling [6, 7].
2Both discretizations of the KPZ equation and microscopic KIR models generically break GI. But even in the few cases explicitly built to preserve GI, the periodic boundary conditions (PBC) customarily adopted in numerical simulations explicitly break down tilting invariance (note in passing that spectral methods rely more than any other on PBC). A recent result supporting this is that cellular automata models for etching have been shown to obey (1) up to d = 6 [13].
3The Gaussian orthogonal ensemble (GOE) if the substrate is static, the Gaussian unitary (GUE) one if it evolves [21, 22].