arxiv:math/ v2 [math.pr] 9 Mar 2007

Transcription

1 arxiv:math/ v2 [math.pr] 9 Mar 2007 Central limit theorems for multiple stochastic integrals and Malliavin calculus D. Nualart and S. Ortiz-Latorre November 2, 2018 Abstract We give a new characterization for the convergence in distribution to a standard normal law of a sequence of multiple stochastic integrals of a fixed order with variance one, in terms of the Malliavin derivatives of the sequence. We also give a new proof of the main theorem in [7] using techniques of Malliavin calculus. Finally, we extend our result to the multidimensional case and prove a wea convergence result for a sequence of square integrable random variables. KEY WORDS: Multiple stochastic integrals. Limit theorems. Gaussian processes. Malliavin calculus. Wea convergence. RUNNING HEAD: Central limit theorems for multiple stochastic integrals MSC2000: 60F05, 60G15, 60H05, 60H07. 1 Introduction Consider a sequence of random variables F belonging to the nth Wiener chaos, n 2, and with unit variance. In [7], Nualart and Peccati have proved that this sequence converges in distribution to a normal N (0,1) law if and only if one of the following two equivalent conditions hold: i) lim E(F 4 ) = 3, ii) lim f l f = 0, for all 1 l n 1, where f is the square integrable ernel associated with the random variable F, and f l f denotes the contraction of l indices of both ernels. In a subsequent paper, Peccati and Tudor [10] gave a multidimensional version of this characterization. Department of Mathematics, University of Kansas, 405 Snow Hall, Lawrence, KS 66045, USA. nualart@math.u.edu Departament de Probabilitat, Lògica i Estadística, Universitat de Barcelona, Gran Via 585, Barcelona, Spain. sortiz@ub.edu 1

2 There have been different extensions and applications of these results. In [3] Hu and Nualart have applied this characterization to establish the wea convergence of the renormalized self-intersection local time of a fractional Brownian motion. In two recent papers, Peccati and Taqqu [9, 10] study the stable convergence of multiple stochastic integrals to a mixture of normal distributions. The aim of this paper is to provide an additional necessary and sufficient condition for the convergence of the sequence F to a normal distribution, in terms of the derivative of F in the sense of Malliavin calculus. This new condition is DF 2 L 2 (Ω) H n. (1) + On the other hand, we give a simple proof of the fact that condition (1) implies the converge in distribution to the normal law based on Malliavin calculus. The main ingredient of the proof is the identity δd = L, where δ, D and L are the basic operators in Malliavin calculus. In this way, we are able to show the result of Nualart and Peccati with the additional equivalent hypotheses (1), without using the Dambis-Dubins-Schwartz characterization of continuous martingales as a Brownian motion with a time change. We also discuss the extension of these results to the multidimensional case. In the last section we study the wea convergence of a sequence of centered square integrable random variables. The result assume conditions on the Malliavin derivatives of the chaotic projections of the sequence {F } N not on the derivatives of the sequence {F } N. Therefore, it can be used with non regular random variables in the Malliavin sense. In [8] Peccati and Taqqu also provide a condition for convergence to the normal law, involving projections of Malliavin derivatives. This type of condition is different from ours, and it is inspired on Clar-Ocone s formula. Condition (1) is a useful tool in establishing the central limit theorem for sequences of random variables defined in terms of a fixed function of a Gaussian process. We apply this approach to derive the wea convergence of the normalized sums of odd powers of the increments of a fractional Brownian motion. In this case, condition (1) follows easily from the ergodic theorem. The paper is organized as follows. In Section 2 we introduce some notation and preliminary results. In section 3 we state and prove the main result of the paper. Section 4 deals with the multidimensional version of the result proved in Section3, andinsection5weapplythepreviousresultstotheweaconvergence of a sequence of centered and square integrable random variables. Finally, in Section 6 we discuss the application of our approach to an example related to the fractional Brownian motion. 2 Preliminaries and notation Let H be a separable Hilbert space. For every n 1 let H n be the nth tensor product of H and denote by H n the nth symmetric tensor product of H, endowed with the modified norm n! H n. Suppose that X = {X(h) : h H} 2

3 is an isonormal Gaussian process on H. This means that X is a centered Gaussian family of random variables indexed by the elements of H, defined on some probability space (Ω,F,P), and such that, for every h,g H, E[X (h)x(g)] = h,g H. We will assume that F is generated by X. Foreveryn 1, leh n bethenthwienerchaosofx,thatis, theclosedlinear subspace of L 2 (Ω,F,P) generated by the random variables {H n (X(h)),h H, h H = 1}, where H n is the nth Hermite polynomial. We denote by H 0 the space of constant random variables. For n 1, the mapping I n (h n ) = n!h n (X(h)) provides a linear isometry between H n and H n. For n = 0, H 0 = R, and I 0 is the identity map. It is well nown (Wiener chaos expansion) that L 2 (Ω,F,P) can be decomposed into the infinite orthogonal sum of the spaces H n. Therefore, any square integrable random variable F L 2 (Ω,F,P) has the following expansion F = I n (f n ), n=0 where f 0 = E[F], and the f n H n are uniquely determined by F. For every n 0 we denote by J n the orthogonal projection on the nth Wiener chaos H n, so I n (f n ) = J n (F). Let {e, 1}be acompleteorthonormalsystemin H. Givenf H n and g H m, for l = 0,...,n m the contraction of f and g of order l is the element of H (n+m 2l) defined by f l g = i 1,...,i l =1 f,e i1 e il H l g,e i1 e il H l. We denote by f l g its symmetrization. Then, f 0 g = f g equals to the tensor product of f and g, and for n = m, f n g = f,g H n. Let us introduce some basic facts on the Malliavin calculus with respect the Gaussian process X. We refer the reader to Nualart [6] for a complete presentation of these notions. Consider the set S of smooth random variables S of the form F = f (X(h 1 ),...,X (h n )), (2) where h 1,...,h n H, f C b (R n ) (the space of bounded functions which have bounded derivatives of all orders) and n N. The derivative operator D on a smooth random variable of the form (2) is defined by DF = n i=1 f x i (X (h 1 ),...,X(h n ))h i, which is an element of L 2 (Ω;H). By iteration one can define D m F which is an element of L 2 (Ω;H m ). For m 1 we denote by D m,2 the completion of S 3

4 with respect to the norm F m,2 given by F 2 m,2 = E[ F 2] + m E[ D i F 2 ]. H i We denote by δ the adjoint of the operator D. That is, δ is an unbounded operator on L 2 (Ω;H) with values in L 2 (Ω), whose domain, denoted by Domδ, is the set of H-valued square integrable random variables u L 2 (Ω;H) such that E[ DF,u H ] c F L2 (Ω), for all F D 1,2. If u belongs to Domδ, then δ(u) is the element of L 2 (Ω) characterized by E[Fδ(u)] = E[ DF,u H ], for any F D 1,2. The operator L defined on the Wiener chaos expansion as L = n=0 nj n is called the infinitesimal generator of the Ornstein-Uhlenbec semigroup. The domain of this operator is the set DomL = {F L 2 (Ω) : i=1 n 2 J n F 2 L 2 (Ω) < + } = D2,2. n=1 The next proposition explains the relationship between the operators D, δ and L. Proposition 1 For F L 2 (Ω) the statement F DomL is equivalent to F DomδD (i.e., F D 1,2 and DF Domδ), and in this case δdf = LF. In the particular case where H = L 2 (A,A,µ), (A,A) is a measurable space, and µ is a σ-finite and non-atomic measure, then H n = L 2 s(a n,a n,µ n ) is the space of symmetric and square integrable functions on A n and for every f H n, I n (f)isthe multiple Wiener-Itô integral (ofordern) off withrespect to X, as defined by Itô in [4]. In this case, F = n=0 I n(f n ) D 1,2 if and only if E[ DF 2 H ] = n f n 2 H n < +, n=1 and its derivative can be identified as the element of L 2 (A Ω) given by D t F = ni n 1 (f n (,t)). (3) n=1 We need the following technical lemma. 4

5 Lemma 2 Consider two random variables F = I n (f), G = I m (g), where n,m 1. Then [ ] n m E DF,DG 2 (n!m!) 2 H = r=1 ((n r)!(m r)!(r 1)!) 2 f r g 2. (4) H (n+m 2r) Proof. Without loss of generality we can assume that H = L 2 (A,A,µ), (A,A) is a measurable space, and µ is a σ-finite and non-atomic measure. In that case, (3) implies that and we have D t F = ni n 1 (f (,t)),d t G = mi m 1 (g(,t)) DF,DG 2 H = nm A I n 1 (f (,t))i m 1 (g(,t))µ(dt). Thans to the multiplication formula for multiple stochastic integrals, see, for instance, Proposition in [6], one obtains DF,DG 2 H = nm A n m 1 r=0 ( )( n 1 m 1 r! r r ) I n+m 2 2r (f (,t) r g(,t))µ(dt). Taing into account the orthogonality between multiple stochastic integrals of different order, we have E[ DF,DG 2 H] n m 1 ( ) 2 ( ) 2 = n 2 m 2 (r!) 2 n 1 m 1 r r r=0 f (,t) r g(,t),f (,s) r g(,s) H (n+m 2 2r)µ(dt)µ(ds). A 2 Notice that A and, as a consequence, Therefore, n m 1 E[ DF,DG 2 H ] = n2 m 2 A f (,t) r g(,t)µ(dt) = f r+1 g, f (,t) r g(,t)µ(dt) = f r+1 g. r=0 which implies the desired result. (r!) 2 ( n 1 r ) 2 ( m 1 r ) 2 f r+1g 2 H (n+m 2 2r), 5

6 3 Main result Fix n 2, n N. Consider a sequence {F } N of square integrable random variables belonging to the nth Wiener chaos. We now that E[ DF 2 H ] = n f 2 H n. (5) The next lemma establish the equivalence between the convergence in L 2 (Ω) of DF 2 H to a constant, and the convergence of E[ DF 4 H ] to the square of the same constant. Lemma 3 Consider a sequence {F = I n (f )} N of square integrable random variables belonging to the nth Wiener chaos such that E [ F 2 ] = f 2 H n + 1. The following conditions are equivalent: 1. lim + E[ DF 4 H ] = n2. 2. lim + DF 2 H = n, in L2 (Ω). Proof. Notice that, using (5), E[( DF 2 H n)2 ] = E[ DF 4 H ] 2nE[ DF 2 H ]+n2 = E[ DF 4 H ] 2n2 f 2 H n +n2, and the result follows easily. Now, we establish the main result of this paper. Theorem 4 Consider a sequence {F = I n (f )} N of square integrable random variables belonging to the nth Wiener chaos such that The following statements are equivalent. E [ F] 2 = f 2 H 1. (6) n + i) As goes to infinity, the sequence {F } N converges in distribution to the normal law N (0,1). ii) lim + E[F 4 ] = 3. iii) For all 1 l n 1, lim + f l f H 2(n l) = 0. iv) DF 2 H L 2 (Ω) + n. 6

7 Proof. We will prove the following implications iv) i) ii) iii) iv). [iv) i)] The sequence of random variables {F } N is tight because it is bounded in L 2 (Ω) by condition (6). Then, by Prohorov s Theorem we have that {F } N is relatively compact, and it suffices to show that the limit of any subsequence converging in distribution is N(0, 1). Suppose that, for a subsequence { l } l N {} N we have F l L l + G. (7) By condition (6) G L 2 (Ω). Therefore, the characteristic function ϕ(t) = E[e itg ] is differentiable and ϕ (t) = ie[ge itg ]. For every N, define ϕ (t) = E[e itf ]. We have ϕ (t) = ie[f e itf ]. By the Continuous Mapping Theorem, (7) implies that F l e itf l L l + Ge itg. (8) The boundedness in L 2 (Ω) plus the convergence in law (8) imply convergence of the expectations. Hence, we obtain ϕ l (t) l + ϕ (t). On the other hand, using the definition of the operator L, Proposition 1 and the definition of the operator δ, we have Therefore, E[F e itf ] = 1 n E[LF e itf ] = 1 n E[ δd(f )e itf ] = 1 n E[ DF,D ( e itf ) H ] = it n E[eitF DF 2 H ]. ϕ l (t) = t n E[eitF l DF l 2 H ]. Furthermore, E[e itf l DF l 2 H ] nϕ(t) DFl E[ 2 ]+n H n E[e itf l ] ϕ(t), which, by the definition of ϕ and hypothesis iv), gives that ϕ l (t) l + tϕ(t). This implies that ϕ(t) satisfies the following differential equation ϕ (t) = tϕ(t) ϕ(0) = 1, 7

8 which is the differential equation satisfied by the characteristic function of the N (0,1). [i) ii)] and [ii) iii)] These implications are proved by Nualart and Peccati in Proposition3, [7]. The proof of the first one is trivial and the proof of the second one involves some combinatorics and it is based on the product formula for multiple stochastic integrals. [iii) iv)]bylemma3itisenoughtoprovethatiii)implieslim + E[ DF 4 H ] = n 2. Using (4) we obtain n 1 E[ DF 4 H ] = (n!) 4 ((n r)!) 4 ((r 1)!) 2 f r f 2 H +n 2 (n!) 2 f 2(n r) 4 H n. r=1 It follows that E[ DF 4 H ] converges to n2 if and only if f l f 2 0, 1 l n 1. H 2(n l) + As f l f 2 H 2(n l) = (2(n l))! f l f 2 H 2(n l) (2(n l))! f l f 2 H 2(n l), we conclude the proof. As a consequence, we obtain. Corollary 5 Fix n 2 and F belonging to the nth Wiener chaos such that E[F 2 ] = 1. Then the distribution of F cannot be normal and E[ DF 4 H ] n2. Proof. If F had a normal distribution or E[ DF 4 H ] = n2, then, according to Theorem 4 and Lemma 3, we would have Var[ DF 2 H ] = 0, but this implies F = 0 or F belonging to the first chaos. 4 Multidimensional case In this section we give a multidimensional version of Theorem 4. For d 2, fix d natural numbers 1 n 1 n d. Consider a sequence of random vectors of the form F = ( F 1,...,Fd ) = ( In1 ( f 1 ),...,Ind ( f d )), (9) where f i H ni, and γ = (γ i,j ) 1 i,j d = ( DF i,df j H) 1 i,j d. The following lemma shows that the convergence of the covariance matrix of F to a diagonal matrix plus the convergencein L 2 (Ω) of the diagonal elements of γ to the constant n i, implies the convergence in L 2 (Ω) of γ to a diagonal matrix. 8

9 Lemma 6 Let {F } N be a sequence of random vectors as (9) such that, for every 1 i,j d, lim + E[Fi Fj ] = δ ij, (10) where δ ij is the Kronecer symbol. We have that implies DF i 2 L 2 (Ω) H n i,1 i d (11) + γ ij L 2 (Ω) ni n j δ ij,1 i,j d. + Proof. We need to show that, for i < j, one has Using (4) we obtain r=1 lim + E[ DFi,DF j 2 H] = 0. n i E[ DF i,dfj 2 H ] = (n i!n j!) 2 f i r=1 ((n i r)!(n j r)!(r 1)!) 2 r f j 2 H (n i +n j 2r) n i (n i!n j!) 2 (n i +n j 2r)! f i ((n i r)!(n j r)!(r 1)!) 2 r f j 2. H (n i +n j 2r) We have reduced the problem to show that lim f i r f j 2 = 0, 1 r n i. H (n i +n j 2r) + The next step is to relate the norm of f i r f j with the norms of fi n i r f i and f j n j r f j. Using the definition of the contractions, we have f i r f j 2 = H fi (n i +n j 2r) ni r f,f i j n j r f j H (2r). (12) Hence, by Cauchy-Schwarz s inequality, we obtain f i r f j 2 f i H (n i +n j 2r) n i r f i H (2r) f j n j r f j H (2r). (13) In the case 1 r n i 1, by assumption (11) and Theorem 4 (implication iv) iii)), the right hand side of equation (13) tends to zero as tends to infinity. In the case r = n i < n j, we have that the right hand side of equation (13) is equal to f i 2 H n i f j n j r f j which tends to zero as tends to infinity, because 2 < +, H n i sup 1 f i 9, H (2r)

10 thans to assumption (10) and, analogously to the previous case, f j n j r f j 0. H (2r) + In the case r = n i = n j, the equality (10) gives ( ) 2 f i r f j 2 = E[F ifj ], H (n i +n j 2r) n i! which tends to zero by assumption (10). Let V d be the set of all (i 1,i 2,i 3,i 4 ) (1,...,d) 4, such that one of the following conditions is satisfied: (a) i 1 i 2 = i 3 = i 4, (b) i 1 i 2 = i 3 i 4 and i 4 i 1, (c) the elements of (i 1,i 2,i 3,i 4 ) are all distinct. The following is a multidimensional version of Theorem 4. Theorem 7 Let {F } N be a sequence of random vectors of the form (9) such that, for every 1 i,j d, lim + E[Fi Fj ] = δ ij, (14) where δ ij is the Kronecer symbol. The following statements are equivalent. i) For every j = 1,...,d, F i variable. converges in distribution to a standard Gaussian ii) For every i = 1,...,d, lim + E[(F i )4 ] = 3. iii) For all 1 i d,1 l n i 1, f i l f i 2 H 2(n i l) iv) For all 1 i d, DF i 2 H L 2 (Ω) + n i. v) For every (i 1,i 2,i 3,i 4 ) V d, ( d lim E + i=1 F i ) 4 = 3d 2, 0. + and 4 lim E[ F i l + ] = 0. l=1 vi) As goes to infinity the sequence {F } N converges in distribution to a d-dimensional standard Gaussian vector N d (0,I d ). 10

11 Proof. As in the one-dimensionalcase, weprovideaproofto the abovetheorem using Malliavin calculus and avoiding the Dambis-Dubins-Schwarz theorem. The equivalences i) ii) iii) iv) follow from Theorem 4. The fact that vi)= v) is easy, and the implication v)= iii) is proved by Peccati and Tudor in [10] using the product formula for multiple stochastic integrals. Hence, it only remains to show the implication iv)= vi). The sequence of random variables {F } N is tight by condition (14). Then, it suffices to show that the limit in distribution of any converging subsequence {F l } l N is N d (0,I d ). For every N, define ϕ (t) = E[e i t,f ]. Let ϕ(t) be the limit of ϕ l (t) as l tends to infinity. As in the proof of Theorem 4 we have for all j = 1,...,d ϕ l ϕ (t) (t). t j l + t j On the other hand, using the definition of the operator L, Proposition 1 and the definition of the operator δ, we have E[F j ei t,f ] = 1 n j E[LF j ei t,f ] = 1 n j E[ δd(f j )ei t,f ] Therefore, = 1 n j E[ DF j,d(ei t,f ) H ] = i n j ϕ l t j (t) = i n j d h=1 d h=1 t h E[e i t,f γ jh ]. t h E[e i t,f l γ jh l ]. (15) Using Lemma 6 and taing the limit of the right-hand side of expression (15) yields ϕ t j (t) = t j ϕ(t), for j = 1,...,d. As a consequence, ϕ is the characteristic function of the law N d (0,I d ). 5 Central limit theorem for square integrable random variables In this section, we will establish a wea convergence result for an arbitrary sequence of centered square integrable random variables. Theorem 8 Let {F } N be a sequence of centered square integrable random variables with the following Wiener chaos expansions Suppose that F = J n (F ). n=1 11

12 (i) lim N + limsup + n=n+1 E[(J nf ) 2 ] = 0, (ii) for every n 1, lim + E[(J n F ) 2 ] = σ 2 n, (iii) n=1 σ2 n = σ2 < +, (iv) for all n 1, D(J n F ) 2 H L 2 (Ω) + nσ2 n. Then, F converges in distribution to the Normal law N(0,σ 2 ) as tends to infinity. Proof. By Theorem 4, conditions (ii) and (iv) imply that for each fixed n 1 the sequence {J n F } N converges in distribution to the normal law N ( 0,σ 2 n). Moreover, by Theorem 7 we have, for each n 1, the following convergence (J 1 F,...,J F ) L + (ξ 1,...,ξ n ), (16) where {ξ n n N are independent centered Gaussian random variables with variances { } σn 2. For every N 1, set n N F N = ξ N = N J n (F ), n=1 N ξ n. n=1 Define also ξ = n=1 ξ n. Let f be a C 1 function such that f and f are bounded by one. Then E[f (F )] E[f (ξ)] E[f (F )] E[f ( F N ) ( ) ] + E[f F N ] E[f(ξ N )] + E[f(ξ N )] E[f (ξ)] ( ) 1/2 E[(J n F ) 2 ] + ( ) E[f F N ] E[f(ξ N )] + E[f(ξ N )] E[f (ξ)]. n=n+1 Taing first the limit as tends to infinity, and then the limit as N tends to infinity, and applying conditions (i), (iii) and 16 we finish the proof. Remar 9 If F D 1,2, and then condition (i) holds. supe( DF 2 H ) <, 12

13 This theorem can be applied to random variables not belonging to D 1,2, and it requires the convergence of the derivatives of the projections. In this sense, it would be interesting to study the relation between the convergence in distribution of a sequence F to a normal law and the convergence in L 2 (Ω) of DF 2 H to a constant. As we have seen, these conditions are equivalent (if F are centered and with lim E[(F 2 )] = σ2 ) for random variables in a fixed chaos. In the general case, we conjecture that this equivalence does not hold. 6 Example Suppose that B H = {Bt H,t 0} is a fractional Brownian motion (fbm) with Hurst parameter H (0,1). That is, B H is a Gaussian stochastic process with zero mean and the covariance function E(B H t BH s ) = 1 2 ( t 2H +s 2H t s 2H). Fix H < 1 2 and an odd integer κ 1. We are interested in the asymptotic behavior of [nt] ( κ, Z (n) t = n κh 1 2 Bj/n (j 1)/n) H BH j=1 as n tends to infinite, where t [0,T]. Set X j = B H j BH j 1. Then, {X j,j 1} is a stationary Gaussian sequence with zero mean, unit variance and correlation ρ H (n) H(2H 1)n 2H 2 as n tends to infinity. We have the following result. Theorem 10 The two-dimensional process (B H,Z (n) ) converges in distribution in the Sorohod space D([0,T]) 2 to (B H,cW), where W is a Brownian motion independent of B H, and c 2 = E[(X 1 X 1+j ) κ ]. j=0 Proof. The proof will be done in two steps. Step 1. We will first show the convergence of finite-dimensional distributions. Let (a,b ], = 1,...,N, be pairwise disjoint intervals conatined in [0,T]. Define the random vectors B = (Bb H 1 Ba H 1,...,Bb H N Ba H N ) and X (n) = (X (n) 1,...,X (n) ), where N X (n) i = n κh 1 2 [na ]<j [nb ] ( B H j/n BH (j 1)/n) κ. We claim that (B,X (n) ) converges in law to (B,V), where B and V are independent and V is a Gaussian random vector with zero mean and independent 13

14 components with variance c 2 (b a ). By the self-similarity of the fbm, it suffices to show the convergence in distribution of (B (n),y (n) ) to (B,V), where B (n) = n H X j, Y (n) = 1 n [na ]<j [nb ] [na ]<j [nb ] with X j = Bj H Bj 1 H and 1 N. We denote by H the Hilbert space defined as the completion of the step functions on [0,T] under the scalar product 1[0,t],1 [0,s] H = E(BH t Bs H ). Using Theorem 8 it suffies to show that: ( ) lim E B (n) n B(n) h = E (( Bb H Ba H )( B H bh Ba H )) h, (17) ( ) lim E B (n) n J my (n) h = 0, (18) ( ) lim E J m Y (n) n J m Y (n) h = δ h (b a )σm, 2 (19) and lim n DJ m Y (n) X κ j, 2 = (b a )mσm 2, (20) H in L 2, for all 1 m κ, and 1,h N. The variances σm 2 must satisfy κ m=1 σ2 m = c2. The projection on the mth Wiener chaos J m Y (n) has the form n c m [na ]<j [nb ] H m(x j ), whereh m (x) denotesthe mth Hermine polynomial. The convergences (17) and (18) are immediate. To prove (19) we write ( ) E J m Y (n) J m Y (n) h = c2 m n [na ]<j [nb ] [na h ]<l [nb h ] δ h c 2 m(b a ) E ( H m (X j )H m (X l ) ) j=0 1 m! ρm H(j), as n tends to infinity. On the other hand, we have DJ m Y (n) 2 = c2 m H n H m 1 (X j )1 (j 1,j] = c2 m n = c2 m n [na ]<j [nb ] [na ]<i,j [nb ] [nb ] i=[na ]+1 2 H H m 1 (X i )H m 1 (X j )ρ H (j i) [nb ] [na ] j=0 H m 1 (X i )H m 1 (X i+j )ρ H (j). 14

15 We claim that the series ξ i = j=0 H m 1(X i )H m 1 (X i+j )ρ H (j), converges almost surely and in L 2, and {ξ i,i 1} is a stationaryergodic [ sequence. The convergence in L 2 follows from the fact that sup j E H m 1 (X i )H m 1 (X i+j ) 2] <, and j=0 ρ H(j) <. On the other hand, the sequence {ξ i,i 1} is ergodic becauuse {X i,i 1} is so. Hence, by the ergodic theorem we have in L 2 lim n DJ m Y (n) 2 = H c2 m(b a ) E(H m 1 (X 1 )H m 1 (X 1+j ))ρ H (j) = c 2 m(b a ) j=0 j=0 1 (m 1)! ρ H(j) m, which implies (20). Step 2. TaingintoaccountthatallL p norms, for1 < p <, areequivalent on a fixed sum of Wiener chaos, in order to show that the sequence Z (n) t is tight in D([0,T]) 2 it suffices to show that E( Z (n) t Z s (n) 2 ) = 1 [nt] n E 2 Xj κ C t s, j=[ns]+1 and this follows easily as above. The convergence of the finite dimensional distributions in the above theorem can also be deduced from general central limit theorems for functionals of Gaussian stationary sequences satisfying the Hermite ran condition (see Brauer and Major [1]). A related result for the function g(x) = x p E( B H 1 p ), where p > 0 and H (0, 3 4 ) was obtained by Corcuera, Nualart and Woerner in [2]. The central limit theorem was proved in this case using the approach of Nualart and Peccati [7] (see [2], Proposition 10). The above result is motivated by the extension of the Itô formula to the fractional Brownian motion in the critical case H = 1 6, using discrete Riemann sums. This problem has been considered by Swanson in [11] in the case of the solution of the one-dimensional stochastic heat equation driven by a space-time white noise{u(t, x), t 0}, which behaves as a fractional Brownian motion with Hurst parameter H = 1 4, for any fiexed x R. In this case, the convergence in law to a Brownian motion is proved for a modified sum of squares of the increments. Acnowledgement 11 David Nualart would lie to than Jason Swanson for very stimulating discussions regarding the Example in Section 6. References [1] P. Breuer and P. Major: Central limit theorems for non-linear functionals of Gaussian fields, J. Multivar. Anal., 13 (1983)

16 [2] J. M. Corcuera, D. Nualart and J. H. C. Woerner: Power variation of some integral fractional processes. Bernoulli 12 (2006) [3] Y. Hu and D. Nualart. Renormalized self-intersection local time for fractional Brownian motion. The Annals of Probability 33 (2005) [4] K. Itô. Multiple Wiener integral. J. Math. Soc. Japan 3 (1951) [5] S. Janson. Gaussian Hilbert Spaces. Cambridge University Press, Cambridge, [6] D. Nualart. The Malliavin Calculus and Related Topics. Springer Verlag, Berlin [7] D. Nualart and G. Peccati. Central limit theorems for sequences of multiple stochastic integrals. The Annals of Probability 33 (2005) [8] G. Peccati and M. Taqqu. Stable convergence of multiple Wiener-Itô integrals. Preprint [9] G. Peccati and M. Taqqu. Stable convergence of generalized stochastic integrals and the principle of conditioning: L 2 theory. Preprint [10] G. Peccati and C. Tudor. Gaussian limits for vector-valued multiple stochastic integrals. Séminaire de Probabilités XXXVIII, pag Lecture Notes in Math., 1857, Springer Verlag, Berlin, [11] J. Swanson: Variations of the solution to a stochastic heat equation. Preprint. 16