## Method of images for one-dimensional discrete symmetric random walk under a reflecting barrier

Kazuhiko Seki ∗

GZR, National Institute of Advanced Industrial Science and Technology (AIST), Onogawa 16-1 AIST West, Ibaraki 305-8569, Japan

(Dated: August 5, 2024)

## Abstract

The transition probability for a one-dimensional discrete symmetric random walk under a reflecting barrier was once given by the method of images. [S. Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943).] However, several inconsistencies have been reported when the method of images is applied in cases where a reflecting barrier is considered, even after the exact solution has been obtained. Here, we explicitly show that the method of images becomes applicable if the image position is shifted.

In a well-known review article, Chandrasekhar presented the transition probability (Green function) constructed by the method of images for a reflecting boundary condition as well as for an absorbing boundary condition. 1 Chandrasekhar's transition probability for a reflecting boundary condition led to confusion; some authors accepted it , whereas others noted several 2 inconsistencies in it. 3-5 An attempt to modify the boundary condition to match Chandrasekhar's transition probability has been reported. 5 However, even when this modification is adopted, the transition probability at the reflecting site with the modified transition rate constant should be set to be one-half of Chandrasekhar's transition probability. 5 van Kampen and Oppenheim derived the exact solution under a reflecting (natural) boundary condition using normal modes; however, the solution involves an infinite summation. 6 Kampen and Oppenheim also introduced an approximation applicable when the lattice spacing is smaller than the distances of interest and noted that the probability density no longer has a zero slope at the site where the natural boundary condition is imposed. 6 Researchers have also obtained the exact solution by solving a recursive equation and have given an intuitive interpretation based on the method of images. 7,8 However, the shift of the image was not mathematically elucidated. Here, we explicitly show that the method of images becomes applicable if the image position is shifted, and derive the transition probability expressed using the shifted image.

As a preliminary step, we first present a known result for an unbiased random walk in an infinite one-dimensional lattice without a boundary. The result helps to construct the transition probability using the method of images by properly locating the image position to satisfy the reflecting (natural) boundary condition. We consider a random walk along the x -axis, with the lattice spacing denoted by b . When the total transition rate constant is denoted by γ , the transition probability for a particle at x at time t ∼ t + dt starting from the initial position x i [ G 0 ( x, x , t i )] in an unbiased random walk in an infinite one-dimensional lattice satisfies

$$\frac { \partial } { \partial t } G _ { 0 } ( x, x _ { i }, t ) = \frac { \gamma } { 2 } G _ { 0 } ( x - b, x _ { i }, t ) - \gamma G _ { 0 } ( x, x _ { i }, t ) + \frac { \gamma } { 2 } G _ { 0 } ( x + b, x _ { i }, t ).$$

We denote the Laplace transform of the time-dependent function f ( ) by t ˆ ( f s ). G 0 ( x, x , t i ) in the Laplace domain can be expressed as 9

$$\hat { G } _ { 0 } ( x, x _ { i }, s ) = \frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \frac { \exp [ - i k ( x - x _ { i } ) ] } { s + \gamma \left [ 1 - \cos \left ( k b \right ) \right ] }.$$

By noting that G 0 ( x, x , t i ) is a function of x -x i alone, we express it as G 0 ( x -x , t i ) ≡ G 0 ( x, x , t i ) when needed. By the inverse Laplace transform of Eq. (2), G 0 ( x -x , t i ) is obtained as

$$G _ { 0 } ( x - x _ { i }, t ) = \frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \cos [ k ( x - x _ { i } ) ] \exp \left [ - \gamma \left ( 1 - \cos \left ( k b \right ) \right ) t \right ]$$

$$= \exp \left ( - \gamma t \right ) I _ { | x - x _ { 1 } | / b } \left ( \gamma t \right ),$$

where the integral representation of the modified Bessel function of the first kind [ I n ( z )] for an integer index is used.

In the Laplace domain, Eq. (1) can be rewritten as

$$( s + \gamma _ { n } ) \, \hat { G } _ { 0 } ( x, x _ { i }, s ) - ( \gamma _ { n } / 2 ) \left [ \hat { G } _ { 0 } ( x - b, x _ { i }, s ) + \hat { G } _ { 0 } ( x + b, x _ { i }, s ) \right ] = \delta _ { x, x _ { i } }.$$

We introduce G x,x ,t b ( i ) ≡ [ ˆ G 0 ( x -b, x , t i ) + ˆ G 0 ( x + b, x , t i ) ] / 2, which satisfies

$$2 \hat { G } _ { b } ( x, x _ { i }, s ) & = \frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \frac { \cos [ k ( x + b - x _ { i } ) ] + \cos [ k ( x - b - x _ { i } ) ] } { s + \gamma _ { n } \, [ 1 - \cos \, ( k b ) ] } \\ & = 2 \frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \frac { \cos [ k ( x - x _ { i } ) ] \cos ( k b ) } { s + \gamma _ { n } \, [ 1 - \cos \, ( k b ) ] },$$

where cos( A -B ) + cos( A + B ) = 2 cos( A ) cos( B ) is used. By introducing Eq. (6) into Eq. (5), we can confirm that Eq. (2) is a solution because we have

$$\frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \cos [ k ( x - x _ { i } ) ] = \delta _ { x, x _ { i } },$$

where ( x -x i ) /b is an integer.

Now, we construct the transition probability using the method of images [ G mrf ( x, x , t i )] to satisfy the reflecting (natural) boundary condition at site 0: 6,10

$$\frac { \partial } { \partial t } G _ { \text{mrf} } ( 0, x _ { i }, t ) = - \frac { \gamma } { 2 } G _ { \text{mrf} } ( 0, x _ { i }, t ) + \frac { \gamma } { 2 } G _ { \text{mrf} } ( b, x _ { i }, t ).$$

The natural boundary condition is known to be equivalent to considering an auxiliary boundary condition at site -b : 6,10

$$\frac { \gamma } { 2 } G _ { \text{mrf} } ( - b, x _ { i }, t ) = \frac { \gamma } { 2 } G _ { \text{mrf} } ( 0, x _ { i }, t ).$$

By adding Eq. (9) with

$$\frac { \partial } { \partial t } G _ { \text{mrf} } ( 0, x _ { i }, t ) = \frac { \gamma } { 2 } G _ { \text{mrf} } ( - b, x _ { i }, t ) - \gamma G _ { \text{mrf} } ( 0, x _ { i }, t ) + \frac { \gamma } { 2 } G _ { \text{mrf} } ( b, x _ { i }, t ),$$

we obtain the natural boundary condition given by Eq. (8). Therefore, imposing the natural boundary condition given by Eq. (8) is equivalent to imposing Eq. (9) at the auxiliary site -b (see Fig. 1). 6,10 Instead of directly considering the natural boundary condition, we consider the auxiliary boundary condition given by Eq. (9): G mrf ( -b, x , t i ) = G mrf (0 , x , t i ).

FIG. 1. Schematics of random walk model with reflecting boundary located at -m . The lattice spacing is denoted by b . The total transition rate is denoted by γ . (a) The natural boundary is located at x = 0. (b) The boundary condition G mrf ( -b, x , t i ) = G mrf (0 , x , t i ) is imposed at the auxiliary site -b . The boundary condition at the auxiliary site [(b)] is equivalent to the natural boundary condition [(a)]. 6,10

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00926v1-figure-28c21a2c-8e8c-4aff-bd81-4d47c6fd80c0.jpg)

By considering the natural boundary located at site 0 and the auxiliary boundary located at site -b , the location of the reflecting boundary might be between site 0 and site -b . We denote the location of the reflecting boundary by -m (see Fig. 1). By considering -m -[ x i - -( m )] = -x i -2 m , we can express the transition probability constructed by the method of images as

$$G _ { \text{mrf} } ( x, x _ { i }, t ) & \equiv G _ { 0 } ( x, x _ { i }, t ) + G _ { 0 } ( x, - x _ { i } - 2 m, t ). & ( 1 1 )$$

In the boundary condition given by Eq. (9), we have G mrf ( -b, x , t i ) = G 0 ( x i + b, t ) + G 0 ( x i -b +2 m,t ) and G mrf (0 , x , t i ) = G 0 ( x , t i ) + G 0 ( x i +2 m,t ). Using Eq. (2) and cos( A ) + cos( B ) = 2 cos[( A + B / ) 2] cos[( A -B / ) 2], we obtain

$$\hat { G } _ { \text{mrf} } ( - b, x _ { i }, s ) = \frac { b } { 2 \pi } \int _ { \stackrel { - \pi / b } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { - } { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stackrel { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { \stack { { \ \stack { \stack { \stack { \ \stack { \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \stack { \ \ \ \ \ \ \ \ \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } { \, \, \,$$

$$\hat { G } _ { \text{mrf} } ( 0, x _ { i }, s ) = \frac { b } { 2 \pi } \int _ { - \pi / b } ^ { \pi / b } d k \, \frac { 2 \cos [ k ( x _ { i } + m ) ] \cos ( k m ) } { s + \gamma \left [ 1 - \cos \left ( k b \right ) \right ] }.$$

We note that Eq. (12) is equivalent to Eq. (13) when m = b/ 2 is set; the boundary condition given by Eq. (9) is therefore satisfied. As a result,

$$G _ { \text{mrf} } ( x, x _ { i }, t ) \equiv G _ { 0 } ( x, x _ { i }, t ) + G _ { 0 } ( x, - x _ { i } - b, t )$$

is the solution of the unbiased lattice random walk with the natural boundary condition at site 0. Using Eq. (4), we can explicitly express G mrf ( x, x , t i ) as

$$G _ { \text{mrt} } ( x, x _ { i }, t ) = \exp \left ( - \gamma t \right ) \left [ I _ { | x - x _ { i } | / b } \left ( \gamma t \right ) + I _ { ( x + x _ { i } + b ) / b } \left ( \gamma t \right ) \right ].$$

For simplicity, we have considered a simple lattice random walk; however, Eq. (14) can be used for ˆ G 0 ( x, x , s i ) for a continuous-time random walk instead of Eq. (2). 9,11 Previously, van Kampen and Oppenheim showed that the probability density no longer has a zero slope at x = -b (nor at x = 0) but at x = -b/ 2. 6 Their conclusion is consistent with Eq. (14). During the final drafting of this article, we became aware that Eq. (14) was proved recently in a very different manner in a different context. 12 Previously, Eq. (15) was derived by solving a recursive equation. 7 An interpretation of Eq. (15) based on the method of images was also given using an intuitive argument for the reflecting boundary. The boundary condition given by G mrf ( -b, x , t i ) = G mrf (0 , x , t i ) implies that a walker just arriving at the auxiliary site is placed back at the physical site located at 0. It also implies that the rational position of the reflecting barrier, although we are considering the discrete lattice, might be the midpoint between the auxiliary site at -1 and the physical site at 0. 8 Here, we derived mathematically that the reflecting barrier should be located at the site -1 2 to satisfy / G mrf ( -b, x , t i ) = G mrf (0 , x , t i ).

What is the boundary condition satisfied by Chandrasekhar's transition probability? By setting m = 0, which is appropriate for Chandrasekhar's transition probability, we find ˆ G mrf ( -b, x , s i ) = k a ( s G ) ˆ mrf (0 , x , s i ) in the Laplace domain with k a ( s ) = ˆ G mrf ( -b, x , s i ) /G ˆ mrf (0 , x , s i ); the value of k a ( s ) satisfyies 0 &lt; k a ( s ) &lt; 1 because of cos( kb ) &lt; 1 in ˆ G mrf ( -b, x , s i ) with m = 0 and indicates a partially absorbing boundary condition in the Laplace domain. Therefore, Chandrasekhar's transition probability is a solution for the partially absorbing boundary condition in the Laplace domain. The solution for the boundary condition given by ˆ G mrf ( -b, x , s i ) = α s G ( ) ˆ mrf (0 , x , s i ) can be constructed by an infinite summation using ˆ G 0 ( x + x i + bj, s ) with j being a non-negative integer for α s ( ) &lt; 1. 7,8 Because the equations

are linear, the transition probability for a perfectly reflecting boundary condition can be expressed using a linear combination of Chandrasekhar's transition probability and the solution for α s ( ) = 1 -k a ( s ). When the continuous limit in space is taken, the solution should converge to a unique solution given by the continuous limit of Eq. (15).

Chandrasekhar also presented the transition probability for a perfectly absorbing boundary condition using the method of images, and the solution should be exact; indeed, G 0 ( | x -x i | , t ) -G 0 ( x + x , t i ) should be zero at x = 0 (the absorbing site). 1,7 Note that only the absorbing site is considered when constructing the transition probability to satisfy the perfectly absorbing boundary condition by the method of images. However, the transition probability for the perfectly reflecting boundary condition is constructed using the method of images on the basis of the slope of the density at the reflecting site, where the definition of the slope requires two sites. This might be the reason for the shift in the transition probability constructed by the method of images when considering the reflecting boundary condition.

## ACKNOWLEDGMENTS

I would like to thank Professor Hisao Hayakawa (Kyoto University) and Professor Yoshihiko Hasegawa (The University of Tokyo) for their interest and for discussion during the workshop on 'Frontiers in Non-equilibrium Physics 2024' (YITP-T-24-01) at Yukawa Institute for Theoretical Physics (Kyoto University, Japan). I would like to acknowledge the warm hospitality during our stay in the YITP.

∗ k-seki@aist.go.jp

1 S. Chandrasekhar, Rev. Mod. Phys. 15 , 1 (1943).

2 F. C. Goodrich, J. Chem. Phys. 22 , 588 (1954).

3 K. Razi Naqvi, K. Mork, and S. Waldenstrøm, Chem. Phys. Lett. 92 , 160 (1982).

4 M. Orlowski, Phys. Status Solidi B 239 , 158 (2003).

5 S. K. Hansen, Phys. Status Solidi B 248 , 2112 (2011).

- 6 N. G. van Kampen and I. Oppenheim, J. Math. Phys. 13 , 842 (1972).
- 7 E. W. Montroll, ENERGETICS IN METALLURGICAL PHENOMENA-Volume-3 , edited by W. M. MUELLER (GORDON &amp; BREACH, NEWYORK, 1967) p. 123.
- 8 J. Schwarz, Martin and D. Poland, J. Chem. Phys. 63 , 557 (1975).
- 9 B. Hughes, Random Walks and Random Environments: Random walks , Oxford science publications No. 1 (Clarendon Press, 1995).
- 10 N. G. van Kampen, Stochastic processes in physics and chemistry , 3rd ed. (North Holland, Amsterdam, 2007).
- 11 E. W. Montroll and G. H. Weiss, J. Math. Phys. 6 , 167 (1965).
- 12 I. Corwin and H. Shen, Commun. Pure Appl. Math. 71 , 2065 (2018).
Kazuhiko Seki
Method of images for one-dimensional discrete symmetric random walk under a reflecting barrier
The transition probability for a one-dimensional discrete symmetric random
walk under a reflecting barrier was once given by the method of images. [S.
Chandrasekhar, Rev. Mod. Phys. 15, 1 (1943).] However, several inconsistencies
have been reported when the method of images is applied in cases where a
reflecting barrier is considered, even after the exact solution has been
obtained. Here, we explicitly show that the method of images becomes applicable
if the image position is shifted.