The Theta Number of Simplicial Complexes
Christine Bachoc, Anna Gundert, Alberto Passuello

TL;DR
This paper generalizes the Lovász theta number from graphs to simplicial complexes using cohomology and Laplacians, providing new bounds on independence numbers through semidefinite programming.
Contribution
It introduces a novel higher-dimensional theta number for simplicial complexes based on cohomology and Laplacians, extending graph bounds to complex structures.
Findings
The higher-dimensional theta number relates to Hoffman's ratio bound and chromatic number.
A hierarchy of semidefinite bounds converges to the independence number.
Analysis of the theta number on dense random complexes shows its effectiveness.
Abstract
We introduce a generalization of the celebrated Lov\'asz theta number of a graph to simplicial complexes of arbitrary dimension. Our generalization takes advantage of real simplicial cohomology theory, in particular combinatorial Laplacians, and provides a semidefinite programming upper bound of the independence number of a simplicial complex. We consider properties of the graph theta number such as the relationship to Hoffman's ratio bound and to the chromatic number and study how they extend to higher dimensions. Like in the case of graphs, the higher dimensional theta number can be extended to a hierarchy of semidefinite programming upper bounds reaching the independence number. We analyze the value of the theta number and of the hierarchy for dense random simplicial complexes.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
The theta number of simplicial complexes
Christine Bachoc
Institut de Mathématiques de Bordeaux, UMR 5251, Université de Bordeaux, 351 Cours de la Libération, 33400 Talence, France.
,
Anna Gundert
Mathematisches Institut, Universität zu Köln, Weyertal 86-90, 50931 Köln, Germany.
and
Alberto Passuello
Institut de Mathématiques de Bordeaux, UMR 5251, Université de Bordeaux, 351 Cours de la Libération, 33400 Talence, France.
Abstract.
We introduce a generalization of the celebrated Lovász theta number of a graph to simplicial complexes of arbitrary dimension. Our generalization takes advantage of real simplicial cohomology theory, in particular combinatorial Laplacians, and provides a semidefinite programming upper bound of the independence number of a simplicial complex. We consider properties of the graph theta number such as the relationship to Hoffman’s ratio bound and to the chromatic number and study how they extend to higher dimensions. Like in the case of graphs, the higher dimensional theta number can be extended to a hierarchy of semidefinite programming upper bounds reaching the independence number. We analyse the value of the theta number and of the hierarchy for dense random simplicial complexes.
1. Introduction
The theta number of a graph was introduced by L. Lovász in his seminal paper [32], in order to provide spectral bounds of the independence number and of the chromatic number of . In modern terms, is the optimal value of a semidefinite program, and as such is computationally easy; in contrast, the independence number and the chromatic number are difficult to compute. These graph invariants satisfy the following inequalities, where denotes the complement of :
[TABLE]
The inequality was one of the main ingredients in Lovász’ proof of the Shannon conjecture on the capacity of the pentagon [32]. More generally, this inequality plays a central role in extremal combinatorics, sometimes in a disguised form: to cite a few, the Delsarte linear programming method in coding theory [8] and recent generalizations of Erdös-Ko-Rado theorems [7, 12, 13] can be interpreted as instances of this inequality. Analogs of the theta number in geometric settings have lead to many advances in packing problems (see [36] and references therein), in particular the very recent solutions to the sphere packing problems in dimensions and [5, 40].
Our aim in this paper is to generalize this graph parameter to higher dimensions, in the framework of simplicial complexes. Let us recall that an (abstract) simplicial complex on a finite set is a family of subsets of called faces that is closed under taking subsets. We refer to Section 1 for basic definitions and results about simplicial complexes. Graphs fit in this framework, being simplicial complexes of dimension . In recent years, considerable work has been devoted to generalizing the classical theory of graphs to this higher-dimensional setting. Much of the efforts have focused on the notion of expansion (see, e.g., [9, 15, 20, 27, 33, 38]), but other natural concepts such as random walks [37], trees [11, 26], planarity [35], girth [10, 34], independence and chromatic numbers [14, 19] have been extended to higher dimensions. Some of these notions were introduced and studied previously in the context of hypergraphs. Pure -dimensional simplicial complexes are essentially -uniform hypergraphs, but the topological point of view brings the machinery of algebraic topology such as homology theory to the subject.
The familiar graph-theoretic notions of independence number and of chromatic number extend in a natural way to this setting: For a -dimensional simplicial complex , an independent set is a set of vertices that does not contain any maximal face of , and the independence number is the maximal cardinality of an independent set. The chromatic number111In the study of hypergraphs, the chromatic number is also known as the weak chromatic number while , the chromatic number of the -skeleton, is known as the strong chromatic number. is the least number of colors needed to color the vertices so that no maximal face of is monochromatic, in other words, it is the smallest number of parts of a partition of the vertices into independent sets.
In order to define the theta number of a pure -dimensional simplicial complex , we will follow an approach that leads in a natural way to the inequality . The main idea is to associate to an independent set a certain matrix, and then to design a semidefinite program that captures as many properties of this matrix as possible. The matrix that we associate to an independent set is (up to a multiplicative factor) a submatrix of the down-Laplacian of the complete complex. In the case of dimension , the down-Laplacian is simply the all-one matrix, and we end up with one of the many formulations of the Lovász theta number.
Our first task will be to compare to the eigenvalue upper bound of proved by Golubev in [19]. This upper bound involves for , the largest eigenvalues of the -th up-Laplacians of and the minimal degrees of the -faces of :
[TABLE]
When every possible -face is contained in at least one -face, i.e., when has a complete -skeleton, this inequality simplifies to
[TABLE]
and can thus be seen as a natural generalization of the celebrated ratio bound for graphs attributed to Hoffman (see, e.g., [4, Theorem 3.5.2]). In that case, we will show that
[TABLE]
therefore provides an upper bound of that is at least as good as (3). In the case of a non-complete -skeleton, Golubev’s bound and turn out to be incomparable, as we will see in examples below.
The theta number of a graph has many very nice properties; some of them, although unfortunately not all of them, can be generalized to higher dimensions. Most of this paper is devoted to determining which of the properties of the graph theta number extend to our notion of the theta number of simplicial complexes.
The relationship to the chromatic number generalizes only partially. Indeed, the inequality immediately leads to the inequality . However, in the case of graphs, the stronger inequality holds. We will see that its natural analog in the setting of -complexes would be that and that this inequality does not hold in general. Instead, we will introduce an ad hoc notion of chromatic number for simplicial complexes, denoted , and show that the inequality holds. While is defined using vertex colorings, the definition of is based on colorings of -faces respecting orientations. Moreover, it is tightly related to a notion of homomorphisms between pure -dimensional simplicial complexes that we introduce and that may be of interest by itself.
A very interesting benefit of the theta number of a graph is that it is possible to expand it into hierarchies of semidefinite upper bounds of the independence number; Lassere’s hierarchy based on polynomial optimization principles is one of the most popular (see [29, 30]). We will see that a similar situation holds in higher dimensions: to a pure -dimensional complex we will associate a sequence for such that
[TABLE]
In order to define , we will proceed in two steps: in a first step, we define a natural sequence for ; in a second step, we modify the definition of slightly in such a way that the sequence of its values decreases.
Our last results concern the theta number of random simplicial complexes from the model proposed by Linial and Meshulam in [31]. This model is a higher-dimensional analog of the Erdős-Rényi model for random graphs and has gained increasing attention in recent years (see [25] for a survey).
We show that is of the order of for probabilities such that for some constant . This result extends the known estimates for the value of the theta number of the random graph .
The paper is organized as follows: Sections 2 and 3 recall basic definitions and properties of simplicial complexes and semidefinite programming. Section 4 recalls properties of the theta number of a graph that serve as a guideline for the theta number of a -dimensional simplicial complex, which is introduced in Section 5. Section 6 computes the theta number of certain basic families of -dimensional simplicial complexes. Section 7 discusses chromatic numbers and Section 8 the hierarchy of theta numbers. The final Section 9 contains the analysis of the theta number of random simplicial complexes.
2. Simplicial complexes
Let be a finite set. We will use the notation for the set of -subsets of . Let us recall that an (abstract) simplicial complex on a vertex set is a family of subsets of (called the faces of ), such that if , then all subsets of also belong to . The dimension of a face is , and we denote by the set of -dimensional faces of , with the convention . Note that we do not require every element in to be a [math]-face of , so can be a proper subset of . The -skeleton of is the simplicial complex .
A simplicial complex is said to be of dimension , if is the maximal dimension of any of its faces. For example, a graph is a simplicial complex of dimension . Going back to the general case, if is of dimension , and if moreover all maximal (with respect to inclusion) faces of are of dimension , then is said to be pure. Unless explicitly mentioned, we will only consider pure complexes.
A basic example of a pure -dimensional simplicial complex is the complete -complex , whose faces are all the subsets of that have at most elements.
We note that in order to define a pure simplicial complex of dimension , it is enough to specify its set of -dimensional faces. In particular, the complementary complex of a pure simplicial complex of dimension , is again a pure simplicial complex of dimension , whose -dimensional faces are those -subsets of that do not belong to (we adopt the convention that the empty complex, whose set of faces is empty, is pure of dimension for all ).
Let be a simplicial complex; we assume that every face of is endowed with an orientation, i.e., a local ordering of its vertices. Then, if and , an oriented incidence number can be defined. Often, the orientation of the faces is induced by a global ordering of the vertex set ; in that case, if where with respect to this ordering,
[TABLE]
The vector space of functions from to is denoted by and its elements are called -dimensional cochains of with coefficients in . The coboundary map is defined for by
[TABLE]
The image of is the subspace of -dimensional coboundaries, and the kernel of is the subspace of -dimensional cocycles. Because the coboundary maps satisfy , we have . The quotient group
[TABLE]
is then called the -th cohomology group of with coefficients in .
Analogously, we can define the homology groups of a simplicial complex. For this, the spaces are endowed with the standard inner product and the boundary map is defined as the adjoint of the coboundary map . We have, for ,
[TABLE]
The spaces of boundaries and of cycles are subspaces of satisfying and thus define the -th reduced homology group of
[TABLE]
Moreover, by duality we have that and . The following diagram summarizes these linear maps for :
[TABLE]
The -th up-Laplacian and -th down-Laplacian of are the following self-adjoint and positive semidefinite operators on :
[TABLE]
By definition, . Furthermore, it is not hard to see that , , , and . For
[TABLE]
we have the Hodge decomposition of into pairwise orthogonal subspaces
[TABLE]
In particular, .
The characteristic functions of faces are called elementary cochains; they form an orthonormal basis of . In order to express the matrices of the Laplacian operators in this basis we introduce the following notation: for , let denote the degree of , i.e., the number of -faces of that contain . For , such that , let
[TABLE]
We note that, if , we can express also as
[TABLE]
For , such that , we set . Then, it is easy to see that
[TABLE]
and
[TABLE]
where we use the same notations for the operators and for their matrices in the basis of elementary cochains.
Example 2.1**.**
In the case of the simplicial complex associated to a graph , defined by , and , we find that is the all-ones matrix and is equal to the combinatorial Laplacian where is the diagonal matrix with the degrees of the vertices as diagonal elements and is the adjacency matrix of the graph.
Example 2.2**.**
For the complete -complex , and for , it is easy to verify that
[TABLE]
Together with the property , we obtain that and that . So is the only non zero eigenvalue of the up and down Laplacians. Computing the traces of these operators gives the multiplicities of this eigenvalue, namely for and for . So we have
[TABLE]
[TABLE]
and, as these dimensions add up to , .
We conclude this section by recalling the definition of the adjacency matrix of a -dimensional simplicial complex : it is the matrix such that where is the diagonal matrix encoding the degrees of the -faces. In other words,
[TABLE]
We note that in dimension this definition coincides with the usual notion of the adjacency matrix of a graph.
3. Semidefinite programming
In this section, we gather basic facts about semidefinite programs. For further information we refer to standard references such as [2], [3], [39].
Semidefinite programs (SDP for short) are special cases of convex optimization programs that admit efficient algorithms, such as algorithms based on the so-called interior point method. They generalize linear programs and have turned out to be very useful for providing polynomial time approximations of hard problems in many areas, especially in combinatorics (see, e.g., [18] and [1, Chapter 6]).
For a matrix we say that is positive semidefinite, denoted by , if is real-valued, symmetric, and if all its eigenvalues are nonnegative. If moreover none of its eigenvalues are equal to zero, is positive definite (). The set of all positive semidefinite matrices is a cone denoted by . The space of real symmetric matrices is endowed with the standard inner product .
Given and symmetric matrices of size , the following optimization problem is a semidefinite program in primal form:
[TABLE]
In other words, this program asks for the supremum of a linear form, where this supremum is taken over the intersection of the cone of positive semidefinite matrices with an affine space.
A feasible solution of this program is a matrix that satisfies the required constraints: and . It is an optimal solution if its objective value is equal to . If there is no feasible solution, we let .
The following dual program is attached to the primal program:
[TABLE]
The terms ’primal’ and ’dual’ do not refer to a specific class of programs: Despite their apparent difference, any of these programs can be put in the form of the other, and, as expected, dualizing twice returns the initial program.
The inequality , referred to as weak duality, always holds, and under some mild conditions even strong duality, i.e., , holds. Strong duality is guaranteed if the SDP satisfies the so-called Slater’s conditions, of which we will use the following version: If an SDP has a strictly feasible primal solution, i.e., if there is a feasible solution of the primal program such that , and a strictly feasible dual solution, i.e., there exists such that , then strong duality holds and, moreover, there are optimal solutions for both the primal and the dual program.
4. The theta number of a graph
In this section, we introduce the theta number of a graph . Our presentation will serve as a guideline for the generalization to higher dimensional simplicial complexes.
Let be an independent set of , i.e., a subset of not containing any edges. The set naturally defines a vector , namely its characteristic vector. We consider the matrix , whose entries are given by:
[TABLE]
The following properties of motivate the definition of : is a positive semidefinite matrix such that if . Furthermore, the cardinality of can be recovered in two different ways from : If and stand as usual for the identity matrix and the all-ones matrix, we have and . So, if we set
[TABLE]
the matrix is feasible for (4) and we get that .
Because (4) is a semidefinite program, its optimal value can be approximated numerically up to arbitrary precision in polynomial time in the size of . If, instead of a sharp numerical value, one aims for a rougher upper bound of , the dual formulation of (4) is often more convenient:
[TABLE]
Here, denotes the largest eigenvalue of .
To illustrate this principle we consider a classical example. For any matrix such that for all , the dual formulation of provides the inequality . A possible choice for is a multiple of the adjacency matrix of , say . The best bound is obtained for minimizing . For -regular graphs, the matrices and commute, so the eigenvalues of are easy to analyze. The optimal choice of then leads to the so-called ratio bound attributed to Hoffman (see, e.g., [4, Theorem 3.5.2]):
[TABLE]
5. The theta number of a simplicial complex
We now move to higher dimensions and define the theta number of a -dimensional simplicial complex . As suggested in the introduction, the down-Laplacian of the complete complex will play the role of the all-ones matrix in (4) and (5). Recall that is the matrix indexed by that is defined by:
[TABLE]
We note that this matrix may not be the down-Laplacian of the complex . Obviously, this is the case if and only if has a complete -skeleton, otherwise the down-Laplacian of is a principal submatrix of . From now on, to avoid confusion, we will denote the matrices associated to by , and reserve the notations , for the complete complex.
Let be an independent set of . Following the same strategy as in the case of graphs, we consider the following matrix , indexed by :
[TABLE]
We have , where as a generalization of the characteristic vector of , we consider the matrix defined as follows:
[TABLE]
where , and is the matrix of the boundary operator with respect to the basis of elementary cochains. The properties of lead to the following definition of :
Definition 5.1**.**
Let be a pure -dimensional complex on , and let be the down Laplacian of the complete complex on . Let:
[TABLE]
Proposition 5.2**.**
We have
[TABLE]
Proof.
Let be an independent set with . As , the matrix is clearly positive semidefinite. We have
[TABLE]
and
[TABLE]
Moreover, from the fact that is an independent set, and from the definition of (7), it is clear that if , or if .
The conditions if are satisfied by the entries of , so the matrix inherits this property.
To sum up, we have proved that the matrix is feasible for . Since its objective value is equal to , we can conclude that . ∎
Now we consider the dual program of (8), in order to obtain another formulation of , similar to (5).
Proposition 5.3**.**
We have
[TABLE]
Proof.
This is just a straightforward rewriting of the dual program. Both programs have the same objective value because Slater’s condition holds: is a strictly feasible solution of (8) and gives rise to a strictly feasible solution of (11). ∎
Remark 5.4**.**
Let us make a few obvious observations about . The first one, is that, as expected, . Indeed, the lower bound follows by taking in (8) while the upper bound follows by taking in (11).
The second observation is that is easy to determine for the empty and the complete -complexes. Indeed, if is the empty -complex, the matrix is feasible for (8) giving that . If is the complete -complex, the semidefinite program (8) has only one feasible solution which is so .
We note that, in these trivial cases, the equality holds.
The benefit of the formulation (11) is that any feasible matrix leads to an upper bound of and therefore to an upper bound of the independence number of . Let us illustrate this principle by showing that we can recover the upper bound proved by Golubev [19] in the case of a -dimensional simplicial complex with complete -skeleton.
We take for some that will be chosen later. Clearly satisfies the conditions required by (11). Then
[TABLE]
We assume that has complete -skeleton, so we have and . Let us denote by the set of non zero eigenvalues of . Then, the eigenvalues of the matrix are: , associated to the eigenspace , and , for , corresponding to eigenvectors in . For , we have and we get:
[TABLE]
We note that, if is regular, i.e., if is a constant number for , then this upper bound is the exact analog of the ratio bound for graphs (6).
We have just seen that, in the case of a -complex with complete -skeleton, is an upper bound of the independence number of which is as least as good as the bound (2). The case of complexes with noncomplete -skeleton turns out to be more tricky; indeed, in some cases provides a good bound of , even a sharp one, and beats the bound (2) given by Golubev, while in other cases, Golubev’s bound is better. We provide examples illustrating this situation in the next section, where we explicitly work out the computation of for certain families of -dimensional complexes. This will also yield counterexamples for certain properties of the theta number related to the chromatic number that we might expect (see Section 7). It will also be interesting to observe the prominent role plaed by the eigenvalues and eigenspaces of the Laplacian operators in these examples .
6. The theta number of certain families of -complexes
6.1. The complete tripartite -complex
To define this complex, we let and partition into three subsets , , of equal size . As -dimensional faces we select all triangles with exactly one vertex in each of these subsets; as -dimensional faces all edges with at most one vertex in each of these subsets. A natural notation for this complex is . It is clear that because is a maximal independent set with vertices. We will show that .
With the notations of (2), , , , and the bound in (2) equals , so this is an example where the theta number beats Golubev’s bound.
We will also show that, for the complementary complex , we have . This complex has a complete -skeleton with and , so Golubev’s bound (2) equals , which is not tight.
Proposition 6.1**.**
We have and .
Proof.
To keep notations light we use the generic notation for throughout the proof. We will verify that , by constructing a suitable matrix feasible for (11). The matrix will be constructed from the projection matrices associated to certain eigenspaces of and .
We denote by the set of edges connecting one vertex in and one vertex in , and similarly for the other kinds of edges. So, . We choose the orientations of the triangular faces and of the edges of following the rule ; this way, for all and .
It turns out that the up-Laplacian has three non zero eigenvalues, , and , respectively with multiplicity , , and . We will need the projection matrices and associated to the eigenvalues and .
The all-one vector is clearly an eigenvector of for the eigenvalue , so . The space is easily seen to be an eigenspace of for the eigenvalue . Similarly, we have two other -dimensional eigenspaces and , and these spaces are pairwise orthogonal. In order to express the projection matrix associated to the sum of these spaces, we introduce the following notation: for , we denote if and both belong to (respectively to , ). Then,
[TABLE]
The down Laplacian has two non zero eigenvalues: with multiplicity and with multiplicity . The vector space is a two-dimensional space of eigenvectors for and for the eigenvalue , and the corresponding projection matrix is given by:
[TABLE]
So far the matrices that we have defined are indexed by . We now will consider matrices indexed by the whole set , therefore we extend the matrices introduced above by adding zero rows and columns for the indices not belonging to (we keep the same notation for the enlarged matrices). We are now ready to define the matrix that will do the job for :
Lemma 6.2**.**
With the previous notations, let
[TABLE]
This matrix satisfies the following properties:
- (1)
* for all * 2. (2)
* for all such that and * 3. (3)
.
Proof.
Properties (1) and (2) follow by direct verification. In order to prove (3), we write where , and , and make the remark that the product of any two of these matrices is zero. Indeed, for and for it follows immediately from the property that the product of up and down Laplacians is zero; for , it is due to the fact that the image of is an eigenspace for the eigenvalue not only for but also for . So, we need to prove that , and are positive semidefinite. For the first two it is obvious because and . So now the only missing piece is a proof that .
For this, we arrange the elements of so that those in come before those in , and we accordingly write by blocks:
[TABLE]
We want to prove that
[TABLE]
By the Schur complement lemma, this is equivalent to . A direct computation shows that , so all boils down to , which is indeed true because is a block-diagonal matrix with three blocks equal to . ∎
Now, we turn our attention to . In order to prove that , we will use the primal formulation (8) and apply a symmetry argument. In the next section we will see a second, simpler, proof, using chromatic numbers, see Example 7.6.
With the previous notations, a feasible matrix must be of the form:
[TABLE]
where is supported on the diagonal and on the triangles that belong to , i.e., the triangles with one vertex in each of , , . It is clear that the automorphism group of permutes transitively the elements of and of , and that, by convexity, (8) has a symmetric solution. So, without loss of generality, we can assume that . Restricting the semidefinite program on this set of matrices leads to a linear program in the variables , , that can be easily solved and leads to the optimal value . We skip the details here.
We note that this approach would not work for because has two orbits: the triangles that are fully contained in one of the subsets , , and the ones that have two vertices in one of these sets and one vertex in another one. ∎
6.2. The complete bipartite -complex
Now and is partitioned in two subsets , , of equal size . As -dimensional faces we select the triangles that meet both sets and , thus having two vertices in one of the parts and the third vertex in the other. We denote this complex by . It is clear that since is an independent set with vertices. This complex has a complete -skeleton and , so the bound (3) equals , showing that and that the theta number agrees with Golubev’s bound.
For the complementary complex , which is nothing else than the disjoint union of two complete complexes , we have . Golubev’s bound is twice the value corresponding to , thus , and it is sharp again. As we will see know, is much larger:
Proposition 6.3**.**
We have and .
Proof.
We let . To compute , we again apply the symmetry principle, like in the case of the complement of the tripartite complex. The automorphism group of has two orbits in : the set of edges contained in or in , having degree , and the set of ’crossing’ edges, with degree . It acts transitively on the -faces. So without loss of generality a feasible matrix of the primal formulation of can be assumed to be
[TABLE]
where and denote the diagonal matrices associated to respectively and . The expressions of and of are linear in the variables , but the condition that is positive semidefinite is slightly more complicated because does not commute with and . In fact, this condition leads to quadratic constraints, as it will become clear if we write the matrices by blocks according to . It is easy to verify that
[TABLE]
and that has two non zero eigenvalues: , with multiplicity and eigenvector the all-one vector, and , with multiplicity . Then, by the Schur complement lemma, the condition
[TABLE]
leads to quadratic inequalities. It is a bit technical but not difficult to see that an optimal solution satisfies , and finally that it is
[TABLE]
leading to the optimal value . ∎
7. Chromatic numbers
Let us first review the case of graphs. For a graph , the clique number and the chromatic number are related by the obvious inequality , and the theta number lies in between these numbers ([32, Lemma 3, Corollary 3]):
[TABLE]
Moreover, the inequality is always at least as strong as the inequality ; indeed, we know that from [32, Corollary 2].
Let us consider the situation for pure -dimensional simplicial complexes. By analogy with graphs, the chromatic number of a complex , is usually defined to be the least number of colors needed to color the vertices of such that no -face is monochromatic. We remark that for the complete -complex , the color classes of an admissible coloring cannot have more than elements, and consequently that . So, for all -dimensional complexes , we have . Given that we have defined a generalization of the theta number to -complexes, that satisfies , it is natural to wonder if the inequality
[TABLE]
is also satisfied. Unfortunately, this is not true in general. Indeed, from the results of Section 6, one can see that (13) is satisfied for the complete tripartite complex and for its complement, but fails for the complete bipartite complex , for which (Proposition 6.3) while .
Let us now see if we can modify the definition of the chromatic number of a simplicial complex, so that it fits better with our theta number. To achieve this, we will adapt the concept of graph homomorphisms to simplicial complexes. Indeed, a nice way to understand the notions of chromatic and clique numbers of graphs is through their connection to graph homomorphisms, as we will recall now.
A homomorphism from a graph to a graph is a mapping from the vertices of to the vertices of that sends an edge of to an edge of . Then, the clique number and the chromatic number have the following interpretations: the clique number is the largest number such that there is a homomorphism from the complete graph to , and similarly is the smallest number such that there is a homomorphism from to . Moreover, one can prove that, if there is a homomorphism from to , then . The combination of these properties immediately leads to (12).
In order to follow a similar approach for simplicial complexes, we introduce an ad-hoc notion of homomorphism.
Definition 7.1**.**
Let and be two pure -dimensional simplicial complexes. A homomorphism from to is a mapping with the following property: There exist orientations of and such that for every , there is such that
- (1)
, 2. (2)
* for all with .*
We note that this definition coincides in dimension with the usual notion of a graph homomorphism as one can always find suitable orientations.
Remark 7.2**.**
In this definition, it is important to understand that a homomorphism may not necessarily be induced by a global mapping between the vertices, i.e., it may be the case that there is no mapping such that for all . As an example consider the -dimensional complex depicted in Figure 1.
Furthermore, condition (2) is not automatically fulfilled. The -dimensional complex depicted in Figure 2 possesses a map satisfying condition (1) but there is no homomorphism from to .
Proposition 7.3**.**
Let and be two pure -dimensional simplicial complexes, and let be a homomorphism from to . Then,
[TABLE]
Proof.
Our strategy will be to start with an optimal solution of the primal formulation (8) of , from which we construct a matrix , feasible for , and having the same objective value as .
So, let be primal optimal for the semidefinite program defining . We remark that, if , then, for all , , and so . As a consequence, by the optimality of , we have .
For , we set
[TABLE]
where the sum is zero if or does not belong to the image of .
We have .
By the property 1) of homomorphisms, if and is not an element of , and if and , then cannot belong to , and so . So, we have that .
Thanks to property 2), if and , the required condition that holds. So, we have proved that is primal feasible for .
It remains to analyze the objective value . We have
[TABLE]
But
[TABLE]
where in the last equality we ignore the terms corresponding to because they are equal to zero, and we apply the property 2). It follows that . ∎
Definition 7.4**.**
Let be a pure -dimensional simplicial complex. Let denote the smallest number such that there exists a homomorphism from to the complete -complex .
It is not hard to see that holds for any pure simplicial complex as a vertex coloring with colors that is a proper graph coloring for gives rise to a homomorphism from to . The complex depicted in Figure 1 serves as an example that the three notions of chromatic numbers considered here differ. It has , and .
Proposition 7.5**.**
We have
[TABLE]
Proof.
If there is then applying (14) leads to (see Remark 5.4). ∎
Example 7.6**.**
Consider the complex defined in Section 6. Clearly, , so we have and hence .
A -dimensional subcomplex of a pure -dimensional simplicial complex is a connected component of if for every -face of any -face of that contains is also in . Note that this condition does not need to hold for lower dimensional simplices, so two distinct connected components can, e.g., share a common vertex. Further observe that the connected components of correspond to the connected components of the graph that has the -faces of as vertices with two vertices forming an edge if the correponding -faces intersect in a common -face.
As different connected components do not share -faces, the inequality can actually be extended to the connected components of .
Proposition 7.7**.**
Let be the collection of connected components of . Then
[TABLE]
It is well-known that a -regular graph has a bipartite connected component if and only if the largest eigenvalue of the Laplacian is . In [23] Horak and Jost present a combinatorial criterion that can be considered as a higher-dimensional analog of this: They show that for a -regular -complex the largest eigenvalue of the Laplacian is if and only if there is a connected component of and an orientation of the -faces of such that for all , . Note that for a connected graph the existence of such an orientation is equivalent to bipartiteness.
If a -dimensional simplicial complex has chromatic number , this guarantees the existence of such an orientation. Hence, we have the following observation.
Proposition 7.8**.**
Let be a -regular -dimensional simplicial complex. If , then the maximal eigenvalue of the up-Laplacian is .
We remark that these results extend to arbitrary complexes for a normalized version of the Laplacian that we do not study here.
8. A hierarchy of semidefinite relaxations for the independence
number of a -simplicial complex
In this section, is again a pure -dimensional simplicial complex. We consider a straightforward generalization of that leads to higher order theta numbers for . We will see that all these numbers provide upper bounds of , until , where . Finally, we will modify this sequence of theta numbers in order to get a decreasing sequence.
It will be convenient to denote by the set of independent sets of dimension . We make the remark that is a simplicial complex, the independence complex of , and that it has complete -skeleton, i.e., . For , the matrices involved in the program defining are indexed by . We define, for :
[TABLE]
and its dual formulation:
[TABLE]
The above definition matches for with that of . Both primal and dual programs are strictly feasible: and respectively give rise to strictly feasible solutions. We note that, if , the feasible matrices of the primal program are diagonal matrices and hence . We have
Proposition 8.1**.**
[TABLE]
Proof.
The same proof as the one of Proposition 5.2 works. For an independent set such that , we define by
[TABLE]
It is then easy to verify, as every subset of an independent set is also an independent set, that is feasible for the primal program (15) and that its objective value is equal to . ∎
However, it is not clear that the sequence is decreasing, because the constraints on the -sets involved in do not occur explicitly in . We now define a variant of that provides a decreasing sequence of upper bounds of .
To start with, we note that, if a matrix is feasible for (15), then the value of for such that only depends on . So, we can associate to a function such that if . If we extend to by , we see that encodes every nonzero entry of . Said differently, we have a one to one correspondence between and the set
[TABLE]
We record for later use that, if corresponds to as above, then
[TABLE]
and
[TABLE]
Now, we introduce, for , a map . It will be more convenient to define on the corresponding functions , in the following way: let
[TABLE]
where
[TABLE]
We are now in the position to define our strengthening of : Let
[TABLE]
Theorem 8.2**.**
The numbers , , satisfy:
- (1)
** 2. (2)
.
Proof.
That is clear since we have only added constraints on in the definition of .
Let be an independent set, with . Let, like in the proof of Proposition 8.1, be defined by:
[TABLE]
The element corresponding to is given by: if , if , and otherwise takes the value [math]. We will need the following lemma:
Lemma 8.3**.**
We have
[TABLE]
for as defined in (20).
Proof.
Let . Let . Every subset of is independent so the number of such that is . So,
[TABLE]
Now let . It is clear that, if is not contained in , . If ,
[TABLE]
∎
Lemma 8.3 shows that is positive semidefinite, and so, iteratively, that is positive semidefinite for every . We conclude that (after a suitable rescaling) is feasible for , and consequently that . We have already remarked that so also .
It remains to prove that the sequence of is decreasing. For this, we start from an optimal solution of , and we show that is feasible for and that .
It is clear that and that is positive semidefinite, as well as for all . That follows easily from (17) and from the definition of . It remains to take care of the objective value. Applying (18),
[TABLE]
where in the sums we restrict to elements in . Taking account of the fact that every subset of an independent set is also an independent set, we obtain
[TABLE]
∎
9. Theta numbers of random complexes
A random model for simplicial complexes of arbitrary fixed dimension was introduced by Linial and Meshulam [31] as a higher dimensional analog of the Erdös-Rényi model for random graphs. It has vertex set , complete -skeleton, and each element of is added as a -dimensional face of independently with probability . Here is a function of , and we let . In this section we analyze the theta number of for ’dense’ complexes, i.e., for in the range .
The study of the theta number of random graphs was initiated by Juhász in [24] who proved that, in the case of constant probability , holds with probability tending to . In subsequent works, the range of probabilities for which Juhász’ result holds was extended, until in [6], Coja-Oghlan was able to cover for some sufficiently large constant .
We will restrict ourselves to the range because we will need the following estimates:
Theorem 9.1** ([16, 22]).**
Let denote the adjacency matrix of . For every there exists such that, if ,
[TABLE]
and
[TABLE]
with probability at least equal to .
With the above, it is rather straightforward to obtain:
Theorem 9.2**.**
For every there exists such that, if ,
[TABLE]
with probability at least equal to .
Indeed, following the method of Juhász, the upper bound is obtained via the dual formulation for the theta number (5) and the matrix , where is the adjacency matrix of , while the lower bound follows from the choice in the primal formulation (4), where , being the adjacency matrix of the complementary graph of .
9.1. The theta number of
We will establish the following similar result for random simplicial complexes :
Theorem 9.3**.**
For every and , there exists such that, if ,
[TABLE]
with probability at least equal to .
For comparison, the independence number of is of the order (see [28]). In the range , the eigenvalues of the adjacency matrix of have been studied in [21]. We will closely follow the methods developed in [21], in particular the role played by the so-called links of , an idea going back to the work of Garland [17]. By definition, for a -dimensional simplicial complex and a -face of , the link is the graph with vertices , and edges . In view of the proof of Theorem 9.3, we will first establish a relationship between the theta number of a simplicial complex and that of its links.
Proposition 9.4**.**
Let be a -dimensional simplicial complex with complete -skeleton. Then
[TABLE]
Proof.
Let . For a matrix , we introduce its localization at denoted and defined by:
[TABLE]
Let denote the diagonal matrix with as diagonal entries. Then we observe that
[TABLE]
and that, if for all such that ,
[TABLE]
Now let be an optimal solution of (8). Taking account of (25) and (26),
[TABLE]
If , we have
[TABLE]
so, since ,
[TABLE]
Now, the crucial observation is that the matrix gives rise to a feasible matrix of the semidefinite program (4) defining the theta number of . Indeed, let be the matrix indexed by and defined by . This matrix inherits some properties of : The matrix is positive semidefinite, the entries of associated to edges of are equal to [math]. With obvious notations, we have and so we obtain
[TABLE]
We have so the announced inequality follows immediately. ∎
Proof of Theorem 9.3.
For the upper bound, we apply Proposition 9.4. The link of a ()-face in a random complex is an Erdös-Renyi random graph on with the same probability . We can thus apply Theorem (9.2) and a union bound to obtain the result. We note that, since the number of such faces is of the order of , for the probability of the bad event to be, say, less than we need to apply Theorem (9.2) for the larger value instead of , explaining the need for an arbitrary large power of in the convergence speed of probabilities.
In order to find a lower bound of , we consider the matrix where denotes the adjacency matrix of the complementary -complex . The feasibility conditions of (8) are fulfilled by except for the normalization condition . We have . Moreover, , so
[TABLE]
The number of -faces of is a random variable binomially distributed in with probability . Hence, by a straightforward application of a Chernoff bound, for every , is at least of the order with probability at least . It remains to upper bound . For this, we apply the localization procedure that we have already encountered in the proof of Proposition 9.4:
[TABLE]
Then, for every , if denotes the vector obtained from by setting to [math] the coordinates of associated to faces not containing ,
[TABLE]
The matrix has the same spectrum as . The latter is identical to the adjacency matrix of the graph on the entries indexed by , and zero elsewhere. So, its non-zero spectrum is that of and hence:
[TABLE]
The links are random graphs so, applying (22) and a union bound, we find that, with probability at least equal to , for a large enough constant ,
[TABLE]
We have obtained the desired upper bound . Putting everything together, we obtain the announced lower bound for . ∎
9.2. The hierarchy of theta numbers of
In this last subsection, we restrict ourselves to the case of random graphs and analyze the hierarchy of theta numbers for constant values of . The restriction to random graphs, i.e., random complexes of dimension , is purely for simplicity. The assumption of constant , however, is essential. Analyzing the complete hierarchy of a random complex for non-constant appears to be a difficult task. It would be interesting to know for which values of the theta number is close to the independence number. Unfortunately, such questions seem to be out of the reach of the methods we apply here.
Theorem 9.5**.**
For every and , there exists such that, if and ,
[TABLE]
with probability at least equal to .
Proof.
We will sometimes use the expression with high probability for an inequality that holds with probability at least for all , with appropriate constants depending on .
For an upper bound of , we apply
[TABLE]
Here, is the graph on with edges if . If is independent, this condition simply means that is an edge of , so is the graph induced by on . If , the number of vertices is itself a random variable. Since , follows a binomial distribution with parameters and . For to be concentrated around its expected value we need for some .
Assuming for some , we have
[TABLE]
because can be viewed as an induced subgraph of . We would like to apply Theorem 9.2. It requires and to be greater that and holds with probability at least . All this will be fine if we assume:
[TABLE]
for a sufficiently large . With a union bound we obtain with high probability:
[TABLE]
For the lower bound, we consider the matrix where is the adjacency matrix of the -skeleton of and we apply (15). We obtain
[TABLE]
In order to estimate we use and remark that has the same non-zero eigenvalues as the adjacency matrix of the graph , itself being the graph induced by on . We have
[TABLE]
so
[TABLE]
Like for the upper bound we have with high probability for some and thus
[TABLE]
for some , under the same conditions on and .
It remains to deal with the ratio . For this we will argue that is almost regular. To be more precise we apply double counting to the set
[TABLE]
The number of -subsets of is so . For a given , the number of containing follows a binomial distribution with parameters and , with expected value . With high probability (requires ) is larger that and so
[TABLE]
Putting everything together and applying another union bound we obtain
[TABLE]
∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] M. Anjos and J.B. Lasserre, Handbook on Semidefinite, Conic and Polynomial Optimization , Springer, 2012
- 2[2] A. Ben-Tal and A. Nemirovski, Lectures on Modern Convex Optimization: Analysis, Algorithms, and Engineering Applications , SIAM, Philadelphia, 2001.
- 3[3] S. Boyd and L. Vandenberghe Convex Optimization , Cambridge University Press, 2004.
- 4[4] A.E. Brouwer and W.H. Haemers, Spectra of Graphs , Springer, New York (2012).
- 5[5] H. Cohn, A. Kumar, S.D. Miller, D. Radchenko, and M.S. Viazovska, The sphere packing problem in dimension 24 , preprint, ar Xiv:1603.06518 [math.NT], 2016, 12pp.
- 6[6] A. Coja-Oghlan, The Lovász number of random graphs , Combinatorics, Probability and Computing 14 (2005), 439-465.
- 7[7] P.E.B. De Corte, D. de Laat, and F. Vallentin, Fourier analysis on finite groups and the Lovász theta-number of Cayley graphs , Experimental Mathematics 23 (2014), 146-152.
- 8[8] P. Delsarte, An algebraic approach to the association schemes of coding theory , Philips Research Repts Suppl. 10, (1973), 1-97.
