Lower Estimates of Transition Densities
Lower Estimates of Transition Densities
been studied for a long time and the ergodic theory of finite-dimensional
diffusion processes is relatively well developed, see, for example, a classical
monograph [21]. In this paper we study the ergodic properties of a class of
ordinary and partial stochastic differential equations that includes stochastic
reaction–diffusion equations in bounded domains. First results on the exis-
tence and uniqueness of invariant measures for stochastic reaction-diffusion
equations were obtained in [12, 26, 42], see also [5], the monographs [6, 10]
and references therein. The rate of convergence to the invariant measure in
infinite dimensions became a subject of interest much later and still is not
well understood. Jacquot and Royer [24] proved exponential ergodicity for
a semilinear parabolic equation with bounded nonlinear drift, Shardlow [39]
applied the theory of Meyn and Tweedie to obtain V -uniform ergodicity for
some semilinear equations in Hilbert spaces. Hairer in [14] proved, under dif-
ferent sets of conditions, uniform exponential ergodicity for equations with
drifts growing faster than linearly. Exponential convergence to equilibrium in
a norm intermediate between the total variation metric and the Wasserstein
metric has been obtained in [31] for the stochastic Navier–Stokes equation.
A closely related problem of asymptotic behavior of the Markov semigroup
Pt φ(x) = Eφ(Xtx ) attracted much attention due to its importance in Mathe-
matical Physics. In particular, exponential convergence of the semigroup in
the spaces Lp (E, µ∗ ), p ∈ [1, ∞), and related questions of the existence of
the spectral gap and logarithmic Sobolev inequality have been studied by
numerous authors, see [1, 2, 3, 7, 8, 11, 18, 22, 43].
The aim of the present paper is to prove V -uniform (exponential) ergodic-
ity with V (x) = |x|E + 1 and, if the drift grows faster than linearly, uniform
exponential ergodicity, for equation (1.1). Our method allows us to find ex-
act bounds on convergence (i.e., to give explicit estimates for the rate of
exponential convergence in the total variation norm or V -variation norm).
In this respect, our results seem to be new even for finite-dimensional SDE’s
(which is also due to our method to estimate the transition density that, to
the best of our knowledge, has not been used in finite dimensions so far). If
the Markov semigroup (Pt ) is symmetric, we obtain explicit lower estimates
for the spectral gap in L2 (E, µ∗ ). Stronger results are obtained in the case
of a drift growing faster than linearly: for a symmetric Markov semigroup,
we show uniform estimates on the spectral gap in the spaces Lp (E, µ∗ ) for
all p ∈ [1, ∞) and in the nonsymmetric case, our estimates remain valid for
p > 1, in particular, in L2 (E, µ∗ ).
Unlike in the aforementioned papers, in the present paper a lower bound
measure and a suitable small set for a skeleton process are found explicitly
in terms of the lower estimates of transition densities and the constants in
an ultimate boundedness condition (or, in particular, a suitable Lyapunov
function). This enables us to apply earlier results on computable bounds for
ERGODICITY FOR STOCHASTIC PDE’S 3
Markov chains, which are expressed in terms of lower bound measures, corre-
sponding small sets and constants from the Lyapunov–Foster geometric drift
condition [33]. The bounds turn out to be uniform with respect to a large
family of drift coefficients, which is important for proving continuous depen-
dence of invariant measures on parameter (cf. Section 8). We also believe
that this uniformity is an important tool for studying the Hamilton–Jacobi–
Bellman equation for the ergodic control problem. On the other hand, the
method employed here has its limitations. Our method strongly relies on
the Girsanov theorem and therefore, we need √ an assumption that F maps
the whole state space into the range of Q. Therefore, any extension to
other types of equations (like stochastic Burgers or Navier–Stokes equations)
would be difficult. Note, however, that in the two recent authors’ papers
[17] and [16] V -uniform ergodicity and spectral gap type results have been
proved for stochastic Burgers, 2D Navier–Stokes and more general reaction–
diffusion equations. Nonetheless, in these papers a different method is used
that allows us neither to find explicit bounds on the convergence constants
nor to show the uniformity of convergence with respect to coefficients.
An important tool for our proofs is a formula for the transition densities
that is derived in this paper. We use this formula to establish suitable lower
estimates on the densities which we believe are of independent interest. They
are obtained by means of the Girsanov theorem and the so called Ornstein–
Uhlenbeck bridge (or pinned Ornstein–Uhlenbeck process). Let us explain
the main idea of this approach.
Let (Ztx ) be an Ornstein–Uhlenbeck process on a separable Hilbert space
H. By this, we mean that (Ztx ) is a solution to a linear stochastic evolution
equation
p
dZtx = AZtx dt + Q dWt ,
(1.2)
Z0x = x ∈ H.
The Ornstein–Uhlenbeck bridge (Zbtx,y ) associated to the Ornstein–Uhlenbeck
process (Ztx ) is informally defined by the formula
P(Ztx ∈ B|Z1x = y) = P(Zbtx,y ∈ B), t < 1,
where x, y ∈ H and B ⊂ H is a Borel set. The importance of various types
of bridge processes for the study of transition densities of finite dimensional
diffusions is well recognised, see, for example, [23]. In infinite-dimensional
framework this concept was developed in [41] in order to study the regularity
of transition semigroups of diffusions on Hilbert spaces. In [28] and [29] an
Ornstein–Uhlenbeck bridge is introduced in order to obtain lower estimates
on the transition kernel of some semilinear stochastic evolution equations.
The basic idea is as follows. Using the equivalence of measures corresponding
4 B. GOLDYS AND B. MASLOWSKI
to Xtx and Ztx and the Girsanov formula, we can write the transition density
of the process Xtx in the form
d(T, x, y) = E(Φ(Z·x )|ZTx = y),
where Φ is a measurable functional defined on trajectories of the Ornstein–
Uhlenbeck process. This form of the density is not suitable for the uniform
estimates that are needed. Therefore, the conditional expectation is trans-
formed into a usual expectation with respect to the measure of the OU
bridge (Zbtx,y ) considered for t ∈ [0, T ]:
d(T, x, y) = EΦ(Zb·x,y ),
which enables us to find the uniform lower estimates. Let us note here a
technical difficulty caused by the fact that we can define the OU bridge for
y in a certain Borel subspace of measure one only, but this turns out to be
sufficient for our needs.
Precise formulations and hypothesis are given in the following Section 2.
In Sections 3 and 4 the properties of the OU bridge, which are needed in
the sequel, are established (some auxiliary results are deferred to the Ap-
pendix). The formula for transition densities is found and the lower estimates
are given in Section 5. These results are applied in Section 6 to establish
our main results, V -uniform ergodicity and uniform exponential ergodicity,
respectively, and find the computable bounds on respective constants. In
Section 7 the corollaries on the Lp (E, µ∗ ) exponential convergence and the
spectral gap are stated. Section 8 is devoted to some extensions and ap-
plications (continuous dependence of invariant measures on a parameter).
Examples (finite-dimensional nonlinear stochastic oscillator and stochastic
parabolic equations) are presented in Section 9.
Hypothesis 2.3.
Z 1
−1/2
(2.7) kQt St Q1/2 kHS dt < ∞,
0
Hypotheses 2.1–2.4 are standing assumptions of the paper and the results
will be enunciated without recalling them again. Obviously, the hypotheses
are used selectively (e.g., Hypothesis 2.4 is not needed for results on the
Ornstein–Uhlenbeck bridge).
Denote by B, P and bB, the Borel σ-algebra of E, the space of probability
measures on E and the space of bounded Borel functions on E, respectively.
Furthermore,
Pt ϕ(x) := Ex ϕ(Xt ), φ ∈ bB, x ∈ E, t ≥ 0,
and
P (t, x, Γ) := Pt 1Γ (x), x ∈ E, Γ ∈ B, t ≥ 0.
Let (Pt∗ ) denote the adjoint Markov semigroup, that is,
Z
(2.10) Pt∗ ν(Γ) := P (t, x, Γ)ν(dx), t ≥ 0, ν ∈ P, Γ ∈ B.
Γ
An invariant measure µ∗
∈ P is defined as a stationary point of the semi-
group (Pt∗ ), that is, Pt∗ µ∗ = µ∗ for each t ≥ 0. Obviously, Pt∗ ν is interpreted
as the probability distribution of Xt if X0 has the initial distribution is ν.
ERGODICITY FOR STOCHASTIC PDE’S 7
−1/2 1/2
Lemma 3.1. The operator Vt = Q1 S1−t Qt is bounded on H and
(3.2) kVt k < 1, t ∈ (0, 1].
Moreover,
(3.3) lim Vt∗ x = lim Vt x = x, x ∈ H.
t→1 t→1
Proof. Estimate (3.2) was proved in [35]. It follows from (3.2) and a
simple identity
∗
Q1 = Q1−t + S1−t Qt S1−t ,
that
1/2 1/2
(3.4) Q1−t = Q1 (I − Vt Vt∗ )Q1 .
To prove (3.3), we will show first that
(3.5) lim hVt x, yi = hx, yi, x, y ∈ H.
t→0
8 B. GOLDYS AND B. MASLOWSKI
−1/2
Indeed, for y ∈ im(Q1 ), we have
1/2 −1/2
lim hVt x, yi = limhS1−t Qt x, Q1 yi = hx, yi.
t→1 t→1
−1/2
For arbitrary y ∈ H, we may find a sequence (yn ) ⊂ im(Q1 ), such that
yn → y in H and then (3.2) yields
hVt x, y − yn i → 0,
uniformly in t ≤ 1, and (3.5) follows. Next, (3.4) yields
1/2 1/2
hQ1−t x, xi = h(I − Vt Vt∗ )Q1 x, Q1 xi, x ∈ H.
1/2
It follows that, for each y ∈ im(Q1 ), we have
lim (|y|2 − |Vt∗ y|2 ) = 0,
t→1
Now, invoking (3.5), we obtain the first part of (3.3). It is enough to prove
the second part of (3.3) for x such that |x| = 1. In this case (3.5) implies
hVt x, xi → 1, and thereby, invoking (3.2),
1 = lim inf hVt x, xi ≤ lim inf |Vt x| ≤ lim sup |Vt x| ≤ 1.
t→1 t→1 t→1
1/2 −1/2
Clearly, Vt∗ = Qt S1−t
∗ Q
1 and the operator
1/2
(3.6) Kt := Qt Vt∗
is of Hilbert–Schmidt type on H. Then the operator
H ∋ x → Kx(t) := Kt x ∈ L2 (0, 1; H)
is bounded.
Let µ denote the probability distribution of the process {Zt , t ∈ [0, 1]}
concentrated on L2 (0, 1; H) and let L : L2 (0, 1; H) → C(0, 1; H) be defined
by the formula
Z t
Lu(t) = St−s Q1/2 u(s) ds.
0
The space im(L) endowed with the norm
kφk = inf{|u| : u ∈ L2 (0, 1; H), Lu = φ}
ERGODICITY FOR STOCHASTIC PDE’S 9
Lemma 3.2. (a) For every t ∈ [0, 1), the operator Bt with the domain
1/2
dom(B) = Q1 (H) extends to a Hilbert–Schmidt operator Bt : H → H and
Z 1
(3.7) kBt kHS dt < ∞.
0
1/2
(b) The operator B with the domain dom(B) = Q1 (H) extends to a
bounded operator B : H → L2 (0, 1; H) and
(3.8) |Bx|L2 (0,1;H) = |x|H , x ∈ H.
(c) We have K = LB and the operator K : H → C(0, 1; E) is γ-radonifying.
−1/2 1/2
Proof. (a) Note first that kQ1 Q1−t k ≤ 1 and by (2.4), the operator
−1/2
Q1−t S1−t is bounded. Therefore, the operator
−1/2 1/2 −1/2 −1/2
(Q1 Q1−t Q1−t S1−t Q1/2 )∗ = Q1/2 S1−t
∗ Q
1
is bounded. Moreover, taking (2.7) and (2.6) into account, we obtain, for a
certain C > 0,
Z 1 Z 1
−1/2
kBt∗ kHS dt ≤ C kQt St Q1/2 kHS dt < ∞
0 0
and (3.7) follows.
For any h ∈ H, we have
Z 1
1/2
|Q1 h|2 = |Q1/2 S1−t
∗
h|2 dt,
0
−1/2
and therefore, for h = Q1 x, we obtain
Z 1
(3.9) |x|2 = |Bx(t)|2 dt.
0
10 B. GOLDYS AND B. MASLOWSKI
1/2
Using the density of Q1 (H) in H, we can extend (3.9) to the whole of H
and (3.8) follows.
1/2
(c) For x ∈ Q1 (H), we have
∗ −1/2
Kt x = Qt S1−t Q1 x
Z t
∗ ∗ −1/2
= St−s QSt−s S1−t Q1 x ds
0
(3.10) Z t
∗ −1/2
= St−s QS1−s Q1 x ds
0
Z t
= St−s Q1/2 Bx(s) ds = L(Bx)(t),
0
for all t ∈ [0, 1]. By (b), this identity extends to all x ∈ H and we find
that K = LB on H. By Hypothesis 2.2, we have µ(C(0, 1; E)) = 1, hence,
L : L2 (0, 1; H) → C(0, 1; E) is γ-radonifying and therefore, K = LB : H →
C(0, 1; E) is γ-radonifying as well.
and
−1/2 c
(3.12) kQt St k ≤ .
tβ
Then Hypothesis 2.3 is satisfied.
Remark 3.4. Conditions (3.11) and (3.12) are well known and often
used in the theory of SPDE’s. Condition (3.11) is a standard assumption
that implies the existence of an H-continuous version of the OU process
(Ztx ), while (3.12) is closely related to the existence and integrability of the
gradient of the OU transition semigroup (cf. [9] for details). Hypothesis 2.3
will be checked in more specific cases in Section 9.
(4.2) b t ) = 0.
lim tr(Q
t→0
To prove that
(4.3) b t ) = 0,
lim tr(Q
t→1
We have also
∞
X 1/2
tr((I − Vt∗ Vt )Q1 ) = tr(Q1 ) − |Vt Q1 ek |2 ,
k=1
and (4.3) follows from (3.3), (3.2) and the dominated convergence. Since the
process (Zbt ) has E-continuous version by Lemma 3.2(b), it follows that
Proof. It was already shown in Proposition 4.1 that the process (Zbt )
has trajectories in C(0, 1; E) and we have
(4.8) k1 = sup |Zbt |E < ∞, P-a.s.
t≤1
Proof. By Hypothesis 2.2 and Theorem 4.3, the processes (Ztx ) and
−1/2
(Zbtx,y ) are concentrated on C(0, 1; E) and (Kt Q1 (S1 x − y)) ∈ E. More-
−1/2
over, the processes (Zbtx,y ) and (Kt Q1 Z1x ) are independent by Proposi-
tion 4.1. Therefore, using well-known properties of conditional expectations,
we obtain
−1/2
E(Φ(Z x )|Z1x = y) = E(Φ(Zb x,y + Kt Q1 (Z1x − y))|Z1x = y)
= EΦ(Zb x,y ), µ1 -a.e.
for any x ∈ E.
Let
Z 1
Yu = S1−s Q1/2 dWs , u ≤ 1,
u
and
−1/2
Hu = Q1−u S1−u Q1/2 , u < 1.
Proof. We have
Z t
−1/2 1/2 −1/2
Kt Q1 Z1 = St−s Q Hs∗ ds Q1 Z1
0
and
Z t
−1/2 −1/2
S1−t Kt Q1 Z1 = S1−s Q1/2 Hs∗ ds Q1 Z1
0
and thereby,
−1/2
Z1 − S1−t Kt Q1 Z1 = Q1−t Q−1
1 Z1 .
= Q1−s Q−1 b
1 Z1 − S1−s Zs .
−1/2
In what follows we will use the notation hαu , hi = hQ1−u Yu , Hu hi.
Therefore,
n
X 1
Ehαnu − αm 2
u , hi = EhYu , ek (1 − u)i2 hek (1 − u), Hu hi2
j=m+1
λ k (1 − u)
(4.16)
n
X
= hek (1 − u), Hu hi2 −→ 0,
n,m→∞
j=m+1
We have
Ehζt − ζr , hihζr , ki
Z r
−1/2
= −EhWt − Wr , hi hQ1−u Yu , Hu ki du
0
Z t
−1/2
− EhWr , ki hQ1−u Yu , Hu hi du
r
Z t Z r
−1/2 −1/2
+E hQ1−u Yu , Hu hi du hQ1−u Yu , Hv ki dv
r 0
= −I1 − I2 + I3 .
We will consider I1 first. Taking into account that the series (4.17) is mean-
square convergent, for each u ∈ (0, 1), we obtain
−1/2
EhWt − Wr , hihQ1−u Yu , Hu ki
∞
X hen (1 − u), Hu ki
= p E(hYu , en (1 − u)ihWt − Wr , hi).
n=1 λn (1 − u)
Next, for u ≤ r,
and therefore,
−1/2
EhWt − Wr , hihQ1−u Yu , Hu ki
∞
X Z
hen (1 − u), Hu ki t
= p hS1−s Q1/2 h, en (1 − u)i ds
n=1 λn (1 − u) r
Z ∞
!
t X hen (1 − u), Hu ki
= p hS1−s Q1/2 h, en (1 − u)i ds
r n=1 λn (1 − u)
Z t
−1/2
= hQ1−u S1−s Q1/2 h, Hu ki ds
r
Z t
−1/2 1/2
= hQ1−u Q1−s Hs h, Hu ki ds
r
and
Z rZ t
−1/2 1/2
(4.19) I1 = hQ1−u Q1−s Hs h, Hu ki ds du,
0 r
−1/2 1/2
where the operator Q1−u Q1−s is bounded. By similar arguments, we find
that I2 = 0 and (4.13) yields
Z rZ t
−1/2 1/2
I3 = hQ1−v Q1−u Hu h, Hv ki du dv
0 r
We have
Z t
−1/2
Ehζt , hi2 = t|h|2 − 2EhWt , hi hQ1−u Yu , Hu hi
0
Z t Z t
−1/2 −1/2
+E hQ1−u Yu , Hu hi du hQ1−v Hv Yv , ki dv
0 0
2
= t|h| − 2J1 + J3 .
Since t < 1 and the functions under the integrals are continuous, we can
change the order of integration in the first integral and obtain
Z tZ t
−1/2 1/2
J3 = hHu h, Q1−u Q1−v Hv hi dv du
0 u
Z tZ t
−1/2 1/2
+ hQ1−v Q1−u Hu h, Hv hi du dv.
0 v
Hence, J3 = 2J1 and (4.20) follows. Combining (4.18) and (4.20), we find
that, for s, t < 1,
Since
for all t ≤ 1. Therefore, (ζt ) is a cylindrical Brownian motion for t ∈ [0, 1]. It
remains to show that, for any t < 1,
Now, we have
Z t
(4.24) EhWt , hihZ1 , ki = hS1−u Q1/2 h, ki du.
0
20 B. GOLDYS AND B. MASLOWSKI
Invoking (4.10) and using the fact that Z1 and (Zbt ) are independent, we
obtain
Z t
E hαu , hihZ1 , ki du
0
Z t
−1/2
= EhQ1−u Yu , Hu hihZ1 , ki du,
0
Z t
−1/2
(4.25) EhQ1−u Q1−u Q−1
1 Z1 , Hu hihZ1 , ki du
0
Z t
1/2
= hQ1−u Hu h, ki du
0
Z t
= hS1−u Q1/2 h, ki du.
0
Comparing (4.24) and (4.25), we obtain (4.23) and the lemma follows.
Then
Z 1
(4.26) E |B1 (s)Zbs | ds < ∞,
0
Z 1
−1/2
(4.27) |B2 (s)x|2 ds = |Q1 S1 x|2 , x ∈ H.
0
−1/2 −1/2
and B3 (s)Z1 has the covariance (Q1 S1−s Q1/2 )∗ (Q1 S1−s Q1/2 ) for each
s ∈ [0, 1).
ERGODICITY FOR STOCHASTIC PDE’S 21
and taking (3.7) into account, we may apply Lemma A.3, which yields the
desired result.
where
Z 1 Z 1
(5.4) ρ(Z x ) := hG(Zsx ), dWs i − 12 |G(Zsx )|2 ds.
0 0
For n ≥ 1, set
F
(x), if |x|E ≤ n,
(5.5) Fn (x) = nx
F , if |x|E > n,
|x|E
and let Gn be defined by Fn (x) := Q1/2 Gn (x). Obviously the approximating
equations
p
dXn (t) = (AXn (t) + Fn (Xn (t))) + Q dWt ,
(5.6)
Xn (0) = x,
have uniquely defined solutions P-a.s. in C(0, 1; E) and denoting by P̃X , P̃Xn
and P̃Z x the distributions in C(0, 1; E) of X, Xn and Z x , respectively, we
have
(5.7) lim P sup |Xn (t) − X(t)|E > 0 = 0.
n→∞ t∈[0,1]
Hence,
(5.8) kP̃Xn (·) − P̃X kvar → 0, m → ∞,
We will now state one of our main results, which provides a formula for
dP ∗ ν
the density dµt t for a given time t > 0 (we may take t = 1). It follows from
the Fubini theorem that the density has the form
Z
dP1∗ ν dP (1, x, ·)
(y) = (y)ν(dx)
dµ1 E dµ1
(5.11) Z
dP (1, x, ·) dµx
= x (y) 1 (y)ν(dx), µ1 -a.e.,
E dµ1 dµ1
provided the product of densities inside the integral on the r.h.s. is (x, y)-
measurable. As mentioned in the preceding proof, the Gaussian measures
µx1 and µ1 are equivalent with the density given by the Cameron–Martin
formula
dµx
g(x, y) := 1 (y)
dµ1
(5.12)
−1/2 −1/2 1 −1/2 2
= exp hQ1 S1 x, Q1 yi − |Q1 S1 x| , x ∈ E,
2
for µ1 -almost all y ∈ E.
Proof. Since both g and h are (x, y)-measurable, taking into account
(5.11) and (5.12), we only have to prove that
dP (1, x, ·)
(5.16) (y) = h(x, y), x ∈ E, µ1 -a.e.
dµx1
Assume at first that the mapping G is bounded and let tki := ki for k ∈
N, i = 0, 1, . . . , k, that is, ∆k = {tk0 , tk1 , . . . , tkk } are equidistant divisions of
the interval [0, 1], tk0 = 0, tkk = 1, tki+1 − tki = 1/k (for brevity, the dependence
of tki on k will be suppressed in the notation). Set, for k ≥ 1,
k−1
X Z 1
(5.17) ρk (Z x ) := hG(Ztxi ), Wti+1 − Wti i − 1
2 |G(Zsx )|2 ds.
i=0 0
for µx1 -almost all y ∈ H. On the other hand, in terms of the cylindrical
Wiener process ζt defined in Lemma 4.7, we have
k−1
X
ρk (Z x ) = hG(Ztxi ), ζti+1 − ζti i
i=0
(5.21) Z Z
ti +1 1
−1/2
+ hG(Ztxi ), Hs∗ Q1−s Ys i ds − 12 |G(Zsx )|2 ds.
ti 0
E(exp ρk (Z x )|Z1x = y)
= E(exp ρk (Z x )|Z1 = y − S1 x)
k−1
X
= E exp hG(Ztxi ), ζti+1 − ζti i
i=0
Z ti+1
− hG(Ztxi ), B1 (s)Zbs − B3 (s)Z1 i ds
ti
Z ! !
1
(5.23) − 1
2 |G(Zsx )|2 ds Z1 = y − S1 x
0
k−1
X
= E exp hG(Zbtx,y
i
), ζti+1 − ζti i
i=0
Z ti+1
− hG(Zbtx,y
i
), B1 (s)Zbs − B3 (s)(y − S1 x)i ds
ti
Z !
1
− 21 |G(Zbsx,y )|2 ds
0
for x ∈ E, y ∈ M, so (5.28), (5.30) and (5.31) yield (5.16) and the proof is
completed.
for y from a set of µ1 -full measure in E. Note that in the proof of Theorem
5.2, we found a set M, µ1 (M) = 1, such that B3 (·)y ∈ L1 (0, 1 : H) and h(x, y)
−1/2
is well defined for y ∈ M1 . Similarly, S1∗ Q1 is a Hilbert–Schmidt operator
−1/2 −1/2
by (2.4), hence, S1∗ Q1 Q1 y ∈ H is well defined for y ∈ M2 , µ1 (M2 ) = 1
and the density g(x, y) is given by the formula (5.12) for y ∈ M2 . We may
take M1 = M ∩ M2 . It follows from Hypothesis 2.4(b) and (4.7) that the
stochastic integral in ρ(Zb x,y ) is a martingale and, hence, for any x ∈ E,
dP (1, x, ·)
(y)
dµ1
Z 1
1
≥ exp − E|G(Zbsx,y )|2 ds
2 0
Z 1
−E |G(Zbsx,y )|(|B1 (s)Zbsx,y | + |B2 (s)x| + |B3 (s)y|) ds
0
−1/2 −1/2 1 −1/2
− |x||S1∗ Q1 Q1 y| − |Q S1 x|2
2 1
Z 1
2
≥ exp −K 1+ E|Zb x,y |2m ds
s E
0
Z 1
−E K(1 + |Zbsx,y |m b
E )(|B1 (s)Zs | + |B2 (s)x| + |B3 (s)y|) ds
0
1 2 −1/2
− c̃|x|E |S1∗ Q−1
1 y| − c̃ kQ1 S1 k2 · |x|2E ,
2
x ∈ E,
for each q < ∞. Therefore, for each η > 0 small enough, there exist constants
c1 (η) > 0 and c3 (η) > 0 and a function Λ = Λη : M1 → R+ such that, x ∈ E
and y ∈ M1 ,
dP (1, x, ·)
(y) ≥ exp −c1 (η) − (K 2 L(2m) + η)|x|2m
E
dµ1
−1/2
− KL(m)kQ1 S1 k · |x|m+1
E
1 2 −1/2
− c3 (η)|x|m+η
E − c̃ kQ1 S1 k2 + η |x|2E − Λ(y) ,
2
and the estimate (5.32) follows.
and
Z t
(5.39) U (y) = sup ∗
St−s QS1−s Q−1
1 y ds ≤ a2 |y|Ẽ , y ∈ Ẽ,
t∈[0,1] 0 E
ERGODICITY FOR STOCHASTIC PDE’S 29
for some a1 , a2 > 0. Then for some constants b1 , b2 , b3 > 0 (dependent only
on A, Q, K and m) and p = max(2, 2m), we have
dP (1, x, ·)
(5.40) (y) ≥ b1 exp{−b2 |x|pE − b3 |y|pẼ }, x ∈ E, y ∈ Ẽ a.e.
dµ1
To see (5.40), we check that, under present conditions (5.37)–(5.39), (5.34) im-
plies, for all for x ∈ E and for µ1 -a.e. y ∈ Ẽ,
dP (1, x, ·)
(y) ≥ exp(−C(1 + |x|2m 2m
E + |y|Ẽ − |x|E + |y|Ẽ
dµ1
+ |x|m+1
E + |y|m
Ẽ
|x|E + |x|m
E |y|Ẽ
+ |x|m+η
E + |y|m+η
Ẽ
+ |x|2E + |y|2Ẽ )),
for arbitrary small η > 0 and a universal C = C(η) < ∞, and (5.40) follows.
It may be of interest to mention some particular cases when the conditions
(5.37)–(5.39) are satisfied. A trivial example is a finite-dimensional one,
H = E = Ẽ = Rd , in which case we obtain
dP (1, x, ·)
(5.41) (y) ≥ b1 exp(−b2 |x|pRd − b3 |y|pRd ), x, y ∈ Rd .
dµ1
Note that in this case the only assumptions in Theorem 5.2, Corollary 5.3
and the present remark are the well posedness and growth conditions in
Hypothesis 2.4 and the strong Feller property of the Ornstein–Uhlenbeck
process (2.7).
Suppose that A = A∗ is strictly negative and define Hλ = dom((−A)λ ), λ ≥
0, with the norm |y|λ = |(−A)λ y|, y ∈ dom((−A)λ ). Let Q = I; then A−1
must be compact and it is easy to check that im(S1 ) ⊂ im(Q1 ) = dom(A),
therefore, (5.37) holds with any Ẽ, Ẽ ֒→ H. Furthermore, we have
for y ∈ Hλ since A−1 Q−1 1 ∈ L(H), thus, (5.38) holds for Ẽ = Hλ with any
λ > 0. If, in addition, kSt kL(H,E) ≤ const · t−σ , t ∈ [0, 1], for σ > 0 such that
σ < λ, then (5.39) holds as well since
Proposition 6.1. Assume that the growth condition (2.11) holds true
and
(6.1) k(p) := sup E|Zt |pE < ∞, p > 0.
t≥0
Then
k2 k(s) + k3
(6.2) Ex |Xt |E ≤ e−k1 t |x|E + + k(1), t ≥ 0.
k1
Suppose that the following stronger version of (2.11) holds: For each x ∈
dom(Ã), there exists x∗ ∈ ∂|x|E such that, for some k1 , k2 , k3 > 0, s > 0, ε >
0, we have
(6.3) hÃx + F (x + y), x∗ iE,E ∗ ≤ −k1 |x|1+ε s
E + k2 |y|E + k3 , y ∈ E.
Then
(6.4) c,
sup sup Ex |Xt |E ≤ M
x∈E t≥1
where
1+ε 1/ε
c = k(1) + max 2(k2 k(s) + k3 ) 1
(6.5) M , +2 .
k1 k1 ε
Proof. Inequality (6.4) has been proven in [14], Proposition 2.1 (see
also a similar result in [20]). The proof of (6.2) follows the lines of similar
proofs based on Yosida approximation techniques (see, e.g., [9]) and we
sketch it only. The process Y x (t) := Xtx − Zt satisfies the equation
Z t
x
(6.6) Y (t) = St x + St−s F (Y x (s) + Zs ) ds, t ≥ 0,
0
and the sequence of approximating processes Yλ (t) is defined by
Z t
(6.7) Yλx (t) = R(λ)St x + R(λ)St−s F (Y x (s) + Zs ) ds, t ≥ 0,
0
where R(λ) := λ(λI − Ã)−1 ∈ L(E) is well defined for λ large enough. It is
well known that
dYλ
(6.8) Yλx → Y x , − ÃYλ − F (Yλ + Z) = σλx → 0, λ → ∞,
dt
in C(0, T ; E) (cf. page 201 of [9]). Since by (2.11)
d− x
(6.9) |Y (t)|E ≤ −k1 |Yλ (t)|E + k2 |Ztx |sE + k3 + |σλ (t)|E ,
dt λ
we obtain
Z t
|Y x (t)|E ≤ e−k1 t |x|E + e−k1 (t−τ ) (k2 |Zτx |sE + k3 ) dτ,
0
ERGODICITY FOR STOCHASTIC PDE’S 31
and thereby,
Z t
(6.10) E|Y x (t)|E ≤ e−k1 t |x|E + ek1 (t−τ ) (k2 k(s) + k3 ) dτ
0
and (6.2) follows.
where Br := {y ∈ E, |y|E < r}, and c1 , c2 , p and Λ are defined in the same
way as in Theorem 5.3. In the following proposition, existence of a universal
small set satisfying a uniform geometric drift condition is shown.
where
Z −1 Z
(6.16) µ̄(Γ) := e−Λ(y) µ1 (dy) e−Λ(y) µ1 (dy), Γ ∈ B(E),
Br Br ∩Γ
is a probability measure. In particular, Br is a small set of the Markov chain
(X̃n ) := (XnT ), with the lower bound measure δµ̄.
32 B. GOLDYS AND B. MASLOWSKI
(b) We have
(6.17) Ex (|XT |E + 1) ≤ 21 (|x|E + 1) + b1Br (x), x ∈ E,
that is, the chain (X̃n ) satisfies the one-step Lyapunov–Foster condition of
geometric drift toward Br , with the constants 12 and b and the Lyapunov
function V (x) = |x|E + 1.
= δµ̄(Γ), Γ ∈ B(E)
and (6.15) follows.
To prove part (b), we use again (6.11) to obtain
Ex (|Xt |E + 1) ≤ k0 |x|E e−k1 t + cb + 1 ≤ 14 |x|E + cb + 1
(6.21) ≤ 12 (|x|E + 1) − 14 |x|E − 12 + cb + 1
≤ 12 (|x|E + 1) + (cb + 12 )1Br (x)
for x ∈ E, t ≥ − k11 log 4k10 , which completes the proof.
In the next theorem our main result on uniform geometric V -ergodicity for
V (x) = |x|E + 1 is stated. It is based on the paper by Meyn and Tweedie [33],
where exact bounds for geometric ergodicity of irreducible Markov chains are
ERGODICITY FOR STOCHASTIC PDE’S 33
found, and Proposition 6.2 above. Following [33], we introduce the constants
b b
v, Mc , γc , λ, b and ξ̄ as follows:
1/2 + γc
v = r + 1, γc = δ−2 (4b + δv), b=
λ < 1,
1 + γc
(6.22)
b 4 − δ2 2
b = v + γc , ξ̄ = 4b
δ5
and
1
Mc = (1 − λ b + bb 2 + ξ̄( bb(1 − λ
b +b b 2 )) > 1.
b)+b
b 2
(1 − λ)
We will show that the Markov chain (X̃n ) has the geometric rate of conver-
gence to the invariant measure with any constant
1
(6.23) ρ∈ 1− ,1 .
Mc
Let bV B denote the Banach space of measurable functions ϕ : E → R such
that
|ϕ(x)|
kφkV = sup < ∞.
x∈E V (x)
where
1
ω=− log ρ > 0 and
T
(6.25)
ρ
M = (1 + γc ) −1 (cb + k0 + 1)e− log ρ
ρ + Mc − 1
and
(6.26) kPt∗ ν − µ∗ kvar ≤ M (Lν + 1)e−ωt , t ≥ 0, ν ∈ P,
R
where Lν = E |x|E ν(dx). The constants ω and M may be chosen the same
for all nonlinear terms F satisfying Hypothesis 2.4(b) with the same con-
stants K and m and (6.11) with the same constants k0 , k1 and cb [or, in
particular, satisfying the growth condition (2.11) with the same k1 , k2 , k3
and s].
34 B. GOLDYS AND B. MASLOWSKI
If the growth of the nonlinear term F is faster than linear, the Markov
process defined by the equation (2.1) may be uniformly ergodic, that is, the
constant Lν in (6.26) may be replaced by another constant independent of
the initial measure ν ∈ P. This has been established earlier in [14] and [20];
however, the lower bound measures are not found there constructively. In
the theorem below explicit bounds are found and, in particular, uniformity
of convergence with respect to coefficients is proven.
ERGODICITY FOR STOCHASTIC PDE’S 35
Theorem 6.4. Assume (6.1) and let the stronger growth condition (6.3)
hold true. Then there exists an invariant measure µ∗ ∈ P and for any ν ∈ P,
(6.30) kPt∗ ν − µ∗ kvar ≤ (1 − δ)−1 e−b
ωt
kν − µ∗ kvar , t ≥ 0,
where ωb = − 12 log(1 − δ) > 0, δ is defined by (6.14) with R = 2M c and r =
∞ and M c is given by (6.5). In particular, the constants on the r.h.s. of
(6.30) are uniform with respect to all nonlinear terms F satisfying the growth
conditions (2.9) and (6.3) with the same constants K, m, k1 , k2 , k3 , s and ε.
≤ e−2nb
ω
kµkvar ≤ (1 − δ)−1 e−(2n+s)b
ω
kµkvar .
Lemma 7.1. Assume that the equation (2.1) has an invariant measure
µ∗ ∈ P. Then the space Lp (E, µ∗ ) is invariant for (Pt∗ ) for each p ∈ [1, ∞).
Moreover, kPt∗ kp→p = 1 and
Z
(7.2) Pt∗ ψ(y) = pt (x, y)ψ(x)µ∗ (dx), ψ ∈ L1 (E, µ∗ ).
E
Putting φ = 1, we obtain
kGt ψk1 = hPt∗ (ψµ∗ ), 1i = hψµ∗ , 1i = kψk1 ,
and therefore,
kGt ψk1 ≤ kψk1 , ψ ∈ L1 (E, µ∗ ), ψ ≥ 0.
Clearly, Gt ψ = Pt∗ (ψµ∗ ). All those arguments extend immediately to an ar-
bitrary ψ ∈ L1 (E, µ∗ ) and therefore, Gt is a contraction on L1 (E, µ∗ ). Other
parts follow easily by a standard density argument.
Theorem 7.2. Assume (6.1) and let the stronger growth condition (6.3)
be satisfied. Then, for each p ∈ (1, ∞), we have
ω
b
(7.3) gap(Lp ) ≥
p
and
(7.4) kPt φ − hµ∗ , φikp ≤ Cp e−(b
ω /p)t
kφkp ,
b > 0 is defined in Theorem 6.4. If, moreover, the semigroup (Pt ) is
where ω
symmetric in L2 (E, µ∗ ), then (7.3) and (7.4) hold for p = 1.
ERGODICITY FOR STOCHASTIC PDE’S 37
Remark 7.3. (1) In Theorem 7.1 and 7.2 the invariant measure µ∗ ,
hence, the space Lp (E, µ∗ ), depends on the coefficients of equation (2.1). It is
interesting to note that the lower bound on gap(Lp ) and Cp are universal for
all systems satisfying Hypothesis 2.4(b) and (6.3) with the same constants.
(2) By Theorem 7.2, the spectral gap exists for all p ∈ (1, ∞). The fact
that this property holds in L1 (E, µ∗ ) is perhaps surprising. Note that it does
not need to hold in general if F = 0. It is known (cf. [13]) that, for a one-
dimensional Ornstein–Uhlenbeck operator LOU 1 considered in L1 (E, µ∗ ), we
have
σ(LOU
1 ) = {λ : Re λ ≤ 0}.
38 B. GOLDYS AND B. MASLOWSKI
where ρb = e−ωt , that is, the skeleton (Yn ) := (Xtn ) is V -uniformly ergodic
with the rate ρb. By [36], Theorem 2.1, it follows that
Z
n 2
(7.13) kPtn ϕkL2 ≤ ρb kϕkL2 , n ∈ N, ϕ ∈ L , ϕ dµ∗ = 0
8. Some extensions.
Theorem 8.1. Let Hypotheses 2.1, 2.2, 2.3 and condition (6.1) be satis-
fied and let F : R+ × E → E be a jointly measurable mapping such that F (t, ·)
is Lipschitz continuous on bounded sets in E and satisfies Hypothesis 2.4(b)
and the growth condition (6.3) with constants independent of t ∈ R+ . Then
∗
kPs,t ∗
ν1 − Ps,t ν2 kvar ≤ (1 − δ)−1 e−b
ω (t−s)
kν1 − ν2 kvar , 0 ≤ s ≤ t, ν1 , ν2 ∈ P,
where ωb = − 12 log(1 − δ) and δ, ω
b depend only on the constants in Hypothe-
sis 2.4(b) and (6.3).
ERGODICITY FOR STOCHASTIC PDE’S 39
Theorem 8.2. Let Hypotheses 2.1–2.4, condition (6.1) and the growth
condition (2.11) hold for equation (8.1) with the constants independent of
α ∈ A, and assume
(8.2) lim Gα (x) = Gα0 (x), x ∈ E,
α→α0
for each t > 0 and x ∈ E. By (5.2), it suffices to show that exp ρα → exp ρα0
in L1 (Ω), where
Z t Z t
ρα (Z x ) = hGα (Zsx ), dWs i − 21 |Gα (Zsx )|2 ds.
0 s
40 B. GOLDYS AND B. MASLOWSKI
9. Examples.
order to apply Theorem 6.4 (uniform ergodicity) and Theorem 7.2 [spectral
gap in Lp (E, µ∗ )], we have to assume the stronger growth condition (6.3).
As a specific example, we consider a nonlinear stochastic oscillator equa-
tion
(9.1) ÿ = f (y, ẏ) + σ ẇt , y(0) = x1 , ẏ(0) = x2 ,
in Rd . We assume that f : Rd × Rd → Rd is a locally Lipschitz function,
x1 , x2 ∈ Rd , σ ∈ L(Rd ) is a regular matrix, and (wt ) is a standard Wiener
process in Rd . Equation (9.1) may be rewritten in the form (2.1) with Xt =
(y(t), ẏ(t)) ∈ R2d = E = H,
0 I 0
A= , F (x) = , x ∈ R2d ,
0 0 f (x)
1/2 0 0
Q = .
0 (σσ ∗ )1/2
According to the Kalman rank condition (see, e.g., [27]), the matrix Qt
is invertible for each t > 0, so the equation with F = 0 is strongly Feller.
Suppose that f has at most polynomial growth, that is,
|f (x)|Rd ≤ K(1 + |x|m
R2d ), x ∈ R2d ,
for some K, m < ∞. Then by Theorem 6.3, the solution is V -uniformly er-
godic (Theorem 6.3), provided the ultimate boundedness condition (6.11)
holds true. For example, we may take d = 1, f (x) = −α2 x2 − α1 x1 for
x = (x1 , x2 ) ∈ R2d , with some α1 , α2 > 0 (a damped linear oscillator). Then
(6.11) holds with constants which may be easily expressed in terms of α1 , α2
and σ and V -uniform ergodicity holds true. Similar results for a more general
version of equation (9.1) can be found in [30].
Note that the semigroup (Pt ) is not symmetric in this case and Corol-
lary 7.4 (on the spectral gap) is not applicable. Indeed, it follows from the
results in [7, 8] that the spectral gap is zero in the present case. This exam-
ple also shows that the assumption of symmetry of Pt in Corollary 7.4 may
not be removed.
where ai , i = 1, . . . , 2n, are in a given bounded set in R2n , a2n+1 ≥ ā, for
a given ā > 0 and n ≥ 0. Then the V -uniform ergodicity (and if n > 0,
uniform exponential ergodicity) holds with constants in (6.24), (6.26) and
(6.30) uniform with respect to f ∈ Ψ. Also, if n > 0, then there is a positive
lower bound on the spectral gap for (Pt ) in Lp (E, µ∗ ), p ∈ (1, ∞), uniform in
f ∈ Ψ. Finally, if b = 0 and Q = I, then this bound holds also for p = 1.
Example 9.4 (The case of Lipschitz drift). Consider (9.2) with the same
differential operator L and initial and boundary conditions in the case when
the noise may degenerate, for simplicity, suppose that c ≤ c0 < 0. For σ ≥
0, let Hσ denote the domain dom((−A)σ ) equipped with the graph norm
|y|σ := |(−A)σ y|. As well known, for σ ∈ (0, 21 ), the norm | · |σ is equivalent
with the norm of Sobolev–Slobodetskii space H 2σ (0, 1),
Z 1Z 1 |y(ξ) − y(η)|2
|y|2H 2σ := |y|2L2 (0,1) + dξ dη, y ∈ Hσ .
0 0 (ξ − η)1+2σ
Assume that f : R → R is Lipschitz continuous. It is easy to check that F :
Hσ → Hσ is continuous and satisfies the growth condition (2.11) if, for some
k̄1 , k̄2 > 0, we have
(f (ξ) − f (η)) sign(ξ − η) ≤ k̄1 |ξ − η| + k̄2 , ξ, η ∈ R.
Assume that Q = (−A)−2∆ for some ∆ ≥ 0. Then setting E = Hσ , we obtain
|G(x)| = |Q−1/2 F (x)| ≤ kQ−1/2 kL(Hσ ,H) |F (x)σ ≤ K(1
c + |x|σ ), x ∈ Hσ ,
APPENDIX
For the reader’s convenience, we collect here some basic facts about mea-
surable linear mapping that are used in the paper. Most of them are well
known.
Let H be a real separable Hilbert space and let µ = N (0, C) be a centered
Gaussian measure on H with the covariance operator C such that im(C) =
H. The space HC = im(C 1/2 ) endowed with the norm |x|C = |C −1/2 x| can
be identified as the reproducing kernel Hilbert space of the measure µ. In the
sequel we will denote by {en : n ≥ 1} the eigenbasis of C and by {cn : n ≥ 1}
the corresponding set of eigenvalues:
Cen = cn en , n ≥ 1.
For any h ∈ H, we define
n
X 1
φn (x) = √ hh, ek ihx, ek i, x ∈ H.
k=1
ck
Proof. Let
Kx = K(t, s)x, x ∈ H.
By assumption, the operator K : H → L2 ((0, a)×(0, a); H)
is Hilbert–Schmidt
for any a < 1 and thereby, by Lemma A.2, there exists the space Ma of full
measure such that
X∞
−1/2 1
KC y= √ hy, ek iKek ,
k=1
ck
where the convergence holds in mean-square and for each y ∈ Ma , in L2 ((0, a)×
(0, a); H). Therefore, KC −1/2 is a measurable
T function of (y, s, t) for s, t ≤ a.
Let an → a be increasing and let M = ∞ n=1 Man . Then KC
−1/2 y is well
defined, for all s, t < 1 is clearly measurable in (y, s, t). Moreover, for each
y ∈ M,
Z an Z an 2
In2 (y) = |K(t, s)C −1/2 y| ds dt
Z 0
an Z 0
an
[3pt] ≤ |K(t, s)C −1/2 y|2 ds dt < ∞.
0 0
Since the sequence In (y) is nondecreasing to a limit I∞ (y) for each y ∈ Ω
and
Z Z 1Z 1
I∞ (y)µ(dy) ≤ kK(t, s)kHS ds dt < ∞,
H 0 0
the lemma follows.
46 B. GOLDYS AND B. MASLOWSKI
REFERENCES
[1] Aida, S. (1998). Uniform positivity improving property, Sobolev inequalities, and
spectral gaps. J. Funct. Anal. 158 152–185. MR1641566
[2] Aida, S. and Shigekawa, I. (1994). Logarithmic Sobolev inequalities and spectral
gaps: Perturbation theory. J. Funct. Anal. 126 448–475. MR1305076
[3] Albeverio, S., Kondratiev, Y. G. and Röckner, M. (1997). Ergodicity of
L2 -semigroups and extremality of Gibbs states. J. Funct. Anal. 144 394–423.
MR1432591
[4] Arima, R. (1964). On general boundary value problem for parabolic equations. J.
Math. Kyoto Univ. 4 207–243. MR0197997
[5] Cerrai, S. (1999). Ergodicity for stochastic reaction–diffusion systems with polyno-
mial coefficients. Stochastics Stochastics Rep. 67 17–51. MR1717811
[6] Cerrai, S. (2001). Second Order PDE ’s in Finite and Infinite Dimensions. A Prob-
abilistic Approach. Springer, Berlin. MR1840644
[7] Chojnowska-Michalik, A. (2001). Transition Semigroups for Stochastic Semilinear
Equations on Hilbert Spaces. Dissertationes Math. 396. MR1841090
[8] Chojnowska-Michalik, A. and Goldys, B. (2002). Symmetric Ornstein–
Uhlenbeck semigroups and their generators. Probab. Theory Related Fields 124
459–486. MR1942319
[9] Da Prato, G. and Zabczyk, J. (1992). Stochastic Equations in Inifnite Dimensions.
Cambridge Univ. Press. MR1207136
[10] Da Prato, G. and Zabczyk, J. (1996). Ergodicity for Infinite Dimensional Systems.
Cambridge Univ. Press. MR1417491
[11] Da Prato, G., Debussche, A. and Goldys, B. (2002). Some properties of invari-
ant measures of non symmetric dissipative stochastic systems. Probab. Theory
Related Fields 123 355–380. MR1918538
[12] Da Prato, G. and Zabczyk, J. (1992). Nonexplosion, boundedness, and ergod-
icity for stochastic semilinear equations. J. Differential Equations 98 181–195.
MR1168978
[13] Davies, E. B. and Simon, B. (1986). L1 -properties of intrinsic Schrödinger operators.
J. Funct. Anal. 65 126–146. MR0819177
[14] Goldys, B. and Maslowski, B. (2001). Uniform exponential ergodicity of stochastic
dissipative systems. Czechoslovak Math. J. 51 745–762. MR1864040
[15] Goldys, B. and Maslowski, B. (2005). On the Ornstein–Uhlenbeck bridge.
Preprint.
[16] Goldys, B. and Maslowski, B. (2004). Exponential ergodicity for stochastic
reaction-diffusion equations. In Stochastic Partial Differential Equations and
Applications VII. Lecture Notes Pure Appl. Math. 245 115–131. Chapman
Hall/CRC Press, Boca Raton, FL. MR2227225
[17] Goldys, B. and Maslowski, B. (2006). Exponential ergodicity for stochastic Burg-
ers and 2D Navier–Stokes equations. J. Funct. Anal. 226 230–255. MR2158741
[18] Gong, F., Röckner, M. and Wu, L. (2001). Poincaré inequality for weighted first
order Sobolev spaces on loop spaces. J. Funct. Anal. 185 527–563. MR1856276
[19] Hairer, M. (2002). Exponential mixing properties of stochastic PDEs through
asymptotic coupling. Probab. Theory Related Fields 124 345–380. MR1939651
[20] Hairer, M. (2002). Exponential mixing for a PDE driven by a degenerate noise.
Nonlinearity 15 271–279. MR1888852
[21] Hasminskii, R. Z. (1980). Stochastic Stability of Differential Equations. Sijthoff &
Noordhoff, Alphen aan den Rijn. MR0600653
ERGODICITY FOR STOCHASTIC PDE’S 47
[43] Wu, L. (2004). Essential spectral radius for Markov semigroups. I. Discrete time
case. Probab. Theory Related Fields 128 255–321. MR2031227