arxiv: v2 [math.pr] 19 Nov 2015

Transcription

1 INFINITESIMAL GENERATORS FOR A CLASS OF POLYNOMIAL PROCESSES arxiv:4.6268v2 [math.pr] 9 Nov 205 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI Abstract. We study the infinitesimal generators of evolutions of linear mappings on the space of polynomials, which correspond to a special class of Markov processes with polynomial regressions called quadratic harnesses. We relate the infinitesimal generator to the unique solution of a certain commutation equation, and we use the commutation equation to find an explicit formula for the infinitesimal generator of free quadratic harnesses. This is an expanded (arxiv version of the paper.. Introduction In this paper we study properties of evolutions of degree-preserving linear mapping on the space of polynomials. We analyze these mappings in a self-contained algebraic language assuming a number of algebraic properties which were abstracted out from some special properties of a family of Markov processes called quadratic harnesses. We therefore begin with a review of the relevant theory of Markov processes that motivates our assumptions. Markov processes that motivate our theory have polynomial regressions with respect to the past σ-fields. Such Markov processes enjoy many interesting properties and appeared in numerous references [2, 3, 0, 2, 5, 6, 7, 22, 23, 24, 25, 26, 27, 28, 29, 30, 3]. As observed by Cuchiero in [4] and Szab lowski [32], transition probabilities P s,t (x,dy of such a process on an infinite state space define the family of linear transformations (P s,t 0 s t that map the linear space P = P(R of all polynomials in variable x into itself. The crucial property that holds in many interesting examples is that transformations P s,t do not increase the degree of a polynomial, that is P s,t : P k P k, where P k denotes a linear space of polynomials of degree k, k = 0,,... This will be the basis for our algebraic approach. Cuchiero in [4], see also [5] introduced the term polynomial process to denote such a process in the time-homogeneous case; we will use this term more broadly to denote the family of operators rather than a Markov process. That is, we adopt the point of view that the linear mappings P s,t of P can be analyzed in abstract without explicit reference to the underlying Markov process and the transition operators. Date: Printed August, 208. (PF27.tex Mathematics Subject Classification. 60J25;46L53. Key words and phrases. Infinitesimal generators; quadratic conditional variances; polynomial processes.

2 2 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI We note that operators P s,t are well-defined whenever the support of X s is infinite. So, strictly speaking, operators P 0,t are well defined only for the so called Markovfamilies that can be started at an infinite number of values of X 0. However, it will turn out that in the cases that we are interested in, even if a Markov process starts with X 0 = 0 we can pass to the limit in P s,t as s 0 and in this way define a unique degree-preserving mapping P 0,t : P P. With the above in mind, we introduce the following definition. Definition.. Suppose P is the linear space of polynomials in one variable x. A polynomial process is a family of linear maps {P s,t : P P, 0 s t} with the following properties: (i for k = 0,,... and 0 s t, (ii P s,t ( =, (iii for 0 s t u P s,t (P k = P k, (. P s,t P t,u = P s,u. Our next task is to abstract out properties of polynomial processes that correspond to a special class of such Markov processes with linear regressions and quadratic conditional variances under the two-sided conditioning. The two-sided linearity of regression, which is sometimes called the harness property, takes the following form: (.2 E(X t X s,x u = u t u s X s + t s u s X u, 0 s < t < u, see e. g., [9, formula (2]. We will also assume that the conditional second moment of X t given X s, X u is a second degree polynomial in the variables X s and X u. The conclusion of [8, Theorem 2.2] says that there are five numerical constants η,θ R, σ,τ 0, and γ +2 στ such that for all 0 s < t < u, (.3 Var(X t X s,x u = (u t(t s u(+σs+τ γs ( +η ux s sx u +θ X u X s u s u s +σ (ux s sx u 2 (u s 2 +τ (X u X s 2 (u s 2 ( γ (X u X s (ux s sx u (u s 2 Such processes, when standardized, are called quadratic harnesses. Typically they are uniquely determined by five constants η, θ, σ, τ, γ from (.3. We will also assume that (X t is a martingale in a natural filtration. This is a consequence of (.2 when E(X t does not depend on t, compare [0, page 47]. The martingale property is easily expressed in the language of polynomial processes, asit just saysthat P s,t (x = x. Somewhat moregenerally, ifthe polynomials m k (x;t in variablex, k 0, aremartingalepolynomialsfor amarkovprocess(x t, that is E(m k (X t ;t X s = m k (X s,s, k 0, then P s,t (m k ( ;t(x = m k (x;s for s < t. It is harder to deduce the properties of P s,t that correspond to (.2 and (.3. We will find it useful to describe the linear mappings of polynomials in the algebraic language, and we will rely on martingale polynomials to complete this task..

3 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 3.. Algebra Q of sequences of polynomials. We consider the linear space Q of all infinite sequences of polynomials in variable x with a nonstandard multiplication defined as follows. For P = (p 0,p,...,p k,... and Q = (q 0,q,...,q k,... in Q, we define their product R = PQ as the sequence of polynomials R = (r 0,r,... Q given by (.4 r k (x = deg(q k j=0 [q k ] j p j (x,k = 0,,... It is easy to check that algebra Q has identity E = (,x,x 2,... In the sequel we will frequently use two special elements of Q (.5 F = (x,x 2,x 3,..., and (.6 D = (0,,x,x 2,... which is the left-inverse of F. It may be useful to note that PF and PD act as shifts: if P = (p 0 (x,p (x,...,p n (x,... then PF = (p (x,p 2 (x,...,p n+ (x,... and PD = (0,p 0 (x,p (x,...,p n (x,... On the other hand, FP is just a multiplication by x, i.e., FP = (xp 0 (x,xp (x,...,xp n (x,... Algebra Q is isomorphic to the algebra of all linear mappings P P under composition: to each P : P P we associate a unique sequence of polynomials P = (p 0,p,...,p k,... in variable x by setting p k = P(x k. Of course, if P does not increase the degree, then p k is of at most degree k. Under this isomorphism, the composition R = P Q of linear operators P and Q on P induces the multiplication operation R = PQ for the corresponding sequences of polynomials P = (p 0,p,...,p k,... and Q = (q 0,q,...,q k,... which was introduced in (.4. It is clear that degree-preserving linear mappings of P are invertible. Proposition.2. If for every n polynomial p n is of degree n then P = (p 0,p,... has multiplicative inverse Q = (q 0,q,... and each polynomial q n is of degree n.

4 4 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI Proof. Write p n (x = n k=0 a n,kx k. The inverse Q = (q 0,q,... is given by the family of polynomials q n of degree n which solve the following recursion: q 0 =, and q n = n x n a n,j q j (x for n. a 0,0 a n,n It is then clear that QP = E. Since both P and Q consist of polynomials of degree n at the n-th position of the sequence, the corresponding linear mappings P, Q on the set of polynomials preserve the degree of a polynomial. Since Q P is the identity on each finitedimensional space P n, we see that P Q is also an identity, so PQ = E. j=0 Inthispaperwestudy P s,t throughthecorrespondingelementsp s,t ofthealgebra Q. We therefore rewrite Definition. in the language of algebra Q. Definition.3 (Equivalent form of Definition.. A polynomial process is a family {P s,t Q : 0 s t} with the following properties: (i for 0 s t and n = 0,,..., the n-th component of P s,t is a polynomial of degree n, (ii P s,t (E FD = E FD, (iii for 0 s t u we assume (.7 P s,t P t,u = P s,u. Since special elements (.5 and (.6 satisfy (,0,0,... = E FD, property (ii is P s,t (,0,0,... = (,0,0,..., see Definition.(ii..2. Martingale polynomials. In this section we re-derive [32, Theorem ]. We note that by Proposition.2 each P s,t is an invertible element of Q. So from (.7 we see that P t,t = E for all t 0. In particular, P 0,t is invertible and M t = P 0,t consists of polynomials m k (x;t in variable x of degree k. From (.7, we get P 0,s P s,t M t = P 0,t M t = E. Multiplying this on the left by M s = P 0,s, we see that (.8 M s = P s,t M t, i.e., M t is a sequence of martingale polynomials for {P s,t } and in addition M 0 = E. Conversely, (.9 P s,t = M s M t. Ref. [32] points out that in general martingale polynomials are not unique. However, any sequence M t = (m 0 (x;t,m (x;t,... of martingale polynomials (with m k (x;t of degree k for k = 0,,..., still determines uniquely P s,t via P s,t = M s M t..3. Quadratic harnesses. Recall that our goal is to abstract out the properties of the (algebraic polynomial process {P s,t } that correspond to relations (.2 and (.3. Since this correspondence is not direct, we first give the self-contained algebraic definition, and then point out the motivation for such definition.

5 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 5 Definition.4. We will say that a polynomial process {P s,t : 0 s t} is a quadratic harness with parameters η, θ, σ, τ, γ if the following three conditions hold: (i (martingale property P s,t (FD F 2 D 2 = FD F 2 D 2, (ii (harness property there exists X Q such that for all t 0 we have (.0 P 0,t F = (F+tXP 0,t, (iii (quadratic harness property element X Q introduced in (.0 satisfies the following quadratic equation (. XF γfx = E+ηF+θX+σF 2 +τx 2. Remark.5. Using martingale polynomials M t = P 0,t re-write assumption (.0 as (.2 FM t = M t (F+tX. from Section.2, we can We now give a brief explanation how the algebraic properties in conditions (i- (iii of Definition.4 are related to the corresponding properties of a (polynomial Markov process..3.. Motivation for the martingale property. For a Markov process (X t with the natural filtration (F t this takes a more familiar form E(X t F s = X s. Translated back into the language of linear operators on polynomials, this is P s,t (x = x. In the language of algebra Q this is P s,t (0,x,0,0,... = (0,x,0,0,... To get the final form of condition (i we note that FD F 2 D 2 = (0,x,0,0, Motivation for harness property (.0. For a Markov processes (X t with martingale polynomials m n (x;t, property (.2 implies that E(X t m n (X t ;t X s = E(X t m n (X u ;u X s = E(E(X t X s,x u m n (X u ;u X s = u t u s X se(m n (X u ;u X s + t s u s E(X um n (X u ;u X s = u t u s X sm n (X s ;s+ t s u s E(X um n (X u ;u X s The resulting identity E(X t m n (X t ;t X s = u t u s X sm n (X s ;s+ t s u s E(X um n (X u ;u X s. in the language of algebra Q with M t = (m 0 (x;t,m (x;t,... becomes (.3 P s,t FM t = u t u s FM s + t s u s P s,ufm u. Sinceeachpolynomialxm n (x;tcanbewrittenasthelinearcombinationofm 0 (x;t,..., m n+ (x;t, onecanfindj t QsuchthatFM t = M t J t. Insertingthisinto(.3, wesee that martingalepropertyeliminates P s,t and after left-multiplication by M s we get (u sj t = (u tj s +(t sj u. In particular, J t depends linearly on t and thus J t = J 0 +t(j J 0.

6 6 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI This shows that for any martingale polynomials harness property implies that there exist Y,X Q such that (.4 FM t = M t (Y+tX. In our special case of M t = P 0,t, we get Y = J 0 = FM 0 = FE = F. This establishes (.2, which of course is equivalent to ( Motivation for quadratic harness property (.. Suppose polynomial process {P s,t } in the sense of Definition.3 arises from a Markov process with polynomial conditional moments which is a harnesses and in addition has quadratic conditional variances (.3. Then from the previous discussion, (.4 holds, and under mild technical assumptions, [8, Theorem 2.3] shows that X, Y satisfy commutation equation which reduces to (. when Y = F. This motivates condition (iii. In fact, it is known that under some additional assumptions on the growth of moments, (.0 and (. imply (.2 and (.3, see [3, Section 4.]; this equivalence is also implicit in the proof of [8, Theorem 2.3] and is explicitly used in [2, page 244]. However, this has no direct bearing on our paper, as in this paper we simply adopt the algebraic Definition Infinitesimal generator. A polynomialprocess{p s,t }with harnessproperty (.0 is uniquely determined by X. Indeed, the n-th element of the sequence on the left hand side of (.0 is the (n + -th polynomial in P 0,t while the n-th element of the sequence on the right hand side of (.0 depends only on the first n polynomials in P 0,t. In fact, one can check that (.5 P 0,t = (F+tX k (E FDD k k=0 Proof. From (.0 we get (.6 P 0,t F(F k D k F k+ D k+ = (F+tXP 0,t (F k D k F k+ D k+ Denote by Y k = P 0,t (F k D k F k+ D k+. Multiplying (.6 by D from the right, we get Y k+ = (F+tXY k D. Since Y 0 = P 0,t (E FD = E FD we get Y k = (F + tx k (E FDD k. The telescoping series gives P 0,t = k=0 Y k. Since the k-th element of F+tX has degree k+, it follows from (.2 that M t is a rational function of t. Formula (.9 shows that the left infinitesimal generator (.7 A t = lim h 0 + h (P t h,t E, t > 0 is a well defined element of Q, and that (.8 A t M t = t M t with differentiation defined componentwise.

7 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 7 Since P s,t is continuous in t, the right infinitesimal generator exists and is given by the same expression. To see this, we compute lim h 0 + h (P t,t+h E on M t. We have lim h 0 + h (P t,t+hm t M t = lim P t,t+h h 0 + h (M t M t+h = lim h 0 +P t,t+h lim h 0 + h (M t M t+h = t M t. The infinitesimal generator A t determines P s,t uniquely; for an algebraic proof see Proposition 2.3. Clearly, M t = E t 0 A sm s ds. The infinitesimal generator A t and its companion operator A t : P P are the main objects of interest in this paper. Remark.6. We note that A t = (0,a (x;t,a 2 (x;t,... always starts with a 0, as from P s,t ( = it follows that A t ( = 0. It is clear that a n (x;t is a polynomial in x of degree at most n. Martingale property implies that a (x;t = 0, as A t (x = 0. The infinitesimal generator considered as an element of Q in the latter case starts with two zeros, A t = (0,0,a 2 (x;t,a 3 (x;t, Basic properties of infinitesimal generators for quadratic harnesses Our first result introduces an auxiliary element H t Q that is related to A t by a commutation equation. Theorem 2.. Suppose that {P s,t Q : 0 s t} is a quadratic harness as in Definition.4 with generator A t. For t > 0, let (2. H t = A t F FA t and denote T t = F th t. Then (2.2 H t T t γt t H t = E+θH t +ηt t +τh 2 t +σt2 t. Proof. Differentiating (.2 and then using (.8 we get Multiplying this from the right by M t by F we get FA t M t = A t M t (F+tX+M t X. FA t = A t F+M t XM t. Comparing this with (2. we see that (2.3 H t = M t XM t. and using (.2 to replace M t (F+tXM t We nowuse this and(.2 in the definition oft t = F th t = M t (F+tXM t th t. Using (2.3 we get T t = M t FM t. To derive equation (2.2 we now multiply (. by M t from the left and by M t from the right. Remark 2.2. As observed in Remark.6 the n-th element of A t is a polynomial of degree at most n. Thus writing H t = (h 0 (x,h (x,..., from (2. we see that h n is of degree at most n+. Martingale property implies that h 0 (x = 0.

8 8 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI We will need several uniqueness results that follow from the more detailed elementwise analysis of the sequences of polynomials. We first show that the generator determines uniquely the polynomial process, at least under the harness (.0 condition. Proposition 2.3. A polynomial process {P s,t : 0 s < t} with harness property (.0 is determined uniquely by its generator A t. Proof. Fix s < t. Combining (.9 with (.2 we get P s,t F FP s,t = M s Mt F FM s M t = M s (Y+tXMt M s (Y+sXM t = (t sm s M t Therefore, from (2.3 we get With s = 0, this implies P s,t F = FP s,t +(t sp s,t H t. M t F = (F th t M t. ( Mt XM t. This equation is similar to equation (.0 and again uniqueness follows from consideration of the degrees of the polynomials. The solution is M t = (F th t k (E FDD k, k=0 compare(.5. Thus M t is determined uniquely, and (.9 shows that operators P s,t are uniquely determined. Since H t is expressed in terms of A t by (2., this ends the proof. Next, weshowthath t = (h 0,h,...isuniquelydeterminedbythecommutation equation(2.2withthe initialcondition h 0 = 0. FromtheproofofProposition2.3 it therefore follows that the entire quadratic harness {P s,t } as well as its generator are also uniquely determined by (2.2. Proposition 2.4. If σ,τ 0, and στ then equation (2.2 has a unique solution among H t Q such that h 0 (x = 0. Proof. Eliminating T t = F th t from (2.2 we can rewrite it into the following equivalent form. (2.4 H t F γfh t = E+θH t +η(f th t +τh 2 t +σ(f th t 2 +( γth 2 t. Write H t = H = (h n (x n=0,,... for a fixed t > 0, with h 0 = 0. We will simultaneously prove that for n, polynomial h n (x is uniquely determined and that its degree is at most n+. The proof is by induction. With h 0 = 0 it is clear that the thesis holds true for n = 0. We therefore assume that h 0,h,...,h n are given polynomials of degrees at most,2,...,n+, respectively.

9 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 9 Looking at the n-th element of (2.4 for n 0, we get the following equation for h n+ (x: ( (2.5 +σt [hn ] n+ (σt 2 +( γt+τ h n+ (x = x n +ηx n+ +σx n+2 +(θ tηh n (x +(γ σtxh n (x+(σt 2 +( γt+τ n [h n ] j h j (x. The degree of the polynomial on the right hand side is at most n + 2, as the highest degree term on the right hand side is (σ +(γ σt[h n ] n+ x n+2. In order to complete the proof, we only need to verify that for all n 0, the coefficient +σt [h n ] n+ (σt 2 +( γt+τ on the left hand side of (2.5 does not vanish. We need to consider separately two cases. Case γ +στ 0: We proceed by contradiction. Suppose that for some n 0 the coefficient at h n+ (x on the left hand side of (2.5 is 0. Since +σt > 0, this implies that σt 2 +( γt+τ cannot be zero. So we get [h n ] n+ = +σt σt 2 +( γt+τ. We now use this value to compute the coefficient at x n+2 on the right hand side of (2.5. We get σ +(γ σt[h n ] n+ = γ +στ σt 2 +( γt+τ 0. Since the left hand side of (2.5 is 0, and the degree of the right hand side of (2.5 is n+2, this is a contradiction. Thisshowsthatthecoefficientofh n+ (x onthelefthandsideof (2.5is non-zero for all n 0. So each polynomial h n+ (x is determined uniquely and has degree at most n+2 for all n 0. Case γ +στ = 0: In this case σt 2 +( γt+τ = (+σt(τ+t. Comparing the coefficients at x n+2 on both sides of (2.5 we get (2.6 (+σt( [h n ] n+ (τ +t[h n+ ] n+2 = σ( [h n ] n+ (τ +t. Since h 0 (x = 0, this gives [h ] 2 = σ/(+σt. So for n = we get [h n ] n+ (τ+t = ( στ/(+σt 0. Dividing both sides of (2.6 by this expression we get recursively[h n ] n+ = σ/(+σt for all n. Therefore, for n, the left hand side of (2.5 simplifies to ( στh n+ (x. Using again στ 0, this shows that polynomial h n+ is determined uniquely and its degree is at most n+2. Of course, (2.5 determines h (x uniquely, too, as h 0 (x = 0. Our main result is the identification of the infinitesimal generator for quadratic harnesses with γ = στ. Such processes were called free quadratic harnesses in [8, Section 4.]. Markov processes with free quadratic harness property were constructed in [9] but the construction required more restrictions on the parameters than what we impose here. In this section we represent the infinitesimal generators using the auxiliary power series ϕ t (D with D defined by (.6. In Section 4 we will use the results of this j=0

10 0 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI section to derive the integral representation (4.3 for the operator A t under a more restricted range of parameters η,θ,σ,τ. For a formal power series ϕ(z = k=0 c kz k we shall write ϕ(d for the series k c kd k. We note that since the sequence D k begins with k zeros, k c kd k is a sequence of finite sums: n ϕ(d = (c 0,c 0 x+c,c 0 x 2 +c x+c 2,..., c j x n j,... So ϕ(d is a well defined element of Q. We will also need D = k=0 Fk D k+ = (0,,2x,3x 2,... which represents the derivative. Theorem 2.5. Fix η,θ R and σ,τ 0 such that στ. Then the infinitesimal generator of the quadratic harness with the above parameters and with γ = στ is given by (2.7 A t = +σt (E+ηF+σF2 D ϕ t (DD, t > 0, where ϕ t (z = k= c k(tz k for small enough z solves the quadratic equation (2.8 (z 2 +ηz +σ(t+τϕ 2 t +((θ tηz 2tσ στ ϕ t +tσ+ = 0 and the solution is chosen so that ϕ t (0 =. (For στ < this solution is written explicitly in formula (4.4 below. We note that for each fixed t equation (2.8 has two real roots for z close enough to 0. As ϕ t (z we choose the smaller root when στ < and the larger root when στ >. For z = 0, equation (2.8 becomes j=0 σ(t+τϕ 2 t (0 (+στ +2tσϕ t(0++tσ = 0, so this procedure ensures that ϕ t (0 =. 3. Proof of Theorem 2.5 The plan of proof is to solve equation (2.2 for H t, and then to use equation (2. to determine A t. 3.. Part I of proof: solution of equation (2.2 when γ = στ. Equation (2.2 takes the form H t F th 2 t +στfh t στth 2 t = E+θH t +η(f th t +τh 2 t +σ(f2 th t F tfh t +t 2 H 2 t. So after simplifications, the equation to solve for the unknown H t is (3. (+σth t F = E+ηF+σF 2 +(θ ηth t σ(t+τfh t +(t+τ(+σth 2 t. Lemma 3.. The solution of (3. with the initial element h 0 = 0 is (3.2 H t = +σt (E+ηF+σF2 ϕ t (DD, where ϕ t satisfies equation (2.8 and ϕ t (0 =.

11 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES Proof. Since t > 0 is fixed, we suppress the dependence on t and we use Remark 2.2 to write H t = H = (0,h (x,... From (3. we read out that h (x = +σt (+ηx+σx2. From Proposition 2.4 we see that (3. has a unique solution. In view of uniqueness, we seek the solution in a special form (3.3 H = +σt (E+ηF+σF2 with c = and c k = c k (t R. Note that H 2 = (+σt 2(E+ηF+σF2 + = k= c k D k k= c k D k (E+ηF+σF 2 c j D j (+σt 2(E+ηF+σF2 η (+σt 2(E+ηF+σF2 j= j= c k D k c j D j k= j= c k D k c j D j σ + (+σt 2(F+ηF2 +σf 3 c j D j σ + (+σt 2(E+ηF+σF2 c k D k 2 c k D j. Inserting this into (3. we get k= j= k=2 j= E+ηF+σF 2 + ( E+ηF+σF 2 c k D k = E+ηF+σF 2 + θ tη +σt (E+ηF+σF2 + t+τ [ (E+ηF+σF 2 +σt k= k= k=2 c k D k σ(t+τ +σt (F+ηF2 +σf 3 j= c k D k k= c k D k c j D j +η(e+ηf+σf 2 c k D k c j D j +σ(f+ηf 2 +σf 3 c j D j +σ(e+ηf+σf 2 c k D k 2 c j D j]. j= The terms with ( F+ηF 2 +σf 3 cancel out, so ( E+ηF+σF 2 factors out. We further restrict our search for the solution by requiring that the remaining factors match, i.e. k= c k+ D k = θ tη +σ k= j= c k D k k= + t+τ ( c k D k c j D j +η c k D k c j D j +σ c k D k 2 c j D j. +σt k= j= k=2 k=2 k= j= j= j=

12 2 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI Collecting the coefficients at the powers of D we get c k+ D k = θ tη c k D k +σt k= + t+τ ( +σt k=2 j= k= k k k c j c k j D k +η c j+ c k j D k +σ c j+2 c k j D k. k= j=0 k= j=0 We now compare the coefficients at the powers of D. Since στ, for k = we get c 2 = θ ηt +σt + t+τ +σt η + t+τ +σt σc 2. So c 2 = (θ +ητ/( στ = β (say. For k 2, we have the recurrence (3.4 c k+ = θ ηt +σt c k + t+τ (k +σt j= k c j c k j +η j=0 k c j+ c k j +σ j=0 c j+2 c k j. We solve this recurrence by the method of generating functions. One can proceed herewith a formalpowerseries, andthen invokeuniqueness toverifythat the power series has positive radius of convergence. Or one can use an a priori bound from Lemma 3.2 below and restrict the argument of the generating function to small enough z. Note that in this argument, t is fixed. Let ϕ(z = k= c kz k. Then ϕ(z = +βz + c k+ z k k=2 = θ ηt t+τ z(ϕ(z + +σt +σt z2 ϕ 2 (z+ηz t+τ +σt (ϕ2 (z +σ t+τ +σt (ϕ(z(ϕ(z βz. This gives quadratic equation (2.8 for ϕ = ϕ t. The following technical lemma assures that the series ϕ(z = k= c kz k converges for all small enough z. Lemma 3.2. Suppose στ and {c k } is a solution of recursion (3.4 with c = for a fixed t. Then for every p > there exists a constant M such that (3.5 c k Mk 2 k p for all k 3. Proof. We will prove the case p = 2 only, as this suffices to justify the convergence of the series. Solving (3.4 for c k+, which appears also in one term on the right hand side of (3.4, we get στ +σt c θ ηt k+ +σt c k + t+τ (k +σt j= k c j c k j + η j=0 k 2 c j+ c k j +σ Since we are not going to keep track of the constants, we simplify this as j=0 c j+2 c k j

13 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 3 c k+ θ ηt στ c k + t+τ (k στ j= k c j c k j + η j=0 k 2 c j+ c k j +σ j=0 c j+2 c k j (n (3.6 c n+ A c n +B k= n c k c n k + k=0 n 2 c k+ c n k + k=0 c k+2 c n k with A = θ ηt στ t+τ and B = στ (+ η +σ. Here, n 2. We now choose M large enough so that (3.5 holds for k = 3,4,5,6, and we also require that (3.7 4A+(24 c 2 +75B M. We now proceed by induction and assume that (3.5 holds for all indices k between 3 and n for some n 6. To complete the induction step we provide bounds for the sums on the right hand side of (3.6. The first sum is handled as follows: n n 3 c k c n k = 2 c c n +2 c 2 c n 2 + c k c n k. k= k=3 We apply the induction bound (3.5 to c 3,...,c n. Noting that c = we get n c k c n k 2 Mn 3 (n 2 +2 c 2 k= 8 Mn 2 (n c 2 M n 2 (n+ 2 +Mn 2 M n 2 M n 4 (n 2 2 +Mn 4 [n/2] k=3 8 Mn 2 (n c 2 (n Mn 2 n 2 n 3 k=3 k 2 (n k 2 k 2 (n k 2 +Mn 2 [n/2] k=3 ++4Mn 2 k 2 n 2 n 3 k=[n/2]+ n 2 k=[n/2]+ k 2 (n k 2 (n k 2 < Mn 2 (n+ 2 ( 8+8 c2 +6π 2 /3 < (6+8 c 2 Mn 2 (n+ 2. Here we used several times the bound (n + /(n 2 2 for n 6, inequality M, and we estimated two finite sums by the infinite series k= /k2 = π 2 /6.

14 4 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI We proceed similarly with the second sum: n n 3 c k+ c n k = 2 c c n +2 c 2 c n + c k+ c n k k=0 M n 3 2 Mn 2 n 2 +2 c 2 (n 2 + 4Mn 3 n 2 k=2 [n/2] k= k 2 + 4Mn 3 n 2 n 3 k=[n/2]+ (n k 2 Mn 2 (n+ 2 ( 8+8 c2 +6π 2 /3 Mn 2 (n+ 2 (6+8 c 2. Finally, we bound the third sum on the right hand side of (3.6. We get n 2 n 3 c k+2 c n k = 2 c 2 c n + c k+2 c n k k=0 2 c 2 Mn 2 n 2 k= + 4Mn 2 n 2 [n/2] k= k 2 + 4Mn 2 n 2 n 3 k=[n/2]+ (n k 2 Mn 2 (n+ 2 ( 8 c2 +6π 2 /3 Mn 2 (n+ 2 (8 c Combining these bounds together and using the induction assumption to the first term on the right hand side of (3.6 we get c n+ Mn 2 (n+ 2 (4A+(24 c 2 +75B In view of (3.7, this shows that c n+ Mn (n+, thus completing the proof of (3.5 2 by induction Part II of proof: solution of equation (2.. We now use (3.2 to determine A t. Since E FD = (,0,0,..., from Remark.6 it follows that A t FD = A t. Therefore, multiplying (2. by D from the right we get Iterating this, we get (3.8 A t = A t = FA t D+H t D. F k H t D k+, k=0 which is well defined as the series consists of finite sums elementwise. Since H t is given by (3.2, we get (3.9 A t = +σt F k (E+ηF+σF 2 ϕ t (DDD k+. k=0 = +σt (E+ηF+σF2 F k D k+ ϕ t (DD. k=0

15 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 5 We now note that (3.0 F k D k+ = D. k=0 (This can be seen either by examining each element of the sequence, or by solving the equation D F FD = E, which is just a product formula for the derivative, by the previous technique. The latter equation is of course of the same form as (2.. Replacing the series in (3.9 by the right hand side of (3.0 we get (2.7. This ends the proof of Theorem Integral representation In this section it is more convenient to use the linear operators on P instead of the sequences of polynomials. We assume that the polynomial process {P s,t : 0 s t} corresponds to a quadratic harness from Definition.4, and as before we use the parameters η,θ,σ,τ and γ to describe the quadratic harness. Infinitesimal generators of several quadratic harnesses, all different than those in Theorem 2.5, have been studied in this language by several authors. For quadratic harnesses with parameters η = σ = 0 and γ = q (,, according to [3], the infinitesimal generator A t acting on a polynomial f is (4. A t (f(x = R x ( f(y f(x y x ν x,t (dy, where ν x,t (dy is a uniquely determined probability measure. By inspecting the recurrences for the orthogonal polynomials {Q n } and {W n }, from [3, Theorem.(ii] one can read out that for q 2 t ( + qτ probability measure ν x,t can be expressed in terms of the transition probabilities P s,t (x,dy of the Markov process by the formula ν x,t (dy = P tq2 (+qτ,t(θ + qx,dy. In this form, the formula coincides with Anshelevich [, Corollary 22] who considered the case η = θ = τ = σ = 0 and γ = q [,]. (However, the domain of the generator in [] is much larger than the polynomials. Earlier results in Refs. [5, page 392], [6, Example 4.9], and [7] dealt with infinitesimal generators for quadratic harnesses such that σ = η = γ = 0. The following result gives an explicit formula for the infinitesimal generator of the evolution corresponding to the free quadratic harness in the integral form similar to (4.. The main new feature is the presence of an extra quadratic factor in front of the integral in expression (4.3 for the infinitesimal generator. Denote (4.2 α = η +θσ ητ +θ, β = στ στ. Theorem 4.. Fix σ,τ 0 such that στ < and η,θ R such that +αβ > 0. Let γ = στ and let ϕ t be a continuous solution of (2.8 with ϕ t (0 =. Then ϕ t is a moment generating function of a unique probability measure ν t and the generator of the quadratic harness with the above parameters on p P is (4.3 A t (p(x = +ηx+σx2 +σt x ( p(y p(x y x ν t (dy, t > 0.

16 6 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI Proof. Since στ <, the solution of quadratic equation (2.8 with ϕ t (0 = is z(tη θ+2tσ+στ + (4.4 ϕ t (z = 2(t+τ(z 2 +zη +σ (z(tη θ+2tσ+στ +2 4(z 2 +zη +σ(+tσ(t+τ 2(t+τ(z 2. +zη +σ (We omit the re-write for the case σ = η = 0. To avoid issues with σ = η = 0, it is better to rewrite (4.4 as ϕ t (z = 2(+tσ A+ A 2 4B. with A = z(tη θ+2tσ+στ + and B = (z 2 +zη+σ(+tσ(t+τ We will identify ν t through its Cauchy-Stieltjes transform G νt (z = z x ν t(dx which in our case will be well defined for all real z large enough. To this end we compute (+στ +2σtz +tη θ (4.5 ϕ t (/z/z = 2(t+τ(σz 2 +ηz + [( στz (α+σβt β ατ] 2 4(+σt(t+τ(+αβ 2(t+τ(σz 2. +ηz + Under our assumptions, the above expression is well defined for real large enough z R. Expression(4.5 coincides with the Cauchy-Stieltjes transform in[2, Proposition 2.3], with their parameters and c SY = στ +σt, α ητ +θ SY = στ, a 2ητ +θστ +θ +t(ηστ +η +2θσ SY = (στ 2 b SY = (σt+(t+τ( η 2 τ +ηθ(στ ++θ 2 σ+( στ 2 ( στ 3. (We added subscript SY to avoid confusion with our use of α in (4.2. Thisshowsthatϕ t (/z/z isacauchy-stieltjestransformofauniquecompactlysupported probability measure ν t. For a more detailed description of measure ν t and explicit formulas for its discrete and absolutely continuous components we refer to [2, Theorem 2.]; see also Remark 4.2 below. It is well known that a Cauchy-Stieltjes transform is an analytic function in the upper complex plane, determines measure uniquely, and if it extends to real z with z large enough then the corresponding moment generating function is well defined for all z small enough and is given by G νt (/z/z = ϕ t (z. This shows that ϕ t (z is the moment generating function of the probability measure ν t.

17 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 7 Next we observe that (3.3 in the operator notation is n H t (x n = +ηx+σx2 c k (tx n k +σt k= Writing c k (t = y k ν t (dy, we therefore get (4.6 H t (f(x = +ηx+σx2 f(y f(x ν t (dy. +σt y x Since the operator version of relation (2. is A t (x n+ = H t (x n + xa(x n, we derive (4.3 from (4.6 by induction on n; for a similar reasoning see [3, Lemma 2.4]. Remark 4.2. Denote by π t,η,θ,σ,τ (dx the univariate law of X t for the free quadratic harness (X t with parameters η,θ,σ,τ as in [9, Section 3]. Then ν t is given by (4.7 ν t (dx = t(t+τ (t2 +θtx+τx 2 π t,η,θ,σ,τ (dx. We read out this answer from [9, Eqtn. (3.4] using the following elementary relation between the Cauchy-Stieltjes transforms: If ν(dx = (ax 2 + bx + cπ(dx and m = xπ(dx then the Cauchy-Stieltjes transforms of π and ν are related by the formula (4.8 G ν (z = (az 2 +bz +cg π (z am az b. In our setting, m = 0, a = τ t(t+τ, b = θ/(t+τ, c = t/(t+τ, and [9, Eqtn. (3.4] gives τz +θt t[(+στ +2σtz +tη θ] (4.9 G π (z = τz 2 + +θtz +t2 2(σz 2 +ηz +(τz 2 +θtz +t 2 t [( στz (α+σβt β ατ] 2 4(+σt(t+τ(+αβ 2(σz 2 +ηz +(τz 2 +θtz +t 2. Inserting this expression into the right hand side of (4.8 we get the right hand side of (4.5. Uniqueness of Cauchy-Stieltjes transform implies ( Some other cases Several other cases can be worked out by a similar technique based on (2.2. Recasting [8, Section IV] in our notation, for m = 0,,... we have (5. D m+ F FD m+ = (m+d m, and if then (5.2 A = H = m= k= c k k! D k c m (m+! Dm+.

18 8 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI The generators for a Lévy processes (ξ t with exponential moments act on polynomials in variable x via κ( x, where tκ(θ = loge(exp(θξ t; this well-known formula appears e.g. [, Section 3]. 5.. Centered Poisson process. For example, the quadratic harness with γ = and η = σ = τ = 0 which corresponds to the Poisson process can be also analyzed by the algebraic technique. Equation (2.4 assumes the form (5.3 HF FH = E+θH We search a solution H of the form c H = k k! D k = ϕ(d. k= Due to (5. equation (5.3 implies ϕ (d = +θϕ(d ( and the solution with ϕ(0 = 0 is ϕ(d = θ e θd. We get It follows from (5.2 that H = θ ( e θd E. Consequently, A = m= θ m (m+! Dm+ = θ 2 m=2 θ m m! Dm. A = θ 2 ( e θd E θd A non-lévy example. Here we consider a quadratic harness with parameters γ =, τ = σ = 0. This quadratic harness appeared under the name quantum Bessel process in [4] (see also [20] for a multidimensional version, and as classical bi-poisson process in []. Then equation (2.4 assumes the form (5.4 HF FH = E+θH+η(F th = E+ηF+(θ tηh. For H = (E+ηF H we obtain (E+ηF( HF F H = (E+ηF(E+(θ tη H. Comparing this with (5.3 we conclude that by uniqueness the solution of (5.4 is H = θ tη (e (E+ηF (θ tηd E. From (3.8 it follows that A =: A(H = (E+ηFA( H. Therefore we obtain A = (θ tη (E+ηF (e (θ tηd E (θ tηd 2.

19 INFINITESIMAL GENERATORS FOR POLYNOMIAL PROCESSES 9 Acknowledgement JW research was supported in part by NCN grant 202/05/B/ST/ WB research was supported in part by the Taft Research Center at the University of Cincinnati. References [] M. Anshelevich, Generators of some non-commutative stochastic processes, Probability Theory Related Fields, 57 (203, pp Arxiv preprint arxiv: [2] D. Bakry and O. Mazet, Characterization of Markov semigroups on R associated to some families of orthogonal polynomials, in Séminaire de Probabilités XXXVII, Springer, 2003, pp [3] P. Barrieu and W. Schoutens, Iterates of the infinitesimal generator and space time harmonic polynomials of a Markov process, Journal of Computational and Applied Mathematics, 86 (2006, pp [4] P. Biane, Quelques proprietes du mouvement brownien non-commutatif., Astérisque, 236 (996, pp [5] P. Biane and R. Speicher, Stochastic calculus with respect to free Brownian motion and analysis on Wigner space, Probability Theory Related Fields, 2 (998, pp [6] M. Bożejko, B. Kümmerer, and R. Speicher, q-gaussian processes: non-commutative and classical aspects, Comm. Math. Phys., 85 (997, pp [7] W. Bryc, Markov processes with free Meixner laws, Stochastic Processes and Applications, 20 (200, pp [8] W. Bryc, W. Matysiak, and J. Weso lowski, Quadratic harnesses, q-commutations, and orthogonal martingale polynomials, Transactions of the American Mathematical Society, 359 (2007, pp [9] W. Bryc, W. Matysiak, and J. Weso lowski, Free quadratic harness, Stochastic Processes and Applications, 2 (20, pp [0] W. Bryc and J. Weso lowski, Conditional moments of q-meixner processes, Probability Theory Related Fields, 3 (2005, pp [] W. Bryc and J. Weso lowski, Classical bi-poisson process: an invertible quadratic harness, Statistics & Probability Letters, 76 (2006, pp [2] W. Bryc and J. Weso lowski, Askey Wilson polynomials, quadratic harnesses and martingales, Annals of Probability, 38 (200, pp [3] W. Bryc and J. Weso lowski, Infinitesimal generators of q-meixner processes, Stochastic Processes and Applications, 24 (204, pp [4] C. Cuchiero, Affine and polynomial processes, PhD thesis, ETH ZURICH, 20. [5] C. Cuchiero, M. Keller-Ressel, and J. Teichmann, Polynomial processes and their applications to mathematical finance, Finance and Stochastics, 6 (202, pp [6] E. Di Nardo, On a representation of time space-harmonic polynomials via symbolic Lévy processes, Sci. Math. Jpn., 76 (203, pp [7] E. Di Nardo and I. Oliva, A new family of time-space harmonic polynomials with respect to Lévy processes, Ann. Mat. Pura Appl. (4, 92 (203, pp [8] P. Feinsilver, Operator calculus, Pacific J. Math., 78 (978, pp [9] R. Mansuy and M. Yor, Harnesses, Lévy bridges and Monsieur Jourdain, Stochastic Process. Appl., 5 (2005, pp [20] W. Matysiak and M. Świeca, Zonal polynomials and a multidimensional quantum Bessel process, Stochastic Processes and their Applications, 25 (205, pp [2] N. Saitoh and H. Yoshida, The infinite divisibility and orthogonal polynomials with a constant recursion formula in free probability theory, Probab. Math. Statist., 2 (200, pp [22] W. Schoutens, Stochastic processes and orthogonal polynomials, Springer Verlag, [23] W. Schoutens and J. L. Teugels, Lévy processes, polynomials and martingales, Stochastic Models, 4 (998, pp [24] A. Sengupta, Time-space harmonic polynomials for continuous-time processes and an extension, Journal of Theoretical Probability, 3 (2000, pp [25] A. Sengupta, Markov processes, time space harmonic functions and polynomials, Statistics & Probability Letters, 78 (2008, pp

20 20 W LODZIMIERZ BRYC AND JACEK WESO LOWSKI [26] A. Sengupta and A. Sarkar, Finitely polynomially determined Lévy processes, Electron. J. Probab, 6 (200, pp. 22. [27] J. L. Solé and F. Utzet, On the orthogonal polynomials associated with a Lévy process, The Annals of Probability, 36 (2008, pp [28] J. L. Solé and F. Utzet, Time space harmonic polynomials relative to a Lévy process, Bernoulli, 4 (2008, pp. 3. [29] P. J. Szab lowski, Lévy processes, martingales, reversed martingales and orthogonal polynomials, arxiv preprint arxiv:22.32, (202. [30] P. J. Szab lowski, On stationary Markov processes with polynomial conditional moments, arxiv preprint arxiv: , (203. [3] P. J. Szab lowski, Markov processes, polynomial martingales and orthogonal polynomials, arxiv preprint arxiv:40.673, (204. [32] P. J. Szab lowski, On Markov processes with polynomials conditional moments, Trans. Amer. Math. Soc. 367 (205, pp (205. Department of Mathematics, University of Cincinnati, PO Box 20025, Cincinnati, OH , USA address: Wlodzimierz.Bryc@UC.edu Faculty of Mathematics and Information Science Warsaw University of Technology pl. Politechniki Warszawa, Poland address: wesolo@alpha.mini.pw.edu.pl