---
title: Planar Corner Growth Model
url: https://www.emergentmind.com/topics/planar-corner-growth-model
type: topic
---

# Planar Corner Growth Model

The planar corner growth model is a directed growth process on the square lattice, usually identified with directed last-passage percolation on \(\mathbb Z^2\). Growth proceeds from a corner of the first quadrant under the rule that a site can be occupied only after its west and south predecessors are available, and the resulting random interface is encoded by maximization over up-right lattice paths with steps \(e_1=(1,0)\) and \(e_2=(0,1)\). In the literature this model appears in both general i.i.d. and exactly solvable exponential or geometric forms, and it serves as a standard \(1+1\)-dimensional KPZ model for the study of deterministic limit shapes, Busemann functions, stationary cocycles, semi-infinite geodesics, and competition interfaces [1510.00859][1709.05771].

## 1. Lattice definition and equivalent formulations

A basic formulation assigns a random weight \(\omega_x\) to each vertex \(x\in\mathbb Z^2\). For coordinatewise ordered points \(x\le y\), the point-to-point last-passage time \(G_{x,y}\) is the maximum total weight collected along an up-right path from \(x\) to \(y\). In one common convention, the terminal weight is excluded from the sum, while the growth cluster is then written as
\[
C_t=\{x\in\mathbb Z_+^2: G_{0,x}+\omega_x\le t\}.
\]
A second convention includes the terminal weight and leads to the recursion
\[
G(i,j)=G(i-1,j)\vee G(i,j-1)+W(i,j),\qquad G(i,0)=G(0,j)=0,
\]
with \(W(i,j)\) the site weight. These are equivalent descriptions of the same directed growth mechanism [1510.00859][1612.08521].

The growth interpretation is that \((i,j)\) becomes infected, red, or occupied only after both \((i-1,j)\) and \((i,j-1)\), together with its own waiting time. The induced time-\(t\) cluster can be written as
\[
Y(t)=\{(i,j)\in\mathbb N^2:G(i,j)\le t\},
\]
so the planar corner growth model is simultaneously a random growth process, a directed polymer optimization problem at zero temperature, and a last-passage percolation model [1612.08521][1910.13378].

Exactly solvable instances arise when the site weights are i.i.d. \(\mathrm{Exp}(1)\) or geometric. More general treatments assume only i.i.d., nondegenerate weights with finite \(p>2\) moment and a lower bound almost surely. In the continuous-weight case, finite maximizing paths are almost surely unique, and the collection of geodesics from the origin forms a geodesic tree [1510.00860][1510.00859].

A recurrent source of confusion is the coexistence of several indexing conventions. Some papers work on \(\mathbb Z_+^2\), others on \(\mathbb N^2\), and some parameterize directions on the simplex \(\mathcal U=\{te_1+(1-t)e_2:0\le t\le 1\}\). These are not different models; they are coordinate and normalization choices adapted to particular questions about shape, stationarity, or geodesics [1510.00859][1812.02689].

## 2. Deterministic asymptotics and shape functions

Under standard ergodicity and moment assumptions, the model has a deterministic limit shape
\[
g(\xi)=\lim_{n\to\infty}\frac1n G_{0,\lfloor n\xi\rfloor},
\]
and, in general i.i.d. settings, \(g\) is continuous, concave, symmetric, and \(1\)-homogeneous. The point-to-line counterpart also has a deterministic limit and is related to \(g\) by convex duality:
\[
g_{\mathrm{pl}}(h)=\sup_{\xi\in\mathcal U}\{g(\xi)+h\cdot\xi\},\qquad
g(\xi)=\inf_{h\in\mathbb R^2}\{g_{\mathrm{pl}}(h)-h\cdot\xi\}.
\]
This dual structure is central to stationary formulations and to the cocycle approach [1510.00859][1404.7786].

In the exponential corner growth model the shape function is explicit:
\[
g(\xi)=\bigl(\sqrt{\xi_1}+\sqrt{\xi_2}\bigr)^2.
\]
The corresponding point-to-line shape is
\[
g_{\mathrm{pl}}(h)=1+\frac{h_1+h_2}{2}+\frac12\sqrt{(h_1-h_2)^2+4}.
\]
These formulas make the exponential model the standard benchmark for exact KPZ analysis [1709.05771].

A key refinement is the stationary exponential model indexed by \(0<\rho<1\). In that construction the bulk weights are \(\mathrm{Exp}(1)\), the horizontal boundary increments are \(\mathrm{Exp}(1-\rho)\), and the vertical boundary increments are \(\mathrm{Exp}(\rho)\). The associated stationary last-passage field has linear shape
\[
g^\rho(s,t)=\frac{s}{1-\rho}+\frac{t}{\rho}.
\]
This linear stationary shape is not the nonstationary limit shape itself; rather, it provides the stationary comparison object from which the nonstationary shape and Busemann structure are extracted [1709.05771].

The exact solvable inhomogeneous models show that strict curvature is not universal. When site weights are exponential with row and column parameters \(a_i+b_j\), or geometric with parameters \(a_i b_j\), the shape is given by explicit variational formulas and can be strictly concave only inside a cone, while becoming linear outside it. This contrasts with the classical expectation for light-tailed i.i.d. models, where the shape is expected to be strictly concave in all directions [1612.08521].

## 3. Stationary cocycles and Busemann functions

A stationary \(L^1\) cocycle is a shift-covariant additive field
\[
B:\Omega\times\mathbb Z^2\times\mathbb Z^2\to\mathbb R
\]
that satisfies stationarity and the cocycle property
\[
B(x,z)=B(x,y)+B(y,z).
\]
Its decisive feature is the recovery identity
\[
\omega_x=\min\{B(x,x+e_1),\,B(x,x+e_2)\},
\]
which makes \(B\) a stationary boundary-gradient object compatible with the original environment. In general i.i.d. planar LPP, such cocycles were constructed from tandem queueing fixed points and used as boundary data for stationary versions of last-passage percolation [1510.00859][1404.7786].

These cocycles solve the variational formulas for the limit shapes. In particular, the point-to-line and point-to-point asymptotics can be written as infima over centered cocycles or cocycles with mean tilt, and the Busemann-type cocycles act as minimizers or correctors. This places the planar corner growth model in a framework analogous to homogenization, with cocycles furnishing stationary gradients of the asymptotic growth profile [1510.00859][1709.05771].

A Busemann function is an almost sure limit of passage-time differences to remote targets:
\[
B^\xi(x,y)=\lim_{n\to\infty}\bigl(G_{x,v_n}-G_{y,v_n}\bigr),
\]
for sequences \(v_n\) tending to infinity in direction \(\xi\). In the general i.i.d. model, existence is proved at exposed differentiability points of the limit shape and on interiors of maximal linear segments with differentiable endpoints. In the exponential model, Busemann functions exist in every direction \(\xi(\alpha)\), and the nearest-neighbor increments have explicit distributions
\[
B^\alpha_{x,x+e_1}\sim \mathrm{Exp}(1-\alpha),\qquad
B^\alpha_{x,x+e_2}\sim \mathrm{Exp}(\alpha),
\]
with mean equal to the shape gradient [1510.00859][1709.05771].

The exactly solvable exponential model admits a stronger all-directions description. The joint distribution of \(\{B^\rho\}_{\rho>1}\) is characterized through multiclass queueing fixed points; for fixed \(x\), the process \(\rho\mapsto B^\rho_{x-e_1,x}\) is càdlàg and equal in distribution to a pure-jump marked point process whose jump times on \((1,\infty)\) form a Poisson point process with intensity \(s^{-1}\,ds\). At each vertex there is a unique random threshold \(\rho^*(x)\) where the preferred predecessor switches, and
\[
\mathbb P\{\rho^*(x)\le \lambda\}=1-\lambda^{-1},\qquad \lambda\ge 1.
\]
These formulas make the directional Busemann field an explicitly tractable stochastic process, not merely a collection of separate one-direction limits [1808.09069].

A further structural sharpening is the strong existence and uniqueness theorem for the tilt-indexed Busemann process. In the i.i.d. planar corner growth model, the Busemann process indexed by tilts in the super-differential of the limit shape can be realized as a measurable function of the underlying i.i.d. weights on the canonical probability space, and any two such realizations coincide almost surely. This removes the need for auxiliary extensions of the probability space and shows that the full tilt-indexed Busemann family is already encoded in the weight field itself [2507.08757].

## 4. Geodesics and the competition interface

Finite maximizing paths are geodesics. Under continuous weights, every vertex in \(\mathbb Z_+^2\) has a unique geodesic from the origin, and these form a geodesic tree. Semi-infinite geodesics are infinite up-right paths whose every finite segment is maximizing. In the general i.i.d. model, stationary cocycles and Busemann functions yield existence of directional semi-infinite geodesics from every starting point, as well as leftmost and rightmost versions when continuity fails [1510.00860].

When the limit shape is differentiable at the relevant directions and the weight distribution is continuous, the directional theory becomes sharp: for each direction there is a unique semi-infinite geodesic from each site, and all such geodesics coalesce. Finite geodesics to targets \(v_n\to\infty\) converge locally to the corresponding Busemann geodesic. In the exponential model this program can be carried out by a comparatively soft argument based on the increment-stationary growth process, coupling, and planar monotonicity, rather than on detailed fluctuation estimates [1510.00860][1812.02689].

The competition interface is the up-right dual path separating the descendants of \(e_1\) and \(e_2\) in the geodesic tree rooted at the origin. In the continuous-weight setting it satisfies a law of large numbers,
\[
\frac{\varphi_n}{n}\to \xi_*,
\]
where \(\xi_*\) is a random direction characterized by the balancing of cocycle increments. The atoms of the law of \(\xi_*\) are exactly the corners of the limit shape, and under differentiability assumptions on the endpoints of linear segments, \(\xi_*\) is almost surely an exposed point of the shape [1510.00860].

The nonexistence of bi-infinite geodesics is one of the main global rigidity results. For a fixed direction, the increment-stationary approach yields absence of bi-infinite directed geodesics in the exponential model. A stronger theorem shows that, with i.i.d. \(\mathrm{Exp}(1)\) weights, there are almost surely no nontrivial bi-infinite geodesics at all; the only bi-infinite geodesics are the trivial axis lines. The proof combines stationary couplings, coarse graining, exit-point estimates, and planarity [1812.02689][1909.06883].

## 5. Fluctuations, coalescence scales, and KPZ geometry

In the stationary exponential corner growth model, the longitudinal fluctuation exponent is \(1/3\) and the transversal exponent is \(2/3\). The variance of \(G^\rho_{0,(m,n)}\) under characteristic scaling is of order \(N^{2/3}\), corresponding to \(N^{1/3}\)-order fluctuations of passage times, while the geodesic deviates from its characteristic line on the \(N^{2/3}\) scale. These are the standard KPZ exponents in \(1+1\) dimensions, proved in this setting by coupling and variance identities rather than solely by asymptotic Fredholm determinant analysis [1709.05771].

In exactly solvable inhomogeneous exponential and geometric models, fluctuations remain KPZ inside the strictly concave cone. For slopes \(r\in(c_1,c_2)\), the rescaled last-passage time converges to the Tracy-Widom GUE law, and the upper tail above the shape satisfies
\[
\mathbb P_{\mathbf a,\mathbf b}\!\left( G(\lfloor nr\rfloor,n)\ge n\gamma_n + ns \right) \le C e^{-c n(s^{3/2}\wedge s)}.
\]
Thus disorder in row and column parameters can alter the macroscopic geometry by creating linear regions, while preserving Tracy-Widom fluctuations in the curved regime [1612.08521].

Coalescence times of semi-infinite geodesics exhibit the KPZ-predicted \(3/2\) exponent. In the exponential model, when two \(\xi[\rho]\)-directed geodesics start \(k\) apart, the natural coalescence distance is of order \(k^{3/2}\). This is quantified by four tail estimates. For small \(\delta\), the probability of slow coalescence outside a box is bounded by
\[
C_1\delta \le
\mathbb P\Big\{{}^{\rho}\!\big(\lfloor\delta N^{2/3}e_1\rfloor,\lfloor\delta N^{2/3}e_2\rfloor\big)\notin 0,v_N\Big\}
\le C_2 |\log\delta|^{2/3}\delta,
\]
and for large \(r\), the fast-coalescence probability inside the box satisfies
\[
e^{-C_1 r^3}\le
\mathbb P\Big\{{}^{\rho}\!\big(\lfloor rN^{2/3}e_1\rfloor,\lfloor rN^{2/3}e_2\rfloor\big)\in 0,v_N\Big\}
\le e^{-C_2 r^3}.
\]
These estimates are derived from stationary LPP, Busemann functions, exit-point control, and Pimentel’s duality [1911.03792].

The resulting picture is that the same directional Busemann field governs both geodesic geometry and fluctuation theory. This suggests a unified structure: shape curvature determines the macroscopic direction, stationary increments produce Busemann limits, and the KPZ exponents dictate the scales on which geodesics separate and merge.

## 6. Inhomogeneous, combinatorial, and higher-dimensional extensions

A major extension replaces i.i.d. weights by row- and column-dependent laws. In one exactly solvable family, the exponential model uses rates \(a_i+b_j\) and the geometric model uses parameters \(a_i b_j\), with \((a_i)\) and \((b_j)\) drawn from ergodic distributions. The shape functions are then
\[
g(s,t)=\inf_{z\in(-\underline\alpha,\underline\beta)}
\left\{
s\,\mathbb E\!\left[\frac{1}{a+z}\right]
+t\,\mathbb E\!\left[\frac{1}{b-z}\right]
\right\}
\]
in the exponential case, and an analogous formula in the geometric case. These models retain enough solvability for full limit-shape theory and Tracy-Widom asymptotics in the curved cone, while allowing linear sectors absent from the homogeneous exponential benchmark [1612.08521].

A related inhomogeneous exponential theory treats additively separable finite-array rates \(\omega^{m,n}(i,j)\sim \mathrm{Exp}(a_m(i)+b_n(j))\). Its first-order asymptotics are given by the explicit centering
\[
\mathcal G_{m,n}^{\mathrm{cent}}
=
\inf_{z\in(-\mathfrak a_m,\mathfrak b_n)}
\left\{
\sum_{i=1}^m \frac{1}{a_m(i)+z}
+
\sum_{j=1}^n \frac{1}{b_n(j)-z}
\right\},
\]
and under vague convergence of the parameter arrays the rescaled cluster converges to a deterministic limit shape. The boundary of that shape can contain flat segments adjacent to the axes, while the cluster can develop visible axis spikes and persistent crevices. A vertical spike occurs iff \(\mathfrak a_s>\mathfrak a\), a horizontal spike iff \(\mathfrak b_s>\mathfrak b\), and flat segments adjacent to the axes are characterized by the finiteness of
\[
\int_{\mathbb R}\frac{\alpha(da)}{(a-\mathfrak a)^2},
\qquad
\int_{\mathbb R}\frac{\beta(db)}{(b-\mathfrak b)^2}.
\]
The same variational structure yields macroscopic flux and particle-profile formulas for step-initial TASEP with particle and hole disorder [1908.09319].

The geometric corner growth model also has a precise combinatorial incarnation in boxed plane partitions and lozenge tilings. For i.i.d. geometric weights
\[
\mathbb P(w_{ij}=k)=(1-q)q^k,
\]
the joint law of a column of last-passage times is identified with a \(g\)-measure on partitions:
\[
\mathbb P\bigl(G(b,c)=\lambda_1,\dots,G(1,c)=\lambda_b\bigr)=P_{g,b,c}(\lambda).
\]
Under the plane-partition and lozenge-tiling bijections, corner counts have generating functions given by normalized Schur polynomials, and asymptotic consequences include limit shapes, Tracy-Widom GUE fluctuations, and GUE-minors limits in fixed-width regimes [1910.13378].

The term “corner growth model” is also used in higher-dimensional settings. In a three-dimensional corner growth model formulated as an SOS process inside a corner, a projection onto the \((1,1,1)\) plane yields an RSOS model on a triangular lattice with periodic boundary conditions, enabling a high-precision velocity estimate
\[
w=0.12606(2),
\]
in disagreement with the conjectured \(w=1/8\). This higher-dimensional result is not a planar corner growth theorem, but it highlights how dimension changes both the natural stationary representation and the asymptotic-shape problem [1212.0946].

Source: https://www.emergentmind.com/topics/planar-corner-growth-model