Pentapartite Entanglement Measures of GHZ and W-Class State in the Noninertial Frame.

We study both pentapartite GHZ and W-class states in the noninertial frame and explore their entanglement properties by carrying out the negativities including 1-4, 2-3, and 1-1 tangles, the whole entanglement measures such as algebraic and geometric averages π5 and Π5, and von Neumann entropy. We illustrate graphically the difference between the pentapartite GHZ and W-class states. We find that all 1-4, 2-3 tangles and the whole entanglements, which are observer dependent, degrade more quickly as the number of accelerated qubits increases. The entanglements of these quantities still exist even at the infinite acceleration limit. We also notice that all 1-1 tangles of pentapartite GHZ state Nαβ=NαIβ=NαIβI=0 where α,β∈(A,B,C,D,E), whereas all 1-1 tangles of the W-class state Nαβ,NαIβ and NαIβI are unequal to zero, e.g., Nαβ=0.12111 but NαIβ and NαIβI disappear at r>0.61548 and r>0.38671, respectively. We notice that the entanglement of the pentapartite GHZ and W-class quantum systems decays faster as the number of accelerated particles increases. Moreover, we also illustrate the difference of von Neumann entropy between them and find that the entropy in the pentapartite W-class state is greater than that of GHZ state. The von Neumann entropy in the pentapartite case is more unstable than those of tripartite and tetrapartite subsystems in the noninertial frame.

1. Introduction

The transfer of quantum states between distant nodes of a quantum network is a basic task for quantum information processing. It is well known that all protocols used for quantum state transmission require entanglement between the sender and the receiver systems. Entanglement, which is at the basis of quantum mechanics and almost every quantum information protocol, has become a very interesting topic, particularly in many-body systems, with the recent development of quantum information technology. The correct understanding of entanglement is of importance due to its special application in many branches such as quantum teleportation, quantum communication and quantum cryptography [1,2,3,4,5,6,7,8,9,10], and quantum algorithms [11,12]. Moreover, it is also helpful in studying quantum communication protocols like quantum key distribution (QKD) [13].

Up to now, the development of quantum computing has required the study of multi-qubit entangled systems, so the entanglement properties of multipartite quantum systems under inertial frames are becoming more and more important. However, with the development of relativistic quantum information science, many authors have paid more and more attention to the development of this field. In order to study the property of quantum entangled state existing in the noninertial frame, we have to employ a relativistic setting [14,15,16]. The relativistic quantum information regarded as a new interesting field has emerged for many years since the relation between quantum information science and relativity theory intermediated by quantum field theory allows us to have a deeper understanding of the fundamental principles in quantum mechanics [17,18]. This also helps us explore how the degree of entanglement would be affected in curved space time, particularly by the acceleration parameter. Until now, quantum information theory has been enriched by the contributions of quantum entanglement made in the noninertial frame [19,20,21]. The properties of multipartite entangled systems are mainly related to the acceleration parameter and the number of particles in the noninertial system, which undoubtedly affect the entanglement degree of the entangled system.

In recent years, many relevant and significant contributions to this field have been made [1,3,14,15,16,22,23,24,25,26,27,28,29,30,31,32,33,34,35]. For example, since several pure multipartite entangled systems was studied [25], where the Unruh effect was discussed [19,20,21], the entanglement has been verified as an observer dependent in the noninertial frame. Compared with the well-known entangled stated-GHZ state [15,36,37,38,39,40,41], the authors paid less attention to the W-class state because its density matrix cannot be written as an X matrix form. Nevertheless, we have employed a special technique to study the density matrix in the non-X matrix form and carried out the tripartite and tetrapartite W-class state cases [42,43,44], except for the generalized GHZ state in the noninertial frame [45,46,47,48]. After studying, we find that the degree of entanglement of the W-class state is more robust than the GHZ and generalized GHZ states. It should be recognized that the entanglement for tripartite and tetrapartite systems still exists even at the infinite acceleration limit. However, Ye and her collaborators carried out the enhancement of multipartite entanglement in an open system in the noninertial frame [49]. Other relevant studies such as the fidelity loss and coherence loss, etc., in the open systems have also been done in Refs. [8,10].

Generally speaking, two main entanglement measures, which are named as negativity and von Neumann entropy, have been used to quantify the genuine entanglement. The negativity is employed to study the multi-tangle and the whole entanglement measures. The von Neumann entropy is concerned with the stability of the entangled system. For example, the three tangle, i.e., 1-2 tangle [50], was proposed to study the entanglement property of the entangled quantum system when tripartite Alice, Bob, and Charlie initially shared an arbitrary fermionic three-qubit pure state as well as the later proposed -tangle [36]. The three tangle has interesting properties, but its analytical calculation becomes a nontrivial task because its calculation depends on the negativities of bipartite and tripartite systems. Similarly, four tangle has been proposed when we studied the tetrapartite systems, which include 1-3 and 2-2 tangles [42,44,45,47,48]. Recently, stimulated by the study of the tripartite and tetrapartite W-class state [42,43,44], Sun and her coauthors studied the entanglement property of a pentapartite W-class state in the noninertial frame and showed how the acceleration parameter and the number of the accelerated qubits affect the entanglement property of the pentapartite W-class entangled system [51]. This was realized by studying the -tangle, including the 1-4, 1-1 tangles and the von Neumann entropy. However, they were not concerned with the 2-3 tangle case due to its complication, which is also an important factor to describe the entanglement property of the entangled system. Moreover, as the W-class and GHZ states are the two most important pure states in quantum information, it is necessary to study their entanglement properties simultaneously and show their difference graphically. To enrich the paper [51], we shall present all 1-4, 2-3 tangles and von Neumann entropy for these two important pure states for completeness, which is the main purpose of this work.

This paper is organized as follows. In Section 2, we briefly review the transformation between Minkowski space and Rindler coordinates. How to construct the density matrix for the simplest case is presented. In Section 3, we study the negativities, including 1-4, 2-3, and 1-1 tangles, whole entanglement measures and and von Neumann entropy, which are illustrated graphically to show their difference. Finally, in Section 4 we summarize our conclusions.

2. Pentapartite Entanglement from One to Five Accelerated Observers

The pentapartite GHZ state that we are going to study in this work is given by where so does the state , while the pentapartite W-class state has the following form [52]:

Here we use the subscripts A, B, C, D, and E to denote Alice, Bob, Charlie, David, and Elly (as we know, Eve is not of confidence), respectively. They initially share a pentapartite GHZ or W-class state in the inertial frame. In this work, we assume that the acceleration of particles always starts from the rightmost one in A(Alice), B(Bob), C(Charlie), D(David), E(Elly) qubits. That is to say, we first suppose that Elly is accelerated in a uniform acceleration but Alice, Bob, Charlie, and David remain stationary, and so on.

For entangled GHZ and W-Class states in the noninertial frame, let us use Rindler coordinates to describe a family of observers with a uniform acceleration and divide Minkowski space-time into two inaccessible regions I and II. The rightward accelerating observers are located in region I and causally disconnected from the analogous counterparts in region II [53,54]. Let us briefly review the connection between the vacuum and excitation states in Minkowski coordinates and those in Rindler coordinates. First, let Alice stay stationary, while Bob moves in a uniform acceleration. We consider Bob accelerated uniformly in the plane. Rindler coordinates are appropriate for describing the viewpoint of an observer moving in a uniform acceleration. Two different sets of Rindler coordinates, which differ from each other by an overall change in sign, are necessary for covering Minkowski space. These sets of coordinates define two Rindler regions disconnected from each other, as shown in Figure 1 [16,55]:

Figure 1

Rindler space time diagram: lines of constant position are hyperbolas and lines of constant proper time for the accelerated observer run through the origin. In present work, we denote regions I and II as Bob and anti-Bob, respectively. The reader can refer to Ref. [55] for more information.

A free Dirac field in dimensional Minkowski space satisfies the Dirac equation , where m is the particle mass and the Dirac gamma matrices. A spinor wave function composed of the complete orthogonal set of fermion and antifermion modes can be expressed as , where and are the creation and annihilation operators for fermions (antifermions) of the momentum k, respectively, satisfying the relation . The quantum field theory for a Rindler observer is constructed by expanding the spinor field in light of a complete set of fermion and antifermion modes in regions I and II as follows:

In a similar way, and are the creation and annihilation operators for fermion (antifermions), respectively, acting on region I (II) for and also satisfy a similar anticommutation relation. The relation between creation and annihilation operators in Minkowski and Rindler space times can be found by the Bogoliubov transformation where with and r is a Bob’s acceleration parameter with the range for . It is known from this equation and its adjoint that Bogoliubov transformation mixes a fermion in region I and antifermions in region II. As a result, it is assumed that the Minkowski particle vacuum state for mode k based on Rindler Fock states is given by where the Rindler region I or II Fock states carry a subscript I and II, respectively, on the kets, but the Minkowski Fock states are indicated by the subscript M on the kets. As what follows, we are only interested in using single mode approximation [15,16,24,56,57,58,59], i.e., and also uniform acceleration ( is considered to relate Minkowski and Unruh modes) for simplicity.

Using the single mode approximation, one can transform Bob’s vacuum state and one-particle state in Minkowski space into Rindler space. Using the creation and annihilation operators on Equation (6) above and using the normalization condition, we can obtain [15,16,24,56,57,58,59] where and () are the mode decomposition of into two causally disconnected regions I and II in Rindler space. It should be pointed out that Bruschi et al. discussed the Unruh effect beyond the single mode approximation [21], in which two complex numbers and (the subindexes L and R corresponding to the Left and Right regions in Rindler diagram, i.e., regions I and II) are used to construct the one-particle state, i.e., . However, in the present case for single mode approximation, one has to satisfy the normalization condition . It is also worth noting that a Minkowski mode that defines the Minkowski vacuum is related to a highly nonmonochromatic Rindler mode rather than a single mode with the same frequency (see Refs. [21,30,60,61] for details). Other relevant contributions [31,59,62,63,64,65] have also been made.

To illustrate how to expand in Rindler coordinates, we are going to give explicit expression when Elly is accelerated, i.e.,

Similarly, we can also obtain the expressions of other cases when the observers Alice, Bob, Charlie, and David are accelerated. Such a procedure also works for the pentapartite W-class case.

After the transformation to the Rindler space, we have to trace out the part of the antiparticle state in region from the density matrix . In this case, when Elly is accelerated, the corresponding density matrix is thus given by

For simplicity, we write out explicitly all nonzero elements for pentapartite GHZ and W-class states in Appendix A. These results will be helpful in calculating the negativity and von Neumann entropy, as shown below.

3. Entanglement Measures: Negativity and von Neumann Entropy

3.1. Negativity

Negativity, which is used to measure the entanglement of multipartite systems, is defined by [66,67,68] where , and represent 1-4, 2-3, and 1-1 tangles, respectively. The expressions , and are the trace norms of the partial transposes of the density matrices. Generally speaking, the trace of any Hermitian operator A is equal to the sum of its eigenvalues [69], , i.e., where represents the negative eigenvalue of the matrix M. It should be pointed out the calculation of these negativities is very complicated and time consuming.

Let us first calculate the negativity 1-4 tangle for the GHZ and W-class states when 1 to 5 observer(s) is (are) accelerated. The explicit expressions of the pentapartite GHZ and W-class states are written out in Appendix B for completeness. (It should be pointed out that the special symbols such as , #, and & appeared in this Appendix B and also in Appendix C are generated systematically by Wolfram Mathematica.) In Figure 2, we plot the negativity 1-4 tangle of pentapartite GHZ and W-class states when only one of five observers is accelerated. It is seen in Figure 2a that the decreases from 1 to (29.29% loss of entanglement), whereas the decreases from 1 to (50% loss of entanglement). This means that entanglement is observer dependent. Furthermore, in the case of the W-class state as shown in Figure 2b, the decreases from 0.8 to 0.7048 (11.891% loss of entanglement), but decreases from 0.8 to 0.29282 (63.4% loss of entanglement).

Figure 2

Panels (a,b) show the variation of 1-4 tangle with the parameter r in the case of pentapartite GHZ and W-class states, respectively, when only one qubit is accelerated.

In Figure 3, we plot the negativity 1-4 tangle when two observers are accelerated. As shown in Figure 3a, the negativity 1-4 tangle in the pentapartite GHZ system decreases from 1 to (50% loss of entanglement), whereas decreases to (60.97% loss of entanglement). In the case of the W-class state as displayed in Figure 3b, the negativity decreases from to 0.5924 (25.94% loss of entanglement), whereas the decreases from 0.8 to 0.24515 (81.86% loss of entanglement). We may conclude that the negativity in the GHZ state decays faster than that of W-class state, but in the GHZ state decays slower than that of W-class state.

Figure 3

Same as Figure 2 but when two qubits are accelerated.

In Figure 4, we plot the negativity 1-4 tangle when three observers are accelerated. We can see in Figure 4a that negativity in the GHZ pentapartite system decreases from 1 to (64.65% loss of entanglement), and decreases to (70.35% loss of entanglement). However, in the case of the W-class state as shown in Figure 4b, decreases from to 0.4529 (43.38% loss of entanglement), whereas decreases to (75.41% loss of entanglement).

Figure 4

Same as Figure 2 and Figure 3 but when three qubits are accelerated.

In Figure 5, we plot the negativity 1-4 tangle when four observers are accelerated. We can see in Figure 5a that negativity in the GHZ pentapartite state decreases from 1 to (77.94% loss of entanglement), and decreases also from 1 to (77.94% loss of entanglement). It should be emphasized that the negativities and at both and are the same, but their explicit expressions given in Appendix B are not the same. Moreover, for the W-class state as shown in Figure 5b, the decreases from 0.8 to 0.1870 (76.62% loss of entanglement), but decreases to a smaller value 0.1436 (82.04% loss of entanglement). We find that the difference of the loss of entanglement for both and of two different pentapartite states is very small.

Figure 5

Same as above but when four qubits are accelerated.

Finally, we find that negativity in the case of GHZ state decreases from 1 to (85.45% loss of entanglement) as seen in Figure 6a, whereas as shown in Figure 6b in the W-class state decreases from 0.8 to 0.0596 (92.54% loss of entanglement). This means that the entanglement of these two pentapartite quantum systems decay most when all particles are accelerated.

Figure 6

Same as above but when all qubits are accelerated.

In Figure 7a,b and Figure 8a,b, we show how the entanglement changes with the number of the accelerated qubits if we only refer to Alice and Elly. It is found that the negativity decreases as the number of accelerated qubits increases, but negativity in the W-class state decreases faster than that in the GHZ state when Elly is taken as a reference (see Figure 8).

Figure 7

Panels (a,b) corresponding to GHZ and W-class states with respect to Alice show the variations of the 1-4 tangle for 1 to 5 arbitrary selected qubits as a function of the acceleration parameter r.

Figure 8

Same as Figure 7 but with respect to Elly.

To calculate the whole entanglement measures such as algebraic and geometric averages and , we have to find the 1-1 tangle of both GHZ and W-class states. In the case of the GHZ state, all 1-1 tangles are equal to zero. In the case of the W-class state, however, some of them that are unequal to zero are expressed as where and . The , , and represent the bipartite subsystems with 0 to 2 accelerated qubits. As shown in Figure 9, it is interesting to see that the entanglement in the 1-1 tangle vanishes at (only one accelerated particle), but vanishes at (two accelerated particles), except for a constant .

Figure 9

Plot of 1-1 tangle for pentapartite W-class state as a function of acceleration parameter r.

We are now in the position to study negativity 2-3 tangle even though it is not required to calculate the whole entanglement measures. However, we want to present them for completeness, as these results have never been presented to our best knowledge. The analytical expressions of both GHZ and W-class states are given in Appendix C. It is found that the negativity of GHZ state at is 1, whereas in W-class state it is . As the acceleration parameter r increases, they all decrease with it, but in W-class state decreases faster than that in the GHZ state. Negativity 2-3 tangle is also dependent on the number of accelerated observers. When only one qubit is accelerated in the GHZ pentapartite state, as displayed in Figure 10a, we notice that are equal to each other, but and are not the same, as shown in Figure 10b. At the infinite acceleration limit, the 2-3 tangle for GHZ state decreases from 1 to ( loss of entanglement) at , but 2-3 tangle and in the case of W-class state, as shown in Figure 10b, decrease from 0.979796 to ( loss of entanglement) and ( loss of entanglement), respectively.

Figure 10

Panels (a,b) show the 2-3 tangle for both GHZ and W-class states, respectively, when only one qubit is accelerated.

When two accelerated observers are considered (see Figure 11), it is found that decreases faster than that of in the case of GHZ state. In the W-class case, the difference between and is almost equal to zero, but the difference between and others such as and is very big. At the infinite acceleration limit, in the case of GHZ state the 2-3 tangles and are equal to 0.5 and 0.3903, respectively, but and in the case of W-class state are equal to 0.6159 and 0.2274, respectively. This implies that the negativity in the W-class state is bigger than that of GHZ state, whereas in the W-class state is smaller than that of GHZ state in the infinite acceleration limit.

Figure 11

Panels (a,b) show the 2-3 tangle for both GHZ and W-class states, respectively, when two qubits are accelerated.

When three observers are accelerated (see Figure 12), we notice that the difference among , and is very small in the GHZ state, as shown in Figure 12a, but in the case of W-class state case, as illustrated in Figure 12b, i.e., their difference is big. In the case of GHZ state when , the 2-3 tangle is 0.353553, but is equal to . However, in the W-class state case, as seen in Figure 12b, we find that their difference is obvious, that is, , and , respectively.

Figure 12

Panels (a,b) show the 2-3 tangle for both GHZ and W-class states, respectively, when three observers are accelerated.

When four qubits are accelerated, in the case of GHZ state, as displayed in Figure 13a, we find that the difference between and is very small, but in the case of W-class state their difference is a little big. At the infinite limit, the 2-3 tangle is equal to in the case of GHZ state, whereas in the W-class state, the 2-3 tangles and , respectively.

Figure 13

Panels (a,b) show the 2-3 tangle for both GHZ and W-class states, respectively, when four qubits are accelerated.

Finally, let us consider the case when all the observers are accelerated simultaneously. We see that all 2-3 tangles of either GHZ state or W-class state are equal to each other, as shown in Figure 14. At , the 2-3 tangles of GHZ and W-class states are equal to and approximately , respectively. The variations of the 2-3 tangles for 1 to 5 arbitrary selected accelerated qubits are displayed in Figure 15. It is found that they all decrease with both the increasing acceleration parameter r and the number of accelerated qubits.

Figure 14

Panels (a,b) show the 2-3 tangle for both GHZ and W-class states, respectively, when all qubits are accelerated.

Figure 15

Panels (a,b) show the 2-3 tangles for both GHZ and W-class states, respectively, when 1 to 5 qubits is (are) accelerated.

3.2. Whole Entanglement Measures

Now, we use the algebraic average -tangle to describe the multipartite entanglement defined by [50,70] from which we are able to calculate the whole residual entanglement -tangle defined by . Moreover, we may use another whole residual entanglement measure named as geometric average [71].

Let us first calculate the whole residual entanglement measures and of both GHZ and W-class states and then plot them. In Figure 16a, it is found that the algebraic average in the GHZ state decreases from 1 to , , , , and for 1 to 5 arbitrary selected qubits, respectively. In Figure 16b, we show the whole residual entanglement in the W-class state. It is interesting to note that decreases from an initial value to ( loss of entanglement), ( loss of entanglement), ( loss of entanglement), ( loss of entanglement), and ( loss of entanglement), respectively.

Figure 16

Panels (a,b) show the whole residual entanglement measure of GHZ and W-class states, respectively, when 1 to 5 observers is (are) accelerated.

Finally, let us show the variation of the geometric average when 1 to 5 accelerated qubits is (are) considered. It is found that is very similar to , i.e., whole residual entanglement and are almost same, as shown in Figure 17. To see this clearly, as shown in Figure 18a,b, there is no difference between residual entanglement and in the GHZ state, but there is a very slight difference in the W-class state.

Figure 17

Same as Figure 16 but for the whole entanglement measures .

Figure 18

Panels (a,b) show the difference between whole entanglement measure when 3 observers are accelerated for the GHZ and W-class states, respectively.

3.3. Entropy

Another useful quantity to measure entanglement is the von Neumann entropy defined by [72,73,74]: where is ith eigenvalue of density matrix . Unlike the negativity, the von Neumann entropy is not required to find the partial transpose of the density matrix except for applying the partial trace to obtain the density matrix of subsystems. We write out explicitly nonzero eigenvalues of GHZ state to calculate entropies in Table 1 but do not list those of W-class state because of complicated expressions.

Table 1

Eigenvalues of GHZ density matrices in the noninertial frame.

Density Matrix	Eigenvalues
ρ_ABCDE_I	λ_31=14(3+cos(2r)) λ_32=sin^2(r)2
ρ_ABCD_IE_I	λ_29=116(11+4cos(2r)+cos(4r)) λ_30=sin^4(r)2 λ_31=λ_32=sin^2(2r)8
ρ_ABC_ID_IE_I	λ_25=164(15cos(2r)+6cos(4r)+cos(6r)+42) λ_26=λ_27=λ_28=λ_29=λ_30=λ_31=12sin^4(r)cos^2(r) λ_32=sin^6(r)2
ρ_AB_IC_ID_IE_I	λ_17=1256(56cos(2r)+28cos(4r)+8cos(6r)+cos(8r)+163) λ_18=λ_19=λ_20=λ_21=12sin^2(r)cos^6(r) λ_22=λ_23=λ_24=λ_25=12sin^6(r)cos^2(r) λ_26=sin^8(r)2 λ_27=λ_28=λ_29=λ_31=λ_32=132sin^4(2r)
ρ_A_IB_IC_ID_IE_I	λ_1=λ_2=11024(382+120cos(4r)+10cos(8r)+                  −2(cos(4r)+7)^2(511cos(4r)+62cos(8r)+cos(12r)+1474)) λ_3=λ_4=λ_5=λ_6=λ_7=λ_8=λ_9=λ_10=λ_11=λ_12=λ_13=λ_14=λ_15=λ_16=λ_17=12sin^4(r)cos^6(r) λ_18=λ_19=λ_20=λ_21=λ_22=λ_23=λ_24=λ_25=λ_26=λ_27=12sin^6(r)cos^4(r) λ_28=λ_29=λ_30=λ_31=λ_32=12sin^8(r)cos^2(r)

Let us show the difference of the von Neumann entropy between GHZ and W-class states. As shown in Figure 19, entropy increases as the number of accelerated observers increases. This means that the system becomes more unstable. In the same condition, it is found that von Neumann entropy of the W-class state increases faster than that of GHZ.

Figure 19

Panels (a,b) show the von Neumann entropy of the GHZ and W-class states when 1, 2, 3, 4, and all observers are accelerated.

4. Conclusions

In this work, we have studied the entanglement measures of pentapartite GHZ and W-class states by investigating the negativity and whole residual entanglement. We have carried out the cases when the 1, 2, 3, 4, or even all observers are accelerated. As we can see, the degree of entanglement will be degraded when the acceleration parameter r increases. However, we have verified again the fact that the degree of entanglement is dependent of the number of the accelerated particles. That is to say, the degree of the entanglement will decrease faster with the number of accelerated particles. The 1-4 and 2-3 tangles for both GHZ and W-class pentapartite states still exist even if the acceleration tends to infinity. Compared between GHZ and W-class states, we find that the degree of the entanglement of W-class state decreases faster than that of GHZ state when the accelerated parameter increases. The maximum values of the GHZ and W-class states are 1 and 0.8 without the acceleration. For 2-3 tangle case, we notice that when only one particle is accelerated, the GHZ entanglement decays faster than the W-class state, but with the increase of accelerated particles, the W-class entanglement decays faster than the GHZ state. However, we also note that the 2-3 tangle corresponding to the cases , and are always the smallest compared to other cases of the same type if we assume that the acceleration of particles always starts from the rightmost two in A(Alice), B(Bob), C(Charlie), D(David), E(Elly) qubits. The whole entanglement measurements show us that entanglement in GHZ state is greater than that of W-class state. However, we find that there is almost no difference between whole residual entanglements and . As far as the von Neumann entropy, compared with the tripartite and tetrapartite entangled systems, the von Neumann entropy of pentapartite system is larger than those of tripartite and tetrapartite cases. As the number of accelerated particles increases, the von Neumann entropy of the system increases accordingly. This implies that the system becomes more and more unstable with the increasing accelerated particles. Before ending this work, we give a useful remark on the difference of the negativity among the tripartite [43], tetrapartite [44], and present pentapartite cases. In the GHZ state case, all 1-1 tangles for them are equal to zero. For the W-class case, however, there only exists a common 1-1 tangle among them. Let us show their difference. For example, the of pentapartite, tetrapartite, and tripartite cases is equal to 0.12111, 0.2071, and 0.412023, respectively, at , but of pentapartite and tetrapartite cases will disappear () at and , respectively, and in the tripartite case at . The of pentapartite, tetrapartite, and tripartite cases is equal to 0.38671, 0.472473, and 0.699185, respectively. This implies that the degree of the entanglement decays faster with the increasing entangled particles so that the system becomes more and more unstable.

Table A1

Nonzero entries for GHZ density matrices.

Density Matrix	Nonzero Entries
ρ_ABCDE_I	[1,1]=[2,2]=12cos^2(r) [32,1]=[1,32]=12cos(r) [32,32]=12
ρ_ABCD_IE_I	[1,1]=12cos^4(r) [32,1]=[1,32]=12cos^2(r) [2,2]=[3,3]=18sin^2(2r) [4,4]=12sin^4(r) [32,32]=12
ρ_ABC_ID_IE_I	[1,1]=12cos^6(r) [1,32]=[32,1]=12cos^3(r) [2,2]=[3,3]=[5,5]=12cos^4(r)sin^2(r) [4,4]=[6,6]=[7,7]=12cos^2(r)sin^4(r) [8,8]=12sin^6(r) [32,32]=12
ρ_AB_IC_ID_IE_I	[1,1]=12cos^8(r) [1,32]=[32,1]=12cos^4(r) [2,2]=[3,3]=[5,5]=[9,9]=12cos^6(r)sin^2(r) [4,4]=[6,6]=[7,7]=[10,10]=[11,11]=[13,13]=132sin^4(2r) [8,8]=[12,12]=[14,14]=[15,15]=12cos^2(r)sin^6(r) [16,16]=12sin^8(r) [32,32]=12
ρ_A_IB_IC_ID_IE_I	[1,1]=12cos^10(r) [1,32]=[32,1]=12cos^5(r) [2,2]=[3,3]=[5,5]=[9,9]=[17,17]=12cos^8(r)sin^2(r) [13,13]=[4,4]=[19,19]=[10,10]=[11,11]=[6,6]=[18,18]=[21,21]=[25,25]=[7,7]=12cos^6(r)sin^4(r) [8,8]=[12,12]=[14,14]=[15,15]=[20,20]=[22,22]=[23,23]=[26,26]=[27,27]=[29,29]=12cos^4(r)sin^6(r) [16,16]=[24,24]=[28,28]=[30,30]=[31,31]=12cos^2(r)sin^8(r) [32,32]=12+12sin^10(r)

Table A2

Nonzero entries for pentapartite W-class state.

Density Matrix	Nonzero Entries
ρ_ABCDE_1	[2,2]=15 [2,3]=[2,5]=[2,9]=[2,17]=[3,2]=[5,2]=[9,2]=[17,2]=15cos(r) [3,3]=[3,5]=[3,9]=[3,17]=[5,3]=[5,5]=[5,9]=[5,17]=[9,3]=[9,5]=[9,9]=[9,17]=[17,3]=[17,5]=[17,9]=[17,17]=15cos^2(r) [4,4],[4,6]=[4,10]=[4,18]=[6,4]=[6,6]=[6,10]=[6,18]=[10,4]=[10,6]=[10,10]=[10,18]=[18,4]=[18,6]=[18,10]=[18,18] =15sin^2(r)
ρ_ABCD_1E_1	[2,2]=[2,3]=[3,2]=[3,3]=15cos^2(r) [2,5]=[2,9]=[2,17]=[3,5]=[3,9]=[3,17]=[5,2]=[5,3]=[9,2]=[9,3]=[17,2]=[17,3]=15cos^3(r) [4,4]=25sin^2(r) [4,6]=[4,7]=[4,10]=[4,11]=[4,18]=[4,19]=[6,4]=[7,4]=[10,4]=[11,4]=[18,4]=[19,4]=15sin^2(r)cos(r) [5,5]=[5,9]=[5,17]=[9,5]=[9,9]=[9,17]=[17,5]=[17,9]=[17,17]=15cos^4(r) [6,6],[6,10]=[6,18]=[7,7]=[7,11]=[7,19]=[10,6]=[10,10]=[10,18]=[11,7]=[11,11]=[11,19]=[18,6]=[18,10]=[18,18]=[19,7]= [19,11]=[19,19]=120sin^2(2r) [8,8]=[8,12]=[8,20]=[12,8]=[12,12]=[12,20]=[20,8]=[20,12]=[20,20]=15sin^4(r)
ρ_ABC_1D_1E_1	[2,2]=[2,3]=[2,5]=[3,2]=[3,3]=[3,5]=[5,2]=[5,3]=[5,5]=15cos^4(r) [2,9]=[2,17]=[3,9]=[3,17]=[5,9]=[5,17]=[9,2]=[9,3]=[9,5]=[17,2]=[17,3]=[17,5]=15cos^5(r) [4,4]=[6,6]=[7,7]=110sin^2(2r) [4,6]=[4,7]=[6,4]=[6,7]=[7,4]=[7,6]=120sin^2(2r) [4,10]=[4,11]=[4,18]=[4,19]=[6,10]=[6,13]=[6,18]=[6,21]=[7,11]=[7,13]=[7,19]=[7,21]=[10,4]=[10,6]=[11,4]=[11,7]= [13,6]=[13,7]=[18,4]=[18,6]=[19,4]=[19,7]=[21,6]=[21,7]=15sin^2(r)cos^3(r) [8,8]=35sin^4(r) [8,12],[8,14]=[8,15]=[8,20]=[8,22]=[8,23]=[12,8]=[14,8]=[15,8]=[20,8]=[22,8]=[23,8]=15sin^4(r)cos(r) [9,9]=[9,17]=[17,9]=[17,17]=15cos^6(r) [10,10],[10,18]=[11,11]=[11,19]=[13,13]=[13,21]=[18,10]=[18,18]=[19,11]=[19,19]=[21,13]=[21,21]=15sin^2(r)cos^4(r) [12,12]=[12,20]=[14,14]=[14,22]=[15,15]=[15,23]=[20,12]=[20,20]=[22,14]=[22,22]=[23,15]=[23,23]=1/5cos(r)^2sin(r)^4 [16,16]=[16,24]=[24,16]=[24,24]=15sin^6(r)
ρ_AB_1C_1D_1E_1	[2,2]=[2,3]=[2,5]=[2,9]=[3,2]=[3,3]=[3,5]=[3,9]=[5,2]=[5,3]=[5,5]=[5,9]=[9,2]=[9,3]=[9,5]=[9,9]=15cos^6(r) [2,17]=[3,17]=[5,17]=[9,17]=[17,2]=[17,3]=[17,5]=[17,9]=15cos^7(r) [4,4]=[6,6]=[7,7]=[10,10]=[11,11]=[13,13]=25sin^2(r)cos^4(r) [4,6]=[4,7]=[4,10]=[4,11]=[6,4]=[6,7]=[6,10]=[6,13]=[7,4]=[7,6]=[7,11]=[7,13]=[10,4]=[10,6]=[10,11]=[10,13]=[11,4] =[11,7]=[11,10]=[11,13]=[13,6]=[13,7]=[13,10]=[13,11]=15sin^2(r)cos^4(r) [4,18]=[4,19]=[6,18]=[6,21]=[7,19]=[7,21]=[10,18]=[10,25]=[11,19]=[11,25]=[13,21]=[13,25]=[18,4]=[18,6]=[18,10]= [19,4]=[19,7]=[19,11]=[21,6]=[21,7]=[21,13]=[25,10]=[25,11]=[25,13]=15sin^2(r)cos^5(r) [8,8]=[12,12]=[14,14]=[15,15]=35sin^4(r)cos^2(r) [8,12]=[8,14]=[8,15]=[12,8]=[12,14]=[12,15]=[14,8]=[14,12]=[14,15]=[15,8]=[15,12]=[15,14]=15sin^4(r)cos^2(r) [8,20]=[8,22]=[8,23]=[12,20]=[12,26]=[12,27]=[14,22]=[14,26]=[14,29]=[15,23]=[15,27]=[15,29]=[20,8]=[20,12]=[22,8]= [22,14]=[23,8]=[23,15]=[26,12]=[26,14]=[27,12]=[27,15]=[29,14]=[29,15]=15sin^4(r)cos^3(r) [16,16]=45sin^6(r) [16,24]=[16,28]=[16,30]=[16,31]=[24,16]=[28,16]=[30,16]=[31,16]=15sin^6(r)cos(r) [17,17]=15cos^8(r) [18,18]=[19,19]=[21,21]=[25,25]=15sin^2(r)cos^6(r) [20,20]=[22,22]=[23,23]=[26,26]=[27,27]=[29,29]=180sin^4(2r) [24,24]=[28,28]=[30,30]=[31,31]=15sin^6(r)cos^2(r) [32,32]=15sin^8(r)
ρ_A_1B_1C_1D_1E_1	[2,2]=[2,3]=[2,5]=[2,9]=[2,17]=[3,2]=[3,3]=[3,5]=[3,9]=[3,17]=[5,2]=[5,3]=[5,5]=[5,9]=[5,17]=[9,2]=[9,3]=[9,5]= [9,9]=[9,17]=[17,2]=[17,3]=[17,5]=[17,9]=[17,17]=15cos^8(r) [4,4]=[6,6]=[7,7]=[10,10]=[11,11]=[13,13]=[18,18]=[19,19]=[21,21]=[25,25]=25sin^2(r)cos^6(r) [4,6]=[4,7]=[4,10]=[4,11]=[4,18]=[4,19]=[6,4]=[6,7]=[6,10]=[6,13]=[6,18]=[6,21]=[7,4]=[7,6]=[7,11]=[7,13]=[7,19]= [7,21]=[10,4]=[10,6]=[10,11]=[10,13]=[10,18]=[10,25]=[11,4]=[11,7]=[11,10]=[11,13]=[11,19]=[11,25]=[13,6]=[13,7]= [13,10]=[13,11]=[13,21]=[13,25]=[18,4]=[18,6]=[18,10]=[18,19]=[18,21]=[18,25]=[19,4]=[19,7]=[19,11]=[19,18]=[19,21]= [19,25]=[21,6]=[21,7]=[21,13]=[21,18]=[21,19]=[21,25]=[25,10]=[25,11]=[25,13]=[25,18]=[25,19]=[25,21]=15sin^2(r)cos^6(r) [8,8],[12,12]=[14,14]=[15,15]=[20,20]=[22,22]=[23,23]=[26,26]=[27,27]=[29,29]=380sin^4(2r) [8,12]=[8,14]=[8,15]=[8,20]=[8,22]=[8,23]=[12,8]=[12,14]=[12,15]=[12,20]=[12,26]=[12,27]=[14,8]=[14,12]=[14,15]= [14,22]=[14,26]=[14,29]=[15,8]=[15,12]=[15,14]=[15,23]=[15,27]=[15,29]=[20,8]=[20,12]=[20,22]=[20,23]=[20,26]=[20,27] =[22,8]=[22,14]=[22,20]=[22,23]=[22,26]=[22,29]=[23,8]=[23,15]=[23,20]=[23,22]=[23,27]=[23,29]=[26,12]=[26,14]= [26,20]=[26,22]=[26,27]=[26,29]=[27,12]=[27,15]=[27,20]=[27,23]=[27,26]=[27,29]=[29,14]=[29,15]=[29,22]=[29,23]= [29,26]=[29,27]=180sin^4(2r) [16,16]=[24,24]=[28,28]=[30,30]=[31,31]=45sin^6(r)cos^2(r) [16,24]=[16,28]=[16,30]=[16,31]=[24,16]=[24,28]=[24,30]=[24,31]=[28,16]=[28,24]=[28,30]=[28,31]=[30,16]=[30,24]= [30,28]=[30,31]=[31,16]=[31,24]=[31,28]=[31,30]=15sin^6(r)cos^2(r) [32,32]=sin^8(r)