, , , January 7, 2005; Revised: April 3, 2005

Hammersley [19] introduced harnesses on

R^{n}

as the probabilistic models of long-range misorientation in the crystalline structure of metals. Several authors studied mathematical aspects of the concept: Mansuy and Yor [24] analyze harnesses on

R_{+}

, Williams [33] analyzes harnesses with arbitrary discrete index set; see also [14] , [34] , [35] , [36] .

The class of random fields on

R_{+}

which we call the quadratic harnesses is related to Hammersley's harnesses by mimicking the relation between the martingale and the quadratic martingale conditions. Such processes have already been studied implicitly by several authors, see the paragraph following Definition 2.2 . In particular, ref. [12] gives a construction of the three parameter family of Markov processes with the quadratic harness property. Some of these processes are related to the free Lévy processes, and some correspond to the non-commutative

q

-Gaussian processes introduced by Frisch and Bourret [17] and studied in [6] .

In full generality, the constructions of the quadratic harnesses on

R_{+}

are not yet well understood. In this paper we concentrate on uniqueness without dealing with the issue of existence. We show that quadratic harnesses are described by five numerical constants, which, under appropriate integrability conditions, determine the process uniquely. We study the related integrability properties of a slightly wider class of processes, improving earlier results of that type [8,Corollary4] , [10,Theorem2] , [32,Theorem2( 1 ∘ )] . We show that martingale polynomials associated with a quadratic harness with finite moments of all orders lead to a

q

-commutation equation

\begin{matrix} [x, y]_{q} = I + τ x^{2} + σ y^{2} + θ x + η y, \end{matrix}

(1.1)

where

x, y

are non-commutative variables,

[x, y]_{q} = x y - q y x,

and

q, η, θ, σ, τ

are the numerical constants that describe the quadratic harness.

This hints on more connections with non-commutative probability,

q

-Fock space constructions, and classical versions of non-commutative processes, see [6] . In this paper we use ( 1.1 ) to prove that quadratic harnesses with finite moments of all orders have orthogonal martingale polynomials, and to derive their three step recurrences. The explicit three terms recurrences that are associated with quadratic harnesses include a four-parameter family of the polynomials in Section 4.2 , and a four-parameter family of

q

-orthogonal polynomials in Section 4.3 . Based on one of the explicit recurrences from this paper, in [9] we extend the construction from [11] to the construction of the quadratic harness which we call the bi-Poisson processes.

Some of our arguments rely on the symmetries of the problem, but perhaps we did not explore the symmetries deeply enough. The family of quadratic harnesses that we study is invariant under the action of the translations and reflections of

R

, and the associated affine Hecke algebra is known to be associated with the Askey-Wilson polynomials, see [25] . Since the explicit recurrences that we found in Section 4 give sub-classes of the reparametrized Askey-Wilson polynomials, it is plausible that the theory of affine root systems [23] might lead to additional progress.

The paper is organized as follows. In Section 2 we define quadratic harnesses and state the main results. Proofs of the main results are given in Section 3 . In Section 4 we derive explicit recurrences for the orthogonal martingale polynomials.

We also introduce an operator technique that resembles umbral calculus as the tool that simplifies the proof of the quadratic harness property.

2 Definitions and main results

2.1 Harnesses and quadratic harnesses

Let

(X_{t})_{t > 0}

be a separable square integrable stochastic process with

σ

-fields:

ℱ_{s, u} = σ {X_{t} : t \in (0, s] \cup [u, \infty)}

. Consider the following two functionals

\begin{matrix} Δ_{t, u} = Δ_{t, u} (X) & = & \frac{X_{u} - X_{t}}{u - t}, \end{matrix}

(2.1)

\begin{matrix} {\tilde{Δ}}_{t, u} = {\tilde{Δ}}_{t, u} (X) & = & \frac{u X_{t} - t X_{u}}{u - t} . \end{matrix}

(2.2)

Both expressions appear in the statement of Theorem 2.2 below and the second one is the time-inverse of the first one. Namely, if

({\tilde{X}}_{t})_{t > 0}

denotes the time-inverse process

(t X_{1 / t})_{t > 0}

, then

{\tilde{Δ}}_{t, u} (X) = Δ_{1 / u, 1 / t} (\tilde{X}), 0 < t < u .

Definition 2.1 ([24] ). An integrable process

(X_{t})_{t > 0}

is a (simple) harness if

E [Δ_{s, t} | ℱ_{r, u}] = Δ_{r, u}

for every

r < s < t < u

A trivial example of a harness is a Gaussian process with covariance ( 2.5 ); additional examples follow Definition 2.2 .

As pointed out in [24] , the harness condition is equivalent to the linearity of regression property

\begin{matrix} E [X_{t} | ℱ_{s, u}] = a t, s, u X_{s} + b t, s, u X_{u}, \end{matrix}

(2.3)

where the coefficients are given by

\begin{matrix} a t, s, u = \frac{u - t}{u - s}, b t, s, u = \frac{t - s}{u - s} . \end{matrix}

(2.4)

It is clear that

a t, s, u + b t, s, u = 1

and

s a t, s, u + u b t, s, u = t

; these two identities will often be used. Since condition ( 2.3 ) is invariant under the time-inversion

(X_{t}) \mapsto (t X_{1 / t})

, another equivalent condition for a harness is

E [{\tilde{Δ}}_{s, t} | ℱ_{r, u}] = {\tilde{Δ}}_{r, u}, r < s < t < u .

The general form of the covariance of a square-integrable harness is as follows.

Proposition 2.1. If

(X_{t})

is a square-integrable centered harness on

R_{+}

, then there are constants

c_{0}, c_{1}, c_{2}, c_{3}

such that for

s \leq t

we have

\begin{matrix} E (X_{t} X_{s}) = c_{0} + c_{1} s + c_{2} t + c_{3} s t . \end{matrix}

(2.5)

Proof. Multiplying ( 2.3 ) by $X_{s}$ and averaging we get $E (X_{s} X_{t}) = a t, s, u E (X_{s}^{2}) + b t, s, u E (X_{s} X_{u})$ . Thus $E (X_{s} X_{t}) = α (s) + β (s) t$ is linear in $t$ and $E (X_{t} X_{u}) = a t, s, u E (X_{s} X_{u}) + b t, s, u E (X_{u}^{2}) = γ (u) + δ (u) t$ is linear in $t$ .
This gives a functional equation $α (s) + β (s) t = γ (t) + δ (t) s$ which, by taking $s = 0$ , gives $γ (t) = α (0) + β (0) t$ . Thus $E (X_{s} X_{t}) = c_{0} + c_{1} s + c_{2} t + c_{3} s t$ . □

For at centered standardized square-integrable

Z

that is independent of a quadratic harness

(X_{t})

, taking

({\tilde{X}}_{t}) = (X_{t} + (t + 1) Z)

we get a quadratic harness with

E ({\tilde{X}}_{s} {\tilde{X}}_{t}) = E (X_{s} X_{t}) + 1 + t + s + t s

. To avoid such non-uniqueness, throughout this paper we assume

\begin{matrix} E (X_{t}) = 0, E (X_{t} X_{s}) = min {t, s} . \end{matrix}

(2.6)

If ( 2.6 ) and ( 2.3 ) hold true and

0 < s \leq t \leq u

, then

\begin{matrix} E [X_{t} | ℱ_{\geq u}] = \frac{t}{u} X_{u} \end{matrix}

(2.7)

and

\begin{matrix} E [X_{t} | ℱ_{\leq s}] = X_{s}, \end{matrix}

(2.8)

see [12,(4)and(5)] . From martingale ( 2.8 ) and reverse martingale ( 2.7 ) condition it follows that the limits

X_{0} : = {lim}_{t ↘ 0} X_{t}

, and

{lim}_{t \to \infty} X_{t} / t

exist with probability one. Under assumption ( 2.6 ) we have

\begin{matrix} {lim}_{t ↘ 0} X_{t} = 0, {lim}_{t \to \infty} X_{t} / t = 0, \end{matrix}

(2.9)

so without loss of generality we may extend

(X_{t})_{t > 0}

to include the value

X_{0} = 0

when convenient. Similarly, we may extend

({\tilde{X}}_{t})_{t > 0} = (t X_{1 / t})_{t > 0}

to include the value

{\tilde{X}}_{0} = 0

corresponding to

t = 0

We now turn to the quadratic harness condition. The familiar quadratic martingale condition associated with the martingale property ( 2.8 ) can be written as

\begin{matrix} E [X_{t}^{2} | ℱ_{\leq s}] = X_{s}^{2} + t - s . \end{matrix}

(2.10)

This suggests that the quadratic harness condition associated with the simple harness property ( 2.3 ) should be written as

\begin{matrix} E [X_{t}^{2} | ℱ_{s, u}] = Q_{t, s, u} (X_{s}, X_{u}), \end{matrix}

(2.11)

where

\begin{matrix} Q_{t, s, u} (x, y) = A t, s, u x^{2} + B t, s, u x y + C t, s, u y^{2} + D t, s, u x + E t, s, u y + F t, s, u \end{matrix}

(2.12)

is a quadratic form in variables

x, y

with time-dependent coefficients; under assumption ( 2.6 ) we trivially have

\begin{matrix} s A t, s, u + s B t, s, u + u C t, s, u + F t, s, u = t . \end{matrix}

(2.13)

Definition 2.2. A square integrable process

(X_{t})_{t > 0}

is a quadratic harness on

R_{+}

if it satisfies conditions ( 2.3 ) and ( 2.11 ).

The well-known examples of quadratic harnesses are the Wiener, Poisson, and Gamma processes. Ref. [32] identifies all quadratic harnesses with covariance ( 2.6 ) when the quadratic form on the right hand side of ( 2.11 ) is such that the corresponding conditional variance is a function of the increments

X_{u} - X_{s}

only; this adds to the already listed examples two additional Lévy processes: the Pascal, and Meixner processes. A related non-commutative form of condition ( 2.11 ) appears in [5] . The main result of [12] asserts that quadratic harnesses with covariance ( 2.6 ) which satisfy the quadratic martingale condition ( 2.10 ) are in fact uniquely determined

q

-Meixner Markov processes. A quadratic harness that does not satisfy condition ( 2.10 ) is analyzed in [11] .

2.2 Five-parameter representation

Generically, the quadratic form on the right hand side of ( 2.11 ) is determined uniquely up to five numerical parameters.

Theorem 2.2. Let

(X_{t})

be a quadratic harness with covariance ( 2.6 ).

Suppose that ( 2.11 ) holds with

F_{t, s, u} \neq 0

, for all

0 < s < t < u

, and that

1, X_{s}, X_{t}, X_{s} X_{t}, X_{s}^{2}, X_{t}^{2}

are linearly independent for all

0 < s < t

. Then there exist

η, θ \in R

σ, τ \geq 0

, and

q \leq 1 + 2 \sqrt{σ τ}

such that

\begin{matrix} Var [X_{t} | ℱ_{s, u}] = \end{matrix}

\begin{matrix} F_{t, s, u} (1 + η {\tilde{Δ}}_{s, u} + θ Δ_{s, u} + σ {\tilde{Δ}}_{s, u}^{2} + τ Δ_{s, u}^{2} - (1 - q) Δ_{s, u} {\tilde{Δ}}_{s, u}) \end{matrix}

(2.14)

for all

0 < s < t < u

, where

\begin{matrix} F_{t, s, u} = \frac{(u - t) (t - s)}{u (1 + σ s) + τ - q s} \end{matrix}

(2.15)

is the normalizing function, and

Δ_{s, u}

{\tilde{Δ}}_{s, u}

are given by ( 2.1 - 2.2 ).

Recall that the conditional variance of

X

with respect to a

σ

-field

ℱ

is defined as

Var [X | ℱ] = E [X^{2} | ℱ] - {(E [X | ℱ])}^{2}

Remark 2.1. Time-inversion

(X_{t}) \mapsto (t X_{1 / t})

preserves the class of quadratic harnesses, modifying the coefficients in ( 2.11 ). More precisely, suppose

(X_{t})

satisfies the assumptions of Theorem 2.2 , and let

{\tilde{X}}_{t} = t X_{1 / t}

be its time inverse.

Then

({\tilde{X}}_{t})

is a quadratic harness with respect to its

σ

-fields

{\tilde{ℱ}}_{s, u} = ℱ_{1 / u, 1 / s}

, and ( 2.14 ) holds with the roles of the parameters

(η, θ)

and

(σ, τ)

switched within each pair:

Var [{\tilde{X}}_{t} | {\tilde{ℱ}}_{s, u}] / {\tilde{F}}_{t, s, u}

= 1 + θ {\tilde{Δ}}_{s, u} (\tilde{X}) + η Δ_{s, u} (\tilde{X}) + τ {\tilde{Δ}}_{s, u}^{2} (\tilde{X}) + σ Δ_{s, u}^{2} (\tilde{X}) - (1 - q) Δ_{s, u} (\tilde{X}) {\tilde{Δ}}_{s, u} (\tilde{X}),

where

{\tilde{F}}_{t, s, u} = \frac{(u - t) (t - s)}{u (1 + τ s) + σ - q s}

Ref. [18] gives criteria for time-inversion invariance of Markov processes.

2.3 Orthogonal martingale polynomials

Suppose that a quadratic harness

(X_{t})

has moments of all orders and martingale polynomials

p_{n} (x; t)

of all degrees

n \geq 0

, that is

\begin{matrix} E [p_{n} (X_{t}; t) | ℱ_{\leq s}] = p_{n} (X_{s}; s), 0 < s < t . \end{matrix}

(2.16)

Clearly,

p_{0} = 1

and

p_{1} (x; t) = x

are natural initial choices, see ( 2.8 ). Since

x p_{n} (x; t)

is a polynomial of degree

n + 1

, it follows that

\begin{matrix} x p_{n} (x; t) = \sum_{k = 0}^{n + 1} C_{k, n} (t) p_{k} (x; t) . \end{matrix}

(2.17)

Theorem 2.3. Suppose a quadratic harness

(X_{t})

with covariance ( 2.6 ) and conditional variance ( 2.14 ) has finite moments of all orders and martingale polynomials

p_{n} (x; t)

. If for each

t > 0

the random variable

X_{t}

has infinite support, then recurrence ( 2.17 ) holds with the infinite matrices

C_{t} : = [\begin{matrix} C_{00} (t) & C_{01} (t) & C_{02} (t) & C_{0, 3} (t) & . . . \\ C_{10} (t) & C_{11} (t) & C_{12} (t) & C_{13} (t) & . . . \\ 0 & C_{21} (t) & C_{22} (t) & C_{23} (t) & . . . \\ 0 & 0 & C_{32} (t) & C_{33} (t) & . . . \\ 0 & 0 & 0 & C_{43} (t) & . . . \\ . . . & . . . & . . . & . . . \end{matrix}]

given by

\begin{matrix} C_{t} = t x + y, t > 0, \end{matrix}

(2.18)

and the infinite matrices

x = C_{1} - C_{0}

y = C_{0}

satisfy equation ( 1.1 ).

From Theorem 2.3 we derive a number of equations that eventually determine the orthogonal martingale polynomials. Namely, it is well known that for orthogonal martingale polynomials

p_{n} (x; t)

recurrence ( 2.17 ) holds with a tri-diagonal matrix

C

, see [13] . Writing ( 2.17 ) as

\begin{matrix} x p_{n} (x; t) = a_{n} (t) p_{n + 1} (x; t) + b_{n} (t) p_{n} (x; t) + c_{n} (t) p_{n - 1} (x; t), n \geq 0, \end{matrix}

(2.19)

from ( 2.18 ) we get

\begin{matrix} a_{n} (t) = σ α_{n + 1} t + β_{n + 1}, b_{n} (t) = γ_{n} t + δ_{n}, c_{n} (t) = ɛ_{n} t + φ_{n}, \end{matrix}

(2.20)

and the coefficients in ( 2.20 ) satisfy a number of equations that result from the

q

-commutation equation ( 1.1 ).

Setting

p_{- 1} (x; t) = 0

p_{0} (x; t) = 1

p_{1} (x; t) = x

, and using ( 2.6 ) we see that the initial values are given by

\begin{matrix} α_{1} = 0, β_{1} = 1, γ_{0} = δ_{0} = 0, ɛ_{1} = 1, φ_{1} = 0 . \end{matrix}

(2.21)

For

n \geq 1

, ( 1.1 ) implies that coefficients

α_{n}, β_{n}, γ_{n}, δ_{n}, φ_{n}, ɛ_{n}

satisfy

\begin{matrix} σ^{2} τ α_{n} α_{n + 1} + σ α_{n} β_{n + 1} q + σ β_{n} β_{n + 1} = σ α_{n + 1} β_{n}, \end{matrix}

(2.22)

\begin{matrix} β_{n + 1} γ_{n + 1} + σ α_{n + 1} δ_{n} \end{matrix}

(2.23)

\begin{matrix} = σ α_{n + 1} (γ_{n} + γ_{n + 1}) τ + (σ α_{n + 1} δ_{n + 1} + β_{n + 1} γ_{n}) q + β_{n + 1} (δ_{n} + δ_{n + 1}) σ \end{matrix}

\begin{matrix} + σ α_{n + 1} θ + β_{n + 1} η, \end{matrix}

\begin{matrix} β_{n + 1} ɛ_{n + 1} + γ_{n} δ_{n} + σ α_{n} φ_{n} \end{matrix}

(2.24)

\begin{matrix} = (σ α_{n + 1} ɛ_{n + 1} + γ_{n}^{2} + σ α_{n} ɛ_{n}) τ + (σ α_{n + 1} φ_{n + 1} + γ_{n} δ_{n} + β_{n} ɛ_{n}) q \end{matrix}

\begin{matrix} + (β_{n + 1} φ_{n + 1} + δ_{n}^{2} + β_{n} φ_{n}) σ + γ_{n} θ + δ_{n} η + 1, \end{matrix}

\begin{matrix} γ_{n - 1} φ_{n} + δ_{n} ɛ_{n} \end{matrix}

(2.25)

\begin{matrix} = (γ_{n - 1} + γ_{n}) ɛ_{n} τ + (γ_{n} φ_{n} + δ_{n - 1} ɛ_{n}) q + (δ_{n - 1} + δ_{n}) φ_{n} σ + ɛ_{n} θ + φ_{n} η, \end{matrix}

\begin{matrix} ɛ_{n} φ_{n + 1} = ɛ_{n} ɛ_{n + 1} τ + ɛ_{n + 1} φ_{n} q + φ_{n} φ_{n + 1} σ . \end{matrix}

(2.26)

a_{n} (t) \neq 0

for all

n

, then the three step recurrence ( 2.19 ) defines a family of polynomials in variable

x

. We now show that these are indeed the martingale orthogonal polynomials for

(X_{t})

Theorem 2.4. Suppose a quadratic harness

(X_{t})

with covariance ( 2.6 ) and conditional variance ( 2.14 ) has finite moments of all orders, parameters

q, η, θ, σ, τ

are such that

(q, σ τ) \neq (- 1, 1)

, and the equations ( 2.21 - 2.26 ) have a solution such that

a_{n} (t) \neq 0

for all

t \geq 0

n \geq 0

Let

{p_{n} (x; t) : n \geq 0, t > 0}

satisfy ( 2.19 ) and ( 2.20 ) for

n = 0, 1, . . .

with

p_{- 1} (x, t) = 0, p_{0} (x, t) = 1

p_{1} (x, t) = x

. Then

{p_{n} (x; t)}

are orthogonal martingale polynomials for

(X_{t})

In Section 4 , we give sufficient conditions in terms of parameters

q, η, θ, σ, τ

for the assumptions of Theorem 2.4 to be satisfied, and we derive explicit three step recurrences for the orthogonal martingale polynomials in several cases of interest.

Theorem 2.4 implies that if quadratic harnesses

(X_{t})_{t > 0}

and

(Y_{t})_{t > 0}

satisfy its assumptions with the same parameters

η, θ, σ, τ, q

, then

E (X_{t_{1}}^{n_{1}} X_{t_{2}}^{n_{2}} . . . X_{t_{k}}^{n_{k}}) = E (Y_{t_{1}}^{n_{1}} Y_{t_{2}}^{n_{2}} . . . Y_{t_{k}}^{n_{k}})

for all

k \geq 1

0 < t_{1} < t_{2} < \cdot \cdot \cdot < t_{k}

, and

n_{1}, . . ., n_{k} \in N

. Under appropriate integrability conditions, this implies that

(X_{t})

is Markov and is uniquely determined by the parameters

η, θ, σ, τ, q

. In particular, this is the case when

{sup}_{n} | b_{n} (t) | < \infty

and

{sup}_{n} | a_{n - 1} (t) c_{n} (t) | < \infty

, as in this case the recurrence ( 2.19 ) corresponds to a compactly supported measure, see [4,Section2] . However, the question of existence of such a process is non-trivial and the constructions are known only in special cases, see [9] , [11] , [12] .

Next, we show that the integrability assumption of Theorem 2.4 is automatically satisfied if

σ τ = 0

. We remark that the result stated below does not use the full power of the quadratic harness condition, and generalizes [10,Theorem2] and [32,Theorem2( 1 ∘ )] . Since

Var [X_{t} | ℱ_{s, u}] = E [X_{t}^{2} | ℱ_{s, u}] - (a_{t, s, u} X_{s} + b_{t, s, u} X_{u})^{2},

using ( 2.9 ) we can pass to the limit in ( 2.14 ) as

u \to \infty

or as

s \to 0

. This gives

\begin{matrix} Var [X_{t} | ℱ_{\leq s}] & = & \frac{t - s}{1 + σ s} (σ X_{s}^{2} + η X_{s} + 1), \end{matrix}

(2.27)

\begin{matrix} Var [X_{t} | ℱ_{\geq u}] & = & \frac{t (u - t)}{u + τ} (τ \frac{X_{u}^{2}}{u^{2}} + θ \frac{X_{u}}{u} + 1) . \end{matrix}

(2.28)

Theorem 2.5. If a square-integrable stochastic process

(X_{t})

with covariance ( 2.6 ) satisfies ( 2.7 ), ( 2.8 ), ( 2.27 ) and ( 2.28 ) with

σ, τ \geq 0

such that

σ τ \leq 1 / 2^{4 r + 10}

for some

r > 2

then

E (| X_{t} |^{r}) < \infty

for all

t > 0

. In particular, if

σ τ = 0

, then

E (| X_{t} |^{r}) < \infty

for all

r, t > 0

3 Proofs of the main results

3.1 Proof of Theorem 2.2

For

0 < x < y < z

define

\begin{matrix} σ_{y, x, z} : = & \frac{A_{y, x, z} + B_{y, x, z} + C_{y, x, z} - 1}{F_{y, x, z}}, \end{matrix}

\begin{matrix} τ_{y, x, z} : = & \frac{x^{2} A_{y, x, z} + x z B_{y, x, z} + z^{2} C_{y, x, z} - y^{2}}{F_{y, x, z}}, \end{matrix}

\begin{matrix} 1 + q_{y, x, z} : = & \frac{B_{y, x, z} (z - x)}{F_{y, x, z}}, \end{matrix}

\begin{matrix} η_{y, x, z} : = & \frac{D_{y, x, z} + E_{y, x, z}}{F_{y, x, z}}, \end{matrix}

\begin{matrix} θ_{y, x, z} : = & \frac{x D_{y, x, z} + z E_{y, x, z}}{F_{y, x, z}} . \end{matrix}

We will show that the left hand sides of the above equations do not depend on the arguments

y, x, z

. To this end we prove three claims. Throughout the proofs of all three Claims,

0 < r < s < t < u

are arbitrary numbers.

Claim 3.1. For all

0 < x < y < z

, and for

f = σ, τ, q, η, θ

we have

\begin{matrix} f_{y, x, z} = f_{\tilde{y}, x, z}, \end{matrix}

(3.1)

provided

x < \tilde{y} < z

Proof of Claim 3.1 . Observe that from ( 2.3 ) and ( 2.11 ) we get

E [X_{s} X_{t} | ℱ_{r, u}] = E [X_{s} E [X_{t} | ℱ_{s, u}] | ℱ_{r, u}]

= a_{t, s, u} E [X_{s}^{2} | ℱ_{r, u}] + b_{t, s, u} X_{u} E [X_{s} | ℱ_{r, u}]

= a_{t, s, u} (A_{s, r, u} X_{r}^{2} + B_{s, r, u} X_{r} X_{u} + C_{s, r, u} X_{u}^{2} + D_{s, r, u} X_{r} + E_{s, r, u} X_{u} + F_{s, r, u})

+ b_{t, s, u} X_{u} (a_{s, r, u} X_{r} + b_{s, r, u} X_{u}) = a_{t, s, u} A_{s, r, u} X_{r}^{2} + (a_{t, s, u} B_{s, r, u} + b_{t, s, u} a_{s, r, u}) X_{r} X_{u}

+ (a_{t, s, u} C_{s, r, u} + b_{t, s, u} b_{s, r, u}) X_{u}^{2} + a_{t, s, u} D_{s, r, u} X_{r} + a_{t, s, u} E_{s, r, u} X_{u} + a_{t, s, u} F_{s, r, u} .

On the other hand,

E [X_{s} X_{t} | ℱ_{r, u}] = E [E [X_{s} | ℱ_{r, t}] X_{t} | ℱ_{r, u}] = a_{s, r, t} X_{r} E [X_{t} | ℱ_{r, u}] + b_{s, r, t} E [X_{t}^{2} | ℱ_{r, u}]

= a_{s, r, t} X_{r} (a_{t, r, u} X_{r} + b_{t, r, u} X_{u})

+ b_{s, r, t} (A_{t, r, u} X_{r}^{2} + B_{t, r, u} X_{r} X_{u} + C_{t, r, u} X_{u}^{2} + D_{t, r, u} X_{r} + E_{t, r, u} X_{u} + F_{t, r, u})

= (a_{s, r, t} a_{t, r, u} + b_{s, r, t} A_{t, r, u}) X_{r}^{2} + (a_{s, r, t} b_{t, r, u} + b_{s, r, t} B_{t, r, u}) X_{r} X_{u}

+ b_{s, r, t} C_{t, r, u} X_{u}^{2} + b_{s, r, t} D_{t, r, u} X_{r} + b_{s, r, t} E_{t, r, u} X_{u} + b_{s, r, t} F_{t, r, u} .

Comparing the coefficients at

X_{r}^{2}

X_{r} X_{u}

X_{u}^{2}

X_{r}

X_{u}

and 1 in the above expressions we get

\begin{matrix} a_{t, s, u} A_{s, r, u} & = & b_{s, r, t} A_{t, r, u} + a_{s, r, t} a_{t, r, u}, \end{matrix}

(3.2)

\begin{matrix} a_{t, s, u} B_{s, r, u} & = & b_{s, r, t} B_{t, r, u}, \end{matrix}

(3.3)

\begin{matrix} b_{t, s, u} b_{s, r, u} + a_{t, s, u} C_{s, r, u} & = & b_{s, r, t} C_{t, r, u}, \end{matrix}

(3.4)

\begin{matrix} a_{t, s, u} D_{s, r, u} & = & b_{s, r, t} D_{t, r, u}, \end{matrix}

(3.5)

\begin{matrix} a_{t, s, u} E_{s, r, u} & = & b_{s, r, t} E_{t, r, u}, \end{matrix}

(3.6)

\begin{matrix} a_{t, s, u} F_{s, r, u} & = & b_{s, r, t} F_{t, r, u} . \end{matrix}

(3.7)

(In the derivation of ( 3.3 ) we used the fact that

b_{t, s, u} a_{s, r, u} = a_{s, r, t} b_{t, r, u}

Adding ( 3.2 ), ( 3.3 ) and ( 3.4 ), we get

a_{t, s, u} (A_{s, r, u} + B_{s, r, u} + C_{s, r, u} - 1) + a_{t, s, u} + b_{t, s, u} b_{s, r, u}

= b_{s, r, t} (A_{t, r, u} + B_{t, r, u} + C_{t, r, u} - 1) + b_{s, r, t} + a_{s, r, t} a_{t, r, u} .

Since a calculation shows that

\begin{matrix} a_{t, s, u} + b_{t, s, u} b_{s, r, u} = b_{s, r, t} + a_{s, r, t} a_{t, r, u}, \end{matrix}

(3.8)

dividing by ( 3.7 ), we get

σ_{s, r, u} = σ_{t, r, u}

which proves ( 3.1 ) when

f = σ

We proceed similarly when

f = τ, q, η, θ

. Adding ( 3.2 ) multiplied by

r^{2}

, ( 3.3 ) multiplied by

r u

( 3.4 ) multiplied by

u^{2}

, and dividing by ( 3.7 ), we obtain

τ_{s, r, u} = τ_{t, r, u}

, after noticing that

s^{2} a_{t, s, u} + u^{2} b_{t, s, u} b_{s, r, u} = t^{2} b_{s, r, t} + r^{2} a_{s, r, t} a_{t, r, u}

Equation ( 3.3 ) multiplied by

(u - r)

and divided by ( 3.7 ) gives

q_{s, r, u} = q_{t, r, u}

. Adding equations ( 3.5 ) and ( 3.6 ) and dividing by ( 3.7 ) gives

η_{s, r, u} = η_{t, r, u}

and similarly multiplying ( 3.5 ) by

r

and ( 3.6 ) by

u

after dividing by ( 3.7 ) gives

θ_{s, r, u} = θ_{t, r, u}

Thus we obtained

f_{s, r, u} = f_{t, r, u}

for

f = σ, τ, q, η, θ

. If

0 < x < y < z

then substitution

r = x, s = y, u = z

yields ( 3.1 ) for

\tilde{y} \in (y, z)

; substitution

r = x, t = y, u = z

gives ( 3.1 ) for

\tilde{y} \in (x, y)

, completing the proof of Claim ( 3.1 ). □

Claim 3.2. For all

0 < x < y < z

, and

f = σ, τ, q, η, θ

we have

\begin{matrix} f_{y, x, z} = f_{y, \tilde{x}, z}, \end{matrix}

(3.9)

provided

0 < \tilde{x} < y

Proof of Claim 3.2 . Now consider the identity

E [X_{t}^{2} | ℱ_{r, u}] = E [E [X_{t}^{2} | ℱ_{s, u}] | ℱ_{r, u}] .

Using ( 2.3 ) and ( 2.11 ) we see that

E [E [X_{t}^{2} | ℱ_{s, u}] | ℱ_{r, u}]

is given by

E [A_{t, s, u} X_{s}^{2} + B_{t, s, u} X_{s} X_{u} + C_{t, s, u} X_{u}^{2} + D_{t, s, u} X_{s} + E_{t, s, u} X_{u} + F_{t, s, u} | ℱ_{r, u}]

= A_{t, s, u} E [X_{s}^{2} | ℱ_{r, u}] + B_{t, s, u} X_{u} E [X_{s} | ℱ_{r, u}] + C_{t, s, u} X_{u}^{2} + D_{t, s, u} E [X_{s} | ℱ_{r, u}]

+ E_{t, s, u} X_{u} + F_{t, s, u} = A_{t, s, u} (A_{s, r, u} X_{r}^{2} + B_{s, r, u} X_{r} X_{u} + C_{s, r, u} X_{u}^{2} + D_{s, r, u} X_{r})

(+ E_{s, r, u} X_{u} + F_{s, r, u}) + B_{t, s, u} X_{u} (a_{s, r, u} X_{r} + b_{s, r, u} X_{u}) + C_{t, s, u} X_{u}^{2}

+ D_{t, s, u} (a_{s, r, u} X_{r} + b_{s, r, u} X_{u}) + E_{t, s, u} X_{u} + F_{t, s, u} = A_{t, s, u} A_{s, r, u} X_{r}^{2}

+ (A_{t, s, u} B_{s, r, u} + B_{t, s, u} a_{s, r, u}) X_{r} X_{u} + (A_{t, s, u} C_{s, r, u} + B_{t, s, u} b_{s, r, u} + C_{t, s, u}) X_{u}^{2}

+ (A_{t, s, u} D_{s, r, u} + D_{t, s, u} a_{s, r, u}) X_{r} + (A_{t, s, u} E_{s, r, u} + D_{t, s, u} b_{s, r, u} + E_{t, s, u}) X_{u}

+ A_{t, s, u} F_{s, r, u} + F_{t, s, u} .

On the other hand,

E [X_{t}^{2} | ℱ_{r, u}] = A_{t, r, u} X_{r}^{2} + B_{t, r, u} X_{r} X_{u} + C_{t, r, u} X_{u}^{2} + C_{t, r, u} X_{r} + D_{t, r, u} X_{u} + F_{t, r, u} .

Comparing the coefficients at

X_{r}^{2}

X_{r} X_{u}

X_{u}^{2}

X_{r}

X_{u}

and

1

, we obtain

\begin{matrix} A_{t, r, u} & = & A_{t, s, u} A_{s, r, u}, \end{matrix}

(3.10)

\begin{matrix} B_{t, r, u} & = & A_{t, s, u} B_{s, r, u} + B_{t, s, u} a_{s, r, u}, \end{matrix}

(3.11)

\begin{matrix} C_{t, r, u} & = & A_{t, s, u} C_{s, r, u} + B_{t, s, u} b_{s, r, u} + C_{t, s, u}, \end{matrix}

(3.12)

\begin{matrix} D_{t, r, u} & = & A_{t, s, u} D_{s, r, u} + D_{t, s, u} a_{s, r, u}, \end{matrix}

(3.13)

\begin{matrix} E_{t, r, u} & = & A_{t, s, u} E_{s, r, u} + D_{t, s, u} b_{s, r, u} + E_{t, s, u}, \end{matrix}

(3.14)

\begin{matrix} F_{t, r, u} & = & A_{t, s, u} F_{s, r, u} + F_{t, s, u} . \end{matrix}

(3.15)

Substituting the right hand sides of equation ( 3.10 - 3.15 ) for

A_{t, r, u}, B_{t, r, u}, . . ., F_{t, r, u}

on the right hand sides of ( 3.2 - 3.7 ) we get

\begin{matrix} A_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = a_{s, r, t} a_{t, r, u}, \end{matrix}

(3.16)

\begin{matrix} B_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = b_{s, r, t} B_{t, s, u} a_{s, r, u}, \end{matrix}

(3.17)

\begin{matrix} C_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = b_{s, r, t} (B_{t, s, u} b_{s, r, u} + C_{t, s, u}) - b_{t, s, u} b_{s, r, u}, \end{matrix}

(3.18)

\begin{matrix} D_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = b_{s, r, t} D_{t, s, u} a_{s, r, u}, \end{matrix}

(3.19)

\begin{matrix} E_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = b_{s, r, t} (D_{t, s, u} b_{s, r, u} + E_{t, s, u}), \end{matrix}

(3.20)

\begin{matrix} F_{s, r, u} (a_{t, s, u} - b_{s, r, t} A_{t, s, u}) & = b_{s, r, t} F_{t, s, u} . \end{matrix}

(3.21)

\begin{matrix}  \end{matrix}

We can now proceed analogously to the proof of Claim 3.1 . Namely, adding ( 3.16 ), ( 3.17 ), ( 3.18 ), again taking into account ( 3.8 ), and dividing by ( 3.21 ), we get

σ_{s, r, u} = σ_{t, s, u}

. Adding ( 3.16 ) multiplied by

r^{2}

, ( 3.17 ) multiplied by

r u

and ( 3.18 ) multiplied by

u^{2}

, and dividing by ( 3.21 ), we obtain

τ_{s, r, u} = τ_{t, s, u}

. Equation ( 3.17 ) multiplied by

(u - r)

and divided by ( 3.21 ) gives

q_{s, r, u} = q_{t, s, u}

. Adding equations ( 3.19 ) and ( 3.20 ) and dividing by ( 3.21 ) gives

η_{s, r, u} = η_{t, s, u}

and similarly multiplying ( 3.19 ) by

r

and ( 3.20 ) by

u

after dividing by ( 3.21 ) gives

θ_{s, r, u} = θ_{t, s, u}

Thus we obtained

f_{s, r, u} = f_{t, s, u}

for

f = σ, τ, q, η, θ

. By Claim 3.1

f_{s, r, u} = f_{t, r, u}

, so

f_{t, r, u} = f_{t, s, u}

, ending the proof of Claim 3.2 . □

Claim 3.3. For all

0 < x < y < z

, and

f = σ, τ, q, η, θ

we have

\begin{matrix} f_{y, x, z} = f_{y, x, \tilde{z}}, \end{matrix}

(3.22)

provided

\tilde{z} > y

Proof of Claim 3.3 . This follows from Claim 3.2 by the time-inversion $(X_{t}) \mapsto (t X_{1 / t})$ .
Alternatively, one can repeat the previous arguments, starting with the identity $E [X_{s}^{2} | ℱ_{r, u}] = E [E [X_{s}^{2} | ℱ_{r, t}] | ℱ_{r, u}] .$ □

Conclusion of proof of Theorem 2.2 . Now, it is easy to deduce that ( 3.1 ), ( 3.9 ) and ( 3.22 ) imply that functions

σ, τ, q, η, θ

are in fact constants. Indeed, given

0 < x_{1} < y_{1} < z_{1}

and

0 < x_{2} < y_{2} < z_{2}

, and

f = σ, τ, q, η, θ

, with

\tilde{x} : = min {x_{1}, x_{2}}

and

\tilde{z} : = max {z_{1}, z_{2}}

, we see that

\begin{matrix} f_{y_{1}, x_{1}, z_{1}} & = f_{y_{1}, x_{1}, \tilde{z}} & by Claimclaim 3 3.3, \end{matrix}

\begin{matrix} = f_{y_{1}, \tilde{x}, \tilde{z}} & by Claimclaim 2 3.2, \end{matrix}

\begin{matrix} = f_{y_{2}, \tilde{x}, \tilde{z}} & by Claimclaim 1 3.1, \end{matrix}

\begin{matrix} = f_{y_{2}, x_{2}, \tilde{z}} & by Claimclaim 2 3.2, \end{matrix}

\begin{matrix} = f_{y_{2}, x_{2}, z_{2}} & by Claimclaim 3 3.3. \end{matrix}

\begin{matrix}  \end{matrix}

\begin{matrix}  \end{matrix}

(If

x_{2} = \tilde{x}

z_{2} = \tilde{z}

, one or both of the last two steps is unnecessary.) So, for all

0 < s < t < u

we have

\begin{matrix} \begin{matrix} (3.23) & σ = \frac{A_{t, s, u} + B_{t, s, u} + C_{t, s, u} - 1}{F_{t, s, u}}, τ = \frac{s^{2} A_{t, s, u} + s u B_{t, s, u} + u^{2} C_{t, s, u} - t^{2}}{F_{t, s, u}}, 1 + q = \frac{B_{t, s, u} (u - s)}{F_{t, s, u}}, η = \frac{D_{t, s, u} + E_{t, s, u}}{F_{t, s, u}}, θ = \frac{s D_{t, s, u} + u E_{t, s, u}}{F_{t, s, u}} . \end{matrix} \end{matrix}

(3.24)

The above equations, along with ( 2.13 ), make a system of linear equations in variables

A_{t, s, u}, \dots, F_{t, s, u}

. This system must be solvable for all

0 < s < t < u

, so its determinant

(u - s)^{3} (u (1 + σ s) + τ - q s) \neq 0

. As

u \to \infty, s \to 0

, the expression

u (1 + σ s) + τ - q s

is positive. Thus

u (1 + σ s) + τ - q s > 0

. Taking the limits

s \to 0, u \to 0

s = 1, u \to \infty

, and

s \to u

, we get

τ \geq 0

σ \geq 0

and

q \leq 1 + σ u + τ / u

respectively. Minimizing the latter over

u > 0

we get

q \leq 1 + 2 \sqrt{σ τ}

The unique solution of the system of equations is given by ( 2.15 ), ( 3.25 - 3.29 ). A calculation shows that ( 2.14 ) holds. □

Remark 3.1. A calculation shows that formula ( 2.14 ) is equivalent to the following formulas, valid for all

0 < s < t < u

\begin{matrix} A_{t, s, u} & = & \frac{(u - t) [u (1 + σ t) + τ - q t]}{(u - s) [u (1 + σ s) + τ - q s]}, \end{matrix}

(3.25)

\begin{matrix} B_{t, s, u} & = & \frac{(u - t) (t - s) (1 + q)}{(u - s) [u (1 + σ s) + τ - q s]}, \end{matrix}

(3.26)

\begin{matrix} C_{t, s, u} & = & \frac{(t - s) [t (1 + σ s) + τ - q s]}{(u - s) [u (1 + σ s) + τ - q s]}, \end{matrix}

(3.27)

\begin{matrix} D_{t, s, u} & = & \frac{(u - t) (t - s) (u η - θ)}{(u - s) [u (1 + σ s) + τ - q s]}, \end{matrix}

(3.28)

\begin{matrix} E_{t, s, u} & = & \frac{(u - t) (t - s) (θ - s η)}{(u - s) [u (1 + σ s) + τ - q s]}, \end{matrix}

(3.29)

and

F_{t, s, u}

is given by ( 2.15 ).

3.2 Proof of Theorem 2.3

For non-commutative variables

x, y

, let

\begin{matrix} Q_{t, s, u} (x, y) = A_{t, s, u} x^{2} + B_{t, s, u} x y + C_{t, s, u} y^{2} + D_{t, s, u} x + E_{t, s, u} y + F_{t, s, u}; \end{matrix}

(3.30)

recall the quadratic form ( 2.12 ) in commuting variables

x, y

. (Note that the order of

x, y

B_{t, s, u}

is changed when compared to its dual version ( 4.35 ).)

Lemma 3.4. Under the assumptions of Theorem 2.3 ,

\begin{matrix} C_{t} = a_{t, s, u} C_{s} + b_{t, s, u} C_{u}, \end{matrix}

(3.31)

\begin{matrix} C_{t}^{2} = Q_{t, s, u} (C_{s}, C_{u}), \end{matrix}

(3.32)

where the coefficients in ( 3.30 ) are given by ( 2.15 ) and ( 3.25 - 3.29 ).

Proof. For a polynomial $φ : R \to R$ consider vectors $v_{s, t} = [E (φ (X_{s}) p_{0} (X_{t}; t)), E (φ (X_{s}) p_{1} (X_{t}; t)), . . ., E (φ (X_{s}) p_{k} (X_{t}; t)), . . .] \in R^{\infty} .$ Since for $s > 0$ random variable $X_{s}$ has infinite support, polynomials $1, X_{s}, X_{s}^{2}, . . .$ are linearly independent, and the corresponding orthogonal polynomials are non-degenerate. This implies that as we change $φ$ for a fixed $s > 0$ , vectors of the form $v_{s, s}$ are dense in $R^{\infty}$ , equipped with the product topology. Indeed, applying the Gram-Schmidt orthogonalization process to $p_{0} (x; s), p_{1} (x; s), . . .$ we get a sequence of orthogonal polynomials $q_{0}, q_{1}, . . .$ such that $E (q_{k} (X_{s})^{2}) \neq 0$ . Therefore, for any $w = [w_{0}, w_{1}, . . .] \in R^{\infty}$ , and $n \geq 0$ we can find numbers $u_{0}, u_{1}, . . ., u_{n}$ such that the first $n$ coordinates of $v_{s, s}$ corresponding to $φ = \sum_{j = 0}^{n} u_{j} q_{j}$ are equal to $w_{0}, w_{1}, . . ., w_{n}$ .
This implies that to verify the identities, it suffices to verify that the identities hold true when multiplied from the left by $v_{s, s}$ .
By the martingale polynomial property, $E (φ (X_{s}) p_{k} (X_{s}; s)) = E (φ (X_{s}) p_{k} (X_{t}; t))$ .
Therefore, from ( 2.17 ) it follows that $v_{s, s} = v_{s, t}$ and $v_{s, s} \times C_{t} = v_{s, t} \times C_{t}$ $= [E (φ (X_{s}) \sum_{j} C_{j, 0} (t) p_{j} (X_{t}; t)), E (φ (X_{s}) \sum_{j} C_{j, 1} (t) p_{j} (X_{t}; t)), . . .]$ $= [E (φ (X_{s}) X_{t} p_{0} (X_{t}; t)), E (φ (X_{s}) X_{t} p_{1} (X_{t}; t)), . . ., E (φ (X_{s}) X_{t} p_{k} (X_{t}; t)), . . .]$ $= [E (φ (X_{s}) X_{t} p_{0} (X_{u}; u)), E (φ (X_{s}) X_{t} p_{1} (X_{u}; u)), . . ., E (φ (X_{s}) X_{t} p_{k} (X_{u}; u)), . . .]$ $= a_{t, s, u} [E (φ (X_{s}) X_{s} p_{k} (X_{u}; u)) : k \geq 0] + b_{t, s, u} [E (φ (X_{s}) X_{u} p_{k} (X_{u}; u)) : k \geq 0] .$ Using the martingale polynomial property again, we see that $v_{s, s} \times C_{t}$ equals to $a_{t, s, u} [E (φ (X_{s}) X_{s} p_{k} (X_{s}; s)) : k \geq 0] + b_{t, s, u} [E (φ (X_{s}) X_{u} p_{k} (X_{u}; u)) : k \geq 0]$ $a_{t, s, u} v_{s, s} \times C_{s} + b_{t, s, u} v_{s, u} \times C_{u} = v_{s, s} \times (a_{t, s, u} C_{s} + b_{t, s, u} C_{u}),$ proving ( 3.31 ). Similar reasoning proves ( 3.32 ): $v_{s, s} \times C_{t}^{2} = v_{s, t} \times C_{t}^{2} = [E (φ (X_{s}) X_{t} p_{k} (X_{t}; t)) : k \geq 0] \times C_{t}$ $= [E (φ (X_{s}) X_{t}^{2} p_{k} (X_{t}; t)) : k \geq 0] = [E (φ (X_{s}) X_{t}^{2} p_{k} (X_{u}; u)) : k \geq 0]$ $= [E (φ (X_{s}) E [X_{t}^{2} | ℱ_{s, u}] p_{k} (X_{u}; u)) : k \geq 0]$ $= [E (φ (X_{s}) Q_{t, s, u} (X_{s}, X_{t}) p_{k} (X_{u}; u)) : k \geq 0]$ $= A_{t, s, u} [E (φ (X_{s}) X_{s}^{2} p_{k} (X_{s}; s)) : k \geq 0] + B_{t, s, u} [E (φ (X_{s}) X_{s} X_{u} p_{k} (X_{u}; u)) : k \geq 0]$ $+ C_{t, s, u} [E (φ (X_{s}) X_{u}^{2} p_{k} (X_{u}; u)) : k \geq 0] + D_{t, s, u} [E (φ (X_{s}) X_{s} p_{k} (X_{s}; s)) : k \geq 0]$ $+ E_{t, s, u} [E (φ (X_{s}) X_{u} p_{k} (X_{u}; u)) : k \geq 0] + F_{t, s, u} [E (φ (X_{s}) p_{k} (X_{s}; s)) : k \geq 0]$ $= A_{t, s, u} v_{s, s} \times C_{s}^{2} + B_{t, s, u} [E (φ (X_{s}) X_{s} p_{k} (X_{s}; s)) : k \geq 0] \times C_{u}$ $+ C_{t, s, u} v_{s, s} \times C_{u}^{2} + D_{t, s, u} v_{s, s} \times C_{s} + E_{t, s, u} v_{s, s} \times C_{u} + F_{t, s, u} v_{s, s}$ $= v_{s, s} \times (A_{t, s, u} C_{s}^{2} + B_{t, s, u} C_{s} C_{u} + C_{t, s, u} C_{u}^{2} + D_{t, s, u} C_{s} + E_{t, s, u} C_{u} + F_{t, s, u} I) .$ □

Lemma 3.5. Suppose a linear operator

C (t) = t x + y

satisfies relation ( 3.32 ) for some quadratic form

Q_{t, s, u}

given by ( 3.30 ). If the coefficients of

Q_{t, s, u}

are given by ( 2.15 ), ( 3.25 - 3.29 ) then ( 1.1 ) holds true.

Proof. Expanding the expressions like

C_{t}^{2} = t^{2} x^{2} + t (x y + y x) + y^{2}

and using ( 3.32 ) we get

\begin{matrix} t^{2} x^{2} + t (x y + y x) + y^{2} = A_{t, s, u} (s^{2} x^{2} + s (x y + y x) + y^{2}) \end{matrix}

(3.33)

\begin{matrix} + B_{t, s, u} (y^{2} + u y x + s x y + s u x^{2}) + C_{t, s, u} (u^{2} x^{2} + u (x y + y x) + y^{2}) \end{matrix}

\begin{matrix} + (D_{t, s, u} + E_{t, s, u}) y + (s D_{t, s, u} + u E_{t, s, u}) x + F_{t, s, u} I . \end{matrix}

Applying the relations ( 2.13 ), ( 3.23 ), and ( 3.24 ) to the coefficients at each of the monomials in ( 3.33 ) we get

x y + y x = τ x^{2} + σ y^{2} + η y + θ x + I + (1 + q) y x,

which is the same as ( 1.1 ). □

3.3 Proof of Theorem 2.4

First we shall show the martingale property ( 2.16 ). To this end, we proceed by induction on

n

. Trivially, 2.16 holds for

n = 0, 1

, see ( 2.8 ). Assume then that

n > 1

and 2.16 holds for

0, 1, \dots, n - 1

and all

0 < s < t

We start from calculating

E [X_{t} p_{n - 1} (X_{u}; u) | ℱ_{\leq s}]

in two ways. On one hand

\begin{matrix} E [X_{t} p_{n - 1} (X_{u}; u) | ℱ_{\leq s}] = E [E [X_{t} | ℱ_{s, u}] p_{n - 1} (X_{u}; u) | ℱ_{\leq s}] = \end{matrix}

\begin{matrix} = & E [(a_{t, s, u} X_{s} + b_{t, s, u} X_{u}) p_{n - 1} (X_{u}; u) | ℱ_{\leq s}] = \end{matrix}

\begin{matrix} = & a_{t, s, u} X_{s} p_{n - 1} (X_{s}; s) + b_{t, s, u} E [X_{u} p_{n - 1} (X_{u}; u) | ℱ_{\leq s}] = \end{matrix}

\begin{matrix} = & a_{t, s, u} (a_{n - 1} (s) p_{n} (X_{s}; s) + b_{n - 1} (s) p_{n - 1} (X_{s}; s) + c_{n - 1} p_{n - 2} (X_{s}; s)) + \end{matrix}

\begin{matrix} + & b_{t, s, u} (a_{n - 1} (u) E [p_{n} (X_{u}; u) | ℱ_{\leq s}] + b_{n - 1} (u) p_{n - 1} (X_{s}; s) + c_{n - 1} (u) p_{n - 2} (X_{s}; s)) . \end{matrix}

On the other hand,

\begin{matrix} E [X_{t} p_{n - 1} (X_{u}; u) | ℱ_{\leq s}] = E [X_{t} E [p_{n - 1} (X_{u}; u) | ℱ_{\leq t}] | ℱ_{\leq s}] = E [X_{t} p_{n - 1} (X_{t}; t) | ℱ_{\leq s}] \end{matrix}

\begin{matrix} = & a_{n - 1} (t) E [p_{n} (X_{t}; t) | ℱ_{\leq s}] + b_{n - 1} (t) p_{n - 1} (X_{s}; s) + c_{n - 1} (t) p_{n - 2} (X_{s}; s) . \end{matrix}

Thus comparing the right hand sides of the above equations we obtain the following equation:

b_{t, s, u} a_{n - 1} (u) E [p_{n} (X_{u}; u) | ℱ_{\leq s}] - a_{n - 1} (t) E [p_{n} (X_{t}; t) | ℱ_{\leq s}] =

= - p_{n} (X_{s}; s) a_{t, s, u} a_{n - 1} (s) + p_{n - 1} (X_{s}; s) (b_{n - 1} (t) - a_{t, s, u} b_{n - 1} (s) - b_{t, s, u} b_{n - 1} (u))

+ p_{n - 2} (X_{s}; s) (c_{n - 1} (t) - a_{t, s, u} c_{n - 1} (s) - b_{t, s, u} c_{n - 1} (u)) .

A trivial verification, using 2.20 , shows that

\begin{matrix} b_{n - 1} (t) - a_{t, s, u} b_{n - 1} (s) - b_{t, s, u} b_{n - 1} (u) & = & 0, \end{matrix}

\begin{matrix} c_{n - 1} (t) - a_{t, s, u} c_{n - 1} (s) - b_{t, s, u} c_{n - 1} (u) & = & 0 . \end{matrix}

Hence denoting

X : = E [p_{n} (X_{u}; u) | ℱ_{\leq s}], Y : = E [p_{n} (X_{t}; t) | ℱ_{\leq s}],

we have

\begin{matrix} b_{t, s, u} a_{n - 1} (u) X - a_{n - 1} (t) Y = - a_{t, s, u} a_{n - 1} (s) p_{n} (X_{s}; s) . \end{matrix}

(3.34)

To obtain a second equation for

X

and

Y

let us consider

E [X_{t}^{2} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}]

On one hand

E [X_{t}^{2} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}] = E [X_{t}^{2} E [p_{n - 2} (X_{u}; u) | ℱ_{\leq t}] | ℱ_{\leq s}] = = E [X_{t}^{2} p_{n - 2} (X_{t}; t) | ℱ_{\leq s}] .

Setting

a_{k} (t) = b_{k} (t) = c_{k} (t) = 0

for

k < 0

, a repeated application of 2.19 gives

(3.35) x 2 p n − 2 ( x ; t ) = p n ( x ; t ) a n − 2 ( t ) a n − 1 ( t ) + p n − 1 ( x ; t ) a n − 2 ( t ) [ b n − 2 ( t ) + b n − 1 ( t ) ] + + p n − 2 ( x ; t ) [ a n − 2 ( t ) c n − 1 ( t ) + b n − 2 2 ( t ) + c n − 2 ( t ) a n − 3 ( t ) ] + + p n − 3 ( x ; t ) c n − 2 ( t ) [ b n − 3 ( t ) + b n − 2 ( t ) ] + p n − 4 ( x ; t ) c n − 3 ( t ) c n − 2 ( t ) ,

we obtain

(3.36) E [ X t 2 p n − 2 ( X u ; u ) | ℱ ≤ s ] = a n − 2 ( t ) a n − 1 ( t ) Y + + p n − 1 ( X s ; s ) a n − 2 ( t ) [ b n − 2 ( t ) + b n − 1 ( t ) ] + + p n − 2 ( X s ; s ) [ a n − 2 ( t ) c n − 1 ( t ) + b n − 2 2 ( t ) + c n − 2 ( t ) a n − 3 ( t ) ] + + p n − 3 ( X s ; s ) c n − 2 ( t ) [ b n − 3 ( t ) + b n − 2 ( t ) ] + + p n − 4 ( X s ; s ) c n − 3 ( t ) c n − 2 ( t ) .

On the other hand, one can rewrite

E [X_{t}^{2} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}]

E [X_{t}^{2} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}] = E [E [X_{t}^{2} | ℱ_{s, u}] p_{n - 2} (X_{u}; u) | ℱ_{\leq s}]

= A_{t, s, u} X_{s}^{2} p_{n - 2} (X_{s}; s) + B_{t, s, u} X_{s} E [X_{u} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}]

+ C_{t, s, u} E [X_{u}^{2} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}] + D_{t, s, u} X_{s} p_{n - 2} (X_{s}; s)

+ E_{t, s, u} E [X_{u} p_{n - 2} (X_{u}; u) | ℱ_{\leq s}] + F_{t, s, u} p_{n - 2} (X_{s}; s) .

After some algebra one gets

(3.37) E [ X t 2 p n − 2 ( X u ; u ) | ℱ ≤ s ] = X C t , s , u a n − 2 ( u ) a n − 1 ( u ) + p n ( X s ; s ) a n − 1 ( s ) R n ( t , s , u ) + p n − 1 ( X s ; s ) S n ( t , s , u ) + p n − 2 ( X s ; s ) T n ( t , s , u ) + + p n − 3 ( X s ; s ) U n ( t , s , u ) + p n − 4 ( X s ; s ) V n ( t , s , u ) ,

where

\begin{matrix} R_{n} (t, s, u) & = & A_{t, s, u} a_{n - 2} (s) + B_{t, s, u} a_{n - 2} (u), \end{matrix}

\begin{matrix} S_{n} (t, s, u) & = & A_{t, s, u} a_{n - 2} (s) [b_{n - 2} (s) + b_{n - 1} (s)] + B_{t, s, u} [a_{n - 2} (u) b_{n - 1} (s) + a_{n - 2} (s) b_{n - 2} (u)] \end{matrix}

\begin{matrix} + C_{t, s, u} a_{n - 2} (u) [b_{n - 2} (u) + b_{n - 1} (u)] + D_{t, s, u} b_{n - 2} (s) + E_{t, s, u} b_{n - 2} (u) + F_{t, s, u}, \end{matrix}

\begin{matrix} T_{n} (t, s, u) & = & A_{t, s, u} [a_{n - 2} (s) c_{n - 1} (s) + b_{n - 2}^{2} (s) + c_{n - 2} (s) a_{n - 3} (s)] \end{matrix}

\begin{matrix} + B_{t, s, u} [a_{n - 2} (u) c_{n - 1} (s) + b_{n - 2} (u) b_{n - 2} (s) + c_{n - 2} (u) a_{n - 3} (s)] \end{matrix}

\begin{matrix} + C_{t, s, u} [a_{n - 2} (u) c_{n - 1} (u) + b_{n - 2}^{2} (u) + c_{n - 2} (u) a_{n - 3} (u)] \end{matrix}

\begin{matrix} + D_{t, s, u} b_{n - 2} (s) + E_{t, s, u} b_{n - 2} (u) + F_{t, s, u}, \end{matrix}

\begin{matrix} U_{n} (t, s, u) & = & A_{t, s, u} c_{n - 2} (s) [b_{n - 3} (s) + b_{n - 2} (s)] + B_{t, s, u} [b_{n - 2} (u) c_{n - 2} (s) + c_{n - 2} (u) b_{n - 3} (s)] \end{matrix}

\begin{matrix} + C_{t, s, u} c_{n - 2} (u) [b_{n - 3} (u) + b_{n - 2} (u)] + D_{t, s, u} c_{n - 2} (s) + E_{t, s, u} c_{n - 2} (u), \end{matrix}

\begin{matrix} V_{n} (t, s, u) & = & A_{t, s, u} c_{n - 2} (s) c_{n - 3} (s) + B_{t, s, u} c_{n - 2} (u) c_{n - 3} (s) + C_{t, s, u} c_{n - 2} (u) c_{n - 3} (u) . \end{matrix}

Comparing the right hand side of ( 3.37 ) with the right hand side of 3.36 we get the second equation with unknowns

X

and

Y

(3.38) C t , s , u a n − 2 ( u ) a n − 1 ( u ) X − a n − 2 ( t ) a n − 1 ( t ) Y = − p n ( X s ; s ) R n ( t , s , u ) + p n − 1 ( X s ; s ) ( a n − 2 ( t ) [ b n − 2 ( t ) + b n − 1 ( t ) ] − S n ( t , s , u ) ) + p n − 2 ( X s ; s ) ( a n − 2 ( t ) c n − 1 ( t ) + b n − 2 2 ( t ) + c n − 2 ( t ) a n − 3 ( t ) − T n ( t , s , u ) ) + p n − 3 ( X s ; s ) ( c n − 2 ( t ) [ b n − 3 ( t ) + b n − 2 ( t ) ] − U n ( t , s , u ) ) + p n − 4 ( X s ; s ) ( c n − 2 ( t ) c n − 3 ( t ) − V n ( t , s , u ) ) .

A calculation based on ( 2.20 ) and ( 3.25 - 3.29 ) gives

a_{n - 2} (t) [b_{n - 2} (t) + b_{n - 1} (t)] - S_{n} (t, s, u) = - \frac{(t - s) (u - t)}{u (1 + σ s) + τ - q s} [θ σ α_{n - 1} + η β_{n - 1} - β_{n - 1} γ_{n - 1} + σ τ α_{n - 1} (γ_{n - 2} + γ_{n - 1}) - - σ α_{n - 1} δ_{n - 2} + σ β_{n - 1} (δ_{n - 2} + δ_{n - 1}) + q (β_{n - 1} γ_{n - 2} + σ α_{n - 1} δ_{n - 1})],

a_{n - 2} (t) c_{n - 1} (t) + b_{n - 2}^{2} (t) + c_{n - 2} (t) a_{n - 3} (t) - T_{n} (t, s, u) = - \frac{(t - s) (u - t)}{u (1 + σ s) + τ - q s} [1 + θ γ_{n - 2} + η δ_{n - 2} - γ_{n - 2} δ_{n - 2} - β_{n - 1} ε_{n - 1} + + τ (γ_{n - 2}^{2} + σ α_{n - 2} ε_{n - 2} + σ α_{n - 1} ɛ_{n - 1}) - σ α_{n - 2} φ_{n - 2} + + q (γ_{n - 2} δ_{n - 2} + β_{n - 2} ε_{n - 2} + σ α_{n - 1} φ_{n - 1}) + σ (δ_{n - 2}^{2} + β_{n - 2} φ_{n - 2} + β_{n - 1} φ_{n - 1})],

c_{n - 2} (t) [b_{n - 3} (t) + b_{n - 2} (t)] - U_{n} (t, s, u) = - \frac{(t - s) (u - t)}{u (1 + σ s) + τ - q s} [θ ɛ_{n - 2} + τ ɛ_{n - 2} (γ_{n - 3} + γ_{n - 2}) - δ_{n - 2} ɛ_{n - 2} + η φ_{n - 2} - φ_{n - 2} γ_{n - 3} + σ (δ_{n - 3} + δ_{n - 2}) φ_{n - 2} + q (δ_{n - 3} ɛ_{n - 2} + γ_{n - 2} φ_{n - 2})],

c_{n - 2} (t) c_{n - 3} (t) - V_{n} (t, s, u) = - \frac{(t - s) (u - t)}{u (1 + σ s) + τ - q s} [ɛ_{n - 3} (τ ɛ_{n - 2} - φ_{n - 2}) + φ_{n - 3} (q ɛ_{n - 2} + σ φ_{n - 2})],

Therefore, by 2.23 - 2.26 the coefficients at

p_{n - 1} (X_{s}; s), \dots, p_{n - 4} (X_{s}; s)

on the right hand side of 3.38 vanish, and 3.38 is equivalent to

(3.39) C t , s , u a n − 2 ( u ) a n − 1 ( u ) X − a n − 2 ( t ) a n − 1 ( t ) Y = − a n − 1 ( s ) p n ( X s ; s ) ( A t , s , u a n − 2 ( s ) + B t , s , u a n − 2 ( u ) ) .

Subtracting from ( 3.39 ) equation ( 3.34 ) multiplied by

a_{n - 2} (t)

, we get

(3.40) a n − 1 ( u ) ( C t , s , u a n − 2 ( u ) − b t , s , u a n − 2 ( t ) ) X = − a n − 1 ( s ) p n ( X s ; s ) ( A t , s , u a n − 2 ( s ) + B t , s , u a n − 2 ( u ) − a t , s , u a n − 2 ( t ) ) .

Now we notice that

C_{t, s, u} a_{n - 2} (u) - b_{t, s, u} a_{n - 2} (t) = = - \frac{(t - s) (u - t)}{(u - s) [u (1 + σ s) + τ - q s]} [(q s - τ) σ α_{n - 1} + (1 + s σ) β_{n - 1}]

see ( 4.2 ) and ( 2.15 ). Similarly,

A_{t, s, u} a_{n - 2} (s) + B_{t, s, u} a_{n - 2} (u) - a_{t, s, u} a_{n - 2} (t) = = \frac{(t - s) (u - t)}{(u - s) [u (1 + σ s) + τ - q s]} [(q u - τ) σ α_{n - 1} + (1 + u σ) β_{n - 1}]

Since ( 2.22 ) implies that

(3.41) a n − 1 ( u ) ( ( q s − τ ) σ α n − 1 + ( 1 + s σ ) β n − 1 ) = a n − 1 ( s ) ( ( q u − τ ) σ α n − 1 + ( 1 + u σ ) β n − 1 ) ,

and

a_{n - 1} (u) \neq 0

by assumption, therefore, ( 3.40 ) becomes

\begin{matrix} ((q s - τ) σ α_{n - 1} + (1 + s σ) β_{n - 1}) X = ((q s - τ) σ α_{n - 1} + (1 + s σ) β_{n - 1}) p_{n} (X_{s}; s) . \end{matrix}

(3.42)

It remains to verify that

\begin{matrix} ((q s - τ) σ α_{n - 1} + (1 + s σ) β_{n - 1}) \neq 0 . \end{matrix}

(3.43)

Since

β_{n - 1} = a_{n - 1} (0) \neq 0

by assumption, this is trivially true if

σ = 0

. Suppose

σ > 0

and ( 3.43 ) isn't true. Since

a_{n - 1} (u) \neq 0

from ( 3.41 ) we see that then the left hand side of ( 3.43 ) must vanish also when

s

is replaced by

u

. Thus

\begin{matrix} σ (α_{n - 1} q + β_{n - 1}) = 0, β_{n - 1} - σ τ α_{n - 1} = 0 . \end{matrix}

(3.44)

Since

β_{n - 1} \neq 0

, the determinant of the system of equations 3.44 is zero,

q + σ τ = 0

, and the second equation gives

σ τ \neq 0

. Since

α_{1} = 0

, from ( 2.22 ) with

q = - σ τ \neq 0

we get

α_{j} = β_{j}

for all

j \geq 1

. As

σ τ \neq 1

, this contradicts 3.44 . Thus ( 3.43 ) holds.

By ( 3.43 ), from ( 3.42 ) we get

X = p_{n} (X_{s}; s),

which means that

{p_{n} (x; t)}_{t > 0}

is a martingale for every

n \geq 0

Now we shall show that polynomials

{p_{n} (x; t)}_{n}

are orthogonal with respect to the distribution of

X_{t}

. Since the polynomials satisfy the three term recurrence 2.19 , it suffices to show that for all

k = 1, 2, \dots

\begin{matrix} E (p_{k} (X_{t}; t)) = 0 . \end{matrix}

(3.45)

Clearly, by the martingale property,

μ_{k} = E (p_{k} (X_{t}; t))

does not depend on

t

We proceed again by induction. From ( 4.2 ) we see that

a_{1} (t) = 1 + σ t

(1 + σ t) p_{2} (x; t) = x^{2} - b_{1} (t) p_{1} (x; t) - t

. Since

p_{1} (x; t) = x

, assumption ( 2.6 ) implies that 3.45 holds for

k = 1, 2

Assume then that for some

n \geq 3

we have

μ_{k} = 0

for all

1 \leq k \leq n - 1

. By 2.19

E (X_{t} p_{n} (X_{t}; t)) = a_{n} (t) E (p_{n + 1} (X_{t}; t)) = a_{n} (t) μ_{n + 1} .

By 2.7 , the martingale property of

{(p_{n} (x; t))}_{n}

, 2.19 and the induction assumption, the left hand side of the above equation transforms into

E (X_{t} p_{n} (X_{t}; t)) = E (\frac{t}{s} E [X_{s} | ℱ_{\geq t}] p_{n} (X_{t}; t)) = \frac{t}{s} E (X_{s} E [p_{n} (X_{t}; t) | ℱ_{\leq s}]) = = \frac{t}{s} E (X_{s} p_{n} (X_{s}; s)) = \frac{t}{s} a_{n} (s) E (p_{n + 1} (X_{s}; s)) = \frac{t}{s} a_{n} (s) μ_{n + 1},

(\frac{a_{n} (s)}{s} - \frac{a_{n} (t)}{t}) μ_{n + 1} = 0 .

Since

\frac{a_{n} (s)}{s} - \frac{a_{n} (t)}{t} = a_{n} (0) \frac{t - s}{t s} \neq 0,

we get

μ_{n + 1} = E (p_{n + 1} (X_{t}; t)) = 0

, and the orthogonality of

{p_{n} (x; t)}_{n}

is proved.

3.4 Proof of Theorem 2.5

Lemma 3.6. Suppose

X, Y

are square-integrable standardized, and there are constants

ɛ, A, B \geq 0, 0 < ρ < 1

such that

\begin{matrix} E [(X - ρ Y)^{2} | Y] & \leq & A + B | Y | + \frac{(1 - ρ^{2}) ɛ^{2}}{1 + ɛ^{2}} Y^{2}, \end{matrix}

(3.46)

\begin{matrix} E [(Y - ρ X)^{2} | X] & \leq & A + B | X | + \frac{(1 - ρ^{2}) ɛ^{2}}{1 + ɛ^{2}} X^{2} . \end{matrix}

(3.47)

2 ɛ < ρ < 1

, then there are constants

C_{1}, C_{2} < \infty

such that for

x > 0

we have

\begin{matrix} Pr (| X | > 2 x / ρ) + Pr (| Y | > 2 x / ρ) \end{matrix}

(3.48)

\begin{matrix} \leq (C_{1} / x^{2} + C_{2} / x + \frac{4 ɛ^{2}}{ρ^{2} + ɛ^{2}} \frac{1 + ρ}{1 - ρ}) (Pr (| X | > x) + Pr (| Y | > x)) . \end{matrix}

Proof. Enlarging

ɛ

if necessary, without loss of generality we may assume that

ɛ > 0

Let

N (x) = P (| X | \geq x) + P (| Y | \geq x)

K = 2 / ρ

. Throughout the proof,

C_{1}, C_{2}

denote positive constants which might differ at each occurrence.

The event

{| X | \geq K x}

, where

x > 0

is fixed, can be decomposed into the sum of two disjoint events

{| X | \geq K x} \cap {| Y | \geq x}

and

{| X | \geq K x} \cap {| Y | < x}

Therefore, denoting

\begin{matrix} P_{1} (x) & = & Pr (| X | \geq x, | Y | \geq x), \end{matrix}

\begin{matrix} P_{2} (x) & = & Pr (| X | \geq K x, | Y | < x), \end{matrix}

\begin{matrix} P_{3} (x) & = & Pr (| Y | \geq K x, | X | < x), \end{matrix}

we get

Pr (| X | \geq K x) \leq P_{1} (x) + P_{2} (x),

and hence by symmetry of the assumptions

\begin{matrix} N (K x) \leq 2 P_{1} (x) + P_{2} (x) + P_{3} (x) . \end{matrix}

(3.49)

To estimate

P_{1} (x)

, we observe the following.

Claim 3.7. If

{| X | \geq x

| Y | \geq x}

, and

0 < a < 1

, then either

(X - ρ Y)^{2} \geq a (1 - ρ)^{2} Y^{2} + (1 - ρ)^{2} (1 - a) x^{2},

(Y - ρ X)^{2} \geq a (1 - ρ)^{2} X^{2} + (1 - ρ)^{2} (1 - a) x^{2} .

Indeed, suppose that both inequalities fail, i.e. on the set

{| X | \geq x

| Y | \geq x}

we have

(X - ρ Y)^{2} < a (1 - ρ)^{2} X^{2} + (1 - ρ)^{2} (1 - a) x^{2}

and

(Y - ρ X)^{2} < a (1 - ρ)^{2} Y^{2} + (1 - ρ)^{2} (1 - a) x^{2}

. Adding the inequalities we obtain

(1 + ρ^{2}) (X^{2} + Y^{2}) - 4 ρ X Y < a (1 - ρ)^{2} (X^{2} + Y^{2}) + 2 (1 - ρ)^{2} (1 - a) x^{2} .

Since

4 ρ X Y \leq 2 ρ (X^{2} + Y^{2})

, this gives

(1 - ρ)^{2} (1 - a) (X^{2} + Y^{2}) < 2 (1 - ρ)^{2} (1 - a) x^{2} .

However, since

| X | \geq x, | Y | \geq x

, we have

(1 - ρ)^{2} (1 - a) (X^{2} + Y^{2}) \geq 2 x^{2} (1 - ρ)^{2} (1 - a),

a contradiction.

Claim 3.7 with

a = 1 / 2

implies

P_{1} (x) \leq Pr (| X - ρ Y | \geq (1 - ρ) \sqrt{Y^{2} + x^{2}} / \sqrt{2}, | Y | \geq x)

+ Pr (| Y - ρ X | \geq (1 - ρ) \sqrt{X^{2} + x^{2}} / \sqrt{2}, | X | \geq x) .

From conditional Chebyshev's inequality and ( 3.46 ) we get

Pr (| X - ρ Y | \geq (1 - ρ) \sqrt{Y^{2} + x^{2}} / \sqrt{2}, | Y | \geq x)

\leq 2 \int_{| Y | \geq x} \frac{A + B | Y | + (1 - ρ^{2}) ɛ^{2} Y^{2} / (1 + ɛ^{2})}{(1 - ρ)^{2} (Y^{2} + x^{2})} d P

\leq \int_{| Y | \geq x} \frac{2 A}{(1 - ρ)^{2} x^{2}} d P + \int_{| Y | \geq x} \frac{2 B | Y |}{(1 - ρ)^{2} | Y | x} d P + \int_{| Y | \geq x} \frac{(1 - ρ^{2}) 2 ɛ^{2} Y^{2} / (1 + ɛ^{2})}{(1 - ρ)^{2} Y^{2}} d P

\leq \frac{C_{1}}{x^{2}} Pr (| Y | \geq x) + \frac{C_{2}}{x} Pr (| Y | \geq x) + \frac{2 ɛ^{2}}{1 + ɛ^{2}} \frac{1 + ρ}{1 - ρ} Pr (| Y | \geq x),

where

C_{1} = 2 A / (1 - ρ)^{2}

C_{2} = 2 B / (1 - ρ)^{2}

. Since the assumptions are symmetric in

X, Y

, this shows that there are constants

C_{1}, C_{2} < \infty

such that

\begin{matrix} P_{1} (x) \leq \frac{C_{1} N (x)}{x^{2}} + \frac{C_{2} N (x)}{x} + \frac{2 ɛ^{2}}{1 + ɛ^{2}} \frac{1 + ρ}{1 - ρ} N (x) . \end{matrix}

(3.50)

To estimate

P_{2} (x)

, we use the trivial estimate

| Y - ρ X | \geq ρ | X | - | Y |

, which shows that the event

{| X | \geq K x, | Y | < x}

implies

| Y - ρ X | \geq ɛ | X | + ((ρ - ɛ) K - 1) x = ɛ | X | + \frac{ρ - 2 ɛ}{ρ} x

Therefore

\begin{matrix} P_{2} (x) \leq & Pr (| Y - ρ X | \geq ɛ | X | + \frac{ρ - 2 ɛ}{ρ} x, | X | \geq K x) \end{matrix}

\begin{matrix} = \int_{| X | \geq K x} Pr (| Y - ρ X | \geq ɛ | X | + \frac{ρ - 2 ɛ}{ρ} x | X) d P \end{matrix}

\begin{matrix} \leq \int_{| X | \geq K x} \frac{A + B | X | + (1 - ρ^{2}) ɛ^{2} X^{2} / (1 + ɛ^{2})}{(ɛ | X | + (ρ - 2 ɛ) x / ρ)^{2}} d P \end{matrix}

\begin{matrix} \leq \int_{| X | \geq K x} \frac{A}{((ρ - 2 ɛ) x / ρ)^{2}} d P + \int_{| X | \geq K x} \frac{B | X |}{ɛ | X | (ρ - 2 ɛ) x / ρ} d P \end{matrix}

\begin{matrix} + \int_{| X | \geq K x} \frac{(1 - ρ^{2}) ɛ^{2} X^{2}}{(1 + ɛ^{2}) ɛ^{2} X^{2}} d P . \end{matrix}

\begin{matrix}  \end{matrix}

Since

N (K x) \leq N (x)

, this shows that

\begin{matrix} P_{2} (x) \leq \frac{C_{1} N (x)}{x^{2}} + \frac{C_{2} N (x)}{x} + \frac{1 - ρ^{2}}{1 + ɛ^{2}} Pr (| X | \geq K x) . \end{matrix}

(3.51)

where

C_{1} = ρ^{2} A / (ρ - 2 ɛ)^{2}

C_{2} = ρ B / (ɛ (ρ - 2 ɛ))

. By symmetry of the assumptions, we also have

\begin{matrix} P_{3} (x) \leq \frac{C_{1} N (x)}{x^{2}} + \frac{C_{2} N (x)}{x} + \frac{1 - ρ^{2}}{1 + ɛ^{2}} Pr (| Y | \geq K x) . \end{matrix}

(3.52)

Combining ( 3.49 ), ( 3.50 ), ( 3.51 ), and ( 3.52 ) we obtain that there are constants

C_{3}, C_{4} > 0

such that

N (K x) \leq \frac{C_{3} N (x)}{x^{2}} + \frac{C_{4} N (x)}{x} + \frac{1 - ρ^{2}}{1 + ɛ^{2}} N (K x) + \frac{4 ɛ^{2}}{1 + ɛ^{2}} \frac{1 + ρ}{1 - ρ} N (x),

which implies ( 3.48 ). □

Lemma 3.8. Suppose

X, Y

are square-integrable standardized,

E (| X |^{p}) + E (| Y |^{p}) < \infty

for some

p > 1

and the assumptions of Lemma 3.6 are satisfied with constants

ɛ > 0, 0 < ρ < 1

such that

\begin{matrix} 2^{p + 3} ɛ^{2} \frac{1 + ρ}{1 - ρ} \leq ρ^{p + 3} . \end{matrix}

(3.53)

Then

E (| X |^{p + 1}) + E (| Y |^{p + 1}) < \infty

Proof. Clearly, ( 3.53 ) implies $2 ɛ < ρ$ . Indeed, $2^{2} ɛ^{2} < 2^{p + 3} ɛ^{2} \frac{1 + ρ}{1 - ρ} \leq ρ^{p + 3} < ρ^{2} .$ We use Lemma 3.6 and we use the notation $N (x) = P (| X | \geq x) + P (| Y | \geq x)$ , $K = 2 / ρ$ introduced in its proof. Fix $M > 0$ . Then, noticing that $K > 1$ from Lemma 3.6 we get $\int_{0}^{M} (p + 1) x^{p} N (x) d x = K^{p + 1} \int_{0}^{M / K} (p + 1) x^{p} N (K x) d x$ $\leq \frac{p + 1}{p - 1} C_{1} \int_{0}^{\infty} (p - 1) x^{p - 2} N (x) d x + \frac{p + 1}{p} C_{2} \int_{0}^{\infty} p x^{p - 1} N (x) d x$ $+ \frac{4 ɛ^{2}}{ρ^{2} + ɛ^{2}} \frac{1 + ρ}{1 - ρ} K^{p + 1} \int_{0}^{M} (p + 1) x^{p} N (x) d x$ $< \frac{p + 1}{p - 1} C_{1} (E (| X |^{p - 1}) + E (| Y |^{p - 1})) + \frac{p + 1}{p} C_{2} (E (| X |^{p}) + E (| Y |^{p}))$ $+ ɛ^{2} 2^{p + 3} / ρ^{p + 3} \frac{1 + ρ}{1 - ρ} \int_{0}^{M} (p + 1) x^{p} N (x) d x .$ Therefore, if ( 3.53 ) holds true, then ${sup}_{M > 0} \int_{0}^{M} (p + 1) x^{p} N (x) d x < \infty$ which implies $E (| X |^{p + 1}) + E (| Y |^{p + 1}) = \int_{0}^{\infty} (p + 1) x^{p} N (x) d x < \infty$ . □

We will also need the following simple observation.

Lemma 3.9. Fix

r \geq 1

. Under the assumptions of Theorem 2.5 , if

E (| X_{t} |^{r}) < \infty

for some

t > 0

, then

E (| X_{t} |^{r}) < \infty

for all

t > 0

Proof. If $E (| X_{t} |^{r}) < \infty$ for some $t > 0$ then $E (| X_{s} |^{r}) < \infty$ for all $s < t$ by ( 2.8 ).
Similarly, $E (| X_{u} |^{r}) < \infty$ for all $u > t$ by ( 2.7 ). □

Conclusion of proof of Theorem 2.5 . Let $r > 2$ be fixed. Since our proof relies only on inequalities, if $σ τ = 0$ we can increase slightly the values of parameters $σ, τ$ on the right hand sides of ( 2.28 ) and ( 2.27 ) to ensure $σ^{'}, τ^{'} > 0$ and $σ^{'} τ^{'} \leq 2^{- 4 r - 10}$ , replacing both equalities by the appropriate inequalities; thus without loss of generality we may assume that $σ τ > 0$ .
Let $p = r - 1 > 1$ . Define $ɛ^{2} = \sqrt{σ τ} / 2$ , $t_{0} = τ / ɛ^{2}, s_{0} = ɛ^{2} / σ$ . Then the correlation coefficient between $X_{s_{0}}$ and $X_{t_{0}}$ is $ρ = \sqrt{s_{0} / t_{0}} = ɛ^{2} / \sqrt{σ τ} = 1 / 2$ . From the inequality form of ( 2.28 ) and ( 2.27 ) we deduce that random variables $X = X_{s_{0}} / \sqrt{s_{0}}$ and $Y = X_{t_{0}} / \sqrt{t_{0}}$ satisfy ( 3.46 ) and ( 3.47 ) with $A = \frac{1 - ρ^{2}}{1 + ɛ^{2}}, B = \frac{1 - ρ^{2}}{1 + ɛ^{2}} max {| η | ɛ / \sqrt{σ}, | θ | \sqrt{τ} / ɛ} .$ Condition ( 3.53 ) is satisfied since with $r = p + 1$ we have $(σ τ)^{1 / 2} \leq 1 / 2^{2 p + 7}$ . Thus $2^{p + 3} ɛ^{2} \frac{1 + ρ}{1 - ρ} = 2^{p + 3} 3 ɛ^{2} < 2^{p + 4} (σ τ)^{1 / 2} \leq 1 / 2^{p + 3} = ρ^{p + 3} .$ Therefore, by Lemma 3.8 , $E (| X_{t_{0}} |^{r}) = E (| X_{t_{0}} |^{p + 1}) < \infty$ and therefore by Lemma 3.9 , $E (| X_{t} |^{r}) < \infty$ for all $t \geq 0$ . □

4 Examples of explicit recurrences

In this section we will say that

(X_{t})

is a quadratic harness with parameters

q, η, θ, σ, τ

if ( 2.3 ), ( 2.6 ), and ( 2.14 ) hold. Our goal is to derive explicit versions of ( 2.19 ) at the expense of additional assumptions on the parameters.

Theorem 4.1. Suppose

(X_{t})

is a quadratic harness with parameters

0 \leq σ τ < 1

- 1 < q \leq 1 - 2 \sqrt{σ τ}

. Moreover, assume that for each

t > 0

random variable

X_{t}

has moments of all orders and infinite support.

Then

(X_{t})

has orthogonal martingale polynomials which are given by recurrence ( 2.19 ) with

p_{0} (x, t) = 1

p_{1} (x, t) = x

, where the coefficients of the recurrence are linear functions of

t

\begin{matrix} a_{n} (t) = σ α_{n + 1} t + β_{n + 1}, b_{n} (t) = γ_{n} t + δ_{n}, c_{n} (t) = (β_{n} t + τ α_{n}) ω_{n}, \end{matrix}

(4.1)

which are determined as follows.

(i) The initial conditions are $α_{1} = 0, β_{1} = 1, γ_{0} = δ_{0} = 0, ω_{1} = 1 .$

(ii) Sequences

(α_{n}), (β_{n})

satisfy

\begin{matrix} [\begin{matrix} α_{n + 1} \\ β_{n + 1} \end{matrix}] = [\begin{matrix} q & 1 \\ - σ τ & 1 \end{matrix}] \times [\begin{matrix} α_{n} \\ β_{n} \end{matrix}], n \geq 1 . \end{matrix}

(4.2)

Moreover, denoting

λ_{n, k} = β_{n} β_{n + k} - σ τ α_{n} α_{n + k}

we have

λ_{n, k} > 0

for all

n \geq 1

k \geq 0

(iii) Setting

λ_{0, 2} = α_{0} = β_{0} = 0

, sequences

(γ_{n}), (δ_{n})

satisfy the following system of linear recurrences

\begin{matrix} γ_{n + 1} & = & \frac{q + σ τ}{λ_{n + 2, 0}} (λ_{n, 2} γ_{n} + (α_{n + 2} β_{n} - β_{n + 2} α_{n}) σ δ_{n}) \end{matrix}

\begin{matrix} + \frac{σ α_{n + 2}}{λ_{n + 2, 0}} (η τ α_{n + 1} + θ β_{n + 1}) + \frac{β_{n + 2}}{λ_{n + 2, 0}} (θ σ α_{n + 1} + η β_{n + 1}), \end{matrix}

(4.3)

\begin{matrix} δ_{n + 1} & = & \frac{q + σ τ}{λ_{n + 2, 0}} (λ_{n, 2} δ_{n} + (α_{n + 2} β_{n} - β_{n + 2} α_{n}) τ γ_{n}) \end{matrix}

\begin{matrix} + \frac{β_{n + 2}}{λ_{n + 2, 0}} (η τ α_{n + 1} + θ β_{n + 1}) + \frac{τ α_{n + 2}}{λ_{n + 2, 0}} (θ σ α_{n + 1} + η β_{n + 1}), n \geq 0 . \end{matrix}

(iv) Sequence

(ω_{n})

satisfies the linear recurrence

\begin{matrix} ω_{n + 1} & = & \frac{(q + σ τ) λ_{n - 1, 1}}{λ_{n + 1, 1}} ω_{n} \end{matrix}

(4.4)

\begin{matrix} + \frac{1 + γ_{n} (τ γ_{n} - δ_{n} + θ) + δ_{n} (q γ_{n} + σ δ_{n} + η)}{λ_{n + 1, 1}}, n \geq 2 \end{matrix}

with the initial term

ω_{2} = (1 + q) \frac{(1 - σ τ)^{2} + (η + θ σ) (θ + η τ)}{(1 - σ τ)^{2} (1 - σ τ (2 + q))} .

Remark 4.1. If

- 1 \leq q < 1 - 2 \sqrt{σ τ}

, then the explicit solution of recurrence ( 4.2 ) is

\begin{matrix} α_{n} = \frac{λ_{+}^{n - 1} - λ_{-}^{n - 1}}{λ_{+} - λ_{-}}, β_{n} = \frac{λ_{+}^{n - 1} (1 - λ_{-}) + λ_{-}^{n - 1} (λ_{+} - 1)}{λ_{+} - λ_{-}}, \end{matrix}

(4.5)

where

λ_{\pm} = \frac{1}{2} (1 + q \pm \sqrt{(1 + q)^{2} - 4 (q + σ τ)}) .

The proof of Theorem 4.1 relies on Theorem 2.4 and the following lemma.

Lemma 4.2. If

σ τ < 1

and

- 1 \leq q \leq 1 - 2 \sqrt{σ τ}

, then the solutions of ( 4.2 ) satisfy

β_{n} > \sqrt{σ τ} α_{n} \geq 0 .

In particular,

β_{n} β_{n + 1} > σ τ α_{n} α_{n + 1}, β_{n}^{2} > σ τ α_{n}^{2}, for all n \geq 1 .

Proof. From the explicit solution in Remark 4.1 , we see that $| λ_{-} | \leq | λ_{+} |$ and hence $α_{n} \geq 0$ for $n > 1$ . We now prove by induction that $β_{n} > \sqrt{σ τ} α_{n}$ for all $n$ . This is trivially true for $n = 1$ , see ( 2.21 ). Suppose the inequality is satisfied for some $n \geq 1$ . Then ( 4.2 ) implies that $β_{n + 1} - \sqrt{σ τ} α_{n + 1} = (1 - \sqrt{σ τ}) β_{n} - \sqrt{σ τ} (q + \sqrt{σ τ}) α_{n} \geq (1 - \sqrt{σ τ}) (β_{n} - \sqrt{σ τ} α_{n}) > 0$ . □

Proof of Theorem 4.1 . By Lemma 4.2 ,

a_{n} (t) = σ α_{n + 1} t + β_{n + 1} \geq β_{n + 1} > 0

for all

t \geq 0

n \geq 0

. Moreover, recurrences ( 4.3 ), ( 4.4 ) are well defined and have unique solution. Therefore, to end the proof we need only to verify that with

\begin{matrix} ɛ_{n} = ω_{n} β_{n}, φ_{n} = τ ω_{n} α_{n}, \end{matrix}

(4.6)

the equations ( 2.21 - 2.26 ) are satisfied. Trivially, ( 2.21 ) holds true. It is easy to verify that the solution of ( 4.2 ) satisfies ( 2.22 ). Using ( 4.6 ) from ( 4.2 ) we also get ( 2.26 ).

It remains to show that equations ( 2.23 ), ( 2.24 ), and ( 2.25 ) hold. We will rewrite these equations using ( 4.2 ) and its equivalent form

\begin{matrix} q β_{n + 1} + σ τ α_{n + 1} = (q + σ τ) β_{n}, β_{n + 1} - α_{n + 1} = - (q + σ τ) α_{n}, n \geq 1 . \end{matrix}

(4.7)

Equation ( 2.23 ) is equivalent to

γ_{n + 1} (β_{n + 1} - σ τ α_{n + 1}) - δ_{n + 1} σ (β_{n + 1} + α_{n + 1} q)

= (σ τ α_{n + 1} + β_{n + 1} q) γ_{n} + δ_{n} σ (β_{n + 1} - α_{n + 1}) + σ α_{n + 1} θ + β_{n + 1} η .

Using ( 4.2 ) and ( 4.7 ), we can rewrite it as

(4.8) γ n + 1 β n + 2 − σ δ n + 1 α n + 2 = ( q + σ τ ) β n γ n − ( q + σ τ ) σ α n δ n + σ α n + 1 θ + β n + 1 η .

Similarly, ( 2.25 ) for our sequences rewrites as

(4.9) ω n ( δ n + 1 ( β n + 1 − σ τ α n + 1 ) − τ γ n + 1 ( β n + 1 + q α n + 1 ) ) = ω n ( ( q β n + 1 + σ τ α n + 1 ) δ n + τ ( β n + 1 − α n + 1 ) γ n + θ β n + 1 + η τ α n + 1 ) .

Both of these equations are satisfied; in fact, ( 4.3 ) was obtained by solving the system of equations ( 4.8 - 4.9 ) when

ω_{n} > 0

To verify that ( 2.24 ) holds, we substitute ( 4.6 ); using ( 4.2 ) and Lemma 4.2 we get

ω_{n + 1} = \frac{β_{n} (q β_{n} + σ τ α_{n}) + σ τ α_{n} (β_{n} - α_{n})}{β_{n + 1} β_{n + 2} - σ τ α_{n + 1} α_{n + 2}} ω_{n} + \frac{1 + γ_{n} (τ γ_{n} - δ_{n} + θ) + δ_{n} (q γ_{n} + σ δ_{n} + η)}{β_{n + 1} β_{n + 2} - σ τ α_{n + 1} α_{n + 2}} .

Since ( 4.7 ) holds, thus ( 2.24 ) is satisfied, too. □

Remark 4.2. From ( 4.3 ) we get

\begin{matrix} γ_{1} = \frac{η + θ σ}{1 - σ τ}, δ_{1} = \frac{η τ + θ}{1 - σ τ}, \end{matrix}

(4.10)

which implies that the initial recurrences are

\begin{matrix} x p_{0} (x; t) & = p_{1} (x; t), \end{matrix}

(4.11)

\begin{matrix} x p_{1} (x; t) & = (1 + σ t) p_{2} (x; t) + \frac{(η + θ σ) t + (η τ + θ)}{1 - σ τ} p_{1} (x; t) + t p_{0} (x; t) . \end{matrix}

(4.12)

\begin{matrix}  \end{matrix}

Thus the first three orthogonal martingale polynomials are

(4.13) p 0 ( x ; t ) = 1 , p 1 ( x ; t ) = x , p 2 ( x ; t ) = 1 1 + σ t x 2 − ( η + θ σ ) t + η τ + θ ( 1 − σ τ ) ( 1 + σ t ) x − t 1 + σ t .

4.1 Free quadratic harnesses

Free harnesses have parameter

\begin{matrix} q = - σ τ . \end{matrix}

(4.14)

The adjective ”free” is motivated by the fact that when

σ τ = 0

this choice is related to free convolutions that arise in free probability, see [12,Section4.3] and [11,Section4] . In general, from Theorem 4.1 it is clear that this choice of

q

significantly simplifies the recurrences. It is not obvious whether this case is further related to free convolutions, see however [22] .

Assuming

σ τ < 1

, ( 4.14 ) implies

- 1 < q < 1 - 2 \sqrt{σ τ}

, so it follows from Theorem 4.1 that the orthogonal martingale polynomials exist. It is easy to check that the solution of ( 4.2 ) is

\begin{matrix} α_{n} = (1 - σ τ)^{n - 2}, β_{n} = (1 - σ τ)^{n - 2}, n \geq 2 . \end{matrix}

(4.15)

For our choice of

q

and

n \geq 0

we have

λ_{n + 2, 0} = β_{n + 2}^{2} - σ τ α_{n + 2}^{2} = (1 - σ τ)^{2 n + 1} .

Using this identity, from ( 4.3 ) we get for

n \geq 2

γ_{n} = \frac{η + 2 θ σ + η σ τ}{{(1 - σ τ)}^{2}},

δ_{n} = \frac{θ + 2 η τ + θ σ τ}{{(1 - σ τ)}^{2}} .

Similarly,

λ_{n + 1, 1} = β_{n + 1} β_{n + 2} - σ τ α_{n + 1} α_{n + 2} = (1 - σ τ)^{2 n}, n \geq 0

Substituting this and

δ_{n}, γ_{n}

into ( 4.4 ), we get

ω_{n} = \frac{(1 - σ τ)^{2} + (η + θ σ) (θ + η τ)}{{(1 - σ τ)}^{2 n}}, n \geq 3 .

We also compute

ω_{2} = \frac{(1 - σ τ)^{2} + (η + θ σ) (θ + η τ)}{{(1 - σ τ)}^{3}} .

Therefore recurrence ( 2.19 ) with initial values ( 2.21 ) gives

x p_{2} (x) = (1 - σ τ) (1 + σ t) p_{3} (x) + \frac{(η + 2 θ σ + η σ τ) t + θ + 2 η τ + θ σ τ}{(1 - σ τ)^{2}} p_{2} (x)

+ \frac{{(1 - σ τ)}^{2} + (η + θ σ) (θ + η τ)}{(1 - σ τ)^{3}} (t + τ) p_{1} (x),

x p_{n} (x) = (1 - σ τ)^{n - 1} (1 + σ t) p_{n + 1} (x) + \frac{(η + 2 θ σ + η σ τ) t + θ + 2 η τ + θ σ τ}{{(1 - σ τ)}^{2}} p_{n} (x)

+ \frac{{(1 - σ τ)}^{2} + (η + θ σ) (θ + η τ)}{(1 - σ τ)^{n + 2}} (t + τ) p_{n - 1} (x)

for

n \geq 2

After renormalizing the

n

-th polynomial by

(1 - σ τ)^{(n - 2) (n - 1) / 2}

, from Theorem 4.1 we get the following.

Proposition 4.3 (free quadratic harnesses). Suppose

(X_{t})

is a quadratic harness with parameters such that

(1 - σ τ)^{2} + (η + θ σ) (θ + η τ) > 0

0 \leq σ τ < 1

, and

q = - σ τ

. If for

t > 0

random variable

X_{t}

has all moments and infinite support, then it has orthogonal martingale polynomials given by the three step recurrences ( 4.11 ), ( 4.12 ), and

\begin{matrix} x p_{2} (x; t) & = & (1 + σ t) p_{3} (x; t) + \frac{(η + 2 θ σ + η σ τ) t + θ + 2 η τ + θ σ τ}{{(1 - σ τ)}^{2}} p_{2} (x; t) \end{matrix}

\begin{matrix} + \frac{{(1 - σ τ)}^{2} + + (η + θ σ) (θ + η τ)}{(1 - σ τ)^{3}} (t + τ) p_{1} (x; t), \end{matrix}

\begin{matrix} x p_{n} (x; t) & = & (1 + σ t) p_{n + 1} (x; t) + \frac{(η + 2 θ σ + η σ τ) t + θ + 2 η τ + θ σ τ}{(1 - σ τ)^{2}} p_{n} (x; t) \end{matrix}

\begin{matrix} + \frac{{(1 - σ τ)}^{2} + (η + θ σ) (θ + η τ)}{(1 - σ τ)^{4}} (t + τ) p_{n - 1} (x; t), n \geq 2 . \end{matrix}

Remark 4.3. The recurrence in Proposition 4.3 is a finite perturbation of the constant coefficient recurrence which was analyzed by many authors, see [29] and the references therein, see also [27] .

Remark 4.4. In [11] we show that the free bi-Poisson process is associated with the generalized free convolution studied in [7] . It is interesting to ask if an analogous situation occurs for the general free harnesses of Proposition 4.3 ; for recent extensions of generalized free convolutions to perturbations of higher-order terms, see [22] .

4.2 Classical quadratic harnesses

The classical quadratic harnesses have parameter

q = 1 - 2 \sqrt{σ τ}

. The adjective ”classical” is motivated by the fact that when

σ τ = 0

, quadratic harnesses with

q = 1

are related to classical probability, see [12,Section4.2] .

Proposition 4.4 (Classical quadratic harnesses). Suppose

(X_{t})

is a quadratic harness with parameters such that

(1 - σ τ)^{2} + (η + θ σ) (θ + η τ) > 0

0 \leq σ τ < 1

, and

q = 1 - 2 \sqrt{σ τ}

. If for

t > 0

random variable

X_{t}

has all moments and infinite support, then it has orthogonal martingale polynomials given by the three step recurrences ( 4.11 ), ( 4.12 ), and

\begin{matrix} x p_{n} (x; t) & = & (1 + (n - 1) ρ + n σ t) p_{n + 1} (x; t) + (γ_{n} t + δ_{n}) p_{n} (x; t) \end{matrix}

(4.16)

\begin{matrix} + ω_{n}^{\circ} ((1 + (n - 2) ρ) t + τ (n - 1)) p_{n - 1} (x; t), n \geq 2 . \end{matrix}

(4.17)

where

ρ = \sqrt{σ τ}

(4.18) γ n = n ( 1 + ( n − 2 ) ρ ) ( 1 − ρ ) 2 ( 1 + ( 2 n − 1 ) ρ ) ( 1 + ( 2 n − 3 ) ρ ) × ( η + ( 2 n − 1 ) θ σ + ( 2 n − 3 ) η ρ + 2 ( n − 1 ) 2 η ρ 2 + ( 2 ( n − 1 ) 2 − 1 ) θ σ ρ ) , n ≥ 1

(4.19) δ n = n ( 1 + ( n − 2 ) ρ ) ( 1 − ρ ) 2 ( 1 + ( 2 n − 1 ) ρ ) ( 1 + ( 2 n − 3 ) ρ ) × ( θ + ( 2 n − 1 ) η τ + ( 2 n − 3 ) θ ρ + 2 ( n − 1 ) 2 θ ρ 2 + ( 2 ( n − 1 ) 2 − 1 ) η τ ρ ) , n ≥ 1

(4.20) ω n ∘ = n ( 1 + ( n − 3 ) ρ ) ( 1 − ρ ) 2 ( 1 + ( 2 n − 2 ) ρ ) ( 1 + ( 2 n − 4 ) ρ ) + n ( n − 1 ) ( 1 + ( n − 2 ) ρ ) ( 1 + ( n − 3 ) ρ ) ( 1 − ρ ) 4 ( 1 + ( 2 n − 2 ) ρ ) ( 1 + ( 2 n − 3 ) ρ ) 2 ( 1 + ( 2 n − 4 ) ρ ) × ( ( 1 + ( n − 2 ) ρ ) θ + ( n − 1 ) η τ ) ( ( 1 + ( n − 2 ) ρ ) η + ( n − 1 ) θ σ ) , n ≥ 2 .

Proof. The assumptions of Theorem 4.1 are satisfied, so the orthogonal martingale polynomials exist, and we only need to solve the recurrences ( 4.2 ), ( 4.3 ), ( 4.4 ), and renormalize the polynomials to simplify the final three step recurrence ( 4.17 ).

From ( 4.2 ) we get

α_{n} = (1 - ρ)^{n - 2} (n - 1), n \geq 1,

β_{n} = (1 - ρ)^{n - 2} (1 + (n - 2) ρ), n \geq 2 .

Thus

λ_{1, k} = β_{k + 1}

k \geq 0

, and a calculation gives

\begin{matrix} λ_{n, k} = (1 - ρ)^{2 n + k - 3} (1 + (2 n + k - 3) ρ), n \geq 2 \end{matrix}

(4.21)

α_{n + 2} β_{n} - β_{n + 2} α_{n} = 2 (1 - ρ)^{2 n - 1}, n \geq 1 .

Since

q + σ τ = (1 - ρ)^{2}

, equations ( 4.3 ) simplify to

(4.22) γ n + 1 = 1 + ( 2 n − 1 ) ρ 1 + ( 2 n + 1 ) ρ γ n + 2 σ 1 + ( 2 n + 1 ) ρ δ n + ( 1 + ( 2 n − 1 ) ρ + 2 n 2 ρ 2 ) η + ( 2 n + 1 + ( 2 n 2 − 1 ) ρ ) θ σ ( 1 − ρ ) 2 ( 1 + ( 2 n + 1 ) ρ ) ,

(4.23) δ n + 1 = 1 + ( 2 n − 1 ) ρ 1 + ( 2 n + 1 ) ρ δ n + 2 τ 1 + ( 2 n + 1 ) ρ γ n + ( 1 + ( 2 n − 1 ) ρ + 2 n 2 ρ 2 ) θ + ( 2 n + 1 + ( 2 n 2 − 1 ) ρ ) η τ ( 1 − ρ ) 2 ( 1 + ( 2 n + 1 ) ρ ) .

Formula ( 4.10 ) shows that ( 4.18 ) and ( 4.19 ) hold true for

n = 1

. Assuming that ( 4.18 ) and ( 4.19 ) hold true for some

n \geq 1

, from equations ( 4.22 ) and ( 4.23 ) a computer-assisted calculation shows that the formulas hold true for

n + 1

as well.

We now show that the solution of ( 4.4 ) is

(4.24) ω n = n ( 1 + ( n − 3 ) ρ ) ( 1 − ρ ) 2 n − 2 ( 1 + ( 2 n − 2 ) ρ ) ( 1 + ( 2 n − 4 ) ρ ) + n ( n − 1 ) ( 1 + ( n − 2 ) ρ ) ( 1 + ( n − 3 ) ρ ) ( 1 − ρ ) 2 n ( 1 + ( 2 n − 2 ) ρ ) ( 1 + ( 2 n − 3 ) ρ ) 2 ( 1 + ( 2 n − 4 ) ρ ) × ( ( 1 + ( n − 2 ) ρ ) θ + ( n − 1 ) η τ ) ( ( 1 + ( n − 2 ) ρ ) η + ( n − 1 ) θ σ ) , n ≥ 1 .

Indeed, a calculation shows that the formula holds true for

n = 1, 2

. Suppose ( 4.24 ) holds for some

n \geq 2

. For

k = 1

, formula ( 4.21 ) holds also for

n = 1

; thus for

n \geq 2

recurrence ( 4.4 ) simplifies to

ω_{n + 1} = \frac{1 + (2 n - 4) ρ}{(1 - ρ)^{2} (1 + 2 n ρ)} ω_{n} + \frac{1 + τ γ_{n}^{2} + σ δ_{n}^{2} + θ γ_{n} + η δ_{n} - 2 ρ δ_{n} γ_{n}}{(1 - ρ)^{2 n} (1 + 2 n ρ)} .

Using ( 4.22 ), ( 4.23 ) and ( 4.24 ), a computer assisted calculation verifies that ( 4.24 ) holds true for

n + 1

Renormalizing the

n

-th polynomial in ( 2.19 ) by the factor

(1 - ρ)^{(n - 2) (n - 1) / 2}

, we get ( 4.17 ) with

ω_{n}

replaced by

ω_{n}^{\circ}

. □

4.3 Orthogonal martingale polynomials when $σ τ = 0$

We use the standard

q

-notation

\begin{matrix} [n]_{q} & = & 1 + q + \cdot \cdot \cdot + q^{n - 1}, \end{matrix}

\begin{matrix} [n]_{q}! & = & [1]_{q} [2]_{q} . . . [n]_{q}, \end{matrix}

with the usual conventions

[0]_{q} = 0, [0]_{q}! = 1

. In this notation, Remark 4.1 gives

α_{n} = [n - 1]_{q}

, and

β_{n} = 1

Passing to the time-inverse

(t X_{1 / t})

if necessary, without loss of generality we may assume that

σ = 0

. In this case, recurrences in Theorem 4.1 have explicit solutions and orthogonal martingale polynomials are monic.

Theorem 4.5. Suppose

(X_{t})

is a quadratic harness with covariance ( 2.6 ) and parameters such that

σ = 0

- 1 < q \leq 1

, and

1 + [n]_{q} η θ + [n]_{q}^{2} τ η^{2} > 0

for all

n

. If for each

t > 0

random variable

X_{t}

has infinite support, then the monic orthogonal martingale polynomials

p_{n} (x; t)

are given by the recurrence

(4.25) x p n ( x ; t ) = p n + 1 ( x ; t ) + ( η t + θ + ( [ n ] q + [ n − 1 ] q ) η τ ) [ n ] q p n ( x ; t ) + ( t + τ [ n − 1 ] q ) ( 1 + [ n − 1 ] q η θ + [ n − 1 ] q 2 τ η 2 ) [ n ] q p n − 1 ( x ; t ) , n ≥ 1 ,

with the initial condition

p_{0} = 1, p_{1} (x) = x

Remark 4.5. It is easy to give simple sufficient conditions on the parameters which imply that

ω_{n} : = 1 + [n - 1]_{q} η θ + [n - 1]_{q}^{2} τ η^{2} > 0

for all

n \geq 1

. Suppose that

- 1 < q \leq 1

(i) If $η θ \geq 0$ then, trivially, $ω_{n} > 0$ .
(ii) If $τ = 0$ , then $ω_{n} > 0$ provided $1 + η θ > max {q, 0}$ . Indeed, for $n \geq 2$ we have
$\begin{matrix} 1 + q \leq [n]_{q} \leq 1 & if - 1 < q < 0, \end{matrix}$ (4.26)

$\begin{matrix} 1 \leq [n]_{q} < 1 / (1 - q) & if 0 \leq q < 1 . \end{matrix}$ (4.27)
For $η θ < 0$ we get $ω_{n} > 0$ from the right hand sides of these equations.
(iii) if $- 1 < q < 1$ and $θ^{2} < 4 τ$ then $ω_{n} > 0$ . Indeed, then the quadratic function $f (x) = 1 + x η θ + x^{2} τ η^{2}$ is non-negative.

(iv) If

τ > 0, η θ < 0

1 + η θ + τ η^{2} > 0

, and

θ^{2} \geq 4 τ

, then the sufficient condition is

\begin{matrix} 1 + min {q, 0} > \frac{| θ | + \sqrt{θ^{2} - 4 τ}}{2 τ | η |} . \end{matrix}

(4.28)

Indeed, inequality

1 + η θ + τ η^{2} > 0

implies

ω_{2} > 0

. Considering separately the cases

θ < 0

and

θ > 0

, the larger root of

1 + x η θ + x^{2} τ η^{2} = 0

\frac{| θ | + \sqrt{θ^{2} - 4 τ}}{2 τ | η |}

. If ( 4.28 ) holds then the left hand sides of inequalities ( 4.26 ), ( 4.27 ) imply

ω_{n} > 0

for

n \geq 3

It might be interesting to point out the explicit recurrence for the case

τ = 0

Corollary 4.6. Suppose

(X_{t})

is a quadratic harness with covariance ( 2.6 ) and parameters such that

τ = 0

- 1 < q < 1

, and

1 + [n]_{q} η θ + [n]_{q}^{2} σ θ^{2} > 0

for all

n

. Suppose that for

t > 0

random variable

X_{t}

has infinite support.

Then as the orthogonal martingale polynomials for

(X_{t})

we can take polynomials

p_{n} (x; t)

given by the recurrence

x p_{n} (x; t) = (1 + t σ [n]_{q}) p_{n + 1} (x; t)

+ (θ + t η + ([n]_{q} + [n - 1]_{q}) θ σ t) p_{n} (x; t) + t (1 + [n - 1]_{q} η θ + [n - 1]_{q}^{2} σ θ^{2}) [n]_{q} p_{n - 1} (x; t) .

Proof of Theorem 4.5 . We use the notation from ( 2.20 ). The integrability assumption of Theorem 4.1 is fulfilled by Theorem 2.5 and the system of recurrences in Theorem 4.1 simplifies. Since

α_{n} = [n - 1]_{q}

β_{n} = 1

, the first equation of recurrence ( 4.3 ) becomes

γ_{n + 1} = η + γ_{n} q

which gives

\begin{matrix} γ_{n} = η [n]_{q} . \end{matrix}

(4.29)

The other recurrences are solved inductively. Suppose that

\begin{matrix} δ_{n - 1} & = & (θ + η τ ([n - 1]_{q} + [n - 2]_{q})) [n - 1]_{q}, \end{matrix}

(4.30)

\begin{matrix} ω_{n} & = & (1 + [n - 1]_{q} η θ + [n - 1]_{q}^{2} τ η^{2}) [n]_{q}, \end{matrix}

(4.31)

hold. The initial conditions say that both formulas are satisfied for

n = 1

, and the second one holds true also for

n = 2

We use ( 4.29 ) and ( 4.30 ) to compute

δ_{n}

from the second equation in ( 4.3 ) as follows

δ_{n} = θ + q δ_{n - 1} + q ([n]_{q} - [n - 2]_{q}) τ γ_{n - 1} + η τ [n - 1]_{q} + η τ [n]_{q} = θ (1 + q [n - 1]_{q})

+ η τ (q ([n - 1]_{q} + [n - 2]_{q}) [n - 1]_{q} + q ([n]_{q} - [n - 2]_{q}) [n - 1]_{q} + [n - 1]_{q} + [n]_{q})

= θ [n]_{q} + η τ (q ([n - 1]_{q} + [n]_{q}) [n - 1]_{q} + [n - 1]_{q} + [n]_{q})

= θ [n]_{q} + η τ ([n - 1]_{q} + [n]_{q}) (1 + q [n - 1]_{q}) = θ [n]_{q} + η τ ([n - 1]_{q} + [n]_{q}) [n]_{q}

Finally, ( 4.4 ) determines

ω_{n + 1}

ω_{n + 1} = q ω_{n} + γ_{n}^{2} τ + q γ_{n} δ_{n} + γ_{n} θ + δ_{n} η + 1 - γ_{n} δ_{n} .

Using the inductive assumption ( 4.31 ), and already established formulas ( 4.29 ) and ( 4.30 ), we get

ω_{n + 1} = (1 + q [n]_{q}) + (1 + q [n - 1]_{q}) [n]_{q} η θ + q [n - 1]_{q}^{2} [n]_{q} τ η^{2} + τ η^{2} [n]_{q}^{2}

+ η (1 + q [n]_{q}) δ_{n} - η [n]_{q} δ_{n} = [n + 1]_{q} + [n]_{q}^{2} η θ + q [n - 1]_{q}^{2} [n]_{q} τ η^{2} + τ η^{2} [n]_{q}^{2} + η q^{n} δ_{n}

= [n + 1]_{q} + [n + 1]_{q} [n]_{q} η θ + τ η^{2} (q [n - 1]_{q} [n]_{q} ([n - 1]_{q} + q^{n - 1}) + [n]_{q}^{2} + q^{n} [n]_{q}^{2}) .

The coefficient at

τ η^{2}

simplifies to

q [n - 1]_{q} [n]_{q}^{2} + [n]_{q}^{2} + q^{n} [n]_{q}^{2} = [n + 1]_{q} [n]_{q}^{2} .

Thus

ω_{n + 1} = (1 + [n]_{q} η θ + [n]_{q}^{2} τ η^{2}) [n + 1]_{q}

Formula ( 4.25 ) comes from ( 2.19 ) after substituting ( 4.29 - 4.31 ) into ( 2.20 ).

□

4.4 Operator solutions

In this section we re-derive the recurrences for some special orthogonal martingale polynomials from Section 4.3 by an operator approach which is related to Lie algebra techniques. This method has a more ad hoc character so we concentrate on two relatively simple cases only.

In the operator approach, we go back directly to Theorem 2.3 . We re-interpret the matrices

x, y

as the linear operators acting on the formal power series in an auxiliary variable

z

. This identifies the martingale polynomial

p_{n} (x; t)

with

z^{n}

, where

z

is an auxiliary variable; similar technique appeared in umbral calculus [28,Ch.1] , in orthogonal polynomials [20] , and in Segal-Bargmann representation [26] .

We seek the solutions of ( 1.1 ) in terms of the

q

-differentiation operator

D q (g) (z) = {\begin{matrix} \frac{g (z) - g (q z)}{(1 - q) z} & q \neq 1, \\ g^{'} (z) & q = 1, \end{matrix}

and the multiplication operator

Z (g) (z) = z g (z),

treated as the linear operators on formal series

g (z)

in variable

z

. Table 1 lists the

q

-commutators of the combinations of these two operators that we need here and in Section 4.5 .

The requirement on the

0

-th column of

C_{t}

reduces to the requirement that on the unit constant function

x 1 = 0

y 1 = z

. When

σ = 0

, we have

a_{n} (t) = 1

so we are looking for the operators that satisfy

\begin{matrix} (t x + y) z^{n} = z^{n + 1} + lower order terms . \end{matrix}

(4.32)


$_{A} \^{B}$	$Z$	$Z D q$	$Z D^{2} q$	$Z^{2} D q$

$D q$	$I$	$D q$	$D^{2} q$	$(1 + q) Z D q - q (1 - q) Z^{2} D^{2} q$

$Z D q$	$Z$	$(1 - q) Z D q + q (1 - q) Z^{2} D^{2} q$		$Z^{2} D q$

$Z D^{2} q$	$(1 + q) Z D q - q (1 - q) Z^{2} D^{2} q$	$Z D^{2} q$

$Z^{2} D q$	$Z^{2}$

Table 1 .

q

-commutators

[A, B]_{q} : = A B - q B A

Example 4.7 (

q

-Meixner processes). The

q

-Meixner process is a quadratic harness with parameters

σ = η = 0

, see [12] . Inspecting Table 1 we verify that

x = D q

and

y = Z (1 + θ D q + τ D^{2} q)

solve ( 1.1 ) when

σ = η = 0

. Indeed,

[x, y]_{q} = [D q, Z]_{q} + θ [D q, Z D q]_{q} + τ [D q, Z D^{2} q]_{q} = I + θ D q + τ D^{2} q .

A calculation shows now that

(t x + y) z^{n} = z^{n + 1} + θ [n]_{q} z^{n} + (t + τ [n - 1]_{q}) [n]_{q} z^{n - 1},

which by the identification of

z^{n}

with

p_{n} (x; t)

implies that the corresponding martingale polynomials satisfy the three step recurrence

\begin{matrix} x p_{n} (x; t) = p_{n + 1} (x; t) + θ [n]_{q} p_{n} (x; t) + (t + τ [n - 1]_{q}) [n]_{q} p_{n - 1} (x; t), n \geq 1 . \end{matrix}

(4.33)

Of course, this is a special case of ( 4.25 ) corresponding to

η = 0

. Feinsilver [16,Section3.4] gives a reparametrization of recurrence ( 4.33 ), and considers

q

-commutator relations

[x, y]_{q} = h

as well as their realizations via operators

D q, Z

without considering ( 1.1 ). Anshelevich [2,Remark6] discusses a reparametrization of recurrence ( 4.33 ) in relation to free Sheffer systems. In [12] recurrence ( 4.33 ) is used as the first step in the construction of the

q

-Meixner Markov processes.

We now use the same method to derive the recurrence for the martingale polynomials of the bi-Poisson process.

Example 4.8 (bi-Poisson process). The bi-Poisson process is a quadratic harness with parameters

σ = τ = 0

. Inspecting Table 1 we verify that

x = D q + η Z (D q + θ D^{2} q)

and

y = Z (1 + θ D q)

solve ( 1.1 ) when

σ = τ = 0

. Indeed,

[x, y]_{q} = [D q, Z]_{q} + θ [D q, Z D q]_{q} + η [Z D q, Z]_{q} + η θ [Z D q, Z D q]_{q}

+ η θ [Z D^{2} q, Z]_{q} + η θ^{2} [Z D^{2} q, Z D q]_{q} = I + θ D q + η Z

+ η θ (1 - q) Z D q + η θ q (1 - q) Z^{2} D^{2} q + η θ (1 + q) Z D q - η θ q (1 - q) Z^{2} D^{2} q + η θ^{2} Z D^{2} q

= I + θ D q + η Z + η θ Z D q + η θ Z D q + η θ^{2} Z D^{2} q

= I + (θ D q + η θ Z D q + η θ^{2} Z D^{2} q) + (η Z + η θ Z D q)

= I + θ x + η y .

Operator

t x + y

satisfies the constraint

(t x + y) 1 = z

. A calculation shows that

(t x + y) z^{n} = z^{n + 1} + (θ + t η) [n]_{q} z^{n} + t (1 + η θ [n - 1]_{q}) [n]_{q} z^{n - 1},

so the constraint ( 4.32 ) holds. By the identification of

z^{n}

with

p_{n} (x; t)

, the corresponding martingale polynomials satisfy the three step recurrence

\begin{matrix} x p_{n} (x; t) = p_{n + 1} (x; t) + (θ + t η) [n]_{q} p_{n} (x; t) + t (1 + η θ [n - 1]_{q}) [n]_{q} p_{n - 1} (x; t), \end{matrix}

(4.34)

n \geq 0

, with

p_{- 1} = 0, p_{0} = 1

. Of course, this is a special case of ( 4.25 ) corresponding to

τ = 0

. But this recurrence was hard to guess without ( 1.1 ), so in [11] it appears for

q = 0

only.

4.5 Dual $q$ -commutation equation. Coherent states

Coherent states and the Segal-Bargmann representation are analytical techniques developed in mathematical physics [26] . The full Segal-Bargmann isomorphism is known to fail even in the relatively simple case of

q

-Brownian motion with

q < 0

, compare [31] . But algebraic duality is available and useful in a more general setting.

Let

\begin{matrix} Q_{t, s, u}^{*} (x, y) = A_{t, s, u} x^{2} + B_{t, s, u} y x + C_{t, s, u} y^{2} + D_{t, s, u} x + E_{t, s, u} y + F_{t, s, u} \end{matrix}

(4.35)

be the quadratic form in the non-commuting variables

x, y

; for background, see e.g. [21,page7] . Note that this is a dual of the quadratic form ( 3.32 ) that appears in the proof of Theorem 2.3 . The following relates ( 4.35 ) to the dual version of the

q

-commutation equation ( 1.1 ).

Proposition 4.9. Let

X_{t}

be a linear function

X_{t} = x + t y

of non-commutative variables

x, y

, and suppose that the coefficients of the quadratic form

Q_{t, s, u}^{*}

are given by ( 2.15 ) and ( 3.25 - 3.29 ). Then the following statements are equivalent.

(i) Operator identity

\begin{matrix} X_{t}^{2} = Q_{t, s, u}^{*} (X_{s}, X_{u}) \end{matrix}

(4.36)

holds for all

s < t < u

;

(ii) The non-commutative variables

x, y

satisfy the equation

\begin{matrix} [x, y]_{q} = σ x^{2} + τ y^{2} + η x + θ y + 1 . \end{matrix}

(4.37)

Proof. Since $X_{u} X_{s} = x^{2} + s u y^{2} + s (x y + y x) + (u - s) y x$ , we have
$X_{t}^{2} - Q_{t, s, u}^{*} (X_{s}, X_{u}) = (1 - A_{t, s, u} - B_{t, s, u} - C_{t, s, u}) x^{2} + (t^{2} - s^{2} A_{t, s, u} - s u B_{t, s, u} - u^{2} C_{t, s, u}) y^{2} + (t - s A_{t, s, u} - s B_{t, s, u} - u C_{t, s, u}) (x y + y x) - (u - s) B_{t, s, u} y x - (D_{t, s, u} + E_{t, s, u}) x - (s D_{t, s, u} + u E_{t, s, u}) y - F_{t, s, u} I .$
Relations ( 2.13 ), ( 3.23 ), ( 3.24 ) are equivalent to $X_{t}^{2} - Q_{t, s, u}^{*} (X_{s}, X_{u}) = - σ F_{t, s, u} x^{2} - τ F_{t, s, u} y^{2}$ $+ F_{t, s, u} (x y + y x) - (1 + q) F_{t, s, u} y x - η F_{t, s, u} x - θ F_{t, s, u} y - F_{t, s, u} I .$ $= F_{t, s, u} ([x, y]_{q} - σ x^{2} - τ y^{2} - η x - θ y - I) = 0 .$ □

4.5.1 Coherent states of $q$ -Meixner process

Continuing Example 4.7 , the analogue of the coherent state in physics is the generating function

\begin{matrix} φ_{t} (z, x) = \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} p_{n} (x; t), \end{matrix}

(4.38)

where

p_{n} (x; t)

are the orthogonal martingale polynomials given by recurrence ( 4.33 ).

Proposition 4.10. The operator

\begin{matrix} X_{t} = (1 + θ Z + τ Z^{2}) D q + t Z \end{matrix}

(4.39)

satisfies ( 4.36 ) with the quadratic form

Q_{t, s, u}^{*}

given by ( 2.15 ), ( 3.25 - 3.29 ) where

σ = η = 0

. Moreover,

\begin{matrix} x φ_{t} (z, x) = (X_{t} φ_{t}) (z, x) . \end{matrix}

(4.40)

Proof. Inspecting Table 1 it is easy to verify that $x = (1 + θ Z + τ Z^{2}) D q$ and $y = Z$ satisfy ( 4.37 ) with $σ = η = 0$ . Indeed, $[x, y]_{q} = [D q, Z]_{q} + θ [Z D q, Z]_{q} + τ [Z^{2} D q, Z]_{q}$ $= I + θ Z + τ Z^{2} = I + θ y + τ y^{2} .$ Therefore, Proposition 4.9 implies ( 4.36 ). The algebraic identity ( 4.40 ) follows from recurrence ( 4.33 ) by the following calculation. $x φ_{t} (z, x) = \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} x p_{n} (x; t)$ $= \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} (p_{n + 1} (x; t) + θ [n] p_{n} (x; t) + τ [n] [n - 1] p_{n - 1} (x; t) + t [n] p_{n - 1} (x; t)) =$ $= D q \sum_{n = 0}^{\infty} \frac{z^{n + 1}}{[n + 1]_{q}!} p_{n + 1} (x; t) + θ z D q \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} p_{n} (x; t)$ $+ τ z^{2} D q \sum_{n = 1}^{\infty} \frac{z^{n - 1}}{[n - 1]_{q}!} p_{n - 1} (x; t) + t z \sum_{n = 1}^{\infty} \frac{z^{n - 1}}{[n - 1]_{q}!} p_{n - 1} (x; t)$ $= (D q + θ Z D q + τ Z^{2} D q + t Z) φ_{t} (z, x) .$ □

Formula ( 4.40 ) is implicit in the usual derivation of the product formula for

φ_{t} (z, x)

, compare [1] . When the parameters

τ, θ

vanish, it appears in [31,Section3] in the context of analyzing ground states for the

q

-deformed Gauss distribution.

Ref. [3] gives a more general analytical scheme which coincides with ( 4.39 ) and ( 4.40 ) when

τ = 0, q = 1

; however, it advocates the normalization by the

L_{2}

-norm of the polynomials, which does not fit all the cases we are interested in.

In [12] we defined the

q

-Meixner process as a Markov process with the initial state

X_{0} = 0

and with the transition probabilities

P_{s, t} (x, d y)

determined as the unique probability measure orthogonalizing the polynomials

Q_{n}

in variable

y

which are given by the three step recurrence

\begin{matrix} y Q_{n} (y | x) \end{matrix}

(4.41)

\begin{matrix} = Q_{n + 1} (y | x, t, s) + (θ [n]_{q} + x q^{n}) Q_{n} (y | x, t, s) \end{matrix}

\begin{matrix} + (t - s q^{n - 1} + τ [n - 1]_{q}) [n]_{q} Q_{n - 1} (y | x, t, s) . \end{matrix}

In that paper, we showed that this Markov process is well defined, that it has orthogonal martingale polynomials

p_{n} (x; t)

given by recurrence ( 4.33 ), and we used this to prove that

(X_{t})

is a quadratic harness with parameters

σ = η = 0

Proposition 4.10 simplifies the verification of the quadratic harness condition in [12,Proposition3.4] , condensing more than three pages of proof into one page.

Proposition 4.11 ([12] ). If

| q | \leq 1

and the polynomials

p_{n} (x; t)

defined by ( 4.33 ) are orthogonal martingale polynomials for a Markov process

(X_{t})

, then

(X_{t})

is a quadratic harness with parameters

σ = η = 0

Proof. For simplicity, we consider only

| q | < 1

, as in this case ( 4.33 ) implies that

| X_{t} | \leq C t

has bounded support. Consider the generating function ( 4.38 ). Since

(X_{t})

is Markov, condition ( 2.3 ) is equivalent to

\begin{matrix} E (φ_{s} (z_{1}, X_{s}) X_{t} φ_{u} (z, X_{u})) \end{matrix}

(4.42)

\begin{matrix} = a t, s, u E (φ_{s} (z_{1}, X_{s}) X_{s} φ_{s} (z, X_{s})) + b t, s, u E (φ_{s} (z_{1}, X_{s}) X_{u} φ_{u} (z, X_{u})) \end{matrix}

holding for all

z_{1}, z

in a neighborhood of

0

(or just as an identity in formal power series in variables

z_{1}, z

We now use the martingale polynomial property which implies that

\begin{matrix} E [φ_{t} (z, X_{t}) | X_{s}] = φ_{s} (z, X_{s}), \end{matrix}

(4.43)

and we use ( 4.40 ) to represent the process through the operator ( 4.39 ) as

X_{t} φ_{t} (z, X_{t}) = X_{t} (φ_{t} (z, X_{t})) .

This gives

E (φ_{s} (z_{1}, X_{s}) X_{t} φ_{u} (z, X_{u})) = E (φ_{s} (z_{1}, X_{s}) X_{t} φ_{t} (z, X_{t}))

= X_{t} E (φ_{s} (z_{1}, X_{s}) φ_{s} (z, X_{s})) = X_{t} G_{s} (z_{1}, z),

where

G_{s} (z_{1}, z) = E (φ_{s} (z_{1}, X_{s}) φ_{s} (z, X_{s})) = \sum_{n = 0}^{\infty} \frac{(z_{1} z)^{n}}{[n]_{q}!^{2}} E (p_{n} (X_{s}; s)^{2})

= \sum_{n = 0}^{\infty} {\frac{(z_{1} z)^{n}}{[n]_{q}!}}^{n} \prod_{k = 1} (t + τ [k - 1]),

and

X_{t}

acts on

G_{s} (z_{1}, z)

as a series in variable

z

. Therefore, equation ( 4.42 ) is equivalent to

X_{t} G_{s} (z_{1}, z) = a t, s, u X_{s} G_{s} (z_{1}, z) + b t, s, u X_{u} G_{s} (z_{1}, z),

and follows from the (trivial) operator identity

X_{t} = a t, s, u X_{s} + b t, s, u X_{u} .

Similarly, condition ( 2.11 ) is equivalent to

\begin{matrix} E (φ_{s} (z_{1}, X_{s}) X_{t}^{2} φ_{u} (z_{2}, X_{u})) = A t, s, u E (φ_{s} (z_{1}, X_{s}) X_{s}^{2} φ_{s} (z_{2}, X_{s})) \end{matrix}

(4.44)

\begin{matrix} + B t, s, u E (φ_{s} (z_{1}, X_{s}) X_{s} X_{u} φ_{u} (z_{2}, X_{u})) + C t, s, u E (φ_{s} (z_{1}, X_{s}) X_{u}^{2} φ_{u} (z_{2}, X_{u})) \end{matrix}

\begin{matrix} + D_{t, s, u} E (φ_{s} (z_{1}, X_{s}) X_{s} φ_{s} (z_{2}, X_{s})) + E_{t, s, u} E (φ_{s} (z_{1}, X_{s}) X_{u} φ_{u} (z_{2}, X_{u})) \end{matrix}

\begin{matrix} + F_{t, s, u} E (φ_{s} (z_{1}, X_{s}) φ_{s} (z_{2}, X_{s})) . \end{matrix}

Notice that for

s \leq u

we have

E (φ_{s} (z_{1}, X_{s}) X_{s} X_{u} φ_{u} (z, X_{u})) = E (φ_{s} (z_{1}, X_{s}) X_{s} X_{u} φ_{u} (z, X_{u}))

= X_{u} E (φ_{s} (z_{1}, X_{s}) X_{s} φ_{s} (z, X_{s})) = X_{u} X_{s} E (φ_{s} (z_{1}, X_{s}) φ_{s} (z, X_{s})) = X_{u} X_{s} G_{s} (z_{1}, z) .

Therefore, equation ( 4.44 ) follows from the operator identity ( 4.36 ) which in expanded form says

X_{t}^{2} = A t, s, u X_{s}^{2} + B t, s, u X_{u} X_{s} + C t, s, u X_{u}^{2} + D t, s, u X_{s} + E t, s, u X_{u} + F t, s, u I,

applied to

G_{s} (z_{1}, z)

treated as a series in variable

z

. □

4.5.2 Coherent states of the bi-Poisson process

We now repeat the methods of Section 4.5.1 to derive new results about the bi-Poisson process from Example 4.8 .

Proposition 4.12. The operator

\begin{matrix} X_{t} = (1 + (θ + η t) Z + t η θ Z^{2}) D q + t Z, \end{matrix}

(4.45)

satisfies ( 4.36 ) for the quadratic form

Q_{t, s, u}^{*}

with the coefficients ( 2.15 ), ( 3.25 - 3.29 ) such that

σ = τ = 0

. Moreover, if

φ_{t} (z, x)

is the generating function ( 4.38 ) of the orthogonal martingale polynomials

p_{n} (x; t)

given by recurrence ( 4.34 ), then

\begin{matrix} x φ_{t} (z, x) = (X_{t} φ_{t}) (z, x) . \end{matrix}

(4.46)

Proof. Inspecting Table 1 , we verify that $x = (1 + θ Z) D q$ and $y = Z + η (Z + θ Z^{2}) D q$ solve ( 4.37 ) with $σ = τ = 0$ . Indeed, $[x, y]_{q} = [D q, Z]_{q} + θ [Z D q, Z]_{q} + η [D q, Z D q]_{q} + η θ [Z D q, Z D q]_{q} + η θ [D q, Z^{2} D q]_{q}$ $+ η θ^{2} [Z D q, Z^{2} D q]_{q} = I + θ Z + η D q + η θ ((1 - q) Z D q + q (1 - q) Z^{2} D^{2} q)$ $+ η θ ((1 + q) Z D q - q (1 - q) Z^{2} D^{2} q) + η θ^{2} Z^{2} D q = I + θ Z + η D q + η θ Z D q$ $+ η θ Z D q + η θ^{2} Z^{2} D q = I + η (I + θ Z) D q + θ (Z + η (Z θ Z^{2}) D q) = I + η x + θ y .$ Therefore, Proposition 4.9 implies ( 4.36 ).
We now derive ( 4.46 ) from ( 4.34 ) by the following calculation. $x φ_{t} (z, x) = \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} x p_{n} (x; t)$ $= \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} (p_{n + 1} (x; t) + (θ + t η) [n]_{q} p_{n} (x; t) + t (1 + η θ [n - 1]_{q}) [n]_{q} p_{n - 1} (x; t)) =$ $= D q \sum_{n = 0}^{\infty} \frac{z^{n + 1}}{[n + 1]_{q}!} p_{n + 1} (x; t) + (θ + t η) z D q \sum_{n = 0}^{\infty} \frac{z^{n}}{[n]_{q}!} p_{n} (x; t)$ $+ t z \sum_{n = 1}^{\infty} \frac{z^{n - 1}}{[n - 1]_{q}!} p_{n - 1} (x; t) + t η θ z^{2} D q \sum_{n = 1}^{\infty} \frac{z^{n - 1}}{[n - 1]_{q}!} p_{n - 1} (x; t)$ $= ((1 + (θ + η t) Z + t η θ Z^{2}) D q + t Z) φ_{t} (z, x) .$ □

Next we give the bi-Poisson version of Proposition 4.11 .

Proposition 4.13. Fix

- 1 \leq q \leq 1

. If polynomials

p_{n} (x; t)

given by ( 4.34 ) are orthogonal martingale polynomials for a Markov process

(X_{t})

, then

(X_{t})

is a bi-Poisson process, i.e., ( 2.3 ) and ( 2.14 ) hold with

σ = τ = 0

, and

1 + η θ \geq max {q, 0}

. Moreover, such Markov process

(X_{t})

is determined uniquely.

Proof. Property ( 2.6 ) follows from the explicit form of the polynomials $p_{0}, p_{1}, p_{2}$ .
Orthogonality gives $E (X_{t}) = E (p_{1} (X_{t}; t) p_{0} (X_{t}; t)) = 0$ and $E (X_{t}^{2}) = E ((p_{2} (X_{t}; t) + (θ + η t) X_{t} + t) p_{0} (X_{t})) = t$ . Martingale polynomial property then implies $E (X_{s} X_{t}) = E (X_{s} E [p_{1} (X_{t}; t) | ℱ_{\leq s}]) = E (X_{s}^{2}) = s$ .
If $- 1 \leq q < 1$ , from ( 4.34 ) it follows that $| X_{t} | \leq C t$ has bounded support; if $q = 1$ then it is not hard to identify the distribution of $| X_{t} |$ for example, from [30,Chapter4] or [13,pages175–181] , verifying that $| X_{t} |$ has a finite exponential moment.
Thus the moment problem has unique solution and the process $(X_{t})$ is determined uniquely. Moreover, polynomials are dense in $L_{2} (X_{s}, X_{u})$ , see [15,Theorem3.1.18] .
Consider the generating function ( 4.38 ) with polynomials $p_{n} (x; t)$ satisfying ( 4.34 ). Notice that property ( 4.43 ) follows again from the martingale polynomial condition. By Proposition 4.12 with $X_{t}$ defined by ( 4.45 ) we have the representation $X_{t} φ_{t} (z, X_{t}) = X_{t} (φ_{t} (z, X_{t})) .$ Since $(X_{t})$ is Markov and polynomials are dense in $L_{2} (X_{s}, X_{u})$ , condition ( 2.3 ) is again equivalent to ( 4.42 ), which we can interpret as the identity between the formal power series in variables $z_{1}, z$ . The latter follows again from the (trivial) operator identity $X_{t} = a t, s, u X_{s} + b t, s, u X_{u}$ , applied to $G_{s} (z_{1}, z) = E (φ_{s} (z_{1}, X_{s}) φ_{s} (z, X_{s})) = \sum_{n = 0}^{\infty} \frac{(z_{1} z)^{n}}{[n]_{q}!^{2}} E (p_{n} (X_{s}; s)^{2}),$ treated as the formal power series in variable $z$ .
Similarly, condition ( 2.11 ) is equivalent to ( 4.44 ), and again $E (φ_{s} (z_{1}, X_{s}) X_{s} X_{u} φ_{u} (z, X_{u})) = X_{u} X_{s} G_{s} (z_{1}, z) .$ Therefore, equation ( 4.44 ) follows from the operator identity ( 4.36 ) applied to $G_{s} (z_{1}, z)$ treated as the formal power series in variable $z$ .
Since the third coefficient in ( 4.34 ) is nonnegative for all $n$ , we get $1 + η θ \geq max {q, 0}$ , compare Remark 4.5 . □

Acknowledgement

We would like to thank M. Bożejko for numerous discussions and hospitality, L. Gallardo for preprint of [18] , M. Ismail, D. Pommeret, M. Yor and V. Zarikjan for helpful discussions. The second and the third authors are very grateful to their hosts for providing excellent research facilities and friendly atmosphere during their visit to the Department of Mathematics, University of Cincinnati, in August and September 2004.

References

W. A. Al-Salam and T. S. Chihara, Convolutions of orthonormal polynomials, SIAM J. Math. Anal., 7 (1976), pp. 16–28.
M. Anshelevich, Free martingale polynomials, Journal of Functional Analysis, (2003), pp. 228–261. axXiv:math.CO/0112194.
N. Asai, I. Kubo, and H.-H. Kuo, Segal-Bargmann transforms of one-mode interacting Fock spaces associated with Gaussian and Poisson measures, Proc. Amer. Math. Soc., 131 (2003), pp. 815–823.
R. Askey and M. Ismail, Recurrence relations, continued fractions, and orthogonal polynomials, Mem. Amer. Math. Soc., 49 (1984), pp. iv+108.
M. Bozejko and W. Bryc, On a class of free Lévy laws related to a regression problem. arxiv.org/abs/math.OA/0410601, 2004.
M. Bozejko, B. Kummerer, and R. Speicher, $q$ -Gaussian processes: non-commutative and classical aspects, Comm. Math. Phys., 185 (1997), pp. 129–154.
M. Bozejko and J. Wysoczanski, Remarks on $t$ -transformations of measures and convolutions, Ann. Inst. H. Poincaré Probab. Statist., 37 (2001), pp. 737–761.
W. Bryc, Some remarks on random vectors with nice enough behaviour of conditional moments, Bull. Polish Acad. Sci., 33 (1985), pp. 677–683.
W. Bryc, W. Matysiak, and J. Wesołowski, Bi-Poisson process II. In preparation, 2005.
W. Bryc and A. Plucinska, A characterization of infinite gaussian sequences by conditional moments, Sankhya A, 47 (1985), pp. 166–173.
W. Bryc and J. Wesołowski, Bi-Poisson process. Submitted. arxiv.org/abs/math.PR/0404241, 2004.
, Conditional moments of $q$ -Meixner processes, Probability Theory Related Fields, 131 (2005), pp. 415–441. arxiv.org/abs/math.PR/0403016.
T. S. Chihara, An introduction to orthogonal polynomials, Gordon and Breach, New York, 1978.
M. Dozzi, Two-parameter harnesses and the Wiener process, Z. Wahrsch. Verw. Gebiete, 56 (1981), pp. 507–514.
C. F. Dunkl and Y. Xu, Orthogonal polynomials of several variables, vol. 81 of Encyclopedia of Mathematics and its Applications, Cambridge University Press, Cambridge, 2001.
P. Feinsilver, Lie algebras and recurrence relations. III. $q$ -analogs and quantized algebras, Acta Appl. Math., 19 (1990), pp. 207–251.
U. Frisch and R. Bourret, Parastochastics., J. Math. Phys., 11 (1970), pp. 364–390.
L. Gallardo and M. Yor, Some new examples of Markov processes which enjoy the time-inversion property, Probability Theory Related Fields, (2005).
J. M. Hammersley, Harnesses, in Proc. Fifth Berkeley Sympos. Mathematical Statistics and Probability (Berkeley, Calif., 1965/66), Vol. III: Physical Sciences, Univ. California Press, Berkeley, Calif., 1967, pp. 89–117.
M. E. H. Ismail and D. Stanton, $q$ -integral and moment representations for $q$ -orthogonal polynomials, Canad. J. Math., 54 (2002), pp. 709–735.
C. Kassel, Quantum groups, vol. 155 of Graduate Texts in Mathematics, Springer-Verlag, New York, 1995.
A. Krystek and Ł. Wojakowski, Associative convolutions arising from conditionally free convolution. Preprint, December 2004.
I. G. Macdonald, Affine Hecke algebras and orthogonal polynomials, vol. 157 of Cambridge Tracts in Mathematics, Cambridge University Press, Cambridge, 2003.
R. Mansuy and M. Yor, Harnesses, Lévy bridges and Monsieur Jourdain, Stochastic Processes and Their Applications, 115 (2005), pp. 329–338.
M. Noumi and J. V. Stokman, Askey-Wilson polynomials: an affine Hecke algebra approach, in Laredo Lectures on Orthogonal Polynomials and Special Functions, Adv. Theory Spec. Funct. Orthogonal Polynomials, Nova Sci. Publ., Hauppauge, NY, 2004, pp. 111–144.
A. Perelomov, Generalized coherent states and their applications, Texts and Monographs in Physics, Springer-Verlag, Berlin, 1986.
A. Ronveaux and W. Van Assche, Upward extension of the Jacobi matrix for orthogonal polynomials, J. Approx. Theory, 86 (1996), pp. 335–357.
G.-C. Rota, Finite Operator Calculus, Academic Press, 1975.
G. Sansigre and G. Valent, A large family of semi-classical polynomials: the perturbed Chebyshev, J. Comput. Appl. Math., 57 (1995), pp. 271–281.
W. Schoutens, Stochastic processes and orthogonal polynomials, vol. 146 of Lecture Notes in Statistics, Springer-Verlag, New York, 2000.
H. van Leeuwen and H. Maassen, A $q$ -deformation of the Gauss distribution, J. Math. Phys., 36 (1995), pp. 4743–4756.
J. Wesołowski, Stochastic processes with linear conditional expectation and quadratic conditional variance, Probab. Math. Statist., 14 (1993), pp. 33–44.
D. Williams, Some basic theorems on harnesses, in Stochastic analysis (a tribute to the memory of Rollo Davidson), Wiley, London, 1973, pp. 349–363.
R. C. Zhang and X. W. Zhuang, Two-parameter harnesses and a characterization of Brownian sheets, Kexue Tongbao (Chinese), 33 (1988), pp. 1694–1697.
Z. G. Zhou, Two-parameter harnesses and the generalized Brownian sheet, Natur. Sci. J. Xiangtan Univ., 14 (1992), pp. 111–115.
X. W. Zhuang, The generalized Brownian sheet and two-parameter harnesses, Fujian Shifan Daxue Xuebao Ziran Kexue Ban, 4 (1988), pp. 1–9.

Department of Mathematics, University of Cincinnati, PO Box 210025, Cincinnati, OH 45221–0025, USA E-mail address : Wlodzimierz.Bryc@UC.edu Faculty of Mathematics and Information Science, Warsaw University of Technology, pl. Politechniki 1, 00-661 Warszawa, Poland E-mail address : matysiak@mini.pw.edu.pl Faculty of Mathematics and Information Science, Warsaw University of Technology, pl. Politechniki 1, 00-661 Warszawa, Poland E-mail address : wesolo@alpha.mini.pw.edu.pl

January 7, 2005; Revised: April 3, 2005

Włodzimierz Bryc, Wojciech Matysiak,

Jacek Wesołowski