Factorized duality, stationary product measures and generating functions

TL;DR

This paper characterizes all factorized duality functions for a broad class of interacting particle and diffusion systems, revealing that only two main families of dualities encompass all cases, using generating functions and stationary measures.

Contribution

It provides a complete classification of factorized duality functions for interacting particle and diffusion systems, extending previous results and unifying duality concepts.

Findings

All factorized duality functions are classified into two main families.

The method unifies duality functions for particle and diffusion systems.

The approach reveals that only these two families cover all cases.

Abstract

We find all factorized duality functions for a class of interacting particle systems. The functions we recover are self-duality functions for interacting particle systems such as zero-range processes, symmetric inclusion and exclusion processes, as well as duality and self-duality functions for their continuous counterparts. The approach is based on, firstly, a general relation between factorized duality functions and stationary product measures and, secondly, an intertwining relation provided by generating functions. For the interacting particle systems, these self-duality and duality functions turn out to be generalizations of those previously obtained in [9] and, more recently, in [8]. Thus, we discover that only these two families of dualities cover all possible cases. Moreover, the same method discloses all self-duality functions for interacting diffusion systems such as the Brownian energy process, where both the process and its dual are in continuous variables.

Equations

$${\hat{\mathbb{E}}}_{\xi }D(\xi (t),\eta )={\mathbb{E}}_{\eta }D(\xi ,\eta (t)).$$

$$(\hat{S}(t){)}_{\text{left}}D=(S(t){)}_{\text{right}}D,$$

$$(\hat{L}{)}_{\text{left}}D=(L{)}_{\text{right}}D.$$

$$D(\xi ,\eta )=\prod _{x\in V}d({\xi }_x,{\eta }_x),\eta \in E^V,\xi \in {\hat{E}}^V.$$

$$d(0,\cdot )\equiv 1.$$

$$Lf(\eta )$$

$$L_{x,y}f(\eta )$$

$${\nu }_{\lambda }(n)$$

$$u(n)$$

$$v(n)$$

$$d(1,n)=a+bn,$$

$$u({\eta }_x)v({\eta }_y)(d(1,{\eta }_x-1)-d(1,{\eta }_x))+u({\eta }_y)v({\eta }_x)(d(1,{\eta }_x+1)-d(1,{\eta }_x))$$

$$=p(x,y)u(1)v(0)(d(1,{\eta }_y)-d(1,{\eta }_x)).$$

$$d(1,n+1)+d(1,n-1)-2d(1,n)=0,$$

$$u({\eta }_x)v({\eta }_y)-u({\eta }_y)v({\eta }_x)=-u(1)v(0)({\eta }_y-{\eta }_x),$$

$$d(0,n)$$

$$d(1,n)$$

$$\mathcal{L}f(z)$$

$${\mathcal{L}}_{x,y}f(z)$$

$$L^{N}f(\frac{1}{N}\eta )$$

$$L_{x,y}^{N}f(\frac{1}{N}\eta )={\eta }_x(\gamma -{\eta }_y)(f(\frac{1}{N}{\eta }^{x,y})-f(\frac{1}{N}\eta ))+{\eta }_y(\gamma -{\eta }_x)(f(\frac{1}{N}{\eta }^{y,x})-f(\frac{1}{N}\eta )),\eta \in {\mathbb{N}}^V.$$

$$\underset{N\to \infty }{\lim }(L^{N}f)(\frac{1}{N}{\eta }^N),$$

$${\mathcal{L}}_{x,y}f(z)$$

$${\nu }_{\lambda }(dz)$$

$$d(k,z)$$

$$d(1,z)$$

$$d(1,z_y)-d(1,z_x)$$

$$\int D(\xi ,\eta ){\mu }_{\lambda }(d\eta )=\prod _{x\in V}\int d({\xi }_x,{\eta }_x){\nu }_{\lambda }(d{\eta }_x)$$

$$\theta (\lambda )$$

$${\hat{\mathbb{E}}}_{\xi }H(\xi (t))=H(\xi ),$$

$${\int }_{E}d({\xi }_x,{\eta }_x){\nu }_{\ast }(d{\eta }_x)={\int }_{E}d({\xi }_x,{\eta }_x){\nu }_{\star }(d{\eta }_x)<\infty ,$$

Full text

For the interacting particle systems, these self-duality and duality functions turn out to be generalizations of those previously obtained in [9] and, more recently, in [8]. Thus, we discover that only these two families of dualities cover all possible cases. Moreover, the same method discloses all self-duality functions for interacting diffusion systems such as the Brownian energy process, where both the process and its dual are in continuous variables.

1 Introduction

Duality and self-duality are very useful and powerful tools that allow to analyze properties of a complicated system in terms of a simpler one. In case of self-duality for particle systems, the dual system is the same and the simplification arises because in the dual one considers only a finite number of particles (e.g. [3]).

Several methods are available to construct dual processes and duality relations. In the context of population dynamics, the starting point is always to consider as dual the backward coalescent process (for a general overview, see [2]). In the context of particle systems, the algebraic method developed in [9] offers a framework to construct duality functions via symmetries of the generator and reversible measures.

However, a complete picture of how to obtain all duality relations is missing. In this latter context of interacting particle systems, natural questions of the same sort are: which particle systems allow self-duality and is it possible to obtain all factorized duality functions for these systems?

One of the useful applications of disposing of all duality functions is that, depending on the target, one can choose appropriate ones: e.g. in the hydrodynamic limit and the study of the structure of the stationary measures, the “classical” duality functions are the appropriate ones (see e.g. [3]), whereas in the study of (stationary and non-stationary) fluctuation fields and associated Boltzmann-Gibbs principles ([12], chapter 11), as well as in the study of speed of relaxation to equilibrium in $L^{2}$ or in the study of perturbation theory around models with duality (cf. [3]), “orthogonal” duality functions turn out to be very useful.

In this paper, we develop an approach to answer the above questions and systematically determine all duality functions and relations for some interacting particle and diffusion systems. In this route, starting from examples, we first investigate a general connection between stationary product measures and factorized duality functions. This shows, in particular, that for infinite systems with factorized self-duality functions, the only stationary measures which are ergodic (w.r.t. either space-translation or time) are in fact product measures.

Then we use this connection between stationary product measures and duality functions to recover all possible factorized duality functions from the stationary product measures. More precisely, we show that, given the first duality function, i.e. the duality function with a single dual particle, all other duality functions are determined. This provides a simple machinery to obtain all duality functions in processes such as Symmetric Exclusion Process (SEP), Symmetric Inclusion Process (SIP) and Independent Random Walkers (IRW). In particular, we recover via this method all orthogonal polynomial duality functions obtained in [8].

Moreover, we prove that in the context of conservative particle systems where the rates for particle hopping depend only on the number of particles of the departure and arrival sites, the processes SEP, SIP and IRW are the only systems which have self-duality with factorized self-duality functions and that the first duality function is necessarily an affine function of the number of particles.

Next, in order to prove that the “possible” duality functions derived via the method described above are indeed duality functions, we develop a method based on generating functions. This method, via an intertwining relation, allows to go from discrete systems (particle hopping dynamics) to continuous systems (such as diffusion processes or deterministic dynamics) and back, and also allows to pass from self-duality to duality and back. The proof of a self-duality relation in a discrete system then reduces to the same property in a continuous system, which is much easier to check directly.

The generating function method also provides new examples of self-duality for processes in the continuum such as the Brownian Energy Process (BEP), which intertwines with the SIP via the generating function. In fact, we show equivalence between self-duality of SIP, duality between SIP and BEP and self-duality of BEP. Finally, this method based on generating functions generalizes the concept of obtaining dualities from symmetries to intertwinings, being a symmetry an intertwining of the generator with itself.

The paper is organized as follows. In Section 2 we introduce the basic definitions of duality and systems considered. Additionally, in Theorem 2.1 we prove which particle systems out of those considered admit factorized self-duality.

In Section 3, we investigate a general relation between factorized duality functions and stationary product measures. We treat separately the finite and infinite contexts in which this relation arises; in the latter case, we exploit this connection to draw some conclusions on the product structure of ergodic measures.

Section 4 is devoted to the derivation of all possible factorized self-duality and duality functions. Here Theorem 2.1 and the relation in the previous section are the two key ingredients.

In Section 5, after an introductory example and a brief introduction on the general connection between duality and intertwining relations, we establish an intertwining between the discrete and the continuum processes. This intertwining relation is then used to produce all the self-duality functions for the Brownian Energy Process.

2 Setting

We start defining what we mean by duality of stochastic processes. Then, we introduce a general class of Markov interacting particle systems with associated interacting diffusion systems.

2.1 Duality

Given two state spaces $\Omega$ and $\hat{\Omega}$ and two stochastic processes $\{\xi(t),\ t\geq 0\}$ and $\{\eta(t),\ t\geq 0\}$ evolving on them, we say that they are dual with duality function $D:\hat{\Omega}\times\Omega\to\mathbb{R}$ (where $D$ is a measurable function) if, for all $t>0$, $\xi\in\hat{\Omega}$ and $\eta\in\Omega$, we have the so-called duality relation

[TABLE]

If the laws of the two processes coincide, we speak about self-duality.

More generally, we say that two semigroups $\{S(t),\ t\geq 0\}$ and $\{\hat{S}(t),\ t\geq 0\}$ are dual with duality function $D$ if, for all $t\geq 0$,

[TABLE]

where “left” (resp. “right”) refers to action on the left (resp. right) variable. Even more generally, we say that two operators $L$ and $\hat{L}$ are dual to each other with duality function $D$ if

[TABLE]

In order not to overload notation, we use the expression $A_{\text{\tiny left}}D(\xi,\eta)$ for $(AD(\cdot,\eta))(\xi)$ and, similarly, $B_{\text{\tiny right}}D(\xi,\eta)=(BD(\xi,\cdot))(\eta)$. We will often write $D(\xi,\eta)$ in place of $(\xi,\eta)\mapsto D(\xi,\eta)$.

2.2 The lattice and the factorization over sites

The underlying geometry of all systems that we will look at consists of a set of sites $V$ either finite or $V=\mathbb{Z}^{d}$. Moreover we are given a family of transition rates $p:V\times V\to\mathbb{R}_{+}$, satisfying the following conditions: for all $x$, $y\in V$,

1. (1)

Vanishing diagonal: $p(x,x)=0\ ,$ 2. (2)

Irreducibility: there exist $x_{1}=x$, $x_{2}$, $\ldots$, $x_{m}=y$ such that $\prod_{l=1}^{m-1}p(x_{l},x_{l+1})>0\ .$

In case of infinite $V$, we further require the following:

1. (3)

Finite-range: there exists $R>0$ such that, for all $x,y\in V$, $p(x,y)=0$ if $|x-y|>R\ ,$ 2. (4)

Uniform bound on total jump rate: $\sup_{x\in V}\sum_{y\in V}p(x,y)<\infty\ .$

Notice that when $p$ is finite-range and translation invariant, then the uniform bound on total jump rate follows automatically.

To each site $x\in V$ we associate a variable $\eta_{x}\in E=\mathbb{N}$, $\{0,\ldots,N\}$ or $\mathbb{R}_{+}$, with the interpretation of either the number of particles or the amount of energy associated to the site $x$. Configurations are denoted by $\eta\in\Omega=E^{V}$.

In all the examples that we will be discussing here, the duality functions factorize over sites, i.e.

[TABLE]

We then call the functions $d(\xi_{x},\eta_{x})$ the singe-site duality functions and further assume

[TABLE]

The above condition (4) is related to the fact that we want to have duality functions which make sense for infinite systems when the dual configuration has a finite total mass. A typical example is when $\eta$, $\xi\in\mathbb{N}^{V}$, where $\eta$ is an infinite configuration while $\xi$ is a finite configuration, so that in the product (3) there are only a finite number of factors different from $d(0,\eta_{x})$. In this sense, the choice $d(0,\cdot)\equiv 1$ is the only sensible one for infinite systems.

When $V$ is finite and $E=\mathbb{N}$ or $\{0,\ldots,N\}$, this condition is not necessary and e.g.  if a reversible product measure $\mu=\otimes_{x\in V}\ \nu$ exists, then the so-called cheap self-duality function $D(\xi,\eta)=\tfrac{1}{\mu(\xi)}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}{\{\xi=\eta\}}=\prod_{x\in V}\frac{1}{\nu(\xi_{x})}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{\xi_{x}=\eta_{x}\}$ does not satisfy (4).

2.3 Interacting particle systems with factorized self-duality

The class of interacting particle systems we consider is described by the following infinitesimal generator acting on local functions $f:\Omega\to\mathbb{R}$ as follows:

[TABLE]

where $L_{x,y}$, the single-edge generator, is defined as

[TABLE]

and $\eta^{x,y}$ denotes the configuration arising from $\eta$ by removing one particle at $x$ and putting it at $y$, i.e.  $\eta^{x,y}_{x}=\eta_{x}-1$, $\eta^{x,y}_{y}=\eta_{y}+1$, while $\eta^{x,y}_{z}=\eta_{z}$ if $z\neq x,y$. Note the conservative nature of the system and the form of the particle jump rates in (6) which depend on the number of particles in the departure and arrival site in a factorized form. Minimal requirements on the functions $u$ and $v$, namely

• (i)

$u(0)=0$, $u(1)=1$ and $u(n)>0$ for all $n>0\ ,$

• (ii)

$v(0)\neq 0$ and $v(N)=0$ if $E=\{0,\ldots,N\}$ and in all other cases v(n)>0\,

guarantee the existence of a one-parameter family of stationary (actually reversible) product measures $\{\otimes\ \nu_{\lambda},\ \lambda>0\}$ with marginals $\nu_{\lambda}$ given by

[TABLE]

for all $\lambda>0$ for which the normalizing constant $Z_{\lambda}<\infty$ and with $\varphi(n)=n!\prod_{m=1}^{n}\frac{v(m-1)}{u(m)}$.

2.4 Basic examples

We recall here the basic examples of self-dual interacting particle systems and corresponding factorized self-duality functions known in literature (cf. e.g. [9]).

• (I)

Independent random walkers (IRW)

• –

$E=\mathbb{N}\ ,$

• –

$u(n)=n\ ,$ $v(n)=1\ ,$

• –

$\nu_{\lambda}\sim\text{Poisson}(\lambda)\ ,$ $\nu_{\lambda}(n)=\tfrac{\lambda^{n}}{n!}e^{-\lambda}\ ,\ \lambda>0\ ,$

• –

$d(k,n)=\frac{n!}{(n-k)!}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{k\leq n\}\ .$

• (II)

Symmetric inclusion process ($\text{SIP}(\alpha)$, $\alpha>0$)

• –

$E=\mathbb{N}\ ,$

• –

$u(n)=n\ ,$ $v(n)=\alpha+n\ ,$

• –

$\nu_{\lambda}\sim\text{Gamma}_{\text{d}}(\alpha,\lambda)\ ,$ $\nu_{\lambda}(n)=\tfrac{\Gamma(\alpha+n)}{\Gamma(\alpha)}\tfrac{\lambda^{n}}{n!}(1-\lambda)^{\alpha}\ ,\ \lambda\in(0,1)\ ,$

• –

$d(k,n)=\tfrac{\Gamma(\alpha)}{\Gamma(\alpha+k)}\tfrac{n!}{(n-k)!}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{k\leq n\}\ .$

• (III)

Symmetric exclusion process ($\text{SEP}(\gamma)$, $\gamma\in\mathbb{N}$)

• –

$E=\{0,\ldots,N\}\ ,$

• –

$u(n)=n\ ,$ $v(n)=\gamma-n\ ,$

• –

$\nu_{\lambda}\sim\text{Binomial}(\gamma,\tfrac{\lambda}{1+\lambda})$, $\nu_{\lambda}(n)=\tfrac{\gamma!}{(\gamma-n)!}\tfrac{\lambda^{n}}{n!}\left(\tfrac{1}{1+\lambda}\right)^{\gamma}\ ,\ \lambda>0\ ,$

• –

$d(k,n)=\tfrac{(\gamma-k)!}{\gamma!}\tfrac{n!}{(n-k)!}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{k\leq n\}\ .$

In the following theorem we show that the only processes with generator of the type (5) which have non-trivial factorized self-duality functions are of one of the types described in the examples above, i.e. IRW, SIP or SEP. Here by “non-trivial” we mean that the first single-site self-duality function $d(1,n)$ is not a constant (as a function of $n$).

THEOREM** 2.1****.**

Assume that the process with generator (6) is self-dual with factorized self-duality function $D(\xi,\eta)=\prod d(\xi_{x},\eta_{x})$ in the form (3) with $d(0,\cdot)\equiv 1$ as in (4). If $d(1,n)$ is not constant as a function of $n$, then

[TABLE]

and the first single-site self-duality function is of the form

[TABLE]

for some $a\in\mathbb{R}$ and $b\neq 0$.

PROOF.  Using the self-duality relation for $\xi_{x}=1$ and no particles elsewhere, together with $u(0)=0$, we obtain the identity

[TABLE]

Setting $\eta_{x}=\eta_{y}=n\geq 1$, this yields anytime $u(n)v(n)\neq 0$

[TABLE]

from which we derive $d(1,n)=a+bn$. Because $d(1,n)$ is not constant as a function of $n$, we must have $b\not=0$. Inserting $d(1,n)=a+bn$ in (10) we obtain

[TABLE]

from which, by setting $\eta_{x}=n$ and $\eta_{y}=0$ we obtain the first in (2.1), while via $\eta_{x}=n$ and $\eta_{y}=1$ we get the second condition. $\;\qed$

REMARK** 2.1****.**

More generally, if we replace (4) with $d(0,n)\neq 0$ in the above statement, we analogously obtain (2.1) and

[TABLE]

for some constants $a,b,c\in\mathbb{R}$, $b,c\neq 0$.

2.5 Interacting diffusion systems as scaling limits

Interacting diffusion systems arise as scaling limits of the particle systems in Section 2.4 above (cf. [9]). More in details, by “scaling limit” we refer to the limit process of the particle systems $\{\tfrac{1}{N}\eta^{N}(t),\ t\geq 0\}_{N\in\mathbb{N}}$, where the initial conditions $\tfrac{1}{N}\eta^{N}(0)=\tfrac{1}{N}(\lfloor z_{x}N\rfloor)_{x\in V}$ converge to some $z\in E^{V}$, with $E=\mathbb{R}_{+}$.

In case of IRW, one obtains a deterministic (hence degenerate diffusion) process $\{z(t),\ t\geq 0\}$ whose evolution is described by a first-order differential operator. In case of $\text{SIP}(\alpha)$, the scaling limit is a proper Markov process of interacting diffusions known as Brownian Energy Process ($\text{BEP}(\alpha)$) (cf. [9]). For the $\text{SEP}(\gamma)$, this limit cannot be taken in the sense of Markov processes, but we can extend the SEP generator to a larger class of functions defined on a larger configuration space and take the many-particle limit. The limiting second-order differential operator is then not a Markov generator, but still a second order differential operator. We will explain this more in detail below.

The limiting differential operators in the case of IRW and $\text{SIP}(\alpha)$ can be described as acting on smooth functions $f:E^{V}\to\mathbb{R}$ as follows:

[TABLE]

with single-edge generators ${\mathscr{L}}_{x,y}$ given, respectively, by

[TABLE]

and

[TABLE]

For the $\text{SEP}(\gamma)$ we proceed as follows. For each $N\in\mathbb{N}$, consider the operator $L^{N}$ working on functions $f:(\mathbb{N}/N)^{V}\to\mathbb{R}$ as

[TABLE]

where

[TABLE]

This operator is not a Markov generator anymore, because the factors $\eta_{x}(\gamma-\eta_{y})$ can become negative. With this operator, we consider the limit

[TABLE]

where $\eta^{N}=(\lfloor Nz_{x}\rfloor)_{x\in V}$ and $f:E^{V}\to\mathbb{R}$ is a smooth function. This then gives the differential operator ${\mathscr{L}}$ which is the analogue of (13) in the context of $\text{SEP}(\gamma)$. This differential operator ${\mathscr{L}}$, with single-edge operators

[TABLE]

does not generate a Markov process but it is still useful because, as we will see in Section 5 below, via generating functions, it is intertwined with the operator (14) for the choice $N=1$.

Naturally, as we can see for the case of $\text{SIP}(\alpha)$ and $\text{BEP}(\alpha)$, when going to the scaling limit, some properties concerning stationary measures and duality pass to the limit. Indeed, $\text{BEP}(\alpha)$ admits a one-parameter family of stationary product measures $\{\otimes\ \nu_{\lambda},\ \lambda>0\}$, where $\nu_{\lambda}\sim\text{Gamma}(\alpha,\lambda)$, namely

[TABLE]

and is dual to $\text{SIP}(\alpha)$ with factorized duality function $D(\xi,z)=\prod_{x\in V}d(\xi_{x},z_{x})$ given by

[TABLE]

After noting that property (4) holds also in this situation, we show that the first single-site duality functions $d(1,x)$ between $\text{SIP}(\alpha)$ and $\text{BEP}(\alpha)$ are affine functions of $z\in\mathbb{R}_{+}$, as we found earlier for single-site self-dualities in Theorem 2.1.

PROPOSITION** 2.1****.**

Assume that $\text{SIP}(\alpha)$ and $\text{BEP}(\alpha)$’s single-edge generators are dual with factorized duality function $D(\xi,z)=\prod d(\xi_{x},z_{x})$ with $d(0,\cdot)\equiv 1$ as in (3)–(4). Then

[TABLE]

for some $a,b\in\mathbb{R}$.

PROOF.  The duality relation for $\xi_{x}=1$ and no particles elsewhere, by using (4), reads

[TABLE]

If we set $z_{x}=z_{y}=z$, then $z^{2}\frac{d^{2}}{dz^{2}}d(1,z)=0$ leads to (17) as unique solution. $\;\qed$

3 Relation between duality function and stationary product measure

In the examples of duality that we have encountered in the previous section, we have a universal relation between the stationary product measures and the duality functions. Given $\{\eta(t),\ t\geq 0\}$, if there is a dual process $\{\xi(t),\ t\geq 0\}$ with factorized duality functions $D(\xi,\eta)=\prod d(\xi_{x},\eta_{x})$ as in (3)–(4) and stationary product measures $\{\mu_{\lambda}=\otimes\ \nu_{\lambda},\ \lambda>0\}$, then there is a relation between these measures and these functions, namely there exists a function $\theta(\lambda)$ such that

[TABLE]

This function $\theta(\lambda)$ is then simply the expectation of the first single-site duality function, i.e.

[TABLE]

In the examples of Section 2.4, we have $\theta(\lambda)=\lambda$ for IRW, $\theta(\lambda)=\tfrac{\lambda}{1-\lambda}$ for $\text{SIP}(\alpha)$ and $\theta(\lambda)=\tfrac{\lambda}{1+\lambda}$ for $\text{SEP}(\gamma)$.

In this section we first investigate under which general conditions this relation holds, and further use it in Sections 3.2.1–3.2.2 as a criterion of characterization of all extremal measures. The reader shall refer to Sections 2.1–2.2 for the general setting in which these results hold.

Later on, we will see that this relation (18) is actually a characterizing property of the factorized duality functions, meaning that all duality functions are determined once the first single-site duality function is fixed.

3.1 Finite case

We start with the simplest situation in which $V$ is a finite set.

First, we assume that the total number of particles/the total energy of the dual process is the only conserved quantity. More precisely, we assume the following property, which we refer to as harmonic triviality of the dual system:

1. [HT]

If $H:\hat{E}^{V}\to\mathbb{R}$ is harmonic, i.e. such that, for all $t>0$,

[TABLE]

then $H(\xi)$ is only a function of $|\xi|:=\sum_{x\in V}\xi_{x}$.

Moreover, let us assume the factorized form of the duality function $D(\xi,\eta)$ as in (3)–(4). Of the single-site functions $d$ we may require the following additional property:

1. [MD]

The function $d$ is measure determining, i.e. for two probability measures $\nu_{\ast}$, $\nu_{\star}$ on $E$ such that for all $x\in V$ and $\xi_{x}\in\hat{E}$

[TABLE]

it follows that $\nu_{\ast}=\nu_{\star}$.

Then we have the following.

THEOREM** 3.1****.**

Assume that $\{\eta(t),\ t\geq 0\}$ and $\{\xi(t),\ t\geq 0\}$ are dual as in (1) with factorized duality function (3) satisfying condition (4). Moreover, assume that [HT] holds and that $\mu$ is a probability measure on $\Omega$. We distinguish two cases:

1. (1)

Interacting particle system case.* If $\hat{E}$ is a subset of $\mathbb{N}$, then we assume that the self-duality functions $D(\xi,\cdot)$ are $\mu$-integrable for all $\xi\in\hat{\Omega}$.* 2. (2)

Interacting diffusion case.* If $\hat{E}=\mathbb{R}_{+}$, then we assume the following integrability condition: for each $\varepsilon>0$, there exists a $\mu$-integrable function $f_{\varepsilon}$ such that*

[TABLE]

Then

1. (a)

$\mu$* is a stationary product measure for the process $\{\eta(t):t\geq 0\}$*

implies

1. (b)

For all $\xi\in\Omega$ and for all $x\in V$, we have

[TABLE]

where $\delta_{x}$ denotes the configuration with a single particle at $x\in V$ and no particles elsewhere. Moreover, if condition [MD] holds, the two statements (a) and (b) are equivalent.

PROOF.  First assume that $\mu$ is a stationary product measure. Define $H(\xi)=\int D(\xi,\eta)\mu(d\eta)$. By $\mu$-integrability in the interacting particle system case (resp. (19) in the interacting diffusion case), self-duality and invariance of $\mu$ we have

[TABLE]

Therefore by [HT] we conclude that $H(\xi)=\psi(|\xi|)$. By using $d(0,\cdot)\equiv 1$ and the factorization of the duality functions, we have that $\psi(0)=1$. For the particle case, we obtain that

[TABLE]

In particular, we obtain that the l.h.s.  does not depend on $x$. Next, for $n\geq 2$, put

[TABLE]

then we have for $x\not=y\in V$, using the factorized duality function and the product form of the measure,

[TABLE]

from which it follows that $\psi(n)=\psi(1)^{n}$. Via an analogous reasoning that uses the factorization of $D(\xi,\eta)$ and the product form of $\mu$, for the diffusion case we obtain, for all $\varepsilon,\rho\geq 0$,

[TABLE]

and hence, by measurability of $\psi(\varepsilon)$, we get $\psi(\varepsilon)=\psi(1)^{\varepsilon}$.

To prove the other implication, put

[TABLE]

We then have by assumption

[TABLE]

and so it follows that $\mu$ is stationary by self-duality, $\mu$-integrability, the conservation of the number of particles and the measure-determining property. Indeed,

[TABLE]

From the factorized form of $D(\xi,\eta)$, (20) implies that for all $x\in V$ and $\xi_{x}\in\hat{E}$

[TABLE]

and also

[TABLE]

therefore $\mu$ is a product measure by the fact that $d$ is measure determining. $\;\qed$

3.2 Infinite case

If $V=\mathbb{Z}^{d}$, then one needs essentially two extra conditions to state an analogous result in which a general relation between duality functions and corresponding stationary measures can be derived.

In this section we will assume that the dual process is a discrete particle system, i.e. $\hat{E}$ is a subset of $\mathbb{N}$, in which the number of particles is conserved. In this case we need an additional property ensuring that for the dynamics of a finite number of particles there are no bounded harmonic functions other than those depending on the total number of particles. Therefore, we introduce the condition of existence of a successful coupling for the discrete dual process with a finite number of particles. This is defined below.

DEFINITION** 3.1****.**

We say that the discrete dual process $\{\xi(t),\ t\geq 0\}$ has the successful coupling property when the following holds: if we start with $n$ particles then there exists a labeling such that for the corresponding labeled process $\{X_{1}(t),\ldots,X_{n}(t),\ t\geq 0\}$ there exists a successful coupling. This means that for every two initial positions ${\bf x}=(x_{1},\ldots,x_{n})$ and ${\bf y}=(y_{1},\ldots,y_{n})$, there exists a coupling with path space measure $\mathbb{P}_{\bf x,y}$ such that the coupling time

[TABLE]

is finite $\mathbb{P}_{\bf x,y}$ almost surely.

Notice that the successful coupling property is the most common way to prove the following property (cf. [13]), which is the analogue of [HT], referred here to as bounded harmonic triviality of the dual process:

1. [BHT]

If $H$ is a bounded harmonic function, then $H(\xi)=\psi(|\xi|)$ for some bounded $\psi:\hat{E}\to\mathbb{R}$.

REMARK** 3.1****.**

The condition of the existence of a successful coupling (and the consequent bounded harmonic triviality) is quite natural in the context of interacting particle systems, where we have that a finite number of walkers behave as independent walkers, except when they are close and interact. Therefore, the successful coupling needed is a variation of the Ornstein coupling of independent walkers, see e.g. [3], [5], [11].

Furthermore, we need a form of uniform $\mu$-integrability of the duality functions which we introduce below and call uniform domination property of $D$ w.r.t. $\mu$ (note the analogy with condition (19)):

1. [UD]

Given $\mu$ a probability measure on $\Omega$, the duality functions $\{D(\xi,\cdot),\ |\xi|=n\}$ are uniformly $\mu$-integrable, i.e. for all $n\in\mathbb{N}$ there exists a function $f_{n}$ such that $f_{n}$ is $\mu$-integrable and such that for all $\eta\in\Omega$

[TABLE]

Under these conditions, the following result holds, whose proof resembles that of Theorem 3.1.

THEOREM** 3.2****.**

Assume as in (1) that $\{\eta(t),\ t\geq 0\}$ is dual to the discrete process $\{\xi(t),\ t\geq 0\}$ with factorized duality function as in (3)–(4). Moreover, assume [BHT] in place of [HT] for the dual process and that $\mu$ is a probability measure on $\Omega$ such that [UD] holds. Then the same conclusions as in Theorem 3.1 follow, where (20) holds for all finite configurations $\xi\in\hat{\Omega}$.

3.2.1 Translation invariant case

In this section, we show that under the assumption of factorized duality functions, minimal ergodicity assumptions on a stationary probability measure $\mu$ on $\Omega$ are needed to ensure (20) and, as a consequence, that $\mu$ is product measure.

Here we restrict to the case $V=\mathbb{Z}^{d}$ because we will use spatial ergodicity.

THEOREM** 3.3****.**

In the setting of Theorem 3.2 with $D(\xi,\eta)$ factorized duality function and $\mu$ probability measure on $\Omega$, if $\mu$ is a translation-invariant and ergodic (under translations) stationary measure for $\{\eta(t),\ t\geq 0\}$, then we have (20) for all finite configurations $\xi$; as a consequence, $\mu$ is a product measure.

PROOF.  To start, let us consider a configuration $\xi=\sum_{i=1}^{n}\delta_{x_{i}}$. By bounded harmonic triviality [BHT], combined with the bound (21) for all such configurations, $\int D(\xi,\eta)\mu(d\eta)$ is only depending on $n$ and, therefore, we can replace $\xi$ by $\sum_{i=1}^{n-1}\delta_{x_{i}}+\delta_{y}$, where $y$ is arbitrary in $\mathbb{Z}^{d}$. Let us call $B_{N}=[-N,N]^{d}\cap\mathbb{Z}^{d}$. Fix $N_{0}$ such that $B_{N_{0}}$ contains all the points $x_{1},\ldots x_{n-1}$. For $y$ outside $B_{N_{0}}$, by the factorization property, $D(\sum_{i=1}^{n-1}\delta_{x_{i}}+\delta_{y},\eta)=D(\sum_{i=1}^{n-1}\delta_{x_{i}},\eta)D(\delta_{y},\eta)$. By the Birkhoff ergodic theorem, we have that

[TABLE]

$\mu$-a.s. as $N\to\infty$. Using this, together with (21), we have

[TABLE]

Iterating this argument gives (20). $\;\qed$

REMARK** 3.2****.**

As follows clearly from the proof, the condition of factorization of the duality function can be replaced by the weaker condition of

[TABLE]

for all $\eta$, $x_{1},\ldots,x_{n}$. We note that this approximate factorization of the duality function leads to (20), though $\mu$ is not necessarily a product measure.

3.2.2 Non-translation invariant case

We continue here with $V=\mathbb{Z}^{d}$ but drop the assumption of translation invariance. Indeed, equality (20) is also valid in contexts where one cannot rely on translation invariance. Examples include spatially inhomogeneous $\text{SIP}(\boldsymbol{\alpha})$ and $\text{SEP}(\boldsymbol{\gamma})$, where the parameters $\boldsymbol{\alpha}=(\alpha_{x})_{x\in V}$ and $\boldsymbol{\gamma}=(\gamma_{x})_{x\in V}$ in Section 2.3 may depend on the site accordingly (cf. e.g. [14]). Also in this inhomogeneous setting the self-duality functions factorize over sites and the stationary measures are in product form, with site-dependent single-site duality functions, resp. site-dependent marginals. We will show that the relation (20) between the self-duality functions and any ergodic stationary measure still holds, and as a consequence this ergodic stationary measures is in fact a product measure. The idea is that the averaging over space w.r.t. $\mu$, used in the proof of Theorem 3.3 above, can be replaced by a time average.

If we start with a single dual particle, the dual process is a continuous-time random walk on $V$, for which we denote by $p(t;x,y)$ the transition probability to go from $x$ to $y$ in time $t>0$. A basic assumption will then be

[TABLE]

for all $x,y\in V$.

THEOREM** 3.4****.**

In the setting of Theorem 3.2 with $D(\xi,\eta)$ duality function and $\mu$ probability measure on $\Omega$, if $\mu$ is an ergodic stationary measure for the process $\{\eta(t),\ t\geq 0\}$ and (23) holds for the dual particle, then we have (20) for all finite configurations $\xi$; as a consequence, $\mu$ is a product measure.

PROOF.  The idea is to replace the spatial average in the proof of Theorem 3.3 by a Cesaro average over time, which we can deal by combining assumption (23) with the assumed temporal ergodicity.

Fix $x_{1},\ldots,x_{n}\in V$, $y\in V$. Define

[TABLE]

It is sufficient to obtain that $H_{n+1}(x_{1},\ldots,x_{n},y)=H_{1}(y)H_{n}(x_{1},\ldots,x_{n})$. We already know by the bounded harmonic triviality that $H_{n}$ only depends on $n$ and not on the given locations $x_{1},\ldots,x_{n}$. Therefore, we have

[TABLE]

By assumption (23), this implies

[TABLE]

where in the last step we used the assumed temporal ergodicity of $\mu$ and Birkhoff ergodic theorem. $\;\qed$

4 From stationary product measures to duality functions.

As we have just illustrated in Sections 3.2.1–3.2.2, relation (18) turns out to be useful in deriving information about the product structure of stationary ergodic measures from the knowledge of factorized duality functions. On the other side, granted some information on the stationary product measures, which follows usually from a simple detailed balance computation, up to which extent does relation (18) say something about the possible factorized duality functions?

In the context of conditions (3)–(4) and in presence of a one-parameter family of stationary product measures $\{\mu_{\lambda}=\otimes\ \nu_{\lambda},\ \lambda>0\}$, relation (18) for $\xi=k\delta_{x}\in\hat{\Omega}$ for some $k\in\mathbb{N}$ reads

[TABLE]

As a consequence, knowing the first single-site duality function $d(1,\cdot)$ and the explicit expression of the marginal $\nu_{\lambda}$ is enough to recover the l.h.s. in (24). However, rather than obtaining $d(k,\cdot)$, at this stage the l.h.s. has still the form of an “integral transform”-type of expression for $d(k,\cdot)$.

In the next two subsections, we show how to recover $d(k,\eta_{x})$ from (24) and the knowledge of $\theta(\lambda)$. This then leads to the characterization of all possible factorized (self-)duality functions.

4.1 Particle systems and orthogonal polynomial self-duality functions

Going back to the interacting particle systems introduced in Section 2.3 with infinitesimal generator (5) and product stationary measures with marginals (7), the integral relation (24) rewrites, for each $k\in\mathbb{N}$ and $\lambda>0$ for which $Z_{\lambda}<\infty$, as

[TABLE]

where $\varphi(n)=n!\prod_{m=1}^{n}\tfrac{v(m-1)}{u(m)}$. Now, if we interpret

[TABLE]

as the Taylor series expansion around $\lambda=0$ of the function $\theta(\lambda)^{k}Z_{\lambda}$, we can re-obtain the explicit formula of $d(k,n)\varphi(n)$ as its $n$-th order derivative evaluated at $\lambda=0$, namely

[TABLE]

and hence, anytime $\varphi(n)>0$,

[TABLE]

Together with the full characterization obtained in Theorem 2.1 of the first single-site self-duality functions $d(1,\cdot)$ - and $\theta(\lambda)$ in turn - we obtain via this procedure a full characterization of all single-site self-duality functions. Beside recovering the “classical” dualities illustrated in Section 2.3, we find the single-site self-duality functions in terms of orthogonal polynomials $\{p_{k}(n),\ k\in\mathbb{N}\}$ of a discrete variable (cf. e.g. [15]) recently discovered via a different approach in [8]. We add the observation that all these new single-site self-duality functions can be obtained from the classical ones via a Gram-Schmidt orthogonalization procedure w.r.t. the correct probability measures on $\mathbb{N}$, namely the marginals of the associated stationary product measures (cf. [8]).

We divide the discussion in three cases, one suitable for processes of IRW-type, the other for SIP and SEP and the last one for the remaining particle systems.

4.1.1 Independent random walkers

We recall that the IRW-case corresponds to the choice of values in (2.1) satisfying the relation $v(1)=v(0)$. If we compute $\theta(\lambda)$ for the general first single-site self-duality function $d(1,n)=a+bn$ obtained in (9), we get

[TABLE]

and, in turn via relation (25), we recover all functions $d(k,\cdot)$, for $k>1$:

[TABLE]

In case $a=0$, $d(k,n)=0$ for $n<k$, while for $n\geq k$, in the summation all terms but the one corresponding to $r=k$ vanish, thus $d(k,n)=\frac{n!}{(n-k)!}b^{k}$. In case $a\neq 0$,

[TABLE]

In conclusion, for the choice $a\cdot b<0$,

[TABLE]

where $\{C_{k}(n;\mu),\ k\in\mathbb{N}\}$ are the Poisson-Charlier polynomials - orthogonal polynomials w.r.t. the Poisson distribution of parameter $\mu>0$ (cf. [15]).

4.1.2 Inclusion and exclusion processes

For SIP and SEP we are in the case $v(1)\neq v(0)$, and hence we abbreviate

[TABLE]

where for $\text{SIP}(\alpha)$ we choose $\sigma=1$ and $\beta=\alpha$, while for $\text{SEP}(\gamma)$ we set $\sigma=-1$ and $\beta=\gamma$. If we compute $\theta(\lambda)$ for $d(1,n)=a+bn$ in (9), we have

[TABLE]

By applying formula (25), we obtain all functions $d(k,n)$ for $k>1$ as follows:

[TABLE]

In case $a=0$, clearly $d(k,n)=0$ for $k<n$, while for $n\geq k$ only the term for $r=k$ is nonzero in the summation:

[TABLE]

In case $a\neq 0$, by using the known relation (cf. [15], p. 51),

[TABLE]

If $\sigma=1$, $\beta=\alpha>0$ and if $a\cdot b<0$, we recognize in (28) the Meixner polynomials as defined in [15], i.e.

[TABLE]

where in our case $\mu=\frac{a}{a-b\alpha}$ and $\{M_{k}(n;\mu,\alpha),\ k\in\mathbb{N}\}$ are the Meixner polynomials - orthogonal polynomials w.r.t. the discrete Gamma distribution of scale parameter $\mu$ and shape parameter $\alpha$.

Furthermore, if $\sigma=-1$, $\beta=\gamma$ with $\gamma\in\mathbb{N}$ and given the additional requirements $a\cdot b<0$ and $-\frac{a}{b}\leq\gamma$, from the expression in (28) we have a representation of the single-site duality functions in terms of the Kravchuk polynomials as defined in [15], i.e.

[TABLE]

where $p=-\frac{a}{b\gamma}$ in our case and $\{K_{k}(n;p,\gamma),\ k\in\mathbb{N}\}$ the Kravchuk polynomials - orthogonal polynomials w.r.t. the Binomial distribution of parameters $\gamma$ and $p$.

As a conclusion of this procedure, we note that all factorized self-duality functions for independent random walkers, inclusion and exclusion processes satisfying (4) are either in the “classical” form of Section 2.3 (case $a=0$) or consist of products of rescaled versions of orthogonal polynomials (case $a\neq 0$). Other factorized self-duality functions for the systems IRW, SIP and SEP do not exist.

REMARK** 4.1****.**

It is interesting to note that, apart from the leading factor $a^{k}$, the remaining polynomials in the expressions of $d(k,n)$ for $a\neq 0$ are “self-dual” in the sense of the orthogonal polynomials literature, i.e. $p_{n}(k)=p_{k}(n)$ in our context (cf. Definition 3.1, [10]). Henceforth, if $d(k,n)$ is interpreted as a countable matrix, the value $a\in\mathbb{R}$ is the only responsible for the asymmetry of $d(k,n)$: upper-triangular for $a=0$ while symmetric for $a=1$.

4.1.3 Trivial factorized self-duality

To conclude, for the sake of completeness, we can implement the same machinery to cover all factorized self-dualities with property (4) for all discrete processes of type (5).

Indeed, from the proof of Theorem 2.1, if the process is neither of the types IRW, SIP and SEP, then the only possible choice is $d(1,n)=a$ for some $a\in\mathbb{R}$, i.e. it is not depending on $n$. From this we get $\theta(\lambda)=a$, and $d(k,n)=a^{k}$ from formula (25). Hence, the self-duality functions must be of the form

[TABLE]

i.e. depending only on the total number of dual particles (and not on the configuration $\eta$). Hence, the duality relation in that case reduces to the trivial relation, for all $t\geq 0$ and $\xi\in\hat{\Omega}$,

[TABLE]

which is just conservation of the number of particles in the dual process. No other self-duality relation with factorized self-duality functions can exist.

4.2 Interacting diffusions and orthogonal polynomial duality functions

As shown in Theorem 3.1, relation (24) still holds whenever the discrete right-variables $n\in\mathbb{N}$ are replaced by continuous variables $z\in\mathbb{R}_{+}$ and sums by integrals. With this observation in mind, we provide a second general method to characterize all factorized duality functions between the continuous process $\text{BEP}(\alpha)$ and its discrete dual $\text{SIP}(\alpha)$.

More precisely, if $d(k,z)$ is a single-site duality function with property (4) between $\text{BEP}(\alpha)$ and $\text{SIP}(\alpha)$, and $\nu_{\lambda}$ is the stationary product measure marginal for $\text{BEP}(\alpha)$ as in (16), then, from the analogue of relation (24) for $k=1$, namely

[TABLE]

we necessarily have by Theorem 3.1 that

[TABLE]

As a consequence, the function $d(k,z)\frac{z^{\alpha-1}}{\Gamma(\alpha)}$ is the inverse Laplace transform of $\theta(\lambda)^{k}\lambda^{-\alpha}$. Given the first single-site duality function $d(1,z)$ in (17), from (29) we obtain

[TABLE]

As a consequence, the r.h.s. in (30) becomes

[TABLE]

and there exist explicit expressions for the inverse Laplace transform of this function. We split the computation in two cases. In case $a=0$, since the inverse Laplace transform of $\lambda^{-\alpha-k}$ is $\frac{z^{\alpha+k-1}}{\Gamma(\alpha+k)}$, we immediately obtain

[TABLE]

i.e. the “classical” single-site duality function as in Section 2.5, up to set $b=\frac{1}{\alpha}$. In case $a\neq 0$, the inverse Laplace transform of (31) is more elaborated:

[TABLE]

As the above expression must equal $d(k,z)\frac{z^{\alpha-1}}{\Gamma(\alpha)}$, it follows that

[TABLE]

As a final consideration, we note that for the choice $a\cdot b<0$,

[TABLE]

where $\mu=-\frac{b\alpha}{a}$ here and $\{L_{k}(z;\alpha-1,\mu),\ k\in\mathbb{N}\}$ are the generalized Laguerre polynomials - orthogonal polynomials w.r.t. to the Gamma distribution of scale parameter $\frac{1}{\mu}$ and shape parameter $\alpha$ as defined in [15].

5 Intertwining and generating functions

In this section, we introduce the generating function method, which allows to go from a self-duality of a discrete process towards duality between a discrete and continuous process, and further towards a self-duality of a continuous process, and back. This then allows e.g. to simplify the proof of a discrete self-duality by lifting it to a continuous self-duality, which is usually easier to verify. The key of this method is an intertwining between discrete multiplication and derivation operators and their continuous analogues, via an appropriate generating function.

To reduce issues of well-definition of operators and their domains, in this section we restrict our discussion to the case of finite vertex set $V$.

5.1 Introductory example

To make the method clear, let us start with a simple example of independent random walkers on a single edge. The generator is

[TABLE]

with $n_{1},n_{2}\in\mathbb{N}$. Define now the (exponential) generating function

[TABLE]

Then it is easy to see that

[TABLE]

where

[TABLE]

Now assume that we have a self-duality function for the particle system, i.e.

[TABLE]

then we have

[TABLE]

where

[TABLE]

In words, a self-duality function of $L$ is “lifted” to a duality function between the independent random walk generator $L$ and its continuous counterpart ${\mathscr{L}}$ by applying the generating function to the $n$-variables. Conversely, given a duality function between the independent random walk generator $L$ and its continuous counterpart ${\mathscr{L}}$, its Taylor coefficients provide a self-duality function of $L$.

We can then also take the generating function w.r.t. the $k$-variables in the function ${\mathcal{D}}$ to produce a self-duality function for ${\mathscr{L}}$, i.e. defining

[TABLE]

we have

[TABLE]

For the classical self-duality function $D(k_{1},k_{2};n_{1},n_{2})=\frac{n_{1}!}{(n_{1}-k_{1})!}\frac{n_{2}!}{(n_{2}-k_{2})!}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{k_{1}\leq n_{1}\}{\mathchoice{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.0mul}{\rm 1\mskip-4.5mul}{\rm 1\mskip-5.0mul}}\{k_{2}\leq n_{2}\}$, we find that

[TABLE]

Beside the factor $e^{z_{1}+z_{2}}$ which depends only on the conserved quantity $z_{1}+z_{2}$, to check the self-duality relation for ${\mathscr{L}}$ w.r.t. the function $e^{v_{1}z_{1}+v_{2}z_{2}}$ is rather straightforward, the computation involving only derivatives of exponentials. By looking at the Taylor coefficients w.r.t. both $v$ and $z$-variables of this self-duality relation, we obtain the self-duality relation for $L$ w.r.t. $D$ where we started from.

In conclusion, all these duality relations turn out to be equivalent, and the proof of self-duality for particle systems requiring rather intricate combinatorial arguments (cf. e.g. [3]) is superfluous once the more direct self-duality for diffusion systems is checked.

5.2 Intertwining and duality

The notion of intertwining between stochastic processes was originally introduced by Yor in [16] in the context of Markov chains and later pursued in [4] and [6] as an abstract framework, in discrete-time and continuous-time respectively, for the problem of Markov functionals, i.e. finding sufficient and necessary conditions under which a random function of a Markov chain is again Markovian.

For later purposes, we adopt a rather general definition of intertwining, in which $\{\eta(t),\ t\geq 0\}$ and $\{\zeta(t),\ t\geq 0\}$ are continuous-time stochastic processes on the Polish spaces $\Omega$ and $\Omega^{\prime}$, respectively, whose expectations read $\mathbb{E}$, $\mathbb{E}^{\prime}$ resp., and ${\mathcal{M}}(\Omega)$ denotes the space of signed measures on $\Omega$. We say that $\{\zeta(t),\ t\geq 0\}$ is intertwined on top of $\{\eta(t),\ t\geq 0\}$ if there exists a mapping $\varLambda:\Omega^{\prime}\to{\mathcal{M}}(\Omega)$ such that, for all $t\geq 0$, $\zeta\in\Omega^{\prime}$ and $f:\Omega\to\mathbb{R}$,

[TABLE]

Working at the abstract level of semigroups, we say that $\{{\mathscr{S}}(t),\ t\geq 0\}$ on a space of functions $f:\Omega^{\prime}\to\mathbb{R}$ denoted by ${\mathcal{F}}(\Omega^{\prime})$, is intertwined on top of $\{S(t),\ t\geq 0\}$, a semigroup on a space of functions $f:\Omega\to\mathbb{R}$ denoted by ${\mathcal{F}}(\Omega)$, with intertwiner $\varLambda$ if $\varLambda$ is a linear operator from ${\mathcal{F}}(\Omega)$ into ${\mathcal{F}}(\Omega^{\prime})$ and if, for all $t\geq 0$ and $f:\Omega\to\mathbb{R}$,

[TABLE]

Similarly, operators ${\mathscr{L}}$ with domain ${\mathcal{D}}({\mathscr{L}})$ and $L$ with domain ${\mathcal{D}}(L)$ are intertwined with intertwiner $\varLambda$ if, for all $f\in{\mathcal{D}}(L)$, $\varLambda f\in{\mathcal{D}}({\mathscr{L}})$ and

[TABLE]

Notice that with a slight abuse of notation we used the same symbol $\varLambda$ for an abstract intertwining operator as for the intertwining mapping. In other words, in case the intertwining mapping as in (33) is given by $\tilde{\varLambda}$, then the corresponding operator is

[TABLE]

An intertwining mapping $\varLambda$ has a probabilistic interpretation if it takes values in the subset of probability measures on $\Omega$. Indeed, in (33) the process $\{\zeta(t),\ t\geq 0\}$ may be viewed as an added structure on top of $\{\eta(t),\ t\geq 0\}$ or, alternatively, the process $\{\eta(t),\ t\geq 0\}$ as a random functional of $\{\zeta(t),\ t\geq 0\}$, in which $\varLambda$ provides this link.

REMARK** 5.1****.**

The connection with duality introduced in Section 2.1 becomes transparant when $\Omega$, $\hat{\Omega}$ and $\Omega^{\prime}$ are finite sets and the operators and functions $L$, $\hat{L}$ ${\mathscr{L}}$, $D$ and $\varLambda$ in (2) and (34) are represented in terms of matrices. There, relations (2) and (34), once rewritten in matrix notation as

[TABLE]

where $L^{\dagger}$ denotes the transpose of $L$, and

[TABLE]

differ essentially only in the terms $L^{\dagger}$ versus $L$ in the r.h.s. of both identities. The presence or absence of transposition can be interpreted as a forward-versus-backward evolution against a forward-versus-forward evolution. More precisely, if $L$, $\hat{L}$ and ${\mathscr{L}}$ are generators of Markov processes $\{\eta(t),t\geq 0\}$, $\{\xi(t),\ t\geq 0\}$ and $\{\zeta(t),\ t\geq 0\}$, respectively, then (35) and (36) relate the evolution of $\{\eta(t),\ t\geq 0\}$ to that of $\{\xi(t),\ t\geq 0\}$, resp. $\{\zeta(t),\ t\geq 0\}$; however, while in (36) the processes run both along the same direction in time, in (35) the processes run along opposite time directions.

Intertwiners as $\varLambda$ in (36) may be also interpreted as natural generalizations of symmetries of generators, indeed (36) with ${\mathscr{L}}=L$ just means that $\Lambda$ commutes with $L$, which is the definition of a symmetry of $L$. As outlined in [9, Theorem 2.6], the knowledge of symmetries of a generator and dualities of this generator leads to the construction of new dualities. The following theorem presents the analogue procedure in presence of intertwiners: a duality and an intertwining lead to a new duality.

THEOREM** 5.1****.**

Let $L$, $\hat{L}$ and ${\mathscr{L}}$ be operators on real-valued functions on $\Omega$, $\hat{\Omega}$ and $\Omega^{\prime}$, respectively. Suppose that there exists an intertwiner $\varLambda$ such that for all $f\in{\mathcal{D}}(L)$, $\varLambda f\in{\mathcal{D}}({\mathscr{L}})$,

[TABLE]

and a duality function $D:\hat{\Omega}\times\Omega\to\mathbb{R}$ for $\hat{L}$ and $L$, namely $D(\xi,\cdot)\in{\mathcal{D}}(L)$ for all $\xi\in\hat{\Omega}$, $D(\cdot,\eta)\in{\mathcal{D}}(\hat{L})$ for all $\eta\in\Omega$ and

[TABLE]

Then, if $\varLambda_{\text{\tiny right}}D(\xi,\cdot)\in{\mathcal{D}}({\mathscr{L}})$ for all $\xi\in\hat{\Omega}$ and $\varLambda_{\text{\tiny right}}D(\cdot,\zeta)\in{\mathcal{D}}(\hat{L})$ for all $\zeta\in\Omega^{\prime}$, $\varLambda_{\text{\tiny right}}D$ is a duality function for $\hat{L}$ and ${\mathscr{L}}$, i.e.

[TABLE]

PROOF.

[TABLE]

Here in the first equality we used that left and right actions commute, in the second equality we used the assumed duality of $\hat{L}$ and $L$, and in the third equality we used the assumed intertwining. $\;\qed$

5.3 Intertwining between continuum and discrete processes

In this section we prove the existence of an intertwining relation between the interacting diffusion processes presented in Section 2.5 and the particle systems of Section 2.3. This intertwining relation provides a second connection, besides the scaling limit procedure (cf. Section 2.5), between continuum and discrete processes, which proves to be better suited for the goal of establishing duality relations among these processes. Indeed, the characterization of all possible factorized self-dualities for particle systems obtained in Section 4 and the intertwining relation below, via the application of Theorem 5.1, produces a characterization of all possible dualities, resp. self-dualities, between the discrete and the continuum processes, resp. of the continuum process.

In the following proposition, we prove the intertwining relation for operators $L^{\sigma,\beta}$ and ${\mathscr{L}}^{\sigma,\beta}$ defined, respectively, on functions $f:E^{V}\to\mathbb{R}$, $E=\mathbb{N}$, as

[TABLE]

where

[TABLE]

and, on real analytic functions $f:\mathbb{R}^{V}\to\mathbb{R}$, as

[TABLE]

where

[TABLE]

Note that $L^{\sigma,\beta}$ in (37) is a special instance of the generator $L^{u,v}$ in (5) with conditions (2.1), while ${\mathscr{L}}^{\sigma,\beta}$ above, matches on a common sub-domain, for particular choices of the parameters $\sigma$ and $\beta$, those in (13).

PROPOSITION** 5.1****.**

Let $G$ be the Poisson probability kernel defined as the operator that maps functions $f:\mathbb{N}\to\mathbb{R}$ into functions $Gf:\mathbb{R}\to\mathbb{R}$ as

[TABLE]

Then, whenever $Gf:\mathbb{R}\to\mathbb{R}$ is a real analytic function, if $G^{\otimes}=\otimes_{x\in V}G^{(x)}$ denotes the tensorized operator mapping functions $f:\mathbb{N}^{V}\to\mathbb{R}$ into functions $f:\mathbb{R}^{V}\to\mathbb{R}$ accordingly, ${\mathscr{L}}^{\sigma,\beta}$ and $L^{\sigma,\beta}$ are intertwined with intertwiner $G^{\otimes}$, namely

[TABLE]

PROOF.  Let us introduce the non-normalized operator

[TABLE]

and the associated tensorized operator $\bar{G}^{\otimes}=\otimes_{x\in V}\bar{G}^{(x)}$. Due to the factorized structure of $L^{\sigma,\beta}$, ${\mathscr{L}}^{\sigma,\beta}$ and $\bar{G}^{\otimes}$, the proof of the intertwining relation (40) with $\bar{G}^{\otimes}$ as an intertwiner reduces to consider and combine the following relations:

[TABLE]

As a first consequence, we have

[TABLE]

We obtain (40) by observing that ($|z|=\sum_{x\in V}z_{x}$)

[TABLE]

and that, for $g(z)=\bar{g}(z)\cdot e^{-|z|}$,

[TABLE]

$\;\qed$

We note that the intertwiner $G^{\otimes}$ has a nice probabilistic interpretation: from an “energy” configuration $z\in\mathbb{R}_{+}^{V}$, the associated particle configurations are generated by placing a number of particles on each site $x\in V$, independently, and distributed according to a Poisson random variable with intensity $z_{x}$.

In the remaining part of this section, under some reasonable regularity assumptions, we are able to invert the intertwining relation (40), namely to find an operator $H^{\otimes}=\otimes_{x\in V}H^{(x)}$ that intertwines $L^{\sigma,\beta}$ and ${\mathscr{L}}^{\sigma,\beta}$, in this order. The natural candidate for $H$ is the “inverse operator” of $G$, whenever this is well-defined. In general, this “inverse intertwiner” lacks any probabilistic interpretation, but indeed establishes a second intertwining relation useful in the next section.

PROPOSITION** 5.2****.**

Let $H$ be the differential operator mapping real analytic functions $g:\mathbb{R}\to\mathbb{R}$ into functions $Hg:\mathbb{N}\to\mathbb{R}$ as

[TABLE]

Then $H$ is the inverse operator of $G$, namely, for all $f:\mathbb{N}\to\mathbb{R}$ such that $Gf:\mathbb{R}\to\mathbb{R}$ is real analytic, we have

[TABLE]

Moreover, the tensorized operator $H^{\otimes}=\otimes_{x\in V}H^{(x)}$ is an intertwiner for $L^{\sigma,\beta}$ and ${\mathscr{L}}^{\sigma,\beta}$, i.e. for all real analytic $g:\mathbb{R}^{V}\to\mathbb{R}$,

[TABLE]

Before giving the proof, we need the following lemma.

LEMMA** 5.1****.**

Let $A$ be the operator acting on functions $f:\mathbb{N}\to\mathbb{R}$ defined as

[TABLE]

Then, the tensorized operator $A^{\otimes}=\otimes_{x\in V}A^{(x)}$ is a symmetry for the generator $L^{\sigma,\beta}$, i.e. for all $f\ :\mathbb{N}^{V}\to\mathbb{R}$

[TABLE]

PROOF.  Instead of going through tedious computations, we exploit the fact that the operator $A^{\otimes}$ has the form

[TABLE]

where $J^{(x)}$ is an operator defined for functions $f:\mathbb{N}^{V}\to\mathbb{R}$ which acts only on the $x$-th variable as

[TABLE]

Since all these operators $\{J^{(x)},\ x\in V\}$ commute over the sites we have

[TABLE]

We conclude the proof by noting that the operator

[TABLE]

is a symmetry for the generator $L^{\sigma,\beta}$, cf. e.g.  [9], [1]. $\;\qed$

PROOF.  [Proposition 5.2] First we compute the following key relations:

[TABLE]

Hence, if we introduce the operator

[TABLE]

and the associated tensorized operator $\bar{H}^{\otimes}=\otimes_{x\in V}\bar{H}^{(x)}$, we obtain

[TABLE]

Now, by using Lemma 5.1 and noting that

[TABLE]

and, by the mixed property of the tensor product,

[TABLE]

we get (42) by applying first (45), then (43) and finally (44):

[TABLE]

$\;\qed$

5.4 Generating functions and duality

As anticipated in the previous section, from the intertwining relation (40) and the functions obtained in Section 4.1, in what follows we find explicitly new duality relations.

Due to the factorized form (3) of the self-duality functions with single-site functions (26) and (28) and the tensor form of the intertwiner $G^{\otimes}$ in (40), the new functions inherit the same factorized form. Moreover, from the definition of $G$ in (39), the whole computation reduces to determine (exponential) generating functions of (26) and (28). To this purpose, some identities for hypergeometric functions are available, cf. e.g.  the tables in [10, Chapter 9]. Moreover, all generating functions obtained satisfy the requirements of analyticity for suitable choices of the parameters $\sigma$, $\beta$, $a$ and $b$ (cf. [10]), hence all operations below make sense.

However, just as the functions found in Section 4.1, the functions here obtained will only be “candidate” (self-)duality functions, since no duality relation as in (1) has been proved, yet. By using the “inverse”  intertwining (42), all these “possible” dualities turn out to be equivalent, i.e. one implies all the others. Thus, in Proposition 5.3 below, we choose to prove directly the self-duality relation for the continuum process, more immediate to verify due to the simpler form of the self-duality functions. Indeed, while the single-site self-duality functions for the $\text{SIP}(\alpha)$ process, for instance, have the generic form of an hypergeometric function

[TABLE]

the single-site duality functions between discrete and continuum processes involve in their expressions hypergeometric functions

[TABLE]

and those for the self-duality of continuum processes are even simpler, namely

[TABLE]

as the number of arguments of the hypergeometric function drops.

The tables below schematically report all single-site (self)-duality functions for the operators $L^{\sigma,\beta}$ in (37) and ${\mathscr{L}}^{\sigma,\beta}$ in (38). Recall that the parameters $a$, $b\in\mathbb{R}$ in (2.1) are properly chosen (cf. Section 4.1).

(I) Independent random walkers (IRW) , \sigma=0\, \beta=1\.

[TABLE]

(II) Symmetric inclusion process ($\text{SIP}(\alpha)$) , \sigma=1\, \beta=\alpha>0\.

[TABLE]

(III) Symmetric exclusion process ($\text{SEP}(\gamma)$) , \sigma=-1\, \beta=\gamma\in\mathbb{N}\.

[TABLE]

More in detail, on the left-most column we place the single-site self-duality functions $d(k,n)$ for the particle systems of Section 2.3: while the top-left functions are those already appearing in [9], [1], cf.  also Section 2.3 and (27) - and, thus, for this reason denoted here as the “classical” ones - the second-to-the-top functions are those derived in Section 4.1 in (26)–(28) and being related to suitable families of orthogonal polynomials. While these two classes of single-site self-duality functions satisfy condition (4) (they are the only ones doing so), the bottom-left single-site self-duality functions correspond to the “cheap” self-duality (cf. end of Section 2.2), namely the detailed-balance condition w.r.t. the measures $\{\otimes_{x\in V}\nu_{\lambda},\ \lambda>0\}$ with marginals (7).

On the mid-column, we find the single-site duality functions between the difference operators $L^{\sigma,\beta}$ and the differential operators ${\mathscr{L}}^{\sigma,\beta}$, obtained from their left-neighbors by a direct application of the operator $G$ in (39) on the $n$-variables. The new functions will depend hence on the two variables $k\in\mathbb{N}$ and $z\in\mathbb{R}_{+}$.

A second application w.r.t. the $k$-variables of the same operator $G$ on the functions just obtained gives us back the right-most column, functions depending now on variables $v,z\in\mathbb{R}_{+}$. These functions represent the single-site self-duality functions for the differential operator ${\mathscr{L}}^{\sigma,\beta}$. As an immediate consequence of Proposition 5.2, we could also proceed from right to left by applying the inverse intertwiner $H$ in (41).

Note that the single-site self-duality functions for ${\mathscr{L}}^{\sigma,\beta}$ on the right-most columns, though they have been derived from different discrete analogues, i.e. classical, orthogonal and cheap single-site functions, within the same table they differ only for a factor which depends only on the conserved quantities $|z|=\sum_{x\in V}z_{x}$ and $|v|=\sum_{x\in V}v_{x}$. Henceforth, when proving the self-duality relation, this extra-factor does not play any role and it is enough to check that the functions

[TABLE]

for constants $c\in\mathbb{R}$, are single-site self-duality functions for the operators ${\mathscr{L}}^{0,1}$, ${\mathscr{L}}^{1,\alpha}$ and ${\mathscr{L}}^{-1,\gamma}$, respectively. This final computation is the content of the next proposition.

PROPOSITION** 5.3****.**

For any constant $c\in\mathbb{R}$, the functions $d(v,z)$ in (46) are single-site self-duality functions for the differential operators ${\mathscr{L}}^{0,1}$, ${\mathscr{L}}^{1,\alpha}$ and ${\mathscr{L}}^{-1,\gamma}$, respectively.

PROOF.  To prove that

[TABLE]

is a single-site self-duality function for the differential operator ${\mathscr{L}}^{0,1}$, we first observe that

[TABLE]

Hence, the self-duality relation for the single-edge generator ${\mathscr{L}}^{1,0}_{x,y}$ rewrites

[TABLE]

which indeed holds.

For the second proof of self-duality for the single-site function

[TABLE]

we use the following shortcut: for $x\in V$,

[TABLE]

Additionally, we recall a formula for the $z_{x}$-derivative of $F_{x}$, namely

[TABLE]

and a recurrence identity

[TABLE]

Hence, the self-duality relation for ${\mathscr{L}}^{1,\alpha}_{x,y}$ w.r.t. the function $F_{x}(\alpha)F_{y}(\alpha)$ rewrites by using (47) as

[TABLE]

By substituting (48), the identity holds.

The proof for the operator ${\mathscr{L}}^{-1,\gamma}$ follows the same lines and we omit it. $\;\qed$

Acknowledgments

F.R. thanks Cristian Giardinà and Gioia Carinci for precious discussions; F.S. thanks Jan Swart for indicating the point of view of an intertwining relation. F.S. acknowledges NWO for financial support via the TOP1 grant 613.001.552.