Lower bounds for weak approximation errors for spatial ... - arXiv

1 downloads 0 Views 468KB Size Report
Jan 16, 2017 - Lemma 2.5 in [4] proves that it holds P-a.s. that ... Observe that Lemma 2.1 implies (i) and (ii). .... Itô's isometry and Lemma 2.6 in [4] imply that.
Lower bounds for weak approximation errors for spatial spectral Galerkin approximations of stochastic wave equations Ladislas Jacobe de Naurois, Arnulf Jentzen, and Timo Welti ETH Zürich, Switzerland

arXiv:1701.04351v1 [math.PR] 16 Jan 2017

January 17, 2017 Abstract Although for a number of semilinear stochastic wave equations existence and uniqueness results for corresponding solution processes are known from the literature, these solution processes are typically not explicitly known and numerical approximation methods are needed in order for mathematical modelling with stochastic wave equations to become relevant for real world applications. This, in turn, requires the numerical analysis of convergence rates for such numerical approximation processes. A recent article by the authors proves upper bounds for weak errors for spatial spectral Galerkin approximations of a class of semilinear stochastic wave equations. The findings there are complemented by the main result of this work, that provides lower bounds for weak errors which show that in the general framework considered the established upper bounds can essentially not be improved.

1

Introduction

In this work we consider numerical approximation processes of solution processes of stochastic wave equations and examine corresponding weak convergence properties. As opposed to strong convergence, weak convergence even in the case of stochastic evolution equations with regular nonlinearities are still only poorly understood (see, e.g., [3, 6, 7, 8, 12] for several weak convergence results for stochastic wave equations and, e.g., the references in Section 1 in [4] for further results on weak convergence in the literature). Therefore, equations available to current numerical analysis are limited to model problems such as the ones considered in the present article that cannot take into account the full complexity of models for evolutionary processes under influence of randomness appearing in real world applications (see, e.g., the references in Section 1 in [4]). The recent article [4] by the authors provides upper bounds for weak errors for spatial spectral Galerkin approximations of a class of semilinear stochastic wave equations, including equations driven by multiplicative noise and, in particular, the hyperbolic Anderson model. The purpose of this work is to show that the weak convergence rates for stochastic wave equations established in Theorem 1.1 in [4] can in the general setting there essentially not be improved. This is achieved by proving lower bounds for weak errors in the case of concrete examples of stochastic wave equations with additive noise and without drift nonlinearity (see Corollary 2.10 below). We argue similarly to the reasoning in Section 6 in Conus et al. [1] and Section 9 in Jentzen & Kurniawan [5]. First results on lower bounds for strong errors for two examples of stochastic heat equations were achieved in Davie & Gaines [2]. Furthermore, lower bounds for strong errors for examples and whole classes of stochastic heat equations have been established in Müller-Gronbach et al. [10] (see also the references therein) and in Müller-Gronbach & Ritter [9], respectively. Results on lower bounds for 1

weak errors in in the case of a few specific examples of stochastic heat equations can be found in Conus et al. [1] and in Jentzen & Kurniawan [5]. Theorem 1.1. For all real numbers η, T ∈ (0, ∞), every R-Hilbert space (H, h·, ·iH , k·kH ), every probability space (Ω, F, P) with a normal filtration (Ft )t∈[0,T ] , every idH -cylindrical (Ω, F, P, (Ft )t∈[0,T ] )Wiener process (Wt )t∈[0,T ] , and every orthonormal basis (en )n∈N : N → H of H there exist an increasing sequence (λn )n∈N : N → (0, D(A) ⊆ H → H with D(A) =  ∞), a linear operator A : P P 2 v ∈ H: |λ he , vi | < ∞ and ∀ v ∈ D(A) : Av = H n∈N n n n∈N −λn hen , viH en , a family of interpolation spaces (Hr , h·, ·iHr , k·kHr ), r ∈ R, associated to −A (cf., e.g., [11, Section 3.7]), a family of R-Hilbert spaces (Hr , h·, ·iHr , k·kHr ), r ∈ R, with ∀ r ∈ R : (Hr , h·, ·iHr , k·kHr ) = S S Hr/2 × Hr/2−1/2 , h·, ·iHr/2 ×Hr/2−1/2 , k·kHr/2 ×Hr/2−1/2 , families of functions PN : r∈R Hr → r∈R Hr , S S N ∈ N ∪ {∞}, and PN : r∈R Hr → r∈R Hr , N ∈ N ∪ {∞}, with ∀ N ∈ N ∪ {∞}, r ∈ R, u ∈  P −r −r Hr , (v, w) ∈ Hr : PN (u) = N h(λ ) e , ui (λ ) e and P (v, w) = (P (v), P (w)) , a linn n H n n N N N r n=1 ear operator A : D(A) ⊆ H0 → H0 with D(A) = H1 and ∀ (v, w) ∈ H1 : A(v, w) = (w, Av), real numbers γ, c ∈ (0, ∞), and functions ξ ∈ L2 (P|F0 ; Hγ ), ϕ ∈ Cb2 (H0 , R), F ∈ Cb2 (H0 , H0 ), B ∈ Cb2 (H0 , HS(H, H0 )) and (Cε )ε∈(0,∞) : (0, ∞) → [0, ∞) with ∀ β ∈ (γ/2, γ] : (−A)−β/2 ∈ HS(H0 ), F ∈ Cb2 (H0 , Hγ ), and B ∈ Cb2 (H0 , L(H, Hγ )) such that (i) it holds that there exist up to modifications unique (Ft )t∈[0,T ] -predictable stochastic processes XN : [0, T ] × Ω → PN (H0 ), N ∈ N ∪ {∞}, which satisfy for all N ∈ N ∪ {∞}, t ∈ [0, T ] that sups∈[0,T ] kXN s kL2 (P;H0 ) < ∞ and P-a.s. that XN t

tA

Z

t (t−s)A

e

= e PN ξ +

PN F(XN s ) ds

0

Z +

t

e(t−s)A PN B(XN s ) dWs

(1.1)

0

(ii) and it holds for all ε ∈ (0, ∞), N ∈ N that     ≤ Cε · (λN )ε−η . c · (λN )−η ≤ E ϕ X∞ − E ϕ XN T T

(1.2)

Here and below we denote for every non-trivial R-Hilbert space (V, h·, ·iV , k·kV ) and every R-Hilbert space (W, h·, ·iW , k·kW ) by Cb2 (V, W ) the set of all globally bounded twice Fréchet differentiable functions from V to W with globally bounded derivatives. Theorem 1.1 is a direct consequence of Theorem 1.1 in [4] (with γ = 2η, β = min{η + ε, 2η}, ρ = 0 in the notation of Theorem 1.1 in [4]) and Corollary 2.10 below (with p = 1/η, δ = 1/2 − η in the notation of Corollary 2.10 below). Inequality (1.2) reveals that the weak convergence rates in Theorem 1.1 in [4] are essentially sharp. More details and further lower bounds for weak approximation errors for stochastic wave equations can be found in Corollary 2.10 below.

2 2.1

Lower bounds for weak errors Setting

Let (H, h·, ·iH , k·kH ) be a separable R-Hilbert space, for every set A let P(A) be the power set of A, let T ∈ (0, ∞), let (Ω, F, P) be a probability space with a normal filtration (Ft )t∈[0,T ] , let (Wt )t∈[0,T ] be an idH -cylindrical (Ω, F, P, (Ft )t∈[0,T ] )-Wiener process, let H ⊆ H be a nonempty orthonormal basis of H, let λ : H → R be a function  with sup ⊆ Ph∈H λh < 0, let2 A : D(A) H → H be the linear operator which satisfies D(A) = v ∈ H : |λ hh, vi | < ∞ and H h∈H h P ∀ v ∈ D(A) : Av = h∈H λh hh, viH h, let (Hr , h·, ·iHr , k·kHr ), r ∈ R, be a family of interpolation spaces associated to −A, let (Hr , h·, ·iHr , k·kHr ), r ∈ R, be the family of R-Hilbert spaces which  satisfies for all r ∈ R that (Hr , h·, ·iHr , k·kHr ) = Hr/2 × Hr/2−1/2 , h·, ·iHr/2 ×Hr/2−1/2 , k·kHr/2 ×Hr/2−1/2 , S S S S let PI : r∈R Hr → r∈R Hr , I ∈ P(H), and PI : r∈R Hr → r∈R Hr , I ∈ P(H), be the functions 2

P which satisfy for all I ∈ P(H),r ∈ R, u ∈ Hr , (v, w) ∈ Hr that PI (u) = h∈I h|λh |−r h, uiHr |λh |−r h and PI (v, w) = PI (v), PI (w) , let A : D(A) ⊆ H0 → H0 be the linear operator which satisfies D(A) = H1 and ∀ (v, w) ∈ H1 : A(v, w) = (w, Av), let µ : H → R be a function which satisfies P |µh |2 < ∞, let B ∈ HS(H, H0 ) be the linear operator which satisfies for all v ∈ H that h∈H |λh | P  Bv = 0, h∈H µh hh, viH h , and let XI = (X I,1 , X I,2 ) : Ω → PI (HR0 ), I ∈ P(H), be random T variables which satisfy for all I ∈ P(H) that it holds P-a.s. that XI = 0 e(T −s)A PI B dWs .

2.2

Lower bounds for the squared norm

Lemma 2.1. Assume the setting in Section 2.1. Then for all I ∈ P(H) it holds P-a.s. that   P    RT 1/2  I,1  µh sin |λ | (T − s) dhh, W i h 1/2 h s H h∈I X 0  . = P  µ |λh |R T (2.1) XI = PI XH =  I,2 1/2 1/2 h X cos |λh | (T − s) dhh, Ws iH |λh | h 1 h∈I

|λh |

/2

0

Proof of Lemma 2.1. Lemma 2.5 in [4] proves that it holds P-a.s. that Z H

T

e

X =

(T −s)A

B dWs =

0

XZ h∈H

T

e(T −s)A Bh dhh, Ws iH

0

!  1/2 −1/2 sin (−A) (T − s) h dhh, W i (−A) s H 0 R  = T 1/2 µ cos (−A) (T − s) h dhh, W i h 0 s H h∈H !  R 1 T /2 µh X sin |λ | (T − s) h dhh, W i 1/2 h s H 0 |λh | =  RT 1 µh 0 cos |λh | /2 (T − s) h dhh, Ws iH h∈H   P    RT 1/2 µh sin |λ | (T − s) dhh, W i h 1 h s H /2 h∈H 0 .  = P  µ |λh |R T  1/2 1/2 h cos |λ | (T − s) dhh, W i |λ | h 1/2 h s H h h∈H 0 X µh

RT

(2.2)

|λh |

Furthermore, Lemma 2.7 in [4] shows for all I ∈ P(H) that it holds P-a.s. that Z

T

H

PI X =

PI e

(T −s)A

Z B dWs =

0

T

e(T −s)A PI B dWs = XI .

(2.3)

0

This and (2.2) complete the proof of Lemma 2.1. Lemma 2.2. Assume the setting in Section 2.1 and let I ∈ P(H). Then (i) it holds that hh, X I,1 iH0 , h ∈ H, is a family of independent centred Gaussian random variables,

1 (ii) it holds that |λh | /2 h, X I,2 H 1 , h ∈ H, is a family of independent centred Gaussian random − /2 variables, and (iii) it holds for all h ∈ H that  1/2 2  sin 2|λ | T |µ | 1 h h Var hh, X I,1 iH0 = 1I (h) T− , 1 |λh | 2 2|λh | /2   1 

 sin 2|λh | /2 T |µh |2 1 1/2 I,2 = 1I (h) Var |λh | h, X T+ , 1 H−1/2 |λh | 2 2|λh | /2  1/2 2  

1 − cos 2|λ | T |µ | 1 h h / 2 . Cov hh, X I,1 iH0 , |λh | h, X I,2 H 1 = 1I (h) 1 − /2 |λh | 4|λh | /2 

3

(2.4) (2.5) (2.6)

Proof of Lemma 2.2. Observe that Lemma 2.1 implies (i) and (ii). It thus remains to prove (iii). Lemma 2.1 implies for all h ∈ H that it holds P-a.s. that Z T  µh 1 I,1 sin |λh | /2 (T − s) dhh, Ws iH , (2.7) hh, X iH0 = 1I (h) 1/2 |λh | 0 Z T

 µh 1/2 1 I,2 |λh | h, X = 1I (h) cos |λh | /2 (T − s) dhh, Ws iH . (2.8) 1/2 H−1/2 |λh | 0 Itô’s isometry hence shows for all h ∈ H that Z   2 |µh |2 T 1 2 I,1 Var hh, X iH0 = E |hh, X iH0 | = 1I (h) sin |λh | /2 (T − s) ds |λh | 0  (2.9) 1 2  sin 2|λh | /2 T |µh | 1 = 1I (h) T− , 1 |λh | 2 2|λh | /2 2 i  h



1/2 1/2 I,2 I,2 = E |λh | h, X Var |λh | h, X H−1/2 H−1/2 Z  2 |µh |2 T 1 = 1I (h) cos |λh | /2 (T − s) ds (2.10) |λh | 0    1 sin 2|λh | /2 T |µh |2 1 = 1I (h) . T+ 1 |λh | 2 2|λh | /2 I,1



Furthermore, observe for all h ∈ H that   h i



1/2 1/2 I,1 I,2 I,1 I,2 Cov hh, X iH0 , |λh | h, X = E hh, X iH0 |λh | h, X H−1/2 H−1/2 Z 2 T   |µh | 1 1 sin |λh | /2 (T − s) cos |λh | /2 (T − s) ds = 1I (h) |λh | 0  2   1/2 |µh |2 sin |λh | T = 1I (h) 1 |λh | 2|λh | /2   1/2 |µh |2 1 − cos 2|λh | T = 1I (h) . 1 |λh | 4|λh | /2

(2.11)

The proof of Lemma 2.2 is thus completed. Corollary 2.3. Assume the setting in Section 2.1 and let I ∈ P(H). Then it holds for all (v, w) ∈ PI (H0 ) that     1 sin 2|λh | /2 T 1 |µh |2 X h hh, viH0 CovOp(X )(v, w) = T− 1/2 0 2 |λh | h∈I 2|λh |     1 1 − cos 2|λh | /2 T

h 1/2 + |λh | h, w H 1 (2.12) 1/2 0 − /2 2|λh |     1 1 − cos 2|λh | /2 T 0 + hh, viH0 1 1 |λh | /2 h 2|λh | /2     1 sin 2|λh | /2 T

0 1/2 + T+ |λh | h, w H 1 ∈ PI (H0 ). 1 1 − /2 |λh | /2 h 2|λh | /2 I

4

Proof of Corollary 2.3. Lemma 2.2 and Lemma 2.1 prove for all x1 = (v1 , w1 ), x2 = (v2 , w2 ) ∈ PI (H0 ) that    hx1 , CovOp(XI )x2 iH0 = Cov hx1 , XI iH0 , hx2 , XI iH0 = E hx1 , XI iH0 hx2 , XI iH0    = E hv1 , X I,1 iH0 + hw1 , X I,2 iH−1/2 hv2 , X I,1 iH0 + hw2 , X I,2 iH−1/2 X  = Var hh, X I,1 iH0 hh, v1 iH0 hh, v2 iH0 h∈H

+

X



1 Cov hh, X I,1 iH0 , |λh | /2 h, X I,2 H

 −1/2

1 hh, v1 iH0 |λh | /2 h, w2 H

−1/2

1 hh, v2 iH0 |λh | /2 h, w1 H

−1/2

h∈H

+

X



1 Cov hh, X I,1 iH0 , |λh | /2 h, X I,2 H

 −1/2

h∈H

+

X

Var



|λh | /2 h, X I,2 1



H−1/2

|λh | /2 h, w1 1



H−1/2

|λh | /2 h, w2 1

H−1/2

h∈H

 =

v1 ,

Xh

 Var hh, X I,1 iH0 hh, v2 iH0

h∈H



+ Cov hh, X

I,1



1 iH0 , |λh | /2 h, X I,2 H

−1/2

  Xh

1 Cov hh, X I,1 iH0 , |λh | /2 h, X I,2 H + w1 ,



1 |λh | /2 h, w2 H 

−1/2

−1/2

i  h

(2.13)

H0

hh, v2 iH0

h∈H



1 + Var |λh | /2 h, X I,2 H



−1/2

1/2

|λh | h, w2

i

H−1/2

 |λh | h 1/2

H−1/2

 Xh  Var hh, X I,1 iH0 hh, v2 iH0 = x1 , h∈H



+ Cov hh, X

I,1



1/2

iH0 , |λh | h, X

I,2



H−1/2

 Xh 

1 Cov hh, X I,1 iH0 , |λh | /2 h, X I,2 H + x1 ,

−1/2



1/2

|λh | h, w2



i h

H−1/2

0

H0

hh, v2 iH0

h∈H



1 + Var |λh | /2 h, X I,2 H



−1/2

1/2

|λh | h, w2

H−1/2

i

 0 . 1 |λh | /2 h H0

This and again Lemma 2.2 complete the proof of Corollary 2.3. Lemma 2.4. Assume the setting in Section 2.1 and let I ∈ P(H). Then it holds for all i ∈ {1, 2} that XI ∈ L2 (P; H0 ) and X |µh |2  I 2  E kX kH0 = T < ∞, |λh | h∈I  1/2 h i 1 X |µ |2  sin 2|λ | T h h 2 < ∞. E kX I,i kH1 i = T+ 1 /2− /2 2 h∈I |λh | (−1)i 2|λh | /2 Proof of Lemma 2.4. Itô’s isometry and Lemma 2.6 in [4] imply that

2   Z T X |µh |2

 I 2  2 (T −s)A

E kX kH0 = E e PI B dWs = T kP Bk = T < ∞. I HS(H,H0 )

|λ h| 0 H0 h∈I

5

(2.14) (2.15)

(2.16)

In addition, Lemma 2.2 shows for all i ∈ {1, 2} that i X h

h i 1 E kX I,i k2H1 i = E |λh | /2− /2 h, X I,i H1

/2−i/2

/2− /2

2 i

h∈H

  1 sin 2|λh | /2 T 1 X |µh |2 = < ∞. T+ 1 2 h∈I |λh | (−1)i 2|λh | /2

(2.17)

The proof of Lemma 2.4 is thus completed. Proposition 2.5. Assume the setting in Section 2.1 and assume inf h∈H |µh | > 0. Then it holds for all I ∈ P(H) \ {H} that X 1       E kXH k2H0 − E kXI k2H0 = E kXH\I k2H0 ≥ T inf |µh |2 > 0. (2.18) h∈H |λh | h∈H\I

Proof of Proposition 2.5. Orthogonality and Lemma 2.1 imply for all I ∈ P(H) \ {H} that         E kXI k2H0 + E kXH\I k2H0 = E kPI XH k2H0 + E kPH\I XH k2H0 (2.19)     = E k(PI + PH\I )XH k2H0 = E kXH k2H0 . This and Lemma 2.4 show for all I ∈ P(H) \ {H} that       E kXH k2H0 − E kXI k2H0 = E kXH\I k2H0 X |µh |2 X 1 =T ≥ T inf |µh |2 > 0. h∈H |λh | |λh | h∈H\I

(2.20)

h∈H\I

The proof of Proposition 2.5 is thus completed. In Corollary 2.7 and Corollary 2.8 below lower bounds on the weak approximation error with the squared norm as test function are presented. Our proofs of Corollary 2.7 and Corollary 2.8 use the following elementary and well-known lemma (cf., e.g., Proposition 6.4 in Conus et al. [1]). Lemma 2.6. Let p ∈ (0, ∞), δ ∈ (−∞, 1/2 − 1/(2p)). Then it holds for all N ∈ N that ∞ X

np(2δ−1) ≥

n=N +1

N p(2δ−1)+1 . [p(1 − 2δ) − 1]2p(1−2δ)−1

(2.21)

Proof of Lemma 2.6. Observe that the assumption that δ ∈ (−∞, 1/2 − 1/(2p)) ensures that p(2δ − 1) ∈ (−∞, −1). This implies for all N ∈ N that Z n+1 Z n+1 ∞ ∞ ∞ X X X p(2δ−1) p(2δ−1) n = n dx ≥ xp(2δ−1) dx n=N +1

n=N +1 ∞

Z =

n

n=N +1

xp(2δ−1) dx = −

N +1

n

(N + 1)p(2δ−1)+1 p(2δ − 1) + 1

(2.22)

N p(2δ−1)+1 ≥ . [p(1 − 2δ) − 1]2p(1−2δ)−1 This completes the proof of Lemma 2.6. Corollary 2.7. Assume the setting in Section 2.1, let c ∈ (0, ∞), p ∈ (1, ∞), let e : N → H be a bijection which satisfies for all n ∈ N that λen = −cnp , and let IN ∈ P(H), N ∈ N, be the sets which satisfy for all N ∈ N that IN = {e1 , e2 , . . . , eN } ⊆ H. Then it holds for all N ∈ N that     T inf h∈H |µh |2 N 1−p E kXH k2H0 − E kXIN k2H0 ≥ . c(p − 1)2p−1 6

(2.23)

Proof of Corollary 2.7. Proposition 2.5 and Lemma 2.6 prove for all N ∈ N that X     E kXH k2H0 − E kXIN k2H0 ≥ T inf |µh |2 h∈H

h∈H\IN

∞ X 1 1 2 −1 = c T inf |µh | h∈H |λh | np n=N +1

T inf h∈H |µh |2 N 1−p ≥ . c(p − 1)2p−1

(2.24)

The proof of Corollary 2.7 is thus completed. Corollary 2.8. Assume the setting in Section 2.1, let c, p ∈ (0, ∞), δ ∈ (−∞, 1/2 − 1/(2p)), let e : N → H be a bijection which satisfies for all n ∈ N that λen = −cnp , let IN ∈ P(H), N ∈ N, be the sets which satisfy for all N ∈ N that IN = {e1 , e2 , . . . , eN } ⊆ H, and assume for all h ∈ H that |µh | = |λh |δ . Then it holds for all N ∈ N that T c2δ−1 N p(2δ−1)+1 . [p(1 − 2δ) − 1]2p(1−2δ)−1

    E kXH k2H0 − E kXIN k2H0 ≥

(2.25)

Proof of Corollary 2.8. Proposition 2.5, Lemma 2.4, and Lemma 2.6 show for all N ∈ N that X X |µh |2     E kXH k2H0 − E kXIN k2H0 = T =T |λh |2δ−1 |λh | h∈H\IN

h∈H\IN

= Tc

2δ−1

∞ X

p(2δ−1)

n

n=N +1

T c2δ−1 N p(2δ−1)+1 . ≥ [p(1 − 2δ) − 1]2p(1−2δ)−1

(2.26)

This completes the proof of Corollary 2.8.

2.3

Lower bounds for the weak error of a particular regular test function

The next proposition, Proposition 2.9 below, follows directly from Lemma 2.2 above and Lemma 9.5 in Jentzen & Kurniawan [5]. Proposition 2.9. Assume the setting in Section 2.1 and let ϕi : H0 → R, i ∈ {1, 2}, be the functions which satisfy for all i ∈ {1, 2}, (v1 , v2 ) ∈ H0 that ϕi (v1 , v2 ) = exp −kvi k2H1 i . Then /2− /2 it holds for all i ∈ {1, 2}, I ∈ P(H) that ϕi ∈ Cb2 (H0 , R) and i i h h I,i 2 H,i 2 E kX kH1 i − E kX kH1 i  /2−h/2 i /2− /2 . (2.27) E[ϕi (XI )] − E[ϕi (XH )] ≥ 2 exp 6 E kX H,i kH1 i /2− /2

Corollary 2.10. Assume the setting in Section 2.1, let c, p ∈ (0, ∞), δ ∈ (−∞, 1/2 − 1/(2p)), let e : N → H be a bijection which satisfies for all n ∈ N that λen = −cnp , let IN ∈ P(H), N ∈ N, be the sets which satisfy for all N ∈ N that IN = {e1 , e2 , . . . , eN } ⊆ H, assume for all h ∈ H that |µh | = |λh |δ , and let ϕi : H0 → R, i ∈ {1, 2}, be the functions which satisfy for all i ∈ {1, 2}, (v1 , v2 ) ∈ H0 that ϕi (v1 , v2 ) = exp −kvi k2H1 i . Then it holds for all i ∈ {1, 2}, N ∈ N that /2− /2 ϕi ∈ Cb2 (H0 , R) and   sin(x) T c2δ−1 2p(2δ−1) N p(2δ−1)+1 IN H E[ϕi (X )] − E[ϕi (X )] ≥ 1 + inf 2δ−1  > 0. (2.28) i x∈[2c1/2 T,∞) (−1) x [p(1 − 2δ) − 1] exp 6p(2δ−1)T c p(2δ−1)+1

7

< 1 prove Proof of Corollary 2.10. Lemma 2.4, Lemma 2.6, and the fact that ∀ x ∈ (0, ∞) : sin(x) x for all i ∈ {1, 2}, N ∈ N that h E kX H,i k2H1

i

h − E kX I,i k2H1

/2−i/2

/2−i/2

i

  1 sin 2|λh | /2 T 1 X |µh |2 = T+ 1 2 |λh | (−1)i 2|λh | /2 h∈H\IN

  1 sin 2|λh | /2 T T X ≥ 1 + inf |λh |2δ−1 1/2 h∈H (−1)i 2|λ | T 2 h h∈H\IN  2δ−1 X  ∞ sin(x) T c np(2δ−1) ≥ 1+ inf i 2 n=N +1 x∈[2c1/2 T,∞) (−1) x   sin(x) T c2δ−1 2p(2δ−1) N p(2δ−1)+1 ≥ 1+ inf > 0. i [p(1 − 2δ) − 1] x∈[2c1/2 T,∞) (−1) x

(2.29)

Furthermore, the assumption that δ ∈ (−∞, 1/2 − 1/(2p)) ensures that p(2δ − 1) ∈ (−∞, −1). Hence, we obtain that Z ∞ ∞ ∞ Z n+1 X X p(2δ−1) p(2δ−1) n ≤1+ x dx = 1 + xp(2δ−1) dx n 1 n=1 n=1 (2.30) 1 p(2δ − 1) =1− = . p(2δ − 1) + 1 p(2δ − 1) + 1 Lemma 2.4 and (2.30) imply for all i ∈ {1, 2} that  h exp −6 E kX H,i k2H1

i

/2−i/2

  ∞ X  H 2  2δ−1 p(2δ−1) ≥ exp −6 E kX kH0 = exp −6T c n n=1

  6p(2δ − 1)T c2δ−1 ≥ exp − > 0. p(2δ − 1) + 1

(2.31)

Combining this and (2.29) with Proposition 2.9 concludes the proof of Corollary 2.10. Corollary 2.11 below specifies Corollary 2.10 to the case where the linear operator A : D(A) ⊆ H → H in the setting in Section 2.1 is the Laplacian with Dirichlet boundary conditions on H. It is an immediate consequence of Corollary 2.10. Corollary 2.11. Assume the setting in Section 2.1, let δ ∈ (−∞, 1/4), let e : N → H be a bijection which satisfies for all n ∈ N that λen = −π 2 n2 , let IN ∈ P(H), N ∈ N, be the sets which satisfy for all N ∈ N that IN = {e1 , e2 , . . . , eN } ⊆ H, assume for all h ∈ H that |µh | = |λh |δ , and let ϕi : H0 → R, i ∈ {1, 2}, be the functions which satisfy for all i ∈ {1, 2}, (v1 , v2 ) ∈ H0 that ϕi (v1 , v2 ) = exp −kvi k2H1 i . Then it holds for all i ∈ {1, 2}, N ∈ N that ϕi ∈ Cb2 (H0 , R) and /2− /2   T (4π 2 )2δ−1 N 4δ−1 sin(x) IN H (2.32) E[ϕi (X )] − E[ϕi (X )] ≥ 1 + inf 4δ−2  > 0. x∈[2πT,∞) (−1)i x [1 − 4δ] exp 12(2δ−1)T π 4δ−1

Acknowledgements This project has been partially supported through the ETH Research Grant ETH-47 15-2 “Mild stochastic calculus and numerical approximations for nonlinear stochastic evolution equations with Lévy noise”.

8

References [1] Conus, D., Jentzen, A., and Kurniawan, R. Weak convergence rates of spectral Galerkin approximations for SPDEs with nonlinear diffusion coefficients. ArXiv e-prints (Aug. 2014). arXiv:1408.1108 [math.PR]. Minor revision requested from Ann. Appl. Probab. [2] Davie, A. M., and Gaines, J. G. Convergence of numerical schemes for the solution of parabolic stochastic partial differential equations. Math. Comp. 70, 233 (2001), 121–134. issn: 0025-5718. doi: 10.1090/S0025-5718-00-01224-2. url: http://dx.doi.org/10. 1090/S0025-5718-00-01224-2. [3] Hausenblas, E. Weak approximation of the stochastic wave equation. J. Comput. Appl. Math. 235, 1 (2010), 33–58. issn: 0377-0427. doi: 10 . 1016 / j . cam . 2010 . 03 . 026. url: http://dx.doi.org/10.1016/j.cam.2010.03.026. [4] Jacobe de Naurois, L., Jentzen, A., and Welti, T. Weak convergence rates for spatial spectral Galerkin approximations of semilinear stochastic wave equations with multiplicative noise. ArXiv e-prints (Aug. 2015). arXiv:1508.05168 [math.PR]. [5] Jentzen, A., and Kurniawan, R. Weak convergence rates for Euler-type approximations of semilinear stochastic evolution equations with nonlinear diffusion coefficients. ArXiv eprints (Jan. 2015). arXiv:1501.03539 [math.PR]. [6] Kovács, M., Lindner, F., and Schilling, R. L. Weak convergence of finite element approximations of linear stochastic evolution equations with additive Lévy noise. ArXiv eprints (Nov. 2014). arXiv:1411.1051 [math.PR]. [7] Kovács, M., Larsson, S., and Lindgren, F. Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise. BIT 52, 1 (2012), 85–108. issn: 0006-3835. doi: 10.1007/s10543-011-0344-2. url: http://dx.doi.org/ 10.1007/s10543-011-0344-2. [8] Kovács, M., Larsson, S., and Lindgren, F. Weak convergence of finite element approximations of linear stochastic evolution equations with additive noise II. Fully discrete schemes. BIT 53, 2 (2013), 497–525. issn: 0006-3835. [9] Müller-Gronbach, T., and Ritter, K. Lower bounds and nonuniform time discretization for approximation of stochastic heat equations. Found. Comput. Math. 7, 2 (2007), 135– 181. issn: 1615-3375. doi: 10.1007/s10208-005-0166-6. url: http://dx.doi.org/10. 1007/s10208-005-0166-6. [10] Müller-Gronbach, T., Ritter, K., and Wagner, T. Optimal pointwise approximation of infinite-dimensional Ornstein-Uhlenbeck processes. Stoch. Dyn. 8, 3 (2008), 519–541. issn: 0219-4937. doi: 10 . 1142 / S0219493708002433. url: http : / / dx . doi . org / 10 . 1142 / S0219493708002433. [11] Sell, G. R., and You, Y. Dynamics of evolutionary equations. Vol. 143. Applied Mathematical Sciences. Springer-Verlag, New York, 2002, xiv+670. isbn: 0-387-98347-3. doi: 10.1007/978-1-4757-5037-9. url: http://dx.doi.org/10.1007/978-1-4757-5037-9. [12] Wang, X. An exponential integrator scheme for time discretization of nonlinear stochastic wave equation. J. Sci. Comput. 64, 1 (2015), 234–263. issn: 0885-7474. doi: 10 . 1007 / s10915-014-9931-0. url: http://dx.doi.org/10.1007/s10915-014-9931-0.

9

Suggest Documents