The Discrete Harmonic Cohomology Module on Networks
Will Dana, David Jekel, Collin Litterell, Austin Stromme November 29, 2015
Contents
1 Graphs and Harmonic Functions 2 1.1 Graphs with Boundary ...... 2 1.2 Harmonic Functions ...... 5 1.3 The Dirichlet Problem and The Kirchhoff Matrix ...... 5 1.4 U as a Functor ...... 6
2 The Harmonic Cohomology Module 9 2.1 Derived Functors ...... 9 2.2 U 1 ...... 9 Km,n,Z 2.3 The Cokernel Interpretation ...... 11
3 Transformations of Boundary Graphs 13 3.1 Layerability ...... 13 3.2 Pasting Networks ...... 15
1 4 Making UΓ,R Nontrivial 18 4.1 Why Z[ζ3]? ...... 18 4.2 Maximal Minors ...... 18 4.3 Bipartite Graphs with 3 Interior Vertices ...... 20
5 A Combinatorial Interpretation of det C 24
6 Miscellaneous Algebraic Lemmas 25 6.1 Factorization ...... 25 6.2 Quotients by Principal Ideals ...... 25 6.3 Multiplying the Conductances by a Constant ...... 27
7 Questions that Remain 28
1 A Terms from Homological Algebra 30 A.1 Exact Sequences ...... 31 A.2 Derived Functors ...... 32 A.3 The Snake Lemma ...... 33
1 Graphs and Harmonic Functions 1.1 Graphs with Boundary A graph G consists of a vertex set V , an edge set E, and a one-to-one assignment of endpoints (an unordered pair of distinct vertices) to each edge. This definition rules out loops from a vertex to itself and multiple edges between two vertices; additionally, we assume in this paper that all graphs are finite. We use the notation v1 ∼ e ∼ v2 if v1 and v2 are the endpoints of e, and say v1 and v2 are incident to the edge e. We also write v1 ∼ v2 if there is an edge between v1 and v2, and we may denote this edge as v1v2. The valence of a vertex is the number of vertices that it is incident to. In this paper, we consider by default only connected graphs: those such that, for any vertices v and w, there is a sequence v, v1, . . . , vn, w in which each pair of consecutive vertices is joined by an edge. 0 0 A morphism of graphs f :(V,E) → (V ,E ) is a pair of functions fV : V → 0 0 V and fE : E → E such that if the endpoints of e are v1 and v2, the endpoints of fE(e) are fV (v1) and fV (v2). A graph with boundary, boundary graph, or bgraph is a graph G = (V,E) together with a partition V = ∂V ∪ int V into boundary vertices ∂V and interior vertices int V , with ∂V nonempty. A morphism of boundary graphs is a graph morphism f with two additional properties [2][p. 4]:
• f sends interior vertices to interior vertices. • For an interior vertex v, if the map is restricted to a map on the neighbors of v, the preimage of every neighbor of f(v) has the same size n, n ≥ 1. (In particular, this implies that the map must be surjective onto the image of the neighbors of the interior vertex.)
A sub-boundary-graph or sub-bgraph of a boundary graph is an ordi- nary subgraph such that the inclusion into the supergraph is a bgraph morphism. In particular, since the inclusion map must be one-to-one, it must be a bijec- tion when restricted to the neighbors of an interior vertex, by the second point above. Let G = (∂V ∪ int V,E) be a graph with boundary and let R be a com- mutative ring with 1. (This theory can be extended with little difficulty to non-commutative rings, but we do not consider any in this paper.) Then G can be given the structure of a network with a function γ : E → R − {0} that assigns a conductance to each edge. In many algebraic contexts, we take the conductances to be units, and that is assumed here unless stated otherwise.
2 Figure 1: There is no bgraph morphism from the 5-star to the 3-star, since it would have to be n-to-1 from 5 neighbors of the first interior vertex to 3 neighbors of the other. However, there is a bgraph morphism from the 6-star to the 3-star, which is 2-to-1 in the neighborhood of the interior vertex.
We refer to the network Γ = (∂V ∪ int V, E, γ). If no conductance function is specified on a network, it is assumed by default that all the conductances are 1. A network morphism is then a morphism of the underlying networks that preserves conductances. There are certain ways of reducing bgraphs to sub-bgraphs that we will want to consider [2, p.25]. A boundary edge is an edge between two bound- ary vertices. A boundary spike is an edge which has an interior vertex as one endpoint and a boundary vertex of valence 1 as the other. An isolated boundary vertex is a boundary vertex that is adjacent to no other vertices. A boundary cutpoint is a boundary vertex which, if removed, disconnects the graph. Often we will want to remove these types of edges/vertices from graphs. Deleting a boundary edge or an isolated boundary vertex means creating a new bgraph by removing that edge or vertex. Deleting (or contracting) a boundary spike means creating a new bgraph by removing the edge and bound- ary vertex associated with the spike and changing the interior vertex into a boundary vertex. Splitting a boundary cutpoint is a slightly more complicated procedure. First, removing the vertex (and its incident edges) creates a disconnected sub- graph. To each connected component, reattach a distinct copy of the deleted boundary vertex with the edges that originally connected it to that component. A graph is layerable if it can be reduced to the empty graph (the graph with no vertices) by a series of deletions of boundary edges, boundary spikes, and isolated boundary vertices. It is quasi-layerable if it can be reduced to the empty graph using these operations along with splitting of boundary cutpoints.1 If one of the connected components obtained by deleting a boundary cutpoint is a single interior vertex, the interior vertex must have originally had degree 1. We refer to this special case as an interior spike. Interior spikes are irrelevant
1The algebraic considerations in this paper suggest that quasi-layerability may be a more natural definition, but for now, we reserve “layerability” to refer to Jekel’s original definition.
3 Figure 2: Reducing a quasi-layerable graph by splitting a boundary cut- point, contracting boundary spikes, deleting boundary edges, contracting fur- ther spikes, deleting further edges, and then deleting the isolated vertices. From the second step on, the graphs are also layerable.
4 to the electrical properties of the network, so for most of this paper we will assume graphs do not have them.
1.2 Harmonic Functions
When the ring R = R and the conductances are restricted to positive values, this structure resembles an electrical network. With this in mind, for some general left R-module M, we can define a potential function on the network Γ to be any function u : V → M. Then for a specific potential, the current along an edge e from vertex v1 to v2 is given by γ(e)(u(v1)−u(v2)). (This linear relationship is derived from Ohm’s law, where our conductance is the reciprocal of resistance.) In investigating electrical networks, we consider assignments of potentials that satisfy Kirchhoff’s Current Law: the net current leaving an interior vertex along all edges is 0. Stated using our notation, X γ(e)(u(v) − u(w)) = 0 ∀v ∈ . e,w:v∼e∼w A function u that satisfies this property is called γ-harmonic, or just har- monic. We denote the set of harmonic functions on the network Γ with values in M by UΓ,M . V Note that UΓ,M ⊆ M has the structure of an R-module. The operator P γ(e)(u(v) − u(w)) which takes a function u to the resulting net e,w:v∼e∼w current at a vertex v is a homomorphism of R-modules, and so its kernel UΓ,M is a submodule. (In the case where R is not commutative, it can only be viewed as a Z-module, but the theory is mostly the same.)
1.3 The Dirichlet Problem and The Kirchhoff Matrix The (discrete) Dirichlet problem asks: given arbitrary potential values at the boundary vertices, is there a harmonic function which takes those values, and is it unique? The question can be simplified significantly with the use of matrices. As mentioned above, the relationship between the potential at a vertex and the net current there is linear, and it can be expressed simply as a matrix obtained from the values of γ. Given a network Γ = (∂V ∪int V, E, γ), label the vertices V = {v1, v2, . . . , vn} and define the n × n Kirchhoff matrix K = {kij} by P γ(e) i = j vi∼e kij = −γ(eij) vi is adjacent to vj through the edge eij 0 otherwise This definition is such that, given a potential function u and the vector ~u = (u(v1), u(v2), . . . , u(vn)), the resulting entry (K~u)i gives the net current at vertex vi.
5 It is standard to number the vertices in such a way that all of the boundary vertices come first. Additionally, K is symmetric, since its definition does not depend on the order of i and j. This allows K to be partitioned into a block matrix
AB K = BT C where: • A is a symmetric |∂V | × |∂V | matrix representing the conductances of edges between boundary vertices; • B is a |∂V |×| int V | matrix representing the conductances of edges between boundary vertices and interior vertices; • C is a symmetric | int V | × | int V | matrix representing the conductances of edges between interior vertices.
Then, partitioning ~u into (~u∂V , ~uint V ) in the same way, the condition of harmonicity can be restated using K: AB ~u∂V ∗ u is harmonic iff T = , B C ~uint V 0 where ∗ can be any value. This is the same as requiring that BT C ~u = 0, so we get a helpfully succinct description of the harmonic functions: T Proposition 1.1. UΓ,M = ker(B ,C)M , where the matrix is interpreted as an R-module homomorphism M V → M int V .
Returning to the Dirichlet problem, if we are given potential values ~u∂V on the boundary, then for the remaining interior values ~uint V to give a harmonic function requires that T B ~u∂V + C~uint V = 0 T C~uint V = −B ~u∂V The Dirichlet problem reduces to considering whether this equation can be solved for ~uint V . This naturally depends on the module M, but in the specific case that the ring R is a field, the problem has a unique solution for any bound- ary values exactly when det C 6= 0. In general, if det C = 0, Γ is said to be Dirichlet singular.
1.4 U as a Functor
As was mentioned above, for an R-module M, UΓ,M can also be interpreted as an R-module. Additionally, given some u ∈ UΓ,M and an R-module homomorphism φ : M → N, we can consider the function φ∗u = φ ◦ u : V → N; then since ! X X γ(e)(φ(u(v)) − φ(u(w)) = φ γ(e)(u(v) − u(w)) , e,w:v∼e∼w e,w:v∼e∼w
6 this new function is harmonic. We then have a homomorphism φ∗ : UΓ,M → UΓ,N . Additionally, for module homomorphisms φ : L → M and ψ : M → N, we have (ψ ◦ φ)∗u = ψ ◦ φ ◦ u = ψ∗φ∗u. As a result, for any fixed network Γ, UΓ,− is a (covariant) functor R-Mod → R-Mod. However, we can also ask what happens when we fix the module and vary the network.
Proposition 1.2. For a fixed R-module M, U−,M is a contravariant functor Network → R-Mod, where for a network morphism f :Γ → Γ0 the image under the functor is the pullback f ∗(u) = u ◦ f.
Proof. First, we must show that, for u ∈ UΓ0,M , u◦f is harmonic on Γ. Consider an arbitrary interior vertex v of Γ; then consider the sum X γ(e)(u(f(v)) − u(f(w))) e,w:v∼e∼w Since v is an interior vertex of Γ, f(v) must be an interior vertex of Γ0. Addi- tionally, f is n-to-1, for some fixed n, on the neighbors of v (thus also the edges incident to v). Combining this with the fact that f preserves conductances, X X γ(e)(u(f(v)) − u(f(w))) = nγ(e0)(u(f(v)) − u(w0)) 0 0 0 0 e,w:v∼e∼w e ,w :f(v)∼e ∼w X = n γ(e0)(u(f(v)) − u(w0)) 0 0 0 0 e ,w :f(v)∼e ∼w = 0.
Having shown this, it follows straightforwardly that U−,M is a contravariant functor, since for network morphisms f :Γ → Γ0 and g :Γ0 → Γ00,(g ◦ f)∗(u) = u ◦ g ◦ f = g∗u ◦ f = f ∗(g∗(u)). In this paper, we are primarily interested in how U acts as a functor on R-modules, with Γ fixed.
Proposition 1.3. UΓ,− is left exact. Proof. Let 0 → L → M → N → 0 be a short exact sequence of R-modules. Then since the submatrix BT C of the Kirchhoff matrix is a linear map, this diagram commutes (where the horizontal maps are simply coordinatewise applications of those in the original sequence):
0 LV M V N V 0
T T T (B C)L (B C)M (B C)N
0 Lint V M int V N int V 0
The rows are still exact and the kernels of the vertical maps are UΓ,L, UΓ,M , and UΓ,N , hence by the Snake Lemma, 0 → UΓ,L → UΓ,M → UΓ,N is exact.
7 However, U is not right exact, because it may not preserve the surjectiv- ity of a map. Consider the Z-modules Q and Q/Z with π : Q Q/Z the natural projection homomorphism and let Γ be the following network (with all conductances 1): 0
0 1/2
0
Displayed on the graph is a potential function taking values in Q/Z. This function is harmonic, since the net current at the left interior vertex is 0+0 = 0 and the net current at the right interior vertex is 1/2 + 1/2 = 0. However, any Q-valued harmonic function that projected to this one would have to be of the form
m1
m3 m4 + 1/2
m2
for integers m1, m2, m3, m4, and for it to be harmonic would require (in Q)
m3 − m1 + m3 − m2 = 0
m4 + 1/2 − m1 + m4 + 1/2 − m2 = 0
or, rearranging into the kind of form that we will use for the rest of the paper,
2m3 = m1 + m2
2m4 + 1 = m1 + m2
This would imply that an even and odd integer are equal, which is a contra- diction. So this particular Q/Z-valued harmonic function is not the image of a Q-valued one under the map induced by Q → Q/Z. We say that this function does not lift. In the remainder of this paper, we will explore the implications of this.
8 2 The Harmonic Cohomology Module 2.1 Derived Functors Given a left exact functor between categories of modules, we can construct derived functors through a process described in the Appendix, which describe the failure of the functor to be exact. In this section, we will take a first look at the derived functors of the functor UΓ,−. In most of this paper, we will be considering for our ring a principal ideal domain R and using R for our module (over itself). In this case, Proposition A.2 gives that we have a simple injective resolution and thus an easy description of the derived functors. Proposition 2.1. For R a PID, F its field of fractions, and π : F → F/R the projection onto the quotient module, consider the induced map π∗ : UΓ,F → UΓ,F/R. Then the first derived functor of UΓ,− at R is given by U 1 ∼ Γ,F/R UΓ,R = im π∗ j ∼ Additionally, UΓ,R = 0 for j ≥ 2. Proof. By Proposition A.2, an injective resolution of R is given by 0 → R → F → F/R → 0 → 0 → ... To obtain the derived functors, we apply the functor to the terms of the sequence, omitting 0, to get the cochain complex π∗ 0 → UΓ,F −→UΓ,F/R → 0 → 0 → ... The first derived functor is then defined to be the quotient of the kernel of the third arrow (all of UΓ,F/R) by the image of the second arrow (that is, 1 im(π∗)). This gives UΓ,R. The higher derived functors are then given by the cohomology at terms which are 0, which is trivially 0.
With this statement, we can now see that the harmonic function shown at the end of the previous section, a Q/Z-valued function on the network Γ which is not in the image of the map π , represents a non-identity element of U 1 . ∗ Γ,Z 1 For a fixed network Γ with conductances in the ring R, the module UΓ,R is of most interest to us, and we refer to it as the discrete harmonic cohomology j module, or just the harmonic cohomology. (Similarly, UΓ,R can be referred to as the jth harmonic cohomology.) In the next section, we will show an example of computing harmonic cohomology from the definition.
2.2 U 1 Km,n,Z As a practice calculation we will compute the first cohomology module for the complete bipartite bgraph Km,n. We say a bgraph is bipartite if there are no interior to interior or boundary to boundary edges. Thus the complete bipartite bgraph Km,n is the bipartite bgraph with m boundary vertices and n interior vertices, and every possible interior to boundary edge. To make the
9 calculation easier, we will prove the following Lemma first. The basic idea is that since C is invertible, given any Q/Z harmonic function we can find a Q harmonic function with the same boundary values; so up to the image of a Q harmonic function, each Q/Z harmonic function is 0 on the boundary. condition We will use the slightly informal notation UΓ,M to denote the collection of all M harmonic functions on Γ that satisfy “condition.” Lemma 2.2. Let R be a PID, F its field of fractions, and π : F V → (F/R)V the natural R-module homomorphism. Suppose Γ is a Dirichlet non-singular network with unit conductances in R. Then
0 on bdry UΓ,F/R U 1 =∼ Γ,R R on bdry π UΓ,F
Proof. Consider the following commutative diagram:
R on bdry ∂V 0 UΓ,F UΓ,F (F/R) 0 π π id
0 on bdry ∂V 0 UΓ,F/R UΓ,F/R (F/R) 0
The maps are the natural ones. We claim the rows are exact. Exactness is easy to see everywhere except at the rightmost points of both rows. To ∂V check exactness there, we first note the map UΓ,F → F is surjective since the Dirichlet problem has a solution, and the projection map F ∂V → (F/R)∂V is ∂V surjective, hence the composition UΓ,F → (F/R) is surjective. Hence exact- ∂V ness of the top row. Now this map UΓ,F → (F/R) is the composition of two ∂V maps UΓ,F → UΓ,F/R → (F/R) , and the whole thing being surjective implies the second map is surjective, which gives surjectivity of the bottom row. Thus we may apply the Snake Lemma, and since the furthest right vertical map has trival kernel and cokernel, we get that
U 0 on bdry 0 → Γ,F/R → U 1 → 0 R on bdry Γ,R π UΓ,F is exact. This gives the result.
Theorem 2.3. Denote by Km,n the network with m boundary vertices, n inte- rior vertices, and every possible interior to boundary edge, all with conductances 1 ∈ Z. Then we have that
U 1 =∼ ( /m)n−1. Km,n,Z Z
10 Proof. Since the bgraph is bipartite, C is invertible over Q, and thus by Lemma 2.2 we can assume all Q/Z harmonic functions are 0 on the boundary. So for each such function, it must be that for all v ∈ int V X mu(v) = u(w) = 0 in Q/Z. w∼v The above equation has solutions u(v) = 0, 1/m, . . . , (m − 1)/m in Q/Z. Since this equation holds for all v ∈ int V and each is independent of the others, it follows that 0 on bdry ∼ n UΓ,F/R = (Z/m) . Now suppose u is a Q harmonic function that has values in R on the boundary, then we again get that it is necessary and sufficient for it to satisfy, for each v ∈ int V , X X mu(v) = u(w) = u(w) in Q. w∼v w∈∂V But since this is in Q, it follows that u is constant on the interior. Applying the same reasoning as before when we project from F to F/R, we then get that u(v) = 0, 1/m, . . . , (m − 1)/m for all v ∈ int V . Thus R on bdry ∼ π UΓ,F = Z/m. Thus by Lemma 2.2 we get the result.
2.3 The Cokernel Interpretation As we saw in the previous section, computing harmonic cohomology from the definition is somewhat unwieldy, and that particular demonstration was depen- dent on the specific family of graphs. Using the snake lemma, we can describe the cohomology in terms of the cokernel of a matrix, allowing us to consider it using the tools of linear algebra and modules. Theorem 2.4. Suppose R is a PID, and F is its field of fractions. Consider some network Γ such that BT C is surjective as a map F V → F int V (in particular, this holds for any Dirichlet nonsingular network.) Then 1 ∼ BT C UΓ,R = coker R T V I where the matrix (B C)R is interpreted as a map R → R . Proof. We have the following commutative diagram, where f, g, h are the natural projections:
(BT C) 0 UΓ,F F V F intV 0 f g h
V intV 0 UΓ,F/R (F/R) (F/R) 0 (BT C)
11 V V Certainly the inclusion maps UΓ,F → F and UΓ,F/R → (F/R) are injec- V tive. The exactness of the upper row at F follows from the definition of UΓ,F as a kernel, and similarly for the lower row at (F/R)V . To finish showing exact- ness of the rows, it then suffices to show that the two maps given by BT C are surjective. Since det(C) is nonzero, it is a unit in F and hence C is invert- ible as a matrix over F . Then for any u ∈ F int V , there is a vector C−1u such T that BT C 0 C−1u = u, so the upper map is surjective. Now, for any u ∈ (F/R)int V , since h is surjective we can pull u back to an element in F int V , which we can pull back to an element in F V and then map to some element v ∈ (F/R)V by g. By commutativity, v is mapped by (BT C) to u. Hence, the lower map is also surjective. Thus, the Snake Lemma yields the following exact sequence
0 → ker f → ker g → ker h → coker f → coker g → coker h → 0
But now observe that both h and g are onto and thus coker g =∼ coker h =∼ 0. It ∼ V ∼ intV 1 is also easy to see that ker g = R and ker h = R . Also, coker f = UΓ,R. Thus we get the following arrows
T (B C)R V intV 1 R R UΓ,R
So from the first isomorphism theorem and exactness we get that
1 ∼ intV intV 1 ∼ intV T UΓ,R = R / ker(R → UΓ,R) = R / im (B C)R. which is the required cokernel. While this formulation does not apply if BT C is not surjective over F , it still covers a wide swath of cases. As a concluding example, we will repeat the calculation of U 1 using its matrix. Km,n,Z First, for Γ = Km,n with conductances 1, the matrix is given by
−1 −1 ... −1 m 0 ... 0 −1 −1 ... −1 0 m . . . 0 ...... ...... −1 −1 ... −1 0 0 . . . m
As long as we remain in Z, performing column reduction operations on the matrix does not change its column space, and thus does not change the cokernel. We can thus subtract the first column from all of the others in the left block to make them 0:
−1 0 ... 0 m 0 ... 0 −1 0 ... 0 0 m . . . 0 ...... ...... −1 0 ... 0 0 0 . . . m
12 In determining the column space, we can ignore 0 columns, so we can freely delete these from the matrix. Then we can add m times the leftmost column to the first column of the right block to get
−1 0 0 ... 0 −1 −m m . . . 0 . . . . . . . . . . −1 −m 0 . . . m
By adding the third, fourth,. . . columns to the second, we cancel out all of its nonzero coefficients, reducing it to 0. Removing this column gives
−1 0 ... 0 −1 m . . . 0 . . . . . . . . −1 0 . . . m