Numerical and analytical results for geometric measure of coherence and geometric measure of entanglement

Quantifying coherence and entanglement is extremely important in quantum information processing. Here, we present numerical and analytical results for the geometric measure of coherence, and also present numerical results for the geometric measure of entanglement. On the one hand, we first provide a semidefinite algorithm to numerically calculate geometric measure of coherence for arbitrary finite-dimensional mixed states. Based on this semidefinite algorithm, we test randomly generated single-qubit states, single-qutrit states, and a special kind of d-dimensional mixed states. Moreover, we also obtain an analytical solution of geometric measure of coherence for a special kind of mixed states. On the other hand, another algorithm is proposed to calculate the geometric measure of entanglement for arbitrary two-qubit and qubit-qutrit states, and some special kinds of higher dimensional mixed states. For other states, the algorithm can get a lower bound of the geometric measure of entanglement. Randomly generated two-qubit states, the isotropic states and the Werner states are tested. Furthermore, we compare our numerical results with some analytical results, which coincide with each other.

Quantum coherence and entanglement are two basic concepts in quantum information theory, which are extensively applied to quantum information processing and quantum computational tasks 1 . Moreover, both quantum coherence and entanglement can be regarded as quantum resources, and they are useful for quantum-enhanced metrology, quantum key distribution and so on [2][3][4][5][6][7][8] . Therefore, characterizing and quantifying coherence and entanglement become significant parts in quantum information theory 9 .
Quantum coherence is defined for a single system, and is widely used in quantum optics in previous studies [10][11][12][13][14][15][16][17] . For any distance measure D between two arbitrary quantum states, a general coherence measure is defined as C D (ρ) = min δ∈I D(ρ, δ) , i.e., the minimum distance from ρ to all possible incoherent states δ ∈ I , where I is the set of all incoherent states [18][19][20][21][22] . From this definition, one can see that C D (ρ) = 0 if and only if ρ is an incoherent state. Many distance-based coherence measures are proposed, such as geometric measure of coherence, relative entropy of coherence and l p norm of coherence. The geometric measure of coherence is defined by using the fidelity between the measured state ρ and its nearest incoherent state 23 . The relative entropy of coherence is another distance-based coherence measure 18,24,25 . Considering the coherence measures based on the matrix norms, the l 1 norm of coherence was introduced and studied in Ref. 18 . Besides, other different coherence measures have also been proposed 25-28 . Furthermore, many experimental results on coherence have been reported [29][30][31][32][33] .
Quantum entanglement is widely regarded as an essential feature of quantum mechanics, and entanglement measures have many applications [34][35][36][37] . A class of entanglement measures are based on the fact that the closer a state is to the set S of separable states, the less entanglement it has 19,22 . According to the distance measure D between quantum states ρ and σ , it is defined as E D (ρ) = min σ ∈S D(ρ, σ ) , i.e., the measure is the minimum distance to all possible separable states [18][19][20][21][22] . One fundamental distance-based entanglement measure is the relative entropy of entanglement 22 , which can be considered as a strong upper bound for entanglement of distillation 38 . Another one is the geometric measure of entanglement (GME) 9,39,40 . Furthermore, the expected Scientific RepoRtS | (2020) 10:12122 | https://doi.org/10.1038/s41598-020-68979-z www.nature.com/scientificreports/ value of entanglement witnesses can be used to estimate the GME [41][42][43] . Other different entanglement measures have been proposed for multipartite systems and mixed states 19,20 . For most quantum states, the analytical solutions of the coherence and entanglement measures are not available, so numerical algorithms must be applied. In some entanglement measures, several numerical algorithms have been used to solve related problems [44][45][46][47] . Moreover, computing many entanglement measures is NP hard for a general state 48,49 , so some upper and lower bounds are proposed to describe entanglement [50][51][52][53][54][55][56][57] and coherence measures 8,[58][59][60][61] . In Refs. 62,63 , a semidefinite program (SDP) was proposed to calculate the fidelity between two states, which inspire us to apply this semidefinite program to numerically obtain the geometric measures of coherence and entanglement. In the following, we will try to provide semidefinite programs, in order to get the numerical results of the geometric measures of coherence and entanglement.
In this work, we first review the definition and properties of fidelity and its semidefinite program. Then we present numerical and analytical results for the geometric measure of coherence. Our algorithm can be used to numerically obtain the geometric measure of coherence for arbitrary finite-dimensional states. We test our semidefinite program for single-qubit states, single-qutrit states, and a special kind of d-dimensional mixed states. For the special kind of d-dimensional mixed states, we also obtain an analytical solution of its geometric measure of coherence. Furthermore, we also propose another algorithm for the geometric measure of entanglement, which can obtain the geometric measure of entanglement for arbitrary two-qubit and qubit-qutrit states, and some special kinds of higher dimensional mixed states. For other states, a lower bound of the geometric measure of entanglement can be acquired by using the algorithm.

Results
A semidefinite program for computing fidelity. We first review the fidelity and its semidefinite program in the following. The fidelity between states ρ and χ is defined as 1 For a pure state |ψ� and an arbitrary state χ , one can get that In Refs. 62,63 , Watrous and colleagues proposed a semidefinite program, whose optimal value equals the fidelity for given positive semidefinite operators, i.e., considering the following optimization problem where L is the collection of all linear mappings. In a X complex Hilbert space, Pos(X) is the set of positive semidefinite operators operating on X . Then the maximum value of 1 2 Tr(X) + 1 2 Tr(X † ) is equal to F(ρ, χ) . X is a randomly generated complex matrix of the same order as χ.
The SDP can not only solve the problem effectively, but also prove the global optimality under weak conditions 64 . This implies that SDP optimization problems can be tackled with standard numerical packages. In this paper, the optimization of the SDP (3) can be solved by using the Matlab parser YALMIP 65 with the solvers, SEDUMI 66 or SDPT3 67,68 . In fact, there exist serval SDP problems in quantum information theory. For example, SDP programs have been used in entanglement detection and quantification [69][70][71][72][73][74] , quantifying quantum resources 75 . Furthermore, the SDP (3) has also been used for calculating the fidelity of quantum channels 76 .

Geometric measures of coherence. In a d-dimension
Hilbert space H with its corresponding reference , a state is incoherent if and only if it is a diagonal density matrix under the reference basis 2,3 . All incoherent states can be represented as 2,3 Thus, the geometric measure of coherence is defined as 23 with the maximum being taken over all possible incoherent states δ ∈ I . Based on the SDP (3), Eqs. (4) and (5), we provide the MATLAB code for the semidefinite program of geometric measure of coherence in Supplemental Material 1.
For an arbitrary single-qubit state ρ , its analytical solutions of C g (ρ) has been derived 23 where ρ 01 is the off-diagonal element of ρ in reference basis. To compare with this analytical solution, we randomly generate 10 5 density matrices and calculate their C g (ρ) by analytical and numerical methods, respectively. The analytical results are calculated based on Eq. (6), and the numerical results are obtained by optimizing the semidefinite program 77 . The maximum deviation between the analytical and numerical results is 3.19 × 10 −9 . For a pure state |ψ� = i i |i� , its geometric measure of coherence is where | i | 2 is the diagonal elements of |ψ��ψ| 3 . However, the corresponding analytical solutions of C g (ρ) are difficult to calculate for general mixed states, so it is necessary to get some lower and upper bounds of C g (ρ) . Here we employ the lower and upper bounds proposed in Ref. 59 and compare them with our numerical results of the optimization program. For a general d × d density matrix ρ , its C g (ρ) satisfies 59 where b ii is from √ ρ = ij b ij |i��j|. Since there is no corresponding analytic solution for general single-qutrit mixed states, we randomly generate 10 5 density matrices to draw their corresponding upper and lower bounds. In Fig. 1, there is a clear dividing line between two bounds indicating that the numerical results obtained by our algorithm coincide with the analytical results from the inequality (8), and points on the upper bound is closer to the dividing line than points on the lower bound for many 3 × 3 density matrices. Now we consider the following mixed state i=0 |i� , I being the d × d identity matrix, and 0 ≤ p ≤ 1 . Since the mixed state ρ is highly symmetric, it will remain unchanged when we exchange its basis order. It limits that the reference incoherent state in geometric measure of coherence must have the same diagonal elements, i.e., the closest incoherent state δ to the density matrix ρ has to be Therefore, we can obtain its analytical solution of geometric measure of coherence.

Proposition 1 For the mixed state
i=0 |i� , its analytical solution of geometric measure of coherence is www.nature.com/scientificreports/ This analytical result is equal to its corresponding upper bound which is the right hand side of the inequality (8). When 2 ≤ d ≤ 20 , we calculate their analytical and numerical results as well as maximum deviation between them. For d = 3 the corresponding graph is drawn and the rest have the similar phenomena like it. In Fig. 2, C g (ρ) and its upper bound coincide for d = 3 , and the maximum deviation between them is 1.51 × 10 −9 . In Fig. 3, the maximum deviation between the numerical and analytical results is about 10 −9 orders of magnitude. Although the average time    The GME is defined as follows 80 where S is the set of all separable states. We replace S with the set P of all PPT states, because S cannot be easily expressed in the semidefinite programs, but P can be expressed since for a given density matrix one can directly calculate its partial transpose 41 . Thus, based on the fact that S is a subset of P 79 , one can obtain a lower bound of E G (ρ) , i.e., where the lower bound The equality in Eq. (17) holds for all two-qubit and qubit-qutrit states 81 , and some special kinds of higher dimensional mixed states. Based on the SDP (3), Eqs. (15) and (16), we provide the MATLAB code for the semidefinite program of ∼ E G (ρ) in Supplemental Material 1. For pure states, the GME is defined as 9 Moreover, it is defined via the convex roof construction for mixed states. If ρ is a two-qubit state, the corresponding expression of E G (ρ) is 9,82,83 The C(ρ) is called concurrence that its expression is where { i } are the square root of eigenvalues of ρ ∼ ρ in descending order and ∼ ρ = (σ y ⊗ σ y )ρ * (σ y ⊗ σ y ) . In order to compare the analytical result Eq. (20) and the numerical result for two-qubit states, we randomly generate 10 5 density matrices and calculate the analytical and numerical results respectively. The maximum difference between them is 1.57 × 10 −9 . Now we apply our semidefinite program to the isotropic states, where the forms of these states are 9 with the maximally entangled state | + � = 1 i=0 |ii� and 0 ≤ F ≤ 1 . The analytical solutions for the GME of these isotropic states were given in Ref. 9 , and the states are separable if and only if F ≤ 1 d 84 . For Eq. (22) when 2 ≤ d ≤ 5 , we calculate that the maximum deviation between the numerical and analytical solution by our semidefinite program and the analytical solution given in 9 , respectively. The results are summarized in Table 1, where ∼ t (s) denotes the average time of each operation. In the example tested above, the semidefinite program always obtain the same value as E G (ρ) within the precision given in Table 1.
Finally, we apply semidefinite program to the Werner states that it can be expressed as a linear combination of two operators of the identity I and the swap F ≡ ij |ij��ji| 9 , i.e., ρ = aI + b F , where a and b are both real coefficients and are limited by Trρ = 1 . When one of the parameters is considered, the states can be expressed as with f ≡ Tr(ρF) . The corresponding analytic solution for the Werner states (23) is 9 (15) (23), so the maximum deviation between the numerical solution and the analytical solution is calculated, respectively. The results are summarized in Table 2, where ∼ t (s) is the average time of each operation.

Discussion
For the geometric measure of entanglement, we compare our algorithm with the algorithm proposed in Ref. 47 . Streltsov and colleagues proposed an algorithm 47 , which can be easily implemented by solving an eigenproblem or finding a singular value decomposition of a matrix. However, their algorithm needs the iteration with many steps, and it may converge to a local minimum which is not the exact value of the geometric measure of entanglement. Our algorithm, which does not need iteration and has no local minimum problems, is based on semidefinite program and easy to implement. Unfortunately, the shortcoming of our algorithm is also obvious. It can be used to calculate the geometric measure of entanglement for arbitrary two-qubit and qubit-qutrit states, and some special kinds of higher dimensional mixed states. But for other states, our algorithm can only get a lower bound of the geometric measure of entanglement.
In this paper, we introduced numerical and analytical results to compute the geometric measures of coherence and the entanglement. In coherence measures, the deviation between the numerical solution and the analytical solution was an order of magnitude of 10 −9 for single-qubit states. Furthermore, we obtained the analytical solution of the geometric measure of coherence C g (ρ) 2 for the special kind of mixed states ρ = p|ψ + ��ψ + | + (1 − p) I d . For randomly generated 3-dimensional density matrices, we have drawn a boundary diagram with a apparently clear boundary line. In entanglement measures, we used PPT states to replace the set of separable states and calculated two-qubit states, the isotropic states and the Werner states by using fidelity and its semidefinite program, and then concluded that their maximum deviation is almost on the order of magnitude of 10 −9 .