Couplings via comparison principle and exponential ergodicity of SPDEs in the hypoelliptic setting

TL;DR

This paper introduces a general framework for analyzing ergodicity of order-preserving Markov semigroups, providing optimal conditions for invariant measures and demonstrating exponential ergodicity and synchronization in hypoelliptic stochastic reaction-diffusion equations.

Contribution

It develops a broad theoretical framework for ergodicity of order-preserving Markov processes and applies it to hypoelliptic SPDEs, refining previous results.

Findings

Established optimal conditions for invariant measure existence and uniqueness.

Proved exponential convergence of transition probabilities.

Demonstrated exponential synchronization-by-noise in hypoelliptic SPDEs.

Abstract

We develop a general framework for studying ergodicity of order-preserving Markov semigroups. We establish natural and in a certain sense optimal conditions for existence and uniqueness of the invariant measure and exponential convergence of transition probabilities of an order-preserving Markov process. As an application, we show exponential ergodicity and exponentially fast synchronization-by-noise of the stochastic reaction-diffusion equation in the hypoelliptic setting. This refines and complements corresponding results of Hairer, Mattingly (2011).

Equations

$$du(t,\xi )=[\Delta u(t,\xi )+f(u(t,\xi ),\xi )]dt+\sum _{k=1}^m{\sigma }_k(\xi )d{W}^k(t),\xi \in {\mathbb{T}}^d,t⩾0,$$

$$\Gamma :=\{(x,y)\in E\times E:x⪯y\}$$

$$W_r(\mu ,\nu ):=\underset{\lambda \in \mathcal{C}(\mu ,\nu )}{\inf }{\int }_{E\times E}r(x,y)\lambda (dx,dy),\mu ,\nu \in \mathcal{P}(E),$$

$${\int }_{E}fd\mu ⩽{\int }_{E}fd\nu .$$

$$P_t(x,\cdot ){⪯}_{st}{P}_t(y,\cdot ).$$

$$P_{t}V(x)⩽V(x)-\gamma {\int }_0^t{P}_{s}V(x)ds+Kt,t⩾0,x\in E;$$

$$P_1(x,A)>\epsilon \text{ and }{P}_1(x,B)>\epsilon .$$

$$W_{d\wedge 1}(P_t(x,\cdot ),P_t(y,\cdot ))⩽C(1+V(x)+V(y){)}^{\theta }(1+M(x){)}^{\theta }\exp (-\lambda t),t⩾0.$$

$$W_{\rho \wedge 1}(P_t(x,\cdot ),\pi )⩽C(1+V(x){)}^{\theta }\exp (-\lambda t),t⩾0.$$

$$\mathrm{Law}(Z^x)=\mathrm{Law}({\tilde{Z}}^x)=P_t(x,\cdot ),\mathrm{Law}(Z^y)=P_t(y,\cdot );$$

$$\mathsf{P}(Z^x⪯Z^y⪯{\tilde{Z}}^x)⩾1-C_1(1+V(x)+V(y))e^{-C_{2}t},$$

$$\mathsf{E}[d(Z^x,Z^y)\wedge 1]⩽1-C_3(1+V(x)+V(y))(1+M(x))e^{-C_{4}t},$$

$$x{⪯}_{inc}y\text{ if and only if }{x}_i⩽y_i\text{ for all i∈N}.$$

$$x{⪯}_{pos}y\text{ if and only if }x(\xi )⩽y(\xi )\text{ for almost all ξ∈D}$$

$$x{⪯}_{Besov}y\text{ if and only if }⟨x,\phi ⟩{⪯}_{pos}⟨y,\phi ⟩\text{ for any non-negative φ∈C∞}$$

$$\mathsf{E}[\rho (X_t,Y_t)\wedge 1]\to 0\text{as t→∞}.$$

$$d_{TV}(\mathrm{Law}(X_t),\mu )⩽r(t),d_{TV}(\mathrm{Law}(Y_t),\mu )⩽r(t),$$

$$Y_n:=X\text{mathbbm}1(X\notin [2^n,2^{n+1}])+(X+1)\text{mathbbm}1(X\in [2^n,2^{n+1}]);X_n=X.$$

$$d_{TV}(\mathrm{Law}(X_n),\mu )=0;d_{TV}(\mathrm{Law}(Y_n),\mu )⩽p_n{2}^{-n}.$$

$$\mathsf{E}[|X_n-Y_n|\wedge 1]=p_n,$$

$$|\phi (x)-\phi (y)|⩽d(x,y{)}^k(\psi (x)+\psi (y)),x,y\in E.$$

$$\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptd(X,Y)\leqslant\varepsilon^{k}\bigl{(}(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(X)^{1/(1-k)}])^{1-k}+(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(Y)^{1/(1-k)}])^{1-k}\bigr{)}.$$

$$\mathsf{E}d(X,Y)⩽\mathsf{E}\phi (Y)-\mathsf{E}\phi (X).$$

$$\mathsf{E}d(X,Y)$$

$$⩽\mathsf{E}[d(\tilde{X},\tilde{Y}{)}^k(\psi (\tilde{X})+\psi (\tilde{Y}))]$$

$$\displaystyle\leqslant K^{k}\varepsilon^{k}\bigl{(}(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(\widetilde{X})^{1/(1-k)}])^{1-k}+(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(\widetilde{Y})^{1/(1-k)}])^{1-k}\bigr{)}$$

$$\displaystyle=K^{k}\varepsilon^{k}\bigl{(}(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(X)^{1/(1-k)}])^{1-k}+(\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\psi(Y)^{1/(1-k)}])^{1-k}\bigr{)}.$$

$$\mathsf{P}(X⪯Y⪯\tilde{X})⩾1-\epsilon$$

$$\mathsf{E}d(X,Y)⩽2{\epsilon }^{1/p}(\mathsf{E}|\phi (X){|}^q{)}^{1/q}+\epsilon .$$

$$f_k(\xi )=1\text{for ξ∈[0,k/n]};f_k(\xi )=0\text{for ξ∈[(k+1)/n,1]}.$$

$$\phi (\check{1})=\phi (f_n)⩾\phi (f_1)+n-1⩾\phi (\check{0})+n,$$

$$\mathsf{E}d(X,Y)$$

Full text

**footnotetext: Email: oleg.butkovskiy@gmail.com. Supported in part by DFG Research Unit FOR 2402 and European Research Council grant 683164.****footnotetext: Email: ms@math.tu-berlin.de. Supported in part by DFG Research Unit FOR 2402.

1 Introduction

The goal of this article is to build a framework for analyzing ergodic properties of order–preserving Markov processes and to provide simple, verifiable, yet general enough sufficient conditions for exponential ergodicity. This framework turns out to be especially powerful for investigating ergodicity of order–preserving stochastic PDEs with highly degenerate additive forcing. Our main example is the stochastic reaction–diffusion equation in the hypoelliptic setting. We show that even if noise comes to the system only from one Brownian motion, then (under certain conditions) this SPDE has a unique invariant measure and its transition probabilities converges to it exponentially in the Wasserstein metric. We also establish exponentially fast synchronization–by-noise of the solutions to this equation. This refines [18, Remark 8.22] and complements [18, Theorem 8.21].

In the mathematical physics literature there is a growing interest in ergodic behavior of nonlinear PDEs forced by smooth in space noise acting only on a few Fourier modes, see, e.g., [29, 30, 16, 5, 23]. Since the noise is smooth in space, it is usually relatively easy to show that these SPDEs have a unique solution and that this solution is a Markov process. On the other hand, since the solution at any fixed time is an infinite–dimensional random variable and the noise acts only onto finitely many degrees of freedom, these processes do not get “enough” noise and hence they are typically only Feller but not strong Feller, see also the discussion in [15, Section 9]. This makes analyzing their ergodic behavior much more challenging.

Indeed, recall that ergodicity of strong Feller processes can be established using the standard classical approach, which combines a local mixing condition on a certain set (the small set condition) and a recurrence condition, see, e.g., [31]. Unfortunately, this method is usually not applicable for Markov processes which are only Feller and not strong Feller because they do not have good mixing properties, see also a detailed discussion in [19, Section 1]). To study ergodicity of these processes three alternative strategies have been suggested recently.

The first approach was developed in [16, 17, 18], see also [15, Section 11]. It introduces the asymptotic strong Feller (ASF) property, which serves as a replacement for the strong Feller property. It is shown there that if a Markov process satisfies ASF as well as certain recurrence and topological irreducibility conditions then the process is exponentially ergodic. Note that for many Markov processes verifying ASF might be rather challenging. In particular, while this method works quite well for stochastic Navier–Stokes (SNS) equations on a torus in the hypoelliptic setting, it is not so clear how to check ASF for the SNS equation on a bounded domain, see the discussion in [11, Section 1].

Another approach establishes exponential ergodicity using generalized couplings [30, 12, 19, 11, 24, 2]. Recall that a coupling is a pair of stochastic processes with given marginal distributions. By contrast, a generalized coupling is a pair of stochastic processes, whose marginals are not necessarily equal to a prescribed pair of probability distributions, but are in a sense close to this pair. Clearly, constructing a generalized coupling is much easier than constructing a coupling. Furthermore, it is shown in the papers mentioned above that existence of a generalized coupling with certain nice properties yields exponential ergodicity. This approach works quite well for a large class of SPDEs in the effectively elliptic setting (that is, when noise acts in a finite but large enough number of directions), but is less useful for studying SPDEs in the hypoelliptic setting.

Finally, the third main approach was introduced in [19]. It utilizes the notion of a $d$–small set (a generalization of a small set), which is particularly well adapted to the study of Markov processes with bad mixing properties. This approach provides another set of sufficient conditions for exponential ergodicity, and it works quite well with stochastic delay equations and nonlinear autoregressions. Unfortunately, verifying this set of conditions for SPDEs is rather difficult.

Our new approach developed in this paper is somehow orthogonal to all of the strategies mentioned above. It is specifically targeted at order preserving Markov semigroups, that is, the semigroups which map increasing bounded functions to increasing bounded functions. On the one hand, this significantly reduces the applicability of this approach; for example, it cannot be used to study stochastic Navier–Stokes equations. On the other hand, for order–preserving Markov processes (e.g., stochastic reaction–diffusion equations) it allows to obtain exponential ergodicity under very weak assumptions; this is rather difficult (or maybe impossible) to achieve with other methods.

The main result (Theorems 2.3 and 2.4) is quite general. It shows that if an order–preserving Markov semigroup satisfies additionally a swap condition (i.e., two Markov processes started with initial conditions $x$, $y$ with $x\preceq y$ can change their order by time $1$ with a small but positive probability), then under a standard Lyapunov–type condition as well as a certain technical assumption the process is exponentially ergodic. We also show that this swap condition cannot be omitted.

We apply the obtained theorems to establish exponential ergodicity of stochastic reaction–diffusion equations on a $d$–dimensional torus $\mathbb{T}^{d}$, $d\in\mathbb{N}$

[TABLE]

where $(W^{1},W^{2},\ldots,W^{m})$ are independent standard Brownian motions; $f$, $\sigma_{k}$ are continuous functions acting from $\mathbb{R}\times\mathbb{T}^{d}\to\mathbb{R}$ and $\mathbb{T}^{d}\to\mathbb{R}$, respectively and satisfying certain conditions. It is clear that if $m=0$ (no noise), then this equation might have multiple invariant measures. On the other hand, if $m=\infty$ (noise acts in every direction), then, under certain additional assumptions on $\sigma$, the process $u$ is strong Feller and it can be shown by the classical methods that it has a unique invariant measure [6, Sections 7 and 11]. Thus, it is natural to ask what the smallest number of directions $m$ that have to be perturbed by noise is, so that equation (1.1) still has a unique invariant measure.

Using the ASF method described above, it was shown in [18, Remark 8.22] that if $f$ is a polynomial and $m=3$, then equation (1.1) has a unique invariant measure and is exponentially ergodic. We extend this result and show exponential ergodicity of $u$ already if $m=1$ (that is, when noise acts only in one direction), see Theorem 4.5 and Remark 4.6; we also do not rely on a specific form of $f$. Note that the convergence to the invariant measure is established in the Wasserstein metric; Theorem 4.8 shows that in the case $m<\infty$ this is optimal. Namely, if no additional assumptions are imposed, then, contrary to the case $m=\infty$, the transition probabilities of Markov process (1.1) might not converge to the invariant measure in the total variation metric. Finally, in Theorem 4.9, we show that any two solutions to (1.1) launched with the same noise from different initial conditions converge to each other exponentially fast (synchronization by noise).

The idea that order–preservation helps to obtain better convergence rates of a Markov process is not new; it was used to study interacting particle systems since 1970s, see [25]. However, it is not clear at all how to extend the techniques used in this book beyond the framework of interacting particle models. Order–preserving Markov processes on a general state spaces were considered in [28, 32]. However, the methods developed there rely on the small set condition. Since in the current paper we study processes on a general state space with bad mixing properties, where this condition might not hold, unfortunately the ideas of [25, 28, 32] cannot be applied in our case.

It is interesting to compare our results with [8]. In that paper the authors consider an order–preserving random dynamical system (RDS) with two additional properties: it has a unique invariant measure and it weakly converges to this measure. It is shown there that this implies that any two trajectories of the RDS converge to each other in probability. By contrast, in the current paper we start with an order–preserving Markov process and prove uniqueness of an invariant measure and convergence of transition probabilities.

Our main tool is a new version of the coupling method specifically tailored for order–preserving Markov processes, see the proof of Theorem 2.3. This is combined with an analysis of the relations between stochastic domination and expected distance of random variables, see Section 3. There we continue the study initiated in [4, Proposition 1] and [8, Proposition 2.4]. Note however, that the methods introduced in [4] and [8] cannot be used to get a quantitative bound even in the case where the Markov process has state space $\mathbb{R}$, see Example 3.2. Therefore we apply a new technique.

While we study in detail only the stochastic reaction–diffusion equations on the torus, the strategy developed in this paper should also work in a very similar way for other order–preserving SPDEs including stochastic reaction–diffusion equation on a bounded domain and stochastic porous medium equations. We would like to mention also that after the preprint of our paper became available, our technique was extended and adapted to study regularization by noise for singular SPDEs [10].

The rest of the paper is organized as follows. We present our main results in Section 2. In Section 3 we investigate the relations between stochastic domination and average distance of random variables. Section 4 is devoted to a detailed study of ergodicity of stochastic reaction–diffusion equations. The proofs of the main results are placed in Section 5.

Convention on constants. Throughout the paper $C$ denotes a positive constant whose value may change from line to line.

Acknowledgments. The authors are deeply indebted to Benjamin Gess and Konstantinos Dareiotis for their help, patience and detailed explanations of some parts of the theory of parabolic PDEs. We also would like to thank Máté Gerencsér, Alessandra Lunardi, Lenya Ryzhik, Vladimír Šverák, and Pavlos Tsatsoulis for useful comments and helpful discussions. We are also grateful to two anonymous reviewers for providing helpful comments on an earlier draft of the manuscript. Part of the work on the project has been done during the visit of OB to the Institute of Science and Technology — Austria (IST) and Max-Planck-Institut für Mathematik in den Naturwissenschaften (Leipzig, Germany). OB is very grateful to these institutions for providing excellent working conditions. OB has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No. 683164).

2 A general framework for ergodicity for order–​ preserving Markov processes

We begin by introducing the basic notation. Let $(E,\rho)$ be a Polish space with partial order $\preceq$ such that the set

[TABLE]

is closed. Let $\mathcal{E}=\mathcal{B}(E)$ be the Borel $\sigma$-field. For sets $A,B\in\mathcal{E}$ we will write $A\preceq B$ if for any $a\in A$, $b\in B$ we have $a\preceq b$. Denote by $\mathcal{P}(E)$ the set of all probability measures on $(E,\mathcal{E})$. Let $\{P_{t}\}_{t\geqslant 0}=\{P_{t}(x,A),x\in E,A\in\mathcal{E}\}_{t\geqslant 0}$ be a Markov transition function over $E$ and denote by $\{\mathsf{P}_{x},x\in E\}$ the corresponding Markov family; that is $\mathsf{P}_{x}$ is the law of the Markov process $\{X_{t},t\geqslant 0\}$ with the given transition function and initial condition $X_{0}=x$. The law of $X$ will be understood in the sense of finite–dimensional distributions; that is, we will not rely on the trajectory–wise properties of $X$.

For a measurable function $r:E\times E\to[0,1]$, we consider the corresponding coupling distance $W_{r}:\mathcal{P}(E)\times\mathcal{P}(E)\to\mathbb{R}_{+}$ given by

[TABLE]

where $\mathcal{C}(\mu,\nu)$ is the set of all couplings between $\mu$ and $\nu$, i.e., probability measures on $(E\times E,\mathcal{E}\otimes\mathcal{E})$ with marginals $\mu$ and $\nu$. If $r$ is a lower semicontinuous metric on $E$, then $W_{r}$ is the usual Kantorovich–Wasserstein distance. If $r$ is the discrete metric, i.e., $r(x,y)=\operatorname{\mathbbm{1}}(x\neq y)$, then $W_{r}$ is the total variation distance, which will be denoted further by $d_{TV}$.

Let us now recall the standard definitions related to the partial order $\preceq$; we refer to, e.g., [26, Chapter IV] for a detailed discussion.

Definition 2.1**.**

1. (i)

A function $f\colon E\to\mathbb{R}$ is called increasing if for any $x,y\in E$ with $x\preceq y$ we have $f(x)\leqslant f(y)$. 2. (ii)

Let $\mu,\nu\in\mathcal{P}(E)$ be two probability measures. We say that $\nu$ stochastically dominates $\mu$ and denote it by $\mu\preceq_{st}\nu$ if for any bounded measurable increasing function $f\colon E\to\mathbb{R}$ we have

[TABLE] 3. (iii)

We say that a Markov transition function $\{P_{t}\}_{t\geqslant 0}$ is order preserving if for any $t>0$ and $x,y\in E$ such that $x\preceq y$ we have

[TABLE]

In other words, a Markov transition function is order preserving if it maps bounded increasing functions to bounded increasing functions. Examples of Markov processes with an order preserving transition function include stochastic–reaction diffusion equations, stochastic porous media equations and others, see, e.g., [8].

Remark 2.2**.**

Strassen’s theorem (see, e.g., [26, Theorem IV.2.4]) provides the following coupling definition of stochastic domination, which is equivalent to the one stated above. We have $\mu\preceq_{st}\nu$ if and only if there exist random elements $X,Y\colon\Omega\to E$ such that $\operatorname{Law}(X)=\mu$, $\operatorname{Law}(Y)=\nu$ and $X\preceq Y$.

Now we are ready to present our main results.

Theorem 2.3**.**

Suppose that there exist measurable functions $V\colon E\to[0,\infty)$, $\varphi\colon E\to\mathbb{R}$, $d\colon E\times E\to[0,\infty)$, such that the following conditions hold:

1. (1)

the Markov transition function $(P_{t})_{t\geqslant 0}$ is order–preserving; 2. (2)

the function $V$ is a Lyapunov function, that is, there exist constants $\gamma,K>0$ such that

[TABLE] 3. (3)

if $x,y\in E$ and $x\preceq y$, then $0\leqslant d(x,y)\leqslant\varphi(y)-\varphi(x)$; 4. (4)

for any $x\in E$ we have $M(x):=\sup_{t\geqslant 0}P_{t}\varphi^{2}(x)<\infty$; 5. (5)

there exist sets $A,B\in\mathcal{E}$ and $\varepsilon>0$ such that $A\preceq B$ and for any $x\in\{V\leqslant 4K/\gamma\}$ we have

[TABLE]

Then for any $\theta>0$ there exist constants $C,\lambda>0$ such that for any $x,y\in E$

[TABLE]

Theorem 2.4**.**

Suppose that all conditions of Theorem 2.3 are satisfied. Assume further that

1. (6)

there exists $\delta\in(0,1]$ such that $\rho\leqslant d^{\delta}$, where $\rho$ is the metric on the Polish space $E$. 2. (7)

there exists $K>0$, $\kappa>0$ such that $M(x)^{\kappa}\leqslant K+KV(x)$ for all $x\in E$.

Then the Markov semigroup has a unique invariant measure $\pi$. Further, for any $\theta>0$ there exist constants $C,\lambda>0$ such that for any $x\in E$

[TABLE]

Sketch of the proof of Theorems 2.3 and 2.4.

Here, for the convenience of the reader, we provide just a very brief roadmap of the proof; a complete proof is given in Section 5.

Fix $x,y\in E$, $t>0$. The proof splits into two independent parts. First, we use a new version of the coupling method and utilize conditions (1), (2), (5) to construct random elements $Z^{x},Z^{y},\widetilde{Z}^{x}$ taking values in $E$ with the following properties:

[TABLE]

for some universal constants $C_{1},C_{2}>0$. Second, using the ideas developed in Section 3, we transform the bound (2) into the following inequality:

[TABLE]

where $C_{3},C_{4}>0$ are again some universal constants. It is at this step, where we are using conditions (3) and (4). Since $\operatorname{Law}(Z^{x})=P_{t}(x,\cdot)$ and $\operatorname{Law}(Z^{y})=P_{t}(y,\cdot)$, inequality (2.7) yields (2.4) and (2.5). ∎

Theorems 2.3 and 2.4 provide general sufficient conditions for an order–preserving Markov process to be ergodic. Condition (3) of Theorem 2.3 is actually a condition on the space $(E,\preceq)$ and $d$ rather than on the Markov semigroup. As shown below it is satisfied in many natural situations, including $E=\mathbb{R}$, $E=L_{p}$, $p\in[1,\infty)$, $E=\mathscr{B}^{\alpha}_{p,\infty}$, $\alpha<0$, $p\in[1,\infty)$; the latter stands for the Besov space of regularity $\alpha$ and integrability $p$, see, e.g., [1, Section 2.7 and Proposition 2.93]. Thus, the only additional assumption for exponential ergodicity (apart from the standard Lyapunov and moment-type conditions) is the swap condition (5). It tells that the state space $E$ contains two sets, one preceding the other, and locally uniformly in the initial condition the Markov process has a small chance to be in either of these sets. The following simple example explains why this condition cannot be dropped.

Example 2.5**.**

Introduce the following trivial order on $E$: for $x,y\in E$ we have $x\preceq y$ if and only if $x=y$. It is clear, that for this order any Markov semigroup is order–preserving. We also see that the set $\Gamma$ defined in (2.1) is closed; furthermore, conditions (3) and (4) of Theorem 2.3 trivially hold with $\varphi\equiv 0$. Thus, any Markov process on $E$, that has a Lyapunov function satisfies conditions (1)–(4) of Theorem 2.3. It is well–known that this is not enough for uniqueness of the invariant measure. Thus, the swap condition (5) cannot be omitted.

We also would like to emphasize that the swap condition (5) is very different in nature from the small set condition or other minorization–type conditions, which are imposed within the classical framework, see e.g., [33]. Indeed, a minorization condition guarantees good mixing properties of transition kernels and, in particular, bounds on the total variation distance between the kernels. On the other hand, the swap condition does not yield such bounds since nothing was assumed about mixing on the sets $A$ and $B$. Lemma 5.9 shows how the swap condition can be verified for the stochastic–reaction diffusion equation.

Now let us provide natural examples of spaces $E$ for which condition (3) of Theorem 2.3 holds.

Example 2.6**.**

Let $S$ be a countable set. Put $E=\{0,1\}^{S}$ equipped with the distance $d(x,y):=\sum_{i=1}^{\infty}2^{-i}|x_{i}-y_{i}|$, where $x,y\in\{0,1\}^{S}$ and $x=(x_{1},x_{2},\dots)$, $y=(y_{1},y_{2},\dots)$. This space often appears within the context of interacting particle systems. Consider the following partial order: if $x,y\in\{0,1\}^{S}$, then

[TABLE]

Then condition (3) holds for the function $\varphi(x):=\sum_{i=1}^{\infty}2^{-i}x_{i}$, $x\in\{0,1\}^{S}$.

Example 2.7**.**

Put $E=\mathbb{R}$ equipped with the standard distance, $d(x,y):=|x-y|$, $x,y\in\mathbb{R}$, and consider the standard order $\leqslant$. Then condition (3) holds for the function $\varphi(x):=x$, $x\in\mathbb{R}$.

Example 2.8**.**

Put $E:=L_{p}(D,\mathbb{R})$, where $p\geqslant 1$ and an arbitrary domain $D\subset\mathbb{R}^{n}$, $n\in\mathbb{N}$. Let $\|\cdot\|_{L_{p}}$ be the standard $L_{p}$ norm in this space. Consider the following partial order

[TABLE]

and let $d(x,y):=\|x-y\|_{L_{p}(D)}^{p}$, $x,y\in L_{p}(D,\mathbb{R})$. Then there exists a function $\varphi$ such that the partial order $\preceq_{pos}$ and function $d$ introduced above satisfy condition (3) of Theorem 2.3. Furthermore, the set $\Gamma$ defined in (2.1) is closed. We postpone the proof of this statement to Section 5.

Example 2.9** ([10]).**

This example is due to [10]. Put $E:=\mathscr{B}^{\alpha}_{p,\infty}(D,\mathbb{R})$, where $p\geqslant 1$, $\alpha<0$ and an arbitrary domain $D\subset\mathbb{R}^{n}$, $n\in\mathbb{N}$. Let $\|\cdot\|_{\mathscr{B}^{\alpha}_{p,\infty}}$ be the standard Besov norm in this space, see, e.g., [1, Definition 2.68]. Consider the following partial order

[TABLE]

and let $d(x,y):=\|x-y\|_{\mathscr{B}^{\alpha}_{p,\infty}}^{p}$, $x,y\in\mathscr{B}^{\alpha}_{p,\infty}(D,\mathbb{R})$. Then, as shown in [10, Lemma A.1], there exists a function $\varphi$ such that the partial order $\preceq_{Besov}$ and function $d$ satisfy condition (3) of Theorem 2.3.

3 Stochastic domination and distance between random variables

In this section we explore the connections between stochastic domination and expected distance of random variables, thus continuing the analysis initiated in [8]. Recall that we are given a Polish space $E$ with metric $\rho$. The following statement played a key role in establishing synchronization–by–noise results in [8].

Proposition 3.1** ([8, Proposition 2.4]).**

Let $X_{t}$, $Y_{t}$ be two stochastic processes taking values in $E$ such that $X_{t}(\omega)\preceq Y_{t}(\omega)$ for all $t\geqslant 0$, $\omega\in\Omega$. Suppose that $X_{t}$ and $Y_{t}$ converge weakly as $t\to\infty$ to a random element with the law $\mu$ . Then,

[TABLE]

It is natural to ask whether the above statement can be quantified. More precisely, assume additionally that

[TABLE]

for some rate function $r$ going to [math] as $t\to\infty$. One can ask whether these bounds guarantee a quantitative estimate on $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\rho(X_{t},Y_{t})\wedge 1]$. Quite surprisingly, the answer to this question is negative: without any additional assumptions $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\rho(X_{t},Y_{t})\wedge 1]$ might tend to [math] very slowly even when $r(t)$ converges exponentially fast. This is illustrated by the following example.

Example 3.2**.**

Let $E=\mathbb{R}$ be equipped with the standard order $\leqslant$. Let $(p_{i})_{i\in\mathbb{Z}_{+}}$ be a sequence of positive numbers summing up to $1$. Let $X$ be a random variable which with probability $p_{i}$ is uniformly distributed on $[2^{i},2^{i+1}]$, $i\in\mathbb{Z}_{+}$. Define for $n\in\mathbb{Z}_{+}$

[TABLE]

Then for each $n$ we clearly have $X_{n}\leqslant Y_{n}$ and it is immediate to see that for $\mu:=\operatorname{Law}(X)$ one has

[TABLE]

On the other hand,

[TABLE]

which is much slower than $p_{n}2^{-n}$.

Thus, to quantify bounds in (3.1) we need to impose extra assumptions. This is where condition (3) of Theorem 2.3 pops up.

Lemma 3.3**.**

Assume that condition (3) of Theorem 2.3 is satisfied for measurable functions $\varphi\colon E\to\mathbb{R}$, $d\colon E\times E\to[0,1]$. Suppose further that there exist a function $\psi\colon E\to\mathbb{R}_{+}$ and $k\in(0,1)$ such that

[TABLE]

Let $X,Y$ be random elements $\Omega\to E$ such that $X\preceq Y$ a.s., $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(X)|<\infty$, $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(Y)|<\infty$ and $W_{d}(\operatorname{Law}(X),\operatorname{Law}(Y))\leqslant\varepsilon$ for some $\varepsilon>0$. Then we have

[TABLE]

Proof.

We begin by observing that since $X\preceq Y$ a.s., condition (3) of Theorem 2.3 yields

[TABLE]

Fix $K>1$ and let $\widetilde{X}$, $\widetilde{Y}$ be random elements such that $\operatorname{Law}(\widetilde{X})=\operatorname{Law}(X)$, $\operatorname{Law}(\widetilde{Y})=\operatorname{Law}(Y)$ and $\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptd(\widetilde{X},\widetilde{Y})\leqslant K\varepsilon$ (they exist since $W_{d}(\operatorname{Law}(X),\operatorname{Law}(Y))\leqslant\varepsilon$). Then we continue (3.4) as follows, using the fact that $d$ is bounded by 1:

[TABLE]

Since $K>1$ was arbitrary, this yields the statement of the lemma. ∎

The lemma above will be very useful for obtaining exponential bounds on the synchronization​ –by–noise. On the other hand, to complete the second part of the proof of Theorem 2.3 (see its sketch above in Section 2), we need to solve the opposite problem. More precisely, Lemma 3.3 considers the case when one random element is less than the other everywhere but their laws are different. The following lemma studies the situation where the laws of the random elements coincide, but one is less than the other with probability smaller than $1$.

Lemma 3.4**.**

Assume that condition (3) of Theorem 2.3 is satisfied for measurable functions $\varphi\colon E\to\mathbb{R}$, $d\colon E\times E\to[0,1]$. Let $X,Y,\widetilde{X}$ be random elements $\Omega\to E$ such that $\operatorname{Law}(X)=\operatorname{Law}(\widetilde{X})$ and

[TABLE]

for some $\varepsilon>0$. Then for any $p,q>1$ such that $1/p+1/q=1$ we have

[TABLE]

Again, we see that the statement of the lemma is satisfied only under some extra conditions (one needs condition (3) of Theorem 2.3 and a reasonable bound on $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(X)|^{q}$). The following two examples show that these extra conditions cannot be dropped.

Example 3.5**.**

Let $E=\mathcal{C}([0,1],\mathbb{R})$ be equipped with the sup norm $\|\cdot\|$ and order $\preceq_{pos}$ (defined as in (2.8)). For $x,y\in E$ put $d(x,y):=\|x-y\|\wedge 1$. Let $n\in\mathbb{N}$. For $k\in\{1,2,\ldots,n\}$ let $f_{k}$ be an element of $\mathcal{C}([0,1],\mathbb{R})$ taking values in $[0,1]$ such that

[TABLE]

Then it is easy to see that $f_{k}\preceq_{pos}f_{l}$ whenever $k\leqslant l\leqslant n$. Let $\eta$ be uniformly distributed on the set $\{1,2,\ldots,n\}$. Put $X:=f_{\eta}$ and $Y=\widetilde{X}:=f_{\eta+1}$, where the summation is taken $\mod n$. Then $X$ and $\widetilde{X}$ are bounded, $\operatorname{Law}(X)=\operatorname{Law}(\widetilde{X})$, $\mathsf{P}(X\preceq Y\preceq\widetilde{X})=1-1/n$, but $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\|X-Y\|\wedge 1]=1$.

It is easy to see that in this example condition (3) of Theorem 2.3 does not hold. Indeed, if a function $\varphi$ satisfies this condition, then for any $n\in\mathbb{N}$, $k\in\{1,2,\ldots,n-1\}$ we have $\varphi(f_{k+1})\geqslant\varphi(f_{k})+d(f_{k},f_{k+1})=\varphi(f_{k})+1$. Therefore

[TABLE]

where $\check{1}$ and $\check{0}$ denote elements of $E$ identically equal to $1$ and [math], respectively. Since $n$ was arbitrary, we see that (3.5) implies that $\varphi(\check{1})$ cannot be finite, which is impossible.

Example 3.6**.**

Let $E=\mathbb{R}$ be equipped with the standard order $\leqslant$. For $x,y\in E$ put $d(x,y):=|x-y|\wedge 1$. Let $n\in\mathbb{N}$ and let $X$ be uniformly distributed on the set $\{1,2,\ldots,n\}$. Let $Y=\widetilde{X}=(X+1)\operatorname{\mathbbm{1}}(X<n)+\operatorname{\mathbbm{1}}(X=n)$. Then condition (3) of Theorem 2.3 holds with $\varphi(x)=x$, $\operatorname{Law}(X)=\operatorname{Law}(\widetilde{X})$, $\mathsf{P}(X\preceq Y\preceq\widetilde{X})=1-1/n$, but $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[|X-Y|\wedge 1]=1$.

Thus, we see that in this example condition (3) of Theorem 2.3 holds, but $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(X)|=\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptX=n/2+1/2$ can be arbitrarily large for large $n$.

Proof of Lemma 3.4.

We begin by observing that the function $\varphi$ is increasing, that is $x\preceq y$ implies $\varphi(x)\leqslant\varphi(y)$ for any $x,y\in E$. Fix $p,q>1$ such that $1/p+1/q=1$. Without loss of generality, we can also assume that $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(X)|^{q}<\infty$ (otherwise the statement of the lemma is trivial).

Let $\Lambda:=\{\omega:X(\omega)\preceq Y(\omega)\preceq\widetilde{X}(\omega)\}$. Then, by assumption, $\mathsf{P}(\Lambda)\geqslant 1-\varepsilon$ and we get (recall that $d\leqslant 1$)

[TABLE]

Since $\operatorname{Law}(X)=\operatorname{Law}(\widetilde{X})$, we see that $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt(\varphi(\widetilde{X})-\varphi(X))=0$. Applying the Cauchy–Schwarz inequality, we finally obtain from (3)

[TABLE]

4 Exponential ergodicity of the stochastic reaction–diffusion equation in the hypoelliptic setting

Fix a filtered probability space $(\Omega,\mathcal{F},(\mathcal{F}_{t})_{t\geqslant 0},\mathsf{P})$. We consider the stochastic reaction–diffusion equation (1.1) evolving on a $d$–dimensional torus $\mathbb{T}^{d}$, $d\in\mathbb{N}$. Note that we have equipped this equation with periodic boundary conditions just to simplify the exposition and to emphasize the main ideas. Everything should work in a similar way for equation (1.1) on a bounded domain equipped with Dirichlet boundary conditions.

We are interested in analytically and probabilistically strong solutions to this equation. Recall the following notions.

Definition 4.1** (see, e.g., [9, Definition 1.7]).**

Let $u_{0}\in L_{2}(\mathbb{T}^{d})$. An $L_{2}(\mathbb{T}^{d})$-valued, continuous in time process $u\colon\Omega\times[0,\infty)\to L_{2}(\mathbb{T}^{d})$ is called an analytically generalized strong solution to (1.1) with initial condition $u_{0}$ if $u(0)=u_{0}$ and for any $\varepsilon>0$, $t>\varepsilon$ the following holds:

[TABLE]

Definition 4.2** (see, e.g., [7, Section 6.1], [9, Definition 1.3]).**

Let $u_{0}\in L_{2}(\mathbb{T}^{d})$. An $L_{2}(\mathbb{T}^{d})$-valued, continuous in time process $u\colon\Omega\times[0,\infty)\to L_{2}(\mathbb{T}^{d})$ is called an analytically strong solution to (1.1) with initial condition $u_{0}$ if conditions (4.1) and (4.2) hold for $\varepsilon=0$ and any $t>0$.

Definition 4.3** (see, e.g., [14, Remark 5.2]).**

A solution to (1.1) (strong or generalized strong) is called a probabilistically strong solution if it is adapted to the filtration $(\mathcal{F}_{t})_{t\geqslant 0}$.

In case of ambiguity, a solution to (1.1) with initial condition $x\in L_{2}(\mathbb{T}^{d})$ will be denoted by $u^{x}$. We refer to [9, 7, 14] for a detailed discussion of relations between different notions of a solution to SPDEs.

We will make the following assumption on $f$ and $\sigma_{i}$.

Assumption A**.**

The function $f$ is jointly continuous and locally Lipschitz in the first variable. Moreover, the following condition holds (dissipativity outside of a compact set): there exist $K_{1},K_{2},K_{3}>0$ such that for any $x,y\in\mathbb{R}$, $\xi\in\mathbb{T}^{d}$

[TABLE]

Suppose that

[TABLE]

As shown in Theorem 4.4 below, the condition on $\sigma$ guarantees local existence of strong solution and condition (4.3) prevents blow-up. Conditions (4.3) and (4.4) are satisfied, for example, for $f(x,\xi):=Kx-x^{3}$, $K\in\mathbb{R}$.

Recall the definition of the order $\preceq_{pos}$ in (2.8).

Theorem 4.4**.**

Suppose that Assumption A holds. Then

1. (i)

for any $u_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$ equation (1.1) has a unique analytically and probabilistically strong solution $u$ with initial condition $u_{0}$; furthermore $u$ has a version which is continuous on $[0,\infty)\times\mathbb{T}^{d}$; 2. (ii)

for any $u_{0}\in L_{2}(\mathbb{T}^{d})$ equation (1.1) has a unique analytically generalized strong and probabilistically strong solution $u$ with initial condition $u_{0}$; furthermore $u$ has a version which is continuous on $(0,\infty)\times\mathbb{T}^{d}$; 3. (iii)

this solution of equation (1.1) is a homogeneous Markov process with state space $L_{2}(\mathbb{T}^{d})$; 4. (iv)

there exists a set $\Omega^{\prime}\subset\Omega$ of full measure such that if $x,y\in L_{2}(\mathbb{T}^{d})$ and $x\preceq_{pos}y$, then

[TABLE] 5. (v)

the corresponding Markov transition function is order preserving with respect to the order $\preceq_{pos}$; 6. (vi)

for any $x\in L_{2}(\mathbb{T}^{d})$, $t\geqslant 0$ we have $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt\|u^{x}(t)\|_{L_{2}(\mathbb{T}^{d})}^{2}<\infty$; moreover there exists a constant $C>0$ such that the following energy estimates hold for any $x\in L_{2}(\mathbb{T}^{d})$, $0\leqslant s\leqslant t$

[TABLE]

where the constants $K_{1}$ and $K_{2}$ were defined in (4.3).

The proof of the theorem is given in Section 5.2. Let us make here a couple of remarks about the theorem. Usually in the literature it is assumed additionally that $f$ is a polynomial or $f$ is globally Lipschitz or that at least some growth bounds on $|f(x,\xi)|$ hold, see, e.g., [7, Section 7], [9, Section 4.2], [18, Section 8.3]. Here we established existence and uniqueness of the solutions to (1.1) without these additional restrictions. The main challenge is that the corresponding Nemytskii operator

[TABLE]

is not $L_{2}(\mathbb{T}^{d})\to L_{2}(\mathbb{T}^{d})$ (or at least $L_{p}(\mathbb{T}^{d})\to L_{2}(\mathbb{T}^{d})$ for some large $p$). Therefore it is difficult to apply here any fixed point principle. Furthermore if $u^{n}$ is a sequence of solutions to (1.1) with smooth initial conditions $u_{0}^{n}$, and $u^{n}$ converges to some $u$ in $L_{2}$, then it is not clear at all that this $u$ is a solution to (1.1). To overcome these obstacles, inspired by some ideas from [9], we have extended the corresponding PDE result [27, Propositions 7.3.1] to $L_{2}$ initial conditions.

The fact that for irregular initial data ($u_{0}\in L_{2}(\mathbb{T}^{d})$) equation (1.1) might not have an analytically strong solution is not surprising. Indeed, even for the standard heat equation $\partial_{t}u=\Delta u$ on $\mathbb{T}^{1}$ we have $\|\Delta u(t)\|_{L_{2}(\mathbb{T}^{d})}\leqslant Ct^{-1}\|u(0)\|_{L_{2}(\mathbb{T}^{d})}$ and thus $\int_{0}^{t}\|\Delta u(s)\|_{L_{2}(\mathbb{T}^{d})}\,ds$ might be infinite.

Now let us present the main result of this section and establish ergodicity of the stochastic reaction–diffusion equation in the hypoelliptic setting. By Theorem 4.4(iii) the solution to (1.1) is a Markov process with the state space $L_{2}(\mathbb{T}^{d})$. Let $(P_{t})_{t\geqslant 0}$ be the family of Markov transition probabilities associated with this process. We introduce the following condition.

Assumption B**.**

There exist $\varepsilon>0$ and $\lambda_{k}\in\mathbb{R}$, $k=1,\ldots,m$ such that

[TABLE]

Theorem 4.5**.**

Suppose that Assumptions A and B hold. Then the stochastic reaction–diffusion equation (1.1) has a unique invariant measure $\pi$. Furthermore, there exist $C>0$, $\lambda>0$ such that for all $x\in L_{2}(\mathbb{T}^{d})$ we have

[TABLE]

Remark 4.6**.**

Assumption B is satisfied for example if $m\geqslant 1$ and $\sigma_{1}(\xi)>\varepsilon$ for all $\xi\in\mathbb{T}^{d}$.

Remark 4.7**.**

Assumption B is imposed only to guarantee the swap condition, see Lemma 5.9. If the solution to (1.1) satisfies (5.43), then only Assumptions A is needed in Theorem 4.5.

As mentioned in the introduction, it is clear that if $m=0$, then equation (1.1) might have multiple invariant measures. On the other hand, if $m=\infty$, equation (1.1) has a unique invariant measure (of course, under certain assumption on the rate of decay of $\|\sigma_{k}\|_{L_{2}}$ as $k\to\infty$) [6, Sections 7 and 11]. The first result showing uniqueness of the invariant measure for finite $m$ (that is when noise acts not in all directions) is due to Hairer [12], who showed it for $m$ large enough. This was later improved by Hairer and Mattingly [18, Remark 8.22], where uniqueness of the invariant measure and exponential ergodicity is proven if $d=1$ and $m=3$. Remark 4.6 refines this result and shows that if the noise acts only in one direction and the corresponding $\sigma_{1}$ is positive, then the solution to equation (1.1) is exponentially ergodic. Note that contrary to [18, Remark 8.22] we also do not rely here on the specific polynomial form of the drift $f$ and do not assume that the space dimension $d\leqslant 3$ as in [18, Assumption RD.2].

The next result shows that, in general, the convergence of transition probabilities to the invariant measure in (4.8) might not take place in the total variation metric, and thus the $W_{\|\cdot\|_{L_{2}}\wedge 1}$ metric can not be replaced in (4.8) by the $d_{TV}$ metric. As discussed in the introduction, this happens due to the fact that the stochastic system do not get “enough” noise. A related result for the case when the drift $f$ is linear can be found in [15, Proposition 9.7].

Theorem 4.8**.**

Suppose that Assumption A holds for a function $f$ which does not depend on space (so $f=f(x)$, $x\in\mathbb{R}$). Then there exists a function $\sigma$ satisfying additionally Assumption B and an initial condition $x\in L_{2}(\mathbb{T}^{d})$, such that

[TABLE]

for any $t\geqslant 0$.

The proof of Theorem 4.8 is given in Section 5.2.

Proof of Theorem 4.5.

Let us verify that all conditions of Theorems 2.3 and 2.4 are satisfied. Let $E=L_{2}(\mathbb{T}^{d})$ with partial order $\preceq_{pos}$ (recall its definition in (2.8)). Set

[TABLE]

The Markov semigroup corresponding to (1.1) is order–preserving with respect to $\preceq_{pos}$ by Theorem 4.4(v) and thus condition (1) of Theorem 2.3 holds. Condition (2) of Theorem 2.3 is satisfied thanks to energy bound (4.6). It follows from Example 2.8 and the definition of $\varphi$ that condition (3) of Theorem 2.3 is met. Condition (4) follows from (4.6) with $M(x):=2\|x\|_{L_{2}(\mathbb{T}^{d})}^{4}+C$, $x\in L_{2}(\mathbb{T}^{d})$. By Lemma 5.9 condition (5) holds. Conditions (6) and (7) hold trivially with $\delta=1/2$ and $\kappa=1/2$, respectively. Thus, all conditions of Theorems 2.3 and 2.4 are satisfied. Therefore the process $u$ has a unique invariant measure and (4.8) holds. ∎

Finally, we establish exponentially fast synchronization by noise for solutions to (1.1): that is, we will prove that any two solutions with initial conditions $x,y\in L_{2}(\mathbb{T}^{d})$ launched with the same noise will converge to each other in probability exponentially fast.

Theorem 4.9**.**

Suppose that Assumptions A and B hold. Then there exist $C>0$, $\lambda>0$ such that for any $x,y\in L_{2}(\mathbb{T}^{d})$ we have

[TABLE]

Note that the fact that $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt[\|u^{x}(t)-u^{y}(t)\|_{L_{2}(\mathbb{T}^{d})}\wedge 1]\to 0$ as $t\to\infty$ follows immediately from Theorem 4.5 and [8, Proposition 2.4]. Unfortunately, as discussed in Section 3, the techniques developed in [8] do not provide a quantitative bound on convergence rate. Therefore to prove this theorem we use the toolkit developed in Section 3.

Proof of Theorem 4.9.

Fix arbitrary $x_{0},y_{0}\in L_{2}(\mathbb{T}^{d})$, $t\geqslant 0$. Define

[TABLE]

Put $X:=u^{z_{-}}(t)$, $Y:=u^{z_{+}}(t)$. Then by the order–preserving property (Theorem 4.4(iv)) we have

[TABLE]

Therefore

[TABLE]

Now let us check that all the conditions of Lemma 3.3 hold. Let $E=L_{2}(\mathbb{T}^{d})$ with partial order $\preceq_{pos}$. Set $d(x,y):=\|x-y\|_{L_{2}}^{2}\wedge 1$, $x,y\in L_{2}(\mathbb{T}^{d})$, and define $\varphi$ as in (4). Then it follows from Example 2.8 that condition (3) of Theorem 2.3 is satisfied. Furthermore, if $d(x,y)=1$, then we have

[TABLE]

If $d(x,y)<1$, then one has

[TABLE]

where in the second inequality we used the following bound

[TABLE]

Combining (4.12) and (4), we get

[TABLE]

where we put $\psi(x):=4\sqrt{2}(1+\|x\|^{2}_{L_{2}})$, $x\in L_{2}(\mathbb{T}^{d})$. Thus, bound (3.2) holds.

It is also clear that by definition we have $X\preceq_{pos}Y$. Further, by Theorem 4.4(vi) we have

[TABLE]

and, similarly, $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt|\varphi(Y)|<\infty$. Finally, by Theorem 4.5 there exist $C,\lambda>0$ such that

[TABLE]

Thus, all the conditions of Lemma 3.3 are met. Therefore inequality (3.3) and energy bound (4.7) yield

[TABLE]

This bound allows us to conclude.

[TABLE]

Taking into account (4.11), we obtain (4.10). ∎

5 Proofs

5.1 Proofs of the results of Section 2

To prove Theorem 2.3 we need a couple of lemmas. The first lemma is quite standard and is in the spirit of [31, Section 15.2]. The lemma deals with the connection between the finiteness of the exponential moment of the first return time of a Markov process to a set and the existence of a Lyapunov function. However we were not able to find in the literature precisely this statement (we found only a number of related ones). Thus, for the convenience of the reader and for the sake of completeness we provide the full proof of the lemma.

Let us consider a measurable space $(\widetilde{E},\widetilde{\mathcal{E}})$ and a Markov kernel $Q$ on it. Let $Z=(Z_{n})_{n\in\mathbb{Z}_{+}}$ be a Markov process with transition kernel $Q$. For a set $A\in\widetilde{\mathcal{E}}$ introduce the first return time to $A$

[TABLE]

Lemma 5.1**.**

Suppose there exists a measurable function $\mathcal{V}\colon\widetilde{E}\to[1,\infty)$ such that for some $\lambda\in(0,1)$, $K>0$

[TABLE]

Fix

[TABLE]

and put $r:=(\lambda+K/M)^{-1}$.

1. (i)

We have

[TABLE] 2. (ii)

Further, suppose that for a set $A\in\widetilde{\mathcal{E}}$ we have for some $\varepsilon>0$

[TABLE]

Then, for any $1<l<r^{\frac{|\log(1-\varepsilon)|\wedge\log M}{2\log M}}$ there exists a constant $C=C(\lambda,l,K,M,\varepsilon)>0$ such that for any $x\in\widetilde{E}$

[TABLE]

Proof.

(i). The proof uses the standard argument. First note, that it follows from the Lyapunov condition (5.1) and the definition of $r$ that

[TABLE]

Put now for $n\in\mathbb{Z}_{+}$

[TABLE]

Then by Dynkin’s formula for discrete time Markov chains (see, e.g., [31, Theorem 11.3.1]) we have for any $x\in\widetilde{E}$

[TABLE]

If $x\in\widetilde{E}\setminus\{\mathcal{V}\leqslant M\}$, then by definition for any $i\in[1,\tau^{(n)}]$ we have $Z_{i-1}\in\widetilde{E}\setminus\{\mathcal{V}\leqslant M\}$. Therefore, (5.7) implies for any $i\in[1,n]$

[TABLE]

Combining this with (5.8), we get

[TABLE]

which by Fatou’s lemma yields (5.3).

If $x\in\{\mathcal{V}\leqslant M\}$, then by the Markov property and bound (5.3),

[TABLE]

which is (5.4).

(ii). Fix $l\in(1,r^{\frac{|\log(1-\varepsilon)|\wedge\log M}{2\log M}})$. For $n\in\mathbb{Z}_{+}$ denote by $\tau^{M}(n)$ the $n$th return time to the set $\{\mathcal{V}\leqslant M\}$:

[TABLE]

Introduce random variables

[TABLE]

We see that the event $\{I_{n}=1\}$ corresponds to the reaching of the desired set $A$ after $n$th visit to the set $\{\mathcal{V}\leqslant M\}$. We have for any $x\in\widetilde{E}$

[TABLE]

Define

[TABLE]

Note that by the strong Markov property and the definition above, for any $x\in\widetilde{E}$

[TABLE]

where in the last inequality we used (5.3) and the fact that $l<r$. Applying the Cauchy–Schwarz inequality and (5.4), we deduce

[TABLE]

where the last inequality follows from the definitions of $l$ and $r$. Thus, (5.11) together with (5.1), (5.1), and the obvious bound $(\mathcal{V}(x)\vee M)\leqslant M\mathcal{V}(x)$ yields

[TABLE]

which implies (5.6). ∎

The next lemma explains why conditions (3) and (4) are imposed in Theorem 2.3.

Lemma 5.2**.**

Assume that condition (3) of Theorem 2.3 is satisfied for measurable functions $\varphi\colon E\to\mathbb{R}$, $d\colon E\times E\to[0,1]$. Let $X,\,Y,\,\widetilde{X},\,\widetilde{Y}$ be random elements $\Omega\to E$ such that $\operatorname{Law}(X)=\operatorname{Law}(\widetilde{X})$, $\operatorname{Law}(Y)=\operatorname{Law}(\widetilde{Y})$ and

[TABLE]

for some $\varepsilon_{1},\varepsilon_{2}>0$. Then

[TABLE]

Proof.

It follows from the gluing lemma (see, e.g., [38, p. 23]) that there exist random elements $Z_{1},Z_{2},Z_{3}$ such that $\operatorname{Law}(Z_{1},Z_{2})=\operatorname{Law}(X,Y)$ and $\operatorname{Law}(Z_{2},Z_{3})=\operatorname{Law}(\widetilde{Y},\widetilde{X})$. Clearly, we also have $\operatorname{Law}(Z_{1})=\operatorname{Law}(Z_{3})$.

Let $\Lambda$ be the set $\{\omega:Z_{1}(\omega)\preceq Z_{2}(\omega)\preceq Z_{3}(\omega)\}$. Then by above

[TABLE]

Thus, we immediately get from Lemma 3.4 that

[TABLE]

Proof of Theorem 2.3.

The proof of the theorem consists of several steps.

Step 1. In this step we fix $x,y\in E$ and $t>0$. Let $\{X^{x}(s),s\geqslant 0\}$ and $\{X^{y}(s),s\geqslant 0\}$ be independent Markov processes with the laws $\mathsf{P}_{x}$ and $\mathsf{P}_{y}$, correspondingly. Let $\mathcal{F}_{s}:=\sigma(X^{x}(r),X^{y}(r),r\leqslant s)$, $s\leqslant t$, be the natural filtration. Introduce stopping times

[TABLE]

Note that by definition these stopping times $\tau_{x\preceq y}$ and $\tau_{y\preceq x}$ take only countably many values.

Let us now extend the state space $E$ and add an additional element denoted by $\diamondsuit$. We assume that $\diamondsuit\preceq\diamondsuit$ and do not impose any further partial order relations between $\diamondsuit$ and other elements of $E$.

Introduce the following random elements:

[TABLE]

We claim now that $\operatorname{Law}(\eta^{x})\preceq_{st}\operatorname{Law}(\eta^{y})$. Indeed, let $f\colon E\cup\{\diamondsuit\}\to\mathbb{R}$ be an arbitrary bounded measurable increasing function. Then

[TABLE]

Note that for any $i=0,..,\lfloor t\rfloor$ we have

[TABLE]

Recall that the kernel $P_{t-i}$ is order preserving and thus for any $z_{1},z_{2}\in E$ such that $z_{1}\preceq z_{2}$ we have

[TABLE]

Since on the set $\{\tau_{x\preceq y}=i\}$ we have $X^{x}(i)\preceq X^{y}(i)$, we can continue (5.1) in the following way:

[TABLE]

Combining this with (5.12), we finally deduce

[TABLE]

Since $f$ was an arbitrary bounded measurable increasing function, we see that indeed $\operatorname{Law}(\eta^{x})\preceq_{st}\operatorname{Law}(\eta^{y})$.

Step 2. We have shown that $\operatorname{Law}(\eta^{x})\preceq_{st}\operatorname{Law}(\eta^{y})$. Therefore, by the Strassen theorem (see, e.g., [26, Theorem IV.2.4]) there exist random variables $\xi^{x}$ and $\xi^{y}$ such that $\operatorname{Law}(\xi^{x})=\operatorname{Law}(\eta^{x})$, $\operatorname{Law}(\xi^{y})=\operatorname{Law}(\eta^{y})$ and $\xi^{x}\preceq\xi^{y}$ a.s.

It follows from the construction, that $\mathsf{P}(\eta^{x}\neq X^{x}(t))\leqslant\mathsf{P}(\tau_{x\preceq y}>t)$ and $\mathsf{P}(\eta^{y}\neq X^{y}(t))\leqslant\mathsf{P}(\tau_{x\preceq y}>t)$. Now we will apply twice the gluing lemma (see, e.g., [38, p. 23]). First, we apply it to the pairs $(X^{x}(t),\eta^{x})$ and $(\xi^{x},\xi^{y})$ and construct two random variables $Y^{x}$ and $Y^{y}$ such that $\operatorname{Law}(Y^{x})=\operatorname{Law}(X^{x}(t))$, $\operatorname{Law}(Y^{y})=\operatorname{Law}(\xi^{y})=\operatorname{Law}(\eta^{y})$ and $\mathsf{P}(Y^{x}\preceq Y^{y})\geqslant 1-\mathsf{P}(\tau_{x\preceq y}>t)$. Then, we apply the gluing lemma to the pairs $(Y^{x},Y^{y})$ and $(\eta^{y},X^{y})$. We deduce that there exist random variables $Z^{x}$ and $Z^{y}$ such that

[TABLE]

In a similar way, we construct another pair of random variables $\widetilde{Z}^{x}$ and $\widetilde{Z}^{y}$ with the following properties

[TABLE]

Step 3. Now it follows from Step 5.1 and Lemma 5.2 that

[TABLE]

Thus everything boils down to bounding the probabilities $\mathsf{P}(t\leqslant\tau_{y\preceq x})$ and \mathsf{P}\bigl{(}t\leqslant\tau_{x\preceq y}). We will do it using Lemma 5.1.

Step 4. First we note that, clearly, $\{(x_{1},x_{2})\in E\times E\colon x_{1}\preceq x_{2}\}\supset A\times B$ (recall the definitions of the sets $A$ and $B$ in condition (5) of the theorem). Therefore,

[TABLE]

Now we apply Lemma 5.1(ii) to the state space $\widetilde{E}:=E\times E$, the kernel $Q$ on it defined by

[TABLE]

the set $A\times B$, and the Lyapunov function $\mathcal{V}(x_{1},x_{2}):=1+V(x_{1})+V(x_{2})$, $x_{1},x_{2}\in E$. It follows from (2.2) and the Gronwall inequality that

[TABLE]

Therefore, condition (5.1) is met. Take $M=4K/\gamma+1$. It is clear that this choice of $M$ satisfies (5.2). Thanks to assumption (2.3) and the definitions of $\mathcal{V}$ and the sets $A$, $B$, we see that condition (5.5) holds for the set $A\times B$ in place of $A$ and the positive constant $M$ chosen above. Therefore all conditions of Lemma 5.1(ii) are satisfied and, thus there exist constants $C>0$, $\lambda>0$ that do not depend on $x,y$ such that

[TABLE]

This, combined with (5.15) and the Chebyshev inequality implies

[TABLE]

Using exactly the same argument, we also get that

[TABLE]

Substituting these bounds into (5.1), we finally deduce

[TABLE]

Taking into account that obviously $W_{d\wedge 1}(P_{t}(x,\cdot),P_{t}(y,\cdot))\leqslant 1$, we obtain (2.4). ∎

Proof of Theorem 2.4.

We begin by observing that by Jensen’s inequality and condition (6) of the theorem we have for any measures $\mu,\nu\in\mathcal{P}(E)$

[TABLE]

We apply now Theorem 2.3 with $\theta:=\frac{\kappa}{\delta(1+\kappa)}$. Taking into account (5.16), we obtain that there exist $C>0$, $\lambda>0$ such that for any $x,y\in E$ we have

[TABLE]

where we have also used condition (7) of the theorem.

From here the proof follows the standard line of argument. Note that by [21, Theorem 17.24] for any fixed $t\geqslant 0$, the mapping $x\mapsto P_{t}(x,\cdot)$ is measurable. Therefore by [38, Theorem 4.8 and Corollary 5.22] the mapping $(x,y)\mapsto W_{\rho\wedge 1}(P_{t}(x,\cdot),P_{t}(y,\cdot))$ is measurable and the function $W_{\rho\wedge 1}$ is convex in both arguments. Therefore we immediately get from (5.17) that for any $\mu,\nu\in\mathcal{P}(E)$

[TABLE]

Now we are ready to prove the existence of an invariant measure for $\{P_{t},t\geqslant 0\}$. Fix any $x\in E$. We make use of (5.18) to obtain that for any $t,s>0$

[TABLE]

Thus the sequence of measures $\{P_{t}(x,\cdot)\}$ is Cauchy in the space $(\mathcal{P}(E),W_{\rho\wedge 1})$. Since this space is complete (see, e.g., [38, Theorem 6.18]), there exists a measure $\pi\in\mathcal{P}(E)$ such that $W_{\rho\wedge 1}(P_{t}(x,\cdot),\pi)\to 0$ as $t\to\infty$. Hence for any $y\in E$ we have

[TABLE]

Therefore for any $t\geqslant 0$ by the monotone convergence theorem we deduce

[TABLE]

where we have used again that the function $W_{\rho\wedge 1}$ is convex. Therefore the measure $\pi$ is invariant for $\{P_{t},t\geqslant 0\}$.

Now let us show the uniqueness of the invariant measure. First, note that if $\pi$ is an arbitrary invariant measure for $\{P_{t},t\geqslant 0\}$, then $\pi(V)<\infty$, see [13, Proposition 4.24]. Thus if $\pi_{1}$ and $\pi_{2}$ are two invariant measures, then by (5.18) we have for any $t\geqslant 0$

[TABLE]

Taking the limit in the right–hand side of the above inequality we get $\pi_{1}=\pi_{2}$.

Finally, bound (2.5) follows directly from (5.18) and the fact that $\pi(V)<\infty$. ∎

Proof of Example 2.8.

First let us show that the set $\Gamma$ defined in (2.1) is closed. Let $(x_{n})_{n\in\mathbb{Z}_{+}}$, $(y_{n})_{n\in\mathbb{Z}_{+}}$ be two sequences of elements of $L_{p}(D,\mathbb{R})$ such that $x_{n}\preceq_{pos}y_{n}$ and $\|x_{n}-x\|_{L_{p}(D)}\to 0$, $\|y_{n}-y\|_{L_{p}(D)}$ as $n\to\infty$. We claim now that $x\preceq_{pos}y$.

Indeed, by passing to an appropriate subsequence $(n_{k})$ we see that $x_{n_{k}}\to x$ almost everywhere (a.e.) and $y_{n_{k}}\to y$ a.e. as $n\to\infty$. Since for each $k$ we have $x_{n_{k}}\leqslant y_{n_{k}}$ a.e., it is now easily seen that $x\leqslant y$ a.e. and thus $x\preceq_{pos}y$. Therefore, the set $\Gamma$ is closed.

Now let us introduce the function $\varphi$. Put

[TABLE]

Let us verify that the partial order $\preceq_{pos}$ and functions $\varphi$ and $d$ satisfy condition (3) of Theorem 2.3. Our original proof of this was a bit long; the following simple proof was suggested to us by the referee, to whom we are very grateful for this.

Fix any $x,y\in L_{p}(D,\mathbb{R})$ such that $x\preceq_{pos}y$. Note that for any real numbers $a<b$, $p\geqslant 1$ one has

[TABLE]

Using this inequality, we get

[TABLE]

and thus condition (3) of Theorem 2.3 holds. ∎

5.2 Proofs of the results of Section 4

We begin with the following auxiliary statement which provides Gronwall-type bounds that will be very useful in the sequel.

Lemma 5.3**.**

Let $C_{1},C_{2},C_{3},C_{4}$ be positive constants. Let $X$ be a continuous $\mathcal{F}_{t}$-adapted process taking values in $[0,\infty)$ such that $X(0)=x$ and

[TABLE]

where $M$ is a continuous local $\mathcal{F}_{t}$-martingale with $M(0)=0$ and

[TABLE]

Then for any $t\geqslant 0$ we have $\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptX(t)<+\infty$ and the following bounds hold:

[TABLE]

where $C=C(C_{1},C_{2},C_{3},C_{4})>0$ is some constant independent of $t$ and $x$.

Proof.

We begin with the proof of (5.21). For $n\in\mathbb{N}$ introduce the stopping time $\tau_{n}:=\inf\{t\geqslant 0:|M(t)|\geqslant n\}$. Then it follows from (5.19) that for any $n\in\mathbb{N}$, $t\geqslant 0$

[TABLE]

Using Fatou’s lemma, we deduce that for any $t\geqslant 0$ we have $\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptX(t)\leqslant x+C_{2}t$ and thus $\hskip 0.58122pt\mathsf{E}\hskip 0.38751ptX(t)<\infty$. This and (5.20) imply that for any $t\geqslant 0$ we have $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt\langle M\rangle_{t}\leqslant C_{4}t+C_{3}xt+C_{2}C_{3}t^{2}/2<\infty$. Therefore $M$ is a martingale and (5.21) follows immediately from (5.19).

To establish (5.22), we introduce a process $Y$, which is a solution to the following equation

[TABLE]

We claim that for any $t\geqslant 0$ we have $X(t)\leqslant Y(t)$. Indeed, assume the contrary and suppose that for some $T_{0}>0$ we have $X(T_{0})>Y(T_{0})$. Then, arguing as in [20, Proof of Proposition 9.2], we introduce

[TABLE]

Since the processes $X$ and $Y$ are continuous we have $S_{0}<T_{0}$ and $X(t)>Y(t)$ for $t\in(S_{0},T_{0})$. Then (5.19) and (5.23) imply

[TABLE]

which contradicts the fact that $X(T_{0})>Y(T_{0})$. Therefore such $T_{0}$ does not exist and $X(t)\leqslant Y(t)$ for any $t\geqslant 0$.

Note now that by Ito’s formula and (5.20)

[TABLE]

where $C_{5}:=(2C_{2}+C_{3})^{2}/(4C_{1})+C_{4}$ and $N$ is a continuous local martingale.

Consider a stopping time $\sigma_{n}:=\inf\{t\geqslant 0:|N(t)|\geqslant n\}$. Then, using (5.2), we get for any $0\leqslant s\leqslant t$

[TABLE]

Therefore, the Gronwall inequality yields

[TABLE]

The proof of (5.22) is completed by passing to the limit in the above inequality using Fatou’s lemma. ∎

In all the remaining theorems and lemmas of this section we will assume that Assumption A is satisfied. To establish Theorem 4.4 we split the process $u$ into the following two parts. The first part, denoted by $w$, is the stochastic convolution, that is a unique analytically and probabilistically strong solution of the following equation

[TABLE]

with the initial condition $w(0)=0$.

Lemma 5.4**.**

The function $w$ is well–defined. Furthermore, there exists a set $\Omega^{\prime}\subset\Omega$ such that $\mathsf{P}(\Omega^{\prime})=1$ and for any $\varepsilon>0$, $T>0$, $\omega\in\Omega^{\prime}$ there exists $C=C(\varepsilon,T,\omega)$ such that

[TABLE]

Proof.

Existence and uniqueness of a strong solution of (5.25) follows from [7, Example 5.39] and assumption (4.5). Bound (5.26) follows from [7, Theorem 5.15(ii)]. ∎

For $T>0$ denote now

[TABLE]

By Lemma 5.4, for any $\omega\in\Omega^{\prime}$, $T>0$ we have $M_{T}(\omega)<\infty$.

The second ingredient of decomposition of $u$ is much more tricky. We fix now arbitrary $\omega\in\Omega^{\prime}$, $T>0$ and consider the (deterministic) PDE

[TABLE]

with the initial condition $v(0)=v_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$. We begin with the case when the initial condition is smooth enough.

Lemma 5.5**.**

For any $v_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$ equation (5.28) has a unique analytically strong solution $v$ with the initial condition $v_{0}$. Furthermore $v\in\mathcal{C}([0,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$ and there exists $C>0$ such that for any $t\in[0,T]$

[TABLE]

where $C_{f}:=\sup_{\begin{subarray}{c}x\geqslant 0\\ \eta\in\mathbb{T}^{d}\end{subarray}}(\operatorname{sign}(x)f(x,\eta))\vee 0$.

Proof.

We begin with the uniqueness part. Suppose that $v$ and $\bar{v}$ are two strong solutions to equation (5.28) with the same initial condition $v_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$. Then, by the chain rule

[TABLE]

where the last inequality follows from assumption (4.4). Since $v(0)=\bar{v}(0)$, an application of the Gronwall lemma immediately implies that $v(t)=\bar{v}(t)$ for any $t\in[0,T]$.

To show existence of a strong solution to (5.28) we fix the initial condition $v_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$ and $T>0$. Note that the function

[TABLE]

is locally Lipschitz in $x$, Hölder in $t$ and continuous in $\xi$ thanks to (5.26) and Assumption A. Therefore it satisfies [27, Inequality (7.3.2)]. Thus, by [27, Propositions 7.3.1.i, 7.3.1.ii, 7.1.10.iii] equation (5.28) has a local solution on some interval $[0,\delta]$, $\delta>0$ and this solution is in $\mathcal{C}([0,\delta],\mathcal{C}^{2}(\mathbb{T}^{d}))$. Let us now show that equation (5.28) has a global solution.

It follows from (4.3) that there exist $x_{-}\leqslant x_{+}\in\mathbb{R}$ such that $f(x,\xi)\geqslant 0$ for any $x\leqslant x_{-}$, $\xi\in\mathbb{T}^{d}$ and $f(x,\xi)\leqslant 0$ for any $x\geqslant x_{+}$, $\xi\in\mathbb{T}^{d}$. Consider the following PDE

[TABLE]

with the initial condition $\psi(0)=(x_{+}+M_{T})\vee\sup_{\eta}v_{0}(\eta)$ (recall the definition of $M_{T}$ in (5.27)). By above, the constant function

[TABLE]

solves this equation on $[0,T]$. On the other hand, since $v(0,\xi)\leqslant\psi(0,\xi)$ and $f(x+w(t,\xi,\omega),\xi)\leqslant f(x+w(t,\xi,\omega),\xi)\vee 0$, the comparison theorem for reaction–diffusion PDEs [36, Theorem 10.1] implies that for every interval $[0,\delta]$ where equation (5.28) has a local solution we have

[TABLE]

By a similar argument,

[TABLE]

Thus, for any interval $[0,\delta]$ where equation (5.28) has a local solution, this solution is uniformly bounded by a constant which does not depend on $\delta$. Therefore, by [27, Proposition 7.3.1.v] equation (5.28) has a global solution on $[0,T]$. By [27, Proposition 7.1.10.iii] this solution is in $\mathcal{C}([0,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$.

Finally, to obtain the desired bound on $\|v(t)\|_{L_{\infty}(\mathbb{T}^{d})}$ in terms of $\|v_{0}\|_{L_{2}(\mathbb{T}^{d})}$ rather than $\|v_{0}\|_{L_{\infty}(\mathbb{T}^{d})}$ we fix $t\in[0,T]$ and consider again PDE (5.32) with the initial condition $\psi(0,\xi):=(v_{0}(\xi)\vee 0)+M_{t}$. By the comparison principle, we have

[TABLE]

Since $\psi(0,\xi)\geqslant M_{t}$ for all $\xi\in\mathbb{T}^{d}$ and the drift is nonnegative, it is immediate to see that $\psi(s,\xi)\geqslant M_{t}$ for all $s\in[0,t]$, $\xi\in\mathbb{T}^{d}$. Therefore, taking into account (5.27), we get that for all $s\in[0,t]$, $\xi\in\mathbb{T}^{d}$

[TABLE]

Introduce now $\psi_{+}$, which is the strong solution of the following PDE

[TABLE]

with the initial condition $\psi_{+}(0):=\psi(0)$. Using (5.34), we see that the comparison principle implies that

[TABLE]

Let $p_{s}$ be the heat kernel on the torus $\mathbb{T}^{d}$. Then it is straightforward to see that for any $\xi\in\mathbb{T}^{d}$

[TABLE]

Combining this with (5.33) and (5.35), we deduce

[TABLE]

By a similar argument, we get

[TABLE]

This yields (5.29). ∎

Now let us move on to less regular initial data.

Lemma 5.6**.**

For any $v_{0}\in L_{\infty}(\mathbb{T}^{d})$ equation (5.28) has a unique analytically generalized strong solution $v$ with initial condition $v_{0}$. Furthermore, we have $v\in\mathcal{C}((0,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$.

Proof.

We begin with the uniqueness part. Suppose that $v$ and $\bar{v}$ are two analytically generalized strong solutions to equation (1.1) with the same initial condition $v_{0}\in L_{\infty}(\mathbb{T}^{d})$. Fix arbitrary $t>0$. Arguing as in the proof of Lemma 5.5, we see that by the Gronwall lemma there exists $C=C(t)>0$ such that for any $\delta\in(0,t)$ we have

[TABLE]

By definition, $v(\delta)\to v(0)$ and $\bar{v}(\delta)\to\bar{v}(0)$ in $L_{2}(\mathbb{T}^{d})$ as $\delta\to 0$. Since $v(0)=\bar{v}(0)$, by passing to the limit as $\delta\to 0$ in (5.36), we deduce that $v(t)=\bar{v}(t)$.

Note again that the function defined in (5.31) is locally Lipschitz in $x$, Hölder in $t$ and continuous in $\xi$. Therefore, by [27, Proposition 7.3.1.i] there exists $\delta>0$ such that on the interval $[0,\delta]$ equation (5.28) has an analytically generalized strong solution $v$ with initial condition $v_{0}$ and $v\in\mathcal{C}((0,\delta],\mathcal{C}^{2}(\mathbb{T}^{d}))$. Since $v(\delta)\in\mathcal{C}^{2}(\mathbb{T}^{d})$, by Lemma 5.5 we can construct an analytically strong solution $v$ on $[\delta,T]$ with the initial condition $v(\delta)$ and $v\in\mathcal{C}([\delta,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$. By gluing these two solutions together we get an analytically generalized strong solution $v\in\mathcal{C}((0,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$. ∎

To consider even less regular initial data (recall that we are interested in the initial conditions from $L_{2}(\mathbb{T}^{d})$) we need the following lemma about approximations of solutions to (5.28).

Lemma 5.7**.**

Let $(v^{n}_{0})_{n\in\mathbb{Z}_{+}}$ be a sequence of $\mathcal{C}^{2}(\mathbb{T}^{d})$ functions converging in $L_{2}$ to $v_{0}\in L_{2}(\mathbb{T}^{d})$. Let $v^{n}$ be the analytically strong solution of (5.28) with the initial condition $v^{n}_{0}$. Then

1. (i)

$v^{n}$* converges in $\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$ to some $v\in\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$;* 2. (ii)

for each $t\in(0,T]$ we have $v(t)\in L_{\infty}(\mathbb{T}^{d})$; 3. (iii)

if $\bar{v}$ is an analytically generalized strong solution of (5.28) with the initial condition $v_{0}$, then for each $t\in[0,T]$ we have $\bar{v}(t,\xi)=v(t,\xi)$ for almost all $\xi\in\mathbb{T}^{d}$. .

Proof.

(i). Let $n,m\in\mathbb{Z}_{+}$. Then arguing as in the derivation of (5.2) we get by the chain rule and assumption (4.4)

[TABLE]

By the Gronwall lemma, this implies that

[TABLE]

Now the statement of the lemma follows immediately from the completeness of the space $\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$ ([21, Theorem I.4.19]).

(ii). Fix $t\in(0,T]$. By Lemma 5.5, we have the following uniform over $n\geqslant N$ bound, where $N$ is large enough:

[TABLE]

Since $v^{n}(t)$ converges to $v(t)$ in $L_{2}$, bound (5.37) implies $v(t)\in L_{\infty}(\mathbb{T}^{d})$.

(iii). Let $\bar{v}$ be an analytically generalized strong solution of (5.28) with the initial condition $v_{0}$. Then arguing as in part (i) of the lemma, we get for any $\delta>0$

[TABLE]

By passing to the limit as $\delta\to 0$ and using the fact that by definition $\bar{v}(\delta)\to v_{0}$ in $L_{2}$, we get

[TABLE]

However, by part (i) of the lemma

[TABLE]

This implies that $\bar{v}(t)=v(t)$ as elements of $L_{2}(\mathbb{T}^{d})$ for any $t\in[0,T]$. ∎

The next lemma establishes existence of solutions to equation (5.28) with $L_{2}(\mathbb{T}^{d})$ initial data. It relies on Lemma 5.7 and extends [27, Propositions 7.3.1].

Lemma 5.8**.**

For any $v_{0}\in L_{2}(\mathbb{T}^{d})$ equation (5.28) has a unique analytically generalized strong solution $v$ with the initial condition $v_{0}$. Furthermore, we have $v\in\mathcal{C}((0,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$.

Proof.

The proof of the uniqueness part is the same as in Lemma 5.6.

Let us show existence of a solution to (5.28). Let $(v^{n}_{0})_{n\in\mathbb{Z}_{+}}$ be a sequence of $\mathcal{C}^{2}(\mathbb{T}^{d})$ functions converging in $L_{2}$ to $v_{0}\in L_{2}(\mathbb{T}^{d})$. Let $v^{n}$ be the analytically strong solution of (5.28) with the initial condition $v^{n}_{0}$ (it exists thanks to Lemma 5.5). By Lemma 5.7(i,ii) there exists $v\in\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$ such that

[TABLE]

We claim now that $v$ is an analytically generalized strong solution to (5.28) with the initial condition $v_{0}$.

Indeed, we have by construction $v(0)=v_{0}$ and $v\in\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$ as a limit of the functions from the space $\mathcal{C}([0,T],L_{2}(\mathbb{T}^{d}))$. Fix any $\varepsilon\in(0,T]$ and let us verify that identity (4.2) holds.

By (5.38), we have $v(\varepsilon/2)\in L_{\infty}$. Therefore, by Lemma 5.8 there exists a process $\bar{v}\in\mathcal{C}((\varepsilon/2,T],\mathcal{C}^{2}(\mathbb{T}^{d}))$, which is an analytically generalized strong solution of (5.28) on interval $[\varepsilon/2,T]$ with the initial condition $v(\varepsilon/2)$. Therefore, by Lemma 5.7(iii) we have $\bar{v}(t)=v(t)$ for any $t>\varepsilon/2$. Thus, identity (4.2) holds and for any $t\geqslant\varepsilon$ the function $v(t)$ has a version which is in $\mathcal{C}^{2}(\mathbb{T}^{d})$. Since $\varepsilon$ was arbitrary, this implies the statement of the lemma. ∎

Now we are ready to present the proof of Theorem 4.4.

Proof of Theorem 4.4.

(i). The proof of the uniqueness part is the same as in Lemma 5.5. To show existence of a strong solution to (1.1) we fix an initial condition $u_{0}\in\mathcal{C}^{2}(\mathbb{T}^{d})$, $T>0$ and first show existence on the time interval $[0,T]$. Let $v$ be an analytically strong solution to (5.28) with the initial condition $v(0)=u_{0}$. Recall the definition of $\Omega^{\prime}$ from Lemma 5.4 and put now

[TABLE]

and $u(t,\xi,\omega)=0$ for $\omega\in\Omega\setminus\Omega^{\prime}$. It follows from Lemmas 5.4 and 5.5 that the function $u$ is an analytically strong solution to (1.1) with the initial condition $u_{0}$. To show the adaptiveness of $u$, we introduce a function $u^{n}$, $n\in\mathbb{Z}_{+}$, which is an analytically and probabilistically strong solution on $[0,T]$ to the following equation

[TABLE]

where for $x\in\mathbb{R}$, $\xi\in\mathbb{T}^{d}$ we put $f^{n}(x,\xi):=f((x\wedge n)\vee(-n),\xi)$. Since $f^{n}$ is uniformly bounded, it follows from [22, Chapter II, Theorem 2.1] that $u^{n}$ is well–defined; thus, identity (5.39) holds on some set $\Omega_{n}$ of full measure. On the other hand, by uniqueness, we have

[TABLE]

and $u^{n}=u$ on $A_{n}$. Since for each $\omega\in\Omega^{\prime}$ the function $u$ is bounded, we see that $\mathsf{P}(A_{n})\to 1$ as $n\to\infty$ and $A_{n}\subset A_{n+1}$. This implies that $u^{n}$ converges to $u$ a.s. as $n\to\infty$. Since each $u^{n}$ is $(\mathcal{F}_{t})$–adapted, their limit $u$ is also $(\mathcal{F}_{t})$–adapted. Therefore $u$ is a probabilistically strong solution of (1.1) on $[0,T]$. Since $T$ was arbitrary, this and uniqueness imply that $u$ is an analytically and probabilistically strong solution of (1.1) on $[0,\infty)$. The continuity of $u$ follows from continuity of $v$ and $w$ (Lemmas 5.4 and 5.5).

(ii). The proof of the uniqueness part is the same as in Lemma 5.6. To show existence fix $T>0$ and put again now

[TABLE]

and $u(t,\xi,\omega)=0$ for $\omega\in\Omega\setminus\Omega^{\prime}$. Here $v$ is an analytically generalized strong solution to (5.28) with the initial condition $v(0)=u_{0}$. It follows from Lemmas 5.4 and 5.6 that the function $u$ is an analytically generalized strong solution to (1.1) with the initial condition $u_{0}$. To show the adaptiveness we consider $(u^{n}_{0})_{n\in\mathbb{Z}_{+}}$, a sequence of $\mathcal{C}^{2}(\mathbb{T}^{d})$ functions converging in $L_{2}$ to $u_{0}\in L_{2}(\mathbb{T}^{d})$. Let $u^{n}$ be the analytically and probabilistically strong solution of (1.1) with the initial condition $u^{n}_{0}$ (it exists thanks to part (i) of the theorem). By Lemma 5.7(iii), we have that for each $t>0$ the function $u^{n}(t)$ converges to $u(t)$ in $L_{2}(\mathbb{T}^{d})$ $\mathsf{P}$-a.s. Since $u^{n}(t)$ is $\mathcal{F}_{t}$–adapted, we see that $u(t)$ is also $(\mathcal{F}_{t})$–adapted. Therefore $u$ is a probabilistically strong solution of (1.1) on $[0,T]$. The proof ends in the same way as in the part (i) of the theorem.

(iii). We begin with the following two observations. First, we note that part (ii) of the theorem implies that the generalized strong solution of equation (1.1) is unique.

Second, for $x\in L_{2}(\mathbb{T}^{d})$ let us denote by $u^{x}$ the generalized strong solution of (1.1) with the initial condition $x$. Then, by the chain rule and Gronwall’s inequality we get that for each $t>0$ there exists $C(t)>0$ such that for any $x,y\in L_{2}(\mathbb{T}^{d})$ we have

[TABLE]

where we also took into account assumption (4.4).

Now the Markov property of $u$ follows from these two properties by repeating literally the argument from [3, Proposition 4.1].

(iv). Take the set $\Omega^{\prime}$ defined in Lemma 5.4. Fix $t>0$, $\omega\in\Omega^{\prime}$, $x,y\in L_{2}(\mathbb{T}^{d})$ such that $x\preceq_{pos}y$. Let $u^{x}$, $u^{y}$ be the generalized strong solutions of (1.1) with the initial conditions $x$ and $y$ respectively. Let us show now that $u^{x}(t,\omega)\preceq_{pos}u^{y}(t,\omega)$ a.s.

We begin with the case when $x,y\in\mathcal{C}^{2}(\mathbb{T}^{d})$. Let $w$ be the stochastic convolution (recall its definition in Lemma 5.4). Let $v^{x}$, $v^{y}$ be analytically strong solutions of (5.28) with the initial conditions $x$ and $y$ respectively. Then by Lemma 5.5 $v^{x}$ and $v^{y}$ are continuous on $[0,t]\times\mathbb{T}^{d}$. Therefore, by comparison principle [36, Theorem 10.1] we have $v^{x}(t,\omega)\preceq_{pos}v^{y}(t,\omega)$ and hence $u^{x}(t,\omega)=v^{x}(t,\omega)+w(t,\omega)\preceq_{pos}v^{y}(t,\omega)+w(t,\omega)=u^{y}(t,\omega)$.

In the general case, we consider $(x^{n})_{n\in\mathbb{Z}_{+}}$ and $(y^{n})_{n\in\mathbb{Z}_{+}}$, two sequences of $\mathcal{C}^{2}(\mathbb{T}^{d})$ functions converging in $L_{2}$ to $x$ and $y$, respectively. By the above, for each $n\in\mathbb{Z}_{+}$ we have

[TABLE]

On the other hand, by part (ii) of the theorem and Lemma 5.7 we have $u^{x^{n}}(t,\omega)\to u^{x}(t,\omega)$, $u^{y^{n}}(t,\omega)\to u^{y}(t,\omega)$ in $L_{2}(\mathbb{T}^{d})$ as $n\to\infty$. This together with (5.40) yields $u^{x}(t)\preceq_{pos}u^{y}(t)$.

(v). Follows immediately from part (iv) of the theorem and the Strassen’s theorem.

(vi) Fix $x\in L_{2}(\mathbb{T}^{d})$. Then $u^{x}$ is an analytically generalized strong solution of (1.1). Therefore by Ito’s lemma, taking into account (4.3), we deduce for any $0<s\leqslant t$

[TABLE]

where the constants $K_{1}$, $K_{2}$ are defined in (4.3) and $M$ is a continuous local martingale with $M(0)=0$ and

[TABLE]

Using the fact that $u^{x}$ is continuous in $L_{2}(\mathbb{T}^{d})$, we can pass to the limit in (5.2) as $s\to 0$ to deduce that (5.2) is also valid for $s=0$.

Therefore all the conditions of Lemma 5.3 are satisfied for the process $X(t):=\|u^{x}(t)\|_{L_{2}(\mathbb{T}^{d})}^{2}$. Thus, $\hskip 0.58122pt\mathsf{E}\hskip 0.38751pt\|u^{x}(t)\|_{L_{2}(\mathbb{T}^{d})}^{2}<\infty$ for any $t>0$, inequality (4.6) follows from (5.21) and inequality (4.7) follows from (5.22). ∎

Before we formulate and prove the final lemma we recall that under Assumption A the solution $u$ with initial condition $u_{0}\in L_{2}(\mathbb{T}^{d})$ satisfies the following mild form of the SPDE (1.1) [7, Theorem 5.4]

[TABLE]

for each $t>0$, $\varepsilon\in(0,t]$ and $\xi\in\mathbb{T}^{d}$, where $p$ denotes the fundamental solution of the heat equation on $\mathbb{T}^{d}$. Letting $\varepsilon\downarrow 0$, we see that almost surely

[TABLE]

where we used the Cauchy-Schwarz inequality to obtain convergence of the first term.

For an arbitrary element $x\in L_{2}(\mathbb{T}^{d})$ let $S_{x}$ and $L_{x}$ be the sets of elements from $L_{2}(\mathbb{T}^{d})$ that are smaller (respectively larger) than $x$, that is

[TABLE]

It is easy to see that the sets $S_{x}$ and $L_{x}$ are closed. For $a\in\mathbb{R}$ put

[TABLE]

Lemma 5.9**.**

Suppose that Assumptions A and B hold. Let $(P_{t})$ be the semigroup associated with equation (1.1). Then for any $M>0$ there exists $\varepsilon>0$ such that

[TABLE]

Proof.

By symmetry it suffices to prove the first claim. We show it in two steps.

Step 1. At this step we will use only Assumption A. We claim that for every $M>0$ and $T\in(0,1]$ there exist $\Gamma\in\mathbb{R}$ and $t_{0}\in(0,T]$ such that

[TABLE]

First we note, that it is sufficient to prove the claim only for large enough $M$. Observe now that by (4.3) there exists $\gamma\in\mathbb{R}$ such that $f(z,\xi)\leqslant 0$ for all $z\geqslant\gamma$ and all $\xi\in\mathbb{T}^{d}$. Fix arbitrary $\beta>\gamma$, and assume that $M$ is large enough so that $\|\check{\beta}\|_{L_{2}}\leqslant M$. By order preservation (Theorem 4.4 (iv)) it suffices to show that for every $M>0$, $T\in(0,1]$ there exist $\Gamma\in\mathbb{R}$ and $t_{0}\in(0,T]$ such that

[TABLE]

By continuity of solutions (Theorem 4.4) there exists $t_{0}\in(0,T]$ such that

[TABLE]

For $x\in L_{2}(\mathbb{T}^{d})$ we define

[TABLE]

Fix now arbitrary $x\in L_{2}(\mathbb{T}^{d})$ satisfying $\check{\beta}\preceq_{pos}x$ and $\|x\|_{L_{2}}\leqslant M$. Then, using order preservation again, we see that

[TABLE]

Choose $\tilde{\Gamma}$ so large that

[TABLE]

Using (5.42), we see that on the set $\tilde{\Omega}\cap\Omega(x)$ (which has measure at least $\kappa/2$) the solution $u^{x}$ satisfies for all $\xi\in\mathbb{T}^{d}$

[TABLE]

where we used the Cauchy–Schwarz inequality in the first inequality and the fact that $f(z,\xi)\leqslant 0$ for $z\geqslant\gamma$ and $\xi\in\mathbb{T}^{d}$. The proof of the claim of Step 1 is complete.

Step 2. This is the only part of the proof of Theorem 4.5, where Assumption B is used. We claim that for every $\Gamma\in\mathbb{R}$ and $\tau>0$ we have

[TABLE]

To show this claim, first of all, we assume, without loss of generality, that $\Gamma>0$. It suffices to show that for some $\tau_{0}=\tau_{0}(\Gamma)>0$ the statement holds for every $\tau\in(0,\tau_{0}]$. Fix numbers $\lambda_{k}$, $k=1,...,m$ as in Assumption B with $\varepsilon=1$ and define $\Lambda(\xi):=\sum_{k=1}^{m}\lambda_{k}\sigma_{k}(\xi)\geqslant 1$ and $\hat{\Lambda}:=\max_{\xi\in\mathbb{T}^{d}}\Lambda(\xi)$. Choose $Q>\hat{\Lambda}(\Gamma+2)+2$ and fix $\tau_{0}\in(0,1]$ such that $\tau_{0}|f(z,\xi)|\leqslant 1$ for all $\xi\in\mathbb{T}^{d}$ and $z\in[\Gamma-Q,\Gamma+2]$ (such $\tau_{0}$ exists since $f$ is continuous). Fix $\tau\in(0,\tau_{0}]$ and $\Theta:=\frac{\Gamma+2}{\tau}$. Then, for

[TABLE]

and $t>0$, the solution $u^{\check{\Gamma}}$ satisfies

[TABLE]

Note that for each $t\geqslant 0$,

[TABLE]

Let

[TABLE]

Then $\mathsf{P}(\hat{\Omega})>0$. Let \psi:=\inf\{s\geqslant 0:\,u^{\check{\Gamma}}(s)\notin[{\mathchoice{{\ooalign{\hbox{\raise 6.8648pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.8648pt\hbox{\displaystyle\widehat{\vrule width=0.0pt,height=6.14998pt\vrule height=0.0pt,width=19.74pt}}}}}\cr\hbox{\displaystyle\Gamma-Q}}}}{{\ooalign{\hbox{\raise 6.8648pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.8648pt\hbox{\textstyle\widehat{\vrule width=0.0pt,height=6.14998pt\vrule height=0.0pt,width=19.74pt}}}}}\cr\hbox{\textstyle\Gamma-Q}}}}{{\ooalign{\hbox{\raise 6.24982pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.24982pt\hbox{\scriptstyle\widehat{\vrule width=0.0pt,height=4.305pt\vrule height=0.0pt,width=13.8179pt}}}}}\cr\hbox{\scriptstyle\Gamma-Q}}}}{{\ooalign{\hbox{\raise 5.83981pt\hbox{\scalebox{1.0}[-1.0]{\lower 5.83981pt\hbox{\scriptscriptstyle\widehat{\vrule width=0.0pt,height=3.07498pt\vrule height=0.0pt,width=9.87pt}}}}}\cr\hbox{\scriptscriptstyle\Gamma-Q}}}}},{\mathchoice{{\ooalign{\hbox{\raise 6.8648pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.8648pt\hbox{\displaystyle\widehat{\vrule width=0.0pt,height=6.14998pt\vrule height=0.0pt,width=21.12505pt}}}}}\cr\hbox{\displaystyle\Gamma+2}}}}{{\ooalign{\hbox{\raise 6.8648pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.8648pt\hbox{\textstyle\widehat{\vrule width=0.0pt,height=6.14998pt\vrule height=0.0pt,width=21.12505pt}}}}}\cr\hbox{\textstyle\Gamma+2}}}}{{\ooalign{\hbox{\raise 6.24982pt\hbox{\scalebox{1.0}[-1.0]{\lower 6.24982pt\hbox{\scriptstyle\widehat{\vrule width=0.0pt,height=4.305pt\vrule height=0.0pt,width=14.78745pt}}}}}\cr\hbox{\scriptstyle\Gamma+2}}}}{{\ooalign{\hbox{\raise 5.83981pt\hbox{\scalebox{1.0}[-1.0]{\lower 5.83981pt\hbox{\scriptscriptstyle\widehat{\vrule width=0.0pt,height=3.07498pt\vrule height=0.0pt,width=10.5625pt}}}}}\cr\hbox{\scriptscriptstyle\Gamma+2}}}}}]\}. If $\psi<\tau$ and $u(\psi,\xi)=\Gamma+2$ for some $\xi\in\mathbb{T}^{d}$, then, on the set $\hat{\Omega}\cap\{\psi<\tau\}$ we have

[TABLE]

so $\psi=0$ which is impossible (since $u$ is continuous). If, on the other hand, $\psi<\tau$ and $u(\psi,\xi)=\Gamma-Q$ for some $\xi\in\mathbb{T}^{d}$, then, on the set $\widetilde{\Omega}\cap\{\psi<\tau\}$, we have

[TABLE]

which is also impossible by the definition of $Q$ and $\Theta$.

Thus, $\psi\geqslant\tau$ almost surely on $\hat{\Omega}$. Therefore, on $\hat{\Omega}$, we have

[TABLE]

for every $\xi\in\mathbb{T}^{d}$, so the proof of Step 2 is complete.

Using order preservation once more, we see that Step 1 and Step 2 imply the assertion in the lemma. ∎

Finally, before we present the proof of Theorem 4.8, we need to introduce some additional notation. If for some $a\in\mathbb{R}$, $x\in L_{2}(\mathbb{T}^{d})$ we have $x(\xi)=a$ for all $\xi\in\mathbb{T}^{d}$, then we will write

[TABLE]

Introduce also the set of all constants in the space $L_{2}(\mathbb{T}^{d})$, that is,

[TABLE]

Let $(P_{t})_{t\geqslant 0}$ be the transition function associated with the solution to (1.1).

Lemma 5.10**.**

Consider equation (1.1) with $m=1$ and $\sigma_{1}\equiv 1$. Suppose that the drift $f$ does not depend on space. This equation has the following properties:

1. (i)

if $u(0)\in A_{c}$, then $u(t)\in A_{c}$ a.s. for any $t\geqslant 0$; 2. (ii)

if $u(0)\notin A_{c}$, then $u(t)\notin A_{c}$ a.s. for any $t\geqslant 0$; 3. (iii)

this equation has a unique invariant measure $\pi$ and $\pi(A_{c})=1$.

Proof.

(i) Suppose that $u(0)\in A_{c}$. Let $a=\Pi[u(0)]$. Consider a stochastic differential equation

[TABLE]

Since the function $f$ is locally Lipschitz and satisfies the one-sided growth condition, equation (5.44) has a unique strong solution [34, Proposition 2.1(b)]. It is immediate to see that the function $u(t,\xi):=X^{a}(t)$ is an analytically and probabilistically strong solution of (1.1) with the initial condition $u(0)$. By Theorem 4.4 this is the unique solution to (1.1). Thus $u(t)\in A_{c}$ for any $t\geqslant 0$.

(ii) Fix the initial condition $x\notin A_{c}$. First, consider the case $x\in\mathcal{C}^{2}(\mathbb{T}^{d})$. Let $a_{1},a_{2}$ be real numbers such that $a_{1}\leqslant x(\xi)\leqslant a_{2}$ for all $\xi\in\mathbb{T}^{d}$. Let $\Omega^{\prime\prime}\subset\Omega$ be the set of full measure where the trajectories of the Brownian motion $W$ are continuous and identity (4.2) holds for $\varepsilon=0$ for the initial conditions $\Pi^{-1}a_{1}$, $\Pi^{-1}a_{2}$, and $x$. Set $\Omega_{0}:=\Omega^{\prime\prime}\cap\Omega^{\prime}$, where $\Omega^{\prime}$ is defined in Theorem 4.4(iv). It follows that $\mathsf{P}(\Omega_{0})=1$.

Assume now the contrary. That is suppose that for some $T>0$, $\omega\in\Omega_{0}$ we have $u^{x}(0)=x$ and $u^{x}(T,\omega)\in A_{c}$. Denote $b:=\Pi[u^{x}(T,\omega)]$. By Theorem 4.4(iv) and part(i) of the theorem, we have for all $\xi\in\mathbb{T}^{d}$

[TABLE]

It follows immediately from the Gronwall lemma, the fact that $f$ is locally Lipschitz, and the comparison principle that solutions to (5.44) are continuous with respect to the initial condition. Therefore there exists some initial condition $a\in(a_{1},a_{2})$ such that $X^{a}(0,\omega)=a$ and $X^{a}(T,\omega)=b$. Denote now $v(t,\xi):=u(T-t,\xi,\omega)-X^{a}(T-t,\omega)$, $t\in[0,T]$, $\xi\in\mathbb{T}^{d}$. Then we have $v(0,\xi)=0$ for any $\xi\in\mathbb{T}^{d}$. On the other hand, for any $t\in[0,T]$, $\xi\in\mathbb{T}^{d}$ we have

[TABLE]

where we also used that the function $f$ is locally Lipschitz and the processes $u(\omega)$ and $X^{a}(\omega)$ are bounded. However, inequality (5.2) together with the fact that $v(0)=0$ implies that $v(t)=0$ for all $t\in[0,T]$ [37, Theorem 3.0.4]. This yields that

[TABLE]

However it was assumed that $u(0)\notin A_{c}$. This contradiction finishes the proof for the case $x\in\mathcal{C}^{2}(\mathbb{T}^{d})$.

In the general case, $x\in L_{2}(\mathbb{T}^{d})$, we will use the Markov property of the solutions to (1.1), that is, Theorem 4.4(iii). We have for any $0<s<t$

[TABLE]

where the first identity is the Markov property, the second identity follows from the fact that $u^{x}(s)\in\mathcal{C}^{2}(\mathbb{T}^{d})$ a.s. by Lemma 5.8, and the third identity follows from the fact that $P_{t-s}(y,A_{c})=0$ if $y\in\mathcal{C}^{2}(\mathbb{T}^{d})\setminus A_{c}$ established above. Since $s$ was arbitrary, we can pass to the limit as $s\to 0$. Using the fact that the process $u^{x}(\omega)$ is by definition continuous in $L_{2}(\mathbb{T}^{d})$ and $x\notin A_{c}$, we get

[TABLE]

which completes the proof.

(iii). It follows from Theorem 4.5 that SPDE (1.1) has a unique invariant measure $\pi$ and that for any $x\in L_{2}(\mathbb{T}^{d})$ one has $P_{t}(x,\cdot)\to\pi$ weakly as $t\to\infty$. Fix $x_{0}\in A_{c}$. Then, by part (i) of the proof and the Portmanteau theorem (see, e.g., [35, Theorem III.1.1(II)]), one has

[TABLE]

where we also used the fact that the set $A_{c}$ is closed in $L_{2}(\mathbb{T}^{d})$. This proves the statement of the lemma. ∎

Proof of Theorem 4.8.

Take $m=1$ and $\sigma_{1}\equiv 1$. Fix any $x\in L_{2}(\mathbb{T}^{d})\setminus A_{c}$. Then for any $t>0$ by Lemma 5.10(ii) one has $P_{t}(x,A_{c})=0$. On the other hand, by Lemma 5.10(iii), $\pi(A_{c})=1$. This implies the statement of the theorem. ∎