enumeration of lattice walks (derivation of formulas)

enumeration of lattice walks (derivation of formulas)

Square lattice

The method of counting walks hinges on the

concept

of

walk

generating function

: a

Laurent polynomial

F

⁢

(

a

,

b

)

in two

variables

a

,

b

such that

[

a

x

⁢

b

y

]

⁢

F

n

equals the number of walks from

(

0

,

0

)

to

(

x

,

y

)

∈

ℤ

×

ℤ

, where

[

a

x

⁢

b

y

]

⁢

P

⁢

(

a

,

b

)

denotes the

coefficient

of

a

x

⁢

b

y

in

P

⁢

(

a

,

b

)

.

The idea is to encode each possible step that can be taken from a vertex to an

adjacent vertex

as a term in

F

⁢

(

a

,

b

)

. In the case of the square lattice, there are four possible steps from a given vertex: left, right, up, or down. They are encoded as

a

,

1

/

a

,

b

and

1

/

b

. So the walk generating

polynomial

is

F

⁢

(

a

,

b

)

=

a

+

1

/

a

+

b

+

1

/

b

. Each term in the expansion of

F

n

corresponds to a walk of length

n

starting at a fixed vertex. For example, a term in

F

4

such as

a

⋅

1

/

b

⋅

1

/

a

⋅

1

/

a

=

a

-

1

⁢

b

-

1

would correspond to the walk of length 4 consisting of the following sequence of steps: right, down,left, left. Moreover, when the term is written as a

monomial

a

x

⁢

b

y

, the exponents

x

,

y

are precisely the coordinates of the end point of the walk (assuming the origin as the starting point of the walk). Thus in this example the walk would end at vertex

(

-

1

,

-

1

)

. When all like terms in the expansion are collected it becomes clear that

[

a

x

⁢

b

y

]

⁢

F

n

counts the number of walks from

(

0

,

0

)

to

(

x

,

y

).

To obtain an explicit formula for the number of walks from the origin to

(

x

,

y

)

we factor

F

=

(

a

+

b

)

⁢

(

1

+

1

/

a

⁢

b

)

, raise to the n-th power, use the binomial expansion, and collect like terms:

F

n

=

(

a

+

b

)

n

⁢

(

1

+

1

/

a

⁢

b

)

n

=

(

∑

i

(

n

i

)

⁢

a

i

⁢

b

n

-

i

)

⁢

(

∑

j

(

n

j

)

⁢

a

-

j

⁢

b

-

j

)

=

∑

(

i

,

j

)

(

n

i

)

⁢

(

n

j

)

⁢

a

i

-

j

⁢

b

n

-

i

-

j

The number of walks

w

⁢

(

x

,

y

,

n

)

of length

n

is the coefficient of

a

x

⁢

b

y

, namely

w

⁢

(

x

,

y

,

n

)

=

∑

i

-

j

=

x

,

n

-

i

-

j

=

y

(

n

i

)

⁢

(

n

j

)

If

n

and

x

+

y

have different parity, there are no integers

i

,

j

satisfying the summation constraints, and thus the coefficient of

a

x

⁢

b

y

is zero, indicating that there are no walks of length

n

from the origin to

(

x

,

y

)

. If, on the other hand,

n

and

x

+

y

have the same parity, the constraints admit only one solution

(

i

,

j

)

, namely

i

=

(

n

+

x

-

y

)

/

2

and

j

=

(

n

-

x

-

y

)

/

2

. Thus, in this case,

w

⁢

(

x

,

y

,

n

)

=

(

n

n

+

x

-

y

2

)

⁢

(

n

n

-

x

-

y

2

)

Triangular lattice

In the case of a triangular lattice, the walk generating function is

F

=

a

+

1

/

a

+

b

+

1

/

b

+

b

/

a

+

a

/

b

, where each term correspond to the possible steps that can be taken from a vertex to an adjacent vertex: right, left, up, down, diagonally up, diagonally down. It can also be written as

F

=

(

1

/

a

+

1

/

b

)

⁢

(

1

+

a

)

⁢

(

1

+

b

)

-

2

. The number of walks of length

n

from the origin to

(

x

,

y

)

is given by the coeeficient of

a

x

⁢

b

y

in

F

n

. Using the binomial expansion, one obtains:

F

n

=

[

(

1

/

a

+

1

/

b

)

⁢

(

1

+

a

)

⁢

(

1

+

b

)

-

2

]

n

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

(

1

/

a

+

1

/

b

)

k

⁢

(

1

+

a

)

k

⁢

(

1

+

b

)

k

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

∑

i

(

k

i

)

⁢

(

1

/

a

)

k

-

i

⁢

(

1

/

b

)

i

⁢

∑

j

(

k

j

)

⁢

a

j

⁢

∑

j

′

(

k

j

′

)

⁢

b

j

′

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

∑

(

i

,

j

,

j

′

)

(

k

i

)

⁢

(

k

j

)

⁢

(

k

j

′

)

⁢

a

j

-

k

+

i

⁢

b

j

′

-

i

The number of walks

(

x

,

y

,

n

)

of length

n

from the origin to

(

x

,

y

)

is the coefficient of

a

x

⁢

b

y

:

w

⁢

(

x

,

y

,

n

)

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

∑

j

-

k

+

i

=

x

,

j

′

-

i

=

y

(

k

i

)

⁢

(

k

j

)

⁢

(

k

j

′

)

Adding the summation constraints for the inner sum, we get

j

+

j

′

=

k

+

x

+

y

, thus

j

′

=

k

+

x

+

y

-

j

and

i

=

j

′

-

y

=

k

+

x

-

j

. Therefore

w

⁢

(

x

,

y

,

n

)

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

∑

j

(

k

j

)

⁢

(

k

j

-

x

-

y

)

⁢

(

k

j

-

x

)

In particular, the number of

closed walks

of length

n

starting and ending at

(

0

,

0

)

is

w

⁢

(

0

,

0

,

n

)

=

∑

k

(

-

2

)

n

-

k

⁢

(

n

k

)

⁢

∑

j

(

k

j

)

3

Honeycomb lattice

First we observe that the

honeycomb

lattice is a

bipartite graph

: a graph with

vertex set

V

=

V

1

∪

V

2

such that every edge is from

V

1

to

V

2

. In the case of the honeycomb lattice,

V

=

ℤ

×

ℤ

=

V

1

∪

V

2

where

V

1

=

{

(

x

,

y

)

∈

V

:

x

+

y

≡

0

(

mod

3

)

}

and

V

2

=

{

(

x

,

y

)

∈

V

:

x

+

y

≡

1

(

mod

3

)

}

. Every edge in the honeycomb lattice joins a vertex in

V

1

to a vertex in

V

2

. Notice that a walk starting at a vertex in

V

1

will end at a vertex in

V

1

(resp.

V

2

) if it has even (resp. odd) length.

From a vertex

(

i

,

j

)

∈

V

1

the possible steps are up, right, or diagonally down. These are encoded by

F

=

a

+

b

+

1

/

a

⁢

b

; on the other hand, if

(

i

,

j

)

∈

V

2

, the possible steps from

(

i

,

j

)

are left, down, or diagonally up, encoded by

G

=

1

/

a

+

1

/

b

+

a

⁢

b

. Since a walk of even length

n

from say

(

0

,

0

)

must alternate between points in

V

1

and points in

V

2

, each such walk corresponds to a term in the expansion of

(

F

⁢

G

)

n

/

2

, and

w

⁢

(

x

,

y

,

n

)

=

[

a

x

⁢

b

y

]

⁢

(

F

⁢

G

)

n

/

2

. Since

F

⁢

G

=

1

+

(

a

+

b

)

⁢

(

a

+

1

/

a

⁢

b

)

⁢

(

b

+

1

/

a

⁢

b

)

, we have

(

F

⁢

G

)

n

/

2

=

∑

k

(

n

k

)

⁢

(

a

+

b

)

k

⁢

(

a

+

1

/

a

⁢

b

)

k

⁢

(

b

+

1

/

a

⁢

b

)

k

=

∑

k

(

n

k

)

⁢

∑

i

(

k

i

)

⁢

a

k

-

i

⁢

b

i

⁢

∑

j

(

k

j

)

⁢

a

k

-

j

⁢

(

1

/

a

⁢

b

)

j

⁢

∑

j

′

(

k

j

′

)

⁢

b

k

-

j

′

⁢

(

1

/

a

⁢

b

)

j

′

=

∑

k

(

n

k

)

⁢

∑

i

,

j

,

j

′

(

k

i

)

⁢

(

k

j

)

⁢

(

k

j

′

)

⁢

a

2

⁢

k

-

i

-

2

⁢

j

-

j

′

⁢

b

k

+

i

-

j

-

2

⁢

j

′

Thus for

n

even,

w

⁢

(

x

,

y

,

n

)

=

[

a

x

⁢

b

y

]

⁢

(

F

⁢

G

)

n

/

2

=

∑

k

(

n

k

)

⁢

∑

2

k

-

i

-

2

j

-

j

′

=

x

,

k

+

i

-

j

-

2

⁢

j

′

=

y

(

k

i

)

⁢

(

k

j

)

⁢

(

k

j

′

)

The constraints under the inner sum imply

j

′

=

k

-

j

-

(

x

+

y

)

/

3

and

i

=

k

-

j

-

x

+

(

x

+

y

)

/

3

. Thus if

x

+

y

≢

0

(

mod

3

)

they are not satisfied by any integers

i

,

j

,

j

′

, and

w

⁢

(

x

,

y

,

n

)

=

0

. This simply reflects the fact that there are no even length walks from a vertex in

V

1

to a vertex in

V

2

. On the other hand, if

x

+

y

≡

0

(

mod

3

)

, using

(

p

q

)

=

(

p

p

-

q

)

we obtain

w

⁢

(

x

,

y

,

n

)

=

∑

k

n

/

2

(

n

/

2

k

)

⁢

∑

j

(

k

j

+

x

-

(

x

+

y

)

/

3

)

⁢

(

k

j

)

⁢

(

k

j

+

(

x

+

y

)

/

3

)

Walks of odd length

n

from

(

0

,

0

)

to

(

x

,

y

)

∈

V

2

can be computed by

w

⁢

(

x

,

y

,

n

)

=

w

⁢

(

x

-

1

,

y

,

n

-

1

)

+

w

⁢

(

x

,

y

+

1

,

n

-

1

)

+

w

⁢

(

x

+

1

,

y

+

1

,

n

-

1

)

, where the terms in the right hand side can be computed by the previous formula.

For

x

=

y

=

0

we obtain the following formula for the number of closed walks of even length

n

starting and ending at

(

0

,

0

)

:

w

⁢

(

0

,

0

,

n

)

=

∑

k

n

/

2

(

n

/

2

k

)

⁢

∑

j

(

k

j

)

3

An alternative expansion yields another formula for

w

⁢

(

x

,

y

,

n

)

that contains a single summation. Writing

F

⁢

G

=

(

1

+

a

2

⁢

b

+

a

⁢

b

2

)

⁢

(

1

+

a

-

2

⁢

b

-

1

+

a

-

1

⁢

b

-

2

)

, we obtain

(

F

⁢

G

)

n

/

2

=

∑

k

(

n

/

2

k

)

⁢

(

a

2

⁢

b

+

a

⁢

b

2

)

k

⁢

∑

k

′

(

n

/

2

k

′

)

⁢

(

a

-

2

⁢

b

-

1

+

a

-

1

⁢

b

-

2

)

k

′

=

∑

k

,

j

(

n

/

2

k

)

⁢

(

k

j

)

⁢

(

a

2

⁢

b

)

k

-

j

⁢

(

a

⁢

b

2

)

j

⁢

∑

k

′

,

j

′

(

n

/

2

k

′

)

⁢

(

k

′

j

′

)

⁢

(

a

-

2

⁢

b

-

1

)

k

′

-

j

′

⁢

(

a

-

1

⁢

b

-

2

)

j

′

=

∑

k

,

j

,

k

′

,

j

′

(

n

/

2

k

)

⁢

(

n

/

2

k

′

)

⁢

(

k

j

)

⁢

(

k

′

j

′

)

⁢

a

2

⁢

k

-

2

⁢

k

′

-

j

+

j

′

⁢

b

k

-

k

′

+

j

-

j

′

The number of walks

w

⁢

(

x

,

y

,

n

)

is the coefficient of

a

x

⁢

b

y

, namely

w

⁢

(

x

,

y

,

n

)

=

∑

2

k

-

2

k

′

-

j

+

j

′

=

x

,

k

-

k

′

+

j

-

j

′

=

y

(

n

/

2

k

)

⁢

(

n

/

2

k

′

)

⁢

(

k

j

)

⁢

(

k

′

j

′

)

The constraints under the summation imply

3

⁢

k

-

3

⁢

k

′

=

x

+

y

and

3

⁢

j

-

3

⁢

j

′

=

2

⁢

y

-

x

, thus

k

′

=

k

-

(

x

+

y

)

/

3

and

j

′

=

j

+

(

x

-

2

⁢

y

)

/

3

. Therefore

w

⁢

(

x

,

y

,

n

)

=

∑

k

(

n

/

2

k

)

⁢

(

n

/

2

k

-

(

x

+

y

)

/

3

)

⁢

∑

j

(

k

k

-

j

)

⁢

(

k

-

(

x

+

y

)

/

3

j

+

(

x

-

2

⁢

y

)

/

3

)

But

∑

j

(

k

k

-

j

)

⁢

(

k

-

(

x

+

y

)

/

3

j

+

(

x

-

2

⁢

y

)

/

3

)

=

(

2

⁢

k

-

(

x

+

y

)

/

3

k

+

(

x

-

2

⁢

y

)

/

3

)

by Vandermonde’s convolution. Therefore

w

⁢

(

x

,

y

,

n

)

=

∑

k

=

0

n

/

2

(

n

/

2

k

)

⁢

(

n

/

2

k

-

(

x

+

y

)

/

3

)

⁢

(

2

⁢

k

-

(

x

+

y

)

/

3

k

+

(

x

-

2

⁢

y

)

/

3

)

Setting

(

x

,

y

)

=

(

0

,

0

)

we obtain the following formula for the number of closed walks starting and ending at

(

0

,

0

)

:

w

⁢

(

0

,

0

,

n

)

=

∑

k

=

0

n

/

2

(

n

/

2

k

)

2

⁢

(

2

⁢

k

k

)

Dice lattice

The dice lattice is also a bipartite graph, with vertex set

V

0

∪

V

1

∪

V

2

, where

V

i

=

{

(

x

,

y

)

∈

ℤ

×

ℤ

:

x

+

y

≡

i

(

mod

3

)

}

, and with edges going from

V

0

to

V

1

∪

V

2

. The possible moves are:

•

From a vertex in

V

0

: east,west,north,south, northeast,southwest, encoded as

a

,

1

/

a

,

b

,

1

/

b

,

a

⁢

b

,

1

/

a

⁢

b

.

•

From a vertex in

V

1

: west, south, northeast, encoded as

1

/

a

,

1

/

b

,

a

⁢

b

.

•

From a vertex in

V

2

: east, north, southwest, encoded as

a

,

b

,

1

/

a

⁢

b

.

Moreover, starting at a vertex in

V

1

, moves

a

,

b

,

1

/

a

⁢

b

take us to a vertex in

V

2

, and moves

1

/

a

,

1

/

b

,

a

⁢

b

take us to a vertex in

V

3

. Of course any move from a vertex in

V

2

∪

V

3

takes us back to a vertex in

V

1

. After these observations and a moment’s reflection, we conclude that for even

n

the number of walks of length

n

starting at

(

0

,

0

)

is the coefficient of

a

x

⁢

b

y

in

[

F

⁢

G

+

G

⁢

F

]

n

/

2

=

2

n

/

2

⁢

(

F

⁢

G

)

n

/

2

where

F

=

a

+

b

+

1

/

a

⁢

b

and

G

=

1

/

a

+

1

/

b

+

a

⁢

b

are defined as in the case of the honeycomb lattice. So the number of walks of even length from

(

0

,

0

)

to

(

x

,

y

)

is

2

n

/

2

times the corresponding number of walks in the honeycomb lattice, and the formulas for the dice lattice follow.

Title

enumeration of lattice walks

(

derivation

of formulas)

Canonical name

EnumerationOfLatticeWalksderivationOfFormulas

Date of creation

2013-03-22 19:21:44

Last modified on

2013-03-22 19:21:44

Owner

csguy (26054)

Last modified by

csguy (26054)

Numerical id

4

Author

csguy (26054)

Entry type

Derivation

Classification

msc 05A15

Generated on Thu Feb 8 20:37:25 2018 by

LaTeXML