Orthogonal Dualities of Markov Processes and Unitary Symmetries
We study self-duality for interacting particle systems, where the particles move as continuous-time random walkers having either exclusion interaction or inclusion interaction. We show that orthogonal self-dualities arise from unitary symmetries of the Markov generator. For these symmetries, we prov...
Збережено в:
| Опубліковано в: : | Symmetry, Integrability and Geometry: Methods and Applications |
|---|---|
| Дата: | 2019 |
| Автори: | , , , , |
| Формат: | Стаття |
| Мова: | Англійська |
| Опубліковано: |
Інститут математики НАН України
2019
|
| Онлайн доступ: | https://nasplib.isofts.kiev.ua/handle/123456789/210242 |
| Теги: |
Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Цитувати: | Orthogonal Dualities of Markov Processes and Unitary Symmetries / G. Carinci, C. Franceschini, C. Giardinà, W. Groenevelt, F. Redig // Symmetry, Integrability and Geometry: Methods and Applications. — 2019. — Т. 15. — Бібліогр.: 32 назв. — англ. |
Репозитарії
Digital Library of Periodicals of National Academy of Sciences of Ukraine| _version_ | 1860249570519285760 |
|---|---|
| author | Carinci, G. Franceschini, C. Giardinà, C. Groenevelt, W. Redig, F. |
| author_facet | Carinci, G. Franceschini, C. Giardinà, C. Groenevelt, W. Redig, F. |
| citation_txt | Orthogonal Dualities of Markov Processes and Unitary Symmetries / G. Carinci, C. Franceschini, C. Giardinà, W. Groenevelt, F. Redig // Symmetry, Integrability and Geometry: Methods and Applications. — 2019. — Т. 15. — Бібліогр.: 32 назв. — англ. |
| collection | DSpace DC |
| container_title | Symmetry, Integrability and Geometry: Methods and Applications |
| description | We study self-duality for interacting particle systems, where the particles move as continuous-time random walkers having either exclusion interaction or inclusion interaction. We show that orthogonal self-dualities arise from unitary symmetries of the Markov generator. For these symmetries, we provide two equivalent expressions that are related by the Baker-Campbell-Hausdorff formula. The first expression is the exponential of an anti-Hermitian operator and thus is unitary by inspection; the second expression is factorized into three terms and is proved to be unitary by using generating functions. The factorized form is also obtained by using an independent approach based on scalar products, which is a new method of independent interest that we introduce to derive (bi)orthogonal duality functions from non-orthogonal duality functions.
|
| first_indexed | 2025-12-07T21:24:52Z |
| format | Article |
| fulltext |
Symmetry, Integrability and Geometry: Methods and Applications SIGMA 15 (2019), 053, 27 pages
Orthogonal Dualities of Markov Processes
and Unitary Symmetries
Gioia CARINCI †, Chiara FRANCESCHINI ‡, Cristian GIARDINÀ §,
Wolter GROENEVELT † and Frank REDIG †
† Technische Universiteit Delft, DIAM, P.O. Box 5031, 2600 GA Delft, The Netherlands
E-mail: G.Carinci@tudelft.nl, W.G.M.Groenevelt@tudelft.nl, F.H.J.Redig@tudelft.nl
‡ Center for Mathematical Analysis Geometry and Dynamical Systems, Instituto
Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal
E-mail: Chiara.Franceschini@tecnico.ulisboa.pt
§ University of Modena and Reggio Emilia, FIM, via G. Campi 213/b, 41125 Modena, Italy
E-mail: Cristian.Giardina@unimore.it
Received December 24, 2018, in final form July 05, 2019; Published online July 12, 2019
https://doi.org/10.3842/SIGMA.2019.053
Abstract. We study self-duality for interacting particle systems, where the particles move
as continuous time random walkers having either exclusion interaction or inclusion interac-
tion. We show that orthogonal self-dualities arise from unitary symmetries of the Markov
generator. For these symmetries we provide two equivalent expressions that are related by
the Baker–Campbell–Hausdorff formula. The first expression is the exponential of an anti
Hermitian operator and thus is unitary by inspection; the second expression is factorized
into three terms and is proved to be unitary by using generating functions. The factorized
form is also obtained by using an independent approach based on scalar products, which
is a new method of independent interest that we introduce to derive (bi)orthogonal duality
functions from non-orthogonal duality functions.
Key words: stochastic duality; interacting particle systems; Lie algebras; orthogonal poly-
nomials
2010 Mathematics Subject Classification: 60J25; 82C22; 22E60
1 Introduction
In a series of previous works, dualities that are orthogonal in an appropriate Hilbert space have
been derived for a class of interacting particle systems with Lie-algebraic structure. This class
includes several well-known processes, for instance the generalized exclusion processes [23, 29],
the inclusion process [16], as well as independent random walkers [11]. These orthogonal dualities
were identified as classical orthogonal polynomials in [12] by using the structural properties of
those polynomials (recurrence relation and raising/lowering operators). In [26] the approach of
generating functions was used instead, by which non-polynomial orthogonal dualities (provided
by some other special functions, e.g., Bessel functions) were also found. Orthogonal duality
functions can also be explained using representation theory: they can be understood as the
intertwiner between two unitarily equivalent representations of a Lie algebra [13, 17].
Often the duality property of a Markov process can be related to the existence of some
(hidden) symmetries of the Markov generator, i.e., operators commuting with the generator of
the Markov process [14, 15]. This occurs for instance when the process has a reversible measure.
In this context detailed balance can be interpreted as a trivial duality, and by acting with
a symmetry of the generator one obtains a non-trivial duality. A natural question that arises is
thus what type of symmetries lead to orthogonal dualities. In this paper, for some particular
mailto:G.Carinci@tudelft.nl
mailto:W.G.M.Groenevelt@tudelft.nl
mailto:F.H.J.Redig@tudelft.nl
mailto:Chiara.Franceschini@tecnico.ulisboa.pt
mailto:Cristian.Giardina@unimore.it
https://doi.org/10.3842/SIGMA.2019.053
2 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
processes, we show that those symmetries have to be unitary and we single out the general
expression they must have.
We expect the association between orthogonal dualities and unitary symmetries to be robust
and apply in great generality to all cases where the duality function is obtained from the action
of a symmetry on the trivial duality. We choose here to focus on a class of processes having
an underlying Lie algebra structure that helps in the explicit characterization of the unitary
symmetry. We thus consider three interacting particle systems (exclusion, inclusion, independent
walkers) for which some orthogonal self-duality function are known and are given by classical
discrete orthogonal polynomials as the Meixner, Krawtchouk and Charlier polynomials. For
this class of processes we provide a full characterization of their unitary symmetries. This result
allows to identify the entire family of orthogonal duality functions, which turns out to be a two-
parameter family. For special values of the parameters we recover the orthogonal polynomial
duality. We expect similar results could be found for higher orthogonal polynomials that would
be associated to other Markov processes and their dualities.
The organization of this paper is as follows. In Section 2 we give an overview of the main
tools required to construct the setting. In Section 2.1 we recall the concept of (self-)duality
between Markov processes and we introduce the notion of equivalence between (self-)duality
functions. In Section 2.2 we introduce three algebras (su(2) algebra, su(1, 1) algebra and the
Heisenberg algebra) and the associated Markov processes that turn out to be interacting particle
systems. In Section 2.3 we recall from [15] a general scheme to construct duality functions for
Markov processes whose generator has an algebraic structure. In this approach there is a one-to-
one correspondence between self-duality functions and symmetries of the Markov generator. In
Section 3.1, by using this connection between duality functions and symmetries we present the
first main result of this paper. Namely, in Theorem 3.1 we provide the expression for the most
general unitary symmetry that will then yield orthogonal duality functions. We also identify the
special values of the parameters appearing in these symmetries for which the duality functions
are orthogonal polynomials. The proof of Theorem 3.1 is contained in Section 3.2. In Section 3.3
we provide a second expression for these unitary symmetries: it is a factorized expression for
function of the algebra generators that we show to be connected to the previous expression via the
Baker, Campbell, Hausdorff formula. In Section 4 we introduce a novel independent procedure
to obtain orthogonal duality functions. This new method relies on the use of a scalar product
in a Hilbert space. In Section 4.1 we prove that the scalar product of two duality functions
is again a new duality function and in Section 4.2 we show that these new duality functions
are biorthogonal by construction. We apply this technique in Section 4.3: for the interacting
particle systems considered in this paper by manipulation of the biorthogonal relation we get
an orthogonal relation.
The literature on stochastic duality for Markov processes is extremely vast. For the reader
convenience we recall [3, 10, 18, 20, 25, 31] for some applications to non-equilibrium statistical
physics, [5, 24] for duality in population models and [2, 8, 9] for the study of singular stochastic
PDE via duality. We also mention the algebraic approach to duality that shows that several
Markov processes dualities in turn derive from algebraic structures, see for instance [6, 7, 15,
22, 28].
The orthogonal dualities that were alluded to at the beginning of this introduction have been
introduced more recently in the literature. One might wonder what the advantages are of having
orthogonality. Assume the process has state space S and invariant measure µ and consider the
process generator as an operator on L2(S , µ). The duality function can be viewed as a family
of functions in the configurations of the original process (labelled by the configurations of the
dual process). Then, when this family happens to be linearly independent and complete so that
it gives a basis, it is natural to ask if/when this can be turned into a family of duality functions
that are orthogonal, thus yielding an orthogonal basis. It is not clear a priori that the natural
Orthogonal Dualities of Markov Processes and Unitary Symmetries 3
orthogonalization Gram–Schmidt procedure conserves the duality property. Thus this has to be
checked independently. In all cases, having an orthogonal basis will be helpful in studying the
contraction properties of the Markov semigroup, and thus quantifying for instance the rate of
relaxation to the invariant measure. Furthermore, in [1] orthogonal duality has been used to
prove a Boltzmann–Gibbs principle where several simplifications occur as a consequence of the
fact that the duality functions constitute an orthogonal basis for the Hilbert space.
2 Preliminaries
We start by recalling the definition of stochastic duality for two processes and introducing
the algebras and the interacting particle system (IPS) of interest. Our goal is to describe
a constructive technique, in which self-duality functions arise from both the symmetric approach
of Section 2.3 as well as from the inner product approach described in Section 4.
2.1 Stochastic duality
The definition of stochastic duality can be formalized for Markov processes as well as their
infinitesimal generators. Although they are not equivalent in general, they become equivalent
under suitable hypothesis regarding the semigroup associated to the generator of the process
discussed in Proposition 1.2 of [19].
Definition 2.1 (Markov duality definitions). Let X = (Xt)t≥0 and Y = (Yt)t≥0 be two con-
tinuous time Markov processes with state spaces S and S dual and generators L and Ldual
respectively. We say that Y is dual to X with duality function D : S ×S dual 7−→ R if
Ex[D(Xt, y)] = Ey[D(x, Yt)],
for all (x, y) ∈ S × S dual and t ≥ 0. If X and Y are two independent copies of the same
process, we say that Y is self-dual with self-duality function D. Duality can also be regarded
at the level of the processes generators. We say that Ldual is dual to L with duality function
D : S ×S dual 7−→ R if
[LD(·, y)](x) =
[
LdualD(x, ·)
]
(y).
If L = Ldual we have self-duality.
Note that self-duality can always be thought as a special case of duality where the dual process
is an independent copy of the first one. The simplification of self-duality for IPS typically arises
from the fact that the computation of correlation functions of the original process reduces to
studying a finite number of variables in the copy process.
Countable state space. If the original process (Xt)t≥0 and the dual process (Yt)t≥0 are
Markov processes with countable state space S and S dual resp., then the duality relation is
equivalent to∑
x′∈S
L(x, x′)D(x′, y) =
∑
y′∈S
Ldual(y, y′)D(x, y′) =
∑
y′∈S
(
Ldual
)T
(y′, y)D(x, y′), (2.1)
where LT denotes the transposition of the generator L. Generators are treated like (eventually
infinite) matrices and in matrix notation the identity (2.1) becomes
LD = D
(
Ldual
)T
. (2.2)
4 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
If Ldual = L we obtain the corresponding identities for self-duality. In this context, the genera-
tor L is given by a matrix known as rate matrix such that
L(x, y) ≥ 0 for x 6= y and
∑
y
L(x, y) = 0.
We say that the process jumps from x to y with rate L(x, y).
Definition 2.2 (duality functions in product form and single site duality functions.). The
duality functions we will present turn out to be of the following product structure
D(x, y) =
∏
i
d(xi, yi)
for (x, y) ∈ S ×S dual. The function inside the product will be regarded as “single site” duality
function and the subscript i removed.
Throughout the paper we will work with duality functions of this structure and so we will
only consider the single site.
Lemma 2.3 (notion of equivalence for duality functions.). If D(x, y) is a duality function
between two processes and the function c : S ×S dual −→ R is constant under the dynamics of
the two processes then Dc(x, y) = c(x, y)D(x, y) is also a duality function. We will refer to D
and Dc as equivalent duality functions.
For example, in the context of the processes we are interested in, we will see that the dynamics
conserves the total number of particles and dual particles, i.e.,
∑
i xi =
∑
i yi is conserved. As
a consequence of this we can always choose a self-duality function up to a multiplicative factor
in terms of the total number of particles. For example, if
D(x, y) =
n∏
i=1
d(xi, yi)
is a self-duality function, then for constants c and b, the function
Db,c(x, y) =
n∏
i=1
bxicyid(xi, yi)
is again a self-duality. This can easily be checked using Definition 2.1. Indeed,
Ex (Db,c(X(t), y)) = Ex
(
n∏
i=1
bXi(t)cyid(Xi(t), yi)
)
= b
∑
i xic
∑
i yiEx (D(X(t), y))
= b
∑
i xic
∑
i yiEx (D(x, Y (t))) = Ey
(
n∏
i=1
bxicYi(t)d(xi, Yi(t))
)
= Ey (Db,c(x, Y (t))) .
Our examples are all such that b = 1 and so we omit it.
2.2 Algebras and IPS
In the next three sections we introduce three algebras with three IPS, each one corresponding
to one of the three algebras. In particular, the probability measure that define the ∗-structure
of the algebra turns out to be the reversible measure of the particle process associated to that
algebra. Here we denote by F (S ) the space of real-valued functions on S , with countable S .
Orthogonal Dualities of Markov Processes and Unitary Symmetries 5
2.2.1 The Lie algebra su(1, 1) and symmetric inclusion process, SIP(k)
Generators of the dual Lie algebra su(1, 1) are K0, K+ and K−. They satisfy[
K0,K±
]
= ∓K± and
[
K+,K−
]
= 2K0.
We shall work in a representation, labeled by k ∈ R+, where the actions of the three generators
on functions f in F (N) is given by
(K+f)(x) := (2k + x)f(x+ 1),
(K−f)(x) := xf(x− 1), (2.3)
(K0f)(x) := (x+ k)f(x),
with f(−1) = 0. We define an inner product on F (N) by
〈f, g〉wp,k =
∑
x
f(x)g(x)wp,k(x), wp,k(x) =
Γ(2k + x)
x!Γ(2k)
px(1− p)2k, (2.4)
where 0 < p < 1, then su(1, 1) acts on the corresponding Hilbert space L2(wp,k) by unbounded
operators with dense domain the set of finitely supported functions on N. The adjoints of the
generators with respect the inner product are given by(
K0
)∗
= K0, (K+)∗ =
1
p
K−, (K−)∗ = pK+. (2.5)
The Casimir element is
Ω = 2
(
K0
)2 −K+K− −K−K+,
which is self-adjoint and commutes with every element of the Lie algebra.
The process associated with this algebra is the symmetric inclusion process SIP(2k), described
below. The inclusion process is introduced first in [14], and also studied further in [15]. The
SIP(2k) is a family of interacting particles processes labeled by the parameter k > 0 and that
can be defined on a generic graph G(V,E). The state space is unbounded so that each site can
have an arbitrary number of particles. The SIP(2k) generator is
LSIP(2k) =
∑
1≤i<l≤|V |
(i,l)∈E
L
SIP(2k)
i,l , (2.6)
L
SIP(2k)
i,l f(x) = xi(2k + xl)
[
f
(
xi,l
)
− f(x)
]
+ xl(2k + xi)
[
f
(
xl,i
)
− f(x)
]
,
where xi,l denotes the particle configuration obtained from the configuration x by moving one
particle from site i to site l, i.e., xi,l = x − δi + δl and so the dynamic conserves the total
number of particles. The generator can be defined on a weighted graph, however for the sake of
simplicity we restrict here to (2.6), since the duality functions will not depend on the weights of
the graph edges.
Clearly, the action of the generator involves only two connected sites and it can be produced
with the representation (2.3) acting on tensor products of F (N) via the expression of the co-
product of the Casimir Ω. Recall that the coproduct is an algebra homomorphism denoted by ∆
and defined by
∆(X) = X1 +X2,
6 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
for a Lie algebra element X. Here the subscript i indicates that the operator acts in the ith
factor of the tensor product. Higher-order coproducts are defined on Lie algebra elements X by
∆n(X) =
n+1∑
i=1
Xi,
which we consider as an operator on F (N)⊗(n+1). One can verify that for the couple of sites
(i, l) the generator of the SIP(2k) on two sites is written in terms of generators of the su(1, 1)
Lie algebra as
L
SIP(2k)
i,l = K+
i K
−
l +K−i K
+
l − 2K0
iK
0
l + 2k2 = −∆(Ω)i,l + 2k2.
This is an operator on F (N)⊗|V |, and the subscript i, l indicates that ∆(Ω) acts on the ith
and lth factor of the tensor product. (Note that ∆(Ω)i,l 6= Ωi + Ωl). Since Ω commutes
with every X ∈ su(1, 1) it follows that L
SIP(2k)
i,l commutes with ∆(X)i,l = Xi + Xl, and hence
LSIP(2k) =
∑
L
SIP(2k)
i,l commutes with ∆|V |−1(X) =
∑
Xi. Last, the reversible measure of the
SIP(2k) process is given by the homogeneous product measure with marginals the Negative
Binomial distributions with parameters 2k > 0 and 0 < p < 1, i.e., with probability mass
function wp,k of equation (2.4).
2.2.2 The Lie algebra su(2) and symmetric exclusion process, SEP(2j)
Generators of the dual su(2) Lie algebra are J0, J+ and J− which satisfy the following commu-
tation relations[
J0, J±
]
= ∓J± and [J+, J−] = −2J0.
We work in a representation of su(2) labeled by j ∈ N/2 on functions f in F ({0, 1, . . . , 2j})
given by
(J+f)(x) := (2j − x)f(x+ 1),
(J−f)(x) := xf(x− 1), (2.7)(
J0f
)
(x) := (x− j)f(x),
where f(−1) = f(2j + 1) = 0. Defining an inner product on F ({0, 1, . . . , 2j}) by
〈f, g〉wp,j =
∑
x
f(x)g(x)wp,j(x), wp,j(x) =
(
2j
x
)(
p
1− p
)x
(1− p)2j , (2.8)
where 0 < p < 1, the generators J0, J+ and J− act on the corresponding Hilbert space L2(wp,j),
with adjoints given by(
J0
)∗
= J0, (J+)∗ =
1− p
p
J−, (J−)∗ =
p
1− p
J+.
The Casimir element is
Ω = 2
(
J0
)2
+ J+J− + J−J+,
which is self-adjoint and commutes with every element in the Lie algebra.
The process associated with this algebra is the exclusion process, defined below. For j = 1/2
the boundary driven simple exclusion process has been studied using duality in [31]. The model
Orthogonal Dualities of Markov Processes and Unitary Symmetries 7
for arbitrary j has been introduced and studied in [29]. The SEP(2j) is a family of interacting
particles processes labeled by the parameter j ∈ N/2 and that can be defined on the same graph
G(V,E), as before. Each site (vertex) of G can have at most 2j particles and the SEP(2j)
generator is
LSEP(2j) =
∑
1≤i<l≤|V |
(i,l)∈E
L
SEP(2j)
i,l ,
L
SEP(2j)
i,l f(x) = xi(2j − xl)
[
f
(
xi,l
)
− f(x)
]
+ (2j − xi)xl
[
f
(
xl,i
)
− f(x)
]
.
As before we can write the generator of the SEP(2j) in two sites using the generators of the
su(2) algebra
L
SEP(2j)
i,l = J+
i J
−
l + J−i J
+
l + 2J0
i J
0
l − 2j2 = ∆(Ω)i,l − 2j2.
Last, the reversible measure of the SEP(2j) process is given by the homogeneous product measure
with marginals the Binomial distribution with parameters 2j > 0 and 0 < p < 1, i.e., with
probability mass function wp,j of equation (2.8).
2.2.3 The Heisenberg algebra and independent random walkers (IRW)
The dual Heisenberg algebra is the Lie algebra with generators a, a† and 1 such that[
a, a†
]
= −1.
The Heisenberg algebra has a representation on F (N) such that(
a†f
)
(x) = f(x+ 1),
(af)(x) = xf(x− 1) (2.9)
and 1 acts as the identity, and where f(−1) = 0. Consider the inner product
〈f, g〉wp =
∑
x
f(x)g(x)wp(x), wp(x) =
px
x!
e−p, (2.10)
where p > 0, then the Heisenberg algebra acts on the corresponding Hilbert space L2(wp) by
unbounded operators with dense domain the set of finitely support functions. The adjoints of a
and a† with respect to the inner product are
a∗ = pa† and
(
a†
)∗
=
1
p
a.
No such element as the Casimir is available for the Heisenberg algebra. The process associated
with this algebra is the process of independent random walkers (IRW), which was first intro-
duced in [30] and is well-known. They are defined in the usual setting, the process consists of
independent particles that perform a symmetric continuous time random walk at rate 1 on the
graph G(V,E). The generator is given by
LIRW =
∑
1≤i<l≤|V |
(i,l)∈E
LIRW
i,l ,
LIRW
i,l f(x) = xi
[
f
(
xi,l
)
− f(x)
]
+ xl
[
f
(
xl,i
)
− f(x)
]
.
8 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
In terms of the generators of the Heisenberg algebra we have
LIRW
i,l = a†ial + aia
†
l − aia
†
l − aia
†
l .
One can verify that LIRW
i,l commutes with ∆(X)i,l for every X in the Heisenberg algebra, so that
LIRW commutes with ∆|V |−1(X). The reversible invariant measure is provided by a homogeneous
product of Poisson distributions with parameter p > 0, i.e., with probability mass function wp
of equation (2.10).
2.3 Self-dualities via symmetries: general approach
and classical self-dualities
A general scheme for constructing self-dualities of continuous time Markov processes whose ge-
nerator has a symmetry, i.e., an operator commuting with its generator, has been first proposed
in [15]. In this section we first recall this approach and then we illustrate it by showing the
symmetry that is associated to classical self-duality functions. By construction, we are guaran-
teed that the functions we find via symmetries are self-dualities, but not orthogonal. However,
orthogonality can be inferred by proving that the symmetry is unitary. A family of unitary
symmetries will be found in Section 3 and, by specializing to some values of the parameters, we
will recover orthogonal dualities in terms of discrete orthogonal polynomials previously found
in [12, 26]. This orthogonality task is also addressed in Section 4 where we show that biortho-
gonality can be achieved by construction.
Recall that, since our processes are defined on a countable state space S , we can work with
the notion of duality in matrix notation, namely equation (2.2).
Definition 2.4. Let A and B be two matrices having the same dimension. We say that A is
a symmetry of B if A commutes with B, i.e.,
[A,B] = AB −BA = 0.
The main idea is that self-duality (in the context of Markov process with countable state
space) can be recovered starting from a trivial duality which is based on the reversible measures
of the processes. Then the action of a symmetry of the model on this trivial self-duality give
rise into a non-trivial one. The following results, whose proof can be found in [15], formalize
this idea.
Theorem 2.5 (symmetries and self-duality). Let d be a self-duality function of the generator L
and let S be a symmetry of L, then D = Sd is again a self-duality function for L.
If there is a description of the process generator in terms of a Lie algebra, then symmetries
can be constructed using this algebraic structure. The two main elements of Theorem 2.5 are
the initial self-duality d and the symmetry operator S. In general, if the process has a reversible
measure the self-duality d can easily be found starting from the reversibility.
Lemma 2.6 (diagonal self-duality and reversibility). If the process associated to generator L
has reversible measure µ, then the diagonal self-duality functions are of the form
d(x, y) =
δx,y
µ(x)
, where x, y ∈ S .
We refer to these diagonal self-duality functions as trivial or “cheap” self-duality functions.
The next lemma summarizes the cheap self-dualities for our three processes: notice that, up to
neglectable factors, they are the inverse of their reversible measure.
Orthogonal Dualities of Markov Processes and Unitary Symmetries 9
Lemma 2.7 (trivial self-duality functions). The processes of interests are self-dual with single
site diagonal self-duality function given by
Dch
p (x, y) =
y!Γ(2k)
Γ(2k + y)
p−yδx,y for the SIP(2k),
(2j − y)!y!
2j!
(
1−p
p
)y
δx,y for the SEP(2j),
y!
py
δx,y for the IRW.
We can now find several self-duality results applying the recipe of Theorem 2.5 using the
trivial self-duality function as starting point. We will give symmetries in terms of the expo-
nential function of Lie algebra elements considered as operators on function spaces, see (2.3),
(2.7) and (2.9). Recall that our representations are defined in terms raising/lowering operators
(i.e., shift by ±1) and diagonal operators. These operators are defined on finitely supported
functions, so that nth powers of the raising and lowering operators will always be 0 for large n.
Consequently, when considering exponential functions of raising and lowering operators acting
on finitely supported functions we do not have to worry about convergence of series, but in other
cases we have to check convergence. The exponential function of a diagonal operator is again
a diagonal operator. Moreover, it will be enough to provide symmetries acting on functions on
one site. Indeed, we have shown that the generator L of the process, which acts on N sites,
commutes with ∆N−1(Xn) for any n ∈ N and any Lie algebra element X. As a consequence L
also commutes with exp
(
∆N−1(X)
)
= exp(X1) exp(X2) · · · exp(XN ).
The following lemma shows how to find the so-called classical self-duality functions which
have a lower triangular structure.
Proposition 2.8 (classical self-duality functions and associated symmetries). The following
results hold.
1. The SIP(2k) is self-dual with single site self-duality function given by
Dcl
p (x, y) := S
(
Dch
p (·, y)
)
(x) =
x!
(x− y)!
Γ(2k)
Γ(2k + y)
p−y1{y≤x}, (2.11)
where S = eK
−
.
2. The SEP(2j) is self-dual with single site self-duality function given by
Dcl
p (x, y) := S
(
Dch
p (·, y)
)
(x) =
x!
(x− y)!
(2j − y)!
2j!
(
1− p
p
)y
1{y≤x},
where S = eJ
−
.
3. The IRW is self-dual with single site self-duality function given by
Dcl
p (x, y) := S
(
Dch
p (·, y)
)
(x) =
x!
(x− y)!
1
py
1{y≤x},
where S = ea.
Proof. We only consider the first item, the proof for the other two is similar. The fact that
Dcl
p (x, y) is a self-duality function is an immediate consequence of Theorem 2.5 since eK
−
com-
mutes with the Casimir Ω. The second equality in (2.11) follows from a straightforward calcu-
lation. Indeed, acting with the symmetry S, we have
Dcl(x, y) = eK
−(
Dch
p (·, y)
)
(x) =
∞∑
i=0
(K−)i
i!
y!Γ(2k)
Γ(2k + y)
(
1
p
)y
δx,y
10 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
=
∞∑
i=0
y!
i!
Γ(2k)
Γ(2k + y)
x!
(x− i)!
(
1
p
)y
1{i≤x}δx−i,y
=
x!
(x− y)!
Γ(2k)
Γ(2k + y)
(
1
p
)y
1{y≤x}. �
By virtue of Lemma 2.3 one can either neglect constants and factors that are constant under
the dynamic of the process or, on the other hand, add convenient choice of these constant factors.
In particular, in Section 4, we will fix the value of these constants in a suitable way.
3 Orthogonal self-dualities and unitary symmetries
In what follows, we will relate the orthogonal polynomials with their hypergeometric functions.
In general, the hypergeometric functions rFs is defined as an infinite series
rFs
(
a1, . . . , ar
b1, . . . , bs
;x
)
=
∞∑
k=0
(a1)k · · · (ar)k
(b1)k · · · (bs)k
xk
k!
,
where (a)k denotes the Pochhammer symbol defined in terms of the Gamma function as
(a)k :=
Γ(a+ k)
Γ(a)
.
Whenever one of the numerator parameters is a negative integer, the hypergeometric func-
tion rFs turns into a finite sum, so it is a polynomial in the other numerator parameters. We
define polynomials as in [21], in particular the following three discrete polynomials: Meixner
polynomials
M(x, y; p) = 2F1
(
−x,−y
2k
; 1− 1
p
)
for x, y ∈ N,
Krawtchouk polynomials
K(x, y; p) = 2F1
(
−x,−y
−2j
;
1
p
)
for x, y = 0, 1, . . . , 2j,
and the Charlier polynomials
C(x, y; p) = 2F0
(
−x,−y
−
;−1
p
)
for x, y ∈ N.
3.1 Main result
In this section we explicitly determine the symmetries S, given in terms of the underlying
Lie algebra generators, which allow to retrieve the orthogonal polynomials. It is important to
mention that, since we start from a (trivial) self-duality which is orthogonal with respect to the
measure w, the operator S that produces the orthogonal self-duality must be unitary. Recall
that a unitary operator in L2(S , w) is a linear operator such that
UU∗ = U∗U = I,
where U∗ is the adjoint of U in L2(S , w). As a consequence of this, we will have that U
preserves the inner product of the Hilbert space L2(S , w) and so the norm of the cheap self-
duality function Dch must be the same of the norm of the orthogonal self-duality function
Dor = SDch∥∥Dor
∥∥2
w
= 〈SDch, SDch〉w =
∥∥Dch
∥∥2
w
.
Orthogonal Dualities of Markov Processes and Unitary Symmetries 11
In the spirit of Proposition 2.8 we list the new orthogonal symmetries for the interacting particles
systems.
Theorem 3.1 (orthogonal self-duality functions and associated symmetries). The following
results hold.
1. For the SIP(2k) we have that
i) The symmetry
Sα,β = exp
(
β
(
−K+ +
1
p
K−
))
exp
(
iαK0
)
(3.1)
extends to a unitary operator for every choice of α, β ∈ R. As a consequence the func-
tions Sα,β
(
Dch
p (x, ·)
)
(y) are orthogonal (single site) self-duality functions in L2(wp,k)
with squared norm
∥∥Dch
p
∥∥2
wp,k
.
ii) Choosing α = α̂ = π and β = β̂ =
√
p arctanh
(√
p
)
we get the Meixner polynomials
up to a constant: Dor
p (x, y) := Sα̂,β̂
(
Dch
p (x, ·)
)
(y) = (p− 1)kM(x, y; p).
2. For the SEP(2j) we have that
i) The symmetry
Sα,β = exp
(
β
(
−J+ +
1− p
p
J−
))
exp
(
iαJ0
)
(3.2)
is unitary for every choice of α, β ∈ R. As a consequence the functions
Sα,β
(
Dch
p (x, ·)
)
(y) are orthogonal (single site) self-duality functions in L2(wp,j) with
squared norm
∥∥Dch
p
∥∥2
wp,j
.
ii) Choosing α = α̂ = π and β = β̂ =
√
p
1−p arctan
(√
p
1−p
)
we get the Krawtchouk
polynomials up to a constant: Dor
p (x, y) := Sα̂,β̂
(
Dch
p (x, ·)
)
(y) = (p− 1)jK(x, y; p).
3. For the IRW we have that
i) The symmetry
Sα,β = exp
(
β
(
− pa† + a
))
exp
(
iαaa†
)
(3.3)
extends to a unitary operator for every choice of α, β ∈ R. As a consequence the func-
tions Sα,β
(
Dch
p (x, ·)
)
(y) are orthogonal (single site) self-duality functions in L2(wp)
with squared norm
∥∥Dch
p
∥∥2
wp
.
ii) Choosing α = α̂ = π and β = β̂ = 1 we get the Charlier polynomials up to a constant:
Dor
p (x, y) := Sα̂,β̂
(
Dch
p (x, ·)
)
(y) = e−
p
2C(x, y; p).
3.2 Proof of the main result
We need the following lemma to introduce the generating function and to compute the action of
the algebra generators in order to prove Theorem 3.1. In particular we only consider the su(1, 1)
algebra and the SIP(2k) process, but for the other two processes the idea is the same.
Definition 3.2 (generating functions). We will always use the definition of generating function
as in [21, formula (9.10.11)], i.e., the generating function G of g(y) is defined as
(Gg)(t) :=
∞∑
y=0
g(y)
Γ(2k + y)
y!Γ(2k)
ty, t ∈ R.
12 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
The generating function of Meixner polynomials M(x, y; p) (see [21]) is
∞∑
y=0
M(x, y; p)
Γ(2k + y)
y!Γ(2k)
ty =
(
1− t
p
)x
(1− t)−2k−x. (3.4)
Lemma 3.3 (intertwining of the su(1, 1) algebra generators). The following results hold
1. GK−g(y) =
(
2kt+ t2 ∂∂t
)
Gg(t) =: K −Gg(t).
2. GK+g(y) =
(
∂
∂t
)
Gg(t) =: K +Gg(t).
3. GK0g(y) =
(
k + t ∂∂t
)
Gg(t) =: K 0Gg(t).
Note that K −, K + and K 0 so defined satisfy the commutation relations of the dual su(1, 1) Lie
algebra.
Proof. We have
GK−g(y) =
∞∑
n=0
yg(y − 1)
Γ(2k + y)
y!Γ(2k)
ty
= 2kt
∞∑
y=0
g(y)
Γ(2k + y)
y!Γ(2k)
ty + t2
∞∑
y=0
g(y)
Γ(2k + y)
y!Γ(2k)
yty−1
=
(
2kt+ t2
∂
∂t
)
Gg(t) = K −Gg(t).
This implicitly defines the operator K − which acts on functions of the t variable as
K − := 2kt+ t2
∂
∂t
.
Similarly,
GK+g(y) =
∞∑
y=0
(2k + y)g(y + 1)
Γ(2k + y)
y!Γ(2k)
ty
=
∞∑
y=0
g(y)
Γ(2k + y)
y!Γ(2k)
yty−1 =
(
∂
∂t
)
Gg(t) = K +Gg(t),
so the operator K + is a first derivative with respect to t, defined as
K +f(t) :=
∂f
∂t
(t).
For K0 we proceed in the same way
GK0g(y) =
∞∑
y=0
(k + y)g(y)
Γ(2k + y)
y!Γ(2k)
ty =
(
k + t
∂
∂t
)
Gg(t) = K 0Gg(t),
and we infer that
K 0f(t) :=
(
k + t
∂
∂t
)
f(t).
Note that for all the above we have called f(t) = (Gg(·))(t). �
Orthogonal Dualities of Markov Processes and Unitary Symmetries 13
Proof of Theorem 3.1. We will only give a proof for the first item as the other two follow
a similar strategy. The first point of the first item regards the unitarity of Sα,β in L2(wp,k),
which is achieved if (Sα,β)∗ = (Sα,β)−1. Using the adjoints (2.5) of K0, K+ and K− we have
that (Sα,β)∗ = exp
(
−iαK0
)
exp
(
β
(
−1
pK
− + K+
))
= (Sα,β)−1 as an operator acting on the
space of finitely supported functions. Since this is dense in L2(wp,k), Sα,β extends to a unitarity
operator. Unitary operators conserve the norm and so the norm of Sα,βD
ch
p (x, y) is the same as
the norm of Dch
p (x, y) in L2(wp,k). In particular, the two squared norms are
∥∥Sα,βDch
p
∥∥2
wp,k
=
∥∥Dch
p
∥∥2
wp,k
=
y!Γ(2k)
Γ(2k + y)
p−y(1− p)2k.
We show now the proof of the second point using a generating function approach. The idea is
to show that the generating function of Dor
p = Sα̂,β̂D
ch
p and the Meixner polynomials are the
same, i.e.,
G
(
Sα̂,β̂
(
Dch
p (x, ·)
))
(y) = G
(
(p− 1)kM(x, ·; p)
)
(y) (3.5)
and so using the generating function of Meixner polynomials in equation (3.4) one has that the
r.h.s. of equation (3.5) is (p − 1)−k(1 − t)−2k−x(1 − t
p
)x
. For the l.h.s. instead of computing
G
(
Sα,βD
ch
p
)
(t) we use Lemma 3.3 to evaluate Sα,β
(
GDch
p
)
(t), here
Sα,β = exp
(
β
(
−K + +
1
p
K −
))
exp
(
iαK 0
)
,
where K +, K − and K 0 are those in Lemma 3.3, which we consider as operators on functions
that are analytic at 0. In other words, we have to find the action of the operator Sα,β on
(
GDch
p
)
(t) =
∞∑
y=0
y!Γ(2k)
Γ(2k + y)
p−yδx,y
Γ(2k + y)
y!Γ(2k)
ty =
(
t
p
)x
.
The action of exp
(
iαK 0
)
on f(t) = Gg(t) is
(
eiα
)K 0
f(t) := G
((
eiα
)K0
g
)
(t) =
∞∑
y=0
Γ(2k + y)
y!Γ(2k)
ty(eiα)y+kg(y) = (eiα)kf
(
eiαt
)
.
Letting α = α̂ = π one has
exp
(
iπK 0
)
f(t) = (−1)kf(−t). (3.6)
To find the action of exp
(
β
(
−K + + 1
pK −
))
we will solve a partial differential equation, whose
solution ψ(t, β) is the action of Sα,β on function f(t). Using Lemma 3.3, this is
ψ(t, β) = e
β
[(
t2
p
−1
)
∂
∂t
+ 2k
p
t
]
f(t) (3.7)
with initial condition ψ(t, 0) = f(t). Here, it is understood that for the operator B :=
(
t2
p −
1
)
∂
∂t + 2k
p t and a function f in its domain, the exponential eβBf is defined as the solution of
the partial differential equation ∂
∂β g(β, t) = Bg(β, t) with initial condition g(0, t) = f(t). Thus,
deriving both sides of (3.7) with respect to β, we get a first-order PDE for ψ:
∂ψ
∂β
−
(
t2
p
− 1
)
∂ψ
∂t
− 2kt
p
ψ = 0. (3.8)
14 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
To solve the PDE we use the method of characteristics: we consider ψ along the characteristic
plane (τ, s), so that along a characteristic curve τ is constant and ψ(t, β) = ψ(t(s), β(s)). We
then have
∂ψ
∂s
=
∂ψ
∂β
∂β
∂s
+
∂ψ
∂t
∂t
∂s
.
Comparing the above with the PDE in equation (3.8) we just have to solve a system of three
first-order ODEs:
∂β
∂s
= 1,
∂t
∂s
=
p− t2
p
,
∂ψ
∂s
=
2kt
p
ψ.
From the first equation we have immediately that β = s, while the second has solution
t(s) =
√
p
tanh
(
s/
√
p
)
+ tanh(c1)
1 + tanh
(
s/
√
p
)
tanh(c1)
.
Using the initial condition t(0) =
√
p tanh(c1) = τ we get c1 = arctanh
(
τ/
√
p
)
and so
t(s) =
√
p
τ/
√
p+ tanh
(
s/
√
p
)
1 + τ/
√
p tanh
(
s/
√
p
) .
Substituting t in the last ODE we find that
ψ(s) =
(
τ sinh
(
s/
√
p
)
+
√
p cosh
(
s/
√
p
))2k
c2.
To find c2 we use the initial condition in the characteristic plane, i.e., ψ(0) = f(τ) = pkc2 so
c2 = f(τ)
pk
and so our solution in the (τ, s) plane is
ψ(τ, s) = f(τ)
(
τ
√
p
sinh
(
s/
√
p
)
+ cosh
(
s/
√
p
))2k
.
In the (t, β) plane this becomes
ψ(t, β) = f
(
√
p
t−√p tanh
(
β/
√
p
)
√
p− t tanh
(
β/
√
p
))(− t
√
p
sinh
(
β/
√
p
)
+ cosh
(
β/
√
p
))−2k
.
Setting β = β̂ = arctanh
(√
p
)√
p the above expression simplifies to
e
β̂
[(
−1+ t2
p
)
∂
∂t
+ 2k
p
t
]
f(t) =
(
1− t√
1− p
)−2k
f
(
t− p
1− t
)
. (3.9)
Equation (3.6) together with (3.9) finally gives
Sα̂,β̂f(t) = (p− 1)k(1− t)−2kf
(
p− t
1− t
)
. (3.10)
Last, we need to set f(t) =
(
GDch
p
)
(t) =
(
t
p
)x
to finally get
Sα̂,β̂
(
GDch
p
)
(t) = (p− 1)k(1− t)−2k−x
(
1− t
p
)x
,
which matches the generating function of the Meixner polynomials. �
In the following section we give a different expression for the three unitary symmetries Sα̂,β̂
of Theorem 3.1.
Orthogonal Dualities of Markov Processes and Unitary Symmetries 15
3.3 Factorized symmetries
We now want to study the unitary symmetries that arise from the previous section. Since we do
not know how to act with these symmetries on functions f(x) ∈ F (N), we wonder if a ‘factorized’
version of Sα̂,β̂ exists, i.e., if we can find a, b and c such that
Sα̂,β̂ = eaK
−
ebK
0
ecK
+
.
The advantage of having a factorized symmetry is that one can directly compute its action
on f(x) (without passing via generating functions), even if, on the other hand, the unitary
property is not an immediate consequence of this form. In the next section we will relate this
factorized form to another symmetry.
Theorem 3.4 (factorized unitary symmetries). The three orthogonal symmetries Sα̂,β̂ can also
be written in a factorized version using the appropriate algebra generators.
1. The action of Sα̂,β̂ in equation (3.1) coincides with the action of eK
−
elog(p−1)K0
epK
+
.
2. The action of Sα̂,β̂ in equation (3.2) coincides with the action of eJ
−
e
log
(
1
p−1
)
J0
e
p
1−pJ
+
.
3. The action of Sβ̂ in equation (3.3) coincides with the action of eae−p/2+iπaa†epa
†
.
Proof. We only show the first item as the other two have similar proofs; to do that we still use
generating functions. To show that
eK
−
(p− 1)K
0
epK
+
g(y) = exp
(
β̂
(
−K+ +
1
p
K−
))
exp
(
iα̂K0
)
g(y), (3.11)
we first consider the generating function G on both sides and then flip the action of G with the
one of the operators to get
eK−(p− 1)K 0
epK +
f(t) = exp
(
β̂
(
−K + +
1
p
K −
))
exp
(
iα̂K 0
)
f(t), (3.12)
where we called f(t) = (Gg)(t) and K −, K 0 and K + are those in Lemma 3.3. The r.h.s. of
equation (3.12) has been evaluated in the proof of Theorem 3.1, equation (3.10) so we just need
to find the action of eK− , (p− 1)K 0
and epK +
. Clearly,
(p− 1)K 0
f(t) = (p− 1)kf(t(p− 1))
since
(p− 1)K 0
f(t) := G
(
(p− 1)K
0
g(y)
)
=
∞∑
y=0
Γ(2k + y)
y!Γ(2k)
ty(p− 1)y+kg(y) = (p− 1)kf(t(p− 1)).
For eK− one can solve the associated PDE as in the proof of Theorem 3.1, or equivalently
considering the limit as p→ 0 on both sides of equation (3.9) and using that lim
p→0
arctanh(
√
p)√
p = 1
leads to
eK−f(t) = (1− t)−2kf
(
t
1− t
)
.
16 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
Last, for epK +
we have that(
epK +
f
)
(t) = ep
∂
∂t f(t) = f(t+ p)
since the action of the first derivative is a shift. Acting on f(t), we have
eK−(p− 1)K 0
epK +
f(t) = eK−(p− 1)K 0
f(t+ p) = (p− 1)keK−f(t(p− 1) + p)
= (p− 1)k(1− t)−2kf
(
t
1− t
(p− 1) + p
)
= (p− 1)k(1− t)−2kf
(
p− t
1− t
)
,
which matches the action of Sα̂,β̂ in equation (3.10). �
Remark 3.5 (Baker–Campbell–Hausdorff formula for dual su(1, 1) algebra). The identity given
in equation (3.11) can also be established as a consequence of the Baker–Campbell–Hausdorff
formula for the su(1, 1) algebra, see [32, formula (24b)] adapted to the dual su(1, 1) algebra. In
formula (24b) one has the following replacement L+ = 1√
pK
−, L− =
√
pK+ and L0 = K0 and
in particular one has to set τ = arctanh(
√
p) and α = π.
The added value of having the factorized version of the symmetry Sα̂,β̂ is that one can
immediately verify its action on the cheap duality Dch
p (x, y): via a straightforward computation
one can produce the orthogonal polynomials of Theorem 3.1, as we show in the proposition
below.
Proposition 3.6 (direct computation of orthogonal polynomials). Acting with the factorized
symmetry on the cheap self-duality function one gets the orthogonal self-duality function. In
particular, for the SIP(2k) this is
eK
−
elog(p−1)K0
epK
+(
Dch
p (x, ·)
)
(y) = Dor
p (x, y).
Proof. The proof follows a straightforward computation, see the appendix. �
4 Orthogonal self-duality via scalar products
In this section we first show how duality and self-duality function emerge as a consequence of
what we call scalar product approach and which is introduced below. We then give some hypothe-
sis to guarantee that such self-duality functions are biorthogonal. To conclude we implement this
recently developed technique to find Meixner polynomials as orthogonal self-duality functions
for the SIP(2k), in a similar way one could find orthogonal self-dualities for SEP(2j) and IRW.
4.1 Scalar product approach
In this section we present a new technique to approach duality: the naive idea is that the scalar
product of two duality functions is still a duality function. We define the scalar product on some
measure space L2(S , µ), in the usual way, i.e.,
〈f, g〉µ =
∑
x∈S
f(x)g(x)µ(x).
We will show that – in the setting of reversible processes – once two duality relations are
available then it is possible to generate new different duality functions starting from the initial
Orthogonal Dualities of Markov Processes and Unitary Symmetries 17
ones. Suppose we have three processes with generators L1, L2 and L3 and state space S1, S2
and S3, respectively. In particular, assume that d1 is a duality function for L1 and L2, while d2
is a duality function for L3 and L2, i.e.,
L1d1(·, y)(x) = L2d1(x, ·)(y) for (x, y) ∈ S1 ×S2 (4.1)
and
L3d2(·, y)(x) = L2d2(x, ·)(y) for (x, y) ∈ S3 ×S2. (4.2)
Then the following proposition holds.
Proposition 4.1 (new duality functions). If µ is a reversible measure for the generator L2 and
if equations (4.1) and (4.2) hold, then the function D : S1 ×S3 → R, given by
D(x, y) = 〈d1(x, ·), d2(y, ·)〉µ
is a duality function for L1 and L3. If L1 = L2 = L3 = L, then D is a new self-duality function
for L.
Proof. For i = 1, 2, 3, Li,xD(x, y) stands for (LiD(·, y))(x) the action of Li on the x variable
of D. Then,
L1,xD(x, y) = 〈L1,xd1(x, ·), d2(y, ·)〉µ =
∑
z
L2,zd1(x, z)d2(y, z)µ(z)
=
∑
z
d1(x, z)L2,zd2(y, z)µ(z) = 〈d1(x, ·), L3,yd2(y, ·)〉µ = L3,yD(x, y),
where we use duality of d1 (resp. d2) in the second (resp. fourth) equality and the self-adjointness
of L2 with respect to µ. �
A first application of the above proposition is shown in the example below, where we re-
cover Laguerre polynomials as duality function between SIP(2k) and BEP(2k), which we do not
introduce here, but it is well explained in [4, Section 2.2].
Example 4.2 (duality via scalar product). A parametrized family of reversible measure for the
SIP(2k) process is
µp(z) =
Γ(2k + z)
Γ(2k)z!
pz, z ∈ N, p ∈ (0, 1)
and the classical self-duality function Dcl
p for SIP(2k) is in equation (2.11), it will be our d1(x, y).
The last ingredient we need is a duality function between BEP(2k) and SIP(2k), a well known
in the literature, see [4, equation (4.9)], is
d2(x, y) =
xyΓ(2k)
Γ(2k + y)
(−1)y, x ∈ R+, y ∈ N.
In particular d2 is the one we need to obtain Laguerre polynomials. Proposition 4.1 assures
us that D(x, y) = 〈d2(x, ·), d2(y, ·)〉µp is a duality function between SIP(2k) and BEP(2k) and
a straightforward computation shows that D is the closed form of the Laguerre polynomials.
Indeed,
D(x, y) = 〈d2(x, ·), d1(y, ·)〉µp =
∞∑
z=0
(−x)zΓ(2k)
Γ(2k + z)
y!
(y − z)!
Γ(2k)
Γ(2k + z)
p−z
Γ(2k + z)
Γ(2k)z!
pz
=
y∑
z=0
(−x)z
z!
y!
(y − z)!
Γ(2k)
Γ(2k + z)
= 1F1
(
−y
2k
;x
)
for y ∈ N and x ∈ R+.
18 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
We can apply Proposition 4.1 for the same generator, to construct the Meixner polynomials
as SIP(2k) self-duality functions.
Example 4.3 (self-duality via scalar product). As for the previous Example 4.2, let µp(z) be the
reversible measure for the SIP(2k) process. Consider now two classical self-duality functions d1
and d2 as in equation (2.11). In particular, we are free to choose them without the constant,
i.e.,
d1(x, y) = d2(x, y) =
x!
(x− y)!
Γ(2k)
Γ(2k + y)
p−y1y≤x.
A simple computation shows that their scalar product in L2(µp) is a Meixner polynomial. Indeed,
D(x, y) = 〈d2(x, ·), d1(y, ·)〉µp
=
∞∑
z=0
x!
(x− z)!
y!
(y − z)!
(
Γ(2k)
Γ(2k + z)
)2
p−2z1z≤x1z≤y ·
Γ(2k + z)
Γ(2k)z!
pz
=
x∧y∑
z=0
1
z!
x!
(x− z)!
y!
(y − z)!
Γ(2k)
Γ(2k + z)
p−z = M(x, y; 1− p) for x, y ∈ N.
The following proposition expands the result of Proposition 4.1 in the context of self-duality.
It turns out that when two self-duality functions, d and D, are in a relation via a scalar product
with a third function F , then, assuming d to be a basis for L2(S , µ), F must also be a self-duality
function.
Proposition 4.4 (basis and self-duality). Assume that {x 7→ d(x, n) |n ∈ S } is a basis of self-
duality functions for L2(S , µ) where µ is a reversible measure for the generator L. Let F =
F (n, z) be a function on S ×S and define D by
D(x, n) := 〈d(x, ·), F (n, ·)〉µ.
If D is self-duality function, so is F .
Proof. Using the short notation we have that
LxD(x, n) = 〈Lxd(x, ·), F (n, ·)〉µ =
∑
z∈S
d(x, z)LzF (n, z)µ(z),
where we used that d is self-duality and that L is self-adjoint with respect to µ. On the other
hand, since D is a self-duality the above quantity must be equal to
LnD(x, n) = 〈d(x, ·), LnF (n, ·)〉µ =
∑
z∈S
d(x, z)LnF (n, x)µ(z).
From the identity LxD(x, n) = LnD(x, n), we have∑
z∈S
d(x, z) [LzF (n, z)− LnF (n, z)]µ(z) = 0
and since d is a basis for L2(S , µ), necessarily LzF (n, z)− LnF (n, z) = 0, i.e., F is also a self-
duality function for L. �
Orthogonal Dualities of Markov Processes and Unitary Symmetries 19
4.2 Biorthogonal self-dualities
How does the orthogonality property play a role? Not all self-duality functions built with this
method turn out to be orthogonal. However, there is a sort of stability with respect to this
orthogonal property in the scalar product construction. More precisely, if we start with two
biorthogonal self-duality functions the scalar product construction yields novel biorthogonal
self-duality functions that may happen to be equal and therefore orthogonal.
To state the next proposition, we will use that the inverse of the reversible measure is a self-
duality function as shown in Lemma 2.6.
Proposition 4.5 (biorthogonal self-duality functions). Let µ1 and µ2 be two reversible measures
and d1, d2 be two self-duality functions for the Markov process with generator L. Suppose that
〈d1(x, ·), d2(·, n)〉µ1 =
δx,n
µ2(n)
and 〈d2(x, ·), d1(·, n)〉µ2 =
δx,n
µ1(n)
. (4.3)
Then the functions
D(x, n) := 〈d1(x, ·), d1(n, ·)〉µ1 , D̃(x, n) := 〈d2(·, x), d2(·, n)〉µ1
are biorthogonal self-duality functions for L, i.e.,〈
D(·,m), D̃(·, n)
〉
µ2
=
δm,n
µ2(m)
.
In particular, if D = D̃ we have the orthogonality relations for D.
Proof. From Proposition 4.1 we have that both D and D̃ are self-duality functions since scalar
product of self-dualities. Assuming now we can interchange the order of summation:〈
D(·,m), D̃(·, n)
〉
µ2
=
∑
x
D(x,m)D̃(x, n)µ2(x)
=
∑
x
(∑
y
d1(x, y)d1(m, y)µ1(y)
)(∑
z
d2(z, x)d2(z, n)µ1(z)
)
µ2(x)
=
∑
y,z
d1(m, y)d2(z, n)µ1(y)µ1(z)
∑
x
d2(z, x)d1(x, y)µ2(x)
=
∑
y,z
µ1(y)µ1(z)d1(m, y)d2(z, n)
δy,z
µ1(y)
=
∑
y
d1(m, y)d2(y, n)µ1(y) =
δm,n
µ2(m)
. �
We now implement this method to get the result below. Here we find Meixner polynomi-
als as biorthogonal self-duality functions and with the aid of some hypergeometric functions
transformation we find an orthogonal duality function.
4.3 From biorthogonal to orthogonality self-duality functions
According to Proposition 4.5 we need two duality functions d1 and d2 satisfying (4.3). For
SIP(2k) recall that
µp(z) :=
Γ(2k + z)
Γ(2k)z!
pz
20 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
is the marginal of the (product) reversible (non normalized) measure, and
Dcl
1
λ
(x, n) :=
x!Γ(2k)λy
(x− y)!Γ(2k + y)
1y≤x
are the (single-site) classical self-duality functions. In the following we denote by 〈·, ·〉p the scalar
product with respect to the non-normalized measure µp. We have the following lemma which
we show for SIP(2k).
Lemma 4.6 (input self-duality functions). For any p, q ∈ R we have〈
Dcl
−p(x, ·), Dcl
−q(·, n)
〉
p
=
δx,n
µq(x)
,
where Dcl are the classical self-duality functions introduced in Proposition 2.8.
Proof. Note that Dcl
1
λ
(x, y) = 0 if y > x. Then
〈
Dcl
1
λ
(x, ·), Dcl
1
α
(·, n)
〉
p
=
x∑
z=n
x!Γ(2k)λz
(x− z)!Γ(2k + z)
z!Γ(2k)αn
(z − n)!Γ(2k + n)
pyΓ(2k + z)
Γ(2k)z!
,
which equals 0 if x < n. Suppose x ≥ n, then shifting the summation index (m = z − n) gives
〈
Dcl
1
λ
(x, ·), Dcl
1
α
(·, n)
〉
p
=
(pλα)nΓ(2k)
Γ(2k + n)
x−n∑
m=0
x!
(x− n−m)!m!
(λp)m.
Now use (A)k+l = (A)k(A + k)l, i.e., Γ(A+k+l)
Γ(A) = Γ(A+k)
Γ(A)
Γ(A+k+l)
Γ(A+k) and the binomial theorem to
obtain
〈
Dcl
1
λ
(x, ·), Dcl
1
α
(·, n)
〉
p
=
x!
(x− n)!
Γ(2k)
Γ(2k + n)
(pλα)n
x−n∑
m=0
(x− n)!
(x− n−m)!m!
(λp)m
=
x!
(x− n)!
Γ(2k)
Γ(2k + n)
(pλα)n(1 + λp)x−n.
Setting λ = −1
p and α = −1
q we get the result, i.e.,
〈
Dcl
−p(x, ·), Dcl
−q(·, n)
〉
p
=
x!Γ(2k)
Γ(2k + x)
(
+
1
q
)x
δx,n =
δx,n
µq(x)
. �
Proposition 4.7 (from biorthogonal to orthogonality self-duality functions). The self-duality
functions
D(x, n) =
〈
Dcl
−q(x, ·), Dcl
−q(n, ·)
〉
q
and
D̃(x, n) =
〈
Dcl
−p(·, x), Dcl
−p(·, n)
〉
q
are biorthogonal with respect to the stationary measure of their associated process. In details,
1. For SIP(2k) we have
D(x, n) = 2F1
(
−x,−n
2k
;
1
q
)
Orthogonal Dualities of Markov Processes and Unitary Symmetries 21
and
D̃(x, n) =
(
q
p(1− q)
)n+x
(1− q)−2k
2F1
(
−x,−n
2k
;
1
q
)
and they are biorthogonal with respect to the measure wp,k. In particular, for the choice
1
q
= 1− 1
p
we have
D(x, n) = M(x, n; p) and D̃(x, n) = (1− p)2kD(x, n)
so that the biorthogonality relation recovers the orthogonality relation of Meixner polyno-
mials.
2. For SEP(2j) we have that
D(x, n) = 2F1
(
−x,−n
−2j
;
q − 1
q
)
and
D̃(x, n) =
(
q(p− 1)
p
)n+x
(1− q)−2j
2F1
(
−x,−n
−2j
;
q − 1
q
)
and they are biorthogonal with respect to the measure wp,j. In particular, for the choice
1
q
= 1− 1
p
we have
D(x, n) = K(x, n; p) and D̃(x, n) = (1− p)2jD(x, n)
so that the biorthogonality relation recovers the orthogonality relation of Krawtchouk poly-
nomials.
3. For IRW we have that
D(x, n) = 2F0
(
−x,−n
−
;
1
q
)
and D̃(x, n) =
(
−q
p
)x+n
eq 2F0
(
−x,−n
−
;
1
q
)
and they are biorthogonal with respect to the measure wp. In particular, for the choice
q = −p
we have
D(x, n) = C(x, n; p) and D̃(x, n) = e−pD(x, n)
so that the biorthogonality relation recovers the orthogonality relation of Charlier polyno-
mials.
22 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
Proof. As always we show the proof for the first item only. Now we apply Proposition 4.5 with
d1(x, n) = Dcl
−q(x, n) and d2(x, n) = Dcl
−p(x, n)
then, from Lemma (4.6), the conditions (4.3) are satisfied for µ1 = µq and µ2 = µp. We have
D(x, n) =
〈
Dcl
−q(x, ·), Dcl
−q(n, ·)
〉
q
= 2F1
(
−x,−n
2k
;
1
q
)
= M
(
x, n;
q
q − 1
)
and
D̃(x, n) =
〈
Dcl
−p(·, x), Dcl
−p(·, n)
〉
q
= (−p)−x−n
∞∑
z=n
(2k)zq
z
z!
z!
(z − x)!(2k)x
z!
(z − n)!(2k)n
=
(−p)−n−xqn
(2k)x(2k)n
∞∑
m=0
(2k)m+n(m+ n)!
m!(m+ n− x)!
qm
=
(−p)−n−xqnn!
(2k)x(n− x)!
∞∑
m=0
(2k + n)m(n+ 1)m
m!(n− x+ 1)m
qm
=
(−p)−n−xqnn!
(2k)x(n− x)!
2F1
(
2k + n, n+ 1
n− x+ 1
; q
)
.
By applying the following 2F1-transformations [21, equations (2.2.6), (2.3.14), (2.2.6)], we obtain
2F1
(
2k + n, n+ 1
n− x+ 1
; q
)
= (1− q)−n−2k
2F1
(
2k + n,−x
n− x+ 1
;
q
q − 1
)
= (1− q)−n−2k (2k)x
(−n)x
2F1
(
−x, 2k + n
2k
;
1
1− q
)
= (1− q)−n−x−2k(−q)x (2k)x
(−n)x
2F1
(
−x,−n
2k
;
1
q
)
.
This gives
D̃(x, n) =
(
− q
p(1− q)
)n+x
(1− q)−2k
2F1
(
−x,−n
2k
;
1
q
)
then, for 1
p = 1− 1
q we get
D̃(x, n) = (1− q)−2k
2F1
(
−x,−n
2k
;
1
q
)
= (1− p)2k
2F1
(
−x,−n
2k
; 1− 1
p
)
= (1− p)2kM(x, n; p). �
Remark 4.8. From the previous proposition we have that we can write the Meixner-duality
function D(x, n) = M(x, n; p) in two forms in terms of scalar product:
D(x, n) =
〈
Dcl
p
1−p
(x, ·), Dcl
p
1−p
(n, ·)
〉
p
p−1
= (1− p)−2k
〈
Dcl
−p(·, x), Dcl
−p(·, n)
〉
p
p−1
. (4.4)
We now write Dcl
p
1−p
(x, y) and Dcl
−p(y, x) as a symmetry acting on the cheap duality. This
allows us to write both expressions for D(x, n) in (4.4) via two symmetries (S1 and S2) acting
on the cheap duality. Before doing that, we will use the following lemma to justify some equality
in the computation below.
Orthogonal Dualities of Markov Processes and Unitary Symmetries 23
Lemma 4.9 (duality of K+ and K− via the cheap duality function). The operators K+ and K−
are dual via
Dch
1 (x, y) =
x!Γ(2k)
Γ(2k + x)
δx,y.
Proof. This follows from(
K+Dch
1 (·, y)
)
(x) =
y!Γ(2k)
Γ(2k + y − 1)
=
(
K−Dch
1 (x, ·)
)
(y). �
As a consequence of the above relation, we have the following corollary.
Corollary 4.10. The operators eK
+
and eK
−
are also in duality via Dch
1 . Moreover, we can
choose parameter α, β on the exponentials and λ on Dch
1 such that the relation(
eαK
−
Dch
1
λ
(·, y)
)
(x) =
(
eβK
+
Dch
1
λ
(x, ·)
)
(y)
is always true for any α, β and λ ∈ R that satisfy α = λβ.
To make notation simpler we write K−1 (resp. K−2 ) for the action of K− on the first (resp.
second) variable and same for K+. Let’s now investigate the two symmetries associated to the
self-duality function in equation (4.4).
Proposition 4.11 (two ways of expressing orthogonal polynomials). Let D be the self-duality
function given by the two scalar products in equation (4.4), then D can be written as a symmetry
acting on Dch. In particular, we have that
D(x, n) = eK
−
1 e
p
p−1
K+
1 Dch
p
p−1
(x, n) (4.5)
for the first scalar product in (4.4), and
D(x, n) = (1− p)−2ke−pK
+
1 e
1
1−pK
−
1 Dch
p(p−1)(x, n) (4.6)
for the second one.
Proof. For the first scalar product in (4.4) we have that, using the first item of Proposition 2.8
D(x, n) =
〈
Dcl
p
1−p
(x, ·), Dcl
p
1−p
(n, ·)
〉
p
p−1
=
〈
eK
−
1 Dch
p
1−p
(x, ·), eK
−
1 Dch
p
1−p
(n, ·)
〉
p
p−1
. (4.7)
The action of eK
−
1 only affect the x variable and it can be placed outside the scalar product,
moreover by Corollary 4.10, the first scalar product in (4.7) is equal to
eK
−
1
〈
Dch
p
1−p
(x, ·), e
p
1−pK
+
2 Dch
p
1−p
(n, ·)
〉
p
p−1
.
The adjoint of K+ with respect to µ p
p−1
is p−1
p K− and the above quantity becomes
eK
−
1
〈
e−K
−
2 Dch
p
1−p
(x, ·), Dch
p
1−p
(n, ·)
〉
p
p−1
.
Using again Corollary 4.10 we finally get
eK
−
1 e
p
p−1
K+
1
〈
Dch
p
1−p
(x, ·), Dch
p
1−p
(n, ·)
〉
p
p−1
= eK
−
1 e
p
p−1
K+
1 Dch
p
p−1
(x, n),
where we computed the scalar product to get the symmetry in (4.5).
24 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
In a similar fashion, for the second scalar product in (4.7) we have that
D(x, n) = (1− p)−2k
〈
Dcl
−p(·, x), Dcl
−p(·, n)
〉
p
p−1
= (1− p)−2k
〈
eK
−
2 Dch
−p(x, ·), eK
−
2 Dch
−p(n, ·)
〉
p
p−1
by Corollary 4.10. Then, by considering the adjoint of K− with respect to µ p
p−1
, the above
expression reads
(1− p)−2ke−pK
+
1
〈
e
p
p−1
K+
2 Dch
−p(x, ·), Dch
−p(n, ·)
〉
p
p−1
and again, by Corollary 4.10 we have the symmetry in (4.6)
(1− p)−2ke−pK
+
1 e
1
1−pK
−
1
〈
Dch
−p(x, ·), Dch
−p(n, ·)
〉
p
p−1
= (1− p)−2ke−pK
+
1 e
1
1−pK
−
1 Dch
p(p−1)(x, n). �
Remark 4.12 (a commutation relation for the exponentials). The expression for D(x, n) in (4.6)
can be written as an action on Dch
p(p−1)(x, n) as follows,
D(x, n) = (1− p)−2ke−pK
+
1 e
1
1−pK
−
1 Dch
p(p−1)(x, n) = e−pK
+
1 e
1
1−pK
−
1 (1− p)−2K0
1Dch
p
p−1
(x, n),
comparing with (4.5) allows us to infer the following relation
eK
−
e
p
p−1
K+
= e−pK
+
e
1
1−pK
−
(1− p)−2K0
. (4.8)
Relation in (4.8) is found in [27, Remark 3.2] adapted to the su(1, 1) Lie algebra generators.
We now do some consideration about the relations among the symmetries in equations (4.5)
and (4.6) with the one obtained in Proposition 3.4. In order to compare their action, we first
realize that
Dch
p
p−1
(x, n) = (p− 1)xDch
p (x, y),
so that their action on Dch
p reads
(p− 1)kD(x, n) = S1D
ch
p (x, y) := eK
−
e
p
p−1
K+
(p− 1)K
0
Dch
p (x, y)
and
(p− 1)kD(x, n) = S2D
ch
p (x, y) := e−pK
+
e
1
1−pK
−
(p− 1)−K
0
Dch
p (x, y).
It is easy to verify that both S1 and S2 are unitary operators on L2(wp,k) since the self-duality
functions (p − 1)kD(x, n) and Dch
p (x, y) have the same norm. One could also check it via the
generating function approach.
Remark 4.13 (equivalence of the symmetries). As a final remark we mention that the symmetry
from Proposition 3.4 and S1 are, as expected, the same. If we consider their action on Dch
p (x, y),
then (
eK
−
e
p
p−1
K+
(p− 1)K
0
Dch
p (·, y)
)
(x) = (p− 1)kM(x, y; p),
while (
eK
−
(p− 1)K
0
epK
+
Dch
p (·, y)
)
(x) = (p− 1)kM(x, y; p)
inferring the commutation relation between K+ and K0:
e
p
p−1
K+
(p− 1)K
0
= (p− 1)K
0
epK
+
.
Orthogonal Dualities of Markov Processes and Unitary Symmetries 25
A Appendix. Proof of Proposition 3.6
We have to show that
eK
−
elog(p−1)K0
epK
+(
Dch
p (x, ·)
)
(y) = Dor
p (x, y),
where
Dch
p (x, y) =
y!Γ(2k)
Γ(2k + y)
p−y and Dor
p (x, y) = (p− 1)k 2F1
(
−x,−y
2k
; 1− 1
p
)
.
We start by acting with the inner operator on the y variable of Dch
p :
eK
−
elog(p−1)K0
epK
+(
Dch
p (x, ·)
)
(y)
= eK
−
elog(p−1)K0
∞∑
i=0
pi
i!
(K+)i
x!Γ(2k)
Γ(2k + x)
(p)−xδx,y
= eK
−
elog(p−1)K0
∞∑
i=0
pi
i!
Γ(2k + y + i)
Γ(2k + y)
x!Γ(2k)
Γ(2k + x)
p.xδx,y+i
since the action of the ith power of K+ is (K+)if(y) = Γ(2k+y+i)
Γ(2k+y) f(y + i). The action of K0 is
diagonal and so we have
eK
−
elog(p−1)K0
px−y
x!
(x− y)!
Γ(2k)
Γ(2k + y)
p−x1{x≥y}
= eK
−
(p− 1)y+k x!
(x− y)!
Γ(2k)
Γ(2k + y)
(
1
p
)y
1{x≥y}.
Finally we use the action of the ith power of K−, i.e., (K−)if(y) = y!
(y−i)!f(y − i),
(p− 1)k
∞∑
i=0
(K−)i
i!
(
p− 1
p
)y x!
(x− y)!
Γ(2k)
Γ(2k + y)
1{x≥y}
= (p− 1)k
∞∑
i=0
1
i!
y!
(y − i)!
(
p− 1
p
)y−i x!
(x− y + i)!
Γ(2k)
Γ(2k + y − i)
1{x≥y−i}1{i≤y}
= (p− 1)k
y∑
i=0∨(y−x)
1
i!
y!
(y − i)!
(
p− 1
p
)y−i x!
(x− y + i)!
Γ(2k)
Γ(2k + y − i)
= (p− 1)k
x∧y∑
s=0
x!
(x− s)!
y!
(y − s)!
1
s!
(
p− 1
p
)s Γ(2k)
Γ(2k + s)
= (p− 1)kM(x, y; p) = Dor
p (x, y).
Where we performed the change of variable y − i = s at the end. Note that, up to the constant
(p−1)k, the last sum is the closed form of Meixner polynomials of parameter p and 2k, variable x
and degree y. Since the result is symmetric in x and y then the action on the x variable would
produce the same result.
Acknowledgements
We thank the anonymous referees for their input which helped to improve the paper.
26 G. Carinci, C. Franceschini, C. Giardiná, W. Groenevelt and F. Redig
References
[1] Ayala M., Carinci G., Redig F., Quantitative Boltzmann–Gibbs principles via orthogonal polynomial duality,
J. Stat. Phys. 171 (2018), 980–999, arXiv:1712.08492.
[2] Borodin A., Corwin I., Gorin V., Stochastic six-vertex model, Duke Math. J. 165 (2016), 563–624,
arXiv:1407.6729.
[3] Borodin A., Corwin I., Sasamoto T., From duality to determinants for q-TASEP and ASEP, Ann. Probab.
42 (2014), 2314–2382, arXiv:1207.5035.
[4] Carinci G., Giardinà C., Giberti C., Redig F., Duality for stochastic models of transport, J. Stat. Phys. 152
(2013), 657–697, arXiv:1212.3154.
[5] Carinci G., Giardinà C., Giberti C., Redig F., Dualities in population genetics: a fresh look with new
dualities, Stochastic Process. Appl. 125 (2015), 941–969, arXiv:1302.3206.
[6] Carinci G., Giardinà C., Redig F., Sasamoto T., Asymmetric stochastic transport models with Uq(su(1, 1))
symmetry, J. Stat. Phys. 163 (2016), 239–279, arXiv:1507.01478.
[7] Carinci G., Giardinà C., Redig F., Sasamoto T., A generalized asymmetric exclusion process with Uq(sl2)
stochastic duality, Probab. Theory Related Fields 166 (2016), 887–933, arXiv:1407.3367.
[8] Corwin I., Ghosal P., Shen H., Tsai L.-C., Stochastic PDE limit of the six vertex model, arXiv:1803.08120.
[9] Corwin I., Petrov L., Stochastic higher spin vertex models on the line, Comm. Math. Phys. 343 (2016),
651–700, arXiv:1502.07374.
[10] Corwin I., Shen H., Tsai L.-C., ASEP(q, j) converges to the KPZ equation, Ann. Inst. Henri Poincaré
Probab. Stat. 54 (2018), 995–1012, arXiv:1602.01908.
[11] De Masi A., Presutti E., Mathematical methods for hydrodynamic limits, Lecture Notes in Mathematics,
Vol. 1501, Springer-Verlag, Berlin, 1991.
[12] Franceschini C., Giardinà C., Stochastic duality and orthogonal polynomials, arXiv:1701.09115.
[13] Franceschini C., Giardinà C., Groenevelt W., Self-duality of Markov processes and intertwining functions,
Math. Phys. Anal. Geom. 21 (2018), Art. 29, 21 pages, arXiv:1801.09433.
[14] Giardinà C., Kurchan J., Redig F., Duality and exact correlations for a model of heat conduction, J. Math.
Phys. 48 (2007), 033301, 15 pages, arXiv:cond-mat/0612198.
[15] Giardinà C., Kurchan J., Redig F., Vafayi K., Duality and hidden symmetries in interacting particle systems,
J. Stat. Phys. 135 (2009), 25–55, arXiv:0810.1202.
[16] Giardinà C., Redig F., Vafayi K., Correlation inequalities for interacting particle systems with duality,
J. Stat. Phys. 141 (2010), 242–263, arXiv:0906.4664.
[17] Groenevelt W., Orthogonal stochastic duality functions from Lie algebra representations, J. Stat. Phys. 174
(2019), 97–119, arXiv:1709.05997.
[18] Imamura T., Sasamoto T., Current moments of 1D ASEP by duality, J. Stat. Phys. 142 (2011), 919–930,
arXiv:1011.4588.
[19] Jansen S., Kurt N., On the notion(s) of duality for Markov processes, Probab. Surv. 11 (2014), 59–120,
arXiv:1210.7193.
[20] Kipnis C., Marchioro C., Presutti E., Heat flow in an exactly solvable model, J. Stat. Phys. 27 (1982),
65–74.
[21] Koekoek R., Lesky P.A., Swarttouw R.F., Hypergeometric orthogonal polynomials and their q-analogues,
Springer Monographs in Mathematics, Springer-Verlag, Berlin, 2010.
[22] Kuan J., An algebraic construction of duality functions for the stochastic Uq
(
A
(1)
n
)
vertex model and its
degenerations, Comm. Math. Phys. 359 (2018), 121–187, arXiv:1701.04468.
[23] Liggett T.M., Interacting particle systems, Classics in Mathematics, Springer-Verlag, Berlin, 2005.
[24] Möhle M., The concept of duality and applications to Markov processes arising in neutral population genetics
models, Bernoulli 5 (1999), 761–777.
[25] Ohkubo J., On dualities for SSEP and ASEP with open boundary conditions, J. Phys. A: Math. Theor. 50
(2017), 095004, 21 pages, arXiv:1606.05447.
[26] Redig F., Sau F., Factorized duality, stationary product measures and generating functions, J. Stat. Phys.
172 (2018), 980–1008, arXiv:1702.07237.
https://doi.org/10.1007/s10955-018-2060-7
https://arxiv.org/abs/1712.08492
https://doi.org/10.1215/00127094-3166843
https://arxiv.org/abs/1407.6729
https://doi.org/10.1214/13-AOP868
https://arxiv.org/abs/1207.5035
https://doi.org/10.1007/s10955-013-0786-9
https://arxiv.org/abs/1212.3154
https://doi.org/10.1016/j.spa.2014.10.009
https://arxiv.org/abs/1302.3206
https://doi.org/10.1007/s10955-016-1473-4
https://arxiv.org/abs/1507.01478
https://doi.org/10.1007/s00440-015-0674-0
https://arxiv.org/abs/1407.3367
https://arxiv.org/abs/1803.08120
https://doi.org/10.1007/s00220-015-2479-5
https://arxiv.org/abs/1502.07374
https://doi.org/10.1214/17-AIHP829
https://doi.org/10.1214/17-AIHP829
https://arxiv.org/abs/1602.01908
https://doi.org/10.1007/BFb0086457
https://arxiv.org/abs/1701.09115
https://doi.org/10.1007/s11040-018-9289-x
https://arxiv.org/abs/1801.09433
https://doi.org/10.1063/1.2711373
https://doi.org/10.1063/1.2711373
https://arxiv.org/abs/cond-mat/0612198
https://doi.org/10.1007/s10955-009-9716-2
https://arxiv.org/abs/0810.1202
https://doi.org/10.1007/s10955-010-0055-0
https://arxiv.org/abs/0906.4664
https://doi.org/10.1007/s10955-018-2178-7
https://arxiv.org/abs/1709.05997
https://doi.org/10.1007/s10955-011-0149-3
https://arxiv.org/abs/1011.4588
https://doi.org/10.1214/12-PS206
https://arxiv.org/abs/1210.7193
https://doi.org/10.1007/BF01011740
https://doi.org/10.1007/978-3-642-05014-5
https://doi.org/10.1007/s00220-018-3108-x
https://arxiv.org/abs/1701.04468
https://doi.org/10.1007/b138374
https://doi.org/10.2307/3318443
https://doi.org/10.1088/1751-8121/aa56f8
https://arxiv.org/abs/1606.05447
https://doi.org/10.1007/s10955-018-2090-1
https://arxiv.org/abs/1702.07237
Orthogonal Dualities of Markov Processes and Unitary Symmetries 27
[27] Rosengren H., A new quantum algebraic interpretation of the Askey–Wilson polynomials, in q-Series from
a Contemporary Perspective (South Hadley, MA, 1998), Contemp. Math., Vol. 254, Amer. Math. Soc.,
Providence, RI, 2000, 371–394.
[28] Schütz G.M., Duality relations for asymmetric exclusion processes, J. Stat. Phys. 86 (1997), 1265–1287.
[29] Schütz G.M., Sandow S., Non-Abelian symmetries of stochastic processes: Derivation of correlation functions
for random-vertex models and disordered-interacting-particle systems, Phys. Rev. E 49 (1994), 2726–2741.
[30] Spitzer F., Interaction of Markov processes, Adv. Math. 5 (1970), 246–290.
[31] Spohn H., Long range correlations for stochastic lattice gases in a nonequilibrium steady state, J. Phys. A:
Math. Gen. 16 (1983), 4275–4291.
[32] Truax D.R., Baker–Campbell–Hausdorff relations and unitarity of SU(2) and SU(1, 1) squeeze operators,
Phys. Rev. D 31 (1985), 1988–1991.
https://doi.org/10.1090/conm/254/03962
https://doi.org/10.1007/BF02183623
https://doi.org/10.1103/PhysRevE.49.2726
https://doi.org/10.1016/0001-8708(70)90034-4
https://doi.org/10.1088/0305-4470/16/18/029
https://doi.org/10.1088/0305-4470/16/18/029
https://doi.org/10.1103/PhysRevD.31.1988
1 Introduction
2 Preliminaries
2.1 Stochastic duality
2.2 Algebras and IPS
2.2.1 The Lie algebra su(1,1) and symmetric inclusion process, SIP(k)
2.2.2 The Lie algebra su(2) and symmetric exclusion process, SEP(2j)
2.2.3 The Heisenberg algebra and independent random walkers (IRW)
2.3 Self-dualities via symmetries: general approach and classical self-dualities
3 Orthogonal self-dualities and unitary symmetries
3.1 Main result
3.2 Proof of the main result
3.3 Factorized symmetries
4 Orthogonal self-duality via scalar products
4.1 Scalar product approach
4.2 Biorthogonal self-dualities
4.3 From biorthogonal to orthogonality self-duality functions
A Appendix. Proof of Proposition 3.6
References
|
| id | nasplib_isofts_kiev_ua-123456789-210242 |
| institution | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| issn | 1815-0659 |
| language | English |
| last_indexed | 2025-12-07T21:24:52Z |
| publishDate | 2019 |
| publisher | Інститут математики НАН України |
| record_format | dspace |
| spelling | Carinci, G. Franceschini, C. Giardinà, C. Groenevelt, W. Redig, F. 2025-12-04T13:09:43Z 2019 Orthogonal Dualities of Markov Processes and Unitary Symmetries / G. Carinci, C. Franceschini, C. Giardinà, W. Groenevelt, F. Redig // Symmetry, Integrability and Geometry: Methods and Applications. — 2019. — Т. 15. — Бібліогр.: 32 назв. — англ. 1815-0659 2010 Mathematics Subject Classification: 60J25; 82C22; 22E60 arXiv: 1812.08553 https://nasplib.isofts.kiev.ua/handle/123456789/210242 https://doi.org/10.3842/SIGMA.2019.053 We study self-duality for interacting particle systems, where the particles move as continuous-time random walkers having either exclusion interaction or inclusion interaction. We show that orthogonal self-dualities arise from unitary symmetries of the Markov generator. For these symmetries, we provide two equivalent expressions that are related by the Baker-Campbell-Hausdorff formula. The first expression is the exponential of an anti-Hermitian operator and thus is unitary by inspection; the second expression is factorized into three terms and is proved to be unitary by using generating functions. The factorized form is also obtained by using an independent approach based on scalar products, which is a new method of independent interest that we introduce to derive (bi)orthogonal duality functions from non-orthogonal duality functions. We thank the anonymous referees for their input, which helped to improve the paper. en Інститут математики НАН України Symmetry, Integrability and Geometry: Methods and Applications Orthogonal Dualities of Markov Processes and Unitary Symmetries Article published earlier |
| spellingShingle | Orthogonal Dualities of Markov Processes and Unitary Symmetries Carinci, G. Franceschini, C. Giardinà, C. Groenevelt, W. Redig, F. |
| title | Orthogonal Dualities of Markov Processes and Unitary Symmetries |
| title_full | Orthogonal Dualities of Markov Processes and Unitary Symmetries |
| title_fullStr | Orthogonal Dualities of Markov Processes and Unitary Symmetries |
| title_full_unstemmed | Orthogonal Dualities of Markov Processes and Unitary Symmetries |
| title_short | Orthogonal Dualities of Markov Processes and Unitary Symmetries |
| title_sort | orthogonal dualities of markov processes and unitary symmetries |
| url | https://nasplib.isofts.kiev.ua/handle/123456789/210242 |
| work_keys_str_mv | AT carincig orthogonaldualitiesofmarkovprocessesandunitarysymmetries AT franceschinic orthogonaldualitiesofmarkovprocessesandunitarysymmetries AT giardinac orthogonaldualitiesofmarkovprocessesandunitarysymmetries AT groeneveltw orthogonaldualitiesofmarkovprocessesandunitarysymmetries AT redigf orthogonaldualitiesofmarkovprocessesandunitarysymmetries |