Connective constant

In mathematics, the connective constant is a numerical quantity associated with self-avoiding walks on a lattice. It is studied in connection with the notion of universality in two-dimensional statistical physics models. While the connective constant depends on the choice of lattice so itself is not universal (similarly to other lattice-dependent quantities such as the critical probability threshold for percolation), it is nonetheless an important quantity that appears in conjectures for universal laws. Furthermore, the mathematical techniques used to understand the connective constant, for example in the recent rigorous proof by Duminil-Copin and Smirnov that the connective constant of the hexagonal lattice has the precise value [math]\displaystyle{ \sqrt{2+\sqrt{2}} }[/math], may provide clues to a possible approach for attacking other important open problems in the study of self-avoiding walks, notably the conjecture that self-avoiding walks converge in the scaling limit to the Schramm–Loewner evolution.

Definition

The connective constant is defined as follows. Let [math]\displaystyle{ c_n }[/math] denote the number of n-step self-avoiding walks starting from a fixed origin point in the lattice. Since every n + m step self avoiding walk can be decomposed into an n-step self-avoiding walk and an m-step self-avoiding walk, it follows that [math]\displaystyle{ c_{n+m} \leq c_n c_m }[/math]. Then by applying Fekete's lemma to the logarithm of the above relation, the limit [math]\displaystyle{ \mu=\lim_{n \rightarrow \infty} c_n^{1/n} }[/math] can be shown to exist. This number [math]\displaystyle{ \mu }[/math] is called the connective constant, and clearly depends on the particular lattice chosen for the walk since [math]\displaystyle{ c_n }[/math] does. The value of [math]\displaystyle{ \mu }[/math] is precisely known only for two lattices, see below. For other lattices, [math]\displaystyle{ \mu }[/math] has only been approximated numerically. It is conjectured that [math]\displaystyle{ c_n \approx \mu^n n^{\gamma -1} }[/math] as n goes to infinity, where [math]\displaystyle{ \mu }[/math] and [math]\displaystyle{ A }[/math], the critical amplitude, depend on the lattice, and the exponent [math]\displaystyle{ \gamma }[/math], which is believed to be universal and dependent on the dimension of the lattice, is conjectured to be [math]\displaystyle{ \gamma=43/32 }[/math].

Known values

Lattice	Connective constant
Hexagonal	[math]\displaystyle{ \sqrt{2 + \sqrt{2}}\simeq 1.85 }[/math]
Triangular	[math]\displaystyle{ 4.15079(4) }[/math]
Square	[math]\displaystyle{ 2.63815853032790(3) }[/math]
Kagomé	[math]\displaystyle{ 2.56062 }[/math]
Manhattan	[math]\displaystyle{ 1.733535(3) }[/math]
L-lattice	[math]\displaystyle{ 1.5657(15) }[/math]
[math]\displaystyle{ (3.12^2) }[/math] lattice	[math]\displaystyle{ 1.7110412... }[/math]
[math]\displaystyle{ (4.8^2) }[/math] lattice	[math]\displaystyle{ 1.80883001(6) }[/math]

These values are taken from the 1998 Jensen–Guttmann paper and a more recent paper by Jacobsen, Scullard and Guttmann. The connective constant of the [math]\displaystyle{ (3.12^2) }[/math] lattice, since each step on the hexagonal lattice corresponds to either two or three steps in it, can be expressed exactly as the largest real root of the polynomial

[math]\displaystyle{ x^{12} - 4x^8 - 8x^7 - 4x^6 + 2x^4 + 8x^3 + 12x^2 + 8x + 2 }[/math]

given the exact expression for the hexagonal lattice connective constant. More information about these lattices can be found in the percolation threshold article.

Duminil-Copin–Smirnov proof

In 2010, Hugo Duminil-Copin and Stanislav Smirnov published the first rigorous proof of the fact that [math]\displaystyle{ \mu=\sqrt{2 + \sqrt{2}} }[/math] for the hexagonal lattice. This had been conjectured by Nienhuis in 1982 as part of a larger study of O(n) models using renormalization techniques. The rigorous proof of this fact came from a program of applying tools from complex analysis to discrete probabilistic models that has also produced impressive results about the Ising model among others. The argument relies on the existence of a parafermionic observable that satisfies half of the discrete Cauchy–Riemann equations for the hexagonal lattice. We modify slightly the definition of a self-avoiding walk by having it start and end on mid-edges between vertices. Let H be the set of all mid-edges of the hexagonal lattice. For a self-avoiding walk [math]\displaystyle{ \gamma }[/math] between two mid-edges [math]\displaystyle{ a }[/math] and [math]\displaystyle{ b }[/math], we define [math]\displaystyle{ \ell(\gamma) }[/math] to be the number of vertices visited and its winding [math]\displaystyle{ W_{\gamma}(a,b) }[/math] as the total rotation of the direction in radians when [math]\displaystyle{ \gamma }[/math] is traversed from [math]\displaystyle{ a }[/math] to [math]\displaystyle{ b }[/math]. The aim of the proof is to show that the partition function

[math]\displaystyle{ Z(x)=\sum_{\gamma: a\to H}x^{\ell(\gamma)}=\sum_{n=0}^{\infty}c_n x^n }[/math]

converges for [math]\displaystyle{ x\lt x_c }[/math] and diverges for [math]\displaystyle{ x\gt x_c }[/math] where the critical parameter is given by [math]\displaystyle{ x_c=1/ \sqrt{2+\sqrt{2}} }[/math]. This immediately implies that [math]\displaystyle{ \mu=\sqrt{2+\sqrt{2}} }[/math].

Given a domain [math]\displaystyle{ \Omega }[/math] in the hexagonal lattice, a starting mid-edge [math]\displaystyle{ a }[/math], and two parameters [math]\displaystyle{ x }[/math] and [math]\displaystyle{ \sigma }[/math], we define the parafermionic observable

[math]\displaystyle{ F(z)=\sum_{\gamma\subset\Omega:a\to z} e^{-i\sigma W_{\gamma}(a,z)}x^{\ell(\gamma)}. }[/math]

If [math]\displaystyle{ x=x_c=1/\sqrt{2 + \sqrt{2}} }[/math] and [math]\displaystyle{ \sigma=5/8 }[/math], then for any vertex [math]\displaystyle{ v }[/math] in [math]\displaystyle{ \Omega }[/math], we have

[math]\displaystyle{ (p-v)F(p) + (q-v)F(q) + (r-v)F(r)=0, }[/math]

where [math]\displaystyle{ p,q,r }[/math] are the mid-edges emanating from [math]\displaystyle{ v }[/math]. This lemma establishes that the parafermionic observable is divergence-free. It has not been shown to be curl-free, but this would solve several open problems (see conjectures). The proof of this lemma is a clever computation that relies heavily on the geometry of the hexagonal lattice.

Next, we focus on a finite trapezoidal domain [math]\displaystyle{ S_{T,L} }[/math] with 2L cells forming the left hand side, T cells across, and upper and lower sides at an angle of [math]\displaystyle{ \pm \pi/3 }[/math]. (Picture needed.) We embed the hexagonal lattice in the complex plane so that the edge lengths are 1 and the mid-edge in the center of the left hand side is positioned at −1/2. Then the vertices in [math]\displaystyle{ S_{T,L} }[/math] are given by

[math]\displaystyle{ V(S_{T,L})=\{ z\in V(\mathbb{H}) : 0 \leq Re(z)\leq \frac{3T+1}{2}, \; |\sqrt{3}Im(z)-Re(z)|\leq 3L\}. }[/math]

We now define partition functions for self-avoiding walks starting at [math]\displaystyle{ a }[/math] and ending on different parts of the boundary. Let [math]\displaystyle{ \alpha }[/math] denote the left hand boundary, [math]\displaystyle{ \beta }[/math] the right hand boundary, [math]\displaystyle{ \epsilon }[/math] the upper boundary, and [math]\displaystyle{ \bar{\epsilon} }[/math] the lower boundary. Let

[math]\displaystyle{ A_{T,L}^x:=\sum_{\gamma \in S_{T,L}:a\to \alpha\setminus\{a\}} x^{\ell(\gamma)},\quad B_{T,L}^x:=\sum_{\gamma \in S_{T,L}:a\to \beta} x^{\ell(\gamma)}, \quad E_{T,L}^x:=\sum_{\gamma \in S_{T,L}:a\to \epsilon \cup \bar{\epsilon}} x^{\ell(\gamma)}. }[/math]

By summing the identity

[math]\displaystyle{ (p-v)F(p) + (q-v)F(q) + (r-v)F(r)=0 }[/math]

over all vertices in [math]\displaystyle{ V(S_{T,L}) }[/math] and noting that the winding is fixed depending on which part of the boundary the path terminates at, we can arrive at the relation

[math]\displaystyle{ 1=\cos(3\pi/8) A_{T,L}^{x_c} + B_{T,L}^{x_c} + \cos(\pi/4) E_{T,L}^{x_c} }[/math]

after another clever computation. Letting [math]\displaystyle{ L\to\infty }[/math], we get a strip domain [math]\displaystyle{ S_T }[/math] and partition functions

[math]\displaystyle{ A_{T}^x:=\sum_{\gamma \in S_{T}:a\to \alpha\setminus\{a\}} x^{\ell(\gamma)},\quad B_{T}^x:=\sum_{\gamma \in S_{T}:a\to \beta} x^{\ell(\gamma)}, \quad E_{T}^x:=\sum_{\gamma \in S_{T}:a\to \epsilon \cup \bar{\epsilon}} x^{\ell(\gamma)}. }[/math]

It was later shown that [math]\displaystyle{ E_{T,L}^{x_c}=0 }[/math], but we do not need this for the proof. We are left with the relation

[math]\displaystyle{ 1=\cos(3\pi/8) A_{T,L}^{x_c} + B_{T,L}^{x_c} }[/math].

From here, we can derive the inequality

[math]\displaystyle{ A_{T+1}^{x_c} - A_{T}^{x_c} \leq x_c (B_{T+1}^{x_c})^2 }[/math]

And arrive by induction at a strictly positive lower bound for [math]\displaystyle{ B_{T}^{x_c} }[/math]. Since [math]\displaystyle{ Z(x_c)\geq\sum_{T\gt 0}B_T^{x_c}=\infty }[/math], we have established that [math]\displaystyle{ \mu\geq\sqrt{2+\sqrt{2}} }[/math].

For the reverse inequality, for an arbitrary self avoiding walk on the honeycomb lattice, we perform a canonical decomposition due to Hammersley and Welsh of the walk into bridges of widths [math]\displaystyle{ T_{-I}\lt \cdots \lt T_{-1} }[/math] and [math]\displaystyle{ T_0\gt \cdots \gt T_j }[/math]. Note that we can bound

[math]\displaystyle{ B_T^x\leq (x/x_c)^T B_T^{x_c}\leq (x/x_c)^T }[/math]

which implies [math]\displaystyle{ \prod_{T\gt 0}(1+B_T^x)\lt \infty }[/math]. Finally, it is possible to bound the partition function by the bridge partition functions

[math]\displaystyle{ Z(x)\leq \sum_{T_{-I} \lt \cdots \lt T_{-1},\; T_0\gt \cdots \gt T_j} 2 \left(\prod_{k=-I}^j B_{T_k}^x\right)=2\left(\prod_{T\gt 0}(1+B_T^x)\right)^2\lt \infty. }[/math]

And so, we have that [math]\displaystyle{ \mu=\sqrt{2+\sqrt{2}} }[/math] as desired.

Conjectures

Nienhuis argued in favor of Flory's prediction that the mean squared displacement of the self-avoiding random walk [math]\displaystyle{ \langle |\gamma(n)|^2 \rangle }[/math] satisfies the scaling relation [math]\displaystyle{ \langle |\gamma(n)|^2 \rangle=\frac{1}{c_n} \sum_{n\;\mathrm{step\; SAW}}|\gamma(n)|^2=n^{2\nu +o(1)} }[/math], with [math]\displaystyle{ \nu=3/4 }[/math]. The scaling exponent [math]\displaystyle{ \nu }[/math] and the universal constant [math]\displaystyle{ 11/32 }[/math] could be computed if the self-avoiding walk possesses a conformally invariant scaling limit, conjectured to be a Schramm–Loewner evolution with [math]\displaystyle{ \kappa=8/3 }[/math].

See also

• Percolation threshold

External links

• Weisstein, Eric W.. "Self-Avoiding Walk Connective Constant". http://mathworld.wolfram.com/Self-AvoidingWalkConnectiveConstant.html.