{'fraction_non_alphanumeric': 0.024967560069801782, 'fraction_numerical': 0.00637612421137411, 'mean_word_length': 4.416091118381194, 'pattern_counts': {'":': 0, '<': 10, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Area law violations in a supersymmetric model

We study the structure of entanglement in a supersymmetric lattice model of fermions on certain types of decorated graphs with quenched disorder. In particular, we construct models with controllable ground state degeneracy protected by supersymmetry and the choice of Hilbert space. We show that in certain special limits these degenerate ground states are associated with local impurities and that there exists a basis of the ground state manifold in which every basis element satisfies a boundary law for entanglement entropy. On the other hand, by considering incoherent mixtures or coherent superpositions of these localized ground states, we can find regions that violate the boundary law for entanglement entropy over a wide range of length scales. More generally, we discuss various desiderata for constructing violations of the boundary law for entanglement entropy and discuss possible relations of our work to recent holographic studies.


I. INTRODUCTION
Entanglement plays a fundamental role in quantum many-body physics. It can be used to identify topological phases of matter not characterized by any broken symmetry [1,2] as well as gapless critical points [3][4][5][6][7]. The structure of many-body entanglement has led to a wide variety of promising new variational states [8][9][10] and can be understood as underlying the functioning of DMRG. Entanglement also provides a new perspective on global issues about the landscape of quantum phases, including questions about the structure of renormalization group flows [11]. Despite these many promising achievements, we are still exploring the basic structure of entanglement and do not yet have a general picture of many-body entanglement. Even a basic question like the scaling of entanglement with system size remains open despite many recent advances. In these notes we describe a large class of models built using the technology of supersymmetry where we can construct states with tunable entanglement properties.
The entanglement entropy, defined as the von Neumann entropy of the reduced density matrix of a spatial subsystem, typically obeys an area law for many-body ground states of local Hamiltonians [12]. That is, S L = −Tr(ρ L log ρ L ) ∼ L d−1 , where L is the linear size of the spatial subregion and d the spatial dimension of the system. There are, however, examples of gapless states that violate the area law logarithmically, i.e., S L ∼ L d−1 log L. The most well-known example being the Landau Fermi liquid ground state [13][14][15][16]. There has also been considerable recent interest in violations of the area law from the point of view of holographic duality [17][18][19]. Under certain constraints, the most important being the absence of long-range interactions, ground state degeneracy, and infinite fine-tuning, it seems that the area law is violated at most logarithmically at the largest scales [20].
When these constraints are relaxed, however, there are constructions that lead to ground states with an anomalous violation of the area law. For instance, when long-range interactions (specifically long range hoppings) are allowed it is possible to construct one dimensional systems for which the entanglement entropy of a pure ground state scales as any function of L less than extensive [21]. Another possibility is to have a large ground state degeneracy, for instance models with an extensive ground state entropy (see e.g. [22][23][24]) are likely to also have an extensive entanglement entropy. Finally, models with entanglement entropy scaling in the whole range between area law and extensive were recently constructed in the context of holography [19]. For the latter models it is hard to say which of the above constraints are violated. Furthermore, it is not well understood in what state the entanglement entropy is computed and if the anomalous scaling persists to all scales in the infrared. However, due to the generality of the construction, the answer to these and related questions may well provide further guidelines as to what type of entanglement properties may be expected.
The way anomalous violations of the area law may arise in systems with a large ground state degeneracy can be understood as follows. Let {|α } be the basis spanning the ground state subspace of the Hilbert space. It follows that a typical ground state is a cohorent state described by the density matrix ρ c = α,β c α c β |α β|. For sufficiently orthogonal states, |α , the reduced density matrix, obtained by tracing out a spatial subregion, will be very similar to the density matrix of a completely mixed state. For systems with a large ground state degeneracy we thus conclude that the entanglement entropy can grow as ∼ log N gs (L) + S α , where N gs (L) is the number of ground states in a system of linear size L and S α is the entanglement entropy of the typical ground state |α . Typically S α obeys the area law. For a system with a ground state degeneracy that is exponential in the volume, the entanglement entropy of the coherent state is thus completely dominated by the first part: log N gs ∼ L d . From this perspective, for a system to have an entanglement entropy that scales ∼ L d−a , with 0 < a < 1 , we need the ground state degeneracy to be exponential in L d−a . To our knowledge all known examples of local Hamiltonians with a large ground state degeneracy have either a = 0 or a = 1 (see e.g. [25]).
A striking example of a model with a large ground state degeneracy is the supersymmetric model for strongly interacting itinerant fermions [26] (for a review see [27]). The Hamiltonian is supersymmetric by construction: it is defined as the anti-commutator of two supercharges. This construction allows for an increased analytic control of a lattice model for strongly interacting fermions. In particular, it is often possible to compute the ground state degeneracy exactly by solving the cohomology problem. From this model a large variety of ground state degeneracies is observed, ranging from a unique ground state to a ground state degeneracy exponential in the number of sites of the system. It is thus natural to ask, as suggested in [18], if there are supersymmetric lattice models with a ground state degeneracy that allows for an anomalous violation of the area law.
Here we construct a class of graphs on which the supersymmetric lattice model has a ground state degeneracy that can be controlled at will. In particular, if the graph has volume L d , with d the dimension of the graph, the ground state degeneracy grows as ∼ 2 t , where t can be any number between 0 and L d .
This work is organized as follows. First, we present a set of criteria that, for our taste, should be obeyed when searching for violations of the boundary law. Second, we present our supersymmetric lattice model and analyze in detail its entanglement properties. Finally, we discuss in the larger context the virtues and defects of our construction.

II. MODEL CRITERIA
It is widely believed that gapped phases in higher dimensions satisfy a boundary law for entanglement entropy, but even this "obvious" statement is unproven outside the one dimensional arena. In searching for violations of the boundary law, we want to set out certain basic rules. These rules represent an attempt to select microscopic models that are realistic and generic for quantum many-body physics.
• [Short-range interactions] We only consider systems without long-range interactions.
Of course, gapless modes are allowed, but we do not "integrate them out" to produce long-range interactions for the remaining degrees of freedom. If this requirement is relaxed, it is known that one can obtain a plethora of entanglement scalings. One cause for concern is the Coulomb interaction which, in insulating phases, is long ranged (although this is quite different from having long range hoppings). In appendix A we give an alternative fermion construction of a state with anomalous entanglement entropy using long-range hoppings.
• [Dimensionality] We only consider systems that can be ascribed an integer spatial dimensionality. We do not consider fractal graphs or other spaces with unusual spatial structure e.g. hyperbolic space. It is not that these spaces are not interesting, only that the notion of boundary law violation becomes less meaningful.
• [Randomness] Our construction will make use of non-translation invariant graphs or quenched randomness which is of a rather peculiar type. While we would like to have models that do not rely on this physics, realistic materials do have impurities, so although the type of randomness we consider is unusual, quenched randomness is present in general.
• [Fine-tuning] We would like to avoid as much as possible fine-tuning of the Hamiltonian. Any finite amount of fine-tuning can be interpreted as some kind of multicritical phenomenon, but with an infinite amount of fine tuning we can also easily construct models with anomalous entanglement e.g. spin chains where every spin in a segment forms a singlet with a spin outside the segment.
Our construction is a supersymmetric lattice model in any dimension with short range interactions and dilute randomness. Except for a vanishing density of impurities, it behaves like a conventional quantum system in integer dimension. Furthermore, although the supersymmetry implies some degree of fine tuning, it is not infinitely fine tuned.

III. THE MODEL
Let us now define the model. The degrees of freedom are spinless fermions living on the square lattice. A fermion at site i is created by the operator c † i with {c i , c † j } = δ ij . The fermions have a hard core, meaning that they are not only forbidden to be on the same site as required by Fermi statistics, but are also forbidden to be on adjacent sites. Their creation operator is is zero if any site next to i is occupied. The Hamiltonian is defined in terms of the supercharges, Q = i d † i and From its definition it follows that the Hamiltonian commutes with the supercharges. Imposing this structure has some immediate consequences: supersymmetric theories are characterized by a positive definite energy spectrum and a twofold degeneracy of each non-zero energy level. The two states with the same energy are called superpartners and are related by the supercharge. In other words, all eigenstates with an energy E s > 0 form doublet representations of the supersymmetry algebra. A doublet consists of two states (|s , Q|s ), such that Q † |s = 0. Finally, all states with zero energy must be singlets: Q|g = Q † |g = 0 and conversely, all singlets must be zero energy states [28].

A. Cohomology
For the supersymmetric models, cohomology has proven to be a very powerful tool to extract information about the zero energy ground state(s) of the models (see for example [23,26,27,[29][30][31]). The key ingredient is the fact that ground states are singlets; they are annihilated both by Q and Q † . This means that a ground state |g is in the kernel of Q: Q|g = 0 and not in the image of Q, because if we could write |g = Q|f , then (|f , |g ), would be a doublet. So the ground states span a subspace H Q of the Hilbert space H of states, such that H Q = ker Q/ImQ. This is precisely the definition of the cohomology of Q. So the ground states of a supersymmetric theory are in one-to-one correspondence with the cohomology of Q. It follows that the solution of the cohomology problem gives the number of zero energy states for each particle number sector. Equivalently, we find that zero energy states are in one-to-one correspondence with the homology elements of Q † .
We compute the cohomology using the 'tic-tac-toe' lemma of [32]. This says that under certain conditions, the cohomology H Q for Q = Q 1 + Q 2 is the same as the cohomology of Q 1 acting on the cohomology of Q 2 . In an equation, H Q = H Q 1 (H Q 2 ) ≡ H 12 , where Q 1 and Q 2 act on different sublattices S 1 and S 2 . We find H 12 by first fixing the configuration on all sites of the sublattice S 1 , and computing the cohomology H Q 2 . Then one computes the cohomology of Q 1 , acting not on the full space of states, but only on the classes in H Q 2 . A sufficient condition for the lemma to hold is that all non-trivial elements of H 12 have the same f 2 (the fermion-number on S 2 ).
Although ground states are in one-to-one correspondence with cohomology elements, the two are not equal unless the cohomology element happens to be a harmonic representative of the cohomology. Harmonic representatives are elements of both the cohomology of Q and the homology of Q † . So they are annihilated by both supercharges, which is precisely the property of a zero energy state. It follows that, although the solution of the cohomology problem gives the number of ground states, it typically does not give the ground states themselves. For a more leisurely introduction to cohomology and an exposition of the relation between the supersymetric model and independence complexes we refer the reader to [27,31].
We now briefly state the cohomology results for the one dimensional chain [26], because we will use these results many times throughout the paper. For the periodic chain of length L = 3n + a the cohomology is trivial for all particle numbers, f , except for f = n, where we have  for a chain of length L = 3n + a and n particles and dim(H Q ) = 0 at all other particle numbers. In particular, we note that the cohomology of an isolated site that can be both empty and occupied (L = 1 with open boundary conditions) is trivial. This is equivalent to the statement that the single site chain has no zero energy states, indeed the empty and occupied state form a doublet of energy E = 1.

B. Decorated graphs with quenched disorder
The graphs we consider in this work are constructed as follows. We start with an original graph, Λ, with N Λ vertices and L Λ links. The only restriction on the graph, Λ, is that it does not contain disconnected subgraphs. Now construct the graph, Λ 3 , from the original graph by adding 2 additional vertices on every link. Then from Λ 3 the decorated graph, Λ 3,t , is constructed by adding t additional sites to the graph and connecting each added site to a unique pair of sites on a link of the original graph, Λ. If we call the collection of sites on the original graph, Λ, the subgraph S 1 and all the added sites subgraph S 2 , we find that S 2 is a collection of L Λ − t open 2-site chains and t periodic 3-site chains. An example is shown in figure 1.

IV. GROUND STATE DEGENERACY
In this section we compute the ground state degeneracy of the supersymmetric model on the decorated graphs with quenched disorder, Λ 3,t . The supersymmetric model on the decorated graphs without disorder, Λ 3 , is discussed in detail in [33]. The cohomology problem is easily solved and the ground state is found to be twofold degenerate [23]. As a warm-up we briefly review this computation. Consider the original graph, Λ, as the subgraph S 1 and the additional sites as the subgraph S 2 . The subgraph S 2 is a collection of two site chains. Remember that an isolated site that can be both empty and occupied has a trivial cohomology. It follows that the non-trivial elements in H Q 2 have the S 1 sites neighboring a two site chain on S 2 either both empty or both occupied. Consequently, there are only two non-trivial elements in H Q 2 , one with S 1 completely empty and one with S 1 completely filled. If we leave S 1 completely empty, we obtain one non-trivial element of H Q 2 in the sector with L Λ fermions, where L Λ denotes the number of links in the original graph Λ. The element with S 1 completely filled, clearly has fermion number N Λ , where N Λ denotes the number of vertices in the original graph Λ. It now quickly follows that these two elements are also in H 12 , since within H Q 2 both states are in the kernel of Q 1 and not in the image of Q 1 . Now remember that for the 'tic-tac-toe' lemma to hold a sufficient condition was that all elements in H 12 have the same number of fermions on sublattice S 2 . Here this condition is clearly not met. One can easily show, however, that the lemma also holds when the elements of H 12 do not differ in their total fermion number by one. We thus conclude that H 12 = H Q provided that N Λ = L Λ ± 1.
As an example consider the square lattice with doubly periodic boundary conditions as the original graph Λ. We find that L Λ = 2N Λ and the total number of sites in Λ 3 is N = 2L Λ + N Λ . Consequently, this lattice has one ground state at 1/5 filling and one at 2/5 filling.
Let us now consider the cohomology problem on the graphs Λ 3,t . As before we define the subgraph S 1 as the original graph Λ and the subgraph S 2 as the rest of the sites. It follows that S 2 is a collection of L Λ − t open 2-site chains and t periodic 3-site chains. Note that a periodic 3-site chain has two zero energy ground states with fermion number one. It follows that the configuration with all S 1 sites empty corresponds to 2 t non-trivial elements in H Q 2 . Furthermore, we note that if two S 1 sites neighboring a 3-site chain are occupied, the 3-site chain effectively becomes an isolated site. Therefore the configuration with all S 1 sites occupied is not an element of H Q 2 . However, if only one of the two S 1 sites neighboring a 3-site chain is occupied, the 3-site chain effectively becomes a 2-site chain, which does have a zero energy state. It follows that other configurations on S 1 may correspond to non-trivial elements in H Q 2 . This complicates the computation of H 12 , however, it is easy to verify that the configuration with all S 1 sites empty is the only allowed configuration in H Q 2 under the following condition. The 3-site chains are distributed over the graph in such a way that any two S 1 sites neighboring a 3-site chain can be connected to each other by drawing a path that contains only 2-site chains on S 2 and S 1 sites. Although it may be possible to solve the cohomology problem in other cases as well, we do not consider these here, since the above constraint is acceptable for the purposes of the present work. Under the above condition, it immediately follows that all 2 t non-trivial elements of H Q 2 are also elements of H Q , since they are all in the kernel of Q 1 within H Q 2 and they all have the same number of fermions on S 2 .
We conclude that the ground state degeneracy of the supersymmetric model on the graph Λ 3,t is 2 t (t > 0) under the following conditions: • The original graph, Λ, does not contain disconnected subgraphs, • The original graph, Λ, has N Λ vertices and L Λ links such that N Λ = L Λ ± 1, • The t periodic 3-site chains are distributed over the graph in such a way that any two S 1 sites neighboring a 3-site chain can be connected to each other by drawing a path that contains only 2-site chains on S 2 and S 1 sites.
It is easily verified that these conditions can be met for t in the entire range from order one to order L d .

V. GROUND STATES FOR LARGE AND SMALL STAGGERING
For the supersymmetric model on Λ 3 graphs in the presence of a staggering parameter, y, it was found that in the limit of small and large staggering the ground states take on a simple form [33]. In particular, the ground state with L Λ fermions is a valence bond solid state in the limit of large staggering and a quantum liquid state in the limit of small staggering. The quantum liquid state can be expressed as a projected valence bond solid state. In this section we will briefly review there results and show that in the presence of quenched disorder similar expressions for the 2 t ground states can be found in the small and large staggering limits.
We introduce the staggering parameter, y, by redefining the supercharge as and similar for Q † , where Q i acts solely on the subgraph S i . The Hamiltonian then reads it is the hopping between staggered and non-staggered sites, and finally, H 2 consist of a diagonal and an off-diagonal part, where the off-diagonal part represents hopping between the added S 2 sites on a link of the graph Λ.
The limit of y → ∞ is easy. To zeroth order in 1/y the Hamiltonian reads H/y 2 ≈ H 2 = {Q † 2 , Q 2 }. The two elements of H Q 2 , one with S 1 completely empty and one with S 1 completely filled, are precisely the ground states of this Hamiltonian. We focus on the state with L Λ particles, that is the state with S 1 completely empty. The S 2 subgraph is a collection of two-site chains. The ground state of a two-site chain contains one particle that resonates between the two sites: 2|∅ . It follows that in the state with S 1 completely empty there is one particle on each two-site chain of S 2 sites that resonates between the two sites. We denote this product state by |VBS inspired by the obvious similarity to a valence bond solid state.
It is clear that this result carries over directly to the decorated graphs with quenched disorder. On the triangles there are now two single particle states that are annihilated by H 2 . In the limit of y → ∞ there is a basis of the 2 t dimensional ground state subspace in terms of valence bond solid states characterized by the internal states of the triangles. Let this basis be |{α } where α labels the state of the triangle on link .
The limit of small staggering is more involved. For y = 0 we have Q = Q 1 and we should solve the cohomology problem by looking at H Q 1 directly. Since S 1 is a collection of disconnected sites, the non-trivial cohomology elements of H Q 1 are all configurations which have at least one neighboring site of each S 1 site occupied. There are many such configurations, with particle numbers in the whole range between N Λ and L Λ . So for y strictly zero, we thus find that there is a large ground state degeneracy. It can be shown, however, that this degeneracy is lifted at second order in degenerate perturbation theory [33]. The effective Hamiltonian at order y 2 that couples the zeroth-order ground states in the sector with L Λ particles can be written as where the sums are over the links and vertices of the original graph Λ, τ runs over the possible single particle configurations on the S 2 sites added to a single link of the original graph Λ.
By the state |1 we denote the vertex configurations in which precisely one neighboring site is occupied. Finally,P projects onto the zeroth-order ground state subspace, that is, it projects out all configurations that contain vertices with all neighboring sites empty. For the decorated graph, Λ 3 , this effective Hamiltonian has precisely one zero energy ground state given byP |VBS 1 . It is instructive for the following to see why the projected VBS state is annihilated by the effective Hamiltonian. First consider the action of the first term and focus on a single link. If both states τ on that link are allowed then the singlet state is annihilated by τ,τ |τ τ |. However, if only one state is allowed, then the action of τ,τ P |τ τ |P is purely diagonal.
It is not hard to see that this happens when the state τ is such that the particle is the only particle neighboring the adjacent S 1 site, that is the neighboring vertex is in state |1 . It follows that the action of the first term in the effective Hamiltonian on the projected VBS state can be written as vertices |1 1|, which precisely cancels the second term.
For the decorated graphs with quenched disorder in the limit of small staggering the above arguments basically go through unaltered, with the exception that the effective Hamiltonian now has 2 t zero energy states that can be written in a maximally localized form asP |α , where |α ∈ |{α } is a VBS state characterized by the states of the triangles. Indeed, one readily verifies that the effective Hamiltonian annihilates all these projected VBS states, the only additional check we need to do is to make sure that these states span the full 2 t dimensional subspace. We check this by verifying thatP has full rank on the subspace spanned by the VBS basis |{α } . It is clear that a linear superposition of the VBS states is annihilated byP if and only if it consists exclusively of configurations that contain vertices with all neighboring sites empty. Since the VBS states differ only in the states on the triangles, we can choose a superposition of VBS states such that on some link the particle sits on the S 2 site that has no neighboring S 1 sites. We could therefore obtain a superposition of VBS states that consists exclusively of configurations where one vertex has all neighboring sites empty, if all links coming out of that vertex contain triangles. This last condition, however, is in contradiction with the condition we imposed on the distribution of triangles over the graph Λ. We thus conclude thatP has full rank on the subspace spanned by the VBS states. It follows that the projected VBS states span the full 2 t dimensional ground state subspace. We conclude that in the limit of small and large staggering there is a maximally localized basis given byP |{α } and |{α } , respectively, that spans the 2 t dimensional ground state subspace of the decorated graphs with quenched disorder, Λ 3,t . For intermediate values of the staggering parameter a simple physical picture of the ground states is lacking.
Finally, before we go on the consider the entanglement entropy of a typical ground state of the model with quenched disorder, we briefly discuss the entanglement properties of the (projected) VBS states. It is clear that each unprojected VBS state satisfies a boundary law for entanglement entropy. The entropy is simply given by the number of valence bonds cut by the entangling surface times the elementary entropy of a single valence bond. The only issue is the projection. We have already shown that the projectorP has no kernel acting on the VBS states. The other important ingredient is to determine the overlap between two different projected VBS states. Given the gapped and short ranged nature of the VBS states, this overlap should be exponentially decaying in distance. In other words, two projected VBS states distinguished by configurations of spatially separated triangles will be orthogonal up to corrections exponential in the distance between the triangles. With this additional fact, we claim the projected VBS states also satisfy a boundary law. This is because after we orthonormalize the projected states we must only mix in exponentially small amounts of states associated with distant triangles, so we do not expect any qualitative change in the degree of entanglement within a given localized state.

VI. SYNTHESIS
We have shown that the supersymmetric model on the graph Λ 3,t has 2 t exact ground states. Once we include the deformation parameter y, we can provide model wavefunctions for these states at small and large y. For large y the ground states are simply VBS states, while for small y the wavefunctions are projected versions of the same VBS states. In either limit there is a basis of the 2 t dimensional ground state subspace in which the basis vectors are maximally localized e.g. associated with local changes in the states of triangles. Let this basis be |{α } where α labels the state of the triangle on link . As demonstrated above, each of these states satisfies a boundary law for entanglement entropy. Except for the complication of the projection, entanglement is due only to the short range valence bonds cut by the entangling surface. Let us now ask what sort of combinations of these local states should be allowed. Let A be a spatial region in the graph not necessarily consisting of sites from only S 1 or S 2 and let B be its complement. From the perspective of cooling to the ground state, a natural choice is the incoherent state This state is already mixed at the level of the total system and hence much of the entropy of ρ A is due to classical effects and not entanglement. Nevertheless, the entropy of A in this state goes like t log 2 + |∂A| where t = t A + t B . We can emulate this effect in a global pure state by entangling the triangles deep inside A with those deep in B. Suppose we begin with a particular product state |α . We allow ourselves to act on this state with single and two triangle unitary transformations, say between any triangle and its neighbors. On general grounds, a circuit of such one and two triangle unitaries consisting of a polynomial in t number of unitaries can produce a state that looks effectively mixed for any A with t A < t/2. This is simply the statement that thermalization can take place in a time polynomial in the system size. However, note that this "thermalization" process has nothing to do with the physical dynamics, it is simply a device to expose what sorts of superpositions in the ground state subspace are reasonable.
We should also point out that states selected from the ground state subspace with respect to the Haar measure on the projective space also have this property of being effectively mixed for t A < t/2 [34]. However, these states are not reachable from the product states |α with a polynomial in t number of one and two triangle unitaries. Hence we count them as unphysical since they could not be prepared in a reasonable time for large t.
The entanglement properties of A thus depend sensitively on t A and the choice of superposition in the ground state subspace. Note that while we assumed that two triangle unitaries between neighboring triangles should be allowed, in fact such neighboring triangles may be quite distant from each other if the triangle density is going to zero. Nevertheless, to obtain boundary law violations we must consider superpositions obtained from such unitaries. Thus fix a graph Λ 3,t with a given distribution of triangles and consider a family of regions A L labeled by linear size L such that t L ∼ L a with d − 1 < a < d. We assume that the total number of triangles and sites is much larger than all region sizes considered. For a state generated by a circuit of fixed depth consisting of one and two triangle unitaries described above, there will be a range of region sizes L min < L < L max for which the leading term in the entanglement entropy of A L is L a . This entropy violates the boundary law and is due to entanglement by construction (the global state is pure). L min is associated with lattice scale effects while L max reflects the fact that a finite depth quantum circuit can only entangle triangles a finite distance away. Once L max is reached, the entanglement of A L due to correlations between triangles will be bounded by L d−1 L max reflecting thermalization out to scale L max . The actual answer may be much less than this since the number of triangles in A is chosen to grow as a subextensive power of L.
Although we focused on a particular realization of the triangle disorder, it is interesting to consider specifying a probability distribution for the disorder. Suppose we introduce a probability measure p(x 1 , ..., x t ) for the positions of the triangles in the lattice. Let us assume this probability distribution is translation invariant so that p(x 1 + a, ..., x t + a) = p(x 1 , ..., x t ) for all a. For concreteness, we imagine a large square lattice base graph with periodic boundary conditions in which triangles are placed. If we further assume, as discussed above, that the states of interest look nearly maximally mixed for subsystems of size less than half the total system size, then the entanglement entropy of a subregion is simply proportional to the number of triangles contained within the subregion. Thus to compute the average entanglement entropy of region R we must average the counting function C R = i χ R (x i ) where χ R (x) is one for x ∈ R and zero otherwise.
The counting function has the property that C AB = C A + C B , where AB stands for the union of A and B, provided A and B are disjoint. Furthermore, the averageC R is clearly translation invariantC R+a =C R because of the translation invariance of p. These two properties are enough to prove thatC R ∝ |R| i.e. the average entanglement entropy satisfies a volume law. This is not terribly surprising given our assumption about the relevant superpositions in the ground state manifold. For example, if the probability distribution p is constant so that p = (1/L d ) t then we immediately findC R = t L d |R|. We thus conclude that to have any chance of obtaining something other than a volume law while maintaining translation invariance for the average entropy we would have to consider some other kind of superpositions in the ground state manifold. Note, furthermore, that our argument also suggests problems quite generally with having an anomalous scaling of the "local ground state degeneracy" with region size in a translation invariant system.

A. Scoring the model
We now briefly review our criteria and evaluate the successes and failures of our model. We also answer a few other potential objections to our construction.
• [Short-range interactions] Our Hamiltonian definitely consists of only short-ranged interactions. However, we did have to consider superpositions of ground states that can only be obtained from unitaries coupling triangles that may be dilute.
• [Dimensionality] When placed on a d-dimensional graph like a hypercubic lattice, our model preserves the dimensionality of the original graph. This is evidenced, for example, by the properties of diffusion on the graph. It is true that a small number of links are "anomalous", but there are not enough of these links to alter the dimensionality of the diffusion problem.
• [Randomness] Quenched randomness played a critical role in our construction. The entanglement entropy depends sensitively on which region we consider, but there are families of regions which show anomalous entanglement entropy for a range of sizes. If we perform a disorder average to restore translational invariance, we recover a volume law for the entanglement entropy.
• [Fine-tuning] Our model is expected to show fine tuning because of the essential role of supersymmetry, however, from a field theoretical point of view we would not expect to have to tune more than a finite number of supersymmetry breaking perturbations to zero in an RG sense. Below we discuss the role of supersymmetry preserving and breaking perturbations on the level of the microscopic model in some detail. There is considerable "tuning" in the sense that we must select special subregions to see the anomalous entanglement properties. It is not completely clear how to quantify the extent to which our construction is generic.
Another possible objection to our model comes from the small and large staggering limits, where it behaves somewhat like a self-contained system coupled to t free spins. However, we emphasize that this comparison is only approximately valid, and in particular, at finite y there are interactions between the spins and the rest of the lattice. The power of supersymmetry is that the exact ground state degeneracy remains even at y of order one where there is no notion of free spins. Thus the model we construct is much less trivial than free spins sprinkled through an otherwise conventional boundary law respecting system. Indeed, more generally, the ground state degeneracy is robust against various supersymmetry preserving perturbations. First of all, the cohomology argument goes through unaltered if we modify the supercharges by including any non-zero site-dependent coefficients in their definition: Q = i λ i d † i and Q † its hermitian conjugate, where λ i are real and positive (phases can be absorbed in the definition of d † i ). So the ground state degeneracy is unchanged for this class of perturbations. Note that the staggering considered above, with y > 0, falls into this class. It may be possible to construct other perturbations of the supercharges that leave the cohomology unchanged. Second of all, perturbations that leave the supersymmetry algebra and the Hilbert space unchanged will preserve the Witten index, as is clear from the definition: where the trace runs over the entire Hilbert space. Its absolute value gives a lower bound to the number of ground states. For the decorated graphs with quenched disorder the Witten index is equal to the number of ground states, since all ground states have the same number of particles. It follows that supersymmetry preserving perturbations of this model that leave the Hilbert space unchanged can increase, but not reduce the ground state degeneracy. Due to the hard-core constraint, one can easily change the Hilbert space, for example by adding an edge, so in some sense the choice of the graph is fine tuning. Finally, we should also comment on supersymmetry breaking perturbations. In the spirit of our disordered model, we consider the possibility of adding supersymmetry breaking perturbations locally. If we consider a dilute system of triangles, then local supersymmetry breaking perturbations should lift the degeneracy associated with triangles near such perturbations. On the other hand, the large gap above the ground state manifold should not be immediately destroyed. Indeed, in the limits of small and large y and dilute number of triangles, the system is roughly equivalent to t free spins in a non-uniform field i.e. each spin has its own magnetic field. However, the case when the triangles are not dilute is potentially much more interesting since the ground state manifold can now support interesting dynamics below the gap of the unperturbed Hamiltonian. It would be interesting to ask what sorts of linear combinations of the VBS states could be selected by such supersymmetry breaking perturbations and if the resulting combinations had any interesting entanglement properties.

B. Comparison to holographic studies
We briefly comment on the potential relations of our work to holographic studies. First, the models we construct have degenerate ground states while holographically the situation is less clear. The AdS 2 black hole certainly has an extensive ground state entropy, but this represents an extreme case in [19]. Second, the ultimate low energy fate of the general class of models considered in [19] is often unknown. As demonstrated by those authors, in some cases the IR fate of the RG flow is known and known to satisfy a boundary law. More generally, it may be that the holographic computations should be interpreted as intermediate scale descriptions, a possibility discussed in [20] where it was termed "IR incompleteness". Certainly we cannot rule out ground states that display anomalous entanglement properties over a wide range of scales so long as they eventually satisfy a boundary law (up to logarithmic violations) at the largest scales. This has been a central theme of our construction e.g. by focusing on the entangling power of finite depth quantum circuits and in the fractal Fermi surface example (see App. A).
Another interesting issue is that it is unclear what state is being used in the holographic calculation. In the case of the extremal black hole, the state is either the completely mixed state or a coherent superposition designed to mimic the mixed state for subsystems less than half the total system size. This is analogous to the gravitational distinction between an eternal black hole and a black hole due to collapse in the finite past. The extremal black hole is a rather special object in this regard. We showed above that both of these types of states can give anomalous entropy in our model, although it should only be called entanglement in the case of the coherent superposition. The only type of state that does not give anomalous entropy in our model is a single localized ground state or non-generic superpositions of such states i.e. with triangles only entangled near the boundary.
One important distinction between our model and the holographic systems comes when considering the thermal entropy to entanglement entropy crossover. We have an exponential number of strictly degenerate ground states separated from other states via a large gap. The temperature dependence of the thermal entropy will thus be of the usual type. However, the thermal physics in the holographic setup looks more like the codimension less than one scenario discussed in [20] and, we suspect, like the thermal physics of the fractal Fermi surface over a wide range of temperatures. Since there are no known examples of field theories with codimension less than one, it seems that the ultimate IR fate of the ground state is to violate the boundary law at most logarithmically. In a sense this may also be the ultimate fate of our model if one assumes that a generic ground state must be connected to the product state via a quantum circuit of finite depth.
Finally, we note that the holographic models are manifestly translationally invariant. For the models we consider, however, it is easily seen that imposing translational invariance on the lattice leads to a ground state degeneracy that is exponential in the volume of the system. We also considered restoring translational invariance for the models with anomalous scaling of the ground state degeneracy by performing a disorder average and we found that this leads to a volume law for the entanglement entropy.

VII. CONCLUSION
The entanglement entropy of a many-body ground state of a local Hamiltonian typically obeys the area law. Violations of the area law may occur in the presence of gapless excitations, long-range interactions, infinite fine-tuning or ground state degeneracy. Gapless excitions, however, are believed to cause at most logarithmic violations of the area law [20]. In this paper we presented a model with short-range interactions on a d dimensional graph that can exhibit anomalous entanglement entropy scaling in a wide range of length scales in the absence of gapless excitations and infinite fine-tuning. The key ingredient to obtain this type of behavior is a controllable ground state degeneracy that results from quenched disorder and is protected by supersymmetry.
We showed that in the limits of small and large staggering there is a basis in which the degenerate ground states take on a simple form. In these limits the degeneracy can be understood in terms of local impurities and we found that these basis elements, which are maximally localized, each obey an area law for the entanglement entropy. We then showed that a generic ground state, which will typically be a superposition of the maximally localized states, can exhibit an anomalous entanglement entropy scaling. More precisely, the ground state degeneracy can be chosen such that the entanglement entropy of a generic ground state grows as L a with d − 1 < a < d in a range of region sizes with linear size L. We carefully addressed the issues of what characterizes a generic ground states and what constraints the range of region sizes in which this anomalous entanglement entropy scaling is observed.
It is clear that this anomalous behavior occurs because the model we presented is not generic. We argue, however, that the ground state degeneracy arises in a reasonably sensible many-body Hamiltonian. We do this by scoring the model against a list of criteria that we impose on microscopic models of quantum many-body systems in order to be realistic and generic. The criteria concern the range of interactions, dimensionality, randomness and fine-tuning. Finally, we discussed the possible relation of this work to holographic models that show anomalous entanglement entropy scaling of the same type. We believe an important question for the holographic models is for which state the holographic entropy computation is done. This is possibly related to the question what ultimately identifies the ground state of our model, or at least its entanglement entropy properties. Related to these issues is the question if there are systems for which anomalous violations of the area law persist on all length scales.
Thus for a range of scales > L > 1/K F (0) we expect to find an anomalous violation of the boundary law for entanglement entropy. This illustrates the power of long range hoppings and fine tuning for producing anomalous entanglement. Note also that while this is not a rigorous argument, we believe it does correctly indicate a violation of the boundary law although we cannot vouch for the details of the scaling form. This kind of fractal Fermi surface in some sense realizes the codimension less than one situation discussed in [20] and effectively seen in holographic studies [19]