Stationary states in infinite volume with non zero current

TL;DR

This paper demonstrates the existence of stationary measures with non-zero current in infinite-volume Ginzburg-Landau stochastic models under specific geometries in three or more dimensions, without external forces.

Contribution

It proves the existence of non-zero current stationary states in infinite-volume stochastic models without external forces, under certain geometric conditions.

Findings

Existence of stationary measures with non-zero current in infinite domains.

Non-zero current states occur without external forces in specific geometries.

Results are valid in three or more dimensions.

Abstract

We study the Ginzburg-Landau stochastic models in infinite domains with some special geometry and prove that without the help of external forces there are stationary measures with non zero current in three or more dimensions.

Equations

$${\mathbb{Z}}_{n,+}^d:=\{x:=(x_1,\dots ,x_d)\in {\mathbb{Z}}^d:x_1\ge n\}.$$

$${\mathbb{Z}}_{-n,-}^d:=\{x:=(x_1,\dots ,x_d)\in {\mathbb{Z}}^d:x_1\le -n\}.$$

$$C_n:=\{x\in {\mathbb{Z}}^d:|x_i|\le n,\forall i=1,\dots d\}.$$

$${\mathcal{X}}_n:={\mathbb{Z}}_{n,+}^d\cup C_n\cup {\mathbb{Z}}_{-n,-}^d.$$

$${\Lambda }_{n,N}={\mathcal{X}}_n\cap [-N,N{]}^d\text{and}{S}_{n,N}={\Lambda }_{n,N+1}\setminus {\Lambda }_{n,N}.$$

$$G\psi (x)=\sum _{\begin{array}{c}y\in {\mathcal{X}}_n\\ y\sim x\end{array}}[\psi (y)-\psi (x)].$$

$$I_{\xi }(\psi )=\sum _{x\in {\Sigma }_{\xi }}{j}_{x-e_1\to x}(\psi )$$

$${\lambda }_x={\lambda }^{-}\cdot p_x^{-}+{\lambda }^{+}\cdot p_x^{+},$$

$$p_x^{\pm }={\mathbb{P}}_x(X(\cdot )\in {\mathbb{Z}}_{\pm n,\pm }^d\text{definitively}).$$

$$G^{N,\Delta}\psi(x)=\left\{\begin{array}[]{ll}\sum_{\begin{subarray}{c}y\in\Lambda_{n,N}\\ {y\sim x}\end{subarray}}[\psi(y)-\psi(x)]+\sum_{\begin{subarray}{c}y\in\Delta\\ {y\sim x}\end{subarray}}[\psi(y)-\psi(x)]&\text{if }x\in{\Lambda_{n,N}},\\ &\\ 0&\text{if }x\in{\Delta}.\end{array}\right.$$

$${\lambda }_x^{(N,\Delta ,\sigma )}=\sum _{y\in \Delta }{\sigma }_y{\mathbb{P}}_x(X^{N,\Delta }(\tau )=y)$$

$$\Delta ={\Delta }_{+}\cup {\Delta }_{-}\text{where}{\Delta }_{\pm }=\{y\in S_{n,N}:y\mp e_1\in {\Lambda }_{n,N}\}$$

$$\Delta =S_{n,N}$$

$$H(\varphi )=\sum _{A\in \mathcal{A}}{V}_A({\varphi }_A),$$

$$H(\varphi )=\sum _x({\varphi }_x^2-1{)}^2+\frac{1}{2}\sum _{x\sim y}({\varphi }_x-{\varphi }_y{)}^2.$$

$$H(\varphi )=\frac{1}{2}\sum _{x\in {\mathcal{X}}_n}{\varphi }_x^2.$$

$$H_{\Lambda }({\varphi }_{\Lambda })=\sum _{A\in \mathcal{A}:A\subset \Lambda }{V}_A({\varphi }_A),$$

$$H_{\Lambda }({\varphi }_{\Lambda })=H_{\Lambda }^0({\varphi }_{\Lambda })+H_{\Lambda }^{\prime }({\varphi }_{\Lambda }),H_{\Lambda }^0({\varphi }_{\Lambda })=a\sum _{x\in \Lambda }{\varphi }_x^2,H_{\Lambda }^{\prime }({\varphi }_{\Lambda })\ge -b|\Lambda |$$

$$|V_{A}(\phi_{A})|+\Big{|}\frac{\partial}{\partial\phi_{x}}V_{A}(\phi_{A})\Big{|}+\Big{|}\frac{\partial^{2}}{\partial\phi_{x}^{2}}V_{A}(\phi_{A})\Big{|}\leq c\sum_{x\in A}\phi_{x}^{2k}$$

$$H_{\Lambda }({\varphi }_{\Lambda }|{\varphi }_{{\Lambda }^c})=H_{\Lambda }({\varphi }_{\Lambda })+\sum _{\begin{array}{c}A:A\cap \Lambda \ne \varnothing \\ A\cap {\Lambda }^c\ne \varnothing \end{array}}{V}_A(\varphi ),\varphi =({\varphi }_{\Lambda },{\varphi }_{{\Lambda }^c})$$

$$H_{\Lambda ,R}({\varphi }_{\Lambda }|{\varphi }_{{\Lambda }^c})=H_{\Lambda }^0({\varphi }_{\Lambda })+g_R(\Vert {\varphi }_{\Lambda }{\Vert }_2^2)H_{\Lambda }^{\prime }({\varphi }_{\Lambda }|{\varphi }_{{\Lambda }^c})$$

$$\Vert {\varphi }_{\Lambda }{\Vert }_2^2=\sum _{x\in \Lambda }{\varphi }_x^2$$

$${\varphi }_x(t)$$

$$w_{x,y}(t)=B_{\{x,y\}}(t)\mathrm{if}x<y,w_{x,y}(t)=-B_{\{y,x\}}(t)\mathrm{if}x>y$$

$$J_{x\to y}=\frac{\partial H}{\partial {\varphi }_x}-\frac{\partial H}{\partial {\varphi }_y}.$$

$${\varphi }_x^{(N,\Delta ,\sigma )}(t)$$

$$L^{n,N,\Delta ,\sigma }=\sum _{\begin{array}{c}x,y\in {\Lambda }_{n,N}\\ \{x,y\}\end{array}}{L}_{x,y}+\sum _{x\in {\Lambda }_{n,N}}\sum _{\begin{array}{c}y\in \Delta \\ y\sim x\end{array}}{\bar{L}}_{x,y},$$

$$L_{x,y}=-\left(\frac{\partial H}{\partial {\varphi }_x}-\frac{\partial H}{\partial {\varphi }_y}\right)\left(\frac{\partial }{\partial {\varphi }_x}-\frac{\partial }{\partial {\varphi }_y}\right)+\frac{1}{\beta }{\left(\frac{\partial }{\partial {\varphi }_x}-\frac{\partial }{\partial {\varphi }_y}\right)}^2$$

$$\bar{L}_{x,y}=\Big{\{}-\left(\frac{\partial H}{\partial\phi_{x}}-{\sigma}_{y}\right)\frac{\partial}{\partial\phi_{x}}+\frac{1}{\beta}\,\frac{\partial^{2}}{\partial\phi^{2}_{x}}\Big{\}}.$$

$${\mu }_{n,N,{\lambda }^{\ast }}(d\varphi |\bar{\varphi })=\frac{1}{{\mathcal{Z}}_{n,N,{\lambda }^{\ast }}(\bar{\varphi })}\cdot e^{-\beta [H_{\Lambda }(\varphi |\bar{\varphi })-{\sum }_{x\in {\Lambda }_{n,N}}{\lambda }_x^{\ast }{\varphi }_x]}d\varphi ,$$

Full text

1. Introduction

Equilibrium statistical mechanics is based on the paradigm of the Boltzmann-Gibbs distribution. This extremely powerful paradigm describes equilibrium thermodynamics and applies to a large class of systems, including phase transitions. By contrast, it does not exist a general and system-independent approach to non-equilibrium statistical mechanics, where instead dynamics plays a key role. The most natural way to create a non-equilibrium state is by putting an extended system in contact with two heat or mass reservoirs at different temperatures or chemical potentials. One could think of a $d$-dimensional box $[-N,N]^{d}$ which identifies the volume of the system and the two reservoirs are attached to the opposite faces along, say, the $x$-direction (for simplicity periodic boundary conditions are chosen in the other directions). Due to the reservoirs, the state has a non-zero current in the $x$-direction. This defines the setting of boundary-driven systems and the stationary measure of those systems is then called a non-equilibrium steady state. Usually one requires that such a state satisfies the macroscopic laws of transport, such as the Fourier’s law, by which the heat current is proportional to the gradient of the temperature, or the Fick’s law, implying proportionality between the mass current and the gradient of the mass density. As a consequence, the current in a large system scales as the inverse of the system length $N$. In particular, an infinite system ($N\to\infty$) has zero current.

The main question in this paper is about the opposite, namely the existence of stationary states of infinite systems having a non zero current. This seems paradoxical because intuition says that some external forces are needed to sustain the current which otherwise would die out. However the fact that the system has “a special geometry” does the trick, as we shall see. Ruelle [18] was the first to give an example of all that by considering a quantum model describing two infinite systems which interact with each other via a third finite system. He proved that indeed, in this setting, there are stationary states with non zero current.

Later on Gallavotti and Presutti [10, 11, 12] studied a similar geometry, namely a finite system in interaction with several distinct infinite systems. The dynamics in [10, 11, 12] is given by the classical Newton equations with Gaussian thermostatic forces added. The focus was however on the existence of dynamics in the infinite-volume and the equivalence between Gaussian thermostats and infinite reservoirs.

We will consider here the analogue of the Ruelle model in stochastic systems, the so called Ginzburg-Landau models. These are lattice systems with unbounded (real valued) spins $\phi_{x}$ called “charges”. The dynamics is stochastic but it conserves the total charge. It is therefore a continuous version of the well-known Kawasaki dynamics in the Ising model. As mentioned, the spatial geometry has an essential role. The crux of our argument is that in the geometrical set up that we consider there may exist non-constant bounded harmonic functions. We will prove that in such a case there are indeed, in $d\geq 3$ dimensions, infinite-volume stationary states with non zero current.

For technical reasons we will prove the statement for super-stable Hamiltonians with non negative, finite range interactions, the class is quite general to include cases where phase transitions are present. We use such assumptions to prove the existence of the infinite volume dynamics, we believe that they could be relaxed but this is not in the spirit of our paper.

In the case of general Hamiltonians we miss the existence of the infinite-volume dynamics but we can prove that the Fick’s law is violated, namely putting the system in contact with two reservoirs which fix the chemical potentials at the right and left faces (as described in the beginning of the introduction) we observe a current which does not decay when the size of the system diverges. See however the remarks after Theorem 3.4.

In the next section we describe the model, in Section 3 we state the main results which are then proved in the successive sections.

2. The model

2.1. The geometrical setup

We consider an infinite system arising from two semi-infinite volumes that are put in contact by means of a channel. For $n\in\mathbb{N}$, we define the $d$-dimensional semi-infinite lattice $\mathbb{Z}^{d}_{n,+}$ as the set of all points to the right of the hyperplane $x_{1}=n$

[TABLE]

Similarly we define the semi-infinite lattice $\mathbb{Z}^{d}_{n,-}$ as the set including all points to the left of the hyperplane $x_{1}=-n$

[TABLE]

Finally the channel $C_{n}$ is defined as the centered squared box of side $2n+1$ connecting the two semi-infinite lattices

[TABLE]

The infinite-volume domain is then obtained as the union

[TABLE]

Often we shall derive results about the infinite volume by first considering a finite volume of linear size $N$ and then studying the limit $N\to\infty$. Thus for all integers $N>n$ we define

[TABLE]

We will use the notation $x\sim y$ to denote nearest neighbor sites in ${\mathcal{X}}_{n}$ and $\{x,y\}$ for the un-oriented bond joining them.

2.2. Harmonic functions

We continue by identifying harmonic functions for our special geometry. Let $\{X(t),\;t\geq 0\}$ be the simple symmetric continuos-time random walk on ${\mathcal{X}}_{n}$ which jumps at rate $1$ to any of its nearest neighbor sites. We denote by $\mathbb{P}_{x}$ the law of this process started from $X(0)=x$. The process is defined by the generator working on functions $\psi:{\mathcal{X}}_{n}\to\mathbb{R}$ as

[TABLE]

We can interpret (2.6) as a conservation law because $\sum_{x\in{\mathcal{X}}_{n}}G\psi(x)=0$ and then $j_{x\to y}(\psi)=\psi(x)-\psi(y)$ can be interpreted as a “current”. When studying the Ginzburg-Landau model we will also have currents and the main point of our analysis will be that there are stationary measures whose average current is equal to $j_{x\to y}(\psi)$ with $\psi$ an harmonic function.

A function $\psi:{\mathcal{X}}_{n}\to\mathbb{R}$ is said to be harmonic if $G\psi(x)=0$ for all $x\in{\mathcal{X}}_{n}$. Harmonic functions are stationary for the evolution defined by (2.6). When studying Fick’s law we will be interested in currents through a section of the channel. Thus, for $|\xi|\leq n$, we consider the total flux $I_{\xi}$ through a section $\Sigma_{\xi}$ in the channel perpendicular to the $x_{1}$-axis, i.e. $\Sigma_{\xi}=\{x\in C_{n}\;:\;x\cdot e_{1}=\xi\}$ where $e_{1}$ denotes the unit vector along the $x_{1}$ axis. We thus define

[TABLE]

The crucial feature of our geometrical setup is that in dimension $d\geq 3$ there are non-constant harmonic functions. We shall say that the random walk $X(\cdot)$ is definitively in a set $A$ if there exists a finite $T>0$ such that for all $t\geq T$ one has $X(t)\in A$.

Definition 2.1** (The harmonic function $\lambda$).**

We fix ${\lambda}^{-},{\lambda}^{+}\in\mathbb{R}$ with $\lambda^{-}<\lambda^{+}$ and define a function $\lambda:{\mathcal{X}}_{n}\to\mathbb{R}$ as

[TABLE]

with

[TABLE]

The following proposition is proved in Appendix A:

Proposition 2.2**.**

The function ${{\lambda}}$ in Definition 2.1 satisfies the following properties.

1. (1)

It is a bounded harmonic function of the process $\{X(t),t\geq 0\}$ with generator $G$. 2. (2)

If the spatial dimension $d\geq 3$ then $p^{+}_{x}+p^{-}_{x}=1$ and ${\lambda}$ is a non-constant function. 3. (3)

The flux $I_{\xi}({\lambda})$ associated to $\lambda$ has the same value for any $|\xi|<n$ and $\frac{I_{\xi}({\lambda})}{n^{d-1}}\leq\frac{c}{n}$ for some $c>0$.

We will also consider harmonic functions in a finite volume $\Lambda_{n,N}\cup\Delta$ with $\Delta\subseteq S_{n,N}$. To this aim we introduce the process $\{X^{N,\Delta}(t),t\geq 0\}$ with generator

[TABLE]

The process $\{X^{N,\Delta}(t),t\geq 0\}$, taking values in $\Lambda_{n,N}\cup\Delta$, is a continuos time random walk that jumps at rate 1 to its nearest neighbors in $\Lambda_{n,N}\cup\Delta$ and is absorbed when it reaches $\Delta$. We call $\tau$ such absorption time.

Definition 2.3** (The harmonic function $\lambda^{(N,\Delta,{\sigma})}$ with boundary condition ${\sigma}$ on $\Delta$).**

We fix $\Delta\subset S_{n,N}$ (see (2.5)) and $\sigma:\Delta\to\mathbb{R}$ and define a function $\lambda^{(N,\Delta,{\sigma})}:\Lambda_{n,N}\cup\Delta\to\mathbb{R}$ as

[TABLE]

*Notice that $\lambda^{(N,\Delta,{\sigma})}_{x}=\sigma_{x}$ for $x\in\Delta$. *

While several results of our paper hold true for a general boundary condition ${\sigma}$ on arbitrary set $\Delta\subseteq S_{n,N}$, two particular cases will be of special interest and are described hereafter.

Hypothesis 2.4** (Special settings).**

(a) Fick’s law.* In this case*

[TABLE]

*and $\sigma_{x}=\lambda^{\pm}$ for $x\in\Delta_{\pm}$. *

(b) The full setting.* In this case*

[TABLE]

and $\sigma_{x}=\lambda_{x}$ for $x\in\Delta$, where ${{\lambda}}$ is the harmonic function in Definition 2.1.

Remark 2.5**.**

Hypothesis 2.4(a) is the natural set-up for the Fick’s law, as discussed in the Introduction. Under Hypothesis 2.4(b) we have that ${\lambda}_{x}^{N,\Delta,\sigma}=\lambda_{x}$ with $x\in\Lambda_{n,N}$ for any integer $N$, see item (3) in the proposition below. This will be used to study the infinite volume dynamics via partial dynamics, that will be defined in Section 2.5.

The following proposition is proved in Appendix A:

Proposition 2.6**.**

The function $\lambda^{(N,\Delta,{\sigma})}$ in Definition 2.3 satisfies the following.

1. (1)

It is an harmonic function with boundary condition $\sigma$ on the set $\Delta$ for the process $\{X^{N,\Delta}(t),t\geq 0\}$ with generator $G^{N,\Delta}$. 2. (2)

Under Hypothesis 2.4(a) we have $\lim_{N\to\infty}\lambda^{(N,\Delta,{\sigma})}_{x}=\lambda_{x}$. 3. (3)

Under Hypothesis 2.4(b) and for any $N\in\mathbb{N}$ we have $\lambda^{(N,\Delta,{\sigma})}_{x}={\lambda}_{x}$ for $x\in\Lambda_{n,N}$.

Remark 2.7**.**

Item (2) of Proposition 2.2 and item (2) in Proposition 2.6 show that in the context of Hypothesis 2.4(a) the current $j_{x\to y}({\lambda}^{N,\Delta,\sigma}):={\lambda}_{x}^{(N,\Delta,\sigma)}-{\lambda}_{y}^{(N,\Delta,\sigma)}$ is not identically zero in the limit $N\to\infty$.

2.3. Hamiltonian

As customary in the theory of lattice systems the energy is given in terms of its potential, thus the formal Hamiltonian is

[TABLE]

where $\mathcal{A}$ is the set of all finite subsets of the lattice ${\mathcal{X}}_{n}$, $\phi_{A}=\{\phi_{x}\}_{x\in A}$ and $V_{A}(\phi_{A})$ are $C^{\infty}$ functions. We may write $V_{A}(\phi)$ for $V_{A}(\phi_{A})$, $\phi_{A}$ in such a case is the restriction of $\phi$ to $A$. To study the infinite volume limit we will restrict to the following case:

Definition 2.8** (Positive interactions.).**

By this we mean Hamiltonians which satisfy the following four conditions.

• •

$V_{A}=0$* if the cardinality $|A|$ of $A$ is $\geq 3$, moreover there is $R$ so that $V_{\{x,y\}}=0$ if $|x-y|>R$.*

• •

$V_{A}=V_{B}$* if $B$ is a translate of $A$.*

• •

$V_{\{x,y\}}(\phi)\geq 0$, $V_{\{x\}}(\phi)\geq a\phi_{x}^{2}-b$, $a>0$.

• •

$V_{\{x,y\}}(\phi)\leq c(\phi_{x}^{2}+\phi_{y}^{2})$**

Remark 2.9**.**

In the first condition we restrict to one and two-body interactions with finite range; in the second one we suppose that the interaction is translational invariant; the third one is special. To understand the origin of the third condition it is convenient to consider the typical two-body interaction, that has the form $-c_{x,y}\phi_{x}\phi_{y}$. In the ferromagnetic case $c_{x,y}>0$ so that we can rewrite it as $\frac{1}{2}c_{x,y}(\phi_{x}-\phi_{y})^{2}-\frac{1}{2}c_{x,y}(\phi_{x}^{2}+\phi_{y}^{2})$. This means that the one body potential at $x$ has an extra term $-\frac{1}{2}\sum_{y}c_{x,y}\phi_{x}^{2}$, the assumption is then that, despite this additional term, the one-body potential is $\geq a\phi_{x}^{2}-b$, $a>0$. Thus the third condition may be seen as a strengthening of the usual super-stability condition for ferromagnetic interactions. The fourth condition is clearly satisfied in the usual case where the two body interaction has the form $-c_{x,y}\phi_{x}\phi_{y}$.

Remark 2.10**.**

The stronger super-stability condition is satisfied in the case of quadratic, ferromagnetic two-body interactions and when the one-body potential grows as $c\phi_{x}^{4}$, $c>0$. A particular case is the Hamiltonian

[TABLE]

which has a phase transition at small temperatures in $\mathbb{Z}^{d},d\geq 2$, as proved by Dinaburg and Sinai [5]. Indeed the one-body potential has a double-well shape with two minima at $\pm 1$ and thus forces the charges to be close to $\pm 1$; the quadratic interaction term forces the charges to be equal. As a consequence, at low temperatures the Gibbs measure concentrates on configurations where the charges are mostly close to $+1$ (or to $-1$).

Remark 2.11**.**

Another Hamiltonian that satisfies the four conditions stated above is the quadratic Hamiltonian

[TABLE]

Here the potentials are only one-body, the interactions are absent. It is however interesting because it has almost explicit solutions obtained by using duality.

We use the assumption of positive interactions to study the infinite-volume dynamics. In finite volumes we can be much more general. In the whole sequel ${\Lambda}$ will denote a bounded set in ${\mathcal{X}}_{n}$ and

[TABLE]

the energy of $\phi_{\Lambda}$ in ${\Lambda}$.

Definition 2.12** (“General” interactions).**

• •

There are integers $K$ and $R$ so that $V_{A}=0$ if the cardinality $|A|$ of $A$ is $\geq K$ or if the diameter of $A$ is $>R$.

• •

$V_{A}=V_{B}$* if $B$ is a translate of $A$.*

• •

There are $a>0$ and $b\geq 0$ so that, for any bounded ${\Lambda}\in{\mathcal{X}}_{n}$,

[TABLE]

• •

There are $k$ and $c$ so that, for any $A\in\mathcal{A}$ and any $x\in A$,

[TABLE]

(2.18) is the usual super-stability condition which states that the energy is the sum of a stable hamiltonian plus a positive quadratic term. The assumption on the derivatives in the last condition will be used when studying the dynamics. To prove the existence of DLR measures in the thermodynamic limit we need more assumptions which are not stated because we will use the above definition only in finite volumes.

When studying dynamics for general hamiltonians we will first introduce a cutoff, use it to prove existence and finally show that it can be removed. We use the following notation: ${\Lambda}$ and $\Delta$ always denote sets in ${\mathcal{X}}_{n}$, their complement being meant as the complement in ${\mathcal{X}}_{n}$. Let ${\Lambda}$ be a bounded set $\phi_{\Lambda}$ and $\phi_{{\Lambda}^{c}}$ configurations in ${\Lambda}$ and its complement, we then set

[TABLE]

We next introduce the cutoff function $g_{R}(\xi)$, $\xi\in\mathbb{R}_{+}$, $R>1$, by setting $g_{R}(\xi)=1$ when $\xi\leq R-1$, $g_{R}(\xi)=0$ when $\xi\geq R$ and $g_{R}(\xi)$ a decreasing $C^{\infty}$ function of $\xi$ in $(R-1,R)$ which interpolates between the values 1 and 0.

Definition 2.13** (“Cutoff Hamiltonians”).**

The general Hamiltonian $H$ with cutoff $R>1$ is:

[TABLE]

where $H^{0}_{\Lambda}(\phi_{\Lambda})$ is defined in (2.18) and

[TABLE]

Thus, when $\|\phi_{\Lambda}\|_{2}^{2}>R$, the Hamiltonian $H_{{\Lambda},R}(\phi_{\Lambda}|\phi_{{\Lambda}^{c}})$ becomes quadratic with no interaction among charges.

2.4. Dynamics.

The stochastic Ginzburg Landau model on $\mathcal{X}_{n}$ describes the time evolution of variables $\phi_{x}(t)$ which represent the amount of “charge” at site $x\in\mathcal{X}_{n}$ at time $t\geq 0$. The evolution is governed by the infinite system of stochastic differential equations

[TABLE]

where the variables $w_{x,y}(t)$ are defined in a space $({\Omega},P)$ as follows. An element ${\omega}\in{\Omega}$ is the collection $\{B_{x,y}(t),t\geq 0\}$ where $x,y$ run over the pairs $x\sim y$ such that $x<y$ in the lexicographic order. $P$ is a product measure such that each $\{B_{x,y}(t),t\geq 0\}$ is a standard Brownian motion. We then set

[TABLE]

Remark 2.14**.**

We will prove an existence theorem of the dynamics for “Positive interactions” (see the previous subsection) and for “General interactions” in the finite-volume case that we will describe below.

Remark 2.15**.**

Restrict the system (2.23) to only two equations, one for $x$ and the other for $y$ with $x\sim y$. By summing the two we see that the total charge $\phi_{x}(t)+\phi_{y}(t)$ is conserved thus the process describes exchanges of charges between the two sites. There is a random white noise term $dw_{x,y}(t)$, to which it is added a drift given by $\{\frac{\partial H}{\partial\phi_{x}}(\phi(t))-\frac{\partial H}{\partial\phi_{y}}(\phi(t))\}$ that we will call the instantaneous expected current from $x$ to $y$, which is thus defined as

[TABLE]

2.5. Partial dynamics.

As mentioned in Remark 2.14 above we will first study a partial dynamics where only finitely-many charges (those contained in a finite volume of linear size $N$) may evolve, while all the others are frozen at their initial values. The infinite volume dynamics will then be obtained in the limit $N\to\infty$.

The partial dynamics in $\Lambda_{n,N}$ freezes all charges outside $\Lambda_{n,N}$. We denote by $\phi$ a configuration in $\Lambda_{n,N}$ , by $\bar{\phi}$ a configuration outside $\Lambda_{n,N}$ and by $(\phi,\bar{\phi})$ a configuration in ${\mathcal{X}}_{n}$. We then write $\phi^{(N,\Delta,{\sigma})}(t)=\{\phi^{(N,\Delta,{\sigma})}_{x}(t|\phi,\bar{\phi},{\omega},\Delta,{\sigma}),x\in{\Lambda}_{n,N}\}$ for the solution (when it exists) of

[TABLE]

We interpret (2.26) by saying that at each bond $\{x,y\}$ with $x\in\Lambda_{n,N}$ and $y\in\Delta$ it is attached a reservoir which exchanges charges at a rate dictated by the chemical potential ${\sigma}_{y}$.

Remark 2.16**.**

Under Hypothesis 2.4(a) we are in the setup of the Fick’s law and the partial dynamics in (2.26) is customary in the analysis of boundary-driven processes, where the boundary processes simulate external reservoirs attached to the right and left faces of the system and generating currents. As we will see in $d\geq 3$ dimensions the currents do not decay as $N\to\infty$ so that the Fick’s law is violated in our geometrical setup. The Hypothesis 2.4(b) is used to study the infinite-volume limit. The choice of these boundary processes is therefore crucial in our analysis and it is at this point that the harmonic function ${\lambda}$ of Definition 2.1 enters into play.

We close this section by observing that the partial dynamics with a cut-off Hamiltonian is a Markov process, as the following proposition precisely states. For a general Hamiltonian $H$ we define the differential operator

[TABLE]

which acts on smooth functions as follows:

[TABLE]

[TABLE]

Proposition 2.17**.**

Let $H$ in (2.26) be a cutoff Hamiltonian (see Definition 2.13). Then, for any $\phi$, $\bar{\phi},\Delta,\sigma$, the equation (2.26) has solution $\phi^{(N,\Delta,\sigma)}(t)=\phi^{(N)}(t|\phi,\bar{\phi},{\omega},\Delta,\sigma)$ for $P$-almost all ${\omega}$. The law of $\{\phi^{(N,\Delta,\sigma)}(t),t\geq 0\}$, defines the transition probability starting from $\phi$ of a Markov diffusion process whose generator is $L^{n,N,\Delta,\sigma}$ in (2.27).

The equations (2.26) with the cutoff Hamiltonian have globally Lipschitz coefficients. The proof of Proposition 2.17 then follows, see for instance the book by Strook and Varadhan [21], and Chapter VII, §2 in [16].

3. Main results

3.1. Finite volumes

We fix arbitrarily $n$ and $N>n$, and shorthand $\phi=\{\phi_{x},x\in{\Lambda}_{n,N}\}$. We also fix $\bar{\phi}$, $\Delta$, $\sigma=\{{\sigma}_{y}\}_{y\in\Delta}$ and shorthand $\lambda^{*}_{x}=\lambda^{(N,\Delta,\sigma)}_{x}$, with $x\in\Lambda_{n,N}\cup\Delta$, see Definition 2.3. Recall that ${\lambda}^{*}_{y}={\sigma}_{y}$ for $y\in\Delta$. Let $\mu_{n,N,{\lambda}^{*}}(d\phi|\bar{\phi})$ be the Boltzmann-Gibbs measure

[TABLE]

where $H_{\Lambda}(\cdot|\bar{\phi})$ is defined in (2.20). The normalizing partition function is

[TABLE]

We will prove in Section 4 the following theorem.

Theorem 3.1**.**

For a general Hamiltonian $H$, let $L^{n,N,\Delta,\sigma}$ be as in (2.27) and $f$ a smooth test function, then

[TABLE]

Morever, if $H$ in (2.26) is a cutoff hamiltonian (in the sense of Definition 2.13), then $\mu_{n,N,{\lambda}^{*}}$ is an invariant measure for the partial dynamics.

We will use Theorem 3.1 to extend the invariance statement to general Hamiltonians. We denote by $P^{n,N,{\lambda}^{*}}(d\phi d{\omega}|\bar{\phi})=\mu_{n,N,{\lambda}^{*}}(d\phi|\bar{\phi})\times P(d{\omega})$ where $P(d{\omega})$ is the law of the Brownian motions $B_{\{x,y\}}(t)$ used to define the dynamics. Furthermore we write $\phi^{(N)}(t)=\{\phi^{(N)}_{x}(t|\phi,\bar{\phi},{\omega},\Delta,{\sigma}),x\in{\Lambda}_{n,N}\}$ for the solution (when it exists) of (2.26) with initial datum $\phi$. We will prove in Section 5 and Appendix B the following theorem.

Theorem 3.2**.**

With the above notation, for any $\bar{\phi}$, $\Delta$ and ${\sigma}$ there is a solution $\phi^{(N)}(t)=\phi^{(N)}(t|\phi,\bar{\phi},{\omega},\Delta,{\sigma}),\;t\geq 0$ of (2.26) for $P^{n,N,{\lambda}^{*}}$-almost all $(\phi,{\omega})$. Moreover for any test function $f$

[TABLE]

Finally, recalling (2.25) for notation,

[TABLE]

Non-validity of Fick’s law. With reference to Hypothesis 2.4(a), and using Propositions 2.2 and 2.6 the above theorem states that, in the limit $N\to\infty$, the current (3.5) is not identically zero, against what stated in the Fick’s law.

To study the infinite-volume dynamics we will use that, under Hypothesis 2.4(b), ${\lambda}^{*}={\lambda}$ and that DLR measures with chemical potential ${\lambda}$ are invariant under the partial dynamics, a statement that we specify next. For this we need more complete notation. We thus write $\phi$ for a configuration on ${\mathcal{X}}_{n}$, $\phi_{{\Lambda}}$ and $\phi_{{\Lambda}^{c}}$ for its restriction to ${\Lambda}_{n,N}$ and ${\Lambda}_{n,N}^{c}$. For $\Delta$ and $\sigma$ as in Hypothesis 2.4(b) we define an evolution on configurations $\phi$ by setting

[TABLE]

whenever the right hand side is well-defined, the definition being non empty because of Theorem 3.2. We will prove in Section 5 the following theorem.

Theorem 3.3**.**

Let $\mu$ be a DLR measure for the formal hamiltonian $H-\sum_{x}{\lambda}_{x}\phi_{x}$ and $\mathcal{P}=\mu\times P$. Then, under Hypothesis 2.4(b), for any $N>n$ and any test function $f$,

[TABLE]

Thus the DLR measures are stationary for all partial dynamics. However the existence of DLR measures for the general Hamiltonians of Definition 2.12 is an assumption, more conditions being needed to ensure their existence, for instance those stated in Definition 2.8 for positive interactions.

3.2. Infinite volume

We restrict here to positive hamiltonians $H$ in the sense of Definition 2.8 and for notational simplicity we consider the specific case of the Dinanburg-Sinai hamiltonian defined in (2.15). Let ${\lambda}$ be the harmonic function of Definition 2.1, $\mu$ a regular DLR measure with formal hamiltonian $H-\sum\phi_{x}{\lambda}_{x}$. By regular we mean that it is supported by configurations $\phi$ such that, for all $x$ large enough, $|\phi_{x}|\leq(\log|x|)^{1/3}$.

We call $\mathcal{P}=\mu\times P$ with $P(d{\omega})$ the law of the Brownians which define the dynamics. We then write $\phi(t)=\phi(t|\phi,{\omega})$ as the solution (when it exists) of (2.23) with initial datum $\phi$. We will prove in Section 6 and Appendix C the following result.

Theorem 3.4**.**

With $\mathcal{P}$-probability 1 there is a solution $\phi(t)=\phi(t|\phi,{\omega})$ of (2.23). For any test function $f$ and any $t>0$

[TABLE]

so that $\mu$ is time-invariant. Finally, recalling (2.25) for notation,

[TABLE]

Theorem 3.4 proves the claim, stated in the introduction, that there are stationary measures in infinite volumes carrying a non zero current. The theorem will be proved by showing that the solution of the partial dynamics converges, as $N\to\infty$, to $\phi(t)$.

Validity of Fick’s law. Let $\xi\in\{-n,\ldots,n\}$ then the stationary current per unit-area through a section $\Sigma_{\xi}$ in the channel is, by (3.9),

[TABLE]

By Proposition 2.2 it follows that ${\mathcal{J}}_{\xi}$ does not depend on $\xi$ and it is bounded by $c/n$. Contrary to what stated after Theorem 3.2 this shows the validity of the Fick’s law if we think of the system as the channel with the semi-spaces $\mathbb{Z}^{d}_{n,\pm}$ as “gigantic” reservoirs. They provide a steady current in the channel but, despite that, they do not change in time: this has evidently to do with the fact that they are infinite, but this is not enough to explain the phenomenon because in $d=2$ the effect is not present.

3.3. The quadratic hamiltonian

The quadratic Hamiltonian in the title is the one defined in (2.16). Being quadratic it may be seen as a cutoff hamiltonian so that the properties stated in Theorem 2.17 apply. In particular, for any $\phi$, $\bar{\phi}$, $\bar{\lambda}$, (2.26) has solution $\phi^{(N)}(t)$ for $P$-almost all ${\omega}$. Moreover the quadratic hamiltonian fits in the class of positive hamiltonians so that Theorem 3.4 applies and the infinite-volume dynamics $\phi(t|\phi,{\omega})$ is well defined with $\mathcal{P}$-probability 1 and the DLR measure $\mu_{\lambda}$ with chemical potential ${\lambda}$ (which is a product measure) is time-invariant.

We have however much more information, in fact, for an Ornestein-Ulhenbeck process it is known that Gaussian measures evolve into Gaussian measures, so that we only need to determine mean and covariance of the process. In our case this can be done using duality. Duality for the quadratic Ginzburg-Landau model follows from the algebraic approach discussed in [2], see [13] for a derivation based on Lie algebra representation theory. For completeness we shall also provide a direct proof in Section 7.

For finite volumes $\Lambda_{n,N}$ duality is stated as follows. Given $\Delta$ and $\sigma$ the duality function is

[TABLE]

where $h_{n}(\xi)$ with $\xi\in\mathbb{R}$ denotes the Hermite polynomial of degree $n$ and $\eta\in\mathbb{N}^{\Lambda_{n,N}}$ with $|\eta|=\sum_{x\in\Lambda_{n,N}}\eta_{x}<\infty$. Duality relates the Ginzburg-Landau evolution of $\{\phi^{N}(t),t\geq 0\}$ to the evolution of the Markov process $\{\eta^{N}(t),t\geq 0\}$ with generator

[TABLE]

where

[TABLE]

with $\eta^{x,y}$ the configuration obtained from $\eta$ by moving a particle from site $x$ to site $y$ and

[TABLE]

Thus the dual process is made of independent particles with absorptions at $\Delta$. We denote by ${\mathcal{E}}_{\eta}$ the expectation with respect to the law of the process $\{\eta^{N}(t),t\ \geq 0\}$ started at $\eta$. Similarly, we denote by ${{E}}_{\phi}$ the expectation with respect to the law of the process $\{\phi^{N}(t),t\geq 0\}$ started at $\phi$. We will prove in Section 7 the following result.

Theorem 3.5**.**

With the above notation we have

[TABLE]

Remark 3.6**.**

Using duality, the mean and covariance of the Gaussian process $\phi^{N}(t)$ can be computed starting the dual process with one and two dual particles. Furthermore duality also implies convergence in the limit $t\to\infty$ to the Gibbs measure $\mu_{n,N,{\lambda}^{*}}(d\phi)$ given by

[TABLE]

where ${\mathcal{Z}}$ is a normalizing constant. Indeed the duality formula (3.15) gives

[TABLE]

Expression (3.16) follows by recalling that, for a Gaussian random variable $Y$ with mean $m$, one has

[TABLE]

Similarly one can check invariance of 3.16. Labelling the particles of $\eta$, we describe $\eta$ as a configuration $X=\{X_{i},i=1,\ldots,|\eta|\}$ where the particles evolve independently and $X_{i}(t)$ is the position of the $i^{th}$ particle at time $t\geq 0$. We have

[TABLE]

where in the last equality it has been used that $\lambda$ is harmonic. Thus

[TABLE]

Remark 3.7**.**

The duality formula can be used to characterize the measure at infinite volume (by taking the $N\to\infty$ limit in (3.15)) and to show existence of the infinite-volume dynamics for general initial conditions $\phi$ which may grow polynomially at infinity.

Remark 3.8**.**

There is a large class of models where duality holds, including both particle systems (symmetric exclusion, Kipnis-Marchioro-Presutti model, independent particles, symmetric inclusion) and several interacting diffusions. We refer to the survey in preparation [2]. Results similar to those of the Ginzburg-Landau model with quadratic Hamiltonian can be obtained in models where duality holds.

4. Proof of Theorem 3.1

Equation (3.3) will be proved via an explicit computation that uses that ${\lambda}^{*}$ is harmonic. This generalizes a previous computation by De Masi, Olla, Presutti in [3]. In this section we shorthand $L=L^{n,N,\Delta,\sigma}$. We have

[TABLE]

where $\langle f,g\rangle_{\mu_{n,N,0}}:=\int f(\phi)g(\phi)\,\mu_{n,N,0}(d\phi)$ and $L^{\dagger}$ denotes the adjoint in $L^{2}(\mu_{n,N,0})$. Hence, to prove (3.3) it is enough to show that

[TABLE]

We compute the adjoint

[TABLE]

As in [3], we have

[TABLE]

Thus we find

[TABLE]

where $x,y\in\Lambda_{n,N}$ and we have defined

[TABLE]

Similarly

[TABLE]

where $x\in\Lambda_{n,N}$ and $y\in\Delta$. Hence (4.2) is equivalent to

[TABLE]

Changing from a sum over bonds to a sum over neighboring sites, we can rewrite this as

[TABLE]

which is clearly satisfied as a consequence of Prop. 2.6. ∎

5. Proof of Theorems 3.2 and 3.3

For brevity we call ${\Lambda}={\Lambda}_{n,N}$, and $\mu(d\phi):=\mu_{n,N,{\lambda}^{*}}(d\phi|\bar{\phi})$. We write $H_{{\Lambda},R}(\phi|\bar{\phi})$ and $\mu^{(R)}$ when we consider the hamiltonian with cutoff $R$, see Definition 2.13. We have already proved that the stochastic differential equations (2.26) with the cutoff hamiltonian have, for any initial datum, global solution with $P$-probability 1, they are denoted here by $\phi^{(R)}(t)$. By what proved in the previous section the Gibbs measure $\mu^{(R)}$ (with the chemical potential ${\lambda}^{*}$) is invariant. We will exploit this to prove a “time super-stability estimate”. We write $\mathcal{P}^{(R)}=\mu^{(R)}\times P$, and for any configuration $\phi$ in ${\Lambda}_{n,N}$,

[TABLE]

Then we have:

Theorem 5.1**.**

Given $T>0$ there are $A>0$ and $B$ (independent of $R$) so that, for all $S>2$,

[TABLE]

Theorem 5.1 will be proved in Appendix B.

We will next prove that we can replace $\mu^{(R)}$ by $\mu$ in (5.1).

Proposition 5.2**.**

Calling $d\mu^{(R)}(\phi)=G^{(R)}(\phi)d\phi$ and $d\mu(\phi)=G(\phi)d\phi$, we have

[TABLE]

where

[TABLE]

Proof. Call $Z^{(R)}$ the partition function then

[TABLE]

so that

[TABLE]

The analogous formula holds for $G(\phi^{\prime})$ so that calling

[TABLE]

one has

[TABLE]

hence (5.2). ∎

Corollary 5.3**.**

There are $A>0$ and $B$ so that, calling $\mathcal{A}=\{\sup_{t\leq T}\|\phi^{(R)}(t)\|_{2}^{2}\geq S\}$,

[TABLE]

Proof. By (5.2) the left hand side of (5.4) is bounded by

[TABLE]

and (5.4) follows from (B.1). ∎

Existence. There are $a^{\prime}>0$ and $b^{\prime}$ so that, for $R>S$,

[TABLE]

having used (5.4) and (5.1). Therefore, calling $\phi(t)$ the solution of (2.26), we have also

[TABLE]

because $\phi^{(R)}(t)=\phi(t)$ in the set

[TABLE]

Thus

[TABLE]

hence the existence of solutions to (2.26) with probability 1.

Time invariance. It is enough to prove that, given any $t>0$,

[TABLE]

for any test function $f$ such that $\sup_{\phi}|f(\phi)|\leq 1$.

Given any ${\epsilon}>0$, let $S$ be such that

[TABLE]

for any $R>S$, then

[TABLE]

For $R>S$

[TABLE]

and

[TABLE]

Since

[TABLE]

we get

[TABLE]

∎

Average current: proof of (3.5). From time-invariance,

[TABLE]

Then (3.5) easily follows, using integration by parts, from the definition of the current (2.25) and the explicit expression for the stationary measure (3.1). ∎

Proof of Theorem 3.3. Fix $N$, we claim that any measure $p$ on $\mathcal{X}_{n}$ of the form $dp(\phi)=d\nu(\phi_{{\Lambda}_{n,N}^{c}})d\mu_{n,N,{\lambda}^{*}}(\phi_{{\Lambda}_{n,N}}|\phi_{{\Lambda}_{n,N}^{c}})$ is invariant under the partial dynamics with $N$. By choosing $\nu$ to be the restriction of $\mu$ to configurations on the complement of ${\Lambda}_{n,N}$ and by using the DLR property, we will then get the invariance statement in the theorem. Let $f(\phi)$ be a smooth test function and let $g_{\phi_{{\Lambda}_{n,N}^{c}}}(\phi_{{\Lambda}_{n,N}}):=f(\phi_{{\Lambda}_{n,N}},\phi_{{\Lambda}_{n,N}^{c}})$. By Theorem 3.2 we get

[TABLE]

∎

6. Proof of Theorem 3.4

We are in the setup of Hypothesis 2.4(b) so that, by item (3) of Proposition 2.6, ${\lambda}_{x}^{N,\Delta,\sigma}={\lambda}_{x}$, throughout the section ${\lambda}$ is the harmonic function in Definition 2.1. We will use the following shorthand notation: given $n$ and $N>n$ we denote by $\phi$ a configuration on ${\Lambda}_{n,N}$, by $\bar{\phi}$ a configuration in the complement of ${\Lambda}_{n,N}$ and by $\mu_{N,\bar{\phi},{\lambda}}(d\phi)$ the Gibbs measure with hamiltonian $H-\sum_{x}\phi_{x}{\lambda}_{x}$ and with boundary condition $\bar{\phi}$.

The starting point is again a time superstability estimate. We can not use the one proved in Appendix B because the parameters in the estimates are volume dependent. In Appendix C we will first prove an equilibrium superstability estimate.

Theorem 6.1**.**

There are $a>0$, $N_{0}>0$ and $b$ so that for all $N>N_{0}$ the following holds. Let $|\bar{\phi}_{x}|\leq(\log|x|)^{1/3}$ for all $x\notin{\Lambda}_{n,N}$, then for any $x_{0}\in{\Lambda}_{n,N/2}$ and $S>0$

[TABLE]

The bound on $\bar{\phi}$ is motivated by Corollary 6.3 stated below.

Definition 6.2**.**

We set:

• •

$M_{{\lambda}}$* is the set of all DLR measures $\mu$ with chemical potential ${\lambda}$ such that*

[TABLE]

• •

The set $\mathcal{G}$ of “good configurations” is:

[TABLE]

Corollary 6.3**.**

With the above notation:

• •

The set $M_{\lambda}$ is non-empty because, if $\phi\in\mathcal{G}$ then, calling $\bar{\phi}_{N}$ the restriction of $\phi$ to the complement of ${\Lambda}_{n,N}$, any weak limit point of $\mu_{N,\bar{\phi}_{N},{\lambda}}$ is in $M_{{\lambda}}$.

• •

If $\mu\in M_{{\lambda}}$ then, for any $a^{\prime}<a$, there is $b^{\prime}$ so that

[TABLE]

and therefore $\mu[\mathcal{G}]=1$.

Proof. The first statement follows from Theorem 6.1. Since \mu\Big{[}|\phi_{x}|\geq S\Big{]}\leq e^{-aS^{4}+b}, then

[TABLE]

which yields (6.3). ∎

We shall next extend the super-stability estimates to the time-dependent case. Given $n$ and $N$ as above and a boundary configuration $\bar{\phi}$ we consider the partial dynamics defined in Section 2.5 with hamiltonian (2.16) and denote by $\{\phi_{x}(t),x\in{\Lambda}_{n,N}\}$ the corresponding process. Recall that the charges outside ${\Lambda}_{n,N}$ are frozen to the initial value $\bar{\phi}$ and that the dynamics does not depend on the chemical potential ${\lambda}$. We denote by $P_{\mu_{N,\bar{\phi},{\lambda}}}=\mu_{N,\bar{\phi},{\lambda}}\times P$ the law of the process $\phi^{(N)}(t|\phi,{\omega})$ when it starts from $\mu_{N,\bar{\phi},{\lambda}}$. In the sequel we fix arbitrarily a positive time $T$ and study the process in the time interval $[0,T]$. Using that $\mu_{N,\bar{\phi},{\lambda}}$ is invariant we will prove in Appendix C the following:

Theorem 6.4**.**

Let $\bar{\phi}$ be a configuration in ${\Lambda}_{n,N}^{c}$ such that $|\bar{\phi}_{x}|\leq(\log|x|)^{1/3}$ for all $x\in{\Lambda}_{n,N}^{c}$. Then, given $T>0$ there are $A>0$ and $B$ (independent of $N$ and $\bar{\phi}$) so that, for all $S>2$,

[TABLE]

The infinite-volume limit. Here we prove Theorem 3.4. We need preliminarily to extend the super-stability estimates from conditional Gibbs measures to DLR measures.

Theorem 6.5**.**

For any $N$, any measure $\mu\in M_{{\lambda}}$ is invariant under the evolution $T_{t}^{(N)}(\phi,{\omega},{\lambda})$, see (3.6). Moreover for all $S>0$,

[TABLE]

where $a^{\prime}$ and $b^{\prime}$ are as in (6.3).

**Proof. ** Invariance has been already proved in Theorem 3.3. We condition the probability on the left hand side of (6.5) to the configuration $\bar{\phi}$ outside ${\Lambda}_{n,N}$. We can use the bound in (6.4) when $\bar{\phi}\in\mathcal{G}_{N}$ and get in this case the bound with $e^{-AS^{4}+B}$. The additional term comes from the contribution of the configurations $\bar{\phi}$ which are not in $\mathcal{G}_{N}$, their probability is bounded using (6.3). ∎

It follows from Theorem 6.5 that:

Corollary 6.6**.**

Let $\mu\in M_{{\lambda}}$ then

[TABLE]

Proposition 6.7**.**

There exist $a^{\prime\prime}>0$ and $b^{\prime\prime}$ so that in $\mathcal{G}^{\prime}_{N}$ we have:

[TABLE]

Proof. We need to bound the differences $|\phi^{(N)}_{x}(t)-\phi^{(2N)}_{x}(t)|$ (recall that $2N$ is the first integer after $N$ in the set $\{2^{n},n\in\mathbb{N}\}$) with $t\leq T$ and $N\geq N_{0}$. We have:

[TABLE]

The contribution of the two-body potential $V$ to $\frac{\partial H}{\partial\phi_{x}}$ is uniformly Lipschitz, the one body term is bounded as follows:

[TABLE]

because $(\log 2N)^{1/3}$ bounds $|\phi_{x}^{(i)}(s)|$, $i=N,2N$.

It then follows:

[TABLE]

We can iterate (6.10) $K$ times with $K$ the largest integer such that $N/8+4K\leq N/4$. After $K$ iterations we get (6.7). ∎

In the theorem below we write $T_{t}^{(N)}(\phi,{\omega}):=\phi^{(N)}(t|\phi,{\omega})$.

Theorem 6.8**.**

For any $(\phi,{\omega})\in\mathcal{G}^{\prime}$ there is $\phi(t|\phi,{\omega})$ which satisfies the infinite volume stochastic differential equations (2.23) and, for any $x$ and any $t\leq T$, $\phi^{(N)}_{x}(t|\phi,{\omega})$ has a limit when $N\to\infty$, that we denote by $\phi_{x}(t|\phi,{\omega})$:

[TABLE]

Moreover, if $\mu\in M_{\lambda}$, then for any test function $f$,

[TABLE]

Proof. Let $N$ be such that $(\phi,{\omega})\in\mathcal{G}^{\prime}_{N}$. Since it satisfies the equations (2.26) then $\phi_{x}^{(N)}(t)$, $t\in[0,T]$, is equi-continuous and bounded and therefore it converges by subsequences to a limit $\phi_{x}(t)$. The limit is independent of the subsequence because $\phi^{(N)}(t)$ is Cauchy by Proposition 6.7.

By the invariance of $\mu$ for the partial dynamics and Theorem 3.3, we have, for all $N$,

[TABLE]

Then (6.12) follows from (6.11) and the Lebesgue dominated convergence theorem. ∎

7. Proof of theorem 3.5

Recalling the notation in Section 3.3, the duality statement (3.15) is a consequence of the following:

Lemma 7.1**.**

For $x,y$ nearest neighbors in $\Lambda_{n,N}$ we have

[TABLE]

For $x\in\Lambda_{n,N}$ and $y\in\Delta$ nearest neighbors we have

[TABLE]

Proof.

To alleviate notation we do not write the argument of the polynomials. We have

[TABLE]

We regroup terms as follows

[TABLE]

and then use the following identities for Hermite polynomials

[TABLE]

[TABLE]

[TABLE]

to find

[TABLE]

Similarly, for the boundaries we have

[TABLE]

This can be rewritten as

[TABLE]

By using the identies (7.1) and (7.3) one arrives to

[TABLE]

∎

Appendix A

Proof of Proposition 2.2.

Item (1).   Introducing the notation $\mathcal{B}^{\pm}=\{X(\cdot)\in\mathbb{Z}^{d}_{\pm n,\pm}\>\text{definitively}\}$, we may write

[TABLE]

By letting the walker $X(\cdot)$ do its first jump to one of its neighbors, and calling $d_{x}$ the number of neighbors of $x$, we can write

[TABLE]

Inserting (A.2) into (A.1) we find

[TABLE]

from which it follows that ${\lambda}$ is an harmonic function.

Item (2).   The proof follows from classical estimates on the recurrence of random walks on $\mathbb{Z}^{d}$. This explains why we need a spatial dimension larger than two. In $d\leq 2$ the random walk $X(\cdot)$ comes back infinitively many times to the channel $C_{n}$ and therefore the only harmonic functions are the constant ones. For completeness we give some details.

To show that $p_{x}^{+}+p_{x}^{-}=1$ it is enough to show that if $d\geq 3$ then the random walk $X(\cdot)$ is definitively in the complement of $C_{n}$. Let

[TABLE]

and define the hitting time of $\{X(t),\;t\geq 0\}$ to $K_{n}$ as

[TABLE]

Denoting by $\mathbb{P}^{X}_{x}$ the law of the $X(\cdot)$ process started from $x$, and defining

[TABLE]

where $S_{n,N-1}$ has been defined in (2.5), we claim that $p_{x}^{+}+p_{x}^{-}=1$ is implied by

[TABLE]

Indeed, we distinguish the following cases:

• •

If $x\in\Lambda_{n,N}$ then with probability 1 the walker $X(\cdot)$ will hit $S_{n,N}^{+}\cap S_{n,N}^{-}$ in a finite time.

• •

If $x\in\mathbb{Z}^{d}_{n,+}\setminus\Lambda_{n,N}$ it can only reach $K_{n}$ after passing through $S_{n,N}^{+}$. The analogous statement holds for $x\in\mathbb{Z}^{d}_{n,-}\setminus\Lambda_{n,N}$.

To prove (A.7) we call $\{Y(t),\;t\geq 0\}$ the usual continous-time random walk on $\mathbb{Z}^{d}$ that jumps with intensity $1$ to any of its nearest neighboring sites. Classical estimates prove that, if $d\geq 3$ then for any compact set $K$,

[TABLE]

where $\mathbb{P}^{Y}_{x}$ denotes the law of the $\{Y(t),\;t\geq 0\}$ process started from $x$. We may couple $X$ and $Y$ in such a way that

[TABLE]

while

[TABLE]

If we call $T$ the first time when $X(t)\in K_{n}$ then $Y(T)\in(K_{n}\cup K_{n-1})$. Therefore the claim (A.7) follows from (A.8) with $K=K_{n}\cup K_{n-1}$.

To prove that $\lambda_{x}$ is non-constant as $x$ varies in ${\mathcal{X}}_{n}$ we observe that $p_{x}^{+}\to 1$ and $p_{x}^{-}\to 0$ when $x_{1}\to+\infty$ and the opposite occurs when $x_{1}\to-\infty$. This sufficies.

Item (3).   To show that the flux is the same on each section it is enough to prove that $I_{\xi}({\lambda})=I_{\xi+1}({\lambda})$ for $|\xi|<n$. For such a $\xi$ we write, recalling (2.6) and that ${\lambda}$ is harmonic,

[TABLE]

We work with the central section inside the channel, i.e. $\xi=0$ and prove that for each $x\in\Sigma_{0}$, we have ${\lambda}_{x-e_{1}}-{\lambda}_{x}\leq c/n$. This in turn follows if we prove that

[TABLE]

where we recall that $p_{z}^{\pm}$ is defined in (2.9) and we take $x^{0}=(x^{0}_{1},..,x^{0}_{d})$, $y^{0}=(y^{0}_{1},..,y^{0}_{d})$ with $x^{0}_{1}=0$, $y^{0}_{1}=-1$ and $x^{0}_{i}=y^{0}_{i}$ for $i=2,\ldots,d$.

Call $\{X(t),t\geq 0\}$ and $\{Y(t),t\geq 0\}$ two copies of the random walk process with generator $G$ in (2.6), starting respectively from $x^{0}$ and $y^{0}$. We will prove that there exists a coupling $Q$ of these two processes so that

[TABLE]

which clearly implies (A.12).

To define the coupling $Q$ it is convenient to realize the process $X(t)$ in terms of its coordinates $X_{i}(t)$. To each $i\in\{1,\ldots,d\}$ we associate an exponential clock which rings with intensity 2, all clocks are independent. When a clock rings we take a variable $\varepsilon$ with values $\pm 1$, all the $\varepsilon$-variables are mutually independent. If the $i$-th clock rings and $\varepsilon$ is the associated variable, then $X_{i}$ tries to jump: $X_{i}\to X_{i}+\varepsilon$, the jump is done if after the jump $X\in\mathcal{X}_{n}$, otherwise it is suppressed.

Definition of Q: The coupling $Q$ is a measure on the sample space $\Omega$ and we will define $X(t)$ and $Y(t)$ on $\Omega$. The elements $\omega\in\Omega$ are of the form $\omega=\{t_{n}^{i,x},\varepsilon_{n}^{i,x},t_{n}^{1,y},\varepsilon_{n}^{1,y}\}$ where $\;i\in\{1,\ldots,d\}$ and $n\in\mathbb{N}$. Under $Q$, the times $t_{n}^{i,x}$ and $t_{n}^{1,y}$ are realizations of Poisson processes of intensity $2$ and the increments $\varepsilon_{n}^{i,x}$ and $\varepsilon_{n}^{1,y}$ are realizations of Bernoulli processes with parameter $1/2$. All these processes are independent of each other. Thus $Q$ is completely defined.

Representation of $X(t)$: We define the processes $\{X_{i}(t)\}_{t\geq 0}$, with $i=1,\ldots,d$, as the collection of walkers that are initialized from $x^{0}_{i}$ and at the times $t_{n}^{i,x}$ jumps by $\varepsilon_{n}^{i,x}$ if the jump is allowed (the walker can not exit ${\mathcal{X}}_{n}$).

Representation of $Y(t)$: We first define the auxiliary processes $\{Y^{\prime}_{i}(t)\}_{t\geq 0}$ where $i=1,\ldots,d$. They start from $y^{0}_{i}$ and they use the variables $\{t_{n}^{i,x},\varepsilon_{n}^{i,x}\}$ for $i=2,\ldots,d$ and the variables $\{t_{n}^{1,y},\varepsilon_{n}^{1,y}\}$ for the first coordinate $\{Y^{\prime}_{1}(t)\}_{t\geq 0}$. To define $\{Y(t),t\geq 0\}$ we introduce the time $\bar{t}$ as

[TABLE]

and define $Y(t)=Y^{\prime}(t)$ for $t\leq\bar{t}$. Then $\{Y(t),t\geq\bar{t}\}$ is constructed by using the variables $\{t_{n}^{i,x},\varepsilon_{n}^{i,x}\}$ with $i=1,\ldots,d$ and starting at time $\bar{t}$ from $Y^{\prime}(\bar{t})$.

Clearly the law of $X(\cdot)$ is $\mathbb{P}_{x^{0}}$ and the law of $Y(\cdot)$ is $\mathbb{P}_{y^{0}}$, thus $Q$ defines the desired coupling. Having defined the coupling $Q$, we now start the analysis of (A.13). To this aim, it is convenient to define two processes $X_{1}^{*}(t)$ and $Y_{1}^{*}(t)$. The process $\{X^{*}_{1}(t),\;t\geq 0\}$ is initialized from [math] and, at times $t_{n}^{1,x}$, jumps by $\varepsilon_{n}^{1,x}$ (with no restrictions). We define similarly the auxiliary process $\{Z_{1}(t),\;t\geq 0\}$: it starts from $-1$ and it uses the variables $\{t_{n}^{1,y},\varepsilon_{n}^{1,y}\}$. To define $\{Y_{1}^{*}(t),t\geq 0\}$ we introduce the time $t^{*}$ as the first time when $Z_{1}(t)=X^{*}_{1}(t)$ and define $Y_{1}^{*}(t)=Z_{1}(t)$ for $t\leq t^{*}$. $\{Y_{1}^{*}(t),t\geq t^{*}\}$ is constructed by using the variables $\{t_{n}^{1,x},\varepsilon_{n}^{1,x}\}$ and starting at time $t^{*}$ from $Y^{\prime}(t^{*})$. As a consequence, $X^{*}_{1}(t)=Y^{*}_{1}(t)$ for $t\geq t^{*}$.

We introduce a stopping time $\tau$ as

[TABLE]

and observe that the following three properties hold true:

• •

[TABLE]

• •

[TABLE]

• •

[TABLE]

Thus, by law of total probability, we may write

[TABLE]

where (A.15) has been used in the second term of the r.h.s. By using (A.16) we have

[TABLE]

By classical estimates for the maximum of a random walk, there exist $a>0$ and $b$ so that

[TABLE]

As a consequence of (A.16) we also have that

[TABLE]

Thus

[TABLE]

for some constant $c^{\prime\prime}$. Choosing $T=n^{2}$ we then get (A.13) because, by (A.17),

[TABLE]

Proof of Proposition 2.6.

Item (1).   We recall that

[TABLE]

Similarly to item (1) of Proposition 2.2, by letting the walker $X^{N,\Delta}(\cdot)$ do its first jump to one of its neighboring sites, and calling $d_{x}$ the number of neighbors of $x\in\Lambda_{n,N}$, we can write

[TABLE]

Inserting (A.25) into (A.24) we find

[TABLE]

from which it follows that $G^{N,\Delta}\lambda^{(N,\Delta,{\sigma})}(x)=0$, i.e. it is an harmonic function.

Item (2).   We are in the setting of Hypothesis 2.4(a). For notational simplicity we assume in this section that the spatial dimension is fixed to $d=3$. For any $x$, let $k$ be an integer so that $x\in\Lambda_{n,k}$. We define

[TABLE]

and we have

[TABLE]

We call

[TABLE]

then

[TABLE]

Similarly, writing $X^{N}$ for $X^{N,\Delta}$, we take $N>k$ and define

[TABLE]

Then, writing ${\lambda}^{N}$ for ${\lambda}^{(N,\Delta,\sigma)}$, we have

[TABLE]

We define

[TABLE]

and get

[TABLE]

Since $X(t)=X^{N}(t)$ for $t\leq\tau_{k}$ then combining (A.30) and (A.34) we find

[TABLE]

By (A.7) ${\epsilon}_{k}\to 0$ as $k\to\infty$ so that we only need to bound ${\epsilon}_{k,N}$. Let

[TABLE]

and

[TABLE]

Then

[TABLE]

It is convenient to change coordinates $x_{1}\to x_{1}-n$. For notational simplicity we rename $N$ instead of $N-n$ and $k$ instead of $k-n$.

The aim is to reduce to an estimate on the simple symmetric random walk in $\mathbb{Z}^{3}$, that we shall call $\{Y(t),t\geq 0\}$. This can be done by generalizing the argument in the proof of item (2) of Prop. 2.2. Thus we introduce a set

[TABLE]

where $D_{0,0}=D$, $D_{m_{2}\pm 1,m_{3}}=D_{m_{2},m_{3}}\pm(2N+1)e_{2}$ which means that $D_{m_{2},m_{3}}$ is translated by $\pm(2N+1)e_{2}$. Analougously $D_{m_{2},m_{3}\pm 1}=D_{m_{2},m_{3}}\pm(2N+1)e_{3}$. We define

[TABLE]

[TABLE]

[TABLE]

and

[TABLE]

Then

[TABLE]

The proof that ${\epsilon}_{k,N}$ vanishes for $k$ and $N$ large follows from:

Proposition A.1**.**

There exist positive constants $c_{1},c_{2},c_{3},a$ such that for all $x\in A_{k}^{+}$

[TABLE]

Proof: The proof follows from classical estimates for random walks, in particular the estimate for the Green function $G(x,z)\leq\frac{c}{|y-z|}$ in 3 dimensions, see [20] Chapter 6, §26. Let

[TABLE]

and

[TABLE]

Then

[TABLE]

We have

[TABLE]

where the first term bounds the contribution of $D_{0,0}=D$ and the second term comes from

[TABLE]

and

[TABLE]

By classical estimates on the displacement of a random walk we obtain

[TABLE]

Hence (A.44) is proved. ∎By letting first $N\to\infty$ and then $k\to\infty$ we obtain ${\lambda}_{x}^{N}\to{\lambda}_{x}$ from (A.35), after recalling that, by (A.7), ${\epsilon}_{k}\to 0$.

Item (3).   We recall that

[TABLE]

so that

[TABLE]

By using that $X^{N,\Delta}(t)=X(t)$ for $t\leq\tau$ we thus find ${\lambda}_{x}={\lambda}_{x}^{N,\Delta,\sigma}$.

Appendix B

The following is a weak (volume-dependent) form of Ruelle’s superstability estimates.

Lemma B.1**.**

There is $c$ (which does depend on ${\Lambda}$ and $\bar{\phi}$, but since they are fixed we may regard $c$ as a constant) so that

[TABLE]

Proof. Let $\Delta$ be the finite set of points in ${\Lambda}^{c}$ interacting with those in ${\Lambda}$. Recalling Definition 2.12 for notation, we write

[TABLE]

Hence, by (2.18),

[TABLE]

The latter term can be regarded as a constant because it only depends on $\bar{\phi}$. Thus $H^{\prime}_{{\Lambda},R}(\phi|\bar{\phi})\geq-g(\|\phi\|_{2}^{2})(|B||{\Lambda}\cup\Delta|+|H^{\prime}_{\Delta}(\bar{\phi}_{\Delta})|)\geq c$ and therefore

[TABLE]

We bound from below the partition function by restricting the integral to $|\phi_{x}|\leq 1$ for all $x\in{\Lambda}$ and we obtain (B.1). ∎

To extend the bound to time intervals we will use the following theorem which will be used again in Appendix C.

Theorem B.2**.**

Let $z(t)$, $t\in[0,T]$, $T>0$ be a process with law $P$. Suppose that for $t\leq T$

[TABLE]

with $M_{t}$ and $N_{t}$ martingales ($N_{0}=0$) and that

[TABLE]

Then

[TABLE]

Proof. Since it is short we give for completeness the proof which can be found in Holley and Strook [14] and in De Masi and Presutti [4]. We write

[TABLE]

By Cauchy-Schwartz

[TABLE]

which is the first term on the right hand side of (B.4).

By Doob’s theorem

[TABLE]

By (B.2)

[TABLE]

which completes the proof of the theorem recalling that $M_{0}=z(0)$. ∎

Proof of Theorem 5.1. Given $S>2$ we define a smooth function $f(\xi)$, $\xi\geq 0$, in such a way that $f(\xi)=1$ for $\xi\geq S$ and $f(\xi)=0$ for $\xi\leq S-1$. When $|\xi|\in[S-1,S]$, $f(\xi)$ is a strictly increasing $C^{\infty}$ function with [math] derivatives at the endpoints. As a consequence $f(\xi)$ is a smooth non decreasing function with derivatives bounded uniformly in $S$, and if the derivatives $f^{\prime}(\xi)\neq 0$ or $f^{\prime\prime}(\xi)\neq 0$ then $\xi\in[S-1,S]$.

Writing $f_{t}$ for $f(\|\phi(t)\|_{2}^{2})$ (5.1) reads as

[TABLE]

which is implied by

[TABLE]

We will bound (B.6) using Theorem B.2 with $z(t)=f_{t}$ and

[TABLE]

where $L$ is the generator $L^{n,N,{\lambda}}$ (with cutoff $R$) of (2.27). The role of the measure $P$ in Theorem B.2 is now taken by $\mathcal{P}^{(R)}$. Since this is time-invariant we get from (B.4)

[TABLE]

$\bullet$  Bound of the term with ${\gamma}_{1}$.

Recalling (2.27), a contribution to $Lf$ comes from the first order derivatives and it is a finite sum of terms of the form

[TABLE]

while the second order derivatives give rise to a sum of terms of the form

[TABLE]

The key point is that both $|f^{\prime}|$ and $|f^{\prime\prime}|$ are bounded by $\leq c\mathbf{1}_{\|\phi\|_{2}^{2}\in[S-1,S]}$ and since $g$ and its derivatives are bounded, the expectation of ${\gamma}_{1}^{2}$ is bounded by the sum of finitely-many terms like

[TABLE]

where $\pi(\phi)$ is a polynomial in $\phi$. By (B.1) this is bounded by $ce^{-(S-1)(a\beta)/2}$, with $c$ depending on $\pi(\phi)$, and in conclusion:

[TABLE]

. $\bullet$  Bound of the term with ${\gamma}_{2}$.

We have

[TABLE]

then, by (B.7), we are reduced to the analysis of terms as those considered for ${\gamma}_{1}$ and we get

[TABLE]

$\bullet$  Bound of the last term in (B.8). We use (B.1) to write

[TABLE]

This concludes the proof. ∎

Appendix C

We start by extending the super-stability estimates in [17], [15] and [1] to the present case, namely for the Hamiltonian (2.15) to which it is added the contribution of a chemical potential ${\lambda}$ which is a harmonic function. This is a special super-stable hamiltonian where the one body term is

[TABLE]

(with $B$ a suitable constant) and where the two-body potential, $|\phi_{x}-\phi_{y}|^{2}$, is nearest neighbor and evidently non negative. We will exploit all that to simplify the proofs given in the general case.

Proof of Theorem 6.1. Dropping the dependence on ${\lambda}$ in the notation we want to bound

[TABLE]

where $H_{{\Lambda}_{n,N}}$ is the energy in ${\Lambda}_{n,N}$ (which includes the chemical potential ${\lambda}$) and $W_{{\Lambda}_{n,N}|{\Lambda}_{n,N}^{c}}$ is the interaction between the charges in ${\Lambda}_{n,N}$ and those in ${\Lambda}_{n,N}^{c}$.

The idea in [17], [15] and [1] is to estimate the integral in (C.2) by introducing a stopping time. To this end we denote by $\Delta_{q}$ the cubes of side $2q+1$ centered at $x_{0}$ taking $q\geq q_{0}$ where $q_{0}$ is such that, for $q\geq q_{0}$,

[TABLE]

We choose $N_{0}$ so that, for $N\geq N_{0}$, $\Delta_{q_{0}}\subset{\Lambda}_{n,N}$. We often write in the sequel for brevity ${\Lambda}={\Lambda}_{n,N}$, ${\Lambda}^{\prime}={\Lambda}_{n,N+1}$ and $Z_{\Lambda}(\bar{\phi})=Z_{{\Lambda}_{n,N}}(\bar{\phi})$.

We partition the configurations on $\mathcal{X}_{\Lambda}$ into the following atoms:

[TABLE]

and, for $q>q_{0}$,

[TABLE]

Thus $q^{\prime}$ stops as soon as $\sum_{x\in\Delta_{q}\cap{\Lambda}}\phi_{x}^{2}\leq|\Delta_{q}|\log q$.

We call $\rho_{q}(\phi_{x_{0}})$ the integral in (C.2) extended to $\mathcal{A}_{q}$ so that

[TABLE]

We split the terms $\rho_{q}(\phi_{x_{0}})$ into three classes.

$\bullet$  $\rho_{q_{0}}(\phi_{x_{0}})$. Here we will prove the bound (C.11) below. We first drop the non negative interaction between the charges in $\Delta_{q_{0}-1}$ and those in the complement getting

[TABLE]

because by the assumption on $q_{0}$ there is no interaction between $\Delta_{q_{0}}$ and the complement of ${\Lambda}$. By (C.1) we have

[TABLE]

We use the term with $\phi_{x}^{2}$ to perform the integrals over the variables $\phi_{x},x\in\Delta_{q_{0}-1}\setminus x_{0}$ so that

[TABLE]

To reconstruct a partition function we write $|\phi_{\Delta_{q_{0}-1}}|\leq 1$ for the set where $|\phi_{x}|\leq 1$ for all $x\in\Delta_{q_{0}-1}$. Then there is $c^{\prime}$ such that

[TABLE]

We claim that $2|\Delta_{q_{0}}|(\log q_{0}+2d)\geq W_{{\Lambda}\setminus\Delta_{q_{0}-1}|\Delta_{q_{0}-1}}$. Proof: let $x\in\Delta_{q_{0}-1}$ and $y\in\Delta_{q_{0}}$, $x\sim y$. We bound $(\phi_{x}-\phi_{y})^{2}\leq 2(\phi_{x}^{2}+\phi_{y}^{2})$. By (C.4) the sum over all such $y$ is bounded by $2|\Delta_{q_{0}}|\log q_{0}$ while the sum over all such $x$ is bounded by $2d(|\Delta_{q_{0}}|-|\Delta_{q_{0}-1}|)$ hence the claim. We then get

[TABLE]

By (C.9) and (C.10) we then finally get:

[TABLE]

$\bullet$  $\rho_{q}(\phi_{x_{0}})$ with $q$ such that $\Delta_{q}\subset{\Lambda}$. With the same procedure we get the analogue of (C.9):

[TABLE]

Finally we need an analogue of (C.10) to reconstruct the partition function. By (C.3) we have

[TABLE]

and proceeding as before we get

[TABLE]

$\bullet$  $\rho_{q}(\phi_{x_{0}})$ with $q$ such that $\Delta_{q}\cap{\Lambda}^{\prime}\neq\emptyset$. We integrate as before over the charges in $\Delta_{q-1}\cap{\Lambda}$ and drop the interaction between ${\Lambda}\setminus\Delta_{q-1}$ and $\Delta_{q-1}\cap{\Lambda}$ as well as the interaction between $\Delta_{q-1}\cap{\Lambda}$ and ${\Lambda}^{c}$. We then get again

[TABLE]

The reconstruction of the partition function is now more complicated because we need to take into account also the interaction between $\Delta_{q-1}\cap{\Lambda}$ and ${\Lambda}^{c}$. We call $B^{\rm in}$ the set of points in ${\Lambda}\setminus\Delta_{q-1}$ which are in a bond with a point in $\Delta_{q-1}\cap{\Lambda}$. $B^{\rm out}$ instead is the set of points in ${\Lambda}^{c}$ which are in a bond with a point in $\Delta_{q-1}\cap{\Lambda}$. Thus we consider the partition function

[TABLE]

$W_{\Delta_{q_{-}1}\cap{\Lambda}|B^{\rm in}}$ is bounded as in (C.13).

[TABLE]

because $y\in{\Lambda}^{\prime}\setminus{\Lambda}$. We have $q\geq N/2$ (because $x_{0}\in{\Lambda}_{n,N/2}$ and the cube of side $2q+1$ and center $x_{0}$ has non-empty intersection with ${\Lambda}_{n,N}^{c}$). Therefore $\log N\leq\log 2q$ so that

[TABLE]

In conclusion we get

[TABLE]

The sum in (C.6) is then bounded as on the right hand side of (6.1) which is therefore proved. ∎

Proof of Theorem 6.4. Let $S>2$ and $f(\xi)$, $\xi\geq 0$, be the same as in Appendix B, thus $f(\xi)$ is a smooth non decreasing function with derivatives bounded uniformly in $S$, and if the derivatives $f^{\prime}(\xi)\neq 0$ or $f^{\prime\prime}(\xi)\neq 0$ then $\xi\in[S-1,S]$. We fix $x\in{\Lambda}_{n,N/4}$ and write $f_{t}$ for $f(|\phi_{x}(t)|)$. It is then enough to prove that

[TABLE]

$\mathcal{E}_{\mu_{N,\bar{\phi},{\lambda}}}$ being the expectation relative to the process with law $P_{\mu_{N,\bar{\phi},{\lambda}}}$.

Since $\phi(t)$ solves the stochastic differential equations (2.26) we have the following martingale decomposition:

[TABLE]

with $M_{t}$ and $N_{t}$ martingales ($N_{0}=0$) and

[TABLE]

see for instance Revuz and Yor, [16]. Since $\mu_{N,\bar{\phi},{\lambda}}$ is time invariant we get, analogously to (B.8),

[TABLE]

Recalling that $x\in{\Lambda}_{n,N/4}$ and that $f^{\prime}$ and $f^{\prime\prime}$ are bounded and equal to 0 unless $|\phi_{x}|\in(S-1,S)$, we have

[TABLE]

Thus, after summing over $|y-x|=1$,

[TABLE]

Moreover

[TABLE]

so that

[TABLE]

Thus

[TABLE]

In conclusion, after using Cauchy-Schwartz and (6.1)

[TABLE]

Since

[TABLE]

$E_{\mu_{N,\bar{\phi},{\lambda}}}[{\gamma}_{2}^{2}]$ is bounded in a similar way as well as $E_{\mu_{N,\bar{\phi},{\lambda}}}[f^{2}]$, we omit the details. ∎