Polyharmonic functions for finite graphs and Markov chains

TL;DR

This paper introduces and analyzes polyharmonic functions on finite graphs and Markov chains, solving boundary value problems and comparing with infinite tree cases, expanding the understanding of discrete potential theory.

Contribution

It defines polyharmonic functions for finite graphs and Markov chains, and solves the Riquier boundary value problem in this discrete setting.

Findings

Unique solutions to boundary problems are established.

Polyharmonic functions are characterized in finite graph and Markov chain contexts.

Connections to infinite tree boundary theory are explored.

Abstract

On a finite graph with a chosen partition of the vertex set into interior and boundary vertices, a $\lambda$-polyharmonic function is a complex function $f$ on the vertex set which satisfies $(\lambda \cdot I - P)^n f(x) = 0$ at each interior vertex. Here, $P$ may be the normalised adjaceny matrix, but more generally, we consider the transition matrix $P$ of an arbitrary Markov chain to which the (oriented) graph structure is adapted. After describing these `global' polyharmonic functions, we turn to solving the Riquier problem, where $n$ boundary functions are preassigned and a corresponding `tower' of $n$ successive Dirichlet type problems are solved. The resulting unique solution will be polyharmonic only at those points which have distance at least $n$ from the boundary. Finally, we compare these results with those concerning infinite trees with the end boundary, as studied by Cohen, Colonnna, Gowrisankaran and Singman, and more recently, by Picardello and Woess.

Equations

$$Pf(x)=\sum _{y}p(x,y)f(y),$$

$$Ph(x)=\lambda h(x)\text{for every}x\in X^o.$$

$$(\lambda \cdot I-P{)}^{n}f=0\text{on }X.$$

$${\Delta }_{\lambda }=\left(\begin{array}{cc}\lambda \cdot I_{X^o}-P_{X^o}& -Q\\ 0& 0\end{array}\right)=\left(\begin{array}{cc}\lambda \cdot I_{X^o}& 0\\ 0& I_{\partial X}\end{array}\right)-P,$$

$${\Delta }_{\lambda }^{n}f=0$$

$$h(x)=\sum _{w\in \partial X}F(x,w)g(w),$$

$$\mathbf{G}(\lambda)=\bigl{(}G(x,y|\lambda)\bigr{)}_{x,y\in X^{o}}=(\lambda\cdot I_{X^{o}}-P_{X^{o}})^{-1}$$

$$\mathbf{F}(\lambda)=\bigl{(}F(x,w|\lambda)\bigr{)}_{x\in X^{o},w\in\partial X}=\mathbf{G}(\lambda)\,Q.$$

$$G(x,y|\lambda )=\sum _{n=0}^{\infty }{p}^{(n)}(x,y)/{\lambda }^{n+1}\text{and}F(x,w|\lambda )=\sum _{n=0}^{\infty }{f}^{(n)}(x,w)/{\lambda }^n,$$

$$\sum _{w\in \partial X}F(x,w)=1\text{for every }x\in X^o.$$

$${\chi }_P(\lambda )=\mathsf{det}(\lambda \cdot I-P)=(\lambda -1{)}^{|\partial X|}{\chi }_{P^o}(\lambda ).$$

$$\bigl{\{}f_{j}^{(k)}:k=1,\dots,\min\{n,\kappa_{j}\}\,,\;j=1,\dots,\mu\bigr{\}}.$$

$$h(x)=\sum _{w\in \partial X}F(x,w|\lambda )g(w),x\in X^o,$$

$$(\lambda \cdot I_{X^o}-P_{X^o})h^o=Qg,$$

$${\Delta }_{\lambda }^n=\left(\begin{array}{cc}(\lambda \cdot I_{X^o}-P_{X^o}{)}^n& -(\lambda \cdot I_{X^o}-P_{X^o}{)}^{n-1}Q\\ 0& 0\end{array}\right).$$

$$(\lambda \cdot I_{X^o}-P_{X^o}{)}^n{h}^o=(\lambda \cdot I_{X^o}-P_{X^o}{)}^{n-1}Qg,$$

$$f_{r}={\Bigl{(}{\displaystyle f_{r}^{o}\atop\displaystyle g_{r}}\Bigr{)}}\,,\quad\Delta_{\lambda}^{o}f_{n}=0\,,\quad\text{and}\quad\Delta_{\lambda}^{o}f_{r}=f_{r+1}^{o}\quad\text{for}\quad r=n\!-\!1,n\!-\!2,\dots,1\,.$$

$$f(x)=\sum_{r=1}^{n}\bigl{[}\mathbf{G}(\lambda)^{r}\,Q\,g_{r}\bigr{]}(x)\,,\;x\in X^{o}\,,$$

$$f_{2}(x)=\sum_{r=2}^{n}\bigl{[}\mathbf{G}(\lambda)^{r-1}\,Q\,g_{r}\bigr{]}(x)\,,\;x\in X^{o}\,,$$

$$\Delta_{\lambda}^{o}f=f_{2}^{o}\,,\quad\text{where}\quad f={\textstyle\Bigl{(}{\displaystyle f^{o}\atop\displaystyle g_{1}}\Bigr{)}}\,$$

$$(\lambda \cdot I_{X^o}-P_{X^o})f^o-Q{g}_1=f_2^o.$$

$${\partial }^{n}X=\{x\in X:p^{(k)}(x,w)>0\text{for some}w\in \partial X\text{and}k\le n-1\},$$

$$(\lambda \cdot I_{X^o}-P_{X^o}){\mathbf{G}}^{(r)}(\lambda )=-k\cdot {\mathbf{G}}^{(r-1)}(\lambda ),$$

$$\mathbf{G}(\lambda {)}^r=\frac{(-1{)}^{r-1}}{(r-1)!}{\mathbf{G}}^{(r-1)}(\lambda ).$$

$${\mathsf{res}}^{\ast }(P_{X^o})=\mathsf{res}(P_{X^o})\setminus \{\lambda :F(o,w|\lambda )=0\text{for some}w\in \partial X\}.$$

$$K^{(X)}(x,w|\lambda )=\frac{F(x,w|\lambda )}{F(o,w|\lambda )},x\in X,w\in \partial X.$$

$$h(x)$$

$$\nu (w)$$

$$f(x)=\sum _{r=1}^n{\int }_{\partial X}{K}_r^{(X)}(x,\cdot |\lambda )d{\nu }_r,\text{where for }w\in \partial X$$

$$K_r^{(X)}(\cdot ,w|\lambda )=\mathbf{G}(\lambda {)}^{r-1}K(\cdot ,w|\lambda )\text{and}{\nu }_r(w)=g_r(w)F(o,w|\lambda ).$$

$$f(x)$$

Full text

1. Introduction

In the setting of the classical Laplacian $\Delta$ on a Euclidean domain, or the Laplace-Beltrami operator on a Riemannian manifold, a polyharmonic function $f$ is one for which $\Delta^{n}f=0$. Their study goes back to work in the $19^{\text{th}}$ century, see e.g. Almansi [1]. A basic reference is the monograph by Aronszajn, Creese and Lipkin [3]. A more recent one is the volume by Gazzola, Grunau and Sweers [8], with a nice introduction to classical problems from elasticity where polyharmonic (in fact biharmonic) functions and $\Delta^{2}$ come up.

While there is a huge body of literature in the smooth case, the literature in the discrete setting is quite restricted: an early reference is Voronkova [14], who analysed the discretized version of $\Delta^{2}f=0$ in a half-strip $[0\,,\infty]\times[0\,,\,H]$. Other quite early references are Yamasaki [16] and Kayano and Yamasaki [10] who investigated the Green kernel for the bi-Laplacian on an infinite network, and a follow-up of this is Venkataraman [13]. Biharmonic Laplacians on trees where also studied by Cohen, Colonna and Singman [6], [7], seemingly without link to [16] and [10]. Prior to that, Cohen, Colonna, Gowrisankaran and Singman [5] were the first to undertake a detailed study of polyharmonic functions on infinite, locally finite trees. In particular, for the standard Laplacian arising from simple random walk on a regular tree, they provided a boundary integral representation which is an analogue of Almansi’s expansion of polyharmonic functions on the unit disk. (To get a flavour of the many close analogies between the potential theory of the unit disk and regular trees, the reader is invited to the introductory sections of Boiko and Woess [4].) Recently, Picardello and Woess [12] extended the study of [5] and proved, among outher, a boundary integral representation of $\lambda$-polyharmonic functions (see below for more details) for arbitrary nearest neighbour transition operators on countable trees, not necessarily required to be locally finite.

In all this work, finite graphs, resp. Markov chains had only marginal appearances: in [16] for the biharmonic Green function of finite subnetworks of an infinite network, and in [5] for finite trees and an associated boundary value problem for biharmonic functions. Anandam [2] also studies polyharmonic functions on finite subtrees of infinite trees.

In the present note, we elaborate a detailed account of the general finite case, in which the mentioned potential theoretic questions turn into issues of linear algebra which can be solved rather easily.

The setting. We start with a finite set $X$, subdivided into the disjoint union of two non-empty subsets $X^{o}$, the interior, and $\partial X$, the boundary. On $X$, we consider a stochastic transition matrix P=\bigl{(}p(x,y)\bigr{)}_{x,y\in X} with the following properties, where $p^{(n)}(x,y)$ denotes the $(x,y)$-entry of the matrix power $P^{n}$.

1. (i)

For all $x\in X^{o}$, there is $w\in\partial X$ such that $p^{(n)}(x,w)>0$ for some $n$. 2. (ii)

For all $w\in\partial X$, we have $p(w,w)=1$, and thus $p(w,x)=0$ for all $x\in X\setminus\{w\}$. 3. (iii)

For all $w\in\partial X$, there is $x\in X^{o}$ such that $p^{(n)}(x,w)>0$ for some $n$.

Thus, $X$ can be given the structure of a digraph, where we have an oriented edge $x\to y$ when $p(x,y)>0$. Then (i) means that the boundary can be reached from any interior point by an oriented path, (ii) means that each boundary point is absorbing, i.e., the only outgoing edge is a loop at that point, and (iii) means that every boundary point is active in the sense that it is reached by some oriented path from an interior point. In probabilistic terms, we have a Markov chain (random process) on $X$, whose evolution is governed by $P\,$: if the current position is $x$ then the next step is from $x$ to $y$ with probability $p(x,y)$.

(1.1) Example.

The most typical situation is the one where we start with a finite resisitive network, that is, a connected, non-oriented graph $(X,E)$ where each edge $e=[x,y]=[y,x]$ carries a positive conductance $a(e)=a(x,y)$. Then we choose our partition $X=X^{o}\cup\partial X$, and we set $m(x)=\sum_{y}a(x,y)$. The transition probabilities become $p(x,y)=a(x,y)/m(x)$, if $x\in X^{o}$ and $y\in X$, while $p(w,w)=1$ for $w\in\partial X$. This defines a reversible Markov chain which is absorbed in $\partial X$, see e.g. Woess [15, Ch. 4]. In particular, setting all $a(x,y)$ equal to $1$, the conductances correspond to the adjacency matrix.

The transition matrix $P$ acts on functions (column vectors) $f:X\to\mathbb{C}$ by

[TABLE]

and the (normalised) graph Laplacian is $I-P$, where $I=I_{X}$ is the identity matrix over $X$. It is typically defined on $X$ without assigning a boundary $\partial X$, but the study undertaken here makes sense only in presence of absorbing points. Note that the more direct analogue of the (negative definite) smooth Laplacian would in reality be $P-I$. More generally, we shall work with suitable variant of $\lambda\cdot I-P$ for $\lambda\in\mathbb{C}$.

A $\lambda$-harmonic function $h:X\to\mathbb{C}$ is one for which

[TABLE]

When $\lambda=1$, we speak of a harmonic function. When speaking of $\lambda$-polyharmonic functions of order $n$, we have two possible approaches: one is to look for functions $f:X\to\mathbb{C}$ which satisfy

[TABLE]

These global $\lambda$-polyharmonic functions can be easily described.

The more interesting version is related with the pre-assignment of boundary values. Let $P_{X^{o}}$ and $Q$ be the restrictions of $P$ to $X^{o}\times X^{o}$ and $X^{o}\times\partial X$, respectively. Then we define the $\lambda$-Laplacian as the matrix given in block-form by

[TABLE]

where the [math]s stand for the zero matrices in the respective dimensions. Here, the identity matrix over $\partial X$ is not multiplied by $\lambda$, so that functions annihilated by $\Delta_{\lambda}$ are $\lambda$-harmonic only in $X^{o}$.

Our main focus is on polyharmonic functions in the sense that they satisfy

[TABLE]

on $X$, or – more reasonably, as we shall see – on the “$n$-th interior” of $X$, i.e., all points in $X^{o}$ from which $\partial X$ cannot be reached in less then $n$ steps. When $\lambda=1$, the two notions (1.3) and (1.5) coincide.

This note is organized as follows. In §2, we first consider ordinary harmonic and polyharmonic functions, that is the case $\lambda=1$. After recalling the well known solution of the Dirichlet problem for harmonic functions with preassigned boundary values (Lemma 2.2), we explain why all global harmonic functions in the sense of (1.3) (with $\lambda=1$) are indeed harmonic (Proposition 2.6). Then we look at all global $\lambda$-polyharmonic functions as in (1.3). In this case, $\lambda$ must belong to the spectrum of $P_{X^{o}}$, and the solutions can be described in terms of a Jordan basis (Proposition 2.7).

In §3, we turn to studying $\Delta_{\lambda}$ and its powers, for $\lambda$ in the resolvent set of $P_{X^{o}}$ (the spectrum being settled in §2). There is a direct analogue to the solution of the Dirichlet problem, and again, any function which satisfies $\Delta_{\lambda}^{n}f=0$ on all of $X$ must be $\lambda$-harmonic (Proposition 3.2). Finally, we give the precise formulation of the Riquier problem, which consists in assigning boundary functions $g_{1}\,,...,g_{n}$ and – loosely spoken – searching for a function $f$ such that the boundary values of $\Delta_{\lambda}^{r-1}f$ coincide with $g_{r}$ for $r=1,\dots,n$. That problem for the special case of finite trees is briefly touched in [5]. Here, we provide the general solution (Theorem 3.4).

Finally, in §4, we undertake a comparison of those results with the case of infinite trees without leaves, which was studied recently in [12] by use of Martin boundary theory.

All results of this note are achieved by applying basic tools from Linear Algebra in the right way. We believe that this material provides a useful basis, firstly as a link to the classical, smooth case (regarding the Laplacian on bounded domains), and secondly, as a basis for handling and understanding polyharmonic functions not only on infinite trees, but also on more general infinite graphs and their boundaries at infinity.

2. The Dirichlet problem, and global $\lambda$-polyharmonic functions

We start with some observations on the case $\lambda=n=1$, that is, ordinary harmonic functions. We start with a simple observation on $\mathsf{spec}(P_{X^{o}})$, the set of eigenvalues of $P_{X^{o}}$.

(2.1) Lemma.

The spectral radius $\rho=\rho(P_{X^{o}})=\max\{|\lambda|:\lambda\in\mathsf{spec}(P_{X^{o}})\}$ satisfies $\rho<1$.

Proof (outline).

Condition (i) on $P$ implies that for each $x\in X^{o}$, there is $n$ such that $\sum_{v\in X^{o}}p^{(n)}(x,v)<1$, that is, $P_{X^{o}}^{n}$ is strictly substochastic in the row of $x$. One easily deduces that there is $m$ such that $P_{X^{o}}^{m}$ is strictly substochastic in every row, which yields the claim. ∎

The following solution of the Dirichlet problem is folklore in the Markov chain community; see e.g. [15, §6.A]. It keeps being “rediscovered” by analysts who deviate into the discrete world, see for example Kiselman [11].

(2.2) Lemma.

For every function $g:\partial X\to\mathbb{C}$ there is a unique harmonic function $h$ on $X$ such that $h|_{\partial X}=g$. It is given by

[TABLE]

where $F(x,w)$ is the probability that the Markov chain starting at $x$ hits $\partial X$ in the point $v$.

We next want to describe the kernel $F(x,w)$ in matrix terminology. Let $\mathsf{res}(P_{X^{o}})=\mathbb{C}\setminus\mathsf{spec}(P_{X^{o}})$ be the resolvent set of $P_{X^{o}}$. For $\lambda\in\mathsf{res}(P_{X^{o}})$, the resolvent is the $X^{o}\times X^{o}$-matrix

[TABLE]

The kernels $G(x,y|\lambda)$ are called Green functions. They are rational functions of $\lambda$. Now we define the $X^{o}\times\partial X$-matrix

[TABLE]

We can extend it to $X\times\partial X$ by setting $F(v,w|\lambda)=\delta_{w}(v)$ for $v,w\in\partial X$. When $\lambda=1$, we just write $G(x,y)$ for $G(x,y|1)$ and $F(x,w)$ for $F(x,w|1)$. For $|\lambda|>\rho$, we can expand

[TABLE]

where the probabilistic meaning is that for the Markov chain starting at $x$, the probability to be at $y$ at time $n$ is $p^{(n)}(x,y)$, while $f^{(n)}(x,w)$ is the probability that the first visit in $w\in\partial X$ occurs at time $n$.

Coming back to the Dirichlet problem, it is a straightforward matrix computation to see that the function $h$, as defined in Lemma 2.2, is harmonic. Its uniqueness follows from invertibility of $(I_{X^{o}}-P_{X^{o}})$. Instead, it may also be instructive to deduce uniqueness from the potential theoretic maximum principle: every real valued harmonic function attains its maximum on $\partial X$, see [15, §6.A].

This also yields one way to see that the Markov chain must hit the boundary almost surely, that is

[TABLE]

Namely, the unique harmonic extension of the constant boundary function $g\equiv 1$ is the constant function $h\equiv 1$ on $X$. Also, the function $x\mapsto F(x,w)$ provides the unique harmonic extension of the boundary function $g=\mathbf{1}_{v}\,$.

(2.5) Corollary.

The geometric and the algebraic multiplicity of the eigenvalue $\lambda=1$ of $P$ coincide and are equal to $|\partial X|$.

Proof.

Lemma 2.2 yields that the geometric multiplicity is $|\partial X|$. The characteristic polynomial of the matrix $P$ is

[TABLE]

By Lemma 2.1, $\chi_{P^{o}}(1)\neq 0$. ∎

Now we can easily describe all free polyharmonic functions of order $n\geq 1$, that is, those which satisfy $(I-P)^{n}f=0$ on $X$.

(2.6) Proposition.

A function $f:X\to\mathbb{C}$ satisfies $(I-P)^{n}f=0$ if and only if $f$ is harmonic.

Proof.

Suppose $n\geq 2$, and let $h=(I-P)^{n-1}f$. Then $h$ is harmonic, and $(I-P)f=h$. Since $(I-P)^{n-1}f=0$ on $\partial X$, the function $h$ solves the Dirichlet problem with boundary values [math]. Therefore $h=0$, that is, $(I-P)^{n-1}f=0$. Proceeding by induction, we obtain that $f$ is harmonic. ∎

Similarly, we can handle the case $(\lambda\cdot I-P)^{n}f=0$, when $\lambda\neq 1$. First of all, when $n\geq 2$ then the function $h=(\lambda\cdot I-P)^{n-1}f$ satisfies $Ph=\lambda\cdot h$. Second, we see that $f=0$ in $\partial X$, so (by abuse of notation) we consider $f$ as a function on $X^{o}$. In other words, $\lambda\in\mathsf{spec}(P_{X^{o}})$.

Let $\kappa=\kappa(\lambda)$ and $\mu=\mu(\lambda)$ be the algebraic and geometric eigenvalue multiplicities of $\lambda$. Let $h_{1}\,,\dots\,,h_{\mu}$ be a basis of $\mathsf{ker}(\lambda\cdot I_{X^{o}}-P_{X^{o}})$. For each $j\in\{1,\dots,\mu\}$, let $\kappa_{j}$ be the length of the associated Jordan chain (= dimension of the associated Jordan block in the Jordan normal form). That is, $\kappa_{1}+\dots+\kappa_{\mu}=\kappa$, and we have functions $f_{j}^{(k)}$, $k=1,\dots,\kappa_{j}$ such that $f_{j}^{(1)}=h_{j}$ and $(\lambda\cdot I_{X^{o}}-P_{X^{o}})f_{j}^{(k)}=f_{j}^{(k-1)}$ for $k\geq 2$. All those functions are extended to $X$ by assigning value [math] on $\partial X$. Then it is clear that $\;\{f_{j}^{(k)}:k=1,\dots,\kappa_{j}\,,\;j=1,\dots,\mu\}\;$ is a basis of the linear space of all global $\lambda$-polyharmonic functions (of arbitrary order). We subsume.

(2.7) Proposition.

With the above notation, for $\lambda\in\mathsf{spec}(P_{X^{o}})$, the space of functions $f:X\to\mathbb{C}$ with $(\lambda\cdot I-P)^{n}f=0$ is spanned by

[TABLE]

(2.8) Corollary.

For a finite network with boundary as in Example (1.1), every global $\lambda$-polyharmonic $h$ function satisfies $Ph=\lambda\cdot h$, and $\lambda\in\mathsf{spec}(P)\subset\mathbb{R}$. Furthermore, $h$ vanishes on $\partial X$ when $\lambda\neq 1$.

Proof.

If we define the diagonal matrix M=\mathsf{diag}\bigl{(}\sqrt{m(x)}\,\bigr{)}_{x\in X^{o}}, then $M\,P_{X^{o}}\,M^{-1}$ is symmetric, so that the spectrum is real and the geometric and algebraic multiplicities of the eigenvalues of $P_{X^{o}}$ coincide. ∎

3. Boundary value problems for $\lambda$-polyharmonic functions

In this section, we assume that $\lambda\in\mathsf{res}(P_{X^{o}})$ and study the operator (resp. matrix) $\Delta_{\lambda}$ of (1.4) and its powers.

Notation: in accordance with the block form used above, for any function $f:X\to\mathbb{C}$ we write f=\Bigl{(}{\displaystyle f^{o}\atop\displaystyle f^{\partial}}\Bigr{)}, where $f^{o}=f|_{X^{o}}$ and $f^{\partial}=f|_{\partial X}\,$. Also, we write $\Delta_{\lambda}^{o}$ for the restriction of the matrix of (1.4) to $X\times X^{o}$, that is, $\Delta_{\lambda}^{o}f=(\Delta_{\lambda}f)^{o}$.

First of all, there is an obvious $\lambda$-variant of the solution of the Dirichlet problem.

(3.1) Lemma.

Let $\lambda\in\mathsf{res}(P_{X_{o}})$. For every function $g:\partial X\to\mathbb{C}$ there is a unique $\lambda$-harmonic function $h$ on $X$ such that $h|_{\partial X}=g$. It is given by

[TABLE]

where $F(x,w|\lambda)$ is defined by (2.4).

Proof.

We write h=\Bigl{(}{\displaystyle h^{o}\atop\displaystyle g}\Bigr{)}, where $h^{o}=h|_{X^{o}}$ and $g$ is the given boundary function. Then the equation $\Delta_{\lambda}=0$ transforms into

[TABLE]

which has the unique solution $h^{o}=\mathbf{G}(\lambda)Qg\,$, as proposed. ∎

Next, we note that

[TABLE]

Thus, if wo look for a solution of $\Delta_{\lambda}^{n}h=0$ then with h=\Bigl{(}{\displaystyle h^{o}\atop\displaystyle g}\Bigr{)} as above, we get the equation

[TABLE]

which has the same solution as in Lemma 3.1. Thus, we have the following general version of Proposition 2.6.

(3.2) Proposition.

A function $f:X\to\mathbb{C}$ satisfies $\Delta_{\lambda}^{n}f=0$ on all of $X$ if and only if $f$ is $\lambda$-harmonic.

For $n\geq 2$, what is more interesting is to assign further boundary conditions. Recall that $\Delta_{\lambda}f$ always vanishes on $\partial X$. The analogue of the Dirichlet problem is the Riquier problem of order $n$. We assign $n$ boundary functions $g_{1}\,,\dots,g_{n}:\partial X\to\mathbb{C}$ and look for a function $f:X\to\mathbb{C}$ such that we have a “tower” of boundary value problems for functions $f_{n}\,,f_{n-1}\,,\dots,f_{1}=f:X\to\mathbb{C}$ as follows:

[TABLE]

(3.4) Theorem.

For $\lambda\in\mathsf{res}(P_{X^{o}})\,$, the unique solution $f=f_{1}$ of (3.3) is given by

[TABLE]

where $\mathbf{G}(\lambda)^{r}$ is the $r$-th matrix power of $\mathbf{G}(\lambda)$.

Proof.

We use induction on $n$. For $n=1$, this is Lemma 3.1. Suppose the statement is true for $n-1$. The function $f_{2}$ is the solution of the Riquier problem of order $n-1$ for the boundary functions $g_{2}\,,\dots,g_{n}$. By the induction hypothesis,

[TABLE]

and this is the unique solution. The last one of the “tower” of equations (3.3) is

[TABLE]

This can be rewritten as

[TABLE]

Inserting the solution for $f_{2}^{o}$ and multiplying by $\mathbf{G}(\lambda)$, we get the solution for $f$, and it is unique. ∎

Note that the solution $f$ does not satisfy (1.5) on all of $X^{o}$. This is due to the fact that our discrete Laplacian is not infinitesmial. Let

[TABLE]

the set of all points in $X$ from which $\partial X$ can be reached in $n-1$ or less steps. Then $\Delta_{\lambda}^{n}f=0$ only on the $n$-th interior $X\setminus\partial^{n}X$, while the values on $\partial^{n}X$ depend on the boundary functions $g_{1}\,,\dots,g_{n}\,$.

The functions $\lambda\mapsto G(x,y|\lambda)$ are rational, and the union of the set of their poles is $\mathsf{spec}(P_{X^{o}})$. For $\lambda\in\mathsf{res}(P_{X^{o}})$, we can differentiate the identity $\lambda\cdot\mathbf{G}(\lambda)-P\,\mathbf{G}(\lambda)=I_{X^{o}}$

$k$ times, and Leibniz’ rule yields

[TABLE]

where $\mathbf{G}^{(r)}(\lambda)$ is the (elementwise) $r$-th derivative of $\mathbf{G}(\lambda)$ with respect to $\lambda$. From this, we get recursively for the matrix powers of $\mathbf{G}(\lambda)$

[TABLE]

We can insert this in the formula of Theorem 3.4 for an alternative form of the solution of the Riquier problem.

4. Comparison with the case of infinite trees; examples

We now want to relate the preceding material, and in particular Theorem 3.4, with the potential theory of countable Markov chains, and more specifically, with Martin boundary theory and $\lambda$-polyharmonic functions on trees, as studied in [12]. We choose and fix a reference point (origin) $o\in X^{o}$ and consider the rational functions $\lambda\mapsto F(o,w|\lambda)$ of (2.4) for $\lambda\in\mathsf{res}(P_{X^{o}})$ and $w\in\partial X$. They have (at most) finitely many zeros. Let

[TABLE]

Every positive real $\lambda>\rho(P)$ belongs to $\mathsf{res}^{*}(P_{X^{o}})$, in particular, $\lambda=1$. For $\lambda\in\mathsf{res}^{*}(P_{X^{o}})$, we define the $\lambda$-Martin kernel

[TABLE]

The function $x\mapsto K^{(X)}(x,w|\lambda)$ is the unique solution of the $\lambda$-Dirichlet problem of Lemma 3.1 with value $1$ at the root $o$ and the boundary function $g_{v}$ proportional to $\delta_{w}\,$, that is $g_{w}(v)=\delta_{w}(v)/F(o,w|\lambda)$. Thus, for a generic boundary function $g:\partial X\to\mathbb{C}$, we can write the solution of the $\lambda$-Dirichlet problem for $x\in X^{o}$ as

[TABLE]

The integral notation indicates that we think of $\nu=\nu_{g}$ as a complex distribution on $\partial X$. In the same way, the solution of the Riquier problem in Theorem 3.4 can be written as

[TABLE]

Now let us look at the case of a nearest neighbour transition operator $P=P_{T}$ on a countable tree $T$ without leaves (i.e., vertices distinct from $o$ which have just one neighbour): there, the geometric boundary is attachted to the tree “at infinity”, and there is no “interior” of $T$ which appears as a subset of the vertex set: the interior is $T$ itself. The Martin kernel $K^{(T)}(x,\xi|\lambda)$ is defined for $x\in T$ and $\xi\in\partial T$, and it satisfies $(\lambda\cdot I-P)K^{(T)}(\cdot,\xi|\lambda)=0$, without any restriction to a sub-matrix such as $P_{X^{o}}\,$. In this setting, [12, Thm. 5.4] says that any $\lambda$-polyharmonic function $f$ of order $n$ on $T$ has a unique representation of the form

[TABLE]

and $\nu_{1}\,,\dots,\nu_{n}$ are distributions on $\partial T$. The normalization is slightly different here from the one chosen in [12], and in particular,

[TABLE]

Let us compare the kernels $K^{(X)}_{r}$ and $K^{(T)}_{r}$. We have

[TABLE]

The only, but crucial difference is that in the first of the two identities, we may multiply from the left by $\mathbf{G}^{(X)}(\lambda)^{r-1}=(\lambda\cdot I_{X^{o}}-P_{X^{o}})^{-(r-1)}$. In the second identity, we may not multiply by $\mathbf{G}^{(T)}(\lambda)^{r-1}$, where $\mathbf{G}^{(T)}(\lambda)=(\lambda\cdot I_{T}-P_{T})^{-1}$ is the resolvent of $P$ as an operator on the Hilbert space $\ell^{2}(T,m)$, with the weights $m(x)$ analogous to Example 1.1 above. Indeed, $K^{(T)}(\cdot,\xi|\lambda)$ does in general not belong to $\ell^{2}(T,m)$.

“Forward only” Laplacians on finite and infinite trees

We now consider a class of examples which constitute the finite analogue of [12, §6]. They were also studied, from the viewpoint of Information Theory, by Hirschler and Woess [9].

In order to carry the above comparison with the infinite case a bit further, we need some more details on the geometry of an infinite tree $T$ with root $o$. We assume that $T$ is locally finite and has no leaves. Each vertex $x\neq o$ has a unique predecessor $x^{-}$, its neighbour which is closer to $o$. For each $x\in T$ there is the unique geodesic path $\pi(o,x)=[o=x_{0}\,,x_{1}\,,\dots,x_{n}=x]$ from $o$ to $x$, where $x_{k}^{-}=x_{k-1}$ for $k=1,\dots,n$. In this case, $|x|=n$ is the length of $x$.

The boundary at infinity $\partial T$ of $T$ consists of all geodesic rays $\xi=[o=x_{0}\,,x_{1}\,,x_{2}\,,\dots]$, where $x_{k}^{-}=x_{k-1}$ for $k\geq 1$. For a vertex $x\in T$, we define the boundary arc

[TABLE]

The collection of all $\partial_{x}T\,$, $x\in T$, is the basis of a topology on $\partial T$, which thus becomes a compact, totally disconnected space, and each boundary arc is open and compact. We now take a Borel probability measure $\mathbb{P}$ on $\partial T$ which is supported by the entire boundary, that is, $\mathbb{P}(\partial_{x}T)>0$ for all $x\in T$. It induces a forward only Markov operator on $T$, as follows.

[TABLE]

Conversely, if we start with transition probabilities $p(x,y)$ such that $p(x,y)>0$ precisely when $y^{-}=x$, then we can construct $\mathbb{P}$ on $\partial_{x}T$ by setting

[TABLE]

This determines $\mathbb{P}$ on the Borel $\sigma$-algebra of $\partial T$.

More generally, a distribution on $\partial T$ is a set function

[TABLE]

If $\nu$ is non-negative real, then it extends uniquely to a Borel measure on $\partial T$. A locally constant function $\varphi$ on $\partial T$ is one such that every $\xi\in\partial T$ has a neighbourhood on which $\varphi$ is constant. Thus, one can write it as a finite linear combination of boundary arcs

[TABLE]

and we can define

[TABLE]

Indeed, in this way, the space of all distributions is the dual of the linear space of all locally constant functions on $\partial T$.

Now take $\lambda\in\mathbb{C}\setminus\{0\}$. Following [12, §6], the $\lambda$-Martin kernel on $T$ is

[TABLE]

For fixed $x$, the function $\xi\mapsto K^{(T)}(x,\xi|\lambda)$ and its derivatives with respect to $\lambda$ are locally constant, whence they can be integrated against distributions on $\partial T$. According to (4.4), we get

[TABLE]

and every $\lambda$-polyharmonic function of order $n$ on $T$ has a unique representation

[TABLE]

where the $\nu_{r}=\nu_{r}^{(T)}$ ($r=1,\dots,n$) are distributions on $\partial T$.

We now consider the finite situation. The graph $X$ under consideration is a finite subtree of $T$ with the same root $o$. The boundary consists of the leaves of the tree:

[TABLE]

We suppose that $\partial X$ is a section of $T$ in the sense of [9]: For every $\xi\in\partial T$, the geodesic ray starting from $o$ that represents $\xi$ intersects $\partial X$ in a unique vertex. (A typical special case is the one where $\partial X=\{x\in T:|x|=L\}$ with $L\in\mathbb{N}$.) For each $x\in X$, we define the finite version of the boundary arc rooted at $x$ as

[TABLE]

In particular, $\partial_{o}X=\partial X$, and $\partial_{w}X=\{w\}$ for $w\in\partial X$.

We consider the restriction to $X$ of the given forward transistion matrix $P_{T}$ on $T$. That is,

[TABLE]

while $p_{X}(x,y)=0$ in all other cases. Exactly as on the whole tree, we have for $x,y\in X$

[TABLE]

The matrix $P_{X^{o}}$ is nilpotent, so that $\mathsf{spec}(P)=\{0,1\}$, and the algebraic multiplicities of those two eigenvalues are $|X^{o}|$ and $|\partial X|$, respectively. For $\lambda\in\mathbb{C}\setminus\{0\}=\mathsf{res}(P_{X^{o}})$ and $x,y\in X^{o}$, we have

[TABLE]

Therefore in this example, the right hand side of (3.6) is obtained by

[TABLE]

if $x\in\pi(o,y)\,$. We note that $\mathsf{res}^{*}(P_{X^{o}})=\mathsf{res}(P_{X^{o}})$ and that $F(o,w|\lambda)=\lambda^{-|w|}\,\mathbb{P}(\partial_{w}T)$ for $w\in\partial X$. We can now compute the kernels $K^{(X)}_{r}$ of (4.3) as follows:

[TABLE]

Then, given boundary functions $g_{1}\,,\dots,g_{n}\,$, the associated solution of the Riquier problem reads

[TABLE]

Now consider (4.6) and the fact that $P_{X^{o}}$ is the restriction of $P_{T}$ to $X^{o}$. In spite of this, when $n\geq 2$ we see that for $w\in\partial X$, the function $x\mapsto K_{n}^{(X)}(x,w|\lambda)$ is not the restriction to $X^{o}$ of $x\mapsto K_{n}^{(T)}(x,\xi|\lambda)$, where $\xi\in\partial_{w}X$. (The value is the same for every such $\xi$, when $x\in X^{o}$.) For a closer look, fix $\xi\in\partial_{w}T$ and let $f(x)=K_{n}^{(T)}(x,\xi|\lambda)$ for $x\in X$. This function solves the Riquier problem on $X$ with boundary functions

[TABLE]

or, equivalently, with boundary measures on $\partial X$

[TABLE]

Indeed, verification of

[TABLE]

leads to known combinatorial identity

[TABLE]

in which $|w|$ and $|x|$ can be arbitrary integers with $|w|>|x|\geq 0$.