SoS and Planted Clique: Tight Analysis of MPW Moments at all Degrees and an Optimal Lower Bound at Degree Four
Samuel B. Hopkins, Pravesh K. Kothari, Aaron Potechin
Introduction
Let be the Erdős-Rényi random graph where each edge is present in with probability independently of others. It is an easy calculation that the largest clique in is of size with high probability. Recovering such a clique using an efficient algorithm has been a long standing open question in theoretical computer science. As early as 1976, Karp [Kar76] suggested the impossibility of finding cliques of size even for any constant in polynomial time. Karp’s conjecture was remarkably prescient and has stood ground after nearly decades of research.
Lack of algorithmic progress on the question motivated Jerrum [Jer92] and Kucera [Kuc95] to consider a relaxed version known as the planted clique problem. In this setting, we are given a graph obtained by planting a clique of size on a graph sampled according to . Information theoretically, the added clique is identifiable as long a . The goal is to recover the added clique via an efficient algorithm for as small an as possible. This variant is also connected to the question of finding large communities in social networks and the problem of signal finding in molecular biology [PS000]. Despite attracting a lot of attention, the best known polynomial time algorithm can only find planted cliques when their size [AKS98, FK03]. The LS+ semi-definite programming hierarchy leads to the state of the art trade off: planted cliques of size can be recovered in time for any .
Recently, this difficulty of finding cliques of size has led to an increasing confidence in planted clique being a candidate for an average case hard problem and has inspired new research directions in cryptography [ABW10], property testing [AAK+07], machine learning [BR13], algorithmic game theory [DBL09, ABC11] and mathematical finance [DBL10].
In this paper we are interested in understanding the Sum of Squares (SoS, also known as Lasserre) semi-definite programming (SDP) hierachy [Las01, Par00] for the planted clique problem. This is a family of algorithms, paramterized by a number called the degree, where the algorithm takes time to execute. The sum of squares hierarchy can be viewed as a common generalization and extension of both linear programming and spectral techniques, and as such has been remarkably successful in combinatorial optimization. In particular it captures the state of the art algorithms for problems such as Sparsest Cut [ARV04], MAX CUT [GW95], Unique Games/Small Set Expansion [ABS10, BRS11, GS12]. Recently, [BBH+12] showed that a polynomial time algorithm from this hierarchy solves all known integrality gap instances of the Unique Games problem, and similar results have been shown for the hard instances of MAX-CUT [DMN13] and Balanced Separator [OZ13]. Very recently, [LRS15] showed that the sum of squares algorithm is in fact optimal amongst all efficient algorithms based on semidefinite programming for a large class of problems that includes constraint satisfaction and the traveling salesman problem. Moreover, Barak, Kelner and Steurer [BKS14, BKS15] used the SoS hierarchy to give improved algorithms to average case problems such as the dictionary learning problem and the planted sparse vector problem that have at least some similarity to the planted clique problem.Thus several researchers have asked whether the SoS hierarchy can yield improved algorithms for this problem as well.
The first published work along these lines was of Meka, Potechin and Wigderson [MPW15] who showed that for every , the degree SoS cannot find planted cliques of size smaller than .We use to denote equality up to factors polylogarithmic in (the size of the graph) and with an arbitrary dependence on the degree parameter . Deshpande and Montanari [DM15] independently proved a tighter lower bound of for the case of degree . In the main result of this paper, we extend the prior works and show that the first non trivial extension of the spectral algorithm, namely the SoS algorithm of degree , cannot find cliques of size , a bound optimal within factors. Our lower bound for degree is obtained by a careful “correction” to the certificate used by [MPW15] and [DM15] in their lower bounds.
A similar result was obtained in an independent work by Raghavendra and Schramm [RS15]. In our second main result, we give a tight analysis of the certificate considered by [MPW15] and [DM15] and show that it yields a lower bound of .
The certificate of [MPW15, DM15] is sufficient to show an lower bound for the degree hierarchy [FK03] (which is a weaker SDP that also runs in time ). However, a generalization of an argument of Kelner (see Section 10) shows that this is not the case for the SoS hierarchy, and our analysis of this certificate is tight. Hence our work shows that to get stronger lower bounds for higher degree SOS it is necessary and sufficient to utilize more complicated constructions of certificates than those used for weaker hierarchies. Whether this additional complexity results in an asymptotically better tradeoff between the running time and clique size remains a tantalizing open question.
Technical Overview
The SoS semidefinite programming hierarchy yields a convex programming relaxation for the planted clique problem. That is, we derive from the input graph a convex program such that if the graph had a clique of size then is feasible. To show that the program fails to solve the planted clique problem with parameter , we show that with high probability there is a solution (known as a certificate) for the program even when is a random graph from (which in particular will not have a clique of size ).
It is known that such an approach cannot work to obtain a lower bound for the SoS program of degree and higher. That is, this natural certificate does not satisfy the conditions of the SoS program. Hence to obtain our lower bound for degree SoS we need to consider a more complicated certificate, that can be thought of as making a global “correction” to the simple certificate of [MPW15, DM15]. For higher degrees, we have not yet been able to analyze the corresponding complex certificate, but we are able to give a tight analysis of the simple certificate, showing that it certifies an lower bound on degree SoS relaxation. The key technical difficulty in both our work and prior ones is analyzing the spectrum of random matrices that have dependent entries. Deshpande and Montanari [DM15] achieved such an analysis by a tour-de-force of case analysis specifically tailored for the degree case. However, the complexity of this argument makes it unwieldy to extend to either the case of the more complex certificate or the case of analyzing the simple certificate at higher degrees. Thus, key to our analysis is a more principled representation-theoretic approach, inspired by Grigoriev [Gri01a], to analyzing the spectrum of these kind of matrices. We hope this approach would be of use in further results for both the planted clique and other problems.
We now give an informal overview of the SoS program for planted clique, the [MPW15, DM15] certificate, our correction to it, and our analysis. See Section 3 and [BS14] for details.
2 The “Simple Moments”
The simplest form of such a pseudo-distribution is to set
3 There is such a thing as too simple
4 Fixing the simple moments
for every clique . Note that when , the correction factor would typically be of the form .While it might seem that there is a chance for these pseudo-expectations to be negative, if then it is exceedingly unlikely that there will exists an such that , and so we ignore this issue in this overview.
If we now focus on the contribution of the terms where the ’s are all distinct, we see that each such set yields the term
5 Analyzing the corrected moments
The above gives some intuition why the corrected moments might be better than the simple moments for one set of polynomials. But a priori it is not at all clear that those polynomials encapsulated all the issues with the simple moments. Moreover, it is also unclear whether or not the correction itself could introduce additional issues, and create new types of negative eigenvectors. Ruling out these two possibilities is the crux of our analysis.
Because of the dependencies in the random matrix , the deviation from expectation varies according to the eigenspace. Thus, in [DM15], the deviation from the expectation is analyzed by first decomposing along the eigenspaces .
A second technical idea is required to carry out this decomposition. Because of the symmetries present in the spaces , this decomposition is very nearly the same as splitting up the matrix in an ostensibly unrelated way. Each entry of for is the indicator for the presence of a clique on . This indicator is just the AND of all the indicators for the presence of the edges . Taking a Fourier decomposition of this suggests a way to decompose as \mathcal{M}=\sum_{\text{subsetsS4nodes}}\mathcal{M}_{S},This is not quite the whole picture, see Section 7.1.1. where the matrix corresponds to the Fourier character . The matrices can be matched up to the subspaces in such a way that those matrices with larger spectral norm (corresponding to larger deviations from expectation) have subspaces with smaller eigenvalues in their kernels!
The foregoing is missing one subtlety. Some monomial matrices do not match nicely to a single subspace. Instead they form cross terms: for example, having in the left kernel but not in the right kernel (not all these matrices need be symmetric). In fact, it is just such a matrix which keeps the simple moments from remaining PSD beyond .
For the four nodes , consider the monomial and the corresponding matrix . The entries do not depend on , and so there are many repeated rows, which creates a much larger spectral norm for than if it had independent entries. [DM15] prove the (tight) bound . At the same time, it turns out only to have in its left kernel, not its right one. Appealing to the above picture, in order to have , we must have .
We make one further observation about the matrix from the previous section: its rows are the tensor squares of the neighborhood indicator vectors from above. Our fix to the simple moments, described above as adjusting individual pseudo-expectations, amounts roughly to adding to the matrix . This carves out of (our worst subspace from an eigenvalue perspective) a new subspaces space with eigenvalues lower-bounded by . Now instead of matching the bad matrix to and as a cross term, we can match it to and as a cross term. Then we only need . With some care in the details, the above picture can be made precise.
However, a crucial point is that the matrix doesn’t satisfy the clique constraints (in that all entries with not a clique should be ). A chunk of our proof goes into analyzing the discrepancy between the matrix and its zeroed out version. Our analysis here requires the use of new combinatorial tools (Section 6.4) combined with the trace moment method.
As we have alluded to already, a key technical step in our proof is to show that certain Fourier-decomposed matrices of the form discussed above have some of the subspaces in their kernels. In the analysis of simple moments for degree , [DM15] use explicit entries for canonical forms of eigenvectors in to accomplish this. However this approach hits analytical roadblocks for the analysis in case of higher degrees. Canonical forms for the eigenvectors are hard to pin down explicitly from the literature in algebraic combinatorics.
A similar approach was utilized heavily by Grigoriev [Gri01a] to prove a sum of squares lower bound for the knapsack problem. While for degree either explicit eigenvectors or our approach will work, although the latter takes some more elbow grease, ours is absolutely vital for our tight analysis of the MPW moments for the higher degrees. We hope that such an approach will be useful for proving better (approaching ) integrality gap for Sum of Squares relaxations of higher degree for Planted Clique and other related problems.
The analysis of the MPW operator at higher degrees also presents other new challenges that do not show up in the special case of degree analyzed in [DM15]. [DM15] deal with the optimization version of the degree SOS program which could be potentially weaker than the one we analyze here (and thus our lower bound is technically stronger). Working with the “optimization” version simplifies the analysis in [DM15] a little bit as the matrix has entries that only have local dependence on the graph . We explicitly work with the feasibility version of the degree 4 SOS program and thus, must deal with the additional complexity of the entries of having a global dependence. As in [MPW15], we deal with this situation by separating into matrices and such that has only local dependence on the graph . [MPW15] deal with by a simple entrywise bound, however, employing such a bound yields no improvement over the bound proved in [MPW15] for us. It turns out that we have to do a fine grained analysis of the matrix itself by a decomposition for such that each piece is essentially only locally dependent on the graph. Once we have such a decomposition, our ideas from the analysis of can be extended to that setting as well.
Finally, our argument for analyzing the spectral norms of each of the pieces of encountered in the decompositions also needs to be much more general than in case of [DM15] to handle higher degrees. For this, we identify a simple combinatorial structure that controls the norm bounds and allows a general hammer for computing the norms of all the matrices that appear in this analysis. Our proofs here are based on the trace power method and build on the combinatorial techniques in [DM15].
6 Preview of Technical Toolkit
In this section we give a preview of the key lemmas that allow us to carry out the analyses described thus far. We have simplified some issues for the sake of exposition; details may be found in Section 6.
We are concerned with the matrices in the aforementioned Fourier decomposition. Let be a bipartite graph on vertices. Let be an matrix with entries
(We are ignoring what happens if .)
This first lemma bounds the spectral norm of such a matrix in terms of the shape of .
We also want to show that these matrices have nontrivial kernels, so we can bound their negative eigenvalues against the parts of the simple moments with larger positive eigenvalues. The following allows us to carry out this matching of Fourier decomposition matrices to eigenspaces of the expectation matrix.
7 Related Work
There’s a large amount of work on understanding Linear and Semidefinite Programming based hierarchies. A detailed survey on the sum of squares hierarchy and references to works related can be found in [BS14]. The earliest works on proving SoS lower bounds were due to Grigoriev [Gri01a, Gri01b] who showed that degree SoS does not beat the random assignment for 3SAT or 3XOR even on random instances from a natural distribution. Some of these lower bounds were rediscovered by Schoenebeck [Sch08]. Lower bounds for SOS essentially rely on gadget reductions from 3SAT or 3XOR and this approach has been understood in some detail [Tul09, BCV+12]. An exception to this methodology is the recent work of Barak et al. in proving SoS lower bounds for pairwise independent CSPs [BCK15]. Even though the lower bounds for CSPs are for random instances, the average-case nature of the problem does not show up as a main analytic issue. There has recently been a surge of interest in understanding the performance of SoS on average-case problems of interest in machine learning, both in proving upper and lower bounds [HSS15, BM15, MW15, GM15, BKS15].
For the planted clique problem, Feige and Krauthgamer gave an analysis of the performance of the LS+ semidefinite heirarchy tight to within constants [FK00, FK03] giving the state of the art algorithm for finding planted cliques in any fixed polynomial time. Other algorithmic techniques not based on convex relaxations have been studied and shown to fail for planted clique beyond , most prominantly Markov Chain Monte Carlo (MCMC) [Jer92]. Recently, Feldman et. al. [FGR+13] showed a lower bound for (a variant of) the planted clique problem in the restricted class of statistical algorithms that generalize MCMC based methods and many other algorithmic techniques. Frieze and Kannan [FK08] proposed an approach for the planted clique problem through optimizing a degree- polynomial related to the random graph. Such polynomials are NP hard to optimize in the worst case but the belief is that the random nature of the polynomials might be helpful. This approach was generalized to higher degree polynomials by Brubaker and Vempala [BV09].
There has also been recent work on variants of the problem that define Gaussian versions of the planted clique and more generally, the hidden submatrix problems showing, for example, strong inditinguishability results about the spectrum of the associated matrices with and without planting [MRZ14]. Finally, the present work builds heavily on independent papers of Meka, Potechin, and Wigderson [MPW15] and Deshpande and Montanari [DM15], which we have already thoroughly discussed.
Section 3 contains preliminaries. Section 4 contains definitions and the necessary background on the simple moments, a.k.a. the MPW operator. Section 5 contains the formal definition of our corrected degree moments. Section 6 lays out the technical framework for the analyses of the corrected degree moments and for the tightened bounds on the MPW operator at higher degrees. Here we define the Fourier decompositions alluded to above and carry out representation-theoretic arguments about their kernels. Section 7 and Section 8 use the tools we have built thus far to prove the main theorems. In Section 9 we prove a technical concentration result for small subgraphs of required for the analysis. In Section 10 we sketch Kelner’s argument showing that our analysis of the MPW moments is nearly optimal.
Preliminaries
We will use the following general notation in the paper.
will denote a draw from unless otherwise stated.
For a square symmetric matrices , we write to mean is positive semidefinite.
For any matrix , denotes its largest singular value, or, equivalently,
For matrices of same dimensions, denotes their entrywise or Hadamard product, i.e., for every .
For a graph and any set of two vertices of , , denotes the indicator of the edge being present in . That is, if is an edge in and otherwise.
For a set of vertices of , , the set of all pairs from .
For a pair of subsets of vertices of , the set of cut edges between and .
For a subspace , denotes the projector to .
Following [BBH+12] and many subsequent papers, we work with SoS using the language of pseudo-expectations.
Let . If a pseudo-expectation satisfying the constraints exists, it can be found in time . If none exists, a certificate of infeasibility of these equations is found instead.
The following observations (actually both the same observation in different forms) will come in handy in our analysis.
by our assumptions on and Cauchy-Schwarz. For each , we know , which implies that the whole expression is nonnegative. ∎
The MPW Operator
For any set , let be the monomial .
Further, we set be the number of -cliques in .
Our analysis improves upon the analysis in [MPW15] and generalizes the improved analysis for the special case of done in [DM15]. By a generalization of the counter example due to Kelner, our analysis can actually be shown to be tight. We defer the details of the counter example to the full version. In the remaining part of this section, we begin the task of proving Theorem 4.4 by introducing certain simplifications and computing the eigen values of the expected value of the matrix under .
Observe that for any , is chosen so as to depend only on the edges with one end point in and the other in . Intuitively, this corresponds to thinking of and as being cliques in by addition of some edges. Moreoever, is chosen so that whenever are actually cliques in . Thus, as noted above, we have the following fact (which is Lemma 5.1 in [MPW15]).
is PSD if is PSD.
Thus, the following lemma completes the proof of Theorem 4.4.
2 The Expectation Matrix
We first describe the entries of the matrix .
Next, we need a basic fact about the (shared) eigenspaces of all the set symmetric matrices, in particular, their number and dimensions which follows from the following well known result from classical theory of Johnson schemes.
For , .
Using a nice basis for the matrices in , one can obtain the following estimates of the eigenvalues of on for each :
Let and . Let be the eigenvalue of on as defined in Fact 4.9. Then,
The Corrected Operator for Degree Four
In this section, we present the pseudodistribution that we will use to show an almost optimal lower bound on the degree SOS algorithm. Our pseudodistribution is obtained by “correcting" the one described in the previous section. The correction itself is inspired by an explicit polynomial described by Kelner who showed that the pseudodistribution from the previous section for degree does not satisfy positivity for .
Let be a real parameter to be chosen later. Let be the linear operator on the linear space of homogeneous, multilinear polynomials of degree such that:
For and , there exists , , such that .
We first use the correction operator to define the corrected moments on all multilinear monomials of degree . For every , , we set:
Then we know there exists satisfying (using Fact 5.2). Thus, the degree moments we defined “think” . We use this relationship to extend the definition to all the monomials. For every , ,
Furthermore, if , then with probability , .
There is and so that for , with probability at least over the draw of , .
We prove the first item; the others are similar.
The proof is by several applications of McDiarmid’s inequality. By a standard Chernoff bound there is a universal constant so that for every and every , with probabilty ,
for some other universal constant . So by McDiarmid’s inequality there is so that with probability ,
By a similar argument there is so that for every , with probability ,
for some other constant . The result follows by a final application of McDiarmid’s inequality (we lose a factor of at this step as opposed to the at previous steps because there are edges to be revealed). ∎
We can now complete the proof of Lemma 5.5 using Lemma 5.7 and Fact 5.2.
In Section 6, we develop some general tools for analyzing the matrices that we encounter before going on to prove Theorem 4.4 and Lemma 5.6.
Tools
We provide background in the required tools from basic representation theory below.
For a finite dimensional complex vector space , let be the set of all linear maps from into . For any finite group and , the pair is said to be a representation of if satisfies, for any ,
where the “” on the LHS corresponds to the group operation and on the RHS, the composition of linear maps on . When the map is clear from the context (as some natural action of the group on ), we abuse notation and just say that is a representation of .
Let be a representation of a group . A subspace is said to be a subrepresentation if for every , for every . That is, is a stable or invariant subspace for all the linear maps , one for each . Observe that in this case, is another representation of . A representation of is said to be irreducible if for any subspace invariant under all the linear maps for , or .
Suppose and are representations of a group . Suppose is a linear map such that for any and ,
Then, for any irreducible representation under , is an irreducible representation in under .
2 Eigenspaces of the Set Symmetric Matrices
The set of all set symmetric matrices is known as the Johnson scheme in algebraic combinatorics. All such matrices commute and thus share eigenspaces. While the matrices in the Johnson scheme are well studied, the description of the eigenspaces in the literature is hard to use for the purpose of our proofs. We thus take a more direct approach and use basic representation theory in what follows to identify a simple symmetry condition on the eigenspaces of set symmetric matrices which will be useful to understand the spectral properties of the matrices we study.
Then, .
3 Kernels of Patterned Matrices
In this section we design some general tools to understand the spectral structure of matrices that have restricted variations around the set symmetric structure discussed in the previous section. The main tool we will use to establish these results is Lemma 6.4 shown in the previous section. Before moving on to this task, we describe a high level overview of what we intend to do. The following paragraph can be skipped to dive directly into the technical details without the loss of continuity.
The study of the eigenspaces of set symmetric matrices lets us completely understand the spectral structure of the expectation matrix . In the next section when we analyze the spectrum of , we will encounter matrices that depend on the underlying graph and thus are not set symmetric. However, if the dependence on the underlying graph is in some sense limited, we hope that some of the nice algebraic properties that set symmetry grants us should perhaps continue to hold. In our case, we will be able to decompose into various pieces and for each of these pieces, the entry at has dependence on the graph based only on the status (edge vs no edge) of a small number of pairs . The goal of this section is to develop tools to understand certain (coarse) spectral properties of such matrices.
and the right automorphism group of as
We are now ready to define patterned matrices:
When , we write for the corresponding patterned matrix.
The following result describes the kernels of certain symmetrized sums of and is the main claim of this section.
We first show that above is well defined in that the definition does not depend on the specific subset used so long as . We adopt the notation (that ignores the “direction” of action) only for the calculations that follow.
We can equivalently write the claim above as:
This completes the proof using Lemma 6.4. ∎
4 Concentration for Locally Random Matrices over G(n,12)𝐺𝑛12G(n,\frac{1}{2})
For , , and a bipartite graph , let be a patterned matrix with . That is,
The next main result of this section considers a different class of matrices that appear in the analysis in Section 8.
Let be a bipartite graph on vertices and suppose is nonempty. Let be a matrix with the entry at , for and are given by
The proofs of both these results are based on the standard idea of analyzing the trace of higher powers of a matrix to prove bounds on its spectral norm. The proof of Lemma 6.11 is similar to the proofs via the trace power method for bounding the norms of matrices as presented in [DM15]. The general format we present here will come in handy for multiple applications to various matrices in Section 7. Lemma 6.12 deals with somewhat more complicated matrices that appear in the analysis of the corrected operator for degree lower bound. Nevertheless, as is common in such proofs, the analysis is based on a combinatorial analysis of the terms that make non zero contribution to the trace powers combined with the simplifying effect of random partitioning based arguments. We describe the details of the proof in the following section.
It is enough to show that occurs as a principal submatrix of . For this, take the submatrix of rows and columns of indexed by tuples in sorted order, i.e., with . ∎
We will use the following lemma to break dependencies in certain random matrices by decomposing them into matrices whose entries, while still dependent, have additional structure.
Suppose when . Let be a sequence of partition of into bins. Each partition induces a matrix based on as follows:
Then, there is a family of partitions such that with .
where is the set of indices so that respect the partition and do not respect any partition for any .
Then, there is a family of partitions of so that with .
We present the proof of 1; the proof of 2 is almost identical. For to be chosen later, we pick partitions uniformly at random and independently so that each is partition of into sets of size each.
Call good at step if for every and if . It is enough to show that after steps the probability that every of size is good at some step .
Fix some with . It is good at step with probability at least . Since the steps are independent, after steps
which is at most for some .
Taking a union bound over all tuples with completes the proof. ∎
The usefulness of the above definition is captured by the following claim:
4.2 Graph-Theoretic Definitions and Lemmas
In this section, we set up some notation and definitions helpful in our proofs of the main results of this section. The next few definitions and notation are generalization of the ones used in [DM15] to general degrees and are useful in the proof of Lemma 6.11.
Let be a graph. A labeled -ribbon is a tuple where is a -ribbon and is a map labeling each vertex of with a vertex in . We require that for an edge in , .
Let be a labeled -ribbon where has vertices. We say is disjoint if for every ,
Let be a labeled -ribbon where has vertices. We say that is contributing if no element of the multiset occurs with odd multiplicity.
The following combinatorial lemma will serve as a tool in the proofs of the main results for this section.
The assumption on implies that contains the cycles
The next few definitions and notation are needed in the proof of Lemma 6.12.
Where is a graph, a labeled fancy -ribbon is a tuple where is a fancy -ribbon and labels each vertex of with a vertex in . We require for any edge that .
When is empty the proof is similar: there are two paths, and . ∎
4.3 Proofs of Lemma 6.11 and Lemma 6.12
By Lemma 6.13 it is enough to prove the analogous claims for the matrix with entries given by
By multiplying by suitable permutation matrices to give , we may assume in the -matching case above that the matching is and in the nonempty graph case that the edge contained is (where we think of the vertex set of as ). Note that .
We apply Lemma 6.14 to obtain a family of matrices for some satisfying . On any entry on which is nonzero it is equal to at that entry, and furthermore for each there is a partition of so that if then , and for all .
Taking a union bound over the matrices , we get that
and so by the triangle inequality applied to , we get
We first handle the case when is non empty. By Lemma 6.13 it is enough to prove the analogous statement for the matrix, also by abuse of notation denoted , which is the sum of matrices (abusing notation again) with entries given by
By multiplication with an appropriate permutation matrix (which cannot change the spectral norm), we may assume that contains the edge . We begin with 2 from Lemma 6.14, whose hypotheses are satisfied by our convention . This gives a family with so that and a corresponding family of partitions . Lemma 6.14 guarantees that for some when and and is zero otherwise.
By a union bound and triangle inequality, we then get
The proof in the case of empty is similar, using Lemma 6.23 in the empty case.
Analyzing Deviations for the Degree-d𝑑d MPW Operator
In this section, we use the tools developed in Section 6 to analyze the spectrum of the deviation matrix and prove Lemma 4.7.
As noted in Section 7, we decompose . For any , depends on a) and b) whether . If depended only on b) above, then it could be decomposed into a sum of patterned matrices defined in Section 6; analyzing these is tractable. Our first step is thus to get rid of the dependence on —the only part depending on the entire graph. We will obtain a matrix that depends only on whether or not (and thus is “locally random" in the sense of [MPW15]).
and for each . We set for every to be
The first is regarding the high level approach. The approach of [DM15] used explicit expressions for a canonical set of eigenvectors in to obtain similar conclusions as us for the case of . This approach gets unwieldy very quickly because the explicit entries of eigenvectors for for are hard to work with [BI84]. We tackle this issue by developing an argument that doesn’t need explicit entries of the eigenvectors. Instead, we use basic representation theory (Section 6.3) to identify a set of symmetries satisfies by vectors in for each and use it obtain the conclusions we require.
Second, [DM15] deal with the optimization version of the degree SOS program which, as noted in the introduction, could be potentially weaker than the one we analyze here (and thus our lower bound is technically stronger). This simplifies the analysis in [DM15] a little bit as the matrix defined above is identically zero for the operator analyzed. We explicitly work with the feasibility version of the degree 4 SOS program and thus, must deal with the additional complexity of handling . It turns out that we have to do a fine grained analysis of the matrix itself. The decomposition we use for is somewhat different from the case of even though, the analysis of each piece of the decomposition proceeds similar to the case of .
Third, for the special case of , essentially the only matrix one has to analyze is the , the matrix obtained by zeroing out all entries in such that : a uniform bound on spectral norm of the remaining component suffices. However, for higher , one has to deal with the “non-disjoint" entries with some care and an argument analogous to the one in [DM15] fails to show PSDness of beyond giving no asymptotic improvement over [MPW15].
Finally, our argument for analyzing the spectral norms of each of the pieces also needs to be much more general than in case of [DM15] to handle higher degrees. For this, we identify a simple combinatorial structure (size of maximum matchings in appropriate bipartite graphs on vertices) that controls the bounds and could also be used to obtain slick proofs of the conclusions required in [DM15] in the context of analyzing for . Our combinatorial argument itself is a generalization of the one given by [DM15] for this case.
With probability at least over the draw of , each satisfies
The next lemma does a (even more) fine grained analysis of the spectrum of :
With probability at least over the draw of , for each :
We can now use Lemma 7.2 and Lemma 7.3 to complete the proof of Lemma 4.7.
For each , we compute and use Lemma 3.4. We write . First, whenever as are projectors to eigenspaces of . Let be the eigenvalues of on eigenspaces . From Lemma 4.10, we have:
In what follows, all our statements hold with probability at least :
Using Lemma 7.2 and Lemma 7.3, for every ,
Then, it is easy to check that for any ,
Next, we bound the cross terms for . Again, using Lemma 7.2 and Lemma 7.3, we have for :
For , it is again easy to check that, for , the above expression is at most .
In the case when one of is at most , we have the bound
In this case, it is easy to check that so long as ,
By an application of Lemma 3.4, the proof is complete. ∎
We first describe the high level idea of the proof.
We start by decomposing where if and otherwise. Notice that each then is obtained by a scaling an appropriate matrix.
Most illuminating is the disjoint case , which is nonzero only at entries with . For any disjoint , depends only whether , which, one could write as an appropriately scaled AND function of the indicators of edges . We can expand this AND function in the monomial (parities of subsets of variables) basis. Each such monomial corresponds to the bipartite graph that contains the pairs that constitute the monomial. This gives a decomposition of into (since the constant term is , being zero mean) components, for each non empty, labeled bipartite graph on .
We can bound the spectral norm of each of the pieces by direct application of tools derived in Section 6.4. The main work in this section goes into showing that depending upon the structure of , an appropriate selection of subspaces lie in left or right kernels of . Thus, for a fixed term , some do not contribute. We identify the maximum spectral norm among contributing terms to obtain the final bound.
To accomplish this goal, we rely heavily on the tools built in Section 6.3 which give us a handle on the symmetries of the eigenspaces . This requires some work based on representation theory of finite groups and is presented in Lemma 6.4 and Lemma 6.9.
The case of for needs even finer decomposition. We decompose each () into matrices that identify the “pattern" of the intersecting vertices. In [MPW15] a similar idea is used to reduce the task of bounding the spectral norm of to a calculation similar to one in the case of . However, unlike [MPW15], we also require properties of the kernels of the components of the decomposition. After restricting to a fixed intersection pattern of vertices, we thus resort to using a generalization of the kernel analysis used for the case. We now proceed with the proof plan as described beginning with the decomposition of each .
1.1 Decomposing L𝐿L
We start by decomposing further as where for any ,
Recall that is the set of all bipartite labeled graphs with left and right vertex sets labeled by and . Recall also from Section 6.3 that for any with , the graph is a copy of on vertex sets where the correspondence between and and and is determined by the sorting map . Finally, recall that for a graph on , we let be the indicator for the presence of edge in , and by convention when for any . For any , define an matrix
We think of as a rescaling and centering of a matrix whose entries are the AND of the indicators for the edges in . Decomposing these ANDs into monomials over those indicators, we see that each monomial corresponds exactly to one bipartite graph , and the centering of corresponds to removing the constant monomial, which corresponds to the empty bipartite graph. Every other monomial recieves equal weight in this expansion, and so from these observations it becomes routine to verify that
Similarly, we further decompose for . Here things are a bit more involved. Let us motivate our decomposition by understanding the structure of the matrix for a little bit. Consider an entry such that . Then, . Thus, the edge structure in the bipartite subgraph on vertex sets and decides the value of for any graph and we can hope to a get a patterned matrix. We now follow this intuition.
We now define a matrix which is nonzero only on entries which intersect in at least places:
Finally, it is again easy to verify that:
1.2 Spectral Analysis of L𝐿L
In order to prove Lemma 7.2, we will first use the decomposition described in the previous section to write as a sum of appropriate patterned matrices. We will then partition the sum into groups, each group corresponding to an equivalence class of (left or right) similar bipartite graphs . We will infer some properties about the kernel and finally use the spectral norm bounds from Section 6.4 to complete the proof.
With probability at least over the draw of ,
Before proving the three claims above, we show how they imply Lemma 7.2:
For the second part, fix an . We first show that some terms in the decomposition in (7.3) do not contribute to .
In the remaining part of this section, we complete the proofs of Claims 7.6 and 7.5.
1.3 Proof of Claims
In this section, we obtain quick proofs of the three claims above using the tools developed in Section 6.
The proof is by appealing to Fact 3.3. Observe that
The next is a direct application of Lemma 6.9.
Finally, we prove Claim 7.5 using Lemma 6.11.
where is some fixed subset of size such that .
2 Proof of Lemma 7.3
We now move on to analyzing the spectrum of the matrix .
The high level plan of the proof is similar to that of Lemma 7.2. We define for each as follows:
We further split where:
First, we observe that . This is because when and are disjoint is exactly and doesn’t depend on the graph. Thus, As before we would like to spot patterned matrices in each to show that appropriate eigenspaces lie in the kernel of . In case of , however, there’s a difference how this needs to be done. This is because each entry of potentially depends on the edges from every vertex in the graph to and . This is unlike the case of where the entry depends only on the edges between and (in fact that’s the reason we separated from in the analysis). Nevertheless, we give a decomposition below that will help us make claims similar to the ones in the case of analyzing in this case too.
Let us first explain the main idea in the decomposition. The entry depends on two events: a) whether and b) the number of subsets of size that form -cliques with . The main observation that motivates our decomposition is the following: in the event that , the deviation in is completely captured (up to low order terms) by just the number of vertices that has an edge to all of in . This allows us to write entries of as a sum of contribution to the deviation due to each vertex separately. For the case of , this argument is in fact exact and there are no low order terms. When , the contributions due to individual vertices contribute the bulk of the deviation and only low order terms remain.
From here on, we are in a situation similar to the one encountered in analyzing in the previous subsection. We show that the components in the decomposition with large spectral norm do not contribute to quadratic forms over eigenspaces with small eigenvalues of the expectation matrix using the idea of patterned matrices from Section 6.3. We show that the remaining components have small spectral norm using the combinatorial techniques combined with the trace moment method developed in Section 6.4.
Using the matrices above, we can show the following approximate factorization for the entries of :
for
Let be the set of vertices in not in so that for all . By definition, if is a clique,
and hence the lemma follows from Theorem 9.4. ∎
We will need another definition before proceeding: For each ,
As in the case of analyzing , we define the filled in versions of as follows:
We start by giving norm bounds on all the matrices involved in the decomposition in Lemma 7.7.
Then, in the sense of Definition 6.15. Thus, by Fact 6.16 Thus, we focus on bounding in the following. We will use the trace moment method for this purpose and the argument is similar to the ones made in Section 6.4 to develop the general purpose spectral concentration results. For this reason, we will be a bit more terse than in the case of the other applications of the trace moment method before. We set up the notation for the general as above and specialize the combinatorial reasoning for each of the specific matrices involved later.
We now investigate when does a term in the expansion above contribute a non-zero value to the LHS.
We can now complete the proof of Lemma 7.3 using Lemma 7.8and Lemma 7.9.
For each , let be the expression given by Lemma 7.7 to approximate the entries of . Then, we have for any :
Analyzing Deviations for the Corrected Degree-4 Operator
The goal of this section is to prove Lemma 5.6, that is, to show that . The proof is organized into main claims that we next present.
(Recall that is the number of -cliques in .) We then set
where is the filled in MPW matrix (Definition 4.5).
Our first claim shows that it is enough to prove PSDness of :
For any there is so that for any with probability if then .
Let be the projector to for every .
For every there is so that for , with probability ,
Next, we analyze the spectrum of . Here at last we see the main technical improvement in these corrected moments—the cross term between and has become a cross term between and and has much-reduced norm.
Next, we bound the spectral norm of . This is a direct corollary of the more general bound in Lemma 7.3, but to have a self-contained proof of the degree-4 case we also give a proof later in this section.
Let be the filled-in matrix for the degree-4 MPW moments with clique size (Definition 4.5). Let be as in Definition 7.1. With probability , .
Let . With probability , .
The proofs of these lemmas follow, but first we complete the proof of Lemma 5.6 and hence of Theorem 5.4.
By Lemma 8.1, it will be enough to exhibit and so that with probability when . Then our final bound will be given by the minimum of and of Lemma 8.1. (Recall that is a parameter inside .) In the following, all that we claim happens with probability at least by a union bound.
First of all, by Lemma 8.5, for every there is so that
We assume is chosen in this way.
For any , by Lemma 8.5 if we choose then . Adding to the previous equation,
By the same reasoning using Lemma 8.4 to add to the previous equation, we have for the same choice of :
So it just remains to add to the left-hand side.
Let be as in Lemma 8.3, which implies that
Choosing we get that this is . So using Lemma 3.5 to add (8.3) and (8.4), we obtain
Together with Lemma 3.5 and (8.5), this implies the lemma, for , as above, and . ∎
We start with the easy parts. Note that , so the bound
2 Lower-Degree Cleanup (Lemma 8.1)
Every with all distinct satisfies
With probability when by Lemma 5.7, we get that . By the same,
which is for and . ∎
We first observe that an eigenvalue lower bound on implies the same on the (principal) submatrix indexed only by cliques (in this case, edges) in . This submatrix is equal to , where
We break Err into two parts so that :
Note that consists of the off-diagonal nonzero entries of Err, while contains the diagonal entries of Err. We start by showing a bound on .
Next we bound . Since it is diagonal, it is enough to give an entrywise bound.
3 Eigenvalue Lower Bound for the Correction
The following is the main claim for this section.
We will need the following graph theoretic machinery for the moment method bound.
The labeled diamond ribbon is contributing if no element of the multiset occurs with odd multiplicity. It is disjoint if the sets
By our disjointness assumption, every element of the multiset must occur with multiplicity at least two and similarly for and . ∎
Note that the matrix has row and column spaces both . Note also that it factors as , where is the matrix whose columns are the vectors . Thus, it will be enough to show that
where here is the identity matrix.
For this, consider the matrix indexed by vertices . It has entries
So, zeroing this matrix on the diagonal, it is enough to prove that
Let . Let for be given by
Then . Note that if or . Thus the obvious generalization of Lemma 6.14 to the two-parameter family applies. This gives us a family of matrices for some and a corresponding family of partitions of .
These will be such that , and
4 Proofs of Remaining Lemmas
Note that for disjoint we have . We bound the maximal sum across any row of . With probability every off-diagonal entry off is at most in absolute value [MPW15, Theorem 10.1]. For each , we then get . At the same time, The diagonal entries are each at most with similar probability, again by [MPW15, Theorem 10.1]. ∎
Now it is enough to bound . This is the deviation introduced by zeroing non-clique entries. Note that each entry of can be decomposed as a function of the underlying graph edge variables:
where we recall that for an edge , the variable is the indicator for that edge. Each of the entries in the sum over nonempty above corresponds to a matrix of the form bounded in Lemma 6.12, where we conclude that each has spectral norm at most with probability . We conclude (also using , see [MPW15, Theorem 10.3]) that as desired. ∎
In this section, we prove the following large-deviation bounds on the number of -cliques a random graph contains and on . Similar results (which are likely sufficient for our needs) appear in the literature; see [Ruc88, Vu01, JLR11] for instance. We provide these proofs for completeness. A coarser concentration result for appears in [MPW15].
For a graph , define to be the number of -cliques in .
Unless otherwise specified, in this section .
This first theorem gives the large deviation bound for the number of cliques of size in .
For all , for all , if then
We also want a large deviations inequality for the number of cliques of size that a clique of size participates in in . Moreover, to carry out the eigenspace splitting arguments needed for Lemma 7.3, we want to know the dependence of this deviation on the number of vertices adjacent to every vertex in the -clique. The following theorem serves both these purposes.
Given any , let be the set of all vertices not in which are adjacent to all vertices in .
There is a universal constant so that for any of size at most , if is a clique in then for any , if then
The key lemma in proving Theorem 9.2 is the following, which bounds how often subsets of of size share (potential) edges.
If and then there are at most multi-sets of subsets of size such that that for all there exists an such that .
Using Lemma 9.5, we can bound the deviation of the number of -cliques in from its expected value; we carry this out now.
Define where if is a clique in and otherwise.
If and then .
The result is trivial for and so we may assume that . Using Corollary 9.8 and Markov’s inequality,
Thus, we just need to give an upper bound on . For all positive even , so this expression is upper bounded by . We now try to minimize over all positive even . Taking the derivative of this expression with respect to yields . Setting this to yields . However, we require to be even so we take to be the smallest positive even integer which is greater than . Now and . Putting everything together, for this , . Plugging this in gives
All that is left is to check that for this to make sure that our application of Corollary 9.8 was valid. Since , this holds, as needed. ∎
Now that we have proven Theorem 9.2, we will derive Theorem 9.4 from Theorem 9.2. The idea is that conditioned on being a clique, by Theorem 9.2, is primarily determined by , which can be easily shown to be tightly concentrated around its expected value. We start with the following lemma
If then for any of size less than , if we first determine all of the edges incident to elements of (which determines ) then if is a clique, when we look at the remainder of the graph, for any ,
so long as the following conditions hold:
To prove Lemma 9.9 we require the following results; proofs of the more elementary ones are deferred to Section 9.1.
For all nonnegative integers and where ,
Note that
This implies that and dividing everything by gives the claimed result. ∎
If then
Applying Proposition 9.10 on and gives
Multiplying this equation by gives the claimed result. ∎
For all nonnegative integers and where , .
For any nonnegative integer and any such that ,
Eventually in the course of proving Lemma 9.9 we will break
into pieces. The following lemmas offer the necessary bounds on each piece.
If then
Applying Proposition 9.13 with and , since ,
Multiplying this equation by and using Proposition 9.12 with gives
If then
Applying Proposition 9.10 on and gives
Multiplying this equation by gives the claimed result. ∎
For all , if then
This lemma follows immediately from applying Theorem 9.2 on the random graph restricted to the vertices . ∎
into four parts and analyze each one separately.
Combining Lemma 9.11, Lemma 9.14, Lemma 9.15, and Lemma 9.16, we have that under the given conditions,
The result now reduces to showing the following equation
which follows from the facts that , , and . ∎
To use Lemma 9.9 to prove Theorem 9.4, we need probabilistic bounds on .
There is a universal constant so that for all ,
We have all we need now to prove Theorem 9.4.
The result is trivial if and follows immediately from Theorem 9.2 if so we may assume that .
In what follows, and denotes universal constants which may vary from line to line. Now recall that by Lemma 9.9, for any ,
so long as the following conditions hold:
Taking , plugging Lemma 9.17 into these equations and using the union bound, we have that
so long as the corresponding conditions hold. Assuming these conditions hold for now, since , , , and ,
as needed. For the first part of Theorem Theorem 9.4, note that this implies that
Taking and ,
as needed. All that is left is to check the conditions for Lemma 9.9, which are as follows.
These conditions are true if . To see this, note that since and
Dividing the first statement by gives the first condition. Using the second and third statements,
Dividing this by gives the second condition, as needed. ∎
For each multi-set of -cliques, define the constraint graph as follows.
Let’s first bound the number of such that is connected.
For any , there are at most multi-sets of -cliques such that is connected.
Since is connected, we can order so that for all there is an such that . Assuming this is the case, we have at most choices for . For each , there are two vertices in which are contained in some where . There are at most choices for this and then there are choices for which two vertices of are contained in . There are at most choices for the other vertices of so for each there are at most choices for . Putting everything together, there are at most
multi-sets of -cliques such that is connected. ∎
Now consider the number of multi-sets of -cliques such that has connected components with sizes . Using Lemma 9.19, there are at most
We now total this up over all possible . For the special case that all connected components of have size , there are at most such . We will show that this term contributes more than all of the other terms combined, which implies that the total number of is at most , as needed.
For a given , . Also, each of the components of must have at least two vertices so to determine the sizes it is sufficient to decide how to distribute the extra vertices among the connected components of . There are at most ways to do this. Thus, the total contribution for terms of a given is at most
Since , the term which comes from contributes more than all the other terms combined, as needed. ∎
Rearranging this equation gives which just says that if we want to pick two elements in we can either pick two elements in , one element from and one element from , or two elements from . ∎
Optimality of MPW Analysis
In this section we sketch an argument due to Kelner showing that the MPW moments are not PSD when .
With high probability, the MW moments are not PSD when . In particular, if then for all , for some appropriately chosen , with high probability,
We will be using the following proposition heavily.
We have the equation that so
For all and all such that ,
Analyzing the third term is trickier. For the third term, taking for each term,
For the second part, for each there are at most such that and . Thus,
From our concentration bounds, with high probability, for all the corresponding term on the right is which is . For all , this is much smaller than . Thus we may ignore the second part.
For the first part, the values are completely independent of the values . When we take the expected value of
over the edges incident to , only the term remains and we obtain . When we take the variance of
over the edges incident to , only the square terms remain so we obtain
From our concentration bounds on the , with high probability this is
Acknowledgements
We thank Jonathan Kelner, Ankur Moitra, Aviad Rubinstein, and David Steurer for many helpful discussions. We especially thank Boaz Barak for introducing the problem to us, for numerous conversations which led to many of the proofs in this paper, and for invaluable help and advice in writing it.
S.B.H. acknowledges the support of an NSF Graduate Research Fellowship under award no. 1144153.