Topological entanglement entropy in d-dimensions for Abelian higher gauge theories

We compute the topological entanglement entropy for a large set of lattice models in d-dimensions. It is well known that many such quantum systems can be constructed out of lattice gauge models. For dimensionality higher than two, there are generalizations going beyond gauge theories, which are called higher gauge theories and rely on higher-order generalizations of groups. Our main concern is a large class of d-dimensional quantum systems derived from Abelian higher gauge theories. In this paper, we derive a general formula for the bipartition entanglement entropy for this class of models, and from it we extract both the area law and the sub-leading terms, which explicitly depend on the topology of the entangling surface. We show that the entanglement entropy SA in a sub-region A is proportional to log GSDA˜\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \left({GSD}_{\tilde{A}}\right) $$\end{document}, where GSDA˜\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {GSD}_{\tilde{A}} $$\end{document} is the ground state degeneracy of a particular restriction of the full model to A. The quantity GSDA˜\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {GSD}_{\tilde{A}} $$\end{document} can be further divided into a contribution that scales with the size of the boundary ∂A and a term which depends on the topology of ∂A. There is also a topological contribution coming from A itself, that may be non-zero when A has a non-trivial homology. We present some examples and discuss how the topology of A affects the topological entropy. Our formalism allows us to do most of the calculation for arbitrary dimension d. The result is in agreement with entanglement calculations for known topological models.


Introduction
The concept of entanglement entropy in quantum many-body systems is increasingly gaining relevance for both the quantum information and the condensed matter theory communities. In the latter case, the interest comes from applying ideas of quantum information that could provide new tools for the study of quantum many-body systems and, in consequence, to deepen the understanding of their quantum phases. In particular, questions about the scaling of this entropy with the system size appear to be relevant as an indicator for quantum entanglement. Of particular interest is the scaling of entanglement entropy for JHEP03(2020)167 ground states of gapped systems, since they often follow an area law [1][2][3]. More precisely, if we consider a distinguished sub-region A of the total system, the scaling of entanglement entropy is linear with the boundary of the region, ∂A. See [2] for a detailed account on the occurrence of area laws for the entanglement entropy of quantum systems.
The growing interest on the study of entanglement entropy in quantum many-body systems arises from different points of view. For example, a source of interest in the scaling area law of entanglement entropy comes from asking whether a quantum many-body system can be simulated by a classical computer. The scaling of entanglement entropy specifies how well a given many-body quantum state can be approximated by a matrix-product state or a PEPS [4]. More importantly for the purposes of this work, the topological entanglement entropy [3,5] arises as an interesting probe for topological order [6][7][8] in quantum states. The entanglement entropy calculated in the ground states of topologically ordered states follows an area law plus a universal correction indicating the presence of long-range entanglement.
Topological phases of matter are usually characterized by exhibiting long-range entanglement and non-local order parameters such as the ground state degeneracy (GSD) and topological spins. In addition, entanglement entropy turns out to be a good measure of the presence of topological order [9,10]. Details about the connection between topological entanglement entropy and topological order are for example exposed in [5,8,11]. For two dimensional topological phases, the scaling of the entanglement entropy presents a constant term [3,5] that corresponds to the topological entropy. This result is examined in detail for the Toric Code in [11][12][13][14]. In this paper, we show that this extends to Abelian higher gauge models in all dimensions. Indeed if the model is topological, it presents a non-zero topological entropy given in terms of the higher cohomology groups of both the bulk and the boundary ∂A of the subregion A.
In general, two dimensional topological order and topological entropy are relatively well understood. The same cannot be said about higher dimensional topological phases. Our interest here is to shed some light on the main features of topological entanglement entropy in dimensions higher than two. To do so, we restrict ourselves to study Abelian lattice models that come from higher gauge theories, such as the ones studied by [15][16][17][18][19][20][21][22][23]. Some of these models are topological. They present a ground state degeneracy that depends on the topology of the underlying manifold, and extended excitations generalizing anyons to higher dimensions. The models presented here can be interpreted as higher gauge generalizations of the Toric Code. This simplification allows us to work in arbitrary dimensions and to understand in detail how these models depend on the topology of the underlying manifold. In the present work we look at the topological entropy and how it depends on the topology of the subregion A.
The definition of entanglement entropy is straightforward: we consider a bipartition of the system into a sub-region A and its complement (B). Let ρ be the density matrix of the ground space state, defined in the whole lattice. Then, ρ A = Tr B (ρ) is the reduced density matrix, obtained by tracing out the contribution from region (B). The entanglement entropy is then defined as the von Neumann entropy of the reduced density matrix, namely, (1.1)

JHEP03(2020)167
In a gapped phase, the entanglement entropy is expected to satisfy an area law as the leading term. The topological information is contained in subleading terms and, in general, it is not easy to extract. Several prescriptions [3,5] were constructed in order to extract the topological correction to the entanglement entropy in two dimensional gapped systems. These prescriptions have been generalized [24,25] for d = 3 and, consequently, used to successfully obtain the entanglement entropy of fracton models [26,27].
In this paper, we study the entanglement entropy of n-dimensional Abelian higher gauge theories, all at once, without the need to adapt the procedure for each dimension. This can be achieved by using the language of homological algebra, in which higher gauge theories are naturally described as shown in [19]. The way the entanglement entropy is obtained, in essence, relies on the fact that the lattice models are constructed as stabilizer codes [28]. Similar calculations for stabilizer codes are performed in [12,[29][30][31][32]. The result obtained for the entanglement entropy relates this quantity to the ground state degeneracy of an associated model, which we write as S A = log GSDÃ , (1.2) where GSDÃ is the ground state degeneracy of a particular restriction of the original model, to be defined in section 4.1. We show that this result allows us to write S A as a sum of two terms. The first one S ∂A is in agreement with the area laws depending on the geometry of the boundary ∂A, while the second term S Topo is explicitly topological, depending on the topology of both the bulk of the region A and its boundary ∂A. This result holds for any higher gauge theory in the sense of [19] and for any dimension d. We explicitly calculate examples where the topological term depends not only on the Betti numbers of the subregion A and its boundary ∂A, but also on more exotic properties. The paper is structured as follows. In section 2 we begin by reviewing Abelian higher gauge theories in detail, this is done by giving explicitly an example. Next, in section 3 we show how these models are easily described in the language of homological algebras, as described in detail in [19]. In section 4 the calculation of the entanglement entropy is performed, and we show how to extract both the area law and the topological entropy from it. In section 5 we apply the results of section 4 to examples in 2D and 3D. We end the paper with some final remarks in section 6.
2 Models from Abelian higher gauge theories. The warm-up example Before we define our class of models in full generality, it is convenient to first present an examples in 3 dimensions. It will be described in the usual way as a many-body system defined on a lattice. The Hilbert space H consists of quantum states attached to elements of the lattice such as vertices, links, and plaquettes. A Hamiltonian H acting on H completes the picture.
In section 3 we will introduce a formalism that allows us to describe all models in this class in a unified way for all dimensions; this is made possible by employing a few constructions coming from homological algebra. As far as the present section, we don't need to be concerned with all the homological details, but we will point out some of the JHEP03(2020)167 chain complexes that will be part of the construction presented in section 3. It is worth mentioning again that the example exhibited in this section is merely illustrative. We choose to present the model for a 3-dimensional cubic lattice to be more intuitive. Through our formalism, to be presented in section 3, these models can be defined and studied in any dimension d for more general cell decompositions, including simplicial complexes.

0, 1, 2-gauge in 3D
This model is a generalization of an ordinary gauge theory. In addition to the states localized on the links and labeled by an Abelian group G 1 , we also have states associated with plaquettes that are labelled by another Abelian group G 2 , and states on vertices associated to the Abelia group G 0 . Here we use the additive notation for the group operation. When G 0 , G 2 are trivial we recover the Quantum Double model defined in 3-dimensions. When G 0 = 0, this class of models corresponds to an Abelian version of the models constructed in [17,18], where 2-groups were considered. We will leave the discussion about the relation between our formalism and the one based on 2-groups (or more generally, n-groups) to section 3. Here we are working with a particular case of Abelian 2-groups, where the only data we need is expressed by the group homomorphism ∂ G 2 : G 2 → G 1 . For simplicity, we consider a 3-dimensional space Σ discretized by a cubic lattice. Let K 0 , K 1 , K 2 and K 3 be the sets of vertices, links, faces and cubes, respectively. For each vertex v ∈ K 0 we have a Hilbert space H v with orthornomal basis {|h , h ∈ G 0 }. For each link l ∈ K 1 we have a Hilbert space H l with orthornomal basis {|g , g ∈ G 1 }. For each plaquette p ∈ K 2 we have a Hilbert space H p with orthornomal basis {|α , α ∈ G 2 }. We call this kind of model a 0, 1, 2-gauge theory to indicate that there are quantum states associated with sets K 0 (vertices), K 1 (links) and K 2 (plaquettes). The total Hilbert space is the tensor product over all local Hilbert spaces, namely It is convenient to represent a basis element of H as functions (f 0 , f 1 , f 2 ), where f i : K i → G i , for i = 0, 1, 2. We will use a notation where each element of the lattice will have a group element attached in order to represent the functions f i . The Hamiltonian has the form and is composed by a sum of commuting projectors labelled by vertices, links, plaquettes and cubes. The vertex operator A v is a sum of gauge transformations. In other words, (2.3) Higher gauge transformations appear in the definition of A l : JHEP03(2020)167 Figure 1. The action of A β l is shown on an initial arbitrary basis state, involving the link l ∈ K 1 and its four adjacent plaquettes.
The elementary transformation A β l acts on the adjacent plaquettes with β ∈ G 2 and on the link with (∂ G 2 β), as shown in figure 1. The 0-holonomy operator, B l , compares the gauge fields of adjacent vertices with the map ∂ G 1 applied to the link degree of freedom, namely The 1-holonomy is not exactly the same as in the Quantum Double models: the operator B p is defined by Note that this operator gives eigenvalue 1 when the holonomy of the plaquette is equal to ∂ G 2 α. To distinguish it from the ordinary holonomy, this operator is usually called fake holonomy. Here we will use the name 1-holonomy instead.
To simplify our notation, we will denote ∂ G 2 α simply by ∂α whenever there is no danger of ambiguity. The 2-holonomy operator B c measures the sum of the degrees of freedom living on the plaquettes that compose the boundary of the cube c, with eigenvalue 1 when the result of the sum equals zero and eigenvalue zero otherwise, as defined in equation (2.6): where o j = {0, 1} takes into account the relative orientation of the faces with respect to the cube, similarly to the 2-holonomy operator. A discussion about the topological properties of this model, such as ground state degeneracy (GSD), will be postponed until we present the general formalism. The data needed to cast this model into the formalism of section 3 is as follows. One needs to specify a pair of chain complexes. The first one represents the discretization of the manifold and is given by

JHEP03(2020)167
where C i is generated by K i and ∂ i is the boundary map, for i = 1, 2, 3. The algebraic data is also a chain complex of the form The groups that label the degrees of freedom are chosen to be G 0 = G 1 = Z 2 = {0, 1} and G 2 = Z 4 = {0, 1, 2, 3}. Moreover, the homomorphism that relates 1, 2-gauge fields is chosen to be such that ∂ G 2 (1) = 1. As in the previous example, we use a graphical description for the basis states as follows: • a dotted line through a plaquette when it holds a |2 p state, and an oriented dashed line for the |1, 3 p states, • for the links we picture a gray surface orthogonal to it for |1 l state, no surface for the |0 l state, • a gray volume for vertices in |1 v state.
The flat configurations of the theory are those that are invariant under the actions of B l , B p and B c for all link l, plaquettes p and cubes c of the lattice. Such configurations consist of dotted or dashed loops, conditions that are enforced by B c . The 1-holonomy operator B p implies that every dashed line encloses a gray surface, which is the boundary of a gray volume. If Σ has the topology of a 3-dimensional ball, there is only one ground state given by For manifolds with exotic topologies, the number of ground states depends on the number of non-equivalent non-contractible closed surfaces one can drawn over Σ. The topological nature of this model will become clearer when we introduce the general formalism in 3.
Now that we are familiar with some examples of Abelian higher gauge theories, we are ready to describe them all at once using a more general mathematical structure. The more general framework that is going to be exhibited in the next section also allows to compute the entanglement entropy, in the more general case, as we will show in section 4.

Review of Abelian higher gauge theories
The model presented on section 2 is an example of what we call Abelian higher gauge theories introduced in [19]. It is useful to describe them using a formalism borrowed from Homological Algebra since it allows us to handle a large class of models of arbitrary dimensions. In this section we recall from [19] only the basic notation and results needed to calculate the entanglement entropy for any such model. We refer to [19] for further details. We have indicated in the last section that the models are parameterized by two chain complexes. The first one is geometrical in nature and accounts for the structure of the JHEP03(2020)167 lattice. As for the second, it is a chain complex of finite Abelian groups encoding the higher gauge group of the model. There will be one model (Hilbert space and Hamiltonian) for any such choice of chain complexes. The choices corresponding to the example was given in section 2.
A simplicial decomposition is a natural choice for lattices of any dimension. Although the formalism can accommodate for any finite cell decomposition we will assume that the lattice K is made of simplices. In other words where K n is the (finite) set of n-dimensional simplices. We would like to point out that there are no further assumptions on K, which makes the formalism very flexible. For instance, K may have a boundary and may not have a uniform dimension It is a standard procedure [33] to associate to K a chain complex which we will denote by (C(K), ∂ C ), where 0 represents the trivial group. We recall that C n is the Abelian group freely generated by K n . In other words, if we write the group operation as an addition operation, c ∈ C n is given by a formal linear combination with n(x) ∈ Z. The homomorphisms ∂ C n : C n → C n−1 are the usual boundary maps. To describe the higher gauge groups that label the degrees of freedom in the simplicial complex, we introduce a chain complex (G, ∂ G ) of finite Abelian groups given by where 0 denotes the trivial group and ∂ G n : G n → G n−1 are group homomorphisms such that ∂ G p • ∂ G p+1 = 0, for any 0 ≤ p ≤ d. Note that the d = 2 case, i.e., the chain can be recast into the language of strict 2-groups or equivalently, crossed modules, which were used in [17] to construct models of topological phases in 3 + 1 dimensions based on higher gauge theories. A crossed module is a quadruple G = (G 1 , G 2 , ∂, ), where G 1 and G 2 are groups, ∂ : G 2 → G 1 is a group homomorphism and : G 1 × G 2 → G 2 is an action of G 1 on G 2 , which satisfies the following conditions: Since we consider only Abelian groups, the chain (3.4) defines a crossed module with trivial action. In fact, the general case where the chain complex (G, ∂ G ) is composed of d + 1 groups can be reformulated in the language of (Abelian) strict (d + 1)-groups. Figure 2. A configuration f ∈ hom(C, G) 0 , consisting on a collection of homomorphisms {f n }. Now, we define a gauge configuration f to be an assignment of a group element g ∈ G n for each element x ∈ K n . In other words, a gauge configuration is a sequence f = {f n } d n=0 of functions such that

JHEP03(2020)167
x → f n (x). (3.8) Strictly speaking, we should call f a higher-gauge configuration. Only in the case when all groups except G 1 are trivial f is a proper gauge configuration, as can be seen from the examples of last section. For simplicity, we will keep using "gauge configuration" to mean a generic f . Because C n is freely generated by K n , each map f n in (3.7) defines a unique group homomorphism f n : C n → G n , which is given by the extension of f n by linearity, i.e., if c ∈ C n as in (3.2), then f n (c) = x∈Kn n(x)f n (x).
We use the same name f n to denote a gauge configuration as in 3.7 and a homomorphism f n : C n → G n since there is a one to one correspondence between them. The set Hom(C n , G n ) of homomorphisms is also an Abelian group if we set (f n +f n )(x) = f n (x) +f n (x), f n ,f n ∈ Hom(C n , G n ).
It is useful to collect all such Abelian groups in a single direct sum. This simple fact allows us to view a gauge configuration f as an element of the direct sum of Abelian groups. Thus, a gauge configuration can be represented by a collection of maps between chain complexes, as depicted by the diagram in figure 2. We would like to point out that figure 2 is not a commuting diagram. When this happens, f ∈ hom(C, G) 0 is called a chain map and, as we will see, the corresponding gauge configuration is gauge equivalent to the trivial one.
We are now in position to define the model by providing a Hilbert space H of states |ψ and a Hamiltonian operator H acting on H. The set hom(C, G) 0 is finite. Let H be the complex vector space generated by the orthonormal basis {|f } f ∈hom(C,G) 0 . In other JHEP03(2020)167 words, a state ψ ∈ H is given by a linear combination where ψ(f ) ∈ C, and the internal product is f |g = δ(f, g). We can explicit the orthonormal basis {|f } as a tensor product over the local degrees of freedom living on the simplexes labeled by the groups.
where the tensor product is made over all the 0-simplexes (vertices v), 1-simplexes (links l), 2-simplexes (plaquette p),. . . through to the d-simplexes associated to a group of the chain (G, ∂ G ). In order to define the Hamiltonian, we will need to introduce more groups other than hom(C, G) 0 given by (3.9). Let us consider the groups hom(C, G) p defined by Hom(C n , G n−p ). (3.12) An element g ∈ hom(C, G) p is a sequence {g n } d n=0 of homomorphisms g n : C n → G n−p . The example of hom(C, G) 1 is shown in figure 3.
An important observation is that the sequence of groups hom(C, G) p can be made into a co-chain complex. This is achieved by considering maps δ p : hom(C, G) p → hom(C, G) p+1 , defined by: with h ∈ hom(C, G) p . In fact, it is straightforward to verify that δ p+1 • δ p = 0, which turns the sequence into a co-chain complex. The expression above shows only the part of the sequence that is relevant for the present application, please refer to [19] for a more detailed account. Associated to this co-chain complex, there are the so-called Brown cohomology groups [34] H p (C, G) = ker(δ p )/Im(δ p−1 ), (3.15) JHEP03 (2020)167 which, as shown in [34], are isomorphic to the direct product of the cohomology groups of the chain complex (C(K), ∂ C ) with coefficients in the homology groups of (G, ∂ G ), i.e., From (3.16), it is clear that the cohomology groups defined in equation (3.15) are topological invariants of the manifold described by the chain complex (C(K), ∂ C ). The Brown cohomology plays a major role in our formalism, as we will see in the following sections. We refer the reader to [19,34] and references therein for a more detailed account on Brown cohomology and on the isomorphism (3.16).
To complete the description, we need to define a chain complex that is the dual of (3.14). This will be done by dualizing the groups G n as follows. Let Hom(G n , U(1)) be the set of homomorphisms a : G n → U(1), for each 0 ≤ n ≤ d. Since each G n is Abelian, this is nothing but the set of irreducible unitary representations of G n , denoted by G n . Let us give G n a structure of an Abelian group. Let a, b ∈ G n and g ∈ G n . Let us write the group operation in G n as a + b and the inverse of a as −a. The group is defined by setting (a + b)(g) = a(g)b(g) and (−a)(g) = (a(g)) −1 . In order to dualize (3.14) we first define the dual hom(C, G) p of (3.12) as As before, an element m ∈ hom(C, G) p is a sequence {m n } d n=0 with m n ∈ Hom(C n , G n−p ). Each m n is completely defined by its values on the generators x ∈ K n . This allows us to introduce a pairing given by Let us define a boundary map δ p : hom(C, G) p → hom(C, G) p−1 given by where m ∈ hom(C, G) p and f ∈ hom(C, G) p−1 . Clearly, δ p • δ p+1 = 0 and thus the chain complex dual to (3.14) that we will need is given by JHEP03(2020)167

Operators and Hamiltonian
The Hamiltonian we presented as an example in section 2 have a similar structure. It is a sum of operators that can be divided in two types. There are higher gauge transformations and diagonal operators measuring higher holonomies. In the general formalism, these two sets of operators in H come from the co-chain complex (3.14) and the chain complex (3.21), respectively. The first set of operators is parametrized by hom(C, G) −1 whereas the second one by hom(C, G) 1 . For t ∈ hom(C, G) −1 and m ∈ hom(C, G) 1 we define: The interpretation of (3.22) and (3.23) can be derived from the special case when the chain complex (3.3) is made of trivial groups except for G 1 . The resulting model has the familiar form of a gauge theory on the lattice. In this case, it follows that a configuration f ∈ hom(C, G) 0 assigns one group element of G 1 for each link of the lattice, as expected in a ordinary gauge theory. Furthermore, t ∈ hom(C, G) −1 gives a group element g(v) ∈ G 1 for each vertex v of the lattice. One can verify that A t performs on each vertex v an ordinary gauge transformation with parameter g(v). As for the general case, f and A t define what we mean by higher gauge configurations and higher gauge transformations. As we pointed out before, we will keep calling them gauge configurations and gauge transformations. Going back to the special case, we need to look at the eigenvalue m, δ 0 f of B m to see what is it measuring. It follows from the definition that δ 0 f ∈ hom(C, G) 1 and m ∈ hom(C, G) 1 , it assigns to each face p of the lattice its holonomy h p and a unitary representation χ p ∈ U(1) respectively. Operators B m are therefore measuring the number p χ p (h p ) depending on all holonomies of the lattice.
The first thing to be noticed is that both A t and B m are not localized as they act on the entire lattice. For the definition of the Hamiltonian, however, we need to define local projectors. This is easily achieved by taking t and m with a local support in the lattice K and averaging over the groups. Definition 3.1 (Localized maps). Let x ∈ K n , g ∈ G n+1 and r ∈Ĝ n−1 . We define the local mapsê[n, x, r] ∈ hom(C, G) 1 and e[n, x, g] ∈ hom(C, G) −1 by where y ∈ K, and f ∈ hom(C, G) p .

Definition 3.2 (Local projector operators)
. Let x ∈ K n , g ∈ G n+1 , r ∈Ĝ n−1 . We define local gauge projector A n,x and local holonomy projector B n,x as: (3.27)

JHEP03(2020)167
The Hamiltonian operator H is defined as It is straightforward to show that A n,x and B n,x are commuting projectors. Furthermore, in the special case where the chain complex (G, ∂ G ) has only G 1 different from the trivial group, we recover the quantum double model with group G 1 . Also, by choosing (G, ∂ G ) we can reproduce the example we have discussed in the last section.
This Hamiltonian is actually frustration free since there is at least one state that gives eigenvalue 1 for all A n,x and B n,x . Let |0 denotes the state labeled by the trivial element of the group hom(C, G) 0 . It corresponds to a configuration that maps all elements of K n to 0 ∈ G n for all n. Let us define One can show that |0 G is non zero and for all 0 ≤ n ≤ d and x ∈ K n . Therefore, a state |ψ is in the ground state H 0 if and only if for all 0 ≤ n ≤ d and x ∈ K n . It is useful to characterize H 0 in another way. Let us consider the following operators: 1. projector A 0 given by that maps any state |f ∈ H into a normalized sum of gauge equivalent states; 2. projector B 0 given by which gives eigenvalue 1 for a state |f ∈ H, only if satisfies f ∈ ker(δ 0 ); in other words, it projects onto the flat holonomy sector of H.
As stated in [19], the projector Π 0 on the ground state subspace H 0 can be written as

JHEP03(2020)167
Furthermore, the dimension of H 0 is determined by the zeroth Brown cohomology group H 0 (C, G) of the cochain complex in eq. (3.14), which explicitly shows how the ground state degeneracy depends on the topology of the underlying manifold. This result can be stated more precisely as follows: Proof. The proof of this theorem can be found in [19] as its main result.

Example
To see how this general framework works, let us review the example given in section 2, this time built from the formalism presented in section 3.

0, 1, 2-gauge
The 0, 1, 2-gauge model presented in section 2.1 comes from the chain complex (C, ∂ G ) given in equation (2.7). We show the relevant maps in figure 4, in particular: • Classical gauge configurations now consider degrees of freedom on plaquettes given by maps f 2 ∈ Hom(C 2 , G 2 ), in addition to the link configurations defined by f 1 ∈ Hom(C 1 , G 1 ), and the configurations for the vertices f 0 ∈ Hom(C 0 , G 0 ). In figure 4 they are represented as straight lines.
• The generalized notion of gauge transformations include 1-gauge transformations given by maps t 0 ∈ Hom(C 0 , G 1 ) and 2-gauge transformation, from every link to its neighbor plaquettes t 1 ∈ Hom(C 1 , G 2 ). In figure 4 they are represented as skewed dotted lines.
• We have as holonomy maps: the 1-holonomy with the functions m 1 ∈ Hom(C 1 ,Ĝ 0 ), the 1-holonomy which is measured by maps m 2 ∈ Hom(C 2 ,Ĝ 1 ) and the 2-holonomy measured by m 3 ∈ Hom(C 3 ,Ĝ 2 ). They are nt represented in the figure, but they would be skewed in the opposite direction of the gauge transformations.
So the Hamiltonian in equation (2.2) is obtained by the decomposition for the maps: (3.37) Now that we have familiarized with the general theory with the help of the example, we are ready to proceed onto the next section, where we present the main results of this JHEP03(2020)167 paper. The calculation of the entanglement entropy for all Abelian higher gauge theories can be carried out in general. The result obtained relates the entanglement entropy of an Abelian higher gauge theory to the GSD of a related theory, as we will precisely see. Furthermore we also compute the topological contribution to the entanglement entropy.

Entanglement entropy in Abelian higher gauge theories
In this section we calculate the entanglement entropy for the class of models defined in [19] and reviewed in section 3. We begin by defining the bipartition of the C(K), ∂ C chain complex into a subcomplex C(K A ), ∂ C A and its complement. We then observe that an associated higher gauge theory can be defined in the subcomplex (C(K A ), ∂ C A ) which will be useful for both the calculation and the interpretation of the results. As usual, we begin by introducing the density matrix ρ in terms of the ground state projector of (3.35). The reduced density matrix ρ A = Tr B (ρ) is then obtained and shown to be best written in terms of the local operators of the higher gauge theory defined in the subcomplex C(K A ), ∂ C A . The entanglement entropy is the von Neumann entropy of the reduced density matrix The result we obtain relates this quantity to a restricted gauge theory in region A. In particular, we show that the entanglement entropy of a higher gauge theory with Hamiltonian as in (3.28) is equal to the logarithm of the ground state degeneracy GSDÃ of a related higher gauge theory restricted to region A, in other words We further analyze the result to extract the topological information. As we will explain in section 4.43, the entropy S A has two terms: where S ∂A scales with the size of the boundary and S Topo is a constant contribution depending on the topology of A and its boundary ∂A.

Bipartition of the geometrical chain complex
We recall from section 3 that the geometrical content of the model is given by the chain complex (C(K), ∂ C ) encoding the lattice. We consider a simplicial chain complex for JHEP03(2020)167 convenience. In order to calculate the entanglement entropy we first need to define the bipartition of the lattice. The system is divided into two regions A and B, where A is the region we have access, as in [27,29,30]. We split the simplicial complex K = d n=0 K n into a subcomplex K A of dimension d by choosing a smaller subset K d,A ⊂ K d together with their subcomplexes of smaller dimension. Then K A is a subcomplex as the boundary maps are well defined, i.e. their image belong to the subcomplex K A . For each 0 ≤ n ≤ d dimension the sets of n-simplices is divided in the form K n = K n,A ∪ K n,B , where K n,A is the set of n-simplices in region A and K n,B the set of simplices in region B. We do this in such a way that K A = d n=0 K n,A is a subcomplex of K. Note that in general the complement K B is not a simplicial complex on its own, as its boundary maps may have image in K A .
Let C n,A be the n-chain group generated by the n-simplices, x ∈ K n,A , of region A. Let also ∂ C n,A : C n,A → C n−1,A be the restriction of the boundary map ∂ C n into the subset Let us apply the construction reviewed in section 3 to the C(K A ), ∂ C A complex together with the same chain complex of Abelian groups in (3.3), namely Homomorphisms between the two chain complexes (4.2) and (4.3) can be constructed giving rise to the groups hom(C A , G) p := n Hom(C n,A , G n−p ).
Elements of such groups are sequences of morphisms f n,A : C n,A → G n−p whose support lies on K A . For example, a gauge configuration on region A is an assignment of a group element g ∈ G n for each element x ∈ K n,A . This is, a collection of maps f A = {f n,A } for n = 0, 1, . . . , d, where: Therefore, gauge configurations in A can be viewed as elements of the group The set of vectors {|f A }, labeled by group elements f ∈ hom(C A , G) 0 , form a basis of the Hilbert space H A . In other words, a state |Ψ A ∈ H A is written as

JHEP03(2020)167
Similarly, we have the group hom(C A , G) −1 whose elements serve as parameters for the higher gauge transformations as well as the dual group hom(C A , G) 1 that parametrizes the higher gauge holonomy operators. More importantly, a higher gauge theory in region A can be defined by considering a chain complex similar to the one in (3.14), that is, where the co-boundary map δ p

Reduced density matrix
As usual, we start by introducing the density matrix ρ of the model with Hamiltonian (3.28), given by where Π 0 : H → H 0 is the ground state projector of eq. (3.35) and GSD stands for the ground state degeneracy of eq. (3.36). We start from a product state, a linear combination of the ground space states, which are independent states. From (3.35) we know that the ground state projector Π 0 can be written in terms of the projectors in (3.33) and (3.34) as (4.6) However, we want to re-parametrize the two sums in the above equation such that they run over non trivial elements only. In other words, we want to factor the redundancies out of the sums, reframing the sum as sum over classes. This can be achieved by looking at the group structure of hom(C, G) −1 and hom(C, G) 1 . Take for instance hom(C, G) −1 whose elements parametrize the higher gauge transformations of the theory. The redundancies in the sum over t ∈ hom(C, G) −1 of (4.6) come from elements that act trivially over quantum states (examples of such elements are shown in section 5). Recall that gauge transformations act on actual states by means of the δ −1 operator. Thus, we can identify the elements of hom(C, G) −1 that act trivially on states: they form a subgroup of hom(C, G) −1 called the kernel and given by ker(δ −1 ) := {t ∈ hom(C, G) −1 | δ −1 (t) = 0}, where 0 ∈ hom(C, G) 0 is the identity element that labels the trivial gauge configuration. Morever, non-trivial gauge transformations are parametrized by elements of hom(C, G) −1 that are not mapped to the identity by δ −1 , they define a subgroup of hom(C, G) 0 known as image and denoted Im(δ −1 ). Both the kernel and the image of the co-boundary map, δ −1 , are related to each other by the first isomorphism theorem [35] which in this case reads hom(C, G) −1 ker(δ −1 ) Im(δ −1 ). (4.7)

JHEP03(2020)167
Elements of the quotient group in the above expression are the cosets of ker(δ −1 ) in hom(C, G) −1 . This is: where the coset [t] = {t + h i , h i ∈ ker(δ −1 )} consists on all elements of hom(C, G) −1 that differ from t by an element in ker(δ −1 ). This is precisely what we need to factor the sums in (4.6). The sum over t ∈ hom(C, G) −1 can be replaced by a sum over the cosets of ker(δ −1 ) in hom(C, G) −1 as follows: where s ∈ [s] is an arbitrary representative of the coset. A similar argument holds for the sum over m ∈ hom(C, G) 1 which allows to factor out the redundancies from the second sum in (4.6). By doing this, we ensure that the sums run over independent group elements only: note that we have used |hom(C, G) −1 | = |ker(δ −1 )||Im(δ −1 )| to simplify the normalization factor of the first sum. A similar identity holds for the second sum. This leaves us with the density matrix of (4.5) written as: We can now proceed to the calculation of the reduced density matrix. Let us consider the bipartition of the geometric chain complex C(K), ∂ C described in section 4.1. This procedure splits the Hilbert space into two subspaces H = H A ⊗ H B . Then, we obtain the reduced density matrix by taking the partial trace over region B, this is: (4.10) To evaluate the partial trace we consider a basis {|f n,B }, where f n,B ∈ hom(C B , G) 0 is the restriction of the tensor product 3.11 over simplexes only belonging to B. Therefore we can write the basis dividing the tensor product: where |f A is the tensor product over the elements belonging to the subcomplex A, the same for complementary B. For simplicity, let us denote its basis as {|b i }, with i = JHEP03(2020)167 1, 2 . . . , dim(H B ). The reduced density matrix is now written: and [m] ∈ hom(C,G) 1 ker(δ 1 ) . Both A t and B m are traceless operators unless they are equal to the identity operator, as we show in proposition A.1 of appendix A. Because of this, the only terms that survive the partial trace are those for which A t and B m act as identity in H B . In other words, the only operators that survive the partial trace are those that act exclusively on H A . In terms of local operators this means that all holonomy operators labeled by simplexes that belong to region A will survive the trace. The case of gauge transformations is more subtle since we need to discard gauge transformations labeled by elements that lie at the boundary of A as well, see appendix A.
To account for such operators we use the restricted gauge theory defined by the complex in (4.4) with a slight modification. We are interested on further restricting such a theory by allowing gauge transformations that act on the interior of K A only. In other words, we want to discard the gauge transformations at the boundary ∂A. To this intent, let us give a more precise notion of interior of A. Let K n,Ã = {x ∈ K n,A | x ∩ ∂A = ∅} be the set of all n-simplices that have no intersection with the boundary of A. Then, the interior of A is the setÃ := d n=0 K n,Ã . Recall from section 3 that the higher gauge transformations of a higher gauge theory are parametrized by elements of the group hom(C, G) −1 . Thus, to account for gauge transformations that act exclusively on the interior of A we just need to consider the subgroup of hom(C A , G) −1 whose support is contained onÃ only. This can be done using the notion of interior of A as follows: consider homomorphisms whose support lie on the interior of A, namely Hom(C n,Ã , G n+1 ). Let then: Hom(C n,Ã , G n+1 ), with elements f ∈ hom(CÃ, G) p consisting on collections of maps f = {f n }: 13) x → f n (x), (4.14) where x ∈ K n,Ã and f n (x) ∈ G n+p . It is straightforward to show that hom(CÃ, G) p is a subgroup of hom(C A , G) p . Moreover, we can define the restriction of the co-coundary operator δ p into the interior of A. This is, δ p A := δ p |Ã, such that the sequence is a co-chain complex, i.e., δ 0 A • δ −1 A = 0. This co-chain complex encodes an Abelian higher gauge theory over K A whose gauge transformations are restricted to act on the interior of A only.

JHEP03(2020)167
We can now return to equation (4.12) and evaluate the partial trace of the density matrix ρ, which yields where the sums now run over independent internal gauge transformations and

non-trivial holonomy values in
Finally, observing that Tr(1 B ) = dim(H B ) we get for the reduced density matrix: The above expression can be further simplified as follows: observe that by applying the first isomorphism theorem [35] on the sequence of eq. (3.14) it is easy to show that the dimension of the Hilbert space factors into dim(H) = hom(C, G) 0 = ker(δ 0 ) Im(δ 0 ) .
Also, in appendix B we show that Im(δ 0 ) = |Im(δ 1 )| this allows us to write: which in turn yields for the reduced density matrix,

Entanglement entropy
Having found the reduced density matrix in (4.19) we are able to calculate its Von Neumann entropy, also known as entanglement entropy. This calculation will require us to evaluate the logarithm of ρ A at some point and this is usually done using a series expansion. In this sense, we will start by calculating the square of ρ A : where in the last equality the factors in the numerator come from rearranging the sums over p ∈ hom(CÃ,G) −1 ker(δ −1 A ) and over q ∈ hom(C A ,G) 1 ker(δ 1 | A ) . This leaves for the square of the density matrix: (4.20)

JHEP03(2020)167
Now we can calculate the logarithm of ρ A by series expansion, which yields: log(ρ A ) = log(λ) λ ρ A . Finally the entanglement entropy is: where we have used Tr(ρ A ) = 1. Let us look at the λ factor more carefully, since it encodes the essential information about the entanglement entropy of the model. By recalling that dim(H A ) = hom(C A , G) 0 = ker(δ 0 A ) Im(δ 0 A ) , we are able to write: Equation (4.22) is already very interesting since it relates 1/λ to the ground state degeneracy (GSD) of the model restricted to H A and for which gauge transformations act in the interior of A only. By replacing this expression into eq. (4.21) we are able to state our first result, that the entanglement entropy is given by: We want to highlight that the only requirement we asked for the bipartition is that the simplicial complex K is divided into a subcomplex K A and its complement. Therefore, this result is very general since it is valid for any higher gauge theory of the type described in sections 2, 3 and constructed in [19] and for any arbitrary dimension.

Topological entanglement entropy
Having a general result for the bipartition entanglement entropy S A , given by equation (4.23), it is natural to ask if it is possible to extract from it both the area law and the sub-leading, possibly topological terms, explicitly exhibiting its dependency on both the geometry and the topology of the region A. Here we give an answer to this question and show how S A depends on A. We demonstrate that the topological contribution to the entanglement entropy comes from topological invariants of the sub-region A and from topological invariants of the entangling surface ∂A. These invariants are related to Brown's cohomology groups [34]. One can express the cohomology groups of Brown as a product of all cohomology groups of the manifold in question with coefficients in the homology groups of the chain complex of groups (4.3) (see [34] for further details). Thus, for higher gauge theories, the topological entanglement entropy depends on the Betti numbers of A and ∂A, but the actual relation is more involved than the case of 1-gauge theories. Let's consider equation (4.23). It depends on the quantity GSDÃ, the ground state degeneracy of a model restricted to the region A without gauge transformations at the boundary ∂A. This quantity can be rewritten in the following way: multiplying and dividing it by |Im(δ −1 A )|, i.e., the order of the group image of δ −1 A , we have

JHEP03(2020)167
where, from theorem 3.3, GSD A = |ker(δ 0 A )|/|Im(δ −1 A )| is the ground state degeneracy of the model restricted to region A. As explained in section 3, GSD A = |H 0 (C A , G)|, where |H 0 (C A , G)| is the order of the zeroth Brown cohomology group of the complex (C A , G). Therefore, GSD A is a topological invariant of A. However, there are more topological contributions to the entanglement entropy coming from the entangling surface ∂A, which can have non-trivial topology. They are hidden in the term |Im(δ −1 A )|/|Im(δ −1 A )|, which essentially counts the number of non-equivalent gauge transformations labelled by simplexes at the boundary ∂A of A. To see this, first we note that, in the same way we constructed a co-chain complex that encodes an Abelian higher gauge theory defined on the region A but with higher gauge transformations restricted to act on the interiorÃ of A, we can also construct a co-chain complex that represents an Abelian higher gauge theory defined on A but with higher gauge transformations that act only on the boundary ∂A of A.
Let K n,∂A = {x ∈ K n,A |x ∩Ã = ∅} be the set of n-simplices that have no intersection with the interior of A. Clearly, K n,Ã ∩ K n,∂A = ∅ and K n,Ã ∪ K n,∂A = K n,A . Higher gauge transformations over A are parameterized by elements of the group hom(C A , G) −1 . Therefore, the transformations that act over the boundary of A are parameterized by the subgroup hom(C ∂A , G) −1 ⊂ hom(C A , G) −1 , whose elements are collections of maps f = {f n } such that, for each 0 ≤ n ≤ d, f n : K n,∂A → G n+1 has support in the boundary of A. Then, we define the restriction δ −1 is a co-chain complex since δ 0 A • δ −1 ∂A = 0. This co-chain complex defines an Abelian higher gauge theory on A with higher gauge transformations labelled by elements of ∂A. From (4.15) and (4.25), we can establish the following result: Proof. See appendix C.
We can then focus in studying the set Im(δ −1 ∂A ). First, note that we can construct, for any −d ≤ p ≤ d, the group hom(C ∂A , G) p = d−1 n=0 Hom(C n,∂A , G n−p ) and the coboundary map δ p ∂A : hom(C ∂A , G) p → hom(C ∂A , G) p+1 , just as we did above. Therefore, we can construct the co-chain complex 27) and, for this co-chain complex, we can also define the Brown cohomology groups H p (C ∂A , G) = ker(δ p ∂A )/Im(δ p−1 ∂A ). Indeed, since δ p ∂A : hom(C ∂A , G) p → hom(C ∂A , G) p+1 is a group homomorphism, the first isomorphism theorem says that hom(C ∂A , G) p /ker(δ p ∂A ) ∼ = Im(δ p ∂A ), (4.30) which implies that Now, from equation (4.28), we have that so, substituting |ker(δ p ∂A )| = |H p (C ∂A , G)||Im(δ p−1 ∂A )| into (4.31), we get (4.29). This equation is a recursion relation that allows us to write |Im(δ p ∂A )| in terms of |Im(δ p−1 ∂A )|. In particular, for p = −1, Applying (4.29) once more, we have This procedure can be continued until we finally reach where α = 1 when p is odd and α = −1 when p is even. Thus, we can write |Im(δ −1 ∂A )| as a product of terms that depend on the geometry of the entangling surface ∂A, i.e., the |hom(C ∂A , G) p | terms, with terms that depend on the topology of ∂A, that is, the We can now go back to the entaglement entropy S A . We have that, from equations (4.23), (4.24) and lemma 4.1,

JHEP03(2020)167
Now, as explained in [19], here we are considering the case where there are finitely many n's such that K n is non-empty, and each of these non-empty sets are also finite. Moreover, all groups G n appearing in the chain complex (G, ∂ G ) are finite. Thus, the groups hom(C, G) p are finite for all p, and their order is given by In particular,   Therefore, the entanglement entropy S A can be written as is the "area law" term, that is, the term which explicitly depends only on the geometry of ∂A, and is the topological entanglement entropy, i.e., the term that explicitly depends on the topology of both A and ∂A. To calculate it, we must be able to compute cohomology groups with coefficients in the homology groups of the chain complex (G, ∂ G ). Cohomology with coefficients is related to the usual integral homology through the universal coefficient theorem (see [33] for a general reference), which states that, for any Abelian group D and any 0 ≤ n ≤ d, JHEP03(2020)167 where H n (C) is the homology group of order n with integer coefficientes of the chain complex C = (C, ∂ C ). The Ext 1 term is related to the torsion part of H n (C), and it is trivial whenever H n (C) is free. Writting H n (C) ∼ = Z βn ⊕ T i , where Z βn corresponds to the free part of H n (C), β n being the Betti number of order n of the chain complex C, and T i is its torsion part, we can see that indeed there is a relation between S Topo and the higher Betti numbers of both the region A and its boundary ∂A. In particular, when we consider manifolds with torsion-free homology groups, the Ext 1 term in (4.46) is trivial and we have as a result

Examples
In this section we calculate the entanglement entropy of the model shown in section 2. We will use this example to demonstrate how the topology of both A and ∂A affect the entanglement entropy. For 1-gauge theories (Abelian Quantum Double Models), we recover well-known general results for any dimension d.

1-gauge theories
We consider the general Abelian 1-gauge theory, also known as the Abelian Quantum Double model, but now in any dimension d. This theory is defined by the chain complexes in figure 5. This example can be seen as a particular case of the d-dimensional version of the model introduced in 2.1, where the only non-trivial group is G 1 .
The geometrical chain complex, generated by the simplicial complex K = d n=0 K n , corresponds to the d-dimensional manifold M in which the theory is defined. Let's consider a bipartition of this manifold into two regions, A and B, where A is a d-dimensional closed immersed submanifold of M , as in the general case discussed in section 4.1. This corresponds to splitting the simplicial complex K as K = d n=0 K n,A ∪ K n,B , where d n=0 K n,A JHEP03(2020)167 is a subcomplex of K. Therefore, following the general procedure shown in section 4.1, we can construct for region A the chain complexes shown in figure 6. From these chain complexes, as was demonstrated in sections 4.1, 4.2 and 4.4, the following cochain complexes are well defined for this theory: We can then apply equation (4.43) to find the entanglement entropy S A of this theory. First, let's compute S ∂A , given by equation (4.44). We have that This quantity depends only on the simplicial complex of the boundary of A and on the higher gauge groups of the theory in question. For 1-gauge, there is only one non-trivial gauge group G 1 , all other groups being equal to {0}, as is shown in figure 6. Therefore, since the order of trivial groups is equal to one, the only term which will survive in the double sum in equation ( Thus, we have indeed an "area law" term, because S ∂A is proportional to the number of vertices in the boundary of A, which is essentially the size of the boundary. We are only left to calculate S Topo . From equation (4.45), we know that We will leave the topology of A and ∂A unspecified for a moment and proceed as far as possible. For a 1-gauge theory, the non-trivial piece of the chain complex (G, ∂ G ) is given by and the only non-trivial homology group related to this chain complex is -25 -

JHEP03(2020)167
Therefore, in the first term of (5.5) the only non-zero contribution to the sum comes from n = 1, while for the second term of (5.5) the only non-zero contribution to the sum comes from n = 0, p = 1. Hence, the topological entanglement entropy in Abelian 1-gauge theories is given by Then, to compute S Topo we need to calculate the order of the cohomology groups H 1 (C A , G 1 ) and H 0 (C ∂A , G 1 ), with coefficients in G 1 . To do so, we employ the universal coefficient theorem for cohomology (4.46), which states that for every 0 ≤ n ≤ d, where H n (C) is the nth homology group with integer coefficients of the chain complex C. Here we only need the n = 0 and n = 1 cases. For n = 0, there is no Ext term [33], so We have also that H 0 (C ∂A ) ∼ = Z β 0 (∂A) , where β 0 (∂A) is the zeroth Betti number of ∂A. Therefore, H 0 (C ∂A , G 1 ) ∼ = Hom(Z β 0 (∂A) , G 1 ) (5.10) and hence Likewise, the Ext term is trivial for the n = 1 case, and we can write So in general, the topological entanglement entropy for Abelian 1-gauge theories reads and it depends on the number of non-contractible curves one can draw over region A, as well as on the number of connected components of ∂A. To proceed further, we need to specify the topology of A. For example, if A has the topology of a d-dimensional ball, its first homology group is trivial and hence we get that S Topo = −β 0 (∂A)log (|G 1 |), i.e., it is only sensitive to the number of connected components of the entangling surface ∂A.
It is important to mention that we could also build theories with degrees of freedom attached to other components of the lattice, such as its 2-simplices (plaquettes). Doing so, we have a particular case of a 2-gauge theory. As we will see in the following discussion, this makes the topological entanglement entropy depend on higher Betti numbers. An example of such model is the 4-dimensional Toric Code [36]. Note that in [36] it is not JHEP03(2020)167 Figure 7. Chain complexes that define Abelian 1-gauge theories in d dimensions with degrees of freedom living at the faces of the lattice. f 2 are the configurations, t 1 are the gauge transformations labeled by links, and there is also m 3 as the 2-holonomies.
called a 2-gauge theory. We do so in order to be consistent with our formalism. Let's then consider an Abelian 2-gauge theory in d-dimensions with degrees of freedom living at the 2-simplices of the lattice. The chain complexes describing this model is presented in figure 7. We can repeat exactly the same steps we did for calculating the entanglement entropy of the 1-gauge model obtained from the chain complexes in figure 5. Hence, we can skip this discussion and right away apply the formula (4.43) to calculate the entanglement entropy of this model. Here, since the only non-trivial group is G 2 , the only non-zero terms in the sum that define S ∂A are the p = 1, n = 1 and p = 2, n = 0 terms, so S ∂A is given by (5.14) Note that this term vanishes for a 3-dimensional region A with periodic boundary conditions. Now, since the only non-trivial piece of the group chain complex shown in figure 7 is the only non-trivial homology group associated to this complex is Thus, the non-zero contribution to the first term of S Topo comes from the term with n = 2, while the non-zero contributions to the second term of S Topo come from the terms with p = 1, n = 1 and p = 2, n = 0. Therefore, the topological entanglement entropy of this model is given by Again, using the universal coefficient theorem (4.46) to calculate the cohomologies, we first note that there is no Ext term for n = 0 and it is trivial for n = 1. Therefore, where we used again that H 0 (C ∂A ) ∼ = Z β0(∂A) . However, applying the theorem to |H 2 (C A , G 2 )|, we have that and the topological entanglement entropy now is related to the number of non-contractible curves and surfaces one can draw over ∂A and A, respectively, and therefore S Topo depends on higher Betti numbers of both A and ∂A.
To give an example where the Ext term is non-trivial, take the 4-dimensional manifold A = RP 3 × [0, 1], i.e., the product of the real projective space with the unit interval. It is a manifold whose boundary is ∂A = RP 3 ∪ RP 3 . The homology groups of RP 3 are otherwise.

(5.19)
Let's compute the second homology group of A. By Künneth's theorem [37], we have that that is, The first homology group of A = RP 3 × [0, 1] can be computed in the same way: that is, Then, from equation (5.18), we have and Z 2 ⊗ Z ∼ = Z 2 [37] and that because Ext 1 (Z, G 2 ) = {0} [33], we have where we used that Ext 1 (Z n , G) ∼ = G/nG, for any n ∈ N and any Abelian group G [33]. Now, the order of the cohomology groups (5.16) and (5.17) of the boundary ∂A = RP 3 ∪ 24) and thus, the topological entanglement entropy is given by i.e., S Topo = log(|G 2 /2G 2 |) + 2log(|G 2 |). (5.26) Expanding the results found in the literature [25,38], here we demonstrated that the topological entanglement entropy depends not only on the topology of the entangling surface ∂A, but also on the topological properties of the bulk region A. Moreover, the Betti numbers of A and ∂A are not the only information needed to obtain the topological entropy. For some models, it may depend also on torsion properties of the sub-region A and its boundary, captured mainly by the Ext 1 functor.

(2D) 0, 1-gauge theories
Let's consider a 0, 1-gauge theory, a particular case of the one exhibited in section 4 with G 2 = 0. For simplicity, we focus in the 2-dimensional case, but the discussion presented here can immediately be extended to any dimension d. The chain complexes defining this theory are shown in figure 8.
Again, we divide the lattice into two regions, A and B, where for A we have a subcomplex d n K n,A . Following the general procedure shown in section 4.1, we can construct for region A the chain complex (C A , ∂ C A ) and find that the entanglement entropy S A is given by equation (4.43). Let's calculate first S ∂A . It is straightforward to see that the JHEP03(2020)167 only non-zero contribution to the sum in (4.44) is given by the n = 0, p = 1 term. So we have that S ∂A = |K 0,∂A |log(|G 1 |), (5.27) and it is again an "area law". Now, consider the (G, ∂ G ) chain complex of this theory: The non-trivial homology groups associated to this chain complex are So, the non-zero contributions to the topological entanglement entropy (4.45) are This result is true for any Abelian finite groups G 0 and G 1 . However, now we have to specify these groups in order to calculate the homology groups (5.29) and (5.30). Thus, let's first consider the model with The topological entanglement entropy of this model is thus Note that this result is equal to the G 1 = Z 2 1-gauge theory. Again, this model exhibits the same topological properties as the Toric Code. We can also consider the model where G 0 = Z 2 and G 1 = Z 2 , with ∂ G 1 : Z 2 → Z 2 being the identity map. In this case, both H 0 (G) and H 1 (G) are equal to the trivial group, and the topological entanglement entropy is equal to zero, confirming the non-topological nature of this model.

(3D) 1, 2-gauge theories
This time we consider the case of the 1, 2-gauge theory (section 2.1 and section 4 with G 0 = 0). Although here we treat the 3-dimensional case, the procedure below can right away be extended to arbitrary dimensions d. The chain complexes that define this theory are reproduced here in figure 9.
Once more, we divide the lattice into two regions, A and B, where region A is such that we have a subcomplex K A . From the general procedure shown in section 4.1, we can construct for region A the chain complex (C A , ∂ C A ) and then find that the entanglement entropy S A is given by equation (4.43). To compute it for this model, let's calculate first JHEP03(2020)167 S ∂A . We see that the only non-zero contributions to the sum in (4.44) are the ones given by the n = 0, p = 1, n = 1, p = 1 and n = 0, p = 2 terms. Therefore Note that, if we were dealing with a two dimensional system with periodic boundary conditions, the terms proportional to log(|G 2 |) would cancel out, as the number of links and vertices is the same, and we would have the same result as the one found in the 1-gauge case. Now, to calculate the topological entanglement entropy S Topo , given by equation (4.45), first we consider the chain complex (G, ∂ G ) of this model: The non-trivial homology groups associated to this chain complex are as follows: Thus, S Topo is given by To proceed further, let's choose the model where G 2 = Z 4 , G 1 = Z 2 and ∂ G 2 : Z 4 → Z 2 such that ∂ G 2 (1) = 1. We can then calculate the homology groups (5.38) and (5.39) to be Therefore, the topological entanglement entropy of this model is Figure 10. Chain complexes for the 1, 2, 3-gauge model.
This formula is equal to the 2-gauge (5.15) case, with degrees of freedom living in the plaquettes of the lattice and G 2 = Z 2 . Note that we never used the fact that we are dealing with a 3-dimensional system to derive equation (5.44), which means that the same result holds for dimension d ≥ 3. Then, if we for example take d = 4, this choice of groups for the 1, 2-gauge model imposes that its long-range entanglement characteristics, detected by the topological entanglement entropy, are the same as the 4D Toric Code one studied in [36]. Other choices of groups may generate 1, 2-gauge models with more unusual behaviors.

(4D) 1, 2, 3-gauge theories
Now we consider the 1, 2, 3-gauge theory in four dimensions. We do this to show that our formalism allows us to readily shift from a 3-dimensional presented in 2.1 case to a 4D, our formalism can be extended to any arbitrary dimension d. The chain complexes that define this theory is shown in figure 10. We again divide the lattice into two regions, A and B, where region A is such that we have a subcomplex d n=0 K n,A . The general procedure shown in section 4.1 allows us to construct for region A the chain complex (C A , ∂ C A ) and thus we can find the entanglement entropy S A using equation (4.43). To calculate it for the 1, 2, 3-gauge model, let's study first the term S ∂A . We see that the only non-zero contributions to the sum in (4.44) are the ones given by the p = 1, n = 0, 1, 2, p = 2, n = 0, 1 and p = 3, n = 0 terms. Therefore that is, To calculate the topological entanglement entropy (4.45), we first consider the chain complex (G, ∂ G ) of the 1, 2, 3-gauge model: The non-trivial homology groups associated to it are the following ones:

JHEP03(2020)167
Hence, S Topo is given by Note that the topological entanglement entropy does not change from the 3D to the 4D case, because the homological groups are the same. To give a more concrete example, let's consider the case where G 1 = G 2 = G 3 = Z 4 and the homomorphisms ∂ G 3 (1) = ∂ G 2 (1) = 1. The lattice is a discretization of a solid ball S 3 . In this case, we can calculate the homology groups (5.47), (5.48) and (5.49). They are Therefore, the topological entanglement entropy of this model is Let's use the universal coefficient theorem to compute these cohomology groups. We have that, as before, |H 0 (C ∂A , Z 2 )| = |Hom(H 0 (C ∂A ), Z 2 )| = 2 β 0 (∂A) , (5.55) Since in this case ∂A = S 2 , we have that β 0 (∂A) = 1, H 1 (C ∂A ) = H 1 (C A ) = {0}. Therefore, Hence, the topological entanglement entropy of this model is given by We see that, although we defined the model over a manifold with trivial topology and the ground state degeneracy of this model does not exhibit a topological dependency, the topological entanglement entropy is different from zero, indicating the presence of longrange entanglement. JHEP03(2020)167

Conclusions
The paper carried out the calculation of the entanglement entropy for all Abelian higher gauge theories in a comprehensive way. Furthermore we could separate the entropy into the topological information and the geometrical one. We started by making a review of the kind of models we treated. Then we described them in very general terms, as introduced in [19]. The calculation followed from the definition of the density matrix ρ as being proportional to the ground state projector, see (4.9). To obtain the reduced density matrix we considered a bipartition of the simplicial complex K into a subcomplex K A and its complement. The partial trace over the unknown region was used to obtain the reduced density matrix ρ A , which included operators that were exclusively supported in K A , see (4.19). From the Von Neumann entropy formula we derived the entanglement entropy and showed that it could be naturally interpreted as the ground state degeneracy of the same model but restricted to the subcomplex K A , see (4.23). In this sense, we mapped the problem of calculating the entanglement entropy of a higher gauge theory to a problem of counting the flat edge states of the theory restricted to region A. Then, we further divided this restricted ground state degeneracy into two contributions, one comming from the bulk region A and the other comming from its boundary ∂A, and we showed that this splitting allows us to write the entanglement entropy as a sum of two terms (4.43): one being the area law, i.e., a term depending only on the geometry of the entangling surface ∂A, and the other being the topological entanglement entropy, a term depending on the topological properties of both A and ∂A.
We demonstrated a formula for the topological entanglement entropy S Topo in terms of the cohomology groups with coefficients in the homology groups of the complex (4.3). The universal coefficient theorem can be applied to give a formula for S Topo in terms of the integral homology groups of the manifold in question, which in turn can be used to express the topological entropy in terms of the Betti numbers of the underlying space and its boundary. However, our equations show that, even in regular 1-gauge theories (Abelian Quantum Double models), S Topo can depend on torsion properties of the manifold.

A Trace of local operators
In this appendix we show how taking the partial trace of the ground state projector, or any product of projection operators of the theory, implies in eq. (4.16).
We begin by writing the density matrix, ρ, using the local decomposition of A 0 and B 0 (see [19] for a detailed account on this). The local decomposition yields  where π : B → B Imφ is the canonical projection sending b ∈ B into its corresponding equivalence class [b] ∈ B Imφ Furthermore, β ∈ Hom B Imφ , U(1) is unique and defined as: notice that β is well defined within equivalence classes since Imφ ⊂ kerβ. To see this, consider b = b ∈ [b], this means that b − b ∈ Imφ ⊂ kerβ, therefore: This is, we have shown that given an irrep β ∈ kerφ then there is a unique morphism β ∈ Hom B Imφ , U(1) . We now need to show that the converse also holds, to this intent, consider β : B Imφ → U(1). Recall that Imφ ⊂ kerβ. Observe also that β is the only map for which the diagram in B.1 commutes.

JHEP03(2020)167
Proof. Consider ker(δ −1 A ), a subgroup of hom(C A , G) −1 . From the definitions of hom(CÃ, G) −1 , δ −1 A , hom(C ∂A , G) −1 and δ −1 ∂A , it is clear that and ker(δ −1 ∂A ) = {f ∈ hom(C ∂A , G) −1 |δ −1 ∂A f = 0} are subgroups of ker(δ −1 A ). We introduce the following equivalence relation on ker(δ −1 That is, two collections of maps f and f in ker(δ −1 A ) are equivalent if they differ by a collection of maps in ker(δ −1 A ), i.e., maps with support in KÃ = ∪ d n=0 K n,Ã that are also sent to the trivial map by the co-boundary operator. This equivalence relation defines the quotient group ker(δ A ) such that f − f = g, so f | ∂A − f | ∂A = g| ∂A = 0, because g ∈ hom(CÃ, G) −1 . Therefore, f | ∂A = f | ∂A . Moreover, φ is an isomorphism. To see this, first take [f ] ∈ ker(φ). So, φ([f ]) = 0 ⇔ f | ∂A = 0, which means that f is a collection of trivial maps. Therefore, ker(φ) = {0} and φ is injective. The map φ is also surjective, because if we take a map g ∈ ker(δ −1 ∂A ), it is a collection of maps which are zero everywhere except in the boundary of A and it can be obtained by applying φ in the class [g] ∈ ker(δ −1 A )/ker(δ −1 A ). Therefore, φ is bijective. Now, let Now, since K n,Ã ∩ K n,∂A = ∅ and K n,Ã ∪ K n,∂A = K n,A for any n = 0, . . . , d, we can write the n-chain group C n,A as a direct sum of subgroups C n,A = C n,Ã ⊕ C n,∂A , i.e., every JHEP03(2020)167 c = x∈K n,A c(x)x ∈ C n,A can be written as c = x∈K n,Ã c(x)x + x∈K n,∂A c(x)x. This implies that any homomorphism f n : C n,A → G n+1 , where G n+1 is some arbitrary finite Abelian group, can be written as f n : C n,Ã ⊕ C n,∂A → G n+1 . Thus, Hom(C n,A , G n+1 ) = Hom(C n,Ã ⊕ C n,∂A , G n+1 ).