4 Loewner order

Definition 4.1 (Loewner order). For A, B matrices, write A≼B if B−A is positive semidefinite. In particular, A≽0 if and only if A is positive semidefinite.

qA(f)=⟨f,Af⟩⟨f,f⟩=0.

A≼B if and only if ∀⁡f≠0,

⟨f,Af⟩⟨f,f⟩≤⟨f,Bf⟩⟨f,f⟩.

⟨f,Af⟩≤⟨f,Bf⟩.

This is indeed an order: A≼B≼C implies A≼C

⟨f,Af⟩≤⟨f,Bf⟩≤⟨f,Cf⟩.

If A≼B, then λk(A)≤λk(B).

λk(A)=min⁡ dim ⁡⁡W=kmax⁡f∈Wf≠0⟨f,Af⟩⟨f,f⟩.

A≼B if and only if A+C≼B+C.

If G is a graph, then LG≼0.

Definition 4.2 (eps-approximation). G is an 𝜀-approximation of H if

(1−𝜀)LH≼LG≼(1+𝜀)LH.

Lemma 4.3. Given the definition ∥M∥=max⁡f≠0∥Mf∥∥f∥ (for M symmetric), we have ∥M∥=max⁡{|λk(M)|}.

Proof. f=∑⁡iαiφi, Mφi=λiφi, ∥f∥2=∑⁡iαi2,

∥Mf∥2=∥∑iαiλiφi∥2=∑iαi2λi2.
∥Mf∥∥f∥=∑iαi2λi2∑iαi2≤max⁡k|λk|.□

Lemma 4.4. Assuming: - G is an 𝜀-approximation of H Then: ∥LG−LH∥≤𝜀.

Proof. −𝜀L~H≼L~G−L~H≼𝜀L~H. λk(L~G−L~H)≤λk(L~H)≤2𝜀. □

Definition ((d,eps)-expander). G=(V,E), |V|=n is a (d,𝜀)-expander if G is an 𝜀-approximation of dnKn.

Equivalent:

(1−𝜀)dnLKn≼LG≼(1+𝜀)dnLKn.

LKn=(n−1)I−AKn=nI−J, where J is the all ones matrix (J(x,y)=1). If f⊥q, then Jf=⟨f,1⟩f=0. In this case, LKnf=nIf−Jf=nf.

So λk(LKn)=n for k≥0.

(1−𝜀)dn(nI−J)≼dI−AG≼(1+𝜀)dn(nI−J)
−𝜀2n(nI−J)≼dI−AG−dn(nI−J)≼𝜀dn(nI−J)
−𝜀(dI−dnJ)≼dnJ−AG≼𝜀(dI−dnJ)

For f⊥1, −𝜀d⟨f,f⟩≤−⟨f,Af⟩≤𝜀d⟨f,f⟩.

So G is a (d,𝜀)-expander if and only if

|λk(AG)|≤𝜀d

for all 1≤k≤n−1.

PIC

Lemma 4.5 (Expander Mixing Lemma). Assuming that:

Then ∀⁡S,T⊆V,
|e(S,T)−dn|S||T||≤𝜀dn|S||T||Sc||Tc|.

Proposition 4.6. Assuming that:

  • G a d-regular graph on n vertices

  • λ,𝜀>0, 𝜀d=λ

Then the following are equivalent:
  • (i) G is a (n,d,λ)-graph
  • (ii) λk(AG)∋[−λ,λ], for 1≤k≤n−1
  • (iii) λk(LG)∈[d−λ,d+λ] for 2≤k≤n
  • (iv) λk(L~G)∈[1−λd,1+λd]=[1−𝜀,1+𝜀] for 1≤k≤n
  • (v) (1−𝜀)dnLKn≼LG≼(1+𝜀)dnLKn
  • (vi) G is (d,𝜀)-expander (also (d,λd)-expander)
  • (vii) ∥LG−dnLKn∥≤𝜀d=λ
  • (viii) ∥AG−dnJ∥≤𝜀d=λ

Lemma 4.7 (Expander Mixing Lemma). Assuming that:

  • G=(V,E) an (n,d,λ)-graph

  • S,T⊆V (and define e(S,T)=∑⁡x∈S∑y∈T𝟙xy∈E)

Then |e(S,T)−dn|S||T||≤λn|S||Sc||T||Tc|≤λ|S||T|

Proof. ⟨𝟙S,LG𝟙T⟩=⟨𝟙S,(dI−AG)𝟙T⟩, ⟨𝟙S,dI𝟙T⟩=d|S∩T|.

⟨𝟙S,AG𝟙T⟩=∑x∑y𝟙S(x)AG(x,y)𝟙T(y)=∑x∑y𝟙xy∈E=e(S,T)

dnLKn=dI−dnJ.

⟨𝟙S,dnJ𝟙T⟩=dn⟨𝟙S,J𝟙T⟩=dn∑x∑yJ(x,y)=dn|S||T|.

|e(S,T)−dn|S||T||=|⟨𝟙S,(dnLKn−LG)𝟙T⟩|≤∥𝟙S∥∥(dnLKn−LG)𝟙T∥≤∥𝟙S∥∥(dnLKn−LG)∥∥𝟙T∥≤λ∥𝟙S∥∥𝟙T∥=λ|S||T| To get the better bound, we should consider functions which are perpendicular to 1: balanced function fS=𝟙S−|S|n1.

LGfS=LG(𝟙S−α)=L𝟙S.

|e(S,T)−dn|S||T||=|⟨fS,(dnLKn−LG)fT⟩|≤λ∥fS∥∥fT∥

∥fS∥2=|S|(1−|S|n)2+(n−|S|)(−|S|n)2=|S|(1−2|S|n+|S|2n2)+(n−|S|)|S|2n2=|S|−2|S|2n+|S|3n2+|S|2n−|S|3n2=|S|−|S|2n=n|S|−|S|2n=|S||Sc|n So
|e(S,T)−dn|S||T||≤|S||Sc||T||Tc|.□

If G is an (n,d,λ)-graph and I⊆V an independent set, then

0=e(I,I)≥dn|I|2−λn|I||Ic|.

λ|I||Ic|≥d|I|3.

|I|≤λd|Ic|=λd(n−|I|).
(1+λd)|I|≤λdn
|I|≤λd(1+λd)n=λd+λn.

Hoffman bound: α(G)≤λd+λn.

Fix d. How small can λ be such that there is an infinite family Gn of (n,d,λ)-graphs?

Note AG2 has d in each entry of the diagonal, so

Tr⁡AG2=dn=∑iλi(AG2)=∑i(λi(AG))2.

So dn≤d2+(n−1)λ2.

So (n−1)λ2≥dn−d2=d(n−d), so

λ2≥d(n−d)n−1=d(n−1−(d−1)n−1)=d(1−d−1n−1).

λ≥d(1−o(1)) as n→∞.

Alon-Boppana Theorem: λ≥2d−1−o(1).

Claim: There exist families of (n,d,λ)-graphs with λ=2d−1. They are called Ramanujan graphs. We will probably not prove existence of these.

Call 𝜀-Ramanujan if λ≥2d−1+𝜀.

Theorem 4.8 (Friedman). Assuming: - 𝜀>0, n→∞ Then:

ℙ(random d-regular graph on n vertices is 𝜀-Ramanujan)→1.

Maxcut in (n,d,λ)-graph:

e(S,Sc)≤dn|S||Sc|+λn|S||Sc|≤(dn+λn)n24=dn4+λn4≤e(G)2+λn4

Diameter, vertex expansion. S⊆V, ∂S={x∈Sc:x∼S}.

|∂S||S|≥e(S,Sc)d|S|=Φ(S)≥λ2(L~G)2≥(1−λd)2

|∂S|≥(1−λ∕d2)|S|.

Exercise: Diameter O(log⁡n).

Vertex expansion

(n,d,λ)-graph G=(V,E). Have for all S,T⊆V,

|e(S,T)−dn|S||T||≤λn|S||Sc||T||Tc|.

Φv(S)=|∂S||S|.

PIC

|∂S|≥e(S,Sc)dΦv(S)≥e(S,Sc)d|S|=Φe(S) If S⊆V, |S|≤n∕2, then ΦV(S)≥Φe(G)≥λ2(L~G)2Φv(S)≥(1−λd)2|S∪∂S|≥(1+12−λ2d)

n→∞, d, λ fixed. λ≤d−δ, δ>0. |S∪∂S|≥(1+κ)|S| for some κ>0. Hence

|B(x,r)|≥(1+κ)r

if |B(x,r−1)|≤n∕2.

diam⁡G=Oλd( log ⁡⁡n) (by considering ball around start and end).

Why aren’t Cayley graphs of abelian groups expanders?

Let Γ be an abelian group and S⊆Γ a set of generators of size d. Let |Γ|=n→∞. Let G=Cay⁡(Γ,S). Then

|B(x,r)|≤(2r+1)d.

This is not exponential in r, so the Cayley graph can’t be an expander.

Theorem 4.9 (Alon-Boppana). Assuming that:

Then as n→∞
λ≥2d−1−O(1).

Proof 1. Pick edge st, pick r∈ℤ≥0.

PIC

φ:V→ℝ,

φ(x)={(d−1)−i∕2if x∈Vi, i≤r0if x∈Vi, i>r
qL(φ)=⟨φ,Lφ⟩⟨φ,φ⟩⟨φ,φ⟩=∑i=0r|Vi|(d−1)i.

⟨φ,Lφ⟩=∑x∼y(φ(x)−φ(y))2=∑i=0r−1e(Vi,Vi+1)(1(d−1)i∕2−1(d−1)(i+1)∕2)2+e(Vr,Vr+1)(d−1)r=∑i=0r−1e(Vi,Vi+1)(d−1)i(1−1d−1)2+e(Vr,Vr+1)(d−1)r e(Vi,Vi+1)≤(d−1)|Vi|. ⟨φ,Lφ⟩≤∑i=0r−1|Vi|(d−1)i(d−1−1)2+|Vr|(d−1)r(d−1)(d−1−1)2=(d−1)−2d−1+1=d−2d−1⟨φ,Lφ⟩≤(d−2d−1)∑i=0r−1|Vi|(d−1)i+(d−1)|Vr|(d−1)r|Vr|≤(d−1)|Vr−1|≤⋯≤(d−1)r−i|Vi|⟹|Vr|(d−1)r≤1r∑i=0r|Vi|(d−1)i⟨φ,Lφ⟩≤(d−2d−1)∑i=0r|Vi|(d−1)i+(2d−1)|Vr|(d−1)r≤(d−2d−1+2d−1−1r+1)⟨φ,φ⟩⟹λ2(L)=min⁡ dim ⁡⁡W=2max⁡f∈Wf≠0⟨f,Lf⟩⟨f,f⟩

Suppose G has 2 deges at distance >2r+2.

PIC

f, f′ as above. ⟨f,f′⟩=0.

QL(αf+βf′)=QL(αf)+QL(βf′).

Let W=span⁡{f,f′}.

λ2(L)≤max⁡(α,β)≠(0,0)QL(αf+βf′)∥αf+βf′∥≤α2QL(f)+β2QL(f′)α2∥f∥+β2∥f′∥≤d−2d−1+2d−1−1r+1□

Corollary 4.10. For all d, there are finitely many (n,d,λ)-graphs with λ<2d−1.

r=clog⁡n.

λn−1(AG)≥2d−1−Od(1log⁡n).

For Alon-Boppana:

Proof 2. Tr⁡A2k=∑⁡xA2k(x,x). Note that

#{closed walks of length 2k in G starting from x}.

is at least

#{closed walks of length 2k in ∏d starting from 0}

(∏⁡d is an infinite d-regular tree). The latter is at least

(d−1)k1k+12kk≈(2d−1)2k+o(1)22k.
≤∑i=1nλi2k≤d2k+(n−1)λ2k.

Exercise: finish details. □

min⁡{|λ1(A)|λn−1(A)}≥⋯

4.1 One-sided expanders

Goal is to find: Gn graphs with n vertices, d-regular such that λn−1(AG)≤λ<d.

Reminder: (Friedman) If G∼Gn,d uniform random d-regular graph then

ℙ(Gn is (n,d,2d−1+𝜀)-graph)→1.
Φ(Gn)≥λ2(L~G)2=λ2(LG)2d=(d−λn−1(AG))2d≥12−𝜀2>0.

Theorem 4.11. Assuming that:

  • d=2l large enough

Then there is 𝜀=𝜀d>0 such that there are d-regular graphs (or multi-graphs) Gn on n vertices with Φ(Gn)≥𝜀 for all sufficiently large n.

Graphs Gn are as follows:

Lemma 4.12. There is c=cd>0 such that

ℙ(Gn is d-regular, i.e. Gn is simple)≥c−o(1).

Lemma 4.13. If Gn is d-regular, then there exists 𝜀=𝜀d>0 such that

ℙ(Φ(Gn)<𝜀)→0.

Proof. If Φ(Gn)<𝜀 then there exists F⊆[n], r=|F|≤n2, there exists F′, F⊆F′ and e(F,F′)=d|F| and |F′|=r+r′, r′=⌊𝜀r⌋.

PIC

ℙ(Φ(Gn)<𝜀)<∑1≤r≤n2∑F⊆F′⊆[n]|F|=r|F′|=r+r′ℙ(πi(F)⊆F′ for all i∈[l]).

ℙ(π1(F)⊆F′)=r+r′rnr.

ℙ≤∑1≤r≤n2n!r!r′!(n−r−r′)!(r+r′rnr)l.

Fact 1: akbk≤(ab)k if a≤b.

Proof: It is equivalent to

ba⋅b(a−1)⋯b(a−k+1)≤ab⋅a(b−1)⋯a(b−k+1).

Compare the product term by term.

Fact 2: n!≥(ne)n. Proof:

en=∑k≥0nkk!≥nnn!.

Using these:

n!r!r′!(n−r−r′)!≤nr+r′r!r′!=nr+r′(r+r′)!r+r′r≤(2e)r+r′(ne+e′)r+r′.

So

ℙ≤∑1≤r≤n2(2e)r+r′(nr+r′)r+r′(r+r′n)lr≤∑1≤r≤n2(2e)2r((1+𝜀)rn)r(l−1−𝜀)

Decompose as ∑⁡1≤r≤n∕2=∑1≤r≤K+∑K<r≤n∕2=S1+S2. Choose 𝜀>0 small so that γ=1+𝜀2 is small. Choose l large so that γl−1−𝜀<12(2e)2.

S2≤∑K<r≤n2(2e)2r((1+𝜀)2)r(l−1−𝜀)≤∑K<r≤n2(12)r≤12K

Now S1.

S1≤∑1≤r≤K(2e)2r((1+𝜀)rn)r(l−1−𝜀)≤(2Kn)l−2∑1≤r≤K(2e)2r≤cK(Kn)l

Now let K=log⁡ log ⁡ log ⁡n, so this is ≤O((log⁡ log ⁡n)n−l( log ⁡ log ⁡ log ⁡n)l)→0. Also, S1≤12K term goes to 0. □

For the lemma about ℙ(d-regular)≥c−o(1):

PIC

These are the bad things. Use Bonferroni inequalities (the partial sum of inclusion exclusion principle inequalities).

(n,d,λ)-graph

n: number of vertices, d-regular, λk(A)∈[−λ,λ], k≠n.

Ramanujan graph: (n,d,2d−1)-graph.

Alon-Boppana: For every 𝜀>0 fixed d≥3, there are finitely many (n,d,2d−1−𝜀)-graphs.

Bipartite: (n,d,λ)-graph, n vertices, d-regular, λk(A)∈[−λ,λ] for k≠n,1.

Theorem 4.14 (Lubotsky-Phillips-Sanak / Margoulis). Assuming that:

  • p prime

  • d=p+1

  • n arbitrarily large

Then there exists a (n,d,2d−1)-graph.

Goal:

Theorem 4.15 (Marcus, Spielman, Srivastava). Assuming that:

  • d≥3

Then there are (n,d,2d−1)-bipartite graphs for arbitrarily large n.

Strategy: Bilu and Linial. Lifts of graphs:

Definition 4.16 (2-lift). A 2-lift of a graph G=(V,E) is a graph Ĝ=(V^,Ê) with

  • x∈V ⟹ x0,x1∈V^.

  • xy∈E ⟹ either x0y0,x1y1∈Ê or x0y1,x1y0∈Ê.

(and no other vertices or edges).

PIC

Definition 4.17 (Signing). S:V×V→ℝ is a signing of G if

S(x,y)={±1if A(x,y)=10if A(x,y)=0

and S(x,y)=S(y,x) (symmetric). So can think of S as a function E→{±1}.

AS+(x,y)=A(x,y)𝟙S(x,y)=1. AS−(x,y)=A(x,y)𝟙S(x,y)=−1. Have A=AS++AS−, S=AS+−AS−.

Lemma 4.18. The eigenvalues of Ĝ are the eigenvalues of AG (old) together with the eigenvalues of S (new).

Proof.

AĜS=(AS+ASS−AS−ASS+)

Let Aφ=λφ. Then

AĜS(φφ)=(AS+φ+AS−φAS−φ+AS+φ)=(AφAφ)=(λφλφ)=λ(φφ).

Let Sη=μη. Now

AĜS(η−η)=(AS+η−AS−ηAS−η−AS+η)=(Sη−Sη)=(μη−μη)=μ(η−η).□

Conjecture 4.19 (Bilu, Linial). If G is d-regular, then there exists a signing S:E(G)→{±1} whose eigenvalues are in [−2d−1,2d−1].

Theorem 4.20 (Bilu, Linial). Can find signings S with eigenvalues λ satisfying |λ|=O(d(log⁡d)3).

Theorem 4.21 (Marcusm, Spielman, Srivastava). Assuming that:

  • G is a d-regular graph

Then there exists a signing with eigenvalues λ with λ≤2d−1.

Theorem 4.21 implies Theorem 4.15.

  • Start with Kd,d.

  • Keep applying Theorem 4.21 to find signing.

  • Build lift with that signing.

  • 2-lift of bipartite graph is bipartite.

  • Spectrum of adjacency matrix of bipartite graph is symmetric around 0. □

Notation. For π∈Sym⁡(X), let |π| be the number of inversions.

Theorem 4.21: S∼Unif⁡({±1}E(G)).

𝔼Sdet⁡(xI−S)=𝔼S(∑π∈Sym⁡(X)(−1)|π|∏y∈V(xI−S)(y,π(y)))=∑k=0nxn−k(−1)k∑T⊆V|T|=k∑π∈Sym⁡(T)𝔼S((−1)|π|∏y∈T(xI−S)(y,π(y)))

For xy∈E: 𝔼S(x,y)2k+1=0, 𝔼S(x,y)2k=1.

PIC

=∑k=0k evennxn−k(−1)k∑M matchingof size ≤k2(−1)k∕2…=∑k≥0xn−2k(−1)kMk(G)=μG(x) where Mk(G) is the number of matchings of size k in G.

μG(x) is the matching polynomial of G.

Heilman-Lieb-Godsil:

Last time, when working on det⁡(xI−As), we showed:

Theorem 4.22 (Godsil-Gutman, 80s).

𝔼S∼{±1}det⁡(xI−As)=μG(x),

where

μG(x)=∑k≥0xn−2k(−1)kmk(G),

where mk(G) is the number of matchings in G with k edges.

Fact: μG(x) is real rooted.

Theorem 4.23 (Heilman-Lieb, 72). Assuming that:

  • G is d-regular

Then
maxroot⁡μG(x)≤2d−1.

If only we could say that maxroot⁡𝔼Sdet⁡(xI−As) is an average of maxroot⁡det⁡(xI−Ah), h∈{±1}E. Hopelessly false:

PIC

Definition 4.24 (Interlacing). Let f be a real rooted polynomial of degree n with roots α1,…,αn and g a real rooted polynomial of degree n−1 with roots β1,…,βn−1. We say g interlaces f if

α1≤β1≤α2≤⋯≤βn−2≤αn−1≤βn−1≤αn.

PIC

Example. If f is real rooted, then so is f′, and also f′ interlaces f.

Definition 4.25 (Common interlacement). We say that real rooted polynomials f, g of degree n have a common interlacement if there is a polynomial h of degree n+1 such that f and g both interlace h. In other words, if the roots of f and g are α1≤⋯αn and β1≤⋯≤βn respectively, then they have a common interlacement if and only if there are some γ0≤⋯≤γn such that

γ0≤α1,β1≤γ1≤α2,β2≤γ2≤⋯≤γn−1≤αn,βn≤γn.

Theorem 4.26 (Fell 80). Assuming that:

  • f, g are real rooted

  • both degree n, monic

Then f and g have a common interlacing if and only if every convex combination of f and g are real rooted.

PIC

PIC

Proof. Assume f and g without repeated or common roots (general case requires more work and details checking).

  • ⇐ Let ht(x)=tf(x)+(t−1)g(x) for t∈(0,1). We know these are real rooted, and waht to show that f, g have a common interlacement. Let λi(t) be the i-th root of ht(x). λi(t) includes (λi(0),λi(1))⊆ℝ.

    Claim: (λI(0),λI(1)) are disjoint.

    Suppose not. Let λk be a root of f and λk∈(λI(0),λI(1)). So there exists t∈(0,1) such that λI(t)=λk.

    ht(λk)=0=tf(λk)+(1−t)g(λk),

    so g(λk)=0, contradiction as we assumed no common roots.

  • ⇒

    PIC

    The dots represent the roots of the polynomial that interlaces them. By monicness, we get the blue and green +s and −s. Then get the orange ones, and use intermediate value theorem. □

Corollary 4.27. Assuming that:

  • f, g are real rooted of degree n

  • have a common interlacing

Then for all t∈[0,1]
min⁡{maxroot⁡(f),maxroot⁡(g)}≤maxroot⁡(tf(x)+(1−t)g(x)).

Approach to Theorem 4.21:

f(x)=𝔼S∈{±1}Edet⁡(xI−As), |E|=m. For h1,…,hk∈{±1}={+,−}, then

fh1,…,hk(x)=𝔼sk+1,…,sm∈{±1}(det⁡(xI−As)|s1=h1,…⁡,sk=hk).

“Method of interlacing families of polynomials”

PIC

Theorem 4.28 (Marcus, Spielman, Srivastava (real rooted)). For every ρ1,…,ρm∈[−1,1] the following polynomial is real rooted

χρ1,…,ρm(x)=∑S∈{±1}mdet⁡(xI−As)∏J:sJ=1(1+ρJ2)∏J:SJ=−1(1−ρJ2).

Theorem 4.28 implies Theorem 4.21:

fh1,…,hk(x)=χh1,…,hk,0,0,…,0(x).
μG(x)=χ0,…,0(x).
tfh1,…,hk,1(x)+(1−t)fh1,…,hk,−1(x)=χh1,…,hk,2t−1,0,…,0(x).

Theorem 4.29 (Marcus, Spielman, Srivastava (matrix)). Assuming that:

  • r1,…,rm∈ℝn are independent random vectors, whose support has ≤2 points

Then
𝔼det⁡(xI−∑j=1mrjrj⊤⁡)

is real rooted.

Theorem 4.29 implies Theorem 4.28:

χρ1,…,ρm(x) is

𝔼det⁡(xI−As).

PIC

𝔼det⁡(xI−As) is real rooted if and only if 𝔼det⁡((x+d)I−As) is real rooted. This equals

𝔼det⁡(xI+(dI−As)).

Note

dI−AS=∑x∼y(ex+S(x,y)ey)(ex+S(x,y)ey)⊤⁡.

Done.

Definition 4.30 (Real stable). A polynomial p(z1,…,zn)∈ℝ[z1,…,zn] is real stable if p(z1,…,zn)≠0 for all (z1,…,zn)∈Hn, H={z∈ℂ|Im⁡z>0}.

Notice that for a single variable, p(z) is real stable if and only if real rooted.

Example. 1−z1z2.

Proposition 4.31 (Stable iff real rooted). p(z1,…,zn) is real stable if and only if for all ℝ≥0n, b∈ℝn, the polynomial t↦p(at+b) is real rooted.

Proof.

  • ⇒ If p not stable, then there exists z0,…,zn∈H with p(z1,…,zn)=0, zj=bj+iaj, aj>0. Then p(ai+b)=0, t=i.
  • ⇐ Suppose ∃⁡a∈ℝ≥0n, b∈ℝn, t∈ℂ∖ℝ with p(at+b)=0. Assume Im⁡t>0. zj+ajt+bj∈H. p(z1,…,zn)=p(at+b)=0.

    PIC

    □

Proposition 4.32 (Stability of real stable). Assuming that:

  • p(z1,…,zn) is real stable

Then the following are also:
  • p(zσ(1),…,zσ(n)), σ∈Sn.

  • p(az1,z2,…,zn), a∈ℝ≥0.

  • p(z2,z2,z3,…,zn).

  • p(c,z2,…,zn), c∈H∪ℝ.

  • z1d1p(−1z1,z2,…,zn).

  • ∂z1p(z1,…,zn).

  • MAP⁡(p(z1,…⁡,zn)). MAP⁡(z1z2+z2z32+2z1z4+z2z32)=z1z2+2z1z4.

  • If p,q are real stable, then so is pq.

Proposition 4.33 (Mother of all stable polynomials). Assuming that:

  • B∈ℝd×d symmetric matrix

  • A1,…,Am∈ℝd×d positive semi definite

Then det⁡(B+z1A1+⋯+znAn) is real stable.

Proof. Let a∈ℝ≥0n, b∈ℝn. Assume positive definite instead of positive semi definite.

det⁡(B+∑j=1n(ajt+bj)Aj)=det⁡((∑j=1najAj)⏟=:At+(∑j=1nbjAj+B)⏟=:C)

A is positive definite, and C is symmetric. So

=det⁡(At+C)=det⁡(A12(tI+A−12CA−12)A12)=det⁡Adet⁡(tI+A−12CA−12)

The roots are eigenvalues of a real symmetric matrix, hence real. □

Proof of Theorem 4.29. Say rj equals rj+ with probability pj+, and equals rj− with probability pj−. We will use Cauchy Binet:

det⁡(AB)=∑det⁡(AS×[n])det⁡(B[n]×C)TODO

p(z1,…,zn)=𝔼det⁡(∑j=1nzjrjrj⊤⁡)=𝔼∑S⊆[m]|S|=ndet⁡(∑j∈Szjrjrj⊤⁡)=𝔼∑S⊆[m]|S|=n(∏j∈Szj)det⁡(∑j∈Srjrj⊤⁡)
det⁡(𝔼∑j=1nzjrjrj⊤⁡)=det⁡(∑j=1nzj𝔼rjrj⊤⁡)

is real stable since 𝔼rjrj⊤⁡≥0 is a covariance matrix.

=det⁡(∑j=1mzj(pj+rj+(rj+)⊤⁡+pj−rj−(rj−)⊤⁡))=∑S⊆[m]|S|=n(∏j∈Szj)∑n∈{±1}Sdet⁡(pjh(j)rjh(j)(rjh(j))⊤⁡)+“terms with squares”

So the earlier thing is real stable. So

𝔼det⁡(∑j=1mzjrjrj⊤⁡)

is real stable. So

𝔼det⁡(∑j=1mzjrjrj⊤⁡+∑j=1mxjejej⊤⁡)

is real stable. Take zj=−1, xj=x. Then

𝔼det⁡(xI−∑j=1mrjrj⊤⁡)

is real stable. □

Theorem 4.34 (Godsil). Assuming that:

  • G is d-regular

  • λ1≤⋯≤λn are the roots of μG(x)

Then
∑k=1nλkl=∑aWal(G),

where Wal(G) is the number of closed walks of length l from a in the path tree Ta(G) of G.

Definition 4.35 (Path tree). Example: Ta(G)

PIC

Proof of ?? implies Theorem 4.23. Wal(G)≤Wal(∏⁡d)≤(2d−1+ol(1))l.

λnl≤∑⁡λkl=∑aWal(G)≤n(2d−1+ol(1))l.

λn≤n1l(2d−1+ol(1)).

Take l→∞. □

Proof of ??. μG′(x)=∑⁡aμG−a(x) Hence

∑a∑k≥0xn−1−2k(−1)kmk(G−a)=∑k≥0(n−2k)xn−1−2k(−1)kmk(G)=μG′(x)

μG(x)=xμG−a(x)−∑b∼aμG−a−b(x)mk(G)=mk(G−a)+∑b∼amk−1(G−a−b)xn−2k(−1)kmk(G)=x⋅xn−1−2k(−1)kmk(G−a)−∑b∼axn−2−2(k−1)(−1)k−1mk−1(G−a−b)∑aμG−a(x)μG(x)=μG′(x)μG(x)=∑j=1n1x−λj=1x∑j=1n11−λjx=1x∑j=1n∑l≥0λjlxl=1x∑l≥0∑jλjlxl Claim: μG−a(x)μG(x)=1x∑⁡l≥0Wal(G)xl. μG−a(x)μG(x)=μG−a(x)xμG−a(x)−∑b∼aμG−a−b(x)=1x11−1x∑b∼aμG−a−b(x)μG−a(x)=1x⋅∑k≥0(1x∑b∼aμG−a−b(x)μG−a(x))k=1x∑k≥01xk(∑b∼a1x∑l≥0Wbl(G−a)xl)k=1x∑k≥01x2k(∑b1∼a∑l1≥0WGl1(G−a)xl1)⋯(∑bk∼a∑lk≥0Wbklk(G−a)xlk)=1x∑l≥0Wal(G)xl

Why? A tree-like walk in Wal(G) that visits a exactly k times is determined by:

  • A sequence b1,…,bk of neighbours of a.

  • A sequence γ1,…,γk of walks in Tbili(G−a) where 2k+l1+⋯+lk=l.

□