Analytic Expression of Quantum Discords in Werner States under LQCC

In this paper, quantum discords in a special kind of states, i.e., Werner states by local quantum operations and classical communication (LQCC) protocols (WLQCC states), are studied. Nineteen parameters to quantify the quantum discords are reduced to four parameters in terms of properties of Werner states and quantum discord. In the case of orthogonal projective measures, analytic expression of quantum discords in WLQCC states is analytically worked out. Some properties of the quantum discord in the WLQCC states are obtained, especially the variation relations between the quantum discords and the parameters characterizing the WLQCC states. By virtue of numerical computations, quantum discords in a Werner state before and after LQCC protocols are compared. It is found that quantum discord in any WLQCC state cannot exceed that in the original Werner state.

Up to now, the form of quantum correlation related to LQCC operations is focused on quantum entanglement. A distinct feature has been found that quantum entanglement cannot be increased through the LQCC channel. It is well known that quantum correlations can be classified into two forms. One form is quantum entanglement, another is quantum correlation different from quantum entanglement (QCDE). The first kind of QCDE, namely quantum discord (QD), was first revealed by Ollivier and Zurek [18] in 2001. Distinctly, it exists even in some separable states, where quantum entanglement does not emerge obviously. This distinct feature attracted many researchers further investigations. Later, other kinds of QCDEs were proposed [19][20][21][22][23][24][25][26][27][28][29][30][31] and their role in quantum information theory is standing out gradually. Moreover, the potential applications of these discovered QCDEs in quantum information processing are gradually exploited, such as quantum computation [32,33], state merging [34], remote state preparation [35], assisted state discrimination [36], quantum correlation swapping [37][38][39], etc. By far, it has been well admitted that QCDE is another form of important quantum resources. Hence, it is interesting to study the behavior of QCDEs in quantum states through LQCC operations. Particularly, similar to the case of quantum entanglement, we wonder whether LQCC can increase QCDEs in quantum states or not.
In this paper, we will study the influences of LQCC on QD in a given quantum state. Considering the difficulties in characterizing the LQCC operation and quantifying QD in a quantum state, we will originally study a comparatively simple case. To be concrete, we will study whether QDs involved in two-qubit Werner states can be increased by the single-state LQCC protocol or not. As for this specific topic, our main motivations are fourfold: (1) study the QD in the concerned states through LQCC operations quantitatively; (2) discover whether LQCC can increase QD or not similar to the quantum entanglement case; (3) if yes or no, seek out the similarities and differences among them qualitatively; (4) try to expose the quantitative features of QDs in the LQCC process. It is worth stressing that so far quantum discord in a generic state is extremely difficult to evaluate. Because of the definition of discord requires the optimization to be executed over all measurements on one subsystem. Until now, the explicit analytical results on QD are known only for a few classes of quantum states [40][41][42][43][44].
The rest of this paper is outlined as follows. In Section 2, WLQCC states are originally introduced and pointedly simplified. In Section 3, the analytical expression of QD in the WLQCC states is derived. In Section 4, some analyses and discussions about the QDs are made and some distinct features are revealed. Finally, a concise summary is given in Section 5.

WLQCC State
Usually, a two-qubit Werner state is written as [45] where I denotes the 4 × 4 identity matrix and |Ψ − = (|01 − |10 )/ √ 2 is the singlet state. The Werner state ρ W ab is characterized by a single real parameter F called fidelity with 1/4 ≤ F ≤. The single-state LQCC protocol maps a Werner state ρ W ab to WLQCC state ab in the form of [46] ab = where A and B are operators act on subsystems a and b, respectively. The success probability of this protocol is Tr(A ⊗ Bρ W ab A † ⊗ B † ). The only restrictions for operators A and B are A † A ≤ I 2 and B † B ≤ I 2 . They can be written as A = U a f a U a and B = V b f b V b , where U a , U a , V b and V b are arbitrary unitary operators, f a and f b are filtering operators taking the forms of f a = ∑ 1 i=0 α i a |i i| and f b = ∑ 1 j=0 α j b |j j|, respectively. Considering physical realizations of filtering operators f a and f b , their eigenvalues must be between zero and one, i.e., 0 ≤ α i a , α j b ≤ 1. Thus any nontrivial LQCC (i.e., any LQCC which is not the zero map) for two-qubit states can be written in the form of [46] where γ is a scale factor in the range 0 < γ ≤ 1, f a = 1 0 0 α a and f b = 1 0 0 α b with 0 ≤ α a , α b ≤ 1. Actually, the scale factor γ can be ignored for that it will vanish in the normalizing process. In spite of the trivial parameter γ, there are still nineteen parameters required to determine a WLQCC state.
According to [47], there exists an equivalence relation V b ρ W ab V † b = V T a ρ W ab V * a by virtue of the inherent permutation symmetry of two-qubit Werner states ρ W ab . Making use of this property, one can further reexpress the state ab as where with W a = U a V T a . Obviously, W a can be seen as an arbitrary unitary operator on subsystem a in the matrix form of here ζ is in the range 0 ≤ ζ ≤ 1 and 0 ≤ µ, ν, θ ≤ 2π. After substituting the matrix forms of ρ W ab , W a , and W † a into Equation (5), one can achieve the matrix form of k ab in computational bases: Furthermore, after the normalizing process, one can get where ξ ab = 6 k ab and Obviously, one can find that the number of parameters required to determine the state k ab is seven, i.e., F, α a , α b , ζ, µ, ν, and θ.

Quantum Discord in the WLQCC State
Quantum discord (QD) is a kind of quantum correlation quantifier proposed in 2001 [18]. In terms of definition, it is easy to prove that QD is invariant under local unitary operations on host qubits. Using this property, according to Equation (4), one can get Hence, to calculate QD one can directly consider the kernel state k ab instead of the state ab .
QD is defined as the discrepancy between two kinds of mutual information expressions which are identical in classical information theory. One kind of mutual information is generalized directly from classical mutual information. It is defined as I(ρ ab ) = S(ρ a ) + S(ρ b ) − S(ρ ab ), where S(·) denotes the von Neumann entropy and ρ a(b) is a marginal state of ρ ab . This kind of mutual information is called quantum mutual information, which is utilized to characterize and quantify the total correlation in state ρ ab . For a WLQCC state ab in Equation (7), its total correlation can be obtained as where and Actually, r k s are the solutions of the quartic equation Three special cases should be clarified here: (1) when F = 1, r 1 = Ω and r 2 = r 3 = r 4 = 0; (2) when α a = 0 and α b = 0,r 1 = Ω and r 2 = r 3 = r 4 = 0; (3) when α a = 1 and α b = 1, Equation (11) is not suitable anymore, and this case will be discussed in the analyses part of this paper. Easily, one can see that I( ab ) is a function of four parameters F, α a , α b and ζ.
Another kind of the mutual information is generalized from the classical analogy which is related to conditional entropy. It is defined as the maximal classical correlation C a (ρ ab ) = S(ρ b ) − min {Ê i a } ∑ i p i a S(ρ b|i ) by optimizing over all possible measurements on part a (or b).
where {Ê i } represents the set of positive-operator-valued measure (POVM) elements such that ∑ iÊ i = I.
To get access to the maximal classical correlation, one must consider all possible POVM for optimization. As a matter of fact, the optimization is quite difficult for a generic state. In this study, von Neumann measurements are employed as an approximation. The measure element performed on subsystem a (or b) of the state ρ ab can be generally expressed asÊ i =ΓP iΓ † , i = 0, 1, whereP i = |i i| is the projector for the subsystem a (or b) along the computational bases |i andΓ ∈ SU(2) is a unitary operator.
1| with parameters τ in the range 0 ≤ τ ≤ 1 and 0 ≤ α, β ≤ π. Note that classical correlation is generally not a symmetric quantity, i.e., C a (ρ ab ) = C b (ρ ab ). C a (ρ ab ) is usually referred to as the "left" classical correlation, while C b (ρ ab ) the 'right' classical correlation.
For our concerned state k ab , it's "left" classical correlation is a function of nine parameters F, α a , α b , ζ, µ, ν, τ, α, and β before optimizing processes: where To get the value of "left" classical correlation for a given Werner derivative (F, α a , α b , ζ, µ, ν, and θ are given), one has to work out the extreme points of ∑ 1 i=0 p i a S( b|i ) at first by solving the following three partial derivative equations: After inspecting the form of ∑ 1 i=0 p i a S( b|i ), we find that parameters α and β induced by measurement are only present in p 0 a in the form of cos(µ + ν − α − β). To simplify partial derivative procedures, we regard the function cos(µ + ν − α − β) as a new parameter X in the scope of [−1, 1]. Thus the latter two equations in Equation (14) are transformed into the following two judgments: Actually, it is still quite difficult to directly solve them. Fortunately, we find that the partial derivative equations hold if τ = 0.5 and X = 0. This implies that the point (τ, X) = (0.5, 0) is an extreme point.
Compare values of this extreme point with other special points, one easily finds that the former value is smaller than the latter ones. This means that this extreme point is at least a local minimal point. We conjecture that it is actually the global minimal point. To verify this conjecture, we have further investigated this feature via numerical calculations and conclusively confirmed that it is exactly the global minimal point. This global minimal point means that to achieve the "left" classical correlation for each state ab , the optimal set of measurement {Ê i a } is determined by parameters τ = 0.5 and α, β, where α, β are related to µ, ν by the equation cos(µ + ν − α − β) = 0.
After substituting τ = 0.5 and X = 0 into Equation (12), we achieve the analytic expression for "left" classical correlation in state ab : where ∆ a = l a + 4α 2 a α 2 b (4F − 1) 2 /Ω 2 . Now let us turn to the "right" classical correlation of ab , it is a function of F, α a , α b , ζ, µ, ν, and θ after optimized on parameters τ, α, and β which induced by the set of measurement elements {Ê i b }: where To get the value of "right" classical correlation for a given Werner derivative, one has to work out the extreme points at first by solving the following three partial derivative equations: Similarly to the solution employed for the optimal procedure in dealing with "left" classical correlation, we regard the function cos(µ − ν + α + β + θ) as a new parameter Y in the scope of [−1, 1]. Thus the latter two equations in Equation (19) are transformed into the following two judgments: With the assistance of numerical calculations, we confirm that the point (τ, Y) = (0.5, 0) is the optimal point for function ∑ 1 i=0 p i b S( k a|i ) for every state k ab . The optimal set of measurement elements {Ê i b } is determined by τ, α, and β, where τ = 0.5 and α, β are determined by the values of µ, ν, and θ.
After substituting values of τ = 0.5 and Y = 0 into Equation (17), we achieve the analytic expression for "right" classical correlation in state k ab : where The quantum discord of k ab is defined as the discrepancy between the total correlation I( k ab ) and the classical correlation C a(b) ( k ab ), i.e., D a(b) ( k ab ) = I( k ab ) − C a(b) ( k ab ). With the analytic expression of "left" and "right" classical correlations, we can achieve the analytical expression for "left" and "right" QD in state ab :

Analyses and Discussions
In the previous section, we have obtained the analytic expressions of the "left" and "right" QDs in the WLQCC states. According to the analytical expression in Equation (25), in this section, we will make some analyses and discussions on these obtained QDs.
Firstly we want to emphasize that, we have developed and used some skills to obtain the analytic expressions of the QD, as is a distinct feature of our present study. In this paper, the concerned state is the WLQCC state ab (see in Equation (2)). Due to the complicated form of LQCC protocol, the number of parameters needed to define the state ab is great. The corresponding number of independent parameters is 19. Hence, it is very difficult to consider the state ab directly. Fortunately, after some analyses, we find one can use some skills to degrade this difficulty. To be concrete, noticing that our main purpose in this paper is to calculate the QD in the WLQCC state ab . Moreover, ab is related to k ab (see in Equation (5)) with local unitary transformations and QD is invariant under the transformations. Hence, in calculating QD in ab one can consider the kernel state k ab instead. The parameters needed in k ab is 7. Hence, the quantity and complexity in calculating QD are reduced to a great extent. See Table 1. Furthermore, through complicated calculation and derivation, one can find that the final number of the independent parameters defining QD in the concerned state k ab is 4. That is to say, the final analytical expression obtained from the QD of the WLQCC states is only a function of 4 independent parameters, i.e., F, α a , α b , and ζ. Using the above skills and derivations, the difficulty is evidently degraded and the analytic expression is obtained.
The second, let us move to investigate the features in the obtained QD. According to the symmetry of the 'left' quantum discord and the "right" one, in the following, we will make some analyses about the variation of the "left" quantum discord D a with F, α a , α b , and ζ. Through analyses, one can get the following distinct features.
(i) D a is an increasing function of F for any given set of (α a , α b , ζ) and reaches its maxima at F = 1 (See Figure 1). This phenomenon can be understood as follows. Inspect the original Werner state in Equation (1) and the WLQCC state in Equation (2). Originally, state ρ W ab consists of two terms, which are mixed with the weights characterized by F. The former represented by I is the maximally mixed state with no quantum correlation. The latter is the maximally entangled state. In this case, a bigger F corresponds to a larger weight of the latter term and accordingly a larger quantum correlation is induced by the mixture of the two terms. As for the WLQCC state ab which are mapped from the original state ρ W ab , quantum correlation in it still increases with the parameter F. This means that the WLQCC operation do not change the dependency relation about QC with F. However, according to [40], one can find that for a given F the quantity of QC in ab is less than that in ρ W ab . This means that the nonzero set of (α a , α b , ζ) can decrease the quantum discord for a given F. That is to say, the LQCC protocol cannot increase the quantum discord in the Werner state ρ W ab . (ii) For given ζ ∈ [0, 0.5], F, and α b , D a monotonically increases with α a and reaches its maxima at α a = 1 (See Figure 2). However, for given ζ ∈ [0.5, 1.0] , F, and α b , D a first increases then decreases with α a , and the transition point is determined by the exact value of ζ. As for the case of α b , one can find that the dependent relation of D a on α b for given ζ, F and α a is similar to that on α a . This feature is related to the LQCC operation and its underlying physics can be explained as follows. The LQCC operation is concerned with two parties, Alice and Bob, who each control one subsystem, and who are restricted to carrying out local quantum operations and classical communication (i.e., LQCC). Specifically, the quantum operations Alice and Bob are allowed to perform are local unitary transformations and local filtrations. As analyzed in the second section in this paper, to calculate the quantum discord in the WLQCC state ab in Equation (2), we have considered the kernel state k ab in Equation (5) instead. From Equation (5), one can see that α a is the parameter characterizing the filtration operator f a , α b is the one for f b , and ζ is the one for the unitary operator W a . (iii) Through analysis, one can find that quantum discord D a is a monotonically increasing function of ζ for given α a , α b , and F. See Figure 3 for example. ζ is the parameter used to characterize the unitary operation W a . In addition, according to the expression of quantum discords in WLQCC states, the parameter ζ vanishes when α a = 1 or α b = 1. Especially, in the situation of α a = α b = 1, LQCC is turned into local unitary operation, which does not change the quantum discords in Werner states. When α a = 0 or α b = 0, the total correlations of WLQCC states are zero. The special case α a = α b = 0 is meaningless for the matrix form of WLQCC states is zero.
Finally, we want to stress that, although the result, i.e., LQCC cannot increase the QD (similar to the behavior of quantum entanglement), has been obtained, it is still an open question whether this conclusion is applicable for other states or in the framework of other QCDEs. In fact, it is easy to verify that this conclusion is not general. In other words, for other general states and for other quantum correlation measures, the property of LQCC cannot increase and the quantum correlation may not exist at all. Nonetheless, for other states with high inherent property of symmetry, such as the higher-dimensional (qutrit) Werner states, we conjecture that the coincidence may remain. We will pay attention to it in the near future.

Summary
To summarize, in this paper, we have derived the analytical expression of QDs in WLOCC states. With the assistance of numerical computations, we find that the QD in a WLQCC state cannot exceed that of the original Werner state. The research in [PRL 81, 3279] exhibits a similar result in the scenario of quantum entanglement: QDs in two-qubit Werner states cannot be increased by single-state LQCC protocols.