Entanglement contour and modular flow from subset entanglement entropies

The Entanglement contour function quantifies the contribution from each degree of freedom in a region A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{A} $$\end{document} to the entanglement entropy SA\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {S}_{\mathcal{A}} $$\end{document}. Recently in [1] the author gave two proposals for the entanglement contour in two-dimensional theories. The first proposal is a fine structure analysis of the entanglement wedge, which applies to holographic theories. The second proposal is a claim that for general two-dimensional theories the partial entanglement entropy is given by a linear combination of entanglement entropies of relevant subsets inside A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \mathcal{A} $$\end{document}. In this paper, we further study the partial entanglement entropy proposal by showing that it satisfies all the rational requirements proposed previously. We also extend the fine structure analysis from vacuum AdS space to BTZ black holes. Furthermore, we give a simple prescription to generate the local modular flows for two-dimensional theories from only the entanglement entropies without refer to the explicit Rindler transformations.


Introduction
The entanglement entropy for quantum many-body systems is defined as the von Neumann entropy of the reduced density matrix. It has become a quite hot topic in the study of modern theoretical physics. In condense matter theories it can be used to distinguish new topological phases and characterize critical points, e.g., [2][3][4][5][6]. In high energy theory, in the context of AdS/CFT [7][8][9] the Ryu-Takayanagi (RT) formula [10,11] (see also [12][13][14][15] for its extension to holographies beyond AdS/CFT) relates quantum entanglement to spacetime geometry, thus entanglement entropy becomes an important tool to study quantum gravity and holography itself.
However, entanglement entropy only contains part of the information in the reduced density matrix. The authors of [16] considered the possibility of the existence of a function f A (i) which captures how much the degrees of freedom in site i contribute to the entanglement entropy S A of the region A. In other words f A (i), which we call the entanglement contour function following [16], characterizes the density function of the entanglement entropy in A. In the continuum limit, we denote this function as f A (x 1 , · · · x d−1 ), where x i with 1 ≤ i ≤ d − 1, are the coordinates that parameterize A, and d is the spacetime dimension of the field theory. By definition it should satisfies the following basic requirements Also in [16], a set of physical requirements for the entanglement contour is proposed. However these requirements cannot uniquely determine the contour function. The fundamental definition of the entanglement contour based on the reduced density matrix is still not JHEP05(2020)018 general divide A into three subsets, A = α L ∪ α ∪ α R . Here α L and α R represent the subsets on the left and right hand side of α. The subsets α L or α R vanishes if α shares one boundary with A, hence the partition only involves two subsets. Also we can choose α = A hence both of α L and α R vanish. The PEE proposal gives an ansatz for the PEE s A (α), which is a linear combination of subset entanglement entropies. Let, for example, S α∪α L denote the entanglement entropy of the subset α ∪ α L , the proposal claims, Consider the case shown in figure 1, if we chose α = A 2 , then we have α L = A 1 and α R = A 3 . The proposal (2.1) gives the s A (A 2 ) in [1] s We can also apply the proposal (2.1) to A 1 and A 3 . For example, if we consider α = A 1 , then we should have α R = A 2 ∪ A 3 while α L vanishes. So according to (2.1) we should have Similarly we can easily get For consistency, one can easily check that Since the PEE proposal (2.1) is a linear combination of subset entanglement entropies, it will be interesting to explore the properties satisfied by the proposal using the inequalities [30,31] obeyed by entanglement entropy. In the following we list some of the important properties satisfied by s A (α): 1. Additivity: splitting an arbitrary subset A i into arbitrary two adjacent parts The upper bound (2.8) indicates that for any site or subset α inside A, s A (α) should not be larger than the entanglement entropy of the subset. This is reasonable when we interpret s A (α) as the contribution from α to the entanglement entropy S A . Because, unlike S α , s A (α) does not count the entanglement between α and other subsets inside A, thus should be no larger.
6. Symmetry: given region A and a partition into the subsets {A i }, we consider a symmetry transformation T of the theory thus T A = A and T A i = A i . Since T is a symmetry, the subsets A i and A i should play the equivalent role in A. As a result we should have In other words the PEE proposal respect the symmetry T .
For holographic CFTs, the proposal furthermore satisfies a lower bound because the mutual information is monogamous [34], i.e. (2.11) The above equality indicates that where the second inequality comes from the Araki-Lieb triangle inequality We can test its rationality using the Bit thread picture [28], which quantifies the quantum entanglement with a set of "Bit threads" with a cross-sectional area of 4 Planck areas and can only end on the boundary or the horizon. In this picture S A is the maximum possible number of threads emanating from A and end on A c . Also the PEE s A (A 2 ) have a quite natural definition as the number of threads emanating from A 2 and end on A c , when JHEP05(2020)018 the number of threads from A to A c is maximal. So far all the configurations that maximize the flux of bit threads on the RT surface E A are considered to be degenerate in the bit thread picture as they give the same S A . However different degenerate configurations will give different s A (A 2 ) thus the entanglement contour can be different. In the bit thread picture s A (A 2 ) has upper and lower bounds. For example, s A (A 2 ) reaches its upper bound when the number of bit threads emanating from A 2 is maximized and they all cross On the other hand when s A (A 1 ∪ A 3 ) reaches its upper bound s A (A 2 ) will reach its lower bound, which is given by (2.13) The above inequality makes sense when A 1 or A 3 vanishes thus the right and side can be positive. It is interesting that, since S A 2 ≥ S A −S A 1 ∪A 3 , the lower bound (2.12) for s A (A 2 ) defined by (2.1) is stronger than the bound (2.13) we get from the bit thread picture. This is not surprising since the bit threads have the largest freedom to move in the bulk thus should give the weakest bound for the PEE. This implies the freedom for the bit threads may be confined in some way. It is also interesting to compare the PEE proposal with mutual information. Assuming A c is the region that purifies A hence A ∪ A c is in a pure state. When the subset A 2 shares one boundary with A, it is easy to see that the s A (A 2 ) calculated by the proposal is just half of the mutual entanglement entropy between A 2 and A c , (2.14) Similarly, in the case of figure 1, we have For holographic CFT, the mutual information is monogamous, so I(A 1 ∪ A 2 : A c ) ≥ I(A 1 : A c ) + I(A 2 : A c ). Using S A = S A c and S A 1 ∪A 3 = S A 2 ∪A c , we can write the lower bound (2.12) as This means the s A (A 2 ) calculated by the proposal is always larger than half of the mutual information between A 2 and A c . Though the PEE proposal (2.1) satisfies all the physical requirements proposed in [16], it is not justified because those requirements are not enough to uniquely determine the PEE. More generally, assuming that the PEE is a linear combination of all the subset entanglement entropies, 1 i.e.

JHEP05(2020)018
Where the coefficients are constants that are independent from the choice of the partition. Imposing the requirements of additivity (2.6) and normalization, one can uniquely determine the coefficients in the ansatz and get (2.1).
3 Fine structure analysis in BTZ black holes

The fine correspondence in BTZ background
In the context of AdS 3 /CFT 2 , a holographic picture for the entanglement contour of a single interval is given in [1]. When the modular flow is local, the boundary modular flow lines in the causal development D A can be defined as the integral curves along the boundary modular flow. Allowing the points on a modular flow line to flow under the bulk modular flow, the trajectory will form a two-dimensional slice in the entanglement wedge. We call these slices the modular planes. From this construction, we can consider the entanglement wedge as a slicing of the modular planes. Each modular plane will intersect with the interval A on a point A(x m ) and its RT surface E A on another point E A (x m ). It was argued in [1] that the cyclic gluing of the point A(x m ) turns on the nonzero contribution to the S A on E A (x m ), thus gives a fine correspondence between the points on A and E A . We can read the entanglement contour function from this fine correspondence. The discussion in [1] is restricted to the pure AdS 3 space. We stress that the prescription can be applied to more general cases with local modular flow. Here we extend the discussion to BTZ black holes.
The key to conduct the fine structure analysis of the entanglement wedge is the bulk and boundary modular flows. In the case of pure AdS 3 , a Rindler transformation, which is a symmetry transformation, can be constructed to map the causal development D A to a "Rindler space" with infinitely far away boundaries. This Rindler transformation can be extend to the bulk [13] which map the corresponding entanglement wedge to a Rindler AdS 3 space, which is the exterior region of a hyperbolic black hole [17], The horizonr = L of the hyperbolic black hole maps to the RT surface in the oringinal AdS 3 , and the Rindler time translation ∂τ maps to the bulk modular flow in the entanglement wedge. So with the inverse bulk Rinlder transformation we can solve the explicit information for bulk and boundary modular flows, and furthermore construct the modular planes in the entanglement wedge [1]. In principle the above prescription can be extend to the BTZ black holes. In other words we can also try to construct a bulk Rindler transformation map the entanglement wedge in BTZ background to the hyperbolic black hole (3.1). However, the construction of the Rindler transformations will be much more complicated. The explicit bulk modular flow in this background is obtained in [18] using other tricks. Constructing modular planes directly in the entanglement wedge following the bulk and boundary modular flows seems to be tough and unwise. However, it will be much easier to do the fine structure analysis in the Rindler bulk space. The images of the static interval A and its RT surface E A (at JHEP05(2020)018  Our goal is to find the fine correspondence between points onÃ andẼÃ and then map this fine correspondence back to the original BTZ black hole.
The bulk and boundary modular flows are just given by ∂τ . It is easy to see that the hyperbolic black hole can be considered as a slicing of the AdS 2 slices with fixedx. The left figure in figure 2 shows a typical AdS 2 slice where the boundary is settled at thẽ r = ∞. Since the boundary modular flow is also the bulk modular flow, the trajectory of the boundary modular flow line (the blue line) is itself rather than a two-dimensional slice, thus the concept of modular plane seems not to exist. This subtlety also exist when we look at the original AdS space. The key to understand this subtlety is to push the boundary modular flow line in the original AdS 3 from r = ∞ to r = r I (r I is large but still finite) first, then take the limit r I → ∞. We argue that the line at r = r I maps to the blue line in the right figure of figure 2. Because both of the larger and large |τ | region in the Rindler bulk map to the large r region in the original AdS 3 . So the blue line will intersect with the horizon (the red line) at some large |τ |. Also let us look at the original AdS 3 , the red line maps to the null geodesics emanating normally from the RT surface. These null geodesics form the null boundary of the entanglement wedge and approach (at large |τ |) the two tips of the causal development D A at the boundary (at large r).
In this case the blue line have an orbit under the modular flow ∂τ , thus the corresponding modular plane is just the AdS 2 slice. Then we take the limit r I → ∞ and the right figure approaches the left figure. In conclusion the modular plane is defined marginally.
Since the AdS 2 slice of the hyperbolic black hole is the modular plane, according to [1] the points onÃ andẼÃ that correspond to each other stay in the same AdS 2 slice. It is easy to see that each pair of points that correspond to each other can be connected by a geodesic that normal to the Rindler horizon, which is just the black line in figure 2 along the radius direction. This fine correspondence still holds when we map back to the original BTZ spacetime. In other words the pair of points on A and E A that correspond to each other should be connected by a geodesic normal to E A . See figure 3.

JHEP05(2020)018
where z = z h is the horizon. We consider a static interval A with length and its middle point settled at x = 0. In the following we will use the conventions in [36] to parameterize the geodesics. Using (x m , z m ) to parameterize E A , we have where z * is the maximum depth of the RT surface. It is related to by Note that when becomes larger than some critical point, the geodesic that compute the entanglement entropy transfer from the connected geodesic to the two component disconnected geodesics [37]. In this paper, we focus on the cases with connected RT surfaces, and leave the disconnected cases for future discussion. It is not hard to find out the slope s m of the vectors normal to E A at any point (x m , z m ), which is given by The plus sign correspond to x m > 0 while the minus sign correspond to x m < 0. Then we consider a static geodesic γ m that passes through (x m , z m ) with a slope s m . We require the JHEP05(2020)018 slope to satisfy (3.6) thus γ m is normal to E A . According to our previous discussion, the partner point of (x m , z m ) on A is the point (x m , 0) where A intersect with γ m . Following [36] γ m can be written as follows: where the coefficients C i are given by . (3.8) Without loss of generality we focus on the x m > 0 part. We find the geodesic γ m intersect with the boundary at the following two points Only the one with the smaller x coordinate belongs to A. One can check that x 2 < x 1 thus x m = x 2 . Using the definitions (3.8) we havē Furthermore one can check that the above fine correspondence (3.10) is consistent with a much simpler equation It is not easy to get (3.12) by directly simplifying (3.10). We get (3.12) by combining the fine structure analysis and the PEE proposal together, which we will show later.
Consider an infinite small segment dx m on A, its partner under the fine correspondence is an infinite small segment on E A , whose length will give the PEE for dx m . More explicitly the contour function f A (x m ) can be read from the following equation

JHEP05(2020)018
We can solve z m as a function of x m from (3.4) and furthermore plug in the fine correspondence (3.10). At last we find the entanglement contour function is given by where we have replacedx m with x.
On the other hand, since the entanglement entropies for arbitrary single intervals are known in the dual CFT 2 with finite temperature, we can use the PEE proposal (2.1) to derive the PEE and furthermore the entanglement contour function based on the proposal. The entanglement entropy for any static interval with length ∆x is given by where β = 2πz h . Now we consider the following partition of A = A 1 ∪ A 3 with We apply (3.15) to the subsets A 1 and A 3 to calculate their entanglement entropies and use the proposal to calculate the PEE. Then we get The proposal furthermore proposes that the entanglement contour on A should be, Since c = 3 2G and β = 2πz h , the above result exactly matches with our previous contour function (3.14) we got from the fine structure analysis.

Correspondence between bulk geodesic chords and boundary PEE
The fine correspondence between the points on the interval A and its corresponding RT surface E A can be extend to the correspondence between any geodesic chords E i to the PEE of a subset A i , and furthermore to a linear combination of subset entanglement entropies according to our proposal (2.1). For any geodesic chord E i we extend it to a whole geodesic E A anchored on the boundary which correspond to a boundary interval A according to the RT proposal. Then A i is the subset of A that correspond to E i according to the fine correspondence. In other words, we have Also the similar construction is conducted for WCFT in [15] in the context of AdS 3 /WCFT correspondence [38,39].

JHEP05(2020)018
Now we use the PEE proposal to determine the fine correspondence in a much simpler way and derive the simpler relation (3.12). Since all the subset entanglement entropies are known, we can use our proposal (2.1) to calculate s A (A i ) and then impose the matching condition (3.19) for the bulk geodesic chords. The strategy is to consider a pair of points A(x m ) and E A (x m ) on A and E A that lead to the partitions A = A 1 ∪ A 3 , E A = E 1 ∪ E 3 . Then these two points are correspond to each other when they satisfy the matching condition (3.19). The right hand side of (3.19) can be easily calculated by integrating the length of the geodesic chords. Integrating ds along the RT surface from The length of E 1 is infinite, so we use (− 2 + δ x , δ z ) to denote the point where we cut off the RT surface. Since δ z is related to the UV cutoff by δ z = , we have δ x = π β coth π β 2 , such that On the other hand the PEE S A (A 1 ) can be calculated using the PEE proposal. The result is given by (3.17) with the x replaced withx m . Then according to (3.19) we match ( This relation is just (3.12) which coincide with the fine correspondence (3.10) we get from the fine structure analysis in the entanglement wedge.
For an arbitrary static geodesic chord E 2 settled on a static geodesic, which is the RT surface of a boundary interval with length and its middle point settled at x = 0.
Assuming the x coordinate of the end points of E 2 are given by x 2 and x 1 with x 2 > x 1 , then the length of E 2 is given by the PEE s A (A 2 ) with the subsets given by This is a finer version of the RT formula and plays a similar role as the kinematic space formalism [40][41][42][43]. Some explicit comparison between the kinematic space picture and the fine structure analysis picture can be found in [25,44].

JHEP05(2020)018 4 Generating local modular flows from entanglement entropy
The entanglement contour contains much more information than entanglement entropy alone. Then it is interesting to ask whether the entanglement contour proposed by (2.1) contains all the information needed to reconstruct the reduced density matrix (or modular Hamiltonian). In this section we show that, for the cases where the modular Hamiltonian is local thus generate a local modular (geometric) flow, the modular flow can be extracted from the PEE via our proposal. This involves an interesting property of the PEE observed in [1]. More explicitly for an arbitrary subset α of an interval A, we have Here A could be any spacelike interval that has the same causal development as A, i.e. D A = D A . And α is the subset of A which intersect with the same class of modular flow lines (or modular planes) as α. In other words the orbits of two end points of α under the modular flow is a pair of modular flow lines. Any spacelike sub-interval α whose end points respectively settle on this pair of modular flow lines will satisfy (4.1). This property is obvious in the fine structure analysis of the entanglement wedge [1] since the cyclic gluing of both α and α turn on the replica story on the same class of modular planes thus turn on contribution from the same bulk geodesic chord. We see no obstacle to extend this property to covariant intervals and furthermore to other holographic cases with local modular flow. The more interesting question is whether this property can be extend to the non-holographic cases. Previously the only way to solve the local modular flow relies on the Rindler transformations, which are symmetry transformations that map the causal development of an interval to a Rindler space with infinitely far away boundaries. The translation along the Rindler time in the Rindler space then maps to the modular flow in the causal development under the inverse Rindler transformations. The Rindler transformations for static intervals in the vacuum state of CFT 2 and static spheres in the vacuum state of higher dimensional CFTs are firstly constructed in [17]. Then in [13,14] a prescription to construct Rindler transformations in 2-dimensional theories using the global symmetries was proposed. Following this prescription the Rindler transformations and also the modular flows for any intervals (beyond the static case) have been constructed for CFT 2 in the vacuum state [13], warped CFT [13,45] and field theories invariant under the BMS 3 transformations (BMS-FTs) [14]. Though the prescription is quite successful in 2-dimensional theories, it is not easy to carry out. Even for the case of CFT 2 in a thermal state, the construction of the Rindler transformations has not been carried out yet. Nevertheless the modular flow in this case is carried out in [18] using other tricks.
In this section, we propose a much simpler prescription to generate the local modular flow for single intervals in 2-dimensional field theories from entanglement entropy. Provided the entanglement entropy for any intervals are known, the PEE s A (α ) for any sub-interval α in the causal development D A can be calculated via the PEE proposal (2.1). Here we need combine the property (4.1) with the PEE proposal. It will be more convenient to consider the partition A = α ∪ α R with α L vanished thus the left end point of α coincide JHEP05(2020)018 with the left end point of A. We denote the other end point of α as O. We only consider A that shares end points with A. The requirement D A = D A is satisfied by this choice of A . We consider a similar partition A = α ∪ α R , where α also shares the left end point with α. We denote the right end point α as O (see figure 4). Then the modular flow line that passes through O should consist of all the points O that satisfy (4.1). In the following we will explicitly use this strategy to generate the modular flows for intervals in CFT 2 , warped CFT and BMSFT.
Modular flows for CFT 2 in a vacuum state or thermal state. In order to describe the modular flow and thermal circle of the state, we need to embed the theory in a coordinate system. It is convenient to consider the CFT 2 duals to the following BTZ background in light-cone coordinates where u and v are the light-cone coordinates, and T u (T v ) is the left (right) moving temperatures. We consider an arbitrary interval with the following end points When T u = T v = 0 the metric (4.2) describes the Poincaré AdS 3 which duals to the vacuum state. The modular flows in the dual CFT is given by [15] Integration along the modular flow will give us the modular flow lines where c 0 is the integration constant that characterizes different modular flow lines.

JHEP05(2020)018
Then we use our prescription to derive the modular flow lines (4.5). The entanglement entropy for any interval is given by where ∆u and ∆v are the coordinate differences between the two end points and u ( v ) is the cutoff on the left (right) moving sector. Assuming all the O that satisfy (4.1) form a line described by the function (u, v(u)). Then according to the proposal (2.1) we have . (4.7) According to (4.1), the PEE s A (α ) should be a constant along the line (u, v(u)). So we should have which reduces in the following equation The solution is given by where c 1 is the integration constant. Choosing c 1 = c 0 , the above solution is exactly (4.5) thus reproduces the modular flow lines. For the cases with a thermal circle, the left and right moving temperatures are usually nonzero. The modular flow is obtained in [18], which is given by The corresponding modular flow lines are given by where c 0 is again an integration constant. Now we reproduce the modular flow lines (4.12) using the PEE proposal and (4.1). The entanglement entropy for an arbitrary interval is given by The equation ∂ u s A (α ) = 0 then reduces to One can check that the solution of the above equation is exactly (4.12) which describes the modular flow lines.
Modular flows for WCFT. The AdS 3 (4.2) with a Dirichlet-Neumann type of boundary conditions [39] duals to a WCFT [38]. Here we focus on the case with T u = 0, T v = 1. The modular flow in this case was carried out in [15] using the Rindler method, The corresponding modular flow lines are given by where c 0 is an integration constant.
On the other hand, the entanglement entropy for an arbitrary interval in WCFT is given by Then the PEE proposal gives 19) and the condition ∂ u s A (α ) = 0 reduces to The solution of the above condition is exactly the modular flow lines (4.17).
For WCFT in general temperatures the entanglement entropy becomes (4.21) Following our prescription we can easily get the corresponding modular flow lines v(u) = c 0 + log sinh luTu 2 − uT u csch luTu JHEP05(2020)018 Figure 5. An example of a more general partition in higher dimensions. The subset A 1 divide the region A = A 1 ∪ A 2 ∪ A 3 ∪ A 4 ∪ A 5 into five subsets.

Discussion
Since we only used the construction of the proposal (2.1) and the general inequalities satisfied by entanglement entropies, the PEE proposal should satisfy the above 6 properties in general two-dimensional theories. One important underlying ingredient of our proposal is that, no matter how we do the partition the subsets that forms A are always in sequential order. This is crucial for the PEE (2.1) to be additive. However, in general the partitions in higher dimensions are more complicated and there is no natural way to define the order for a general partition, For example, see the case in figure 5. In these cases we need to impose further constrains on the entanglement entropies to make the PEE additive.
Since the entanglement contour function respect the symmetries, for highly symmetric configurations 2 the contour functions may only depend on one coordinate. For example consider A to be a spherical shell in a theory that is rotational symmetric with respect to the center of the shell, the rotation symmetry implies that the contour function only depend on the radius coordinate. Note that the contour function does not depend on how we do the partition, so it is enough to only consider partitions that respect the symmetry, i.e. with all the subsets being concentric spherical shells (see the left figure of figure 6). Since a natural order can be defined along the radius coordinate, the PEE proposal in these configurations should also satisfy all the six requirements [16]. Similarly we can consider configurations with translation symmetries and with the partition shown in the right figure of figure 6. We call these kind of configurations the quasi-one-dimensional configurations. A more explicit discussion on the entanglement contour for quasi-one dimensional cases in higher dimensions can be found in [29].
For holographic theories the fine structure analysis indicates that the PEE is invariant under the modular flow of A in the sense of (4.1). However this is not obvious at all from the field theory side. As we have seen using the PEE proposal and the property (4.1) we can generate the right modular flows in a much easier way. These are highly non-trivial results and gives strong support for the PEE proposal. It will be interesting to test the property (4.1) in more general theories that are not holographic. JHEP05(2020)018 Figure 6. Quasi-one-dimensional configurations with rotation and translation symmetries. Here Combine our proposal with the holographic picture, we get the correspondence between bulk geodesic chords and boundary PEEs, which can be considered as a finer version of the RT formula or its analogues [12][13][14][15] in holographic theories beyond AdS/CFT. One can consult [44] for an interesting application of this correspondence to interpret the bulk volume with boundary entanglement entropies. Since the minimal entanglement wedge cross sections are geodesic chords, the PEE can also be closely related to the quantities that are assumed to be dual to the minimal cross section, which include the entanglement of purification [52,53] and the logarithmic negativity [54].
On the other way around, if we can calculate the entanglement contour using the fine correspondence in holography [1,25] or other methods (for example the Gaussian formula [16,[19][20][21][22]) for lattice models , we can calculate the entanglement entropies of the subsets with less symmetries via the PEE proposal [29].
As we have noted the study of entanglement contour is still at an initial stage. Many foundational questions remain to be answered. For example, in general can the entanglement contour be uniquely determined? How can we generalize the PEE proposal to disconnected regions? Does the entanglement contour exist in cases with non-local modular Hamiltonian? If it exists then how can we extend the above discussions to the non-local cases? 3 We leave all these questions for future investigation.
Note: during the accomplishment of this paper, [25] appears which has some overlap with our paper. The discussion on uniqueness of the entanglement contour and the justification of the PEE proposal in Poincaré invariant field theories can be found in a later paper [56].

JHEP05(2020)018
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.