Characterizing the width of entanglement

The size of controllable quantum systems has grown in recent times. Therefore, the spatial degree of freedom becomes more and more important in experimental quantum systems. However, the investigation of entanglement in many-body systems mainly concentrated on the number of entangled particles and ignored the spatial degree of freedom so far. As a consequence, a general concept together with experimentally realizable criteria has been missing to describe the spatial distribution of entanglement. We close this gap by introducing the concept of entanglement width as measure of the spatial distribution of entanglement in many-body systems. We develop criteria to detect the width of entanglement based solely on global observables. As a result, our entanglement criteria can be applied easily to many-body systems since single-particle addressing is not necessary.


Introduction
To experimentally implement quantum technologies [1,2], such as information processing, simulation [3] or metrology [4], characterizing and understanding multipartite entanglement [5] is important. By characterizing multipartite entanglement one is able to understand experimental setups, identify experimental limitations, and investigate possible noise sources in a better way. The size of controllable quantum systems has grown in recent times and large arrays of atoms [6] or clouds of macroscopic singlet states [7] have been produced. Also ideas of coupling several ion traps to build a quantum computer exist [8]. However, for large quantum systems the spatial degree of freedom becomes more and more important since for large systems external fields cannot be approximated by constant fields anymore. Furthermore, entanglement can be easily protected against constant fields but is very vulnerable to spatially varying fields. As a result, the spatial distribution of entanglement compared to the spatial distribution of external fields is decisive for the temporal evolution of spatially extended quantum systems. This affects also the performance of quantum metrology schemes for gradient fields. For example, the precision of measuring a gradient with the help of pairs of singlet states aligned in a chain of fixed length = · L N d scales as N 2 if arbitrary widths of entanglement are allowed, but it is constant with respect to N if only entanglement of neighboring particles is allowed (see section 2). The spatial distribution of entanglement plays also an important role in the investigation of quantum phase transitions [9][10][11][12][13][14] and makes a distinction between different given ground states of the generalized Heisenberg spin-chain possible [15,16]. Entanglement of multipartite system can be characterized with different quantities, such as the entanglement depth or k-producibility [17,18], which is defined as the minimum number of entangled particles necessary to create a given state. In systems without any spatial ordering, entanglement depth is a powerful variable to characterize the entanglement properties of this system. However, in systems with spatial order such as spin chains, or in the presence of gradient fields, the dynamics of a system may also depend on whether entanglement exists only between neighboring particles or between distant particles.
The investigations done so far concentrated on entanglement depth or k-producibility [17][18][19] or required addressability of single subsystems [9][10][11][12][13]. Our criteria, developed in this paper, are based solely on global observables. Therefore, they open the possibility to study correlation propagation and other physical characteristics of many-body systems without the necessity of addressing single subsystems. To develop criteria for the width of entanglement, we follow the route of detecting multipartite entanglement with the help of squeezing parameters [17,20] or minimal energy [18,19]. However, instead of determining limits depending on the number of entangled particles, we estimate limits for given spatial entanglement distributions.
The paper is organized as follows: first, we introduce the concept of entanglement width to characterize the spatial distribution of entanglement. Then, we give an example of how the width of entanglement influences the time evolution of a quantum system before we present methods to characterize the width of entanglement with the help of global observables. We conclude by demonstrating how quantum phase transitions may manifest themselves in the width of entanglement.

The concept of entanglement width
To define the width of entanglement there has to be some spatial ordering or structure of the involved partitions/particles. This can be e.g. the arrangement of ions in a linear chain (figure 1(a)), cold atoms in a lattice ( figure 1(b)) or a general graph where interactions are only allowed between certain particles (figure 1(c)). The width of entanglement w of a pure state y Yñ = ñ | ⨂ | j j , given by the product of one or several entangled states y ñ | j , is defined as the maximal distance w of two entangled particles within the states y ñ | j (see figure 1). Here, the distance is defined by counting the particles following the given structure starting and ending at two entangled particles. A completely separable state exhibits an entanglement width of w=1. The entanglement width of a mixed state is defined by the minimum with w over all decomposition By definition, the entanglement depth is a lower bound of the entanglement width. However, the entanglement width does not make any statement about the entanglement depth. For example, the width of entanglement w=6 in figure 1 (a) stays the same, no matter if all particles are entangled with each other or only the two outer ones (1 and 6), whereas the entanglement depth changes from k=6 to k=2. Depending on the given structure, we might define the width for several directions as demonstrated in figure 1(b) (total width w x =3 in the horizontal directions and w y =2 in the vertical direction) or along allowed interactions visualized by a graph as demonstrated in figure 1(c) (total width w = 4). Furthermore, states with equal entanglement depth but different entanglement width can lead to dramatic different effects as we demonstrate in the following example.
Motivated by [15] consider a chain of particles j at the positions = + · x x j d j 0 , with equal spacing d between two particles. Furthermore, consider the observable j j x j y j z denoting the Pauli matrices acting on particle j and λ being a parameter determining the observable. Such an observable can be created, e.g., by a time evolution : whereas the entanglement depth disregards any spatial ordering, the definition of entanglement width requires the particles to be spatially ordered, e.g. in a spin chain or a grid. The entanglement depth of the state Yñ | in (a) is given by k=3 (since maximally three particles are entangled). This is a lower bound on the entanglement width in (b), which equals w=6 (since entanglement occurs over a distance of six particles in the chain). In the case of the 2d grid in (c) we distinguish between the width in the horizontal direction (w x =3) and vertical direction (w y =2). In a general graph (d) where edges denote possible interaction between two particles there might exist different ways to connect two particles. In this case, we take always the shortest way to determine the entanglement width.
of spins in a gradient field which rotates the spins depending on their position. Here, λ is given by the gradient of the field. Another approach to create an observable similar to  J is given by a standing light wave coupled to cold atoms in a lattice [15]. In this way, it is possible to distinguish between four given ground states with different spatial entanglement configuration which are important for condensed matter and high energy physics. Whereas Eckart et al needed the previous knowledge that their state is given by exactly one out of four states, the following method to investigate the width of entanglement does not need this strict restriction.  J is given by the global spin operator, if l d is an integer and e.g.
where particles are combined to pairs forming together the singlet state y ñ = ñ -ñ -| (| | ) 01 10 2 . However, for  l Ï d the variance depends on the spatial distribution of the entangled pairs ( j, k ) as can be seen in figure 2 for N=16 and = - . For example for l = ( ) d N 1 2 the 'hugging' configuration with w=N, where particle j is entangled with particle + -N j 1 (solid green line), reaches the minimal variance of However, if all particles with odd number j are entangled with their right neighbor = + k j 1 (dashed blue line), the width of entanglement is given by w=2, which leads to the variance As a consequence, the variance of  J depends crucially on the spatial distribution of entanglement. Also the quantum Fisher information (QFI), which is an important measure for quantum metrology and entanglement [21], is strongly influenced by the spatial distribution of entanglement. For example, the time evolution given in equation (3) leads to a QFI given by = F N hug 2 for the hugging configuration and = · L d N fixed, whereas it is constant in N for the right neighbor configuration. As a consequence, the two here considered states exhibit different behavior although their are equal in their entanglement depth.

First criterion for entanglement width
In this section we demonstrate that the observable D  J , defined in equation (2), is able to distinguish between short-range and long-range entanglement. We will estimate the minimal variance D  ( ) J 2 for states with nearestneighbor entanglement and demonstrate that states possessing long range entanglement are able to violate these bounds. Since D  ( ) J 2 is a concave function, it reaches its minimum for pure states. Furthermore, since nearestneighbor entanglement implies the entanglement of at most two particles, we find for pure states . The product state (non-encircled particles, w = 1) is chosen in such a way, that it minimizes the variance. Although, both entangled states possess the same entanglement depth k=2, they exhibit quite different behavior due to their different spatial distribution of entanglement.
2 with the two-particle variance . The minimum of this two-particle variance is given by (for a proof see appendix A). The remaining task is to optimize over all possible combinations of entangled pairs. This is a classical finite optimization task and for many parameters λ easy to estimate. For example for l = ( ) the optimal pairing is given by 1, and we find e > -2 3 for all pairs of a a , j k . Therefore, the lower bound is exactly given by equation (5). A simple lower bound for the optimal pairing is given by In this way, we can find easily lower bounds not only for nearest-neighbor entanglement but also for general maximal widths w. For example, in figure 3 we show for N=16 the variance D  ( ) J hug 2 which can beat the limits for w=2 and w=4 for l » ( ) d N 7 2 . In general, the criterion introduced here is not directly a criterion of the entanglement width but distinguishes states where entangled particles are equally coupled to the observable (e =  D =  J 1 0 ) from states with non-equal coupling of entangled particles (the smaller ε the larger D  J ). As a consequence, the observable  J defined in equation (2) cannot only be defined with the help of the sine function but also with any for some noneneighboring particles j k , . As a consequence, our criterion detects only states with similar symmetries than the used function f (x). In general, small deviations from this symmetry are also allowed as long as they are small

Second criterion for entanglement width
Another method to distinguish between short-range and long-range entanglement is to use several observables with different correlations. For example the Hamiltonian includes correlations between all spins with equal weight. Both observables can be used to detect multipartite entanglement [18,21,22]. However, they can only together distinguish between nearest-neighbor and long-range entanglement.
The difference between the energy and the total spin may be a good indicator for entanglement beyond nearest-neighbors, which implies an entanglement width of  A quantum state which has obviously entanglement beyond nearest-neighbor entanglement is the state where always two non-nearest-neighbor spins form together the Bell state y ñ = ñ -ñ -| (| | ) 01 10 2 . As a consequence, we find c = - sin . The violation decreases with N 1 . However, note that a maximal violation proportional to N 1 appears also in other entanglement criteria, see e.g. [24]. The reason may be, that our criterion is highly symmetric and the minimal overlap between maximal entangled symmetric states and product states decreases faster than N 1 [25], which is a consequence of the de Finetti theorem.
As an application for this criterion, we consider spin chains with nearest-neighbor and next-nearest This model is often called the J 1 -J 2 -model. It is used to understand phenomena in magnetic materials such as Tomonaga-Luttinger-liquids states and spin-Peierls states [26,27]. Here, the interactions between nearest-and next-nearest neighbors compete with each other, leading to frustrated spins.
To estimate the quality of our criterion, we compare it with a criterion for entanglement depth. As an example we investigate the entanglement of the ground state of a spin chain with N=8 spins and an interaction Hamiltonian given by H 2 defined in equation (14) for different values of α 1 .
In figure 4 we compare the expectation value of the energy H 2 (blue crosses) with the minimal achievable energy for states with entanglement depth  k 1 (blue dashed dotted line) and  k 2 (blue dashed line) as well as the correlation function χ (red * ) and the limit for states with entanglement width  w 2 (red line). As can be seen in figure 4 the ground state exhibits entanglement depth of  k 3 for  a 0.3 and  a 1.1. On the other hand, we detect entanglement width of  w 3 for  a 0.86 using χ. As a consequence, there exist a regime (black square) where we detect already an entanglement width of > w 2 with the help of the correlation function χ, but no multipartite entanglement with the help of global measurements. Furthermore, the correlation function χ shows two sudden jumps at a » 0.5 and a » 0.7 which are indicators that for  ¥ N phase transitions may occur at these points [12]. Contrary to that, we find no hint for phase transitions in the expectation value á ñ H 2 . The behavior found with our entanglement criterion fits well to existing results in the literature. As shown by Majumdar and Ghosh [28,29], the ground state of H 2 is given by for a = 1 2. Therefore, the ground state at this point shows bipartite entanglement between nearest neighbors as also indicated by our results. Furthermore, the relation between entanglement and phase transitions for the Hamiltonian H 2 has been investigated with the help of the concurrence [12] and generalized geometric measures [13]. Similar to our results, Gu et al found with the help of the concurrence that the entanglement changes from nearest neighbor to next-nearest neighbor entanglement around a » 0.75 for N=8 spins [12]. Biswas et al showed that phase transitions which occur only for  ¥ N manifest themselves in the generalized geometric measure for entanglement also for finite N around a » 0.7 [13]. Our entanglement criterion is easily accessible experimentally in contrast to the methods used by Biswas et al and Gu et al, which need the exact knowledge of the state and addressability of single particles.

Conclusion
In summary, we have introduced the concept of entanglement width and demonstrated that states with equal entanglement depth but different entanglement width behave differently. We developed criteria based solely on global observables which are able to distinguish between different values of entanglement width. With the help of these entanglement criteria we investigated the ground state of a spin chain described by H 2 and showed that we are able to detect long-range entanglement before we detected multipartite entanglement.
Furthermore, phase transitions as predicted by other theories for spin-chains described by H 2 manifest themselves in our concept of entanglement width whereas they were not visible in our measure for entanglement depth. Although, we just started to investigate the width of entanglement it turned out to be an important indicator to investigate the behavior of quantum systems in the presence of space dependent interactions. In this way, we could used it as indicator for many-body phenomena such as phase transitions. It would be very interesting to further investigate this possibility.
In addition, detecting the width of entanglement with global observables is not limited to the two methods demonstrated in this paper. Another example is e.g. given by measuring the variance of the total spin with and without a linear dependency on the position and a consecutive comparison. Here, the variance of the total spin without external field gives information about the entanglement, whereas the position dependent observable measures the spatial distribution of correlations, classical as well as quantum (for more details see appendix D). We are certain that many more criteria detecting the width of entanglement will be developed in the future.
Since the two vectors always lie in a plane, we assume s á ñ = 0 y for both particles. This implies phase relations between the parameters a−d which implies that the maximum of g can be reached by choosing all parameters to be real. As a consequence we get e y e e e e ñ = + and the circulant correlation matrix C. As a result, the minimum of c N cl for only nearest-neighbor entangled states is bounded from below by the minimal eigenvalue of C. The correlation matrix C is a circulant matrix [16], and the eigenvalues of C are given by we have found an eigenvector of C with the minimal eigenvalue and the length N.
We note two interesting facts: (i) the eigenvectors of circulant matrices like C are independent of the coefficient c k , only the eigenvalues depend on c k . (ii) Although the vectors  x j live in a three-dimensional space, the vectors  x j forming the eigenvectors are only two-dimensional. Since the eigenvectors are independent of the coefficient c k the scheme presented here can be generalized to all correlation functions which are circulant, that is = -c c j k j k , and the minimum of χ is always given by l · ( ) N min . As a consequence, the scheme presented here can be used to find the minimal energy of arbitrary circulant Hamiltonian.

Appendix C. Minimal energy H 2 for k=2
To derive the minimal energy á ñ H 2 achievable by states with entanglement depth  k 2 we divide the spin chain into blocks, where the ions within a block may be entangled but ions belonging to different blocks must be separable analog to [11,12]. Using the methods of [12], the minimal achievable energy is bounded by where C j describes the optimization over each single block. For single-particle blocks C j is given by For two-particle blocks we obtain with the help of lemma A1 of [12] a a a = - As a consequence, for states with entanglement depth  k 2 we derive the limit  á ñ - Appendix D. Detecting the width of entanglement with the help of a gradient Another way to detect the entanglement width is to change the weight of different interactions through an external space-dependent field, e.g. by measuring the collective spin in the presence of a magnetic gradient. In the following, we assume for simplicity an even number of spins and that the magnetic field is zero in the middle of the chain. Therefore, the interaction of the magnetic field with particle k can be described by As depicted in figure D1. This bound is tight, since it is reached at the point á ñ = á ñ = . Therefore, all states with expectation values á ñ B 2 exceeding this bound must posses entanglement between these two subgroups. If we allow entanglement only between nearest neighbors, then compared to á ñ B 2 cl only the correlation between the two middle particles within the chain are allowed to change from two to six. As a consequence, all states with á ñ > á ñ + B B 4 2 2 cl posses entanglement beyond nearest neighbors.
Furthermore, the behavior of á ñ B 2 for N=4 leads us to the conjecture, that the largest possible value of á ñ B 2 for a given á ñ J 2 with only entanglement within within the subgroups ¼ For N=4 the conjecture coincides with the inequality given in equation (D2). Furthermore, first tests with random states confirm this bound.