Eigenvalue distribution in random trees Stephan Wagner (joint work - - PowerPoint PPT Presentation

eigenvalue distribution in random trees
SMART_READER_LITE
LIVE PREVIEW

Eigenvalue distribution in random trees Stephan Wagner (joint work - - PowerPoint PPT Presentation

Eigenvalue distribution in random trees Stephan Wagner (joint work with Kenneth Dadedzi) Naples, 13 September 2018 Tree eigenvalues To every graph, and in particular every tree, we can associate an adjacency matrix, a Laplacian matrix, and


slide-1
SLIDE 1

Eigenvalue distribution in random trees

Stephan Wagner (joint work with Kenneth Dadedzi) Naples, 13 September 2018

slide-2
SLIDE 2

Tree eigenvalues

To every graph, and in particular every tree, we can associate an adjacency matrix, a Laplacian matrix, and several other interesting matrices. v1 v2 v3 v4 v5 A =       1 1 1 1 1 1 1 1       L =       1 −1 −1 2 −1 −1 3 −1 −1 −1 1 −1 1      

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

2 / 33

slide-3
SLIDE 3

Tree eigenvalues

We will be interested in the eigenvalues of these matrices. v1 v2 v3 v4 v5 A =       1 1 1 1 1 1 1 1       L =       1 −1 −1 2 −1 −1 3 −1 −1 −1 1 −1 1       Eigenvalues of A: −1.84776, −0.765367, 0, 0.765367, 1.84776 Eigenvalues of L: 0, 0.518806, 1, 2.31111, 4.17009

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

3 / 33

slide-4
SLIDE 4

A large random tree

A tree with 100 vertices (left) and the distribution of the eigenvalues of its adjacency matrix (right).

  • 2
  • 1

1 2

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

4 / 33

slide-5
SLIDE 5

A large random tree

A tree with 100 vertices (left) and the distribution of the eigenvalues of its Laplacian matrix (right).

1 2 3 4 5 6

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

5 / 33

slide-6
SLIDE 6

The spectral distribution

Definition

Let T be a fixed tree with n vertices. Picking one of the n eigenvalues (counted with multiplicity) uniformly at random, we obtain a random variable XT : P(XT = x) = multiplicity of x as an eigenvalue of T n . The distribution of this random variable is the spectral distribution of T.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

6 / 33

slide-7
SLIDE 7

The spectral distribution

Definition

Let T be a fixed tree with n vertices. Picking one of the n eigenvalues (counted with multiplicity) uniformly at random, we obtain a random variable XT : P(XT = x) = multiplicity of x as an eigenvalue of T n . The distribution of this random variable is the spectral distribution of T. What can we say about the spectral distribution of a large random tree T?

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

6 / 33

slide-8
SLIDE 8

Uniform models

The simplest type of model uses the uniform distribution on the set of trees of a given order within a specified family (e.g. the family of all labelled trees, all unlabelled trees or all binary trees).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

7 / 33

slide-9
SLIDE 9

Uniform models

The simplest type of model uses the uniform distribution on the set of trees of a given order within a specified family (e.g. the family of all labelled trees, all unlabelled trees or all binary trees). The analysis of such models often involves exact counting and generating functions.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

7 / 33

slide-10
SLIDE 10

Uniform models

The simplest type of model uses the uniform distribution on the set of trees of a given order within a specified family (e.g. the family of all labelled trees, all unlabelled trees or all binary trees). The analysis of such models often involves exact counting and generating functions. In particular, this is the case for simply generated families of trees.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

7 / 33

slide-11
SLIDE 11

Simply generated trees

On the set of all rooted ordered (plane) trees, we impose a weight function by first specifying a sequence 1 = w0, w1, w2, . . . and then setting w(T) =

  • i≥0

wKi(T)

i

, where Ki(T) is the number of vertices of outdegree i in T. Then we pick a tree of given order n at random, with probabilities proportional to the

  • weights. For instance,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

8 / 33

slide-12
SLIDE 12

Simply generated trees

On the set of all rooted ordered (plane) trees, we impose a weight function by first specifying a sequence 1 = w0, w1, w2, . . . and then setting w(T) =

  • i≥0

wKi(T)

i

, where Ki(T) is the number of vertices of outdegree i in T. Then we pick a tree of given order n at random, with probabilities proportional to the

  • weights. For instance,

w0 = w1 = w2 = · · · = 1 generates random plane trees,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

8 / 33

slide-13
SLIDE 13

Simply generated trees

On the set of all rooted ordered (plane) trees, we impose a weight function by first specifying a sequence 1 = w0, w1, w2, . . . and then setting w(T) =

  • i≥0

wKi(T)

i

, where Ki(T) is the number of vertices of outdegree i in T. Then we pick a tree of given order n at random, with probabilities proportional to the

  • weights. For instance,

w0 = w1 = w2 = · · · = 1 generates random plane trees, w0 = w2 = 1 (and wi = 0 otherwise) generates random binary trees,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

8 / 33

slide-14
SLIDE 14

Simply generated trees

On the set of all rooted ordered (plane) trees, we impose a weight function by first specifying a sequence 1 = w0, w1, w2, . . . and then setting w(T) =

  • i≥0

wKi(T)

i

, where Ki(T) is the number of vertices of outdegree i in T. Then we pick a tree of given order n at random, with probabilities proportional to the

  • weights. For instance,

w0 = w1 = w2 = · · · = 1 generates random plane trees, w0 = w2 = 1 (and wi = 0 otherwise) generates random binary trees, wi = 1

i! generates random rooted labelled trees.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

8 / 33

slide-15
SLIDE 15

Branching processes

A classical branching model to generate random trees is the Galton-Watson tree model: fix a probability distribution on the set {0, 1, 2, . . .}.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

9 / 33

slide-16
SLIDE 16

Branching processes

A classical branching model to generate random trees is the Galton-Watson tree model: fix a probability distribution on the set {0, 1, 2, . . .}. Start with a single vertex, the root.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

9 / 33

slide-17
SLIDE 17

Branching processes

A classical branching model to generate random trees is the Galton-Watson tree model: fix a probability distribution on the set {0, 1, 2, . . .}. Start with a single vertex, the root. At time t, all vertices at level t (i.e., distance t from the root) produce a number of children, independently at random according to the fixed distribution (some of the vertices might therefore not have children at all).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

9 / 33

slide-18
SLIDE 18

Branching processes

A classical branching model to generate random trees is the Galton-Watson tree model: fix a probability distribution on the set {0, 1, 2, . . .}. Start with a single vertex, the root. At time t, all vertices at level t (i.e., distance t from the root) produce a number of children, independently at random according to the fixed distribution (some of the vertices might therefore not have children at all). A random Galton-Watson tree of order n is obtained by conditioning the process.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

9 / 33

slide-19
SLIDE 19

Branching processes

A classical branching model to generate random trees is the Galton-Watson tree model: fix a probability distribution on the set {0, 1, 2, . . .}. Start with a single vertex, the root. At time t, all vertices at level t (i.e., distance t from the root) produce a number of children, independently at random according to the fixed distribution (some of the vertices might therefore not have children at all). A random Galton-Watson tree of order n is obtained by conditioning the process. Simply generated trees and Galton-Watson trees are essentially equivalent. For example, a geometric distribution for branching will result in a random plane tree, a Poisson distribution in a random rooted labelled tree.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

9 / 33

slide-20
SLIDE 20

Simply generated and Galton-Watson trees

An example: Consider the Galton-Watson process based on a geometric distribution with P(X = k) = pqk (p = 1 − q).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

10 / 33

slide-21
SLIDE 21

Simply generated and Galton-Watson trees

An example: Consider the Galton-Watson process based on a geometric distribution with P(X = k) = pqk (p = 1 − q). The tree above has probability p7(pq)2(pq2)2(pq3)2 = p13q12,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

10 / 33

slide-22
SLIDE 22

Simply generated and Galton-Watson trees

An example: Consider the Galton-Watson process based on a geometric distribution with P(X = k) = pqk (p = 1 − q). The tree above has probability p7(pq)2(pq2)2(pq3)2 = p13q12, as does every tree with 13 vertices.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

10 / 33

slide-23
SLIDE 23

Random increasing trees

Another random model that produces very different shapes uses the following simple process, which generates random recursive trees:

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

11 / 33

slide-24
SLIDE 24

Random increasing trees

Another random model that produces very different shapes uses the following simple process, which generates random recursive trees: Start with the root, which is labelled 1.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

11 / 33

slide-25
SLIDE 25

Random increasing trees

Another random model that produces very different shapes uses the following simple process, which generates random recursive trees: Start with the root, which is labelled 1. The n-th vertex is attached to one of the previous vertices, uniformly at random.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

11 / 33

slide-26
SLIDE 26

Random increasing trees

Another random model that produces very different shapes uses the following simple process, which generates random recursive trees: Start with the root, which is labelled 1. The n-th vertex is attached to one of the previous vertices, uniformly at random. In this way, the labels along any path that starts at the root are increasing. Clearly, there are (n − 1)! possible recursive trees of order n, and there are indeed interesting connections to permutations. The model can be modified by not choosing a parent uniformly at random, but depending on the current outdegrees (to generate, for example, binary increasing trees).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

11 / 33

slide-27
SLIDE 27

Random increasing trees

An example of an increasing tree: 7 10 13 8 11 12 3 9 6 2 4 5 1

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

12 / 33

slide-28
SLIDE 28

Convergence of the distribution

Theorem (Bhamidi/Evans/Sen, 2012)

Under fairly general assumptions on the random tree model (for example, uniformly random labelled trees and random recursive trees are covered), there exists a (model dependent) deterministic probability distribution S such that the spectral distribution of random n-vertex trees converges in distribution to S in the topology of weak convergence of probability measures on R.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

13 / 33

slide-29
SLIDE 29

Convergence of the distribution

Theorem (Bhamidi/Evans/Sen, 2012)

Under fairly general assumptions on the random tree model (for example, uniformly random labelled trees and random recursive trees are covered), there exists a (model dependent) deterministic probability distribution S such that the spectral distribution of random n-vertex trees converges in distribution to S in the topology of weak convergence of probability measures on R. Not much is known about the limiting distribution S, not even in special

  • cases. For example, it is not known whether the distribution is purely

discrete or contains a continuous component.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

13 / 33

slide-30
SLIDE 30

Single eigenvalues

In the following, we will focus on simply generated trees, or equivalently Galton-Watson trees. By Tn, we denote a random tree with n vertices from some given simply generated family (e.g. rooted labelled trees, binary trees).

Theorem

Let α be a fixed real number. The measure of {α} in the spectral distribution of Tn (i.e., the proportion of α among the eigenvalues of Tn) converges in probability to a constant cα ≥ 0 (that depends not only on α, but also the family of trees). If α is an eigenvalue of some tree in the respective family, then cα > 0.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

14 / 33

slide-31
SLIDE 31

Single eigenvalues

Theorem

Let α be a fixed real number. The measure of {α} in the spectral distribution of Tn (i.e., the proportion of α among the eigenvalues of Tn) converges in probability to a constant cα ≥ 0 (that depends not only on α, but also the family of trees). If α is an eigenvalue of some tree in the respective family, then cα > 0.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

15 / 33

slide-32
SLIDE 32

Single eigenvalues

Theorem

Let α be a fixed real number. The measure of {α} in the spectral distribution of Tn (i.e., the proportion of α among the eigenvalues of Tn) converges in probability to a constant cα ≥ 0 (that depends not only on α, but also the family of trees). If α is an eigenvalue of some tree in the respective family, then cα > 0. An analogous statement holds for Laplacian eigenvalues; in this case, cα > 0 provided that α is an eigenvalue of an augmented Laplacian matrix

  • f a tree in the respected family (obtained by adding 1 to the diagonal

entry corresponding to the root in the Laplacian matrix).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

15 / 33

slide-33
SLIDE 33

Additive functionals

Definition

An invariant F(T) defined for rooted trees T is called additive if it satisfies a recursion of the form F(T) =

k

  • i=1

F(Ti) + f(T), where T1, . . . , Tk are the branches of the tree and f(T) is a so-called toll function.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

16 / 33

slide-34
SLIDE 34

Additive functionals

Definition

An invariant F(T) defined for rooted trees T is called additive if it satisfies a recursion of the form F(T) =

k

  • i=1

F(Ti) + f(T), where T1, . . . , Tk are the branches of the tree and f(T) is a so-called toll function.

Definition

A fringe subtree of a rooted tree is a subtree consisting of a vertex v and all its descendants. The vertex v is the natural root of this subtree.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

16 / 33

slide-35
SLIDE 35

Additive functionals

For a vertex v of a rooted tree T, let Tv denote the fringe subtree rooted at v. If F is additive with toll function f, then we have F(T) =

  • v∈V (T)

f(Tv). This is easily proved by induction.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

17 / 33

slide-36
SLIDE 36

Additive functionals

For a vertex v of a rooted tree T, let Tv denote the fringe subtree rooted at v. If F is additive with toll function f, then we have F(T) =

  • v∈V (T)

f(Tv). This is easily proved by induction. Thus, if we set f(T) =

  • 1

T is isomorphic to S,

  • therwise,

then the corresponding additive functional F counts the number of

  • ccurrences of S as a fringe subtree.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

17 / 33

slide-37
SLIDE 37

Additive functionals

For us, the relevant functional is Nα(T), the multiplicity of α as an eigenvalue of T. Letting T1, . . . , Tk be the branches of T, the sum

k

  • i=1

Nα(Ti) is the multiplicity of α as an eigenvalue of the disjoint union k

i=1 Ti.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

18 / 33

slide-38
SLIDE 38

Additive functionals

For us, the relevant functional is Nα(T), the multiplicity of α as an eigenvalue of T. Letting T1, . . . , Tk be the branches of T, the sum

k

  • i=1

Nα(Ti) is the multiplicity of α as an eigenvalue of the disjoint union k

i=1 Ti.

This is also the subgraph of T obtained by removing the root. So Cauchy’s interlacing theorem immediately shows that nα(T) = Nα(T) −

k

  • i=1

Nα(Ti) ∈ {−1, 0, 1}. In other words, the toll function only takes the values 0 and ±1.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

18 / 33

slide-39
SLIDE 39

The convergence result

Convergence of the proportion Nα(Tn)/n is now a consequence of the following general theorem due to Janson:

Theorem (Janson 2012)

If the toll function f is bounded, then the associated additive functional F satisfies F(Tn) n

p

→ E(T ), where T denotes an unconditioned Galton-Watson tree.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

19 / 33

slide-40
SLIDE 40

The convergence result

Convergence of the proportion Nα(Tn)/n is now a consequence of the following general theorem due to Janson:

Theorem (Janson 2012)

If the toll function f is bounded, then the associated additive functional F satisfies F(Tn) n

p

→ E(T ), where T denotes an unconditioned Galton-Watson tree. The constant E(T ) is clearly positive in the special case “occurrences of a fixed fringe subtree”. This will be exploited in the following construction.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

19 / 33

slide-41
SLIDE 41

Forcing subtrees

Let α be an eigenvalue of some rooted tree R, and consider a tree S consisting of a root to which two (or possibly more) copies of R are attached.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

20 / 33

slide-42
SLIDE 42

Forcing subtrees

Let α be an eigenvalue of some rooted tree R, and consider a tree S consisting of a root to which two (or possibly more) copies of R are attached. Now suppose S occurs k times as a fringe subtree in a larger tree T.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

20 / 33

slide-43
SLIDE 43

Forcing subtrees

Let α be an eigenvalue of some rooted tree R, and consider a tree S consisting of a root to which two (or possibly more) copies of R are attached. Now suppose S occurs k times as a fringe subtree in a larger tree T. When the roots of these k fringe subtrees are removed, at least 2k of the remaining components are isomorphic to R. Thus the multiplicity

  • f α as an eigenvalue of the resulting forest is at least 2k.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

20 / 33

slide-44
SLIDE 44

Forcing subtrees

Let α be an eigenvalue of some rooted tree R, and consider a tree S consisting of a root to which two (or possibly more) copies of R are attached. Now suppose S occurs k times as a fringe subtree in a larger tree T. When the roots of these k fringe subtrees are removed, at least 2k of the remaining components are isomorphic to R. Thus the multiplicity

  • f α as an eigenvalue of the resulting forest is at least 2k.

By the interlacing property, the multiplicity of α as an eigenvalue of T is at least k.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

20 / 33

slide-45
SLIDE 45

Forcing subtrees

Let α be an eigenvalue of some rooted tree R, and consider a tree S consisting of a root to which two (or possibly more) copies of R are attached. Now suppose S occurs k times as a fringe subtree in a larger tree T. When the roots of these k fringe subtrees are removed, at least 2k of the remaining components are isomorphic to R. Thus the multiplicity

  • f α as an eigenvalue of the resulting forest is at least 2k.

By the interlacing property, the multiplicity of α as an eigenvalue of T is at least k. Thus the multiplicity of α is always greater than or equal to the number of

  • ccurrences of S as a fringe subtree. Positivity of cα follows.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

20 / 33

slide-46
SLIDE 46

How to compute cα

By the aforementioned result of Janson, we have cα = E(nα(T )) = P(nα(T ) = 1) − P(nα(T ) = 1). An analytic expression can be obtained by means of generating functions.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

21 / 33

slide-47
SLIDE 47

Generating functions

Recall that a family of simply generated trees is characterised by a weight function on rooted ordered trees: w(T) =

  • i≥0

wKi(T)

i

, for a given weight sequence 1 = w0, w1, w2, . . ..

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

22 / 33

slide-48
SLIDE 48

Generating functions

Recall that a family of simply generated trees is characterised by a weight function on rooted ordered trees: w(T) =

  • i≥0

wKi(T)

i

, for a given weight sequence 1 = w0, w1, w2, . . .. This can be translated to an identity for the (weighted) generating function T(x) =

  • T

w(T)x|T|, the sum being over all rooted ordered trees T. The coefficient of xn is the total weight of all trees with n vertices.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

22 / 33

slide-49
SLIDE 49

Generating functions

We have T(x) =

  • T

w(T)x|T|

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-50
SLIDE 50

Generating functions

We have T(x) =

  • T

w(T)x|T| =

  • k≥0
  • T1

· · ·

  • Tk

wkw(T1) · · · w(Tk)x|T1|+···+|Tk|+1

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-51
SLIDE 51

Generating functions

We have T(x) =

  • T

w(T)x|T| =

  • k≥0
  • T1

· · ·

  • Tk

wkw(T1) · · · w(Tk)x|T1|+···+|Tk|+1 = x

  • k≥0

wk

T1

w(T1)x|T1|k

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-52
SLIDE 52

Generating functions

We have T(x) =

  • T

w(T)x|T| =

  • k≥0
  • T1

· · ·

  • Tk

wkw(T1) · · · w(Tk)x|T1|+···+|Tk|+1 = x

  • k≥0

wk

T1

w(T1)x|T1|k = x

  • k≥0

wkT(x)k

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-53
SLIDE 53

Generating functions

We have T(x) =

  • T

w(T)x|T| =

  • k≥0
  • T1

· · ·

  • Tk

wkw(T1) · · · w(Tk)x|T1|+···+|Tk|+1 = x

  • k≥0

wk

T1

w(T1)x|T1|k = x

  • k≥0

wkT(x)k = xΦ(T(x)).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-54
SLIDE 54

Generating functions

We have T(x) =

  • T

w(T)x|T| =

  • k≥0
  • T1

· · ·

  • Tk

wkw(T1) · · · w(Tk)x|T1|+···+|Tk|+1 = x

  • k≥0

wk

T1

w(T1)x|T1|k = x

  • k≥0

wkT(x)k = xΦ(T(x)). Here, Φ(t) =

k≥0 wktk is the weight generating function associated with

the specific family. For example, Φ(t) = et for rooted labelled trees, Φ(t) =

1 1−t for rooted ordered (plane) trees, or Φ(t) = 1 + t2 for binary

trees.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

23 / 33

slide-55
SLIDE 55

Generating functions

If there exists a positive real solution τ to the equation Φ(t) = tΦ′(t), then

  • ne can show that the generating function T(x) has a square-root

singularity at ρ = τ/Φ(τ) = 1/Φ′(τ), and a lot of information (in particular, an asymptotic formula for the number/total weight of n-vertex trees) can be inferred.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

24 / 33

slide-56
SLIDE 56

Generating functions

If there exists a positive real solution τ to the equation Φ(t) = tΦ′(t), then

  • ne can show that the generating function T(x) has a square-root

singularity at ρ = τ/Φ(τ) = 1/Φ′(τ), and a lot of information (in particular, an asymptotic formula for the number/total weight of n-vertex trees) can be inferred. For us, it is mostly important that the constant cα can be expressed as cα = 1 τ

  • T

w(T)nα(T)ρ|T| = 1 τ

  • T:nα(T)=1

w(T)ρ|T| −

  • T:nα(T)=−1

w(T)ρ|T| .

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

24 / 33

slide-57
SLIDE 57

Categories of trees

In order to evaluate the expression, we will categorise rooted trees. For a rooted tree T with root r, let Ψ(T, z) = det(zI − A(T)) be the characteristic polynomial of (the adjacency matrix of) T, and let Ψr(T, z) be the characteristic polynomial of (the adjacency matrix of) the forest that results when the root is removed.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

25 / 33

slide-58
SLIDE 58

Categories of trees

In order to evaluate the expression, we will categorise rooted trees. For a rooted tree T with root r, let Ψ(T, z) = det(zI − A(T)) be the characteristic polynomial of (the adjacency matrix of) T, and let Ψr(T, z) be the characteristic polynomial of (the adjacency matrix of) the forest that results when the root is removed. Consider the quotient Q(T, z) = Ψr(T, z) Ψ(T, z) .

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

25 / 33

slide-59
SLIDE 59

Categories of trees

In order to evaluate the expression, we will categorise rooted trees. For a rooted tree T with root r, let Ψ(T, z) = det(zI − A(T)) be the characteristic polynomial of (the adjacency matrix of) T, and let Ψr(T, z) be the characteristic polynomial of (the adjacency matrix of) the forest that results when the root is removed. Consider the quotient Q(T, z) = Ψr(T, z) Ψ(T, z) . We note that nα(T) = 1 if and only if Q(T, z) has a pole at z = α,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

25 / 33

slide-60
SLIDE 60

Categories of trees

In order to evaluate the expression, we will categorise rooted trees. For a rooted tree T with root r, let Ψ(T, z) = det(zI − A(T)) be the characteristic polynomial of (the adjacency matrix of) T, and let Ψr(T, z) be the characteristic polynomial of (the adjacency matrix of) the forest that results when the root is removed. Consider the quotient Q(T, z) = Ψr(T, z) Ψ(T, z) . We note that nα(T) = 1 if and only if Q(T, z) has a pole at z = α, nα(T) = −1 if and only if Q(T, z) has a zero at z = α,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

25 / 33

slide-61
SLIDE 61

Categories of trees

In order to evaluate the expression, we will categorise rooted trees. For a rooted tree T with root r, let Ψ(T, z) = det(zI − A(T)) be the characteristic polynomial of (the adjacency matrix of) T, and let Ψr(T, z) be the characteristic polynomial of (the adjacency matrix of) the forest that results when the root is removed. Consider the quotient Q(T, z) = Ψr(T, z) Ψ(T, z) . We note that nα(T) = 1 if and only if Q(T, z) has a pole at z = α, nα(T) = −1 if and only if Q(T, z) has a zero at z = α, nα(T) = 0 if and only if limz→α Q(T, z) is finite and not zero.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

25 / 33

slide-62
SLIDE 62

Recursions

Let T be a rooted tree with root r, and let T1, . . . , Tk be its branches with roots r1, . . . , rk. It is not too hard to show that Ψr(T, z) Ψ(T, z) = 1 z − k

j=1 Ψrj (Tj,z) Ψr(Tj,z)

,

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

26 / 33

slide-63
SLIDE 63

Recursions

Let T be a rooted tree with root r, and let T1, . . . , Tk be its branches with roots r1, . . . , rk. It is not too hard to show that Ψr(T, z) Ψ(T, z) = 1 z − k

j=1 Ψrj (Tj,z) Ψr(Tj,z)

,

  • r

Q(T, z) = 1 z − k

j=1 Q(Tj, z)

.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

26 / 33

slide-64
SLIDE 64

Recursions

Let T be a rooted tree with root r, and let T1, . . . , Tk be its branches with roots r1, . . . , rk. It is not too hard to show that Ψr(T, z) Ψ(T, z) = 1 z − k

j=1 Ψrj (Tj,z) Ψr(Tj,z)

,

  • r

Q(T, z) = 1 z − k

j=1 Q(Tj, z)

. We find that Q(T, z) has a zero at z = α if and only if at least one of the Q(Tj, z) has a pole (one direction is clear, for the converse one needs to ensure that there cannot be cancellations in the denominator), Q(T, z) has a pole at z = α if and only if k

j=1 Q(Tj, α) = α,

Q(T, α) = β if and only if k

j=1 Q(Tj, α) = α − 1/β.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

26 / 33

slide-65
SLIDE 65

Functional equations

Thus we have a recursive description of trees according to their “type” (zero, pole or other value). For β ∈ R ∪ {∞}, we set Hβ(x) =

  • T:Q(T,α)=β

w(T)x|T|.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

27 / 33

slide-66
SLIDE 66

Functional equations

Thus we have a recursive description of trees according to their “type” (zero, pole or other value). For β ∈ R ∪ {∞}, we set Hβ(x) =

  • T:Q(T,α)=β

w(T)x|T|. Clearly, we have T(x) =

  • β∈R∪{∞}

Hβ(x).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

27 / 33

slide-67
SLIDE 67

Functional equations

Thus we have a recursive description of trees according to their “type” (zero, pole or other value). For β ∈ R ∪ {∞}, we set Hβ(x) =

  • T:Q(T,α)=β

w(T)x|T|. Clearly, we have T(x) =

  • β∈R∪{∞}

Hβ(x). Moreover, cα = 1 τ

  • T:nα(T)=1

w(T)ρ|T| −

  • T:nα(T)=−1

w(T)ρ|T| = 1 τ (H∞(ρ) − H0(ρ)).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

27 / 33

slide-68
SLIDE 68

Functional equations

The recursive characterisation of Q(T, z) can be turned into the following system of functional equations: H0(x) = xΦ(T(x)) − xΦ(T(x) − H∞(x)), Hβ(x) = [yα−1/β]Φ

γ∈R

Hγ(x)yγ . In the second equation, y is just a formal variable, and there are actually infinitely many equations (one for each β).

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

28 / 33

slide-69
SLIDE 69

Numerical values

The infinite system of equations can actually be solved numerically, and yields results such as the following:

Theorem

The value of the constant c1 (limiting proportion of the eigenvalue 1) is ≈ 0.0213 for random (rooted) labelled trees, ≈ 0.0255 for random plane trees, ≈ 0.0141 for random pruned binary trees.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

29 / 33

slide-70
SLIDE 70

Numerical values

The infinite system of equations can actually be solved numerically, and yields results such as the following:

Theorem

The value of the constant c1 (limiting proportion of the eigenvalue 1) is ≈ 0.0213 for random (rooted) labelled trees, ≈ 0.0255 for random plane trees, ≈ 0.0141 for random pruned binary trees. In other words, the eigenvalue 1 makes up about 2.13% (2.55%, 1.41%) of the spectrum of large random labelled (plane, pruned binary) trees, respectively.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

29 / 33

slide-71
SLIDE 71

Central limit theorem

Making use of ideas of Drmota, Gittenberger and Morgenbesser on infinite systems of functional equations, one also obtains a central limit theorem, i.e. a result of the following form:

Theorem

For a given simply generated family of trees and a real number α that

  • ccurs as eigenvalue of at least one tree in the family, there exist constants

cα > 0 and bα > 0 such that mean and variance of Nα(Tn) are equal to µn = cαn + o(n) and σ2

n = bαn + o(n) respectively, and the normalised

random variable Nα(Tn) − µn

  • σ2

n

converges weakly to a standard normal distribution.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

30 / 33

slide-72
SLIDE 72

The eigenvalue 0

The eigenvalue 0 is special in some sense: n0(T) = ±1 for all trees T (there are no trees with n0(T) = 0), so there are only two possibilities: Q(T, z) has either a pole or a zero at z = 0.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

31 / 33

slide-73
SLIDE 73

The eigenvalue 0

The eigenvalue 0 is special in some sense: n0(T) = ±1 for all trees T (there are no trees with n0(T) = 0), so there are only two possibilities: Q(T, z) has either a pole or a zero at z = 0. There is also a combinatorial interpretation: since the characteristic polynomial of the adjacency matrix coincides with the matching polynomial for all forests, N0(T) = n − 2µ(T), where n is the number of vertices and µ(T) the matching number.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

31 / 33

slide-74
SLIDE 74

The eigenvalue 0

The system of functional equations reduces to a 2 × 2-system, and explicit values of c0 can be found for some families of trees: c0 = 2W(1) − 1 ≈ 0.134 for random (rooted) labelled trees (where W is the Lambert W-function), c0 = √ 5 − 2 ≈ 0.236 for random plane trees, c0 = 7 − 4 √ 3 ≈ 0.0717 for random pruned binary trees.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

32 / 33

slide-75
SLIDE 75

More trees and spectra

Similar results can also be obtained for other classes of trees and other types of spectra, for instance: a central limit theorem for the multiplicity of Laplacian eigenvalues in simply generated trees (based on very similar functional equations); a central limit theorem for the multiplicity Nα in recursive trees; the constant c0 in this case is 1 1 + log t 1 − log t dt ≈ 0.193; a “forcing subtrees” type result for the distance spectrum, showing that the proportion of certain eigenvalues in the distance spectrum of large random rooted trees can be bounded below by positive constants.

Eigenvalue distribution in random trees

  • S. Wagner, Stellenbosch University

33 / 33