A quantitative witness for Greenberger-Horne-Zeilinger entanglement.

Along with the vast progress in experimental quantum technologies there is an increasing demand for the quantification of entanglement between three or more quantum systems. Theory still does not provide adequate tools for this purpose. The objective is, besides the quest for exact results, to develop operational methods that allow for efficient entanglement quantification. Here we put forward an analytical approach that serves both these goals. We provide a simple procedure to quantify Greenberger-Horne-Zeilinger-type multipartite entanglement in arbitrary three-qubit states. For two qubits this method is equivalent to Wootters' seminal result for the concurrence. It establishes a close link between entanglement quantification and entanglement detection by witnesses, and can be generalised both to higher dimensions and to more than three parties.

It is a fundamental strength of physics as a science that most of its basic concepts have quantifiability built into their definition. Just think of, e.g., length, time, or electrical current. Their quantifiability allows to measure and compare them in different contexts, and to build mathematical theories with them1. There is no doubt that entanglement is a key concept in quantum theory, but it seems to resist in a wondrous way that universal principle of quantification. The reason for this is, in the first place, that entanglement comes in many different disguises related to its resource character, i.e., what one would like to do with it. In principle, there are numerous task-specific entanglement measures23. However, most of them cannot be calculated easily (nor measured or estimated) for generic mixed quantum states, and therefore it is difficult to use them.

There are notable exceptions, the concurrence4 and the negativity for bipartite systems5. These measures have already provided deep insight into the nature of entanglement, but they also have their shortcomings. The concurrence is strictly applicable only to two-qubit systems while for the negativities it is not known how to distinguish entanglement classes. The generalisations of the concurrence (such as the residual tangle6) do quantify task-specific entanglement even for multipartite systems but again it is not known how to estimate them for general mixed quantum states.

There is another difficulty. An N-qubit density matrix ρ is characterised mathematically by 22 − 1 real parameters. Reducing it to its so-called normal form7—which contains the essential entanglement information—removes 6N parameters. The entanglement measure is determined by the remaining exponentially many parameters which need to be processed to calculate the precise value. Even an operational method similar to that of Wootters-Uhlmann89 would quickly reach its limits with increasing N. Therefore it is desirable to develop methods which provide useful approximate answers even for larger systems. If one asks for mere entanglement detection, witnesses10 are such a tool because here the number of required parameters (both for measurement and processing) can be reduced substantially. There are also estimates of entanglement measures using witness operators1112 which, however, have not yet produced practical methods for entanglement quantification.

Here we develop an easy-to-handle quantitative witness for Greenberger-Horne-Zeilinger (GHZ) entanglement13 in arbitrary three-qubit states. It yields the exact three-tangle for the family of GHZ-symmetric states14, and those states which are locally equivalent to them. For all other states, the method gives an optimised lower bound to the three-tangle. Due to this feature we call the approach a witness.

We start by defining the GHZ symmetry14 and stating our central result. Then we prove the validity of the statement for two qubits. We obtain a method equivalent to that of Wootters-Uhlmann, i.e., it gives the exact concurrence for arbitrary density matrices. Subsequently we explain the extension of the approach to arbitrary three-qubit states.

Results

The procedure

The N-qubit GHZ state in the computational basis is defined as . It is invariant under: (i) Qubit permutation. (ii) Simultaneous spin flips i.e., application of . (iii) Correlated local z rotations: where σ, σ, σ are Pauli matrices. An N-qubit state is called GHZ symmetric and denoted by ρS if it remains invariant under the operations (i)–(iii). An arbitrary N-qubit state ρ can be symmetrized by the operation where the integral denotes averaging over the GHZ symmetry group including permutations and spin flips. Notably, the GHZ-symmetric N-qubit states form a convex subset of the space of all N-qubit states.

Observation:

If an appropriate entanglement measure µ is known exactly for GHZ-symmetric N-qubit states ρS, it can be employed to quantify GHZ-type entanglement in arbitrary N-qubit states ρ. Here, µ(ψ) is a positive -invariant function of homogeneous degree 2 in the coefficients of a pure quantum state ψ, and µ(ρ) is its convex-roof extension15. The estimate for µ(ρ) is found in the following sequence of steps:

(1) Given a state ρ, derive a normal form ρNF(ρ), i.e., apply local filtering operations so that all local density matrices are proportional to the identity7 (see Section Methods). If ρNF(ρ) = 0 the procedure terminates here, and µ(ρ) = 0.

(2) Renormalise ρNF/tr ρNF and transform it using local unitaries to obtain the state according to appropriate criteria (see below) so that the entanglement of ρS(ρNF/tr ρNF) is enhanced.

(3) Project the state onto the GHZ-symmetric states . The estimate for µ(ρ) is obtained after renormalisation

Two qubits

For two qubits the entanglement measure under consideration is the concurrence C(ρ) (Refs. 4,8). From the symmetrization ρS(ρ) of an arbitrary two-qubit state ρ we find (for details see Supplementary Information): In the symmetrization entanglement may be lost, as illustrated by the state for which inequality (3) gives the poor estimate C(Ψ−) ≥ 0. Therefore, the optimisation steps (1) and (2) are necessary to avoid unwanted entanglement loss in the symmetrization (3). The goal is to augment the right-hand side of inequality (3) up to the point that equality is reached. We will show now that for two qubits this can indeed be achieved.

It is fundamental that the maximum of an -invariant function µ(ρ) under general local operations can be reached by applying the optimal transformation where A = A1 A and is an invertible local operation7. Consider first the normal form ρNF(ρ) which is obtained from ρ by iterating determinant-one local operations7 (see also Methods). Such operations (represented by matrices) describe stochastic local operations and classical communication (SLOCC). Consequently, the normal form is locally equivalent to the original state ρ, that is, it lies in the entanglement class of ρ. Note that the iteration leading to the normal form minimises the trace of the state. Subsequent renormalisation increases the absolute values of all matrix elements in equation (3). Here, the correct rescaling of the mixed-state entanglement measure is crucial. This is why homogeneity degree 2 of µ(ψ) is required1617.

Hence, transforming ρ to its normal form increases the moduli of ρ00,00, ρ00,11, ρ11,00, ρ11,11 (and also the concurrence) as much as possible for a state that is SLOCC equivalent with ρ. The sum of the off-diagonal matrix elements in equation (3) reaches its maximum if ρ00,11 is real and positive. As this can be achieved by a z rotation on one qubit we may consider it part of finding the normal form and drop the absolute value bars in equation (3). Then, the sum of matrix elements equals, up to a factor 1/2, the fidelity of ρNF/tr ρNF with the Bell state . The question is how large this fidelity may become.

To find the answer we transform ρNF/tr ρNF to a Bell-diagonal form using local unitaries (this is always possible71819). If then we apply another operation to maximise (see Supplementary Information). The result is a Bell-diagonal with maximum real off-diagonal element (please note that denotes a normalised state, whereas ρNF is not normalised). However, Bell-diagonal two-qubit density matrices with this property can be made GHZ symmetric without losing entanglement4 (see also Supplementary Information).

Hence, our optimised symmetrization procedure (1)–(3) leads to the exact concurrence for arbitrary two-qubit states ρ. In passing, we have demonstrated that the concurrence is related via C(ρ) = max(0, 2f − 1) · tr ρNF to the maximum fidelity that can be achieved by applying invertible local operations to ρ.

Three qubits

For three qubits, the GHZ-symmetric states are described by two parameters14 and therefore form a two-dimensional submanifold in the space of all three-qubit density matrices. It turns out that it has the shape of a flat isosceles triangle, see Fig. 1. A convenient parametrisation is as it makes the Hilbert-Schmidt metric in the space of density matrices conincide with the Euclidean metric. This way geometrical intuition can be applied to understand the properties of this set of states. All entanglement-related properties of GHZ-symmetric states are symmetric under sign change x −x as this is achieved by applying σ to one of the qubits.

Figure 1

The triangle of GHZ-symmetric three-qubit states.

The upper corners correspond to and the lower corner to ρS(001), cf. Ref. 14. The grey area shows GHZ-class states (τ3 > 0) whereas the yellow area comprises states with vanishing τ3 (“W”). The border between GHZ-class and W-class states is the GHZ/W line, equation (10) (red solid line). We also show a state ρS(x0, y0) together with the point (, ) that is required to determine the three-tangle τ3(x0, y0), equation (6).

The GHZ-class entanglement of three-qubit states is quantified by the three-tangle τ3 (Refs. 6,17, see also Methods). For GHZ-symmetric three-qubit states ρS(x0, y0) the exact solution for the three-tangle20 (see also Methods) is where x0 ≥ 0 and (, ) are the coordinates of the intersection of the GHZ/W line with the direction that contains both GHZ+ and ρS(x0, y0) (cf. Fig. 1). The grey surfaces in Fig. 2 illustrate this solution.

Figure 2

Illustration of the procedure for finding the three-tangle of a general mixed three-qubit state ρ.

In the xy plane, there is the triangle of GHZ-symmetric states while on the vertical axis, the three-tangle for each GHZ-symmetric state (cf. equation (6)) is shown. Simple projection generates a non-optimal GHZ-symmetric state. The optimisation steps (1), (2) move the symmetrization image to with enhanced three-tangle.

Now we turn to constructing a quantitative witness for the three-tangle of arbitrary three-qubit states by using the solution in equation (6). As before, the main idea is that an arbitrary state can be symmetrized according to equation (2) and thus is projected into the GHZ-symmetric states. Again, we assume ρ000,111 real and nonnegative, so that x(ρ) ≥ 0. From Figs. 1 and 2 it appears evident that the entanglement of the symmetrization image ρS(ρ) can be improved by moving its point (x(ρ), y(ρ)) closer to GHZ+. More precisely, the entanglement measure is enhanced upon increasing one of the coordinates without decreasing the other (cf. equations (3) and (6)).

In this spirit, finding the normal form in step (1) is appropriate as it yields the largest possible three-tangle for a state ρNF/tr ρNF locally equivalent to the original ρ (cf. Ref. 7). As the normal form is unique only up to local unitaries it does not automatically give the state with minimum entanglement loss in the symmetrization. Therefore, the unitary optimisation step (2) is required to generate the best coordinates.

In the symmetrization the information contained in various matrix elements is lost. For two qubits, however, the concurrence of the optimised Bell-diagonal states depends only on and the loss of in the symmetrization does not harm. In contrast, the three-qubit normal form depends on 45 parameters. We may not expect that τ3(ρ) depends only on two of them and, hence, entanglement loss in the symmetrization (3) is inevitable (cf. Supplementary Information). Consequently, steps (1)–(3) lead to a lower bound for the three-tangle that coincides with the exact τ3(ρ) at least for those states which are locally equivalent to a GHZ-symmetric state. The most straightforward optimisation criterion in step (2) is to maximise . Alternative criteria which generally do not give the best τ3(ρ) but can be handled more easily (possibly analytically) are maximum fidelity , minimum Hilbert-Schmidt distance of from GHZ+, or maximum Re .

Discussion

Evidently this approach can be generalised. Therefore we conclude with a discussion of some of its universal features. The essential ingredients are an exact solution of the entanglement measure for a sufficiently general family of states with suitable symmetry, and the entanglement optimisation for a given arbitrary state ρ via general local operations. The former determines the border where the entanglement vanishes. The latter ensures an appropriate fidelity of the image ρS(ρ) with the maximally entangled state. This reveals a remarkable relation between entanglement quantification through SL(2, ) invariants and the standard entanglement witnesses which we briefly explain in the following.

A well-known witness for two-qubit entanglement is . It detects the entanglement of an arbitrary normalised two-qubit state ρ2qb if On the other hand, from our concurrence result we see, by dropping the optimisation over SLOCC A = A1 A2, that is a (non-optimised) quantitative witness for two-qubit entanglement. In other words, yields one of the many possible lower bounds to the exact result. Analogously it is straightforward to establish the relation between the standard GHZ witness and the non-optimal quantitative witness . The latter represents a linear lower bound to the three-tangle obtained via the optimisation steps (1)–(3) (see Supplementary Information).

Finally we mention that our approach can be used without optimisation, i.e., either without step (1), or (2), or both. This renders the witness less reliable but more efficient. At best it requires only four matrix elements (for any N). We note that, if we apply the witness to a tomography outcome the measurement effort can be reduced by using the prior knowledge of the state and choosing the local measurement directions such that the fidelity with the expected GHZ state is measured directly. This implements optimisation step (2) right in the measurement.

Methods

Normal form of an N-qubit state

The normal form of a multipartite quantum state is a fundamental concept that was introduced by Verstraete et al.7. It applies to arbitrary (finite-dimensional) multi-qudit states. Here we focus on N-qubit states only.

In the normal form of an N-qubit state ρ, all local density matrices are proportional to the identity. Therefore the normal form is unique up to local unitaries. Remarkably, the normal form can be obtained by applying an appropriate local filtering operation where . Therefore ρNF is locally equivalent to the original state ρ. The normal form ρNF is peculiar since it has the minimal norm of all states in the orbit of ρ generated by local filtering operations. Practically, the normal form can be found by a simple iteration procedure described in Ref. 7. It is worth noticing that GHZ-symmetric states – which play a central role in our discussion – are naturally given in their normal form.

Three-tangle of three-qubit GHZ-symmetric states

The pure-state entanglement monotone that needs to be considered for three-qubit states is the three-tangle τ3(ψ), i.e., the square root of the residual tangle introduced by Coffman et al.6: Here ψ with are the components of a pure three-qubit state in the computational basis. The three-tangle becomes an entanglement measure also for mixed states via the convex-roof extension15 i.e., the minimum average three-tangle taken over all possible pure-state decompositions {p, ψ}. In general it is difficult to carry out the minimisation procedure in equation (9), but there exist various approaches for special families of states161720212223. For GHZ-symmetric three-qubit states, the convex roof of the three-tangle can be calculated exactly (see equation (6)). This solution is shown in Fig. 2 and can be understood as follows. The border between the W and the GHZ states is the GHZ/W line which has the parametrised form14 with −1 ≤ v ≤ 1. The solution for the convex roof is obtained by connecting each point of the GHZ/W line (x, y, τ3 = 0) with the closest of the points (, , τ3 = 1). That is, the three-tangle is nothing but a linear interpolation between the points of the border between GHZ and W states, and the maximally entangled states GHZ±.

Author Contributions

The authors contributed equally to this work.