Local Times of Gaussian Processes Joachim Lebovits

To cite this version:

Joachim Lebovits. Local Times of Gaussian Processes: with respect to Gaussian Processes Part II. 2017. ￿hal-01489778￿

HAL Id: hal-01489778 https://hal.inria.fr/hal-01489778 Preprint submitted on 14 Mar 2017

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Local Times of Gaussian Processes Stochastic Calculus with respect to Gaussian Processes: Part II

Joachim Lebovits∗† March 14, 2017

Abstract The aim of this work is to define and perform a study of local times of all Gaussian processes that have an integral representation over a real interval (that maybe infinite). Very rich, this class of Gaussian processes, contains Volterra processes (and thus fractional ), multifractional Brownian motions as well as processes, the regularity of which varies along the time. Using the -based anticipative stochastic calculus with respect to Gaussian processes developed in [Leb17], we first establish a Tanaka formula. This allows us to define both weighted and non-weighted local times and finally to provide occupation time formulas for both these local times. A complete comparison of the Tanaka formula as well as the results on Gaussian local times we present here, is made with the ones proposed in [MV05, LN12, SV14]. Keywords: Gaussian processes and local times, stochastic analysis & White noise theory, Wick-Itô integrals, Tanaka formulas, fractional and multifractional Brownian motions. AMS Subject Classification: 60G15; 60H40; 60H05; 60G22 1 Introduction The purpose of this paper is to define and perform a study of local times of all Gaussian processes that have an integral representation over a real interval (that maybe infinite). Theis class of Gaussian G processes, denoted , is defined as being the set of Gaussian processes G := (Gt)t R which can be written under the form: ∈ Gt = gt(u) dBu, (1.1) ZR where R denotes the set of real numbers, R denotes a closed interval of R (that may be equal to

R), B := (Bu)u R is Brownian motion on R and (gt)t R is a family of a measurable square integrable functions on R.∈ The set G obviously contains Volterra processes,∈ as well as Gaussian Fredholm processes1. In order to perform this study of Gaussian local times we will take advantage of the White Noise-based anticipative stochastic calculus with respect to Gaussian processes in G , developed in [Leb17]. Our main results are: A Tanaka formula, that reads, for every T > 0: • T T Gt c = c + sign(Gt c) d⋄Gt + δ c (Gt) dRt, | − | | | Z0 − Z0 { } 2 where the equality holds in L (Ω), where t Rt denotes the variance function of G, which will be supposed to be a continuous function, of bounded7→ variations; the meaning of the different terms will be explained below. Two occupation time formulas that read, for every positive real-valued Borel function Φ: • T (G) Φ(Gs(ω)) ds = ℓT (y)(ω) Φ(y) dy; Z0 ZR · T L (G) Φ(Gs(ω)) dRs = T (y)(ω) Φ(y) dy, Z0 ZR · ∗MAThematic Center of Heidelberg, University of Heidelberg, INF 294, 69120 Heidelberg, Germany. †Laboratoire Analyse, Géométrie et Applications, C.N.R.S. (UMR 7539), Université Paris 13, Sorbonne Paris Cité, 99 avenue Jean-Baptiste Clément 93430, Villetaneuse, France. Email address: [email protected] 1See [Leb17, Introduction] for a precise definition of these processes.

1 (G) L (G) where ℓT (y) (resp. T (y)) denotes the non weighted (resp. the weighted) local time of G (which will be both defined precisely), at point y, up to time T . We will moreover establish that these two latter random variables are L2 random variables. If this work is self-contained, it would be helpful for many readers to refer to [Leb17] to get a precise idea of what White Noise-based stochastic calculus wrt Gaussian processes is. Outline of the paper The remaining of this paper is organized as follows. In Section 2, we recall basic

facts about white noise theory and about the family of operators (MH )H (0,1), which is instrumental for our running example. Section 3 is devoted to the obtention of an Itô formula∈ for generalized functionals of G.This result is a preparatory work for Section 4. More precisely, a Tanaka formula is given in Subsection 4.1, while the definition of weighted and non weighted local times of elements of G , as well as their integral representation is given in Subsection 4.2. An occupation time formula for both, weighted and non weighted, local times of G is established in Subsection 5.1 and constitutes the second main result of this paper. We also prove in Section 5.1 that both, weighted and non weighted local times of G, belong, as two parameters processes, to L2(R Ω, λ µ). Finally, a complete comparison of the Tanaka formula as well as the results on Gaussian× local times⊗ we present here, with the ones proposed in [MV05, LN12, SV14] is made in Subesction 4.1 for Tanaka formula and at the end of Subsection 5.2 for the other results.

2 Background on White Noise Theory and White Noise-based stochastic calculus

Introduced by T. Hida in [Hid75], White Noise Theory is, roughly speaking, the stochastic analogous of deterministic generalized functions (also known as tempered distributions). The idea is to realize nonlinear functional on a Hilbert space as functions of white noise (which is defined as being the time derivative of Brownian motion). We recall in this section the minimum standard set-up for classical white-noise theory. Readers interested in more details may refer to [HKPS93, Kuo96] and [Si12].

2.1 The spaces of stochastic test functions and stochastic distributions

Define N (resp. N∗) the set of non negative integers (resp. positive integers). Let S (R) be the Schwartz space endowed with its usual topology. Denote S ′(R) the space of tempered distributions, which is the dual space of S (R), and F or (F ) the Fourier transform of any element F of S ′(R). For every p F p positive real p, denote L (R) the set of measurable functions f such that R f(u) du < + . When 1 b ixξ | | ∞ f belongs to L (R), f is defined on R by setting f(ξ) := R e− f(x) dx. DefineR the measurable space (Ω, ) by setting Ω := S ′(R) and := (S ′(R)), whereR denotes the σ-algebra of Borel sets. The b b Bochner-MinlosF theorem ensures thatF thereB exists a unique pBrobability measure µ on (Ω, ) such that, for every f in S (R), the map < .,f >: (Ω, ) R defined by < .,f > (ω) =< ω,fF > (where < ω,f > is by definition ω(f), i.e. the actionF of ω→on f) is a centred Gaussian random variable with 2 S variance equal to f L2(R) under µ. The map f <.,f > being an isometry from ( (R),<,>L2(R)) to 2 k k 7→ 2 (L (Ω, , µ),<,>L2(Ω, ,µ)), it may be extended to L (R). One may thus consider the centred Gaussian F F 2 random variable <.,f >, for any f in L (R). In particular, let t be in R, the indicator function 1[0,t] is defined by setting: 1 (s):=1 if 0 s t, 1 (s) := 1 if t s 0 and 1 (s) := 0 otherwise. [0,t] ≤ ≤ [0,t] − ≤ ≤ [0,t] Then the process (Bt)t R, where Bt(ω) := B(t,ω) := < ω, 1[0,t] > is a standard Brownian motion with respect to µ. It then admits∈ a continuous version which will be denoted B. Define, for f in L2(R), e e e I1(f)(ω) := <ω,f >. Then I1(f)(ω)= R f(s) dBs(ω) µ a.s., where R f(s) dBs denotes the Wiener n − 1/2 2 dn 2 R 1/4 n − x /2R x integral of f. For every n in N, let en(x) := ( 1) π− (2 n!) e dxn (e− ) be the n-th Hermite − S function. It is well known (see [Tha93]) that (ek)k N is a family of functions of (R) that forms an orthonormal basis of L2(R, dt). The following properties∈ about the Hermite functions (the proof of which can be found in [Tha93]) will be useful. Theorem 2.1. There exist positive constants C and γ such that, for every k in N,

1/12 γx2 ek(x) C (k + 1)− 1 x 2√k+1 + e− 1 x >2√k+1 . | |≤ · {| |≤ } · {| | } 

2 2 + 2p 2 Let ( ) be the family norms defined by f := ∞ (2k + 2) < f,e > 2 R , for all (p,f) in p p Z p k=0 k L ( ) | | ∈ | | 2 Z L2(R). The operator A, defined on S (R), by settingP A := d + x2 + 1, admits the sequence × − dx2 (en)n N as eigenfunctions and the sequence (2n + 2)n N as eigenvalues. Define, for p in N, the spaces S ∈ 2 S ∈ 2 p(R) := f L (R), f p < + and p(R) as being the completion of L (R) with respect to the norm { ∈. We summarize| | here∞} the minimum− background on White Noise Theory, written e.g. | | p in [LLVH14,− p. 692-693]. More precisely, let (L2) denote the space L2(Ω, , µ), where is the σ- field G G 2 generated by (<.,f >)f L2(R). According to Wiener-Itô’s theorem, for every random variable Φ in (L ) ∈ 2 n there exists a unique sequence (fn)n N of functions in L (R ) such that Φ can be decomposed as + 2 n ∈ 2 n Φ = n=0∞ In(fn), where L (R ) denotes the set of all symmetricb functions f in L (R ) and In(f) denotesP the n th multiple Wiener-Itô integral of f with the convention that I0(f0) = f0 for constants − + b + n 2 f0. For any Φ := n=0∞ In(fn) satisfying the condition n=0∞ n! A⊗ fn 0 < + , define the element 2 + n n | | ∞ Γ(A)(Φ) of (L ) byP Γ(A)(Φ) := ∞ I (A⊗ f ), wherePA⊗ denotes the n th tensor power of the n=0 n n − operator A (see [Jan97, AppendixP E] for more details about tensor products of operators). The operator 1 Γ(A) is densely defined on (L2). It is invertible and its inverse Γ(A)− is bounded. We note, for ϕ in 2 2 2 n (L ), ϕ := ϕ 2 . For n in N, let Dom(Γ(A) ) be the domain of the n-th iteration of Γ(A). Define k k0 k k(L ) the family of norms ( ) by: p p Z kk ∈ p p 2 p Φ := Γ(A) Φ = Γ(A) (Φ) 2 , p Z, Φ (L ) Dom(Γ(A) ). k kp k k0 k k(L ) ∀ ∈ ∀ ∈ ∩ 2 p 2 For p in N, define ( p) := Φ (L ): Γ(A) (Φ) exists and belongs to (L ) and define ( p) as the S 2{ ∈ } S− completion of the space (L ) with respect to the norm p. As in [Kuo96], we let ( ) denote the kk− S projective limit of the sequence (( p))p N and ( )∗ the inductive limit of the sequence (( p))p N. S ∈ S S− ∈ This means in particular that ( ) (L)2 ( )∗ and that ( )∗ is the dual space of ( ). Moreover, S ⊂ ⊂ S S S ( ) is called the space of stochastic test functions while ( )∗ the Hida distribution space. We will note S S <<, >> the duality bracket between ( )∗ and ( ). If φ, Φ belong to (L2), then we have the equality S S S S << Φ, ϕ >>= < Φ,ϕ>(L2) = E[Φ ϕ]. Besides, denote <,> the duality bracket between ′(R) and (R) + and recall that every tempered distribution F can be written as F = n=0∞ < F,en >en, where the convergence holds in S ′(R). The next proposition, that will be usedP extensively in the sequel, is a consequence of the definition of ( ) and ( )∗. S S S + Proposition 2.2. Let F be in in ′(R). Define <.,F >:= n=0∞ < F,en ><.,en >. Then there exists 2 p0 in N such that that <.,F > belongs to ( p ), and henceP to ( )∗. Moreover we have <.,F > = − 0 p0 2 + S S N k k− F p . Conversely, define Φ := n=0∞ bn < .,en >, where (bn)n N belongs to R . Then Φ belongs to | |− 0 ∈ P + 2 2p0 ( )∗ if and only if there exists an integer p0 in N such that n=0∞ bn (2n + 2)− < + . In this S + ∞ ∞ S S latter case F := n=0 bnen belongs to p0 (R) and then to P′(R). It moreover verifies the equality 2 + 2p 2 − F = ∞ bP2 (2n + 2)− 0 = Φ . | | p0 n=0 n k k p0 − P − 2.2 ( )∗-process, ( )∗-derivative and ( )∗-integral S S S Let (R, (R),m) be a sigma-finite measure space. Through this section, I denotes an element of (R). B B A measurable function Φ : I ( )∗ is called a stochastic distribution process, or an ( )∗-process. An → S 1 S ( )∗-process Φ is said to be differentiable at t0 I if lim r− (Φt +r Φt ) exists in ( )∗. We note S ∈ r 0 0 − 0 S dΦt → 0 the ( )∗- derivative at t of the stochastic distribution process Φ. It is also possible to define an dt S 0 ( )∗-valued integral in the following way (one may refer to [Kuo96, p.245-246] or [HP57, Def. 3.7.1 p.77] forS more details).

Theorem-Definition 2.1 (integral in ( )∗). Assume that Φ : I ( )∗ is weakly in L1(I,m), i.e. assume that for all ϕ in ( ), the mapping uS << Φ , ϕ>>, from I to →R, belongsS to L1(I,m). Then there S 7→ u exists an unique element in ( )∗, noted Φ m(du), such that, for all ϕ in ( ), S I u S R << Φ(u) m(du), ϕ >> = << Φu,ϕ >>m(du). ZI ZI

We say in this case that Φ is ( )∗-integrable on I (with respect to the measure m), in the Pettis sense. In the sequel, unless otherwiseS specified, we will always refer to the integral in Pettis’ sense

3 2.3 S-transform and Wick product For f in L2(R), define the Wick exponential of <.,f >, noted : e<.,f> :, as the (L2) random variable <.,f> 1 f 2 equal to e − 2 | |0 . The S-transform of an element Φ of ( ∗), noted S(Φ), is defined as the function <.,η> S from S (R) to R given by S(Φ)(η) := <<Φ, : e :>> for any η in S (R). For any (Φ, Ψ) ( )∗ ( )∗, ∈ S × S there exists a unique element of ( )∗, called the Wick product of Φ and Ψ, and noted Φ Ψ, such that S(Φ Ψ)(η) = S(Φ)(η) S(Ψ)(η)S for every η in S (R). Note that, when Φ belongs to (⋄L2), SΦ(η) is ⋄ <.,η> 1 η 2 <.,η> nothing but E[Φ : e :] = e− 2 | |0 E[Φ e ]. The following result will be intensively used in the sequel. 2 η max{p;q} Lemma 2.3. [LLV14, Lemma 2.3.] The inequality S(X Y )(η) X p Y q e| | is valid for 2 | ⋄ |≤k k− k k− any (p, q) in N and (X, Y ) in ( p) ( q). S− × S− Some useful properties of S transforms are listed in the proposition below. The proof of the results stated in this proposition can be found in [Kuo96, Chap 5]. Proposition 2.4 (Some properties of S transforms). When Φ is deterministic then Φ Ψ=ΦΨ, for all + + ⋄ Ψ in ( )∗. Moreover, let Φ= k=0∞ ak<., ek > and Ψ= n=0∞ In(fn) be in ( )∗.Then their S-transform S S + S + n is given, for every η in (R),P by S(Φ)(η)= k=0∞ ak < η,eP k >L2(R) and S(Ψ)(η)= k=0∞ . Finally, for every (f,η,ξ) in L2(R) S (R)P R, we have the equality: P × × iξ<.,f> 1 ( η 2+2iξ ξ2 f 2) S(e )(η)= e 2 | |0 − | |0 . (2.1)

One may refer to [Jan97, Chap.3 and 16] for more details about Wick product. The following results on the S-transform will be used extensively in the sequel. See [Kuo96, p.39] and [HKPS93, p.280-281] for proofs. Denote (A; B) the set of B-valued functions defined on A. F Lemma 2.5. The S-transform verifies the following properties:

(i) The map S :Φ S(Φ), from ( )∗into (S (R); R), is injective. 7→ S F (ii) Let Φ : I ( )∗ be an ( )∗process. If Φ is ( )∗-integrable over I wrt m, then one has, for all η S → S S S in (R), S( I Φ(u) m(du))(η)= I S(Φ(u))(η) m(du). R R (iii) Let Φ : I ( )∗ be an ( )∗-process differentiable at t I. Then, for every η in S (R) the map u [SΦ(u→)](ηS) is differentiableS at t and verifies S[ dΦ (t)](∈ η)= d S[Φ(t)](η) . 7→ dt dt   The next theorems provide a criterion for integrability in ( )∗, in term of S-transform. S Theorem 2.6. [Kuo96, Theorem 13.5] Let Φ : I ( )∗ be a stochastic distribution such that, for all η in S (R), the real- valued map t S[Φ(t)](η) is→ measurableS and such that there exist a natural integer 7→1 a η 2 p, a real a and a function L in L (I,m) such that S(Φ(t))(η) L(t) e | |p , for all η of S (R) and for | |≤ almost every t of I. Then Φ is ( )∗- integrable over I, wrt to m. S We end this section with the following theorems that will be useful in the next section.

Theorem 2.7. [Ben03, Theorem 2.17] For any differentiable map F : I S ′(R), the element <., F (t)> is a differentiable stochastic distribution process which satisfies the equality:→ d dF < ., F (t) > = < ., (t) > dt dt Theorem 2.8. [CG62, Lemmas 1 and 2 p.73-74] Let I be an interval of R, t F (t) be a map from I 7→ into S p(R), t ϕ(t, .) be a map from I into S (R) and t0 I. If both maps t F (t) and t ϕ(t, .) are continuous− (respectively7→ differentiable) at t , then the function∈ t < F (t), ϕ7→(t, .) > is continuous7→ 0 7→ (respectively differentiable) at t0. In this latter case we have

d < F (t), ϕ(t, .) > dF d (t )=< (t ), ϕ(t, .) > + < F (t), [t ϕ(t, .)](t ) > . dt 0 dt 0 dt 7→ 0

4 2.4 Stochastic integral with respect to R Denote G := (Gt)t R the process defined, for every t in , by Gt :=< .,gt >, where (gt)t R is a ∈ 2 ∈ family of functions of L (R). Denote (t,s) Rt,s the covariance function of G. We hence have Rt,s := 7→ R2 E[Gt Gs]= L2(R), for every (s,t) in . We will note in the sequel Rt instead of Rt,t. G is a a.s. Gaussian process which fulfills For the sake of notational simplicity we can and will assume that G0 = 0. Moreover, when the Gaussian process G will admit a continuous modification, we will systematically use it and still call it G.

White Noise derivative of G R S Define the map g : ′(R) by setting g(t) := gt. When g is differentiable at point t, one denotes gt′ its derivative. Denote λ the→ Lebesgue measure on R and define L1 (R) := f : R R is measurable ; f loc { → ∈ L1((a,b)), for all finite interval (a,b) s.t. [a,b] R . In this section and in the next one, we will make the following assumption: ⊂ } e (A ) a = b (i) The map g is continuous on R, A a = b (i) (ii) The map g is differentiable λ almost everywhere on R, A  − A(ii)  a = b 1 R ( iiiiii)) (iii) There exists q in N∗ such that t gt′ q belongs to L ( ),  7→ | |− loc A(iii)  2 a = b (iv) For every (a,b) in R such that a b, one has, in S ′(R), the equality: A a = b≤ (iv) pfojerpfoerofjerof b e   gb ga = gu′ du. (Ea,b)  − Za Define the set R by setting R := t R; g is differentiable at point t . For the sake of notational D D { ∈ } simplicity we will write A(i), (resp. A(ii), A(iii) or A(iv)), in the sequel, when one wants to refer to (i) (resp. to (ii), (iii) or (iv)) of Assumption2 (A ). Remark 1. A first consequence of Assumption (A ) is that g is “weakly” locally absolutely continuous on R; this means that the map t is absolutely continuous on every finite interval [a,b] of R, for every η in S (R). 7→

Theorem 2.9 (Gaussian White Noise [Leb17, Theorem-Definition 3.1]). Define W (G) := (W (G)) t t RD (G) (G) ∈ by setting: W t := < .,gt′ >, where the equality holds in ( )∗. Then W is is the ( )∗-derivative of S (G) S (G) the process (G ) . We will sometimes note dGt instead of W . Moreover the map t W t t RD dt t t p ∈ 7→ k k− is continuous if and only if t gt′ p is continuous. 7→ | |− Generalized functionals of G In order to get a Tanaka formula in Section 4, we define here generalized functionals of G, by using [Kuo96, 1 Section 7.1]. We identify, here and in the sequel, any function f of Lloc(R) with its associated tempered distribution, denoted Tf , when it exists. In particular, one notes in this case: = R f(t) φ(t) dt, for every φ in S (R). In this latter case we say that the tempered distribution T := TfR is of function type. Define the sets := t R; R =0 and c := t R; R > 0 . ZR { ∈ t } ZR { ∈ t } Theorem-Definition 2.2. Let F be a tempered distribution. For every t in c , define ZR + ∞ 1 1 k F (Gt) := < F,ξt,k > Ik gt⊗ , √2πR k! Rk t kX=0 t 

1/4 1/2 k/2 x2 where, for every (x, k) in R N, ξt,k(x) := π (k!) R exp ek(x/(√2Rt)). Then for all real × t {− 4Rt } t, F (Gt) is a Hida distribution, called generalized functional of Gt.

As stated in [Ben03], when F = f is of function type, F (Gt) coincides with f(Gt).

2 An easy way to see if Assumption (A ) holds is to check if whether the following condition holds. (i) The map g is continuous on R and differentiable on every finite interval (a, b) s.t. [a, b] ⊂ R, D ( ) N∗ ′ 1 R (ii) There exists q in s.t. t 7→ |gt| q belongs to L ( ).  − loc

e 5 S-Transform of G and W (G) The following theorem makes explicit the S-transforms of G, of the Gaussian white noise W (G) and of generalized functionals of G. Denote γ the heat kernel density on R R i.e. + × 1 x2 γ(t, x) := exp − if t = 0 and 0 if t =0. (2.3) √2πt { 2t } 6 Theorem 2.10. [Leb17, Theorem 3.6] For every η in S (R) one has the following equalities: R (i) S(Gt)(η)= L2(R), for every t in , (G) d R (ii) S(Wt )(η)= = dt [L2(R)], for every t in D; S c (iii) For p N, F p(R), and t in R, S(F (Gt))(η)= F,γ (Rt, . ) . ∈ ∈ − Z h − i Furthermore, there exists a constant Dp, independent of F,t and η, such that:

c 2 2p 2p 1/2 2 2 t R, S(F (Gt))(η) max Rt− ; Rt Rt− Dp F p exp η p . ∀ ∈ Z | | ≤ { } | |− {| | } The Wick-Itô integral with respect to Gaussian processes We still assume in this section that Assumption (A ) holds and still denote I a Borel set of R.

Definition 1 (Wick-Itô integral wrt Gaussian process). Let X : R ( )∗ be a process such that the (G) R → S R process t Xt Wt is ( )∗-integrable on . The process X is then said to be dG-integrable on (or integrable7→ on R⋄), wrt theS Gaussian process G. The Wick-Itô integral of X wrt G, on R, is defined by setting:

(G) Xs d⋄Gs := Xs Ws ds. (2.4) ZR ZR ⋄

For any I in (R), define X d⋄G := 1 (s) X d⋄G . B I s s R I s s R R The Wick-Itô integral of an ( ∗)-valued process, wrt G is then an element of ( )∗. It is easy to see that Wick-Itô integration wrt G, isS linear and, assuming it belongs to (L2), centered.S c S Theorem 2.11. [Leb17, Theorem 3.13] Let p N, [a,b] be an interval of R and let F : [a,b] p(R) − p 1/4 p∈ 1/4 1 Z → be a continuous map. If t max Rt− − ; Rt − belongs to L ([a,b]) (resp. there exists an integer q 7→ p 1/4{ p 1/4 } 1 such that t gt′ q max Rt− − ; Rt − belongs to L ([a,b])), then the stochastic distribution process 7→ | |− { } F (t, G ) is ( )∗-integrable (resp. dG-integrable) on [a,b] (resp. on (a,b)). t S 2.5 Gaussian Processes in G of “reference” To see the generality of the results on local times we present here, we will consider in this paper, classical Gaussian processes, made with elements of G . These processes are: Brownian motion and , fractional and multifractional Brownian motions as well as Vγ - processes.

Fractional, Multifractional Brownian motions and Vγ - processes Readers interested in an exhaustive presentation of fBm or mBm may refer to [Nua06] for fBm and to [LLVH14] for mBm, as well as to the references therein. Recall that fBm, which was Introduced in [Kol40] and popularized in [MVN68], is a centered Gaussian process, the covariance function of which is denoted RH and is given by: 1 R (t,s) := ( t 2H + s 2H t s 2H ), H 2 | | | | − | − | where H belongs to (0, 1), and is usually called the Hurst exponent. When H = 1/2, fBm reduces to standard Brownian motion. MBm, which is a Gaussian extension of fBm, was introduced in [PLV95] and in [BJR97] in order to match any prescribed non-constant deterministic local regularity and to decouple this property from long range dependence. A mBm on R, with functional parameter h : R h h h → (0, 1), is a Gaussian process B := (Bt )t R defined, for all real t, by Bt := B(t,h(t)), where B := ∈ (B(t, H))(t,H) R (0,1) is fractional Brownian field on R (0, 1) (which means that B is a Gaussian field, ∈ × × such that, for every H in (0, 1), the process (B(t, H))t R is a fBm with Hurst parameter H). In other words, a mBm is simply a “path” traced on a fractional Brownia∈ n field. Note also that when h is constant,

6 mBm reduces to fBm. For any deterministic function γ : R+ R, Vγ - processes are defined in [MV05, γ γ → Prop. 1] as being the processes, denoted B := (Bt )t [0,T ], by setting: t ∈ γ e e Bt := ε(t s) dWs; t [0,T ], (2.5) Z0 − ∀ ∈ 2 e 2 2 with γ : R+ R such that γ is of class C everywhere in R+ except in 0; and such that (γ )′ is non → 2 increasing. The map ε : R∗ R is defined by setting ε := (γ ) . Subset of G , the set V contains + → ′ γ Gaussian processes, that can be more irregular than any fBm.p However it does not contain fBm (nor mBm) since Vγ only contains processes the regularity of which remains constant along the time).

Operators (MH )H (0,1) ∈ The operator MH will be useful in the sequel, not only to provide one with a representation of fBm and of mBm under the form (1.1), but also to verify that Assumption (A ) made page 5, hold for both fBm and mBm. Let H belongs to (0, 1); following [EVdH03] and [LLV14, Section 2.2], define the 2 √2π 1/2 H L (R)-valued operator M , in the Fourier domain by: M (u)(y) := y − u(y), y R∗, H H cH 1 | | ∀ ∈ 2π 2 where cx is defined, for every x in (0, 1) by cx := Γ(2x+1) sin(dπx) . This operator is wellb defined on the 2 S 1 homogeneous Sobolev space LH (R) := u ′(R) : u = Tf ;f Lloc(R) and u H < + , where { ∈ ∈ k2 k ∞} the norm H derives from the inner product denoted , H , which is defined on LH (R) by: u, v H := 1 1k·k2H h·b ·i 2 2 h i c2 R ξ − u (ξ) v (ξ) dξ. MH being an isometry from (LH (R), H ) into (L (R), L2(R)), it is H | | k·k k·k clearR that, for every (H,t,s) in (0, 1) R2, = R (t,s). We will say that b b × H [0,t] H [0,s] L2(R) H an mBm is normalized when its covariance function, denoted Rh, verifies the equality:

c2 ht,s 1 2ht,s 2ht,s 2ht,s Rh(t,s)= t + s t s , (2.6) c(h(t)) c(h(s)) 2 | | | | − | − |   h(t)+h(s) where ht,s := 2 and cx has been above, right after MH (u)(y). Example 2.12 (Gaussian Processes in G of “reference”). LetdH be real in (0, 1) and h : R (0, 1) be a deterministic measurable function. Define the processes → B := < ., 1 >; t R ; B := < ., 1 t 1 >; t [0, 1] ; { [0,t] ∈ } { [0,t] − · [0,1] ∈ } H h B := < .,MH (1[0,t]) >; t R ; Bb := < .,Mh(t)(1[0,t]) >; t R ; { ∈ } { ∈ } γ γ B := < ., 1 ε(t .)) > if t R∗ & B := 0 . { [0,t) · − ∈ + 0 } Wee know, thanks to Section 2.1, that B is a Browniane motion on R, that B is Brownian bridge on [0, 1], that BH is a fBm of Hurst index H, that Bh is a normalized mBm of functional parameter h and that γ b B is a Vγ - process. H h(t) h Ae word on notation: B. or B. will always denote an fBm with Hurst index H or h(t), while B. will stand for an mBm.

3 Itô Formulas for generalized functionals of G

The main result of this section is an Itô Formula, in (L2)-sense, for generalized functionals of G (Theorem 3.4). This latter will not only provide us with a Tanaka formula, but will also give us a precise reason to define the weighted local time of G, the way we do it, in Section 4. To do so we first need to establish an Itô formula in ( )∗. This is done in Subsection 3.1. Let [a,b] be an intervalS of R and j : [a,b] R be a function of bounded variation. Denote α the → j signed measure such that j(t) = αj ([a,t]), for every t in [a,b]. For any function f : [a,b] R, denote b b → a f(s) dj(s) or a f(s) dαj (s) the Lebesgue-Stieljes integral of f with respect to j, assuming it exists. In 1 1 thisR latter case, weR will write that f L (I, dj(t)) or L (I, αj ). In the particular case where the function f is continuous on [a,b], the Lebesgue-Stieljes∈ integral of f exists and is also equal to the Riemann-Stieljes integral of f, which is denoted and defined by:

b n (n) (R.S.) f(s) dj(s) := lim f(ξi ) (j(xi) j(xi 1)), (3.1) Z π 0 − − a → Xi=1

7 (n) where the convergence holds uniformly on all finite partitions Pπ := a := x0 x1 xn := b (n) { ≤ ≤···≤ } of [a,b] such that max (xi xi 1) π and such that ξi belongs to [xi 1, xi]. The following result, will 1 i n − − ≤ − be used extensively≤ in≤ the sequel of this section.

Lemma 3.1. [Leb17, Lemma 4.1] Let [a,b] be a finite interval of R, I (resp. J) an interval of R+ (resp. of R) and let L : [a,b] I J be a C1-function. Let f : [a,b] I and j : [a,b] J be two continuous functions of bounded variation× × on [a,b]. Then one has the following→ equality: → b ∂L b ∂L L(b,f(b), j(b)) L(a,f(a), j(a)) = (s,f(s), j(s)) ds + (s,f(s), j(s)) df(s) − Za ∂u1 Za ∂u2 b ∂L + (s,f(s), j(s)) dj(s). (3.2) Za ∂u3 In the remaining of this paper, and unless otherwise specify, the measure m denotes a measure, that may be σ-finite or signed. Lemma 3.2. Let T > 0 and v : [0,T ] R R be a continuous function such that there exists a couple 2 × → λT y (CT , λT ) of R R+∗ such that max v(t,y) CT e for all real y. Define; for every a > λT , the × t [0,T ]| | ≤ ∈ map J : R (0, 1/4a) R R by setting: v + × × →

Jv(t,u1,u2) := v(t, x) γ(u1, x u2) dx. (3.3) ZR · −

Then Jv is well defined. Moreover lim Jv(t,u1,u2)= v(t0,l0), (t0,l0) in [0,T ] R. + (t,u1 ,u2) (t0,0 ,l0) ∀ × → 3.1 Itô Formula in ( )∗ for generalized functionals of G on an interval of c S ZR Before establishing an Itô formula in ( )∗, let us first fix some notations. For a tempered distribution G (n) S th (1) and a positive integer n, let G denote the n distributional derivative of G. We also write G′ := G . Hence, by definition, the equality < G′,ϕ>= < G, ϕ′ > holds for all ϕ in S (R). For a map t F (t) S ∂nF − (n) th S7→′ from [a,b] to p(R) we will note ∂xn (t) the quantity (F (t)) , that is the n derivative in (R), − ∂nF S of the tempered distribution F (t). Hence we may consider the map t n (t) from [a,b] to ′(R). 7→ ∂x ∂F F (t0+r) F (t0) Moreover for any t0 in [a,b], we will note (t0) the quantity lim − , when it exists in S p(R), ∂t r 0 r − ∂F → for a certain integer p. When it exists, ∂t (t0) is a tempered distribution, which is said to be the derivative of the distribution F (t), with respect to t at point t = t0. In line with Subsection 2.4, we then define, ∂nF (n) for t0 in [a,b] and a positive integer n, the following quantities: ∂xn (t0, Gt0 ) := (F (t0)) (Gt0 ) and ∂F ∂F ∂t (t0, Gt0 ) := ∂t (t0) (Gt0 ).  c 1 S Theorem 3.3. Let [a,b] be an interval of R, p N, and let F be an element of C ([a,b], p(R)) − ∂F ∂2F Z S∈ A such that both maps ∂x and ∂x2 , from [a,b] into p(R), are continuous. Assume moreover that ( ) holds and that the map t R is both continuous and− of bounded variations on [a,b]. Then the following 7→ t equality holds in ( )∗: S b b b ∂F ∂F 1 ∂2F F (b, Gb) F (a, Ga)= ∂t (s, Gs) ds + ∂x (s, Gs) d⋄Gs + 2 ∂x2 (s, Gs) dRs. (3.4) − Za Za Za

Proof. Assumption A(iii) and the fact the map t Rt is upper and lower bounded on [a,b] allow one ∂F ∂2F 7→ to apply Theorem 2.11 to ∂t (s, Gs), ∂x2 (s, Gs) (with the Borel measure defined by t Rt), and to ∂F 7→ ∂x (s, Gs). This and [Leb17, Remark 6] entail that all the integrals on the right side of (3.4) exist. Thanks to Lemma 2.5 it is then sufficient to show the equality of the S-transforms of both sides of (3.4). Using Theorem 2.10 one then just has to establish, for every η in S (R), the following equality: b b ∂F ∂F Ia,b(η)= < ∂t (s),γ(Rs, . ) > ds + < ∂x (s),γ(Rs, . ) >< gs′ ,η >ds Za − Za − b ∂2F + < ∂x2 (s),γ(Rs, . ) > dRs, (3.5) Za −

8 where we set I (η) :=< F (b),γ(R , . ) > < F (a),γ(R , . ) >. Define the map a,b b − b − a − a L : [a,b] R+ R by setting: L(u1,u2,u3) :=< F (u1),γ(u2, . u3) >. The assumptions made on F × × 1 ∂γ −1 ∂2γ entail that L is a C -function. Moreover, using the equality ∂t = 2 ∂x2 , one gets:

∂F ∂γ < ∂x (s),γ(Rs, . ) > = < F (s), ∂x (Rs, . ) >, 2 − − < ∂ F (s),γ(R , . ) > = 1 < F (s), ∂γ (R , . ) > . ∂x2 s − s 2 ∂t s − s Since I (η)= L(b, R ,γ()) L(a, R ,γ()), (3.5) also reads as (3.2) with f(s) := R a,b b b − a a s and j(s) := jη(s) :=. According to Remark 1, jη is absolutely continuous on [0,T ]. Lemma 3.1 therefore applies and achieves the proof.

Remark 2. As Theorem 3.4 below will show, one can extend (3.4) to the case where [a,b] R = , if one makes additional assumptions on R and F . However, we will not extend (3.4) directly∩ since Z 6 one∅ does not need such a general formula to get Theorem 3.4, and hence a Tanaka formula.

3.2 Itô Formula in (L2) for certain generalized functionals of G on any interval of R To get a Tanaka formula, one needs both: to generalize (3.4) to the case where [a,b] contains elements of , and to specify (3.4) to the case where F (t, x) = F (x). However, Theorem 3.3 does not extend ZR immediately to the case where [a,b] intersects the set R because the generalized functionals of G are not defined in this situation. As in [Ben03], we now extendZ our previous results to deal with this difficulty. Having in mind the definition of distributions of function type, given right above Theorem-Definition 2.2, one has the following result.

Theorem 3.4. Let [a,b] be an interval of R, F be a distribution of function type that is continuous on a,b := [a,b] . Assume that the first distributional derivative of F , noted F ′, is also of function type ZR ∩ ZR and that F and F ′ both have polynomial growth. Assume moreover that (A ) holds, that the map t Rt is both continuous and of bounded variations on [a,b] and that the following assumptions are satisfied:7→ (i) λ( a,b)= α ( a,b)=0, ZR R ZR (ii) t R−1/2 belongs to L1([a,b], dR ), 7→ t t 1/2 1 (iii) q N such that t gt′ q Rt− L ([a,b], dt), ∃ ∈ 7→ | |− ∈

where λ denote the Lebesgue measure on R and where the signed measure αR has been defined at the beginning of Section 3. Then, the following equality holds in (L2): b 1 b F (Gb) F (Ga)= F ′(Gt) d⋄Gt + F ′′(Gt) dRt. (3.6) − Za 2 Za

(G) Proof. Let us first show that t F ′(Gt) Wt is ( )∗- integrable on [a,b] and that t F ′′(Gt) is 7→ ⋄ a,bS 7→ ( )∗-integrable on [a,b] with respect to α . Denote Γ := t [a,b]; R > 0 . S R R { ∈ t } (G) a) t F ′(G ) W is ( )∗-integrable on (a,b). 7→ t ⋄ t S 2 Denote Ra,b := sup Rt and D := sup ek(u) : k N,u R . We have the following upper bound, t [a,b] { | | ∈ ∈ } ∈ a,b valid for all k in N and t in ΓRb , 2 2 ′ 2 1/2 k d 1/2 2 < F ′, ξ > = < F,ξ > = π k! R < F, [e (x/√2R )exp x /4R ] > , | t,k | | t,k | t | dx k t {− t} | ′ d where ξt,k denotes dx [x ξt,k(x)] (which has been defined right after Theorem-Definition 2.2). The right hand side of the previous7→ equality can be decomposed into two parts. Since we get the upper bound of these two parts similarly, we only show here how to get the upper bond of 1/2 1/2 2 < F, (2 R )− e′ (x/√2R )exp x /4R > =: J . | t k t {− t} | t,k

9 Denote K (resp. N) the real (resp. the integer) such that F (x) K(1 + x N ), for every real x. Using | |≤ | | k k+1 the relatione ek′ (x)= ek 1(x) ek+1(x) (see [Kuo96, p.354]),e valid for all positive integer k, and 2 − − 2 the polynomial growthq of F , one immediatelyq gets:

2 1/2 x /4Rt Jt,k (2 Rt)− F (x) 2(k + 1) D e− dx ≤ ZR | | p N b N u2/4 (1 + R) DK (1 + u ) e− du √k +1. ≤ p ZR | | b e It is now easy to obtain the next inquality, for a real D′ which does not depend on k nor t:

a,b 2 k (t, k) Γ N, < F ′, ξ > D′(k + 1)! R . (3.7) ∀ ∈ R × | t,k | ≤ t 1 Moreover, since the norm operator of A− is equal to 1/2 (see [Kuo96, p.17]), we get for a real D that a,b does not depend on t, using Theorem-Definition 2.2 and (3.7), for all t in ΓR :

+∞ 2 1 1 2k 2 1 2k 1 − F ′(Gt) 1 = 2πR (k!) Rt− < F ′, ξt,k > A− gt 0 D Rt− . (3.8) k k− t | | ≤ · kP=0 In view of Lemma 2.3 and of (3.8), it is clear that there exists a real K, such that, for every η in S (R) (G) (G) η 2 1/2 η 2 a,b R q − q and t in ΓR D, S(F ′(Gt) Wt )(η) F ′(Gt) 1 Wt q e| | K gt′ q Rt e| | . One can ∩ | ⋄ |≤kb k− k k− ≤ | |− then apply Thm. 2.6 and conclude that F ′(G ) d⋄G ( )∗. a t t ∈ S R b) t F ′′(G ) is ( )∗-integrable on [a,b], with respect to α . 7→ t S R Since F ′ also has polynomial growth we get, by replacing F by F ′ in (3.8), using Lemma 2.3, S(F ′′(Gt))(η) 2 1/2 η a,b | |≤ − 1 K Rt e| | , for every t in ΓR , Thanks to Assumption (i), we know that the right hand side of the previ- 1 a,b b ous inequality belongs to L ([a,b] , dR ). Theorem 2.6 then applies and ensures that F ′′(G ) dR \ZR t a s s exists in ( )∗. To establish (3.6), it is then sufficient to show that the following equalityR holds in (L2): S 1 F (Gb) F (Ga)= F ′(Gt) d⋄Gt + F ′′(Gt) dRt. (3.9) − Z a,b 2 Z a,b ΓR ΓR

a,b Having in mind (iii) of Theorem 2.10, we proceed as in the proof of [Leb17, Theorem 4.4]. Since ΓR is an open set of [a,b], it can be written under the form: Γa,b = [a,c) (c , d ) (d, b], (3.10) R ⊔ i i ⊔ iGN ∈ where all the intervals in (3.10) are disjoint and where, by convention, (x, y) = (x, y] = , for every ∅ reals x and y such that x y. Note moreover that every element of ci, di, i N (as well as c, if ≥ a,b { ∈ } [a,c) = , and d, if (d, b] = ) belongs to R . Exactly as in the proof of [Leb17, Theorem 4.4], one needs to distinguish6 ∅ between two6 ∅ cases: Z

2 a,b First case: (c′, d′) [a,b] with c′ < d′ s.t. [a,c′) and (d′,b] are both subsets of Γ . ∃ ∈ R 2 a,b Second case: ∄ (c′, d′) [a,b] with c′ < d′ s.t. [a,c′) and (d′,b] are both subsets of Γ . ∈ R In both cases, the same reasoning as the one we used in the first case of the proof of [Leb17, Theorem 4.4] (with Theorem 3.3 instead of Lemma 3.1 and Lemma 3.2 with v(t, .)= F ) applies and allows us to conclude. Remark 3. As it has been noticed in [Ben03, Remark 4.6] (in the case of fBm), Equality (3.6) holds 2 2 2 in (L ) in the sense that both sides belong to (L ). If F ′′(Gs) belongs to (L ) and under the additional v v 2 assumption that F ′′(G ) dR < + , then F ′′(G ) dR is an (L )-valued Bochner integral (see u k s k0 s ∞ u s s Appendix A.1 forR the definition of Bochner integral)R and all members of both sides of (3.6) belong to (L2). Note moreover that the assumption t R is bounded on [a,b] is sufficient but not necessary; however, if 7→ t one wants to relieve assumption on R, one will have to weigh down the ones on F , F ′ or on the process G. Once again, all the Gaussian processes of “reference” fulfill assumptions of Theorem 3.4. The only non obvious case is the one of Vγ - processes. In this latter case, t Rt is absolutely continuous on [0,b] and therefore condition (iii) of Theorem 3.4 is fulfilled. Besides, because7→ R is: absolutely continuous on [0,b],

10 non decreasing and such that limu 0+ Ru′ =+ , it is easy to show, using [Leb17, (B.6)], that condition (i) of Theorem 3.4 entails condition→ (ii). Finally,∞ one has the equality:

(2e)−1 (2e)−1 (2e)−1 −1/2 −1/2 1/2 1/2 I := Ru dRu = Ru Ru′ du =2 (Ru )′ du =2 R1/2e. (3.11) Z0 Z0 · Z0 · 4 Tanaka formula & Local times

The main results of this section are a Tanaka formula (Theorem 4.1) and two occupation time formulas (Theorem 5.1). In Subsection 4.1 we perform a complete comparison between our Tanaka formula and the ones provided so far, for general Gaussian processes (namely [LN12, Theorem 4.3] and, to a lesser extent, in [MV05, Theorem 36]). In particular, it will be shown that the Tanaka formula we provide here not only fully extends the one provided in [MV05, Theorem 36] but is also less restrictive, in many places, than the one provided by [LN12, Theorem 4.3]. We then apply our Tanaka formula to the Gaussian processes of “reference” and compare the resulting Tanaka formulas hence obtained with the Tanaka formulas established so far, for these particular Gaussian processes. In Subsection 4.2 we provide a definition and a basic study of both (non weighted) and weighted local times of elements G in G . After having provided an integral representation of each of these two processes, we establish, under classical assumptions, the belonging of both these local times, as two parameters processes, to L2(R Ω, λ µ). The proof, that (non weighted) local time of G belongs to L2(R Ω, λ µ) we give here,× is not⊗ only new and based on White noise analysis; it also offers the advantage× to be⊗ easily adaptable to the weighted local time of G. It is then shown in what extent the results we provide in this Subsection are new, and we also give a complete comparison of our results to the one obtained in [LN12], for weighted local times (that are the only local times considered in this latter reference). Finally, we apply our result to the Gaussian processes of “reference” and show how we extend the results known for them, so far. In addition to the references [GH80, MR06] and [RY99, Chapter VI], that provide a complete overview on local times, a summary on fractional Brownian motion’s local time (that includes the concept of weighted fractional local time) can be found in [Cou07, Section 7.8] (see also references therein, in particular [CNT01]). Note that a study of fBm’s local time in the white noise theory’s framework can be found in [HØS05]. The non-weighted local times of particular mBms have been studied in [BDG06, BDG07] in the case of an harmonizable mBm while the weighted local time of a Volterra-type mBm has been studied in [BDM10]. This approach has been extended to multifractional Brownian sheets in [MWX08]. Through this section, let T denotes a positive real and recall that T := [0,T ]. ZR ZR ∩ 4.1 Tanaka Formula A consequence of Theorem 3.4 is the following Tanaka formula. Theorem 4.1 (Tanaka formula for G). Let T > 0 be such that [0,T ] R and c be real number. Assume ⊂ that Assumption (A ) holds and that the map t Rt is both continuous and of bounded variations on [0,T ] and such that: 7→ (i) λ( T )= α ( T )=0, ZR R ZR (ii) t R−1/2 L1([0,T ], dR ), 7→ t ∈ t 1/2 1 (iii) q N such that t gt′ q Rt− L ([0,T ], dt). ∃ ∈ 7→ | |− ∈ Then, the following equality holds in (L2): T T Gt c = c + sign(Gt c) d⋄Gt + δ c (Gt) dRt, (4.1) | − | | | Z0 − Z0 { } where the function sign is defined on R by sign(x) := 1R∗ (x) 1R (x). + − − Proof. This is a direct application of Theorem 3.4 with F : x x c . 7→ | − | Remark 4. Thanks to (3.11), we clearly see that condition (ii) is satisfied if t Rt is absolutely continuous. 7→

11 Comparison with Tanaka formula for general and particular Gaussian processes As Equality (4.10) will show, one can write (4.1) with the (weighted) local time instead of the last term on the right hand side. Because the approach used in [LN12] is more intrinsic (since all the assumptions are made exclusively on the variance function R), it is not easy to compare conditions (H1), (H3a) & (H4), required by [LN12, Theorem 4.3], to conditions (i) to (iii) of Theorem 4.1. However, one can note the following facts. It is clear that one does not need here to make any assumption on the covariance function (t,s) R but 7→ t,s only on the variance function t R . In [LN12] (t,s) R is assumed to be continuous on [0,T ]2, 7→ t 7→ t,s while only t Rt is assumed to be continuous on [0,T ] here. Moreover, we do not need to assume that 7→ T t Rt is increasing, nor such that R = 0 . Besides, the fact that t Rt is of bounded variations on [07→,T ] is sufficient here and one thusZ does{ not} need to make Assumptio7→n (H3a). Finally we do not need any assumption of the type of [LN12, (H4)] to get a Tanaka formula. Besides, in view of results of [Leb17, Theorem 5.6] and (3.11), it is clear that Theorem 4.1 applies to all Vγ - processes and, therefore, that Equality (4.1) (or if one prefers Equality (4.10) below) fully generalizes [MV05, Theorem 36] (with F : x (x c)+ instead of x x c ). Note moreover that Theorem 4.17→ also generalizes− the classic7→al | Tanaka− | formula for Brownian motion (see [RY99, Theorem 1.5]), the one given for fBm in [Ben03, Corollary 4.8] for any Hurst index H in (0, 1), the one provided in [BDM10, Theorem 12] for Volterra type mBm (it is assumed in there that h takes values in (1/2, 1)) and the one given in [LLV14, Theorem 6.1] for normalized mBm.

4.2 Local times of Gaussian processes in G In all this subsection one does not assume anymore that Assumptions (A ) holds. G is just assumed to be of the form (1.1) and such that the map g : t gt, defined at the beginning of Subsection 2.4, is continuous on R. One assumes, through this section,7→ that [0,T ] R. Let us first define the local time of G. ⊂ Definition 2. (non weighted) local time of G The (non weighted) local time of G at any point a R, up to time T , denoted ℓ(G)(a), is defined by: ∈ T (G) 1 ℓT (a) := lim λ( s [0,T ]; Gs (a ε,a + ε) ), ε 0+ 2ε { ∈ ∈ − } → where λ denotes the Lebesgue measure on R and where the limit holds in ( )∗, when it exists. S Having in mind the definition of generalized functional of G (given in Theorem-Definition 2.2) one has the following result. Proposition 4.2. [Integral representation on [0,T ] of non weighted local time] T Assume that the variance map s Rs is continuous on [0,T ] and such that the Lebesgue measure of R is equal to 0, then: 7→ Z

1. The map s δa(Gs) is ( )∗-integrable on [0,T ] for every a R∗. The map s δ0(Gs) is 7→ S −1/2 1 ∈ 7→ ( )∗-integrable on [0,T ] if s R belongs to L ([0,T ], ds). S 7→ s 2. The following equality holds in ( )∗, for every a R∗: S ∈ T (G) ℓT (a)= δa(Gs) ds, (4.2) Z0 Equality (4.2) still holds for a =0 if s R−1/2 L1([0,T ], ds). 7→ s ∈ T Remark 5. The reason why we have to assume that the Lebesgue measure of R is equal to 0 lies in the S T Z fact that, when F belongs to ′(R), F (Gt) is not defined if t belongs to R (Theorem-Definition 2.2). Z T However, for any a in R∗, if δ (G ) has a limit in ( )∗, as s tends to s , one may extend δ (G ) a s S 0 ∈ ZR a s at point s0 by setting δa(Gs0 ):=0. We can not say anything in general, when a =0. Before giving the proof of Proposition 4.2, let us state the following result, the proof of which is obvious and therefore left to the reader, since t R is a continuous function. 7→ t

12 Lemma 4.3. One has the following result: a T η S (R), a R∗, ε (0, | | ), γ (0,ε), such that: (s,y) D( ,γ) D(a,γ), ∀ ∈ ∀ ∈ ∀ ∈ 2 ∃ ∈ ∀ ∈ ZR × (y )2 < a2/16, (4.3) − − s −

T T where we have set D(a,γ) := (a γ,a + γ), D( R ,γ) := u [0,T ]; d(u, R ) < γ and where d(A, x) denote the distance between the subset− A of R andZ the real{x.∈ Z } Proof of Proposition 4.2. 1. Thanks to [Kuo96, Theorem 7.3], we know that s δa(Gs) is an ( )∗-process, the S-transform of which is equal, for every (a,η,s) R S (R) Γ , to:7→ S ∈ × × R 2 1 (a ) S(δa(Gs))(η)= exp − − , (4.4) √2πRs 2Rs  −1/2 where ΓR := u [0,T ]; Ru > 0 . When a = 0, the right hand side of (4.4) is bounded by Rs , which belongs, by assumption,{ ∈ to L1([0},T ], dt). Theorem 2.6 then applies and establishes 1. in the case where a = 0. Let a be in R∗. The weak measurability of t δa(Gt) is clear. In view of [Kuo96, p. 246], one just has 7→1 to show that t S(δa(Gt))(η) belongs to L ([0,T ], dt), for every η in S (R). Let η be in S (R) and a 7→ | | T ε be in (0, 2 ). In view of (4.4) and thanks to Lemma 4.3, one knows that, for every s0 in R , there 2 Z −1/2 a /32Rs exists a neighborhood s0 of s0 on which S(δa(Gs))(η) Rs e− =: la,η(s). Since la,η belongs to L1([0,T ], dt) for everyV η in S (R), one| has established|≤ 1. in the case a = 0 and therefore ends the proof of 1. 6 2. Let ε> 0, and a be a real number. Define T a( )∗ S Iε(a) := λ( s [0,T ]; Gs D(a,ε) ); and Φε(a) := 1(a ε,a+ε)(y) δy(Gs) ds dy. { ∈ ∈ } ZR Z − 0 

A simple application of Theorem 2.6 and Fubini’s theorem shows that Φε(a) is well defined, for every a R. Besides, since I (a)= T 1 (G ) ds, one has: ∈ ε 0 D(a,ε) s R T Iε(a)= 1(a ε,a+ε)(y) δy(Gs) dy ds, (4.5) Z0 ZR −  where we used the following lemma, the proof of which can be found in Appendix A.2.

2 Lemma 4.4. Let j S ′(R) be of function type and let g be in L (R). Define Ψ : R ( )∗ by setting ∈ j → S Ψ (y)= j(y) δ (<.,g >). If Ψ is ( )∗-integrable on R, then one has the equality j y j S ( )∗ j(y) δy(<.,g >) dy S= j(<.,g >). (4.6) ZR Moreover, if j(<.,g >) belongs to (L2), then (4.6) is also true in (L2).

A direct application of Fubini’s Theorem and (4.5) shows the equality of the S-transform of Φε(a) and I (a). The equality in ( )∗results of the injectivety of the S-transform. It then remains to show that ε S Iε(a) S lim exists in ( )∗. Besides, let (η,a) be in (R) R. Let ε be in (0, 1) such that D(a,ε) R+∗ or ε 0 2ε → S × ⊂ D(a,ε) R∗ , if a R∗. Using (ii) of Lemma 2.5 and (4.4) we get: ⊂ − ∈ T T Iε(a) 1 1 1 2 1 S( )(η)= exp − (y ) ds dy =: ρη(y,s) ds dy. 2ε 2εZ  Z √2πRs { 2Rs − }  2εZ Z D(a,ε) 0 D(a,ε) 0  (4.7) We will use the following lemma, the proof of which can be found in Appendix A.2. S T Lemma 4.5. For every η in (R), define the map Gη : R∗ R, by setting Gη(y) := 0 ρη(y,s) ds. −1/2 → 1 T Then Gη is continuous on R∗. Moreover, if s Rs belongs to L ([0,T ] R , ds), then oneR can extend T 7→ \Z Gη to R by setting Gη(0) := 0 ρη(0,s) ds. Furthermore, this extension is continuous on R. R

13 T The equality Gη(a) = S( 0 δa(Gs) ds)(η) and the continuity of Gη, on D(a,ε), allows us to write 1 T Iε(a) T lim Gη(y) dy = SR( δa(Gs) ds)(η). Thanks to (4.7), one finally gets lim S( )(η)= S( δa(Gs) ds)(η). ε 0 2ε D(a,ε) 0 ε 0 2ε 0 Moreover,→ R using the same computationsR as in the proof of Lemma 4.5, it is easy→ to show that, for everyR a in R∗, η η S (R), y D(a, a /2), G (y) M e| |0 , ∀ ∈ ∀ ∈ | | | η |≤ a 1 −1/2 2 Iε(a) η − 0 where Ma := T (8(1 + e) a + sup Rs exp( a /32Rs) ). One then gets: S( 2ε )(η) Ma e| | . | | { − (S)∗ } | | ≤ Iε(a) T  [Kuo96, Theorem 8.6] then allows us to write lim 2ε = δa(Gs)ds, which ends the proof. ε→0 0 R (G) Remark 6. We deduce from the proof of Proposition 4.2 that ℓT (a) exists for every a in R∗; and also −1/2 1 for a =0, if s Rs belongs to L ([0,T ], ds). Besides, the example of the Gaussian process (Gs)s R, ∈ 7→ 1 −1/2 1 defined by setting gs := k 0 s(1 + k)− ek, shows that the condition s Rs belongs to L ([0,T ], ds) ≥ 7→ is required in order thatP the map s δ (G ) is ( )∗-integrable on [0,T ]. 7→ 0 s S Theorem-Definition 4.1. [Weighted local time of G and its integral representation in ( )∗] S T Assume that the map s Rs is continuous and of bounded variations on [0,T ] and such that αR( R ) is equal to 0. Then: 7→ Z

1. The map s δa(Gs) is ( )∗- integrable on [0,T ] with respect to the measure αR, for every a R∗. 7→ S ∈ −1/2 The map s δ0(Gs) is ( )∗-integrable on [0,T ] with respect to the measure αR, if s Rs 17→ S 7→ belongs to L ([0,T ], dRs).

2. For every a R, when s δa(Gs) is ( )∗-integrable on [0,T ] with respect to the measure αR, one ∈ 7→ S L (G) can define the weighted local time of G at point a, up to time T , denoted T (a), as being the ( )∗- process defined by setting: S T L (G) T (a) := δa(Gs) dRs, (4.8) Z0 where the equality holds in ( )∗. S Proof. Since the proof of Point 1 follows exactly the same steps as proof of Proposition 4.2 (one just has to replace the Lebesgue measure by αR), it is left to the reader.

Assume that s Rs is absolutely continuous on [0,T ] and such that its derivative, denoted R′, is differentiable a.e.7→on [0,T ]. If the following conditions are fulfilled: (i) s R′ δ (G ) is ( )∗-integrable on [0,T ], 7→ s a s S (G) (ii) s R′′ ℓs (a) is ( )∗-integrable on [0,T ], 7→ s S

(G) (G) ′ (G) ′ (iii) s R′ ℓs (a) has a limit in 0 and T (denoted respectively ℓ (a) R and ℓ (a) R ), 7→ s s s s=0 s s s=T then an integration by parts yields: T (G) (G) (G) (G) L (a)= ℓ (a) R′ ℓ (a)R′ ℓ (a) R′′ ds, (4.9) T s s s=T − s s s=0 − Z s s 0 where R′′ denotes the second derivative of R. Remark 7. 1. In view of the definition of weighted local time, it is clear that (4.1) now reads: T L (G) Gt c = c + sign(Gt c) d⋄Gt + T (c). (4.10) | − | | | Z0 − As we stated at the end of Subsection 4.1, one recovers hence the Tanaka formula provided in [MV05, Theorem 36] and in [LN12, Theorem 4.3]. 2. The definition of weighted local time of G (Equality (4.8)) generalizes the one given, for fBm, in [Cou07, Definition 14]. Note moreover that Equality (4.9) above generalizes, to all Gaussian processes which belong to G , the relation between weighted and non-weighted local times given, for fBm, in [Cou07, Definition 14].

14 5 Occupation times formulas and Comparison to previous re- sults 5.1 Occupation times formulas The following theorem establishes that, under suitable assumptions, both weighted and non-weighted local times are (L2) random variables. It also provides an occupation time formula for both ℓ(G) and (G) L . Denote Mb(R) the set of positive Borel functions defined on R and recall that λ denotes the Lebesgue measure on R.

Theorem 5.1. Assume that the function s Rs is continuous on [0,T]. 7→ T 2 (i) Assume that λ ( T )=0 and that s R−1/2 L1([0,T ], ds). If E eiξGs ds dξ < + , R ZR 7→ s ∈ R 0 ∞ (G) 2  R R  then the map a ℓT (a) belongs to L (λ µ). Moreover one has the following equality, valid for µ-a.e. ω in Ω, 7→ ⊗ T M (G) Φ b(R), Φ(Gs(ω)) ds = ℓT (y)(ω) Φ(y) dy. (5.1) ∀ ∈ Z0 ZR ·

T (ii) Assume that αR( R )=0. Assume moreover that s Rs is of bounded variations on [0,T ], and Z T7→ 2 (G) such that s R−1/2 L1([0,T ], dR ). If E eiξGs dR dξ < + , then a L (a) 7→ s ∈ s R 0 s ∞ 7→ T belongs to L2(λ µ). Moreover one has the following R R equality, valid for µ-a.e. ω in Ω, ⊗ T M L (G) Φ b(R), Φ(Gs(ω)) dRs = T (y)(ω) Φ(y) dy. (5.2) ∀ ∈ Z0 ZR ·

(G) L (G) 2 Remark 8. 1. The first consequence of Theorem 5.1 is that both ℓT (a) and T (a) are (L ) random variables for every real a. In a forthcoming paper we discuss the joint continuity of processes (T,a) 7→ ℓ(G)(a) and (T,a) L (G)(a). In particular, under assumptions of local non-determinism3, it is shown T 7→ T (G) L (G) that both map a ℓT (a) and a T (a) are continuous (the continuity at point 0 requires that 1/2 7→ 1 7→ 1 t Rt− belongs to L ([0,T ], dt) or to L ([0,T ], dRt)). 2.7→Since E[X]= S(X)(0), for every (L2) random variable, it is easy to get, for almost every real a, the equalities: T T (G) 1 a2 (G) 1 a2 E[ℓ (a)] = exp ds;& E[L (a)] = exp dRs, T √2π R 2 Rs T √2π R 2 Rs Z0 s {− } Z0 s {− } when assumptions of Theorem 5.1 are fulfilled, −1/2 1 3. As we stated right after the proof of Theorem-Definition 4.1, the map t Rt belongs to L ([0,T ], dRt) 7→ T 2 if t R is absolutely continuous. In this latter case (5.2) holds under the single condition E eiξGs dR dξ < 7→ t R 0 s + , which is equivalent to T T (∆(t,s))−1/2dR dR < + , where ∆(t,s) := E[(G G R)2 ]R. As it is  ∞ 0 0 s t ∞ t − s shown in [GH80, (21.9) & (21.10)]R R any of these two latter conditions is necessary and sufficient to ensure the existence of a weighted local time, that is an (L2) random variable. This shows that the assumptions we made on the variance function R are almost minimal. 4. The condition t R−1/2 L1([0,T ], dt) L1([0,T ], dR ) is verified for all the Gaussian processes of 7→ t ∈ ∩ t “reference” (see Proposition 5.3 below for the case of Vγ - processes). Let us give the proof of Theorem 5.1 before discussing the results provided therein and before making the comparison with the results obtained in [LN12, §4].

(ω) Proof. Let us start with the proof of (i). Define, for every (ω,a)inΩ R, the map denoted fT : R R (ω×) (ω) → and, for every positive integer n, the family of maps denoted fT,n : R R by setting: fT (ξ) := 2 → 1 T iξGs(ω) (ω) ξ /n (ω) 2π 0 e ds and fT,n(ξ) := e− fT (ξ). Recalling the notations about Fourier transform defined in SectionR 2.1, define also: K(ω)(a) := f (ω) (a); & K(ω)(a) := f (ω) (a). (5.3) T F T T,n F T,n   3see [GH80, §3.24] and references therein for an overview of this notion

15 (ω) (ω) 2 (L2) (ω) (ω) It is clear that lim K = K , in L (R), since lim f = f . Denote KT,n(a) (resp. KT (a)) n + T,n T n + T,n T → ∞ (ω) (ω) → ∞ the random variable ω KT,n(a) (resp. ω KT (a)). 7→ (G) 7→ Step 1: lim KT,n(a)= ℓ (a), where the equality holds in ( )∗. n + T → ∞ S (ω) 1 2 (ω) (ω) iaξ Let ω in Ω. Since f belongs to L (R) L (R), one can write K (a)= f (ξ) e− dξ. For every T,n ∩ T,n R T,n η in S (R), using Lemma 2.5, we then have: R

T 2 1 iaξ ξ /n iξGs S(KT,n(a))(η)= e− − S(e )(η) ds dξ. ZR 2π  Z0  Using Equality (2.1) and Fubini’s theorem, that both obviously apply here, we get:

T T 1 2 2 1 1 ξ ( +Rs)+2iξ(a <η,gs>) 1 S(K (a))(η)= e− 2 n − dξ ds =: J (s) ds. T,n √2π √2π √2π n,η Z0  ZR   Z0

Using the change of variable u := ξ √ε + R , where ε := 2 , in J yields: · n s n n n,η (a <η,g >)2 T exp − s S(K (a))(η)= 1 {− 2(εn+Rs) } ds. T,n √2π √ε +R Z0 n s Using the same arguments that the one we used in the proof of Lemma 4.5 allows us to apply the Lebesgue dominated convergence. We therefore get, using (4.4) and then (4.2),

T T 2 1 (a <η,gs>) lim S(KT,n(a))(η)= exp − ds = S( δa(Gs) ds)(η), n + 2Rs → ∞ Z0 √2πRs {− } Z0 (G) = S(ℓT (a))(η). (5.4) The domination used in the proof of Lemma 4.5 applies here also and prove the existence of a function L, which belongs to in L1([0,T ]), such that:

T η 2 (p,n,η,a) N N∗ S (R) R, S(KT,n(a))(η) e| |p L(u) du. (5.5) ∀ ∈ × × × | |≤ · Z0 Formulas (5.4) and (5.5) allow us finally to apply [Kuo96, Theorem 8.6.] and we then establish that (G) lim KT,n(a)= ℓ (a), in ( )∗. n + T → ∞ S Step 2: Let us show that there exists a strictly increasing function ϕ : N∗ N∗ such that, for almost 2 → every real a, lim KT,ϕ(n)(a)= KT (a), in (L ). n + → ∞ 2 From the first step we know that, for µ - almost every ω, lim K(ω) K(ω) = 0. From the other n + k T,n − T kL2(R) 2 → ∞ 2 hand, one gets sup K(ω) K(ω) . Since, by assumption, ω K(ω) is an (L2) T,n L2(R) T L2(R) T L2(R) n N∗k k ≤ k k 7→ k k random variable,∈ the Lebesgue dominated convergence theorem applies. Thus we get the convergence (ω) (ω) 2 2 lim Ω KT,n KT 2 R dω = 0; that may also be written lim E (KT,n(a) KT (a)) da = n + k − kL ( ) n + ZR − → ∞ R → ∞   0. This shows that there exists a strictly increasing function ϕ : N∗ N∗ such that for almost every 2 → real a, lim KT,ϕ(n)(a) = KT (a), where the limit holds in (L ). Finally, results of steps 1 and 2 yield, n + → ∞ (G) (G) for almost every real a, KT (a) = ℓT (a). We have then proved that the map a ℓT (a) belongs to L2(λ µ). 7→ ⊗ Step 3: Proof of Equality (5.1): An application of the monotone class theorem allows us to establish (5.1) only for every Φ in S (R). q q For any Φ in S (R), denote Φ the function defined by setting Φ(ξ) := Φ( ξ) and, for any F in S ′(R), q q q − denote F the tempered distribution defined by < F , Φ >=< F, Φ >, for every Φ in S (R). For every ω (ω) (ω) T (ω) in Ω, define θ : S (R) R by setting θ (Φ) := Φ(G (ω)) ds. Since θ belongs to S ′(R), one T → T 0 s T can compute its Fourier transform, denoted (θ(ω)).R We obtain, for almost every real u: F T

16 T } (ω) iuGs(ω) (ω) (θT )(u)= e− ds =2π fT (u) (5.6) F Z0 } (ω) 1 (ω) (ω) 1 (ω) (ω) Besides, since θT = 2π ( (θT )) we get, using (5.6), θT = 2π (2πfT )= KT ; and hence, for all Φ in S (R), F F F (ω) (ω) θT (Φ) =. (5.7)

(G) (ω) Thanks to Step 2, we know that, for almost every real a, KT (a)= ℓT (a). The definition of θT as well as (5.7) then yield: T (G) Φ(Gs(ω)) ds = ℓT (a) Φ(a) da. Z0 ZR

(ω) 1 T iξGs(ω) The proof of (ii) is obvious in view of the one of (i), since one just has to take fT (ξ) := 2π 0 e dRs (ω) T R and θT (Φ) := 0 Φ(Gs(ω)) dRs. It is therefore left to the reader. R 5.2 Comparison with the existing results on Local times of general and particular Gaussian processes

Let us first recall that, for any Borel set of R+, the occupation measure of G on , denoted ν , is B B B defined for all in (R) by setting ν ( ) := λ( s ; Gs ). When ν is absolutely continuous with D B B D { ∈B ∈ D} B respect to the Lebesgue measure on R, one defines the local time of G on , noted LG( , x); x R , as being the Radon-Nikodym derivative of ν . Folllowing [GH80, Section 2]B we say in{ thisB case that∈ G }is (LT). When = [0,T ], we write L (T, x); Bx R instead of L ([0,T ], x); x R . B { G ∈ } { G ∈ } Remark 9. 1. In Theorem 5.1, instead of giving our proof of (i), one may think to use [GH80, Theorem (21.9)]. In particular, we know thanks to [GH80, (21.10)] (or the equivalent condition [GH80, (22.3)]), 2 that G is (LT) and that LG(T,a); a R belongs to L (λ µ). However we do not know, at this { (G∈) } ⊗ stage, that LG(T,a); a R and ℓT (a); a R are the same process (since one does not know, { (G) ∈ } { ∈ } 2 in particular that ℓT (a)) is a random variable, which belongs to (L ), for almost every real a). Note moreover that [GH80, (21.10) or (22.3)] can not be used in the proof of (ii) of Theorem 5.1 either; at least without any additional assumption on the variance function t Rt. In contrast, the proof we gave for (i) can be adapted for establishing proof of (ii). 7→ 2. Note that (5.2) allows one to get the occupation times formula given in [RY99, Corollary p. 224] in the case where G is a continuous Gaussian semimartingale4. It is also clear that (5.2) generalizes, to Gaussian non-, the occupation times formulas given in [RY99, Corollary p. 224]. In addition to the general references cited in the introduction of this section, [LN12, Section 4], is the only one recent paper in which local times of general Gaussian processes are studied at such a level of generality. The authors consider therein only what we have called here weighted local time, which they denote by Lt(x) and define, when it exists, as the density of the occupation measure:

t mt( ) := 1B(Gs) dRs, B (R). B Z0 ∀ ∈B

By the very definition of Lt(x), they obtained the occupation formula (5.2). They then show that Lt(x) is an (L2) random variable, by giving its chaos decomposition. In comparison with [LN12, Section 4], and as we already stated in Section 4.1, one first does not need 2 to assume that (t,s) Rt,s is continuous on [0,T ] , but only that t Rt and t gt are continuous 7→ 7→ 7→ T on [0,T ]. Besides, the fact that t Rt is of bounded variations on [0,T ] and such that αR( R )=0 is sufficient here and one thus does7→ not need to make Assumption (H3a) of [LN12] (i.e. toZ assume T that t Rt is increasing and such that R = 0 ). Finally, and as we stated in Point 3 of Remark 7→ T 2 Z { } 8, the assumption E eiξGs dR dξ <+ is necessary and sufficient to ensure the existence of a R 0 s ∞ weighted local time, that R isR an (L2) random  variable. It is then less restrictive than [LN12, (H4)]. Finally,

Assumption [LN12, (H4)] is not required here but one needs to assume that s R−1/2 L1([0,T ], dR ). 7→ s ∈ s 4See [Leb17, Proposition 5.7] to get more informations on the structure of Gaussian .

17 Let us now see the case of the Gaussian processes of “reference”.

(G) Example 5.2. 1. If G is a Brownian motion or a Brownian bridge, then ℓT (a), given in Definition 2, always exists and belongs to L2(R Ω, λ µ), as a two parameters process. Moreover, assumptions × ⊗ (B) L (B) of Proposition 4.2, Theorem-Definition 4.1 and Theorem 5.1 are all fulfilled. Of course, ℓT and T both coincide and we recover classical results for Brownian motion local time (see [RY99, Chapter 6] for example) in all its generality. 2. (i) If G is a fBm BH , then assumptions of Proposition 4.2, of Theorem-Definition 4.1 and Theorem (G) 5.1 are all fulfilled, whatever the value of H in (0, 1) is. Thus ℓT (a), given in Definition 2, always exists (G) 2 and ℓT (a)(ω); (a,ω) belongs to L (R Ω, λ µ). We therefore recover the results for non weighted local{ time given in [HØ02],} as well as the× occupation⊗ times formulas given, for both weighted and non weighted fractional local times, in [Cou07, Proposition 30 & Definition 14]. In particular, if we except the occupation time formula, the results we provide here on weighted local time when G = BH (i.e. on L (BH ) T ), are new. (ii) To our best knowledge, there was no proof of existence of (weighted or non-weighted) local time of general mBm (i.e. in sense of [LLVH14, Definition 1.2], recalled in our introduction) in the literature so far. However, several results exist for certain classes of mBm. More precisely, in [BDG06, BDG07] the authors are interested in the mBm, which corresponds to Y(1,0) in [ST06]. They establish, on each time interval [a,b] (0, + ), the existence of a square integrable local time (i.e. in L2(R Ω, λ µ)), assuming sup h(t); t ⊂0 < β∞ 1, where h is β-Hölder continuous. In [BDM10], a Volterra-type× ⊗ representation of mBm{ is studied.≥ } It∧ is established, using the sufficient condition given in [GH80, Theorem 22.1], the existence of a square integrable local time, (assuming h is a differentiable and (1/2, 1)-valued function). As a consequence, the authors get (5.1). Starting from the Itô formula [AMN01, Theorem 2], they also get (5.2). All these results are recovered by Theorem 5.1. In particular, the results we provide on weighted local time when G = Bh, are new. 3 Finally, we provide new results on weighted and non weighted local times of Vγ - processes in the Proposition 5.3 below.

γ Proposition 5.3. Let B be a Vγ - process. Then both weighted and non-weighted local times of G (denoted respectively ℓ(γ) and L (γ)) exist and belong to L2(R Ω, λ µ), as two parameters processes. T T × ⊗ Moreover, for µ-a.e. ω in Ω and all Φ in Mb(R), one has the following equalities: T (γ) Φ(Gs(ω)) ds = ℓT (y)(ω) Φ(y) dy. (5.8) Z0 ZR T L (γ) Φ(Gs(ω)) dRs = T (y)(ω) Φ(y) dy. (5.9) Z0 ZR

Proof. It is sufficient to apply Proposition 4.2, Theorem-Definition 4.1 and Theorem 5.1. Since t Rt is absolutely continuous (as it has been shown right above (3.11)), one sees, in view of Remark7→ 8, 2 2 T iξBγ T iξBγ that one just has to prove that: E e s dR dξ < + & E e s ds dξ < R 0 s ∞ R 0 + . However, thanks to [GH80, (21.10) R R & (22.3)], one knows that this is R equivalent R to show that: T∞ T 1/2 T T 1/2 γ 2 − − γ 0 0 (∆(t,s)) dt ds < + & 0 0 (∆(t,s)) dRs dRt < + , where ∆(t,s) := E[(Bt Bs ) ]. ∞ ∞ − 2 ThanksR R to [MV05, Proposition 1.1],R weR know that: R t s ∆(t,s) 2 R t s , for every (t,s) in [0,T ] . We then use the same arguments as the one given| right− | ≤ above (3.≤11),| as− well| as the fact that R is increasing, to get:

T T T t 2 1/2 1/2 T ΘT := (∆(s,u))− du ds 2 Ru− du dt 8(1 + RT )+ , Z Z ≤ Z Z ≤ √α 0 0 0 0  where α := inf u (0,T ], R′ 1 (where we set inf =+ ). This achieves the proof. { ∈ u ≤ } ∅ ∞ Acknowledgments

I want to express my deep gratitude to Jacques Lévy Véhel for his advices and for the very stimulating discussions we had about this work. I also want to thank Professor T. Hida for his warm welcome at the University of Nagoya,

18 where a part of this paper was written, as well as Professor L. Chen and the Institute for Mathematical Sciences of Singapore (NUS), where another part of this paper was written. This work is dedicated to the memory of Professor Marc Yor. A Appendix Appendix A.1 Bochner Integral A.1 Bochner Integral The following notions about Bochner integral come from [HP57, p.72, 80 and 82] and [Kuo96, p.247].

Definition 3. Bochner integral [Kuo96, p.247] Let I be a Borelian subset of R endowed with the Lebesgue measure. One says that Φ : I ( )∗ is Bochner integrable on I if it satisfies the two following conditions: → S 1 Φ is weakly measurable on I i.e u << Φu, ϕ >> is measurable on I for every ϕ in ( ). 7→ S 1 2 p N such that Φu ( −p) for almost every u I and u Φu belongs to L (I). ∃ ∈ ∈ S ∈ 7→ k k−p I ds The Bochner-integral of Φ on is denoted I Φs .

∗ R Proposition A.1. If Φ : I ( ) is Bochner-integrable on I then there exists an integer p such that Φs ds → S I −p ≤ Φs ds. Moreover Φ is also Pettis-integrable on I and both integrals coincide on I. R I k k−p R Remark 10. The previous proposition shows that there is no risk of confusion by using the same notation for both Bochner and Pettis integrals.

(n) ∗ (n) Theorem A.2. Let p N and (Φ ) be a sequence of processes from I to ( ) such that Φu ( −p) for ∈ n∈N S ∈ S almost every u I and for every n. Assume moreover that Φ(n) is Bochner-integrable on I, for every n, and ∈ (m) (n) ∗ that lim Φs Φs ds = 0. Then there exists an ( ) -process (almost surely ( −p)-valued), (n,m)→(+∞,+∞) I − −p S S R denoted Φ, defined and Bochner-integrable on I, such that (n) lim Φs Φs −p ds =0 (A.1) n→+∞ ZI k − k ∗ Furthermore, if there exists an ( ) -process, denoted Ψ, which verifies (A.1), then Ψs = Φs for a.e. s in I. Finally (n) S ∗ one has lim Φs ds = Φs ds, where the equality and the limit both hold in ( ) . n→+∞ I I S A.2 BochnerR IntegralR A.2 Proof of Lemmas 4.4 and 4.5 Proof of Lemma 4.4. The existence of the left hand side of (4.6) results from (4.4) and Theorem 2.6. Using (ii) of Lemma 2.5 and again (4.4), we get, for every η in S (R), 1 2 S j y δ <.,g > dy η j y S δ <.,g > η dy j y −(y−) dy ( ( ) y( ) )( )= ( ) ( y( )( ) = ( ) exp 2|g|2 ZR ZR √2π g ZR 0 | |0  1 2 j y −(y−) dy S j <.,g > η , = ( ) exp 2|g|2 = ( ( ))( ) √2π g ZR 0 | |0  where the last equality results from [Kuo96, Theorem 7.3]. The conclusion comes from the injectivity of S- transform. 

Proof of Lemma 4.5. Let (a, η) in R∗ S (R) and ε be in (0, |a| ). Among the assumptions that are required to apply the theorem × 2 of continuity under sign , only the domination of ρη is not obvious and will therefore be stablished here. Let T −1/2 γ be chosen in (0,ε) suchR that (4.3) holds. Denote Γγ := u [0,T ]; d(u, R ) γ , mγ := sup Ru ; u { ∈ Z a≥2 } { ∈ T −1/2 − 32R [0,T ], s.t. d(u, R ) γ and, for every r in (0,T ], denote Mr := sup Rs e s ; s [0,r] . Using Lemma Z ≥ } T { ∈ } 4.3, one gets, for every (s,y) in D( ,γ) D(a,γ), ZR ×

ρη(y,s) = ρη(y,s) 1Γγ (s)+ ρη(y,s) 1D(ZT ,γ)(s) mγ 1Γγ (s)+ Mγ 1D(ZT ,γ)(s) =: f(s). | | | | | | R ≤ R Since f belongs to L1([0,T ],ds), the domination is complete. The case a = 0 is simpler. Indeed, one has:

∗ −1/2 (ε, η) R S (R), (y,s) D(0,ε) [0,T ], ρη(y,s) R . ∀ ∈ × ∀ ∈ × | |≤ s −1/2 1 T Since s Rs belongs to L ([0,T ] ,ds) by assumption, the domination is now complete and the theorem 7→ \ZR of continuity under sign applies here and ends the proof.  R 19 References

[AMN01] E. Alòs, O. Mazet, and D. Nualart. Stochastic calculus with respect to Gaussian processes. Ann. Probab., 29(2):766–801, 2001. [BDG06] B. Boufoussi, M. Dozzi, and R. Guerbaz. On the local time of multifractional Brownian motion. Stochastics, 78(1):33–49, 2006. [BDG07] B. Boufoussi, M. Dozzi, and R. Guerbaz. Sample path properties of the local time of multifractional Brownian motion. Bernoulli, 13(3):849–867, 2007. [BDM10] B. Boufoussi, M. Dozzi, and R. Marty. Local time and Tanaka formula for a Volterra-type multifractional Gaussian process. Bernoulli, 16(4):1294–1311, 2010. [Ben03] C. Bender. An Itô formula for generalized functionals of a fractional Brownian motion with arbitrary Hurst parameter. Stochastic Processes and their Applications, 104:81–106, 2003. [BJR97] A. Benassi, S. Jaffard, and D. Roux. Elliptic Gaussian random processes. Rev. Mat. Iberoamericana, 13(1):19– 90, 1997. [CG62] I.M. Chilov and G.E. Gelfand. Les distributions, volume 2. Dunod, 1962. [CNT01] L. Coutin, D. Nualart, and C. A. Tudor. Tanaka formula for the fractional Brownian motion. . Appl., 94(2):301–315, 2001. [Cou07] L. Coutin. An Introduction to (Stochastic) Calculus with Respect to Fractional Brownian Motion. Séminaire de Probabilités XL, 1899:3–65, 2007. [EVdH03] R.J. Elliott and J. Van der Hoek. A general fractional white noise theory and applications to finance. Mathe- matical Finance, 13(2):301–330, 2003. [GH80] D. Geman and J. Horowitz. Occupation densities. Ann. Probab., 8(1):1–67, 1980. [Hid75] T. Hida. Analysis of Brownian functionals. Carleton Univ., Ottawa, Ont., 1975. Carleton Mathematical Lecture Notes, No. 13. [HKPS93] T. Hida, H. Kuo, J. Potthoff, and L. Streit. White Noise. An Infinite Dimensional Calculus, volume 253. Kluwer academic publishers, 1993. [HØ02] Y. Hu and B. Øksendal. Chaos expansion of local time of fractional Brownian motions. Stochastic Anal. Appl., 20(4):815–837, 2002. [HØS05] Y. Hu, B. Øksendal, and D.M Salopek. Weighted local time for fractional Brownian motion and applications to finance. Stoch. Anal. Appl., 23(1):15–30, 2005. [HP57] E. Hille and R.S. Phillips. Functional Analysis and Semi-Groups, volume 31. American Mathematical Society, 1957. [Jan97] S. Janson. Gaussian Hilbert Spaces, volume 129 of Cambridge Tracts in . Cambridge University Press, Cambridge, 1997. [Kol40] A. Kolmogorov. Wienersche Spiralen und einige andere interessante Kurven in Hilbertsche Raum. C. R. (Dokl.) Acad. Sci. URSS, 26:115–118, 1940. [Kuo96] H.H. Kuo. White Noise Distribution Theory. CRC-Press, 1996. [Leb17] J. Lebovits. Stochastic Calculus with respect to Gaussian Processes, Part I, Preprint. 2017. Available at the adresse https://docs.google.com/viewer?a=v&pid=sites&srcid=ZGVmYXVsdGRvbWFpbnxsZWJvdml0c2pvYWNoaW18Z3g6NDEwYTU1NzBhODYzYmM5Yw . [LLV14] J. Lebovits and J. Lévy Véhel. White noise-based stochastic calculus with respect to multifractional Brownian motion. Stochastics An International Journal of Probability and Stochastic Processes, 86(1):87–124, 2014. [LLVH14] J. Lebovits, J. Lévy Véhel, and E. Herbin. Stochastic integration with respect to multifractional Brownian motion via tangent fractional Brownian motions. Stochastic Process. Appl., 124(1):678–708, 2014. [LN12] P. Lei and D. Nualart. Stochastic calculus for Gaussian processes and application to hitting times. Commun. Stoch. Anal., 6(3):379–402, 2012. [MR06] M.B. Marcus and J. Rosen. Markov processes, Gaussian processes, and local times, volume 100 of Cambridge Studies in Advanced Mathematics. Cambridge University Press, Cambridge, 2006. [MV05] O. Mocioalca and F. Viens. Skorohod integration and stochastic calculus beyond the fractional Brownian scale. J. Funct. Anal., 222(2):385–434, 2005. [MVN68] B. Mandelbrot and J.W. Van Ness. Fractional Brownian motions, fractional noises and applications. SIAM Rev., 10:422–437, 1968. [MWX08] M. Meerschaert, D. Wu, and Y. Xiao. Local times of multifractional Brownian sheets. Bernoulli, 14(3):865–898, 2008. [Nua06] D. Nualart. The and Related Topics. Springer, 2006. [PLV95] R. Peltier and J. Lévy Véhel. Multifractional Brownian motion: definition and preliminary results, 1995. rapport de recherche de l’INRIA, n0 2645. [RY99] D. Revuz and M. Yor. Continuous martingales and Brownian motion, volume 293 of Grundlehren der Math- ematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, third edition, 1999.

20 [Si12] Si Si. Introduction to Hida distributions. World Scientific Publishing Co. Pte. Ltd., Hackensack, NJ, 2012. [ST06] S. Stoev and M. Taqqu. How rich is the class of multifractional Brownian motions? Stochastic Processes and their Applications, 116:200–221, 2006. [SV14] T. Sottinen and L.. Viitasaari. Stochastic Analysis of Gaussian Processes via Fredholm Representation. 2014. Preprint, arXiv:1410.2230. [Tha93] S. Thangavelu. Lectures of Hermite and Laguerre expansions. Princeton University Press, 1993.

21