Two-qubit trace-norm geometric discord: the complete solution

TL;DR

This paper provides a comprehensive solution for calculating the trace-norm geometric discord in any two-qubit state, simplifying the problem to analyzing critical points of a specific mapping.

Contribution

It introduces a complete method for determining trace-norm geometric discord in two-qubit states, extending previous partial solutions.

Findings

Explicit formulas for trace-norm geometric discord for various two-qubit states

Reduction of the problem to critical point analysis of a specific mapping

Comparison with known solutions demonstrating the method's effectiveness

Abstract

We present the complete solution of the problem of determination of trace-norm geometric discord for arbitrary two-qubit state. Final answer is achieved due to effective reduction of the problem to the study of critical points of certain mapping depending on projectors. Our results are illustrated on various, also new, families of two-qubit states and compared to already known special solutions.

Equations

$$⟨n,\sigma ⟩=n_1{\sigma }_1+n_2{\sigma }_2+n_3{\sigma }_3,\sigma =({\sigma }_1,{\sigma }_2,{\sigma }_3),n=(n_1,n_2,n_2)$$

$$⟨n,\sigma ⟩⟨m,\sigma ⟩=⟨n,m⟩{\text{openone}}_2+i⟨n\times m,\sigma ⟩$$

$$\varrho =\frac{1}{2}\left({\text{openone}}_2+⟨n,\sigma ⟩\right),n\in {\mathbb{R}}^3$$

$$\varrho =\frac{1}{4}\left({\text{openone}}_2\otimes {\text{openone}}_2+⟨x,\sigma ⟩\otimes {\text{openone}}_2+{\text{openone}}_2\otimes ⟨y,\sigma ⟩+\sum _{j,k=1}^3{K}_{jk}{\sigma }_j\otimes {\sigma }_k\right)$$

$$x_j=\mathrm{tr}(\varrho {\sigma }_j\otimes {\text{openone}}_2),y_j=\mathrm{tr}(\varrho {\text{openone}}_2\otimes {\sigma }_j),K_{jk}=\mathrm{tr}(\varrho {\sigma }_j\otimes {\sigma }_k)$$

$$⟨Um,\sigma ⟩=\hat{U}⟨m,\sigma ⟩{\hat{U}}^{\ast }$$

$$\varrho \to {\hat{U}}^{\ast }\otimes {\hat{V}}^{\ast }\varrho \hat{U}\otimes \hat{V},\hat{U},\hat{V}\in \mathrm{SU}(2)$$

$$x\to x^{\prime }=Ux,y\to y^{\prime }=Vy,K\to K^{\prime }=UK{V}^{\mathrm{T}}$$

$$D_1(\varrho )=\underset{{\mathbb{P}}_{\mathcal{A}}}{\min }||\varrho -{\mathbb{P}}_{\mathcal{A}}(\varrho )|{|}_1$$

$$P_1+P_2={\text{openone}}_2,P_j{P}_k={\delta }_{jk}{P}_k$$

$$\mathbb{P}(A)=P_{1}A{P}_1+P_{2}A{P}_2$$

$$P_k=u{P}_k^0{u}^{\ast }\text{for some}u\in \mathrm{SU}(2)$$

$$P_1^0=\left(\begin{array}{cc}1& 0\\ 0& 0\end{array}\right),P_2^0=\left(\begin{array}{cc}0& 0\\ 0& 1\end{array}\right)$$

$$⟨\mathcal{P}m,\sigma ⟩=\mathbb{P}(⟨m,\sigma ⟩),m\in {\mathbb{R}}^3$$

$${\mathcal{P}}_0=\mathrm{diag}(0,0,1)$$

$$\mathcal{P}=V{\mathcal{P}}_0{V}^{\mathrm{T}},V\in \mathrm{SO}(3)$$

$${\mathcal{M}}_0={\text{openone}}_3-{\mathcal{P}}_0,\mathcal{M}={\text{openone}}_3-\mathcal{P}$$

$$\mathcal{M}=V{\mathcal{M}}_0{V}^{\mathrm{T}},V\in \mathrm{SO}(3)$$

$$\dim \mathrm{Ran}{\mathcal{M}}_0=\dim \mathrm{Ran}\mathcal{M}=2$$

$$S(\mathcal{M})=\varrho -{\mathbb{P}}_{\mathcal{A}}(\varrho )=\frac{1}{4}\left(⟨\mathcal{M}x,\sigma ⟩\otimes {\text{openone}}_2+\sum _{k=1}^3⟨\mathcal{M}K{e}_k,\sigma ⟩\otimes ⟨e_k,\sigma ⟩\right)$$

$$D_1(\varrho )=\underset{\mathcal{M}}{\min }\mathrm{tr}|S(\mathcal{M})|=\underset{\mathcal{M}}{\min }\mathrm{tr}\sqrt{Q(\mathcal{M})}$$

$$\underset{\mathcal{M}}{\min }\mathrm{tr}\sqrt{Q(\mathcal{M})}=\underset{\mathcal{M}}{\min }\mathrm{tr}\sqrt{Q^{\prime }(\mathcal{M})}$$

$$||S(\mathcal{M})|{|}_1=\frac{1}{\sqrt{2}}\sqrt{||\mathcal{M}x|{|}^2+\mathrm{tr}(\mathcal{M}K{K}^{\mathrm{T}})+\sqrt{{\left[||\mathcal{M}x|{|}^2+\mathrm{tr}(\mathcal{M}K{K}^{\mathrm{T}})\right]}^2-4\left[||K^{\mathrm{T}}\mathcal{M}x|{|}^2+\mathrm{tr}(E^{\mathrm{T}}E-\mathcal{M}{E}^{\mathrm{T}}E)\right]}}$$

$$Q(\mathcal{M})=\frac{1}{16}\left((\mathrm{tr}(\mathcal{M}K{K}^{\mathrm{T}})+⟨\mathcal{M}x,x⟩){\text{openone}}_2\otimes {\text{openone}}_2+{\text{openone}}_2\otimes 2⟨K^{\mathrm{T}}\mathcal{M}x,\sigma ⟩+\sum _{j,k}\mathrm{tr}(K^{\mathrm{T}}\mathcal{M}{F}_j\mathcal{M}K{F}_k){\sigma }_j\otimes {\sigma }_k\right)$$

$$(F_j{)}_{kl}=-{\epsilon }_{jkl},j,k,l=1,2,3$$

$$U^{\mathrm{T}}{F}_{j}U=\sum _k{U}_{jk}{F}_k$$

$$x_0=U_{0}x\text{and}{I}_0=U_{0}K{V}_0^{\mathrm{T}}$$

$$I_0=\mathrm{diag}(i_1,i_2,i_3)$$

$$E_0=\mathrm{diag}({ϵ}_1,{ϵ}_2,{ϵ}_3)$$

$$\mathrm{tr}(I_0\mathcal{M}{F}_j\mathcal{M}{I}_0{F}_k)=-2{ϵ}_k{V}_{j3}{V}_{k3}$$

Full text

qubits, quantum correlations, geometric quantum discord, trace norm

I Introduction

It is a common belief that question of the explicit quantitative characterization of correlations for a two-qubit systems can always be answered, but actually in many cases we land on a uncharted territory, where the explicit answer is not known Modi ; Adesso ; HHWPZF . One of these situations is the strict solution of the question of trace-norm quantum geometric discord $D_{1}$ for an arbitrary state of the two-qubit system. The two-qubit system is the simplest nontrivial compound system which we have at our disposal to study quantum correlations and all relevant tools invented to measure them Modi . The geometric quantum discord measured by means of the distance induced by the trace-norm seems to be one of the solid candidates. Contrary to the geometric discord defined by means of the Hilbert-Schmidt distance, it is well defined and satisfies requirements of the bona fide measure of quantum correlations. As it is well known, the Hilbert-Schmidt distance yields to an anomalous behavior of related to it geometric discord $D_{2}$, but it has an attractive property, that it can be calculated relatively easy. The procedure for calculating geometric discord involves the minimization which for generic states for the trace-distance is rather challenging. As we have shown for qutrits and higher dimensional bipartite systems can hardly be operationally/analytically performed, except for selected families of states LFJ .

In the present work we show how explicitly the $D_{1}$ can be determined, we reveal the geometry of the problem of minimization, what in turn deepens the understanding of the intrinsic geometry of the two-qubit quantum state space. For the Hilbert-Schmidt distance quantum discord $D_{2}$ such a question has been answered in a compact form by Dakić et al. in 2010 GQD , namely the explicit formula for the $D_{2}$ has been obtained. The ease of computation does not heal the nonphysical behavior of the $D_{2}$ under the local evolution of the system lj-hs and its non-contractibility under completely positive trace preserving mappings.

The trace-norm measurement-induced geometric discord is well defined correlation measure. For two qubits it is also equal to the standard geometric discord defined as the distance from a given state to the set of classical-quantum states Na ; RSI and it is more convenient to use quantity based on the disturbance induced by measurement. Such quantity will be simply called geometric quantum discord. However, its calculation is still rather challenging, and by now no explicit answer for arbitrary two-qubit state has been known. According to our best knowledge supported by the recent review on the quantum correlations and geometric quantum discord HHWPZF , prior to the present work only few answers were found for selected families of states:

• (a)

for the Bell diagonal states (or states with maximally mixed marginals) Paula ,

• (b)

for the X-shaped states cicca ,

• (c)

for the states with the correlation matrix having only one non-zero singular value and arbitrary Bell vectors of the marginals – such class contains in particular so called quantum-classical states cicca .

The objective of this study is to provide the complete solution of the problem of giving the explicit value of geometric quantum discord $D_{1}(\rho)$ for arbitrary two-qubit state $\rho$. We find the answer to the minimization problem which has nice geometric meaning. As an illustration, we analyse previously known partial solutions for mentioned above families of states and indicate mechanism making that these solutions were so much simpler.

The paper is organized as follows. In Section 2, we briefly review basic notions and fix the relevant notation. Then, in Section 3, we prove the theorem giving characterization of the central formula for the trace-norm of the disturbance of the system by a local projective measurements. To guide the reader gently to the full solution, in Section 4, we compute the $D_{1}(\rho)$ for the relatively large family of states for which the minimization procedure is not necessary at all or is straightforward. The Section 5 is central for the final result, there we study the critical points of the trace-norm disturbance mapping and obtain that there are two possibilities: singular critical points and smooth critical points. As a preliminary application we apply new results to the states with maximally mixed marginals. In the Section 6 we prove the main result: the determination of the trace norm geometric discord in general case. Finally, the Section 7 will allow reader to make contact with known solutions and go beyond the family of the X- shaped states. We finish the paper with some conclusions and comments on the main findings and obstructions to their extensions to higher dimensional systems.

II Preliminary notions

In the description of two-qubit systems we follow general formulation of $d$ - level quantum systems introduced in LFJ , specified to the case $d=2$.

II.1 Two qubits

Qubit is a $2$- level quantum system. The corresponding Hilbert space equals to $\mathbb{C}^{2}$ and the observables are given by hermitian elements of full matrix algebra $\mathbb{M}_{2}(\mathbb{C})$. As a basis in $\mathbb{M}_{2}(\mathbb{C})$ we use standard Pauli matrices $\sigma_{1},\,\sigma_{2},\,\sigma_{3}$ and the identity matrix $\openone_{2}$. In the following we will use the notation

[TABLE]

Then from the properties of Pauli matrices it follows that

[TABLE]

The set $\mathcal{E}_{2}$ of all states of $2$ - level system can be parametrized as follows

[TABLE]

where $||n||\leq 1$. So $\mathcal{E}_{2}$ is given by the unit ball in $\mathbb{R}^{3}$ and the pure states correspond to the unit sphere $||n||=1$.

Consider now two qubits ${\mathcal{A}}$ and ${\mathcal{B}}$. It is convenient to parametrize the set of states of composite system as follows

[TABLE]

where $x,\,y\in\mathbb{R}^{3}$ and $K=(K_{jk})$ is the correlation matrix. Notice that

[TABLE]

The parametrization (II.3) is chosen in such a way that the marginals $\mathrm{tr}_{{\mathcal{A}}}\varrho$ and $\mathrm{tr}_{{\mathcal{B}}}\varrho$ are given by the vectors $x$ and $y$ as in (II.2).

Let us discuss now the adjoint representation of the group $\mathrm{SU(2)}$. Let $\hat{U}\in\mathrm{SU(2)}$ and define $3\times 3$ matrix $U$ by

[TABLE]

The matrix $U$ is real and orthogonal. In this way to each $\hat{U}\in\mathrm{SU(2)}$ there corresponds $U\in\mathrm{SO(3)}$. In contrast to general case, when $d=2$ such obtained group is exactly equal to the group $\mathrm{SO(3)}$. Consider now the local transformations of the state of two qubits

[TABLE]

The corresponding vectors $x,\,y$ and the correlation matrix $K$ transform as follows

[TABLE]

where $U,\,V$ are the adjoint representations of $\hat{U},\,\hat{V}$ respectively. Since in the case of qubits we can use the full group of orthogonal transformations to diagonalize correlation matrix, any two-qubit state is locally equivalent to the state with diagonal $K$.

II.2 Trace-norm geometric discord

Let $\varrho$ be a state of bipartite system ${\mathcal{A}}{\mathcal{B}}$. When we perform local measurement on the subsystem ${\mathcal{A}}$, the state $\varrho$ may be disturbed due to such measurement. The trace-norm (one-sided) measurement induced geometric discord is defined as the minimal disturbance induced by projective measurement $\mathbb{P}_{{\mathcal{A}}}$ on subsystem ${\mathcal{A}}$, computed using the trace distance in the set of states. It can be compared with the standard geometric discord equal to the distance from a given state to the set of classical - quantum states GQD . At it was already stated in the Introduction, in the case of qubits these two notions coincide, and it is more convenient to use the quantity based on the disturbance induced by the measurement which will be simply called trace - norm geometric discord. The formal definition is as follows Paula

[TABLE]

where $||A||_{1}=\mathrm{tr}\,\,|A|$.

In the case of qubits, the local projective measurement $\mathbb{P}_{{\mathcal{A}}}$ is given by the one - dimensional projectors $P_{1},\,P_{2}$ on $\mathbb{C}^{2}$, such that

[TABLE]

and $\mathbb{P}_{{\mathcal{A}}}=\mathbb{P}\otimes\mathrm{id}$, where

[TABLE]

One - dimensional projectors $P_{k}$ can be always chosen as

[TABLE]

where

[TABLE]

Define now a real orthogonal projector $\mathcal{P}$ on $\mathbb{R}^{3}$

[TABLE]

If $\mathcal{P}_{0}$ denotes such projector given by (II.7), where we take $P_{1}^{0}$ and $P_{2}^{0}$, then

[TABLE]

and

[TABLE]

Define also orthogonal complements to $\mathcal{P}_{0}$ and $\mathcal{P}$

[TABLE]

Obviously $\mathcal{M}_{0}=\mathrm{diag}(1,1,0)$,

[TABLE]

and

[TABLE]

Notice that in this case the projectors $\mathcal{M}$ run over the whole set of projectors with dimension $2$.

Now we compute the disturbance of the state (II.3) caused by measurement $\mathbb{P}_{{\mathcal{A}}}$. We have

[TABLE]

where $e_{k},\,k=1,2,3$ are the vector of the canonical basis of $\mathbb{R}^{3}$. So

[TABLE]

where $Q(\mathcal{M})=S(\mathcal{M})S(\mathcal{M})^{\ast}$ and the minimum is taken over all projectors $\mathcal{M}$ on two dimensional subspaces of $\mathbb{R}^{3}$.

Remark II.1

Notice that the quantity $D_{1}(\varrho)$ does not depend on the choice of the state within the class of locally equivalent states. In other words

[TABLE]

where $Q^{\prime}$ is obtained from $Q$ by taking vector $x^{\prime}=Ux$ and matrix $K^{\prime}=UKV^{\mathrm{T}}$ instead of $x$ and $K$.

III Trace norm of disturbance $S(\mathcal{M})$

In this section we find the elegant formula for the trace norm of the disturbance $S(\mathcal{M})$ (see cicca for the other version of such formula).

Theorem III.1

The trace norm of the disturbance $S(\mathcal{M})$ of the state (II.3) is given by the formula

[TABLE]

where $||\cdot||$ denotes the euclidian norm in $\mathbb{R}^{3}$ and $E=\mathrm{adj}\,K$ is the adjunct matrix of the correlation matrix $K$ (i. e. the transpose of its cofactor matrix).

Proof: We start the proof of this theorem with the formula for $Q(\mathcal{M})=S(\mathcal{M})S(\mathcal{M})^{\star}$. By a direct computation we obtain

[TABLE]

where $F_{k}$ are the generators of $\mathrm{SO(3)}$, given by

[TABLE]

which satisfy covariance relations

[TABLE]

To further simplify (III.2) consider locally equivalent state for which the correlation matrix is diagonal. So, for a proper transformations $U_{0},\,V_{0}\in\mathrm{SO(3)}$, we have

[TABLE]

and $I_{0}$ is diagonal

[TABLE]

Let $E_{0}=\mathrm{adj}(I_{0})$ be the adjunct matrix of $I_{0}$. Then

[TABLE]

where $\epsilon_{1}=i_{2}i_{3},\,\epsilon_{2}=i_{1}i_{3},\epsilon_{3}=i_{1}i_{2}$. First we show the identity

[TABLE]

where $\mathcal{M}=V\mathcal{M}_{0}V^{\mathrm{T}},\,V\in\mathrm{SO(3)}$. The formula (III.4) follows since

[TABLE]

and

[TABLE]

For any orthogonal transformation $V$ define the mapping $\tau_{V}\,:\,\mathbb{M}_{2}(\mathbb{C})\to\mathbb{M}_{2}(\mathbb{C})$ as follows

[TABLE]

Now using the identity (III.4) we can transform the formula (III.2) to obtain

[TABLE]

Observe that

[TABLE]

where vector $v$ is given by the third column of the matrix $V$. Similarly

[TABLE]

Let us now introduce the matrix

[TABLE]

This matrix has a block diagonal form with blocks

[TABLE]

Put

[TABLE]

then

[TABLE]

The spectral analysis of $R_{0}(\mathcal{M})$ is easy since characteristic polynomials for the blocks are given by

[TABLE]

where the $+(-)$ sign corresponds to the lower (upper) block. One can also check that the vectors $I_{0}\mathcal{M}x_{0}$ and $E_{0}v$ are orthogonal, so

[TABLE]

and the eigenvalues of $R_{0}(\mathcal{M})$ are doubly degenerate

[TABLE]

Thus we get

[TABLE]

where

[TABLE]

The formula (III.8) can be further simplified and we arrive at

[TABLE]

Now applying the inverse to the transformation (III.3), we finally obtain the formula (III.1). $\Box$

IV Solving the minimization problem in the simple case

To compute the geometric discord of a given state, we have to find minimum of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$. In general it is a hard problem and its general solution will be described in the next sections. Now, let us consider a class of two-qubit states for which the solution of this problem is straightforward. This class (denoted by $\mathcal{E}_{0}$) contains the states (II.3) with arbitrary Bloch vectors $x,\,y$ and the correlation matrix of the form

[TABLE]

where $V_{0}$ is some orthogonal matrix and $t$ is a real parameter. From general properties of the two-qubit states it follows that $t$ belongs to the interval $|t|\leq 1$, but the actual value of $t$ depends on the choice of the matrix $V_{0}$ and vectors $x,\,y$. To apply the formula (III.1), notice that

[TABLE]

and

[TABLE]

so

[TABLE]

moreover

[TABLE]

Now the formula (III.1) gives

[TABLE]

so

[TABLE]

and minimal value is achieved for such $\mathcal{M}$ which projects on $x^{\perp}$ in $\mathbb{R}^{3}$. Thus we have

Theorem IV.1

For every state $\varrho\in\mathcal{E}_{0}$, geometric discord is given by

[TABLE]

Characterization of the class $\mathcal{E}_{0}$ is not an easy task. Partial information we can obtain considering locally equivalent states with diagonal correlation matrix. Therefore let us consider the correlation matrix $K$ of the form

[TABLE]

where $I=\mathrm{diag}(i_{1},\,i_{2},\,i_{3}),\;i_{1},i_{2},i_{3}=\pm 1$. In the following we restrict our analysis to $I_{1}=\openone_{3}$ and $I_{2}=\mathrm{diag}(1,-1,1)$. In the first case we obtain the family of states

[TABLE]

where $a,b\in\mathbb{R},\,w,z\in\mathbb{C}$. Similarly, in the second case, we have

[TABLE]

Notice that the relation between parameters $a,\,b,\,w,\,z$ and Bloch vectors $x,\,y$ is as follows

[TABLE]

and

[TABLE]

To obtain more specific information, we must restrict the number of parameters. The simplest is the case when Bloch vectors equal to zero vector. Then (IV.1) contains the one - parameter family of Werner states

[TABLE]

On the other hand, (IV.2) contain the family of isotropic states

[TABLE]

Isotropic states can be obtained as a mixture of maximally mixed state $\varrho_{\infty}=\frac{1}{4}\openone_{4}$ and maximally entangled state $\Psi_{+}$ given by

[TABLE]

i.e.

[TABLE]

On the other hand, the Werner states $\tilde{\varrho}_{\mathrm{W}}$ with reversed parametrization ($t\to-t$) can be obtained as he mixture of $\varrho_{\infty}$ and the state $\Phi_{-}$, where

[TABLE]

So both states are locally equivalent.

More interesting is the case of states with non-trivial Bloch vectors. For simplicity we consider the vectors of the form

[TABLE]

In this case the properties of the states (IV.1) and (IV.2) differ significantly. In particular, $\varrho_{1}$ is positive - definite when

[TABLE]

On the other hand, the positivity region for $\varrho_{2}$ is given by the conditions

[TABLE]

V Critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$

V.1 Formulation of the problem

In this section we start to analyze the real problem of finding minimum of the trace norm $||S(\mathcal{M})||_{1}$ in general case. As a first step we will study critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$. Since the projector $\mathcal{P}=\openone_{3}-\mathcal{M}$ is one dimensional, there is a unit vector $v\in\mathbb{R}^{3}$ such that $\mathcal{P}=P_{v}$. The vector $v$ is given by the third column of the matrix $V\in\mathrm{SO(3)}$ relating $\mathcal{P}$ and $\mathcal{P}_{0}$ (II.9). To study critical points of trace norm of disturbance it is useful to consider the auxiliary function $\mathbf{g}$ defined on the unit sphere $\mathbb{S}^{2}\subset\mathbb{R}^{3}$ with values in $\mathbb{R}^{2}$. The function is defined as follows

[TABLE]

where

[TABLE]

and

[TABLE]

Notice that using the functions (V.2) and (V.3), the formula (III.1) can be rewritten as

[TABLE]

In the following we will consider the mapping $\mathbf{g}$ as a function of two-dimensional projectors $\mathcal{M}$ or vectors $v\in\mathbb{S}^{2}$.

The formulas (V.2) and (V.3) can be further simplified if we introduce the following operators

[TABLE]

and

[TABLE]

where $x$ is a Bloch vector. Now we obtain

[TABLE]

and

[TABLE]

From the equation (V.6) we obtain in particular that

[TABLE]

where $\lambda_{1}\leq\lambda_{2}\leq\lambda_{3}$ are the eigenvalues of real non-negative matrix $L_{+}$. The relation (V.8) allows to obtain the first general result concerning the value of geometric discord $D_{1}(\varrho)$. Namely we have

Theorem V.1

For any two-qubit state (II.3) we have

[TABLE]

where $\lambda_{1}\leq\lambda_{2}\leq\lambda_{3}$ are the eigenvalues of $L_{+}=KK^{\mathrm{T}}+W_{x}$ .

From the relation (V.8) we can also derive well known exact formula for Hilbert - Schmidt norm geometric discord $D_{2}(\varrho)$ GQD . By definition

[TABLE]

Using (III.2) we obtain

[TABLE]

so

[TABLE]

Obviously, for any two-qubit state

[TABLE]

The right hand side of (V.9) in some cases gives not only the lower bound but the exact value of discord $D_{1}$. It can happen for such states $\varrho$ for which $D_{1}(\varrho)=\sqrt{D_{2}(\varrho)}$. Below we show it for an explicit family of states. Consider the states MM

[TABLE]

For the family (V.10), we have

[TABLE]

and

[TABLE]

so

[TABLE]

The order of eigenvalues of (V.11) depends on $\theta$. For $\theta\in[0,\pi/4]$

[TABLE]

and for such $\theta$

[TABLE]

On the other hand, one can check that

[TABLE]

so in this case we obtain only the lower bound. Now for $\theta\in(\pi/4,\,\pi/2]$

[TABLE]

and for such $\theta$

[TABLE]

which equals to the value of $D_{1}(\varrho_{\theta})$.

V.2 Singular critical points

We see that the problem of finding minimum of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ is relatively easy when it is achieved on the set

[TABLE]

Each critical point $v_{0}\in\mathbb{D}$ is called a singular critical point in contrast to the smooth critical point i.e. critical point belonging to the set $\mathbb{S}^{2}\setminus\mathbb{D}$. When the minimum is achieved at some singular critical point, the minimal value of trace norm $||S(\mathcal{M})||_{1}$ is given by the minimum of the function $g_{1}$. To study the properties of the set $\mathbb{D}$, consider another representations of the functions $g_{1}$ and $g_{2}$. Let $L_{-}=KK^{\mathrm{T}}-W_{x}$. One can check that

[TABLE]

and

[TABLE]

where we dropped the term proportional to $1-||v||^{2}$. Using (V.13) and (V.14) one can show that the set (V.12) can be discrete or it can be equal to the whole sphere $\mathbb{S}^{2}$. It follows from the following Lemma, which can be proved by a direct computation.

Lemma V.1

Let $\widetilde{\lambda}_{1},\,\widetilde{\lambda}_{2},\,\widetilde{\lambda}_{3}$ be the eigenvalues of the matrix $L_{-}$. The vectors $v\in\mathbb{D}$ satisfy

[TABLE]

Remark V.1

One can check that in the discrete case, the set $\mathbb{D}$ can contain only one or two projectors, if we identify the vectors $v$ and $-v$.

For the further analysis it is crucial to have a convenient criterion that enables to decide whether the value of $D_{1}(\varrho)$ is achieved at a critical point belonging to the set $\mathbb{D}$ or $\mathbb{S}^{2}\setminus\mathbb{D}$. To obtain such criterion, consider the function

[TABLE]

defined on the region

[TABLE]

For any $t_{\ast}>0$ consider the function

[TABLE]

And finally, define the family of regions $\Omega_{t_{\ast}},\,t_{\ast}\geq 0$ such that $\Omega_{0}=\Omega$ and for $t_{\ast}>0$

[TABLE]

Suppose now that the function $g_{1}$ restricted to $\mathbb{D}$ achieves at some vector $v_{\ast}\in\mathbb{D}$ its absolute minimum. Let $\mathrm{Ran}\,\mathbf{g}$ be the range of the function $\mathbf{g}$. The values of $D_{1}(\varrho)$ is achieved at $v_{\ast}\in\mathbb{D}$ if and only if

[TABLE]

So we obtain

Proposition V.1

Let the function $g_{1}$ restricted to $\mathbb{D}$ achieves its absolute minimum at the vector $v_{\ast}\in\mathbb{D}$. The value of $D_{1}(\varrho)$ is achieved at this vector, if and only if the following conditions are satisfied

[TABLE]

for all $v\in\mathbb{S}^{2}$, and

[TABLE]

for all $v\in\mathbb{S}^{2}$ such that $g_{1}(v)\leq g_{1}(v_{\ast})$.

Remark V.2

Notice that if the conditions (V.16) and (V.17) are not satisfied, then the value $g_{1}(v_{\ast})$ gives the upper bound for $D_{1}(\varrho)$ i.e.

[TABLE]

Now we determine the minimal value $g_{1}(v_{\ast})$. To this end, consider the following cases:

1. $\widetilde{\lambda}_{1}>\widetilde{\lambda}_{2}>\widetilde{\lambda}_{3}$. The solutions of (V.15) are given by

[TABLE]

and the expression for $g_{1}(v),\;v\in\mathbb{D}$ reads

[TABLE]

Note that $\widetilde{\lambda}_{2}\geq 0$. The equality (V.19) is valid in any orthonormal basis. But if we choose the basis of ordered and normalized eigenvectors of $L_{-}$, then

[TABLE]

2. $\widetilde{\lambda}_{1}=\widetilde{\lambda}_{2}>\widetilde{\lambda}_{3}$. The equation (V.15) gives $v_{1}=v_{2}=0$ and $v_{3}^{2}=1$. So (V.19) is also valid and we obtain

[TABLE]

3. $\widetilde{\lambda}_{1}=\widetilde{\lambda}_{2}=\widetilde{\lambda}_{3}$. In this case $\mathbb{D}=\mathbb{S}^{2}$ and the formula (V.19) is valid.

So we have the following

Proposition V.2

For all vectors $v\in\mathbb{D}$

[TABLE]

where $\mathrm{int}$ denotes the intermediate value and $\widetilde{\lambda}_{1},\widetilde{\lambda}_{2},\widetilde{\lambda}_{3}$ are the eigenvalues of the matrix $L_{-}$.

Remark V.3

Notice that if the Bloch vector $x$ is zero vector, then the function $g_{1}$ is constant on the set $\mathbb{D}$ and equals to $2\,\mathrm{int}\,\{\lambda_{1},\lambda_{2},\lambda_{3}\}$, where $\lambda_{1},\lambda_{2},\lambda_{2}$ are the eigenvalues of the matrix $KK^{\mathrm{T}}$.

V.3 Smooth critical points

Smooth critical points $v_{0}\in\mathbb{S}^{2}\setminus\mathbb{D}$ satisfy

[TABLE]

Applying differential analysis to the set $\mathbb{S}^{2}\setminus\mathbb{D}$ one can find all smooth critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$. Let

[TABLE]

with the constrain $||v||=1$ and where $\lambda$ is a Lagrange multiplier. For the smooth critical points one obtains

[TABLE]

where $2\omega=\lambda\,\sqrt{g_{1}^{2}-g_{2}}$. Moreover $4\,\mu=2\,f^{2}$ i.e.

[TABLE]

Remark V.4

The condition (V.21) in particular means that smooth critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ are critical points of the function $\mathbf{g}$ defined by (V.1).

To find critical points of $\mathbf{g}$ we look for the possible solutions of the equation

[TABLE]

A vector $v\in\mathbb{S}^{2}$ is a solution of (V.23) if there exists a real number $\mu$ such that $v$ is a solution of the above eigenvector problem. Applying the above results one can prove the following Theorem:

Theorem V.2

Any smooth critical point of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ satisfies (V.23) i.e. is a critical point of the function $\mathbf{g}$.

Remark V.5

One can also show that critical points contained in the set $\mathbb{D}$ satisfy (V.23), so the thesis of Theorem V.2 is true for all critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$.

The condition in equation (V.23) can be rewritten in the following way. Notice that the two-dimensional projector $\mathcal{M}$ is a critical point of the function $\mathbf{g}$ if and only if it commutes with the matrix

[TABLE]

i.e. $v$ is the solution of (V.23) if and only if

[TABLE]

On the other hand, the information about the Lagrange multiplier $\omega$ can be recovered from the equation (V.25) that leads to

[TABLE]

V.4 First application: states with maximally mixed marginals

General analysis of the equations (V.23) or (V.25) will be presented in the next section, here we consider the first application to the states with maximally mixed marginals (MMM states) i.e. such states $\varrho$ that

[TABLE]

In this case the Bloch vectors $x$ and $y$ are zero vectors. Let $K$ be the correlation matrix of MMM state $\varrho$. Using the local transformation one can bring the matrix $KK^{\mathrm{T}}$ to the diagonal form $\mathrm{diag}\,(i_{1}^{2},\,i_{2}^{2},\,i_{3}^{2})$ with $i_{1}^{2}\geq i_{2}^{2}\geq i_{3}^{2}$. Now the equation (V.15) for the set $\mathbb{D}$ has a form

[TABLE]

Moreover

[TABLE]

so the general form of $P_{v}$ satisfying (V.26) depends on the degeneracy of the eigenvalues of $G_{\mu}$.

Consider first the case when $i_{1}^{2}>i_{2}^{2}>i_{3}^{2}$. Suppose that $\mu\notin\{i_{1}^{2},\,i_{2}^{2},\,i_{3}^{2}\}$, then the eigenvalues of $G_{\mu}$ are non-degenerate and the only possible solutions of (V.26) correspond to the vectors

[TABLE]

which, according to (V.27) cannot belong to the set $\mathbb{D}$. So they can be smooth critical points of $\mathcal{M}\to||S(\mathcal{M})||_{1}$. The condition (V.22) enforces $\mu=i_{1}^{2}$ or $\mu=i_{2}^{2}$, which is in a contradiction with the assumption that $\mu\notin\{i_{1}^{2},\,i_{2}^{2},\,i_{3}^{2}\}$. So we conclude that the square of the discord can be achieved in the set $\mu\in\{i_{1}^{2},\,i_{2}^{2}\}$, since $i_{3}^{2}$ can be excluded because of the lower bound given by Theorem V.1

[TABLE]

Since $i_{2}^{2}<i_{1}^{2}$, we get

[TABLE]

Similarly one also obtains the formula (V.29) in the case when $i_{1}^{2}=i_{2}^{2}>i_{3}^{2}$ or $i_{1}^{2}>i_{2}^{2}=i_{3}^{2}$. The case $i_{1}^{2}=i_{2}^{2}=i_{3}^{2}$ is trivial: the set $\mathbb{D}$ is equal to the whole sphere $\mathbb{S}^{2}$ and the function (V.22) is constant and we again obtain (V.29). Thus we confirm the known result (see e.g. Paula )

Conclusion V.1

For any MMM state trace - norm geometric discord is equal to the intermediate value in the set of singular values of the corresponding correlation matrix.

VI Main result: determination of quantum discord in general case

In this section we use the equation (V.25) to solve the main problem, namely to find the value of geometric discord in general case. It is convenient to consider the matrix (V.24) with respect to ordered basis of orthonormal eigenvectors of the operator $L_{-}$ (we will assume that $\widetilde{\lambda}_{1}\geq\widetilde{\lambda}_{2}\geq\widetilde{\lambda}_{3}$). When we also neglect the parts which commute with the projectors $P_{v}$, we obtain the matrix

[TABLE]

where

[TABLE]

and

[TABLE]

In the following we will also use the matrices $\widetilde{\Lambda}=\Lambda-(\widetilde{\lambda}_{1}-\widetilde{\lambda}_{2})\,\openone_{3}$ and

[TABLE]

Notice that if $L_{-}$ has a degenerate spectrum, the form of the operator $\widetilde{G}_{\mu}$ is particulary simple and the desired minimal value of $\mathcal{M}\to||S(\mathcal{M})||_{1}$ can be found in a straightforward way.

VI.1 Non-degenerate spectrum of $L_{-}$

Assume that $\widetilde{\lambda}_{1}>\widetilde{\lambda}_{2}>\widetilde{\lambda}_{3}$. First we find the location of critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$. Notice that in the following all components of the vector $v$ and the Bloch vector $x$ are computed with respect to the ordered basis of orthonormal eigenvectors of $L_{-}$.

Proposition VI.1

Assume that the matrix $L_{-}=KK^{\mathrm{T}}-W_{x}$ has non - degenerate eigenvalues. The critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ ($\mathcal{M}=\openone_{3}-P_{v}$) are contained in the subset of vectors $v\in\mathbb{S}^{2}$ satisfying the condition

[TABLE]

i.e. in the big circles in $\mathbb{S}^{2}$, where the components of $v$ are given with respect to the ordered basis of eigenvectors of $L_{-}$. In particular the value of $D_{1}(\varrho)$ is achieved only on such vectors.

Proof: Suppose that critical points do not satisfy the condition $v_{1}v_{2}v_{3}=0$ and moreover let $v_{3}x_{1}-v_{1}x_{3}\neq 0$. The condition $v_{1}v_{2}v_{3}\neq 0$ implies that the critical points are outside the set $\mathbb{D}$. Notice that $\widetilde{\mu}+\langle\,{x},\,{v}\,\rangle^{2}\geq 0$. Since (V.25) is satisfied, the following strict inequality is true

[TABLE]

By a direct computation one can find the representation of $\widetilde{\mu}$

[TABLE]

along the critical points of the function $\mathbf{g}$, described by the condition

[TABLE]

or

[TABLE]

To find the critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ we use the condition (V.22), which now reads

[TABLE]

To obtain the result, let us assume that the equation (VI.1) has a non-empty set of solutions. Then we get that

[TABLE]

But this representation can not be valid outside the set $\mathbb{D}$ and we obtain the contradiction with our assumption. Let now the equation (VI.2) has the non-empty set of solutions. Then we get

[TABLE]

and after simplifications we obtain

[TABLE]

But this is impossible, as this leads to the condition $v_{2}=0$. Finally, let $v_{3}x_{1}-v_{1}x_{3}=0$. Then the equation (V.25) gives

[TABLE]

This condition together with (V.22) leads again to $v_{2}=0$. $\Box$

The assumption of non-degeneracy of the spectrum of $L_{-}$ is relevant in the above result. On the other hand, the proof suggests that the location of critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ given by Proposition VI.1 is not optimal. Indeed, if for example $x_{1}\neq 0$, the critical points cannot satisfy $v_{1}=0$ if $v_{2}v_{3}\neq 0$. Similarly, if $x_{3}\neq 0$, the critical points never lay on the big circle $v_{3}=0$ with $v_{1}v_{2}\neq 0$, unless $\widetilde{\lambda}_{1}+\widetilde{\lambda}_{3}>2\,\widetilde{\lambda}_{2}$ and the following condition is satisfied

[TABLE]

The additional informations of this kind can simplify the determination of quantum discord in some situations. Nevertheless, we will not discussing this improvements here and focus on the direct method of finding the absolute minimum of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$.

Now we can go to the solution of the main problem. This is a problem of minimizing the quadratic form $\mu$, restricted to the big circles. The calculations are elementary but tedious, so we do not present all details.

1. Let us start with $v_{1}=0$. One can find

[TABLE]

and

[TABLE]

So

[TABLE]

2. Let $v_{3}=0$. Then

[TABLE]

Inserting the minimum $\widetilde{\mu}_{\mathrm{min}}$ one gets

[TABLE]

3. Let finally $v_{2}=0$. This case is not so straightforward. One needs to introduce the additional parameter. We start with the formula

[TABLE]

To describe the absolute minimum of the function (VI.3), we need to consider auxiliary functions. The complexity of formulae below mainly follows from the presence of the absolute value in (VI.3). Let us define functions $p(\theta)$ and $r(\theta)$, as follows. The common domain of $p(\theta)$ and $r(\theta)$ is the closed interval $[-\pi,\pi]$. To define these functions we need the gap angle $\varphi_{\mathrm{gap}}\in(0,\pi/2)$, determined by the equation

[TABLE]

Moreover, let $\sigma\,:\,[-\pi,\pi]\to[-1,1]$ be any function such that

[TABLE]

and $\mathcal{O}_{\mathrm{gap}}$ be the open set

[TABLE]

Now we are ready to define

[TABLE]

and

[TABLE]

Finally, let $\theta_{\ast}$ and $\theta_{\ast\ast}$ be the angles defined below in two steps. First, let $x_{1}x_{3}\neq 0$. Define

[TABLE]

where

[TABLE]

and

[TABLE]

with

[TABLE]

If $x_{1}x_{3}=0$ we pass to the one - sided limits $x_{1}\to 0^{\pm}$ or $x_{3}\to 0^{\pm}$ (or both) in the formulas (VI.4) and (VI.5). The left and right limits can only differ by sign, but this has no importance as the functions $p(\theta)$ and $r(\theta)$ give the same result in these situations.

To obtain the upper bound for quantum discord, define

[TABLE]

and

[TABLE]

Then one can show that

[TABLE]

The above analysis leads to the following result.

Theorem VI.1

In the case of non-degenerate spectrum of the matrix $KK^{\mathrm{T}}-W_{x}$, the value of trace - norm geometric discord is given by

[TABLE]

where

[TABLE]

with $\mu_{\ast},\;\mu_{\ast\ast}$ given by (VI.6) and (VI.7) respectively.

VI.2 Degenerate spectrum of $L_{-}$

When there is any degeneracy in the spectrum of $L_{-}$, the analysis simplifies significantly. From the above discussion we obtain:

Theorem VI.2

1. Let $\widetilde{\lambda}_{1}=\widetilde{\lambda}_{2}\geq\widetilde{\lambda}_{3}$, then

[TABLE]

2. Let $\widetilde{\lambda}_{1}>\widetilde{\lambda}_{2}=\widetilde{\lambda}_{3}$, then

[TABLE]

VII Examples

In this section we present the applications of Theorems VI.1 and VI.2 to the cases already known in the literature as well as to the new examples of states. But we start with considering again, from the wider perspective, the results obtained in Sections IV and V.D.

VII.1 The class $\mathcal{E}_{0}$

In this case the correlation matrix is of the form $K=t\,V_{0}$, where $V_{0}$ is orthogonal and the Bloch vectors $x$ and $y$ are arbitrary. Thus

[TABLE]

and

[TABLE]

So we can apply Theorem VI.2 to obtain

[TABLE]

VII.2 The class of MMM states

In this case Bloch vectors $x$ and $y$ are zero vectors, so

[TABLE]

Let $\lambda_{1}\geq\lambda_{2}\geq\lambda_{3}$ be the eigenvalues of $KK^{\mathrm{T}}$. When the eigenvalues are degenerate, we apply Theorem VI.2 and obtain $D_{1}(\varrho)=\sqrt{\lambda_{2}}$. In the non-degenerate case we must apply Theorem VI.1 and the computation is more involved. Observe that in the case of MMM states

[TABLE]

where obviously $\mu_{\ast\ast}=\lambda_{2}$. To show that $\mu_{\ast}$ also equals to $\lambda_{2}$ one checks that $\theta_{\ast}=0$. Then

[TABLE]

so by (VI.6), $\mu_{\ast}=\lambda_{2}$. Thus

[TABLE]

VII.3 Pure states

It is well known that up to the local equivalence, every pure state has the form

[TABLE]

The corresponding correlation matrix and Bloch vector are given by

[TABLE]

So

[TABLE]

an we can apply Theorem VI.2 to obtain

[TABLE]

On the other hand, the measure of entanglement defined by negativity, gives the same value.

VII.4 X-states

Two - qubit X - state has the X - shaped form

[TABLE]

We can always take off - diagonal matrix elements $\varrho_{14}$ and $\varrho_{23}$ to be positive cicca . Then the correlation matrix $K$ is diagonal and the Bloch vector $x$ has only third component non-zero. So the matrix $KK^{\mathrm{T}}$ is diagonal and commutes with the operator $W_{x}$. It follows that the matrix $L_{-}$ has also diagonal form and the basis of eigenvector of $L_{-}$ coincides with the canonical basis. Since the spectrum of $L_{-}$ is in general non-degenerate, to find the value of $D_{1}(\varrho_{X})$ we apply the Theorem VI.1. Let us start with computing $\mu_{\ast}$ and $\mu_{\ast\ast}$. Notice that this case when $x_{1}x_{3}=0$ as $x_{1}=x_{2}=0$. Let

[TABLE]

then

[TABLE]

We can choose the sign of $x_{1}$ such that $\sin\theta>0$, so the angles are in the interval $[0,\pi/2]$. Computing this limit we obtain:

$\bullet$ for $\widetilde{\lambda}_{1}>\widetilde{\lambda}_{3}+||x||^{2},\;\theta_{\ast}=0$, then $N_{\ast}=\widetilde{\lambda}_{1}-\widetilde{\lambda}_{3}-||x||^{2},\;\sigma(\theta_{\ast})=\pm 1$

$\bullet$ for $\widetilde{\lambda}_{1}=\widetilde{\lambda}_{3}+||x||^{2},\;\theta_{\ast}=\pi/4$, then $N_{\ast}=0$

$\bullet$ for $\widetilde{\lambda}_{1}<\widetilde{\lambda}_{3}+||x||^{2},\;\theta_{\ast}=\pi/2$, then $N_{\ast}=\widetilde{\lambda}_{3}+||x||^{2}-\widetilde{\lambda}_{1},\;\sigma(\theta_{\ast})=\pm 1$.

Let $\theta_{\ast}=0$, then $p(0)=\cos^{2}(0\pm\varphi_{\mathrm{gap}})=\cos^{2}(\varphi_{\mathrm{gap}})$, so

[TABLE]

If $\theta_{\ast}=\pi/4$, then $\widetilde{\lambda}_{1}=\widetilde{\lambda}_{3}+||x||^{2}$ and

[TABLE]

Finally, if $\theta_{\ast}=\pi/2$, then $p(\theta_{\ast})=1$ and again

[TABLE]

Summarizing

[TABLE]

Now consider $\mu_{\ast\ast}$. Observe that $\cos\theta_{\ast\ast}=\lim\limits_{x_{1}\to 0}\cos\theta=0$. So $\theta_{\ast\ast}=\pi/2$ and

[TABLE]

Thus

[TABLE]

Next we see that $d_{1}=\widetilde{\lambda}_{1}$ and $d_{2}=\widetilde{\lambda}_{2}+||x||^{2}$. Finally one can check that

[TABLE]

and

[TABLE]

so

[TABLE]

and this formula is in the full accordance with the results presented in cicca .

VII.5 The class containing quantum-classical states

Now we apply our results to the class of states with arbitrary Bloch vector $x$ and correlation matrix having only one non-zero singular value. Such class in particular contains quantum - classical states cicca . We can always assume that

[TABLE]

One can check that possible critical points lay on the big circle $v_{2}=0$ or are given by $(0,\pm 1,0)$. By a direct inspection one gets that the minimal value can be achieved only on the circle $v_{2}=0$. So we are left with one dimensional minimization problem. This can be done by a direct analysis (see cicca ), but we apply our general result. To find value of $D_{1}(\rho)$ we need only $\mu_{\ast}$ and $\mu_{\ast\ast}$ which take the form

[TABLE]

with $\theta_{\ast}\in(\varphi_{\mathrm{gap}},\varphi_{\mathrm{gap}})$ and $\theta_{\ast\ast}\in(-\pi/2,-\varphi_{\mathrm{gap}})\cup(\varphi_{\mathrm{gap}},\pi/2)$. After a direct but quite lengthy computations one finds that the absolute minimum is given by $\mu_{\ast}$ and the value of discord reads

[TABLE]

Notice that this formula is in a full agreement with the result obtained in cicca .

VII.6 Beyond the X-states

Consider the following two-parameter family of states

[TABLE]

where for $\widetilde{\gamma}=\sqrt{1+16\gamma^{2}}$ we define

[TABLE]

One can check that (VII.3) defines a state if $|\gamma|\leq\sqrt{3}$ and $a$ is the function of $\gamma$, known only numerically. The corresponding correlation matrix and Bloch vector are as follows

[TABLE]

so the parameter $a$ is equal to the norm $||x||$ of the Bloch vector. Interesting property of the above correlation matrix is that

[TABLE]

so the matrix $L_{-}$ is simple

[TABLE]

and the corresponding ordered eigenvalues are given in the explicit way

[TABLE]

Since the spectrum of $L_{-}$ is non-degenerate, to find the value of $D_{1}(\varrho_{\gamma,a})$, we must apply the Theorem VI.1. First notice that the Bloch vector in the basis of eigenvectors of $L_{-}$ reads

[TABLE]

Then

[TABLE]

Since $N_{\ast\ast}=||x||^{2}$

[TABLE]

and

[TABLE]

So $\theta_{\ast\ast}\in(0,\pi/2)$ and $\theta_{\ast\ast}=\varphi_{\mathrm{gap}}$. It follows that $r(\theta_{\ast\ast})=1$ and

[TABLE]

Next, since $N_{\ast}=||x||^{2}\sqrt{1+16\,\gamma^{2}}$ we obtain that

[TABLE]

So applying (VI.6) we obtain

[TABLE]

On the other hand

[TABLE]

and

[TABLE]

So

[TABLE]

and

[TABLE]

VIII Summary

Let us now summarize the results reported in the article. First, we have obtained new compact and elegant formula for the trace norm of the disturbance $S(\mathcal{M})$ of the general two-qubit state. We want to emphasize, that to obtain this result we have used the qubit structure of subsystems in the crucial way. In the general qudit case we were able only to find simplified formula for the square of the modulus of disturbance and only for the special classes of states LFJ . It seems that the computation of trace norm of $S(\mathcal{M})$ in general case is a highly non-trivial or even not possible at all. Then, we have applied the formula for $||S(\mathcal{M})||_{1}$ to compute the value of geometric discord in the situation when the minimization procedure is not needed at all or is straightforward. This is the case of the class of states with arbitrary Bloch vectors and the correlation matrix proportional to some orthogonal matrix in three dimensions i.e. $K=t\,V_{0}$. In particular, for such states with zero Bloch vectors, $||S(\mathcal{M})||_{1}=|t|$ and $D_{1}(\varrho)=|t|$. The presence of non-zero Bloch vectors only slightly changes the trace norm and the minimization can be done in elementary way. This again must be contrasted with the general qudit case, where the similar results are valid only for vanishing Bloch vectors and the classes of Werner and isotropic states LFJ . As a first step towards the computation of $D_{1}(\varrho)$ in general case, we considered critical points of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$. To solve this problem, we introduce the function $\mathbf{g}$ defined on the unit sphere $\mathbb{S}^{2}\subset\mathbb{R}^{3}$ with values in $\mathbb{R}^{2}$. Any critical point of the mapping $\mathcal{M}\to||S(\mathcal{M})||_{1}$ is a critical point of the function $\mathbf{g}$ and the properties of this function can be used to isolate singular and smooth critical points. It is remarkable that the information about critical points is encoded in the spectrum of the matrix $L_{-}=KK^{\mathrm{T}}-W_{x}$. One of the most important results of our study says that when this spectrum is non-degenerate, the critical points are located on big circles in $\mathbb{S}^{2}$. It means that the problem of finding geometric discord reduces to minimization of the quadratic form restricted to big circles. The complete solution of this problem is rather involved and it simplifies significantly, when any degeneracy of the spectrum of $L_{-}$ is present.

We apply the Theorems VI.1 and VI.2 to the cases already discussed in the literature (Bell diagonal, X-shaped and quantum-classical states) as well as new examples beyond the class of X-states. It is also worth to emphasize that in particular cases, methods based on partial results can be more efficient then the direct application of the general result, what is explicitly shown by the discussions of MMM states in Sections V.D and VII.B.

Having the complete answer for the two-qubit system we are able to estimate possible obstacles in finding general explicit formula for the trace-distance quantum discord for qudits. As we have already stressed, the geometry of the two-qubit quantum state space is exceptional and the set of geometrical tools which is at our disposal is very reach. Let us recall here the most effective one, that any correlation matrix can be diagonalized by means of local transformations, due to the accordance of the dimensions of relevant spaces and groups i.e. $\mathrm{SO(3)}\simeq\mathrm{Adj(SU(2))}$. For the next higher level system of two-qutrits $d=3$ we have that $\mathrm{Adj(SU(3)))}$ is a small subgroup of $\mathrm{SO(8)}$, too small to allow diagonalization of a generic two-qutrit correlation matrix by means of local transformations. For higher dimensional cases this effects only aggravates, what means that using local orbit of the set of diagonal correlation matrices we can cover only a part of the quantum state space. The solution of the minimization problem for some families of two-qudit states is given in our previous works LFJ ; jfl-pla .

The strict analytical general formulas characterising quantum correlations valid for all or sufficiently wide class of quantum states of compound systems, in general are not known and only selected cases are described in the literature. Such explicit knowledge would be of value that hardly can be overestimated, as it gives insight into the structure of the quantum states, but on the other hand due to the complexity of this space, computation of physically important measures of quantum correlations is so hard task. We hope that results obtained in the present paper will fill the gap in the knowledge on measures of quantum correlations of compound systems for the most fundamental physical system of two qubits.