ypical behavior of the harmonic measure in critical Galton Watson trees with infinite variance offspring distribution Shen LIN LPMA, Université Pierre et Marie Curie, Paris, France E-mail: shen.lin.math@gmail.com arxiv:603.0200v2 math.pr] 25 Feb 207 Abstract We study the typical behavior of the harmonic measure in large critical Galton Watson trees whose offspring distribution is in the domain of attraction of a stable distribution with index, 2]. Let µ n denote the hitting distribution of height n by simple random wal on the critical Galton Watson tree conditioned on non-extinction at generation n. We extend the results of ] to prove that, with high probability, the mass of the harmonic measure µ n carried by a random vertex uniformly chosen from height n is approximately equal to n λ, where the constant λ > depends only on the index. In the analogous continuous model, this constant λ turns out to be the typical local dimension of the continuous harmonic measure. Using an explicit formula for λ, we are able to show that λ decreases with respect to, 2], and it goes to infinity at the same speed as 2 when approaches. Keywords. size-biased Galton Watson tree, reduced tree, harmonic measure, uniform measure, simple random wal and Brownian motion on trees. AMS 200 Classification Numbers. 60J80, 60G50, 60K37. Introduction In this paper, we continue our previous wor ] by extending its results to the critical Galton Watson trees whose offspring distribution may have infinite variance. o be more precise, let ρ be a non-degenerate probability measure on Z + with mean one, and we assume throughout this paper that ρ is in the domain of attraction of a stable distribution of index, 2], which means that ρr = r + r L r for any r 0,, 0 where the function Lx is slowly varying as x 0 +. he finite variance condition for ρ is sufficient for the previous statement to hold with = 2. When, 2, by classical results of 8, Chapters XIII and XVII], the condition is satisfied if and only if the tail probability ρ = ρx, + x varies regularly with exponent as x +. Supported in part by the grant ANR-4-CE25-004 ANR GRAAL
Under the probability measure P, for every integer n 0, we let n be a Galton Watson tree with offspring distribution ρ, conditioned on non-extinction at generation n. Conditionally given the tree n, we consider a simple random wal on n starting from the root. he probability distribution of its first hitting point of generation n will be called the harmonic measure µ n, which is supported on the set n n consisting of all vertices of n at generation n. Let q n > 0 be the probability that a critical Galton Watson tree 0 survives up to generation n. It has been shown by Slac 6] that, as n, the probability q n decreases as n up to multiplication by a slowly varying function, and q n # n n converges in distribution to a non-trivial limit distribution on R + that we will specify shortly. he main result of the present wor generalizes heorem of ] from the finite variance case with = 2 to all, 2]. heorem. Let Ω n be a random vertex uniformly chosen from n n. If the offspring distribution ρ has mean one and belongs to the domain of attraction of a stable distribution of index, 2], there exists a constant λ >, which only depends on, such that for every δ > 0, we have lim P n λ δ µ n Ω n n λ+δ =. 2 n Roughly speaing, the previous theorem asserts that if we loo at a typical vertex at level n of the conditional Galton Watson tree n, then this random vertex carries with high probability a mass of order n λ given by the harmonic measure µ n. Since λ >, we see particularly that the harmonic measure on the conditional Galton Watson tree is not uniformly spread and exhibits a fractal behavior. he fractal nature of the harmonic measure in critical Galton Watson tree has been initially investigated by Curien and Le Gall 5] from another perspective. o compare with heorem, their main result in the finite variance case has also been generalized to the present setting as heorem. in 0] : under the same assumption of heorem, we have the existence of another constant β 0, only depending on, such that for every δ > 0, we have the convergence in P-probability {v µ n n n : n β δ µ n v n β+δ} P. 3 n In other words, with high probability the mass given by the harmonic measure µ n to a random vertex at level n drawn with respect to the same harmonic measure is of order n β with β <. We can thus say that the harmonic measure µ n is mainly supported on a subset of size approximately equal to n β, which is much smaller than the whole size of n n. Notice that the family β, < 2 is shown in 0, heorem.3] to be bounded from above in R +. As the hitting distribution µ n of generation n by simple random wal on n is unaffected if we remove the branches of n that do not reach height n, we can simplify the model by considering merely simple random wal on the reduced tree n associated with n, which consists of all vertices of n that have at least one descendant at generation n. When the critical offspring distribution ρ has infinite variance, the study of scaling limits of n goes bac to Vatutin 7] and Yaymiv 8]. If we scale the graph distance by the factor n, as n the discrete reduced tree n n converges to a random continuous tree that we now describe. For every, 2], we define the -offspring distribution θ as follows. If = 2, we let θ 2 be the Dirac measure at 2. If < 2, θ is the probability measure on Z + given by θ 0 = θ = 0, θ = Γ! Γ2 = 2 3! 2, 2,
where Γ denotes the usual Gamma function of Euler. he mean of θ is equal to. o construct the random tree, we begin with the genealogical tree 0 of a family of particles that evolves in continuous time according to the following rules. At time 0 there is a single particle, this particle lives for a random time uniformly distributed over 0, ], then dies and gives birth to a random number of new particles. his number of offspring is distributed according to θ and is assumed to be independent of the lifetime of the initial particle. Inductively, all the particles evolve independently of each other, and each new particle appeared at time t 0, dies at a time uniformly distributed over t, ] and gives birth to an independent number of new particles according to the supercritical offspring distribution θ. After taing completion of 0 with respect to its natural intrinsic metric d, we obtain a random compact rooted R-tree that will be called the reduced stable tree of parameter. Its boundary is composed of all points of at height. We refer to Section 2. for a precise definition of. he description above maes clear the recursive structure of : let U be a uniform random variable over 0, ] and let N N be a random variable of law θ. We tae i i to be a sequence of independent copies of, and we assume the independence between U, N and i, i. If we attach to the top of a single line segment of length U the roots of the rescaled trees U,..., U N, the resulting tree rooted at the origin of the initial line segment will have the same distribution as. See Fig. for an illustration. U U 2 U N Height U Height 0 Figure : Recursive structure of the reduced stable tree As the continuous analog of simple random wal, Brownian motion on starting from the root can be defined up until its first hitting time of. It behaves lie linear Brownian motion as long as it stays inside a line segment of. It is reflected at the root of and when it arrives at a branching point, it chooses one of the adjacent line segments with equal probabilities. We define the continuous harmonic measure µ on as the quenched distribution of the first hitting point of by Brownian motion. he behavior of µ is closely related to that of the discrete harmonic measure µ n. In order to state a result analogous with heorem in the continuous setting, we must first mae sense of a typical point chosen from the boundary. o this end, we introduce another noncompact random rooted R-tree Γ endowed with the same branching structure as, such that each point of Γ at height y 0, corresponds bijectively to a point of 0 at height e y 0,. It is easy to verify that the resulting new tree Γ is a continuous-time Galton Watson tree of branching rate with the supercritical offspring distribution θ. In particular, 3
Γ 2 is the Yule tree which describes the genealogy of the classical Yule process. By definition, the boundary Γ of Γ is the set of all infinite geodesics in Γ starting from the root these will be called geodesic rays. Since and Γ share the same branching structure, both and Γ can be canonically identified with a common random subset of N N. For every r > 0, we write Γ r ] = exp r E#Γ r for the level set of Γ at height r, and we now that. By a martingale argument, one can define W := lim e r #Γ r r. As θ log <, the Kesten Stigum theorem for continuous-time Galton Watson trees, see e.g. 2, heorem III.7.2] implies that the previous convergence holds in the L -sense, EW ] =, and W > 0 almost surely. It follows from heorem III.8.3 in 2] that E e uw] = u + u for any u 0,. Moreover, according to heorem of 6], q n # n n converges in distribution to this martingale limit W as n. For every x Γ, we let Hx denote the height of x in Γ, and we write Γ x] for the tree of descendants of x in Γ, viewed as an infinite random R-tree rooted at x. For every r > 0, we write Γ r x] for the level set at height r of the tree Γ x]. If one thins of Γ x] as a subtree of Γ, the set Γ r x] consists of all the points of Γ at height r + Hx that are descendants of x. We similarly define W x] := lim r e r #Γ r which has the same distribution as W. he uniform measure ω on Γ is defined as the unique probability measure on Γ satisfying that, for every x Γ and for every geodesic ray v Γ passing through x, ω Bv, Hx = exp x], Hx W x] W, where Bv, Hx stands for the set of all geodesic rays in Γ that coincide with v up to height Hx. In existing literature, we also call ω the branching measure on the boundary of Γ. Since can be identified with Γ as explained above, we let ω be the random probability measure on induced by ω, which will be referred to as the uniform measure on. he following result is an extension of heorem 2 in ]. heorem 2. For every, 2], with the same constant λ as in heorem, we have P- a.s. ω dv-a.e. log µ B d v, r lim r 0 log r log ω B d v, r lim r 0 log r = λ, 4 =, 5 where B d v, r stands for the closed ball of radius r centered at v in the metric space, d. 4
By Duquesne and Le Gall 7, heorem 5.5], the Hausdorff measure of with respect to d is P-a.s. equal to. According to Lemma 4. in 2], assertion 5 of the preceding theorem implies that P-a.s. the uniform measure ω has the same Hausdorff dimension as the whole boundary of the reduced stable tree. Meanwhile, taing account of 4, we may interpret λ as the local dimension of the harmonic measure µ at a typical point of the boundary. In contrast, the Hausdorff dimension of the harmonic measure µ is P-a.s. equal to the constant β appearing in 3, since heorem.2 of 0] shows that P-a.s. µ dv-a.e., Since β < < λ, P-a.s. the set satisfies that µ B = whereas ω B = 0. log µ B d v, r lim = β. r 0 log r B = { v log µ B d v, r } : lim = β r 0 log r Corollary 3. For every, 2], P-a.s. the two measures µ and ω on the boundary of are mutually singular. When approaches, the Hausdorff dimension β of the harmonic measure µ remains bounded see heorem.3 in 0]. Naively, it would seem that the typical local dimension C λ of µ would behave asymptotically lie for some positive constant C. However, the explosion of λ turns out to be much faster than when decreases to. Proposition 4. We have that 0 < lim inf 2 λ lim sup 2 λ <. Our proof of Proposition 4 relies on the fact that the constant λ appearing in heorems and 2 can be expressed explicitly with help of the conductance of. Informally, if we thin of the random tree as an electric networ of resistors with unit resistance per unit length, the effective conductance C between the root and the boundary is a random variable larger than. From a probabilistic point of view, it is the mass under the Brownian excursion measure for the excursion paths away from the root that hit height. Following the definition of and its electric networ interpretation, the distribution of C satisfies the recursive distributional equation C d == U + U C + C 2 + + C N, 6 where C i i are i.i.d. copies of C, the integer-valued random variable N is distributed according to θ, and U is uniformly distributed over 0, ]. All these random variables are supposed to be independent. Proposition 5. For any, 2], the constant λ appearing in heorems and 2 is given by Moreover, λ is decreasing for all, 2]. λ = E W C ]. 7 5
he surprisingly simple formula 7 provides a unified expression of λ in terms of W and C. We will see in Section 2.5 that the product W C has the same first moment as the conductance Ĉ of a size-biased version of the reduced stable tree. Note that for < 2, the -offspring distribution θ is decreasing for the usual stochastic partial order. his property allows one to construct simultaneously all reduced stable trees,, 2] as a nested family, so that 2 for all < 2 2 see Section 2.4 in 0]. While the monotonicity of the typical local dimension λ is affirmed by the previous result, one question still open is whether the Hausdorff dimension β of the continuous harmonic measure µ is also decreasing with respect to. he rest of this paper is divided into three parts. he continuous model of Brownian motion on is studied in Section 2, where we prove heorem 2, Propositions 4 and 5. hen we set up notation and terminology for the discrete setting in Section 3, while Section 4 is devoted to proving heorem. he arguments are basically adapted from those used in the finite variance case. However, we emphasize that several modifications are indispensable to circumvent the problem of infinite variance. Instead of repeating some highly analogous reasoning, we refer the reader to ] for details that we omit. Acnowledgments. he author is grateful to an anonymous referee for many comments that greatly improved the paper. 2 he continuous setting 2. he reduced stable tree Let us begin with a formal definition of the reduced stable tree of parameter, 2]. We set V = N n, n=0 where by convention N = {, 2,...} and N 0 = { }. If v = v,..., v n V, we set v = n the generation of the vertex v in particular, = 0, and if n, we define the parent of v as v = v,..., v n and then say that v is a child of v. For two elements v = v,..., v n and v = v,..., v m belonging to V, their concatenation is vv := v,..., v n, v,..., v m. he notions of a descendant and an ancestor of an element of V are defined in the obvious way, with the convention that every v V is both an ancestor and a descendant of itself. An infinite subset Π of V is an infinite discrete rooted ordered tree without leaves if there exists a collection of positive integers v = v Π N for every v V such that Π = { } {v,..., v n V : v j v,...,v j for every j n}. Recall that the generating function of the -offspring distribution θ is given see for instance 6, p.74] as θ r = r + r, r 0, ]. 8 0 For fixed, 2], we introduce a collection K v v V of independent random variables distributed according to θ under the probability measure P, and define a random infinite discrete tree Π := { } {v,..., v n V : v j K v,..., v j for every j n}. 6
We point out that Π is just a supercritical Galton Watson tree with offspring distribution θ. In particular, Π 2 is an infinite binary tree. Let U v v V be another collection, independent of K v v V, consisting of independent real random variables uniformly distributed over 0, ] under the same probability measure P. We set now Y = U and then by induction, for every v Π different from the root, Y v = Y v + U v Y v. Note that a.s. 0 Y v < for every v Π. Consider the set 0 := { } 0, Y ] v Π \{ } here is a straightforward way to define a metric d on 0 R-tree and, for every x = v, r 0 of the tree 0 when < 2. {v} Y v, Y v ]., so that 0, d is a non-compact, we have d, 0, x = r. See Figure 2 for an illustration Height 2 2 22 23 Y Y 2 3 32 Y 3 2 3 Y Height 0 Figure 2: he random tree 0 when < 2 he reduced stable tree is the completion of 0 with respect to the metric d. It is immediate to see that, d is a compact R-tree, and = 0 where := {x : d, 0, x = } is the boundary of, which can be identified with a random subset of N N. he point, 0 is called the root of. For every x, we set its height Hx = d, 0, x. A genealogical order on can be defined by setting x y if and only if x belongs to the geodesic path from the root to y. Conditionally on, we can define Brownian motion B t t 0 on starting from the root and up to its first hitting time of. See e.g. 5, Section 2.] for the details of this construction. he harmonic measure µ is the distribution of B, which is a random probability measure on N N. 7
2.2 he continuous-time Galton Watson tree o introduce a new tree sharing the same branching structure as, we start with the same infinite discrete tree Π introduced in Section 2.. Consider now a collection V v v V of independent real random variables exponentially distributed with mean under the probability measure P. We set Z = V and then by induction, for every v Π different from the root, Z v = Z v + V v. he continuous-time Galton Watson tree hereafter to be called CGW tree for short of stable index is the set Γ := { } 0, Z ] v Π \{ } {v} Z v, Z v ], which is equipped with the metric d defined in the same way as d in the preceding subsection. For this metric, Γ is a.s. an infinite R-tree. For every x = v, r Γ, we eep the notation Hx = r = d, 0, x for the height of x. Observe that if U is uniformly distributed over 0, ], the random variable log U is exponentially distributed with mean. Hence we may and will suppose that the collection V v v V is constructed from the collection U v v V in the previous subsection via the formula V v = log U v for every v V. A homeomorphism from 0 onto Γ is given by the mapping Ψ defined as Ψv, r := v, log r for every v, r 0. By stochastic analysis, we can write for every t 0,, t ΨB t = W HB s 2 ds 0 where W t t 0 is Brownian motion with constant drift /2 towards infinity on the CGW tree Γ this process is defined in a similar way as Brownian motion on, except that it behaves lie Brownian motion with drift /2 on every line segment of the tree. It is easy to see that the Brownian motion W is transient. From now on, when we spea about Brownian motion on the CGW tree or on other similar infinite trees, we will always mean Brownian motion with drift /2 towards infinity. By definition, the boundary of Γ is the set of all geodesic rays in Γ starting from the root, 0, and it can be canonically embedded into N N. Due to the transience of Brownian motion on Γ, there is an a.s. unique geodesic ray denoted by W that is visited by W t t 0 at arbitrarily large times. he distribution of the exit ray W yields a random probability measure ν on N N. hans to 9, we have in fact ν = µ, provided we thin of both µ and ν as probability measures on N N. Infinite continuous trees. We let be the set of all pairs Π, z v v Π that satisfy the following conditions: i Π is an infinite discrete tree without leaves in the sense of Section 2.; ii z v 0, for all v Π ; iii z v < z v for every v Π\{ } ; iv for every v Π := {v, v 2,... N N : v, v 2,..., v n Π, n }, we have lim z n v,...,v n =. 9 8
If Π, z v v Π, we consider the associated tree := { } 0, z ] v Π\{ } {v} z v, z v ], equipped with the natural distance defined as above. he set Π can be identified with the boundary of the tree, which is defined as the collection of all geodesic rays in. We eep the notation Hx = r for the height of a point x = v, r. he genealogical order on is defined as previously and again is denoted by. If v = v, v 2,... Π and x = v, r, we write x v if v = v, v 2,..., v for some integer 0. When we say that we consider a tree, it means that we are given a pair Π, z v v Π satisfying the above properties, with being the associated tree. Clearly, for every, 2], the CGW tree Γ is a random element in, and we write Θ d for its distribution. Let us fix = Π, z v v Π. Under our previous notation, the root has offspring. We denote by, 2,..., the subtrees of rooted at the first branching point, z. o be more precise, for every i, we define the shifted discrete tree Πi] = {v V : iv Π}, and i is the infinite continuous tree corresponding to the pair Πi], z iv z v Πi]. Under Θ d, the offspring number is distributed according to θ. he branching property of the CGW tree means that under the conditional law Θ d, the subtrees,..., are i.i.d. following the same law Θ. If r > 0, the level set of at height r is r = {x : Hx = r}. If x r, let x] denote the subtree of descendants of x in. o define it formally, we write v x for the unique element of V such that x = v x, r, then x] is the tree associated with the pair Πv x ], z vxv r v Πvx]. As usual, r x] stands for the level set at height r of the tree x]. As we have seen in the Introduction, the martingale limit W = lim e r #r r exists Θ d -a.s., and W Θ d =. For every x, we similarly set W x] = lim r e r #r x]. If v is a geodesic ray passing through x, let Bv, Hx denote the set of geodesic rays in that coincide with v up to height Hx. hen Θ d -a.s., the uniform measure ω on is the probability measure on characterized by ω Bv, Hx = exp Hx W x] W for every x and v such that x v. For a fixed infinite continuous tree, we define the harmonic measure µ on as the distribution of the exit ray chosen by Brownian motion on with drift /2 towards infinity. 2.3 he invariant measure and the size-biased CGW tree We write := N N for the set of all pairs consisting of a tree and a distinguished geodesic ray v. We define a transformation S on by shifting, v = v, v 2,... at the first branching point of, that 9
is, S, v = v, ṽ, where ṽ = v 2, v 3,... and v is the tree rooted at the first branching point of that is chosen by v. he next result extends Proposition 5 in 0] from the Yule tree Γ 2 to the CGW tree Γ. Proposition 6. For all, 2], the probability measure W Θ d ω dv on is invariant under S. Proof. Fix a bounded measurable function F on. Under Θ d, the height z of the first branching point is exponentially distributed with mean, while the offspring number is distributed according to θ. Recall that θ = W = i= and e z W i. 0 Using these remars and the branching property of the CGW tree, we see that F S, vw Θ d ω dv = = = θ =2 i= θ e z e z dz =2 0 F, uw Θ d ω du, which shows the required invariance. F i, ue z W i Θ d = ω i du F, uw Θ d ω du o fully understand this invariant measure, we construct a size-biased version Γ of the CGW tree Γ, following the popularized idea of size-biasing a Galton Watson tree see e.g. 4, 3]. Let us start with the root. It lives for a random exponential lifetime with parameter, then it dies and simultaneously gives birth to a random number N of children, with N distributed as the size-biased offspring number N. We pic one of these children uniformly at random, say v {, 2,..., N }. We give independently the other children independent descendant trees distributed as Γ, whereas v lives for another independent exponential lifetime with parameter and then reproduces a random number N 2 of children with N 2 an independent copy of N. Again, we choose one of the children of v uniformly at random, call it v 2 {, 2,..., N 2 }, and give the others independent CGW descendant trees distributed as Γ. Repeating this procedure independently for infinitely many times, we obtain a random infinite tree Γ which will be called the size-biased CGW tree with parameter, 2]. A formal definition of Γ with the distinguished geodesic ray v = v, v 2,... as a random element of can be given similarly as in Section 2.6 of ]. Notice that the successive heights Z v of the vertices on the spine v are distributed as a homogeneous Poisson process on 0, with intensity. Keeping Proposition 6 in mind, we show by the next result why it is convenient to thin of Γ when studying a random geodesic ray in Γ sampled according to ω. Lemma 7. he pair Γ, v follows the distribution W Θ d ω dv. 0
2.4 he size-biased reduced tree Recall the bijection Ψ: v, r 0 v, log r Γ introduced in Section 2.2. If we apply its inverse Ψ v, s = v, e s to the size-biased CGW tree Γ and then tae the natural compactification, we obtain a random compact rooted tree called the size-biased reduced tree of parameter, 2]. A Poissonian description of is given as follows. First, under the mapping Ψ the geodesic ray v = v, v 2,... in Γ corresponds to a distinguished point at height in that we will still denote by v. Along the ancestral line of v in, we eep the same notation v for the branching points, with their respective heights Y v := exp Z v,, increasing to as. It is immediate to see that Y v is distributed on 0, as a Poisson process with intensity measure x x 0, dx. Moreover, if we set Y v0 = 0, the sequence V := Y v Y v Y v,, consists of i.i.d. random variables with the density function x over 0, ]. For any r > 0, we write r 0 for the same random tree as 0 with the distance d multiplied by the factor r, which means that the infinite discrete tree Π associated to 0 remains the same while all the uniform random variables U v, v V involved in the definition of 0 are replaced by ru v, v V. Conditionally on Y v, independently to each branching point v of height Y v on the ancestral line of v in, we graft a random number J of independent trees, all distributed according to Y v 0. he random number J is independent of those trees grafted at v, and it has the same distribution as N. As in the case of the size-biased CGW tree Γ, we assume the independence of J and the trees at grafted v among all different levels. v Height Y v v 3 J 3 = 2 v 2 J 2 = Y v2 v J = 3 Y v Height 0 Figure 3: Schematic representation of the size-biased reduced tree when < 2 o finally get the size-biased reduced tree, we tae the completion of the tree obtained above with respect to its natural distance see Fig. 3. his compactification is equivalent to
adding all the points on the boundary of the grafted trees. All these points together with the distinguished one v form the boundary, which is defined as the set of all points in that are at height. Remar. he -stable continuum random tree arises as the scaling limit of rescaled discrete Galton Watson trees with a critical offspring distribution in the domain of attraction of a stable distribution of index, 2] see e.g. the monograph 6] of Duquesne and Le Gall. Under the same assumption, if we condition the critical Galton Watson tree to survive up to generation n and rescale it by the factor n, after letting n, we get as limit an -stable tree conditioned to reach height see heorem 4. in 7], from which the reduced stable tree can be obtained as a pruned tree by eeping only the points at height less than that possess a descendant that reaches height. he probabilistic structure of the size-biased reduced tree can also be derived from a spinal decomposition heorem 4.5 in 7] of the conditional -stable tree along the ancestral line of a typical point at height. Since this point of view is not needed for proving our main results, we omit the details. 2.5 he continuous conductance Recall that the random variable C is defined, from an electric networ point of view, as the effective conductance between the root and the boundary in. his is equivalent to defining C as the mass assigned by the excursion measure away from the root for Brownian motion with no drift on, to the set of trajectories that reach height before coming bac to the root. his probabilistic definition of the conductance has the advantage of being easily generalized to a more general setting. If is an infinite continuous tree, we define its conductance C between the root and as the mass assigned by the excursion measure away from the root for Brownian motion with drift /2 on, to the set of trajectories that go to infinity without returning to the root. By the correspondence between the CGW tree Γ and the reduced stable tree, we can verify that CΓ has the same law as C. See 5, Section 2.3] for details. In Section 2.3 of 0], it is shown that for, 2] the law of C is characterized by the distributional identity 6 in the class of all probability measures on, ]. Now let us tae into account the martingale limit W. From 0 and the branching property of Γ, it follows that the joint distribution of W, C satisfies the distributional equation W, C d == U W + + W N, U + U C + + C N, in which W i, C i i are i.i.d. copies of W, C, the integer-valued random variable N is distributed according to θ, and U is uniformly distributed over 0, ]. All these random variables are supposed to be independent. Under P, we define for every, 2] a positive random variable Ĉ distributed as C under the probability measure W Θ d. If the size-biased reduced tree is viewed as an electric networ of ideal resistors with unit resistance per unit length, we will see that its effective conductance between the root and has the same law as Ĉ. For this reason, we sometimes call Ĉ the size-biased conductance. Lemma 8. he random variable Ĉ satisfies the distributional identity Ĉ d V == V + Ĉ + C 2 + + C, 2 N 2
where in the right-hand side the random variable V has density function x over 0, ], the integer-valued random variable N has the size-biased distribution of θ, and C i i 2 are i.i.d. copies of C. All these random variables V, N, Ĉ and C i i 2 are independent. Proof. he law of Ĉ is determined by EgĈ ] = EW gc ] for every nonnegative measurable function g on R +. Using under the same notation, we have ] E gĉ = E = E U W + + W N g U + N U W g U + Defining N as the size-biased version of N, we see that ] E gĉ = EN ] E U W g U + U C C + + C N U + + C N ] U C + + C N. ] ] Recall that EN ] = = E U ]. he statement of the lemma thus follows if we let V be a random variable independent of N, Ĉ and C i i 2, satisfying EfV ] = E for every nonnegative measurable function f. U fu ] he law γ of the conductance C has been studied at length in Proposition 2. of 0]. We now discuss some similar properties of Ĉ. For every v 0,, n 2, x, and c i i 2, N, we define Gv, n, x, c i i 2 := so that 2 can be reformulated as v + v x + c 2 + + c n. 3, 4 d Ĉ == GV, N, Ĉ, C i i 2 5 where V, N, Ĉ, C i i 2 are as in 2. Let M be the set of all probability measures on, ] and let Φ : M M map a probability distribution σ to Φ σ = Law GV, N, X, C i i 2 where X is distributed according to σ and independent of V, N, C i i 2. Proposition 9. We fix the parameter, 2]. he distributional equation 2 characterizes the law γ of Ĉ in the sense that, γ is the unique fixed point of the mapping Φ on M, and for every σ M, the -th iterate Φ σ converges to γ wealy as. 3
2 he law γ has a continuous density over,. 3 For any monotone continuously differentiable function g :, R +, we have ] EĈ Ĉ g Ĉ + ] E gĉ = E gĉ + C 2 + + C where Ĉ, N, and C i i 2 are as in 5. 4 We define, for all l 0, the Laplace transforms ϕ l := Eexp l C /2] and ϕ l = Eexp l Ĉ /2] := e lx/2 γ dx. ] N, 6 hen ϕ solves the linear differential equation 2l ψ l + l ψ l ϕ l ψl = 0. 7 5 If < < 2, only the first moment of Ĉ is finite, whereas all moments of Ĉ2 are finite. Proof. he results for = 2 are stated without proof in ]. he arguments given below is similar in spirit to that of Proposition 2. in 0], but due to the fact that N has an infinite mean for < 2, sharper estimates are necessary in the present setting. First of all, the stochastic partial order on M is defined by saying that σ σ if and only if there exists a coupling X, Y of σ and σ such that a.s. X Y. It is clear that for any, 2], the mapping Φ is increasing for the stochastic partial order. We endow the set M of all probability measures on, ] that have a finite first moment with the -Wasserstein metric d σ, σ := inf { E X Y ] : X, Y coupling of σ, σ }. he metric space M, d is Polish and its topology is finer than the wea topology on M. Let us show that Φ maps M into M for any, 2]. For any a, b > 0 and 0 < r <, the weighted harmonic-geometric means inequality says that r a + r a r b r. b We fix some r 2, and apply the preceding inequality to obtain V r r r r V + X + C 2 + + C X + C 2 + + C r. V V N N It suffices to show that the right-hand side of the last display has a finite first moment provided that X has a finite mean. Notice that by 3, ] ] E V r V r = E U U r U r is finite. On the other hand, since a.s. X and C i for all i 2, X E + C 2 + + C ] r N X + C 2 + + C ] N E r N = θ r EX] + E C ] ] 2 + + E C. 8 2 4
It is shown in 0, Proposition 2.] that C has a finite mean. Meanwhile, we now by Stirling s formula that θ = O + when. As r > 2, it follows that 2 r θ <. he sum in 8 is finite and therefore GV, N, X, C i i 2 has a finite first moment. Next, we observe that Φ is strictly contractant on M, d. o see this, let X, Y be a coupling between σ, σ M under the probability measure P, and assume that V, N, C i i 2 and X, Y are independent under P. hen GV, N, X, C i i 2 and GV, N, Y, C i i 2 provide a coupling of Φ σ and Φ σ. Using the fact that a.s. X, Y, C i, we have GV, N, X, C i i 2 GV, N, Y, C i i 2 V V = V + X + C 2 + + C V + Y + C N 2 + + C N V X Y = V X + C 2 + + C + V V Y + C N 2 + + C + V N V + N X Y. V 2 By definition, ] V E + N V 2 = 2 ] V θ E + V 2 = ] U θ E + U 2 2 ] U θ E + U 2 2 where U is uniformly distributed over 0, ]. Notice that for for any integer 2, ] U E + U 2 = + log 2 <. hus, taing expectation and minimizing over the couplings between σ and σ, we get ] d Φ σ, Φ σ V E + N d σ, σ V 2 N N log N ] + E N 2 d σ, σ =: c d σ, σ with c <. So for, 2], the mapping Φ is contractant on M and by completeness it has a unique fixed point γ in M. Furthermore, for every σ M, we have Φ σ γ for the metric d, hence also wealy, as. However, γ is also a fixed point of Φ according to 5. he equality γ = γ will follow if we can verify that γ is the unique fixed point of Φ in M. o this end, it will be enough to show that Φ σ γ as, for every σ M. Let us apply Φ to the Dirac measure δ at infinity to see that Ṽ Φ δ = Law, V Φ 2 δ = Law V + Ṽ 5 + C 2 + + C N,
where Ṽ is a copy of V being independent of V, N and C i i 2 under P. We claim that Φ 2 δ has a finite first moment. In fact, applying 3 yields E V + Ṽ V ] + C 2 + + C = E N E We integrate with respect to U to obtain E U + Ṽ U U + Ṽ = 2 θ E U + U + Ṽ U U + C 2 + + C N Ṽ + C 2 + + C N ] U + C 2 + + C U ] Ṽ log + C 2 + + C ] + C 2 + + C = E Ṽ + C 2 + + C Ṽ 2 E log + C 2 + + C ]. Using the inequality loga + b a + log b for a, b, we see that Ṽ E log + C 2 + + C ] ] /2 E Ṽ + E log C 2 + + C ]. ] he random variable Ṽ /2 has clearly a finite mean. Meanwhile, for any r 0,, we have log a r ar for every a, and thus E log C 2 + + C ] r E C 2 + + C r ] r E C If we tae r <, by the same arguments following 8 we deduce that 2 θ E log C 2 + + C ] r 2 2 + + C ] r.. ] θ E C 2 + + C ] r <. 9 herefore, we get Φ 2 δ M. By monotonicity, Φ 2 σ M for every σ M, and from the preceding results we have Φ σ γ for every σ M. his implies that γ = γ is the unique fixed point of Φ in M. 2 For every t, we set F t := PĈ t, and F t := PĈ + C 2 + + C t for every integer 2. It follows from 2 that, for every t >, V F t = P V + Ĉ + C 2 + + C = = =2 t /t θ dv v 0 t t dx 6 x N x 2 F t t vt vt E N F N x ], 20
where N is distributed according to θ. Since obtain from the last display that F t = for every t, 2] and 2, we where A 0 := 2 t F t = t 2 dx x A 0 x 2, t, 2], 2 E N N F x ], 2. From 20 and 2, we observe that F is continuous on,. Recall that C has a continuous density see 0, Proposition 2.]. It follows that all functions F, 2 are continuous on,. By dominated convergence the function x E N F N x ] is also continuous on,. Using 20 and 2 again, we see that F is continuously differentiable on,. Consequently γ has a continuous density on,. 3 o derive the identity 6, we first differentiate 20 with respect to t to see that tt d F t dt F t = E N N F t ] holds for t,. Let g :, R + be a monotone continuously differentiable function. We multiply both sides of the last display by g t and integrate for t from to to obtain E Ĉ Ĉ g Ĉ ] + E gĉ ] = E N gĉ + C 2 + + C N ], which gives 6. 4 aing gx = exp xl/2 for l > 0 in 6, we obtain 7 by using the generating function of N given in 8. 5 he first moment of Ĉ is finite because γ M. aing gx = x in 6 yields E Ĉ 2] = E Ĉ] + E N ] E C ]. For every, 2, as E N ] = the second moment of Ĉ is infinite. Finally, by Proposition 6 in 5], all moments of C 2 are finite. Using E Ĉ2 Ĉ2 g Ĉ2 ] + 2 E gĉ2 ] = 2 E gĉ2 + C 2 ] with gx = x n, n, one can easily show by induction that all moments of Ĉ2 are finite. 2.6 Proof of heorem 2 he following proof of heorem 2 for all, 2] proceeds in much the same way as in the case = 2. We sip some details, for which we refer the reader to Section 2.5 in ]. First, we have seen in Proposition 6 that the probability measure W Θ d ω dv on is invariant under the shift S. aing into account that W > 0, Θ d -a.s., we can verify, in a similar fashion as in 0, Proposition 2.6], that the shift S acting on the probability space, W Θ d ω dv is ergodic. 7
Let H n, v denote the height of the n-th branching point on the geodesic ray v. For every n, it holds H n = n i=0 H S i where S i stands for the i-th iterate of the shift S. It follows from the ergodic theorem that W Θ d ω n H n n One calculates as in the proof of Proposition 6 to get H, v W Θ d ω dv = dv-a.s. and thus Θ d ω dv-a.s. H, v W Θ d ω dv. 22 z W Θ d =. 23 For a fixed geodesic ray v = v, v 2,..., we write x n,v = v,..., v n, H n+, v for the n + -st branching point on v, and we set F n, v := log ω {u : x n,v u}. Using the ergodic theorem again, we have Θ d ω dv-a.s., n F n n where the limit is identified with F, v W Θ d ω dv, z W Θ d =. 24 he calculation can be done analogously as in display 8 of ]. We only point out that we need W log W Θ d <, which is true because θ log 2 < cf. heorem I.0.2 in 2]. By considering the ratio F n /H n and taing n, we derive that Θ d -a.s. ω dv-a.e., lim r r log ω Bv, r =, from which the convergence 5 readily follows. hen we turn to the harmonic measure µ and set G n, v := log µ {u : x n,v u}. Using the flow property of µ see Lemma 2.3 in 0] and the ergodic theorem, we obtain the Θ d ω dv-almost sure convergence n G n n We can evaluate the preceding limit as i= = e z W i log G, v W Θ d ω dv. 25 W log C i C + + C Θ d C C + + C Θ d. Putting 25 together with 22 and 23, we see that Θ d -a.s. ω dv-a.e., lim r r log µ Bv, r = W log C C + + C Θ d. 8
which implies that the convergence 4 holds with the limit λ := E N W log C + + C ] N C, 26 where in the right-hand side we eep the same notation as in. Finally, it remains to chec that λ defined by 26 satisfies < λ <. Since C, an upper bound for λ is E N W log + C 2 + + C ] N, which equals EN log+c 2 + +C N ] by independence and the fact that EW ] =. For any fixed r 0,, there exists a positive constant M < such that log + x M + x r for every x > 0. For similar reasons as in 9, we have λ θ E log + C 2 2 + + C ] M + E C ] r θ +r <. 2 On the other hand, from 23 and the display following 25 we also now that λ = e z W i log C + + C C i= Θ d. i By concavity of the logarithm, i= W i log C + + C C i i= W i log W + + W W. i his inequality is strict if and only if for some i {,..., }, W i W + + W C i C + + C. Since the latter property holds with positive probability under Θ d, λ > = e z W i log W + + W W Θ d i i= F, v W Θ d ω dv. By 24, the right-hand side of the last display is equal to. herefore, we have λ >. 2.7 Proof of Proposition 5 According to its definition 26, λ = Ĉ + C E 2 + + C ] N log. 27 Ĉ 9
Since we can identify the latter with EĈ ] by taing gx = logx in 6, we get the identity λ = EĈ ] = EW C ]. Let us show the monotonicity of λ to complete the proof of Proposition 5. Indeed, we will prove that the law γ of the size-biased conductance decreases for, 2] in the sense of the usual stochastic partial order. Observe that the function Gv, n, x, c i i 2 defined in 4 is decreasing with respect to v, increasing with respect to n, increasing with respect to c i for every i 2. Since PV > x = x for every x 0, ], the random variable V is increasing with respect to, 2] for the stochastic order. Meanwhile, we deduce from 8 that Putting it into the identity E r N ] = r r r, r 0, ]. 28 r P N = = and recalling that a.s. N 2, we get =0 r P N > = Er ˆN ] r r 2 P N = r 2, =3 which implies that for every integer 3, P N 2 3 + = 2! = θ. By taing the derivative with respect to, we see that the size-biased offspring number N is decreasing with respect to, 2] for the stochastic order. Recall that the conductance C is also decreasing with respect to, 2] for the stochastic order see Section 2.4 of 0]. Now we tae < 2 2. According to the previous discussion, for any fixed x, G V 2, N 2, x, C 2 i i 2 G V, N, x, C i i 2, which implies that Φ 2 σ Φ σ for any probability distribution σ on, ]. By Proposition 9, for any, 2], the law γ of Ĉ has no atom and the iterates of Φ σ converge wealy to γ. As a result, Ĉ 2 is stochastically dominated by Ĉ. herefore, we have shown that the law of the size-biased conductance Ĉ is decreasing with respect to, 2] for the stochastic order. In particular, the first moment EĈ ] decreases with respect to, 2] and so does λ. 20
2.8 he asymptotic behavior of λ as We give here the proof of Proposition 4. First of all, let us collect some facts about the Riemann zeta function ζs = n /ns. It is well-nown that ζs has a simple pole at s = with residue. hen for the derivative ζ s = logn n n, we have s lim 2 ζ =. 29 he following inequality for Gamma function ratios see for example, page 4] is also needed. For every integer 2 and every 0 x, x Γ + x Γ x, One can thus find two postive constants C, C 2 such that for any, 3/2 and any 2, θ = Starting from 27 and the fact that Ĉ, C, we have λ E log + C 2 + + C ] N Γ Γ2 Γ + C +, C 2 +. 30 E log C + C 2 + + C ], N where in the last expectation we introduced another copy C of C independent of all other random variables. Let U be a sequence of independent random variables uniformly distributed over 0, ]. Since each C is dominated by U for the stochastic partial order because of 6, E log C + C 2 + + C ] E log U + U2 + + U ]. In 0, page ], we have already shown that from which it follows that E log U + U2 + + U ] hus, together with 30 and 29, we get λ θ 2 + log. 2 2 + log, lim sup 2 λ lim 2 4C 2 log = 4C 2 <. 2 For a lower bound, let us write Since Ĉ, C, Ĉ λ = E log + C 2 + + C ] E log N Ĉ]. 3 Ĉ E log + C 2 + + C ] E log N ] = N 2 θ log. 2
Still using 30 and 29, we see that and hence lim inf lim inf 2 θ log C, 2 Ĉ E log + C 2 + + C ] C > 0. 32 N On the other hand, let V 0, be a random variable of density function x. Due to 2, the size-biased conductance Ĉ is dominated by V for the stochastic order, and thus Elog Ĉ ] Elog V ]. However, we define the analytic function Hx := 0 t x dt, for x, t that interpolates the harmonic numbers Hn = n = /, n. An integration by parts shows that H = ] logx x dx = E log V. 0 Since Hx logx converges to the famous Euler-Mascheroni constant as x goes to +, we have ] lim E log Ĉ] = lim E log V = 0. Combining this with 3 and 32, we obtain lim inf he proof of Proposition 4 is therefore completed. 3 he discrete setting 3. Galton Watson trees 2 λ C > 0. We briefly introduce the general notion of discrete rooted ordered trees, which extends the case of infinite trees without leaves that we have seen in Section 2.. A discrete rooted ordered tree t is by definition a subset of V such that the following holds: i t ; ii If u = u,..., u n t\{ }, then ū = u,..., u n t ; iii For every u = u,..., u n t, there exists an integer u t 0 such that, for every j N, u,..., u n, j t if and only if j u t. he quantity u t above is called the number of children of u in t. A vertex with no child is called a leaf. From now on, we will just say tree instead of discrete rooted ordered tree for short. We write for the non-strict genealogical order on t. A tree t is always viewed as a graph whose vertices are the elements of t and whose edges are the pairs {ū, u} for all u t\{ }. Since the lexicographical order on the vertices is not really used in our arguments, we will not pay much attention to this order structure. 22
For an infinite tree t, we say it has a single infinite line of descent if there exists a unique sequence of positive integers u n n such that u, u 2,..., u n t for all n. he height of a tree t is written as ht := sup{ u : u t}. he set of all vertices of t at generation n is denoted by t n := {u t: u = n}. If u t, the tree of descendants of u is tu] := {w V : uw t}. Let t be a tree of height larger than n, and consider a simple random wal X = X 0 on t starting from the root, which is defined under the probability measure P t. We write τ n := inf{ 0: X = n} for the first hitting time of generation n by X, and we define the discrete harmonic measure µ t n supported on t n as the law of X τn under P t. Critical Galton Watson trees. For every integer n 0, we let n be a discrete Galton Watson tree with critical offspring distribution ρ, conditioned on non-extinction at generation n, viewed as a random subset of V. We suppose that the random trees n, n 0 are defined under the probability measure P. We let n be the reduced tree associated with n, consisting of all vertices of n that have at least one descendant at generation n. Note that u n for every u n. When taling about a reduced tree, we always implicitly assume that the correct relabeling of the vertices has been done to mae it a tree without changing the genealogical order. Set q n := P h 0 n the non-extinction probability up to generation n for a Galton Watson tree of offspring distribution ρ. If L is the slowly varying function appearing in, it has been established in 6, Lemma 2] that qn Lq n n as n. 33 By the asymptotic inversion property of slowly varying functions see e.g. 3, Section.5.7], it follows that q n n ln as n, for some other function l slowly varying at. Size-biased Galton Watson tree. We denote by a size-biased Galton Watson tree with offspring distribution ρ. It is defined similarly as the size-biased CGW tree Γ. We refer to Lyons, Pementle and Peres 3] for the details. he unique infinite line of descent v, v 2,... in is called its spine. If N denotes a random variable distributed according to the size-biased distribution of ρ, that is, P N = = ρ for every 0, then the offspring numbers of v are i.i.d. copies of N. Let ] n be the finite tree obtained from by eeping only its first n generations, i.e., ] n := {u : u n}. As shown in 3], the law of the random tree ] n is that of n biased by # n n. Moreover, conditionally given the first n levels of, the vertex v n on the spine is uniformly distributed on the n-th level of. For every integer n, let ] n be the finite tree obtained from by erasing the infinite tree of descendants of the vertex v n. By convention, the vertex v n itself is ept in ] n. We emphasize that in general ] n ] n, since the height of ] n can be strictly larger than n. We let ] n be the reduced tree associated with the tree ] n up to generation n, which consists of all vertices of ] n that have at least one descendant at generation n. Notice that ] n is also the reduced tree associated with ] n up to generation n. 23
3.2 Convergence of discrete conductances ae a tree t of height larger than n and consider the new tree t obtained by adding to the graph t an edge between its root and an extra vertex. We define under the probability measure P t a simple random wal X on t starting from the root. Let τ be the first hitting time of by X, and for every integer i n, let τ i be the first hitting time of generation i of the tree t by X. We write C i t := P t τ i < τ. his notation is justified by the fact that C i t can be interpreted as the effective conductance between and the i-th generation of t in the graph t, see e.g. 4, Chapter 2]. Here are some basic facts concerning the conductance of the reduced Galton Watson tree n. Proposition 0. Suppose that the critical offspring distribution ρ is in the domain of attraction of a stable distribution of index, 2]. hen it holds that n C n n d n C. Recall that C is the conductance of the continuous reduced tree. Given the convergence result of rescaled discrete trees n n towards Proposition 3.2 in 0], the previous proposition can be shown in the same way as Proposition in ]. he following moment estimate for the conductance C n n is Lemma 3.9 in 0]. Lemma. For every r 0,, there exists a constant K = Kr, ρ depending on r and the offspring distribution ρ such that, for every integer n, E n Cn n r ] K. 3.3 Bacward size-biased Galton Watson tree In the proof of heorem, we will need a rear-view variant of the size-biased Galton Watson tree, which can be traced bac to the inflated Galton Watson tree introduced by Peres and Zeitouni in 5]. Let us explain succinctly how this bacward size-biased Galton Watson tree is constructed. More details can be found in Section 3.4 of ]. First, the random tree has a unique infinite ray of vertices u 0, u, u 2,..., referred to as its spine. We fix a genealogical order on the spine by declaring that, for every n 0, the vertex u n is at generation n of. Under this rule, u 0 is a child of u, which is a child of u 2, etc. o describe the finite subtrees in branching off every node of the spine, we tae an i.i.d. sequence N n n of random variables following the size-biased distribution of ρ. Independently to each vertex u n, n, we graft N n independent copies of an ordinary ρ-galton Watson descendant tree, so that u n has exactly N n children: one of them, u n, is on the spine while the rest of them are the roots of independent ρ-galton Watson trees. Last but not least, we define the genealogical order on by eeping the genealogical orders inherited from the grafted Galton Watson trees and combining them with the genealogical order on the spine. For instance, u 2 is an ancestor of any vertex in one of the subtrees grafted at u. he notion of generation for every vertex in can also be defined in a consistent manner: for any vertex v not on the spine, there is a unique vertex u m on the spine such that v belongs to a finite subtree grafted at u m, then we say that the generation of v in is equal to m + plus the initial generation of v inside the corresponding grafted tree. See Figure 4 for an illustration. 24