Field-Theoretic Thermodynamic Uncertainty Relation -- General formulation exemplified with the Kardar-Parisi-Zhang equation

We introduce a field-theoretic thermodynamic uncertainty relation as an extension of the one derived so far for a Markovian dynamics on a discrete set of states and for overdamped Langevin equations. We first formulate a framework which describes quantities like current, entropy production and diffusivity in the case of a generic field theory. We will then apply this general setting to the one-dimensional Kardar-Parisi-Zhang equation, a paradigmatic example of a non-linear field-theoretic Langevin equation. In particular, we will treat the dimensionless Kardar-Parisi-Zhang equation with an effective coupling parameter measuring the strength of the non-linearity. It will be shown that the field-theoretic thermodynamic uncertainty relation holds up to second order in a perturbation expansion with respect to a small effective coupling constant.


Introduction
The thermodynamic uncertainty relation (TUR) in a non-equilibrium steady state (NESS) provides a bound on the entropy production in terms of mean and variance of an arbitrary current [1]. Specifically, in the NESS, after a time t a fluctuating integrated current X(t) has a mean X(t) = j t, and a diffusivity D = lim t→∞ (X(t) − j t) 2 /(2 t). With the entropy production rate σ the expectation of the total entropy production in the NESS is given by σ t. These quantities satisfy the universal thermodynamic uncertainty relation i.e. σ is bounded from below by j 2 /D. The TUR has been proven for a Markovian dynamics on a general network by Gingrich et al. [2,3] and further investigated for a number of different settings, both in the classical (see, e.g., [4,5,6,7,8,9,10,11,12,13,14,15]) and the quantum domain (see, e.g., [16,17,18,19,20,21,22]). It has led to a deeper understanding of systems far from equilibrium as it introduces a lower bound on the dissipation given the knowledge of the occurring fluctuations. Such a relation is of interest for the modeling and analysis of e.g. biomolecular processes, which may often be described as a Markov network (see e.g. [23,24,25]).
Of particular interest is the work by Gingrich et al. [8], where the authors extend the relation from mesoscopic Markov jump processes to overdamped Langevin equations. Here a temporal coarse-graining procedure is described, which allows the formulation of a discrete Markov jump process in terms of an overdamped Langevin equation for the mesoscopic states of the model. These authors observe that for purely dissipative dynamics the TUR is saturated. An additional spatial coarse-graining performed in [8] results in a macroscopic description, where it is found that the tightness of the resulting uncertainty relation increases with the strength of the Gaussian potential wells (see [8], fig. 9).
In this work, we present a field-theoretic equivalent to the TUR. Such a thermodynamic uncertainty relation for general field-theoretic Langevin equations may prove helpful in further understanding complex dynamics like turbulence for fluid flow or non-linear growth processes, described by the stochastic Navier-Stokes equation (e.g. [26]) or the Kardar-Parisi-Zhang equation [27], respectively. Both are prominent representatives of field-theoretic Langevin equations. For the latter, we highlight the recent progress concerning a study of the inward growth of interfaces in liquid crystal turbulence as an experimental realization. On the theory side, analytic results on the effect of aging of two-time correlation functions for the interface growth were found [28]. Furthermore we refer the reader to three review articles [29,30,31] concerning the latest developments around the Kardar-Parisi-Zhang universality class. The paper is organized as follows. In order to state a field-theoretic version of the thermodynamic uncertainty relation, we translate in section 2 the notion of current, diffusivity and entropy production known from the setting of coupled Langevin equations to their respective equivalents for general field-theoretic Langevin equations. As an illustration of the generalizations introduced in section 2, we will then study the one-dimensional Kardar-Parisi-Zhang (KPZ) equation as a paradigmatic example of such a field-theoretic Langevin equation. As the calculation of the current, diffusivity and entropy production in the NESS requires a solution to the KPZ equation, we will use spectral theory and construct an approximate solution in the weak-coupling regime of the KPZ equation in section 3. With this approximation, we will then derive in section 4 the thermodynamic uncertainty relation to quadratic order in the coupling parameter.

Thermodynamic Uncertainty Relation for a Field Theory
In this section, we will present a generalization of the thermodynamic uncertainty relation introduced in [1] to a field theory. Consider a generic field theory of the form Here Φ γ (r, t) is a scalar field or the γ-th component of a vector field (γ ∈ [1, n]; n ∈ N) with r ∈ Ω ⊂ R d , F γ [{Φ µ (r, t)}] represents a (possibly nonlinear) functional of Φ µ and η γ (r, t) denotes Gaussian noise, which is white in time, and with K(r − r ′ ) as spatial noise correlations. Prominent examples of (1) are the stochastic Navier-Stokes equation for turbulent flow (see e.g. [26]) or the Kardar-Parisi-Zhang equation for non-linear growth processes [27] to name only two. The latter will be treated in the subsequent sections within the framework established in the following. Let us begin with the introduction of some notions. A natural choice of a local fluctuating current j(r, t) is with Φ(r, t) = (Φ 1 (r, t), . . . , Φ n (r, t)) ⊤ . The local current j(r, t) is fluctuating around its mean, i.e. j(r, t) = j(r, t) + δj(r, t), with δj(r, t) denoting the fluctuations. Given that the system (1) possesses a NESS, the long-time behavior of the local current (2) can be described as with where · denotes averages with respect to the noise history. As the thermodynamic uncertainty relation in a Markovian network is formulated for some form of integrated currents, we define in analogy the projection of the local current onto an arbitrarily directed weight function g(r) The integral in (6) represents the usual L 2 -product of the two vector fields j(r, t) and g(r) with j(r, t) · g(r) = k j k (r, t)g k (r) as the scalar product between j and g. With this projected current j g (t), we associate a fluctuating 'output' Hence j g (t) = ∂ t Ψ g (t) and in the NESS The fluctuating output Ψ g (t) provides us with the means to define a measure of the precision of the system output, namely the squared variational coefficient ǫ 2 , as If the system is in its non-equilibrium steady state, we can rewrite (9) as Let us now connect the variance of the output Ψ g (t) to the Green-Kubo diffusivity given by Using (6) and (2), it is straightforward to verify that Thus, By dividing both sides of (12) by 2t and taking the limit of t → ∞ it is found in analogy to [32], that with D g from (11) and therefore Since in the NESS Ψ g (t) is stochastically independent of the initial configuration Ψ g (0), we can simplify the expression for the diffusivity in the NESS according to With the result of (14) and ǫ 2 from (9), an alternative formulation of the precision in a NESS is We proceed with expressing the total entropy production ∆s tot . The total entropy production is given by the sum of the entropy dissipated into the medium along a single trajectory, ∆s m , and the stochastic entropy, ∆s, of such a trajectory; see e.g. [33]. The medium entropy is given by, Here p[Φ(r, t)|Φ(r, t 0 )] denotes the functional probability density of the entire vector field Φ(r, t), i.e. the field configuration after some time t has elapsed since a starting-time t 0 < t, conditioned on an initial value Φ(r, t 0 ), i.e. a certain field configuration at the starting time t 0 . In contrast, p[ Φ(r, t)| Φ(r, t 0 )] is the conditioned probability density of the time reversed process, i.e. starting in the final configuration at time t 0 and ending up in the original one at time t. For the sake of simplicity, we will write in the following For the system (1), the action functional (see e.g. [33,34,35,36,37,38,39,40] and references therein) is given by where K −1 (r − r ′ ) is the inverse of the noise correlation kernel K(r − r ′ ) from (1). The two integral kernels fulfill Note, that we do not explicitly state an expression for the Jacobian ensuing from the variable transformation η(r, t) → Φ(r, t) in (18). This is justified as we will use the action functional to derive a general expression for the medium entropy where it turns out that the Jacobian does not contribute (s.f. [33,41]). Inserting (17), (18) into (16) and noticing that only the time-antisymmetric part of the action functional (18) and its time-reversed counterpart survives, leads to (see also [33,41,42]) The stochastic entropy change ∆s for the same trajectory, is given by (see also [41]) Thus, the total entropy production ∆s tot reads With (22) we may also define the rate of total entropy production σ in a NESS according to The expressions stated in (9) and (22) provide us with the necessary ingredients to formulate the field-theoretic thermodynamic uncertainty relation as with σ from (23), D g from (13) and J g from (8). The higher the precision, i.e. the smaller ǫ 2 , the more entropy ∆s tot is generated, i.e. the higher the thermodynamic cost. Or, in other words, in order to sustain a certain NESS current J g , a minimal entropy production rate σ ≥ J 2 g /D g is required.

Theoretical Background
Within this section we will lay the groundwork for the calculation of the quantities entering the TUR for the KPZ equation. The main focus thereby is on the perturbative solution of the KPZ equation in the weak-coupling regime and the discussion of issues with diverging terms due to a lack of regularity.

The KPZ Equation in Spectral Form
Consider the one-dimensional KPZ equation [27] on the interval [0, b], b > 0, with Gaussian white noise η(x, t) subject to periodic boundary conditions and, for simplicity, vanishing initial condition h(x, 0) = 0, x ∈ [0, b] (i.e. the growth process starts with a flat profile). HereL = ν∂ 2 x is a differential diffusion operator, ∆ 0 a constant noise strength, and λ the coupling constant of the non-linearity. A Fourier-expansion of the height field h(x, t) and the stochastic driving force η(x, t) reads The set of {φ k (x)} is given by and thus h k (t), η k (t) ∈ C in (26). A similar proceeding for the case of the Edwards-Wilkinson equation was used in [43,44,45,46]. Inserting (26) into (25) For the {φ k (x)} the relation φ l (x)φ m (x) = φ l+m (x)/ √ b holds and thus the double-sum in the Fourier expansion of the KPZ equation can be rewritten in convolution form setting k = l + m. This yields which implies ordinary differential equations for the Fourier-coefficients h k (t), The above ODEs (30) are readily 'solved' by the variation of constants formula, which leads for flat initial condition h k (0) ≡ 0 to k ∈ Z. Note, that the assumption of flat initial conditions is not in conflict with (21) as in the NESS, in which the relevant quantities will be evaluated, the probability density becomes stationary. With (31), a non-linear integral equation for the k-th Fourier coefficient has been derived. In subsection 3.4, the solution to (31) will be constructed by means of an expansion in a small coupling parameter λ. We close this section with the following general remarks.
(i) Equation (31) has been derived on a purely formal level. In particular, the integral dt ′ e µ k (t−t ′ ) η k (t ′ ) has to be given a meaning. In a strict mathematical formulation, this integral has to be written as which is called a stochastic convolution (see e.g. [47,48,49,50]). This has its origin in the fact that the noise η(x, t) in (25) is mathematically speaking a generalized time-derivative of a Wiener process W (x, t) (see also subsection 3.2, (35)). In this spirit, (31) with the first integral on the right hand side replaced by (32) may be called the mild form of the KPZ equation (in its spectral representation) and h(x, t) = k∈Z h k (t)φ k (x), h k (t) solution of equation (31), is then called a mild solution of the KPZ equation. In mathematical literature, proofs of existence and uniqueness of such a mild solution can be found for various assumptions on the regularity of the noise (see e.g. [49,51,52] and references therein). An assumption will be adopted (see subsection 3.2), which guarantees the existence of h(x, t) L2([0,b]) , i.e. the norm on the Hilbert space of square-integrable functions L 2 . This norm, or respectively the corresponding L 2 -product, denoted in the following by (·, ·) 0 , of h with any L 2 -function g, i.e. (h, g) 0 , will be used in subsection 4.1 and subsection 4.3 to calculate the necessary contributions to a field-theoretic thermodynamic uncertainty relation. Furthermore, with this assumption on the noise, it is shown in Appendix C for the mild solution that almost surely h( . This justifies the choice H = L 2 ([0, b]) in the following calculations.
(ii) The Fourier expansion applied above can be understood in a more general sense. For the case of periodic boundary conditions, the differential operator L possesses the eigenfunctions {φ k (x)} and corresponding eigenvalues {µ k } from (27) and (28), respectively. It is well-known that the set {φ k (x)} constitutes a complete orthonormal system in the Hilbert space L 2 (0, b) of all square-integrable functions on (0, b). Thus the Fourier-expansion performed above can also be interpreted as an expansion in the eigenfunctions of the operatorL.
(iii) With this interpretation, (31) also holds for a 'hyperdiffusive' version of the KPZ equation in which the operatorL is replaced byL p ≡ (−1) p+1 ∂ 2p x , with p ∈ N and adjusted eigenvalues {µ p k }. This may be used to introduce a higher regularity to the KPZ equation.
(iv) Besides the complex Fourier expansion in (26) [49]) and will be used in the next section. The relationship between h k (t) and h k (t) reads with h k (t) as the complex conjugate.

A Closer Look at the Noise
In the following discussion of the noise it is instructive to pretend, for the time being, that the noise is spatially colored with noise correlator K(x − x ′ ) instead of assuming directly spatially white noise. The noise η(x, t) is given by a generalized time-derivative of a Wiener process Such a Wiener process W (x, t) can be written as (e.g. [47,49]) Here {α k } ∈ R are arbitrary expansion coefficients that may be used to introduce a spatial regularization of the Wiener process, {β k (t)} ∈ R are stochastically independent standard Brownian motions and {γ k (x)} from (33). A well-known result for the two-point correlation function of two stochastically independent Brownian motions β k (t) reads [47] In the following it will be shown that the noise η defined by (35) and (36) possesses the autocorrelation results in the one assumed in (25). Furthermore, an explicit expression of the kernel K(x − x ′ ) by means of the Fourier coefficients {α k } of W (x, t) from (36) will be given. To this end, first an expression for the two-point correlation function of the Wiener process itself can be derived according to To represent the noise structure dictated by (25), the expression in (39) has to be an even, translationally invariant function in space. Thus, the following relation has to be fulfilled Then the two-point correlation function of the Wiener process is given by (27), equation (41) implies for the two-point correlation function of the Fourier coefficients W k (t) This result leads immediately to For the relation between (41) and the noise from (38), we differentiate (41) with respect to t and t ′ yielding The following identification can be made which structurally represents the standard implicit assumption that K(x − x ′ ) is translationally invariant, positive definite and even. Note, that the regularity of the noise-kernel K(|x − x ′ |) is given by the behavior of the set of {α k } for k → ∞, where {α k } are the dimensionless Fourier coefficients of the underlying Wiener process from (36) for all k. For the case of α k = 1 ∀ k ∈ Z, spatially white noise is obtained. Thus, the derivation via the Wiener process has indeed led to a translationally invariant real-valued two-point correlation function for η(x, t), given by (38), (45), which describes white in time and spatially colored Gaussian noise. In the following, we will use (45) to approximate spatially white noise to meet the required form in (25). Now the assumption mentioned in the remarks in subsection 3.1 can be made more precise. In the following it will be assumed that (see Appendix C) This assumption excludes white noise for k ∈ N, but via the introduction of a cutoff parameter Λ ∈ N, Λ ≫ 1 arbitrarily large but finite, for the range of k, white noise is accessible, i.e. for k ∈ R with Note that for the linear case, i.e. the Edwards-Wilkinson model, the authors of [44] also introduce a cutoff, albeit in a slightly different manner. Such a cutoff amounts to an orthogonal projection of the full eigenfunction expansion of (25) to a finite-dimensional subspace spanned by the eigenfunctions . Mathematically, this projection may be represented by a linear projection operator P Λ , which maps the Hilbert space (29). This mapping, however, causes a problem in the non-linear term of (29), where by mode coupling the k-th Fourier mode (−Λ ≤ k ≤ Λ) is influenced also by modes with |l| > Λ. This issue can be resolved by choosing Λ large enough, for modes with h l (t) ∼ exp[µ l t], µ l from (28), (61), |l| > Λ will be damped out rapidly so that the bias introduced by limiting l to the interval R is small. Note that the restriction to h ∈ span{φ −Λ , . . . , φ Λ } also implies the introduction of restricted summation boundaries in the convolution term in (31), namely This restriction to finitely many Fourier modes is not as harsh as it might seem, since for very large wavenumbers the dynamics of the KPZ equation is governed by the Edwards-Wilkinson equation, which, due to its equilibrium behavior, does not contribute to the thermodynamic uncertainty relation (24) (see e.g. [53,54,55,56]).
With the cutoff Λ, condition (46) is of course fulfilled for α k = 1 ∀ k ∈ R and α k = 0 ∀k / ∈ R. Inserting this choice of α k into (45) yields (49) Also, the choice of α k = 1 ∀ k ∈ R implies for the correlation function of the Fourier coefficients η k (t) from (43) To end this section, a noise operatorK describing spatial noise correlations will be introduced asK with kernel K(x − x ′ ) from (49) and its inverseK −1 given bŷ where its kernel reads Field-Theoretic Thermodynamic Uncertainty Relation 13

Dimensionless Form of the KPZ Equation
Before the KPZ equation is analyzed further, it is prudent to relate all physical quantities to suitable reference values so that the scaled quantities are dimensionless and that the equation is characterized by only one dimensionless parameter. In anticipation of the calculations below, we choose this parameter to represent a dimensionless effective coupling parameter λ eff , that replaces the coupling constant λ from (25). To this end the following characteristic scales are introduced, Here H is a characteristic scale for the height field (not to be confused with the notation for the Hilbert space), N a scale for the noise field, b is the characteristic length scale in space and T the time scale of the system. Choosing the three respective scales according to leads to the dimensionless KPZ equation on the interval x ∈ [0, 1] Here, the effective dimensionless coupling constant is given by and The effective coupling constant λ eff is found in various works concerning the KPZ-Burgers equation; see e.g. [57,58,59,37]. In the following sections we will perform all calculations for the dimensionless KPZ equation. This requires one simple adjustment in the linear differential operatorL on x s ∈ [0, 1], which is now given bŷ with eigenvalues to the orthonormal eigenfunctions Furthermore, the noise correlation function in Fourier space from (50) now reads The scaling also effects the noise operators defined in (51), (52) at the end of subsection 3.2. The scaled ones read Note that for the sake of simplicity the subscript s will be dropped in the calculations below where all quantities are understood as the scaled ones.

Expansion in a Small Coupling Constant
Returning to the nonlinear integral equation of the k-th Fourier coefficient of the heights field, h k (t) from (31), now in its dimensionless form and with the restricted spectral range given by (48) and all quantities dimensionless, an approximate solution will be constructed. Note, that the summation of the discrete convolution in (66) with h (0) Thus every h (n) k , n > 1, can be expressed in terms of h (0) m , m ∈ R, i.e. the stochastic convolution according to (32), which is known to be Gaussian. In the following calculations multipoint correlation functions have to be evaluated, which can be simplified by Wick's theorem, where a recurring term reads h (0) . It is thus helpful to determine this correlation function in general once and use this result later on. With (63) and k, l ∈ Z (and therefore also for k, l ∈ R) it follows that: Since for the auxiliary expression Π k,l the symmetries hold, it is found that

Thermodynamic Uncertainty Relation for the KPZ Equation
In this section we will show that the thermodynamic uncertainty relation from (24) holds for the KPZ equation driven by Gaussian white noise in the weakcoupling regime. In particular, the small-λ eff expansion from subsection 3.4 will be employed. To recapitulate, the two ingredients needed for the thermodynamic uncertainty relation are (i) the long time behavior of the squared variation coefficient or precision ǫ 2 of Ψ g (t) from (9); (ii) the expectation value of the total entropy production in the steady state, ∆s tot from (22).

Expectation and Variance for the Height Field
With (7) adapted to the KPZ equation, namely with g(x) as any real-valued L 2 -function fulfilling 1 0 dxg(x) = 0, i.e. g(x) possessing non-zero mean, we rewrite the variance as As is shown below, ǫ 2 can be evaluated for arbitrary time t > 0. However, the final interest is on the non-equilibrium steady state of the system. Therefore, the long-time asymptotics will be studied.

Evaluation of Expectation and Variance
In the small λ eff expansion, the expectation of the output Ψ g (t) from (74), with h(x, t) solution of the dimensionless KPZ equation (55) to (57) reads: where g k and g k are the k-th Fourier coefficient of the weight function g(x) and its complex conjugate, respectively. Here the result from (68) is used as well as the fact that odd moments of Gaussian random variables vanish identically.
Replacing h Note, that in the case of k = 0 the second term in the last line of (77) is evaluated in the limit µ k → 0, which yields t. Since the interest is on the steady state behavior, the long-time asymptotics of the two expressions in (77) above is studied. So, eq. (76) yields where g k = g −k ∀k as g(x) ∈ R. Using the explicit form of µ k from (61), the expression in (78) can be simplified according to with Λ from (47). Equivalently, the steady state current from (8) reads The first term of the variance as defined in (75) reads in the small-λ eff expansion where moments proportional to λ eff (and λ 3 eff ) vanish due to (68) and (69)  In Appendix A, we present the rather technical derivation of Subtraction of (78) squared from (82) leads to Again, with µ k from (61), the above expression in (83) can be reduced to Here H (2) Λ = Λ l=1 1/l 2 is the so-called generalized harmonic number, which converges to the Riemann zeta-function ζ(2) for Λ → ∞. Using (13), eq. (84) yields the diffusivity D g , With (84) and (79) squared, the first constituent of the thermodynamic uncertainty relation, ǫ 2 = Var[Ψ g (t)]/ Ψ g (t) 2 from (9), is given for large times by Note, since ǫ 2 ≈ 4/(λ 2 eff t), the long time asymptotics of the second term has to scale as ∆s tot ∼ λ 2 eff t for the uncertainty relation to hold. Note further, that the result for the precision of the projected output Ψ g (t) in the NESS is independent of the choice of g(x).

Alternative Formulation of the Precision
Before we continue with the calculation of the total entropy production, we would like to mention an intriguing observation. From the field-theoretic point of view, it seems natural to define the precision ǫ 2 as This is due to the fact that the height field h(x, t) is at every time instance an element of the Hilbert-space L 2 ([0, 1]) as mentioned in subsection 3.1. Hence, the difference between h(x, t) and its expectation is measured by its L 2 -norm. Also the expectation squared is in this framework given by the L 2norm squared. At a cursory glance, the definitions in (87) and (9) seem to be incompatible. However, for the case of the above calculations of ǫ 2 for the onedimensional KPZ equation, it holds up to O(λ 3 eff ) in perturbation expansion that Thus, with (88), it is obvious that in terms of the perturbation expansion both definitions of the precision, as in (9) and (87), respectively, are equivalent. Equation (88) can be verified by direct calculation along the same lines as above in this section. By studying these calculations it is found perturbatively that the height field h(x, t) is spatially homogeneous, which is reflected by h k (t)h l (t) ∼ δ k,−l (see (73)) for the correlation of its Fourier-coefficients. Further, the long-time behavior is solely determined by the largest eigenvalue of the differential diffusion operatorL = ∂ 2 x , namely by µ 0 = 0 (see e.g. (78) and (83), the essential quantities for deriving (88)). In the following, we would like to give some reasoning why the above two statements should also hold for a broad class of field-theoretic Langevin equations as in (1). For simplicity, we restrict ourselves in (1) to the case of one-dimensional scalar fields Φ(x, t) and F [Φ(x, t)] =LΦ(x, t) +N [Φ(x, t)]. HereL denotes a linear differential operator andN a non-linear (e.g. quadratic) operator.L should be selfadjoint and possess a pure point spectrum with all eigenvalues µ k ≤ 0 (e.g.L = (−1) p+1 ∂ 2p x , p ∈ N, i.e. an arbitrary diffusion operator subject to periodic boundary conditions). For this class of operatorsL there exists a complete orthonormal system of corresponding eigenfunctions {φ k } in L 2 (Ω). If it is further known, that the solution Φ(x, t) of (1) belongs at every time t to L 2 (Ω), we can calculate e.g. the second moment of the projected output Ψ g (t) according to (Ψ g (t)) 2 = ( Ω dx Φ(x, t)g(x)) 2 , where g(x) ∈ L 2 (Ω) as well. As is the case in e.g. equation (81), the second moment is determined by the Fourier-coefficients Φ k (t) of Φ(x, t) and g k of g(x), namely Like the KPZ equation, (1) is driven by spatially homogeneous Gaussian white noise η(x, t) with two-point correlations of the Fourier-coefficients η k (t) given by η k (t)η l (t) ∼ δ k,−l . Therefore, we expect the solution to (1) subject to periodic boundary conditions to be spatially homogeneous as well, at least in the steady state, which implies see e.g. [60,61]. Hence, with (90), the expression in (89) becomes Comparing (91) to Φ(x, t) 2 0 , which is given by we find in the NESS provided that the long-time behavior is dominated by the Fourier-mode with largest eigenvalue, i.e. k = 0 with µ 0 = 0. Under the same condition, the first moment of the projected output reads in the NESS and thus Similarly, Note, that g 0 and Φ 0 (t) have to be real throughout the argument (which is indeed the case for expansions with respect to the eigenfunctions of the general diffusion operatorsL from above). Hence, under the assumption that the prior mentioned requirements are met, which, of course, would have to be checked for every individual system (as was done in this section for the KPZ equation), the asymptotic equivalence in (93) and (97) validates the statement in (88) (and therefore, in the NESS, also (87)) for a whole class of one-dimensional scalar SPDEs from (1).

Total Entropy Production for the KPZ Equation
The total entropy production for the KPZ equation is obtained by inserting This stationary solution is the same as the one for the linear case, namely for the Edwards-Wilkinson model. Note that in (100) we denote by · 2 0 the standard L 2 -norm.
Stationary Total Entropy Production With (22), the total entropy production in the NESS for the KPZ equation reads Using ḣ , ∂ 2 x h 0 , and the initial condition h(x, 0) = 0, the first term in (101) vanishes and thus (102) For Gaussian white noise, the expectation value of (102) is given by For a derivation of this result see Appendix B. Note that (103) and its derivation remains true for h ∈ span{φ −Λ , . . . , φ Λ }. More generally, the expectation of the total entropy production may also be written as withK −1 from (65).
Evaluating the Expectation of the Stationary Total Entropy Production Above, an expression for the stationary total entropy production ∆s tot and its expectation value were derived (see eq. (103)). Inserting the Fourier representation from (26) and (62) with R k from (48). As (105) above is already of order λ 2 eff , it suffices to expand the Fourier coefficients h i (t ′ ) to zeroth order, which yields with h (0) i (t ′ ) given by (68). Via a Wick contraction and using (73), the fourpoint correlation function in (106) reads Inserting (107) into (106) leads to the following form of the total entropy production in the NESS, Note that the long time behavior of ∆s tot is indeed of the form required, i.e. ∆s tot ∼ λ 2 eff t (see remark after (86)), for the uncertainty relation to hold. With µ k from (61), the expression for the total entropy production from (108) reads Thus, with (23) and (109), the total entropy production rate reads With (86) and (109), or, equivalently, (80), (85) and (110), the constituents of the thermodynamic uncertainty relation are known. Hence, the product entering the TUR from (24) for the KPZ equation reads Here, we deliberately refrain from writing ∆s tot ǫ 2 = 5 − 1/Λ as this would somewhat mask the physics causing this result. This point will be discussed further in the following.

Edwards-Wilkinson Model for a Constant Driving Force
To give an interpretation of the two terms in (108) and consequently in (111), we believe it instructive to briefly calculate the precision and total entropy production for the case of the one-dimensional Edwards-Wilkinson model modified by an additional constant non-random driving 'force' v 0 and subject to periodic boundary conditions. To be specific, we consider already in dimensionless form and with space-time white noise η. We denote (112) in the sequel with FEW for 'forced Edwards-Wilkinson equation'. Following the same procedure as described in section 3, we find the following integral equation for the k-th Fourier coefficient of the height field in FEW, where again a flat initial configuration was assumed and µ k = −4π 2 k 2 as above. With (113), we get immediately in the NESS and thus Ψ g (t) 2 = g 2 0 v 2 0 t 2 as well as Thus, As already discussed above in section 3, the Fokker-Planck equation corresponding to (112) has the stationary solution p s [h] = exp − dx (∂ x h) 2 and thus, with (22) and (113), the total entropy production reads in the NESS With (115) and (116), the TUR product for (112) is given by i.e. the thermodynamic uncertainty relation is indeed saturated for the Edwards-Wilkinson equation subject to a constant driving 'force' v 0 . For the sake of completeness we state the expressions for the current, diffusivity and rate of entropy production in the non-equilibrium steady state, namely With the calculations for FEW, we can now give an interpretation of the two terms in (109) and (111). The first term in squared brackets in (109) originates from the first term of (108), where the latter represents the action of all higher-order Fourier modes on the mode k = 0 (see (107)). To illustrate this point further, observe that, in the NESS, we get according to (76) to (80) for the current: and from the calculation above we see that it contains only the impact of Fourier modes l = 0 on the mode k = 0, which belongs to the constant eigenfunction φ 0 (x) = 1. In other words, the modes l = 0 act like a constant external excitation, just in the same manner as v 0 acts for FEW in (114). Comparing (119) to (114), we may set and get J g = g 0 v 0 in both cases. Following now the calculations for FEW, we would expect from (116) which is in fact exactly the first term in the squared brackets from (108) and (109), respectively. Since with (120) also the expression for ǫ 2 from (115) coincides with the first summand on the r.h.s. of (86), it is clear that both cases result in the saturated TUR. This explains the value 2 on the r.h.s. of (111).
Turning to the second term of (109), we see that it stems from the second term in (108). In contrast to the first term in (108), the second one does not only measure the effect of the modes on the k = 0 mode but also on all other modes k = 0. It further features interactions of the k and l modes among each other via mode coupling. Hence, the mode coupling seems responsible for the larger constant on the right hand side of (111), since by neglecting the mode coupling term in (109), the thermodynamic uncertainty relation was saturated also for the KPZ equation up to O(λ 2 eff ). To conclude this brief discussion, we give the respective relations of the KPZ current (80), diffusivity (85) and total entropy production rate (110) to FEW, namely with J FEW g , D FEW g and σ FEW from (118). We see that the additional mode coupling term in KPZ leads to corrections in D KPZ and σ KPZ of at least second order in λ eff . For the case of λ eff → 0 the KPZ equation becomes the standard Edwards-Wilkinson equation (EW), namely ∂ t h(x, t) = ∂ 2 x h(x, t) + η(x, t), which possesses a genuine equilibrium steady state. Therefore, for the standard EW we have J EW g = 0, σ EW = 0 and D EW g = g 2 0 /2. From (122) it follows that for λ eff → 0, (J g , σ, D g ) KPZ → (J g , σ, D g ) FEW and from (118), (120) that (J g , σ, D g ) FEW → (J g , σ, D g ) EW = (0, 0, g 2 0 /2). Hence, the non-zero expressions for J KPZ g and σ KPZ result solely from the KPZ non-linearity. The impact of the latter on the k = 0 Fourier mode (i.e. the spatially constant mode) results in contributions to J KPZ g and σ KPZ that can be modeled exactly by FEW, the Edwards-Wilkinson equation driven by a constant force v 0 from (112).

Conclusion
We have introduced an analog of the TUR [1,2] in a general field-theoretic setting (see (24)) and shown its validity for the Kardar-Parisi-Zhang equation up to second order of perturbation. To ensure convergence of the quantities entering the thermodynamic uncertainty relation for the case of Gaussian spacetime white noise, we had to introduce an arbitrarily large but finite cutoff Λ of the corresponding Fourier spectrum. While this cutoff solves the issue of divergences, it naturally leads to subtleties in treating the non-linearity, as its Fourier spectrum is affected also by modes that are beyond the considered spectral range. In order to minimize the resulting bias, the cutoff has to be chosen large enough such as to guarantee the dominance of the diffusive term over the non-linear term. To circumvent the introduction of a cutoff to ensure convergence, a possible solution may be to induce a higher regularity by treating spatially colored noise instead of Gaussian white noise and/or choosing a higher order diffusion operatorL (see e.g. [65,66]) . This is currently under investigation.
As is obvious from (111), the field-theoretic version of the TUR for the KPZ equation displays a greater constant than the one in [1]. This is due to the mode-coupling of the fields as a consequence of the KPZ non-linearity. To illustrate this point, we also treated the Edwards-Wilkinson equation in section 4.3, driven out of equilibrium by a constant velocity v 0 , see (112). By identifying v 0 with the influence of higher-order Fourier modes on the mode k = 0, we may interpret the first term in (108) as the contribution from the forced Edwards-Wilkinson equation, for which the TUR with constant equal to 2 is saturated (see (117)), an observation which is in accordance with findings in [8] for finite dimensional driven diffusive systems. The second term in (108) is the contribution to the entropy production made up by the interaction between Fourier modes of arbitrary order, which is due to the mode coupling generated by the KPZ non-linearity. It is this additional entropy production that weakens the dissipation bound in the TUR. Note, that also the first term in (108) is due to the mode-coupling, however is special in thus far that it measures only the impact of the other modes on the zeroth k-mode and does not include a response of the mode k = 0.
Regarding future research, an intriguing topic is the question as to whether the findings in [8] concerning conditions for the saturation of the dissipation bound in the TUR for an overdamped two-dimensional Langevin equation can be recovered in the present field-theoretic setting. Furthermore, it would be of great interest to employ the developed framework to other field-theoretic Langevin equations in order to observe the resulting dissipation bounds in the corresponding TURs. Of special interest in this context is the stochastic Burgers equation, especially, if excited by a noise term suitable for generating genuine turbulent response (see [67]). A comparison of the predictions made in the present paper to numerical simulations of the KPZ equation seems to be another intriguing task. Besides numerical calculations, it would also be of great interest to test our predictions via experimental realizations of KPZ interfaces. Lastly, the formulation of a genuine non-perturbative, analytic formalism would also be of utmost interest.
A Evaluation of (81) Using (73), the first term in (81) reads Note that the case of k = 0 is treated like in (77). The second term in (81) is given by where we used Wick's-theorem, (73) and (69). Note that the two Kronecker-deltas in the last term of (124) can also be written as δ 0,k δ 0,l δ k,−l , such that the whole expression is multiplied by δ k,−l . Again with Wick's-theorem, (73) and (70) we can calculate the third and forth term of (81) accordingly and find As can be seen from (123) to (125), all four terms in (81) contain a δ k,−l and thus (81) reduces to The first term of (126) is readily evaluated with (123) as Oliver Niggemann, Udo Seifert The second term of (126) reads with (124): dr Πn,n(r, r) Hence, with h Here the choice of the minimum of (t ′ ∧ r) is arbitrary, since for (t ′ ∧ r) = r the other case is obtained by simply interchanging r ↔ t ′ under the integral and vice versa; thus the results for both choices are equivalent. In the following (t ′ ∧ r) = r is chosen. Hence, the integral expression in (129) can be evaluated as Thus, with (129) where we changed m → l. To save computational effort, rewrite the last two terms of (126) in the following way Hence, it suffices to calculate one of the two expectation values. With (125) we see that where we substituted m → k − m and used the symmetry of Π k,l (t, t ′ ) from (72). Note that for k = 0, the above expression in (133) vanishes. Thus in the following calculations k = 0 is assumed. In this setting, (133) reads with (73) − 16π 4 e µ k t t 0 dt ′ e −µ k t ′ l∈R k \{0,k} kl 2 (k − l)e µ k−l t ′ t ′ 0 dr e −µ k−l r Π k,k (t, r)Π l,l (t ′ , r) = −16π 4 l∈R k \{0,k} kl 2 (k − l) 4µ k µ l 1 2µ k (µ l + µ k + µ k−l ) for k = 0 and t ≫ 1, where we changed summation index m → l. Thus, with the results from (127) (135)

B Expectation of the Total Entropy Production
The Fokker-Planck equation for the KPZ equation from (55)  Here it is understood that · now denotes the expectation value with regard to the stationary distribution p s [h]. The total stationary entropy production ∆stot is given by (see (102)) Hence its expectation value reads which is evaluated with the aid of (144): (147)

C Regularity Results for the one-dimensional KPZ Equation
Dealing with the one-dimensional KPZ equation allows us to make use of the equivalence to the stochastic Burgers equation and adapt the regularity results for the latter from [49,52,68,69,51]. In subsection 3.1 and subsection 3.2, we found that our operatorsL andK share the same set of eigenfunctions, which simplifies the results obtained by the authors of [52,51] to the following. Under the assumption that and therefore we get the following result for the 1d-KPZ equation: Here As spatial white noise is excluded by (151), we introduced a finite cutoff Λ of the Fourier spectrum instead of using our approximation to the KPZ solution as a spectral Galerkin scheme and letting Λ → ∞. However, for the KPZ equation driven by spatially colored noise satisfying (151) or even an adapted version of (151) to a higher order diffusion operator as defined in subsection 3.1 (see e.g. [65,66]), in future work we want to derive a TUR taking the full Fourier spectrum into account. We would like to conclude with the following remark. Since a couple of years, there exists a complete existence and regularity theory for the KPZ equation driven by space-time white noise introduced by Hairer [70] (see also [71,72] and for further reading on the so-called regularity structures developed in [70] see [73]). In [70] it is shown that the solutions of the KPZ equation with mollified noise converge after a suitable renormalization to the solution of the renormalized KPZ equation with space-time white noise, when removing the regularization. It is due to this renormalization procedure (where a divergent quantity needs to be subtracted) and the poor regularity of the solution, that at present it is not obvious to us how the method developed in [70] can be of use for constructing a TUR.