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

TL;DR

This paper develops numerical algorithms and analytical solutions to quantify geometric measures of coherence and entanglement in quantum states, enabling better analysis of quantum information resources.

Contribution

It introduces semidefinite programming algorithms for calculating geometric coherence and entanglement, along with analytical solutions for specific mixed states.

Findings

Algorithms successfully compute geometric measures for various states.

Numerical results match analytical solutions where available.

Lower bounds are provided for complex states.

Abstract

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.

Figures

Tables

Table 1. Table 1: The maximum deviation between the numerical solution and the analytical solution with Δ ​ E G ​ ( ρ ) = E G ​ ( ρ ) − E ∼ G ​ ( ρ ) Δ subscript 𝐸 𝐺 𝜌 subscript 𝐸 𝐺 𝜌 subscript similar-to 𝐸 𝐺 𝜌 \Delta E_{G}(\rho)=E_{G}(\rho)-\overset{\sim}{E}_{G}(\rho) , and the average time t ∼ ​ ( s ) similar-to 𝑡 𝑠 \overset{\sim}{t}(s) of each operation for the isotropic states ( 22 ).

d	2	3	4	5
$\Delta E_G\left(\rho \right)$	$1.03\times {10}^{-9}$	$2.02\times {10}^{-9}$	$3.21\times {10}^{-9}$	$5.05\times {10}^{-9}$
$\stackrel{\sim }{t}\left(s\right)$	$0.29$	$0.49$	$2.23$	$27.35$

Table 2. Table 2: The maximum deviation between the numerical solution and the analytical solution with Δ ​ E G ​ ( ρ ) = E G ​ ( ρ ) − E ∼ G ​ ( ρ ) Δ subscript 𝐸 𝐺 𝜌 subscript 𝐸 𝐺 𝜌 subscript similar-to 𝐸 𝐺 𝜌 \Delta E_{G}(\rho)=E_{G}(\rho)-\overset{\sim}{E}_{G}(\rho) , and the average time t ∼ ​ ( s ) similar-to 𝑡 𝑠 \overset{\sim}{t}(s) of each operation for the Werner states ( 23 ).

d	2	3	4	5
$\Delta E_G\left(\rho \right)$	$5.00\times {10}^{-10}$	$2.26\times {10}^{-9}$	$3.17\times {10}^{-9}$	$5.57\times {10}^{-9}$
$\stackrel{\sim }{t}\left(s\right)$	$0.32$	$0.46$	$1.47$	$14.83$

Equations

$$F(\rho ,\chi )\equiv \mathrm{Tr}\sqrt{{\rho }^{\frac{1}{2}}\chi {\rho }^{\frac{1}{2}}}.$$

$$\begin{array}{r}F(|\psi ⟩,\chi )=\sqrt{⟨\psi |\chi |\psi ⟩}.\end{array}$$

$$\mathrm{maximize}:$$

$$\mathrm{subject}\mathrm{to}:$$

$$X\in L(\mathbb{X}),$$

$$\rho ,\chi \in Pos(\mathbb{X}),$$

$$\delta =\sum _{i=0}^{d-1}{p}_i|i⟩⟨i|.$$

$$C_g(\rho )=1-[\underset{\delta \in \mathbb{I}}{\max }F(\rho ,\delta ){]}^2,$$

$$C_g(\rho )=\frac{1}{2}(1-\sqrt{1-4|{\rho }_{01}{|}^2}),$$

$$C_g(|\psi ⟩)=1-\underset{i}{\max }\{|{\lambda }_i{|}^2\},$$

$$1-\frac{1}{d}-\frac{d-1}{d}\sqrt{1-\frac{d}{d-1}(\mathrm{Tr}{\rho }^2-\sum _i{\rho }_{ii}^2)}\le C_g(\rho )\le \min \{1-\underset{i}{\max }\{{\rho }_{ii}\},1-\sum _i{b}_{ii}^2\},$$

$$\rho =p|{\psi }^{+}⟩⟨{\psi }^{+}|+(1-p)\frac{\text{mathbbm}1}{d},$$

$$\delta =\sum _{i=0}^{d-1}\frac{1}{d}|i⟩⟨i|.$$

$$C_g(\rho )=1-\frac{1}{d^2}[(d-1)\sqrt{1-p}+\sqrt{1+(d-1)p}{]}^2.$$

$$|{\psi }^{AB}⟩=|{\psi }^A⟩\otimes |{\psi }^B⟩.$$

$${\rho }^{AB}=\sum _i{p}_i{\rho }_i^A\otimes {\rho }_i^B,$$

$${\rho }_{A|B}=\sum _{ij,kl}{\rho }_{ij,kl}|i⟩⟨j{|}_A\otimes |k⟩⟨l{|}_B,$$

$${\rho }_{A|B}^{T_A}=\sum _{ij,kl}{\rho }_{ij,kl}|j⟩⟨i{|}_A\otimes |k⟩⟨l{|}_B\ge 0.$$

$$E_G(\rho )=1-[\underset{\sigma \in \mathbb{S}}{\max }F(\rho ,\sigma ){]}^2,$$

$$E_G(\rho )\ge {\stackrel{\sim }{E}}_G(\rho ),$$

$${\stackrel{\sim }{E}}_G(\rho )=1-[\underset{\tilde{\sigma }\in \mathbb{P}}{\max }F(\rho ,\tilde{\sigma }){]}^2.$$

$$E_G(|\psi ⟩)=1-\underset{|\varphi ⟩\in \mathbb{S}}{\max }|⟨\psi |\varphi ⟩{|}^2,$$

$$E_G(\rho )=\frac{1}{2}(1-\sqrt{1-C(\rho {)}^2}).$$

$$C(\rho )=\max \{0,{\lambda }_1-{\lambda }_2-{\lambda }_3-{\lambda }_4\},$$

$$\rho =\frac{1-F}{d^2-1}(\text{mathbbm}1-|{\Phi }^{+}⟩⟨{\Phi }^{+}|)+F|{\Phi }^{+}⟩⟨{\Phi }^{+}|,$$

$$\rho =\frac{d^2-fd}{d^4-d^2}\text{mathbbm}1\otimes \text{mathbbm}1+\frac{f{d}^2-d}{d^4-d^2}\hat{F},$$

$$E_G(\rho )=\frac{1}{2}(1-\sqrt{1-f^2}),$$

Full text

Introduction

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}(\rho)=\min_{\delta\in\mathbb{I}}D(\rho,\delta)$, i.e., the minimum distance from $\rho$ to all possible incoherent states $\delta\in\mathbb{I}$, where $\mathbb{I}$ is the set of all incoherent states [18, 19, 20, 21, 22]. From this definition, one can see that $C_{D}(\rho)=0$ if and only if $\rho$ 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 $\rho$ 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 [26, 25, 27, 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 $\mathbb{S}$ of separable states, the less entanglement it has [19, 22]. According to the distance measure $D$ between quantum states $\rho$ and $\sigma$, it is defined as $E_{D}(\rho)=\min_{\sigma\in\mathbb{S}}D(\rho,\sigma)$, 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) [39, 40, 41]. Furthermore, the expected value of entanglement witnesses can be used to estimate the GME [42, 43, 44]. Other different entanglement measures have been proposed for multipartite systems and mixed states [19, 20].

Most quantum states have no corresponding analytical solutions for the coherence and entanglement measures, so numerical algorithms must be applied. In some entanglement measures, several numerical algorithms have been used to solve related problems [45, 46, 47, 48]. Moreover, computing many entanglement measures is NP hard for a general state [49, 50], so some upper and lower bounds are proposed to describe entanglement [51, 52, 53, 54, 55, 56, 57, 58] and coherence measures [59, 60, 61, 62, 63]. In Refs. [64, 65], a semidefinite program (SDP) was proposed to calculate the fidelity between two states, so the geometric measures of coherence and entanglement probably can be numerically obtained based on the semidefinite program. 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, the algorithm can get a lower bound of the geometric measure of entanglement.

Results

A semidefinite program for computing fidelity. We first review the fidelity and its semidefinite program in the following. The fidelity between states $\rho$ and $\chi$ is defined as [1]

[TABLE]

For a pure state $|\psi\rangle$ and an arbitrary state $\chi$, one can get that

[TABLE]

In Refs. [64, 65], 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

[TABLE]

where $L$ is the collection of all linear mappings. In a $\mathbb{X}$ complex Hilbert space, $Pos(\mathbb{X})$ is the set of positive semidefinite operators operating on $\mathbb{X}$. Then the maximum value of $\frac{1}{2}{\rm Tr}(X)+\frac{1}{2}{\rm Tr}(X^{\dagger})$ is equal to $F(\rho,\chi)$. $X$ is a randomly generated complex matrix of the same order as $\chi$.

The SDP can not only solve the problem effectively, but also prove the global optimality under weak conditions [66]. This implies that SDP optimization problems can be tackled with standard numerical packages. In this paper, the optimization of the SDP (Results) can be solved by using the Matlab parser YALMIP [67] with the solvers, SEDUMI [68] or SDPT3 [69, 70]. In fact, there exist serval SDP problems in quantum information theory. For example, SDP programs have been used in entanglement detection and quantification [71, 72, 73, 74, 75, 76], quantifying quantum resources [77]. Furthermore, the SDP (Results) has also been used for calculating the fidelity of quantum channels [78].

Geometric measures of coherence. In a $d$-dimension Hilbert space $\mathcal{H}$ with its corresponding reference basis $\{|i\rangle\}_{i=0}^{d-1}$, 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]

[TABLE]

Thus, the geometric measure of coherence is defined as [23]

[TABLE]

with the maximum being taken over all possible incoherent states $\delta\in\mathbb{I}$. Based on the SDP (Results), Eqs. (4) and (5), we provide the MATLAB code for the semidefinite program of geometric measure of coherence in Supplemental Material.

For an arbitrary single-qubit state $\rho$, its analytical solutions of $C_{g}(\rho)$ has been derived [23]

[TABLE]

where $\rho_{01}$ is the off-diagonal element of $\rho$ in reference basis. To compare with this analytical solution, we randomly generate $10^{5}$ density matrices and calculate their $C_{g}(\rho)$ 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 [79]. The maximum deviation between the analytical and numerical results is $3.19\times 10^{-9}$.

For a pure state $|\psi\rangle=\sum_{i}\lambda_{i}|i\rangle$, its geometric measure of coherence is

[TABLE]

where $|\lambda_{i}|^{2}$ is the diagonal elements of $|\psi\rangle\langle\psi|$ [3]. However, the corresponding analytical solutions of $C_{g}(\rho)$ are difficult to calculate for general mixed states, so it is necessary to get some lower and upper bounds of $C_{g}(\rho)$. Here we employ the lower and upper bounds proposed in Ref. [61] and compare them with our numerical results of the optimization program. For a general $d\times d$ density matrix $\rho$, its $C_{g}(\rho)$ satisfies [61]

[TABLE]

where $b_{ii}$ is from $\sqrt{\rho}=\sum_{ij}b_{ij}{|i\rangle\langle 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\times 3$ density matrices.

Now we consider the following mixed state

[TABLE]

with $|\psi^{+}\rangle=\frac{1}{\sqrt{d}}\sum_{i=0}^{d-1}|i\rangle$, $\mathbbm{1}$ being the $d\times d$ identity matrix, and $0\leq p\leq 1$. Since the mixed state $\rho$ 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 $\delta$ to the density matrix $\rho$ has to be

[TABLE]

Therefore, we can obtain its analytical solution of geometric measure of coherence.

Proposition 1. For the mixed state $\rho=p|\psi^{+}\rangle\langle\psi^{+}|+(1-p)\frac{\mathbbm{1}}{d}$ with $|\psi^{+}\rangle=\frac{1}{\sqrt{d}}\sum_{i=0}^{d-1}|i\rangle$, its analytical solution of geometric measure of coherence is

[TABLE]

This analytical result is equal to its corresponding upper bound which is the right hand side of the inequality (8). When $2\leq d\leq 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}(\rho)$ and its upper bound coincide for $d=3$, and the maximum deviation between them is $1.51\times 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 $\overset{\sim}{t}(s)$ of each operation increases exponentially, it is within an acceptable range in the low dimensional case.

Geometric measures of entanglement. A separable bipartite pure state can be written in the following product form

[TABLE]

For mixed states, if it can be represented as convex weights $p_{i}$ and product states $\rho_{i}^{A}\otimes\rho_{i}^{B}$ [80]

[TABLE]

then $\rho^{AB}$ is called separable. For a bipartite state, if there is no negative eigenvalues after the partial transposition of subsystem $A$, this bipartite state is called the PPT states [81], i.e., a bipartite state

[TABLE]

is PPT, when its partial transposition with respect to the subsystem $A$ satisfys

[TABLE]

The GME is defined as follows [82]

[TABLE]

where $\mathbb{S}$ is the set of all separable states. We replace $\mathbb{S}$ with the set $\mathbb{P}$ of all PPT states, because $\mathbb{S}$ cannot be easily expressed in the semidefinite programs, but $\mathbb{P}$ can be expressed since for a given density matrix one can directly calculate its partial transpose [42]. Thus, based on the fact that $\mathbb{S}$ is a subset of $\mathbb{P}$ [81], one can obtain a lower bound of $E_{G}(\rho)$, i.e.,

[TABLE]

where the lower bound $\overset{\sim}{E}_{G}(\rho)$ is defined by

[TABLE]

The equality in Eq. (17) holds for all two-qubit and qubit-qutrit states [83], and some special kinds of higher dimensional mixed states. Based on the SDP (Results), Eqs. (15) and (16), we provide the MATLAB code for the semidefinite program of $\overset{\sim}{E}_{G}(\rho)$ in Supplemental Material.

For pure states, the GME is defined as [9]

[TABLE]

Moreover, it is defined via the convex roof construction for mixed states. If $\rho$ is a two-qubit state, the corresponding expression of $E_{G}(\rho)$ is [9, 84, 85]

[TABLE]

The $C(\rho)$ is called concurrence that its expression is

[TABLE]

where $\{\lambda_{i}\}$ are the square root of eigenvalues of $\rho\overset{\sim}{\rho}$ in descending order and $\overset{\sim}{\rho}=(\sigma_{y}\otimes\sigma_{y})\rho^{\ast}(\sigma_{y}\otimes\sigma_{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\times 10^{-9}$.

Now we apply our semidefinite program to the isotropic states, where the forms of these states are [9]

[TABLE]

with the maximally entangled state $|\Phi^{+}\rangle=\frac{1}{\sqrt{d}}\sum_{i=0}^{d-1}|ii\rangle$ and $0\leq F\leq 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\leq\frac{1}{d}$ [86]. For Eq. (22) when $2\leq d\leq 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 $\overset{\sim}{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}(\rho)$ 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 $\mathbbm{1}$ and the swap $\hat{F}\equiv\sum_{ij}|ij\rangle\langle ji|$[9], i.e., $\rho=a\mathbbm{1}+b\hat{F}$, where $a$ and $b$ are both real coefficients and are limited by ${\rm Tr}{\rho}=1$. When one of the parameters is considered, the states can be expressed as

[TABLE]

with $f\equiv{\rm Tr}(\rho\hat{F})$. The corresponding analytic solution for the Werner states (23) is [9]

[TABLE]

where $f\leq 0$ and [math] otherwise. For $2\leq d\leq 5$, we apply our semidefinite program to states (23), so the maximum deviation between the numerical solution and the analytical solution is calculated, respectively. The results are summarized in Table 2, where $\overset{\sim}{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. [48]. Streltsov and colleagues proposed an algorithm [48], 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}(\rho)=1-\frac{1}{d^{2}}[(d-1)\sqrt{1-p}+\sqrt{1+(d-1)p}]^{2}$ for the special kind of mixed states $\rho=p|\psi^{+}\rangle\langle\psi^{+}|+(1-p)\frac{\mathbbm{1}}{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}$.

Acknowledgements

C.Z. gratefully acknowledges Otfried Gühne and Haidong Yuan for helpful discussions. This work is supported by the National Natural Science Foundation of China (Grant No. 11734015), and K.C.Wong Magna Fund in Ningbo University.

Author contributions statement

Z.Z., Y.-L.D. and C.Z. wrote the main manuscript text and Y.D. prepared figures 1-3. All authors reviewed the manuscript.

Additional information

Competing financial interests: The authors declare no competing interests.

Corresponding authors

Correspondence and requests for materials should be addressed to Y.-L.D. (email:yldong@suda.edu.cn) or C.Z. (email:chengjie.zhang@gmail.com).

Supplemental Material

Here we provide some detailed calculations of the single-qubit states, randomly generate $3$-dimensional density matrices, and a special kind of $d$-dimensional density matrices for coherence measures. The corresponding MATLAB code for the semidefinite program of geometric measure of coherence is

⬇

1ction [C, D]=coherence(rho)

2

3 d=length(rho);

4

5 delt=sdpvar(d,d, 'diagonal', 'real');

6

7 X=sdpvar(d,d, 'full', 'complex');

8

9 constr=[[rho X; ctranspose(X) delt]$\geq$0, delt$\geq$0, trace(delt)==1];

10

11 result=solvesdp(constr, -trace(X)-trace(ctranspose(X)), sdpsettings('verbose', 1));

12

13%check for errors

14 if (result. problem $\neq$ 0 )

15 disp(result. info);

16 end

17

18%return the maximum fidelity

19 F=double((trace(X)+trace(ctranspose(X)))/2);

20

21 D=double(delt);

22

23 C=1-F^2;

24

25end

We have used the parser YALMIP [1] with the solvers, SEDUMI [2] or SDPT3 [3].

We also provide some detailed calculations of the two-qubit states, the isotropic states and the Werner states for entanglement measures. The corresponding MATLAB code for the semidefinite program of geometric measure of entanglement is

⬇

1ction [E, E2]=GME(rho)

2

3%d(i) is the dimension of the matrix

4 d1=2;

5 d2=2;

6

7 sigma=sdpvar(d1d2, d1d2, 'hermitian', 'complex');

8

9 X=sdpvar(d1d2, d1d2, 'full', 'complex');

10

11%performing partial transpose for sigma and the pt subroutine comes from pptmixer.

12 constr=[[rho X; ctranspose(X) sigma]$\geq$0, sigma$\geq$0, pt(sigma, [1,0], [d1, d2])$\geq$0, trace(sigma)==1];

13

14 result=solvesdp(constr, -trace(X)-trace(ctranspose(X)), sdpsettings('verbose', 1));

15

16%check for errors

17 if (result. problem $\neq$ 0)

18 disp(result. info);

19 end

20

21%return the maximum fidelity

22 F=double((trace(X)+trace(ctranspose(X)))/2);

23

24 D=double(sigma);

25

26 E=1-F^2;

27

28end

We have used the pt subroutine from the program PPTmixer, where the PPTmixer was presented in [4] (or [5]).