This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Log-concavity in planar random walks

Swee Hong Chan Department of Mathematics, UCLA, Los Angeles, CA 90095. [email protected] Igor Pak Department of Mathematics, UCLA, Los Angeles, CA 90095. [email protected]  and  Greta Panova Department of Mathematics, USC, Los Angeles, CA 90089. [email protected]
Abstract.

We prove log-concavity of exit probabilities of lattice random walks in certain planar regions.

1. Introduction

In the study of random walks, there is a fundamental idea based on its Markovian property: when some “life event” happens to the walk, the future trajectory of the walk can be changed, and this transformation can be exploited to obtain both quantitative and qualitative conclusions on the random walk distributions. These “life events” can be rather mundane, for example the first time when the walk returns to the starting point, crosses with some other random walk, enters an obstacle, etc. On the other hand, the conclusions can be quite remarkable and include the classical reflection principle, and the Karlin–McGregor formula also known as the Lindström–Gessel–Viennot lemma in the discrete setting, see §\S5.3.

In this paper we use a variation on this approach for random walks in simply connected planar regions. The conclusion is qualitative in some sense that at the end we prove log-concavity of a certain natural exit probability distribution. We chose to present a discrete version of the result rather than the (somewhat cleaner) continuous version as the former is more powerful and amenable to generalizations (see Section 4), while the latter follows easily by taking limits.

Theorem 1 (Special case of Theorem 7).

Let  Γ2\Gamma\subset\mathbb{Z}^{2}  be a simply connected region in the plane whose boundary  Γ=αη+ηβ\partial\Gamma=\alpha\cup\eta_{+}\cup\eta_{-}\cup\beta  is comprised of two vertical intervals  α,β\alpha,\hskip 0.85355pt\beta, and two  xx-monotone lattice paths  η+,η\eta_{+},\hskip 0.85355pt\eta_{-} , see Figure 1. We assume that  α{x=0}\alpha\subset\{x=0\}  and  β{x=m}\beta\subset\{x=m\}  for some m>0m>0.

Let  {Xt}\{X_{t}\}  be the nearest neighbor lattice random walk which starts at the origin  X0=OαX_{0}=O\in\alpha, and is absorbed whenever  XtX_{t}  tries to exit the region Γ\hskip 0.85355pt\Gamma. Denote by  TT  the first time tt such that  XtβX_{t}\in\beta, and let  p(k)p(k)  be the probability that XT=(m,k)X_{T}=(m,k). Then  {p(k)}\{p(k)\}  is log-concave:

p(k)2p(k+1)p(k1)for all k,  such that (m,k±1)β.p(k)^{2}\,\geq\,p(k+1)\hskip 1.70709ptp(k-1)\qquad\text{for all \ $k\in\mathbb{Z}$, \ such that \ \ $(m,k\pm 1)\hskip 0.85355pt\in\beta$}.
Refer to caption
Figure 1. Random walk XtX_{t} in a region Γ\Gamma as in the theorem.

In particular, the theorem implies that the sequence  {p(k)}\{p(k)\}  is unimodal (see e.g. [Brä]). This also points to the difficulty of proving the result by a direct calculation, as in general there is no natural point at which the probability maximizes. We refer to  p(k)p(k)  as exit probabilities, since one can think of them as probabilities the walks exits the region through different points on the interval β\beta.

Perhaps surprisingly, Theorem 1 is a byproduct of our recent work [CPP2] on the Stanley inequality for the case of posets of width two; in fact, we obtain the inequality as a corollary of the theorem (see §\S5.1).


2. Proof of Theorem 1

We start with the following combinatorial result. Throughout this section, a lattice path is a path in 2\mathbb{Z}^{2}, possibly self-intersecting along vertices and edges.

Lemma 2.

In the conditions of Theorem 1, let  A,BαA,B\in\alpha  and  C,C,D,DβC,C^{\prime},D,D^{\prime}\in\beta  be six points on the boundary of the region Γ\Gamma, such that

A=(0,a),B=(0,b),C=(m,c),C=(m,cr),D=(m,d+r),D=(m,d),A=(0,a),\ B=(0,b),\ C=(m,c),\ C^{\prime}=(m,c-r),\ D^{\prime}=(m,d+r),\ D=(m,d),

and suppose  (ab)(cd)r(a-b)\leq(c-d)-r, where  r>0r>0. Then there is an injection

Φ:{(ξAC,ξBD)}{(ξAC,ξBD)},\Phi:\,\bigl{\{}(\xi_{AC},\xi_{BD})\bigr{\}}\,\longrightarrow\,\bigl{\{}(\xi_{AC^{\prime}},\xi_{BD^{\prime}})\bigr{\}},

where

\circ  ξAC:AC\xi_{AC}:A\to C,  ξBD:BD\xi_{BD}:B\to D,  ξAC:AC\xi_{AC^{\prime}}:A\to C^{\prime}, and  ξBD:BD\xi_{BD^{\prime}}:B\to D^{\prime}  are lattice paths which lie inside Γ\Gamma,

\circ  the sum of numbers of horizontal edges in  ξAC\xi_{AC}  and  ξBD\xi_{BD}  which project onto  [i,i+1][i,i+1],  is equal to

the sum of numbers of horizontal edges in  ξAC\xi_{AC^{\prime}}  and  ξBD\xi_{BD^{\prime}}  which project onto  [i,i+1][i,i+1], for all  0im10\leq i\leq m-1,

\circ  the sum of numbers of vertical edges in  ξAC\xi_{AC}  and  ξBD\xi_{BD}  which project onto  jj,  is equal to

the sum of numbers of vertical edges in  ξAC\xi_{AC^{\prime}}  and  ξBD\xi_{BD^{\prime}}  which project onto  jj, for all  0jm0\leq j\leq m.

Here by project we mean the vertical projection onto the xx-axis. By adding the edge and vertex count equalities above over all ii and jj, we conclude that injection Φ\hskip 0.85355pt\Phi  preserves the total length of these paths:

|ξAC|+|ξBD|=|ξAC|+|ξBD|.\bigl{|}\xi_{AC}\bigr{|}\hskip 1.70709pt+\hskip 1.70709pt\bigl{|}\xi_{BD}\bigr{|}\,=\,\bigl{|}\xi_{AC^{\prime}}\bigr{|}\hskip 1.70709pt+\hskip 1.70709pt\bigl{|}\xi_{BD^{\prime}}\bigr{|}\hskip 0.85355pt.
Proof of Lemma 2.

The proof is an explicit construction of the injection Φ\Phi, illustrated in Figure 2. See also Figure 3 in the next section for a detailed construction of each step.

(0)(0)  We start with the lattice paths  ξAC\xi_{AC}  and  ξBD\xi_{BD}  drawn inside the region Γ\Gamma. Note that by the definition of the points in the lemma, we have  |CC|=|DD|=r|CC^{\prime}|=|D^{\prime}D|=r.

(1)(1)  Let  :=b+cdr\ell:=b+c-d-r  and let  B:=(0,)B^{\prime}:=(0,\ell). By the assumption in the lemma, BB^{\prime} lies above AA on the line spanned by α\alpha. Denote by  η^\widehat{\eta}_{-}  the lattice path  BDB\to D  formed by following interval α\alpha down from BB, then path  η\eta_{-}  and then interval β\beta up to DD. Let  χ\chi  be the lattice path obtained by shifting  η\eta_{-}  up at distance (b)(\ell-b). Similarly, let  χ^:BC\widehat{\chi}:B^{\prime}\to C^{\prime}  be the lattice path obtained by shifting  η^\widehat{\eta}_{-}  up at distance (b)(\ell-b).

Note that the path  ξAC\xi_{AC}  starts below  χ^\widehat{\chi}  and ends above  χ^\widehat{\chi}. Thus  ξAC\xi_{AC}  intersects  χ^\widehat{\chi}  at least once, where the intersection points could be multiple and include AA. Order these points of intersection according to the order in which they appear on ξAC\hskip 0.85355pt\xi_{AC}, and denote by EE the last such point of intersection. Finally, denote by  ξEC\xi_{EC}  the last part of the path  ξAC\xi_{AC}  between EE and CC, and note that  ξEC\xi_{EC}  lies above the xx-monotone lattice path χ^\widehat{\chi}.

(2)(2)  Denote by  η^+\widehat{\eta}_{+}  the lattice path  ACA\to C, starting at AA, following α\alpha up, then  η+\eta_{+} and ending by following β\beta down to CC. Let  ζBC:BC\zeta_{B^{\prime}C^{\prime}}:B^{\prime}\to C^{\prime}  be the lattice path obtained by shifting  ξBD\xi_{BD}  up at distance b\ell-b. By the same argument as above, path  ζBC\zeta_{B^{\prime}C^{\prime}}  start above and ends below  η^+\widehat{\eta}_{+}. Denote by FF the last point of intersection of  η^+\widehat{\eta}_{+}  and  ζBC\zeta_{B^{\prime}C^{\prime}}  according to the order on which they appear on ζBC\zeta_{B^{\prime}C^{\prime}} . Finally, denote by  ζFC{\zeta}_{FC^{\prime}}  the last part of the path  ζBC\zeta_{B^{\prime}C^{\prime}}  between FF and CC, and note that  ζFC\zeta_{FC^{\prime}}  lies below the xx-monotone lattice path η+\eta_{+} .

(3)  Observe that  ζBC\zeta_{B^{\prime}C^{\prime}}  lies above  χ^\widehat{\chi}  since shifting both paths down gives  ξBD\xi_{BD}  lying above  η\eta_{-}, respectively. Also, the path  ξEC\xi_{EC}  lies below  η+\eta_{+}  and above  χ^\widehat{\chi}  by definition. Since CC^{\prime} is below CC, we have that EE and CC lie on different sides of ζFC\zeta_{FC^{\prime}}. Thus lattice paths  ξEC\xi_{EC}  and  ζFC\zeta_{FC^{\prime}}  must intersect in the connected component Λ\Lambda of the region between  η+\eta_{+}  and χ^\widehat{\chi} that contains interval [CC]β[CC^{\prime}]\subset\beta. There could be many such intersections, of course, including multiple intersections when the paths form loops.

Lemma 3 (Fomin [Fom, Thm 6.1]).

Let  γ1:EC\gamma_{1}:E\to C  and  γ2:FC\gamma_{2}:F\to C^{\prime}  be two intersecting paths between boundary points of the simply connected region  ΛΓ\Lambda\subseteq\Gamma. Let  Gγ1γ2G\in\gamma_{1}\cap\gamma_{2}  be an intersection point, and suppose

γ1:=Eγ1Gγ1′′Candγ2:=Fγ2Gγ2′′C,\gamma_{1}\hskip 1.70709pt:=\hskip 1.70709ptE\hskip 1.70709pt\to_{\gamma_{1}^{\prime}}\hskip 1.70709ptG\hskip 1.70709pt\to_{\gamma_{1}^{\prime\prime}}\hskip 1.70709ptC\qquad\text{and}\qquad\gamma_{2}\hskip 1.70709pt:=\hskip 1.70709ptF\hskip 1.70709pt\to_{\gamma_{2}^{\prime}}\hskip 1.70709ptG\hskip 1.70709pt\to_{\gamma_{2}^{\prime\prime}}\hskip 1.70709ptC\hskip 0.85355pt,

by which we mean that  GG  separates  γi\gamma_{i}  into two paths:  γi\gamma_{i}^{\prime}  and  γi′′\gamma_{i}^{\prime\prime}, where  i{1,2}i\in\{1,2\}. Then there is a well defined key intersection point  GG  as above, such that the map  {(γ1,γ2)}{(γ1,γ2)}\bigl{\{}(\gamma_{1},\gamma_{2})\bigr{\}}\to\bigl{\{}(\gamma_{1}^{\ast},\gamma_{2}^{\ast})\bigr{\}}  is an injection, where  the pair of paths  (γ1,γ2)(\gamma_{1}^{\ast},\gamma_{2}^{\ast})  is obtained from  (γ1,γ2)(\gamma_{1},\gamma_{2})  by a swap at GG :

γ1:=Eγ1Gγ2′′Candγ2:=Fγ2Gγ1′′C.\gamma_{1}^{\ast}\hskip 1.70709pt:=\hskip 1.70709ptE\hskip 1.70709pt\to_{\gamma_{1}^{\prime}}\hskip 1.70709ptG\hskip 1.70709pt\to_{\gamma_{2}^{\prime\prime}}\hskip 1.70709ptC^{\prime}\qquad\text{and}\qquad\gamma_{2}^{\ast}\hskip 1.70709pt:=\hskip 1.70709ptF\hskip 1.70709pt\to_{\gamma_{2}^{\prime}}\hskip 1.70709ptG\hskip 1.70709pt\to_{\gamma_{1}^{\prime\prime}}\hskip 1.70709ptC\hskip 0.85355pt.

The lemma is a special case of the (much more general) result by Fomin; below we include a proof sketch for completeness. Denote by GG the  key intersection  of paths  ξEC\xi_{EC}  and  ζFC\zeta_{FC^{\prime}}  defined by the lemma. Finally, denote by  ξGC\xi_{GC}  the last part of the path  ξEC\xi_{EC}  between GG and CC, and by  ζGC{\zeta}_{GC^{\prime}}  the last part of the path  ζFC\zeta_{FC^{\prime}}  between GG and CC.

Refer to caption
Figure 2. Construction of paths  ζGC\zeta_{GC^{\prime}}  and  μGD\mu_{G^{\prime}D^{\prime}}  in the proof of the lemma.

(4)  Denote by  ξAG:AG\xi_{AG}:A\to G  the lattice path  ξAC\xi_{AC}  without the last part  ξGC\xi_{GC}. In the notation of the lemma, define  ξAC\xi_{AC^{\prime}}  to be the path  ξAG\xi_{AG}  followed by  ζGC{\zeta}_{GC^{\prime}}.

(5)  Let GG^{\prime} be the point on ξBD\xi_{BD}  obtained by shifting GG down at distance \ell, and denote by  ξGD\xi_{G^{\prime}D}  the last part of the path  ξBD\xi_{BD}  between GG and DD. Similarly, let  μGD\mu_{G^{\prime}D^{\prime}}  be the path obtained shifting down at distance \ell the path  ζGC{\zeta}_{GC^{\prime}}\hskip 0.85355pt. Denote by  ξBG:BG\xi_{BG^{\prime}}:B\to G^{\prime}  the lattice path  ξBD\xi_{BD}  without the last part  ξGD\xi_{G^{\prime}D}. In the notation of the lemma, define  ξBD\xi_{BD^{\prime}}  to be the path  ξBG\xi_{BG^{\prime}}  followed by  μGD\mu_{G^{\prime}D^{\prime}}.

Claim:The map  Φ:(ξAC,ξBD)(ξAC,ξBD)\Phi:\hskip 0.85355pt\big{(}\xi_{AC},\xi_{BD}\big{)}\hskip 0.85355pt\to\hskip 0.85355pt\big{(}\xi_{AC^{\prime}},\xi_{BD^{\prime}}\big{)}  constructed above is an injection.

To prove the claim, we consider the inverse of Φ1\Phi^{-1}. Start with a pair of lattice paths  (ξAC,ξBD)(\xi_{AC^{\prime}},\xi_{BD^{\prime}})  and follow the steps as above after relabeling  CCC\leftrightarrow C^{\prime}, DDD\leftrightarrow D^{\prime}. This construction will not work at all cases as the shifted paths are no longer guaranteed to intersect because of topological considerations. However, when  (ξAC,ξBD)=Φ(ξAC,ξBD)(\xi_{AC^{\prime}},\xi_{BD^{\prime}})=\Phi(\xi_{AC},\xi_{BD}), the construction will work for the same reason and produce the pair of lattice paths  (ξAC,ξBD)(\xi_{AC},\xi_{BD})  as in the steps (0)(0)(5)(5). Here we are using the key intersection point in step (3)(3) to ensure the construction is well defined and can be inverted at this step. We see that Φ\Phi is an involution between pairs of lattice paths {(ξAC,ξBD)}\bigl{\{}(\xi_{AC},\xi_{BD})\bigr{\}} and the pairs of paths {(ξAC,ξBD)}\bigl{\{}(\xi_{AC^{\prime}},\xi_{BD^{\prime}})\bigr{\}}, which intersect as in (3)(3) after the translations in (1)(1) and (2)(2). The details are straightforward.

Finally, the projection conditions on Φ\Phi as in the lemma are straightforward. Indeed, we effectively swap parts of lattice paths:  ξGC\xi_{GC}  with  ζGC{\zeta}_{GC^{\prime}}, and  ξGD\xi_{G^{\prime}D}  with  μGD\mu_{G^{\prime}D^{\prime}}. Since path  ξGC\xi_{GC}  is shifted path  μGD\mu_{G^{\prime}D^{\prime}}, and path  ξGD\xi_{G^{\prime}D}  is shifted path  ζGC{\zeta}_{GC^{\prime}}, this implies both conditions. ∎

Proof of Theorem 1.

In the notation of Lemma 2, set  a=b=0a=b=0, so both points  A=BA=B  lie at the origin. Further, set  r=1r=1, c=k+1c=k+1 and d=k1d=k-1, so  C=(m,k+1)C=(m,k+1),  C=D=(m,k)C^{\prime}=D^{\prime}=(m,k)  and  D=(m,k1)D=(m,k-1). In this case, the injection  Φ\Phi  shows that the number of pairs of lattice paths  O(m,k+1)O\to(m,k+1)  and  O(m,k1)O\to(m,k-1), is less of equal than the number of lattice paths  O(m,k)O\to(m,k), squared.

We are not done, however, as our paths overcount the paths implied by the probabilities in the theorem, since we consider only the first time by the lattice random walk XtX_{t} is at β\beta. Recall that Φ\Phi preserves the number of horizontal edges which project onto point mm and onto [m1,m][m-1,m]. The assumption in the theorem that TT is the first time the walk is at β\beta can be translated as having exactly one of point (m,)(m,\ast) and exactly one edge  (m1,)(m,)(m-1,\ast)\to(m,\ast)  corresponding to the last step of the walk XtX_{t}. Therefore, this property is also preserved under Φ\Phi. This completes the proof. ∎

Remark 4.

The level of generality in Lemma 2 may seem like an overkill as we only use a special case in the proof of the theorem. In reality, explaining the special case needed for Theorem 1 is no easier than the general case in the lemma. In fact, setting  A=BA=B  only makes it more difficult to keep track of the paths. Furthermore, other properties in the lemma are used heavily in the next section.

Sketch of proof of Lemma 3.

Let  Ω2\Omega\subset\mathbb{Z}^{2}  be a simply connected region and let  γ:PQ\gamma:P\to Q be a lattice walk in Ω\Omega, where  P,QΩP,Q\in\Omega. Define a loop-erased walkLE(γ){\text{\rm LE}}(\gamma)  by removing cycles as they appear in γ\gamma. Formally, take the first self-intersection point XX, where  γ:PγXγ′′Xγ′′′Q\gamma:P\to_{\gamma^{\prime}}X\to_{\gamma^{\prime\prime}}X\to_{\gamma^{\prime\prime\prime}}Q, and remove part γ′′\gamma^{\prime\prime}. Repeat this until the resulting path  LE(γ){\text{\rm LE}}(\gamma)  has no self-intersections.

Suppose  P1,Q1,Q2,P2ΩP_{1},Q_{1},Q_{2},P_{2}\in\partial\Omega  are points on the boundary (oriented clockwise) and in this order. For a pair of paths  (γ1,γ2)(\gamma_{1},\gamma_{2}), γ1:P1Q2\gamma_{1}:P_{1}\to Q_{2}, γ2:P2Q1\gamma_{2}:P_{2}\to Q_{1}, note that  γ1\gamma_{1}  and  γ2\gamma_{2}  intersect by planarity, and so do  LE(γ1){\text{\rm LE}}(\gamma_{1})  and  γ2\gamma_{2}. Let  XX  be the intersection point of  γ2\gamma_{2}  and  LE(γ1){\text{\rm LE}}(\gamma_{1})  which lies closest to P1P_{1} along the path  LE(γ1){\text{\rm LE}}(\gamma_{1}). Note that point XX can appear multiple times on  γ1LE(γ1)\gamma_{1}\supseteq{\text{\rm LE}}(\gamma_{1}).

Consider the point  Xγ1X\in\gamma_{1}  such that the edge  YXY\to X  in γ1\gamma_{1} is not deleted in LE(γ1){\text{\rm LE}}(\gamma_{1}). Similarly, choose the first XX on γ2\gamma_{2}. This defines the key intersection of paths  γ1\gamma_{1}  and  γ2\gamma_{2}. As in the lemma, swap the future of these paths to obtain paths  (γ1,γ2)(\gamma_{1}^{\ast},\gamma_{2}^{\ast}). To see that the map  (γ1,γ2)(γ1,γ2)(\gamma_{1},\gamma_{2})\hskip 0.85355pt\to\hskip 0.85355pt(\gamma_{1}^{\ast},\gamma_{2}^{\ast})  in an injection, note that it is invertible on all pairs  (γ1,γ2)(\gamma_{1}^{\ast},\gamma_{2}^{\ast}) such that  LE(γ1){\text{\rm LE}}(\gamma_{1}^{\ast})  intersects γ2\gamma_{2}^{\ast}. We omit the details. ∎

Remark 5.

There is a much larger probabilistic context in which the loop-erased random walk plays a prominent role, see e.g. [LL, §\S11] and §\S5.3.


3. Large example and subtleties in the proof

The construction in the proof above may seem excessively complicated at first, given that the map  Φ\Phi  is easy to define in the example in Figure 2. Indeed, in that case one can simply shift  ξBD\xi_{BD}  up  to define path  ζBD\zeta_{B^{\prime}D^{\prime}}, find the last (only in this case) intersection point GG with  ξAC\xi_{AC}, and swap the futures of these paths as we do in (4)(4) and (5)(5). Voila!

Unfortunately, this simplistic approach does not work for multiple reasons. Let’s count them here:

(i)  Path  ζGC\zeta_{GC^{\prime}}  does not have to be inside Γ\Gamma. This is why we defined point FF in (2)(2).

(ii)  Path  ζGD\zeta_{G^{\prime}D^{\prime}}  does not have to be inside Γ\Gamma. This is why we defined path χ\chi and point EE in (1)(1).

(iii)  Paths  χ\chi  and  ξAC\xi_{AC}  do not have to intersect at all. This is why we defined path   χ^\widehat{\chi}  in (1)(1).

(iv)  Paths  ζBC\zeta_{B^{\prime}C^{\prime}}  and  ξAC\xi_{AC}  do have to intersect for geometric reason, since  (ab)(cd)r(a-b)\leq(c-d)-r. On the other hand, paths  ξEC\xi_{EC}  and  ζFC\zeta_{FC^{\prime}}  intersect for topological reasons. This is why in (3)(3) we consider a connected component of Λ\Lambda between paths  η+\eta_{+}  and  χ^\widehat{\chi}. Note that the latter can in fact intersect, possibly multiple times.

(v)  Paths  ζBC\zeta_{B^{\prime}C^{\prime}} and  ξAC\xi_{AC}  can have multiple loops intersecting each other in multiple ways. There is no natural way to define the “last intersection” which would be easy to reverse. This is why in (3)(3) we invoked Lemma 3, whose proof requires loop-erased walk and symmetry breaking.111In the first draft of the paper we were not aware of the issue  (v), only to discover it when lecturing on the result.

To help the reader understand the issues  (i),  (ii) and  (iv), consider a large example in Figure 3. Here paths  ξAC\xi_{AC}  and  χ\chi  intersect multiple times in step (1)(1) defining EE. Then, in step (2)(2), path  ζBC\zeta_{B^{\prime}C^{\prime}}  intersects the boundary  η+\eta_{+}  multiple times. Note that FF is defined as the last intersection along path  ζBC\zeta_{B^{\prime}C^{\prime}}, not along  η+\eta_{+} . The same property applies to EE, even if in the example it is the last intersection on both paths.

What exactly goes wrong in  (v)  in the definition of GG is rather subtle and we leave this as an exercise to the reader. Note that there is no issue similar to  (v)  in the definition of points EE and FF, since the boundaries  η±\eta_{\pm}  are xx-monotone.

Remark 6.

Note that we explicitly use xx-monotonicity of the boundary paths  η±\eta_{\pm} . In §\S4.5 below, we address what happens when the boundary is not xx-monotone.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 3. Steps of the construction of injection Φ\Phi in the proof of the lemma.


4. Generalizations and applications

4.1. General transition probabilities

In the notation of the introduction, consider a more general random walk  XtX_{t}  which moves to neighbors with general (not necessarily uniform) transition probabilities:

P[(i,j)(i±1,j)]=π±(i,j),P[(i,j)(i,j±1)]=ω±(i,j),{\text{\bf P}}\bigl{[}(i,j)\to(i\pm 1,j)\bigr{]}\hskip 1.70709pt=\hskip 1.70709pt\pi_{\pm}(i,j),\qquad{\text{\bf P}}\bigl{[}(i,j)\to(i,j\pm 1)\bigr{]}\hskip 1.70709pt=\hskip 1.70709pt\omega_{\pm}(i,j),

with obvious constraints  π±(i,j),ω±(i,j)0\pi_{\pm}(i,j),\hskip 0.85355pt\omega_{\pm}(i,j)\geq 0, and

π+(i,j)+π(i,j)+ω+(i,j)+ω(i,j)=1,\pi_{+}(i,j)\hskip 1.70709pt+\hskip 1.70709pt\pi_{-}(i,j)\hskip 1.70709pt+\hskip 1.70709pt\omega_{+}(i,j)\hskip 1.70709pt+\hskip 1.70709pt\omega_{-}(i,j)\,=1,

for all  (i,j)Γ(i,j)\in\Gamma. We say these transition probabilities are yy-invariant if they are translation invariant with respect to vertical shifts:

π±(i,j)=π±(i,j),ω±(i,j)=ω±(i,j)for all i,j and j.\pi_{\pm}(i,j)\hskip 1.70709pt=\pi_{\pm}(i,j^{\prime}),\quad\omega_{\pm}(i,j)\hskip 1.70709pt=\hskip 1.70709pt\omega_{\pm}(i,j^{\prime})\quad\text{for all \, $i,\hskip 0.85355ptj$ \hskip 1.70709ptand \hskip 1.70709pt$j^{\prime}$.}
Theorem 7.

Let  Γ2\Gamma\subset\mathbb{Z}^{2}  be the lattice region as in Theorem 1. Let  {Xt}\{X_{t}\}  be the lattice random walk which starts at the origin  X0=OαX_{0}=O\in\alpha, moves according to yy-invariant transition probabilities  π±(i,j)\pi_{\pm}(i,j),  ω±(i,j)\omega_{\pm}(i,j)  as above, and is absorbed whenever  XtX_{t}  tries to exit the region Γ\hskip 0.85355pt\Gamma. Denote by  TT  the first time tt such that  XtβX_{t}\in\beta, and let  p(k)p(k)  be the probability that XT=(m,k)X_{T}=(m,k). Then  {p(k)}\{p(k)\}  is log-concave:

p(k)2p(k+1)p(k1)for all k,  such that (m,k±1)β.p(k)^{2}\,\geq\,p(k+1)\hskip 1.70709ptp(k-1)\qquad\text{for all \ $k\in\mathbb{Z}$, \ such that \ \ $(m,k\pm 1)\hskip 0.85355pt\in\beta$}.

The proof of the theorem follows verbatim the proof of Theorem 1 in the previous section, since yy-invariance is exactly the property preserved by the injection Φ\Phi in Lemma 2.  \square

4.2. Monotone walks

Let  {Xt}\{X_{t}\}  be a random walk as above with yy-invariant transition probabilities. We say that the walk is monotone if  π(i,j)=ω(i,j)=0\pi_{-}(i,j)=\omega_{-}(i,j)=0  for all  (i,j)Γ(i,j)\in\Gamma.

Corollary 8.

Let  Γ2\Gamma\subset\mathbb{Z}^{2}  be the lattice region as in Theorem 1, and let  γ={x=}Γ\gamma=\{x=\ell\}\cap\Gamma  be a vertical interval,  0<<m0<\ell<m. Fix points  A=(0,a)αA=(0,a)\in\alpha  and  B=(m,b)βB=(m,b)\in\beta. Let  {Xt}\{X_{t}\}  be a monotone lattice random walk which starts at point  X0=AαX_{0}=A\in\alpha, moves according to yy-invariant transition probabilities  π+(i,j)\pi_{+}(i,j),  ω+(i,j)\omega_{+}(i,j)  as above, is absorbed whenever  XtX_{t}  tries to exit the region Γ\hskip 0.85355pt\Gamma, and arrives at BB at time  u=(m+ba)u=(m+b-a):  Xu=BX_{u}=B. Denote by  TT  the first time tt such that  XtγX_{t}\in\gamma, and let  q(k)q(k)  be the probability that XT=(,k)X_{T}=(\ell,k). Then  {q(k)}\{q(k)\}  is log-concave:

q(k)2q(k+1)q(k1)for all k,  such that (,k±1)γ.q(k)^{2}\,\geq\,q(k+1)\hskip 1.70709ptq(k-1)\qquad\text{for all \ $k\in\mathbb{Z}$, \ such that \ \ $(\ell,k\pm 1)\hskip 0.85355pt\in\gamma$}.
Proof.

Define two regions:  Γ1:={(i,j)Γ:0i}\Gamma_{1}:=\{(i,j)\in\Gamma~{}:~{}0\leq i\leq\ell\}  and  Γ2:={(i,j)Γ:1im}\Gamma_{2}:=\{(i,j)\in\Gamma~{}:~{}\ell-1\leq i\leq m\}. Note that the regions are overlapping along interval γ\gamma and interval  γ:={(1,j)Γ}\gamma^{\prime}:=\{(\ell-1,j)\in\Gamma\}, see Figure 4. Suppose  (1,k)(,k)(\ell-1,k)\to(\ell,k)  is the unique edge of the lattice path  ABA\to B  which projects onto  [1,][\ell-1,\ell]. In the notation of Theorem 7, we have

q(k)=p1(k)p2(k),q(k)\hskip 1.70709pt=\hskip 1.70709ptp_{1}(k)\hskip 1.70709ptp_{2}(k),

where  p1(k)p_{1}(k)  and  p2(k)p_{2}(k)  are the exit probabilities in the region Γ1\Gamma_{1} and in the region Γ2\Gamma_{2} rotated 180180^{\circ}. Since  {pi(k)}\{p_{i}(k)\}  are log-concave by Theorem 7, so is  {q(k)}\{q(k)\}, as desired. ∎

Refer to caption
Figure 4. Left:  a monotone walk XtX_{t} crossing the vertical γ\gamma and γ\gamma^{\prime} (red lines) at height kk. Middle:  steps in the triangular lattice. Right:  the square–octagon lattice.

4.3. General lattices

One can further generalize random walks from 2\mathbb{Z}^{2} to general lattices. For example, we can include steps  ±(1,1)\pm(1,1)  and yy-invariant transition probabilities

P[(i,j)(i±1,j±1)]=ν±(i,j),{\text{\bf P}}\bigl{[}(i,j)\to(i\pm 1,j\pm 1)\bigr{]}\hskip 1.70709pt=\hskip 1.70709pt\nu_{\pm}(i,j),

such that  ν±(i,j)=ν±(i,j)\nu_{\pm}(i,j)\hskip 0.85355pt=\hskip 0.85355pt\nu_{\pm}(i,j^{\prime}), for all i,ji,j and jj^{\prime}. One can view this result as a random walk on the triangular lattice instead, see Figure 4. Both the statement and the proof of Theorem 7 extend verbatim once the reader observes that all lattice paths which must intersect for topological reasons do in fact intersect at lattice points.

Similarly, one can use this approach and general transition probabilities to set some of them zero and obtain random walks on other lattices. For example, one can obtain the square–octagon lattice as in the figure by restricting the walks to vertices of the lattice. Theorem 7 applies to this case then. We omit the details.

4.4. Dyck and Schröder paths

In the context of enumerative combinatorics, it is natural to consider lattice paths with steps  (1,1)(1,1)  and  (1,1)(1,-1). Such paths are called Dyck paths. When step  (2,0)(2,0)  is added, such paths are called Schröder paths.

Fix two points  A=(0,0)A=(0,0)  and B=(m,b)2B=(m,b)\in\mathbb{Z}^{2}  and two nonintersecting Dyck paths  η±:AB\eta_{\pm}:A\to B. Note that  m+b0(mod2)m+b\equiv 0\pmod{2}, since otherwise there are no such paths. Denote by Γ\Gamma the region between these paths. Let 0<<m0<\ell<m, and denote by  N(k)\textrm{N}(k)  the number of Dyck paths   ζ:AB\zeta:A\to B  which lie inside Γ\Gamma and contain point (,k)(\ell,k).

Corollary 9.

In the notation above, the sequence  {N(k)}\{\textrm{\em N}(k)\}  is log-concave:

N(k)2N(k+2)N(k2)for all k,  such that (,k±2)γ.\textrm{\em N}(k)^{2}\,\geq\,\textrm{\em N}(k+2)\hskip 1.70709pt\textrm{\em N}(k-2)\qquad\text{for all \ $k\in\mathbb{Z}$, \ such that \ \ $(\ell,k\pm 2)\hskip 0.85355pt\in\gamma$}.
Proof sketch.

The proof follows verbatim the proof of Corollary 8 via two observations. First, the vertical translation and topological properties used in the proof of Lemma 2 work with diagonal steps. Second, the intersection points of the paths are at the ends of the steps, not midway, because the Dyck paths here have endpoints on the underlying grid spanned by (1,1)(1,1) and (1,1)(1,-1) which is invariant under the (0,2)(0,2) translation. ∎

Similarly, fix two nonintersecting Schröder paths  η±:AB\eta_{\pm}:A\to B, and denote by Γ\Gamma the region between these paths. Let 0<<m0<\ell<m, and denote by  F(k)\textrm{F}(k)  the number of Schröder paths   ζ:AB\zeta:A\to B  which lie inside Γ\Gamma and contain point (,k)(\ell,k). The same argument as above gives the following.

Corollary 10.

In the notation above, the sequence  {F(k)}\{\textrm{\em F}(k)\}  is log-concave:

F(k)2F(k+2)F(k2)for all k,  such that (,k±2)γ.\textrm{\em F}(k)^{2}\,\geq\,\textrm{\em F}(k+2)\hskip 1.70709pt\textrm{\em F}(k-2)\qquad\text{for all \ $k\in\mathbb{Z}$, \ such that \ \ $(\ell,k\pm 2)\hskip 0.85355pt\in\gamma$}.
Remark 11.

Note that this result does not directly apply to the, otherwise similar, Motzkin paths, with steps (1,1)(1,1), (1,1)(1,-1) and (1,0)(1,0). The reason is that the intersection points used in the injection Φ\Phi in Lemma 2 might no longer be on the underlying lattice points and appear in the middle of the steps, e.g. at (12,12)(\frac{1}{2},\frac{1}{2}).

Example 12.

Take  A=(0,0)A=(0,0),  B=(2n,0)B=(2n,0), so  m=2nm=2n. Fix maximal and minimal Dyck paths

η+:\displaystyle\eta_{+}: (0,0)(1,1)(n,n)(n+1,n1)(2n,0),\displaystyle(0,0)\hskip 1.70709pt\to\hskip 1.70709pt(1,1)\hskip 1.70709pt\to\hskip 1.70709pt\ldots\hskip 1.70709pt\to\hskip 1.70709pt(n,n)\hskip 1.70709pt\to\hskip 1.70709pt(n+1,n-1)\hskip 1.70709pt\to\hskip 1.70709pt\ldots\hskip 1.70709pt\to\hskip 1.70709pt(2n,0),
η:\displaystyle\eta_{-}: (0,0)(1,1)(n,n)(n+1,n+1)(2n,0).\displaystyle(0,0)\hskip 1.70709pt\to\hskip 1.70709pt(1,-1)\hskip 1.70709pt\to\hskip 1.70709pt\ldots\hskip 1.70709pt\to\hskip 1.70709pt(n,-n)\hskip 1.70709pt\to\hskip 1.70709pt(n+1,-n+1)\hskip 1.70709pt\to\hskip 1.70709pt\ldots\hskip 1.70709pt\to\hskip 1.70709pt(2n,0).

Set  :=n\ell:=n, which makes the picture symmetric. Then Corollary 9 implies log-concavity of binomial coefficients{(nk),0kn}\bigl{\{}\binom{n}{k},\hskip 0.85355pt0\leq k\leq n\}. On the other hand, for the Schröder paths ABA\to B, Corollary 10 implies log-concavity of Delannoy numbers{D(k,nk),0kn}\bigl{\{}{\text{\rm D}}(k,n-k),\hskip 0.85355pt0\leq k\leq n\}, see [OEIS, A008288], a new enumerative result, see §\S5.5.

Finally, let  η+\eta_{+}  be as above and let

η:(0,0)(1,1)(2,0)(3,1)(2n,0).\eta_{-}:\hskip 1.70709pt(0,0)\hskip 1.70709pt\to\hskip 1.70709pt(1,1)\hskip 1.70709pt\to\hskip 1.70709pt(2,0)\hskip 1.70709pt\to\hskip 1.70709pt(3,1)\hskip 1.70709pt\to\hskip 1.70709pt\ldots\hskip 1.70709pt\to\hskip 1.70709pt(2n,0).

Then Corollary 9 implies log-concavity of ballot numbers{B(k,nk),n/2kn}\bigl{\{}{\text{\rm B}}(k,n-k),\hskip 0.85355ptn/2\leq k\leq n\}, where

B(k,nk)=2kn+1k+1(nk).{\text{\rm B}}(k,n-k)\,=\,\hskip 1.70709pt\frac{2k-n+1}{k+1}\hskip 0.85355pt\binom{n}{k}.

4.5. Boundary matters

First, let us note that both Theorem 1 and Theorem 7 easily extend to regions without either or both of the boundaries  η±\eta_{\pm} . In this case the vertical boundaries α,β\alpha,\beta become either rays of lines, and the region Γ\Gamma is an infinite strip in one or two directions, see Figure 5.

Corollary 13.

In the notation of Theorem 1, let Γ\hskip 0.85355pt\Gamma be an infinite strip with one or two ends. Then the exit probability distribution  {p(k)}\{p(k)\}  is log-concave.

The proof follows immediately from Theorem 1 by taking a sequence  {Γn}\{\Gamma_{n}\}  of regions where the boundary goes to infinity, i.e.  ΓnΓ\Gamma_{n}\to\Gamma, and noting that log-concavity is preserved in the limit. Alternatively, one can easily modify the proof of Lemma 2 to work for unbounded regions; in fact the construction simplifies in that case. We omit the details.

One can also ask whether the xx-monotonicity assumption in the theorem can be dropped. Note that the proof of Lemma 2 breaks in step (2)(2) as the shifted path  ζBC\zeta_{B^{\prime}C^{\prime}}  no longer has to lie inside Γ\Gamma, see Figure 5. Although we do not believe that Theorem 1 extends to non-monotone boundaries, it would be interesting to find a formal counterexample.

Refer to caption
Figure 5. Infinite regions with one and two ends, the issue with non-monotone boundary in step (2)(2) of the proof of the lemma, and a ladder graph GG.
Example 14.

In the notation above, let m=1m=1 and consider an infinite strip between two lines which forms a ladder graph GG as in Figure 5. When restricted to GG, the nearest neighbor random walk moves along α\alpha with equal probability  13\frac{1}{3}  of going up or down, until it eventually moves to the right, at which point it stops. In this case the exit probabilities can be calculated explicitly:

p(±2r)=n=0132n+1(2nnr),p(2r+1)=p(2r1)=n=0132(n+1)(2n+1nr),p(\pm 2r)\,=\,\sum_{n=0}^{\infty}\hskip 1.70709pt\frac{1}{3^{2n+1}}\hskip 1.70709pt\binom{2n}{n-r}\,,\qquad p(2r+1)\,=\,p(-2r-1)\,=\,\sum_{n=0}^{\infty}\hskip 1.70709pt\frac{1}{3^{2(n+1)}}\hskip 1.70709pt\binom{2n+1}{n-r}\,,

for all  r0r\geq 0. A direct calculation gives:

p(±k)=1ϕ2k5for all k0 ,whereϕ=5+12is the golden ratio.p(\pm k)\,=\,\frac{1}{\phi^{2k}\hskip 0.85355pt\sqrt{5}}\quad\text{for all \, $k\geq 0$}\hskip 0.85355pt,\quad\text{where}\ \ \phi\hskip 1.70709pt=\hskip 1.70709pt\frac{\sqrt{5}+1}{2}\ \ \text{is the \emph{\color[rgb]{0.0,0,0.7}\definecolor[named]{pgfstrokecolor}{rgb}{0.0,0,0.7}golden ratio}}.

Thus, log-concavity is an equality at all k0k\neq 0. We leave it as an exercise to the reader to give a direct bijective proof of this fact. Note that these equalities disappear for  m2m\geq 2, cf. §\S5.7.


5. Final remarks

5.1.

Log-concavity of the number of monotone lattice paths as in Corollary 8 is equivalent to the Stanley inequality for posets of width two, as noted in [CFG, GYY]. For general posets, the Stanley inequality was proved in [Sta1]. An explicit injection in the width two case was given in [CFG] and generalized by the authors [CPP1, CPP2]. The construction in [CPP2] was the basis of this paper.

5.2.

Except for the special case of monotone lattice paths and monotone boundary discussed above, we are not aware of the problem even being considered before. The generality of our results is then rather surprising given that even simple special cases appear to be new (see below).

5.3.

The reflection principle is due to Mirimanoff (1923), and is often misattributed to André, see [Ren]. It is described in numerous textbooks, both classical [Fel, Spi] and modern [LL, MP]. For the Karlin–McGregor formula (1959) and its generalizations, including the Brownian motion version of Fomin’s result (Lemma 3), see e.g. [LL, Ch. 9]. For the Lindström–Gessel–Viennot lemma and applications to enumeration of lattice paths, see the original paper [GV] and the extensive treatment in [GJ, §\S5.4]. It is also related to a large body of work on tilings in the context of integrable probability, see [Gor]. For the algebraic combinatorics context of Fomin’s result in connection with total positivity, see [Pos, §\S5].

5.4.

Note also that the log-concavity of exit probabilities does not seem to be a consequence of any standard non-combinatorial approaches. For example, the real-rootedness fails already for Delannoy numbers  {D(k,8k),0k8}\bigl{\{}{\text{\rm D}}(k,8-k),\hskip 0.85355pt0\leq k\leq 8\}, see §\S4.4. We refer to [Bre, Sta2] for surveys of classical methods on unimodality and log-concavity, and to [Pak] for a short popular introduction to combinatorial methods. See also surveys [Brä, Huh] for more recent results and advanced algebraic and analytic tools.

5.5.

In the context of Example 12, Dyck, Schröder and Motzkin paths play a fundamental role in enumerative combinatorics in connection with the Catalan numbers [OEIS, A000108], Schröder numbers [OEIS, A006318] and Motzkin numbers [OEIS, A001006], respectively. Ballot numbers and Delannoy numbers appear in exactly the same context. We refer to [Sta3, Ch. 5] for numerous properties of these numbers.

While binomial coefficients and ballot numbers are trivially log-concave via explicit formulas, the log-concavity of Delannoy numbers appears to be new. Non-real-rootedness in this case suggests that already this special case is rather nontrivial. It would be interesting to see if log-concavity of Delannoy numbers can be established directly, in the style of basic combinatorial proofs in [Sag].

5.6.

There is a large literature on exact and asymptotic counting of various walks in the quarter plane with small steps, see e.g. [Bou, BM]. Most notably, both Kreweras walks (1965) and Gessel walks (2000) fit our framework, while some others do not. It would be interesting to further explore this connection.

5.7.

In the context of §\S4.5 and Example 14, consider a simple random walk constrained to a strip  0xm0\leq x\leq m, reflected at  x=0x=0,  and with no top/bottom boundaries. This special case is especially elegant. The exit probabilities  p(mt)p(mt), as mm\to\infty, converge to the hyperbolic secant distribution, which is log-concave in tt.222See the MathOverflow answer:  https://mathoverflow.net/a/395065/4040. This is in sharp contrast with the case of a simple random walk which starts at the origin, but is not constrained to be in the  x0x\geq 0  halfplane. Denote by  q(k)q(k)  the hitting probabilities of the point (k,m)(k,m) on the line x=mx=m. In this case it is well known that hitting probabilities  q(mt)q(mt), as  mm\to\infty, converge to the Cauchy distribution, see e.g. [Spi, p. 156], which is not log-concave (in tt).


Acknowledgements

We are grateful to Robin Pemantle for interesting discussions on the subject, and to Timothy Budd for his help resolving the  MathOverflow  question. Nikita Gladkov pointed out a subtle issue in step (3)(3) of the construction in the previous version of this paper, necessitating the addition of Lemma 3. We are also thankful to Per Alexandersson and Cyril Banderier for helpful remarks on the paper. We also thank the anonymous referees for useful comments to improve presentation. The last two authors were partially supported by the NSF.


References

  • [Bou] M. Bousquet-Mélou, Counting walks in the quarter plane, in Trends Math., Birkhäuser, Basel, 2002, 49–67.
  • [BM] M. Bousquet-Mélou and M. Mishna, Walks with small steps in the quarter plane, in Contemp. Math. 520, AMS, Providence, RI, 2010, 1–39.
  • [Brä] P. Brändén, Unimodality, log-concavity, real-rootedness and beyond, in Handbook of enumerative combinatorics, CRC Press, Boca Raton, FL, 2015, 437–483.
  • [Bre] F. Brenti, Unimodal, log-concave and Pólya frequency sequences in combinatorics, Mem. AMS 81 (1989), no. 413, 106 pp.
  • [CPP1] S. H. Chan, I. Pak and G. Panova, The cross–product conjecture for width two posets, preprint (2021), 30 pp;  arXiv: 2104.09009.
  • [CPP2] S. H. Chan, I. Pak and G. Panova, Extensions of the Kahn–Saks inequality for posets of width two, preprint (2021), 24 pp.;  arXiv:2106.07133.
  • [CFG] F. R. K. Chung, P. C. Fishburn and R. L. Graham, On unimodality for linear extensions of partial orders, SIAM J. Algebraic Discrete Methods 1 (1980), 405–410.
  • [Fel] W. Feller, An introduction to probability theory and its applications, vol. I (Third ed.), John Wiley, New York, 1968, 509 pp.
  • [Fom] S. Fomin, Loop-erased walks and total positivity, Trans. AMS 353 (2001), 3563–3583.
  • [GV] I. M. Gessel and X. Viennot, Binomial determinants, paths, and hook length formulae, Adv. Math. 58 (1985), 300–321.
  • [Gor] V. Gorin, Lectures on random tilings, monograph draft, Nov. 25, 2019, 191 pp.; https://tinyurl.com/w22x6qq .
  • [GJ] I. P. Goulden and D. M. Jackson, Combinatorial enumeration, Wiley, New York, 1983, 569 pp.
  • [GYY] R. L. Graham, A. C. Yao and F. F. Yao, Some monotonicity properties of partial orders, SIAM J. Algebraic Discrete Methods 1 (1980), 251–258.
  • [Huh] J. Huh, Combinatorial applications of the Hodge–Riemann relations, in Proc. ICM Rio de Janeiro, Vol. IV, World Sci., Hackensack, NJ, 2018, 3093–3111.
  • [LL] G. F. Lawler and V. Limic, Random walk: a modern introduction, Cambridge Univ. Press, Cambridge, UK, 2010, 364 pp.
  • [MP] P. Mörters and Y. Peres, Brownian motion, Cambridge Univ. Press, Cambridge, UK, 2010, 403 pp.
  • [Pak] I. Pak, Combinatorial inequalities, Notices AMS 66 (2019), 1109–1112; an expanded version of the paper is available at  https://tinyurl.com/py8sv5v6.
  • [Pos] A. Postnikov, Total positivity, Grassmannians, and networks, preprint (2006), 79 pp.;  arXiv:math/0609764.
  • [Ren] M. Renault, Lost (and found) in translation, Amer. Math. Monthly 115 (2008), 358–363.
  • [Sag] B. E. Sagan, Inductive and injective proofs of log concavity results, Discrete Math. 68 (1988), 281–292.
  • [OEIS] N. J. A. Sloane, The online encyclopedia of integer sequences, oeis.org.
  • [Spi] F. Spitzer, Principles of random walk (Second ed.), Springer, New York, 1976, 408 pp.
  • [Sta1] R. P. Stanley, Two combinatorial applications of the Aleksandrov--Fenchel inequalities, J. Combin. Theory, Ser. A 31 (1981), 56--65.
  • [Sta2] R. P. Stanley, Log-concave and unimodal sequences in algebra, combinatorics, and geometry, in Graph theory and its applications, New York Acad. Sci., New York, 1989, 500--535.
  • [Sta3] R. P. Stanley, Enumerative Combinatorics, vol. 1 (Second ed.) and vol. 2, Cambridge Univ. Press, 2012 and 1999.