The critical Ising model on trees, concave recursions and nonlinear capacity

Robin Pemantle * * Department of Mathematics, University of Pennsylvania, 209 South 33rd Street, Philadelphia, PA 19104 pemantle@math.upenn.edu Research supported in part by NSF grant # DMS 0103635. and Yuval Peres † † Departments of Statistics and Mathematics, University of California, Berkeley, CA 94720. peres@stat.berkeley.edu Research supported in part by NSF Grants DMS-0104073 and #DMS-0244479.

November 27, 2006

Abstract
We consider the Ising model on a general tree under various boundary conditions: all plus, free and spin-glass. In each case, we determine when the root is influenced by the boundary values in the limit as the boundary recedes to infinity. We obtain exact capacity criteria that govern behavior at critical temperatures. For plus boundary conditions, an L 3   capacity arises. In particular, on a spherically symmetric tree that has n α b n   vertices at level n   (up to bounded factors), we prove that there is a unique Gibbs measure for the ferromagnetic Ising model if and only if α 1 / 2   . Our proofs are based on a new link between nonlinear recursions on trees and L p   capacities.
Keywords: Ising model, reconstruction, capacity, nonlinear potential theory, trees, iteration, spin-glass, recursion. Subject classification: Primary: 60K35; Secondary: 31C45.

1 Introduction

Let T   be a finite rooted tree. Let | v |   denote the distance from a vertex v V ( T )   to the root o   and write v w   if v   is the parent of w   , i.e., the neighbor of w   closer to the root than w   . Consider the space Ω = Ω ( T ) = { + 1 , 1 } V ( T )   of configurations on the vertices of T   . For each w 0   there is a unique edge v w   with v w   ; let J w = J ( v w )   be a positive number, so that { J w : 0 w V ( T ) }   is a fixed set of interaction strengths on the edges of T   . We assume throughout that the interaction strengths are bounded:
0 < J m i n J v J m a x v V ( T ) , v o . (1.1)
This assumption loses little generality; see the end of Section 4. Fix an inverse temperature β   and define the weight of a configuration η Ω   to be the product over all pairs of neighboring vertices:
W ( η ) = v w exp ( β J w η ( v ) η ( w ) ) .   The Ising model under various boundary conditions can be obtained by restricting to suitable subsets of Ω   and assigning probabilities proportional to W   . Our aim in this paper is to pinpoint the locations of the phase transitions that occur in these models as N   . In each case the critical temperature for phase transitions to occur is known. We refine these results by giving sharp criteria for the existence of a phase transition in terms of capacities.

2 Main results

Let T   be any tree, rooted at a vertex o   , and let T   denote the set of maximal paths oriented away from the root; these are either infinite or end at a leaf of T   . For finite trees, we may identify T   with the set of leaves in T   different from o   . Let { R ( e ) : e E ( T ) }   be a set of resistances (nonnegative numbers) assigned to the edges of T   . Let μ   be a flow on T   , that is, a nonnegative function on E ( T )   such that the at every vertex (except for the root and the leaves) inflow equals outflow: whenever v w   ( v   is a parent of w   ) and w   is not a leaf, we have μ ( v w ) = y : w y μ ( w y )   . Denote | μ | : = y : o y μ ( o y )   . Fix p > 1   and set s = p 1   . For y T   define
V μ ( y ) : = e y ( μ ( e ) R ( e ) ) s ; (2.2)
V ( μ ) : = sup { V μ ( y ) : y T } ; (2.3)
c a p p ( T ) : = sup { | μ | : V ( μ ) = 1 } . (2.4)
These capacities have been studied on more general networks as part of discrete nonlinear potential theory; see, e.g., Murakami and Yamasaki (1992), Soardi (1993, 1994) and the references therein. However, all the properties of c a p p   that we will use follow readily from the definition. We note that c a p 2 ( T )   reduces to the electrical conductance between o   and T   .
Let P   denote the probability measure on Ω   proportional to W   :
P ( η ) = W ( η ) ξ Ω W ( ξ ) .   This is a ferromagnetic Ising model with no external field and free boundary conditions.
There is another construction of the measure P   as an error-transmission model. To the edge leading to a vertex v   from its parent, assign the positive bias
θ v = e β J v e β J v e β J v + e β J v = tanh ( β J v ) . (2.5)
Conditional on the sign η ( o )   at the root, let the sign at each other vertex v   be determined recursively, by copying the sign at the parent with probability ( 1 + θ v ) / 2   and reversing sign with probability ( 1 θ v ) / 2   . When J v = J   does not depend on v   , we write θ   for the common bias. Now suppose T   is an infinite, locally finite tree, rooted at o, and let T ( N )   be the induced finite subgraph of T   with vertices { v V ( T ) : | v | N }   . Letting P ( N )   be the free-boundary Ising measure on T ( N )   , we ask about P ( N ) ( η ( 0 ) = + 1 | η ( v ) : v T ( N ) )   .
In particular, this converges in probability to 1 / 2   if and only if the free boundary Gibbs measure on T   is extremal, see Georgii (1988). The question of extremality of the Gibbs measure with free boundary on regular trees was settled by Bleher et al (1995), see also Ioffe (1996a) for an elegant alternative proof.
The same question for general trees was solved by Ioffe (1996b) and Evans et al (2000), where the critical value is computed for an arbitrary tree. However, the question of extremality at the critical temperature was left open. In this paper we settle the critical case by showing that zero L 2   capacity (with respect to certain resistances) implies extremality.
We prove
Theorem 2.1 Let T   be an infinite, locally finite tree, rooted at o, with no leaves except possibly at o and interaction strengths J v   satisfying ( 1.1 ). For vertices y , w   , write y w   if y   is on the path from o   to w   . Assign to each edge e = v w   with v w   , the resistance
R w : = R ( e ) : = o < y w ( tanh β J y ) 2 (2.6)
Then the Gibbs measure for free boundary is extremal if and only if c a p 2 ( T ) = 0   .
One direction of this theorem (that extremality implies zero capacity) was already proved in Evans et al (2000).
Let Ω + = Ω + ( T ) Ω ( T )   be the configurations with η ( v ) = + 1   for v T   . Then the probability measure P +   on Ω +   defined by P + ( η ) = W ( η ) ξ Ω + W ( ξ )   is the Ising model with plus boundary conditions and no external field. The critical value of the interaction strength here has long been known for regular trees (see Preston 1974, 1976). Lyons (1991) computes the critical temperature for general trees and allows the interaction strengths to vary as well. We refine the known results by determining what happens at criticality. The sharp criterion turns out to involve an “ L 3   -capacity”. We prove
Theorem 2.2 Let T   be any infinite, locally finite tree rooted at 0 and having no leaves except possibly 0. Let { J v }   be bounded interaction strengths, i.e., satisfying ( 1.1 ), and assign resistance R v = o < y v ( tanh ( β J y ) ) 1   to the edge between v   and its parent. Then the decreasing limit lim N P ( N , + ) ( η ( 0 ) = + 1 )   is equal to 1 / 2   if and only if c a p 3 ( T ) = 0   .
Here P ( N , + )   is the measure on configurations on the first N   levels of T   with plus boundary conditions imposed at level N   .
For ease of reading, we state the result more explicitly in the special case of spherically symmetric trees, and when the interaction strength is constant.
Corollary 2.3 Under the hypotheses of Theorem  2.2 , assume spherical symmetry as well: θ v = θ | v |   and deg ( v ) = d | v |   depend only on | v |   . Then there are multiple Gibbs states if and only if
n 1 n i = 1 ( d i θ i ) 2 < . (2.7)
In particular, for a spherically symmetric tree T   , suppose that the level cardinalities satisfy
| T n | θ n n α . (2.8)
Then there is a unique Gibbs state for the Ising model at criticality if and only if α 1 / 2   .
Note that for T   satisfying ( 2.8 ), endowed with edge resistances θ n   at level n   , the standard L 2   capacity of T   is zero as long as α 1   .
Corollary 2.4 Suppose that J v J   is constant, and let θ : = tanh ( β J )   . Then phase transition occurs with plus boundary conditions if and only if c a p 3 ( T ) > 0   with resistances θ n   at distance n   from the root. If T   is spherically symmetric, this is equivalent to n 1 θ 2 n | T n | 2 < .  
Define a measure P s g   by making the signs η ( v )   for v T   i.i.d. fair coin flips, and requiring that the measure be proportional to W   conditionally on the values on T   :
P s g ( η ) = 2 | T | W ( η ) ξ | T = η | T W ( ξ ) .   This is equivalent to the following spin-glass model considered by Chayes et al (1986):
the Hamiltonian has no external field and interactions of a fixed magnitude, with the signs of the interactions determined by i.i.d. fair coin flips; the boundary conditions are all plus.
The question is whether, conditional upon the signs of the interactions, the sign at the root is influenced at all by the boundary values in the limit as N   . A critical interaction strength is given in Chayes et al (1986) for regular trees; we improve this to the case of general trees and settle what happens at the critical case. The result is a standard (i.e., L 2   ) capacity criterion, exactly equal to the criterion for the case of a free boundary.
Theorem 2.5 Let T   be an infinite, locally finite tree, rooted at o, with no leaves (except possibly o) and interaction strengths J v   satisfying ( 1.1 ). Assign resistances R v = 0 < y v ( tanh ( β J y ) ) 1   as in ( 2.6 ). Then P s g ( η ( o ) = 1 | η | T ) 1 / 2   in probability under the spin-glass measure if and only if c a p 2 ( T ) = 0   .
Let x v   denote the likelihood ratio of having spin 1 versus 1   at v   , given the boundary.
The method in in the plus boundary case is to show that { x v : v V ( T ) }   satisfy a recursion of the form
x v = v w f w ( x w ) . (2.9)
This reduces the problem to the question of whether, on a given infinite tree, this recursion has a nonzero solution. We give a general solution to this problem, recursively establishing a set of inequalities relating solutions and sub-solutions of these equations to generalized capacities. In the cases of free and spin-glass boundary conditions, the log likelihood ratios are random variables { X v : v V ( T ) }   and we obtain versions of ( 2.9 ) for certain moments { m v }   of { X v }   .
The rest of the paper is organized as follows. The next section focuses entirely on the deterministic aspect of the problem, namely, when the recursion ( 2.9 ) has a nontrivial solution or sub-solution. The theorems in this section are broad enough to handle the recursions arising from the three types of boundary conditions in the Ising model. Then we spend one section on each of the three models and conclude with some questions.

3 Recursions on trees and potential theory

Let T   be any locally finite rooted tree and let { f v : v V ( T ) }   be a collection of nonnegative functions indexed by the vertices of T   . We are interested in whether the simultaneous inequalities
x v v w f w ( x w ) (3.10)
have any nonzero solutions. A special case of interest is when f v f   does not depend on v   . Our characterization is in terms of generalized capacities, which we defined in ( 2.4 ).
Fix p > 1   and let s = p 1   . In addition to ( 2.2 ) and ( 2.4 ), we define
V ( T ) : = inf { V ( μ ) : | μ | = 1 } ; (3.11)
( μ ) : = e E ( T ) μ ( e ) p R ( e ) s . (3.12)
We quote several easy, and well known, consequences of the definitions of energy and capacity.
  • ( i )   The supremum in the definition ( 2.4 ) of c a p p   is achieved if non-empty.
    Furthermore, inf { ( μ ) : | μ | = 1 } = V ( T ) = c a p p ( T ) s   and the infimum is achieved if non-empty.
    ( i i )   Joining several trees at the root sums their capacities.
    ( i i i )   Multiplying all resistances by α   decreases capacity by a factor of α   .
    ( i v )   A single edge of resistance R   in series with a tree T   has capacity c a p p ( T ) 1 + R s ( c a p p ( T ) s ) 1 / s .  
For instance, to see ( i v )   , observe that there is a one to one correspondence between flows μ   from the root to the boundary in T   and flows μ R   in the enhanced tree, such that | μ R | = | μ |   and V ( μ R ) = V ( μ ) + R s | μ | s   . These yield the following lemma, which is all we need below. Recall that T ( v )   denotes the subtree from v   .
Lemma 3.1 Fix p > 1   and s = p 1   . For any vertex v   , define φ ( v ) : = R v c a p p ( T ( v ) ) ,   where R o   =1 by convention. (In particular, φ ( v ) = R v   if v   is a leaf.) Then for any vertex v   , φ ( v ) = v w ( R v / R w ) φ ( w ) ( 1 + φ ( w ) s ) 1 / s .  
Proof: If w o   , let T ( w )   be the tree rooted at the parent of w   consisting of T ( w )   plus the edge between w   and its parent. Then
φ ( v ) = R v c a p p ( T ( v ) ) = v w R v c a p ( T ( w ) )
= v w ( R v / R w ) R w c a p p ( T ( w ) ) ( 1 + R ( w ) s c a p p ( T ( w ) ) s ) 1 / s
which gives the desired expression.   We now relate these computations to the system ( 3.10 ). In the following theorem, f ( )   denotes liminf x f ( x )   and s   denotes p 1   .
Theorem 3.2 Suppose that there exist κ 1 > 0   , p = 1 + s > 1   and a collection of positive constants { a v : v V ( T ) }   such that for every v V ( T )   and x 0   ,
f v ( x ) a v x ( 1 + ( κ 1 x ) s ) 1 / s . (3.13)
Then the solution to the system x v = v w f w ( x w )   with boundary conditions x w =   when w   is a leaf satisfies
x 0 c a p p ( T ) κ 1 , (3.14)
where the resistances are given by R v = 0 y v a y 1   . Similarly, if
a v x ( 1 + ( κ 2 x ) s ) 1 / s f v ( x ) (3.15)
then
c a p p ( T ) κ 2 x 0 . (3.16)
Proof: Let g ( v ) = R v c a p p ( T ( v ) ) / κ 2   , with g ( v ) =   if v 0   is a leaf. We show by induction that g ( v ) x v   for all v   . If v   is a leaf, this is true by definition. Assume v   is not a leaf and, by induction, that g ( w ) x w   for all v w   . Applying the previous lemma gives g ( v ) = v w ( R v / R w ) g ( w ) ( 1 + ( κ 2 g ( w ) ) s ) 1 / s .   By monotonicity of x / ( 1 + ( c x ) s ) 1 / s   , and the induction hypothesis, this is at most v w ( R v / R w ) x w ( 1 + ( κ 2 x w ) s ) 1 / s ,   which is equal to v w a w x w ( 1 + ( κ 2 x w ) s ) 1 / s .   This is at most v w f w ( x w )   by assumption, finishing the induction. An exactly analogous induction establishes that x v G ( v )   for all v   , where G ( v ) = R v c a p p ( T ( v ) ) / κ 1   . Setting v = 0   now recovers the statement of the theorem.   With regard to sub-solutions, that is, to the system of inequalities ( 3.10 ), we have the following immediate corollary, used in Section 4 to analyze Ising models with free boundaries.
Corollary 3.3 Under the hypothesis ( 3.13 ), any solution to x v v w f w ( x w )   satisfies x 0 c a p p ( T ) / κ 1   .  
Although the finitary result in Theorem  3.2 is the most useful, the following corollary for infinite trees is more elegant. The corollary follows directly from the fact that c a p p ( T )   is the decreasing limit of c a p p ( T ( N ) )   , so we omit the details.
Corollary 3.4 Let T   be infinite and locally finite, having no leaves except possibly the root. If c a p p ( T ) = 0   then there are no nontrivial solutions to ( 3.10 ) on T   . If c a p p ( T ) > 0   then any such sub-solution x v v w   satisfies x 0 c a p p ( T ) / κ 1   . Furthermore, there is a sub-solution with the property that x 0 c a p p ( T ) / κ 2   .
To see the value in what we have proved, we turn to some special cases.
Corollary 3.5 Suppose f   satisfies
  • ( i )   f ( x ) = a x Θ ( x p )   near 0 for some p > 1   ; ( i i )   f ( x ) / x   is bounded away from a   except near zero; ( i i i )   f ( x )   is bounded away from zero except near zero.
Then there is a nontrivial sub-solution x v x w f ( x w )   on the vertices of T   if and only if c a p p ( T ) > 0   with resistances a n   at distance n   from the root.  
Remark: This condition is known to be related to many other conditions of interest. For example, it is sufficient that T   have branching number greater than a 1   and necessary that the branching number be at least a 1   . For p = 2   , it is equivalent to the transience of a homesick random walk with parameter a 1   (see Lyons 1992).
Corollary 3.6 Suppose that T   , an infinite, locally finite, leafless tree, is spherically symmetric, meaning that the degree of v   depends only on v   . Suppose f v = f | v |   depends only on | v |   as well. Assume the inequalities ( 3.13 ) and ( 3.15 ). Then there is a nonzero solution to x v = v w f w ( x w )   if and only if n = 1 n j = 1 1 ( a j d j ) s < .  
We use this in the next section with s = 2   to obtain an exact summability criterion for phase transition of the Ising model on an arbitrary spherically symmetric tree. This refinesthe work of Lyons (1989), who computed the critical value but did not settle the behavior at criticality.

4 Plus boundary conditions

In this section T   is an infinite tree with no leaves except possibly the root and T ( N )   denotes the truncation to distance at most N   from the root. We fix interaction strengths { J w : 0 w V ( T ) }   satisfying ( 1.1 ), set θ v = tanh ( β J v ) )   , and consider the family of measures P ( N , + )   of configurations of zeros and ones on T ( N )   with plus boundary conditions.
The goal is to determine whether P ( N , + ) ( η ( 0 ) = + 1 )   converges to 1 / 2   or is bounded below by 1 / 2 + ε   as N   . This is accomplished in the following theorem, already stated in the introduction.
Theorem  2.2  Let T   be any infinite, locally finite tree rooted at 0 and having no leaves except possibly 0. Let { J v }   be bounded interaction strengths, i.e., satisfying ( 1.1 ), and assign resistances R v = 0 < y v ( tanh ( β J y ) ) 1   as in ( 2.6 ). Then the decreasing limit lim N P ( N , + ) ( η ( 0 ) = + 1 )   is equal to 1 / 2   if and only if c a p 3 ( T ) = 0   .
The key to the proof of Theorem  2.2 and to the main results in each of the next two sections is the following recursive likelihood computation. Recall that T ( N )   is the truncation of T   to N   generations. For any tree denote by T ( v )   the subtree rooted at v   , so that for | v | N   , the tree T ( N ) ( v )   has vertex set { w V ( T ) : v w , | w | N }   .
If ξ : T { ± 1 }   is a boundary configuration, let P ξ   denote the Ising measure with boundary condition ξ   and let P v ( N , ξ )   denote the Ising measure on T ( N ) ( v )   whose boundary condition is ξ | T ( N ) ( v )   .
Lemma 4.1 For each v o   let θ v = tanh ( β J v )   . Let x v ( N ) = x v ( N , ξ ) = log [ P v ( N , ξ ) ( η ( v ) = + 1 ) P v ( N , ξ ) ( η ( v ) = 1 ) ]   be the log-likelihood ratio at the root given the boundary. Then for | v | < N   , x v ( N ) = v w f w ( x w ( N ) ) ,   where
f v ( x ) : = log [ cosh ( x / 2 ) + θ v sinh ( x / 2 ) cosh ( x / 2 ) θ v sinh ( x / 2 ) ] . (4.17)
Proof of Lemma  4.1 : Let η   be a configuration on T ( N ) ( v )   . If | v | < N   then for each child w   of v   , let η w   be the restriction of η   to the subtree T ( N ) ( w )   . We may then write W ( η ) = v w W ( η w ) exp ( η ( v ) η ( w ) β J w ) .   Let S   be the set of children of v   . Summing over all vectors { ε w : w S }   of ± 1   's, and within that, over all η   coinciding with ε   on S   , and writing Z   for the normalizing factor (the sum of weights of all configurations), we obtain
P v ( N , ξ ) ( η ( v ) = + 1 ) = Z 1 ε , η v w [ W ( η w ) exp ( η ( w ) β J w ) ]
= Z 1 v w η w : T ( N ) ( w ) { ± 1 } W ( η w ) exp ( η ( w ) β J w )
= Z 1 v w [ e β J w P w ( N , ξ ) ( η ( w ) = 1 ) + e β J w P w ( N , ξ ) ( η ( w ) = 1 ) ] .
Similarly, P v ( N , ξ ) ( η ( v ) = 1 ) = Z 1 v w [ e β J w P w ( N , ξ ) ( η ( w ) = + 1 ) + e β J w P w ( N , ξ ) ( η ( w ) = 1 ) ] .   Divide these expressions, and then divide top and bottom by v w P w ( N , ξ ) ( η ( w ) = 1 ) ,   to obtain
P v ( N , ξ ) ( η ( v ) = + 1 ) P v ( N , ξ ) ( η ( v ) = 1 ) = v w e ( β J w + x w ( N ) ) + e β J w e ( β J w + x w ( N ) ) + e β J w
= v w cosh ( β J w ) ( e x w ( N ) + 1 ) + sinh ( β J w ) ( e x w ( N ) 1 ) cosh ( β J w ) ( e x w ( N ) + 1 ) sinh ( β J w ) ( e x w ( N ) 1 ) .
Next, divide numerator and denominator by cosh ( β J w )   and recall that tanh ( β J w ) = θ   .
It follows that the log of the likelihood ratio above satisfy
x v ( N ) = v w log e x w ( N ) + 1 + θ ( e x w ( N ) 1 ) e x w ( N ) + 1 θ ( e x w ( N ) 1 ) . (4.18)
Finally, divide numerator and denominator by e x w ( N ) / 2   to complete the proof.   Proof of Theorem  2.2 : Specialize to plus boundary conditions. Thus we write P v ( N , + )   for P v ( N , ξ )   where ξ + 1   . Let x v ( N ) = x v ( N , + ) = log ( P v ( N , + ) ( η ( v ) = + 1 ) P v ( N , + ) ( η ( v ) = 1 ) )   be the log-likelihood ratio of plus to minus at the root of the subtree T ( N ) ( v )   . Lemma  4.1 shows that x v ( N ) = v w f w ( x w ( N ) ) ,   with f v   as in equation ( 4.17 ). Write θ   for θ v   and differentiate f v   :
f v ( x ) = θ cosh 2 ( x / 2 ) θ 2 sinh 2 ( x / 2 ) = θ 1 + ( 1 θ 2 ) sinh 2 ( x / 2 ) , x > 0 . (4.19)
It is also clear that the denominator in ( 4.19 ) is increasing in x   and has the expansion 1 + ( 1 θ 2 ) x 2 / 4 + O ( x 3 )   near 0. Thus f v ( x )   is a concave function of x > 0   . Inverting and integrating, we also see that the Taylor expansion of f v   near 0 has the form
f θ ( x ) = f v ( x ) = θ x θ ( 1 θ 2 ) 12 x 3 + O ( x 4 ) . (4.20)
Another consequence of ( 4.19 ) is that f v ( x )   is an odd function of x   . Claim:there are continuous functions κ 2 ( θ ) κ 1 ( θ ) > 0   such that
θ x ( 1 + κ 2 x 2 ) 1 / 2 f v ( x ) θ x ( 1 + κ 1 x 2 ) 1 / 2 . (4.21)
Indeed, if we establish this, then the boundedness assumption ( 1.1 ) implies that constants κ 1   and κ 2   may be chosen so that ( 4.21 ) holds for all v   with θ v   in place of θ   . It then follows from Theorem  3.2 that x 0 ( N )   is bounded between c a p 3 ( T ( N ) ) / κ 2   and c a p 3 ( T ( N ) ) / κ 1   for all N   . Taking decreasing limits finishes the proof of the theorem. It remains to prove the claim.
The choice of κ 2   is easy. When J v   is bounded away from 0 and infinity, θ v   and 1 θ v   are bounded away from 0, so the x 3   term in ( 4.20 ) is bounded away from zero. Thus δ ( θ ) > 0   and κ 2 ( θ )   may be chosen so that ( 4.21 ) holds for x δ   . Since the O ( x 4 )   remainder term varies continuously with p   , the choices of δ   and κ 2   can be made continuously in p   as well.
Now note that sup x ( x / ( 1 + ( κ x ) 2 ) 1 / 2 ) < κ 1   . Let f θ   denote f v   for v   with θ v = θ   . Letting M ( θ ) = f θ ( δ ( θ ) )   and choose κ 2 ( θ ) κ 2 ( θ )   so that M ( θ ) > κ 2 ( θ ) 1   gives a choice of κ 2 ( θ )   varying continuously with p   and satisfying ( 4.21 ).
To choose κ 1   , we may again use the expansion near zero to get a κ 1 * ( θ )   that works for x δ   . We It is shown later that each f θ   is concave (Lemma  6.3 below). Also, f θ   is bounded by log 1 + θ 1 θ   . Thus f θ ( θ ε ) x M   for some ε , M > 0   and all x δ   . On the other hand, θ x / ( 1 + ( κ x ) 2 ) 1 / 2   converges to θ x   uniformly on the interval [ δ , M / θ ]   as κ 0   . Thus κ 1 ( θ ) κ 1 * ( θ )   may be chosen so that θ x / ( 1 + ( κ 1 x ) 2 ) 1 / 2 ( θ ε ) x M   for all x δ   . Since the right-hand side is an upper bound for f θ   on [ δ , )   , and since the choice of κ 1   is evidently continuous in θ   , this finishes the proof of the theorem.   We conclude this section with a discussion of the boundedness condition ( 1.1 ). Given any tree T   with associated interactions { J ( e ) : e E ( T ) }   a new tree T   may be constructed by subdividing edges of T   according to the following scheme. Fix an ε > 0   . Replace each edge e   with θ e < ε   by a series of n   edges e ( 1 ) , , e ( n )   with θ e ( j ) = θ e 1 / n   , where n   is the least integer making θ e 1 / n   greater than ε   .
From the error propagation description of the Ising measure, we see that the measure on { ± 1 } V ( T )   gotten by restricting the Ising measure on T   to the vertices of T   coincides with the Ising measure on T   . Distances in T   no longer coincide with distances in T   , but it is easy to see that the various definitions of phase transition in this article are unchanged if limits on T   are taken with respect to distances in T   . The associated resistor network to T   may be described as follows. Each edge not subdivided retains the same resistance. A subdivided edge with resistance R ( e ) = A / θ e   is replaced by n   edges in series, of resistances A θ e j / n   for j = 1 , , n   . Since θ e 1 / n < ε 1 / 2   , the total resistance is less than 1 / ( 1 ε )   times the greatest resistance, which is A / θ e   . Thus the resistance of the new network is equal to the old resistance up to a bounded factor, and hence has capacity within a bounded factor of the original capacity. We conclude that no generality is lost by assuming J v   to be bounded away from zero.
There is some generality lost in assuming J v   to be bounded above, but for good reason, as shown by the following example. Let T   be a spherically symmetric tree with | T n | n α 2 n   for some α > 1 / 2   . As seen in Corollary  2.4 , there is a phase transition on T   with constant interaction strength satisfying θ = 1 / 2   . Now replace each edge in generation n   by n   edges having θ e = 2 1 / n   . The resistance of each new series of edges in generation n   is of order n   times the old resistance, so when α 3 / 2   , the new tree has zero capacity.Thus the capacity criterion breaks down when the interaction strengths are allowed to have θ v 1   , i.e., J v   .

5 Free boundary conditions

The question we ask in this section is: if you generate a configuration on T ( N )   from the free boundary measure, then look only at the boundary, do you have non-vanishing information about the root as N   ? To formalize this, let ξ   be the random boundary configuration induced by the free measure P ( N )   on configurations on all of T ( N )   . In the notation of Lemma  4.1 , let X v ( N ) : = x v ( N , ξ )   be the log-likelihood ratio of plus to minus at v   given the boundary. It is important to keep track of the spaces on which these are defined, so we re-iterate:
  • x v ( N )   is a map from configurations on T ( N ) ( v )   to real numbers. By natural projection, it may be thought of as a map from configurations on T   to reals.
    X v ( N )   is a random variable on Ω ( T v ( N ) )   , which we may view naturally as a projection of Ω ( T ( N ) )   , or by extension on Ω ( T )   .
    X v ( N )   is measurable with respect to v ( N )   , which we define to be the σ   -field generated by boundary values { η ( y ) : y T v ( N ) }   .
We want to know whether the P ( N )   law of X 0 ( N )   (the free law) converges weakly to a point mass 0 as N   . Evans et al (2000) showed, in the case of constant interaction strength J   , that X 0 ( N )   does not go to zero when T   has positive L 2   capacity with resistances θ 2 n   , where θ = tanh ( J / t )   . As mentioned in the introduction, they (as well as Ioffe (1996) have results in the other direction which leave the critical case open. We sharpen this by showing that zero capacity implies X 0 ( N ) D 0   .
Theorem  2.1  Let T   be an infinite, locally finite tree, rooted at o, with no leaves except possibly at o and interaction strengths J v   satisfying ( 1.1 ) and set θ v = tanh ( β J v ) )   . Supposethat c a p 2 ( T ) = 0   with resistances assigned as in ( 2.6 ) with k = 2   . Then X 0 ( N )   converges in distribution to 0.
Proof of Theorem  2.1 : By Lemma  4.1 , when | v | < N   ,
X v ( N ) = v w f w ( X w ( N ) ) , (5.22)
holds for every sample, with f w   as in equation ( 4.17 ). To make sense of this functional recursion, we will derive from it a system of real inequalities:
m v ( N ) v θ w 2 m w ( N ) 1 + κ m w ( N ) . (5.23)
The quantity m v ( N )   will be an expectation of X v ( N )   but it is not obvious what measure should be used to take the expectation. Define the measures Q v N +   (respectively Q v N   ) on the σ   -field v ( N )   of boundary values by letting Q v N + ( ξ ) : = P v ( N ) ( η : η | T v ( N ) = ξ | η ( v ) = + 1 )   be the conditional distribution of the free boundary given a plus at v   (respectively, given a minus at v   ). Define m v ( N ) : = X v ( N ) d Q v N + = X v ( N ) d Q v N .   The properties of the measures Q v N ±   summarized in the following lemmas make these appropriate for the study of the free boundary.
Lemma 5.1 For any v   with | v | < N   , Q v N + = v w [ ( 1 + θ w ) 2 Q w N + + ( 1 θ w ) 2 Q w N ] .   In particular, the projection of Q v N +   onto boundary configurations on T ( N ) ( w )   is ( 1 + θ w ) 2 Q w N + + ( 1 θ w ) 2 Q w N .    
Lemma 5.2 For any odd function φ   , φ ( X v ( N ) ) d Q v N + = φ ( | X v ( N ) | ) tanh ( | X v ( N ) | / 2 ) d P v ( N ) .  
Lemma 5.3 There is a positive, continuous function κ   such that when f θ   is defined as in ( 4.17 ) with θ = θ v   , then f θ ( X v ( N ) ) d Q v N + θ X v ( N ) d Q v N + 1 + κ ( θ ) X v ( N ) d Q v N + .  
To finish the proof from these lemmas, use ( 5.22 ) and Lemma  5.1 to evaluate
m v ( N ) = v w f w ( X w ( N ) ) d Q v N +
= v w f w ( X w ( N ) ) d ( ( 1 + θ w ) Q w N + + ( 1 θ w ) Q w N )
= v w θ w f w ( X w ( N ) ) d Q w N + . (5.24)
Apply Lemma  5.3 to see that this is at most v w θ w 2 m v ( N ) 1 + κ ( θ v ) m v ( N ) .   By continuity of κ ( θ )   and the boundedness assumption ( 1.1 ), we arrive at ( 5.23 ). Theorem  3.2 now applies to show that m 0 ( N ) c a p 2 ( T ( N ) ) κ   with resistances as in the hypothesis of the theorem. Hence c a p 2 ( T ) = 0   implies m 0 ( N ) 0   as N   . Finally, by Lemma  5.2 , this implies X 0 ( N ) D 0   as N   , finishing the proof.   It remains to prove the lemmas. Lemma  5.1 is immediate from the Markov property.
Proof of Lemma  5.2 : We first compare Q v N +   to the boundary measure induced by the free measure P v ( N )   . We claim that
d Q v N + d P v ( N ) = 1 + tanh ( X v ( N ) / 2 ) . (5.25)
Indeed, from Bayes' rule, one gets d Q v N + d P v ( N ) = P v ( N ) ( η ( v ) = + 1 | v ( N ) ) P v ( N ) ( η ( v ) = + 1 ) .   The denominator is 1 / 2   by symmetry, while the numerator is exp ( X v ( N ) ) / ( 1 + exp ( X v ( N ) ) ) = ( 1 + tanh ( X v ( N ) / 2 ) ) / 2   by definition of X v ( N )   . This proves the claim. Now if φ   is any odd function, then φ ( x ) = ( φ ( x ) φ ( x ) ) / 2   , and thus
φ ( X v ( N ) ) d Q v N + = 1 2 ( φ ( X v ( N ) ) φ ( X v ( N ) ) ) d Q v N +
= ( φ ( X v ( N ) ) φ ( X v ( N ) ) ) e X v ( N ) / 2 e X v ( N ) / 2 + e X v ( N ) / 2 d P v ( N )
= φ ( X v ( N ) ) e X v ( N ) / 2 e X v ( N ) / 2 e X v ( N ) / 2 + e X v ( N ) / 2 d P v ( N ) ,
yielding the desired conclusion.   Proof of Lemma  5.3 : Abbreviate the notation by writing X   for X v ( N )   , E   for integration against P v ( N )   and E +   for integration against Q v N +   . First, for any c > 0   , the product E + f θ ( X ) ( 1 + c E + X ) = E + f θ ( X ) + c ( E + f θ ( X ) ) ( E + X )   is equal, by Lemma  5.2 , to the sum E [ f θ ( | X | ) tanh | X / 2 | ] + E [ f θ ( | X | ) tanh | X / 2 | ) ] E [ c | X | tanh | X / 2 | ] .   Since the functions f θ ( x ) tanh ( x / 2 )   and c x tanh ( x / 2 )   are both nondecreasing on [ 0 , )   , they are positively correlated functions of | X |   (under P v ( N )   or any other law), and hence
( E + f θ ( X ) ) ( 1 + c E + X ) E f θ ( | X | ) tanh | X / 2 | + E ( c | X | f θ ( | X | ) tanh 2 | X / 2 | )
= E [ f θ ( | X | ) tanh | X / 2 | ( 1 + c | X | tanh | X / 2 | ) ] .
Recall that tanh ( x ) = x Θ ( x 3 )   . Refer to the Taylor expansion for f θ = f v   in equation ( 4.20 ) to see that for κ ( θ )   sufficiently small, there is a range x [ 0 , δ ]   for which
f θ ( x ) ( 1 + κ ( θ ) x tanh ( x / 2 ) ) < θ x . (5.26)
Since f θ   is itself bounded and less than θ x ε ( θ ) x   on [ δ , )   , we may choose κ ( θ )   smaller if necessary so that ( 5.26 ) holds for all x 0   . Clearly the choice of κ   can be made continuously in θ   . It follows that ( E + f θ ( X ) ) ( 1 + κ ( θ ) E + X ) E θ | X | tanh | X / 2 | = θ E + X   by Lemma  5.2 . Dividing by ( 1 + κ ( θ ) E + X )   proves the lemma.  

6 Spin-glasses

Let P v ( N , 0 )   denote the spin-glass measure P s g   on configurations on the tree T ( N ) ( v )   (see Section 1 for definitions). Our object in this section is to determine when P 0 ( N , s g ) ( η ( 0 ) = + 1 | ( N ) )   converges in distribution to 0, where ( N ) = 0 ( N )   is the σ   -field generated by boundary values on T ( N )   . By the Markov random field property (or by the definitions of P   and P s g   ), the measures P ( N )   and P ( N , s g )   agree when conditioned on the boundary, so the functions X v ( N )   of the previous section compute conditional probabilities with respect to P ( N , s g )   . Thus our task is to see when X 0 ( N ) D 0   under the laws P ( N , s g )   .
Theorem  2.5  Let T   be an infinite, locally finite tree, rooted at 0, with no leaves except possibly at 0 and interaction strengths J v   satisfying ( 1.1 ) and set θ v = tanh ( β J v )   . Then X 0 ( N ) D 0   under the spin-glass measure if and only if c a p 2 ( T ) = 0   with resistances R v = y v θ y 2   as assigned in ( 2.6 ).
Proof: The structure of the proof is similar to that of Theorem  2.1 . We begin with equation ( 5.22 ):
X v ( N ) = v w f w ( X w ( N ) ) .   This time we have unconditional independence of boundary values instead of conditional independence, so instead of m v ( N )   we work with U v ( N ) : = ( X v ( N ) ) 2   and the mean u v ( N ) : = ( X v ( N ) ) 2 d P v ( N , s g ) ,   where the integrating measure in this case is just i.i.d. fair coin-flips on the boundary of T ( N ) ( v )   . In place of Lemma  5.1 we have the observation that the random variables X w ( N )   have mean zero and are independent as w   ranges over the children of a fixed v   . Lemmas  5.2 and  5.3 are replaced by the following two lemmas. Define g v ( x ) : = ( f v ( x ) ) 2 .  
Lemma 6.1 For all v   and all N > | v |   , E ( U v ( N ) ) 2 3 ( E U v ( N ) ) 2 .  
Lemma 6.2 There are continuous functions κ 2 ( c , θ v ) κ 1 ( c , θ v ) > 0   such that for any random variable V   satisfying E V 2 c ( E V ) 2   , one has
h 2 ( E V ) E g v ( V ) h 1 ( E V ) , (6.27)
with h i ( x ) = θ 2 v x / ( 1 + κ i ( c , θ v ) x )   .
From these two lemmas the proof is finished as follows. Let E   denote expectation with respect to i.i.d. unbiased (spin-glass) boundary conditions. Since each f v   is an odd function, the quantities f ( X w ( N ) )   are independent mean-zero as w   varies over the children of v   , which gives rise to the recursive formula
u v ( N ) = E ( X v ( N ) ) 2
= E ( v w f v ( X w ( N ) ) ) 2
= v w E f v ( X w ( N ) ) 2
= v w E g v ( U w ( N ) ) .
Apply Lemma  6.2 with V = U v ( N )   and c = 3   (obtaining the hypothesis from Lemma  6.1 ), to get v w h 2 ( u w ( N ) ) u v ( N ) v w h 1 ( u w ( N ) ) .   By continuity and the boundedness assumption ( 1.1 ), we may take κ i   in the definition of h i   to be constants independent of v   . By Theorem  3.2 we see that lim N u 0 ( N )   is estimated up to a constant factor by c a p 2 ( T )   with resistances as stated in the hypothesis of the theorem.
Since X 0 ( N )   has mean zero and is bounded by 0 v log [ ( 1 + θ v ) / ( 1 θ v ) ]   , it follows that the random variables X 0 ( N )   converge in distribution to 0 if and only if their variances u 0 ( N )   go to zero. This completes the proof of Theorem  2.5 .   It remains to prove Lemmas  6.1 and  6.2 . Before proving Lemma  6.1 , we record some preliminary facts.
Lemma 6.3 Suppose f   is a differentiable, weakly increasing and concave function on [ 0 , )   , with f ( 0 ) = 0   . Then x 2 f x   is concave.
Proof: Let φ ( x ) = f ( x 0 ) + ( x x 0 ) f ( x 0 )   be the tangent line for f   at x 0   . Concavity implies that φ ( x ) f ( x )   for all x 0   and that h ( x 0 ) f ( x 0 ) / x 0   . Thus φ ( x ) = a x + b   with b 0   , whence x 2 φ x   is a concave support function, lying above x 2 f x   withequality at x 0 2   . We conclude that x 2 f x   is the minimum of a family of concave functions.
 
Lemma 6.4 Let g : [ 0 , ) [ 0 , )   be concave with g ( 0 ) = 0   , and let Y   be a nonnegative random variable with positive finite variance. Then
E [ g ( Y ) 2 ] [ E g ( Y ) ] 2 E Y 2 ( E Y ) 2 . (6.28)
Proof: Let Z = Y / E Y   and h ( z ) = g ( z E Y ) / E ( g ( Y ) )   . Then E Z = E h ( Z ) = 1   , so there must exist z 1 , z 2 > 0   such that h ( z 1 ) z 1   and h ( z 2 ) z 2   . We also may assume that h ( z )   is not identically equal to z   , and thus by concavity there is a unique fixed point x > 0   for which h ( x ) = x   . For any z 0   , | h ( z ) x | | z x | ,   and therefore, E [ h ( Z ) 2 ] = E ( h ( Z ) x ) 2 + 2 x x 2 E ( Z x ) 2 + 2 x x 2 = E Z 2 ,   proving the lemma.  
Lemma 6.5 For any non-negative random variable X L 4   , and any concave function f   with f ( 0 ) = 0   , E f 4 ( X ) ( E f 2 ( X ) ) 2 E X 4 ( E X 2 ) 2 .  
Proof: by Lemma  6.3 , the function g : = x 2 f x   is concave. Applying Lemma  6.4 to the function g   and the random variable Y = X 2 L 2   gives E f 4 ( X ) ( E f 2 ( X ) ) 2 = E g 2 ( X ) [ E g ( X ) ] 2 E X 2 ( E X ) 2 = E X 4 ( E X 2 ) 2 ,   proving the lemma.   Proof of Lemma  6.1 : Recall the definitions of U v ( N )   and u v ( N )   and define the fourth moment s v ( N )   :
  • U v ( N ) = ( X v ( N ) ) 2   ; u v ( N ) = E U v ( N )   ; s v ( N ) = E ( U v ( N ) ) 2 = E ( X v ( N ) ) 4   .
For any v   , the random variables { f w ( X w ( N ) ) : v w }   are independent with mean zero, so any monomial of these will have mean zero unless all exponents are even. The basic recursion ( 5.22 ) yields
u v ( N ) = E ( v w f w ( X w ( N ) ) ) 2
= v w E f w ( X w ( N ) ) 2 .
Hence
( u v ( N ) ) 2 = v w ( E f w ( X w ( N ) ) 2 ) 2 + v { w , w } 2 E f w ( X w ( N ) ) 2 E f w ( X w ( N ) ) 2 . (6.29)
The fourth power expands similarly:
s v ( N ) = E ( v w f w ( X w ( N ) ) ) 4
= v w E f w ( X w ( N ) ) 4 + v { w , w } 6 E f w ( X w ( N ) ) 2 E f w ( X w ( N ) ) 2 . (6.30)
It is required to show that s v ( N ) 3 ( u v ( N ) ) 2   .
Proceed by induction on N | v |   . First suppose N | v | = 1   and that v   has d   children.
Then X v ( N )   is the sum of d   independent mean-zero random variables, each equal to ± log ( p / ( 1 p ) )   . In this case, s v ( N ) / ( u v ( N ) ) 2 = 3 2 / d < 3   . Now suppose N | v | > 1   . By induction, s w 3 u w 2   for each child w   of v   . Applying Lemma  6.5 , we see that for each such w   , E f w ( X w ( N ) ) 4 3 ( E f w ( X w ( N ) ) 2 ) 2 .   Plugging this into equation ( 6.30 ) and comparing with equation ( 6.29 ) shows that s v ( N ) 3 ( u v ( N ) ) 2   , completing the induction.   Proof of Lemma  6.2 : We observed in the proof of Lemma  6.5 that g v   is concave; it is bounded as well. For the upper bound, first note that g v ( x ) h ( x ) : = θ 2 v x 1 + κ ( θ v ) x   for some κ   varying continuously with θ v   . The proof of this is same as the proof of ( 4.21 ), using the Taylor expansion g v ( x ) = θ v x ( 1 + θ v ) ) θ v ( 1 θ v ) x 2 / 6 + O ( x 3 )   together with boundedness and concavity of g v   . Jensen's inequality gives E g v ( V ) E h ( V ) h ( E V )   which proves the upper bound with κ 1 = κ   .
For the lower bound, since g v ( x ) = θ v x O ( x 2 )   near 0, we have g v ( x ) θ v x λ x 2   for some λ   and all x   in some interval [ 0 , δ ]   . Choosing λ   larger if necessary, we can ensure that g v ( x ) θ v x λ x 2   for all x 0   . Hence E g v ( V ) θ 2 v E V c λ ( E V ) 2 .   Choose δ ( θ v ) > 0   so that the right-hand side is positive for x ( 0 , δ ( θ v ) )   . Choose κ 2 ( θ v )   so that θ 2 v x 1 + κ 2 ( θ v ) x [ θ 2 v x λ x 2 ] g v ( δ / 2 ) 4 c .   This satisfies ( 6.27 ) when E V δ   . But when E V > δ   , then hypothesis on V   implies that P ( V > δ / 2 ) 1 / ( 4 c )   and therefore that E g v ( V ) g v ( δ / 2 ) / ( 4 c )   . Hence ( 6.27 ) is valid for all x 0   . Together with the evident continuous dependence of κ i   on θ v   , this proves the lemma.  

7 Concluding remarks

Although we have in general no explicit probabilistic interpretation of L p   capacities, in the case of integer values of p   there is a more probabilistic formulation. Positive L p   capacity is equivalent to the existence of a probability measure μ   on T   such that p   independent paths picked from μ   will coincide along a path of finite average resistance. This corresponds to the representation of L p   -energy as a p   -fold integral over T   .
Finally, we remark that other statistical mechanical models lead to recursions similar to ( 5.22 ) but with functions f v   that are not necessarily concave. The Potts model with 1 < q < 2   is essentially similar to the Ising model, but when q > 2   , the functions f v   are not concave and qualitatively different behavior arises. See Häggström (1996) for a discussion of this as pertains to the random cluster model, and Pemantle and Steif (1999) for the Heisenberg and other continuous-state models on general trees.
Acknowledgement. Much of the research presented here was performed at the Mittag Leffler Institute. We are grateful to E. B. Dynkin for telling us about the importance of L p   capacities in connection with superprocesses.
References

  1. Bleher, P. M., Ruiz, J. and Zagrebnov V. A. (1995) On the purity of limiting Gibbs state for the Ising model on the Bethe lattice, J. Stat. Phys 79, 473–482.
  2. Chayes, J. T., Chayes, L., Sethna, J. and Thouless, D. (1986). Comm. Math. Phys. 106 41 89.
  3. Evans, W., Kenyon, C., Peres, Y. and Schulman, L. J. (2000). Broadcasting on trees and the Ising Model, Ann. Appl. Prob. 10, 410–433.
  4. H. O. Georgii (1988). Gibbs Measures and Phase Transitions. W. de Gruyter, Berlin.
  5. Häggström, O. (1996). The random-cluster model on a homogeneous tree. Probab. Th. Related Fields 104, 231–253.
  6. Ioffe, D. (1996a). A note on the extremality of the disordered state for the Ising model on the Bethe lattice. Lett. Math. Phys. 37, 137–143.
  7. Ioffe, D. (1996b) Extremality of the disordered state for the Ising model on general trees. Trees (Versailles, 1995), 3–14, Progr. Probab. 40, Birkhäuser, Basel.
  8. Lyons, R. (1989). The Ising model and percolation on trees and tree-like graphs. Comm. Math. Phys. 125 337 353.
  9. Lyons, R. (1992). Random walks, capacity, and percolation on trees. Ann. Prob. 20 2043 2088. Murakami, A. and Yamasaki, M. (1992). Nonlinear potentials on an infinite network. Mem. Fac. Sci. Shimane Univ. 26, 15–28
  10. Pemantle, R. and Steif, J. E. (1999). Robust phase transitions for Heisenberg and other models on general trees. Ann. Probab. 27, 876–912.
  11. Preston, C. J. (1974). Gibbs states on countable sets. Cambridge University Press: Cambridge.
  12. Preston, C. J. (1976). Random fields. Lecture notes in Mathematics, vol. 534. Springer: Berlin.
  13. Soardi, P. M. (1993). Morphisms and currents in infinite nonlinear resistive networks. Potential Anal. 2, 315–347.
  14. Soardi, P. M. (1994) Potential Theory on Infinite Networks. Lect. Notes Math. 1590, Springer, Berlin.