This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Topological entanglement entropy meets holographic entropy inequalities

Joydeep Naskar [email protected] Department of Physics, Northeastern University, Boston, MA 02115, USA The NSF AI Institute for Artificial Intelligence and Fundamental Interactions, Cambridge, MA, U.S.A.    Sai Satyam Samal [email protected] Department of Physics and Astronomy, Purdue University, West Lafayette, IN 47907, USA
(January 17, 2025)
Abstract

Topological entanglement entropy (TEE) is an efficient way to detect topological order in the ground state of gapped Hamiltonians. The seminal work of Kitaev and Preskill Kitaev and Preskill (2006) and simultaneously by Levin and Wen Levin and Wen (2006) proposed information quantities that can probe the TEE. In the present work, we explain why the subtraction schemes in the proposed information quantities Levin and Wen (2006); Kitaev and Preskill (2006) work for the computation of TEE and generalize them for arbitrary number of subregions by explicitly noting the necessary conditions for an information quantity to capture TEE. Our conditions differentiate the probes defined by Kitaev-Preskill and Levin-Wen into separate classes. While there are infinitely many possible probes of TEE, we focus particularly on the cyclic quantities 𝒬2n+1\mathcal{Q}_{2n+1} and multi-information n\mathcal{I}_{n}. We also show that the holographic entropy inequalities are satisfied by the quantum entanglement entropy of the non-degenerate ground state of a topologically ordered two-dimensional medium with a mass gap.

I Introduction

Understanding topological phase transitions and detecting the topological order have been one of the most exciting research directions since the discovery of the fractional quantum hall (FQH) effect Tsui et al. (1982). A topological phase transition falls beyond the Landau-Ginzburg paradigm and hence requires a new approach Wen (2007). A topological phase is protected by a finite energy gap at zero temperature Wen (2007) and closing the gap by continuously tuning the parameters in the Hamiltonian would lead to a different topological phase Alicea (2012). These exotic phases are associated with the presence of gapped excitations called anyons Leinaas and Myrheim (1977); Wilczek (1982); Arovas et al. (1984); Khare (2005). These quasi-particles have braiding statistics that is intermediate between bosons and fermions Kitaev (2006); Nayak et al. (2008); Rao (2016) and have recently been experimentally verified in FQH liquids Nakamura et al. (2020); Bartolomei et al. (2020); Ruelle et al. (2023).

However, we must point out that looking at anyonic excitations is not the only way to characterize the topological phase. Ground state degeneracy also provides another quantum number that characterizes different topological phases. It was shown that the ground state degeneracy of FQH states (filling factor ν=1m\nu=\frac{1}{m}) goes as mgm^{g} on a Reimann manifold with gg genus Wen and Niu (1990) and is non-degenerate on a sphere Haldane (1983). Interestingly, the Von-Neumann (VN) entanglement entropy of the ground state also gives information about the topological phase Kitaev and Preskill (2006); Levin and Wen (2006). The VN entropy of gapped ground state in the thermodynamic limit goes as the area of the system Srednicki (1993); Eisert et al. (2010) with correction, γ\gamma, that is universal and independent of the area of the sub-system,

S(ρ)=αLγ+.\displaystyle S(\rho)=\alpha L-\gamma+\cdots\,. (1)

The ellipsis terms vanish as the system size goes to infinity. The universal constant, γ\gamma, captures the non-local feature of the ground state and encodes information about the topological phase Kitaev and Preskill (2006).

Intuitively, in order to probe the topological entanglement entropy (TEE), one can use subtraction scheme similar to Ref. Kitaev and Preskill (2006); Levin and Wen (2006) which is invariant under any smooth deformations. It cancels out all the boundary contributions in Eq. (1) and leaves only the non-local contributions (proportional to γ\gamma) to the TEE. This motivates us to ask the following question: are there other information quantities that one can write in order to probe the TEE and thereby detecting topological order and if yes, what are the general properties of such information quantities? By borrowing ideas from holography, we propose new information quantities that can probe the TEE and also comment on the necessary conditions which must be satisfied for being a topological information quantity. We will be restricting ourselves to case where the ground state is non-degenerate or equivalently our system is on a disk in two spatial dimensions. Since the ground state of the gapped Hamiltonian can be described by an effective topological quantum field theory (TQFT). We will be presenting a TQFT calculation on a disk-like geometry along the lines of Ref. Kitaev and Preskill (2006) and present these new information quantities that can probe the TEE.

I.1 Definition of topological entanglement entropy

In the literature, there is no consensus on the precise definition of TEE. We will review the three commonly used definitions as follows:

  • The first definition comes from the area law, Eq. (1). The additive constant term γ<0-\gamma<0 is defined as the TEE, and is believed to capture some global features of entanglement in the ground state. We will use this quantity γ\gamma as the definition of TEE, i.e, Stopo=γS_{\text{topo}}=\gamma. Our sense of ‘topological’ means that StopoS_{\text{topo}} depends only on global features and not the system’s dimensions, when talking about area law. In the context of anyons, we mean that StopoS_{\text{topo}} only depends on the total quantum dimension 𝒟\mathcal{D} and γ=log(𝒟)\gamma=\log{\mathcal{D}}. It is in some sense a strong assumption that the ground state is described by a TQFT. Nevertheless, we will proceed with this.

  • The second definition comes from the work of Levin and Wen Levin and Wen (2006) where TEE is defined as a tripartite information quantity called the conditional mutual information, Eq. (2),

    StopoLW=\displaystyle S_{\text{topo}}^{\text{LW}}= S(AB)+S(BC)S(B)S(ABC),\displaystyle S(AB)+S(BC)-S(B)-S(ABC)\,, (2)

    computed on a geometry as shown in Fig. 1b. In the Eq. (2), S(X)S(X) is the VN entropy associated with region XX. However, on changing this particular geometry (for example, permuting labels of regions or choosing a disk) renders this quantity useless. We would like to point out that the information quantity used in this case i.e., conditional mutual information is related to the entropy inequality called the strong sub-additivity (SSA) Lieb and Ruskai (1973).

  • The work of Kitaev and Preskill Kitaev and Preskill (2006) motivated a third definition of the TEE involving a different tripartite information quantity related to an entropy inequality called monogamy of mutual information (MMI) Hayden et al. (2013) as follows,

    StopoKP=\displaystyle S_{\text{topo}}^{\text{KP}}= S(AB)+S(AC)+S(BC)\displaystyle S(AB)+S(AC)+S(BC) (3)
    S(A)S(B)S(C)S(ABC),\displaystyle-S(A)-S(B)-S(C)-S(ABC),

    where we choose the sign convention to set Stopo0S_{\text{topo}}\geq 0. This was first evaluated on a disk (see Fig. 1a and a non-simply connected geometry, Fig. 1b. It was found to be invariant under permutation of labels of regions. The topological property of StopoKPS_{\text{topo}}^{\text{KP}} holds for different geometries, however, the value StopoKPS_{\text{topo}}^{\text{KP}} depends on the topology (more specifically, on the number of connected components of the geometry). We call such quantities fixed-topology TEE probes. The Ref. Bao et al. (2015a) gives an explicit formula of StopoKPS_{\text{topo}}^{\text{KP}} in terms of the zeroth Betti number.

    Refer to caption
    Figure 1: (a) Geometry used by Kitaev and Preskill Kitaev and Preskill (2006) for computing the TEE. The information quantity used here comes from the entropy inequality called the monogamy of mutual information (MMI) Hayden et al. (2013). (b) Geometry used by Levin and Wen Levin and Wen (2006) for obtaining the TEE. In this case the information quantity used is known as conditional mutual information and it is related to entropy inequality known as strong sub-additvity (SSA) Lieb and Ruskai (1973).

I.2 Summary of results

A holography-inspired approach Bao et al. (2015a) showed that the cyclic family of holographic entropy inequalities (HEI) Bao et al. (2015b) gives the TEE, γ\gamma. This was demonstrated by explicitly by applying the area law, Eq. (1). In this current work, we extend the previously done analysis Bao et al. (2015a, b) in the following sense,

  • We don’t explicitly use the area law, Eq. (1) and instead perform a TQFT calculation on a 2d disk-like geometry. We demonstrate that the TEE depends only on the total quantum dimension 𝒟\mathcal{D}, in agreement with the results of Ref. Kitaev and Preskill (2006); Levin and Wen (2006); Bao et al. (2015a).

  • We show that all HEI Bao et al. (2024a) are satisfied by gapped Hamiltonians with unique ground state, Eq. (1) and, especially the facet inequalities, can be used to evaluate the TEE, γ\gamma.

  • We further show that a certain class of entanglement entropy inequalities that need not hold true for holographic quantum states, could hold true for the 2d-disk geometry Kitaev and Preskill (2006). Hence, these entropy relations can also be used as probes for TEE.

  • We contrast the definitions of TEE in Eq. (2) and Eq. (3) and classify all possible definitions of StopoS_{\text{topo}} into two classes, suggesting necessary conditions on the subtraction scheme Kim et al. (2023a); Levin (2024) of entanglement entropies.

The organization of this paper is as follows: In section, Sec. II, we begin by reviewing the necessary TQFT tools for our discussion. Our main results are in section, Sec. III, where we consider the HEIs in the context of TEE and use these results in subsection IV.1 to discuss entanglement-based probes of topological phases. Finally, we discuss our results and outline some future directions in section, Sec. IV. Throughout this paper, we will neglect the spurious TEE which often arise in the symmetry protected topological phases Zou and Haah (2016); Santos et al. (2018); Williamson et al. (2019); Stephen et al. (2019); Kato and Brandão (2020); Kim et al. (2023b). However, we will bring them into the discussion towards the end in the section, Sec. IV.

II Necessary TQFT techniques

In this section, we review the calculation for TEE (StopoS_{\text{topo}}) for a 2d disk-like geometry with 33 subregions, Fig. 1Kitaev and Preskill (2006). The procedure involves stitching together a time-reversal (TR) copy of the system with itself and making punctures at each intersection point Kitaev and Preskill (2006). Each puncture can host anyon of a particular type. As a result of making punctures, the VN entropy of the subregion AA (SAS_{A}), (see Fig. 1a) is proportional to the VN entropy of the sphere with 33 punctures (𝒮3\mathcal{S}_{3}), i.e. 2SA=𝒮32S_{A}=\mathcal{S}_{3}. The factor of 22 comes from the fact that we have doubled the system by attaching a TR copy of the system with itself. Note that, in this procedure we don’t create any anyons and therefore the net anyonic charge is zero Kitaev and Preskill (2006). Hence, the VN entropy of the sphere with 33 punctures (𝒮3\mathcal{S}_{3}) is equal to the VN entropy of system with 33 anyonic excitation such that the the total charge vanishes.

Now, we proceed with the calculations for the VN entropy for a sphere with punctures. For concreteness, let us assume that we have a sphere with 44 punctures and later we give a general formula for the sphere with nn punctures. Each puncture hosts anyons and let us denote the anyons by aa, bb, cc and dd with quantum dimension did_{i} where i=a,b,c,di=a,b,c,d. Since, we have the constraint that the total anyonic charge in the system is zero, this implies that the VN entropy 𝒮4\mathcal{S}_{4} is given as,

𝒮4=abcdμ=1Nabcd𝟙Pabcd𝟙Nabcd𝟙log(Pabcd𝟙Nabcd𝟙missing),\displaystyle\mathcal{S}_{4}=-\sum_{abcd}\sum_{\mu=1}^{N_{abcd}^{\mathds{1}}}\frac{P_{abcd}^{\mathds{1}}}{N_{abcd}^{\mathds{1}}}\log\bigg(\frac{P_{abcd}^{\mathds{1}}}{N^{\mathds{1}}_{abcd}}\bigg{missing})\,, (4)

where Pabcd𝟙P_{abcd}^{\mathds{1}} refers to the probability of fusing aa, bb, cc and dd to identity (𝟙\mathds{1}) and Nabcd𝟙N_{abcd}^{\mathds{1}} denotes the number of ways all four fuses to identity (𝟙\mathds{1}). Here, identity (𝟙\mathds{1}) refers to vacuum or an anyon with trivial anyonic charge. By using the fact that the probability of fusing anyons of type aa and bb to give cc is equal to Pabc=NabcdcdadbP_{ab}^{c}=\frac{N_{ab}^{c}d_{c}}{d_{a}d_{b}} Preskill (2004), one obtains the probability of fusing aa, bb, cc and dd to identity as follows,

Pabcd𝟙=Nabcd𝟙dadbdcdd𝒟6,\displaystyle P_{abcd}^{\mathds{1}}=\frac{N^{\mathds{1}}_{abcd}d_{a}d_{b}d_{c}d_{d}}{\mathcal{D}^{6}}\,, (5)

where 𝒟=idi2\mathcal{D}=\sqrt{\sum_{i}d_{i}^{2}} is the total quantum dimension. Substituting Eq. (5) into Eq. (4) and using the expression for the fusion probability, PabcP_{ab}^{c} one obtains the final expression for the VN entropy as,

𝒮4=6logD4aPalogda,\displaystyle\mathcal{S}_{4}=6\log D-4\sum_{a}P_{a}\log d_{a}\,, (6)

where Pa=da2𝒟2P_{a}=\frac{d_{a}^{2}}{\mathcal{D}^{2}} is interpreted as the probability of finding an anyon of type aa in a gas of anyons in thermal equilibrium where effectively no more fusions between any pair of anyons are taking place Preskill (2004). One can generalize the obtained formula for VN entropy of a sphere with 44-punctures to sphere with nn-punctures,

𝒮n=2(n1)log𝒟n𝒦,\displaystyle\mathcal{S}_{n}=2(n-1)\log\mathcal{D}-n\mathcal{K}\,, (7)

where we have defined 𝒦=aPalogda\mathcal{K}=\sum_{a}P_{a}\log d_{a} which we will use in the subsequent calculations.

Now, we can proceed with the calculation for the StopoS_{\text{topo}} in the geometry proposed in Ref. Kitaev and Preskill (2006) with three regions AA, BB and CC. One can probe the TEE (StopoS_{\text{topo}}) by computing the multi-information (we will give a standard definition later), Stopo=I(A:B:C)3=2S3+32S4=log𝒟S_{\text{topo}}=-I(A:B:C)\equiv-\mathcal{I}_{3}=-2S_{3}+\frac{3}{2}S_{4}=\log\mathcal{D} Kitaev and Preskill (2006). Hence, TEE only depends on the total quantum dimension. We would like to generalize the calculation for a disk with an arbitrary number of subregions. We begin by considering a disk with 55-subregions and we will calculate 5I(A:B:C:D:E)\mathcal{I}_{5}\equiv-I(A:B:C:D:E). Heuristically, one can verify 5\mathcal{I}_{5} also successfully probes the TEE by plugging Eq. (1) into the formula for 5\mathcal{I}_{5} and observing that all the boundary terms cancels exactly. On a physical ground we expect that we get the same StopoS_{\text{topo}} as we obtained in the case with 33-subregions by computing 3\mathcal{I}_{3}. This is because the underlying topology is unchanged and increasing the number of subregions should not affect the TEE.

However, we would like to do a TQFT calculation to show that the we obtain the same TEE in the case of 55-subregions. Let us begin by labeling the 55-subregions as AA, BB, CC, DD and EE. We will now point out a subtle issue while probing the TEE using 5\mathcal{I}_{5}. The expression for 5\mathcal{I}_{5} involves VN entropy of two types of regions, which we call contagious (e.g. AB) and non-contagious regions (e.g. AC). Contagious regions share a boundary while non-contagious regions have no boundary in common. In the case where we have 33-subregions, there is no possibility of non-contagious regions and hence both SABS_{AB} and SACS_{AC} were identified with 𝒮4\mathcal{S}_{4}. With 55-subregions, it turns out that one cannot simply identify SACS_{AC} with S5S_{5} and SABS_{AB} with S4S_{4} and carry out the calculations. The reason is, making a puncture at point where 33-subregions meet is different from the point where 55-subregions meet. We can circumvent the problem by smoothly deforming the geometry such that we only have intersection of at most 33-subregions, as shown in Fig. 2b.

Refer to caption
Figure 2: Geometry considered for computing the topological entanglement entropy (TEE) with 55 subregions. By stitching a time-reversal (TR) copy of the system at spatial infinity and making punctures at the intersection points of subregions, one can map the Von Neumann (VN) entropy of a subregion to the VN entropy of sphere with nn number of punctures Kitaev and Preskill (2006). (a) Disk partitioned into 55 subregions with all the subregions meeting at a given point. In order to compute the TEE, one first needs to smoothly deform it to a new geometry which have intersection of at most 33 subregions. Note that this smooth deformation doesn’t give rise to any topological phase transition and the system remains in the same topological phase. (b) The new deformed geometry with intersection of at most 33 subregions. By introducing punctures at the intersections of the subregions after gluing with TR copy, one finds the TEE, Stopo=log𝒟S_{\text{topo}}=\log\mathcal{D}.

A smooth deformation doesn’t introduce any topological phase transition and system remains in the same topological phase (protected by a gap). This smooth deformation is equivalent to choosing a new set of subregions and doesn’t affect the underlying Hamiltonian. By considering small deformations to the geometry, Fig. 2b similar to Ref. Kitaev and Preskill (2006) and using the Eq. (1) one finds that 5\mathcal{I}_{5} probes the StopoS_{\text{topo}} successfully. This enables us to compute the TEE using the five-party multi-information, 5\mathcal{I}_{5}. A careful TQFT calculation reveals that, Stopo=5=S3S432S8+S7+2S6=log𝒟S_{\text{topo}}=-\mathcal{I}_{5}=-S_{3}-S_{4}-\frac{3}{2}S_{8}+S_{7}+2S_{6}=\log\mathcal{D} Note (1). This is the same as we obtained using disk with 33-subregions and computing 3\mathcal{I}_{3}, which we expect as the topology is the same.

III Holographic entropy inequalities in TQFT

We will briefly review the idea of area law in holography to motivate the holographic entropy inequalities. The AdSd+1/CFTdAdS_{d+1}/CFT_{d} correspondence Maldacena (1998) is a duality between a quantum theory of gravity in d+1d+1 dimensional Anti-de Sitter (AdS)(AdS)-spacetime and a conformal field theory (CFT) on its dd dimensional boundary. The entanglement entropy S(A)S(A) of a subregion AA on the boundary is related to the area of the minimal surface ΓA\Gamma_{A} in the bulk homologous to AA, and is given by the Ryu-Takayanagi formula Ryu and Takayanagi (2006a, b),

S(A)=|ΓA|4GN,S(A)=\frac{|\Gamma_{A}|}{4G_{N}}\,, (8)

to the leading order, where GNG_{N} is Newton’s constant.

We will refer to the quantum states on the boundary CFT that have a semi-classical dual in the bulk and obey the RT formula, Eq. (8) as holographic states. These states are a proper subset of all quantum states on the boundary CFT. The goal of the holographic entropy cone program Bao et al. (2015b); Czech et al. (2024); He et al. (2019); Hernández Cuenca (2019); He et al. (2020); Hernández-Cuenca et al. (2023); Bao et al. (2024b); Bao and Naskar (2024); Bao et al. (2024a) is to non-trivially constrain these multipartite holographic states. This program has led to the discovery of many holographic entropy inequalities and here we will study those inequalities in the context of topological entanglement entropy.

III.1 The holographic entropy cone

Consider a five-party information quantity, 𝒬5\mathcal{Q}_{5}, which is given by Eq. (9). The quantity 𝒬5\mathcal{Q}_{5} along with other information quantities formed using the cyclic family of holographic inequalities were shown to give the TEE Bao et al. (2015a) by assuming the area law, Eq. (1). Here, we verify the claim by calculating 𝒬5\mathcal{Q}_{5} explicitly using TQFT, Eqs. (10)–(11).

𝒬5\displaystyle\mathcal{Q}_{5} =S(ABC)+S(BCD)+S(CDE)+S(DEA)+S(EAB)S(AB)S(BC)S(CD)S(DE)S(EA)\displaystyle=S(ABC)+S(BCD)+S(CDE)+S(DEA)+S(EAB)-S(AB)-S(BC)-S(CD)-S(DE)-S(EA)
S(ABCDE),\displaystyle\ \ \ -S(ABCDE)\,, (9)
=12[3S7+2S63S63S5],\displaystyle=\frac{1}{2}\big{[}3S_{7}+2S_{6}-3S_{6}-3S_{5}\big{]}\,, (10)
=12[36log𝒟21𝒦+20log𝒟12𝒦30log𝒟+18𝒦24log𝒟+15𝒦]=log𝒟.\displaystyle=\frac{1}{2}\big{[}36\log\mathcal{D}-21\mathcal{K}+20\log\mathcal{D}-12\mathcal{K}-30\log\mathcal{D}+18\mathcal{K}-24\log\mathcal{D}+15\mathcal{K}\big{]}=\log\mathcal{D}\,. (11)

In fact, it is straightforward to do a similar proof for the family of cyclic inequalities, see Sec. S2 of SM 111See Supplementary Material for details.. Having illustrated an example, we will now generalize our claims about HEI. We will begin by giving some definitions.

Definition III.1 (Holographic entropy inequality).

A holographic entropy inequality (HEI) is an information quantity 𝒬\mathcal{Q} on nn-subregions, such that all nn-party holographic states satisfy

𝒬=i=12n1aiSi0,\mathcal{Q}=\sum_{i=1}^{2^{n}-1}a_{i}S_{i}\geq 0, (12)

where SiS_{i} are the entanglement entropy of subregions accompanied by integer coefficients aia_{i}.

Given nn parties, one can construct an entropy basis of all possible (2n12^{n}-1) subregion entanglement entropies and order it lexicographically. For example, if n=3n=3 and the subregions are {A,B,C}\{A,B,C\}, then the entropy basis is {SA,SB,SC,SAB,SAC,SBC,SABC}\{S_{A},S_{B},S_{C},S_{AB},S_{AC},S_{BC},S_{ABC}\}. An entropy vector is a vector in this basis. We will call it a holographic entropy vector when the system is a holographic system, i.e, the nn subregions are chosen on a boundary X\partial X of a bulk manifold XX. The holographic entropy cone is the space of all allowed holographic entropy vectors. It is a convex, rational polyhedral cone, which can be characterized by the facets of the cone. Thus, knowing all the facet HEIs is sufficient for most purposes and any other non-facet HEI will be redundant.

Definition III.2 (Facet HEI).

A holographic entropy inequality 𝒬\mathcal{Q} involving nn-parties is a facet of the HEC iff there exists a codimension-1 set of linearly independent holographic entropy vectors which saturate 𝒬\mathcal{Q}. We will refer to such HEIs as ‘facet HEIs’.

We will now define the nn-partite information quantity n\mathcal{I}_{n} introduced already, which we will call multi-information. Given nn regions AiA_{i}, where i=1,,ni=1,\dots,n, we define n\mathcal{I}_{n} to be,

n=\displaystyle\mathcal{I}_{n}= inS(Ai)i<jnS(AiAj)+\displaystyle\sum_{i}^{n}S(A_{i})-\sum_{i<j}^{n}S(A_{i}A_{j})+\cdots
+(1)n+1S(A1An),\displaystyle+(-1)^{n+1}S(A_{1}\cdots A_{n})\,, (13)

which for up to n=3n=3 are given by (using shorthand S(Ai)SAiS(A_{i})\equiv S_{A_{i}}),

1=SA1;2=SA1+SA2SA1A2,3=SA1+SA2+SA3SA1A2SA1A3SA2A3+SA1A2A3.\begin{split}\mathcal{I}_{1}=&S_{A_{1}}\ ;\ \mathcal{I}_{2}=S_{A_{1}}+S_{A_{2}}-S_{A_{1}A_{2}}\,,\\ \mathcal{I}_{3}=&S_{A_{1}}+S_{A_{2}}+S_{A_{3}}-S_{A_{1}A_{2}}-S_{A_{1}A_{3}}-S_{A_{2}A_{3}}\\ &+S_{A_{1}A_{2}A_{3}}\,.\end{split} (14)

Note that n\mathcal{I}_{n} is generally sign indefinite for n>3n>3 in holographic systems and the results from Kitaev and Preskill (2006) indicates that 3\mathcal{I}_{3} is sign-definite. Here, we find that n\mathcal{I}_{n} is sign definite and equals log(𝒟)-\log{\mathcal{D}}, for all n3n\geq 3. In the Sec. S1 of SM Note (1) we explicitly calculate 3\mathcal{I}_{3}, 4\mathcal{I}_{4} and 5\mathcal{I}_{5}. We introduce another quantity m,n\mathcal{I}_{m,n}, which is defined same as m\mathcal{I}_{m} but for mnm\leq n, where nn is the total number of subregions, such that n,nn\mathcal{I}_{n,n}\equiv\mathcal{I}_{n}. For example, consider a system with 55 subregions (A1A_{1}, A2A_{2}, A3A_{3}, A4A_{4} and A5A_{5}) and arbitrarily choose 33 subregions, say, A1A_{1}, A2A_{2} and A3A_{3} then 3,5=3\mathcal{I}_{3,5}=\mathcal{I}_{3}, Eq. (14). The set of multi-information {i,n}i=3n\{\mathcal{I}_{i,n}\}_{i=3}^{n} forms a basis for the facet HEI He et al. (2019) involving nn-parties (excluding subaddivity). With the above definitions, we can now formulate our following conjectures and propositions.

Conjecture III.1.

n\mathcal{I}_{n} is topological for n3n\geq 3 on a 2d2d disk geometry divided into nn subregions.

Conjecture III.2.

All m,n\mathcal{I}_{m,n} are topological for 3mn3\leq m\leq n on a 2d2d disk geometry divided into nn subregions.

See the supplementary material for more discussions around the conjectures III.1 and III.2.

Proposition III.1.

All facet HEIs of the HEC (n3)(n\geq 3) are topological, when evaluated on a disk using the quantum entanglement entropy in the unique ground state of a topologically ordered two-dimensional medium with a mass gap.

Proof.

From conjecture III.2, we have that all m,n\mathcal{I}_{m,n} are topological. All facet HEIs can be expressed in the \mathcal{I}-basis He et al. (2019). Since the basis elements m,n\mathcal{I}_{m,n} is topological. Therefore, all facet HEIs are topological. ∎

Proposition III.2.

The sum of coefficients {ai}\{a_{i}\} of all facet HEIs, except subadditivity(SA) is negative.

This is proved in Ref. Bao et al. (2024a). For more details, see sec. S3 of SM.

Proposition III.3.

All facet HEIs of the HEC (n3)(n\geq 3) are satisfied by the quantum entanglement entropy in the unique ground state of a topologically ordered two-dimensional medium with a mass gap.

Proof.

From proposition III.1, we have all facet HEIs are topological. We will now prove that they are also non-negative. As follows from sections II and S1, we have γ=log(𝒟)\gamma=\log{\mathcal{D}} on the 2d2d-disk geometry.

Since we have a connected geometry, each term contributes a γ-\gamma. This is because each term has the form given by Eq. (7), and due to topological nature of HEIs shown in proposition III.1, the nn-dependence of the HEI cancels out, leaving behind 2γ-2\gamma, which we divide by 22 to account for the TR copy.

By using Lemma III.2 we have the sum of coefficients is negative for all facet HEIs of the HEC (n3)n\geq 3) Bao et al. (2024a).

Therefore, the HEIs are strictly positive in this geometry. ∎

We will take a step further and discuss the case when these inequalities are saturated. Now consider the 2d planar geometry where each region is disconnected. In this case, each kk-party entanglement term contributes kk number of γ-\gamma terms as there are kk boundaries. We know from He et al. (2020) that all facet HEIs of the HEC are balanced, i.e, each party appears equal number of times with positive and negative coefficients. As a result, the topological terms cancel and the HEIs are strictly zero.

We have an alternative proof for the same based on conjecture III.3 stated after lemma III.1. We give a third proof of this argument in sec. S3 of SM exploiting the structural forms of HEIs. We have verified numerically that this holds for all known superbalanced facet HEIs.

III.2 What about non-facet inequalities?

We will now look at the true HEIs that are not facets of the HEC. For example, consider the following information quantity,

𝒬6,1\displaystyle\mathcal{Q}_{6,1} =S(AEF)+S(BEF)+S(ADE)+S(ADF)+S(BDE)+S(BDF)+S(ABCD)+S(ABCE)\displaystyle=S(AEF)+S(BEF)+S(ADE)+S(ADF)+S(BDE)+S(BDF)+S(ABCD)+S(ABCE) (15)
+S(ABCF)+S(C)S(ABCEF)S(ABCDF)S(ABCDE)S(AD)S(AE)S(AF)\displaystyle\ \ \ +S(ABCF)+S(C)-S(ABCEF)-S(ABCDF)-S(ABCDE)-S(AD)-S(AE)-S(AF) (16)
S(BD)S(BE)S(BF)S(CDEF),\displaystyle\ \ \ -S(BD)-S(BE)-S(BF)-S(CDEF)\,, (17)
=12[4𝒮94𝒮7+𝒮10𝒮6]=12log𝒟6𝒦(Not topological).\displaystyle=\frac{1}{2}\big{[}4\mathcal{S}_{9}-4\mathcal{S}_{7}+\mathcal{S}_{10}-\mathcal{S}_{6}\big{]}=12\log\mathcal{D}-6\mathcal{K}\ \ \ \textbf{(Not topological)}\,. (18)

Calculating the information quantity 𝒬6,1\mathcal{Q}_{6,1} on the disk-like geometry, Fig. 3 using Eq. (7), we find that 𝒬6,1\mathcal{Q}_{6,1} is not topological, Eq. (18). However, it is possible to find some fixed configuration(s) of geometry where Q6,1Q_{6,1} is topological and gives the TEE. We classify all such information quantities that are topological only under specific choice of geometries as fixed-geometry TEE probes (in contrast to fixed-topology TEE probes). The definition of StopoS_{\text{topo}} using strong sub-additivity in Levin and Wen (2006) is another example of this class. More explicitly, we can check that on a 2d2d-disk, we have

StopoLW\displaystyle S^{\text{LW}}_{\text{topo}} =S(AB)+S(BC)S(B)S(ABC)\displaystyle=S(AB)+S(BC)-S(B)-S(ABC) (19)
=12[2𝒮4𝒮3𝒮3]=2log𝒟𝒦,\displaystyle=\frac{1}{2}\big{[}2\mathcal{S}_{4}-\mathcal{S}_{3}-\mathcal{S}_{3}\big{]}=2\log\mathcal{D}-\mathcal{K}, (20)

whereas the same quantity evaluated on a torus-like geometry with the subregion BB disconnected gives the TEE (see Fig. 1b).

Now we will look at another example where the given inequality does not hold for all holographic states, but the information quantity defined using it, 𝒬6,2\mathcal{Q}_{6,2} turns out to probe the TEE, Eq. (21). The factor of 22 in front of log(𝒟)\log{\mathcal{D}} comes from the sum of coefficients being 2-2 [we will illuminate this point in III.4]

𝒬6,2=\displaystyle\mathcal{Q}_{6,2}= S(ADE)+S(ADF)+S(AEF)+S(BDE)+S(BDF)+S(BEF)+S(CDE)+S(CDF)+S(CEF)\displaystyle S(ADE)+S(ADF)+S(AEF)+S(BDE)+S(BDF)+S(BEF)+S(CDE)+S(CDF)+S(CEF) (21)
+S(ABC)S(AD)S(AE)S(AF)S(BD)S(BE)S(BF)S(CD)S(CE)S(CF)\displaystyle+S(ABC)-S(AD)-S(AE)-S(AF)-S(BD)-S(BE)-S(BF)-S(CD)-S(CE)-S(CF) (22)
2S(DEF)S(ABCDEF),\displaystyle-2S(DEF)-S(ABCDEF)\,, (23)
=\displaystyle= 12[4𝒮93𝒮64𝒮7+𝒮10],\displaystyle\frac{1}{2}\big{[}4\mathcal{S}_{9}-3\mathcal{S}_{6}-4\mathcal{S}_{7}+\mathcal{S}_{10}\big{]}\,, (24)
=\displaystyle= 12[4log𝒟]=2log𝒟(Topological).\displaystyle\frac{1}{2}\big{[}4\log\mathcal{D}\big{]}=2\log\mathcal{D}\ \ \ \textbf{(Topological)}\,. (25)
Refer to caption
Figure 3: A disk with 66-subregions obtained after smooth deformations which is used to obtain the value of the information quantity 𝒬6,1\mathcal{Q}_{6,1}, Eq. (15) and 𝒬6,2\mathcal{Q}_{6,2}, Eq. (21). It turns out that the information quantity 𝒬6,1\mathcal{Q}_{6,1} can not be used to probe the TEE however the latter quantity, 𝒬6,2\mathcal{Q}_{6,2} can probe the TEE.

This is surprising as it naively seems to suggest that the space of allowed entropy vector space of the unique ground state of a gapped TQFT is strictly contained within the HEC. However, the space of allowed entropy vector space in the case of degenerate ground states is expected to be larger than the HEC (and thus those states can violate the HEIs). Alternatively, it might be the case, that inclusion of the spurious entanglement entropy may violate these inequalities that seem to naively hold. We postpone these interesting questions for our future work and instead generalize the above observations below. We will need a few more definitions before we do that.

Definition III.3 (kk-balanced HEIs).

A HEI 𝒬\mathcal{Q} is said to be kk-balanced, iff every kk-clustering of single subregions are balanced, i.e, they appear equal number of times with positive and negative coefficients.

For example, consider the monogamy of mutual information (MMI), Eq. (3) involving three subregions. Each subregion A,B,CA,B,C are balanced, implying 11-balanced, every 22-clusters AB,AC,BCAB,AC,BC are also balanced, implying 22-balanced. However, the 33-cluster ABCABC is not balanced and so on. Note that 11-balance and 22-balance can be alternatively referred as balance and superbalance, respectively. We will make the following conjecture III.3 based on empirical observations.

Conjecture III.3.

Any 22-balanced information quantity 𝒬\mathcal{Q} is topological with indefinite sign.

Now we will state our general proposition applicable to all facet and non-facet inequalities.

Proposition III.4.

Any information quantity 𝒬i\mathcal{Q}_{i} whose sum of coefficients is negative (positive) and satisfies 22-balance, is both topological and sign-definite, on a 2d disk-like geometry.

Proof.

By conjecture III.3, 𝒬i\mathcal{Q}_{i} is topological. If the sum of coefficients is negative (positive), then, each term contributes a γ-\gamma as a result of connected geometry and therefore 𝒬i\mathcal{Q}_{i} is sign-definite. ∎

This explains why facet HEIs (n3)(n\geq 3) are both valid and topological. We will make this more concrete below.

Lemma III.1.

All facet HEIs, except subadditivity (SA) are 22-balanced.

The proof is given in Ref. He et al. (2020).

Since, all facet HEIs (except SA) satisfy lemma III.1 and propositions III.2,III.4, therefore, all facet HEIs are both valid and topological.

One may question the validity of the conjecture III.3. Empirically, we have checked all known (over 2000) holographic entropy inequalities (n3)(n\geq 3) in the 2d disk geometry and found them to be both valid and topological. We have also tested them over many superbalanced information quantities that are holographically false to gather empirical evidence in favour of proposition III.4 and have not found any counterexample.

We will try to reason the validity of conjecture III.3 from the perspective of entanglement entropy in quantum field theories Casini and Huerta (2004, 2023); Hubeny et al. (2019). In a local QFT, the entanglement entropy of any finite region diverges. Balanced information quantities cancel these divergences, however, it treats the purifying region OO as a separate entity. Superbalanced quantities on the other hand treat the purifier OO at the same footing as other regions. This give rise to a scheme-independent cancellation of divergences. It has been argued that superbalanced quantities characterize certain topological properties of the configurations Hubeny et al. (2019). We believe that this picture is captured clearly in the TQFT setting where the total quantum dimension is related to the TEE. In the following section, we will suggest how our results are useful in the context of entanglement-based probes of topological phases.

IV Discussions

IV.1 Entanglement-based probes of topological phases

Topological phase is associated with the existence of anyons. The TEE captures this anyonic content exactly because TEE is a function of the total quantum dimension of a given gapped Hamiltonian. Therefore calculating TEE of a system in a unique ground state is one of the possible pathways to detect topological order. Coming up with information quantities that can capture this long-range correlations, equivalently, the topological order is a non-trivial task. Our analysis sheds light on how one can write down some information quantities that can probe the TEE on 2d2d-disk-like geometry. We showed that an information quantity which is superbalanced (see definition III.3), can be used as a probe for the TEE. Hence, using the proposed prescription one can write down information quantities that can successfully probe the TEE and hence detect any topological order. Thus, a general formula for calculating the TEE is given by a superbalanced information quantity 𝒬\mathcal{Q} can be dubbed as,

𝒬=cγ=clog(𝒟),\mathcal{Q}=-c\gamma=-c\log{\mathcal{D}}, (26)

where cc is the sum of coefficients in 𝒬\mathcal{Q}. If one is interested in obtaining StopoS_{\text{topo}}, then

Stopo=1c𝒬=γ,S_{\text{topo}}=\frac{-1}{c}\mathcal{Q}=\gamma, (27)

where the formula Eq. (26) (and Eq. (27)) is applicable for the 2d2d disk geometry. By changing the topology, especially, connected components of the system, the formula needs to be appropriately modified taking the topology into account.

As shown in the Sec. S1 and Sec. S2 of SM Note (1), the multi-information n\mathcal{I}_{n} and cyclic entropy inequalities 𝒬2n+1\mathcal{Q}_{2n+1} are good candidates for probing TEE . While the cyclic inequalities are motivated from holography, the multi-information doesn’t have a supporting holographic counterpart. Both of these quantities are generalization of the MMI used in Kitaev and Preskill (2006), which is given by 𝒬3\mathcal{Q}_{3} and 3\mathcal{I}_{3}, respectively. The advantage of generalizing these quantities is that if the physical system in consideration is naturally divided into more than three regions, these quantities give a direct approach to calculate the TEE. Another advantage of using other facet HEI probes is that, it might be the case that for a non-trivial topology, one has the pre-factor vanishing for some inequalities having the same sum of coefficients, say cc, in those cases, one has to choose a different probe, having a sum of coefficients, cc^{\prime} (and ccc^{\prime}\neq c).

To reiterate our classification of TEE probes, we have defined all superbalanced information quantities as fixed-topology TEE probes while all balanced information quantities are termed as fixed-geometry TEE probes, for which the 2d disk geometry is not suitable.

IV.2 Lessons from holography

We will take inspiration from holography to comment on multi-information n\mathcal{I}_{n} being the potential signal of multipartite entanglement in gapped Hamiltonians. Recent works in holography Balasubramanian et al. (2024); Gadde et al. (2024) propose various information quantities as measures (and signals) of multipartite entanglement entropy, each with their pros and cons. In particular, holographically, n\mathcal{I}_{n} is zero on separable quantum states and sign-indefinite for entangled states. Unfortunately, it is zero, for some states with true multiparite entanglement. Therefore, it is a good but not the best signal for probing multipartite entanglement in holography. In gapped Hamiltonians with unique ground states (i.e. topologies with 2d-disk), n=γ\mathcal{I}_{n}=-\gamma (for n3n\geq 3), suggestive of nn-party global correlations in the ground state. Of course, if one picks a geometry in which all subregions are disconnected, then n\mathcal{I}_{n} vanishes identically. So, we already know that n\mathcal{I}_{n} is a good signal. It would be interesting to import the alternative definitions of signals of multipartite entanglement from holography and probe them in topological systems. We will report results in this direction in future work.

IV.3 Including spurious entanglement entropies

It was shown previously that, by including spurious contributions to entanglement entropy, the equality γ=log(𝒟)\gamma=\log{\mathcal{D}} ceases to hold Zou and Haah (2016); Santos et al. (2018); Williamson et al. (2019); Stephen et al. (2019); Kato and Brandão (2020); Kim et al. (2023b). One limitation of our proof method is that we have first proved that the information quantities formed using HEIs are topological, i.e., they satisfy the equality γ=log(𝒟)\gamma=\log{\mathcal{D}}, and then use it to show that the HEIs are valid. One may argue that if instead one has an inequality γlog(𝒟)\gamma\geq\log{\mathcal{D}}, are the facet HEIs still valid? Since the sum of coefficients for a facet HEI is negative, if there are physical processes that produce a positive contribution to the entanglement entropy (such as spurious contributions) for every or some subregions, depending on the relative magnitude of those contributions, it is a priori not clear if these inequalities continue to hold and one might expect to see violations similar to the case of quantum extremal surfaces in holography Akers et al. (2021). Perhaps, these exotic systems will shed more light on the nature of the HEC, if they cease to hold true for superbalanced information quantities with negative sum of coefficients that are not valid for holographic states, but continue to hold for facet HEIs. This is a very interesting question, and we will take a deep dive on this question in a future work. The answer will also shed light on the space of quantum states of gapped Hamiltonians with unique ground states.

IV.4 On changing the geometry

Another limitation is that we have used the geometry of a 2d-disk for our proof. Does the proof hold when one considers other topologies? It is clear that the geometry where each region is disconnected from each other, all superbalanced quantities will identically evaluate to zero. However, we would like to believe that the facet HEIs would always give a non-negative value of the entanglement entropy for the subtraction scheme. It would be interesting if there is a geometry that contradicts this. For the case of connected geometries, where some HEIs may yield zero, others may not. In summary, we would like to conjecture that the full holographic entropy cone, inclusive of all facet HEIs is a probe of global correlations in gapped Hamiltonians. We also think that one may be able to come up with geometries that violate some or all of the superbalanced information quantities with negative sum of coefficients that do not hold holographically.

IV.5 Beyond gapped Hamiltonians

Lastly, we would like to learn beyond holography and gapped Hamiltonians, and ask which inequalities hold true, in general, for all quantum states. All quantum entanglement entropy inequalities (QEEI) must be obeyed by the GHZ states. This puts a lower bound on the sum of the coefficients of the QEEIs Bao et al. (2024a) to be zero. On the other hand, the 2d-disk geometry considered in our analysis puts the upper bound on the sum of coefficients to zero for the validity of such inequalities. Therefore, one only needs to look at those inequalities whose sum of coefficients are identically zero. We believe that these inequalities form the boundary that separates inequalities that are valid for topological systems from those that do not. Although one can easily find geometries that violate superbalanced inequalities that do not hold in holography whose sum of coefficients is zero 222We thank Michael Levin for working out this example with us., we leave it as a future exercise, to ascertain if the same is true for inequalities that hold holographically. This may shed light on whether the space of ground states of gapped Hamiltonians is closed or open.

V Acknowledgements

We thank Ning Bao for suggesting this project direction. J.N. would like to thank Ning Bao, Keiichiro Furuya, Eugene Tang and particularly Michael Levin for useful discussions. S.S.S. would like thank Colleen Delaney, Jukka I. Väyrynen, Bowen Yan, Ayush Raj and Soham Ray for discussions. We thank Ning Bao, Eugene Tang and Jukka I. Väyrynen for comments on the draft. J.N. is supported by the NSF under Cooperative Agreement PHY2019786 and the Graduate Assistantship by the Departmenet of Physics, Northeastern University.

References

Supplementary Materials

In the following supplementary materials, we are show the explicit calculations of some results used in the article. In the section, Sec. S1, we review the computation of the topological entanglement entropy (TEE) on disk with a given number of subregions by calculating the multi-information, Eq. (III.1). We show the explicit calculations with 33, 44 and 55 subregions and verified it numerically upto 2020 subregions. Additionally, we also show that one can probe the TEE by considering m,n\mathcal{I}_{m,n}, with mnm\leq n, where nn is the number of subregions. In the section, Sec. S2, we give a TQFT proof that the information quantity associated with cyclic-inequalities can also probe the TEE which was previously done by assuming area law Bao et al. (2015a). In the section sec. S3 we present an alternate proof that the holographic inequalities inequalities are satisfied in the non-degenerate ground states of TQFT.

S1 Multi-information as a probe for topological entanglement entropy

In this section we are going to review the calculation for the topological entanglement entropy along the lines of Ref. Kitaev and Preskill (2006). Here, we will show that the StopoS_{\text{topo}} can be probed by 5\mathcal{I}_{5} in a disk-like geometry with 55 subregions, Fig. 2b. The expression for 5\mathcal{I}_{5} is given as follows,

Stopo(5 subregions)5=\displaystyle S_{\text{topo}}\ (\text{$5$ subregions})\equiv-\mathcal{I}_{5}= SASBSCSDSE\displaystyle-S_{A}-S_{B}-S_{C}-S_{D}-S_{E} (S1)
+SAB+SAC+SAD+SAE+SBC+SBD+SBE+SCD+SCE+SDE\displaystyle+S_{AB}+S_{AC}+S_{AD}+S_{AE}+S_{BC}+S_{BD}+S_{BE}+S_{CD}+S_{CE}+S_{DE} (S2)
SABCSABDSABESACDSACESADESBCD\displaystyle-S_{ABC}-S_{ABD}-S_{ABE}-S_{ACD}-S_{ACE}-S_{ADE}-S_{BCD} (S3)
SBCESBDESCDE\displaystyle-S_{BCE}-S_{BDE}-S_{CDE} (S4)
+SABCD+SABCE+SACDE+SABDE+SBCDE\displaystyle+S_{ABCD}+S_{ABCE}+S_{ACDE}+S_{ABDE}+S_{BCDE} (S5)
SABCDE.\displaystyle-S_{ABCDE}\,. (S6)

It is straightforward to show that in this case, all the boundary terms cancels exactly and under small and smooth deformations Kitaev and Preskill (2006) Stopo5S_{\text{topo}}\equiv\mathcal{I}_{5} probes the topological entanglement entropy, Eq. (1). The procedure to calculate 5\mathcal{I}_{5} involves gluing together a time-reversal copy of the system at spatial infinity and making punctures at the intersection of the subregions. The resultant geometry is shown in the Fig. S1a.

Refer to caption
Figure S1: Geometry used to calculate the topological entanglement entropy (StopoS_{\text{topo}}). (a) Stitching a time reversal copy at spatial infinity and making punctures at the intersections of subregions. (b) subregion AA is now identified with a sphere with 33-punctures whose Von-neumann entropy is calculated using TQFT and is equal to SA=𝒮3=4log𝒟3𝒦S_{A}=\mathcal{S}_{3}=4\log\mathcal{D}-3\mathcal{K}. (c)subregion CDECDE is now identified with sphere with 66 punctures and the Von-Neumann entropy is given as SCDE=𝒮6=10log𝒟6𝒦S_{CDE}=\mathcal{S}_{6}=10\log\mathcal{D}-6\mathcal{K}. Note that in this case only punctures at boundaries are contributing Kitaev and Preskill (2006)

If we look at any subregion, then we find that it is a sphere with nn number of punctures, e.g., region AA is identified with sphere with 33-punctures and the region CDECDE is identified with a sphere with 66-punctures. The Von-Neumann entropy of any subregion can be computed using TQFT and is given by Eq. (7). This implies, we have,

Stopo(5 subregions)5=\displaystyle S_{\text{topo}}\ (\text{$5$ subregions})\equiv-\mathcal{I}_{5}= 12[𝒮3𝒮4𝒮4𝒮3𝒮5\displaystyle\frac{1}{2}\big{[}-\mathcal{S}_{3}-\mathcal{S}_{4}-\mathcal{S}_{4}-\mathcal{S}_{3}-\mathcal{S}_{5} (S7)
+𝒮5+𝒮7+𝒮6+𝒮6+𝒮6+𝒮7+𝒮7+𝒮5+𝒮7+𝒮6\displaystyle+\mathcal{S}_{5}+\mathcal{S}_{7}+\mathcal{S}_{6}+\mathcal{S}_{6}+\mathcal{S}_{6}+\mathcal{S}_{7}+\mathcal{S}_{7}+\mathcal{S}_{5}+\mathcal{S}_{7}+\mathcal{S}_{6} (S8)
𝒮7𝒮8𝒮6𝒮8𝒮8𝒮7𝒮7𝒮7𝒮8𝒮6\displaystyle-\mathcal{S}_{7}-\mathcal{S}_{8}-\mathcal{S}_{6}-\mathcal{S}_{8}-\mathcal{S}_{8}-\mathcal{S}_{7}-\mathcal{S}_{7}-\mathcal{S}_{7}-\mathcal{S}_{8}-\mathcal{S}_{6} (S9)
+𝒮8+𝒮6+𝒮7+𝒮7+𝒮6𝒮5]=log𝒟,\displaystyle+\mathcal{S}_{8}+\mathcal{S}_{6}+\mathcal{S}_{7}+\mathcal{S}_{7}+\mathcal{S}_{6}-\mathcal{S}_{5}\big{]}=\log\mathcal{D}\,, (S10)

by using Eq. (7). This shows that 5\mathcal{I}_{5} also probes the topological entanglement entropy of a disk with 55 subregions. Similar set of calculations can be performed on a disk with 33 and 44 subregions,

Stopo(3 subregions)3=\displaystyle S_{\text{topo}}\ (\text{$3$ subregions})\equiv-\mathcal{I}_{3}= SASBSC+SAB+SBC+SCASABC=12(4𝒮3+3𝒮4),\displaystyle-S_{A}-S_{B}-S_{C}+S_{AB}+S_{BC}+S_{CA}-S_{ABC}=\frac{1}{2}\big{(}-4\mathcal{S}_{3}+3\mathcal{S}_{4}\big{)}\,, (S11)
=\displaystyle= 12(16log𝒟+12𝒦+18log𝒟12𝒦)=log𝒟,\displaystyle\frac{1}{2}\big{(}-16\log\mathcal{D}+12\mathcal{K}+18\log\mathcal{D}-12\mathcal{K}\big{)}=\log\mathcal{D}\,, (S12)
Stopo(4 subregions)4=\displaystyle S_{\text{topo}}\ (\text{$4$ subregions})\equiv-\mathcal{I}_{4}= SASBSCSD+SAB+SAC+SAD+SBC+SBD+SCD\displaystyle-S_{A}-S_{B}-S_{C}-S_{D}+S_{AB}+S_{AC}+S_{AD}+S_{BC}+S_{BD}+S_{CD} (S13)
SABCSABDSBCDSACD+SABCD,\displaystyle-S_{ABC}-S_{ABD}-S_{BCD}-S_{ACD}+S_{ABCD}\,, (S14)
=\displaystyle= 12(2𝒮3𝒮4+2𝒮5)=log𝒟.\displaystyle\frac{1}{2}\big{(}-2\mathcal{S}_{3}-\mathcal{S}_{4}+2\mathcal{S}_{5}\big{)}=\log\mathcal{D}\,. (S15)

Note that all the quantities 3\mathcal{I}_{3}, 4\mathcal{I}_{4}, and 5\mathcal{I}_{5} give the same value for TEE, i.e. Stopo=log𝒟S_{\text{topo}}=\log\mathcal{D}. It is expected because the underlying topology does not change if we increase the number of subregions, therefore, we obtained the same TEE. We verified it numerically for up to n=20n=20 subregions. Next, we shall consider the case where we calculate 3\mathcal{I}_{3} for three arbitrarily chosen subregions out of 55 available subregions on a disk, Fig. S1a. Consider a system consisting of AA, CC and EE, then we can probe the topological entanglement entropy as follows,

Stopo(5 subregions)=3,5(5 subregions)\displaystyle S_{\text{topo}}\ (\text{$5$ subregions})=-\mathcal{I}_{3,5}\ (\text{$5$ subregions}) =12[SA+SC+SESACSAESEC+SAEC],\displaystyle=-\frac{1}{2}\big{[}S_{A}+S_{C}+S_{E}-S_{AC}-S_{AE}-S_{EC}+S_{AEC}\big{]}\,, (S16)
=12[𝒮3+𝒮4+𝒮5𝒮7𝒮6𝒮7+𝒮8],\displaystyle=-\frac{1}{2}\big{[}\mathcal{S}_{3}+\mathcal{S}_{4}+\mathcal{S}_{5}-\mathcal{S}_{7}-\mathcal{S}_{6}-\mathcal{S}_{7}+\mathcal{S}_{8}\big{]}\,, (S17)
=log𝒟.\displaystyle=\log\mathcal{D}\,. (S18)

Now, let us consider a more general case where we have (2n+1)(2n+1)-subregions (odd number of subregions) on a disk-like geometry, Fig. S2a,c. Now, we shall probe the topological entanglement entropy by choosing 33 subregions that are non-contagious, say A2n+1A_{2n+1}, A2A_{2} and A4A_{4}. In this case, we can write the following,

Stopo(2n+1 subregions)=3,2n+1\displaystyle S_{\text{topo}}\ (\text{$2n+1$ subregions})=-\mathcal{I}_{3,2n+1} =SA2n+1+SA2+SA4SA2n+1A2SA2n+1A4SA2A4+SA2n+1A2A4,\displaystyle=S_{A_{2n+1}}+S_{A_{2}}+S_{A_{4}}-S_{A_{2n+1}A_{2}}-S_{A_{2n+1}A_{4}}-S_{A_{2}A_{4}}+S_{A_{2n+1}A_{2}A_{4}}\,, (S19)
=12[𝒮2n+1+𝒮4+𝒮4𝒮2n+3𝒮2n+3𝒮8+𝒮2n+5],\displaystyle=-\frac{1}{2}\big{[}\mathcal{S}_{2n+1}+\mathcal{S}_{4}+\mathcal{S}_{4}-\mathcal{S}_{2n+3}-\mathcal{S}_{2n+3}-\mathcal{S}_{8}+\mathcal{S}_{2n+5}\big{]}\,, (S20)
=log𝒟.\displaystyle=\log\mathcal{D}\,. (S21)

where we have used the general formula for the VN entropy of a sphere with nn-number of punctures, Eq. (7). In all of the examples we demonstrated that the there are many ways to probe the topological entanglement entropy and all of them yields Stopo=log𝒟S_{\text{topo}}=\log\mathcal{D}. Since the underlying geometry is the same, therefore the TEE computed is also same.

S2 Proof for the general case of cyclic inequalities

In this section we are going to look at the general case of cyclic inequalities for (2n+1)(2n+1)-subregions (odd number of subregions) on a disk-like geometry, Fig. S2. Expression for the cyclic inequality can be written as follows,

i=12n+1Sai+i=12n+1Sai+SA,\sum_{i=1}^{2n+1}S_{a_{i}^{+}}\geq\sum_{i=1}^{2n+1}S_{a_{i}^{-}}+S_{A}\,, (S22)

where we have defined A={a1,,a2n+1}A=\{a_{1},\dots,a_{2n+1}\} to be the set of 2n+12n+1 subregions, SA=Sa1a2n+1S_{A}=S_{a_{1}\dots a_{2n+1}} and

ai(k)=aiai+k1,ai±:=ai(2n+1±12)=ain±1.a^{(k)}_{i}=a_{i}\cdots a_{i+k-1},\;a^{\pm}_{i}:=a_{i}^{(\frac{2n+1\pm 1}{2})}=a_{i}^{n\pm 1}\,. (S23)

The corresponding information quantity turns out to be able to probe the topological entanglement entropy, Stopo=i=12n+1Sai+i=12n+1SaiSAS_{\text{topo}}=\sum_{i=1}^{2n+1}S_{a_{i}^{+}}-\sum_{i=1}^{2n+1}S_{a_{i}^{-}}-S_{A}. This was proven in Ref. Bao et al. (2015a) by assuming the area law Eq. (1). We quote their result on the derived formula for γ\gamma below,

Stopo=γ{i=12n+1(b0[(AiAi+n]b0[(Ai+1Ai+n])+b0[(A1A2n+1)]},S_{\text{topo}}=-\gamma\left\{\sum_{i=1}^{2n+1}\left(b_{0}[\partial(A_{i}\dots A_{i+n}]-b_{0}[\partial(A_{i+1}\dots A_{i+n}]\right)+b_{0}[\partial(A_{1}\dots A_{2n+1})]\right\}, (S24)

where b0[A]b_{0}[\partial A] denotes the zeroth Betti number (the number of connected components) of the boundary of a region A, and all indices are taken modulo (2n+1)(2n+1). Now we will give a TQFT proof of the same for a fixed topology, i.e., the 2d disk. The first step involves a continuous deformation of the geometry where we have intersection of at most 33 subregions. Second step involves gluing together a time-reversal copy of the system and making punctures at each intersection point, as shown in the Fig. S2.

Refer to caption
Figure S2: (a) Extracting the topological entanglement entropy of disk geometry with with 2n+12n+1 subregions. (b) First step involves deforming the geometry where we have at most intersection of 33 subregions. (c) Gluing a time-reversal copy of the system and making punctures at the intersection of 33 subregions. As a result of this Von-Neumann entropy of a subregion is equal to Von-Neumann entropy of sphere with punctures where each puncture host a anyon of a particular type.

The final step involves identifying the Von Neumann entropy of the subregions in terms of a sphere with particular number of punctures. This is done as follows,

i=12n+1Sai+\displaystyle\sum_{i=1}^{2n+1}S_{a_{i}}^{+} =Sa1++i=2n1Sai++San++San+1++San+2++i=n+32n1Sai++S2n++Sa2n+1+,\displaystyle=S_{a_{1}}^{+}+\sum_{i=2}^{n-1}S_{a_{i}}^{+}+S_{a_{n}}^{+}+S_{a_{n+1}}^{+}+S_{a_{n+2}}^{+}+\sum_{i=n+3}^{2n-1}S_{a_{i}}^{+}+S_{2n}^{+}+S_{a_{2n+1}}^{+}\,, (S25)
=Sa1++(n12+1)Sai++San++San+1++San+2++(2n1n3+1)San+3+++S2n++Sa2n+1+,\displaystyle=S_{a_{1}}^{+}+(n-1-2+1)S_{a_{i}}^{+}+S_{a_{n}}^{+}+S_{a_{n+1}}^{+}+S_{a_{n+2}}^{+}+(2n-1-n-3+1)S_{a_{n+3}}^{+}++S_{2n}^{+}+S_{a_{2n+1}}^{+}\,, (S26)
=Sa1++(n2)Sa2++San++San+1++San+2++(n3)San+3++S2n++Sa2n+1+,\displaystyle=S_{a_{1}}^{+}+(n-2)S_{a_{2}}^{+}+S_{a_{n}}^{+}+S_{a_{n+1}}^{+}+S_{a_{n+2}}^{+}+(n-3)S_{a_{n+3}}^{+}+S_{2n}^{+}+S_{a_{2n+1}}^{+}\,, (S27)
=12[𝒮2n+3+(n2)𝒮2n+4+𝒮2n+3+𝒮2n+2+𝒮2n+3+(n3)𝒮2n+3+𝒮2n+3+𝒮2n+2],\displaystyle=\frac{1}{2}\big{[}\mathcal{S}_{2n+3}+(n-2)\mathcal{S}_{2n+4}+\mathcal{S}_{2n+3}+\mathcal{S}_{2n+2}+\mathcal{S}_{2n+3}+(n-3)\mathcal{S}_{2n+3}+\mathcal{S}_{2n+3}+\mathcal{S}_{2n+2}\big{]}\,, (S28)
=12[2𝒮2n+2+(n+1)𝒮2n+3+(n2)𝒮2n+4].\displaystyle=\frac{1}{2}\big{[}2\mathcal{S}_{2n+2}+(n+1)\mathcal{S}_{2n+3}+(n-2)\mathcal{S}_{2n+4}\big{]}\,. (S29)

The factor of half comes from the fact that we have doubled the system by gluing it with its time-reversal copy. Similar set of arguments and calculations gives,

i=12n+1Sai\displaystyle\sum_{i=1}^{2n+1}S_{a_{i}}^{-} =Sa1+i=2nSai+San+1+San+2+San+3+i=n+42n1Sai+S2n+Sa2n+1,\displaystyle=S_{a_{1}}^{-}+\sum_{i=2}^{n}S_{a_{i}}^{-}+S_{a_{n+1}}^{-}+S_{a_{n+2}}^{-}+S_{a_{n+3}}^{-}+\sum_{i=n+4}^{2n-1}S_{a_{i}}^{-}+S_{2n}^{-}+S_{a_{2n+1}}^{-}\,, (S30)
=12[𝒮2n+1+(n1)𝒮2n+2+𝒮2n+1+𝒮2n+2+𝒮2n+3+(n4)𝒮2n+3+𝒮2n+3+𝒮2n+2],\displaystyle=\frac{1}{2}\big{[}\mathcal{S}_{2n+1}+(n-1)\mathcal{S}_{2n+2}+\mathcal{S}_{2n+1}+\mathcal{S}_{2n+2}+\mathcal{S}_{2n+3}+(n-4)\mathcal{S}_{2n+3}+\mathcal{S}_{2n+3}+\mathcal{S}_{2n+2}\big{]}\,, (S31)
=12[2𝒮2n+1+(n+1)𝒮2n+2+(n2)𝒮2n+3].\displaystyle=\frac{1}{2}\big{[}2\mathcal{S}_{2n+1}+(n+1)\mathcal{S}_{2n+2}+(n-2)\mathcal{S}_{2n+3}\big{]}\,. (S32)

The last term SAS_{A} is given as,

SA=12𝒮2n+1.\displaystyle S_{A}=\frac{1}{2}\mathcal{S}_{2n+1}\,. (S33)

This implies, we have,

i=12n+1Sai+i=12n+1SaiSA=\displaystyle\sum_{i=1}^{2n+1}S_{a_{i}}^{+}-\sum_{i=1}^{2n+1}S_{a_{i}}^{-}-S_{A}= 12[2𝒮2n+2+(n+1)𝒮2n+3+(n2)𝒮2n+4]\displaystyle\ \frac{1}{2}\big{[}2\mathcal{S}_{2n+2}+(n+1)\mathcal{S}_{2n+3}+(n-2)\mathcal{S}_{2n+4}\big{]} (S34)
12[2𝒮2n+1+(n+1)𝒮2n+2+(n2)𝒮2n+3]12𝒮2n+1,\displaystyle-\frac{1}{2}\big{[}2\mathcal{S}_{2n+1}+(n+1)\mathcal{S}_{2n+2}+(n-2)\mathcal{S}_{2n+3}\big{]}-\frac{1}{2}\mathcal{S}_{2n+1}\,, (S35)
=\displaystyle= log𝒟.\displaystyle\log\mathcal{D}\,. (S36)

where we have plugged in Eq. (7) to obtain the final result.

S3 Proof of holographic entropy inequalities in non-degenerate ground states of TQFT

Now we will introduce the tripartite form of the facet HEIs given in Hernández-Cuenca et al. (2023).

Definition S3.1 (Tripartite form Hernández-Cuenca et al. (2023)).

An information quantity 𝒬\mathcal{Q} is said to be in the tripartite form if it is expressed as

𝒬=iI3(Xi:Yi:Zi|Wi)\mathcal{Q}=\sum_{i}-I_{3}(X_{i}:Y_{i}:Z_{i}|W_{i}) (S37)

where the arguments Xi,Yi,Zi,Wi[N]X_{i},Y_{i},Z_{i},Wi\subset[N] are disjoint subsystems, the sum runs over any finite number of terms, and we allow for the conditioning to trivialize, Wi=W_{i}=\emptyset, in which case I3(Xi:Yi:Zi|)=I3(Xi:Yi:Zi)I_{3}(X_{i}:Y_{i}:Z_{i}|\emptyset)=I_{3}(X_{i}:Y_{i}:Z_{i}) and, they are defined to be

I3(Xi:Yi:Zi|Wi)=I3(Xi:Yi:ZiWi)I3(Xi:Yi:Wi)I_{3}(X_{i}:Y_{i}:Z_{i}|W_{i})=I_{3}(X_{i}:Y_{i}:Z_{i}W_{i})-I_{3}(X_{i}:Y_{i}:W_{i}) (S38)

and,

I3(Xi:Yi:Zi)=Xi+Yi+ZiXiYiXiZiYiZi+XiYiZiI_{3}(X_{i}:Y_{i}:Z_{i})=X_{i}+Y_{i}+Z_{i}-X_{i}Y_{i}-X_{i}Z_{i}-Y_{i}Z_{i}+X_{i}Y_{i}Z_{i} (S39)

We denote IpCqI^{p}C^{q} for a 𝒬\mathcal{Q} that has pp number of I3(Xi:Yi:Zi)-I_{3}(X_{i}:Y_{i}:Z_{i}) and qq number of I3(Xi:Yi:Zi|Wi)-I_{3}(X_{i}:Y_{i}:Z_{i}|W_{i}) terms in the sum (S37).

We borrow the following conjecture from Hernández-Cuenca et al. (2023), which holds true for all known superbalanced facet HEIs.

Conjecture S3.1.

All facet inequalities (except SA) are expressible in the IpCqI^{p}C^{q} form with p1p\geq 1 and q0q\geq 0.

The definition S3.1, taken together with the conjecture S3.1, gives a straightforward proof for proposition III.2 used in the text.

Proposition S3.1.

All facet HEIs of the HEC (n3)(n\geq 3) are satisfied by the quantum entanglement entropy in the unique ground state of a topologically ordered two-dimensional medium with a mass gap.

Proof.

From proposition III.1, we have all facet HEIs are topological. We will now prove that they are also non-negative. As follows from sections II and S1, we have γ=log(𝒟)\gamma=\log{\mathcal{D}} on the 2d2d disk geometry.

The conditional multi-information terms (denoted as CC) in a facet HEI expressed in triparite form (Eq. S37) is a difference of two multi-information, both of which are topological and therefore, cancels each other. The only contribution comes from the multi-information terms (denoted as II). Since, II is topological and I0I\leq 0 (by extension of conjecture III.2), and by conjecture S3.1, p1p\geq 1, thus 𝒬0\mathcal{Q}\geq 0.

Therefore, the facet HEIs are non-negative on a 2d2d disk geometry. ∎