A characterization of the infinitely divisible

## squared Gaussian processes

### November 27, 2006

Abstract
We show that, up to multiplication by constants, a Gaussian process has an infinitely divisible square if and only if its covariance is the Green function of a transient Markov process.
Short title : Infinitely divisible squared Gaussian. Keywords : Gaussian processes, infinite divisibility, Markov processes, local time.
AMS 2000 subject classification : 60E07, 60G15, 60J25, 60J55.
1 Introduction The question of the infinite divisibility of squared Gaussian vectors is an old problem which was first raised by Paul Lévy in 1948 [11. Given a centered Gaussian vector $\left({\phi }_{1},\cdots ,{\phi }_{p}\right)$  when can the vector $\left({\phi }_{1}^{2},\cdots ,{\phi }_{p}^{2}\right)$  be written as a sum of $n$  i.i.d. $p$  vectors for every $n\in \mathbb{N}$  ? Many authors have worked on this problem. We refer the readers to Paranjape [15, Griffiths [9, Evans [7, Moran and Vere-Jones [14 [18and the references therein for more on this problem.
In 1984, Griffiths [10has established a characterization of the $p$  -dimensional centered Gaussian vectors with an infinitely divisible square. His criterion is difficult to use since it requires the computation of the signs of all the cofactors of the covariance matrix. Indeed, except for the Brownian motion (and more generally Gaussian Markov processes), there were no examples, in the literature, of processes satisfying this remarkable property nor examples of processes lacking it.
We have recently shown [8that the family of fractional Brownian motions provide examples of both kinds. A fractional Brownian motion is a real-valued centered Gaussian process with a covariance given by $g\left(x,y\right)=|x{|}^{\beta }+|y{|}^{\beta }-|x-y{|}^{\beta }$  where the index $\beta$  is in $\left(0,2\right)$  . We proved that when $\beta$  is in $\left(0,1\right]$  then the square of this process is infinitely divisible, and when $\beta$  is in $\left(1,2\right)$  it is not. The critical index $1$  corresponds to the Brownian motion. In both cases we have used the criterion of Griffiths. We have shown in [8that if the Green function of a transient Markov process is symmetric, then it is the covariance of a Gaussian process with an infinitely divisible square. In particular when the index $\beta$  is in $\left(0,1\right]$  , the covariance of the corresponding fractional Brownian motion can be interpreted as the Green function of the symmetric stable process with index $\left(\beta +1\right)$  , killed at its hitting time of $0$  . To treat the second case, we have shown directly that the condition of Griffiths is not satisfied.
In view of these examples, a natural question arises. Is the representation of the covariance function of a centered Gaussian process as the Green function of a symmetric Markov process, a necessary condition for the infinite divisibility of its square ? We show in this paper that, up to a multiplication by a constant function, the answer is affirmative. The result is presented in Section 3 in the form of a necessary and sufficient condition . The proof is based this time on a criterion for infinite divisibility established by Bapat [1, which we recall in Section 2. Although, at first sight, the verification of this criterion is as difficult as that of the equivalent criterion of Griffiths (one has to check here too the sign of each cofactor of the covariance matrix) it allows to significantly shorten the arguments.
Gaussian processes with a covariance equal to a Green function, play an important role in the study of Markov processes. Indeed, the Isomorphism Theorem of Dynkin [6provides an indentity in law connecting each of these Gaussian processes to the local time process of the corresponding symmetric Markov process. This identity has been exploited to study properties of the local time process of the Markov process, using similar properties of the Gaussian processes. We choose to mention here only two references, but many more papers on the subject are quoted in [2for example. In [12Marcus and Rosen have studied sample path properties of the local time process using similar properties of the Gaussian process. In [2Bass, Eisenbaum, and Shi used the Isomorphism Theorem to establish the transience of the most visited sites of a symmetric stable process. The question of characterising the Gaussian processes with a covariance that is equal to a Green function of a symmetric Markov process has been open since the Isomorphism Theorem of Dynkin was first proved. The importance of an answer to this question is twofold; it will give a powerful tool for the study of these Gaussian processes as well as that of their associated Markov processes. The results of Section 3 provide an answer to this question.
In Section 5, we show that the Brownian sheet does not have an infinitely divisible square. Recalling that linear combinations of the squared components of centered Gaussian vectors are infinitely divisible, we find this example and that of the fractional Brownian motion with index in $\left(1,2\right)$  , somewhat counterintuitive.
2 The criterion of Bapat Let $A=\left({A}_{ij}{\right)}_{1\le i,j\le p}$  be a $p×p$  matrix. We write $A\ge 0$  if ${A}_{ij}\ge 0$  for all $i$  , $j$  .
Definition 2.1 A matrix $A=\left({A}_{ij}{\right)}_{1\le i,j\le p}$  is said to be an $M$  -matrix if (i) ${A}_{ij}\le 0$  for all $i\ne j$  (ii) $A$  is nonsingular and ${A}^{-1}\ge 0$  .
We refer the reader to [3for the theory of $M$  -matrices.
Definition 2.2 A diagonal matrix is called a signature matrix if all its entries are either $1$  or $-1$  .
The Laplace transform, $\psi \left({t}_{1},\cdots ,{t}_{p}\right)$  of the square of a $p-$  dimensional Gaussian vector is given by $\psi \left(t\right)=\left[det\left(I+GT\right){\right]}^{-1/2}$  where $t=\left({t}_{1}\cdots ,{t}_{p}\right)$  , $G$  is a positive definite $p×p$  -matrix, $T$  is the diagonal matrix with entries ${T}_{ii}={t}_{i}$  , and $I$  is the $p×p$  -identity matrix.
Bapat [1obtained the following characterization of the matrices $G$  for which the Laplace transform $\psi$  is infinitely divisible, that is, for which ${\psi }^{\delta }\left(t\right)$  is a Laplace transform for any $\delta >0$  .
Theorem A. The Laplace transform $\psi$  is infinitely divisible if, and only if, there exists a signature matrix $S$  such that $S{G}^{-1}S$  is an $M$  -matrix.
Remark 2.3 It is elementary to check, using Theorem A, that a centered $3$  -dimensional Gaussian vector $\left({\eta }_{1},{\eta }_{2},{\eta }_{3}\right)$  such that $IE\left({\eta }_{1}{\eta }_{2}\right)IE\left({\eta }_{2}{\eta }_{3}\right)IE\left({\eta }_{3}{\eta }_{1}\right)<0$  can not have an infinitely divisible square.
3 A necessary and sufficient condition for infinite divisibility In a previous work [8on squared Gaussian processes, we have established the following result. Theorem B. Let $g$  be the Green function of a strongly symmetric transient Borel right Markov process with state space $\left(E,\mathcal{ℰ}\right)$  . Let $\eta$  be a centered Gaussian process with covariance $g$  . Then the process ${\eta }^{2}$  is infinitely divisible.
The set $E$  (by the assumption of a right process) is a Borel subset of a compact metric space and $\mathcal{ℰ}$  is the $\sigma$  -field of its Borel sets. Note that if $\left({\eta }_{x}^{2}\left(x\right),x\in E\right)$  is an infinitely divisible squared Gaussian process, then for any $\mathcal{ℰ}$  -measurable real valued function $d$  , the process $\left({d}^{2}\left(x\right){\eta }_{x}^{2}\left(x\right),x\in E\right)$  is also infinitely divisible. We therefore have Corollary 3.1. Let $g$  be the Green function of a symmetric transient Markov process on a state space $\left(E,\mathcal{ℰ}\right)$  . Then for any $\mathcal{ℰ}$  -measurable real valued function $d$  on $E$  , there exists a centered Gaussian process $\eta$  with covariance equal to $\left(d\left(x\right)g\left(x,y\right)d\left(y\right),\left(x,y\right)\in E×E\right)$  . The process ${\eta }^{2}$  is infinitely divisible.
We shall first treat the case when $E$  is a finite set. Theorem 3.2 completes Corollary 3.1, for this case, by showing that its sufficient condition on the covariance, is also necessary for the infinite divisibility of ${\eta }^{2}$  .
Theorem 3.2 Let $\eta$  be a centered Gaussian vector, indexed by a finite set $E$  , with a positive definite covariance function $\left(G\left(x,y\right),\left(x,y\right)\in E×E\right)$  .
The vector ${\eta }^{2}$  is infinitely divisible if, and only if, there exists a real valued function $d$  on $E$  such that for any $x,y\in E$  , $G\left(x,y\right)=d\left(x\right)g\left(x,y\right)d\left(y\right),$  where the function $g$  is the Green function of a transient symmetric Markov process.
Remark 3.3 Note that the Green function $g$  of a symmetric Markov process $X$  is always positive definite. Indeed, it is semi-positive definite (see for example Marcus and Rosen [12or Eisenbaum [8) and it has been shown in [8section II, that for any ${x}_{1},{x}_{2},...{x}_{n}$  in the state space of $X$  , the matrix $\left(g\left({x}_{i},{x}_{j}\right){\right)}_{1\le i,j\le n}$  is invertible. Let $\eta$  be a centered Gaussian process indexed by an infinite set $E$  . Then it has an infinitely divisible square if for every finite subset $F$  of $E$  the covariance of $\left({\eta }_{x},x\in F\right)$  satisfies the condition of Theorem 3.2. This does not guarantee, a-priori, that for some deterministic function $d$  , the covariance of $\left(d\left(x\right){\eta }_{x},x\in E\right)$  is the Green function of a transient symmetric Markov process. Restricting our attention to $E=\mathbb{R}$  , and under the additional continuity assumption on the covariance function, the following theorem establishes the necessity of that condition.
Theorem 3.4 Let $\eta$  be a centered Gaussian process, indexed by $\mathbb{R}$  , with a positive definite covariance function $\left(G\left(x,y\right),\left(x,y\right)\in {\mathbb{R}}^{2}\right)$  . Assume that $G$  is jointly continuous. If the process ${\eta }^{2}$  is infinitely divisible , then there exists a measurable real valued positive function $d$  on $\mathbb{R}$  such that for any $x,y\in \mathbb{R}$  $G\left(x,y\right)=d\left(x\right)g\left(x,y\right)d\left(y\right),$  where the function $g$  is the Green function of a symmetric transient Markov process.
Theorem 3.4 is proved in Section 4. Although Theorem 3.2 was a good hint to anticipate Theorem 3.4, we could not use it directly to prove it.
Our argument is based on an explicit construction of the function $d$  and of the Green function $g$  of the Markov process that will be associated with the Gaussian process $\eta$  . As a first step of this construction we show that the covariance $G$  has to be positive. The proof of Theorem 3.4 remains valid, with simple changes, when the index set $\mathbb{R}$  is replaced by a separable locally compact metric space. Moreover, we will see in Section 4, that as a by-product of the proofs of Theorem 3.2 and Theorem 3.4, we obtain the following characterization of the associated Gaussian processes, which makes use of the definition below.
Definition 3.5 A $p×p$  matrix $A$  is said to be diagonally dominant if for every $i=1,...,p$  : ${\sum }_{j=1}^{p}{A}_{ij}\ge 0.$  Theorem 3.6 (i) A positive definite matrix $G$  is the Green function of a finite state space transient Markov process, if, and only if, ${G}^{-1}$  is a diagonally dominant $M$  -matrix.
(ii) Let $E$  be a separable locally compact metric space. Let $\left({\eta }_{x},x\in E\right)$  be a Gaussian process with a continuous positive definite covariance $\left(G\left(x,y\right),\left(x,y\right)\in {E}^{2}\right)$  . Then $G$  is the Green function of a transient Markov process on $E$  , if and only if, for every ${x}_{1},{x}_{2},...{x}_{p}$  in $E$  , the inverse of the matrix $\left(G\left({x}_{i},{x}_{j}\right){\right)}_{1\le i,j\le p}$  is a diagonally dominant $M$  -matrix.
4 Proofs Proof of Theorem 3.2 Let $E$  be the finite set $\left\{{x}_{1},\cdots ,{x}_{p}\right\}$  . Theorem A guarantees the existence of a signature matrix $S$  with diagonal ${S}_{i},i=1,\cdots ,p$  such that the covariance matrix of $\left({S}_{i}{\eta }_{{x}_{i}},i=1,\cdots p\right)$  is positive. We shall therefore assume, from the onset, that the covariance $G$  is positive.
We will actually prove that the covariance $G$  satisfies
 $\begin{array}{c}G\left(x,y\right)=\frac{D\left(y\right)}{D\left(x\right)}g\left(x,y\right),\end{array}$ (1)
where the function $g$  is the Green function of a transient Markov process and $D$  is a strictly positive function.
This will be sufficient to establish Theorem 3.2. Indeed, if $G$  satisfies (1) then we have $g\left(x,y\right)=\frac{D\left(x\right)}{D\left(y\right)}G\left(x,y\right).$  Let $U$  be the corresponding potential with density $g\left(x,y\right)$  with respect to a reference measure $\mu$  . That is, $Uf\left(x\right)=\int g\left(x,y\right)f\left(y\right)\mu \left(dy\right).$  We then set, $m\left(dy\right)=\mu \left(dy\right)/{D}^{2}\left(y\right)$  . With respect to $m$  , the kernel $U$  has densities $\stackrel{~}{g}\left(x,y\right)=D\left(x\right)G\left(x,y\right)D\left(y\right)$  . Consequently $\stackrel{~}{g}$  is a symmetric Green function. Setting $d=1/D$  , one obtains the necessity of the condition of Theorem 3.2.
To prove (1), let $G$  be the covariance matrix of a $p$  -dimensional centered Gaussian vector with an infinitely divisible square. By Theorem A, ${G}^{-1}$  is an $M$  -matrix. This implies (see Berman and Plemmons [3) that
 $\begin{array}{c}{G}^{-1}=cI-B,\end{array}$ (2)
where $B\ge 0$  and $c$  is strictly greater than the absolute value of any eigenvalue of $B$  . Further by [3( Chapter 6, p.137 M36), since ${G}^{-1}$  is an $M$  -matrix there exists a diagonal matrix $D$  such that ${D}_{ii}>0$  for all $i$  and $D{G}^{-1}{D}^{-1}$  has strictly positive row sum. This means that for any $i$  ${\sum }_{j=1}^{p}\left(D{G}^{-1}{D}^{-1}{\right)}_{ij}>0.$  According Definition 3.5, $D{G}^{-1}{D}^{-1}$  is diagonally dominant. Set $T=D\left(\frac{1}{c}B\right){D}^{-1}$  . Then $D{G}^{-1}{D}^{-1}=D\left(cI-B\right){D}^{-1}=c\left(I-T\right).$  Note that for any $i$  ${\sum }_{j=1}^{p}{T}_{ij}<1.$  The matrix $T$  is therefore the transition matrix of a transient Markov chain $\left({X}_{n}{\right)}_{n\in \mathbb{N}}$  with state space $E=\left\{{x}_{1},{x}_{2},...,{x}_{p}\right\}$  satisfying ${T}_{ij}=I{P}_{{x}_{i}}\left({X}_{1}={x}_{j}\right)$  , and ${\sum }_{j=1}^{p}I{P}_{{x}_{i}}\left({X}_{1}={x}_{j}\right)=1-I{P}_{{x}_{i}}\left({X}_{1}=\Delta \right)$  , where $\Delta$  denotes a cemetery point. Let ${\ell }_{\infty }^{{x}_{j}}$  denote the total number of visits of $X$  to the state ${x}_{j}$  . The Green function of $X$  is defined by $g\left({x}_{i},{x}_{j}\right)=I{E}_{{x}_{i}}\left({\ell }_{\infty }^{{x}_{j}}\right).$  It can be computed as follows $I{E}_{{x}_{i}}\left({\ell }_{\infty }^{{x}_{j}}\right)=I{E}_{{x}_{i}}\left({\sum }_{n=0}^{\infty }{1}_{\left\{{X}_{n}={x}_{j}\right\}}\right)={\sum }_{n=0}^{\infty }{T}_{ij}^{n}=\left(I-T{\right)}_{ij}^{-1},$  which is defined and is finite, since, by (2) and the discussion following it, the spectral radius of $T$  is stricly smaller than $1$  . Hence we can write $cDG{D}^{-1}=g;$  that is, for every $x$  , $y$  in $E$  $cG\left(x,y\right)=\frac{D\left(y\right)}{D\left(x\right)}g\left(x,y\right).$  We shall now use the well-known method to turn a finite state space Markov chain into a continuous time Markov process by subordination to a Poisson process with rate c (see for example Cinlar [5). By its construction, this Markov process is transient with potential density ( Green function) equal to $g$  and (1) is established. $\square$  Remark 4.1 Let $Y$  be the Markov process with potential density $g$  . Equation (1) looks as if $G$  is the potential density of an h-path transform of $Y$  . This is really the case if $D$  is excessive for $Y$  . Unfortunately this is not true in general. Consequently we see that the collection of covariance functions that correspond to Gaussian processes with an infinitely divisible square is somewhat richer than the set of Green functions of symmetric Markov processes. This remark remains true also when the index set $E$  of the Gaussian process is infinite.
To select covariance matrices that correspond to Green functions of Markov processes we have the condition given by Theorem 3.6 (i). Indeed, assume that ${G}^{-1}$  is a diagonally dominant M-matrix. Keeping the notations of the proof of Theorem 3.2, we can choose $D=I$  and obtain $T=\frac{1}{c}B$  . For any $i$  ${\sum }_{j=1}^{p}{T}_{ij}\le 1$  . Since the spectral radius of $T$  is strictly smaller then $1$  , for at least one $i$  the above inequality must be strict. Therefore $T$  is the transition matrix of a transient Markov chain and we can conclude that $G$  is the Green function of a transient Markov process with finite state space. To see that the condition is also necessary, consider the Green function $g$  of a transient symmetric Markov process $X$  with a finite state space $E=\left\{{x}_{1},{x}_{2},...,{x}_{p}\right\}$  . Then the inverse of the matrix $G=\left(g\left({x}_{i},{x}_{j}\right){\right)}_{1\le i,j\le p}$  is a diagonally dominant M-matrix. Indeed ${G}_{ij}={\lambda }_{j}\left(I-T{\right)}_{ij}^{-1}$  , where ${\lambda }_{j}$  is the expected value of the exponential sojourn times in state $j$  , $T$  is the transition matrix of the Markov chain $\left(X\left({S}_{n}\right){\right)}_{n\in \mathbb{N}}$  , and ${S}_{n}$  is the $n$  -th jump time of $X$  Proof of Theorem 3.4 We first show that $G$  has to be positive. For a fixed $n$  , we define the function ${d}_{n}$  on $\mathbb{R}$  by ${d}_{n}\left(x\right)=\frac{k}{{2}^{n}}if\frac{k}{{2}^{n}}\le x<\frac{k+1}{{2}^{n}}.$  Let $K$  be the compact set $\left[a,b\right]$  with $a  , and let ${d}_{n}\left(K\right)=\left\{{d}_{n}\left(x\right):x\in K\right\}$  .
The set ${d}_{n}\left(K\right)$  is finite set. Since the process $\left({\eta }_{{d}_{n}\left(x\right)}^{2},x\in K\right)$  is infinitely divisible, thanks to Theorem A there exists a signature function ${s}_{n}$  (ie a function taking values in $\left\{-1,+1\right\}$  ) on ${d}_{n}\left(K\right)$  such that for every $x,y\in K$  ${s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right)\ge 0.$  Note that
 $\begin{array}{c}{s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right)=|G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right)|\end{array}$ (3)
Since $G$  is continuous, we obtain by letting $n$  tend to $\infty$  that ${lim}_{n}{s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)\right)$  exists for all $\left(x,y\right)$  for which $G\left(x,y\right)\ne 0$  , and for such $\left(x,y\right)$  $\left({lim}_{n}{s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)\right)G\left(x,y\right)=|G\left(x,y\right)|.$  For $\left(x,y\right)$  in ${K}^{2}$  such that $G\left(x,y\right)=0$  , there exists a finite sequence ${a}_{1},{a}_{2},...,{a}_{p}$  such that $G\left(x,{a}_{1}\right)G\left({a}_{1},{a}_{2}\right)...G\left({a}_{p-1},{a}_{p}\right)G\left({a}_{p},y\right)\ne 0$  . Indeed, for $z\in K$  , set $C\left(z\right)=\left\{u\in \mathbb{R}:G\left(z,u\right)\ne 0\right\}$  . For each $z$  , $C\left(z\right)$  is an open set and ${\cup }_{z\in K}C\left(z\right)$  is a covering of the compact set $K$  . Thus there exists a finite subcovering $C\left({z}_{1}\right),C\left({z}_{2}\right),...C\left({z}_{m}\right)$  of the set $K$  . The sets of the covering are not disjoint. Let $C\left({z}_{{i}_{1}}\right)$  be one of these sets. $C\left({z}_{{i}_{1}}\right)$  is a union of disjoint intervals. If $C\left({z}_{{i}_{1}}\right)$  doesn't cover $K$  , then there exists ${z}_{{i}_{2}}$  in $\left\{{z}_{1},{z}_{2},...,{z}_{m}\right\}$  such that $C\left({z}_{{i}_{2}}\right)$  covers some of the endpoints of $C\left({z}_{{i}_{1}}\right)$  that are in $K$  . If $C\left({z}_{{i}_{1}}\right)\cup C\left({z}_{{i}_{2}}\right)$  do not cover $K$  there exists ${z}_{{i}_{3}}$  in $\left\{{z}_{1},{z}_{2},...,{z}_{m}\right\}$  such that $C\left({z}_{{i}_{3}}\right)$  covers some of the endpoints of $C\left({z}_{{i}_{1}}\right)\cup C\left({z}_{{i}_{2}}\right)$  that are in $K$  , and we may continue on until we exhaust all the finite covering above. Then we just have to make use of the sequence $\left({z}_{{i}_{1}},{z}_{{i}_{2}},...,{z}_{{i}_{m}}\right)$  to construct $\left({a}_{1},{a}_{2},...,{a}_{p}\right)$  ( $p\le m$  ) connecting $x$  to $y$  such that : $G\left(x,{a}_{1}\right)G\left({a}_{1},{a}_{2}\right)...G\left({a}_{p-1},{a}_{p}\right)G\left({a}_{p},y\right)\ne 0$  . Since ${s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}^{2}\left({d}_{n}\left({a}_{1}\right)\right){s}_{n}^{2}\left({d}_{n}\left({a}_{2}\right)\right)...{s}_{n}^{2}\left({d}_{n}\left({a}_{p}\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)={s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)$  we obtain, using ( 3 ), the existence of ${lim}_{n}{s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right)$  for all $x$  , $y$  . Set $H\left(x,y\right)={lim}_{n}{s}_{n}\left({d}_{n}\left(x\right)\right){s}_{n}\left({d}_{n}\left(y\right)\right).$  The function $H$  is symmetric and by ( 3 ) $H\left(x,y\right)=sign\left(G\left(x,y\right)\right)$  and by its definition for all $x,y,z\in K$  $H\left(x,y\right)=H\left(x,z\right)H\left(z,y\right).$  Hence there exists a signature function ${S}_{K}$  on $K$  (take $H\left(.,{z}_{0}\right)$  for a fixed ${z}_{0}$  in $K$  ) such that for any $x,y\in K$  ${S}_{K}\left(x\right){S}_{K}\left(y\right)G\left(x,y\right)=|G\left(x,y\right)|$  Denote by ${S}_{n}$  the function ${S}_{\left[-n,n\right]}$  , then repeating the above argument and letting $n$  tend to $\infty$  , we finally obtain the existence of $S$  satisfying $S\left(x\right)S\left(y\right)G\left(x,y\right)=|G\left(x,y\right)|.$  If $S$  is not identically equal to $1$  , then there exists a point $x$  at which $S\left(x\right)=1$  and $limin{f}_{y\to x}S\left(y\right)=-1$  . By continuity this means that $G\left(x,x\right)$  has to be equal to $0$  (because $G\left(x,{y}_{n}\right)<0$  for a sequence $\left({y}_{n}{\right)}_{n\ge 0}$  that converges to $x$  ).
This is excluded since $G$  is positive definite. Consequently $G$  is positive.
We define the measure $m$  on $\mathbb{R}$  by $m\left(dy\right)=\left(1\wedge \frac{1}{\sqrt{G\left(y,y\right)}}\right){e}^{-|y|}dy.$  Note that the measure $m$  is finite and $\int \sqrt{G\left(y,y\right)}m\left(dy\right)<\infty$  . Making use of the covariance inequality ( $G\left(x,y\right)\le \sqrt{G\left(x,x\right)G\left(y,y\right)}$  for $x$  , $y$  in $\mathbb{R}$  ), and the dominated convergence Theorem, we see that the function
 $\begin{array}{c}\chi \left(x\right)=\int G\left(x,y\right)m\left(dy\right),\end{array}$ (4)
is continuous.
We now consider the fixed compact set $K=\left[a,b\right]$  with $a  . For any integer $n$  , we keep the definition of ${d}_{n}$  introduced at the beginning of the proof. We set ${I}_{n}=\left\{k\in \mathbb{N}:-n{2}^{n}+1\le k\le n{2}^{n}-1and\frac{k}{{2}^{n}}\in {d}_{n}\left(K\right)\right\}.$  Define ${G}_{n}$  on the set $\left\{\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right):k,\ell \in {I}_{n}\right\}$  by ${G}_{n}\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right)=G\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right){\int }_{\frac{\ell }{{2}^{n}}}^{\frac{\ell +1}{{2}^{n}}}m\left(dy\right).$  Since the process ${\eta }^{2}$  is infinitely divisible and the positive measure $m$  has a support equal to $\mathbb{R}$  , ${G}_{n}^{-1}$  is an M-matrix ( ${G}_{n}$  need not be symmetric). Hence we can write ${G}_{n}^{-1}={c}_{n}I-{B}_{n},$  where ${B}_{n}\ge 0$  and ${c}_{n}$  is stricly greater than the absolute value of any eigenvalue of ${B}_{n}$  .
Define ${\chi }_{n}\left(\frac{k}{{2}^{n}}\right)={\sum }_{\ell \in {I}_{n}}{G}_{n}\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right)$  and let ${D}_{n}$  be the diagonal matrix $diag\left({\chi }_{n}\left(\frac{k}{{2}^{n}}\right),k\in {I}_{n}\right)$  , Then, ${D}_{n}{e}_{n}={G}_{n}{e}_{n}$  , where ${e}_{n}$  is equal to $\left(1,1,...1\right)$  . Consequently, for any $k\in {I}_{n}$  ${\sum }_{\ell \in {I}_{n}}\left[{D}_{n}^{-1}{G}_{n}^{-1}{D}_{n}{\right]}_{k,\ell }=\left({\chi }_{n}\left(\frac{k}{{2}^{n}}\right){\right)}^{-1}>0.$  Setting ${T}_{n}=\frac{1}{{c}_{n}}{D}_{n}^{-1}{B}_{n}{D}_{n}$  , we have
 $\begin{array}{c}{D}_{n}^{-1}{G}_{n}^{-1}{D}_{n}={c}_{n}\left(I-{T}_{n}\right).\end{array}$ (5)
and for every $k$  , ${\sum }_{\ell \in {I}_{n}}{T}_{n}\left(k,\ell \right)<1$  . Consequently, the matrix ${T}_{n}$  is the transition matrix of a transient Markov chain.
Rewriting (5) it follows that
 $\begin{array}{c}{G}_{n}={D}_{n}{O}_{n}\left({D}_{n}{\right)}^{-1},\end{array}$ (6)
with ${O}_{n}=\frac{1}{{c}_{n}}\left(I-{T}_{n}{\right)}^{-1}$  .
Let $A$  be a square matrix of size $|{I}_{n}|$  defined by $A=\left(A\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right){\right)}_{k,l\in {I}_{n}}$  .
We associate with $A$  an operator on the set of the continuous functions with compact support. We denote this operator $\mathbf{A}$  and define it as follows. Let $f$  be a continuous function on $\mathbb{R}$  with a compact support then the function $\mathbf{A}f$  is given by $\mathbf{A}f\left(x\right)={\sum }_{\ell \in {I}_{n}}A\left({d}_{n}\left(x\right),\frac{\ell }{{2}^{n}}\right)f\left(\frac{\ell }{{2}^{n}}\right)\forall x\in \mathbb{R}$  with the convention that $A\left(\frac{k}{{2}^{n}},\frac{\ell }{{2}^{n}}\right)=0$  when $k$  is outside of ${I}_{n}$  .
That way we associate with ${D}_{n}$  (resp. ${D}_{n}^{-1}$  , ${G}_{n}$  , ${O}_{n}$  ) the operator ${\mathbf{D}}_{\mathbf{n}}$  (resp.
${\mathbf{D}}_{\mathbf{n}}^{-1}$  , ${\mathbf{G}}_{\mathbf{n}}$  , ${\mathbf{O}}_{\mathbf{n}}$  ). Note that we have for every function $f$  and every $x$  in $\mathbb{R}$  ,
 $\begin{array}{ccc}{\mathbf{D}}_{\mathbf{n}}f\left(x\right)={\chi }_{n}\left({d}_{n}\left(x\right)\right)f\left({d}_{n}\left(x\right)\right)& & \end{array}$
 $\begin{array}{ccc}\left({\mathbf{D}}_{\mathbf{n}}{\right)}^{-1}f\left(x\right)=\left({\chi }_{n}\left({d}_{n}\left(x\right)\right){\right)}^{-1}f\left({d}_{n}\left(x\right)\right)& & \end{array}$
 $\begin{array}{ccc}{\mathbf{G}}_{\mathbf{n}}f\left(x\right)={\sum }_{\ell \in {I}_{n}}G\left({d}_{n}\left(x\right),\frac{\ell }{{2}^{n}}\right)f\left(\frac{\ell }{{2}^{n}}\right){\int }_{\frac{\ell }{{2}^{n}}}^{\frac{\ell +1}{{2}^{n}}}m\left(dy\right).& & \end{array}$
By (5), we know that the sequence $\left({\mathbf{O}}_{\mathbf{n}}{\right)}_{n\ge 0}$  satisfies
 $\begin{array}{ccc}{\mathbf{O}}_{\mathbf{n}}={\mathbf{D}}_{\mathbf{n}}^{-1}{\mathbf{G}}_{\mathbf{n}}{\mathbf{D}}_{\mathbf{n}}.& & \end{array}$
Moreover, for each $n$  , ${\mathbf{O}}_{\mathbf{n}}$  is the Green operator of a finite state space Markov process.
Let $O$  be the operator defined on ${\mathcal{C}}_{K}$  ( the continuous functions with support in $K$  ), by
 $\begin{array}{ccc}Of\left(x\right)={D}^{-1}\mathcal{G}Df\left(x\right)x\in K,& & \end{array}$
where
 $\begin{array}{ccc}Df\left(x\right)={\chi }_{{}_{K}}\left(x\right)f\left(x\right),& & \end{array}$
 $\begin{array}{ccc}{D}^{-1}f\left(x\right)=\left({\chi }_{{}_{K}}\left(x\right){\right)}^{-1}f\left(x\right)& & \end{array}$
and
 $\begin{array}{ccc}\mathcal{G}f\left(x\right)={\int }_{K}G\left(x,y\right)f\left(y\right)m\left(dy\right)& & \end{array}$
with
 $\begin{array}{ccc}{\chi }_{{}_{K}}\left(x\right)={\int }_{K}G\left(x,y\right)m\left(dy\right).& & \end{array}$
Note that by the continuity of $G$  and the fact that $G\left(x,x\right)>0$  and $m$  has $\mathbb{R}$  as its support, ${\chi }_{{}_{K}}$  is continuous, strictly positive on $K$  and thus bounded below by a strictly positive constant.
Lemma 4.2 For every function $f\in {\mathcal{C}}_{K}$  , the sequence $\left({\mathbf{O}}_{\mathbf{n}}f{\right)}_{n\ge 0}$  converges to $Of$  uniformly on $K$  Proof of Lemma 4.2 For any continuous function $f$  on ${\mathcal{C}}_{K}$  and any $x\in K$  , we have ${\mathbf{O}}_{\mathbf{n}}f\left(x\right)=\frac{1}{{\chi }_{n}\left({d}_{n}\left(x\right)\right)}{\int }_{K}G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right){\chi }_{n}\left({d}_{n}\left(y\right)\right)f\left({d}_{n}\left(y\right)\right)m\left(dy\right),$  and $Of\left(x\right)=\frac{1}{{\chi }_{{}_{K}}\left(x\right)}{\int }_{K}G\left(x,y\right){\chi }_{{}_{K}}\left(y\right)f\left(y\right)m\left(dy\right).$  Note that ${\chi }_{n}\left({d}_{n}\left(x\right)\right)={\int }_{K}G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right)m\left(dy\right).$  Since $G$  is uniformly continuous on compacts and $m$  is finite, the sequence $\left({\chi }_{n}\left({d}_{n}\left(x\right)\right){\right)}_{n\in \mathbb{N}}$  converges uniformly on $K$  to ${\chi }_{{}_{K}}\left(x\right)$  . Hence the sequence $G\left({d}_{n}\left(x\right),{d}_{n}\left(y\right)\right){\chi }_{n}\left({d}_{n}\left(y\right)\right)f\left({d}_{n}\left(y\right)\right)$  converges uniformly on $K×K$  to $G\left(x,y\right){\chi }_{{}_{K}}\left(y\right)f\left(y\right)$  . $\square$  We would like to show that $O$  is a Green operator. For this we shall use the following lemma which in essence is due to Hunt (see Meyer [13(chapter X, p.255) or Taylor [16). Lemma C Let $V$  be a kernel on a measurable space $\left(E,\mathcal{ℰ}\right)$  such that $V$  satisfies the complete maximum principal and the function $V1$  is bounded.
There exists then a subMarkovian resolvent $\left({V}_{p}\right)$  such that ${V}_{0}=V$  .
As Meyer has noted in [13(p.253), if we assume that $V$  is continuous then it is sufficient to verify the complete maximum principle for continuous functions with compact support only. More precisely, we have to verify that for any $a\ge 0$  , any two positive continuous functions with compact support, $\left(f,h\right)$  , if for all $x$  in $\left\{h>0\right\}$  $a+Vf\left(x\right)\ge Vh\left(x\right),$  then the inequality remains valid for all $x$  in $E$  .
Lemma 4.3 The kernel $O$  satisfies the complete maximum principle on $K$  .
Proof of Lemma 4.3 First note that $O$  maps continuous functions with support in $K$  to continuous functions and is therefore, a continuous kernel on $K$  . Let $f$  and $h$  be in ${\mathcal{C}}_{K}^{+}$  and $a\ge 0$  and suppose that
 $\begin{array}{c}a+Of\left(x\right)\ge Oh\left(x\right)forallx\in \left\{h>0\right\}.\end{array}$ (7)
Recall that $\left\{h>0\right\}$  is contained in $K$  . We need to show that ( 7 ) is satisfied by all $x\in K$  . Suppose that this is not true. Then there exists a constant $b>0$  and a point ${x}_{0}$  in $K$  so that $a+Of\left({x}_{0}\right)  Let $\epsilon >0$  and $N$  be such that for any $n>N$  $|{\mathbf{O}}_{\mathbf{n}}f\left(x\right)-Of\left(x\right)|<\epsilon \forall x\in K$  and $|{\mathbf{O}}_{\mathbf{n}}h\left(x\right)-Oh\left(x\right)|<\epsilon \forall x\in K.$  Such $N$  exists by Lemma 4.2 and the fact that $f,h\in {\mathcal{C}}_{K}^{+}$  .
By ( 7 ), for every $n>N$  and $x$  in $\left\{h>0\right\}$  ,
 $\begin{array}{c}a+2\epsilon +{\mathbf{O}}_{\mathbf{n}}f\left(x\right)\ge {\mathbf{O}}_{\mathbf{n}}h\left(x\right).\end{array}$ (8)
Recall now that for each $n$  , ${O}_{n}$  is a potential of a Markov process on ${I}_{n}$  . It follows from J.C. Taylor [17, that ${O}_{n}$  satisfies the complete maximum principle on ${I}_{n}$  and, by its definition, ${\mathbf{O}}_{\mathbf{n}}$  satisfies the complete maximum principle on $K$  . Consequently ( 8 ) is valid for all $x$  in $K$  . In particular, we have $a+2\epsilon +{\mathbf{O}}_{\mathbf{n}}f\left({x}_{0}\right)\ge {\mathbf{O}}_{\mathbf{n}}h\left({x}_{0}\right).$  On the other hand we have $a-2\epsilon +{\mathbf{O}}_{\mathbf{n}}f\left({x}_{0}\right)<{\mathbf{O}}_{\mathbf{n}}h\left({x}_{0}\right)-b.$  Choosing $\epsilon   leads to the desired contradiction. $\square$  In order to simplify the rest of the proof of Theorem 3.4 we use now, instead of $O$  , the operator ${V}_{K}$  defined by ${V}_{K}f\left(x\right)=\frac{1}{{\chi }_{{}_{K}}\left(x\right)}{\int }_{K}G\left(x,y\right)f\left(y\right)m\left(dy\right).$  First we note that ${V}_{K}1=1$  . Further since ${V}_{K}f=O\left(f/{\chi }_{{}_{K}}\right)$  , ${V}_{K}$  satisfies the complete maximum principle on $K$  , by Lemma 4.3.
Let ${K}_{n}$  be the compact set $\left[-n,n\right]$  for $n\in \mathbb{N}$  , and denote by ${V}_{{}_{{K}_{n}}}$  the corresponding operator as defined above for $K.$  Further, arguing as for ${\chi }_{{}_{K}}$  , one can show that $\chi$  defined in (4) is strictly positive.
Lemma 4.4 Let $V$  be the operator defined on bounded Borel functions by
 $\begin{array}{c}Vf\left(x\right)=\frac{1}{\chi \left(x\right)}{\int }_{\mathbb{R}}G\left(x,y\right)f\left(y\right)m\left(dy\right).\end{array}$ (9)
Then there exists a sub-Markovian resolvent $\left({V}_{p}\right)$  such that ${V}_{0}=V$  .
Proof of Lemma 4.4 Note that $V1=1$  , and that $V$  maps continuous functions with compact support to continuous functions on $\mathbb{R}$  . Let $f$  and $h$  be two positive continuous functions with compact supports and $a\ge 0$  and suppose that
 $\begin{array}{c}a+Vf\left(x\right)\ge Vh\left(x\right)forallx\in \left\{h>0\right\},\end{array}$ (10)
we need to show that this is satisfied for all $x$  in $\mathbb{R}$  .
We denote by $H$  a compact set that contains both the compact supports of $f$  and of $h$  . There exists ${n}_{0}$  such that for $n>{n}_{0}$  , H is contained in ${K}_{n}$  . Hence ( 10 ) can be written as $a\frac{\chi \left(x\right)}{{\chi }_{{}_{{K}_{n}}}\left(x\right)}+{V}_{{}_{{K}_{n}}}f\left(x\right)\ge {V}_{{}_{{K}_{n}}}h\left(x\right)forallx\in \left\{h>0\right\}.$  We have shown, when defining $m$  , that $\chi$  is a continuous function. Moreover, for any $n\in \mathbb{N}$  , the function ${\chi }_{{}_{{K}_{n}}}$  is continuous and the sequence $\left({\chi }_{{}_{{K}_{n}}}\right)$  is increasing and converges simply to $\chi$  . Consequently, by Dini's Theorem, $\left({\chi }_{{}_{{K}_{n}}}\right)$  converges uniformly to $\chi$  and $\chi \left(x\right)$  is strictly positive. Since the sequence $\left({\chi }_{{}_{{K}_{n}}}{\right)}_{n>{n}_{0}}$  is bounded below by a strictly positive constant on $H$  , $\chi /{\chi }_{{}_{{K}_{n}}}$  converges to $1$  uniformly on $H$  . Hence for every $\epsilon >0$  , there exists $N$  such that for every $n>N$  , $a\left(1+\epsilon \right)+{V}_{{}_{{K}_{n}}}f\left(x\right)\ge {V}_{{}_{{K}_{n}}}h\left(x\right)forallx\in \left\{h>0\right\}.$  Since the operator ${V}_{{}_{{K}_{n}}}$  satisfies the complete maximum principle on ${K}_{n}$  , the above inequality is still true for $x$  in ${K}_{n}$  . We now multiply each side of the inequality by ${\chi }_{{}_{{K}_{n}}}\left(x\right)/\chi \left(x\right)$  to obtain $a\left(1+\epsilon \right){\chi }_{{}_{{K}_{n}}}\left(x\right)/\chi \left(x\right)+Vf\left(x\right)\ge Vh\left(x\right)forallx\in {K}_{n}.$  Since ${\chi }_{{}_{{K}_{n}}}\left(x\right)/\chi \left(x\right)\le 1$  we finally get $a\left(1+\epsilon \right)+Vf\left(x\right)\ge Vh\left(x\right)forallx\in {K}_{n},$  and letting $n$  tend to $\infty$  $a\left(1+\epsilon \right)+Vf\left(x\right)\ge Vh\left(x\right)forallx\in \mathbb{R}.$  Since this is true for any $\epsilon >0$  , ( 10 ) is established for all $x$  in $\mathbb{R}$  .
The sub-Markovian resolvent $\left({V}_{p}\right)$  allows one to construct a semi-group $\left({P}_{t}\right)$  with $\left({V}_{p}\right)$  as its resolvent (see e.g. Bliedtner and Hansen [4Corollary 8.3 of chapter II) and with it a transient Markov process on $E=\mathbb{R}$  with $0$  -potential equal to $V$  . On the other hand, ( 9 ) implies that $m$  is a reference measure for the Markov process with potential V and that the potential density $h\left(x,y\right)$  of $V$  with respect to $m$  is equal to $h\left(x,y\right)=\frac{1}{\chi \left(x\right)}G\left(x,y\right).$  Setting then $\mu \left(dy\right)=\chi \left(y\right)m\left(dy\right)$  , we see that the operator $V$  admits the symmetric densities $\left(\frac{1}{\chi \left(x\right)}G\left(x,y\right)\frac{1}{\chi \left(y\right)},x,y\in \mathbb{R}\right)$  with respect to $\mu$  . $\square$  Remark 4.5 Assume that $G$  is such that for any $p$  , the inverse of the covariance matrix of $\left({\eta }_{{x}_{1}},{\eta }_{{x}_{2}},...,{\eta }_{{x}_{p}}\right)$  is a diagonally dominant M-matrix. One can then follow the steps of the proof of Theorem 3.4, without the need to define the matrices ${D}_{n}$  (we check similarly as in Remark 4.1, that the matrix ${T}_{n}=\frac{1}{{c}_{n}}{B}_{n}$  is a transition matrix) and conclude that there exists a symmetric Markov process with potential densites equal to $\left(G\left(x,y\right),\left(x,y\right)\in {\mathbb{R}}^{2}\right)$  . This together with Theorem 3.6 (i) lead to Theorem 3.6 (ii).
Finally, every centered Gaussian process with infinitely divisible square is equal to a deterministic function times a Gaussian process that is associated with a Markov process. The Isomorphism Theorem of Dynkin can hence be easily adapted to incorporate the deterministic functions, so that it can be used in studying any Gaussian processes with infinitely divisible square.
5 The case of the Brownian sheet A Brownian sheet is a centered Gaussian process $\left({W}_{x,s},x\ge 0,s\ge 0\right)$  with a covariance given by $IE\left({W}_{x,s}{W}_{y,t}\right)=\left(x\wedge y\right)\left(s\wedge t\right)$  Remember that up to a multiplicative constant $\left(x\wedge y;x,y\in {\mathbb{R}}_{+}\right)$  is the Green function of the linear Brownian motion killed at its first hitting time of $0$  . Hence this covariance is the product of two Green functions and one may ask whether $\left({W}_{x,s}^{2},x\ge 0,s\ge 0\right)$  is infinitely divisible. The answer is given in the following proposition.
Proposition 5.1 (i) For every $\left({x}_{i},{s}_{i}{\right)}_{1\le i\le 3}$  of ${\mathbb{R}}_{+}^{6}$  the vector $\left({W}_{{x}_{i},{s}_{i}}^{2},1\le i\le 3\right)$  is infinitely divisible.
(ii) The process $\left({W}_{x,s}^{2},x\ge 0,s\ge 0\right)$  is not infinitely divisible.
Proof. Using Griffiths criterion(or Bapat's criterion) it is easy to check that a sufficient condition for the infinite divisibility of the square of a three-dimensional Gaussian vector, indexed by $\left\{1,2,3\right\}$  and with a covariance $g$  , is
 $\begin{array}{c}g\left(i,j\right)g\left(k,k\right)\ge g\left(i,k\right)g\left(j,k\right),\end{array}$ (11)
for any choice of $i,j,k$  in $\left\{1,2,3\right\}$  .
Since any Green function satisfies ( 11 ), so does the covariance function of the Brownian sheet, as a product of two Green functions.
By Bapat's criterion, we know that a covariance matrix $G$  such that $G\ge 0$  , corresponds to an infinitely divisible square Gaussian vector if ${G}^{-1}$  is an M-matrix. We choose $\left({x}_{i},{s}_{i}{\right)}_{1\le i\le 4}$  such that $0<{x}_{1}<{x}_{3}<{x}_{2}<{x}_{4}$  and $0<{s}_{4}<{s}_{1}<{s}_{3}<{s}_{2}$  . Let $G$  be the matrix $\left(\left({x}_{i}\wedge {x}_{j}\right)\left({s}_{i}\wedge {s}_{j}\right){\right)}_{1\le i,j\le 4}$  . We compute the coefficient ${G}_{1,2}^{-1}$  ${G}_{1,2}^{-1}={x}_{1}{x}_{3}{s}_{4}\left({x}_{2}-{x}_{3}\right)\left({s}_{3}-{s}_{1}\right){s}_{4}>0.$  Hence ${G}^{-1}$  is not an M-matrix and the corresponding Gaussian vector does not have an infinitely divisible square. $\square$  Acknowledgements We would like to thank Gennady Samorodnitsky and Jay Rosen for helpful discussions and comments which helped us to improve a previous version of this work.
References

1. Bapat R. B. (1989). Infinite divisibility of multivariate gamma distribution and M-matrices.Sankhya, 51, 73-78.
2. Bass R., Eisenbaum N. and Shi Z. (2000). The most visited sites of symmetric stable processes.Proba. Theory and Relat. Fields 116 (3),391-404.
3. Berman A. and Plemmons R.J. (1979) Nonnegative Matrices in the Mathematical Sciences. Academic Press, Inc. New York..
4. Bliedtner J. and Hansen W. (1986).Potential Theory. Universitext Springer.
5. Cinlar E. (1975).Introduction to Stochastic Processes. Prentice-Hall, Inc. Englewood Cliffs, N.J.
6. Dynkin E. B. (1983). Local times and quantum fields. Seminar on Stochastic Processes 82, 69-84, Birkhauser Boston.
7. Evans S. N. (1991). Association and infinite divisibility for the Wishart distribution and its diagonal marginals.Jour. Multivar. Anal. 36 ,199-203.
8. Eisenbaum N. (2003). On the infinite divisibility of squared Gaussian processes. Probab. Theory and Relat. fields 125, 381-392.
9. Griffiths R. C. (1970). Infinite divisible multivariate gamma distributions. Sankhya, Series A,32 ,393-404.
10. Griffiths R. C. (1984). Characterization of infinitely divisible multivariate gamma distributions.Jour. Multivar. Anal. 15 ,12-20.
11. Lévy P. (1948).The arithmetical character of the Wishart distribution.Proc. Camb. Phil. Soc. 44, 295-297.
12. Marcus M.B. and Rosen J. (1992). Sample path properties of the local times of strongly symmetric Markov processes via Gaussian processes. (Special Invited paper).Ann. Probab. 20 ,1603-1684.
13. Meyer P. A. (1966). Probabilités et Potentiel. ed. Hermann, Paris.
14. Moran P.A.P. and Vere-Jones D. (1969). The infinite divisibility of multi-gamma distributions. Sankhya, Series A, 31, 191-194.
15. Paranjape S.R. (1978). Simpler proofs for the infinite divisibility of multivariate gamma distributions. Sankhya, Series A, 40 , 393-398.
16. Taylor J. C. (1972). On the existence of sub-Markovian resolvents. Inventiones Math. 17, 85-93.
17. Taylor J. C. (1975). A characterization of the kernel ${lim}_{\lambda \to 0}{V}_{\lambda }$  for sub-markovian resolvents $\left({V}_{\lambda }\right)$  .Annals of Probab. Vol.3,No.2 ,355-357.
18. Vere-Jones D. (1967). The infinite divisibility of a bivariate gamma distribution.Sankhlya, Series A,29,421-422.

 Nathalie Eisenbaum Haya Kaspi Laboratoire de Probabilités Industrial Engineering and Management Université Paris VI CNRS Technion 4 place Jussieu Haifa, Israel 32000 75252 Paris Cedex 05, France iehaya@tx.technion.ac.il nae@ccr.jussieu.fr