The KPZ Equation of Kinetic Interface Roughening: A Variational Perspective

TL;DR

This paper explores the KPZ equation, a model for interface roughening, using a variational approach to better understand its behavior across different dimensions.

Contribution

The paper introduces a variational approach for analyzing the KPZ equation, enabling analytical progress in any substrate dimensionality.

Findings

The non-equilibrium potential (NEP) shows distinct scaling behavior with the nonlinearity parameter λ in different dimensions.

Numerical evidence shows that the NEP's asymptotic behavior strongly depends on initial conditions in both KPZ and GB models.

The variational approach provides insights into the stochastic thermodynamics of interface roughening processes.

Abstract

Interfaces of rather different natures—as, e.g., bacterial colony or forest fire boundaries, or semiconductor layers grown by different methods (MBE, sputtering, etc.)—are self-affine fractals, and feature scaling with universal exponents (depending on the substrate’s dimensionality d and global topology, as well as on the driving randomness’ spatial and temporal correlations but not on the underlying mechanisms). Adding lateral growth as an essential (non-equilibrium) ingredient to the known equilibrium ones (randomness and interface relaxation), the Kardar–Parisi–Zhang (KPZ) equation succeeded in finding (via the dynamic renormalization group) the correct exponents for flat d=1 substrates and (spatially and temporally) uncorrelated randomness. It is this interplay which gives rise to the unique, non-Gaussian scaling properties characteristic of the specific, universal type of non-equilibrium roughening. Later on, the asymptotic statistics of process h(x) fluctuations in the scaling regime was also analytically found for d=1 substrates. For d>1 substrates, however, one has to rely on numerical simulations. Here we review a variational approach that allows for analytical progress regardless of substrate dimensionality. After reviewing our previous numerical results in d=1, 2, and 3 on the time evolution of one of the functionals—which we call the non-equilibrium potential (NEP)—as well as its scaling behavior with the nonlinearity parameter λ, we discuss the stochastic thermodynamics of the roughening process and the memory of process h(x) in KPZ and in the related Golubović–Bruinsma (GB) model, providing numerical evidence for the significant dependence on initial conditions of the NEP’s asymptotic behavior in both models. Finally, we highlight some open questions.

Figures

Full text

1. Introduction

Originally formulated in 1986 to describe kinetic interface roughening, the nonlinear stochastic partial differential equation known as the Kardar–Parisi–Zhang (KPZ) equation readily became a paradigm of non-equilibrium dynamic critical scaling systems, defining as its canonical member the rather comprehensive KPZ universality class. This foundational work [1] established a profound connection between the dynamics of interface growth and the broader field of non-equilibrium physics. The mathematical formulation of the KPZ equation initially posed significant challenges, as its strong nonlinearity exceeded the scope of existing methods for studying stochastic partial differential equations. A major breakthrough in addressing the “well-definedness” of this equation occurred in 2013 [2]. The resolution of this fundamental mathematical problem was essential for advancing the analytical understanding of complex stochastic systems and extracting deeper physical insights.

Universality is a cornerstone concept in statistical physics, asserting that diverse physical systems and mathematical models can exhibit identical statistical behaviors in their long-time and large-scale limits, irrespective of their specific microscopic details. The KPZ equation describes and connects a broad spectrum of significant stochastic growth-like processes across physics, chemistry, and biology, spanning both classical and quantum systems. This wide applicability means that the KPZ universality class offers a powerful coarse-grained predictive framework, enabling the inference of macroscopic behaviors across diverse systems without requiring detailed knowledge of their microscopic interactions. Researchers can infer macroscopic behaviors and statistical properties across vastly different physical, chemical, and biological systems without needing to delve into their specific microscopic interactions. This “surprising and broad universality” of the one-dimensional KPZ equation (even appearing in deterministic Hamiltonian systems) further elevates its status from a mere model to a fundamental principle that unifies the study of non-equilibrium phenomena, suggesting that common mathematical structures dictate emergent behavior irrespective of the underlying physical substrate.

The study of the KPZ equation has revealed deep and unexpected mathematical connections that have been pivotal in its analytical understanding and in broadening its perceived universality, due to its relation with both the Burgers equation of turbulence and the restricted partition function of directed polymers in random media (DPRM). In particular, the KPZ equation and its associated universality class have remarkably extended their reach into the domain of biological growth phenomena, providing a framework for understanding complex pattern formation and spatial dynamics in living systems. The KPZ equation is being increasingly applied to movement ecology due to its ability to model complex stochastic processes that resemble ecological patterns. Since the scaling behavior of the surface of Escherichia coli and Bacillus subtilis colonies growing on agar plates was established to be self-affine and consistent with models in the so-called “KPZ universality class” [3], there have been numerous applications of the KPZ equation in movement ecology [4,5], among which we might cite the following:

Modeling nonlinearities in ecological systems. The KPZ equation’s inclusion of nonlinear growth terms is particularly useful for modeling ecological processes where interactions between individuals or between species and their environments lead to complex, emergent movement patterns. For example, the propagation of invasive species, herd dynamics, or predator–prey interactions often exhibit such nonlinearities, and KPZ can capture these effects efficiently. A foundational exploration of movement patterns, emphasizing the need for nonlinear models like the KPZ to understand large-scale spatial patterns, can be found in [6].Spatial pattern formation in ecosystems. KPZ-like models can explain spatial clustering and aggregation, phenomena common in ecosystems due to habitat fragmentation or social behaviors in species. In this sense, the KPZ equation aligns well with empirical observations in movement ecology, where spatially heterogeneous environments cause uneven dispersal patterns. Book [7] elaborates on spatial pattern formation and how stochastic models, including those similar to KPZ, apply to understanding the distribution and aggregation of organisms.Capturing anomalous diffusion in animal movement. Traditional diffusion models often fail to capture the complexity of animal movement, which frequently exhibits subdiffusive or superdiffusive behavior due to factors like resource hotspots or behavioral adaptations. The KPZ equation, with its growth and roughening terms, provides a more flexible framework for modeling these deviations from classical diffusion. The article [8] reviews the inadequacy of simple diffusion models and how alternative stochastic approaches like KPZ are better suited to describe animal search behaviors that diverge from Brownian assumptions.Modeling stochastic processes in animal movement. The KPZ equation is an excellent tool for capturing random yet structured movement patterns, such as those seen in animal migrations or dispersal behaviors. Unlike simple random walks, the KPZ framework can describe anisotropic movement, a key aspect of ecological processes where environmental and physiological conditions lead to non-uniform patterns. The article [9] discusses the role of stochastic models, like KPZ, in modeling animal movement patterns beyond simple Brownian motion.Interfacing with statistical mechanics for predictive ecology. KPZ’s roots in statistical mechanics make it a natural fit for ecological models that require predictions based on system dynamics under uncertainty. By using KPZ, researchers can leverage known statistical properties to predict critical transitions, tipping points, or responses to environmental changes. The article [10] discusses the interface between movement ecology and statistical mechanics, highlighting how frameworks like KPZ facilitate a predictive approach to ecological modeling.In sum, the KPZ equation offers a robust framework for addressing the complex, nonlinear, and stochastic characteristics of movement patterns in ecology. Its applicability extends beyond traditional diffusion models, making it valuable for studying the spatial and temporal heterogeneity observed in animal populations. By leveraging the KPZ equation, ecologists can create models that more accurately reflect real-world movement dynamics, supporting better predictive and management tools in ecological research [11,12,13,14]. The KPZ universality class is characterized by definite values of either the roughness exponent α or the dynamical exponent z, and by the relation $[eqn]$ . If L denotes the substrate’s length along any dimension, then the global width

$$[eqn]$$

scales as $[eqn]$ . Here h is the local height (flat substrates) or radius (curved substrates), and $[eqn]$ denote spatial average and $[eqn]$ the ensemble one. Moreover, the correlation length $[eqn]$ on the interface scales as $[eqn]$ . For finite L, the scaling regime ends at saturation, when $[eqn]$ . As $[eqn]$ , it is a dynamical critical phenomenon. Putting it all together, it turns out that in the scaling regime $[eqn]$ , where $[eqn]$ is called the growth exponent.

Several kinetic models (such as the Eden model for bacterial colonies) yielded values consistent with the set $[eqn]$ , $[eqn]$ (hence $[eqn]$ ) for flat 1d substrates. The two necessary ingredients can be abstracted from the kinetic models: driving environmental fluctuations and interface relaxation. Note that the requirement of space- and time-continuous models to respect basic symmetries (t-, $[eqn]$ -, and h-translations plus $[eqn]$ -rotation and inversion) leaves us only with terms of the form $[eqn]$ (i.e., i-laplacians) or at most $[eqn]$ if the equilibrium condition is relaxed. Among them, $[eqn]$ and $[eqn]$ are the most relevant operators in the sense of the renormalization group. In 1982, Edwards and Wilkinson (EW) [15] proposed (for granular aggregates on arbitrary-dimensional flat substrates) the following extended Langevin equation:

$$[eqn]$$

as the continuum limit of those kinetic models. Here $[eqn]$ is a diffusion coefficient accounting for surface relaxation, and F is a particle flux. The stochastic field $[eqn]$ , with $[eqn]$ and

$$[eqn]$$

is spatially uncorrelated Gaussian white noise accounting for driving environmental fluctuations. However, something was still missing since in 1d, the EW model yielded the correct value of $[eqn]$ but not that of $[eqn]$ . Moreover for $[eqn]$ , the system remained in equilibrium.

Scarcely four years after, the missing relevant factor was identified as lateral growth, due to the fact that particles aggregate normally to the surface [1]. (The new term, gauged by parameter $[eqn]$ , is the lowest order in the expansion of the curvature.) This makes the linear stochastic partial differential Equation (1) into a nonlinear one, the Kardar–Parisi–Zhang (KPZ) equation.

$$[eqn]$$

(For readers less familiar with the KPZ equation, we note that it describes the temporal evolution of an interface $[eqn]$ under the combined effects of surface relaxation $[eqn]$ , nonlinear lateral growth $[eqn]$ , and stochastic noise $[eqn]$ . The richness of the emergent collective behavior from this simple yet non-trivial equation makes it a paradigm of non-equilibrium systems). Now for $[eqn]$ , the system remains out of equilibrium and the interface moves at a rate depending on $[eqn]$ . Once they recognized that Equation (2) can be assimilated to a Burgers equation for $[eqn]$ , they realized that $[eqn]$ must remain invariant under dynamic renormalization, which led to the aforementioned relation $[eqn]$ . Fortunately, for flat 1d substrates, the dynamic renormalization group fixed point (DRGFP) can be reached perturbatively. This allowed them to obtain analytically the aforementioned exponents $[eqn]$ , $[eqn]$ , $[eqn]$ .

Of course, $[eqn]$ becomes the most relevant operator for $[eqn]$ . Golubović and Bruinsma [16] studied the effect of this on the lower critical dimension and the exponent $[eqn]$ , and concluded that it leads to a different universality class than the KPZ one, showing that the $[eqn]$ term does not suffice (either by itself or through its interplay with the noise) to determine the universality class, but it is their interplay with interface relaxation that does.

Keeping two terms for interface relaxation (i.e., taking curvature effects into account), the growth rate equation reads

$$[eqn]$$

In Section 2, we review a variational formulation of the KPZ equation and related ones and point out some issues like the necessity of Galilean invariance in discrete versions, the Markov character of process $[eqn]$ , and the existence of steady-state distribution (before reaching saturation!) in $[eqn]$ substrates. In Section 3, as a prolegomenon to Section 4, we recall previous numerical results on the time evolution of the NEP in $[eqn]$ , 2, and 3, as well as its scaling behavior with the nonlinearity parameter $[eqn]$ . In Section 4 (always from a variational perspective), we show the intimate connection between the NEP and entropy production (both from a path-integral formulation or starting directly from the KPZ equation) and address the memory of process $[eqn]$ with respect to initial conditions, both in KPZ and GB. In Section 5, we highlight some open questions that can be approached from a variational perspective.

2. The Variational Formulation

2.1. Model A—Like Variational Formulation

Although describing a non-equilibrium process, Equation (2) originated in purely kinetic models. For some time, it was believed that its drift $[eqn]$ cannot be derived from an effective free energy [17].

Using the Hopf–Cole transformation

$$[eqn]$$

(mapping h onto a field $[eqn]$ obeying a linear equation, but with multiplicative noise), it was proved in [18] that for any d, Equation (2) can be cast as model A (namely, a relaxational system without conservation laws, in the nomenclature of [19])

$$[eqn]$$

with functional $[eqn]$ given by Equation (5), and functional $[eqn]$ (known as the motility) by

$$[eqn]$$ $$[eqn]$$

The Lyapunov property $[eqn]$ is easy to check, and the minimum of $[eqn]$ is reached by constant functions. In addition, the effective potential $[eqn]$ —the restriction of $[eqn]$ to $[eqn]$ —is minimum at $[eqn]$ for $[eqn]$ (Figure 1). Hence, $[eqn]$ is a Lyapunov functional for the KPZ problem, and the deterministic KPZ equation displays simple dynamics: the asymptotic stability of constant solutions indicates that arbitrary initial conditions approach constant profiles for long times. While the deterministic case is simple, the stochastic dynamics remain highly non-trivial, producing self-affine fractal profiles. Notably, the Lyapunov functional offers little intuition about these stochastic features.

As shown in [18], $[eqn]$ and $[eqn]$ (together with the density $[eqn]$ , defined by $[eqn]$ ) fulfill the relations postulated in [20] on the grounds of global shift invariance.

It is worth noting that several studies [21,22,23,24] have shown that in oscillator systems subjected to noise, synchronization is achieved when the coupling is strong enough. This dynamic process can be related to a roughness process and exhibit, in one dimension, scale invariance within the KPZ universality class. Since Kuramoto-type models can be obtained from energy functionals, this provides further support for the idea of the existence of such a functional for KPZ.

2.2. Discrete Version

Discretization is a necessity if the substrate dimensionality is $[eqn]$ , where the DRGFP is not reachable perturbatively and one must resort to numerical integration to find the critical exponents. In the past, a criterion to judge the quality of a discretization of Equation (2) was that it should preserve the two main symmetries of Equation (2): Galilean invariance (GI) and a fluctuation–dissipation relation peculiar of 1d (1dFDR). Note that although there is an explicit FDR for diffusion phenomena, there is in principle no such thing for growth processes. Recently, however, a kind of FDR—similar to the diffusive one—was found for growth phenomena [25,26]. In turn, Galilean invariance implies that the physics of interface growth remains the same when viewed from a reference frame moving laterally at constant speed, provided the surface height is appropriately tilted. This symmetry protects parameter $[eqn]$ from being renormalized and imposes the constraint $[eqn]$ .

Consistency in employing the Hopf–Cole transformation constraints the discrete version of $[eqn]$ to have the form $[eqn]$ . In [27], scaling was obtained to a very good approximation (Figure 2) using this discretization, which breaks GI and 1dFDR explicitly.

2.3. Memory of Process h(x,t)

From the Family–Vicsék scaling of roughening processes [29], it turns out that

$$[eqn]$$

Over the last 15 years or so, numerous works by mathematicians and physicists [2,30,31,32,33,34,35] have shown that the asymptotic (in the scaling regime) probability distribution $[eqn]$ (for any x) not only has skewness, but shows a remarkable dependence on the initial conditions. In other words, process $[eqn]$ has memory.

Equation (2) is an extended Langevin equation; thus, the joint process  $[eqn]$ obeys a Fokker–Planck equation and is Markovian. Now $[eqn]$ obeys a marginal distribution and is generically non-Markovian (unless a conservation law, like, e.g., $[eqn]$ in the EW model, forces it to be Markovian).

Besides the model A variational formulation of Equations (4) and (5), Equation (2) can also be cast (again, for any d) as [18]

$$[eqn]$$

with

$$[eqn]$$

Again, $[eqn]$ and its effective potential is minimum at $[eqn]$ for $[eqn]$ (Figure 3). But at variance with $[eqn]$ , $[eqn]$ is not bounded from below, so it is not a Lyapunov functional for $[eqn]$ . However, $[eqn]$  does provides intuition on the stochastic dynamics (for instance, it resembles an extended escape problem).

In particular, Equation (8) can be transformed into

$$[eqn]$$

which shows explicitly that $[eqn]$ depends on the whole trajectory $[eqn]$ , from 0 to t.

As in the case of reaction–diffusion systems, Equation (2) looks formally as a stochastically forced gradient flow (in the sense that the coupling terms have the same form for all sites), but it is not strictly gradient. Equation (7) is regarded as a gradient flow in the sense that its motility is $[eqn]$ .

2.4. On the Existence of a Stationary State for Process h(x,t)

The EW model of Equation (1) is an equilibrium system and as such, its joint probability distribution functional (pdf’l) $[eqn]$ is granted a stationary state $[eqn]$ . According to conservation law, $[eqn]$ , $[eqn]$ is also stationary for any x. In the KPZ model, $[eqn]$ also obeys a Fokker–Planck equation, but since its collective mode grows steadily in the scaling regime [see Equation (6)], process $[eqn]$ cannot be stationary in that regime—what is stationary is $[eqn]$ , as already argued. In principle, there is no theoretical reason why—even in 1d and with PBC—Equation (2) should attain a steady state $[eqn]$ . Yet less that $[eqn]$ be the same as that of the (equilibrium) EW model. (Moreover, recall that most results indicating strong dependence on initial conditions have been obtained for 1d substrates.)

It is just a fact that, by imposing $[eqn]$ in the functional Fokker–Planck equation for $[eqn]$

$$[eqn]$$

it admits as a stationary solution for 1d substrates and periodic boundary conditions (PBCs) the exponential $[eqn]$ of the EW Ginzburg–Landau functional

$$[eqn]$$

This is what is meant by the abovementioned 1dFDR [17].

The novelty is that Equation (10) also admits as a stationary solution (for any substrate dimensionality and independently of BC) the exponential

$$[eqn]$$

with $[eqn]$ given by Equation (8). Hence $[eqn]$ is by definition a non-equilibrium potential (NEP) [36]. Clearly, this is a formal solution since $[eqn]$ is not normalizable, but $[eqn]$ is also formal. In fact, using the short-time propagator $[eqn]$ obtained in a path-integral representation, namely

$$[eqn]$$

it becomes clear that

$$[eqn]$$

does not coincide with $[eqn]$ , indicating that $[eqn]$ is at most marginally stable.

3. KPZ NEP’s Time Behavior and Its λ-Dependence

We recall our numerical results on the time behavior of the NEP in dimensions $[eqn]$ [37]. These results (i) provide the first systematic characterization of the dimensional dependence of NEP, (ii) quantify memory effects in NEP that had not been reported, and (iii) establish explicit connections with stochastic thermodynamics.

Figure 4, Figure 5 and Figure 6 display the behavior of

$$[eqn]$$

that is, the average of $[eqn]$ over realizations of process $[eqn]$ , in substrate dimensionalities $[eqn]$ , 2, and 3, respectively. In all cases, $[eqn]$ asymptotically. The NEP also serves as a criterion for the EW–KPZ crossover, which can be estimated from the maximum of the curves for $[eqn]$ .

Figure 7 and Figure 8 display the mean $[eqn]$ -dependence over 100 samples of the NEP’s asymptotic slope and of the time of occurrence of the NEP’s maximum, respectively, for $[eqn]$ , 2, and 3. Whereas the first magnitude shows no dependence whatsoever on d, a mild $[eqn]$ -dependence is seen in the second one.

We feel that, while RG fails for $[eqn]$ , the variational approach could offer a qualitative or complementary insight for obtaining the critical exponents.

4. Novel Results

4.1. Stochastic Thermodynamics

From a path-integral approach [38,39], we proceed as follows: The probability distribution functional (pdf) of a given history (involving h at all sites) is

$$[eqn]$$

where superscript F indicates increasing time sequence, $[eqn]$ is a normalization factor, and $[eqn]$ is the forward stochastic action functional. The forward stochastic Lagrangian, given by Equation (12)

$$[eqn]$$

is still a functional of $[eqn]$ . We use calligraphic types for functionals of the configuration alone, and blackboard types for functionals of the whole history. The conditional pdf of tracing the same history back is

$$[eqn]$$

with the backward stochastic Lagrangian given by Equation (13).

$$[eqn]$$

With the entropy production rate being $[eqn]$ , one can readily obtain integral fluctuation theorems and thermodynamic uncertainty relations (TURs) [40,41].

Whereas the foregoing approach gives us a hint of the relevance of the NEP $[eqn]$ to the entropy production rate, the following approach starting directly from the KPZ equation enhances this relevance: We start expressing the KPZ equation in a semi-rigorous fashion (i.e., as one would write it in Itô interpretation if the noise was multiplicative), but explicitly taking into account that the KPZ equation is ill-posed when not regularized [2,42].

$$[eqn]$$

Here $[eqn]$ is the increment of a functional Wiener process, $[eqn]$ , $[eqn]$ , and we assume that the KPZ noise is spatially correlated over a small regularization length-scale ℓ. Although such regularization is usually omitted when introducing the KPZ equation, it appears in explicit calculations (in connection with a UV cut-off, e.g., a lattice).

For a KPZ equation with such a noise, the stochastic entropy change of the medium associated to a particular realization of the evolution of the KPZ front is given as [43]

$$[eqn]$$

where the stochastic integral in (15) is to be understood in the Stratonovich sense, i.e., the following:

$$[eqn]$$

where

$$[eqn]$$

with $[eqn]$ , and $[eqn]$ , $[eqn]$ . It is now clear that with such Stratonovich interpretation of (15), one has (compare the result with Equation (14))

$$[eqn]$$

In particular, considering formally the case of a spatially uncorrelated noise, $[eqn]$ , one has

$$[eqn]$$

which clearly shows that if the potential $[eqn]$ is understood as a Stratonovich integral, then $[eqn]$ . Also notice that only within the Stratonovich interpretation, one has

$$[eqn]$$

As a plot twist, Equation (15) illuminates the “correct” way to define the potential, in the case in which the noise is correlated, for the thermodynamic interpretation to hold.

4.2. Memory of Initial Conditions in KPZ

The dependence of the asymptotic h-statistics on initial conditions is an established fact [30,31,32,33,44]. A numerical test, undertaken with initial conditions of the form

$$[eqn]$$

shows that the NEP also displays memory effects.

Within the studied interval, the slopes show a strong dependence on the initial conditions, and in some cases, the curves also display concavity (Figure 9).

4.3. The Related Golubović–Bruinsma Model

As indicated in the Introduction, Golubović and Bruisnma [16] have analyzed an extension of the KPZ model. The equation for the the growth rate they have found is as follows:

$$[eqn]$$

As in the KPZ case, this equation can be derived from a NEP, given by

$$[eqn]$$

such that

$$[eqn]$$

Equation (23) expresses the Golubovic–Bruinsma model as a stochastic gradient flow driven by the functional derivative of the non-equilibrium potential $[eqn]$ . This formulation highlights the shared variational structure between this model and the standard KPZ equation, allowing for a unified analysis of their dynamical and statistical properties within the same theoretical framework.

From this form of the potential, we can immediately predict that, for plane initial conditions, the long time behavior will be dominated by a linear dependence on t, with a slope depending on $[eqn]$ . We can also anticipate the existence of memory effects.

In the following figures, we show the time behavior of the indicated functional as well as, changing initial conditions, we prove the existence of memory effects. In all these figures, the parameters we have used are $[eqn]$ .

From Figure 10, it is clear that the linear contributions have a similar behavior than in the EW case, while the nonlinear contribution, dominating at long time, reaches a linear dependence on t. The slope should scale with $[eqn]$ , which is an aspect to be discussed in a forthcoming work.

In Figure 11, the memory effects are clearly shown when changing the initial conditions. A detailed study will be the subject of a forthcoming work.

Considering the associated Fokker–Planck equation, and considering the stationary case we observe that, independently of the bc, the exponential

$$[eqn]$$

fullfils the stationary condition

$$[eqn]$$

with $[eqn]$ given by Equation (22). Clearly, as in the KPZ case, this is a formal solution since $[eqn]$  cannot be normalized.

An interesting aspect was studied in [45] and corresponds to the large deviations of the interface height in this model. Such a study follows [46,47] on similar aspects for the KPZ equation.

4.4. Isolating the Collective Mode

This prima facie promising approach was also worked out in [40,41]. By writing $[eqn]$ , with

$$[eqn]$$

and being $[eqn]$ , Equation (2) becomes

$$[eqn]$$

Equation (2) yields $[eqn]$ , with

$$[eqn]$$

parametrically depending on $[eqn]$ and F. Note the term $[eqn]$ , whose variance is depressed from $[eqn]$ by $[eqn]$ . Hence Equation (24) can be written in the form

$$[eqn]$$

Unfortunately, the inherent constraint $[eqn]$ hinders the usefulness of this otherwise promising approach.

5. Conclusions and Outlook

Here we highlighted, from a variational approach, several aspects of the KPZ equation as a theory of rough interface growth. We have addressed several misconceptios about KPZ—such as the alleged absence of a free-energy-type functional (the NEP), the supposed role of Galilean invariance in fixing the universality class, the Markovian character of the process, and the conjecture linking the stationary EW pdf to $[eqn]$ KPZ with periodic boundaries—and dismantled them.

In addition, we presented numerical results of the NEP (from which one can extract information on the $[eqn]$ dependence of its long-time behavior and the crossover time from EW to KPZ regime), and have shown numerically that also the NEP displays memory of initial conditions (a fact established in several theoretical works). We have also referred to our work on the stochastic thermodynamic description of KPZ based on a variational approach, allowing us to numerically obtain several known fluctuation theorems and TUR for entropy production.

We have also shown the existence of a NEP in a related model, the Golubovic–Bruinsma one. We have shown the time dependence of such a functional, as well as the existence of memory effects. We expect to extend this approach to other related systems. We also expect to exploit our functional approach to investigate large deviation functions as discussed in [45,46,47].

We are currently working on the possibility of obtaining information about the critical exponents and the existence of an upper critical dimension, using the time behavior of the NEP. While RG fails for $[eqn]$ , the variational approach offers a qualitative or complementary insight, even if exact exponents are still an open challenge.