This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Conformally rigid graphs

Stefan Steinerberger Department of Mathematics, University of Washington, Seattle, WA 98195, USA steinerb@uw.edu  and  Rekha R. Thomas Department of Mathematics, University of Washington, Seattle, WA 98195, USA rrthomas@uw.edu
Abstract.

Given a finite, simple, connected graph G=(V,E)G=(V,E) with |V|=n|V|=n, we consider the associated graph Laplacian matrix L=DAL=D-A with eigenvalues 0=λ1<λ2λn0=\lambda_{1}<\lambda_{2}\leq\dots\leq\lambda_{n}. One can also consider the same graph equipped with positive edge weights w:E>0w:E\rightarrow\mathbb{R}_{>0} normalized to eEwe=|E|\sum_{e\in E}w_{e}=|E| and the associated weighted Laplacian matrix LwL_{w}. We say that GG is conformally rigid if constant edge-weights maximize the second eigenvalue λ2(w)\lambda_{2}(w) of LwL_{w} over all ww, and minimize λn(w)\lambda_{n}(w^{\prime}) of LwL_{w^{\prime}} over all ww^{\prime}, i.e., for all w,ww,w^{\prime},

λ2(w)λ2(𝟙)λn(𝟙)λn(w).\lambda_{2}(w)\leq\lambda_{2}(\mathbbm{1})\leq\lambda_{n}(\mathbbm{1})\leq\lambda_{n}(w^{\prime}).

Conformal rigidity requires an extraordinary amount of structure in GG. Every edge-transitive graph is conformally rigid. We prove that every distance-regular graph, and hence every strongly-regular graph, is conformally rigid. Certain special graph embeddings can be used to characterize conformal rigidity. Cayley graphs can be conformally rigid but need not be, we prove a sufficient criterion. We also find a small set of conformally rigid graphs that do not belong into any of the above categories; these include the Hoffman graph, the crossing number graph 6B and others. Conformal rigidity can be certified via semidefinite programming, we provide explicit examples.

1. Introduction

1.1. Introduction

The purpose of this paper is to introduce a property, which we call conformal rigidity, that a finite, simple graph G=(V,E)G=(V,E) may or may not have. It is the discrete analogue of a property of interest in spectral geometry and can only be expected in the presence of structure. Roughly put, it says that among all possible ways one could weigh the edges, constant weight edges are already, in a suitable sense, the canonical choice. We start by defining things more formally. Let G=(V,E)G=(V,E) be a finite, simple, connected graph on nn vertices. A classical object associated to GG is its graph Laplacian L=DAL=D-A, where Dn×nD\in\mathbb{R}^{n\times n} is the diagonal matrix collecting the degrees of the vertices and A{0,1}n×nA\in\left\{0,1\right\}^{n\times n} is the adjacency matrix of GG. The matrix LL is positive semidefinite (psd) and has eigenvalues

0=λ1<λ2λn0=\lambda_{1}<\lambda_{2}\leq\dots\leq\lambda_{n}

which are fundamental objects in spectral graph theory. The theory generalizes to the case where edges are weighted wij>0w_{ij}>0; the weighted Laplacian Lwn×nL_{w}\in{\mathbb{R}}^{n\times n} is the psd matrix defined as:

(4) (Lw)ij:={wij if (i,j)Ejiwij in the (i,i) entry0 otherwise \displaystyle(L_{w})_{ij}:=\left\{\begin{array}[]{rl}-w_{ij}&\textup{ if }(i,j)\in E\\ \sum_{j\sim i}w_{ij}&\textup{ in the }(i,i)\textup{ entry}\\ 0&\textup{ otherwise }\end{array}\right.

with eigenvalues 0=λ1(w)<λ2(w)λn(w)0=\lambda_{1}(w)<\lambda_{2}(w)\leq\dots\leq\lambda_{n}(w). We use jij\sim i to denote that (i,j)E(i,j)\in E. We only work with undirected graphs, so there is no difference between jij\sim i and iji\sim j. The constant vector 𝟙:=(1,1,,1)\mathbbm{1}:=(1,1,\ldots,1) spans the eigenspace of λ1(w)=0\lambda_{1}(w)=0. By the Rayleigh-Ritz Theorem [15, §4.2], we can write λ2(w)\lambda_{2}(w) as

λ2(w)=minvVf(v)=0(i,j)Ewij(f(i)f(j))2vVf(v)2\lambda_{2}(w)=\min_{\sum_{v\in V}f(v)=0}\frac{\sum_{(i,j)\in E}w_{ij}(f(i)-f(j))^{2}}{\sum_{v\in V}f(v)^{2}}

where the minimum ranges over all f:Vf:V\rightarrow\mathbb{R} that are not identically 0. The condition vVf(v)=0\sum_{v\in V}f(v)=0 forces orthogonality to the constant vector 𝟙\mathbbm{1}. Likewise, and without orthogonality conditions, we have

λn(w)=maxf:V(i,j)Ewij(f(i)f(j))2vVf(v)2.\lambda_{n}(w)=\max_{f:V\rightarrow\mathbb{R}}\frac{\sum_{(i,j)\in E}w_{ij}(f(i)-f(j))^{2}}{\sum_{v\in V}f(v)^{2}}.

It is well known that GG is connected if and only if λ2(w)>0\lambda_{2}(w)>0. In fact, λ2(w)\lambda_{2}(w) measures the connectivity of a graph: larger values correspond to a more connected graph. The largest eigenvalue λn(w)\lambda_{n}(w) is more intricate but corresponds roughly to the largest possible oscillation that can occur among the functions of the graph.

Refer to captionRefer to captionRefer to caption
Figure 1. Three conformally rigid graphs. Left: Hoffman graph. Middle: complement of the Shrikhande graph. Right: CNG 6B.

It is a natural question how, for any given graph G=(V,E)G=(V,E), one should choose weights ww to maximize λ2(w)\lambda_{2}(w) or to minimize λn(w)\lambda_{n}(w), respectively. Since the eigenvalues scale linearly with ww, we require both wij0w_{ij}\geq 0 and the normalization

(i,j)Ewij=|E|\sum_{(i,j)\in E}w_{ij}=|E|

which we do throughout the paper. Maximizing λ2(w)\lambda_{2}(w) can be understood in several ways, including creating the ‘densest’ network, while minimizing λn(w)\lambda_{n}(w) is akin to a smoothness condition on the weights. It is a way of ‘regularizing’ the graph, finding its most symmetric representation. We are interested in graphs for which these optimization procedures cannot improve on the unweighted case.

Definition 1.1.

We say a graph G=(V,E)G=(V,E) is conformally rigid if, for all edge weights w,w:E0w,w^{\prime}:E\rightarrow\mathbb{R}_{\geq 0} normalized to (i,j)Ewij=|E|\sum_{(i,j)\in E}w_{ij}=|E| and (i,j)Ewij=|E|\sum_{(i,j)\in E}w_{ij}^{\prime}=|E|,

λ2(w)λ2(𝟙)λn(𝟙)λn(w).\lambda_{2}(w)\leq\lambda_{2}(\mathbbm{1})\leq\lambda_{n}(\mathbbm{1})\leq\lambda_{n}(w^{\prime}).
λ2(w)\lambda_{2}(w)λ2(𝟙)\lambda_{2}(\mathbbm{1})λn(𝟙)\lambda_{n}(\mathbbm{1})λn(w)\lambda_{n}(w^{\prime})
Figure 2. If a graph is conformally rigid, then the the interval spanned by the spectrum of the Laplacian, for any assignment of weights (fixing the total sum), always contains [λ2(𝟙),λn(𝟙)][\lambda_{2}(\mathbbm{1}),\lambda_{n}(\mathbbm{1})].

This corresponds, in a suitable sense, to graphs which are already so symmetric that varying edge weights cannot further increase λ2\lambda_{2} or decrease λn\lambda_{n}. The definition can be phrased in terms of the spectrum of the Laplacian matrix as in Fig. 2. An alternative viewpoint is as follows: by changing the weights, we change the eigenvalues. As mentioned already, there is a linear scaling, doubling all the weights doubles all the eigenvalues, which makes it natural to normalize the sum of the weights. Conformal rigidity is concerned with the setting where the largest attainable second eigenvalue over all possible weights is simply the eigenvalue of the unweighted graph itself, meaning

supw>0λ2(w)(i,j)Ewij=supw>0inff:VvVf(v)=0(i,j)Ewij(f(i)f(j))2((i,j)Ewij)(vVf(v)2)=λ2(𝟙)\sup_{w>0}\frac{\lambda_{2}(w)}{\sum_{(i,j)\in E}w_{ij}}=\sup_{w>0}\inf_{f:V\rightarrow\mathbb{R}\atop\sum_{v\in V}f(v)=0}\frac{\sum_{(i,j)\in E}w_{ij}(f(i)-f(j))^{2}}{\left(\sum_{(i,j)\in E}w_{ij}\right)\left(\sum_{v\in V}f(v)^{2}\right)}=\lambda_{2}(\mathbbm{1})

and where, simultaneously, the same is true for the largest eigenvalue λn(w)\lambda_{n}(w^{\prime})

infw>0λn(w)(i,j)Ewij=infw>0supf:V(i,j)Ewij(f(i)f(j))2((i,j)Ewij)(vVf(v)2)=λn(𝟙).\inf_{w^{\prime}>0}\frac{\lambda_{n}(w^{\prime})}{\sum_{(i,j)\in E}w^{\prime}_{ij}}=\inf_{w^{\prime}>0}\sup_{f:V\rightarrow\mathbb{R}}\frac{\sum_{(i,j)\in E}w^{\prime}_{ij}(f(i)-f(j))^{2}}{\left(\sum_{(i,j)\in E}w^{\prime}_{ij}\right)\left(\sum_{v\in V}f(v)^{2}\right)}=\lambda_{n}(\mathbbm{1}).

These equations show the implicit mathematical structure of the question: there is an entire additional optimization problem sitting on top of the Rayleigh-Ritz formulations of λ2\lambda_{2} and λn\lambda_{n} which are optimization problems in themselves.

Let us quickly consider some examples. It may not come as a surprise that the complete graph KnK_{n} is conformally rigid. An elementary proof can be found in § 5.1.

Proposition 1.2.

The complete graph KnK_{n} is conformally rigid. Moreover, the only ww for which λ2(w)=λ2(𝟙)\lambda_{2}(w)=\lambda_{2}(\mathbbm{1}) and λn(𝟙)=λn(w),\lambda_{n}(\mathbbm{1})=\lambda_{n}(w^{\prime}), is w=w=𝟙w=w^{\prime}=\mathbbm{1}.

Conformal rigidity does not imply that 𝟙\mathbbm{1} is the unique weight vector that maximizes λ2(w)\lambda_{2}(w) and minimizes λn(w)\lambda_{n}(w^{\prime}): a conformally rigid graph GG may have non-constant weights w,ww,w^{\prime} with λ2(w)=λ2(𝟙)\lambda_{2}(w)=\lambda_{2}(\mathbbm{1}), an example being a circulant graph in Example 2.4, or λn(w)=λn(𝟙)\lambda_{n}(w^{\prime})=\lambda_{n}(\mathbbm{1}) with an example being the Haar Graph 565, §6.3.

1.2. Motivation

Our definition of conformal rigidity in graphs has a number of different motivations, we describe three below.

1.2.1. Conformal Optimization.

The first motivation comes from (continuous) spectral geometry. Let (M,g)(M,g) be a smooth, compact Riemannian manifold of dimension n2n\geq 2. The Laplace-Beltrami operator Δ-\Delta has a discrete (infinite) sequence of eigenvalues 0=λ1<λ20=\lambda_{1}<\lambda_{2}\leq\dots and it is a classical problem to understand which manifold minimizes or maximizes λk\lambda_{k} under suitable conditions. Let us now consider changes of the underlying metric gg and introduce the conformal class of the metric

[g]={h=e2ug,uC(M)}.[g]=\left\{h=e^{2u}g,u\in C^{\infty}(M)\right\}.

This corresponds to local rescalings of the metric that preserve angles: one can blow up certain parts of the manifold and one can shrink other parts but one cannot change the local conformal structure. We can now define

Λk(M,[g])=suph[g]λk(h)volh(M)2n.\Lambda_{k}(M,[g])=\sup_{h\in[g]}\lambda_{k}(h)\cdot\mbox{vol}_{h}(M)^{\frac{2}{n}}.

The volume term counters the fact that this procedure may increase or decrease volume (and Laplacian eigenvalues are sensitive to that). Alternatively, one could only consider metrics for which volh(M)2n\mbox{vol}_{h}(M)^{\frac{2}{n}} is of a fixed size; this would correspond to our edge normalization (i,j)Ewij=|E|\sum_{(i,j)\in E}w_{ij}=|E|. Korevaar [20] proved that this supremum is finite for each value of kk. Very much as in our setting, metrics that attain the supremum would correspond to the ‘most symmetric’ realization of the manifold in its conformal class. As a motivating result, we recall a classic result of Hersch [18] and El Soufi - Ilias [7] showing that

Λ2(𝕊n,[dx])=nωn2n,\Lambda_{2}(\mathbb{S}^{n},[dx])=n\cdot\omega_{n}^{\frac{2}{n}},

where ωn\omega_{n} is the volume of the standard sphere and equality is attained if and only if the sphere is round. Phrased differently, the round sphere is already the most symmetric realization. Conformally rigid graphs are, in this sense, the discrete analogues of spheres in this context: their unweighted version is already spectrally extremal among all possible ways of distributing weights over the edges.

1.2.2. Frame Bounds.

Another reason for why we are interested in conformally rigid graphs is the following basic fact corresponding to the spectrum of the complete graph KnK_{n}: given any list of nn numbers x1,,xnx_{1},\dots,x_{n}\in\mathbb{R}, we have

i=1nxi=01i<jn(xixj)2=ni=1nxi2.\sum_{i=1}^{n}x_{i}=0\implies\sum_{1\leq i<j\leq n}(x_{i}-x_{j})^{2}=n\sum_{i=1}^{n}x_{i}^{2}.

One way of phrasing this result is as follows: it is possible deduce the 2\ell^{2}-norm of (x1,,xn)(x_{1},\dots,x_{n}) just from knowing all the pairwise differences xixjx_{i}-x_{j} (up to, of course, additive constants which cannot be recovered). There is a natural question: if we no longer have access to all pairwise differences but only some (encoded by a graph), how accurately can we recover the norm? This is determined by the spectrum

i=1nxi=0λ2i=1nxi2(i,j)E(xixj)2λni=1nxi2,\sum_{i=1}^{n}x_{i}=0\implies\qquad\lambda_{2}\sum_{i=1}^{n}x_{i}^{2}\leq\sum_{(i,j)\in E}(x_{i}-x_{j})^{2}\leq\lambda_{n}\sum_{i=1}^{n}x_{i}^{2},

and can be interpreted as frame bounds, with λn/λ2\lambda_{n}/\lambda_{2} being the condition number. If we are restricted to working with a fixed set of differences, a fixed EE, then it is a natural question whether one can improve the accuracy of recovery by adding weights. Conformally rigid graphs have the property that allowing weights can neither increase λ2\lambda_{2} nor decrease λn\lambda_{n}, the graphs are already optimal in this regard. One way of stating this precisely is via the following Proposition which follows immediately from the previously mentioned sup/inf characterization.

Proposition 1.3.

Let G=(V,E)G=(V,E) be conformally rigid. Then, for all possible (positive) edge weights wijw_{ij} there always exists a function f:Vf:V\rightarrow\mathbb{R} with mean value 0, meaning vVf(v)=0\sum_{v\in V}f(v)=0, which oscillates little in the sense of

(i,j)Ewij(f(i)f(j))2λ2(𝟙)|E|((i,j)Ewij)(vVf(v)2).\sum_{(i,j)\in E}w_{ij}(f(i)-f(j))^{2}\leq\frac{\lambda_{2}(\mathbbm{1})}{|E|}\left(\sum_{(i,j)\in E}w_{ij}\right)\left(\sum_{v\in V}f(v)^{2}\right).

Moreover, there exists a function f:Vf:V\rightarrow\mathbb{R} which oscillates a lot insofar as

(i,j)Ewij(f(i)f(j))2λn(𝟙)|E|((i,j)Ewij)(vVf(v)2).\sum_{(i,j)\in E}w_{ij}(f(i)-f(j))^{2}\geq\frac{\lambda_{n}(\mathbbm{1})}{|E|}\left(\sum_{(i,j)\in E}w_{ij}\right)\left(\sum_{v\in V}f(v)^{2}\right).

Note that when all the weights are wij=1w_{ij}=1, then the Proposition is merely a form of the classical Minimax formulation for eigenvalues. The new ingredient, given by conformal rigidity, is that we allow for general weights.

1.2.3. Sparsification.

A final motivation comes from a recently proposed framework for Graph Sparsification due to Babecki and the authors [1]. In this framework, the hypercube graph cannot be further sparsified. We were wondering how general this result is which is what lead us to the notion of conformally rigid graphs. If (up to scaling) 𝟙\mathbbm{1} is the unique weight that maximizes λ2(w)\lambda_{2}(w), then GG cannot be sparsified in any setting of [1]. This led us to wonder if all Cayley graphs have this property. As we will see in this paper, not all Cayley graphs have this property. This line of investigation led us to the notion of conformal rigidity. It is an interesting question to understand which conformally rigid graphs have unique maximizers for one or both of λ2(w)\lambda_{2}(w) and λn(w)\lambda_{n}(w).

1.3. Related results

In this section we collect together various results in the literature that are close to the notion of conformal rigidity.

1.3.1. The optimization perspective

The problem of maximizing λ2(w)\lambda_{2}(w) subject to wij=|E|\sum w_{ij}=|E| appears to have been considered first by Fiedler [9, 10] and references therein. Fiedler calls the maximum value of λ2(w)\lambda_{2}(w) the absolute connectivity of the graph. Finding weights that maximize λ2(w)\lambda_{2}(w) subject to a general weighted inequality of the form dij2wij1\sum d_{ij}^{2}w_{ij}\leq 1 underlies the problem of finding the fastest mixing Markov process on GG. This problem was studied by Sun-Boyd-Xiao-Diaconis [27] where they modeled it as a semidefinite program which is a type of convex optimization problem. We will expand on this in § 3. Boyd-Diaconis-Parrilo-Xiao [3] exploited the symmetries of the graph to speed up computations. We will adopt the methods in both [3] and [27] to certify the conformal rigidity of graphs. The semidefinite program viewpoint of maximizing λ2(w)\lambda_{2}(w) was also studied in [12] where they connect it to the problem of separators in GG. Both [12] and [27] interpret this problem as a graph embedding problem. The problem of minimizing λn(w)\lambda_{n}(w) is also a semidefinite program with interpretations in terms of graph embeddings and partitions. The difference λn(w)λ2(w)\lambda_{n}(w)-\lambda_{2}(w) of the Laplacian LwL_{w} is called the spectral width of LwL_{w}. The problem of minimizing the spectral width of GG is

(5) min{λn(w)λ2(w):wij=|E|,w0},\displaystyle\min\left\{\lambda_{n}(w)-\lambda_{2}(w)\,:\,\sum w_{ij}=|E|,\,\,\,w\geq 0\right\},

and is also a semidefinite program with connections to graph embeddings and separators [14]. The spectral width bounds the uniform sparsest cut in GG since

1nλ2(w)w(S,VS)|S||VS|1nλn(w)\frac{1}{n}\lambda_{2}(w)\leq\frac{w(S,V\setminus S)}{|S||V\setminus S|}\leq\frac{1}{n}\lambda_{n}(w)

where w(S,VS)w(S,V\setminus S) is the weight of the cut in GG given by the partition V=S(VS)V=S\cup(V\setminus S).

Remark. The spectral width minimization problem looks similar to conformal rigidity and so we pause to point out an important difference. In our setup we are solving

(6) min{λn(w)λ2(w): 1w=𝟙w=|E|,w,w0}.\displaystyle\min\left\{\lambda_{n}(w)-\lambda_{2}(w^{\prime})\,:\,\mathbbm{1}^{\top}w=\mathbbm{1}^{\top}w^{\prime}=|E|,\,\,\,w,w^{\prime}\geq 0\right\}.

where ww does not have to equal ww^{\prime}. If GG is conformally rigid then 𝟙=w=w\mathbbm{1}=w=w^{\prime} is an optimal solution to (6) and the optimal value of (6) is the smallest spectral width of GG. However, if 𝟙\mathbbm{1} is an optimal solution of (5), GG may not be conformally rigid since (6) allows for different weights in λ2(w)\lambda_{2}(w) and λn(w)\lambda_{n}(w^{\prime}) potentially making their difference smaller than λn(𝟙)λ2(𝟙)\lambda_{n}(\mathbbm{1})-\lambda_{2}(\mathbbm{1}).

1.3.2. The discrete Wirtinger inequality

Another topic that is related to, and partially inspired our results, is the discrete Wirtinger inequality proved by Fan, Taussky and Todd [8] in 1955 (see also an earlier paper of Schoenberg [26]).

Theorem 1.4 (Fan-Taussky-Todd [8]).

If x1,,xn+1x_{1},\dots,x_{n+1} are real numbers so that x1=xn+1x_{1}=x_{n+1} and i=1nxi=0\sum_{i=1}^{n}x_{i}=0, then

i=1n(xixi+1)24sin2(πn)i=1nxi2\sum_{i=1}^{n}(x_{i}-x_{i+1})^{2}\geq 4\sin^{2}\left(\frac{\pi}{n}\right)\sum_{i=1}^{n}x_{i}^{2}

with equality unless xi=Acos(2πi/n)+Bsin(2πi/n)x_{i}=A\cos(2\pi i/n)+B\sin(2\pi i/n).

In modern language, this inequality is a way of describing the second eigenvalue λ2\lambda_{2} of the graph Laplacian on the cycle graph CnC_{n} with nn vertices. It is clear that this argument, interpreting the quadratic form as that of the Laplacian of CnC_{n} and the constant as the spectral gap, will generalize to weights, but the arising matrix computations are not trivial [19, 21, 30]. We present a somewhat dual statement of the usual discrete Wirtinger inequality with weights.

Corollary 1.5.

The cycle graph CnC_{n} is conformally rigid. In particular, if w1,wnw_{1},\dots w_{n} are positive edge weights, then there always exist real numbers x1,,xn+1x_{1},\dots,x_{n+1} so that x1=xn+1x_{1}=x_{n+1} as well as i=1nxi=0\sum_{i=1}^{n}x_{i}=0 and

i=1nwi(xixi+1)24sin2(πn)(1ni=1nwi)i=1nxi2\sum_{i=1}^{n}w_{i}(x_{i}-x_{i+1})^{2}\leq 4\sin^{2}\left(\frac{\pi}{n}\right)\left(\frac{1}{n}\sum_{i=1}^{n}w_{i}\right)\sum_{i=1}^{n}x_{i}^{2}

One way of phrasing this Corollary is that the largest constant in the discrete Wirtinger inequality is attained if all the weights are the same: one cannot hope to get an improved discrete Wirtinger inequality by a clever choice of weights. Note that this is simply Proposition 1.3 applied to the cycle graph. One also obtains a dual result (from the second part of Proposition 1.3): for any choice of weights there exists a function that oscillates at least as rapidly as the most oscillating function for equal weights, however, this case seems to be less relevant in applications.

1.4. Organization of the paper

In § 2 we state our main results beginning with a summary. In § 3 we explain how semidefinite programming can be used to both test and certify conformal rigidity. We also explain how the automorphisms of the graph can sometimes drastically reduce the computational effort and also provide structural insight. § 4 discusses the connection between conformal rigidity and graph embeddings into Euclidean space. Section 5 contains the proofs of all our main results. Section 6 surveys various ways of proving that an explicitly given graph is conformally rigid via semidefinite programming certificates.

2. Main Results

2.1. Summary

Our main result is that the class of conformally rigid graphs form an interesting and highly nontrivial collection of graphs.

  1. (1)

    All edge-transitive graphs are conformally rigid and this is true for a simple reason (see §2.1). The smallest example of a conformally rigid graph that is not edge-transitive (that we know of) is the Hoffman graph on 16 vertices.

  2. (2)

    All distance-regular graphs and, as a special case, all strongly-regular graphs are conformally rigid.

  3. (3)

    Some Cayley graphs are conformally rigid and others are not. We provide a sufficient condition for Cayley graphs to be conformally rigid (Theorem 2.3). The smallest Cayley graph that we found that is conformally rigid but not edge-transitive is a circulant graph on 18 vertices (see §2.2).

  4. (4)

    As an application of the criterion for Cayley graphs, we found (empirically) that there seems to be a large number of circulant graphs that are conformally rigid (see §2.3). They seem to have a rich structure. There are also infinite families of circulants that are not conformally rigid.

  5. (5)

    Conformal rigidity can be characterized in terms of the existence of special spectral embeddings of the graph on the eigenspaces of λ2(𝟙)\lambda_{2}(\mathbbm{1}) and λn(𝟙)\lambda_{n}(\mathbbm{1}). These embeddings are edge-isometric, meaning that that all edge lengths must be equal (see § 4).

  6. (6)

    For any given graph G=(V,E)G=(V,E) it can be decided whether it is conformally rigid with a finite amount of computation by phrasing it as a semidefinite program. We explain how this is done and use it to verify conformal rigidity of a number of seemingly isolated examples that do not seem to fall into any of the above categories

    • the Hoffman graph on 16 vertices.

    • the crossing number graph 6B on 20 vertices [24].

    • the distance-2 graph of the Klein graph on 24 vertices.

    • the (20,8)(20,8) accordion graph [16] on 40 vertices.

We summarize our results in the figure below. A solid arrow implies containment while a dashed arrow implies partial containment.

distance-regulardistance-transitivestrongly regularconformally rigidarc-transitivevertex- and edge-transitiveedge-transitiveCayleyTheorem 2.3Proposition 2.1Theorem 2.2
Figure 3. Summary of our main results.

Throughout this paper we assume that the vertex set of G=(V,E)G=(V,E) is enumerated by the first nn integers, V=[n]V=[n], and that all graphs are connected, undirected and simple. Also, for simplicity, we sometimes write λ2\lambda_{2} for λ2(𝟙)\lambda_{2}(\mathbbm{1}) and λn\lambda_{n} for λn(𝟙)\lambda_{n}(\mathbbm{1}).

2.2. Edge-transitive graphs are conformally rigid

Recall that an automorphism of G=(V,E)G=(V,E) is a permutation π\pi of its vertices so that (i,j)E(i,j)\in E if and only if (π(i),π(j))E(\pi(i),\pi(j))\in E. The set of all automorphisms of GG, denoted Aut(G)\textup{Aut}(G) and called the automorphism group of GG, is a subgroup of SnS_{n}. The graph GG is vertex-transitive if all vertices of GG lie in a single orbit of Aut(G)\textup{Aut}(G). The action of Aut(G)\textup{Aut}(G) on VV induces an action on EE via (i,j)(π(i),π(j))(i,j)\mapsto(\pi(i),\pi(j)). The graph GG is edge-transitive if all edges of GG lie in a single orbit of Aut(G)\textup{Aut}(G). Complete graphs are edge-transitive and so are many other common families of graphs such as cycles, complete bipartite graphs and edge graphs of hypercubes. All distance-transitive graphs are edge-transitive (see Fig. 3). The following is a first observation.

Proposition 2.1.

Every edge-transitive graph is conformally rigid.

The proof of this Proposition is a simple corollary of the technique of symmetry reduction (see for example [3]) which we recall in § 3. All proofs of results in this section can be found in § 5. One way of interpreting Proposition 2.1 is as follows: edge-transitivity is a strong symmetry property that essentially says that all edges in GG are the “same” which points toward weighting them equally. Therefore, in some sense, the most interesting conformally rigid graphs are the ones that are not edge-transitive since they are conformally rigid for a different reason. In the rest of the paper, we will focus on such graphs.

2.3. Distance-regular graphs

Recall that a finite graph G=(V,E)G=(V,E) is distance-regular if for any two vertices u,vu,v, the number of vertices that is distance ii from uu and distance jj from vv only depends on i,ji,j and the distance d(u,v)d(u,v).

Theorem 2.2.

Every distance-regular graph is conformally rigid.

Note that, in particular, every strongly-regular graph is distance-regular. Strongly-regular graphs are distance-regular graphs of diameter two (where diameter is the largest distance between any pair of points). For example, the Shrikhande graph and its complement (see Fig. 1) are strongly-regular. The proof of Theorem 2.2 uses structural theory of distance-regular graphs (see [11]) and the duality theory of semidefinite programs which are developed in § 3 and § 4.

Refer to captionRefer to captionRefer to caption
Figure 4. Three conformally rigid graphs that are distance-regular but neither strongly-regular nor edge-transitive. Left: a bipartite (0,2)(0,2)-graph. Middle: a generalized quadrangle graph. Right: Doob graph D(1,1)D(1,1).

2.4. Cayley Graphs

A natural class of graphs with a great deal of symmetry are Cayley graphs. The Cayley graph Γ(G,S)\Gamma(G,S) has vertex set V=GV=G where GG a finite group and the edge set E={(g,sg):gG,sS}E=\left\{(g,s\circ g):g\in G,s\in S\right\} where SGS\subset G is a generating set. Since we only deal with undirected graphs, we want the edge set to be symmetric which is achieved by requiring a symmetry in the generating set SS, meaning sSs1Ss\in S\implies s^{-1}\in S. We always assume this.

Cayley graphs need not be conformally rigid. The smallest example we found is shown on the left in Fig. 5. This unweighted graph has λ2=2\lambda_{2}=2 and λ6=5\lambda_{6}=5. As it turns out, both of these eigenvalues can be improved by changing the edge weights: by increasing the weights on the three edges connecting the two triangles (while decreasing the weight on the other edges to maintain the normalization), we can increase λ2\lambda_{2} and attain λ2(w)=18/7>2\lambda_{2}(w)=18/7>2. Likewise, by setting the three edge weights connecting the two triangles to 0 (and increasing the remaining weights uniformly) we can decrease the largest eigenvalue to λ6(w)=4.5<5\lambda_{6}(w)=4.5<5 (and thus, in particular, can achieve values smaller than 5 for small positive weights by continuity of the eigenvalues). On the other hand, the circulant graph shown on the right of Fig. 5 is conformally rigid but not edge-transitive. Its conformal rigidity is implied by the following Theorem.

Refer to captionRefer to caption
Figure 5. Left: The triangular prism graph is Cayley and not conformally rigid. Right: a Cayley graph on 18\mathbb{Z}_{18} (generated by S={5,1,1,5}S=\left\{-5,-1,1,5\right\}) that is conformally rigid and not edge-transitive.
Theorem 2.3.

Let Γ(G,S)\Gamma(G,S) be a finite Cayley graph. If there exist eigenvectors φ2,φn\varphi_{2},\varphi_{n} corresponding to λ2,λn\lambda_{2},\lambda_{n} such that for φ{φ2,φn}\varphi\in\left\{\varphi_{2},\varphi_{n}\right\} and sSs\in S

(7) gGφ(g)φ(gs)is independent ofs,\displaystyle\sum_{g\in G}\varphi(g)\varphi(g\circ s)\qquad\mbox{is independent of}~{}s,

then Γ(G,S)\Gamma(G,S) is conformally rigid.

Comments. Several comments are in order.

  1. (1)

    The condition in (7) is simple to check for any given eigenvector φ\varphi: one has to compute |S||S| different sums and check whether the values all coincide. This is illustrated below where we use Theorem 2.3 to prove that the circulant shown in Fig. 5 is conformally rigid.

  2. (2)

    Recall that SS is symmetric (sSsSs\in S\implies-s\in S). By changing variables one sees that gΓφ(g)φ(gs)=gΓφ(g)φ(g(s)),\sum_{g\in\Gamma}\varphi(g)\varphi(g\circ s)=\sum_{g\in\Gamma}\varphi(g)\varphi(g\circ(-s)), and hence it suffices to check |S|/2|S|/2 sums in (7).

  3. (3)

    The difficulty in applying Theorem 2.3 is that one has to find suitable eigenvectors. What we have observed in practice is that when the eigenspace corresponding to, say, λ2\lambda_{2}, has multiplicity larger than 1, then the property will typically only be true for a single element in the space.

  4. (4)

    The result is sufficient. It is not clear to us whether it is necessary. In small-scale numerical experiments (carried out exclusively in the setting of circulants) we found that every time a circulant turned out to be conformally rigid, there indeed seemed to be a suitable choice of eigenvectors φ2\varphi_{2} and φn\varphi_{n} which made the result applicable. It would be interesting to understand this better.

Example 2.4.

For the circulant graph on the right in Fig. 5, V=18={0,1,2,,17}V={\mathbb{Z}}_{18}=\{0,1,2,\ldots,17\} and S={1,1,5,5}S=\{1,-1,5,-5\} we choose

λ2=2,φ2=(1,1,0,1,1,0,1,1,0,1,1,0,1,1,0,1,1,0)\lambda_{2}=2,\,\,\,\varphi_{2}=(-1,-1,0,1,1,0,-1,-1,0,1,1,0,-1,-1,0,1,1,0)

and

λ18=8,φ18=(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1).\lambda_{18}=8,\,\,\,\varphi_{18}=(-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1,-1,1).

A simple computation then shows that for each s{1,1,5,5}s\in\{1,-1,5,-5\},

i18φ2(i)φ2(i+s)=6 and i18φ18(i)φ18(i+s)=18.\sum_{i\in{\mathbb{Z}}_{18}}\varphi_{2}(i)\varphi_{2}(i+s)=6\,\,\,\textup{ and }\sum_{i\in{\mathbb{Z}}_{18}}\varphi_{18}(i)\varphi_{18}(i+s)=-18.

Theorem 2.3 then implies the conformal rigidity of this circulant. ∎

2.5. Circulants

Denote a circulant with nn vertices and generators SS by Cn(S)C_{n}(S). It is typical to list only one of ss or s-s in SS, with the assumption that the other one is included. The set of vertices of Cn(S)C_{n}(S) is V=n={0,1,2,,n1}V={\mathbb{Z}}_{n}=\left\{0,1,2,\dots,n-1\right\} and its edges are (v,v+s)(v,v+s), as we vary over all vVv\in V and sSs\in S (mod nn). Since circulants are Cayley graphs, Theorem 2.3 applies to them. We start with some examples.

  1. (1)

    The cycle graph CnC_{n} and the complete graph KnK_{n} are circulants. These are edge-transitive and hence conformally rigid by Proposition 2.1.

  2. (2)

    Paley graphs of prime order are circulants, edge-transitive and thus also conformally rigid.

  3. (3)

    Cocktail Party graphs are circulants, distance-transitive and hence, edge-transitive, and therefore conformally rigid.

  4. (4)

    The smallest circulant on two generators which is not edge-transitive but nonetheless conformally rigid is C18({1,5})C_{18}(\left\{1,5\right\}) (see Fig. 5).

Circulants are rich and rather nontrivial. In particular, understanding when a circulant is edge-transitive appears to be a nontrivial problem, no definite characterization appears to be known. There are infinitely many circulants that are not conformally rigid.

Proposition 2.5.

For all n7n\geq 7, the circulant Cn({1,2})C_{n}(\left\{1,2\right\}) is not conformally rigid.

The restriction n7n\geq 7 is necessary, C5({1,2})=K5C_{5}(\left\{1,2\right\})=K_{5} is conformally rigid and so is C6({1,2})C_{6}(\left\{1,2\right\}) (also known as the octahedral graph). The list of examples in (1)-(4) contains several infinite families of conformally rigid circulants, however, all of them are edge-transitive. It is a natural question whether there are infinitely many conformally rigid circulants that are not edge-transitive. We were unable to answer this question. Numerically, they appear to be fairly frequent. As an example, the complete list of conformally rigid circulants on n=21n=21 vertices with 2 generators is listed below (with graphs that are not edge transitive underlined):

C21({1,6})¯,C21({1,8}),C21({2,5}),C21({2,9})¯,C21({3,4})¯,C21({3,10})¯\underline{C_{21}(\left\{1,6\right\})},C_{21}(\left\{1,8\right\}),C_{21}(\left\{2,5\right\}),\underline{C_{21}(\left\{2,9\right\})},\underline{C_{21}(\left\{3,4\right\})},\underline{C_{21}(\left\{3,10\right\})}
C21({4,10}),C21({5,9})¯,C21({6,8})¯,C21({10,11}).C_{21}(\left\{4,10\right\}),\underline{C_{21}(\left\{5,9\right\})},\underline{C_{21}(\left\{6,8\right\})},C_{21}(\left\{10,11\right\}).
Figure 6. Conformally rigid circulants on n=21n=21 vertices with two generators. Underlined circulants are not edge-transitive.

Questions. We conclude with a number of questions.

  1. (1)

    Is it possible to characterize conformally rigid circulants Cn({a,b})C_{n}(\left\{a,b\right\}) in terms of n,a,bn,a,b? While there seem to be a great many, we were unable to detect a simple pattern or even identify a single infinite family.

  2. (2)

    Somewhat unrelated, but of independent interest in this context: is it possible to characterize when Cn({a,b})C_{n}(\left\{a,b\right\}) is edge-transitive in terms of n,a,bn,a,b? Some partial results can be found in [2, 5, 6, 25, 28, 29].

  3. (3)

    Theorem 2.3 provides a rigorous criterion that can be used to prove conformal rigidity of Cayley graphs and, in particular, circulants. Does the Theorem identify all conformally rigid circulants? Numerically this seems to be true for small conformally rigid circulants.

  4. (4)

    Are there infinitely many conformally rigid circulants Cn({a,b})C_{n}(\left\{a,b\right\}) that are not edge-transitive? The smallest one appears to be C18({1,5})C_{18}(\left\{1,5\right\}) described in Example 2.4, but there appear to be many more after that.

  5. (5)

    All these questions would be also interesting in the context of circulants with three or more generators.

2.6. Sporadic outliers

At this point, we have identified a number of reasons why a graph may be conformally rigid; most of the examples can be identified as belonging to a number of different groups, these being

  1. (1)

    edge-transitive graphs,

  2. (2)

    distance-regular graphs, in particular, strongly-regular graphs,

  3. (3)

    and certain Cayley graphs, in particular, certain circulant graphs.

The purpose of this section is to describe some isolated sporadic examples of graphs that do not belong into any of these groups.

Theorem 2.6.

The following graphs are conformally rigid, not edge-transitive and not distance-regular:

  1. (1)

    the Hoffman graph on 16 vertices (see Fig. 1)

  2. (2)

    the crossing number graph 6B on 20 vertices [24] (see Fig. 1)

  3. (3)

    the Haar graph 565 on 20 vertices (see Fig. 11)

  4. (4)

    the distance-2 graph of the Klein graph on 24 vertices (see Fig. 7)

  5. (5)

    the (7,1)(7,1)-bipartite (0,2)(0,2)-graph on n=48n=48 vertices (see Fig. 7)

  6. (6)

    the (7,2)(7,2)-bipartite (0,2)(0,2)-graph on n=48n=48 vertices

  7. (7)

    the (20,8)(20,8) accordion graph [16] on 40 vertices (see Fig. 7)

  8. (8)

    the non-Cayley vertex-transitive graph (24,23) on n=24n=24 vertices.

Refer to captionRefer to captionRefer to caption
Figure 7. Sporadic outliers: conformally rigid graphs that are neither edge-transitive nor distance-regular. Left: distance-2 graph of the Klein graph. Middle: a bipartite (0,2)(0,2)-graph on n=48n=48 vertices. Right: the (20,8)(20,8) accordion graph [16].

The graphs in this list are mostly well-known; the less well known examples (3), (5), (6), (8) are implemented in Mathematica 13.3 as (”Haar”, 565), (”ZeroTwoBipartite”, 7, 1),(”ZeroTwoBipartite”, 7, 2) and (”NoncayleyTransitive”, 24, 23). The graphs in Theorem 2.6 were all proven to be conformally rigid by identification of a suitable dual semidefinite programming certificate (see § 6). However, we also found one other example, the non-Cayley vertex-transitive graph (20,10) on n=20n=20 vertices, that is not edge-transitive, not distance-regular, but very likely conformally rigid (beyond any reasonable numerical doubt); however, we were unable to find a certificate that would rigorously prove its conformal rigidity.

3. Semidefinite programming certificates for conformal rigidity

3.1. Introduction.

The conformal rigidity of a graph can be checked and certified using semidefinite programming which is a branch of convex optimization [4]. This connection provides both theoretical and computational tools as we will see in § 4 and § 6. We first explain the necessary background.

Assume that G=([n],E)G=([n],E) is an undirected connected graph with positive edge weights wijw_{ij}. Let 𝒮+n\mathcal{S}_{+}^{n} be the cone of n×nn\times n real symmetric psd matrices. We use X0X\succeq 0 to denote that X𝒮+nX\in\mathcal{S}^{n}_{+}. For A,B𝒮+nA,B\in\mathcal{S}^{n}_{+}, ABA\succeq B stands for AB0A-B\succeq 0.

Recall from § 1.3.1 that the problems of maximizing λ2(w)\lambda_{2}(w) and minimizing λn(w)\lambda_{n}(w) with the sum of edge weights being bounded, have been studied in a number of papers such as [9, 12, 13, 27]. We follow the exposition in [27] where one observes

(8) λ2:= max {λ2(w):wij=|E|,w0}.\displaystyle\lambda_{2}^{\ast}:=\textup{ max }\left\{\lambda_{2}(w)\,:\,\sum w_{ij}=|E|,\,\,w\geq 0\right\}.

is a convex optimization problem that is equivalent to

(9) p:= min{ijEwij:λ2(w)1,w0}.\displaystyle p^{\ast}:=\textup{ min}\left\{\sum_{ij\in E}w_{ij}\,:\,\lambda_{2}(w)\geq 1,\,\,w\geq 0\right\}.

Problem (8) has an optimal solution since its objective function is continuous and its feasible region is compact. Therefore, (9) also has an optimal solution. If ww^{*} is optimal for (9), then λ2(w)=1\lambda_{2}(w^{*})=1 and hence there is no harm in optimizing over λ2(w)1\lambda_{2}(w)\geq 1. The two problems are off by a scaling; λ2=|E|/p\lambda_{2}^{*}=|E|/p^{*} and ww^{*} is optimal for (9) if and only if (|E|w)/p(|E|w^{*})/p^{*} is optimal for (8). Using II to denote the identity matrix and JJ to denote the matrix filled with 1’s, the constraint λ2(w)1\lambda_{2}(w)\geq 1 can be modeled by LwI(1/n)JL_{w}\succeq I-(1/n)J (the eigenvalues of I(1/n)JI-(1/n)J are 1,1,,1,01,1,\ldots,1,0) and (9) is the semidefinite program (SDP):

(10) p:= minijEwijs.t. LwI1nJw0.\displaystyle\begin{split}p^{\ast}:=&\textup{ min}\sum_{ij\in E}w_{ij}\\ \textup{s.t. }&L_{w}\succeq I-\frac{1}{n}J\\ &w\geq 0.\end{split}

SDPs are optimization problems in the space of symmetric matrices where we are optimizing a linear function over a (convex) region described by linear matrix inequalities, or equivalently, an affine slice of the cone of psd matrices, see [4]. Calling (10) the primal SDP, its dual SDP is derived in [27] as

(11) d:= max Trace Xs.t. Xii+Xjj2Xij1ijE𝟙X𝟙=0,X0.\displaystyle\begin{split}d^{\ast}:=&\textup{ max }\textup{Trace }X\\ \textup{s.t. }&X_{ii}+X_{jj}-2X_{ij}\leq 1\,\,\forall ij\in E\\ &\mathbbm{1}^{\top}X\mathbbm{1}=0,X\succeq 0.\end{split}

Strong duality holds between the primal and dual SDPs and so p=dp^{*}=d^{*}. A pair of primal/dual feasible solutions ww^{*} and XX^{*} are optimal for their respective problems if and only if they satisfy the complementary slackness conditions

(12) (1(Xii+Xjj2Xij))wij=0ijE\displaystyle(1-(X_{ii}^{*}+X_{jj}^{*}-2X_{ij}^{*}))w_{ij}^{*}=0\,\,\forall ij\in E
(13) LwX=X\displaystyle L_{w^{*}}X^{*}=X^{*}

Analogous to the above, the problem of minimizing λn(w)\lambda_{n}(w) can be written as

(14) λn:= min {λn(w):wij=|E|,w0}.\displaystyle\lambda_{n}^{\ast}:=\textup{ min }\left\{\lambda_{n}(w)\,:\,\sum w_{ij}=|E|,\,\,w\geq 0\right\}.

which is equivalent to

(15) q:= max{ijEwij:λn(w)1,w0}.\displaystyle q^{\ast}:=\textup{ max}\left\{\sum_{ij\in E}w_{ij}\,:\,\lambda_{n}(w)\leq 1,\,\,w\geq 0\right\}.

As before, if ww^{*} is optimal for (15) then λn(w)=1\lambda_{n}(w^{*})=1, λn=|E|/q\lambda_{n}^{*}=|E|/q^{*} and ww^{*} is optimal for (15) if and only if (|E|w)/q(|E|w^{*})/q^{*} is optimal for (14). Problem (15) is the SDP:

(16) q:= maxijEwijs.t. LwI1nJw0\displaystyle\begin{split}q^{\ast}:=&\textup{ max}\sum_{ij\in E}w_{ij}\\ &\textup{s.t. }L_{w}\preceq I-\frac{1}{n}J\\ &w\geq 0\end{split}

whose dual is

(17) r:= min Trace Ys.t. Yii+Yjj2Yij1ijE𝟙Y𝟙=0,Y0.\displaystyle\begin{split}r^{\ast}:=&\textup{ min }\textup{Trace }Y\\ &\textup{s.t. }Y_{ii}+Y_{jj}-2Y_{ij}\geq 1\,\,\forall ij\in E\\ &\mathbbm{1}^{\top}Y\mathbbm{1}=0,\,\,Y\succeq 0.\end{split}

Again, strong duality holds so that q=rq^{*}=r^{*}, and a pair of primal/dual feasible solutions (w,Y)(w^{*},Y^{*}) are optimal for their respective problems if and only if they satisfy complementary slackness:

(18) (Yii+Yjj2Yij1)wij=0ijE\displaystyle(Y^{*}_{ii}+Y^{*}_{jj}-2Y^{*}_{ij}-1)w^{*}_{ij}=0\,\,\forall ij\in E
(19) LwY=Y\displaystyle L_{w^{*}}Y^{*}=Y^{*}

We can now test for conformal rigidity via the above optimization problems.

Observation 3.1.

A graph GG is conformally rigid if and only if

  1. (1)

    𝟙/λ2(𝟙)\mathbbm{1}/\lambda_{2}(\mathbbm{1}) is an optimal solution to (9) (or (10)) with p=|E|/λ2(𝟙)p^{*}=|E|/\lambda_{2}(\mathbbm{1}), and

  2. (2)

    𝟙/λn(𝟙)\mathbbm{1}/\lambda_{n}(\mathbbm{1}) is an optimal solution to (15) (or (16)) with q=|E|/λn(𝟙)q^{*}=|E|/\lambda_{n}(\mathbbm{1}).

Under strong duality, the pair of primal and dual SDPs will have the same optimal value, which allows the following re-phrasing:

Observation 3.2.

A graph GG is conformally rigid if and only if

  1. (1)

    The SDP (11) has an optimal solution XX^{*} with Trace(X)=|E|/λ2(𝟙)\textup{Trace}(X^{*})=|E|/\lambda_{2}(\mathbbm{1}), and

  2. (2)

    the SDP (17) has an optimal solution YY^{*} with Trace (Y)=|E|/λn(𝟙)\textup{Trace }(Y^{*})=|E|/\lambda_{n}(\mathbbm{1}).

Usually, these observations cannot be used directly to certify conformal rigidity since SDPs are solved using numerical algorithms making it difficult to know precisely what their optimal values and solutions are. In § 6 we illustrate several methods to overcome these difficulties. A more computationally friendly certificate for conformal rigidity can be obtained from the complementary slackness conditions. Indeed, when GG is conformally rigid, the optimal weights are 𝟙/λ2(𝟙)\mathbbm{1}/\lambda_{2}(\mathbbm{1}) and 𝟙/λn(𝟙)\mathbbm{1}/\lambda_{n}(\mathbbm{1}) and so complementary slackness implies the following.

Proposition 3.3.

A graph GG is conformally rigid if and only if there is a feasible solution XX for (11) and a feasible solution YY for (17) such that

(20) Xii+Xjj2Xij=1ijE,LX=λ2(𝟙)X\displaystyle X_{ii}+X_{jj}-2X_{ij}=1\quad\,\,\forall ij\in E,\,\,\,LX=\lambda_{2}(\mathbbm{1})X
(21) Yii+Yjj2Yij=1ijE,LY=λn(𝟙)Y.\displaystyle Y_{ii}+Y_{jj}-2Y_{ij}=1\quad\,\,\forall ij\in E,\,\,\,LY=\lambda_{n}(\mathbbm{1})Y.

The condition LX=λ2(𝟙)XLX=\lambda_{2}(\mathbbm{1})X is saying that the columns of XX are in the eigenspace of LL for λ2(𝟙)\lambda_{2}(\mathbbm{1}). Therefore, the rank of XX is bounded above by the multiplicity of λ2(𝟙)\lambda_{2}(\mathbbm{1}) as an eigenvalue of LL. Analogously, the columns of YY must come from the eigenspace of LL for λn(𝟙)\lambda_{n}(\mathbbm{1}) and so the rank of YY cannot exceed the multiplicity of λn(𝟙)\lambda_{n}(\mathbbm{1}). In particular, if λ2(𝟙)\lambda_{2}(\mathbbm{1}) (respectively, λn(𝟙)\lambda_{n}(\mathbbm{1})) has multiplicity 11, then XX (respectively, YY) can be taken to be the rank one matrix

(22) puu( respectively, quu)\displaystyle p^{*}uu^{\top}\textup{}(\textup{ respectively, }q^{*}uu^{\top})

where uu is a normalized eigenvector of λ2(𝟙)\lambda_{2}(\mathbbm{1}) (respectively, λn(𝟙)\lambda_{n}(\mathbbm{1})). These observations were noted in [27]. A very short explanation why the notion of conformal rigidity simplifies when the eigenspace associated to λ2\lambda_{2} (or λn\lambda_{n}) has multiplicity 1 is as follows: in that case, we have

λ2(𝟙)=minvVf(v)=0(i,j)E(f(i)f(j))2vVf(v)2\lambda_{2}(\mathbbm{1})=\min_{\sum_{v\in V}f(v)=0}\frac{\sum_{(i,j)\in E}(f(i)-f(j))^{2}}{\sum_{v\in V}f(v)^{2}}

We note that introducing weights wijw_{ij} (normalized to (i,j)Ewij=|E|\sum_{(i,j)\in E}w_{ij}=|E|) and then continuously varying the weights leads to a continuously perturbation of the matrix. If λ2\lambda_{2} has multiplicity 1, then both the eigenvalue λ2(w)\lambda_{2}(w) and the associated eigenvector depend continuously on the change in the weights (and, in particular, do not change much for small changes in the weight). Suppose now that, for the eigenvector in the unweighted case, not all terms in the sum (f(i)f(j))2(f(i)-f(j))^{2} are of the same size: then we could consider a small perturbation of the weights, putting more emphasis on larger summands and less on smaller ones and arrive at a weight for which λ2(w)>λ2(𝟙)\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}) which contradicts conformal rigidity.

We now do a small example to illustrate the above ideas, one in which all claims can be checked by hand.

Example 3.4.

Let GG be the 44-cycle with edges 12,23,34,1412,23,34,14 which is edge-transitive and hence conformally rigid by Proposition 2.1. The Laplacian of GG has eigenvalues: 0,2,2,40,2,2,4 with eigenvectors:

(1,1,1,1),(1,0,1,0),(0,1,0,1)λ=2,(1,1,1,1)λ=4.(1,1,1,1),\underbrace{(-1,0,1,0),(0,-1,0,1)}_{\lambda=2},\underbrace{(-1,1,-1,1)}_{\lambda=4}.

The weights

w=(1/2,1/2,1/2,1/2), and w=(1/4,1/4,1/4,1/4)w^{*}=(1/2,1/2,1/2,1/2),\,\,\,\textup{ and }\,\,\,w^{\prime*}=(1/4,1/4,1/4,1/4)

are optimal solutions to (9) and (15). To verify these claims, we construct certificates X,YX,Y as in Proposition 3.3. Take

X=12(1010010110100101),Y=14(1111111111111111).X=\frac{1}{2}\begin{pmatrix}1&0&-1&0\\ 0&1&0&-1\\ -1&0&1&0\\ 0&-1&0&1\end{pmatrix},\,\,\,\,Y=\frac{1}{4}\begin{pmatrix}1&-1&1&-1\\ -1&1&-1&1\\ 1&-1&1&-1\\ -1&1&-1&1\end{pmatrix}.

Check that XX is a feasible solution of (11) that satisfies (20); the eigenvalues of XX are 0,0,1,10,0,1,1 which shows that X0X\succeq 0. The analogous checks can be made for YY, and (X,Y)(X,Y) is a pair of SDP certificates for the conformal rigidity of GG.

Note that Trace X=2=𝟙w\textup{Trace }X=2=\mathbbm{1}^{\top}w^{*} and Trace Y=1=𝟙w\textup{Trace }Y=1=\mathbbm{1}^{\top}w^{\prime*}. Since the multiplicity of λ2\lambda_{2} is 22, rank(X)2\operatorname{rank}(X)\leq 2. In fact, rank(X)=2\operatorname{rank}(X)=2. Since the multiplicity of λ4=4\lambda_{4}=4 is 11, we expect YY to be a scaling of uuuu^{\top} where u=(1,1,1,1)u=(-1,1,-1,1) is the eigenvector of 44 shown above, and this is indeed the case. ∎

3.2. Symmetry reduction

The symmetries of GG allow the SDPs (10) and (16) to be reduced to have less variables, allowing them to be more easily solvable and helping to understand their structure. We briefly explain this well-known technique in our setting and use it to prove the conformal rigidity of some graphs where direct methods were not successful. This exposition closely follows that in [3].

Recall that a permutation πAut(G)\pi\in\textup{Aut}(G) acts on vertices, edges and weights of GG via iπ(i),ijπ(i)π(j),wijwπ(i)π(j)i\mapsto\pi(i),\,\,ij\mapsto\pi(i)\pi(j),\,\,w_{ij}\mapsto w_{\pi(i)\pi(j)}. Let PP denote the n×nn\times n permutation matrix corresponding to π\pi with Pij=1P_{ij}=1 if π(j)=i\pi(j)=i and 0 otherwise. If we denote the weight obtained by the action of π\pi on ww by πw\pi\cdot w, then Lπw=PLwPL_{\pi\cdot w}=PL_{w}P^{\top}.

Now consider the SDP (10) and its feasible region

(23) :={wE:LwI1nJ,w0}.\displaystyle\mathcal{F}:=\{w\in{\mathbb{R}}^{E}\,:\,L_{w}\succeq I-\frac{1}{n}J,\,\,w\geq 0\}.

If ww\in\mathcal{F} and πAut(G)\pi\in\textup{Aut}(G), then πw0\pi\cdot w\geq 0 and

Lπw=PLwPP(I1nJ)P=I1nJ.L_{\pi\cdot w}=PL_{w}P^{\top}\succeq P(I-\frac{1}{n}J)P^{\top}=I-\frac{1}{n}J.

Therefore, \mathcal{F} is invariant under the action of Aut(G)\textup{Aut}(G).

By the invariance of \mathcal{F}, if ww^{*} is an optimal solution of (10), then πw\pi\cdot w^{*}\in\mathcal{F} for all πAut(G)\pi\in\textup{Aut}(G) and therefore, by the convexity of \mathcal{F}, we also have that

(24) w¯:=1|Aut(G)|ππw.\displaystyle\overline{w^{\ast}}:=\frac{1}{|\textup{Aut}(G)|}\sum_{\pi}\pi\cdot w^{\ast}\in\mathcal{F}.

Further, w¯\overline{w^{\ast}} is also an optimal solution to (10) since 𝟙(πw)=𝟙w=p\mathbbm{1}^{\top}(\pi\cdot w^{*})=\mathbbm{1}^{\top}w^{*}=p^{*}. This optimal solution w¯\overline{w^{\ast}} is fixed under the action of Aut(G)\textup{Aut}(G). Thus we may restrict the optimization in (10) to the fixed point subset, ¯\overline{\mathcal{F}}\subset\mathcal{F}.

Let O(i),i=1,,kO(i),i=1,\ldots,k be the orbits of edges in EE under the action of Aut(G)\textup{Aut}(G). Any fixed point w¯w\in\overline{\mathcal{F}} must have the same value in all coordinates lying in an edge orbit O(i)O(i). Let wiw_{i} be the (variable) weight on the edges in O(i)O(i). Then (10) is equivalent to the symmetry reduced SDP:

(25) p:= mini=1k|O(i)|wis.t. i=1kwiLiI1nJwi0i\displaystyle\begin{split}p^{\ast}:=&\textup{ min}\sum_{i=1}^{k}|O(i)|w_{i}\\ \textup{s.t. }&\sum_{i=1}^{k}w_{i}L_{i}\succeq I-\frac{1}{n}J\\ &w_{i}\geq 0\,\,\forall i\end{split}

where LiL_{i} is the Laplacian of the graph Gi=([n],O(i))G_{i}=([n],O(i)). The symmetry reduction for the SDP (16) is analogous and yields the SDP:

(26) q:= maxi=1k|O(i)|wis.t. i=1kwiLiI1nJwi0i\displaystyle\begin{split}q^{\ast}:=&\textup{ max}\sum_{i=1}^{k}|O(i)|w_{i}\\ \textup{s.t. }&\sum_{i=1}^{k}w_{i}L_{i}\preceq I-\frac{1}{n}J\\ &w_{i}\geq 0\,\,\forall i\end{split}

The dual SDPs (11) and (17) can also be symmetry reduced in a similar way.

4. Certificates and Spectral Embeddings

It was noted in [12] and [27] that the dual certificates XX and YY in Proposition 3.3 provide embeddings of GG. In the context of conformal rigidity, this interpretation can be developed further, and it will be our main tool to prove that all distance-regular graphs are conformally rigid.

Definition 4.1.

Let G=([n],E)G=([n],E) be a connected graph and let λ>0\lambda>0 be an eigenvalue of its Laplacian LL, of multiplicity mm. Let λm\mathcal{E}_{\lambda}\cong{\mathbb{R}}^{m} be the eigenspace of λ\lambda and Pn×kP\in{\mathbb{R}}^{n\times k} a matrix whose columns lie in λ\mathcal{E}_{\lambda}. The collection of vectors 𝒫={pi,i=1,,n}k\mathcal{P}=\{p_{i},\,\,i=1,\ldots,n\}\subset{\mathbb{R}}^{k}, such that p1,pnp_{1}^{\top},\ldots p_{n}^{\top} are the rows of PP, is called an embedding of GG on λ\mathcal{E}_{\lambda}.

Indeed, the configuration 𝒫\mathcal{P} provides a realization/embedding of GG in k{\mathbb{R}}^{k} by assigning the vector pip_{i} to vertex ii, and connecting pip_{i} and pjp_{j} by an edge for all ijEij\in E. The embedding is centered in the sense that pi=0\sum p_{i}=0 since the columns of PP are orthogonal to 𝟙\mathbbm{1}. Spectral embeddings as above were defined by Hall [17] and several variants have been studied. Typically, there are more requirements on the matrix PP such as the columns of PP should form an orthonormal basis of λ\mathcal{E}_{\lambda}. When PP is chosen to have this extra property, we will write UU in the place of PP and 𝒰\mathcal{U} in the place of 𝒫\mathcal{P}. The embedding is spherical if pi=α>0\|p_{i}\|=\alpha>0 for all i=1,,ni=1,\ldots,n. We will see embeddings with these additional properties for distance-regular graphs.

Definition 4.2.

An embedding 𝒫\mathcal{P} of GG on λ\mathcal{E}_{\lambda} is said to be edge-isometric if there exists a constant c>0c>0 such that pipj=c\|p_{i}-p_{j}\|=c for all ijEij\in E.

The following is an adaptation of an observation in [12] and [27] in the context of conformal rigidity.

Proposition 4.3.

A connected graph G=([n],E)G=([n],E) is conformally rigid if and only if GG has an edge-isometric embedding 𝒫\mathcal{P} on λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})} and an edge-isometric embedding 𝒬\mathcal{Q} on λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})}.

Proof.

Suppose GG is conformally rigid. Then there exists matrices XX and YY satisfying the conditions of Proposition 3.3. Since XX and YY are psd, they admit a Gram decomposition X=PPX=PP^{\top} and Y=QQY=QQ^{\top}. The columns of PP lie in λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})} since the columns of XX do. Similarly, the columns of QQ lie in λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})} as the columns of YY do. Therefore, 𝒫\mathcal{P} and 𝒬\mathcal{Q}, consisting of the rows of PP and QQ respectively, are embeddings of GG on λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})} and λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})} respectively. The condition

1=Xii+Xjj2Xij=pipj2ijE1=X_{ii}+X_{jj}-2X_{ij}=\|p_{i}-p_{j}\|^{2}\,\,\,\forall ij\in E

ensures that 𝒫\mathcal{P} is edge-isometric. Similarly for 𝒬\mathcal{Q}.

Conversely, suppose we have embeddings as in the statement of the Proposition. Let 𝒫\mathcal{P} be the edge-isometric embedding on λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})}. Suppose pipj=c>0\|p_{i}-p_{j}\|=c>0 for all ijEij\in E and the columns of PP lie in λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})}. Define X:=c2PP𝒮+nX:=c^{-2}PP^{\top}\in\mathcal{S}^{n}_{+}. Then 𝟙X𝟙=0\mathbbm{1}^{\top}X\mathbbm{1}=0 since the columns of PP are eigenvectors of λ2(𝟙)\lambda_{2}(\mathbbm{1}) and are hence orthogonal to 𝟙\mathbbm{1}. Since XX and PP have the same column space, LX=λ2(𝟙)XLX=\lambda_{2}(\mathbbm{1})X. If ijEij\in E, then

Xii+Xjj2Xij=1c2(pi2+pj22pipj)=1c2(pipj2)=1.X_{ii}+X_{jj}-2X_{ij}=\frac{1}{c^{2}}(\|p_{i}\|^{2}+\|p_{j}\|^{2}-2p_{i}^{\top}p_{j})=\frac{1}{c^{2}}(\|p_{i}-p_{j}\|^{2})=1.

Therefore, XX satisfies the conditions in Proposition 3.3. By the same argument we get a certificate YY needed in Proposition 3.3 from the edge-isometric embedding of GG on λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})}. Together, XX and YY certify that GG is conformally rigid. ∎

Remark 4.4.

The certificate X=c2PPX=c^{-2}PP^{\top} in the proof of Proposition 4.3 is an optimal solution to (11). Therefore, Trace X=c2Trace PP=|E|/λ2(𝟙)\textup{Trace }X=c^{-2}\textup{Trace }PP^{\top}=|E|/\lambda_{2}(\mathbbm{1}) which implies that c2=(λ2(𝟙)pi2)/|E|c^{2}=(\lambda_{2}(\mathbbm{1})\sum\|p_{i}\|^{2})/|E|.

Similarly, if 𝒬={q1,,qn}\mathcal{Q}=\{q_{1},\ldots,q_{n}\} is the edge-isometric embedding on λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})} that gave rise to the certificate YY, then Trace Y=(c)2qi2=|E|/λn(𝟙)\textup{Trace }Y=(c^{\prime})^{-2}\sum\|q_{i}\|^{2}=|E|/\lambda_{n}(\mathbbm{1}) which implies that c2=(λn(𝟙)qi2)/|E|c^{\prime 2}=(\lambda_{n}(\mathbbm{1})\sum\|q_{i}\|^{2})/|E|.

Next, we define a more ambitious type of embedding towards a characterization of when GG and its complement are both conformally rigid.

Definition 4.5.

An embedding 𝒫\mathcal{P} of G=([n],E)G=([n],E) on λ\mathcal{E}_{\lambda} as in Definition 4.2 is edge-nonedge-isometric if there are two nonzero constants α,β>0\alpha,\beta>0 such that

  1. (1)

    pipj=αijE\|p_{i}-p_{j}\|=\alpha\,\,\,\forall ij\in E, and

  2. (2)

    pipj=βijE\|p_{i}-p_{j}\|=\beta\,\,\,\forall ij\not\in E.

Conformal rigidity of GG does not imply that GG has an embedding on λ2\mathcal{E}_{\lambda_{2}} or λn\mathcal{E}_{\lambda_{n}} that is edge-nonedge-isometric.

Example 4.6.

The Hoffman graph on 1616 vertices in Fig 1 is conformally rigid, but its complement is not. Its Laplacian eigenvalues, with multiplicities, are:

(λ,mult(λ)):(8,1),(6,4),(4,6),(2,4),(0,1).(\lambda,\textup{mult}(\lambda)):(8,1),(6,4),(4,6),(2,4),(0,1).

The eigenspace 8\mathcal{E}_{8} is spanned by the vector

v=(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)v^{\top}=(-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1)

and therefore, consider the ordered configuration

𝒬={1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}.\mathcal{Q}=\left\{-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1\right\}.

Once can check that qiqj=2\|q_{i}-q_{j}\|=2 for all ijEij\in E, while if ijEij\not\in E, then qiqj=0,2\|q_{i}-q_{j}\|=0,2. For instance, (1,9),(1,10)E(1,9),(1,10)\in E and (1,2),(1,16)E(1,2),(1,16)\not\in E. Therefore, there are two different distances between non-adjacent vertices of GG in the embedding provided by 𝒬\mathcal{Q}. Since the eigenspace 8\mathcal{E}_{8} is 11-dimensional, any other embedding on this eigenspace must be a scaling of the configuration 𝒬\mathcal{Q}. The property of there being two different distances between non-adjacent vertices is invariant under scaling. ∎

In the rest of this section, assume that G=([n],E)G=([n],E) is a regular graph of valency dd. Then its complement Gc=([n],Ec)G^{c}=([n],E^{c}) is a (nd1)(n-d-1)-regular graph on [n][n]. Since the adjacency matrices are related as AGc=JAGIA_{G^{c}}=J-A_{G}-I where JJ is the matrix of all ones, we have that LGc=(nd)IJ+AGL_{G^{c}}=(n-d)I-J+A_{G}. If GG and GcG^{c} are both connected, then 0 is an eigenvalue of both LGL_{G} and LGcL_{G^{c}}, and λ>0\lambda>0 is an eigenvalue of LGL_{G} if and only if nλn-\lambda is an eigenvalue of LGcL_{G^{c}}, and both eigenvalues have the same eigenspace which we will write as λ=nλc\mathcal{E}_{\lambda}=\mathcal{E}^{c}_{n-\lambda}.

Conformal rigidity is not closed under graph complementation. A simple example is the 66-cycle C6C_{6} which is conformally rigid as it is edge-transitive, but its complement is not, see Example 6.4. Proposition 3.3 provides the following certification for the conformal rigidity of both GG and GcG^{c}.

Corollary 4.7.

Let G=([n],E)G=([n],E) be a connected regular graph whose complement GcG^{c} is also connected. Then GG and GcG^{c} are both conformally rigid if and only if there exists X,Y,X,Y𝒮+nX,Y,X^{\prime},Y^{\prime}\in\mathcal{S}^{n}_{+} such that

(27) Xii+Xjj2Xij=1ijE,cols(X)λ2(𝟙),   1X𝟙=0\displaystyle X_{ii}+X_{jj}-2X_{ij}=1\forall ij\in E,\,\,\,\textup{cols}(X)\subset\mathcal{E}_{\lambda_{2}(\mathbbm{1})},\,\,\,\mathbbm{1}^{\top}X\mathbbm{1}=0
(28) Yii+Yjj2Yij=1ijE,cols(Y)λn(𝟙)   1Y𝟙=0\displaystyle Y_{ii}+Y_{jj}-2Y_{ij}=1\forall ij\in E,\,\,\,\textup{cols}(Y)\subset\mathcal{E}_{\lambda_{n}(\mathbbm{1})}\,\,\,\mathbbm{1}^{\top}Y\mathbbm{1}=0
(29) Xii+Xjj2Xij=1ijE,cols(X)λn(𝟙),   1X𝟙=0\displaystyle X^{\prime}_{ii}+X^{\prime}_{jj}-2X^{\prime}_{ij}=1\forall ij\not\in E,\,\,\,\textup{cols}(X^{\prime})\subset\mathcal{E}_{\lambda_{n}(\mathbbm{1})},\,\,\,\mathbbm{1}^{\top}X^{\prime}\mathbbm{1}=0
(30) Yii+Yjj2Yij=1ijE,cols(Y)λ2(𝟙),   1Y𝟙=0.\displaystyle Y^{\prime}_{ii}+Y^{\prime}_{jj}-2Y^{\prime}_{ij}=1\forall ij\not\in E,\,\,\,\textup{cols}(Y^{\prime})\subset\mathcal{E}_{\lambda_{2}(\mathbbm{1})},\,\,\,\mathbbm{1}^{\top}Y^{\prime}\mathbbm{1}=0.

It would be particularly convenient if we could take X=YX^{\prime}=Y and Y=XY^{\prime}=X. Here is an embedding version of this wish, that follows from Proposition 4.3.

Corollary 4.8.

Let GG be a connected regular graph such that GcG^{c} is also connected. If GG has an embedding 𝒫\mathcal{P} on λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})} and an embedding 𝒬\mathcal{Q} on λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})} that are both edge-nonedge-isometric then both GG and GcG^{c} are conformally rigid.

A class of graphs that have edge-nonedge-isometric embeddings on λ2(𝟙)\mathcal{E}_{\lambda_{2}(\mathbbm{1})} and λn(𝟙)\mathcal{E}_{\lambda_{n}(\mathbbm{1})} are the connected strongly-regular graphs. If GG is a connected strongly-regular graph then there are only two distances among pairs of vertices in GG; d(i,j)=1d(i,j)=1 if ijEij\in E and d(i,j)=2d(i,j)=2 if ijEij\not\in E. Also GG has only two nonzero eigenvalues which are typically called rr and ss. Complements of strongly regular graphs are also strongly-regular and suppose r=nsr^{\prime}=n-s and s=nrs^{\prime}=n-r are the nonzero eigenvalues of GcG^{c}. We will see in Corollary 5.4 that GG has a spherical edge-isometric embedding 𝒰={u1,,un}\mathcal{U}=\{u_{1},\ldots,u_{n}\} on r\mathcal{E}_{r} where uiu_{i}^{\top} form the rows of a matrix UU whose columns form an orthonormal basis of r\mathcal{E}_{r}. Since r=sc\mathcal{E}_{r}=\mathcal{E}^{c}_{s^{\prime}}, 𝒰\mathcal{U} is also a spherical embedding of GcG^{c} on sc\mathcal{E}^{c}_{s^{\prime}} and by Corollary 5.4, it is an edge-isometric embedding of GcG^{c}. Since the edges of GcG^{c} are the nonedges of GG and vice-versa, 𝒰\mathcal{U} is an edge-nonedge-isometric embedding of GG on r\mathcal{E}_{r} (and of GcG^{c} on sc\mathcal{E}^{c}_{s^{\prime}}). Similarly, there is an edge-nonedge-isometric spherical embedding 𝒱\mathcal{V} of GG on s=rc\mathcal{E}_{s}=\mathcal{E}^{c}_{r^{\prime}} (and GcG^{c} on rc\mathcal{E}^{c}_{r^{\prime}}). The embeddings 𝒰\mathcal{U} and 𝒱\mathcal{V} are spherical 22-designs [11, Corollary 6.2]. Lemma 6.3 in [11] says that if 𝒰\mathcal{U} is a spherical 22-design such that uiuj{α,β}u_{i}^{\top}u_{j}\in\{\alpha,\beta\} whenever uiuju_{i}\neq u_{j}, and GG is the graph with vertex set 𝒰\mathcal{U} and two vectors adjacent if their inner product is α\alpha, then GG is strongly-regular.

Example 4.9.

Let GG be the Petersen graph which is strongly-regular. It is easy to see that GG and GcG^{c} are both conformally rigid since they are both edge-transitive. The eigenvalues of GG and GcG^{c} with their multiplicities are:

G:(5,4),(2,5),(0,1) and Gc:(8,5),(5,4),(0,1).G:(5,4),(2,5),(0,1)\,\,\,\textup{ and }\,\,\,G^{c}:(8,5),(5,4),(0,1).

Take 𝒰\mathcal{U} to be the columns of

U=(121200121200001231231301231231300012312312313001231300161616131313161323013213213213213213213213213212).{\small U^{\top}=\left(\begin{array}[]{cccccccccc}\frac{1}{2}&-\frac{1}{2}&0&0&-\frac{1}{2}&\frac{1}{2}&0&0&0&0\\ -\frac{1}{2\sqrt{3}}&\frac{1}{2\sqrt{3}}&-\frac{1}{\sqrt{3}}&0&-\frac{1}{2\sqrt{3}}&\frac{1}{2\sqrt{3}}&\frac{1}{\sqrt{3}}&0&0&0\\ -\frac{1}{2\sqrt{3}}&-\frac{1}{2\sqrt{3}}&\frac{1}{2\sqrt{3}}&-\frac{1}{\sqrt{3}}&0&0&\frac{1}{2\sqrt{3}}&\frac{1}{\sqrt{3}}&0&0\\ -\frac{1}{6}&-\frac{1}{6}&-\frac{1}{6}&\frac{1}{3}&-\frac{1}{3}&-\frac{1}{3}&-\frac{1}{6}&\frac{1}{3}&\frac{2}{3}&0\\ -\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&\frac{1}{3\sqrt{2}}&\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&\frac{1}{3\sqrt{2}}&\frac{1}{\sqrt{2}}\\ \end{array}\right)}.

Then 𝒰\mathcal{U} is an edge-nonedge-isometric embedding of GG on 2\mathcal{E}_{2} and of GcG^{c} on 8c=2\mathcal{E}^{c}_{8}=\mathcal{E}_{2}. Indeed, uiuj=2/3\|u_{i}-u_{j}\|=\sqrt{2/3} for all ijEij\in E and uiuj=2/3\|u_{i}-u_{j}\|=2/\sqrt{3} for all ijEij\not\in E.

Take 𝒱\mathcal{V} to be the columns of

V=(16160161616000161601616160001616132132231321321320232313213102325131013101310232525232513101310).{\small V^{\top}=\left(\begin{array}[]{cccccccccc}\frac{1}{\sqrt{6}}&\frac{1}{\sqrt{6}}&0&-\frac{1}{\sqrt{6}}&-\frac{1}{\sqrt{6}}&-\frac{1}{\sqrt{6}}&0&0&0&\frac{1}{\sqrt{6}}\\ \frac{1}{\sqrt{6}}&0&-\frac{1}{\sqrt{6}}&-\frac{1}{\sqrt{6}}&\frac{1}{\sqrt{6}}&0&0&0&\frac{1}{\sqrt{6}}&-\frac{1}{\sqrt{6}}\\ \frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&-\frac{\sqrt{2}}{3}&\frac{1}{3\sqrt{2}}&\frac{1}{3\sqrt{2}}&-\frac{1}{3\sqrt{2}}&0&\frac{\sqrt{2}}{3}&-\frac{\sqrt{2}}{3}&\frac{1}{3\sqrt{2}}\\ \frac{1}{3\sqrt{10}}&-\frac{2}{3}\sqrt{\frac{2}{5}}&\frac{1}{3\sqrt{10}}&\frac{1}{3\sqrt{10}}&\frac{1}{3\sqrt{10}}&-\frac{2}{3}\sqrt{\frac{2}{5}}&\sqrt{\frac{2}{5}}&-\frac{2}{3}\sqrt{\frac{2}{5}}&\frac{1}{3\sqrt{10}}&\frac{1}{3\sqrt{10}}\\ \end{array}\right)}.

Then 𝒱\mathcal{V} is an edge-nonedge-isometric embedding of GG on 5\mathcal{E}_{5} and of GcG^{c} on 5c=5\mathcal{E}^{c}_{5}=\mathcal{E}_{5}. In this case, vivj=2/3\|v_{i}-v_{j}\|=2/\sqrt{3} for all ijEij\in E and vivj=2/3\|v_{i}-v_{j}\|=\sqrt{2/3} for all ijEij\not\in E. The certificates are

X=Y=32UU and Y=X=34VV.X=Y^{\prime}=\frac{3}{2}UU^{\top}\,\,\textup{ and }\,\,Y=X^{\prime}=\frac{3}{4}VV^{\top}.

The scalings 3/23/2 and 3/43/4 are determined from the edge lengths in the two embeddings as in Remark 4.4. The columns of UU (resp. VV) form an orthonormal basis of 2=8c\mathcal{E}_{2}=\mathcal{E}^{c}_{8} (resp. 5=5c\mathcal{E}_{5}=\mathcal{E}^{c}_{5}). Since ui=1/2\|u_{i}\|=1/\sqrt{2} and vi=2/5\|v_{i}\|=\sqrt{2/5} for all i=1,,10i=1,\ldots,10, the embeddings 𝒰\mathcal{U} and 𝒱\mathcal{V} are both spherical. ∎

5. Proofs

5.1. The complete graph

We give a simple proof that the complete graph is conformally rigid. This proof is the only proof where we show that the constant weights are the only choice of weights maximizing λ2\lambda_{2} and minimizing λn\lambda_{n}; it is structurally quite different from the other proofs in the paper and makes heavy use of the particular structure of the complete graph.

Proof of Proposition 1.2.

Consider a weighted complete graph on nn vertices with Laplacian LwL_{w} such that ijEwij=|E|=n(n1)/2\sum_{ij\in E}w_{ij}=|E|=n(n-1)/2. If the smallest nonzero eigenvalue of LwL_{w} is nn, then all the nonzero eigenvalues of LwL_{w} are nn since there is always one eigenvalue which is 0 and

tr(Lw)=i=1njiwij=n(n1)=i=1nλi(w)(n1)λ2(w).\mbox{tr}(L_{w})=\sum_{i=1}^{n}\sum_{j\neq i}w_{ij}=n(n-1)=\sum_{i=1}^{n}\lambda_{i}(w)\geq(n-1)\lambda_{2}(w).

Likewise, if the largest nonzero eigenvalue is nn, then all the nonzero eigenvalues are nn since there is always one eigenvalue which is 0 and

tr(Lw)=i=1njiwij=n(n1)=i=1nλi(w)(n1)λn(w).\mbox{tr}(L_{w})=\sum_{i=1}^{n}\sum_{j\neq i}w_{ij}=n(n-1)=\sum_{i=1}^{n}\lambda_{i}(w)\leq(n-1)\lambda_{n}(w).

In either case, there are only two eigenspaces; span(𝟙)\textup{span}(\mathbbm{1}) with eigenvalue 0 and 𝟙\mathbbm{1}^{\perp} with eigenvalue nn. If f:Vf:V\rightarrow\mathbb{R} has mean value 0 (i.e., f𝟙f\in\mathbbm{1}^{\perp}), then

Lwf=nf,L_{w}f=nf,

and if ff is constant, then Lwf=0L_{w}f=0. The matrix Jn×nJ\in\mathbb{R}^{n\times n} of all ones, plays the inverse role: if f:Vf:V\rightarrow\mathbb{R} has mean value 0, then Jf=0nJf=0\in\mathbb{R}^{n}, and if ff is constant, then Jf=nfJf=nf. Therefore

Lw+J=nIL_{w}+J=n\cdot I

which forces w=𝟙w=\mathbbm{1}. ∎

5.2. Edge-transitive graphs

Proof of Proposition 2.1.

If GG is edge-transitive, then the SDP (10) is the one-variable symmetry reduced problem (25):

(31) min {|E|α:αLI1nJ,α0}.\displaystyle\textup{min }\left\{|E|\alpha\,:\,\alpha L\succeq I-\frac{1}{n}J,\,\,\alpha\geq 0\right\}.

The solution α=1/λ2(𝟙)\alpha=1/\lambda_{2}(\mathbbm{1}) is optimal for this problem and hence, w=𝟙/λ2(𝟙)w^{*}=\mathbbm{1}/\lambda_{2}(\mathbbm{1}) is optimal for (10). The argument for λn(w)\lambda_{n}(w) is analogous.

5.3. Distance-Regular Graphs

This section provides some background material and then proves Theorem 2.2. We start by recalling the basic setup. We use d(i,j)d(i,j) to denote the length of the shortest path between vertices ii and jj in a connected graph G=([n],E)G=([n],E). The diameter of GG is the largest distance between any pair of vertices, formally diam(G)=maxi,j[n]d(i,j)\textup{diam}(G)=\max_{i,j\in[n]}d(i,j).

Definition 5.1.

A connected graph GG is called distance-regular if for any two vertices ii and jj

|{u[n]:d(u,i)=p,d(u,j)=q}||\{u\in[n]\,:\,d(u,i)=p,\,\,d(u,j)=q\}|

only depends on d(i,j)d(i,j) and p,qp,q\in\mathbb{N}.

Distance-regular graphs are always regular and the number of distinct distances d(i,j)d(i,j) among the vertices equals diam(G)+1\textup{diam}(G)+1 (since d(i,i)=0d(i,i)=0). We also formally recall the definition of strongly-regular graphs.

Definition 5.2.

A regular graph G=([n],E)G=([n],E) that is not complete or empty is strongly-regular if any two adjacent vertices i,ji,j have aa common neighbors, and any two non-adjacent vertices have bb common neighbors.

A connected strongly-regular graph is a distance-regular graph of diameter 22. We will now quickly survey known results about distance-regular graphs and refer to [11, Chapter 13] for additional results and details.

Let G=([n],E)G=([n],E) be a distance-regular graph with Laplacian LL, and λ\lambda an eigenvalue of LL with multiplicity mm. Consider an embedding 𝒰={u1,,un}m\mathcal{U}=\{u_{1},\ldots,u_{n}\}\subset{\mathbb{R}}^{m} of GG on λ\mathcal{E}_{\lambda} as in Definition 4.1, with the extra requirement that the columns of Un×mU\in{\mathbb{R}}^{n\times m} form an orthonormal basis of λ\mathcal{E}_{\lambda}. The projection matrix UUUU^{\top} (from nλ{\mathbb{R}}^{n}\rightarrow\mathcal{E}_{\lambda}) is the Gram matrix of 𝒰\mathcal{U} with ijij-entry equal to uiuju_{i}^{\top}u_{j}. When GG is distance-regular, the embedding 𝒰\mathcal{U} is very rigid as formalized in the following Lemma.

Lemma 5.3.

[11, Lemma 1.2] Let GG be a distance-regular graph and 𝒰\mathcal{U} be an embedding of GG on the eigenspace λ\mathcal{E}_{\lambda} as described above. If ii and jj are two vertices of GG, then the inner product uiuju_{i}^{\top}u_{j} is determined by d(i,j)d(i,j).

We can use this Lemma to show that, for distance-regular graphs, 𝒰\mathcal{U} is an edge-isometric embedding.

Corollary 5.4.

If GG is a connected distance-regular graph and λ>0\lambda>0, then 𝒰\mathcal{U} is an edge-isometric spherical embedding of GG on λ\mathcal{E}_{\lambda}.

Proof.

By construction, 𝒰\mathcal{U} is an embedding of GG on λ\mathcal{E}_{\lambda}. An immediate consequence of Lemma 5.3 is that, since d(i,i)=0d(i,i)=0, ui2=uiui\|u_{i}\|^{2}=u_{i}^{\top}u_{i} is independent of ii, and hence all uiu_{i} have the same 2\ell^{2}-norm. Since the columns of UU are orthonormal, the square of its Frobenius norm is

UF2=i,juij2=mand thusui=mn.\|U\|_{F}^{2}=\sum_{i,j}u_{ij}^{2}=m\qquad\mbox{and thus}\qquad\|u_{i}\|=\sqrt{\frac{m}{n}}.

In particular, all the points in {u1,,un}m\left\{u_{1},\dots,u_{n}\right\}\subset\mathbb{R}^{m} are on a sphere centered at the origin. It remains to show that if (i,j)E(i,j)\in E, that uiuj\|u_{i}-u_{j}\| is a positive constant independent of i,ji,j. If (i,j)E(i,j)\in E, then d(i,j)=1d(i,j)=1. Applying Lemma 5.3 one more time, we see that there exists a constant c2c_{2} such that

(i,j)Euiuj=c2.(i,j)\in E\implies u_{i}^{\top}u_{j}=c_{2}.

Thus, whenever (i,j)E(i,j)\in E, we have

uiuj2\displaystyle\|u_{i}-u_{j}\|^{2} =ui2+uj22uiuj=2mn2c2\displaystyle=\|u_{i}\|^{2}+\|u_{j}\|^{2}-2u_{i}^{\top}u_{j}=\frac{2m}{n}-2c_{2}

which is indeed a nonnegative constant independent of ii and jj. It remains to show that 2m/n2c22m/n-2c_{2} is a positive constant. We argue by contradiction. If 2m/n2c22m/n-2c_{2} was 0, then we would have ui=uju_{i}=u_{j} for all (i,j)E(i,j)\in E. Since the graph is connected, this means that ui=uju_{i}=u_{j} for all 1i,jn1\leq i,j\leq n. This contradicts the fact that the columns of UU are an orthonormal set as soon as there are at least two columns. So, suppose Un×1U\in\mathbb{R}^{n\times 1} which happens if the dimension of the eigenspace λ\mathcal{E}_{\lambda} is 1. Appealing again to Lemma 5.3, we see that all entries of UU are ±1/n\pm\sqrt{1/n}. Since the the column UU is an eigenvector of the graph Laplacian and λ>0\lambda>0, orthogonality to the constant vector implies that both numbers ±1/n\pm\sqrt{1/n} have to appear (and, in fact, have to appear the same number of times), which forces nn to be even. This suggests the decomposition V=[n]=ABV=[n]=A\cup B where

A={i[n]:ui=1n}andB={i[n]:ui=1n}.A=\left\{i\in[n]:u_{i}=-\frac{1}{\sqrt{n}}\right\}\qquad\mbox{and}\qquad B=\left\{i\in[n]:u_{i}=\frac{1}{\sqrt{n}}\right\}.

Moreover, appealing once more to Lemma 5.3, if (i,j)E(i,j)\in E, then uiuj=c3u_{i}u_{j}=c_{3} for some constant c3c_{3} that is independent of ii and jj. This constant has to be either 1/n1/n or 1/n-1/n. If it were 1/n1/n, then this would mean that there is not a single edge between AA and BB which would mean that the graph is not connected which is a contradiction. Therefore, c3=1/nc_{3}=-1/n. This, in turn, means that all edges run between AA and BB which shows that the embedding is indeed edge-isometric. ∎

The embedding 𝒰\mathcal{U} is in fact a spherical 22-design [11, Corollary 6.2].

Proof of Theorem 2.2.

Proposition 4.3 shows that conformal rigidity is equivalent to having suitable edge-isometric embeddings of the graph on λ2\mathcal{E}_{\lambda_{2}} and λn\mathcal{E}_{\lambda_{n}} which is guaranteed by the preceding Corollary. ∎

We note that this entire approach is specific to distance-regular graphs: one does not in general obtain edge-isometric spherical embeddings via this approach (though, some of the ideas can be used to produce certificates at a greater level of generality, see §6.5). For instance, the claw graph on 44 vertices is edge-transitive and hence conformally rigid. Its Laplacian has two nonzero eigenvalues and the embedding 𝒰\mathcal{U} constructed from an orthonormal basis of either eigenspace is edge-isometric, but not spherical.

5.4. Cayley graphs

Proof of Theorem 2.3.

Suppose there exists an assignment of weights such that λ2(w)>λ2(𝟙)\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}). This means that for every function f:Gf:G\rightarrow\mathbb{R} with mean value 0 and normalized in 2(Γ)\ell^{2}(\Gamma), we have, for non-constant weights wg,sw_{g,s} on the edge connecting gg and gsg\circ s,

sSgGwg,s(f(g)f(gs))2λ2(w)>λ2(𝟙).\sum_{s\in S}\sum_{g\in G}w_{g,s}(f(g)-f(g\circ s))^{2}\geq\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}).

We use φ\varphi to denote an eigenvector corresponding to λ2(𝟙)\lambda_{2}(\mathbbm{1}) when wg,s=1w_{g,s}=1 (and we also assume it to be normalized in 2\ell^{2}). We introduce

φa(v)=φ(av).\varphi_{a}(v)=\varphi(a\circ v).

Then, from the above inequality we have

aGsSgGwg,s(φa(g)φa(gs))2(#G)λ2(w).\sum_{a\in G}\sum_{s\in S}\sum_{g\in G}w_{g,s}(\varphi_{a}(g)-\varphi_{a}(g\circ s))^{2}\geq(\#G)\cdot\lambda_{2}(w).

We can rewrite the left hand side as

aGsSgGwg,s(φa(g)φa(gs))2\displaystyle\sum_{a\in G}\sum_{s\in S}\sum_{g\in G}w_{g,s}(\varphi_{a}(g)-\varphi_{a}(g\circ s))^{2} =sSgGaGwg,s(φa(g)φa(gs))2\displaystyle=\sum_{s\in S}\sum_{g\in G}\sum_{a\in G}w_{g,s}(\varphi_{a}(g)-\varphi_{a}(g\circ s))^{2}
=sSgGaGwa1g,s(φ(g)φ(gs))2\displaystyle=\sum_{s\in S}\sum_{g\in G}\sum_{a\in G}w_{a^{-1}g,s}(\varphi(g)-\varphi(g\circ s))^{2}
=sSgG(aGwa1g,s)(φ(g)φ(gs))2.\displaystyle=\sum_{s\in S}\sum_{g\in G}\left(\sum_{a\in G}w_{a^{-1}g,s}\right)(\varphi(g)-\varphi(g\circ s))^{2}.

This means that we can replace the weight wg,sw_{g,s} by the average in the entire orbit

wg,s1#GaGwag,sw_{g,s}\rightarrow\frac{1}{\#G}\sum_{a\in G}w_{ag,s}

and still obtain a function with the desired property of exceeding λ2(w)>λ2(𝟙)\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}). Thus we may assume that wg,sw_{g,s} only depends on ss, and

sSwsgG(φ(g)φ(gs))2λ2(w)>λ2(𝟙).\sum_{s\in S}w_{s}\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2}\geq\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}).

Squaring out, we see that (invoking the 2\ell^{2}-normalization of φ\varphi)

gG(φ(g)φ(gs))2\displaystyle\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2} =gGφ(g)2+φ(gs)22φ(g)φ(gs)\displaystyle=\sum_{g\in G}\varphi(g)^{2}+\varphi(g\circ s)^{2}-2\varphi(g)\varphi(g\circ s)
=22gGφ(g)φ(gs).\displaystyle=2-2\sum_{g\in G}\varphi(g)\varphi(g\circ s).

Now suppose this last sum is independent of ss, i.e., for some c>0c>0 and all sSs\in S

c=gG(φ(g)φ(gs))2.c=\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2}.

Therefore

sSwsgG(φ(g)φ(gs))2=csSws.\sum_{s\in S}w_{s}\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2}=c\sum_{s\in S}w_{s}.

A sum of finitely many real numbers does not change if we replace each summand by the global average and thus

sSwsgG(φ(g)φ(gs))2\displaystyle\sum_{s\in S}w_{s}\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2} =csSws=csS(1|S|tSwt)\displaystyle=c\sum_{s\in S}w_{s}=c\sum_{s\in S}\left(\frac{1}{|S|}\sum_{t\in S}w_{t}\right)
=sS(1|S|tSwt)gG(φ(g)φ(gs))2\displaystyle=\sum_{s\in S}\left(\frac{1}{|S|}\sum_{t\in S}w_{t}\right)\sum_{g\in G}(\varphi(g)-\varphi(g\circ s))^{2}
=sSgG(1|S|tSwt)(φ(g)φ(gs))2.\displaystyle=\sum_{s\in S}\sum_{g\in G}\left(\frac{1}{|S|}\sum_{t\in S}w_{t}\right)(\varphi(g)-\varphi(g\circ s))^{2}.

However, at this point the weights are independent of ss and thus constant. This is a contradiction to λ2(w)>λ2(𝟙)\lambda_{2}(w)>\lambda_{2}(\mathbbm{1}). The argument for λn\lambda_{n} is the same. ∎

5.5. Proof of Proposition 2.5

The proof makes use of the fact that eigenvalues and eigenvectors of circulant graphs are well understood. In particular, in the case when all the weights are equal,

λ2(Cn({1,2}))=42cos(2πn)2cos(4πn).\lambda_{2}(C_{n}\left(\left\{1,2\right\}\right))=4-2\cos\left(\frac{2\pi}{n}\right)-2\cos\left(\frac{4\pi}{n}\right).

One way of seeing this is to realize that the graph Laplacian can be written as a circulant matrix

L=(c0cn1cn2c1c1c0cn1c2c2c1c0c3cn1cn2cn3c0)L=\begin{pmatrix}c_{0}&c_{n-1}&c_{n-2}&\dots&c_{1}\\ c_{1}&c_{0}&c_{n-1}&\dots&c_{2}\\ c_{2}&c_{1}&c_{0}&\dots&c_{3}\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ c_{n-1}&c_{n-2}&c_{n-3}&\dots&c_{0}\\ \end{pmatrix}

where c0=4c_{0}=4 and c1=c2=cn1=cn=1c_{1}=c_{2}=c_{n-1}=c_{n}=-1 while all other cic_{i} are 0. It is known that the eigenvalues of such a circulant matrix are given by

λj=c0+c1ωj+c2ω2j++cn1ω(n1)jwhereω=exp(2πin),\lambda_{j}=c_{0}+c_{1}\omega^{j}+c_{2}\omega^{2j}+\dots+c_{n-1}\omega^{(n-1)j}\qquad\mbox{where}~{}\omega=\exp\left(\frac{2\pi i}{n}\right),

where jj runs from 0jn10\leq j\leq n-1. In our case, this expression simplifies to

λj\displaystyle\lambda_{j} =4ωjω2jω(n1)jω(n2)j\displaystyle=4-\omega^{j}-\omega^{2j}-\omega^{(n-1)j}-\omega^{(n-2)j}
=4cos(2πjn)cos(4πjn)cos(2π(n1)jn)cos(2π(n2)jn)\displaystyle=4-\cos\left(\frac{2\pi j}{n}\right)-\cos\left(\frac{4\pi j}{n}\right)-\cos\left(\frac{2\pi(n-1)j}{n}\right)-\cos\left(\frac{2\pi(n-2)j}{n}\right)
=42cos(2πjn)2cos(4πjn).\displaystyle=4-2\cos\left(\frac{2\pi j}{n}\right)-2\cos\left(\frac{4\pi j}{n}\right).

The case j=0j=0 corresponds to the eigenvalue λ1=0\lambda_{1}=0. To make the cosine as large as possible, we need the arguments in the cosine to be as close as possible to a multiple of 2π2\pi and this leads to j=1j=1 and j=n1j=n-1. We will now investigate a different choice of weights on the circulant that maintains the circulant structure. For a ε>0\varepsilon>0 to be thought of as very small, we set c0=4c_{0}=4 as well as

c1=cn1=1+εandc2=cn2=1ε.c_{1}=c_{n-1}=-1+\varepsilon\qquad\mbox{and}\qquad c_{2}=c_{n-2}=-1-\varepsilon.

Then

(32) λj=4(2ε)cos(2πjn)(2+ε)cos(4πjn).\displaystyle\lambda_{j}=4-(2-\varepsilon)\cos\left(\frac{2\pi j}{n}\right)-(2+\varepsilon)\cos\left(\frac{4\pi j}{n}\right).

Our goal will be to show that for a small value of ε\varepsilon we can ensure that λ2\lambda_{2} increases. However, it is clear that some restriction on nn will be required. When n=5n=5, then C5({1,2})=K5C_{5}\left(\left\{1,2\right\}\right)=K_{5} is conformally rigid. Moreover, C6({1,2})C_{6}\left(\left\{1,2\right\}\right) is the octahedral graph which is conformally rigid (it is edge-transitive). In both of these cases, the procedure cannot work and it is easy to see from the spectrum why: K5K_{5} has eigenvalues (0,5,5,5,5)(0,5,5,5,5) and a perturbation induced by ε\varepsilon affects all eigenvalues and some will increase while others decrease. Something similar is happening when n=6n=6, the octahedral graph has spectrum (0,4,4,4,6,6)(0,4,4,4,6,6), in particular, the multiplicity of the second eigenvalue is 3, this allows for it to be conformally rigid.

Lemma 5.5.

For every n7n\geq 7 there exists ε0>0\varepsilon_{0}>0 so that for all ε(0,ε0)\varepsilon\in(0,\varepsilon_{0}), the multiplicity of λ2\lambda_{2} is 2.

Proof.

It suffices to prove the result when ε=0\varepsilon=0 and then appeal to continuity of the expression in ε\varepsilon. Note first that the expression (32) is the same for jj and j-j and the expression is nn-periodic and thus it is the same for jj and njn-j, so it suffices to understand the behavior when jn/2j\leq n/2. We introduce the function

f(x)=42cos(2πx)2cos(4πx)f(x)=4-2\cos\left(2\pi x\right)-2\cos\left(4\pi x\right)

and it suffices to understand this function for 0x1/20\leq x\leq 1/2 (see Fig. 8). A quick computation shows that

f(16)=f(12)=4.f\left(\frac{1}{6}\right)=f\left(\frac{1}{2}\right)=4.

This appearance of 1/61/6 corresponds exactly to the conformal rididity of the octahedral graph. The function is monotonically increasing on (0,1/6)(0,1/6) and thus, once n7n\geq 7, we have

min1jn/2f(jn)=f(1/n)\min_{1\leq j\leq n/2}f\left(\frac{j}{n}\right)=f(1/n)

which is attained only when j=1j=1 and the remaining case arises, from symmetry, when j=n1j=n-1. This proves the result for ε=0\varepsilon=0, the existence of an entire range (0,ε0)(0,\varepsilon_{0}) follows from continuity. ∎

Refer to captionRefer to caption
Figure 8. The octahedral graph C6({1,2})C_{6}\left(\left\{1,2\right\}\right) and the function f(x)f(x) for 0x1/20\leq x\leq 1/2 with level set f(1/6)=f(1/2)=4f(1/6)=f(1/2)=4.

Having proved this Lemma, the argument now follows relatively easily. We have once n7n\geq 7, for an entire range 0εε00\leq\varepsilon\leq\varepsilon_{0} that

λ2(ε)=4(2ε)cos(2πn)(2+ε)cos(4πn).\lambda_{2}(\varepsilon)=4-(2-\varepsilon)\cos\left(\frac{2\pi}{n}\right)-(2+\varepsilon)\cos\left(\frac{4\pi}{n}\right).

This expression is increasing and thus there is a way of changing the weights to increase λ2\lambda_{2} which shows that Cn({1,2})C_{n}\left(\left\{1,2\right\}\right) is not conformally rigid for n7n\geq 7.

6. Certifying Conformal Rigidity via SDPs

In this section we certify the conformal rigidity of several graphs mentioned in § 2.6 via semidefinite programming (see § 3). The solutions returned by an SDP solver are typically numerical and only approximately true. Often with some massaging they can be turned into verifiable proofs of conformal rigidity. We explain our methods, which rely on the theory explained in § 3.

6.1. The Rationalization Trick.

Solving the SDP (10) numerically on a computer, we obtain an approximate (primal) solution, an approximate (dual) solution, and a duality gap. In practice, we considered a duality gap of <107<10^{-7} as reasonable numerical evidence that the graph under consideration is indeed conformally rigid. However, this is not a proof. If the graph is conformally rigid, then the constant assignment of weights

wij=1λ2(G)andwij=1λn(G)w_{ij}=\frac{1}{\lambda_{2}(G)}\qquad\mbox{and}\qquad w_{ij}=\frac{1}{\lambda_{n}(G)}

is an optimal (primal) solution for the SDPs (10) and (16) respectively. Their optimality can be certified by optimal solutions XX and YY to the dual SDPs (11) and (17), via strong duality. We first obtain approximate numerical dual optimal solutions XX and YY using the SDP solver MOSEK [23] which we called through Mathematica [22]. If there is an exact solution comprised of rational numbers, then approximating the numerical entries by nearby rational numbers can lead to a closed form expression for true optimal solutions. We found this method to be remarkably effective in cases where λ2(𝟙)\lambda_{2}(\mathbbm{1}) and λn(𝟙)\lambda_{n}(\mathbbm{1}) were both rational numbers.

We first illustrate the procedure for the Hoffman graph (see Fig. 9). The solver produced a dual solution XX, a matrix with entries like -0.434896 and 0.565104. A little bit of an inspection shows that only five different numbers truly appear and that the numbers are exceedingly close to

the rational numbers{359384,167384,25384,217384,409384}.\mbox{the rational numbers}\qquad\left\{-\frac{359}{384},-\frac{167}{384},\frac{25}{384},\frac{217}{384},\frac{409}{384}\right\}.

Replacing the numerical entries in the matrix by these entries (see Fig. 9), we can quickly check (symbolically) that this is indeed a dual optimal XX. The same procedure can be repeated for λn\lambda_{n} where the numerical dual solution only consists of the numbers 0.25-0.25 and 0.250.25 arranged in a block structure (see Fig. 9). Again, once a candidate has been identified, verifying that it is an optimal YY is easy.

Refer to captionRefer to captionRefer to caption
Figure 9. Left: the Hoffman graph. Middle: the dual optimal solution XX for λ2\lambda_{2}. Right: the dual optimal solution YY for λn\lambda_{n}.

It should be noted that the two aspects that may cause complications are as follows: it is not a priori clear whether there is a dual optimum that is only comprised of rational numbers and, even if there is, it is not a priori clear how large the denominators are (equivalent to asking what type of accuracy one requires of the numerical approximation). On the other hand, once a candidate for an optimal solution has been identified, it is easy to prove that it is indeed one.

Refer to captionRefer to captionRefer to caption
Figure 10. Left: the CrossingNumberGraph6B. Middle: a numerical certificate for λ2\lambda_{2}. Right: the certificate for λn\lambda_{n}.

Another example is, CrossingNumberGraph6B, where the rationalization method by itself is not successful (see Fig. 10). The Laplacian has eigenvalues λ2=1\lambda_{2}=1 (with multiplicity 3) and λn=6\lambda_{n}=6 (with multiplicity 1). While the rationalization method works for λn\lambda_{n} (where the solution is particularly easy), it does not work for λ2\lambda_{2}. In the next section we will use a more sophisticated method to obtain an optimal solution XX. It was mentioned in [27] that if XX is feasible for the dual SDP (11), then so is (I(1/n)J)X(I(1/n)J)(I-(1/n)J)X(I-(1/n)J). Same applies to YY. For many graphs, applying this transformation to the dual solutions returned by the SDP solver was enough to make them nice and verifiable.

6.2. The Projection Method.

There is an alternate method that makes use of equations (13) and (19); the rows (and columns) of the dual optima X,YX,Y are linear combinations of eigenvectors of the Laplacian matrix corresponding to λ2\lambda_{2} and λn\lambda_{n}, respectively. This suggests a natural idea.

  1. (1)

    Compute a numerical approximation for the dual optimum.

  2. (2)

    Project each row into the appropriate eigenspace of the Laplacian matrix.

  3. (3)

    Convert the coefficients into reasonable nearby numbers (either rational numbers or algebraic numbers of low degree, like 2\sqrt{2}).

  4. (4)

    Use this to recover an actual solution.

We illustrate on the solution XX (for λ2\lambda_{2}) for CrossingNumberGraph6B. The first row of a numerical candidate starts as

(2.50568,1.67231,1.67231,0.00569952,)(2.50568,1.67231,1.67231,0.00569952,\dots)

and one might be tempted to think of these numbers as 5/2,5/3,5/3,0,5/2,5/3,5/3,0,\dots but this does not work. However, there exists a particularly nice basis for the smallest eigenspace of λ2\lambda_{2} given by the vectors

u1\displaystyle u_{1} =(0,0,4,3,1,3,1,1,1,2,2,1,1,1,1,4,3,3,0,0)\displaystyle=(0,0,4,-3,1,-3,1,1,1,-2,2,-1,-1,-1,-1,-4,3,3,0,0)
u2\displaystyle u_{2} =(0,3,1,0,2,0,2,1,1,2,2,1,1,2,2,1,0,0,3,0)\displaystyle=(0,-3,-1,0,2,0,2,-1,-1,2,-2,1,1,-2,-2,1,0,0,3,0)
u3\displaystyle u_{3} =(3,0,1,0,1,0,1,2,2,2,2,2,2,1,1,1,0,0,0,3).\displaystyle=(-3,0,-1,0,-1,0,-1,2,2,2,-2,-2,-2,1,1,1,0,0,0,3).

Note that these vectors are not orthogonal but they are in closed form. Computing the inner product of the first row of the numerical XX with these three vectors, we get the numbers (19.9994,19.9994,34.9999)(19.9994,-19.9994,-34.9999) and it’s pretty clear what these are supposed to be. Correcting them and solving the linear system, we obtain a corrected matrix whose first row starts with

73,32,32,16,16,\frac{7}{3},\frac{3}{2},\frac{3}{2},-\frac{1}{6},-\frac{1}{6},\dots

which isn’t even particularly close to the original numerical values. This then turns out to indeed be a correct dual optimal solution XX.

This method also worked for the (20,8)(20,8)-accordion graph (see Fig. 7); the solution XX for λn\lambda_{n} is trivial, the solution YY for λ2\lambda_{2} is not. The numerical solution returned by the solver is not very insightful, however, projecting into the eigenspace corresponding to λ2\lambda_{2} is very helpful because it only has multiplicity 2: the new coordinates happened to either be 0 or, in case they were both nonzero, happen to have a ratio of (1+5)/2(1+\sqrt{5})/2 which helped identify the correct solution.

If the multiplicity of either λ2\lambda_{2} or λn\lambda_{n} is one, then there is a particularly easy dual optimum constructed from an eigenvector of that eigenvalue which is a special case of the projection method. For the Haar Graph 565 shown in Fig. 11 with 2020 vertices, we used the projection method to find XX; this graph has λ2=2.76393\lambda_{2}=2.76393 with multiplicity 44. However, λ20=10\lambda_{20}=10 with multiplicity 11, and so we can use (22) to find an optimal Y𝒮+20Y\in\mathcal{S}^{20}_{+} for (16) which has optimal value q=5q^{*}=5 (see the next subsection). An eigenvector of λ20\lambda_{20} is

u=(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)u=(-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1,1,1)

which has norm 20\sqrt{20}. Therefore Y=5(1/20)uuY=5(1/20)uu^{\top} is a dual optimal solution for (16) with trace equal to 5=q5=q^{*}.

6.3. Symmetry Reduction

Conformal rigidity of the Haar Graph 565 is especially easy via symmetry reduction. Haar Graph 565 has 5050 edges which fall into 22 orbits of size 4040 and 1010 under the automorphism group of the graph. This Haar graph and its edge orbit graphs are shown in Fig. 11.

Refer to captionRefer to captionRefer to caption
Figure 11. Haar 565 (left) and its decomposition (middle and right) under the automorphism group of the graph.

The symmetry reduced SDP (25) has 22 variables and the constraint

w1L1+w2L2I120J.w_{1}L_{1}+w_{2}L_{2}\succeq I-\frac{1}{20}J.

The eigenvalues of w1L1+w2L2I+120Jw_{1}L_{1}+w_{2}L_{2}-I+\frac{1}{20}J are linear in w1,w2w_{1},w_{2} as shown below and can be computed explicitly in Mathematica:

(33) 0,1+4w1,1+(5±5)w1,1+(3±5)w1+2w2,1+4w1+2w2,1+8w1+2w2\displaystyle\begin{split}0,\,-1+4w_{1},\,-1+(5\pm\sqrt{5})w_{1},\,-1+(3\pm\sqrt{5})w_{1}+2w_{2},\\ -1+4w_{1}+2w_{2},-1+8w_{1}+2w_{2}\end{split}

Therefore, problem (10) is to minimize 40w1+10w240w_{1}+10w_{2} subject to all eigenvalues in (33) being nonnegative and w1,w20w_{1},w_{2}\geq 0. This is a 22-variable linear program with an optimal solution

w1=w2=1λ2(𝟙)0.361803.w_{1}=w_{2}=\frac{1}{\lambda_{2}(\mathbbm{1})}\sim 0.361803.

Therefore, 𝟙/λ2(𝟙)\mathbbm{1}/\lambda_{2}(\mathbbm{1}) is an optimal solution of (10). To solve (16), we maximize 40w1+10w240w_{1}+10w_{2} subject to all the eigenvalues in (33) being nonpositive and w1,w20w_{1},w_{2}\geq 0. This is again a 22-variable linear program with an optimal solution w1=w2=110w_{1}=w_{2}=\frac{1}{10} and optimal value q=5q^{*}=5. Therefore, 𝟙/λ20(𝟙)\mathbbm{1}/\lambda_{20}(\mathbbm{1}) is an optimal solution of (16) and Haar Graph 565 is conformally rigid. The solver we used returned the optimal solution w1=0,w2=1/2w_{1}=0,w_{2}=1/2 which is not useful to prove conformal rigidity but is useful to verify the optimality of w1=w2=1/10w_{1}=w_{2}=1/10. It also provides an example of a graph where a multiple of 𝟙\mathbbm{1} is not the only weight that minimizes λn(w)\lambda_{n}(w).

A more involved example is the distance-22 graph of the Klein graph (see Fig. 7) which has n=24n=24 vertices and 168168 edges that come in 22 orbits each of size 8484. In this case, λ2(𝟙)=11.3542\lambda_{2}(\mathbbm{1})=11.3542 and λn(𝟙)=16.6458\lambda_{n}(\mathbbm{1})=16.6458, each with multiplicity 88. The eigenvalues of w1L1+w2L2I+124Jw_{1}L_{1}+w_{2}L_{2}-I+\frac{1}{24}J are:

(34) 0,1+8w1+8w2,1+7w1+7w2±7(w12w1w2+w22)\displaystyle\begin{split}0,\,\,-1+8w_{1}+8w_{2},\,\,-1+7w_{1}+7w_{2}\pm\sqrt{7(w_{1}^{2}-w_{1}w_{2}+w_{2}^{2})}\\ \end{split}

Therefore, Problem (10) is to minimize 84w1+84w284w_{1}+84w_{2} subject to the nonnegativity of the eigenvalues and w1,w20w_{1},w_{2}\geq 0. This is a 22-dimensional convex optimization problem with optimal solution w1=w2=1/λ2(𝟙)0.0880728w_{1}=w_{2}=1/\lambda_{2}(\mathbbm{1})\sim 0.0880728. The feasible region and the optimal position of the objective function are shown in Fig. 12.

Refer to captionRefer to caption
Figure 12. Optimizing λ2(w)\lambda_{2}(w) (left) and λn(w)\lambda_{n}(w) (right) on the distance-22 graph of the Klein graph, after symmetry reduction.

Similarly, Problem (16) is the 22-dimensional convex optimization problem of maximizing 84w1+84w284w_{1}+84w_{2} subject to the nonpositivity of the eigenvalues and w1,w20w_{1},w_{2}\geq 0. It has optimal solution w1=w2=1/λn(𝟙)0.0600754w_{1}=w_{2}=1/\lambda_{n}(\mathbbm{1})\sim 0.0600754. The feasible region and the optimal position of the objective function are shown in Fig. 12.

6.4. Complementary Slackness

Proposition 3.3 provides certificates XX and YY for the conformal rigidity of a graph. This approach can be especially effective. The conditions LX=λ2XLX=\lambda_{2}X and LY=λnYLY=\lambda_{n}Y allows us to parametrize XX and YY and reduce the number of variables. We explain the parametrization of XX; the parametrization of YY is similar.

Lemma 6.1.

Set U:=[u1uk]n×kU:=[u_{1}\ldots u_{k}]\in{\mathbb{R}}^{n\times k} where {u1,,uk}\{u_{1},\ldots,u_{k}\} is a set of kk vectors in n\mathbb{R}^{n} that form a basis for the eigenspace λ2\mathcal{E}_{\lambda_{2}} of LL. Then X𝒮+nX\in\mathcal{S}^{n}_{+} and LX=λ2XLX=\lambda_{2}X if and only if S𝒮+kS\in\mathcal{S}^{k}_{+} and X=USUX=USU^{\top}.

Proof.

If X=USUX=USU^{\top} and S0S\succeq 0, then X0X\succeq 0 since for all vnv\in{\mathbb{R}}^{n}

vXv=vUSUv=(Uv)S(Uv)0.v^{\top}Xv=v^{\top}USU^{\top}v=(U^{\top}v)^{\top}S(U^{\top}v)\geq 0.

The columns of XX are combinations of the columns of UU and so, LX=λ2XLX=\lambda_{2}X.

Now suppose X0X\succeq 0 and LX=λ2XLX=\lambda_{2}X. Then rank(X)k\operatorname{rank}(X)\leq k where k=mult(λ2)k=\textup{mult}(\lambda_{2}). Since X0X\succeq 0, XX has a reduced singular value decomposition of the form X=ADAX=ADA^{\top} where the columns of AA are orthonormal and Dk×kD\in{\mathbb{R}}^{k\times k} is a diagonal matrix with nonnegative diagonal entries. The column space of XX is the column space of AA which means that A=UBA=UB for some matrix Bk×kB\in{\mathbb{R}}^{k\times k}. Therefore,

X=UBD(UB)=U(BDB)U=USU.X=UBD(UB)^{\top}=U(BDB^{\top})U^{\top}=USU^{\top}.

The matrix S=BDBS=BDB^{\top} is psd since DD is psd. ∎

Note that in Lemma 6.1 we do not require that the columns of UU form an orthonormal basis of λ2\mathcal{E}_{\lambda_{2}}, just that it is a basis. Using Lemma 6.1, we can rephrase Proposition 3.3 as follows.

Corollary 6.2.

Let {u1,,uk}\{u_{1},\ldots,u_{k}\} be a basis of the eigenspace of LL with eigenvalue λ2\lambda_{2} and U=[u1uk]n×kU=[u_{1}\cdots u_{k}]\in{\mathbb{R}}^{n\times k}. Let {v1,,v}\{v_{1},\ldots,v_{\ell}\} be a basis of the eigenspace of LL with eigenvalue λn\lambda_{n} and V=[v1v]n×V=[v_{1}\cdots v_{\ell}]\in{\mathbb{R}}^{n\times\ell}. Then GG is conformally rigid if and only if there exists

  1. (1)

    a matrix X=USUX=USU^{\top} such that

    (35) S𝒮+k,Xii+Xjj2Xij=1ijE,\displaystyle S\in\mathcal{S}^{k}_{+},\,\,\,\,X_{ii}+X_{jj}-2X_{ij}=1\,\,\forall ij\in E,
  2. (2)

    a matrix Y=VTVY=VTV^{\top} such that

    (36) T𝒮+,Yii+Yjj2Yij=1ijE.\displaystyle T\in\mathcal{S}^{\ell}_{+},\,\,\,\,Y_{ii}+Y_{jj}-2Y_{ij}=1\,\,\forall ij\in E.
Proof.

By construction, 𝟙X𝟙=𝟙USU𝟙=0\mathbbm{1}^{\top}X\mathbbm{1}=\mathbbm{1}^{\top}USU^{\top}\mathbbm{1}=0 since the columns of UU are orthogonal to 𝟙\mathbbm{1}. Similarly for YY. All other constraints in Proposition 3.3 are satisfied by XX (and YY) under the given parametrization. ∎

Example 6.3.

We can use the above method to prove the conformal rigidity of the Hoffman graph (see Fig. 1) which has 1616 vertices and 3232 edges. Its Laplacian eigenvalues with multiplicities are:

(λ,mult(λ)):(8,1),(6,4),(4,6),(2,4),(0,1).(\lambda,\,\textup{mult}(\lambda)):\,\,(8,1),(6,4),(4,6),(2,4),(0,1).

The rows of the following matrix span 2\mathcal{E}_{2}. Hence, UU^{\top} is

(1111111100022000111111110020020011111111020000201111111120000002)\small{\left(\begin{array}[]{cccccccccccccccc}-1&1&-1&-1&-1&1&1&1&0&0&0&-2&2&0&0&0\\ -1&-1&1&1&-1&1&1&-1&0&0&-2&0&0&2&0&0\\ -1&-1&1&-1&1&1&-1&1&0&-2&0&0&0&0&2&0\\ -1&-1&-1&1&1&-1&1&1&-2&0&0&0&0&0&0&2\\ \end{array}\right)}

Setting SS to be a symbolic symmetric matrix of size 4×44\times 4, and X=USUX=USU^{\top} we can solve the linear equations Xii+Xjj2Xij=1X_{ii}+X_{jj}-2X_{ij}=1 for all ijEij\in E. Plugging the solution back into SS we get the matrix:

S=(s110000s220000s3300001s11s22s33)S=\left(\begin{array}[]{cccc}s_{11}&0&0&0\\ 0&s_{22}&0&0\\ 0&0&s_{33}&0\\ 0&0&0&1-s_{11}-s_{22}-s_{33}\\ \end{array}\right)

which is psd if and only if all its diagonal entries are nonnegative. Therefore, there are infinitely many choices for the certificate XX. Picking s11=s22=s33=1/4s_{11}=s_{22}=s_{33}=1/4, we get the certificate XX in Fig. 13. The eigenspace of 8\mathcal{E}_{8} is spanned by

(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),(-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1),

and running the same procedure, we get the unique certificate YY in Fig. 13.

Refer to captionRefer to caption
Figure 13. Certificates for the conformal rigidity of the Hoffman graph. Left: a certificate XX for λ2\lambda_{2}, Right: the certificate YY for λn\lambda_{n}.

We use Corollary 4.7 and the above parametrizations to show that while C6C_{6} is conformally rigid, its complement is not.

Example 6.4.

The cycle C6C_{6} is conformally rigid since it is edge-transitive. Let G=C6cG=C_{6}^{c} which is not edge-transitive and has 99 edges. Its Laplacian has eigenvalues 0,2,3,3,5,50,2,3,3,5,5. The eigenspace 5c\mathcal{E}^{c}_{5} is spanned by the orthonormal vectors

u=(1/(23),1/3,1/(23),1/(23),1/3,1/(23),\displaystyle u=(-1/(2\sqrt{3}),-1/\sqrt{3},-1/(2\sqrt{3}),1/(2\sqrt{3}),1/\sqrt{3},1/(2\sqrt{3}),
v=(1/2,0,1/2,1/2,0,1/2).\displaystyle v=(1/2,0,-1/2,-1/2,0,1/2).

Setting

T=(abbc),T=\begin{pmatrix}a&b\\ b&c\end{pmatrix},

consider a potential certificate Y=[uv]T[u,v]Y^{\prime}=[u\,v]\,T\,[u,v]^{\top}. If C6cC_{6}^{c} were conformally rigid, we would have Trace Y=95=Trace (T)=a+c\textup{Trace }Y^{\prime}=\frac{9}{5}=\textup{Trace }(T)=a+c. Therefore, we can set c=9/5ac=9/5-a in YY^{\prime}. Then imposing the linear conditions, we get from the edge (1,3)(1,3) of C6cC_{6}^{c} that a+3b=7/10a+\sqrt{3}b=7/10 and from the edge (2,5)(2,5) of C6cC_{6}^{c} that a+3b=3/5a+\sqrt{3}b=3/5 which is a contradiction. So no such YY^{\prime} exists and C6cC_{6}^{c} is not conformally rigid.∎

6.5. Using UUUU^{\top}

We conclude by noting a particularly simple method that often ends up producing certificates XX and YY (and provably does so for distance-regular graphs). For λ2\lambda_{2}, the procedure works as follows: one collects an orthonormal basis of the eigenspace λ2\lambda_{2} corresponding to the Laplacian LL in a matrix UU. A rescaling of the matrix UUcUUUU^{\top}\rightarrow cUU^{\top} to ensure that the trace is |E|/λ2|E|/\lambda_{2} often produces certificates that seem to work: in particular, running a check on the graphs built into the internal database of Mathematica shows that this approach leads to a valid certificate (for λ2\lambda_{2}) for all conformally rigid graphs having less than 20 vertices. A case where this does not work is the Haar graph 565 on 20 vertices (see Fig. 11). However, there are other examples such as the non-Cayley vertex-transitive graph (24,23) on n=24n=24 vertices, where all the other approaches failed while cUUcUU^{\top} indeed produces a valid certificate for λ2\lambda_{2}.

Acknowledgment. SS was supported by the NSF (DMS-2123224). RT acknowledges helpful discussions with Sameer Agarwal.

References

  • [1] C. Babecki, S. Steinerberger and R. Thomas, Spectrahedral Geometry of Graph Sparsifiers, SIAM J. Disc. Math, 39(1) (2025), 449–483.
  • [2] J. Berggren, An algebraic characterization of symmetric graphs with a prime number of vertices, Bull. Austral. Math. Soc. 7 (1972), 131–134.
  • [3] S. Boyd, P. Diaconis, P. Parrilo, and L. Xiao, Fastest mixing Markov chain on graphs with symmetries. SIAM Journal on Optimization, 20(2) (2009), 792–819.
  • [4] S.P. Boyd and L. Vandenberghe. Convex Optimization. Cambridge University Press, 2004.
  • [5] C. Y. Chao, On the classification of symmetric graphs with a prime number of vertices, Trans. Amer. Math. Soc. 158 (1971), 247–256.
  • [6] Y. Cheng and J. Oxley, On weakly symmetric graphs of order twice a prime. J. Combin. Theory Ser. B 42 (1987), 196–211.
  • [7] A. El Soufi and S. Ilias, Immersions minimales, premiere valeur propre du laplacien et volume conforme, Math. Ann., 275:2 (1986), 257–267.
  • [8] K. Fan, O. Taussky, J. Todd, Discrete analogs of inequalities of Wirtinger. Monatsh. Math.59 (1955), 73–90.
  • [9] M. Fiedler, Absolute algebraic connectivity of trees. Linear and Multilinear Algebra, 26 (1990), 85-106.
  • [10] M. Fiedler, Some minimax problems for graphs. Discrete Mathematics 121.1-3 (1993), 65-74.
  • [11] C. Godsil. Algebraic Combinatorics. Chapman & Hall, Inc., 1993.
  • [12] F. Göring, C. Helmberg, and M. Wappler. Embedded in the shadow of the separator. SIAM Journal on Optimization 19.1 (2008), 472-501.
  • [13] F. Göring, C. Helmberg, and S. Reiss. Graph realizations associated with minimizing the maximum eigenvalue of the Laplacian. Mathematical Programming 131 (2012), 95-111.
  • [14] F. Göring, C. Helmberg, and S. Reiss. On minimizing the spectral width of graph laplacians and associated graph realizations. SIAM Journal on Optimization 23.2 (2013), 834-856.
  • [15] R. Horn and C. Johnson, Matrix Analysis. Cambridge University Press, Cambridge, 1985.
  • [16] J. B. Gauci and J. P. Zerafa, Accordion graphs: Hamiltonicity, matchings and isomorphism with quartic circulants. Discrete Applied Mathematics, 321 (2022), 126-137.
  • [17] K.M. Hall, An rr-dimensional quadratic placement algorithm. Management Science, 17 (1970), 219-229.
  • [18] J. Hersch, Quatre proprietes isoperimetriques de membranes spheriques homogenes, C. R. Acad. Sci. Paris Ser. A-B, 270 (1970), A1645–A1648.
  • [19] I. Izmestiev, A general discrete Wirtinger inequality and spectra of discrete Laplacians. arXiv preprint arXiv:1502.03186.
  • [20] N. Korevaar, Upper bounds for eigenvalues of conformal metrics, Journal of Differential Geometry 37 (1993), 73–93.
  • [21] G. Milovanovic and I. Milovanovic, On discrete inequalities of Wirtinger’s type. J. Math. Anal. Appl. 88 (1982), 378–387.
  • [22] Wolfram Mathematica, Wolfram Research, https://www.wolfram.com/mathematica/
  • [23] MOSEK, https://www.mosek.com/
  • [24] E. Pegg and G. Exoo, Crossing number graphs. The Mathematica Journal 11 (2009), 161-170.
  • [25] C. Praeger, R. Wangand M. Y. Xu, Symmetric graphs of order a product of two distinct primes. J. Combin. Theory Ser. B 58 (1993), 299–318.
  • [26] I. J. Schoenberg, The finite Fourier series and elementary geometry. American Mathematical Monthly (1950), 390-404.
  • [27] J. Sun, S. Boyd, L. Xiao and P. Diaconis, The fastest mixing Markov process on a graph and a connection to a maximum variance unfolding problem. SIAM Review 48 (2006), 681-699.
  • [28] R. J. Wang and M. Y. Xu, A classification of symmetric graphs of order 3p. J. Combin. Theory Ser. B 58 (1993), 197–216
  • [29] H. Zhang, On edge transitive circulant graphs. Tokyo J. of Mathematics, 19 (1996), p. 51-55.
  • [30] X. Zhang, A refinement of the discrete Wirtinger inequality. Journal of mathematical analysis and applications, 200 (1996), 687-697.