Processing math: 100%
Press "Enter" to skip to content

Month: November 2010

Intertwining and commutation relations for birth-death processes

I have posted today on arXiv a paper entitled Intertwining and commutation relations for birth-death processes, joint work with Aldéric Joulin.

Given a birth-death process on N with semigroup (Pt)t0 and a discrete gradient u depending on a positive weight u, we establish intertwining relations of the form uPt=Qtu, where (Qt)t0 is the Feynman-Kac semigroup with potential Vu of another birth-death process. We provide applications when Vu is positive and uniformly bounded from below, including Lipschitz contraction and Wasserstein curvature, various functional inequalities, and stochastic orderings. The proofs are remarkably simple and rely on interpolation, commutation, and convexity.

Let us give the main ingredient. We consider a birth-death process (Xt)t0 on the state space N:={0,1,2,}, i.e. a Markov process with transition probabilities given by

Pxt(y)=Px(Xt=y)=(λxt)1y=x+1+(νxt)1y=x1+(1(λx+νx)t)1y=x+tε(t).

The transition rates λ and ν are respectively called the birth and death rates of the process (Xt)t0. We assume that the process is irreducible, positive recurrent, and non-explosive. This holds when the rates satisfy to λ>0 on N and ν>0 on N and ν0=0 and

x=1λ0λ1λx1ν1ν2νx<andx=1(1λx+νxλxλx1++νxν1λxλ1λ0)=.

The unique stationary distribution μ of the process is reversible and is given by

μ(x)=μ(0)xy=1λy1νy, xNwithμ(0):=(1+x=1λ0λ1λx1ν1ν2νx)1.     (1)

Let us denote by F (respectively F+) the space of real-valued (respectively positive) functions f on N, and let bF be the subspace of bounded functions. The associated semigroup (Pt)t0 is defined for any function fbFF+ and xN as

Ptf(x)=Ex[f(Xt)]=y=0f(y)Pxt(y).

This family of operators is positivity preserving and contractive on Lp(μ), p[1,]. Moreover, the semigroup is also symmetric in L2(μ) since λxμ(x)=ν1+xμ(1+x) for any xN (detailed balance equation). The generator L of the process is given for any fF and xN by

Lf(x)=λx(f(x+1)f(x))+νx(f(x1)f(x))=λxf(x)+νxf(x),

where and are respectively the forward and backward discrete gradients on N:

f(x):=f(x+1)f(x)andf(x):=f(x1)f(x).

Our approach is inspired from the remarkable properties of two special birth-death processes: the M/M/1 and the M/M/ queues. The M/M/ queue has rates λx=λ and νx=νx for positive constants λ and ν. It is positive recurrent and its stationary distribution is the Poisson measure μρ with mean ρ=λ/μ. If Bx,p stands for the binomial distribution of size xN and parameter p[0,1], the M/M/ process satisfies for every xN and t0 to the Mehler type formula

L(Xt|X0=x)=Bx,eνtμρ(1eνt).     (2)

The M/M/1 queueing process has rates λx=λ and νx=ν1N{0} where 0<λ<ν are constants. It is a positive recurrent random walk on N reflected at 0. Its stationary distribution μ is the geometric measure with parameter ρ:=λ/ν given by μ(x)=(1ρ)ρx for all xN. A remarkable common property shared by the M/M/1 and M/M/ processes is the intertwining relation

L=LV     (3)

where LV=LV is the discrete Schrödinger operator with potential V given by

  • V(x):=ν in the case of the M/M/ queue
  • V(x):=ν1{0}(x) for the M/M/1 queue.

The operator LV is the generator of a Feynman-Kac semigroup (PVt)t0 given by

PVtf(x)=Ex[f(Xt)exp(t0V(Xs)ds)].

The intertwining relation (3) is the infinitesimal version at time t=0 of the semigroup intertwining

Ptf(x)=PVtf(x)=Ex[f(Xt)exp(t0V(Xs)ds)].     (4)

Conversely, one may deduce (4) from (3) by using a semigroup interpolation. Namely, if we consider

s[0,t]J(s):=PVsPtsf

with V as above, then (4) rewrites as J(0)=J(t) and (4) follows from (3) since

J(s)=PVs(LVPtsfLPtsf)=0.

Let us fix some uF+. The u-modification of the original process (Xt)t0 is a birth-death process (Xu,t)t0 with semigroup (Pu,t)t0 and generator Lu given by

Luf(x)=λuxf(x)+νuxf(x),

where the birth and death rates are respectively given by

λux:=ux+1uxλx+1andνux:=ux1uxνx.

One can check that the measure λu2μ is reversible for (Xu,t)t0. As consequence, the process (Xu,t)t0 is positive recurrent if and only if λu2 is μ-integrable. We define the discrete gradient u and the potential Vu by

u:=(1/u)andVu(x):=νx+1νux+λxλux.

Let φ:RR+ be a smooth convex function such that for some constant c>0, and for all rR,

φ(r)rcφ(r).     (5)

In particular, φ vanishes at 0, is non-increasing on (,0) and non-decreasing on (0,). Moreover, the behavior at infinity is at least polynomial of degree c. Note that one can easily find a sequence of such functions converging pointwise to the absolute value ||.

Theorem 1 (Intertwining and sub-commutation) Assume that for every xN and t0, we have

Ex[exp(t0Vu(Xu,s)ds)]<.

Then for every fbF, xN, t0,

uPtf(x)=PVuu,tuf(x)=Ex[uf(Xu,t)exp(t0Vu(Xu,s)ds)].     (6)

Moreover, if Vu0 then for every fbF, xN, t0,

φ(uPtf)(x)Ex[φ(uf)(Xu,t)exp(t0cVu(Xu,s)ds)].     (7)

Proof: Let us prove (7). If we define

s[0,t]J(s):=PcVuu,sφ(uPtsf)

then (7) rewrites as J(0)J(t). Hence it suffices to show that J is non-decreasing. We have the intertwining relation

uL=LVuuu,     (8)

where Lu is the generator of the u-modification process (Xu,t)t0 and where

LVuu:=LuVu.

Now

J(s)=PcVuu,s(T)whereT=LcVuuφ(uPtsf)φ(uPtsf)uLPtsf.

Letting gu=uPtsf, we obtain, by using (8),

T=LcVuuφ(gu)φ(gu)LVuugu

and thus

T=λu(φ(gu)φ(gu)gu)+νu(φ(gu)φ(gu)gu)+Vu(φ(gu)gucφ(gu)).

Now (5) and Vu0 give T0. Since the Feynman-Kac semigroup (PcVuu,t)t0 is positivity preserving, we get (7). The proof of (6) is similar but simpler (T is identically zero). ◻

The identity (6) implies a propagation of monotonicity: if f is non-increasing then Ptf is also non-increasing.

Actually, the intertwining relations above have their counterpart in continuous state space. Let A be the generator of a one-dimensional real-valued diffusion (Xt)t0 of the type

Af=σ2f+bf,

where f and the two functions σ,b are sufficiently smooth. Given a smooth positive function a on R, the gradient of interest is af=af. Denote (Pt)t0 the associated diffusion semigroup. Then it is not hard to adapt to the continuous case the argument of theorem~1 to show that the following intertwining relation holds:

aPtf(x)=Ex[af(Xa,t)exp(t0Va(Xa,s)ds)].

Here (Xa,t)t0 is a new diffusion process with generator

Aaf=σ2f+baf

and drift ba and potential Va given by

ba:=2σσ+b2σ2aaandVa:=σ2aab+aaba.

In particular, if the weight a=σ, where σ is assumed to be positive, then the two processes above have the same distribution and by Jensen's inequality, we obtain

|σPtf(x)|Ex[|σf(Xt)|exp(t0(σσb+bσσ)(Xs)ds)].

Hence under the assumption that there exists a constant ρ such that

infσσb+bσσρ,

then we get |σPtf|eρtPt|σf|. This type of sub-commutation relation is at the heart of the Bakry-Emery calculus for diffusions.

Leave a Comment
Syntax · Style · .