arXiv:0806.3881v3 [math.OA] 15 May 2009 2 1 prtr lcrclrssac ewr,e network, resistance electrical operator, upre yteUiest fIw eateto Mathemati of Department Iowa of University the by supported 73,6J0 08,6D5 66,8C0 11.Secondary 81Q10. 82C10, 46L60, 60D05, 60J85, 60J10, 47B32, esnegfroant ogrneodr atc models lattice order, m long-range banded ferromagnet, walk, Heisenberg random theory, graph spectral trip trees, Gel’fand system, embedding, Neumann von self-adjoint, tially partial space self-adjoint, path essentially orthogonality, paths, space, infinite boundary, t Martin Fourier kernel, Poisson analysis, harmonic function, harmonic geometry, e od n phrases. and words Key Date 2000 h oko EJwsprilyspotdb S rn DMS-045 grant NSF by supported partially was PETJ of work The PRTRTER FEETIA EITNENETWORKS RESISTANCE ELECTRICAL OF THEORY OPERATOR : ahmtc ujc Classification. Subject Mathematics a ,21.Cenversion. Clean 2019. 4, May eiae oBbPwr,i eonto fhspoern sp pioneering his of recognition in Powers, Bob to Dedicated AL .T OGNE N RNP .PEARSE J. P. ERIN AND JORGENSEN T. E. PALLE iihe om rp nry iceeptnilter,g theory, potential discrete energy, graph form, Dirichlet ff cierssac,rssac erc nqeeso curre of uniqueness metric, resistance resistance, ective rmr:0C0 57,3C0 23,4E2 62,47B25, 46F25, 46E22, 42C30, 31C20, 05C75, 05C50, Primary: smty ieroeaos none ieroeaos e operators, linear unbounded operators, linear isometry, 1 M-tts erdcn enl,fae,self-similar. frames, kernels, reproducing KMS-states, , e igdHletsae inrtasom dynamical transform, Wiener space, Hilbert rigged le, asom onayrpeetto,Mro process, Markov representation, boundary ransform, esr,wihe rps onayter,Hilbert theory, boundary graphs, weighted measure, tie,rssac om,rssac erc isotropic metric, resistance forms, resistance atrices, sNFVGEgatDMS-0602242. grant VIGRE NSF cs 13,4C5 73,5C3 22,82C41. 82C22, 52C23, 47B39, 42C25, 31C35, : 51 h oko PPwspartially was EPJP of work The 7581. irit. ahLpain transfer Laplacian, raph t,di nts, ff usion ssen- 2
Abstract. An electrical resistance network (ERN) is a weighted graph (G, c). The con- ductance function cxy weights the edges, which are then interpreted as resistors of possibly varying strengths. The effective resistance metric R(x, y) is the natural notion of distance between two vertices x, y in the ERN. The space of functions of finite energy (modulo constants) is a Hilbert space with inner product , which we call the energy space . The evaluation functionals on give E HE HE rise to a reproducing kernel vx for the space. Once a reference vertex o is fixed, these { } functions vx satisfy ∆vx = δx δo, where ∆ is the network Laplacian. This kernel yields − a detailed description of the structure of = in arm, where in is the closure HE F ⊕ H F of the space of finitely supported functions and arm is the closed subspace of harmonic H functions. The energy splits accordingly into a “finite part” expressed as a sum taken over E the vertices, and an “infinite part” expressed as a limit of sums. Intuitively, the latter part corresponds to an integral over some sort of boundary bdG, which is developed explicitly in 8. The kernel vx also allows us to recover easily many known (and sometimes difficult) § { } results about . As does not come naturally equipped with a natural o.n.b., we HE HE provide candidates for frames (and dual frames) when working with an infinite ERN. In particular, the presence of nonconstant harmonic functions of finite energy leads to different plausible definitions of the effective resistance metric on infinite networks. We characterize the free resistance RF (x, y) and the wired resistance RW (x, y) in terms of Neumann or Dirichlet boundary conditions on a certain operator. (In the literature, these correspond to the limit current and minimal current, resp.) We develop a library of equivalent formulations for each version. Also, we introduce the “trace resistance” RS (x, y), computed in terms of the trace of the Dirichlet form to finite subnetworks. This E provides a finite approximation which is more accurate from a probabilistic perspective, and gives a probabilistic explanation of the discrepancy between RF and RW . For R = RF or R = RW , the effective resistance is shown to be negative semidefinite, so that it induces an inner product on a Hilbert space into which it naturally embeds. We show that for (G, RF ), the resulting Hilbert space is and for (G, RW ) it is in. Under HE F the free embedding, each vertex x is mapped to the element vx of the energy kernel; under the wired embedding it is mapped to the projection fx of vx to in. This establishes as F HE the natural Hilbert space in which to study effective resistance. We obtain an analytic boundary representation for elements of arm in a sense anal- H ogous to that of Poisson or Martin boundary theory. We construct a Gel’fand triple S ⊆ S ′ and obtain a probability measure P and an isometric embedding of into H2E ⊆ HE L (S ′, P). This gives a concrete representation of the boundary in terms of the measures (1 + vx )dP S / in, where xn is a sequence tending to infinity. n ∈ ′ F { } The spectral representation for the graph Laplacian ∆ on is drastically different HE from the corresponding representation on ℓ2. Since the ambient Hilbert space is defined HE by the energy form, many interesting phenomena arise which are not present in ℓ2; we highlight many examples and explain why this occurs. In particular, we show how the deficiency indices of ∆ as an operator on indicate the presence of nontrivial boundary HE of an ERN, and why the ℓ2 operator theory of ∆ does not see this. Along the way, we prove that ∆ is always essentially self-adjoint on the ℓ2 space of functions on an ERN, and examine conditions for the network Laplacian and its associated transfer operator to be bounded, compact, essential self-adjoint, etc. We consider two approaches to measures on spaces of infinite paths in an ERN. One arises from considering the transition probabilities of a random walk as determined directly by the network, i.e., p(x, y) = cxy/ y x cxy. The other applies only to transient networks, ∼ and arises from considering the transitionP probabilities induced by a unit flow to infinity. The latter leads to the notion of forward-harmonic functions, for which we also provide a characterization in terms of a boundary representation. Using our results we establish precise bounds on correlations in the Heisenberg model for quantum spin observables, and we improve earlier results of R. T. Powers. Our focus is on the quantum spin model on the rank-3 lattice, i.e., the ERN with Z3 as vertices and with edges between nearest neighbors. This is known as the problem of long-range order in the physics literature, and refers to KMS states on the C∗-algebra of the model. OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 3
1 Contents
2 1. Introduction 5 3 1.1. Outline 5 4 1.2. What this paper is about 13 5 1.3. What this paper is not about 16 6 1.4. General remarks 20 7 1.5. Acknowledgements 20 8 2. Electrical resistance networks 22 9 2.1. The electrical resistance network model 22 10 2.2. The energy 24 11 3. Currents and potentials on electrical resistance networks 27 12 3.1. Currents on electrical resistance networks 27 13 3.2. Potential functions and their relationship to current flows. 28 14 3.3. The compatibility problem 31 15 4. The energy Hilbert space 35 16 4.1. The evaluation operator Lx and the reproducing kernel vx 37 17 4.2. The finitely supported functions and the harmonic functions 38 18 4.3. Relating the energy form to the Laplacian 40 19 5. The resistance metric 50 20 5.1. Resistance metric on finite networks 50 21 5.2. Resistance metric on infinite networks 53 22 5.3. Harmonic resistance 64 23 5.4. Comparison of resistance metric to other metrics 65 24 5.5. Generalized resistance metrics 66 25 6. von Neumann construction of the energy space 68 HE 26 6.1. von Neumann’s embedding theorem 68 27 6.2. as an invariant of G 69 HE 28 7. The boundary and boundary representation 72 29 7.1. Motivation and outline 72 30 7.2. Gel’fand triples and duality 74 31 7.3. The boundary as equivalence classes of paths 82 32 7.4. The structure of ′ 84 SG 33 8. The Laplacian on 86 HE 34 8.1. Properties of ∆ on 86 HE 35 8.2. The defect space of ∆ 90 V 36 8.3. Dual frames and the energy kernel 91 2 37 9. The ℓ theory of ∆ and the transfer operator 97 2 38 9.1. The Laplacian on ℓ (1) 97 39 9.2. The transfer operator 100 2 40 9.3. The Laplacian and transfer operator on ℓ (c) 104 41 10. The dissipation space D and its relation to 109 H HE 42 10.1. The structure of 110 HD 43 10.2. The divergence operator 113 44 10.3. Analogy with calculus and complex variables 115 45 10.4. Solving potential-theoretic problems with operators 116 46 11. Probabilistic interpretations 118 47 11.1. The path space of a general random walk 118 48 11.2. The forward-harmonic functions 122 49 12. Examples and applications 128 50 12.1. Finite graphs 128 51 12.2. Infinite graphs 132 52 13. Infinite trees 134 4 PALLEE.T.JORGENSENANDERINP.J.PEARSE
14. Lattice networks 140 1 14.1. Simple lattice networks 140 2 14.2. Noncompactness of the transfer operator 148 3 14.3. Non-simple integer lattice networks 150 4 14.4. Defect spaces 154 5 15. Application to magnetism and long-range order 157 6 2 15.1. Kolmogorov construction of L (Ω, P) 157 7 15.2. The GNS construction 159 8 15.3. Magnetism and long-range order in electrical resistance networks 161 9 15.4. KMS states 163 10 16. Future directions 166 11 Appendix A. Some functional analysis 168 12 A.1. von Neumann’s embedding theorem 168 13 Appendix B. Some operator theory 170 14 B.1. Projections and closed subspaces 170 15 B.2. Partial isometries 171 16 B.3. Self-adjointness for unbounded operators 171 17 B.4. Banded matrices 175 18 Appendix C. Navigation aids for operators and spaces 178 19 C.1. A road map of spaces 178 20 C.2. A summary of the operators on various Hilbert spaces 178 21 References 179 22 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 5
1 1. Introduction
2 This article is dedicated to the construction of an operator-theoretic context for study of 3 certain potential-theoretic function spaces on graphs, and an investigation of the resulting 4 structures. The primary object of study is an electrical resistance network (ERN), a graph 5 with weighted edges. Our foundation is the effective resistance metric as the intrinsic 6 notion of distance, and we approach the analysis of the ERN by studying the space of 7 functions on the vertices which have finite (Dirichlet) energy. There is a large existing 8 literature on this subject, but ours is unique in several respects, most of which are due to 9 the following.
10 We use the effective resistance metric to find canonical Hilbert spaces of functions • 11 associated with the ERN. 12 We adhere to the intuition arising from the metaphor of electrical resistance net- • 13 works, including Kirchhoff’s Law and Ohm’s Law. 14 We apply the results of our Hilbert space construction to the isotropic Heisenberg • 15 ferromagnet and prove a theorem regarding long-range order in quantum statistical 16 mechanics for certain lattice networks. 17 It is known (see [LP09] and the references therein) that the resistance metric is • 18 unique for finite graphs and not unique for certain infinite networks. We are able 19 to clarify and explain the difference in terms of certain Hilbert space structures, 20 and also in terms of Dirichlet vs. Neumann boundary conditions for a certain 21 operator. Additionally, we introduce trace resistance, and harmonic resistance and 22 relate these to the aforementioned.
23 By using the intrinsic inner product (associated to the effective resistance) we are able to 24 obtain results which are more physically realistic than many found elsewhere in the litera- 2 25 ture. This inner product is quite different than the standard ℓ inner product for functions 26 defined on the vertices of a graph, and holds many surprises. Many of our results apply 27 much more generally than those already present in the literature. The next section elabo- 28 rates on these rather vague remarks and highlights the advantages and differences inherent 29 in our approach, in a variety of circumstances. 30 This work is uniquely interdisciplinary, and as a consequence, we have made effort to 31 address the union (as opposed to the intersection) of several disparate audiences: graph 32 theory, resistance networks, spectral geometry, fractal geometry, physics, probability, un- 33 bounded operators in Hilbert space, C*-algebras, and others. It is inevitable that parts of 34 the background material there will be unknown to some readers and so we have included 35 the appendices to mediate this. After presenting our results at various talks, we felt that the 36 inclusion of this material would be appreciated by most.
37 1.1. Outline. 38 39 2 — Electrical resistance networks. We introduce the ERN as a connected simple § 0 1 40 graph G = G , G equipped with a positive weight function c on the edges. The edges 1 0 { 0 } 41 G G G are ordered pairs of vertices, so c is required to be symmetric. Hence, ⊆ × 1 42 each edge (x, y) G is interpreted as a conductor with conductance cxy (or a resistor ∈1 43 with resistance c−xy . Heuristically, smaller conductances (or larger resistances) correspond 44 to larger distances; see the discussion of 5 just below. We make frequent use of the § 45 weight that c defines on the vertices via c(x) = y x cxy, where y x indicates that 1 ∼ ∼ 46 (x, y) G . The graphs we are most interested in areP infinite graphs, but we do not make ∈ 47 any general assumptions of regularity, group structure, etc. We require that c(x) is finite at 6 PALLEE.T.JORGENSENANDERINP.J.PEARSE
0 each x G , but we do not generally require that the degree of a vertex be finite, nor that 1 ∈ c(x) be bounded. 2 In the “cohomological” tradition of von Neumann, Birkhhoff, Koopman, and others 3 [vN32c, Koo36, Koo57], we study the ERN by analyzing spaces of functions defined on 4 it. These are constructed rigourously as Hilbert spaces in 6.1; in the meantime we collect 5 0 § some results about functions u, v : G R defined on the vertices. The network Laplacian 6 → (or discrete Laplace operator) operates on such a function by taking v(x) to a weighted 7 average of its values at neighbouring points in the graph, i.e., 8
v(x) v(y) (∆v)(x) := c (v(x) v(y)) = − , (1.1) xy − c 1 Xy x Xy x −xy ∼ ∼ 1 where x y indicates that (x, y) G . (The rightmost expression in formula (1.1) is written 9 ∼ ∈ so as to resemble the familiar difference quotients from calculus.) This is the usual second- 10 difference operator of numerical analysis, when adapted to a network. There is a large 11 literature on discrete harmonic analysis (basically, the study of the graph/network Lapla- 12 cian) which include various probabilistic, combinatoric, and spectral approaches. It would 13 be difficult to give a reasonably complete account, but the reader may find an enjoyable ap- 14 proach to the probabilistic perspective in [Spi76,Tel06], the combinatoric in [ABR07], the 15 analytic in [Fab06], and the spectral in [Chu01, GILb]. More sources are peppered about 16 the relevant sections below. Our formulation (1.1) differs from the stochastic formulation 17 often found in the literature, but the two may easily be reconciled; see (2.6). 18 Together with its associated quadratic form, the bilinear (Dirichlet) energy form 19
1 (u, v) := c (u(x) u(y))(v(x) v(y)) (1.2) E 2 xy − − xXG0 Xy x ∈ ∼ 0 acts on functions u, v : G R and plays a central role in the (harmonic) analysis on 20 → (G, c). (There is also the dissipation functional D, a twin of which acts on functions 21 1 E defined on the edges G and is introduced in the following section.) The first space of 22 functions we study on the ERN is the domain of the energy, that is, 23
0 . dom := u : G R .. (u) < . (1.3) E { → E ∞} In 6.1, we construct a Hilbert space from the resistance metric (andshowit to be a canoni- 24 § cal invariantfor (G, c) in 6.2), thereby recoveringthe familiar result that dom is a Hilbert 25 § E space with inner product . (Actually, this is not quite true, as is only a quasinorm; see 26 E E the discussion of 6.1 just below for a more accurate description.) 27 § For finite graphs, we prove the simple and folkloric key identity which relates the energy 28 and the Laplacian: 29
(u, v) = u, ∆v = ∆u, v , u, v dom , (1.4) E h i1 h i1 ∈ E 2 where u, ∆v 1 = x G0 u(x)∆v(x) indicates the standard ℓ inner product. The formula 30 h i ∈ (1.4) is extended toP infinite networks in Theorem 4.41 (see (1.9) for a preliminary discus- 31 sion), where a third term appears. Indeed, understanding the mysterious third term is the 32 motivation for most of this investigation. 33
3 — Currents and potentials on electrical resistance networks. We collect several 34 § well-known and folkloric results, and reprove some variants of these results in the present 35 context. Currents are introduced as skew-symmetric functions on the edges; the intuition 36 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 7
1 is that I(x, y) = I(y, x) > 0 indicates electrical current flowing from x to y. In marked − 2 contrast to common tradition in geometric analysis [ABR07, PS07], we do not fix an ori- 3 entation. For us, an orientation is a choice of one of (x, y), (y, x) for each edge, and hence { } 4 just a notation to be redefined as convenient. In particular, any nonvanishing current de- 5 fines an orientation; one makes the choice so that I is a positive function. At this point we 6 give the definition of the dissipation, an inner product defined for functions on the edges, 7 and its associated quadratic form:
1 1 2 D(I) = c− I(x, y) . (1.5) 2 xy (x,Xy) G1 ∈ 8 Most of our results in this section are groundwork for the sections to follow; several 9 results are folkloric or obtained elsewhere in the literature. We include items which relate 10 directly to results in later sections; the reader seeking a more well-rounded background is 11 directed to [LPW08, LP09, Soa94, CdV98, Bol98] and the excellent elementary introduc- 12 tion [DS84]. After establishing the Hilbert space framework of 4, we exploit the close § 13 relationship between the two functionals and D, and use operators to translate a prob- E 14 lem from the domain of one functional to the domain of the other. We also introduce 15 Kirchhoff’s Law and Ohm’s Law, and in 3.3 we discuss the related compatibility prob- § 16 lem: every function on the vertices induces a function on the edges via Ohm’s Law, but not 17 every function on the edges comes from a function on the vertices. This is related to the 18 fact that most currents are not “efficient” in a sense which can be made clear variationally 19 (cf. Theorem 3.26) and which is important in the definition of effective resistance metric 20 in Theorem 5.2. We recover the well-known fact that the dissipation of an induced current 21 is equal to the energy of the function inducing it in Lemma 3.16; this is formalized as an 22 isometric operator in Theorem 10.12. We show that the equation
∆v = δα δω (1.6) − 23 always has a solution; we call such a function a dipole. In(1.6) and everywhere else, we 0 24 use the notation δ to indicate a Dirac mass at x G , that is, x ∈ 1, y = x, δx = δx(y) := (1.7) 0, else. 25 Proving the existence of dipoles allows us to fill gaps in [Pow76a,Pow76b] (see 1.2.1 just § 26 below) and extend the definition of effective resistance metric in Theorem 5.2 to infinite 27 dimensions. 28 As is discussed at length in Remark 3.11, the study of dipoles, monopoles, and harmonic 29 functions is a recurring theme of this paper:
∆v = δα δω, ∆w = δω, ∆h = 0. − − 0 30 As mentioned above, for any network G and any vertices x, y G , there is a dipole in ∈ 31 dom . However, dom does not always contain monopoles or nonconstant harmonic E E 32 functions; the existence of monopoles is equivalent to transience of the network [Lyo83]; 33 we give a new criterion for transience in Lemma 4.51. In Theorem 14.5, we show that d 34 the integer lattice networks (Z , 1) support monopoles iff d 3, but in Theorem 14.17 we d ≥ 35 show all harmonic functions on (Z , 1) are linear and hence do not have finite energy. (Both 36 of these results are well known; the first is a famous theorem of Polya — we include them 37 for the novelty of method of proof.) In contrast, the binary tree in Example 13.4 support 8 PALLEE.T.JORGENSENANDERINP.J.PEARSE monopoles and nontrivial harmonic functions, both of finite energy (any network support- 1 ing nontrivial harmonic functions also supports monopoles, cf. [Soa94, Thm. 1.33]). It is 2 apparent that monopoles and nontrivial harmonic functions are sensitive to the asymptotic 3 geometry of (G, c). 4
4 — The energy Hilbert space . We use the natural Hilbert space structure on the 5 § HE space of finite-energy functions (with inner product given by ) to reinterpret previous 6 E results as claims about certain operators, and thereby clarify and generalize results from 7 2– 3. This is the energy space . 8 § § HE We construct a reproducing kernel for from first principles (i.e., via Riesz’s Lemma) 9 0 HE in 4.1. If o G is any fixed reference point, define vx to be the vector in which 10 § ∈ HE corresponds (via Hilbert space duality) to the evaluation functional Lx: 11
L u := u(x) u(o). x − Then the functions v form a reproducing kernel, and v is a solution of the discrete 12 { x} x Dirichlet problem ∆v = δ δ . Although these functions are linearly independent, they 13 x x − o are usually neither an orthonormal basis (onb) nor a frame. However, the span of v is 14 { x} dense in dom and appears naturally when the energy Hilbert space is constructed from 15 E the resistance metric by von Neumann’s method; cf. 6.1. Note that the Dirac masses 16 § δ 0 , which are the usual candidates for an onb, are not orthogonal with respect to the 17 { x}G energy inner product (1.2); cf. (2.11). In fact, Theorem 4.47 shows that δ 0 may not 18 { x}G even be dense in the energy Hilbert space! Thus, v is the only canonical choice for a 19 { x} representing set for functions of finite energy. 20 In 4.2 we use the Hilbert space structure of to better understand the role of the 21 § HE nontrivial harmonic functions. In particular, Lemma 4.22 shows that we may decompose 22 into the functions of finite support ( in) and the harmonic functions of finite energy 23 HE F ( arm): 24 H = in arm. (1.8) HE F ⊕H In 4.3, we prove a discrete version of the Gauss-Green formula (Theorem 4.41) which 25 § appears to be absent from the literature: 26
(u, v) = u(x)∆v(x) + u(x) ∂v (x), u , v (1.9) E ∂Ò ∀ ∈HE ∈ MP xXG0 xXbd G ∈ ∈ ∂v where (x) denotes the normal derivative of v, and is a space containing span vx ; see 27 ∂Ò MP { } 4.3 for precise definitions. For the moment, both the boundary and the normal derivatives 28 § are understood as limits (and hence vanish trivially for finite graphs); we will be able to 29 define these objects more concretely via techniques of Gel’fand in 7. 30 § It turns out that the boundary term (that is, the rightmost sumin(1.9)) vanishes unless 31 the network supports nontrivial harmonic functions (that is, nonconstant harmonic func- 32 tions of finite energy). More precisely, in Theorem 4.47 we prove that there exist u, v 33 ∂v ∈HE for which u , 0 if and only if the network is transient. That is, the random walk 34
bd G ∂Ò on the networkP with transition probabilities p(x, y) = cxy/c(x) is transient. We also give 35 several other equivalent conditions for transience, in 4.3.1. 36 § It is easy to prove (see Corollary 4.65) that nontrivial harmonic functions cannot lie in 37 2 0 2 0 ℓ (G ). This is why we do not require u, v ℓ (G ) in general, and why we stringently 38 ∈ avoid including such a requirement in the definition of the domain of the Laplacian. Such 39 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 9
1 a restriction would remove the nontrivial harmonic functions from the scope of our anal- 2 ysis, and we will see that they are at the core of some of the most interesting phenomena 3 appearing on an infinite ERN.
4 5 — Effective resistance metric. The effective resistance metric R is foundational to § 5 our study, instead of the shortest-path metric more commonly used as graph distance. The 6 shortest-path metric on a weighted graph is usually defined to be the sum of the resistances 7 in any shortest path between two points. The effective resistance metric is also defined via 8 c, but in a more complicated way. The crucial difference is that the effective resistance 9 metric reflects both the topology of the graph and the weighting c; two points are closer to- 10 gether when there is more connectivity (more paths and/or paths with greater conductance) 11 between them. The effective resistance metric is a much more accurate way to measure 12 distance when travel from point x to point y can be accomplished simultaneously through 13 many paths, for example, flow of electrical current, fluid diffusion through porous media, 14 or data transfer over the internet. 15 In 5.1, we give a multifarious definition of the effective resistance metric R, which § 16 may be physically characterized as the voltage drop between two vertices when electrical 17 leads with a fixed current are applied to the two vertices. Most of these formulations 18 appear elsewhere in the literature, but some appear to be specific to the physics literature, 19 some to probability, and some to analysis. We collect them and prove their equivalence in 20 Theorem 5.2, including a couple new formulations that will be useful in later sections. 21 It is somewhat surprising that when these formulas are extended to an infinite network 22 in the most natural way, they are no longer equivalent. (Note that each of the six formulas 23 has both a free and wired version, but some appear much less natural in one version than F 24 in the other.) Some of the formulas lead to the “free resistance” R and others lead to the W 25 “wired resistance” R ; here we follow the terminology of [LP09]. In 5.2, we precisely § 26 characterize the types of extensions that lead to each, and explain this phenomenon in 27 terms of projections in Hilbert space, Dirichlet vs. Neumann boundary conditions, and via 28 probabilistic interpretation. Additionally, we discuss the “trace resistance” given in terms 29 of the trace of the Dirichlet form , and we study the “harmonic resistance” which is the F W E 30 difference between R and R and is not typically a metric.
31 6 — Construction of the energy space . In 6.1, we use a theorem of von Neumann § HE § 32 to give an isometric imbedding of the metric space (G, R) into ; cf. Theorem 6.1. For F W HE 33 infinite networks, (G, R ) embedsinto and(G, R ) embeds into in. In 6.2 we discuss HE F § 34 how this enables one to interpret as an invariant of the original ERN. HE
35 7 — The boundary bd G and boundary representation. We study the boundary bd G in § 36 terms of the Laplacian by reinterpreting the boundary term of (1.9) as an integral over a 37 space which contains . This gives a representation of bd G as a measure space whose HE 38 structure is well-studied. 39 In Theorem 7.1 of 7.1, we observe that an important consequence of (1.9) is the fol- § 40 lowing boundary representation for the harmonic functions:
u(x) = u ∂hx + u(o), (1.10)
∂Ò Xbd G
41 for u arm, where hx = P armvx is the projection of vx to arm;see(1.8). This formula ∈H H H 42 is in the spirit of Choquet theory and the Poisson integral formula and is closely related to 43 Martin boundary theory. 10 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Unfortunately, the sum in (1.10) is only understood in a limiting sense and so provides 1 limited insight into the nature of bd G. This motivates the development of a more concrete 2 expression. We use a self-adjoint extension ∆ of ∆ to construct a Gel’fand triple G 3 ∗ V S ⊆ P = ∆ 4 G′ and a Gaussian probability measure . Here, G : dom( ∗ ∞) is a suitable dense HE ⊆ S S V (Schwartz) space of “test functions” on the ERN, and ′ is the corresponding dual space 5 SG of “distributions” (or “generalized functions”). This enables us to identify bd G as a subset 6 of , and in Corollary 7.27, we rewrite (1.10) more concretely as 7 SG′
u(x) = u(ξ)hx(ξ) dP(ξ) + u(o), (1.11) Z SG′ again for u arm and with hx = P armvx. Thus we study the metric/measure structure 8 ∈ H H of G by examining an associated Hilbert space of random variables. This is motivated in 9 part by Kolmogorov’s pioneering work on stochastic processes (see 15.1) as well as on 10 § a powerful refinement of Minlos. The latter is in the context of the Gel’fand triples men- 11 tioned just above; see [Nel64] and 7.2 below. Further applications to harmonic analysis 12 § and to physics are given in 11– 15.3. 13 § §
8 — The Laplacian on . We study the operator theory of the Laplacian in some 14 § HE detail in 8.1, examining the various domains and self-adjoint extensions. We identify one 15 § domain for the Laplacian which allows for the choice of a particular self-adjoint extension 16 for the constructions in 7. Also, we give technical conditions which must be considered 17 § when the graph contains vertices of infinite degree and/or the conductance functions c(x) 18 0 is unbounded on G . This results in an extension of the Royden decomposition to = 19 HE in in arm, where in is the -closure of span δ δ and in is the orthogonal 20 F 1 ⊕F 2 ⊕H F 2 E { x − o} F 1 complement of in within in. Example 14.35 shows a case where in is not dense in 21 F 2 F F 2 in. 22 F In 8.2, we study the defect space of ∆ , that is, the space spanned by solutions to 23 § V ∆u = u. In 8.2.1, we relate the boundary term of (1.9) to the the boundary form 24 − § 1 βbd(u, v) := ∆∗ u, v u, ∆∗ v , u, v dom(∆∗ ) (1.12) 2i h V iE − h V iE ∈ V of classical functional analysis; cf. [DS88, XII.4.4]. This gives a way to detect whether 25 § or not a given network has a boundary by examining the deficiency indices of ∆. In Theo- 26 rem 8.19, weshowthatif ∆ fails to be essentially self-adjoint, then arm , 0 . In general, 27 H { } the converse does not hold: Corollary 9.28 shows that ∆ has no defect when deg(x) < 28 ∞ and c(x) is bounded. (Thus, any homogeneous tree of degree 3 or higher with constant 29 conductances provides a counterexample to the converse.) 30 In 8.3, we study the relation between the reproducing kernel v and the spectral 31 § { x} properties of ∆ and its self-adjoint extensions. In particular, we examine the necessary 32 conditions for vx to be a frame for , and the relation between vx and δx. 33 { } HE 2 9 — The ℓ theory of ∆ and T. We consider some results for ∆ and T as operators on 34 2 § 2 ℓ (1), where the inner product is given by u, ∆v := u(x)∆v(x) and on ℓ (c), where the 35 h i1 inner product is given by u, ∆v c := c(x)u(x)∆v(x).P 36 h i 2 We prove that the Laplacian is essentiallyP self-adjoint on ℓ (1) under very mild hy- 37 potheses in 9.1. The subsequent spectral representation allows us to give a precise char- 38 § acterization of the domain of the energy functional in this context. In 9.2, we examine 39 E § boundedness and compactness of ∆ and T in terms of the decay properties of c. The space 40 2 ℓ (c) considered in 9.3 is essentially a technical tool; it allows for a proof that the terms 41 § of the Discrete Gauss-Green formula are absolutely convergent and hence independent of 42 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 11
1 any exhaustion. However, it is also interesting in its own right, and we show an interesting 1 2 connection with the probabilistic Laplacian c− ∆. Results from this section imply that ∆ is 2 3 also essentially self-adjoint on , subject to the same mild hypotheses as the ℓ (1) case. HE 4 The energy Hilbert space contains much different information about a given infinite HE 2 5 graph system (G, c) than does the more familiar ℓ sequence space, even when appropriate 6 weights are assigned. In the language of Markov processes, is better adapted to the 2 HE 7 study of (G, c) than ℓ . One reason for this is that is intimately connected with the HE 8 resistance metric R.
9 10 — and D. The dissipation space D is the Hilbert space of functions on the § HE H H 10 edges when equipped with the dissipation inner product. We solve problems in discrete 11 potential theory with the use of the drop operator d (and its adjoint d∗), where
dv(x, y) := c (v(x) v(y)). (1.13) xy − 12 The drop operator d is, of course, just an implementation of Ohm’s Law, and can be in- 13 terpreted as a weighted boundary operator in the sense of homology theory. The drop 14 operator appears elsewhere in the literature, sometimes without the weighting cxy; see 15 [Chu01, Tel06, Woe00]. However, we use the adjoint of this operator with respect to the 2 16 energy inner product, instead of the ℓ inner product used by others. This approach appears 17 to be new, and it turns out to be more compatible with physical interpretation. For exam- 2 18 ple, the displayed equation preceding [Woe00, (2.2)] shows that the ℓ adjoint of the drop 19 operator is incompatible with Kirchhoff’s node law. Since the resistance metric may be 20 defined in terms of currents obeying Kirchhoff’s laws, we elect to make this break with the 21 existing literature. Additionally, this strategy will allow us to solve the compatibility prob- 22 lem described in 3.3 in terms of a useful minimizing projection operator P , discussed § d 23 in detail in 10.4. Furthermore, we believe our formulation is more closely related to the § 24 (co)homology of the ERN as a result. 25 We decompose into the direct sum of the range of d and the currents which are HD 26 sums of characteristic functions of cycles
= ran d [χ ], (1.14) HD ⊕ ϑ 27 where ϑ is a cycle, i.e., a path in the graph which ends where it begins. In(1.14) and else- 28 where, we indicate the closed linear span of a set by [χ ] := cl span χ . From (1.8) (and ϑ { ϑ} 29 the fact that d is an isometry), it is clear that the first summand of (1.14) can be further 30 decomposed into weighted edge neighbourhoods dδx and the image of harmonic functions 31 under d in Theorem 10.8. After a first draught of this paper was complete, we discovered 32 that the same approach is taken in [LP09]. One of us (PJ) recalls conversations with Raul 33 Bott concerning an analogous Hilbert space operator theoretic approach to electrical net- 34 works; apparently attempted in the 1950s in the engineering literature. We could not find 35 details in any journals; the closest we could come is the fascinating paper [BD49] by Bott 36 et al. A further early source of influence is Norbert Wiener’s paper [WR46]. 37 In 10.4, we describe how d solves the compatibility problem and may be used to solve § ∗ 38 a large class of problems in discrete potential theory. Also, we discuss the analogy with 39 complex analysis.
40 11 — Probabilistic interpretations. In [LP09,DS84,Tel06,Woe00] and elsewhere, the § 41 randomwalk on an ERN is defined by the transition probabilities p(x, y) := cxy/c(x). In this 1 42 context, the probabilistic transition operator is P = c− T and one uses the stochastically 1 43 renormalized Laplacian ∆c := c− ∆, where c is understood as a multiplication operator; 12 PALLEE.T.JORGENSENANDERINP.J.PEARSE see Definition 2.3. This approach also arises in the discussion of trace resistance in 5.2.3 1 § and allows one to construct currents on the graph as the average motion of a random walk. 2 As an alternative to the approach described above, we discuss a probabilistic interpre- 3 tation slightly different from those typically found in the literature: we begin with a volt- 4 age potential as an initial condition, and consider the induced current I. The components 5 of such a flow are called current paths and provide a way to interpret potential-theoretic 6 problems in a probabilistic setting. We study the random walks where the transition prob- 7 ability is given by I(x, y)/ z x I(x, z). We consider the harmonic functions in this con- 8 ∼ text, which we call forward-harmonicP functions, and the associated forward-Laplacian of 9 Definition 11.19. We give a complete characterization of forward-harmonic functions as 10 cocycles, following [Jor06]. 11
12 — Examples. We collect an array of examples that illustrate the various phenomena 12 § encountered in the theory and work out many concrete examples. Some elementary finite 13 examples are given in 12.1 to give the reader an idea of the basics of electrical resistance 14 § network theory. In 12.2 we move on to infinite graphs. 15 § 13 — Trees. When the electrical resistance network has a tree as an underlying net- 16 § work, the resistance distance coincides with shortest-path metric, as there is always exactly 17 one path between any two vertices; cf. Lemma 5.46 and the preceding discussion. When 18 the tree has exponential growth, as in the case of homogeneous trees of degree 3, one 19 ≥ can always construct nontrivial harmonic functions, and monopoles of finite energy. In 20 fact, there is a very rich family of each, and this property makes this class of examples a 21 fertile testing ground for many of our theorems and definitions. In particular, these exam- 22 ples highlight the relevance and distinctions between the boundary (as we construct it), the 23 Cauchy completion, and the graph ends of [PW90, Woe00]. In particular, they enable one 24 to see how adjusting decay conditions on c affects these things. 25
d 14 — Integer lattices. The lattice electrical resistance networks (Z , c) have vertices at 26 § d the points of R which have integer coordinates, and edges between every pair of vertices 27 (x, y) with x y = 1. The case for c = 1 is amenable to Fourier analysis, and in 14.1 we 28 | − | § obtain explicit formulas for many expressions: 29
Lemma 14.4 gives a formula for the potential configuration functions v . 30 • { x} Theorem 14.7 gives a formula for the resistance distance R(x, y). 31 • Theorem 14.9 gives a formula for the resistance distance to infinity in the sense 32 • R(x, ) = limy R(x, y). 33 ∞ →∞ d Theorem 14.5 gives a formula for the solution w of ∆w = δ on Z ; it is readily 34 • − o seen that this w has finite energy (i.e., is a monopole) iff d 3. 35 ≥ In [P´ol21], P´olya proved that the random walk on this graph is transient if and only if 36 d 3; see [DS84] for a nice exposition. We offer a new characterization of this dichotomy 37 ≥ d (there exist monopoles on Z if and only if d 3) which we recover in this section via a 38 ≥ new (and completely constructive) proof. In Remark 14.21 we describe how in the infinite 39 integer lattices, functions in may be approximated by functions of finite support. 40 HE d 15 — Magnetism. The integer lattice networks (Z , 1) investigated in 14 comprise the 41 § § framework of infinite models in thermodynamics and in quantum statistical mechanics. In 42 15.3 we employ these formulas in the refinement of an application to the theory of the 43 § (isotropic Heisenberg) model of ferromagnetism as studied by R. T. Powers. In addition 44 to providing an encapsulated version of the Heisenberg model, we give a commutative 45 analogue of the model, extend certain results of Powers from [Pow75, Pow76a, Pow76b, 46 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 13
1 Pow78, Pow79], and discuss the application of the resistance metric to the theory of ferro- 2 magnetism and “long-range order”. This problem was raised initially by R. T. Powers, and 3 may be viewed as a noncommutative version of Hilbert spaces of random variables. 4 Ferromagnetism in quantum statistical mechanics involves algebras of noncommutative 5 observables and may be described with the use of states on C∗-algebras. As outlined in 6 the cited references, the motivation for these models draw on thermodynamics; hence the 7 notions of equilibrium states (formalized as KMS states, see 15.4). These KMS states § 8 are states in the C∗-algebraic sense (that is, positive linear functionals with norm 1), and 9 they are indexed by absolute temperature. Physicists interpret such objects as representing 10 equilibria of infinite systems. 11 In the present case, we consider spin observables arranged in a lattice of a certain rank, 12 d = 1, 2, 3,..., and with nearest-neighbor interaction. Rigourous mathematical formula- 13 tion of phase transitions appears to be a hopeless task with current mathematical technol- 14 ogy. As an alternative avenue of enquiry, much work has been conducted on the issue of 15 long-range order, i.e., the correlations between observables at distant lattice points. These 16 correlations are measured relative to states on the C∗-algebra; in this case in the KMS states 17 for a fixed value of temperature. 18 While we shall refer to the literature, e.g. [BR79, Rue69] for formal definitions of key 19 terms from the C∗-algebraic formalism of quantum spin models, physics, and KMS states, 20 we include a minimal amount of background and terminology from the physics literature.
21 16 — Future Directions. We conclude with a brief discussion of several projects which § 22 have arisen from work on the present paper, as well as some promising new directions that 23 we have not yet had time to pursue.
24 Appendices. We give some background material from functional analysis in Appen- 25 dix A, and operator theory in Appendix B. In Appendix C, we include some diagrams to 26 help clarify the properties of the many operators and spaces we discuss, and the relations 27 between them.
28 1.2. What this paper is about. Theeffective resistance metric provides the foundation for 29 our investigations because it is the natural and intrinsic metric for an electrical resistance 30 network, as the work of Kigami has shown; see [Kig01] and the extensive list of references 31 by the same author therein. Moreover, the close relationship between diffusion geometry + 32 (i.e., geometry of the resistance metric) [MM08, SMC08, CKL 08, CM06] and random 33 walks on graphs leads us to expect/hope there will be many applications of our results to 34 several other subjects, in addition to fractals: models in quantum statistical mechanics, 35 analysis of energy forms, interplay between self-similar measures and associated energy 36 forms, certain discrete models arising in the study of quasicrystals (e.g., [BM00, BM01]), 37 and multiwavelets (e.g. [BJMP05, DJ06, DJ07, Jor06]), among others. A general theme 38 of these areas is that the underlying space is not sufficiently regular to support a group 39 structure, yet is “locally” regular enough to allow analysis via probabilistic techniques. 40 Consequently, the analysis of functions on such spaces is closely tied to Dirichlet energy 41 forms and the graph Laplacian operator associated to the graph. This appears prominently 42 in the context of the present paper as follows:
43 (1) The embedding of the metric space ((G, c), R) into the Hilbert space of func- HE 44 tions of finite energy,in such a way that the original metric may be recovered from 45 the norm, i.e., 14 PALLEE.T.JORGENSENANDERINP.J.PEARSE
2 R(x, y) = (vx vy) = vx vy , E − k − kE where vx is the image of x under the embedding. 1 ∈HE (2) The relation of the energy form to the graph Laplacian via the equation 2
(u, v) = u(x)∆v(x) + u(x) ∂v (x), (1.15) E ∂Ò xXG0 xXbd G ∈ ∈ introduced just above in the discussion of 2. Each summation on the right hand 3 § side of (1.15) is more subtle than it appears. These details for the first sum are 4 given in Theorem 4.41, and the details for the second sum are the focus of almost 5 all of 7. 6 § (3) The presence of nonconstant harmonic functions of finite energy. These are pre- 7 W cisely the objects which support the boundary term in (1.15) and imply R (x, y) < 8 F R (x, y). They are also responsible for the boundary described in 7. 9 § (4) The solvability of the Dirichlet problem ∆w = δy, where δy is a Dirac mass 10 0 − at the vertex y G . The existence of finite-energy solutions w is equivalent 11 ∈ to the transience of the random walk on the network. Such functions are called 12 monopoles and (via Ohm’s law) they induce a unit flow to infinity as discussed 13 in [DS84, LP09, LPW08]. 14
Remark 1.1. In addition to uses in graph theory and electrical networks, the discrete Lapla- 15 cian ∆ has other uses in numerical analysis: many problemsin PDE theory lend themselves 16 to discretizations in terms of subdivisions or grids of refinements in continuous domains. 17 A key tool in applying numerical analysis to solving partial differential equations is dis- 18 cretization, and use of repeated differences; especially for using the discrete ∆ in approxi- 19 mating differential operators, and PDOs. See e.g., [AH05]. 20 One picks a grid size δ and then proceeds in steps: 21
(1) Start with a partial differential operator, then study an associated discretized oper- 22 d ator with the use of repeated differences on the δ-lattice in R . 23 (2) Solve the discretized problem for h fixed. 24 (3) As δ tends to zero, numerical analysts evaluate the resulting approximation limits, 25 and they bound the error terms. 26
When discretization is applied to the Laplace operator in d continuous variables, the 27 d result is our ∆ for the network (Z , c); see 14 for details and examples. However, when 28 § the same procedure is applied to a continuous Laplace operator on a Riemannian manifold, 29 the discretized ∆ will be the network Laplacian on a suitable infinite network (G, c) which 30 d in general may have a much wilder geometry than Z . 31 This yields numerical algorithms for the solution of partial differential equations, and 32 in the case of second order PDEs, the discretized operator is the discrete Laplacian studied 33 in this investigation. 34
1.2.1. Motivation and applications. Applications to infinite networks of resistors serve as 35 motivations, but our theorems have a wider scope, have other applications; and are, we 36 believe, of independent mathematical interest. Our interest originates primarily from two 37 sources. 38
(1) A series of paperswritten by Bob Powers in the 1970swhich he introduced infinite 39 systems of resistors into the resolution of an important question from quantum 40 statistical mechanics in [Pow75, Pow76a, Pow76b, Pow78, Pow79]. 41 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 15
1 (2) The pioneering work of Jun Kigami on the analysis of PCF self-similar fractals, 2 viewing these objects as rescaled limits of networks; see [Kig01].
3 Indeed, our larger goal is the cross-pollination of these two areas, and we hope that the 4 results of this paper may be applicable to analysis on fractal spaces. A first step in this 5 direction is given in Theorem 16.3. To this end, a little more discussion of each of the 6 above two subjects is in order. 7 Powers was interested in magnetism and the appearance of “long-range order”, which 8 is the common parlance for correlation between spins of distant particles; see 15 for a § 9 larger discussion. Consequently, he was most interested in graphs like the integer lattice d 10 Z (with edges between vertices of distance 1, and all resistances equal to 1), or other 11 regular graphs that might model the atoms in a solid. Powers established a formulation of 12 resistance metric that we adopt and extend in 5, where we also show it to be equivalent to § 13 Kigami’s formulation(s). Also, the proofs of Powers’ originalresults on effective resistance 14 metric contain a couple of gaps that we fill. In particular, Powers does not seem to have be 15 aware of the possibility of nontrivial harmonic functions until [Pow78], where he mentions 16 them for the first time. It is clear that he realized several immediate implications of the 17 existence of such functions, but there more subtle (and just as important!) phenomena that 18 are difficult to see without the clarity provided by Hilbert space geometry. 19 Powers studied an infinite graph G by working with an exhaustion, that is, a nested se- 20 quence of finite graphs G1 G2 ... Gk G = k Gk. For example, Gk might be all d ⊆ ⊆ ⊆ ⊆ 21 the vertices of Z lying inside the ball of Euclidean radiusS k, and the edges between them. 22 Powers used this approach to obtain certain inequalities for the resistance metric, express- 23 ing the consequences of deleting small subsets of edges from the network. Although he 24 makes no reference to it, this approach is very analogous to Rayleigh’s “short-cut” meth- 25 ods, as it is called in [DS84]. 26 Powers’ use of an expanding sequence of graphs may be thought of as a “limit in the 27 large” in contrast to the techniques introduces by Kigami, which may be considered “limits 28 in the small”. Self-similarity and scale renormalization are the hallmarks of the theory of 29 fractal analysis as pursued by Kigami, Strichartz and others (see [HKK02, Kig01, Kig03, 30 Hut81, Str06, BHS05, Bea91, JP94, Jor04], for example) but these ideas do not enter into 31 Powers’ study of resistors. One aim of the present work is the development of a Hilbert 32 space framework suitable for the study of limits of networks defined by a recursive al- 33 gorithm which introduces new vertices at each step and rescales the edges via a suitable 34 contractive rescaling. As is known from, for example [JP94, Jor06, Str98a, Str06, Tep98], 35 there is a spectral duality between “fractals in the large”, and “fractals in the small”.
36 1.2.2. The significance of Hilbert spaces. A main theme in this paper is the use of Hilbert 37 space technology in understanding metrics, potential theory, and optimization on infinite 38 graphs, especially through finite-dimensional approximation. We emphasize those aspects 39 that are intrinsic to infinite electrical resistance networks, and our focus is on analytic as- 40 pects of graphs; as opposed to the combinatorial and algebraic sides of the subject, etc. 41 Those of our results stated directly in the framework of graphs may be viewed as discrete 42 analysis, yet the continuum enters via spectral theory for operators and the computation of 43 probability of sets of infinite paths. In fact, we will display a rich variety of possible spec- 44 tral types, considering the spectrum as a set (with multiplicities), as well as the associated 45 spectral measures, and representations/resolutions. 46 Related issues for Hilbert space completions form a recurrent theme throughout our pa- 47 per. Given an electrical resistance network, we primarily study three spaces of functions 16 PALLEE.T.JORGENSENANDERINP.J.PEARSE
2 0 naturally associated with it: , D, and to a lesser extent ℓ (G ). Our harmonic anal- 1 HE H ysis of functions on G is studied via operators between the respective Hilbert spaces as 2 discussed in 10 and the Hilbert space completions of these three classes are used in an 3 § essential way. In particular, we obtain the boundary of the graph (a necessary ingredient 4 of (1.15) and the key to several mysteries) by analyzing the finite energy functions on G 5 which cannot be approximated by functions of finite support. However, this metric space 6 is naturally embedded inside the Hilbert space , which is already complete by defini- 7 HE tion/construction. Consequently, the Hilbert space framework allows us to identify certain 8 vectors as corresponding to the boundary of (G, c), and thus obtain a concrete understand- 9 ing of the boundary. 10 However, the explicit representations of vectors in a Hilbert space completion (i.e., the 11 completion of a pre-Hilbert space) may be less than transparent; see [Yoo07]. In fact, 12 this difficulty is quite typical when Hilbert space completions are used in mathematical 13 physics problems. For example, in [JO00´ ,Jor00], one beginswith a certain space of smooth 14 d functions defined on a subset of R , with certain support restrictions. In relativistic physics, 15 one must deal with reflections, and there will be a separate positive definite quadratic 16 form on each side of the “mirror”. As a result, one ends up with two startlingly different 17 2 Hilbert space completions: a familiar L -space of functions on one side, and a space of 18 distributions on the other. In [JO00´ ,Jor00], one obtains holomorphic functions on one side 19 of the mirror, and the space of distributions on the other side is spanned by the derivatives 20 of the Dirac mass, each taken at the same specific point x0. 21
1.3. What this paper is not about. Many of the topics discussed in this paper may appear 22 to have been previously discussed elsewhere in the literature, but there are certain important 23 subtleties which actually make our results quite different. This section is intended to clarify 24 some of these. It is the opinionof the authorsthat most interesting results of this paper arise 25 primarily from three things: 26
(1) differences between finite approximations to infinite networks, and how & when 27 these differences vanish in the limit, and 28 (2) the phenomena that result when one works with a quadratic form whose kernel 29 contains the constant functions, and 30 (3) the boundary (which is not a subset of the vertices) that naturally arises when 31 a network supports nonconstant harmonic functions of finite energy, and how it 32 explains other topics mentioned above. 33
In classical potential theory, working modulo constant functions amounts to working with 34 2 the class of functions satisfying f < , but abandoning the ℓ requirement f < . 35 k ′k2 ∞ k k2 ∞ This has some interesting consequences, and the nontrivial harmonic functions play an 36 especially important role; see Remark 4.66. What would one hope to gain by removing the 37 2 ℓ condition? 38
(1) From the natural embedding of the metric space (G, R) into the Hilbert space 39 HE of functions of finite energy given by x vx, the functions vx are not generally in 40 2 7→ ℓ . See Figure 10 of Example 14.16 for an illustration. 41 2 (2) The resistance metric does not behave nicely with respect to ℓ conditions. Several 42 formulations of the resistance distance R(x, y) involve optimizing over collections 43 2 of functions which are not necessarily contained in ℓ , even for many simple ex- 44 amples. 45 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 17
2 0 1 (3) Corollary 4.65 states that nontrivial harmonic functions cannot lie in ℓ (G ). Con- 2 2 sequently, imposing an ℓ hypothesis removes the most interesting phenomena 3 from the scope of study; see Remark 4.66.
4 The infinite trees studied in Examples 13.2–13.6 also provide examples of these situations.
5 While there already is a large literature on electrical networks and on graphs (see 6 e.g., [CW92, CW07, DK88, Dod06, DS84, Pow76b, CdV04, CR06, Chu07, FK07], and the 7 preprint [Str08] which we received after the first version of this paper was completed), we 8 believe that our present operator/spectral theoretic approach suggests new questions and 9 new theorems, and allows many problems to be solved in greater generality. 10 The literature on analysis on graphs breaks down into a variety of overlapping subar- 11 eas, including: combinatorial aspects, systems of resistors on infinite networks, random- 12 walk models, operator algebraic models [DJ08,Rae05], probability on graphs (e.g., infinite 13 particle models in physics [Pow79]), Brownian motion on self-similar fractals [Hut81], 14 Laplace operators on graphs, finite element-approximations in numerical analysis [BS08]; 15 and more recently, use in internet-search algorithms [FK07]. Even just the study of Laplace 16 operators on graphs subdivides further, due to recently discovered connections between 17 graphsand fractals generated by an iterated functionssystem (IFS); see e.g.,[Kig03,Str06]. 18 Other major related areas include discrete Schr¨odinger operators in physics, informa- 19 tion theory, potential theory, uses of the graphs in scaling-analysis of fractals (constructed 20 from infinite graphs), probability and heat equations on infinite graphs, graph C∗-algebras, 21 groupoids, Perron-Frobenius transfer operators (especially as used in models for the in- 22 ternet); multiscale analysis, renormalization, and operator theory of boundaries of infinite 23 graphs (more current and joint research between the co-authors.) The motivating appli- 24 cations from [Pow75, Pow76a, Pow76b, Pow78, Pow79] include the operator algebra of 25 electrical networks of resistors (lattice models, C∗-algebras, and their representations), and 26 more specifically, KMS-states from statistical mechanics. While working and presenting 27 our results, we learned of even more such related research directions from experts working 28 in these fields, and we are thankful to them all for taking the time to explain some aspects 29 of them to us. 30 The main point here is that the related literature is vast but our approach appears to be 31 entirely novel and our results, while reminiscent of classical theory, are also new. We now 32 elucidate certain specific differences.
33 1.3.1. Spectral theory. The spectral theory for networks contrasts sharply with that for 34 fractals, as is seen by considering the measures involved; they do not begin to become 35 similar until one considers limits of networks. The spectrum of discrete Laplacians on in- 36 finite networks is typically continuous (lattices or trees provide examples, and are worked 37 explicitly in 12). By contrast, in the analysis on fractals program of Kigami, Strichartz, § 38 and others, the Laplace operator has pure point spectrum; see [Tep98] in particular. The 39 measures used in the analysis of networks are weighted counting measures, while the mea- 40 sures used in fractal analysis are based on the self-similar measures introduced by Hutchin- 41 son [Hut81]. There is an associated and analogous entropy measure in the study of Julia 42 sets; cf. [Bea91] and the recent work on Laplacians in [RT08]. 43 Our approach differs from the extensive literature on spectral graph theory (see [Chu01] 44 for an excellent introduction, and an extensive list of further references) due to the fact 2 45 that we eschew the ℓ basis for our investigations. We primarily study ∆ as an operator 46 on , and with respect to the energy inner product. The corresponding spectral theory is HE 2 47 radically different from the spectral theory of ∆ in ℓ . Most other work in spectral graph 18 PALLEE.T.JORGENSENANDERINP.J.PEARSE
2 theory takes place in ℓ , even implicitly when working with finite graphs: the adjoint of 1 2 the drop operator (see Definition 10.2) is taken with respect to the ℓ inner product and 2 consequently violates Kirchhoff’s laws. In fact, the discussion preceding [Woe00, (2.2)] 3 shows how this version of the adjoint is incompatible with Kirhhoff’s Law as mentioned in 4 the summary of 10 just above. Additionally, [Chu01] and others work with the spectrally 5 § 1/2 1/2 renormalized Laplacian ∆s := c− ∆c− . However, ∆s is a bounded Hermitian operator 6 (with spectrum contained in [0, 2]) and so is unsuitable for our investigationsof bdG based 7 on defect indices, etc. 8 As we have only encountered relatively few instances where the complete details are 9 worked out for spectral representations in the framework of discrete analysis, we have at- 10 tempted to provide several explicit examples. These are likely folkloric, as the geometric 11 possibilities of graphs are vast, and so is the associated range of spectral configurations. A 12 list of recent and past papers of relevance includes [Str08, Car72, Car73a, Car73b, CR06, 13 Chu07,CdV99,CdV04,Jor83], and Wigner’s original paper on the semicircle law [Wig55]. 14 The present investigation also led to a spectral analysis of the binary tree from the perspec- 15 tive of dipoles in [DJ08]; this study discovered that the spectrum of ∆ on the binary tree is 16 also given by Wigner’s semicircle law. 17 There is also a literature on infinite/transfinite networks and generalized Kirchhoff laws 18 using nonstandard analysis, etc., see [Zem91, Zem97]. However, this context allows for 19 edges with resistance 0, which we do not allow (for physical as well as theoretical rea- 20 sons). One can neglect the resistance of wires in most engineering applications, but not 21 when considering infinite networks (the epsilons add up!). The resulting theory therefore 22 divergesrapidly from the observationsof the present paper; according to our definitions, all 23 networks support currents satisfying Kirchhoff’s law, and in particular, all induced currents 24 satisfy Kirchhoff’s law. 25
1.3.2. Operator algebras. There are also recent papers in the literature which also ex- 26 amine graphs with tools from operator algebras and infinite determinants. The papers 27 [GILb, GILc, GILa] by Guido et al are motivated by questions for fractals and study the 28 detection of periods in infinite graphs with the use of the Ihara zeta function, a variant of 29 the Riemann zeta function. There are also related papers with applications to the operator 30 algebra of groupoids [Cho08,FMY05], and the papers [BM00,BM01] which apply infinite 31 graphs to the study of quasi-periodicity in solid state physics. However, the focus in these 32 papers is quite different from ours, as are the questions asked and the methods employed. 33 While periods and quasi-periods in graphs play a role in our present results, they enter 34 our picture in quite different ways, for example via spectra and metrics that we compute 35 from energy forms and associated Laplace operators. There does not seem to be a direct 36 comparison between our results and those of Guido et al. 37
1.3.3. Boundaries of graphs. There is also no shortage of papers studying boundaries of 38 infinite graphs: [PW90, Saw97, Woe00] discuss the Martin boundary, [PW90, Woe00] also 39 describe the more geometrically constructed “graph ends”, and [Car72, Car73a, Car73b] 40 use unitary representations. There are also related resultsin[CdV99,CdV04] While there 41 are connections to our study, the scope is different. 42 Martin boundary theory is really motivated by constructing a boundary for a Markov 43 process, and the geometry/topology of the boundary is rather abstract and a bit nebulous. 44 Additionally, one needs a Green’s function, and it must satisfy certain hypotheses before 45 the construction can proceed. Furthermore, the focus of Martin boundarytheory is the non- 46 negative harmonic functions. Our boundary construction is more general in that it applies 47 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 19
1 to any electrical network as in Definition 2.7 and it remains correct for all harmonic func- 2 tions of finite energy, including constant functions and harmonic functions which change 3 sign. However, it is also more restrictive in the sense that an electrical resistance network 4 may support functions which are bounded below but do not have finite energy. 5 We should also point out that our boundaryconstruction is related to, but different from, 6 the “graph ends” introduced by Freudenthal and others. The ends of a graph are the natural 7 discrete analogue of the ends of a minimal surface (usually assumed to be embedded in 3 8 R ), a notion which is closely related to the conformaltype of the surface. Starting with the 9 central book [Woe00] by Wolfgang Woess, the following references will provide the reader 10 with an introductionto the study of harmonic functions on infinite networks and the ends of 11 graphs and groups: [Woe86, Woe87, Woe89], and [Woe95] on Martin boundaries, [PW90] 12 on ends, [Woe96] on Dirichlet problems, [Woe97] on random walk. A comparison of the 13 examples in 14 and 13 illustrates that varying the resistances produces dramatic changes § § 14 in the topology of the boundary. 15 Our boundary essentially consists of infinite paths which can be distinguished by har- 16 monic functions; see 7.3 for details. It follows that transient networks with no nontrivial § 17 harmonic functions have exactly one boundary point (corresponding to the unique mono- d 18 pole). In particular, the integer lattices (Z , 1) have precisely 1 boundary point for d 3, 2 ≥ 19 and have 0 boundary points for d = 1, 2. The Martin boundary of (Z , 1) consists of two 2 20 points; similarly, (Z , 1) has two graph ends; cf. [PW90].
21 1.3.4. Measures and measure constructions. A reader glancing at our paper will notice a 22 number of incarnations of measures on infinite sample spaces: it may be a suitable space 23 of paths ( 11.1– 11.2 and 15.1) or an analogue of the Schwartz space of tempered distri- § § § 24 butions (section 7.2). The latter case relies on a construction of “Gel’fand triples” from § 25 mathematical physics. The reader may wonder why they face yet another measure con- 26 struction, but each construction is dictated by the problems we solve. Taking limits of finite 27 subsystems is a universal weapon used with great success in a variety of applications; we 28 use it here in the study of resistance distances on infinite graphs ( 5.2); boundaries, bound- § 29 ary representations for harmonic functions ( 7.2, 8.3, and 11.1– 11.2); and equilibrium § § § § 30 states and phase-transition problems in physics ( 15.1– 15.2). 2 § § 31 (1) as an L space. The central Hilbert space in this study, the energy space , HE HE 32 appears with a canonical reproducing kernel, but without any canonical basis, and there is 2 33 no obvious way to see as an L (X, µ) for some X and µ. Therefore, a major motivation HE 2 34 for our measure constructions is just to be able to work with as an L space. In 15.1, 2 HE § 35 we use a construction from probability to write = L (Ω, µ) in a way that makes the HE 36 energy kernel vx x G0 into a system of (commuting) random variables. Here, Ω is an { } ∈ 37 infinite Cartesian product of a chosen compact space S; one copy of S for each point 0 38 x G . In 15.2, we use a non-commutative version of this probability technology: rather ∈ § 39 than Cartesian products, we will use infinite tensor products of C∗-algebras , one for 0 A 40 each x G . The motivation here is an application to a problem in quantum statistical ∈ 41 mechanics. The “states” on the C∗-algebra of all observables are the quantum mechanical 42 analogues of probability measures in classical problems. Heuristically, the reader may 43 wish to think of them as non-commutative measures; see e.g., [BR97]. 44 (2) Boundary integral representation of harmonic functions. As it sometimes happens, 45 the path to bd G is somewhat circuitous: we begin with the discovery of an integral over the 46 boundary,which leads us to understand functions on the boundary, which in turn points the 47 way to a proper definition of the boundaryitself. A closely related motivation for a measure 48 is the formulation of an integral representation of harmonic functions u : ∈HE 20 PALLEE.T.JORGENSENANDERINP.J.PEARSE
u(x) = u(ξ)hx(ξ) dP(ξ) + u(o). (1.16) Z SG′ where hx = P armvx. Thus the focus of 7.2 is a formalization of the imprecise “Riemann 1 H ∂hx § sums” u(x) = u + u(o) of 4.3 as an integral of a bona fide measure. To carry this 2 bd G ∂Ò § out, we constructP a Gel’fand triple G G′ , where G is a dense subspace of and 3 S ⊆HE ⊆ S S HE G’ is its dual, but with respect to a strictly finer topology. We are then able to produce a 4 S 2 Gaussian probability measure P on G′ and isometrically embed into L ( G′ , P). In fact, 5 2 S HE S L ( G′ , P) is the second quantization of . However, the focus here is not on realizing 6 S 2 HE as an L space (or subspace), but in obtaining the boundary integral representation of 7 HE harmonic functions as in (1.16). Our aim is then to build formulasthat allow us to compute 8 values of harmonic functions u from an integral representation which yields u(x) as 9 ∈HE an integral over bd G ′ . Note that this integration in (1.16) is with respect to a measure 10 ⊆ SG depending on x just as in the Poisson and Martin representations. 11 (3) Concrete representation of the boundary. We would like to realize bd G as a measure 12 space defined on a set of well-understood elements; this is the focus of the constructions 13 in 7. The goal is a measure on the space of all infinite paths in G which yields the 14 § boundary bd G in such a way that G bd G is a compactification of G which is compatible 15 ∪ with the energy form and the Laplace operator ∆, and hence also the natural resistance 16 E metric on (G, c). This type of construction has been carried out with great success for the 17 case of bounded harmonic functions (e.g., Poisson representation and the Fatou-Primalov 18 theorem) and for nonnegative harmonic functions (e.g., Martin boundary theory), but our 19 scope of enquiry is the harmonic functions of finite energy. Finally, we would like to use 20 this Gaussian measure on to clarify bd G as a subspace of . Such a relationship is 21 SG′ SG′ a natural expectation, as the analogous thing occurs in the work of Poisson, Choquet, and 22 Martin. 23
1.4. General remarks. 24
Remark 1.2. Since we aim for several different audiences (operator theory, analysis of 25 fractals, mathematical physics, etc), we have included more details in our exposition and 26 proofs than would otherwise be typical for a paper with a narrow focus. 27
Remark 1.3. Throughout the introductory discussion of electrical resistance networks in 28 2– 5, we discuss collections of real-valued functions on the vertices or edges of the graph 29 § § G. Such objects are most natural for the heuristics of the physical model, and addition- 30 ally allow for induced orientation/order and make certain probabilistic arguments possible. 31 However, in the latter portions of this paper, we need to incorporate complex-valued func- 32 tions into the discussion in order to make full use of spectral theory and other methods. 33
Remark 1.4. For the aid of the reader, we have included a list of symbolsand abbreviations 34 used in this document. Wherever possible, we have attempted to ensure that each symbol 35 has only one meaning. In cases of overlap, the context should make things clear. In 36 Appendix C, we also include some diagrams which we hope clarify the properties of the 37 many operators and spaces we discuss, and the relations between them. 38
1.5. Acknowledgements. 39 While working on the project, the co-authors have benefitted from interaction with col- 40 leagues and students. We thank everyone for generously suggesting improvements as our 41 paper progressed. The authors are grateful for stimulating comments, helpful advice, and 42 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 21
1 valuable references from John Benedetto, Il-Woo Cho, Raul Curto, Dorin Dutkay, Alexan- 2 der Grigor’yan, Dirk Hundertmark, Richard Kadison, Keri Kornelson, Michel Lapidus, 3 Russell Lyons, Diego Moreira, Peter M¨orters, Paul Muhly, Bob Powers, Marc Rieffel, 4 Karen Shuman, Sergei Silvestrov, Jon Simon, Myung-Sin Song, Bob Strichartz, Andras 5 Telcs, Sasha Teplyaev, Ivan Veselic, Lihe Wang, Wolfgang Woess, and Qi Zhang. The 6 authors are particularly grateful to Russell Lyons for several key references and examples, 7 and to Jun Kigami for several illuminating conversations and for suggesting the approach 8 in (5.41). Initially, the first named author (PJ) learned of discrete potential theory from 9 Robert T. Powers at the University of Pennsylvania in the 1970s, but interest in the subject 10 has grown exponentially since. 22 PALLEE.T.JORGENSENANDERINP.J.PEARSE
2. Electrical resistance networks 1
In this section, we introduce the mathematical model of an electrical resistance network 2 (ERN) as a graph G whose edges are understood as resistors and whose vertices are the 3 nodes at which these resistors are connected. Thus, the resistance data is specified by a 4 function Ω, so that Ω(x, y) is the resistance of the edge (resistor) between the vertices x 5 and y. With the network data (G, Ω) fixed, we begin the study of functions defined on the 6 vertices. We define many basic terms and conceptsused throughoutthe paper, including the 7 Dirichlet energy form and the Laplace operator ∆. Additionally, we prove a key identity 8 E relating to ∆ for finite graphs: Lemma 2.13. In Theorem 4.41, this will be extended to 9 E infinite graphs, in which case it is a discrete analogue of the familiar Gauss-Green identity 10 from vector calculus. The appearance of a somewhat mysterious boundary term in the 11 Theorem 4.41 prompts several questions which are discussed in Remark 4.7. Answering 12 these questions comprises a large part of the sequel; cf. 7. In fact, Theorem 4.41 provides 13 § much of the motivation for energy-centric approach we pursue throughout our study; the 14 reader may wish to look ahead to Remark 4.66 for a preview. 15
2.1. The electrical resistance network model. This section contains the basic definitions 16 used throughout the sequel. 17
0 1 0 Definition 2.1. A graph G = G , G is given by the set of vertices G and the set of edges 18 1 0 0 { } 1 G G G . Two vertices are neighbours (or are adjacent) iff there is an edge (x, y) G 19 ⊆ × ∈ 1 connecting them, and this is denoted x y. This relation is symmetric, as (y, x) G 20 1 ∼ 0 ∈ whenever (x, y) G . The set of neighboursof x G is 21 ∈ ∈ 0 . G(x) = y G .. y x . (2.1) { ∈ ∼ }
In our context, the set of edges of G will be determined by the conductance function, so 22 that all graph data is implicitly provided by c. 23
Definition 2.2. The conductance cxy is a symmetric function 24
c : G0 G0 [0, ), (2.2) × → ∞ in the sense that cxy = cyx. It is our convention that x / y if and only if cxy = 0; that is, 25 1 ∼ there is an edge (x, y) G if and only if 0 < c(x, y) < . 26 ∈ ∞ Conductance is the reciprocal of resistance, and this is the origin of the name “resistance 27 1 network”. It is important to note that c−xy gives the resistance between adjacent vertices; 28 1 this feature distinguishes c−xy from the effective resistance R(x, y) discussed later, for which 29 x and y need not be adjacent. 30
0 Definition 2.3. The conductances define a measure or weighting on G by 31
c(x) := cxy. (2.3) Xy x ∼ 0 Whenever G is connected, it follows that c(x) > 0, for all x G . The notation c will also 32 ∈ be used, on occasion, to indicate the multiplication operator (cv)(x) := c(x)v(x). 33
0 0 Definition 2.4. A path γ from α G to ω G is a sequence of adjacent vertices 34 ∈ ∈ (α = x0, x1, x2,..., xn = ω), i.e., xi xi 1 for i = 1,..., n. The path is simple if any vertex 35 ∼ − appears at most once (so that a path is simply connected). 36 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 23
0 1 Definition 2.5. A graph G is connected iff for any pair of vertices α, ω G , there exists ∈ 2 a finite path γ from α to ω.
3 Remark 2.6. Note that for resistors connected in series, the resistances just add, so this 4 condition implies there is a path of finite resistance between any two points. We emphasize 5 that all graphs and subgraphs considered in this study are connected.
6 At this point, the reader may wish to peruse some of the examples of 12. § 7 Definition 2.7. An electrical resistance network (ERN) is a connected graph (G, c) whose 0 8 conductance function satisfies c(x) < for every x G . We interpret the edges as being ∞ ∈ 9 defined by the conductance: x y iff c > 0. ∼ xy 10 Note that c need not be bounded in Definition 2.7. Also, we will typically assume an 11 ERN to be simple in the sense that there are no self-loops, and there is at most one edge 12 from x to y. This is mostly for convenience: basic electrical theory says that two conductors 1 2 13 cxy and cxy connected in parallel can be replaced by a single conductor with conductance 1 2 14 cxy = cxy + cxy. Also, electric current will never flow along a conductor connecting a node 15 to itself. Nonetheless, such self-loops may be useful for technical considerations: one can 16 remove the periodicity of a random walk by allowing self-loops. This can allow one to 17 obtain a “lazy walk” which is ergodic, and hence amenable to application of tools like the 18 Perron-Frobenius Theorem. See, for example, [LPW08, LP09]. 19 We will be interested in certain operators that act on functions defined on electrical 20 resistance networks.
21 Definition 2.8. The Laplacian on G is the linear difference operator which acts on a func- 0 22 tion v : G R by → (∆v)(x) := c (v(x) v(y)). (2.4) xy − Xy x 23 ∼ 0 24 A function v : G R is called harmonic iff ∆v 0. → ≡ 25 Definition 2.9. The transfer operator on G is the linear operatorT which acts on a function 0 26 v : G R by →
(T v)(x) := cxyv(y). (2.5) Xy x ∼ 27 Hence, the Laplacian may be written ∆= c T, where (cv)(x) := c(x)v(x). − 28 We won’t worry about the domain of ∆ or T until 8. For now, consider both of these 0 § 29 operators as defined on any function v : G R. The reader familiar with the literature → 30 will note that the definitions of the Laplacian and transfer operator given here are normal- 31 ized differently than may be found elsewhere in the literature. For example, [DS84] and 32 other probabilistic references use
1 1 ∆ := c− ∆= 1 P, so (∆ v)(x) := c (v(x) v(y)), (2.6) c − c c(x) xy − Xy x ∼ 1 33 where P := c− T is the probabilistic transition operator corresponding to the transition 34 probabilities p(x, y) = cxy/c(x). For another example, [Chu01] and other spectral-theoretic 35 references use 24 PALLEE.T.JORGENSENANDERINP.J.PEARSE
1/2 1/2 1/2 1/2 cxyv(y) ∆s := c− ∆c− = 1 c− T c− , so (∆sv)(x) := v(x) . (2.7) − − Xy x c(y) ∼ p 1 However, these renormalized version are much more awkward to work with in the 2 present context; especially when dealing with the inner product and kernels of the Hilbert 3 spaces we shall study. Not only are (2.4) and (2.5) are better suited to the electrical resis- 4 tance network framework (as will be evinced by the operator theory developed in 4 and 5 § succeeding sections) but both ∆c and ∆s are bounded operators, and hence do not allow for 6 the delicate spectral analysis carried out in 7– sec:Lap-on-HE. 7 § §
2.2. The energy. In this section we study the relation between the energy and Laplacian 8 E ∆ on finite networks, as expressed in Lemma 2.13. This formula will be used prolifically, 9 as it also holds on infinite networks in many circumstances. In fact, a noticeable portion of 10 4 is devoted to determining when this is so. 11 § Definition 2.10. The graph energy of an electrical resistance network is the quadratic form 12 0 defined for functions u : G R by 13 → 1 (u) := c (u(x) u(y))2. (2.8) E 2 xy − x,Xy G0 ∈ There is also the associated bilinear energy form 14
1 (u, v) := c (u(x) u(y))(v(x) v(y)). (2.9) E 2 xy − − x,Xy G0 ∈ For both (2.8) and(2.9), note that cxy = 0 for vertices which are not neighbours, and hence 15 1 only pairs for which x y contribute to the sum; the normalizing factor of corresponds 16 ∼ 2 to the idea that each edge should only be counted once. The domain of the energy is 17
0 . dom = u : G R .. (u) < . (2.10) E { → E ∞}
The close relationship between the energy and the conductances is highlighted by the 18 simple identities 19
(δ ) = c(x), and (δ ,δ ) = c , (2.11) E x E x y − xy 0 where δ is a (unit)Dirac mass at x G . Theeasy proofis left as an exercise. A significant 20 x ∈ upshot of (2.11) is that the Dirac masses are not orthogonal with respect to energy. 21
Remark 2.11. It is immediatefrom (2.8) that (u) = 0ifandonlyif u is a constant function. 22 E The energy form is positive semidefinite, but if we work modulo constant functions, it 23 becomes positive definite and hence an inner product. We formalize this in Definition 4.1 24 and again in 6.1. In classical potential theory (or Sobolev theory), this would amount to 25 § working with the class of functions satisfying f < , but abandoning the requirement 26 k ′k2 ∞ that f < . As a result of this, the nontrivial harmonic functions play an especially 27 k k2 ∞ important role in this paper. In particular, it is precisely the presence of nontrivial harmonic 28 functions which prevents the functions of finite support from being dense in the space of 29 functions of finite energy; see 4.2. 30 § OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 25
1 Traditionally (e.g., [Kat95,FOT94¯ ]) the study of quadratic forms would combine (u, v) E 2 and u, v 2 . In our context, this is counterproductive, and would eclipse some of our h iℓ 3 most interesting results. Some of our most intriguing questions for elements v 2 0 ∈ HE 4 involve boundary considerations, and in these cases v is not in ℓ (G ) (Corollary 4.65). 5 One example of this arises in the discrete Gauss-Green formula (Theorem 4.41); another 6 arises in study of forward-harmonic functions in 11.2. § 7 The following proposition may be found in [Str06, 1.3]or[Kig01, Ch. 2], for example. § 8 Proposition 2.12. The following properties are readily verified:
9 (1) (u, u) = (u). E E 1 10 (2) (Polarization) (u, v) = [ (u + v) (u v)]. E 4 E −E − 11 (3) (Markov property) ([u]) (u), where [u] is any contraction of u. E ≤E 12 For example, let [u] := min 1, max 0, u . The following result relates the Laplacian to { { }} 13 the graph energy on finite networks, and can be interpreted as a relation between dom 2 0 E 14 and ℓ (G ).
15 Lemma 2.13. Let G be a finite electrical resistance network. For u, v dom , ∈ E (u, v) = u(x)∆v(x) = v(x)∆u(x). (2.12) E xXG0 xXG0 ∈ ∈ 16 Proof. Direct computation yields
1 (u, v) = c u(x)v(x) u(x)v(y) u(y)v(x) + u(y)v(y) E 2 xy − − x,Xy G0 ∈ 1 1 = c(x)u(x)v(x) + c(y)u(y)v(y) 2 2 xXG0 yXG0 ∈ ∈ 1 1 u(x)T v(x) u(y)T v(y) − 2 − 2 xXG0 yXG0 ∈ n ∈ = c(x)u(x)v(x) u(x)T v(x) − xXG0 xX,y G0 ∈ ∈ = u(x) (c(x)v(x) T v(x)) − xXG0 ∈ = u(x)∆v(x). (2.13) xXG0 ∈ 17 Of course, the computation is identical for x G0 v(x)∆u(x). ∈ P 18 We include the following well-known result for completeness.
19 Corollary 2.14. On a finite electrical resistance network, all harmonic functions of finite 20 energy are constant.
21 Proof. If h is harmonic, then (h) = x G0 h(x)∆h(x) = 0. See Remark 2.11. E ∈ P 22 Connectedness is implicit in the calculations of both Lemma 2.13 and Corollary 2.14; 23 recall that all electrical resistance networks considered in this work are connected. We 24 will extend Lemma 2.13 to infinite graphs in Theorem 4.41, where the formula is more 25 complicated: 26 PALLEE.T.JORGENSENANDERINP.J.PEARSE
(u, v) = u(x)∆v(x) + “boundary term” . E { } xXG0 ∈ It is shown in Theorem 4.47 that the presence of the boundary term corresponds to the 1 existence of nontrivial harmonic functions, in contrast to Corollary 2.14. In fact, one can 2 interpret Corollary 2.14 as the reason why the boundary term alluded to above vanishes on 3 finite networks. We study the interplay between and ∆ further in 9.2– 9.3. 4 E § § OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 27
1 3. Currents and potentials on electrical resistance networks
2 The potential theory for an electrical resistance network is studied via an experiment 3 in which 1 amp of current is passed through the network, inserted into one vertex and 4 extracted at some other vertex. The voltage drop measured between the two nodes is the 5 effective resistance between them, see 5. § 6 When the voltages are fixed at certain vertices, it induces a current in the network in 7 accordance with the laws of Kirchhoff and Ohm. This induced current is introduced for- 8 mally in Definition 3.17. Induced currents are important for studying flows of minimal 9 dissipation, and will also be useful in the study of forward-harmonic functions in 11.2. If § 10 a voltage drop of 1 volt is imposed between two vertices, the effective resistance between 11 these two vertices is the reciprocal of the dissipation of the induced current. 12 In Theorem 3.28 we show that there always exists an harmonic function satisfying the 13 boundary conditions implied by the above described experiment, in order to fill a gap 14 in [Pow76b]. In Theorem 3.26 and Theorem 3.26 it is shown that these harmonic functions 15 correspond to currents which minimize energy dissipation.
16 3.1. Currents on electrical resistance networks. 0 0 17 Definition 3.1. A current is a skew-symmetric function I : G G R. × → 18 Definition 3.2. An orientation is a subset of the edges which includes exactly one of each 19 pair (x, y), (y, x) . For a given current I, one may pick an orientation by requiring that { } 20 I(x, y) > 0 on every edge for which I is nonzero, and arbitrarily choosing (x, y) or(y, x) 21 outside the support of I. We refer to this as an orientation induced by the current; this will 22 be used extensively in 11.2 to study the forward-harmonic functions. § 0 23 The energyis a functionaldefined on functions v : G R which give voltages between → 24 different vertices in the network. The associated notion defined on the edges of the network 25 is the dissipation of a current.
26 Definition 3.3. The dissipation of a current may be thought of as the energy lost as a 1 27 current flows through an electrical resistance network. More precisely, for I, I , I : G 1 2 → 28 R,
1 1 2 D(I) := c− I(x, y) . (3.1) 2 xy (x,Xy) G1 ∈ 29 The associated bilinear form is the dissipation form:
1 1 D(I , I ) := c− I (x, y)I (x, y). (3.2) 1 2 2 xy 1 2 (x,Xy) G1 ∈ 1 30 Again, the normalizing factor of 2 corresponds to the idea that each edge only contributes 31 once to the sum. The domain of the dissipation is
. dom D := I .. D(I) < . (3.3) { ∞} 32 Remark 3.4. When an orientation for G is chosen,it is easy to see that dom D is a Hilbert O 2 33 space under the inner product (3.2). Indeed, dom D = ℓ ( , c). O 34 Definition 3.5. A cycle ϑ is a set of n edges corresponding to a sequence of vertices 0 1 35 (x , x , x ,..., x = x ) G , for which (x , x + ) G for each k. Denote the set of 0 1 2 n 0 ⊆ k k 1 ⊆ 36 cycles in G by . L 28 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Definition 3.6. For physical realism, we often require that a current flow satisfy Kirch- 1 hoff’s node law, i.e., that the total current flowing into a vertex must equal the total current 2 flowing out of a vertex: 3
I(x, y) = 0, x G0. (3.4) ∀ ∈ Xy x ∼ This is indeed the version of Kirchhoff’s law you would find in a physics textbook; with 4 our convention I(x, y) > 0 indicates that the current flows from x to y. 5 However, if we are performing the experiment described above, then there are boundary 6 conditions at α, ω to take into account, and Kirchhoff’s node law takes the nonhomoge- 7 neous form 8
1, x = α,
I(x, y) = δα δω = 1, x = ω, (3.5) − Xy x − ∼ 0, else, 0 where δx is the usual Dirac mass at x G . 9 ∈ Definition 3.7. A current flow from α to ω is a current I dom D that satisfies (3.5). The 10 ∈ set of all current flows is denoted (α, ω). 11 F We usually use α to denote the beginning of a flow and ω to denote its end. Shortly, we 12 will see that the currents corresponding to potentials are precisely the current flows. 13
Remark 3.8. Although trivial, it is important to note that the characteristic function of a 14 1 current path χ : G 0, 1 trivially satisfies (3.5). Also, the characteristic function 15 γ → { } of a cycle satisfies (3.4) in much the same way. As a consequence, if I (α, ω), then 16 ∈ F I + tχ (α, ω) for any t R by a brief computation. In other words, perturbation 17 ϑ ∈ F ∈ on a cycle preserves the Kirchhoff condition. However, the dissipation will vary because 18 χ D( ϑ) > 0. 19
3.2. Potential functions and their relationship to current flows. From the proceeding 20 0 section, it is clear that a special role is played by functions v : G R which satisfy 21 → the equation ∆v = δα δω. Such a function is the solution to a discrete Dirichlet problem, 22 − where the “boundary”has been chosen to be α and ω (not to be confused with the boundary 23 term discussed in Remark 4.7). 24
Definition 3.9. A dipole is a function v dom which satisfies 25 ∈ E
∆v = δα δω (3.6) 0 − for some vertices α, ω G . The collection of all such functions is denoted (α, ω). 26 ∈ P Note that when G is finite, (α, ω) contains only a single element. This follows from 27 P Corollary 2.14 because the difference of any two solutions to (3.6) is harmonic. 28
Remark 3.10. The definition of a monopole that we give here is a heuristic definition; we 29 0 give the precise definition in Definition 4.34. A monopole at ω is a function w : G R 30 → which satisfies 31
∆w = kδω, w dom , k C. (3.7) ∈ E ∈ In the sequel, we are primarily concerned with monopoles wo, where o = ω is some fixed 32 vertex which acts as a point of reference or “origin”. Also, we typically take k = 1, as the 33 − induced current of such a monopole is a unit flow to infinity in the language of [DS84]. 34 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 29
1 Remark 3.11. The study of dipoles, monopoles, and harmonic functions is a recurring 2 theme of this paper:
∆v = δα δω, ∆w = δω, ∆h = 0. − − 3 In Theorem 3.28, we will show that (α, ω), is nonempty for any α and ω, on any network P 4 (G, c); the existence of monopoles and nontrivial harmonic functionsis a much more subtle 5 issue. 6 In Corollary 4.21,weoffer a more refinedproofof the existence of dipoles, using Hilbert 7 space techniques. Perhaps a more interesting question is when (α, ω) contains more than P 8 element; the linearity of ∆ shows immediately that any two dipoles in (α, ω) differ by a P 9 harmonicfunction. We have shown that when a connected graph is finite the only harmonic 10 functions are constant (Corollary 2.14), and therefore (α, ω) consists only of a single P 11 function, up to the addition of a constant. The situation for monopoles is similar, as the 12 difference of two monopoles at ω is also a harmonic function. 13 Not all electrical resistance networks support monopoles; the current induced by a 14 monopole is a finite flow to infinity and hence indicates that the random walk on the net- 15 work is transient, by [Lyo83]. See also [DS84, LP09] for terminology and proofs. It is 16 well-known that for a reversible Markov chain, if the random walk started at one vertex is 17 transient, then it is transient when started at any vertex. We give a very brief proof of this 18 in Lemma 3.29; and a new criterion for transience in Lemma 4.51. 19 On some networks, a monopole can be understood as the limit of a sequence of dipoles 20 v where ∆v = δ δ and x . In such a situation, a monopole can be considered xn xn xn − o n →∞ 21 as a dipole where one of the Dirac masses “sits at ”. However, this is not possible on ∞ 22 all networks, as is illustrated by the binary tree in Example 13.4. Again, the linearity of ∆ 23 shows immediately that any two monopolesat ω differ by a harmonicfunction. When these 24 monopoles correspond to a “distribution of dipoles at infinity” (i.e., a limit of sums axvx 25 where the v ’s are dipoles with x in the limit), the addition of a harmonic function x → ∞ P 26 transforms the distribution at infinity. It will take some work to make these ideas precise; 27 for now the reader can consider this remark simply as a preview of coming attractions. The 28 presence of monopoles is also extremely closely related to the existence of “long-range 3 29 order”, and the theoretical foundation of magnetism in R ; see 15.3. § 30 Furthermore, it is possible for an electrical resistance network to support monopoles but d 31 not nontrivial harmonic functions. In 14, we show that the integer lattice networks (Z , 1) § 32 support monopoles (Theorem 14.5). However, all harmonic functions are linear and hence 33 do not have finite energy; cf. Theorem 14.17. Both of these results are well-known in the 34 literature in different contexts, and/or with different terminology.
35 Lemma 3.12. The dipoles (α, ω) and the current flows (α, ω) are convex sets. Fur- P F 36 thermore, if v (α, ω), then v + h (α, ω) for any harmonic function h; similarly, if ∈ P ∈ P 37 I (α, ω), then I + J (α, ω) for any function J satisfying (3.4). ∈F ∈F 38 Proof. If v (α, ω), c 0 and c = 1, then the linearity of ∆ gives i ∈P i ≥ i P
∆ civi = ci∆vi = ci(δα δω) = δα δω. X X X − − 39 The computation for the other parts is similar.
40 Theorem 3.13. obtains its minimum for some unique v (α, ω), and D obtains its E ∈ P 41 minimum for some unique I (α, ω). ∈F 30 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Proof. Each of these is a quadratic form on a convex set, by Lemma 3.12, so the result 1 is an immediate application of [Rud87, Thm. 4.10] or [Nel69], e.g. To underscore the 2 . uniqueness, suppose that (v ) = (v ). Then with ε := inf (v) .. v (α, ω) , the 3 E 1 E 2 {E ∈ P } parallelogram law gives 4
(v v ) = 2 (v ) + 2 (v ) 4ε2 = 0, E 1 − 2 E 1 E 2 − since (v ) = ε because v were chosen to be minimal. 5 E i i Definition 3.14. Ohm’s Law (V = RI) appears in the present context as 6 1 v(x) v(y) = I(x, y). (3.8) − cxy
Remark 3.15. It will shortly become evident (if it isn’t already) that current flows satis- 7 fying Kirchhoff’s law correspond to harmonic functions via Ohm’s law and that current 8 flows satisfying the nonhomogeneous Kirchhoff’s law (3.5) correspond to dipoles, that is, 9 solutions of the Dirichlet problem (3.6) with Neumann boundary conditions. To make this 10 precise, we need the notion of induced current given in Definition 3.17 and justified by 11 Lemma 3.16. 12 0 Lemma 3.16. Every function v : G R induces a unique current via I(x, y) := c (v(x) 13 → xy − v(y)), and the dissipation of this current is the energy of v: 14
D(I) = (v). (3.9) E Moreover, if v (α, ω), then I (α, ω). 15 ∈P ∈F Proof. It is clear that Ohm’s Law defines a current. The equality (3.9) is a very brief 16 calculation and follows straight from the definitions; see (2.9) and (3.1). A proofof (3.9) 17 is also given in [DS84]. 18 If v (α, ω), then ∆v = δα δω and 19 ∈P −
(δα δω)(x) = (∆v)(x) = c (v(x) v(y)) = I(x, y) (3.10) − xy − Xy x Xy x ∼ ∼ verifies the nonhomogeous form of Kirchhoff’s law. 20
Definition 3.17. Given v (α, ω), the induced current is defined via Ohm’s Law as in 21 ∈ P the statement of Lemma 3.16. That is, 22
I(x, y) := c (v(x) v(y)). (3.11) xy − Remark 3.18. Note that (3.9) holds when the current I is induced by v. It makes no sense 23 to attempt to apply the same equality to a general current: there may be NO associated 24 potential because of the compatibility problem described just below. Nonetheless, Theo- 25 rem 10.27 provides a way to give the identity analogous to (3.9) for general currents by 26 using the adjoint of the operator implicit in (3.11). 27
Remark 3.19. If ∆v = δα δω has a solution v , then any other solution is of the form 28 − 0 v = v0 + h where h is harmonic, by linearity of ∆. So to minimize energy, one must 29 consider such perturbations: 30
d [ (v + th)] = = 0 (v , h) = 0. dt E 0 t 0 ⇐⇒ E 0 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 31
x y
Figure 1. A Dirac mass on an edge of Z2.
1 Conversely, if (v, h) = 0, then E
(v + th) = (v) + 2t (v, h) + t2 (h) (v), E E E E ≥E 2 shows that energy is minimized for t = 0. In particular, energy is minimized for v which 3 contains no harmonic component. In Lemma 4.22 this important principle is restated in 4 the language of Hilbert spaces: energy is minimized for the v which is orthogonal to the 5 space of harmonic functions with respect to . E 6 Analogous remarks hold for I which minimizes D(I). However, note that Kirchhoff’s 7 Law is blind to conductances and so I (α, ω) does not imply that D(I) is minimal. In ∈ F 8 the next section, we show that induced currents are minimal with respect to D when they 9 are induced by a minimal potential v.
10 3.3. The compatibility problem. The converse to Lemma 3.16 is not always true, but a 11 partial converse is given by Theorem 3.26. Given an electrical resistance network (G, c), 12 one can always attempt to construct a Ohm’s function by fixing the value v(x0) at some 0 13 point x G , and then applying Ohm’s law to determine the value of v for other vertices 0 ∈ 14 x x . However, this attempt can fail if the network contains a cycle (see Example 12.2 for ∼ 0 15 an example) because the existence of a cycle is equivalent to the existence of two distinct 16 paths from one point to another. This phenomenonis worked out in detail for a simple case 17 in Example 3.20. 18 In general, it may happen that there are two different paths from x0 to y0, and the net 19 voltage drop v(x ) v(y ) computedalong these two paths is not equal. Such a phenomenon 0 − 0 20 makes it impossible to define v. Note that Kirchhoff’s law does not forbid this, because 21 (3.4)–(3.5) is expressed without reference to the conductances c. We refer to this as the 22 compatibility problem: a general current function may not correspond to a potential, 23 even though every potential induces a well defined current flow (see Lemma 3.16). In this 24 section we provide the following answer: for any current, there exists a unique associated 25 current which does correspond to a potential.
26 Example 3.20 (The Dirac mass on an edge). Consider a Dirac mass on an edge of the 2 27 network (Z , 1) as depicted in Figure 1. We use such a current here to illustrate the com- 28 patibility problem. To find a potential corresponding to this current, consider the following 29 dilemma: I(x, y) = 1 and I 0 elsewhere corresponds to a potential (up to a constant) ≡ 30 which has v(x) = 1 and v(y) = 0, as in Figure 2. Since I(x, w) = 0, we have v(w) = v(x), 31 and since I(y, z) = 0, we have v(z) = v(y). But then v(z) = 1 , 0 = v(z), contradicting the < 32 fact that I(w, z) = 0! ւ χ 33 Definition 3.21. A current I satisfies the cycle condition iff D(I, ϑ) = 0 for every cycle 34 ϑ . (We call χ a cycle.) ∈ L ϑ 32 PALLEE.T.JORGENSENANDERINP.J.PEARSE
x w y z
Figure 2. A failed attempt at constructing a potential to match Figure 1.
Remark 3.22. From the preceding discussion, it is clear that for a current satisfying the 1 cycle condition, voltage drop between vertices x and y may be measured by summing the 2 currents along any single path from x to y, and the result will be independentof which path 3 was chosen. In the Hilbert space interpretation of 10 the cycle condition is restated as “I 4 § is orthogonal to cycles”. The next two results must be folklore (perhaps dating back to the 5 th 19 century?) but we include them for their relevance in 10, especially the Hilbert space 6 § decomposition of Theorem 10.8 (see also Figure 4). While writing a second draft of this 7 document, the authors discovered a similar treatment in [LP09, 9]. 8 § Lemma 3.23. I is an induced current if and only if I satisfies the cycle condition. 9
Proof. ( ) If I is induced by v, then for any ϑ , the sum 10 ⇒ ∈ L 1 I(x, y) = (v(x) v(y)) = 0, (3.12) c − (xX,y) ϑ xy (xX,y) ϑ ∈ ∈ χ since every term v(xi) appears twice, once positive and once negative, whence D(I, ϑ) = 0. 11 ( ) Conversely, to prove that there is such a v, we must show that v(x ) v(y ) is 12 ⇐ 0 − 0 independence of the path from x0 to y0 used to compute it. In a direct analogy to basic 13 vector calculus, this is equivalent to the fact that the net voltage drop around any closed 14 cycle is 0. 15
1 (v(x) v(y)) = I(x, y) = D(I, χ ) = 0, − c ϑ (xX,y) ϑ (xX,y) ϑ xy ∈ ∈ 0 Now define v by fixing v(x0) for some point x0 G , and then coherently use v(x) v(y) = 16 1 ∈ − I(x, y) to compute v at any other point. 17 cxy
The presence of cycles is not always obvious! As an exercise, we invite the reader to 18 determine the cycles involved in Example 3.20. 19
Lemma 3.24 (Resurrection of Kirchhoff’s Law). Let I be the current induced by v. Then 20 v (α, ω) if and only if I satisfies the nonhomogeneous Kirchhoff’s law. 21 ∈P Proof. ( ) Computing directly, 22 ⇒
I(x, y) = c (v(x) v(y)) = ∆v(x) = δα δω. (3.13) xy − − Xy x Xy x ∼ ∼ ( ) Conversely, to show ∆v = δα δω, 23 ⇐ −
∆v(x) = c (v(x) v(y)) defof ∆ xy − Xy x ∼ OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 33
= I(x, y)(3.8) Xy x ∼ = δα δω, I (α, ω). − ∈F
1 Corollary 3.25. Let I be the current induced by v. Then v is harmonic if and only if I 2 satisfies the homogeneous Kirchhoff’s law.
3 Proof. Mutatis mutandis, this is the same as the proof of Lemma 3.24.
4 Theorem 3.26. I minimizes D on (α, ω) if and only if I is induced by a potential v that F 5 lies in (α, ω). Moreover, v also minimizes over (α, ω). P E P 6 Proof. ( ) Since I minimizes D, we have ⇒ d [D(I + tJ)] = 0, (3.14) dt t=0 7 for any current J satisfying the homogeneous Kirchhoff’s law. From Remark 3.8, this χ 8 applies in particular to J = ϑ, where ϑ is any cycle in . d L 9 Note that D(I, χ ) = 0 iff D(I + tχ ) = 0. To see this, replace I by I + tχ in (3.1), ϑ dt ϑ t=0 ϑ χ 10 differentiate D(I + t ϑ) term-by-termh withi respect to t and evaluates at t = 0 to obtain that 11 (3.14) is equivalent to
1 1 χ I(x, y) ϑ(x, y) = I(x, y) = 0, ϑ . cxy cxy ∀ ∈ L (x,Xy) G1 (xX,y) ϑ ∈ ∈ 12 By Lemma 3.23, this shows that I is induced by some v; and by Lemma 3.24, we know 13 v (α, ω). From Lemma 3.16, it is clear that the v must also be the energy-minimizing ∈ P 14 element of (α, ω). P 15 ( ) Since I is induced by v (α, ω), the only thing we need to check is that I is min- ⇐ ∈P 16 imal with respect to any harmonic current (i.e. a current induced by a harmonic function); 17 this follows from Lemma 3.23 and the first part of the proof. If h is any harmonic function 0 18 on G , denote the induced current by H as before. Then Lemma 3.16 gives
d d [D(I + tH)] = [ (v + th)] = 0, dt t=0 dt E t=0 19 by the minimality of v.
20 Remark 3.27. Part of the motivation for Theorem 3.26 is to fix an error in [Pow76b]. The 21 author was not apparently aware of the possibility of nontrivial harmonic functions, and 22 hence did not see the need for taking the element of (α, ω) with minimal energy. This P 23 becomes especially important in Theorem 5.2. 24 Theorem 3.26 is generalized in Theorem 10.27 where we exploit certain operators to 25 obtain, for any given current I, an associated minimal current. This minimal current is 26 induced by a potential, even if the original is not, and provides a resolution to the compat- 27 ibility problem described at the beginning of 3.3, just above. § 28 In 10.4.1 we revisit this scenario and show how the minimal current may be obtained § 29 by the simple application of a certain operator, once it has been properly interpreted in 30 terms of Hilbert space theory. See Theorem 10.27 and its corollaries in particular.
31 Theorem 3.28. (α, ω) is never empty. P 34 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Proof. It is clear that (α, ω) , ∅ because one always has the characteristic function of 1 F a current path from α to ω (since we are assuming the underlying graph is connected); 2 see Definition 11.12 and Remark 3.8. From Theorem 3.13 one sees that there is always a 3 flow which minimizes dissipation. By Theorem 3.26, this minimal flow is induced by an 4 element of (α, ω). 5 P The following result is well-known in probability (see, e.g., [Str05]), but we include it 6 here for completeness and the novel method of proof. 7 0 Corollary 3.29. If the random walk on (G, c) is transient when started from y G , then 8 0 ∈ it is transient when started from any x G . 9 ∈ Proof. By [Lyo83], the hypothesis means there is a monopole w dom with ∆w = δ . 10 y ∈ E y But then by Theorem 3.28 and the linearity of ∆, v + wx is a monopole at x, for any 11 v (x, y). 12 ∈P Remark 3.30. Theorem 3.28 fills a gap in [Pow76b]. A key point is that the finite dissipa- 13 tion of the flow ensures the finite energy of the inducing voltage function, by Lemma 3.16. 14 A different proof of Theorem 3.28 is obtained in Corollary 4.21 by the application of 15 Hilbert space techniques. 16 Theorem 3.28 also follows from results of [Soa94, III.4] since the difference of two 17 § Dirac masses corresponds to a “balanced” flow, i.e., the same amount of current flows in 18 as flows out. 19 2 0 Remark 3.31. There are examples for which the elements of (α, ω) do not lie in ℓ (G ); 20 P see Figure 10 of Example 14.2. 21 0 Proposition 3.32. IfG is finite and v (α, ω), then v(ω) v(x) v(α) for all x G . 22 ∈P ≤ ≤ ∈ Proof. This is immediate from the maximum principle for harmonic functions on the finite 23 0 set G with boundary α, ω . See[LP09, 2.1], for example, or [LPW08]. 24 { } § Remark 3.33. In 5, we will see that Proposition 3.32 extends to a more general result: if v 25 § is the unique element of (α, ω) of minimal energy, then the same conclusion follows. One 26 P way to see this is to define u(x) = Px[τα < τω] (i.e., the probability that the random walk 27 W started at x reaches α before ω). By Theorem 5.18, v is defined by v(x) = u(x)R (α, ω), 28 W where R (α, ω) is the (wired) effective resistance between α and ω. 29 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 35
1 4. The energy Hilbert space
2 In this section, we study the Hilbert space of voltage functions, where the inner HE 3 product is given by the energy form. A key feature of is the presence of a family of HE 4 dipoles vx x G0 indexed by the vertices. In Corollary 4.13 we show how this family of { } ∈ 5 dipoles forms a reproducing kernel for in the sense of Aronszajn. HE 6 In Theorem 4.41, we establish a discrete version of the Gauss-Green Formula which 7 extends Lemma 2.13 to the case of infinite graphs. The appearance of a somewhat mys- 8 terious boundary term prompts several questions which are discussed in Remark 4.7. An- 9 swering these questions comprises a large part of the sequel; cf. 7. We are able to prove § 10 in Lemma 4.68 that this boundary term vanishes for finitely supported functions on G, and 2 0 11 in Corollary 4.65 that nontrivial harmonic functions cannot be in ℓ (G ). This is discussed 12 further in Remark 4.66 and provides the motivation for energy-centric approach we pursue 13 throughout our study. 14 The energy Hilbert space will facilitate our study of the resistance metric R in 5. HE § 15 In particular, it provides an explanation for an issue stemming from the “nonuniqueness of 16 currents” in certain infinite networks; see [LP09,Tho90]. This disparity leads to differences 17 between two apparently natural extensions of the effective resistance to infinite networks, 18 which are greatly clarified by the geometry of Hilbert space. Also, presents an analytic HE 19 formulation of the type problem for random walks on an electrical resistance network: 20 transience of the random walk is equivalent to the existence of monopoles, that is, finite- 21 energy solutions to a certain Dirichlet problem. In fact, this approach will readily allow us 22 to obtain explicit formulas for effective resistance on integer lattice networks in 14, with § 23 applications to a physics problem of [Pow76b] in 15. § 24 Definition 4.1. The energy form is symmetric and positive definite, and its kernel is E 25 the set of constant functions on G. Let 1 denote the constant function with value 1. Then 26 dom /R1 is a vector space with inner product and corresponding norm given by E u, v := (u, v) and u := (u, u)1/2. (4.1) h iE E k kE E 27 Upon completion with respect to this inner norm, we obtain the energy Hilbert space . 0 HE 28 Fix a reference vertex o G to act as an “origin”. It will readily be seen that all results ∈ 29 are independent of this choice.
30 Remark 4.2. In 6, we provide an alternative construction of via techniques of von § HE 31 Neumann and Schoenberg. This provides for a more explicit description of the structure 32 of and its relation to the metric geometry of (G, R), and shows that is the natural HE HE 33 Hilbert space in which to embed (G, R). However, this must be postponed until after the 34 introduction of the effective resistance metric.
35 Remark 4.3 (Three warnings about ). HE 36 (1) has no canonical o.n.b.; the usual candidates δx are not orthogonal and typi- HE { } 37 cally not even dense, as we discuss further below. 38 (2) Multiplication operators are not Hermitian; see Lemma 4.4 and Remark 7.20. 39 (3) Pointwise identities should not be confused with Hilbert space identities; see Re- 40 mark 4.19 and Lemma 4.33.
41 To elaborate on the last point, note that elements of are technically equivalence classes HE 42 of functionswhich differonlybya constant; thisis whatis meantbythe notationdom /R1. E 43 In other words, if v1 = v2 + k for k C, then v1 = v2 in . When working with represen- ∈ HE 44 tatives, we typically abuse notation and use u to denote the equivalence class of u. Often, 36 PALLEE.T.JORGENSENANDERINP.J.PEARSE we choose u so that u(o) = 0 (occasionally without warning). A different but no less useful 1 choice is to pick k so that v = 0 outside a finite set when v is a function of finite support 2 (see Definition 4.16). 3
Lemma 4.4. If ϕ and Mϕ denotes the multiplication operator Mϕ : u ϕu, then 4 ∈ HE 7→ Mϕ is Hermitian if and only if ϕ = 0 in . 5 HE Proof. From the formula (2.9), 6
1 Mϕu, v = cxy(ϕ(x)u(x)v(x) ϕ(x)u(x)v(y) ϕ(y)u(y)v(x) + ϕ(y)u(y)v(y)). h iE 2 − − x,Xy G0 ∈ By comparison with the corresponding expression, this is equal to u, Mϕv iff (ϕ(y) 7 h iE − ϕ(x))u(y)v(x) = (ϕ(y) ϕ(x))u(x)v(y). However, since we are free to vary u and v, it must 8 − be the case that ϕ is constant. 9
Definition 4.5. An exhaustion of G is an increasing sequence of finite and connected sub- 10 graphs G , so that G G + and G = G . 11 { k} k ⊆ k 1 k S Definition 4.6. The notation 12
:= lim (4.2) k xXG0 →∞ xXGk ∈ ∈ is used whenever the limit is independent of the choice of exhaustion G of G. We 13 { k} typically justify this independence by proving the sum to be absolutely convergent. 14
Remark 4.7. One of the main results in this section is a discrete version of the Gauss-Green 15 theorem presented in Theorem 4.41: 16
u, v = u(x)∆v(x) + u(x) ∂v (x), u, v . (4.3) h iE ∂Ò ∈HE xXG0 xXbd G ∈ ∈ This differs from the literature, where it is common to find (u, v) = u, ∆v 2 given as a 17 E h iℓ definition (of or of ∆, depending on the context), e.g. [Kig01, Str08]. After reading a 18 E preliminary version of this paper, a reader pointed out to us that a similar formula appears 19 in [DK88, Prop 1.3]; however, these authors apparently do not consider the extension of 20 this formula to infinite networks. 21 ∂v We refer to u as the “boundary term” by analogy with classical PDE theory. 22 bd G ∂Ò This terminologyP should not be confused with the notion of boundary that arises in the dis- 23 cussion of the discrete Dirichlet problem. In particular, the boundary discussed in [Kig03] 24 0 and [Kig08] refers to a subset of G . By contrast, when discussing an infinite network G, 25 our boundary bdG is never contained in G. Green’s identity follows immediately from 26 (4.3) in the form 27
(u(x)∆v(x) v(x)∆u(x)) = v(x) ∂u (x) u(x) ∂v (x) . (4.4) Ò − ∂Ò − ∂ xXG0 xXbd G ∈ ∈ Note that our definition of the Laplace operator is the negative of that often found in the 28 PDE literature, where one will find Green’s identity written 29
(u∆v v∆u) = (u ∂ v v ∂ u). Ò ∂Ò ∂ ZΩ − Z∂Ω − OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 37
1 As the boundary term may be difficult to contend with, it is extremely useful to know 2 when it vanishes. We have several results concerning this:
3 (i) Lemma 4.68 shows the boundary term vanishes when either argument of u, v has h iE 4 finite support, 5 (ii) Lemma 4.47 gives necessary and sufficient conditions on the electrical resistance 6 network for the boundary term to vanish for any u, v , ∈HE 7 (iii) Lemma 4.63 show the boundary term vanishes when both arguments of u, v and 2 h iE 8 their Laplacians lie in ℓ .
9 In fact, Lemma 4.47 expresses the fact that it is precisely the presence of monopoles that 10 prevents the boundaryterm from vanishing. An example with nonvanishing boundary term 11 is given in Example 13.5.
12 4.1. The evaluation operator Lx and the reproducing kernel vx. 0 13 Definition 4.8. For any vertex x G , define the linear evaluation operator Lx on by ∈ HE L u := u(x) u(o). (4.5) x − 14
0 1/2 15 Lemma 4.9. For any x G , one has L u k (u) , where k depends only on x. ∈ | x |≤ E 16 Proof. Since G is connected, we can choose an n-step path from x to o and apply the 17 Schwarz inequality to obtain
n 2 2 2 cxi,xi 1 Lxu = u(x) u(o) = − (u(xi) u(xi 1)) | | | − | c − − Xi=1 r xi,xi 1 − n n 1 2 cxi,xi 1 (u(xi) u(xi 1)) − ≤ cxi,xi 1 − − Xi=1 − Xi=1 k2 (u), ≤ E 1/2 n 1 18 for k = i=1 . cxi,xi 1 P − 19 Definition 4.10. Let vx be defined to be the unique element of for which HE
vx, u = u(x) u(o), for every u . (4.6) h iE − ∈HE 20 This is justified by Lemma 4.9 and the Riesz Representation Theorem.
21 Definition 4.11. Let be a Hilbert space of functions on X. An operator S on is said H H 22 to have a reproducing kernel kx x X iff { } ∈ ⊆H
(Sv)(x) = kx, v , x X, v . (4.7) h iH ∀ ∈ ∀ ∈H 23 If S is projection to a subspace L , then one says k is a reproducing kernel for L. If ⊆H { x} 24 S = I, then is a reproducing kernel Hilbert space with kernel k. H 25 Theorem 4.12 (Aronszajn’s Theorem [Aro50]). Let f be a reproducing kernel for . { x} H 26 Define a sesquilinear form on the set of all finite linear combinations of these elements by
ξx fx, ηy fy := ξxηx fx(y). (4.8) *Xx Xy + Xx 38 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Then the completion of this set under the form (4.8) is again . 1 H Corollary 4.13. vx x G0 is a reproducing kernel for . Thus, span vx is dense in . 2 { } ∈ HE { } HE Proof. Choosing representatives with vx(o) = 0, it is trivial to check that vx, vy = vx(y) = 3 h iE vy(x) and then apply Aronszajn’s Theorem. 4
We note that there is a rich modern literature dealing with reproducing kernels and 5 their manifold application to both continuous analysis problems (see e.g., [AD06, AL08, 6 AAL08, BV03, Zha09]), and infinite discrete stochastic models. One of the differences 7 between these studies and our present work is the approach we take in Definition 4.9, i.e., 8 the use of “relative” reproducing kernels. 9
Definition 4.14. The family of functions vx x G0 is called the energy kernel. Note that vo 10 { } ∈ corresponds to a constant function, since vo, u = 0 for every u . Therefore, will 11 h iE ∈ HE 0 often ignore or omit this term and sometime write summations over the set G o . 12 \ { } Remark 4.15. Definition 4.14 is justified by Corollary 4.13. In this paper, the functions vx 13 will play a role analogous to fundamental solutions in PDE theory; see 10.3. 14 § The functions vx are R-valued. This can be seen by first constructing the energy kernel 15 for the Hilbert space of R-valued functions on G, and then using the decomposition of 16
a C-valued function u = u1 + iu2 into its real and imaginary parts. Alternatively, see 17 Lemma 4.29. 18 Reproducing kernels will help with many calculations and explain several of the re- 19 lationships that appear in the study of electrical resistance networks. They also extend 20 the analogy with complex function theory discussed in 10.3. The reader may find the 21 § references [Aro50, Yoo07, Jor83] to provide helpful background on reproducing kernels. 22
4.2. The finitely supported functions and the harmonic functions. 23
Definition 4.16. Let δx be the Dirac mass at x, i.e., the class of which contains the 24 HE characteristic function of the singleton x . Then the set of “finitely supported functions” 25 { } in is 26 HE
. 0 span δ = u dom .. u(x) = k for some k, for all but finitely many x G , (4.9) { x} { ∈ E ∈ } . Then spt v H means that v span δ .. x H , and v has finite support iff H is finite. 27 ⊆ ∈ { x ∈ } Define in to be the closure of span δ with respect to . 28 F { x} E Definition 4.17. The set of harmonic functions of finite energy is denoted 29
0 . arm = arm(G ) := v .. ∆v 0 . (4.10) H H { ∈HE ≡ }
Lemma 4.18. The Dirac masses δx x G0 form a reproducing kernel for ∆. That is, for any 30 0 { } ∈ x G , one has δx, u = ∆u(x). 31 ∈ h iE Proof. Compute δx, u = (δx, u) directly from formula (2.9). 32 h iE E Remark 4.19. Note that one can take the definition of the Laplacian to be the operator A 33 defined via the equation 34
δx, u = Au(x). h iE OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 39
1 This point of view is helpful, especially when distinguishing between identities in Hilbert 2 space and pointwise equations. For example, if h arm, then ∆h and the constant ∈ H 3 function 1 are identified in because u, ∆h = u, 1 = 0, for any u . However, HE h iE h iE ∈HE 4 one should not consider a (pointwise) solution of ∆u(x) = 1 to be a harmonic function.
0 5 Lemma 4.20. For any x G , ∆v = δ δ . ∈ x x − o 6 Proof. Using Lemma 4.18, ∆vx(y) = δy, vx = δy(x) δy(o) = (δx δo)(y). h iE − − 7 By applying Lemma 4.20 to vα vω, we see: − 8 Corollary 4.21. The space of dipoles (α, ω) is nonempty. P 9 Lemma 4.18 is extremely important. Since in is the closure of span δ , it implies F { x} 10 that the finitely supported functions and the harmonic functions are orthogonal. This result 11 is called the “Royden Decomposition” in [Soa94, VI] and also appears elsewhere, e.g., § 12 [LP09, 9.3]. § 13 Theorem 4.22. = in arm. HE F ⊕H 14 Proof. For all v , Lemma 4.18 gives δx, v = ∆v(x). Since in = span δx , this ∈ HE h iE F { } 15 equality shows v in whenever v is harmonic. Conversely, if δx, v = 0 for every x, ⊥ F h iE 16 then v must be harmonic. Recall that constants functions are 0 in . HE 17 Corollary 4.23. span δx is dense in iff arm = 0. { } HE H 18 Remark 4.24. Corollary 4.23 is immediate from Theorem 4.22, but we wish to emphasize 19 the point, as it is not the usual case elsewhere in the literature. Part of the importance of 20 the energy kernel v arises from the fact that the Dirac masses are generally inadequate { x} 21 as a representing set for . This leads to unusual consequences, e.g., one may have HE
u , u(x)δx, in . HE xXG0 ∈ 22 More precisely, u x G u(x)δx maynottendto 0 as k , for some exhaustion Gk . k − ∈ k kE →∞ { } P 23 Definition 4.25. Let fx = P invx denote the image of vx under the projection to in. Simi- F F 24 larly, let hx = P armvx denote the image of vx under the projection to arm. H H 25 For future reference, we state the following immediate consequence of orthogonality.
26 Lemma 4.26. With fx = P invx, fx x G0 is a reproducing kernel for in, but fx arm. F { } ∈ F ⊥ H 27 Similarly, with hx = P armvx, hx x G0 is a reproducing kernel for arm, but hx in. H { } ∈ H ⊥F 28 Remark 4.27. The role of vx in with respect , is directly analogous to role of 2 HE 2h· ·iE 29 the Dirac mass δx in ℓ with respect to the usual ℓ inner product. This analogy will be 0 30 developed further when we show that v is the image of x G under a certain isometric x ∈ 31 embedding into , in 6. It is obvious that δx , and the following result shows that HE § ∈HE 32 δ is always in span v when deg(y) < . However, it is not true that v is always in y { x} ∞ y 33 span δ , or even in its closure. This is discussed further in 6. { x} § 0 34 Lemma 4.28. For any x G , δx = c(x)vx y x cxyvy. ∈ − ∼ P 35 Proof. Lemma 4.18 implies δx, u = c(x)vx y x cxyvy, u for every u , so apply 0 h iE h − ∼ iE ∈HE 36 this to u = vz, z G . Since δx, vx , it mustP also be that y x cxyvy . ∈ ∈HE ∼ ∈HE P 40 PALLEE.T.JORGENSENANDERINP.J.PEARSE
0 4.2.1. Real and complex-valued functions on G . While we will need complex-valued 1 functions for some later results concerning spectral theory, it will usually suffice to consider 2 R-valued functions elsewhere. 3
Lemma 4.29. The reproducing kernels vx, fx, hx are all R-valued functions. 4
Proof. Computing directly, 5
1 vz, u = (vz(x) vz(y))(u(x) u(y)) = vz, u . h iE 2 − − h iE x,Xy G0 ∈ Then applying the reproducing kernel property, 6
vz, u = u(x) u(o) = u(x) u(o) = vz, u . h iE − − h iE Thus vz, u = vz, u for every u arm, and vz must be R-valued. The same compu- 7 h iE h iE ∈ H tation applies to fz and hz. 8
Definition 4.30. A sequence of functions un converges pointwise in iff k C 9 {0 }⊆HE HE ∃ ∈ such that u (x) u(x) k, for each x G . 10 n − → ∈ Lemma 4.31. If un converges to u in , then un converges to u pointwise in . 11 { } E { } HE Proof. Define wn := un u so that wn 0. Then 12 − k kE → n wn(x) wn(o) = vx, wn vx wn →∞ 0, | − | |h iE| ≤ k kE · k kE −−−−−−→ so that lim wn exists pointwise and is a constant function. 13
4.3. Relating the energy form to the Laplacian. Before completing the extension of 14 Lemma 2.13 to infinite networks, we need some definitions. 15
Definition 4.32. Let := span vx x G0 denote the vector space of finite linear combina- 16 V { } ∈ tions of dipoles, and let ∆ be the closure of the Laplacian when taken to have the dense 17 V domain . 18 V Since ∆ agrees with ∆ pointwise, we may suppress reference to the domain for ease 19 V of notation. When given a pointwise identity ∆u = v, there is an associated identity in , 20 HE but the next lemma shows that one must use the adjoint. 21
Lemma 4.33. For u, v , ∆u = v pointwise if and only if v = ∆∗ u in . 22 ∈HE V HE Proof. We show that u dom ∆∗ for simplicity, so let ϕ span vx be given by ϕ = 23 n ∈ V ∈ { } 24 i=1 aivxi . Then P n n ∆ϕ, u = ai δx δo, u = ai(∆u(xi) ∆u(o)) Lemma 4.18 h iE h i − iE − Xi=1 Xi=1 n = a (v(x ) v(o)) ∆u(x) = v(x) i i − Xi=1 n = ai vx , v (4.6) h i iE Xi=1 = ϕ, v , h iE OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 41
1 which gives the estimate ∆ϕ, u = ϕ, v ϕ v , by Schwarz, which means |h iE| |h iE| ≤ k kEk kE 2 u dom ∆∗ . The converse is trivial. ∈ V 0 3 Definition 4.34. A monopole at x G is an element wx which satisfies ∆wx = δx ∈ ∈ HE 4 pointwise. By Lemma 4.33, this is equivalent to
wx, ∆u = δx, u , for all u dom ∆ . (4.11) h iE h iE ∈ V 5 Let wo always denote the unique energy-minimizing monopole at the origin. When this 6 function exists (i.e., when wo ), it lies in in. In this case, we indicate the distin- v ∈ HE f F 7 guished monopoles wx := vx + wo and w x := fx + wo, where fx = P invx. It may happen v f F 8 that wx = w x; cf. Theorem 14.5. 9 The space of monopoles at x is denoted , and MPx
v f span wx, w x x G0 , wo := { } ∈ ∃ ∈HE (4.12) MP span vx , else. { } 10 Note that a monopole need not be in dom ∆ ; see Example 13.8 or Example 14.39. V 11 However, it is always the case that wx dom ∆∗ , which is the content of the following ∈ V 12 lemma.
13 Remark 4.35. Observe that combining (4.11) with Lemma 4.18 immediately gives
wx, ∆u = ∆u(x), (4.13) h iE 14 so that a collection of monopoles wx x G0 is a reproducing kernel for ran ∆ . Note that the { } ∈ V 15 expression ∆u(x) is defined in terms of differences, so the right-hand side is well-defined 16 even without reference to another vertex, i.e., it makes sense independent of any choice of 17 representative.
18 Remark 4.36. The presence of monopoles in is equivalent to the transience of the HE 19 underlying network, that is, the transience of the simple random walk on the network with 20 transition probabilities p(x, y) = cxy/c(x). To see this, note that if wx is a monopole, 21 then the current induced by wx is a unit flow to infinity with finite energy. It was proved 22 in [Lyo83] that the network is transient if and only if there exists a unit current flow to 23 infinity; see also [LP09, Thm. 2.10]. As mentioned in Corollary 3.29, the existence of a 24 monopole at one vertex is equivalent to the existence of a monopole at every vertex.
25 Lemma 4.37. When the network is transient, contains the spaces span v , span f , MP { x} { x} 26 and span hx , where fx = P invx and hx = P armvx. { } F H v f 27 Proof. The first two are obvious, since v = w w and f = w w by Definition 4.34. x x − o x x − o 28 For the harmonics, note that these same identities give
wv w = v = f + h = w f w + h , x − o x x x x − o x v f v f 29 which implies that h = w w . (Of course, w = w when arm = 0.) x x − x x x H 30 Corollary 4.38. arm , 0 iff there is more than one monopole at x. H 31 Proof. As usual, if this is true for any x, it is true for all. Suppose contains a monopole v v HE 32 wx , wx. Then h := wx wx is a nonzero harmonic function in . − HE 42 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Definition 4.39. If H is a subgraph of G, then the boundary of H is 1
. ∁ bd H := x H .. y H , y x . (4.14) { ∈ ∃ ∈ ∼ } The interior of a subgraph H consists of the vertices in H whose neighbours also lie in H: 2
. int H := x H .. y x = y H = H bd H. (4.15) { ∈ ∼ ⇒ ∈ } \ For vertices in the boundary of a subgraph, the normal derivative of v is 3
∂v (x) := cxy(v(x) v(y)), for x bd H. (4.16) ∂Ò − ∈ Xy H ∈ Thus, the normal derivative of v is computed like ∆v(x), except that the sum extends only 4 over the neighbours of x which lie in H. 5
Definition 4.39 will be used primarily for subgraphs that form an exhaustion of G, in 6 the sense of Definition 4.5: an increasing sequence of finite and connected subgraphs G , 7 { k} so that Gk Gk+1 and G = Gk. Also, recall that bd G := limk bd G from Defini- 8 ⊆ →∞ k tion 4.40. S P P 9
Definition 4.40. A boundary sum (or boundary term) is computed in terms of an exhaus- 10 tion G by 11 { k} := lim , (4.17) k Xbd G →∞ bdXGk whenever the limit is independent of the choice of exhaustion, as in Definition 4.6. The 12 boundary bd G is examined more closely as an object in its own right in 7. 13 § The key point of the following result is that for u, v in the specified set, the two sums 14 are both finite. The decomposition is true for all u, v by taking limits of (4.19), but 15 ∈ HE is clearly meaningless if it takes the form . 16 ∞−∞ Theorem 4.41 (Discrete Gauss-Green Formula). If u and v , then 17 ∈HE ∈ MP u, v = u(x)∆v(x) + u(x) ∂v (x). (4.18) h iE ∂Ò xXG0 xXbd G ∈ ∈
Proof. It suffices to work with R-valued functions and then complexify afterwards. By the 18 same computation as in Lemma 2.13, we have 19
1 ∂v cxy(u(x) u(y))(v(x) v(y)) = u(x)∆v(x) + u(x) (x). (4.19) 2 − − ∂Ò xX,y Gk x Xint G x Xbd G ∈ ∈ k ∈ k Taking limits of both sides as k gives (4.18). It remainsto see that one of the sums 20 →∞ on the right-hand side is finite (and hence that both are). For this part, we work just with u 21 and polarize afterwards. Note that if v = wz is a monopole, then 22
u(x)∆v(x) = u(x)δz(x) = u(z). xXG0 xXG0 ∈ ∈ This is obviously independent of exhaustion, and immediately extends to v . 23 ∈ MP OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 43
1 Remark 4.42. It is clear that (4.18) remains true much more generally than under the 2 specified conditions. Clearly, the formula holds whenever x G0 u(x)∆v(x) < . Un- ∈ | | ∞ 3 fortunately, given any hypotheses more specific than this, theP limitless variety of infinite 4 networks almost always allow one to construct a counterexample; i.e. one cannot give a 5 conditionfor which the formula is true for all u , for all networks. To see this, suppose ∈HE 6 that v = i∞=1 aiwxi with each wxi a monopole at the vertex xi. Then P ∞ u(x)∆v(x) = aiu(xi), xXG0 Xi=1 ∈ 7 and one would need to provide a condition on sequences a that would ensure ∞ a u(x ) { i} i=1 i i 8 is absolutely convergent for all u . Such a hypothesis is not likely to be usefulP (if it is ∈HE 9 even possible to construct) and would depend heavily on the network under investigation. 10 Nonetheless, the formula remains true in many specific contexts. For example, it is 11 clearly valid whenever v is a dipole, including all those in the energy kernel. We will also 12 see that it holds for the projectionsof v to in and to arm. Consequently, for v which are x F H 13 limits of elements in , we often use this result in combination with ad hoc arguments. MP ∂u
14 Lemma 4.43. For all u dom ∆ , 0 ∆u = . Thus, the Discrete Gauss-Green G bd G Ò ∈ V − ∂ 15 formula (4.18) is independent of representatives.P P
16 Proof. On each (finite) Gk in any given exhaustion,
∂u ∆u(x) + (x) = cxy(u(x) u(y)) = 0, ∂Ò − x Xint G x Xbd G xX,y Gk ∈ k ∈ k ∈ 17 since each edge appears twice in the sum; once with each sign (orientation). For the second 18 claim, we apply the formulaof the first to see that the result remains true when u is replaced 19 by u + k:
✟✟ ✟✟
(u + k)∆v + (u + k) ∂v = u∆v + u ∂v + k ∆v✟+ ∂v .
Ò Ò ∂Ò ∂ ✟ ∂ XG0 Xbd G XG0 Xbd G XG✟0 Xbd G ✟ 20 4.3.1. More about monopoles and the space . MP 21 Theorem 4.44 ([Soa94, Thm. 1.33]). Let u be a nonnegative function on a recurrent 22 network. Then u is superharmonic if and only if u is constant.
23 Corollary 4.45. If arm , 0, then there is a monopole in . H HE 24 Proof. If h arm and h , 0, then h = h h with h arm and h 0, by [Soa94, ∈ H 1 − 2 i ∈ H i ≥ 25 Thm. 3.72]. Since the hi cannot both be 0, Theorem 4.44 implies the network is transient. 26 Then by [Lyo83, Thm. 1], the network supports a monopole.
27 Definition 4.46. The phrase “the boundary term is nonvanishing” indicates that (4.18) 28 holds with nonzero boundary sum when applied to u, v , for every representative of u h iE 29 except the one specified by u(x) = u, w , for w x. h iE ∈ MP 30 Recall from Remark 4.36 that the network is transient iff there are monopoles in . HE 31 Theorem 4.47. The network is transient if and only if the boundary term is nonvanishing. 32 Moreover, the boundary term vanishes for the elements of ran ∆ . V 44 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Proof. ( ) If the network is transient, then as explained in Remark 4.36, thereisa w 1 ⇒ ∈HE with ∆w = δz. Now let wz := P inw so that for any u dom ∆ ,(4.18) 2 F ∈ V
∂wz u, wz = u(z) + u . h iE ∂Ò Xbd G
∂wz It is immediate that u = 0 if and only if the computation is done with the repre- 3 bd G ∂Ò sentative of u specifiedP by u(z) = u, wz . 4 h iE 0 ( ) Suppose that there does not exist w with ∆w = δz, for any z G . Then 5 ⇐ ∈ HE ∈ = span v as discussed in Definition 4.34. Therefore, it suffices to show that 6 MP { x}
u, vx = u∆vx, h iE xXG0 ∈ but this is clear because both sides are equal to u(x) u(o)by(4.6) and Lemma 4.20. 7 − For the final claim, note that if u ran ∆ , then (4.13) gives 8 ∈ V
∂wx ∂wx
u(x) = u, wx = u∆wx + u = u(x) + u , Ò h iE ∂Ò ∂ XG0 Xbd G Xbd G so that the boundary term must vanish. 9
Remark 4.48. It follows from Theorem 4.47 that a monopole wx cannot lie in ran ∆ . 10 V However, one can have wx ran ∆∗ , as in Example 14.39. 11 ∈ V Lemma 4.49. The network is transient if and only if there is a sequence εk with εk 0 12 1 { } → and supk (εk + ∆)− δx B < . 13 k kE ≤ ∞ 1 Proof. For both directions of the proof, we let fk := (εk + ∆)− δx. 14 1 ( ) Let ∆ be any self-adjoint extension of ∆ , and let E(dλ) be the corresponding 15 ⇒ ∗ V V projection-valued measure. Then 16
1 = + ∆ 1 = ∞ Rεu (ε ∗ )− u E(dλ)u, (4.20) V Z0 ε + λ 1 where we use the notation Rε := (ε + ∆ )− for the resolvent. Note that ∆ Rε (∆ Rε)∗ = 17 ∗ V ∗ V ⊆ ∗ V ∆ = ∆ ∆ ∆ ∆ ∆ 18 ∗ ∗ Rε∗ ∗ Rε. On the other hand, ∗ ∗ and therefore Rε ∗ Rε ∗ . Combining V V V ⊆ V V ⊆ V ∆ ∆ = ∆ 19 these gives ∗ Rε Rε ∗ . Now we apply this and (4.20) to u ∗w to get V ⊆ V λ = + ∆ 1 = + ∆ 1∆ = ∆ + ∆ 1 = ∞ fk (εk ∗ )− δx (εk ∗ )− ∗ w ∗ (εk ∗ )− w E(dλ)w. V V V V V Z0 εk + λ Note that Rε is bounded, and so w dom Rε automatically. This integral implies 20 ∈
2 2 ∞ λ 2 ∞ 2 2 fk E(dλ)w E(dλ)w = w . k kE ≤ Z0 εk + λ! k kE ≤ Z0 k kE k kE 1 Thus we have supk (εk + ∆)− δx = sup fk B = w < . 21 k kE k kE ≤ k kE ∞ ( ) We show the existence of a monopole at x. Since ε f + ∆ f = δ , the bound 22 ⇐ k k k x sup fk B implies that 23 k kE ≤ 1For concreteness, one may take the Friedrichs extension, see (B.9) but this is not necessary. See also Defini- tion 7.6 and 8.1 in this regard. § OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 45
∆ fk δx = εk fk εk B 0. k − kE k k≤ → 1 Let w be a weak- limit of fk . Then for ϕ dom ∆ , ∗ { } ∈ V
∆ϕ, w = lim ∆ϕ, fk = lim ϕ, ∆ fk = lim ϕ, δx εk fk = ϕ, δx , h iE k h iE k h iE k h − iE h iE →∞ →∞ →∞ 2 so that w is a monopole at x.
3 Lemma 4.50. On any network, ran ∆ in and hence arm ker ∆∗ . V ⊆F H ⊆ V 4 Proof. If v V, then clearly ∆ v in. To close the operator, we consider sequences ∈ V ∈ F 5 u V which are Cauchy in , and for which ∆u is also Cauchy in , and then include { n}⊆ E { n} E 6 u := lim un in dom ∆ by defining ∆ u := lim ∆ un. Since fn := ∆ un has finite support V V V V 7 for each n, the -limit of f must lie in in. E { n} F 8 The second claim follows upon taking orthogonal complements; note that in is closed, cℓ F 9 so we actually have (ran ∆ ) in. Alternatively, it can be proven directly, as follows. V ⊆ F 10 To see arm dom ∆∗ , we need h, ∆ v C v . However, this is trivially true H ⊆ V h V iE ≤ k kE 11 because h, ∆ v = 0 for any h arm, by Theorem 4.22. To see arm ker ∆∗ , h V iE ∈ H H ⊆ V 12 compute the value of ∆∗ h via ∆∗ h, v = h, ∆ v = 0. V h V iE h V iE cℓ 13 Theorem 4.51. The random walk on (G, c) is transient if and only if (ran ∆∗ ) = in. V F 14 Proof. ( ) If the network is transient, we actually have a monopole at every vertex by ⇒ 15 Lemma 3.29. Then any u span δx is in ran ∆∗ because the monopole wx is in dom ∆ , ∈ { } V V 16 and so in ran ∆∗ . The other inclusion is Lemma 4.50. F ⊆ V 0 17 ( ) If δx ran ∆ for some x G , then ∆ w = δx for w dom ∆ dom and so ⇐ ∈ V ∈ V ∈ V ⊆ E 18 w is a monopole. Then the induced current dw is a unit flow to infinity, and the network is 19 transient, again by [Lyo83].
20 The next result has a similar flavour, but instead concerns ran(I + ∆ ), as this is the V 21 orthogonal complement of the defect space; see Definition 8.13. Example 13.9 illustrates 22 an application of Lemma 4.52.
23 Lemma 4.52. Let wz be a monopole at z. Then wz ran(I + ∆ ) if and only if there ∈HE 0 ∈ V 24 is a function u satisfying u + ∆u = 0 on G z . ∈HE \ { }
25 Proof. ( ) Since wz = v + ∆v for some v , set u := v wz. Then for x , z, it is easy ⇒ ∈HE − 26 to check u(x) + ∆u(x) = 0. ( ) Set v := w + au for a := 1/(u(z) + ∆u(z)). Then it is easy ⇐ z − 27 to check that v(x) + ∆v(x) = wz(x).
28 Definition 4.53. For an infinite graph G, we say u(x) vanishes at iff for any exhaustion ∞ 29 Gk , one can always find k and a constant C such that u(x) C < ε for all x < Gk. { } k − k∞ 30 One can always choose the representative of u so that C = 0, but this may not be ∈ HE 31 compatible with the choice u(o) = 0.
32 Definition 4.54. Say γ = (x0, x1, x2,... )isa path to iff xi xi 1 for each i, and for any ∞ ∼ − 33 exhaustion G of G, { k}
k, N such that n N = x < G . (4.21) ∀ ∃ ≥ ⇒ n k 46 PALLEE.T.JORGENSENANDERINP.J.PEARSE
The next two results are almost converse to each other, although the exact converse of 1 Lemma 4.55 is false; see Figure 10 of Example 14.2. Lemma 4.55 is related to [Soa94, 2 Thm. 3.86], in which the result is stated as holding almost everywhere with respect to the 3 notion of extremal length. 4
Lemma 4.55. If u and u vanishes at , then u in. 5 ∈HE ∞ ∈F Proof. Let u = f + h vanish at . This implies that for any exhaustion Gk and any 6 ∈HE ∞ { } ε> 0, there is a k and C for which h(x) C < ε outside Gk. A harmonic function can 7 k − k∞ only obtain its maximum on the boundary, unless it is constant, so in particular, ε bounds 8 h(x) C on all of Gk. Letting ε 0, h tends to a constant function and u = f . 9 k − k∞ → 0 Lemma 4.56. If h arm is nonconstant, then from any x G , there is a path to 10 ∈ H 0 ∈ infinity γ = (x0, x1,... ), with h(x j) < h(x j+1) for all j = 0, 1, 2,.... 11
cxy Proof. Since h(x) = y x c(x) h(y) supy x h(y) and h is nonconstant, we can always find 12 ∼ ≤ ∼ y x for which h(yP1) > h(x0). This follows from the maximal principle for harmonic 13 ∼ functions; cf. [LP09, 2.1], [LPW08, Ex. 1.12], or [Soa94, Thm. 1.35]. Thus, one can 14 § inductively construct a sequence which defines the desired path γ. Note that γ is infinite, 15 0 so the condition h(x j) < h(x j+1) eventually forces it to leave any finite subset of G , so 16 Definition 4.54 is satisfied. 17
The next result is the contrapositive of Lemma 4.55, but it is instructive to prove it 18 directly. 19
Lemma 4.57. Ifh arm is nonconstant, then h has at least two different limiting values 20 ∈H at . 21 ∞ 0 Proof. Choose x G for which hx = P armvx is nonconstant. Then Lemma 4.56 22 ∈ H ∈ HE gives a path to infinity γ1 along which hx is strictly increasing. Since the reasoning of 23 Lemma 4.56 works just as well with the inequalities reversed, we also get γ to along 24 2 ∞ which hx is strictly decreasing. This gives two different limiting values of hx, and hence hx 25 cannot vanish at . 26 ∞ Remark 4.58 (The harmonicregion of a network). Inorderto see whyit is necessaryto take 27 the step of choosing x for which hx is nonconstant, consider a network built by conjoining 28 a copy of the integer lattice (Z, 1) to the binary tree ( , 1), by identifying their origins. 29 T Such a network has a portion (Z) which does not support nontrivial harmonic functions, 30 and so h must be constant for any x in this region. Since h is a reproducing kernel for 31 x { x} arm, however, h cannot be constant for x in the tree portion of the network. This idea 32 H x motivates the proof of Lemma 4.57 and the definition of boundary points given in 7.3. 33 § Roughly speaking, the “harmonic portion” of a network consists of those points which can 34 be distinguished by a harmonic function, and the boundary consists of those paths to 35 ∞ which can be distinguished by harmonic functions. We will also see in Remark 5.43 that 36 the fact that a network may have nonharmonic regions (like Z in this example) contributes 37 ha of the failure of the harmonic resistance R to be a metric. 38
4.3.2. Special applications of the Discrete Gauss-Green formula. In this subsection, we 39 use Lemma 2.13 to infinite networks to establish that ∆ is Hermitian when its domain is 40 correctly chosen (Corollary 4.61), and that Lemma 2.13 remains correct on infinite net- 41 works for vectors in span v (Theorem 4.68). 42 { x} Lemma 4.59. If v span vx , then u, v = x G0 u(x)∆v(x). 43 ∈ { } h iE ∈ P OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 47
1 Proof. It suffices to consider v = vx, whence
u(y)∆vx(y) = u(y)(δx δo)(y) = u(x) u(o) = u, vx , − − h iE XG0 XG0
2 by Lemma 4.20 and the reproducing property of Corollary 4.13.
3 Theorem 4.60. For u, v span v , ∈ { x}
u, ∆v = ∆u(x)∆v(x). (4.22) h iE xXG0 ∈ 4 Furthermore, x G0 ∆u(x) = 0 for u span vx . ∈ ∈ { } P 5 Proof. Let u span v be given by the finite sum u = ξ v . Since v is a constant, we ∈ { x} x x x o 6 may assume the sum does not include o. Then P
∆u(y) = ξ ∆v (y) = ξ (δ δ )(y) = ξ . (4.23) x x x x − o y Xx Xx
7 Now we have
u, ∆u = ξxξy vx, ∆vy = ξxξy vx,δy δo . h iE h iE h − iE Xx,y Xx,y
8 Since it is easy to compute vx,δy δo = δxy + 1 (Kronecker’s delta), we have h − iE
2 2 2 2 u, ∆u = ξxξy(δxy + 1) = ξx + ξx = ∆u(x) + ∆u(x) , (4.24) h iE | | | | Xx,y Xx Xx Xx Xx
9 by (4.23). Since u span vx , ∆u span δ x δo (see (4.23)), so that u, ∆u < and ∈ { } ∈ { − } h iE ∞ 10 (4.24) is convergent. Therefore, x ∆u(x) is absolutely convergent, hence independent of 11 exhaustion. Since P
∆v (x) = 1 1 = 0 y − xXG0 ∈ 12 by Lemma 4.20, it follows that x ∆u(x) = 0, and the second sum in (4.24) vanishes. Then 13 (4.22) follows by polarizing. P
14 Corollary 4.61. The Laplacian ∆ is Hermitian and even semibounded on dom ∆ (see V V 15 Definition B.10) with
0 ∆u(x) 2 u, ∆u < . (4.25) ≤ | | ≤ h iE ∞ xXG0 ∈
16 Proof. For u, v span v , two applications of Lemma 4.60 yield ∈ { x}
∆u, v = ∆u(x)∆v(x) = ∆u(x)∆v(x) = ∆v, u . h iE h iE xXG0 xXG0 ∈ ∈ 17 This property is clearly preserved under closure of the operator. 48 PALLEE.T.JORGENSENANDERINP.J.PEARSE
Now let u dom ∆ and choose un V with limn un u = limn ∆un 1 ∈ V { } ⊆ →∞ k − kE →∞ k − ∆u = 0. Then Fatou’s lemma [Mal95, Thm. I.7.7] yields 2 kE
2 2 ∆u(x) = lim ∆un(x) lim un, ∆un = u, ∆u , (4.26) | | | | ≤ n h iE h iE xXG0 xXG0 →∞ ∈ ∈ which gives the central inequality in (4.25) and hence semiboundedness. 3
1 2 Remark 4.62. The notation u ℓ means x G0 u(x) < and the notation u ℓ means 4 2 ∈ ∈ | | ∞ 2∈ x G0 u(x) < . When discussing an elementP u of , we say u lies in ℓ if it has a 5 ∈ | | ∞ 2 HE representative which does, i.e., if u + k ℓ for some k C. This constant is clearly 6 P ∈ ∈ necessarily unique on an infinite network, if it exists. 7
The next result is a partial converse to Theorem 4.41. 8
2 Lemma 4.63. If u, v, ∆u, ∆v ℓ , then u, v = x G0 u(x)∆v(x),and u, v dom . 9 ∈ h iE ∈ ∈ E 2 1 P Proof. If u, ∆v ℓ , then u∆v ℓ , and the following sum is absolutely convergent: 10 ∈ ∈ 1 1 u(x)∆v(x) = u(x)∆v(x) + u(y)∆v(y) 2 2 xXG0 xXG0 yXG0 ∈ ∈ ∈ 1 1 = c u(x)(v(x) v(y)) c u(y)(v(x) v(y)) 2 xy − − 2 xy − xXG0 Xy x yXG0 Xx y ∈ ∼ ∈ ∼ 1 = c (u(x) u(y))(v(x) v(y)), 2 xy − − xXG0 Xy x ∈ ∼ which is (2.9). Absolute convergence justifies the rearrangement in the last equality; the 11 1 rest is merely algebra. Substituting u in for v in the identity just established, u∆u ℓ 12 ∈ shows u dom , and similarly for v. 13 ∈ E Remark 4.64. All that is required for the computation in the proof of Lemma 4.63 is that 14 1 2 u∆v ℓ , which is certainly implied by u, ∆v ℓ . However, this would not be sufficient 15 ∈ ∈ to show u or v lies in dom . 16 E We will see in Theorem 4.55 that if h arm is nonconstant, then h + k is bounded 17 ∈ H away from 0 on an infinite set of vertices, for any choice of constant k. So the next result 18 should not be surprising. 19
2 Corollary 4.65. If h is a nontrivial harmonic function, then h cannot lie in ℓ . 20 ∈HE 2 Proof. If h ℓ , then (h) = x G0 h(x)∆h(x) = x G0 h(x) 0 = 0 by Lemma 4.63. But 21 ∈ E ∈< ∈ · since h is nonconstant, (h) >P0! P 22 E ւ 2 Remark 4.66 (Restricting to ℓ misses the most interesting bit). When studying the graph 23 2 . 2 Laplacian, some authors define dom ∆ = v ℓ .. ∆v ℓ . Our philosophy is that dom 24 { ∈ ∈ }0 E is the most natural context for the study of functions on G , and this is motivated in detail 25 in 6.1. Some of the most interesting phenomena in dom are due to the presence of 26 § E nontrivial harmonic functions, as we show in this section and the examples of 13– 14. 27 § § Consequently, Corollary 4.65 shows why one loses some of the most interesting aspects of 28 2 the theory by only studying those v which lie in ℓ . Example 13.2 illustrates the situation 29 of Corollary 4.65 on a tree network. In general, if a at least two connected components of 30 2 G o are infinite, then v < ℓ for vertices x in these components. 31 \ { } x OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 49
1 4.3.3. The Discrete Gauss-Green formula for vertices of infinite degree. If there are ver- 2 tices of infinite degree in the network, then it does not necessary follow that span δ { x} ⊆ 3 span v , or that span δ . However, we do have the following version of Theo- { x} { x} ⊆ MP 4 rem 4.41. When all vertices have finite degree, Theorem 4.68 follows from Theorem 4.41 5 by Lemma 4.28.
6 Definition 4.67. Let := span δx x G0 denote the vector space of finite linear combina- F { } ∈ 7 tions of Dirac masses, and let ∆ be the closure of the Laplacian when taken to have the F 8 domain . F 9 Note that is a dense domain only when arm = 0, by Corollary 4.23. Again, since ∆ F H 10 agrees with ∆ pointwise, we may suppress reference to the domain for ease of notation. F 11 The next result extends Lemma 2.13 to infinite networks.
12 Theorem 4.68. If uorvliesin dom ∆ , then u, v = x G0 u(x)∆v(x). F h iE ∈ P 13 Proof. First, suppose u dom ∆ and choose a sequence un span δx with un u ∈ F { }⊆ { } k − kE → 14 0. From Lemma 4.18, one has δx, v = ∆v(x), and hence un, v = x G0 un(x)∆v(x) h iE h iE ∈ 15 holds for each n. Define M := sup un , and note that M < , sinceP this sequence is {k kE} ∞ 16 convergent (to u ). Moreover, un, v M v by the Schwarz inequality. Since un k kE |h iE| ≤ · k kE 17 converges pointwise to u in by Lemma 4.31, this bound will allow us to apply Fatou’s HE 18 Lemma (as stated in [Mal95, Lemma 7.7], for example), as follows:
u, v = lim un, v hypothesis E n E h i →∞h i = lim un(x)∆v(x) un span δx n ∈ { } →∞ xXG0 ∈ = u(x)∆v(x). xXG0 ∈ 0 19 Note that the sum over G is absolutely convergent, as required by Definition 4.5. 20 Now suppose that v dom ∆ and observe that this implies v in also. By Theo- ∈ F ∈ F 21 rem 4.22, one can decompose u = f + h where f = P inu and h = P armu, and then F H
u, v = f, v + h, v = f, v , h iE h iE h iE h iE 22 since h is orthogonal to v. Now apply the previous argument to f, v . h iE 50 PALLEE.T.JORGENSENANDERINP.J.PEARSE
5. The resistance metric 1
We now introduce the natural notion of distance on (G, c): the resistance metric R. 2 While not as intuitive as the more common shortest-path metric, it reflects the topology of 3 the graph more accurately and is often more useful for modeling and practical applications. 4 The effective resistance is intimately related to the random walk on(G, c), the Laplacian, 5 and the Dirichlet energy form [LP09, LPW08, Soa94, Kig01, Str06, DS84]. 6 In 5.1, we give several formulationsof this metric (Theorem 5.2), each with its ownad- 7 § vantages. Many of these are familiar from the literature: (5.1) from [Pow76b] and[Per99, 8 8], (5.2) from [DS84],(5.3) from [DS84, Pow76b],(5.4)–(5.5) from [Kig01, Str06]. 9 § In 5.2, we extend these formulations to infinite networks. Due to the possible presence 10 § of nontrivial harmonic functions, some care must be taken when adjusting these formu- 11 lations. It turns out that there are two canonical extensions of the resistance metric to 12 infinite networks which are distinct precisely when arm , 0 (cf. [LP09] and the refer- 13 H ences therein): the “free” resistance and the “wired” resistance. We are able to clarify and 14 explain the difference in terms of the reproducing kernels for and for in, and also in 15 HE F terms of Dirichlet vs. Neumann boundary conditions; see Remark 5.19. We also explain 16 the discrepancy in terms of projections in and attempt to relate this to conditioning of 17 HE the random walk on the network; see Remark 5.22 and Remark 5.34. Additionally, we 18 introduce trace resistance and harmonic resistance and relate these to the free and wired 19 resistances. (Note: unlike the others, harmonic resistance is not a metric.) In the limit, the 20 trace resistance coincides with the free resistance. 21
5.1. Resistance metric on finite networks. We make the standing assumption that the 22 network is finite in 5.1. However, the results actually remain true on any network for 23 § which arm = 0. 24 H Definition 5.1. If one amp of current is inserted into the electrical resistance network at x 25 and withdrawn at y, then the (effective) resistance R(x, y) is the voltage drop between the 26 vertices x and y. 27
Theorem 5.2. The resistance R(x, y) has the following equivalent formulations: 28
.. R(x, y) = dist∆(x, y) := v(x) v(y) . ∆v = δx δy (5.1) { −. − } = dist (x, y) := (v) .. ∆v = δx δy (5.2) E {E . − } = dist (x, y) := min D(I) .. I (x, y) (5.3) D { ∈F } .. = distR(x, y) := 1/ minv dom (v) . v(x) = 1, v(y) = 0 (5.4) ∈ E{E } . 2 = distκ(x, y) := minv dom κ 0 .. v(x) v(y) κ (v) (5.5) ∈ E{ ≥ | − | ≤ E } 2 .. = dists(x, y) := sup v(x) v(y) . (v) 1 . (5.6) v dom {| − | E ≤ } ∈ E
Proof. (5.1) (5.2). We may choose v satisfying ∆v = δ δ by Theorem 3.28. Then 29 ⇐⇒ x − y
(v) = v(z)∆v(z) = v(z)(δ (z) δ (z)) = v(x) v(y), (5.7) E x − y − zXG0 zXG0 ∈ ∈ where first equality is justified by Theorem 2.13. 30 (5.2) (5.3). Note that every v (x, y) corresponds to an element I (x, y) via 31 ⇐⇒ ∈ P ∈ F Ohm’s Law by Lemma 3.16, and (v) = D(I) by the same lemma. Also, this current flow 32 E is minimal by Theorem 3.26. 33 OPERATORTHEORYOFELECTRICALRESISTANCENETWORKS 51
1 (5.2) (5.4). Suppose that ∆v = δ δ . Since (v + k) = (v) and ∆(v + k) = ∆v ⇐⇒ x − y E E 2 for any constant k, we may adjust v by a constant so that v(y) = 0. Define
v v(x) u := − v(x) v(y) − 3 so that u(x) = 0 and u(y) = 1. Observe that (5.1) gives (v) = v(x) v(y), whence E −
2 1 (u) = (v)/(v(x) v(y)) = (v(x) v(y))− min (u). E E − − ≥ E 1 4 This shows (v) [min (u)]− and hence dist distR. E ≤ E E ≤ 5 For the other inequality, suppose u minimizes (u), subject to u(x) u(y) = 0. Then by E − 6 Theorem 2.13 and the same variational argument as described in Remark 3.19, we have
(ρ, u) = ρ(z)∆u(z) = 0, E zXG0 ∈ 7 for every function ρ for which ρ(x) = ρ(y) = 0. It follows that ∆u(z) = 0 for z , x, y, and 8 hence ∆u = ξδ +ηδ . Observe that (u) = ( u) = (1 u), and so the same result follows x y E E − E − 9 from minimizing with respect to the conditions u(y) = 1 and u(x) = 0. This symmetry E 1 10 forces η = ξ and we have ∆u = ξδ ξδ . Now for v = u one has ∆v = δ δ , and so − x − y ξ x − y
(u) = ξ2 (v) = ξ2(v(x) v(y)) = ξ(u(x) u(y)) = ξ, E E − − 1 11 where the second equality follows by (5.1). Then ξ = (v) = (u), whence dist distR. E E E ≥ 12 (5.4) (5.5). Starting with (5.5), it is clear that ⇐⇒
2 .. v(x) v(y) distκ(x, y) = inf κ 0 . | −(v) | κ, v dom { ≥ E ≤ ∈ E} v(x) v(y) 2 = sup | −(v) | , v dom , v nonconstant . { E ∈ E } v 13 Given a nonconstant v dom , one can substitute u := v(x) v(y) into the previous line to ∈ E | − | 14 obtain
✘✘ u(x) u(y) 2✘v(x✘) v(y) 2 dist (x, y) = sup | − | | ✘−✘ | , v dom , v nonconstant κ (u)✘v(x✘) v(y) 2 { E | − | ∈ E } 1 = sup (u) , u dom , u(x) u(y) = 1 { E ∈ E | − | } = 1/ inf (u), u dom , u(x) u(y) = 1 . {E ∈ E | − | } 15 Since we can always add a constant to u and multiply by 1 without changing the energy, ± 16 this is equivalent to letting u rangeover the subset of dom for which u(x) = 1 and u(y) = 0 E 17 and we have (5.4). 18 (5.5) (5.6). It is immediate that (5.5) is equivalent to ⇐⇒
2 v(x) v(y) . sup | − | .. (v) < . ( (v) E ∞) E 2 2 1/2 19 For any v dom , define w := v/ √ (v) so that w(x) w(y) = v(x) v(y) (v)− with ∈ E E2 | − | | − | E 20 (w) = 1. Clearly then w(x) w(y) dist (x, y). The other inequality is similar. E | − | ≤ s 52 PALLEE.T.JORGENSENANDERINP.J.PEARSE
x x x