arxiv: v1 [math.pr] 26 Jul 2014

Transcription

1 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION arxiv:47.79v [math.pr] 26 Jul 24 W LODEK BRYC Abstract. For a parameter < q <, we use the Jackson q-integral to define integration with respect to the so called q-brownian motion. Our main results are the q-analogs of the L 2 -isometry and of the Ito formula for polynomial integrands. We also indicate how the L 2 -isometry extends the integral to more general functions. Note This is the expanded version of the manuscript with additional material which is marked by the gray background like in this note. Please note that the numbering of formulas in this version does not match the shorter version!!!. Introduction In this paper we use the Jackson q-integral to define an operation of integration with respect to a Markov process which first arose in non-commutative probability and which we shall call the q-brownian motion. The non-commutative q-brownian motion was introduced in [4] who put the formal approach of [] on firm mathematical footing. Stochastic integration with respect to the non-commutative q-brownian was developed in Ref. [9]. According to [3, Corollary 4.5], there exists a unique classical Markov process (B t [, with the same univariate distributions and the same transition operators as the non-commutative q-brownian motion. Markov process (B t is a martingale and which converges in distribution to the Brownian motion as q. With some abuse of terminology we call (B t the q-brownian motion and we review its basic properties in Section 2. Our definition of integration with respect to (B t uses the orthogonal martingale polynomials and mimics the well-known properties of the Ito integral. The integral of an instantaneous function f(b s,s of the process is denoted by f(b s,sdlb, where to avoid confusion with the Ito integral with respect to the martingale (B t we use dl instead of d. Date: Created: Dec. 23. Revised: May 24 and July 24. Printed July 3, 28. (File: q-stoch-int-23-v3.tex. Key words and phrases. q-brownian motion; stochastic quantum calculus.

2 2 W LODEK BRYC We define this integral for < q < and for polynomials f(x,t in variable x with coefficients that are bounded for small enough t. Our definition uses the Jackson integral [5] which is reviewed in Section 3. This approach naturally leads to the q-analog of the L 2 -isometry (. E( f(b s,sdlb 2 ( f(b E s,s 2 d q s, with the Jackson q-integral appearing on the right hand side. In Theorem 4.3 we establish this formulafor polynomial f(x,t and in Corollary4.4 we use it to extend the integral to more general functions. Proposition 4.5 shows that the integral is well defined for analytic functions f(x that do not depend on t. In Section 5 we use Corollary 4.4 to exhibit a solution of the linear q-equation dlz azdlb. Our second main result is a version of Ito s formula which takes the form (.2 f(b t,t f(, ( (s x f(b s,sdlb + (D q,s f(b qs,sd q s+ ( (s x f(b qs,sd q s. In Theorem 4.6, this formula is established for polynomial f(x,t and q (,. In Remark 4.9 we point out that we expect this formula to hold in more generality. As in the Ito formula, the operators D q,s, (s x and (s x should be interpreted as acting on the appropriate variable of the function f(x,s, which is then evaluated at x B s or at x B qs. The time-variable operator D q,s f is the q-derivative with respect to variable s: (.3 (D q,s f(x,s f(x,s f(x,qs. ( qs This expression is well defined for s > which is all that is needed in (.2 when q (,. The operators (s x and x (s act only on the space variable x but they depend on the time variable s, and are defined as the singular integrals with respect to two time-dependent probability kernels: (.4 ( (s x f(x,s (.5 ( (s x f(x,s R f(y,s f(x,s ν x,s (dy y x (y xf(z,s+(x zf(y,s+(z yf(x,s R (x y(y z(z x 2 µ x,s (dy,dz. Probability measures ν x,s (dy and µ x,s (dy,dz can be expressed in terms of the transition probabilities P s,t (x,dy of the Markov process (B t t [, as follows: and ν x,s (dy P sq 2,s(qx,dy µ x,s (dy,dz P qs,s (x,dyν y,s (dz P qs,s (x,dyp q 2 s,s(qy,dz. Transition probabilities P s,t (x,dy appear in (2.2 below. The same probability kernel ν x,s (dy appears also in [, Proposition 2].

3 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 3 As q probability measures ν x,s (dy and µ x,s (dy,dz convergeto degenerated measures. Using Taylor expansion one can check that for smooth enough functions ( (s x f(x,s converges to f/ x and ( (s x f(x,s converges to 2 2 f/ x 2 as q., so (.2 is a q-analog of the Ito formula. We note that there seems to be some interest in q-analogs of the Ito formula. In particular, Ref. [, 2] discuss formal q-versions of the Ito formula and its applications to financial mathematics... Limits as q. We conclude this introduction with the verification that as q, formulas (. and (.2 coincide with the usual formulas of Ito calculus, at least on polynomials. Formula (. becomes the familiar ( 2 E f(w s,sdw E(f(W s,s 2 ds, as the Jackson q-integral becomes the Riemann integral in the limit. As q formula (.2 becomes ( t ( f(w t,t x f (W s,sdw+ s f (W s,sds+ ( 2 2 x 2f (W s,sds, but this is less obvious. Clearly, D q,s f becomes D s in the limit. Since P q 2 s,s(qx,dy δ x (dy in distribution as q and f(y f(x y x f (x + 2 f (x(y x +..., we see that (s x becomes x, when acting on smooth enough functions of variable x. Next, we expand the numerator of the integrand in (.4 into the Taylor series at x (here we take f to be a polynomial in variable x. We get (.6 (y xf(z+(x zf(y+(z yf(x (y x [f(x+f (x(z x+ 2 ] f (x(z x 2 +(x z [f(x+f (x(y x+ 2 ] f (x(y x 2 +(z yf(x +R(x,y,x 2 (x y(y z(z xf (x+r(x,y,z.

4 4 W LODEK BRYC After inserting (.6 into (.5, the first term gives 2 x, since µ 2 x,s (dy,dz is a probability measure. The reminder term R(x, y, z consists of the higher order terms. The coefficient at f (k (x/k! in R is (y x(z x k +(x z(y x k (y x(z x ( (z x k (y x k k 2 (x y(y z(z x (z x j (y x k 2 j, andit isclearthat fork 3wehavelim q (z x j (y x k 2 j µ x,s (dy,dz j, since µ x,s (dy,dz D δ x (dyδ x (dz in distribution as q q-brownian motion With B, the univariate distribution of B t is a q-gaussian distribution (2. γ t;q (dy q 2π 4t ( qy 2 k ((+q k 2 ( q y2 t qk ( q k+ dy, supported on the interval y 2 t/ q. Here, parameter t is the variance. (We note that there are several other distributions that are also called q-gaussian, see the introduction to [8]. The transitionprobabilitiesp s,t (x,dy arenon-homogeneousand, as noted in [5], arewelldefined forallx R. The explicit formwillnot beneeded in thispaper, but forcompletenesswenotethat for x 2 s/ q ands < t transitionprobability P s,t (x,dy has density supported on the interval y 2 t/ q which can be written as an infinite product (2.2 P s,t (x,dy q 2π 4t ( qy 2 k k (t sq k ( q k+( t(+q k 2 ( qy 2 q k (t sq 2k 2 ( qq k (t+sq 2k xy +( q(sy 2 +tx 2 q 2kdy. With B, we have γ t;q (dy P,t (,dy and the separable version of the process satisfies B t 2 t/ q for all t almost surely. Formulas (2. and (2.2 are taken from [7, Section 4.], but they are well known, see [3, Theorem 4.6] and [, 6, 2]. We note that for x > 2 s/ q transition probability P s,t (x,dy has an additional discrete component; for q this discrete component is explicitly written out in [2, Section 5.3]. As q, the transition probabilities P s,t (x;dy converge in variation norm to the (Gaussian transition probabilities of the Wiener process. This can be deduced from the convergence of densities established in the appendix of [4]. Our definition of the integral relies on the so called continuous q-hermite polynomials {h n (x;t : n,,...}. These are monic polynomials in variable x that are defined by the three step recurrence (2.3 xh n (x;t h n+ (x;t+t[n]h n (x;t,

5 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 5 where [n] [n] q + q + + q n, h (x;t, h (x;t x. (These are renormalized monic versions of the standard continuous q-hermite polynomials as given in [3, Section 3.] or in [7, Section 3.26]. Polynomials{h n (x;t} playaspecialrolebecausetheyareorthogonalmartingale polynomials for (B. The orthogonality is (2.4 h n (x;th m (x;tγ t;q (dx δ m,n [n]!t n, where [n]! [n] q! [] q [2] q...[n] q. The martingale property (2.5 h n (x;s h n (y;tp s,t (x,dy holds for all real x and s < t. When q recursion (2.3 becomes the recursion for the Hermite polynomials, which are the martingale orthogonal polynomials for the Brownian motion [9]. It is easy to check from (2.3 that (2.6 h n (x;t t n/2 h n (x/ t;. From (2.6 it is clear that if x 2 t/ q then (2.7 h n (x;t C n t n/2 for some constant C n. Explicit sharp bound with C n ( q n/2 n deduced from [3, (3..]. 3. The Jackson q-integral k [n k ] can be For reader s convenience we review basic facts about the Jackson integral [5]. In addition to introducing the notation, we explicitly spell out the technical assumptions that suffice for our purposes. In particular, since we will integrate functions defined on [, only, we assume that q (,. We follow [6] quite closely, but essentially the same material can be found in numerous other sources such as [3, Section.4]. Suppose a : [, R is bounded in aneighborhoodof, andlet b : [, R be such that (3. b(t b( Ct δ for some C < and δ > in a neighborhood of. Definition 3.. Fort > andq (,, the Jacksonq-integralofa(s with respect to b(s is defined as (3.2 a(sd q b(s : a(q k t ( b(q k t b(q k+ t. k When b(t t, formula (3.2 takes the following form that goes back to [5] (3.3 a(sd q s ( qt q k a(q k t. Recalling (.3, it is clear that (3.2 is a(s(d q,sb(sd q s. We will also need the q-integration by parts formula [5, 5]. This formula is derived in [6, page 83] under the assumption of differentiability of functions a(t, b(t. Ref. [3, Theorem.4.] gives a q-integration by parts formula under k

6 6 W LODEK BRYC more general assumptions, but since the conclusion is stated in a different form than what we need, we give a short proof. Proposition 3.2. For any pair of functions a,b that satisfy bound (3. near t we have (3.4 a(sd q b(s a(tb(t a(b( b(qsd q a(s. Proof. Underourassumptionsona,b,theseriesdefining a(sd qb(sand b(qsd qa(s converge, so we can pass to the limit as N in the telescoping sum identity N a(q n t ( b(q n t b(q n+ t N + b(q n+ t ( a(q n t a(q n+ t n n a(tb(t a(q N+ tb(q N+ t. We also note the q-anti-differentiation formula, which is a special case a(s of formula (3.4 (3.5 b(t b( (4. (D q,s b(sd q s. 4. Integration of polynomials with respect to B t Our definition is designed to imply the relation h n (B s,sdlb [n+] q h n+ (B t ;t and relies on expansion into polynomials {h n (x;t}. Let q (, and suppose that f(x,t d m a m(tx m is a polynomial of degree d in variable x with coefficients a m (t that depend on t. Then, for t >, we can expand f into the q-hermite polynomials (2.3, d b m (t (4.2 f(x,t [m]! h m(x;t. m Lemma 4.. If all the coefficients a m (t of the polynomial f(x,t are bounded in a neighborhood of then the coefficients b m (t in expansion (4.2 are bounded in the same neighborhood of. If all the coefficients a m (t satisfy estimate (3. then coefficients b m (t satisfy estimate (3.. Proof. Since h,...,h m are monic, they span the polynomials of degree m so that by orthogonality (2.4 we have x j h m (x,tγ t;q (dx for j < m. For the same reason, x m h m (x,tγ t;q (dx h 2 m(x,tγ t;q (dx t m [m]!. Thus the coefficients are b m (t t m f(x,th m (x;tγ t;q (dx t m t m/2 d jm d jm a j (t a j (t x j h m (x;tγ t;q (dx x j h m (x/ t;γ t;q (dx,

7 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 7 where in the last line we used (2.6. Since random variable B t / t has the same distribution as B, changing the variable of integration to y x/ t we get b m (t d a j (tt (j m/2 y j h m (y;γ ;q (dy jm [m]!a m (t+ d jm+ a j (tt (j m/2 y j h m (y;γ ;q (dy. Thus b m ( [m!]a m ( and for < t < we get b m (t b m ( [m]! a m (t a m ( +C t. Definition 4.2. Let f(x,t be a polynomial in variable x with coefficients bounded in a neighborhood of t. We define the dlb integral as the sum of the Jackson q-integrals with respect to random functions h m+ (B s ;s as follows: (4.3 f(b s,sdlb d m [m+]! b m (sd q h m+ (B s ;s. Since B t 2 t/ q for all t almost surely, from Lemma 4. and inequality (2.7 we see that the Jackson q-integrals on the right hand side of (4.3 are well defined. As a special case we get the following formula for the integrals of the deterministic functions. Example 4.. For non-random integrands, we have f(sdlb f(sd q B s. Indeed, since h the expansion (4.2 is f(sh (x;s and since h (x;s x, we get d q h (B s ;s d q B s. By computing the fourth moments one can check that the distribution of random variable f(sd qb s in general is not γ σ 2,q. When the coefficients of f(x,t satisfy estimate (3., we can use (3.4 to write the dlb integral as (4.4 f(b s,sdlb d m b m (t [m+]! h m+(b t ;t d m [m+]! Taking b m (t or we get formula (4.. Formula (4.4 leads to a couple of explicit examples. h m+ (B qs ;qsd q b m (s. Example 4.2. Noting that x h, x 2 h 2 (x;t+t, and h 3 (x;t x 3 (2+qxt, we get the following q-analogs of the formulas that are well known for the Ito integrals. B s dlb ( (B t 2 t, +q

8 8 W LODEK BRYC (B s 2 dlb (+q(+q +q 2 (B t 3 2+q (+q(+q +q 2 tb t + sd q B s. We remark that a similar definition could be introduced for a wider class of the q-meixner processes from [5, 7], compare [8, Theorem 7]. 4.. The L 2 -isometry. Our first main result is the q-analog of the Ito isometry. Theorem 4.3. If < q < and f(x,t is a polynomial in x with coefficients bounded in a neighborhood of t then (. holds. Proof. With (4.2, by orthogonality (2.4 for the q-hermite polynomials, we have So the right hand side of (. is ( ( 2 E f(b s,s d b 2 m(ss m. [m]! m (4.5 d m [m]! b 2 m(ss m d q (s d m [m+]! b 2 m(sd q (s m+. On the other hand, expanding the right hand side of (4.3 we get E ( ( f(b s,sdlb 2 d m,n (( E h n+ (B ;tq k h tq k n+ (B ;tq k+ tq k+ [m+]![n+]! ( k,j b m (q k tb n (q j t h m+ (B tq j ;tq j h m+ (B Since h n (B t ;t is a martingale in the natural filtration, for j < k we have ( E h m+ (B tq ;tq j h j m+ (B tq ;tq j+ B,B. j+ tq k tq k+ tq ;tq j+. j+ Similarly, ( E h n+ (B ;tq k h tq k n+ (B ;tq k+ B tq k+ tq,b j tq j+ for j > k. So the only contributing terms are j k and E ( ( f(b s,sdlb 2 d m,n (( E h n+ (B ;tq k h tq k n+ (B ;tq k+ tq k+ b m (q k tb n (q k t [m+]![n+]! k ( h m+ (B tq k ;tq k h m+ (B tq k+ ;tq k+. Noting that (( ( E h n+ (B ;tq k h tq k n+ (B ;tq k+ tq k+ h m+ (B ;tq k h tq k m+ (B ;tq k+ tq k+

9 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 9 for m n, we see that ( ( 2 E f(b s,sdlb d m [m+]! 2 k ( ( 2 b 2 n(q k te h m+(b ;tq k h tq k m+ (B ;tq k+ tq k+ d m [m+]! b 2 n(q k tt m+( q k(m+ q (k+(m+. k Here we used the fact that for s < t, by another application of the martingale property and (2.4, we have ( ( 2 ( ( E h k (B t ;t h k (B s ;s E h 2 k (B t ;t E h 2 k (B s ;s [k]!(t k s k. This shows that the right hand side of (4.3 matches (4.5. Suppose that for every t > the function f(x,t as a function of x is square integrable with respect to γ t;q (dx so that we have an L 2 (γ t;q (dx-convergent expansion (4.6 f(x,t n b n (t [n]! h n(x;t. Corollary 4.4. If each of the coefficients b n (t in (4.6 is bounded in a neighborhood of and the series t b 2 [n]! n(ss n d q s n converges, then the sequence of random variables n b k (s (4.7 [k]! h k(b s ;sdlb k converges in mean square as n. (4.8 When the mean-square limit (4.7 exists, it is natural to denote it as f(b s,sdlb. We remark that the verification of the assumptions of Corollary 4.4 may not be easy. In Section 5 we consider function ( f(x,t ( qq k x+tq 2k k and we rely on explicit expansion(4.6 to showthat the integral(4.8 is well defined for t < /( q. The following result deals with analytic functions that do not depend on variable t. Proposition 4.5. Suppose f(x k a kx k with infinite radius of convergence. Then f(b s dlb is well defined for all t >.

10 W LODEK BRYC Proof. Since the support of γ t;q (dx is compact and f is bounded on compacts, it is clear that we can expand f(x into (4.6. We will show that the coefficients b n (t of the expansion are bounded in a neighborhood of. (The proof shows also that b n (t satisfy estimate (3.; the later is not needed here, but this property is assumed in Theorem 4.6 below. As in the proof of Lemma 4., we have b n (t t n f(xh n (x;tγ t;q (dx t n/2 h n (y; a k t k/2 y k γ ;q (dy. Since on the support of γ ;q (dy the series in the integrand converges absolutely and is dominated by a constant that does not depend on y, we can switch the order of summation and integration so that b n (t t n/2 a k t k/2 h n (y;y k γ ;q (dy k kn k t n/2 a k t k/2 h n (y;y k γ ;q (dy, where in the last line we used the orthogonality of h n (y; and y k for k < n. Choose R > such that R 2 > 4t/( q. Expressing the coefficients a k by the contour integrals we get b n (t 2πi t n/2 kn z R 2πi f(z z k+dztk/2 z R f(z kn h n (y;y k γ ;q (dy t k/2 n/2 z k+ t k/2 h n (y;y k γ ;q (dydz. Since y 2/ q on the support of γ ;q (dy, from (2.4 we get ( k ( k h n (y;y k γ ;q (dy 2 2 h n (y; γ ;q (dy [n]!. q q Summing the resulting geometric series under the dz-integral, we get ( n (4.9 b n (t max f(z 2 [n]! z R R R q q R q 2 t. This shows that b n (t is bounded in any neighborhood of t. Next, we use (3.3 and inequality (4.9 to prove the convergence of the series. b 2 ( qt n+ n n [n]! (ssn d q s b 2 n n [n]! (qk tq (n+k n k ( n 4t n ( n C R 2 ( q q k(n+ 4t q C ( q R 2 ( q q n+ n k n ( n 4t C R 2 <, ( q n with constant C tr 2 ( qmax z R f(z 2 /(R q 2 t 2.

11 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION We now apply Corollary 4.4 to infer that f(b s dlb is well defined. Since R was arbitrary, the conclusion holds for all t > The q-ito formula. Our second main result is the q-version of the Ito formula. Theorem 4.6. Let < q <. Suppose f(x,t is a polynomial in x with coefficients that satisfy estimate (3. in a neighborhood of t. Then (.2 holds Proof of Theorem 4.6. By linearity, it is enough to consider a single term in (4.4. We first consider the constant term, f(x,t b(t. Then x f as the expression f(x,s f(y,s under the integral in the definition (.4 vanishes. Similarly, x f, asthe expression(z xf(y,s+(x yf(z,s+(y zf(x,s in the definition (.5 vanishes. So (.2 in this case reduces to the q-integral identity b(t b( (D q,sb(sd q s which was already recalled in (3.5. Now we consider a non-constant term of degree m+. We write (4.4 as (4. b(th m+ (B t ;t [m+] q We first note the following identity. b(sh m (B s ;sdlb + h m+ (B qs ;qs(d q,sb(sd q s. Lemma 4.7. For a fixed s >, we have (s x (h m+ (x;s [m+]h m (x;s. Proof. This is a special case of [5, Lemma 2.3], applied to τ θ. Thus, noting that (s x acts only on variable x, we have (4. b(th m+ (B t ;t x (b(sh m+(x;s (s xb s dlb + h m+ (B qs ;qs(d q,sb(sd q s. Next, we consider the second term on the right hand side of (4.. Using the q-product identity (4.2 D q,s (b(sh(s b(sd q,s (h(s+h(qsd q,s (b(s, we write h m+ (x;qs(d q,s b(s D q,s (b(sh m+ (x,s b(sd q,s (h m+ (x,s. This recovers the second term in (.2: (4. becomes (4.3 b(th m+ (B t ;t + (s D q,s (b(sh m+ (x;s x (b(sh m+ (x;s xb qs d q s+ D (s xb s dlb x (b(sh m+ (x;s xb qs where D x (s is a linear operator on polynomials in variable x, whose action on monomials x n n j a j(sh j (x;s is defined by D (s x (xn n a j (sd q,s (h j (x;s. j d q s,

12 2 W LODEK BRYC It remains to identify D (s x with (s x. To do so we use the approach from [5]. We fix s > and introduce an auxiliary linear operator A that acts on polynomials p(x in variable x by the formula (4.4 A(p(x D (s x (xp(x xd (s x (p(x. Lemma 4.8. (4.5 A(h m (x;s [m] q h m (x;qs. Proof. By definition, A(h m (x;s D (s x (xh m (x;s xd (s x (h m (x;s, and we evaluate each term separately. From recursion (2.3 we get (4.6 D (s x (xh m (x;s D q,s (h m+ (x;s s[m] q D q,s (h m (x;s. Next, we have xd (s x (h m(x;s xd q,s (h m (x;s D q,s (xh m (x;s D q,s (h m+ (x;s+s[m]h m (x;s Since D q,s (s, applying (4.2 we get D q,s (h m+ (x;s [m] q D q,s (sh m (x;s. xd (s x (h m(x;s D q,s (h m+ (x;s s[m] q D q,s (h m (x;s [m] q h m (x;qs. To end the proof, we subtract this from (4.6. Using martingale property (2.5 and Lemma 4.7, we rewrite (4.5 as A(h m (x;s [m] q h m (y;sp qs,s (x,dy hm (z;s h m (y;s ν s,y (dzp qs,s (x,dy. z y So by linearity, for any polynomial p in variable x we have p(z p(y (Ap(x µ x,s (dy,dz. z y Since(4.4givesD x (s (x n A(x n +xd x (s (x n andd x (s (x D q,s (h (x;s this determines D x (s on monomials: n D x (s (xn x k A(x n k k ( n k By the geometric sum formula, the integrand is n x kzn k y n k z y k x kzn k y n k µ x,s (dy,dz. z y (x zyn +(y xz n +(z yx n. (z x(x y(y z This shows that D (s x (s x on polynomials. With this identification, (4.3 concludes the proof of Theorem 4.6. Remark 4.9. Operators (s x and (s x are well defined on functions with bounded second derivatives in variable x, so we expect that (.2 holds in more generality.

13 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 3 For an analytic function f(x one can use contour integration to verify that each term in (.2 can be approximated uniformly on compacts by the same expression applied to a polynomial. Thus formula (.2 for polynomials yields f(b t f(+ ( (s x f(b s,sdlb + ( (s x f(b qs,sd q s. Proof of Remark 4.9. Using Cauchy formula with R large enough, we have (4.7 g(x,t : ( (t x f(x,t f(ζ 2πi (ζ x(ζ y dζdν x,t(dy. ζ R With p n (x n k a kx k f(x, we see that g(x,t p n (ζ 2πi (ζ x(ζ y dζdν x,t(dy ζ R + 2πi ζ R f(ζ p n (ζ (ζ x(ζ y dζdν x,t(dy can be approximated uniformly in x from the support of γ t;q by polynomials. In particular, expanding p n (x into the q-hermite polynomials, from Lemma 4.7 we deduce that g(x,t k b n (t [n]! (t x (h n(x,t k b n+ (t h n (x;t [n]! and the convergence is uniform over x from the support of γ t. The estimates from proof of Proposition 4.5 show that ( (s x f(b s,sdlb is well defined. Next we note that ( (s x f(x,t f(ζ 2πi (ζ x(ζ y(ζ z dζdµ x,s(dy,dz, ζ R so ( (s x f(x,t is also a limit of ( (s x p n (x,t, and the convergence is uniform over x from the support of γ t and over all s [,t]. By Theorem 4.6, formula (.2 is valid for polynomials p n (x and the contour integral identities show that we can pass to the limit as n. 5. A q-analogue of the stochastic exponential In this section we consider a more general function f(x,t for which the meansquare limit (4.8 is applicable. Our goal is to exhibit a solution to the differential equation dlz azdlb which we interpret in integral form as (5. Z t c+a Z s dlb with deterministic parameters a,c R. In view of the fact that we can integrate only instantaneous functions of B t, we seek a solution in this form, and we seek the series expansion. The solution parallels the development in Ito theory and relies on the identity (4. which is a special case of our definition (4.4. Thus, as

14 4 W LODEK BRYC a mean-square expansion the solution of (5. is (5.2 Z t c n a n h n (B t ;t [n] q! c k ( ( qaq k B t +a 2 tq 2k. Formal calculations give c+a Z s dlb c+ca n a n h n+ (B t ;t Z t, [n+] q! as h (x;t. Note that the product expression on the right hand side of (5.2 is well defined for all t, as with probability one B t < 2 t/ q. However, the solution of (5. lives only on the finite interval t < a 2 ( q. This is because for < q < the product j ( qj converges, so the series n a 2n [n]! s n d q s n a 2n t n+ [n+]! n ( q n a 2n t n+ n+ j2 ( qj converges when ( qa 2 t < and hence Corollary 4.4 can be applied only to this range of t. This is the same range of t where (Z t is a martingale [2, Corollary 4]. 6. Deterministic integrands In the Ito-Wiener theory, the distribution of the integral of a deterministic function b(s is Gaussian with the variance b2 (sds. In this section we show that for < q < the distribution of the random variable b(sdlb depends on b(s and q in a more complicated way. This gives rise to a plethora of q- Gaussian laws that all converge to the Gaussian law as q. To illustrate this phenomenon, we consider Z s r dlb q kr( B B q k q k+ k for r. Then E(Z and from (. we see that (6. E(Z 2 s 2r d q s q q 2r+ [2r+] q. To show that the distribution of Z depends on r and q in a complicated way, we compute the fourth moment. Recall that the 4-th moment of γ t;q is t 2 (2 + q, where t is the variance. We will show that this pattern fails for Z.

15 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 5 Noting that by martingale property E(B u B t B s for s < t < u, to compute E(Z 4 we need only the formulas ( (6.2 E (B t B s 4 (t s((q +2t 3qs for t < s and ( (6.3 E (B t 2 B t 2 (B u 2 B u 2 (u 2 u (t 2 t ( (6.4 E (B t 2 B t (B u 2 B u 3 ( q(u 2 u (t 2 t for t < t 2 u < u 2. The calculations use the first four martingale polynomials: the first three were ( already listed in Example 4.2 and the fourth one is h 4 (x;t x 4 q 2 +2q +3 tx 2 + ( q 2 +q + t 2. We use them to compute recurrently conditional moments. This gives martingale property that we already used, the quadratic martingale property E((B t 2 B s (B s 2 +t s, and a bit more cumbersome formulas (6.5 E((B t 3 B s (B s 3 +(t s(2+qb s (6.6 E((B t 4 B s (B s 4 +(t s(3+2q+q 2 (B s 2 +(t s((2+qt (+q+q 2 s After a calculation these formulas give E ( (B t B s 4 B s (t s (( q 2 (B s 2 +(q +2t ( q 2 +q + s. To derive (6.2 we take the expected value of this expression. To derive(6.3, we use martingale property and quadratic martingale property to compute E ( (B u 2 B u 2 ( B t 2 E (B u 2 B t 2 2 ( B t 2 E (B u B t 2 2 B t 2 u2 u. To derive (6.4 we use again the formulas for conditional moments to compute E ( (B u 3 B u ( q(u2 u B u. This gives E u 2 B ( (B t 2 B t (B u 2 B u 3 ( q(u 2 u ( E(B t 2 B u E(B t B u.

16 6 W LODEK BRYC We now compute E(Z 4. Denoting ξ k B B, we see that q k q k+ k k E(Z 4 q 4rk Eξk 4 +6 q 2r(j+k E(ξjξ 2 k+4 2 q r(3j+k E(ξ k ξj, 3 k k j k j as all other terms vanish. From (6.2 we see that E((B B 4 ( tq k tq k+ q 2 t 2 q 2k (2+3q, so q 4rk Eξk 4 ( q2 (2+3q q 4r+2. k From (6.3 we get k q 2r(j+k E(ξj 2 ξ2 k k ( q2 q (2r+(j+k k j k j ( q 2 q 2r+ ( q 2r+ 2 (+q 2r+. From (6.4 we get k k q r(3j+k E(ξ k ξj 3 ( q 3 q (3r+j+(r+k ( q 3 q r+ ( q r+ ( q 4r+2. k j This gives j E(Z 4 ( q2( 2+3q 6q r+ +q r+2 +4q 2r+ 3q 2r+2 q 3r+3 ( q r+ ( q 2r+ 2. (+q 2r+ So using (6. we get E(Z 4 (E(Z q 6qr+ +q r+2 +4q 2r+ 3q 2r+2 q 3r+3 ( q r+ (+q 2r+. Acknowledgements. I would like to thank Magda Peligrad for positive feedback and encouragement. Referee s comments helped to improve the presentation. References [] Michael Anshelevich. Generators of some non-commutative stochastic processes. Probability Theory and Related Fields, 57(3-4:777 85, 23. [2] Philippe Biane. Processes with free increments. Math. Z., 227(:43 74, 998. [3] Marek Bożejko, Burkhard Kümmerer, and Roland Speicher. q-gaussian processes: noncommutative and classical aspects. Communications in Mathematical Physics, 85(:29 54, 997. [4] Marek Bożejko and Roland Speicher. An example of a generalized Brownian motion. Comm. Math. Phys., 37(3:59 53, 99. [5] Wlodek Bryc and Jacek Weso lowski. Infinitesimal generators of q-meixner processes. Stoch. Proc. Appl., 24:95 926, 24. [6] W lodzimierz Bryc, Pawe l Szab lowski, and Wojciech Matysiak. Probabilistic aspects of Al-Salam Chihara polynomials. Proc. Amer. Math. Soc., 33:27 34, 25. arxiv.org/abs/math.ca/3455. [7] W lodzimierz Bryc and Jacek Weso lowski. Conditional moments of q-meixner processes. Probab. Theory Related Fields, 3:45 44, 25. arxiv.org/abs/math.pr/436. [8] Rafael Díaz and Eddy Pariguan. On the Gaussian q-distribution. Journal of Mathematical Analysis and Applications, 358(: 9, 29.

17 ON INTEGRATION WITH RESPECT TO THE q-brownian MOTION 7 [9] Catherine Donati-Martin. Stochastic integration with respect to q Brownian motion. Probability theory and related fields, 25(:77 95, 23. [] U Frisch and R Bourret. Parastochastics. Journal of Mathematical Physics, (2:364 39, 97. [] Emmanuel Haven. Quantum calculus (q-calculus and option pricing: a brief introduction. In Sofge D. Lawless W. van Rijsbergen K. Bruza, P., editor, Quantum Interaction, volume 5494 of Lect. Notes in Artificial Intelligence, pages Springer, 29. [2] Emmanuel Haven. Itô s lemma with quantum calculus (q-calculus: Some implications. Foundations of Physics, 4(3: , 2. [3] Mourad E. H. Ismail. Classical and quantum orthogonal polynomials in one variable, volume 98 of Encyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, 25. With two chapters by Walter Van Assche, With a foreword by Richard A. Askey. [4] Mourad EH Ismail and Dennis Stanton. On the Askey-Wilson and Rogers polynomials. Canad. J. Math, 4(5:25 45, 988. [5] F. H. Jackson. On q-definite integrals. Quarterly Journal of Pure Applied Mathematics, :93 23, 9. [6] Victor G Kac and Pokman Cheung. Quantum calculus. Springer, 22. [7] Roelof Koekoek and Rene F. Swarttouw. The Askey scheme of hypergeometric orthogonal polynomials and its q-analogue, 998. Delft University of Technology Report no. 98-7, koekoek/askey.html. [8] Wim Schoutens. Lévy-Sheffer and IID-Sheffer polynomials with applications to stochastic integrals. Journal of Computational and Applied Mathematics, 99: , 998. [9] Wim Schoutens. Stochastic processes and orthogonal polynomials, volume 46 of Lecture Notes in Statistics. Springer-Verlag, New York, 2. [2] Pawel J Szab lowski. q-wiener and (α,q-ornstein Uhlenbeck processes. a generalization of known processes. Theory of Probability & Its Applications, 56(4: , 22. Department of Mathematical Sciences, University of Cincinnati, PO Box 225, Cincinnati, OH , USA address: Wlodzimierz.Bryc@UC.edu