# Measuring geometric quantum discord using one bit of quantum information

###### Abstract

We describe an efficient DQC1-algorithm to quantify the amount of Geometric Quantum Discord present in the output state of a DQC1 computation. DQC1 is a model of computation that utilizes separable states to solve a problem with no known efficient classical algorithm and is known to contain quantum correlations as measured by the discord. For the general case of a -qubit DQC1-state we provide an analytical expression for the Geometric Quantum Discord and find that its typical (and maximum) value decreases exponentially with . This is in contrast to the standard Quantum Discord whose value for typical DQC1-states is known to be independent of . We experimentally demonstrate the proposed algorithm on a four-qubit liquid-state nuclear magnetic resonance quantum information processor. In the special case of a two-qubit DQC1 model, we also provide an expression for the Quantum Discord that only requires the outcome of the DQC1 algorithm.

###### pacs:

03.67.-a, 03.67.Lx, 76.60.-k.Since its inception in 2001, Quantum Discord (QD) Henderson and Vedral (2001); Ollivier and Zurek (2001) and related measures have been used to quantify the amount of non-classical correlations in a physical system. In addition to quantifying the most well known quantum correlations of entanglement, they also capture the quantum correlations that exist in separable states. Discord measures are based on the premise that if a measurement on one part of a bipartite state disturbs the total state, there must be correlations stronger than what is found in the classical world.

Quantum discord has been studied extensively over the past few years, from operational definitions Cavalcanti et al. (2011); *Madhok:2011bs; *zurek2003quantum, to witnesses Bylicka and Chruscinski (2010); *Laine:2010mi; *Ma:2011qa; *Maziero:2010kl; *rahimi2009single; *Smirne:2011pi; *Yu:2011fu; Dakic et al. (2010); Passante et al. (2011) and analytical expressions for specific sets of states Ali (2010); Ali et al. (2010); Chen et al. (2011); Girolami and Adesso (2011a); Lu et al. (2011); Luo (2008). However, it is generally difficult to calculate as it requires both full state knowledge and an optimization over all projective measurements of a subsystem. Even in the simple case of two qubits, a general closed form expression does not exist. The measure of Geometric Quantum Discord (GQD) was introduced Dakic et al. (2010) as a simple geometric measure of the distance from a given state to the closest classical (zero discord) state.

DQC1 is a model of mixed state quantum computation Knill and Laflamme (1998) that contains limited entanglement Datta et al. (2005) yet is thought to outperform classical methods. It does generate non-classical correlations as measured by the quantum discord Datta et al. (2008); Lanyon et al. (2008); Passante et al. (2011), which, in hind sight, is not surprising since almost all quantum states have non-zero discord Ferraro et al. (2010). Quantum discord has indeed been witnessed in the final states of DQC1 computations Lanyon et al. (2008), even at very small values in highly mixed quantum states Passante et al. (2011). To further explore the role of quantum correlations in the DQC1 model and beyond, their quantification becomes a necessary pursuit.

In this work, we derive an analytical expression for the geometric discord in the final state of a DQC1 computation of arbitrary dimension, in terms of quantities that can be efficiently estimated on a DQC1 computer. We then demonstrate its experimental evaluation in a four-qubit implementation of a DQC1 algorithm. Additionally, we derive a simple expression for the quantum discord in the special case of a two-qubit DQC1-state.

## I Introduction

Quantum discord is the most well-known measure of non-classical correlations and is thought to differentiate quantum and classical systems Henderson and Vedral (2001); Ollivier and Zurek (2001). It is defined as the minimum difference between two classically equivalent formulations of the mutual information between subsystems and , and , where is the Shannon entropy when is described by a classical probability distribution, and is the von Neumann entropy when describes a quantum system. In the quantum case the conditional entropy depends on the measurement basis, , where is a complete set of orthonormal projectors on such that , is the probability of observing outcome on system , and is the state of system conditional on the measurement of system returning measurement outcome . The Quantum Discord, , thus reduces to

(1) |

where is the reduced density matrix of system . In order to calculate this quantity, full state knowledge is required in addition to a minimization over all possible projective measurements on subsystem .

While there is no analytical expression for the quantum discord of a general state, the discord for various sets of two-qubit states have been found Luo (2008); Ali et al. (2010); Chen et al. (2011); Lu et al. (2011). In addition, the quantum discord for all two-qubit states has been reduced to the problem of solving a set of transcendental equations Girolami and Adesso (2011a). One of the results we report here is an expression for the QD in the final state of a two-qubit DQC1 algorithm that can be calculated with only the outcome of the DQC1 algorithm (trace of the unitary) and the initial polarization of the top register.

Motivated by the difficultly in computing the QD, Dakic, Vedral and Brukner Dakic et al. (2010) proposed the GQD, which is defined as

(2) |

where is the set of all zero discord states (.) These can be written as . The quantity is the square of the Hilbert-Schmidt norm of Hermitian operators. Note that the GQD, like the original quantum discord, is not symmetric in the subsystems, and is defined as where .

Shortly after the proposal of this measure of quantum correlations, Luo and Fu Luo and Fu (2010) showed that Equation (2) is equivalent to the minimization

(3) |

where is a projective measurement on system , and . The minimization can be performed analytically for arbitrary -dimensional Dakic et al. (2010) and -dimensional Vinjanampathy and Rau (2011) states. The resulting expressions are in terms of quantities that are not efficiently experimentally accessible. However, a method for measuring a tight lower bound on the GQD has been recently proposed Girolami and Adesso (2011b), requiring only a constant number of measurements on up to four copies of the state. In this report we provide an analytical expression for the GQD of -dimensional DQC1-states, and describe a DQC1-algorithm to efficiently estimate it.

The DQC1 model of computation is an example of a mixed state quantum information processor Knill and Laflamme (1998). As shown in Figure 1, this model has access to only one qubit with non-zero polarization, accompanied by a register of maximally mixed qubits. The computation consists of two gates: a Hadamard on the top, single qubit register, followed by an -qubit unitary on the bottom register that is controlled by the state of the top qubit. The final state of this simple algorithm is what we call the DQC1-state and can be written as

Expectation value measurements of the Pauli matrices and yield the real and imaginary values of the trace of the unitary: and . Since calculating the trace of a unitary does not have a known efficient classical algorithm, a DQC1 computer is believed to be more powerful than its classical counterpart.

Despite its apparent improvement over classical methods, DQC1 is known to have zero bipartite entanglement between the top and bottom registers and can only have a small amount of entanglement across any other bipartite splitting Datta et al. (2005). It can, however, have quantum correlations characterized by the quantum discord Datta et al. (2008). This has been confirmed in two experiments to date: measurement of the discord in a two-qubit optics setup using full state tomography Lanyon et al. (2008), and a four-qubit NMR implementation that witnesses discord with a small number of experiments Passante et al. (2011).

## Ii Analytical expression for GQD

In order to find an expression for the geometric discord in a DQC1-state, let us write Eqn. (3) as

(4) |

The first term, the purity of the state of the total system, is invariant under unitary transformations. Therefore it can be calculated for the initial state, and only depends on the initial polarization of the top register:

(5) |

In order to calculate the other two terms, consider parameterizing the measurement on the top qubit as , where and ; and are real and ; and . Thus, the state after measurement is

(6) |

For the DQC1-state, the terms and are equivalent, and evaluate to

where are the eigenphases of (i.e. ), and we have used . Thus, we are left with the task of minimizing

over the measurement parameters and . Examining first and second order derivatives of with respect to and , we find the optimal measurement parameters to be

(7) | |||||

(8) | |||||

(9) |

which reduces the GQD to

(10) |

where , and can be evaluated with a DQC1-algorithm using back-to-back applications of the control-. The depth of the circuit that evaluates is at most double the one that evaluates , and therefore has the same efficiency of evaluation as the standard DQC1 algorithm. Worth noting here that indicates that the optimal measurement is always in the transverse (X-Y) plane of the Bloch sphere, independent of the implemented unitary. Also note that Eqn. (9) allows for the ability to experimentally evaluate the optimal measurement parameter using the same data collected for the evaluation of .

The GQD does not depend on the eigenvectors of the unitary, but rather the distribution of eigenphases. This gives rise to classes of unitaries that generate the same GQD in a DQC1 circuit. For instance, reproducing a result from Ref. Dakic et al. (2010), the set of unitaries that generate zero GQD (and hence, zero QD) will have , implying for some . This is possible if and only if the eigenphases of are . That is to say , where is a binary observable ().

For a unitary drawn randomly according to the Haar measure, its eigenphases are randomly distributed over the unit circle Diaconis (2003). Thus, for large , approaches zero and approaches its maximum value of , i.e. the GQD of a DQC1-state decreases exponentially with the number of qubits. This is in contrast to the average QD of a DQC1-state, which was shown in Ref. Datta et al. (2008) to be independent of the number of qubits for large .

## Iii Experimental measurement of the geometric discord

We measure the quantum correlations, as quantified by the GQD, in a liquid-state nuclear magnetic resonance implementation of the DQC1 algorithm. The qubits in NMR are ensembles of spin-1/2 nuclei. In this experiment, we use the four carbon-13 nuclei in the molecule trans-crotonic acid (molecular information and Hamiltonian parameters can be found in Ref. Passante et al. (2009).) The experiment is implemented in a 16.7 T magnetic field Bruker Avance spectrometer, where, for carbon-13 nuclei () at room temperature, the initial polarization of the thermal state is . The spins are manipulated with radio frequency pulses numerically generated using the GRAPE algorithm Khaneja et al. (2005); Ryan et al. (2008). The pulses are designed to have a fidelity of no less than 0.998 and are adjusted for non-linearities in the pulse generation and transmission by placing a pickup coil at the location of the sample and running a feedback loop to iteratively adjust the pulse shape for optimal transmission.

For the particular instance of the DQC1 circuit in this experiment, we choose to implement one of the unitary matrices used for approximating the Jones polynomial for braids with four strands Passante et al. (2009). This is a problem that completely encapsulates the power of the DQC1 model and is of interest in many fields of physics and math. The unitary has the form , where and . The DQC1-state, after the application of the control-, is known to have nonzero discord as witnessed in a recent experiment Passante et al. (2011). Using Eqn. (10), we calculate the expected value of the GQD to be . As shown in Figure 2, this is achieved for an optimal measurement with and rad. Also shown in Figure 2(ii) is the experimental value of GQD in the DQC1-state, which is found to be by measuring the outcome of a DQC1 circuit with back-to-back applications of control-. For our experimental polarization (), the GQD in the final state is .

## Iv The quantum discord of two-qubit DQC1-states

In the case where the bottom register is a single qubit, we report an analytical expression for the QD of the DQC1-state, thereby contributing to the sets of two-qubit states for which the QD has been analytically solved. Commencing from Eqn. (1), a minimization must be performed over the conditional entropy term, , which, in the case of a two-qubit DQC1-state, reduces to

where , , and and parameterize the measurement on the top qubit (see the GQD case above). The minimization with respect to the measurement parameters and results in an optimal measurement characterized by

This allows us to write the analytical expression for the QD of a two-qubit DQC1-state as

(11) |

where is the binary entropy, and , which can be measured directly with a DQC1-algorithm.

Hence, the quantum discord of a two-qubit DQC1-state is zero for or , and maximum for . Recall, from Eqn. (10), that the GQD is zero for and maximum for . Consider a general single-qubit unitary, expressed as , where is a rotation by about the unit vector in three dimensions , are the Pauli matrices, and a global phase. It follows from that and . Thus, the QD and GQD are simultaneously maximum for , and simultaneously zero for , independent of the rotation axis . Examples of unitaries that can be written as rotations include bit-flip (or NOT) and Hadamard.

## V Conclusion

In this article we provided analytical expressions for the quantum correlations present in a DQC1-state in terms of experimentally accessible quantities. The geometric quantum discord can be computed for -dimensional DQC1-states by implementing a DQC1 algorithm with a second application of the controlled unitary. The algorithm was experimentally demonstrated for a four-qubit liquid-state NMR implementation of the DQC1 model. We also showed, that for the special case of a two-qubit system, the quantum discord can be calculated using outcome of the DQC1-circuit (the trace of the unitary).

For a typical DQC1-state, where the unitary is chosen uniformly at random according to the Haar measure, we found that the geometric measure of discord scales with , in contrast to the quantum discord, which is known to be independent of (for large ). This suggests that the geometric quantum discord fails to completely quantify the correlations defined in the original entropic measure of quantum discord. Finally, the work we have presented here further supports the suspicion that the apparent speedup exhibited by (the dynamics of) DQC1 is not necessarily captured by the geometric measure of quantum discord at the conclusion of the algorithm.

###### Acknowledgements.

G.P. would like to thank M. Ditty for his technical expertise with the spectrometer. This work was funded by NSERC, QuantumWorks, and CIFAR.## References

- Henderson and Vedral (2001) L. Henderson and V. Vedral, J. Phys. A: Math. Gen. 34, 6899 (2001).
- Ollivier and Zurek (2001) H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 (2001).
- Cavalcanti et al. (2011) D. Cavalcanti, L. Aolita, S. Boixo, K. Modi, M. Piani, and A. Winter, Phys. Rev. A 83, 032324 (2011).
- Madhok and Datta (2011) V. Madhok and A. Datta, Phys. Rev. A 83, 032323 (2011).
- Zurek (2003) W.H. Zurek, Physical Review A 67, 012320 (2003).
- Bylicka and Chruscinski (2010) B. Bylicka and D. Chruscinski, Phys. Rev. A. 81, 062102 (2010).
- Laine et al. (2010) E.-M. Laine, J. Pillo, and H.-P. Breuer, EPL 92, 60010 (2010).
- Ma et al. (2011) Z.-H. Ma, Z.-H. Chen, and J.-L. Chen, Arxiv preprint arXiv:1104.0299 (2011).
- Maziero and Serra (2010) J. Maziero and R. Serra, Arxiv preprint arXiv:1012.3075 (2010).
- Rahimi and SaiToh (2010) R. Rahimi and A. SaiToh, Phys. Rev. A. 82, 022314 (2010).
- Smirne et al. (2011) A. Smirne, D. Brivio, S. Cialdi, B. Vacchini, and M. Paris, Arxiv preprint arXiv:1105.0174 (2011).
- Yu et al. (2011) S. Yu, C. Zhang, Q. Chen, and C. Oh, Arxiv preprint arXiv:1102.4710 (2011).
- Dakic et al. (2010) B. Dakic, V. Vedral, and C. Brukner, Phys. Rev. Lett. 105, 190502 (2010).
- Passante et al. (2011) G. Passante, O. Moussa, D.A. Trottier, and R. Laflamme, Phys. Rev. A 84, 044302 (2011).
- Ali (2010) M. Ali, J. Phys. A: Math. Theor. 43, 495303 (2010).
- Ali et al. (2010) M. Ali, A.R.P. Rau, and G. Alber, Phys. Rev. A. 81, 042105 (2010).
- Chen et al. (2011) Q. Chen, C. Zhang, S. Yu, X. Yi, and C. Oh, Arxiv preprint arXiv:1102.0181 (2011).
- Girolami and Adesso (2011a) D. Girolami and G. Adesso, Phys. Rev. A. 83, 052108 (2011a).
- Lu et al. (2011) X.-M. Lu, J. Ma, Z. Xi, and X. Wang, Phys. Rev. A. 83, 012327 (2011).
- Luo (2008) S. Luo, Phys. Rev. A. 77, 042303 (2008).
- Knill and Laflamme (1998) E. Knill and R. Laflamme, Phys. Rev. Lett. 81, 5672 (1998).
- Datta et al. (2005) A. Datta, S. T. Flammia, and C. M. Caves, Phys. Rev. A 72, 042316 (2005).
- Datta et al. (2008) A. Datta, A. Shaji, and C. M. Caves, Phys. Rev. Lett. 100, 050502 (2008), arXiv:0709.0548.
- Lanyon et al. (2008) B. P. Lanyon, M. Barbieri, M. P. Almeida, and A. G. White, Phys. Rev. Lett. 101, 200501 (2008), arXiv:quant-ph/0807.0668.
- Ferraro et al. (2010) A. Ferraro, L. Aolita, D. Cavalcanti, F.M. Cucchietti, and A. Acín, Phys. Rev. A 81, 052318 (2010).
- Luo and Fu (2010) S. Luo and S. Fu, Phys. Rev. A. 82, 034302 (2010).
- Vinjanampathy and Rau (2011) S. Vinjanampathy and A. Rau, “Calculation of quantum discord for qubit-qudit or n qubits,” (2011), arXiv:1106.4488.
- Girolami and Adesso (2011b) D. Girolami and G. Adesso, “Observable measure of quantum correlations,” (2011b), arXiv:1110.5083.
- Diaconis (2003) R. Diaconis, Bull. Amer. Math. Soc. 40, 155 (2003).
- Passante et al. (2009) G. Passante, O. Moussa, C.A. Ryan, and R. Laflamme, Physical Review Letters 103, 250501 (2009).
- Khaneja et al. (2005) N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbruggen, and S. Glaser, J. Magn. Reson. 172, 296 (2005).
- Ryan et al. (2008) C. A. Ryan, C. Negrevergne, M. Laforest, E. Knill, and R. Laflamme, Phys. Rev. A 78, 012328 (2008).