A transfer principle and applications to eigenvalue estimates for graphs Omid Amini, David Cohen-Steiner

To cite this version:

Omid Amini, David Cohen-Steiner. A transfer principle and applications to eigenvalue estimates for graphs. Commentarii Mathematici Helvetici, European Mathematical Society, 2018. ￿hal-01109634v2￿

HAL Id: hal-01109634 https://hal.inria.fr/hal-01109634v2 Submitted on 29 Mar 2019

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS

OMID AMINI AND DAVID COHEN-STEINER

Abstract. In this paper, we prove a variant of the Burger-Brooks transfer principle which, combined with recent eigenvalue bounds for surfaces, allows to obtain upper bounds on the eigenvalues of graphs as a function of their genus. More precisely, we show the existence of nr a universal constants C such that the k-th eigenvalue λk of the normalized Laplacian of a graph G of (geometric) genus g on n vertices satisfies d (g + k) λnr(G) ≤ C max , k n where dmax denotes the maximum valence of vertices of the graph. Our result is tight up to a change in the value of the constant C, and improves recent results of Kelner, Lee, Price and Teng on bounded genus graphs. To show that the transfer theorem might be of independent interest, we relate eigenvalues of the Laplacian on a metric graph to the eigenvalues of its simple graph models, and discuss an application to the mesh partitioning problem, extending results of Miller-Teng-Thurston- Vavasis and Spielman-Teng to arbitrary meshes.

1. Introduction The spectrum of the Laplacian of a finite graph reflects information about the structural properties of the graph and has been successfully used in a large variety of applications to other domains. In particular, the eigenvalues of a bounded degree graph provide information on the existence of good clusterings of that graph, see [1] for clustering in two classes and [20, 22] for k-way clusterings, whose optimal quality is shown to relate to the k-th eigenvalue. In particular, upper bounds on the eigenvalues of a class of graphs directly translate into effi- cient clustering algorithms with quality guarantees. This motivated a series of work, starting with Spielman and Teng [27], who gave an O(1/n) bound for the Fiedler value of a bounded degree planar graph on n vertices, using a suitably centered circle packing representation of the graph. Kelner extended this result to an O((g + 1)/n) bound for (geometric) genus g graphs [16]. The argument uses Riemann-Roch theorem to find a circle packing representa- tion of the graph. Recently, Kelner, Lee, Price and Teng proved an O((g + 1) log(g + 1)2k/n) upper bound for the k-th eigenvalue [17], using a multicommodity flow problem to suitably uniformize the graph metric. The study of the spectrum of a finite graph is in many ways related to the spectral theory of Riemannian manifolds, and results in geometric analysis have been a source of inspiration to state and prove corresponding results concerning finite graphs. In particular, eigenvalue bounds for surfaces have a somewhat parallel history. Hersch [13] first proved an O(1/vol(M)) 2 bound for the Neumann value of the sphere S equipped with a Riemannian metric. Yang and Yau [28] then showed that for genus g surfaces an O((g +1)/vol(M)) bound holds, and Li and Yau improved the latter result by replacing the genus with the finer conformal invariant they defined [21]. It is interesting to notice that these proofs are quite similar at a high level to 1 2 OMID AMINI AND DAVID COHEN-STEINER the ones later used in the graph setting. Conformal uniformization was used in place of circle packing representations, but the very same topological argument for centering the packing in the discrete case was used in the manifold case as well. For higher eigenvalues, Korevaar [19] established an O((g + 1)k/vol(M)) for genus g surfaces, and Hassannezhad [15] improved this to O((g + k)/vol(M)), by combining the two methods of constructing disjoint capacitors of Grigor’yan, Netrusov and Yau [14], and Colbois and Maerten [8]. While traditionally bounds on graph eigenvalues are used to prove bounds for Riemannian manifolds [3, 4, 5, 6, 9, 10], it is intriguing to see that the spectral theory of Riemannian manifolds has not been much used so far to provide information on the spectral properties of general finite graphs. Our aim in this paper is to show how eigenvalue bounds for surfaces combined with ba- sic spectral theory of (singular) surfaces, and a suitable transfer principle, allows to obtain eigenvalue estimates for graphs in terms of their geometric genus. In this way, we are able to extend the above mentioned result of [15, 19, 28] for surfaces to the graph setting using a suitable variant of the Burger-Brooks transfer method, c.f. Theorem 1.2. Our results are tight and improve the recent results of Kelner, Lee, Price and Teng [17] on bounded genus graphs. In addition to providing a uniform arguably more conceptual proof of the results of [17, 16, 27], we hope that our method makes the above mentioned existing similarities between the methods used in the spectral theory of surfaces and graphs more transparent. The transfer principle proved in this paper may be of independent interest. In fact, we shall show it can be used to provide uniform upper and lower bounds on the eigenvalues of metric graphs in terms of the eigenvalues of their simple graph models. Furthermore, it allows to generalize to completely arbitrary meshes the mesh partitioning results of Miller, Teng, Thurston, and Vavasis [24] and Spielman and Teng [27].

1.1. Statement of the main theorem on eigenvalues of bounded genus graphs. Let G = (V,E) be a finite simple graph, that we assume connected all through the paper. For two vertices u, v ∈ V , we write u ∼ v if the two vertices u and v are connected by an edge G in G. The valence of a vertex v of G is denoted by dv , or simply dv if there is no risk of confusion and the graph G is understood from the context. We denote by dmax the maximum degree of vertices of the graph, and by n the of vertices. The geometric genus of G is by definition the minimum integer g such that G can be embedded with no crossing on the compact orientable surface of genus g. Denote by C(G) the vector space of all real valued functions f defined on the set of vertices of G. The (discrete) Laplacian ∆ and the normalized Laplacian L of G are defined as follows: the Laplacian ∆ : C(G) → C(G) is the linear operator which sends a function f ∈ C(G) to ∆(f) ∈ C(G) defined by X ∀ v ∈ V (G), ∆(f)(v) = f(v) − f(u). u:u∼v Let S be the linear operator on C(G) whose matrix in the standard basis of C(G) is diagonal with entries the valences of the vertices of G, i.e., for any f ∈ C(G)

∀ v ∈ V (G),S(f)(v) = dvf(v). The normalized Laplacian is the operator S−1/2∆S−1/2. We denote by λ0(G) = 0 < λ1(G) ≤ λ2(G) ≤ · · · ≤ λn−1(G) A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 3 the set of eigenvalues of ∆, which we call the standard spectrum of G, and by nr nr nr λ0 (G) = 0 < λ1 (G) ≤ · · · ≤ λn−1(G) the set of all eigenvalues of the normalized Laplacian L, which we call the normalized spec- trum. The standard and normalized spectrum of G are easily seen to satisfy the inequalities nr nr dmin λk (G) ≤ λk(G) ≤ dmax λk (G) for any k. In this paper we prove the following theorem. Theorem 1.1. There exists a universal constant C such that the eigenvalues of the normalized Laplacian of any graph G on n vertices satisfy: d (g + k) ∀k ∈ , λnr(G) ≤ C max , N k n

where dmax and g are the maximum valence and the geometric genus of G, respectively. The linear dependance in the maximum degree is clearly optimal, as can be seen by con- sidering star graphs, which have lower bounded Fiedler value. The above result also implies a similar bound for the eigenvalues of the standard Laplacian, at the expense of an extra dmax factor. We note that Kelner, Lee, Price and Teng [17] give a similar bound for the standard spectrum with a linear rather than quadratic dependence in dmax. However, their bound has a gk log(g + 1)2 dependence instead of our (g + k) dependence. In addition to simplifying and improving the result of [17] for bounded genus graphs, we note that the dependence in g and k of our estimate is tight, at least when g is sufficiently high, see Remark 2.6. Informally, the improvement over [17] means that the asymptotic behavior of graphs’ eigen- values do not depend on the (geometric) genus of the graph. This fact, which may be seen as a one-sided discrete form of Weyl’s law for surfaces, is consistent with the intuition that at a small scale, bounded genus graphs behave like planar graphs. Finally, we note that the result in [17] also applies to graphs in any fixed proper minor-closed family (where the genus g is replaced with a parameter h depending on the family), while the stronger bounds of Theorem 1.1 cannot be extended to minor-closed classes, as we show by explicit examples in Remark 2.7.

1.2. Two-fold covers and their associated discrete Laplacians. Let M be a measured topological space, and denote by µ the measure on M.A 2-fold cover of M is a finite collection U = (Uv)v∈V , for a finite index set V , of open subsets Uv of non-zero measure such that almost every point in M is covered by exactly two subsets. To any 2-fold cover of a measured space we associate a discrete Laplacian as follows: We first form a graph G = (V,E) on the set of vertices V and with edges {u, v} ∈ E for two vertices u, v such that µ(Uv ∩ Uu) 6= 0. We define a weight function ω : E → R which to ω any edge e = {u, v} of G, associates the weight ω(e) = µ(Uu ∩ Uv). The weighted valence dv of a vertex v of G is defined by

ω X dv = µ(Uu ∩ Uv). u:u∼v

The discrete Laplacian associated to the 2-fold cover U denoted by LU is the normalized graph Laplacian associated to the weighted graph (G, ω). This is defined from the weighted Laplacian by normalizing using the weighted valence (as in the previous section). Formally, 4 OMID AMINI AND DAVID COHEN-STEINER define the weighted Laplacian ∆U : C(G) → C(G) by X   ∀v ∈ V, ∆U (f)(v) = f(v) − f(u) w({u, v}), u:u∼v ω for any f ∈ C(G). Let SU be the diagonal operator with entries the weighted valence dv of vertices v ∈ V , i.e., for any f ∈ C(G), ω ∀ v ∈ V,SU (f)(v) = dv f(v). 1 1 − 2 − 2 Then we let LU := SU ∆U SU . Denote by λk(LU ) the k-th smallest eigenvalue of LU . When (M, µ) carries a natural notion of Laplacian, it is possible to relate the eigenvalues of the Laplacian on M to the eigenvalues of LU for any 2-fold cover U. More precisely, let the measured space (M, µ) belong to any of the following three classes:

(C 1) a smooth manifold with a smooth Riemannian metric g, and µ the measure associated to the metric g; (C 2) a compact smooth surface with a conformal class of smooth Riemannian metrics g, and µ a Radon measure absolutely continuous with respect to µg, c.f. Section 2.1; (C 3) a metric graph with µ the Lebesgue measure. In any of the above cases, we can define a Laplacian on (M, µ) c.f. Section 2.1 and Section 3, and we denote by λk(M, µ), or simply λk(M) if there is no risk of confusion, the eigenvalues for the corresponding Laplacian. Our transfer principle is stated as follows. Theorem 1.2. Let (M, µ) be a measured space as in (C 1), (C 2), or (C 3) above. Assume all the elements in a 2-fold cover U of M have Neumann value at least η. Then for all positive integers k we have: λ (M) λ (L ) ≤ 2 k . k U η The main difference with the classical versions of the transfer principle [5, 3, 23] is that we discretize the continuous Laplacian as a weighted normalized graph Laplacian instead of a combinatorial one, which allows for a closer connection between the two. Our variant here uses a different notion of graph approximation that involves particular weights. In addition, the above mentioned results take as input a partition of M, while our theorem is expressed in terms of two-fold covers, which adds more flexibility. In order to prove Theorem 1.1, we apply the above theorem in the case where (M, µ) is a measured surface equipped with a conformal class of smooth Riemannian metrics g. This version seems to be required to get our Theorem 1.1 on eigenvalues of bounded (geometric) genus graphs.

1.3. Organization of the paper. The necessarily background on Laplacian eigenvalues in measured surfaces is recalled in Section 2. The proof of Theorem 1.2 for measured surfaces (Case (C 2) among the above three cases) is given in Section 2. The proof in the two other cases is similar and is thus omitted. Section 2 contains also the proof of Theorem 1.1. In Section 3, we apply Theorem 1.2 in the case of metric graphs with Lebesgue measure (C 3)), to obtain a uniform quantitative complement to a theorem of X. Faber [11] on the spectral convergence of finite graphs to metric graphs. Moreover, we give in Section 4 an algorithmic A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 5 application of the above theorem to mesh partitioning in numerical analysis, generalizing the results of Miller-Teng-Thurston-Vavasis [24] and Spielman-Teng [27] to anisotropic meshes.

2. Eigenvalues of bounded genus graphs In this section we give the proofs of Theorem 1.2 and Theorem 1.1. We start by recall- ing the variational approach to study eigenvalue problems for surfaces with measures [18], which provides a setting to study eigenvalue problems for singular surfaces. This makes the statement of Theorem 1.2 precise in the case of a measured metric surface. 2.1. Eigenvalues on measured surfaces. Let M be a smooth compact surface, possibly with boundary, which we suppose equipped with a smooth Riemannian metric g. Denote by µg the induced volume form on M. Let µ be a Radon measure on M which we suppose absolutely ∞ 2 continuous with respect to the measure µg. For a C -smooth function f ∈ L (M, µ), the

Rayleigh quotient RMg (f, µ) is defined by R 2 M |∇g f| dµg RMg (f, µ) := R 2 . M f dµ The eigenvalues of the measured metric surface (Mg, µ) are defined by the variational formula:

(1) λk(Mg, µ) := inf sup RMg (f, µ), Λk+1 ∗ f∈Λk+1 2 where Λk+1 ⊂ L (M, µ) varies over subspaces of dimension k + 1 which consist only of ∗ smooth functions on M, and Λk+1 = Λk+1 \{0}. Note that in the case µ = µg, we recover the usual variational characterization of the eigenvalues of the Laplacian ∆g associated to the Riemannian surface Mg. To see the point of introducing this formalism, assume that the two metrics g and h on M are conformally equivalent. From the conformal invariance of the Dirichlet integral, we see that RMg (f, µ) = RMh (f, µ). In particular, letting µ = µh, we see that the spectra of the metric h within the conformal class of Mg coincides with the spectra of measured surface (Mg, µ) for an appropriate Radon measure µ. Now, if h is a metric with conical singularities, it is a classical fact that M is conformally equivalent to a constant curvature metric g, the conformal factor being square integrable with respect to the corresponding area form µg. Thus the framework of measured metric surfaces allows in particular to define spectra of surfaces with conical singularities as the spectra of the measured metric surface (Mg, µh) [18]. Let U be an open subset of M, and denote by U the topological closure of U in M. The Neu- R 2 R 2 mann value λ(U) of U is defined as the infimum of the Rayleigh ratio |∇gf| dµg/ f dµ U R U over all smooth functions f on U which extend continuously to U and satisfy U fdµ = 0, R 2 ||∇gf|| dµg λ(U) := inf U . R R 2 f: U fdµ=0 U f dµ

2.2. Proof of Theorem 1.2. We suppose M, µ, g, and µg as above, and consider a 2-fold cover U = (Uv)v∈V of M. Denote by η the minimum of λ(Uv) for v ∈ V . Let G be the associated weighted graph with vertex set V and weight matrix WU = [ω({u, v})]u,v, where ω({u, v}) = µ(Uu ∩ Uv) for u 6= v. Let LU be the matrix of the associated normalized graph −1/2 −1/2 Laplacian. We have LU = I − SU WU SU , where the matrix SU is diagonal with entries ω P given by the weighted valences of the vertices dv = u:u∼v ω({u, v}). 6 OMID AMINI AND DAVID COHEN-STEINER

Proof of Theorem 1.2. Let v ∈ V and f any smooth function on M. By restricting f to Uv and substracting the mean over Uv, we get: Z Z Z 2 1 2 ||∇gf|| dµg ≥ λ(Uv) f − fdµ dµ Uv Uv µ(Uv) Uv  Z 1 Z  ≥ η f 2 dµ − f dµ2 , Uv µ(Uv) Uv

Summing the last inequalities over v ∈ V yields: Z Z 2 X 1 2 (2) ||∇ f||2dµ ≥ 2||f||2 − fdµ , η g g 2 µ(U ) M v v Uv 2 where the L norm ||.||2 is with respect to the measure µ. Denote by 1Uv the −1/2 2 function of the open set Uv, and let φv = µ(Uv) 1Uv . Define Φ : L (M) → C(G), by Z Φ(f)(v) := fφv, M on any vertex v of G. We see that the quadratic form in f in the right hand side of Equation (2) 2 2 is given by 2||f||2 − ||Φf||2.  Let  > 0, and denote by Λk+1 a (k + 1)-dimensional space of smooth functions on M such that for any f ∈ Λk+1 \{0}, we have R 2 M ||∇gf|| dµg R 2 ≤ (1 + )λk(M). M f dµ Note that by the variational characterization of the eigenvalues (see e.g. (1) and (6)), such a  space exists. For any f ∈ Λk+1, by inequality (2), we have: Z 2(1 + )λk(M) 2 2 2 2 2 ||f||2 ≥ ||∇gf|| dµg ≥ 2||f||2 − ||Φf||2. η η M That is:  (1 + )λ (M) ||Φf||2 ≥ 2 1 − k ||f||2. 2 η 2 Let Φ∗ denote the adjoint of the operator Φ : L2(M) → C(G). From the variational char- acterization of the eigenvalues, this implies that the compact self-adjoint operator Φ∗Φ on 2 L (M) has at least k + 1 eigenvalues greater than or equal to 2(1 − (1 + )λk(M)/η). We can assume that this latter quantity is positive, otherwise there is nothing to prove since all the eigenvalues of the normalized Laplacian are at most 2. Since the non zero eigenvalues of Φ∗Φ are the same as the non zero eigenvalues of ΦΦ∗, we thus deduce that 2(1 + )λ (M) (3) λ (2I − ΦΦ∗) ≤ k . k η To conclude the proof, it suffices to notice that Z µ(U ∩ U ) [ΦΦ∗] = φ φ = u v . u,v u v 1/2 M (µ(Uu)µ(Uv)) A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 7

Because U is a 2-fold cover, the v-th entry of the diagonal matrix DU is equal to µ(Uv). We ∗ thus easily check that 2I − ΦΦ = LU . Therefore, inequality (3) gives 2(1 + )λ (M) λ (L ) ≤ k . k U η Since this holds for any  > 0, the theorem follows. 

2.3. Proof of Theorem 1.1. We first give some background about graphs embedded in a surface, and refer to [25] for more details. We assume that all surfaces are compact, orientable and without boundary. An embedding of a graph G in a surface M is a drawing of G on M so that all vertices of G are distinct on M, and every edge of G form a simple arc on M connecting its two endpoint vertices. Interior of edges and vertices are assumed to be pairwise disjoint. A face of and embedding, or simply a face of G if the embedding is clear from the context, is a connected component of the complementary of G in M. 2 An embedding is called cellular if every face is homeomorphic to an open disk in R . The genus g(G) is the minimum integer g such that G has an embedding in a surface M of genus g. The following result will allow us to suppose that a graph G with a given genus g(G) is embedded in a cellular way. Proposition 2.1 ([25, Proposition 3.4.1]). Every embedding of a connected graph G in a surface of genus g(G) is cellular. Suppose from now on that the connected graph G is embedded in a cellular way in a surface 2 M of genus g, so that every face F is homeomorphic to an open disk DF in R . The boundary 1 of the face F in M is the image of the boundary ∂DF ' S under a continuous map, which is locally a homeomorphism away from the preimage of the vertices. We denote by F the set of all faces of G. For any face F ∈ F, we define a boundary walk of F to be any walk in the graph G consisting of vertices and edges as they are encountered when walking along the whole boundary of F , following the circle ∂DF , and starting at some vertex. Note that some edges may appear more than once in a boundary walk. The degree of a face F ∈ F, denoted deg(F ), is the number of edges on any boundary walk of F .

We define a new multigraph G = (V , E ) embedded in M, and containing G as an induced subgraph, by coning over boundary walks of faces as follows. The vertex set V of G consists  of the vertices in G and a new vertex vF for each face F of F, i.e., V = V t vF F ∈F . For each face F of F, let v1, . . . , vdeg(F ) be the vertices of G which appear in this order in a boundary walk of F . Note that a vertex might appear more than once. The edge set E of G consists of the edges in E, and new edges {vF , vi}, for i = 1,..., deg(F ). The embedding of G in M is obtained in the following natural way: each face F ∈ F is homeomorphic to a 2 disk DF in R , and under this homeomorphism, the vertices v1, . . . , vdeg(F ) in the boundary walk of F , appear in this cyclic order on the boundary of DF . Choose the center of DF as the image of vF and the rays from vF to vi as the image of the edges {vF , vi}. We refer to all the new edges {vF , vi} added in the process as cone edges of G . Call an embedding of a graph in M a weak triangulation if the degree of any face of the embedding is three. We use this terminology since it can happen that two different faces of the embedding share more than one edge, in which case we do not have a triangulation. The embedded (multi)graph G constructed above has the following properties. 8 OMID AMINI AND DAVID COHEN-STEINER

Proposition 2.2. The embedding of G in M is a weak triangulation, and each face of G is G G incident to exactly one edge of G. Moreover, for any vertex v of G, we have dv = 2dv , where G G dv and dv denote the valence of v in G and G, respectively. Proof. By definition of the embedding, each face of G consists of two cone edges and an edge of G, which proves the first assertion. To prove the second statement, let F1,...,Fk ∈ F be all the faces of G which are incident to the vertex v ∈ V . For each i = 1, . . . , k, the number of edges {vFi , v} in E is half the number of edges of G in a boundary walk of Fi. Each edge e ∈ E incident to v appears precisely twice in the union of the edges of the boundary walks of F1,...,Fk. This shows that the total number of edges of E of the form {vFi , v}, for G i = 1, . . . , k, is equal to dv , which proves the claim.  Definition 2.3 (Open star). For each vertex v of G, we define the open star of v in M with respect to the embedding of G , or simply the open star of v, denoted by Sv, as the interior of the union of all the faces of G which contain v in their boundaries. Let now (M, g, µ) be any surface as in Section 2.1, so the measure µ is absolutely continuous with respect to the volume form µg of the smooth Riemannian metric g. Proposition 2.4. The open stars of vertices of G form a 2-fold cover of M. Proof. By Proposition 2.2, the boundary walk of each face of G is a triangle which has exactly two vertices in G. It follows that every point of M \ G appears in exactly two open stars, which proves the claim by absolute continuity of µ with respect to µg. 

We now introduce a metric h on M with conical singularities (and will later assume µ = µh). For reasons that will soon become clearer, we set the length of each edge of E ⊂ E to be −1 equal to one, and the length of each cone edge in E to be cos(π/(2dmax)) /2. We equip M with the natural metric h such the triangles have the Euclidean metric induced by their edge lengths. Note that for any triangle T of G , the angle of T at any of its vertices that belongs to the graph G is equal to π/(2dmax). The metric h has only conical singularities, and we denote by λk(M) = λk(Mh), the eigenvalues of the surface M as defined in Section 2.1. Thus, λk(M) is the k-th eigenvalue of the measured metric surface (Mg, µ = µh), where g is a metric of constant curvature in the conformal class of h. Using Theorem 1.2 we can relate the eigenvalues of M to those of G. Denote by U the 2-fold cover of M given by the open stars Sv of vertices of V , c.f. Proposition 2.4. The intersection of two distinct open stars Su and Sv, for two vertices u and v of G, has non zero measure if and only if u and v are neighbors in G. Moreover, all non-empty intersections have the same measure, equal to the area of two triangles in M. Therefore, the normalized Laplacian LU equals the normalized Laplacian of G. Hence, in order to apply the transfer result Theorem 1.2, we only need to lower bound the Neumann value of the open stars Sv of the vertices of G. We do so by again applying the transfer result to a specific 2-fold cover of each open star Sv, for v ∈ V . Thanks to the choice of edge lengths, the vertices of G in (M, h) have non negative curvature. It follows that by cutting Sv along an arbitrary cone edge, we can unfold Sv to 0 the plane without overlap. Denote by Sv the unfolded star of v as isometrically embedded in 2 R . The cutting operation can only decrease the Neumann value so it is sufficient to bound 0 2 from below the Neumann value of the unfolded open subset Sv of R . 0 0 We call a kite in Sv the union of two triangles in Sv which share an edge of G. So for any edge e = {v, u} ∈ E, there is a kite Ke, and the union of the kites Ke for e incident to v is A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 9

0 equal to the planar set Sv. For any edge e ∈ E incident to v, the kite Ke has two diagonals op 0 composed of the edge e and the diagonal opposite to e, that we denote by diage . Cut Sv op 0 along all the opposite diagonals diage for e ∈ E incident to v. This cuts Sv into an open region Pv with polygonal boundary containing the vertex v, together with one triangle Te for each kite Ke for e ∈ E incident to v. Define    Uv := Ke e={v,u}∈E ∪ Te e={v,u}∈E ∪ Pv . 0 The cover Uv of Sv is a 2-fold cover.

Proposition 2.5. Any X ∈ Uv is a convex set of diameter at most two.

Proof. The triangle Te, for e = {u, v} ∈ E, is obviously convex of diameter one, and so is the kite Ke. As for the region Pv, to prove the convexity of P , it will be enough to show that the angle of P at v is at most π. As we previously observed, by the choice of the edge lengths, all the triangles of G has angle π/2dmax at any of their vertices which belong to G. The number 0 G G of triangles of Sv is at most 2dv , since dv ≤ dmax, it follows that the angle at v of P is at most π, and the convexity follows. The claim on the diameter follows from the fact that all the edges of G have length equal to one. 

Since any element of Uv is a planar convex sets of diameter at most two, there exists a universal constant C1 > 0 such that for any X ∈ Uv, we have λ(X) ≥ C1 for λ(X) the 0 Neumann value of X [26]. For the 2-fold cover of Sv, the non-zero element µh(X ∩ Y ) for X 6= Y ∈ Uv have the same value, equal to the area of a triangle in M. Therefore, the normalized Laplacian LUv equals the normalized Laplacian of the graph whose edges are the 1 pairs X 6= Y ∈ Uv whose intersection has positive measure. This graph Sd is obtained from the star graph with d edges by inserting a new vertex in the middle of each edge. (Recall that a star graph with d edges has a central vertex connected to d other vertices.) The 1 Neumann value of Sd is lower bounded by an absolute constant C2 > 0 independent of d. Hence, applying the transfer theorem 1.2 to the 2-fold cover Uv, it follows that there exists a 0 universal constant C3 = C1.C2/2 > 0 such that the Neumann value of Sv is bounded from 0 below by C3, i.e., λ(Sv) ≥ C3. This gives λ(Sv) ≥ C3. We get from these observations, and Theorem 1.2 applied to the 2-fold cover of U, that

nr λk(M) (4) λk = λk(LU ) ≤ 2 . C3 A result of Hassannezhad [15] states now that there is a universal constant A such that for each k:

(5) λk(M)µ(M) ≤ A(g + k). Note that this result is not explicitly stated in the framework of measured metric surfaces in [15], however the proof given in [15] works also in this setting and gives the above statement. Putting Equations (4) and (5) together, and observing that µ(M) ≥ C4n/dmax for some constant C4, we conclude that for C = 2A/(C3C4), we have d (g + k) λnr(G) ≤ C max , n which is the statement of Theorem 1.1.  10 OMID AMINI AND DAVID COHEN-STEINER

Remark 2.6. It is shown in [7] that for large g, there are area one and genus g Riemannian surfaces S with 4π λ (S) ≥ (g − 1) + 8π(k − 1) −  k 5 for any  > 0. Now, the classical Brooks-Burger method implies the existence of a bounded degree genus g graph G with n vertices such that λk(G) ≥ Cλk(S)/n. Hence, at least for large enough n and g, there are graphs whose eigenvalues match the behaviour of the estimate in Theorem 1.1. Remark 2.7. The following example shows that the strong estimates as in Theorem 1.1 cannot hold for more general classes of graphs closed under taking minor. Recall that the Cartesian product G1G2 of two graphs G1 = (V1,E1) and G2 = (V2,E2) has vertex set V1 × V2 and there is an edge between (v1, v2) and (u1, u2) in V1 × V2 if either u1 = v1 and {u2, v2} ∈ E2, or u2 = v2 and {u1, v1} ∈ E1. The Laplacian eigenvalues of G1G2 are of the form λi(G1) + λj(G2) for i = 1,..., |V1| and j = 1,..., |V2|. Let d be a fixed large enough integer, and for any ` ∈ N, consider the Cartesian product C2lG of a cycle C2` of length 2l with a d-regular graph G on t vertices, for an integer t ∈ N. For any fixed t ∈ N, we get in this way a family of graphs by varying ` and G. All these graph are of treewidth bounded by some f(t) for a (linear) function f of t. Bounded treewidth graphs form a minor-closed family, so all these graphs belong to a fixed proper minor-closed th family Ft. For G a random d-regular graph on t vertices, and for the l eigenvalue of tl C2l G ∈ Ft, for l ∈ , we have λl(C2` G) = Ω( ) with high probability as t tends to  N  |C2lG| infinity. This shows that there do not exist in general constants h = h(Ft) and C = C(Ft) 2 associated to Ft ensuring that the inequality λk(G) ≤ C dmax (gt + k)/n hold for any graph G ∈ Ft on n vertices, and for any k ∈ N (unlike what happens for the class of bounded genus graphs). In particular, the strong estimates as in Theorem 1.1 cannot hold for general minor-closed classes of graphs.

3. Eigenvalues of the Laplacian on metric graphs We briefly review the basic definitions concerning the spectral theory of metric graphs, and refer e.g. to [2, 29] for more details. Let G = (V,E) be a finite connected graph and let ` : E → R>0 be a (length) function on the edges of G. The length of e is denoted by `e. We define the metric realization of (G, `) as follows: for each edge e = uv of G take a closed interval Ie ⊂ R of length `e, and a surjection πe : ∂Ie → {u, v} (which identifies the two extremities of Ie with the vertices of G in e). Define the topological space (with the quotient ) G  Γ := (V t Ie)/ x = πe(x) ∀e ∈ E & x ∈ ∂Ie . e The space Γ has a natural metric, the shortest path metric induced by piecewise isometric paths between points, see e.g. [2]. We call a metric graph any Γ isometric to a metric realization of a pair (G, `), as above. The pair (G, `) is called a model of Γ; when G is a simple graph, the model is called simple. Note that there are plenty of models for a metric graph Γ, e.g. any finite subset of points of Γ can be part of a simple model of Γ. 1 For any point p ∈ Γ, we denote by Tp Γ the set of unit tangent vectors to Γ at p. For an 1 ~ ~ interval I = [a, b], in R, we define Ta I = {1}, with 1 the unit vector in R. For a metric graph Γ and a point p ∈ Γ, let (G, `) be a simple model of Γ with p ∈ V (G), and let e1, . . . , ed be the A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 11

1 edges of G incident to v. Define Tp Γ as the set of all unit tangent vectors at p of the intervals 1 Iej , as above. Let ~u ∈ Tp Γ be a unit tangent vector, and let I = Ie be the corresponding interval (corresponding to the edge e of a simple graph model (G, `)). For  > 0 sufficiently small, we denote by p + ~u the unique point in I at distance  from p on I. A function f :Γ → R is piecewise smooth if there exists a simple graph model G = (V,E) of Γ such 2 that the restriction of f to the intervals Ie, for e ∈ E, are of class C . The space of piecewise smooth function on Γ is denoted by S(Γ). Let f :Γ → R be a piecewise smooth function on 1 a metric graph Γ. Let p ∈ Γ and ~u ∈ Tx Γ a unit tangent vector to at x. The (outgoing) slope of f along ~u denoted by d~u(f) is defined by f(p + ~u) − f(p) d~u(f) := lim . →0+ 

For a point p ∈ Γ, we define σp as the sum of the slopes of f along unit tangents: X σp := d~uf(p), 1 ~u3Tp Γ

Note that for all but at most a finite number of points p ∈ Γ, we have σp = 0. A metric graph Γ has a natural Lebesgue measure denoted by dx. The Laplacian of Γ is the (measure valued) operator ∆ on Γ which to a function f ∈ S(Γ) associates the measure 00 X ∆(f) := −f dx − σpδp. p∈Γ Define the Zhang space Zh(Γ) as the space of all functions f ∈ S(Γ) such that f 00 ∈ 1 L (Γ, dx). The inner product ( , ) and the Dirichlet pairing ( , )Dir on Zh(Γ) are defined by Z ∀f, g ∈ Zh(Γ), (f, g) := fg dx, and Γ Z Z Z 0 0 0 0 (f, g)Dir := f∆(g) = g∆(f) = f g dx = (f , g ). Γ Γ Γ A function f in Zh(Γ) is an eigenfunction of the Laplacian on Γ with eigenvalue λ if for any function g ∈ Zh(Γ), we have (f, g)Dir = λ(f, g). The eigenvalues of ∆ are all nonnegative and, assuming Γ is connected, they form a discrete subset 0 = λ0(Γ) < λ1(Γ) < λ2(Γ) < ··· < λn(Γ) < . . . of R. In addition, λk(Γ) has the following (usual) variational characterization: (f, f)Dir (6) λk(Γ) = inf sup Λk+1⊂Zh(Γ) f∈Λk+1 (f, f) dim(Λk+1)=k+1 Definition 3.1 (Dilation of a metric graph). Let Γ be a metric graph with a simple graph model (G, `), and β ∈ R>0. The metric graph βΓ is defined as the metric realization of the pair (G, β`). The following proposition is straightforward, see e.g. [2]. Proposition 3.2. Let Γ be a metric graph and β > 0 a real. For any integer k ≥ 0, we have 1 λk(βΓ) = β2 λk(Γ).

By a metric star S we mean the metric realization of a pair (Sd, `) with Sd = K1,d a star graph of arbitrary valence d, and ` a length function on E(Sd). For such a metric star, define

`max(S ) := maxe∈E(Sd) `(e). 12 OMID AMINI AND DAVID COHEN-STEINER

Lemma 3.3. For any metric star S , we have: π2 λ1(S ) ≥ 2 . 4`max(S )

Proof. Assume that S is the metric realization of a pair (Sd, `) with d ∈ N. We adapt the argument in [12, Example 3] to the case where the branches of S have non-necessary equal lengths. Let us parametrize each each edge e of Sd with the interval [0, `e] starting from the leaf vertex towards the central vertex of Sd. In this parametrization, an eigenfunction√ φ of the Laplacian, with corresponding eigenvalue λ, must be of the form ae cos( λxe), where xe is the length parameter of the edge e in S , for e ∈ E(Sd). This follows in particular from the fact that the slope of an eigenfunction must be zero at leaves.√ Now let a be the value of the eigenfunction φ at the center of S . If a = 0, we get that λ`e ∈ π/2 + N, for any edge e, which implies the claim. If a is non zero, then we use the fact that at the center of S , the sum of the (out-going) slopes of φ along the branches must be zero [2, Proposition 15.1], which gives X √ ae sin( λ`e) = 0.

e∈E(Sd) √ Since ae cos( λ`e) = a for any edge e of Sd, this implies X √ tan( λ`e) = 0,

e∈E(Sd) and so, again, at least one of the arguments in the tangents must be at least π/2, and the lemma follows. 

For a simple graph G and a vertex v ∈ V (G), we denote by ΣG(v) the star subgraph of G with central vertex v and with the edge set all the incident edges to v. Let (G, `) be a simple graph model of a metric graph Γ. For any v ∈ V (G), we define the metric star with center v (with respect to G) of Γ denoted by SG(v), or simply Sv if there is no risk of confusion, as the subset of Γ isometric to the metric realization of ΣG(v) with length function given by `. Denote by `max,G the maximum length of edges in G, and note that `max,G = maxv∈V (G){`max(Sv)}. Given a simple graph model (G, `) of a metric graph Γ, the family of all the metric stars nr Sv, for v ∈ V (G), forms a 2-fold cover S of Γ. Denote by λk (G, `) the k-th eigenvalue of LS . Lemma 3.3 together with Theorem 1.2 yields the following bound:

Theorem 3.4. Let Γ be a metric graph with a simple graph model (G, `). For any k ∈ N, we have 2 π nr λk(Γ) ≥ 2 λk (G, `). 8`max,G We now show that under certain natural conditions, it is possible to achieve eigenvalue upper bounds closely matching the lower bounds of the above corollary. For a simple graph model (G, `) of Γ denote by `min,G the minimum length of edges e in E(G). Definition 3.5. A simple graph model of a metric graph Γ is called length-balanced if for any edge e ∈ E(G), we have `e ≤ 2`min,G. We have the following theorem. A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 13

Theorem 3.6. There are absolute constants c1, c2 > 0 such that for any length-balanced simple graph model (G, `) of Γ on n vertices, and for any non-negative integer k ≤ n − 1, we have c2 2 nr 2 `min,G λk(Γ) ≤ λk (G, `) ≤ c1`min,G λk(Γ). dmax Before giving the proof, we state an interesting corollary of the above theorem. We first need the following definition.

Definition 3.7. Let Γ be a metric graph. Define `min as the supremum of `min,G over all length-balanced simple graph models (G, `) of Γ. It is easy to see that there is a length-balanced simple graph model G of Γ such that `min = `min,G. For such a simple graph model (G, `), define the model (Gk, `) as the k-th subdivision of G where each edge e is subdivided into k edges of equal lengths `e/k. Note that Gk is length-balanced, has at least k + 1 vertices, and has minimum edge length equal to `min/k. Thus as a consequence of Theorem 3.6, we get

Corollary 3.8. With the notations as above, there are absolute constants c1 and c2 such that for any metric graph Γ, we have

c2 2 2 nr 2 `minλk(Γ) ≤ k λk (Gk, `) ≤ c1`minλk(Γ). dmax Our results, especially corollary 3.8 above, should be viewed as a quantitative complement to a theorem of X. Faber [11] on the spectral convergence of finite graphs to metric graphs, in the sense that they provide uniform upper and lower bounds on the eigenvalues of Γ in terms of eigenvalues of simple graph models of Γ.

1 nr nr Proof of Theorem 3.6. First note that since λk(βΓ) = β2 λk(Γ) and since λk (G, β`) = λk (G, `), by the very definition, it will be enough to prove the theorem for `min = 1. The right hand side inequality follows from Theorem 3.4, and the well-balanced property of the simple graph model G of Γ. We now prove the other inequality, namely the existence nr of c2 such that for any k ≤ n − 1, c2 λk(Γ) ≤ dmaxλk (G, `) (still under the assumption that `min = 1 and the length-balanced property of the model (G, `)). Since the lengths of all edges are between 1 and 2, we get λnr(G, `) ≥ 1 λ (G). Indeed, letting g = D1/2f, we have the k 2dmax k S following expression for the Rayleigh quotient P 2 P 2 (g, LS g) e={u,v}∈E `(e)(f(u) − f(v)) 1 e={u,v}∈E(f(u) − f(v)) = P ` 2 ≥ P 2 (g, g) v dvf(v) 2dmax v f(v) ` P (where dv = u:u∼v `({u, v})), which using the variational characterization of the eigenvalues 0 proves the claim. So it will be enough to show the existence of a constant c2 such that 0 c2λk(Γ) ≤ λk(G).

Consider Wk+1 the vector space of dimension k+1 generated by the first k+1 eigenfunctions g0, . . . , gk ∈ C(G) associated to λi(G), for i = 0, . . . , k. Note that in particular X (g(u) − g(v))2 λ (G) ≥ k P g(v)2 u,v∈V (G) v u∼v 14 OMID AMINI AND DAVID COHEN-STEINER for any g ∈ Wk+1 \{0}. We construct an injective linear map Ψ : C(G) → Zh(Γ) such that for any g ∈ C(G) \{0}, we have 2 (Ψ(g), Ψ(g))Dir X (g(u) − g(v)) ≤ 8 . (Ψ(g), Ψ(g)) P g(v)2 u,v∈V (G) v u∼v

Applying the variational characterization of λk(Γ), given in Equation (6), to the test space Ψ(Wk), for k ≤ n − 1, will then give the result. Consider an edge e = {u, v} of G, and denote by u and v the two points at distance 1 e e 4du and 1 from u and v on e, respectively, where d and d denote the valence of the vertices 4dv u v 1 u and v in G, respectively. Note that the length of each segment [ue, ve] in Γ is at least 2 . For any vertex v of G, denote by Bv the union of all segments [v, ve] on the edges e adjacent to v in G (i.e., B is the ball of radius 1 around v in Γ). For any function g ∈ C(G), defined v 4dv on the set of vertices of G, let Ψ(g) be the function on Γ which takes value equal to g(v) on each ball Bv, and which is affine linear of slope (g(v) − g(u))/`([ue, ve]) on each segment [ue, ve], for any edge e ∈ E(G). Obviously, Ψ is an injective linear map from C(G) to Zh(Γ). Let now g ∈ C(G) \{0} and denote f = Ψ(g). We have Z 02 X 1 2 X 2 (f, f)Dir = f dx = (g(u) − g(v)) ≤ 2 (g(u) − g(v)) . `([ue, ve]) Γ e={u,v}∈E(G) {u,v}∈E(G)

Denote by B the union ∪v∈V (G)Bv. Since each ball Bv has total length equal to 1/4, we have Z Z 1 X f 2dx ≥ f 2dx = g(v)2. 4 Γ B v∈V (G) It thus follows from the two above estimates that for any g ∈ C(G) \{0}, we have 2 (Ψ(g), Ψ(g))Dir X (g(u) − g(v)) ≤ 8 , (Ψ(g), Ψ(g)) P g(v)2 u,v∈V (G) v u∼v and the theorem follows.  4. Anisotropic mesh partitioning In this final section we discuss a practical application of our transfer theorem to the mesh partitioning problem in scientific computing. Parallelizing finite elements computations re- quires to split the base mesh in such a way that communication between different pieces is minimized. This is naturally formalized as a (possibly multi-way) sparsest cut problem, which we may want to solve using spectral clustering. Guarantees for such methods in this setting were proved by Miller-Teng-Thurston-Vavasis and Spielman-Teng [24, 27]. More precisely, these papers show that spectral partitioning provides good cuts for meshes in d-dimensional Euclidean space provided that all d-simplices in the mesh are well-shaped, i.e. not too far from being equilateral. It is not hard to design a 2-fold cover of a general mesh such that our transfer result provides guarantees for spectral clustering applied to anisotropic meshes. Specifically, let T d be a triangulation of a domain D ⊂ R . Performing a barycentric subdivision of all d-simplices 0 gives a triangulation T . For a d-simplex σ of T , let now Uσ be the interior of the union of 0 σ with the d + 1 d-simplices of T that share a facet with σ. The collection of Uσ forms a A TRANSFER PRINCIPLE AND APPLICATIONS TO EIGENVALUE ESTIMATES FOR GRAPHS 15

2-fold cover U of the domain, and the corresponding Laplacian LU is defined using weights

wσ1,σ2 that are proportional to the sum of the volumes of σ1 and σ2. Hence, assuming that neighboring d-simplices in T have volumes within a ratio of κ > 1, we see that the eigenvalues of LU and those of the normalized Laplacian of the dual graph of T are also within a ratio of κ.

−1 −1 −2 Proposition 4.1. The Neumann value of Uσ is at least C κ  for some universal con- stant C > 0, where  is the maximum diameter of simplices in T .

0 Proof. Let τi, i = 1 . . . d + 1, be the d-simplices in T that share a facet with σ, and σi be the 0 d-simplex in T that is included in σ and shares a facet with τi. The interiors of σ, τi, and of τi ∪ σi form a 2-fold cover of Uσ. The entries of the corresponding Laplacian are within a factor κ of the those of the normalized Laplacian of the intersection graph of the elements of the cover, which is a once subdivided star graph. Such a star graph has Fiedler value lower bounded by a constant. Now each element in the cover is a convex set with diameter at most 2, so by [26] their Neumann value is lower bounded by a constant times −2. The claim then follows from theorem 1.2.  Therefore, Theorem 1.2 applied to the cover U yields that the Fiedler value of the dual 2 2 graph of T is at most 2Cκ λ1(D) . By Cheeger’s inequality, a suitable spectral partitioning 0p 0 algorithm gives a balanced cut of size at most κC λ1(D)/, for some constant C . We note that if d-simplices in T are nearly equilateral, then  ' (vol(D)/n)1/d, where n is the number of simplices in T . Hence in this case we recover the n1/d behaviour proved in [24, 27] for the size of the cut, since the assumption that simplices are well-shaped implies an upper bound on κ. However, the methods used in those works do not seem to apply to the case of general anisotropic meshes. Acknowledgments: This work was started during a stay of the first named author at Labora- toire J. A. Dieudonn´eat Universit´ede Nice Sophia-Antipolis, and pursued during another visit at INRIA Sophia-Antipolis. He thanks both these institutions, specially Philippe Maisonobe and Jean-Daniel Boissonnat, for their support. This research has been partially supported by the European Research Council under Ad- vanced Grant 339025 GUDHI (Geometry Understanding in High Dimensions).

References

[1] N. Alon and V.D. Milman, λ1, Isoperimetric inequalities for graphs, and superconcentrators, Journal of Combinatorial Theory, Series B 38 (1985), 73–88. [2] M. Baker and R. Rumely, Harmonic analysis on metrized graphs, Canadian Journal of Mathematics 59 (2007), 225–275. [3] R. Brooks, The spectral geometry of a tower of coverings, Journal of Differential Geometry 23 (1986), 97–107. [4] R. Brooks and E. Makover, Random construction of Riemann surfaces, J. Differential Geom. 68 (2004), 121–157. [5] M. Burger, Petites valeurs propres du Laplacien et topologie de Fell, th`esede doctorat, 1986, Econom Druck AG (Basel). [6] P. Buser, Cubic graphs and the first eigenvalue of a Riemann surface, Math. Z. 162 (1978), no. 1, 87–99. [7] B. Colbois and A. El Soufi, Extremal eigenvalues of the Laplacian in a conformal class of metrics : the “conformal spectrum”, Annals of Global Analysis and Geometry 24 (2003), 337–349. [8] B. Colbois and D. Maerten, Eigenvalue estimates for the Neumann problem of a bounded domain, J. Geom. Anal 18 (2008), 1022–1032. 16 OMID AMINI AND DAVID COHEN-STEINER

[9] B. Colbois and Y. Colin de Verdi`ere, Sur la multiplicit´ede la premi`ere valeur propre d’une surface de Riemann `acourbure constante, Commentarii Mathematici Helvetici 63 (1988), 194–208. [10] Y. Colin de Verdi`ere, Spectres de vari´et´esRiemanniennes et spectres de graphes, Proceedings of the International Congress of Mathematicians. Vol. 1. 1986. [11] X. Faber, Spectral convergence of the discrete Laplacian on models of a metrized graph, New York J. Math.12 (2006), 97–121. [12] L. Friedlander, Extremal properties of eigenvalues for a metric graph, Annales de l’Institut Fourier 55 (2005), 199–211. [13] J. Hersch, Quatre propri´et´esisop´erim´etriquesde membranes sph´eriqueshomog`enes, Comptes-rendus de l’acad´emiedes sciences s´erieA-B, 270:A1645-A1648, 1970. [14] A. Grigor’yan, Y. Netrusov, and S.T. Yau, Eigenvalues of elliptic operators and geometric applications, Surveys in Differential Geometry IX (2004), 147–217. [15] A. Hassannezhad, Conformal upper bounds for the eigenvalues of the Laplacian and Steklov problem, Journal of Functional Analysis 261 (2011),3419–3436. [16] J. Kelner, Spectral partitioning, eigenvalue bounds, and circle packings for graphs of bounded genus, SIAM Journal on Computing 35 (2006), 882-902. [17] J. Kelner, J.R. Lee, G. Price, S.H. Teng, Metric uniformization and spectral bounds for graphs, Geometric and Functional Analysis 21 (2011), 1117–1143. [18] G. Kokarev, Variational aspects of Laplace eigenvalues on Riemannian surfaces, Advances in Mathematics 258 (2014), 191–239. [19] N. Korevaar, Upper bounds for eigenvalues of conformal metrics, Journal of Differential Geometry 37 (1993), 73–93. [20] J.R. Lee, S.O. Gharan, L. Trevisan, Multi-way spectral partitioning and higher-order Cheeger inequalities, In Proceedings of the forty-fourth annual ACM Symposium on Theory of Computing (STOC ’12). ACM, New York, NY, USA, 1117-1130. [21] P. Li, S.T. Yau, A new conformal invariant and its application to the Willmore conjecture and the first eigenvalue of compact surfaces, Inventiones mathematicae 69 (1982), 269–291. [22] A. Louis, P. Raghavendra, P. Tetali, S. Vempala, Many sparse cuts via higher eigenvalues, In Proceedings of the forty-fourth annual ACM Symposium on Theory of Computing (STOC ’12). ACM, New York, NY, USA, 1131-1140. [23] T. Mantuano, Discretization of compact Riemannian manifolds applied to the spectrum of the Laplacian, Annals of Global Analysis and Geometry 27 (2005), 33–46. [24] G.L. Miller, S.H. Teng, W. Thurston, S.A. Vavasis, Geometric separators for finite element meshes, SIAM J. Scientific Computing 19 (1998), 364–386. [25] B. Mohar, C. Thomassen, Graphs on surfaces, Johns Hopkins University Press, 2001. [26] L. E. Payne, H.F. Weinberger, An optimal Poincar´einequality for convex domains, Arch. Rational. Mech. Anal. 5 (1960), 286–292. [27] D. Spielman, S.H. Teng, Spectral partitioning works: Planar graphs and finite element meshes, Linear Algebra and its Applications 421 (2007), 284–305. [28] P. Yang, S.T. Yau, Eigenvalues of the Laplacian of compact Riemann surfaces and minimal submanifolds, Ann. Scuola Norm. Sup. Pisa Cl. Sci. 4 (1980), 55–63. [29] S. Zhang, Admissible pairing on a curve. Inventiones mathematicae 112 (1993), 171–193.

CNRS - Departement´ de mathematiques´ et applications, Ecole´ Normale Superieure,´ Paris E-mail address: [email protected]

INRIA, 2004 Route des Lucioles, BP93, Sophia-Antipolis, France E-mail address: [email protected]