Research Article | Open Access
Computing Coherence Vectors and Correlation Matrices with Application to Quantum Discord Quantification
Coherence vectors and correlation matrices are important functions frequently used in physics. The numerical calculation of these functions directly from their definitions, which involves Kronecker products and matrix multiplications, may seem to be a reasonable option. Notwithstanding, as we demonstrate in this paper, some algebraic manipulations before programming can reduce considerably their computational complexity. Besides, we provide Fortran code to generate generalized Gell-Mann matrices and to compute the optimized and unoptimized versions of associated Bloch’s vectors and correlation matrix in the case of bipartite quantum systems. As a code test and application example, we consider the calculation of Hilbert-Schmidt quantum discords.
Correlation functions are fundamental objects for statistical analysis and are thus ubiquitous in most kinds of scientific inquiries and their applications [1, 2]. In physics, correlation functions have an important role for research in areas such as quantum optics and open systems [3, 4], phase transitions and condensed matter physics [5, 6], and quantum field theory and nuclear and particle physics . Another area in which correlation functions are omnipresent is quantum information science (QIS), an interdisciplinary field that extends the applicability of the classical theories of information, computation, and computational complexity [8–10].
Investigations about the quantum correlations in physical systems have been one of the main catalyzers for developments in QIS [11–15]. There are several guises of quantum correlations, and quantum discord stands among the most promising quantum resources for fueling the quantum advantage [16–26]. When computing or witnessing quantum discord, or other kinds of correlation or quantumness quantifiers, we are frequently faced with the need for calculating coherence vectors and correlation matrices [27–41]. And it is the main aim of this paper to provide formulas for these functions which are amenable for more efficient numerical calculations when compared with the direct implementation of their definitions.
In order to define coherence vectors and correlation matrices, let us consider a composite bipartite system with Hilbert space . Hereafter, the corresponding dimensions are denoted by for . In addition, let , withbeing a basis for the special unitary group . Any density operator describing the state of the system can be written in the local basis as follows:whereand is the identity operator in . One can readily verify that the components of the coherence (or Bloch’s) vectors and and of the correlation matrix are given byIt is worth mentioning that the mean value of any observable in , for , can be obtained using these quantities.
In https://github.com/jonasmaziero/LibForQ.git, we provide Fortran code to compute the coherence vectors, correlation matrices, and quantum discord quantifiers we deal with here. Besides these functions, there are other tools therein that may be of interest to the reader. The instructions on how to use the software are provided in the readme file. Related to the content of this section, the subroutine bloch_vector_gellmann_unopt(, , ) returns the coherence vectors or and the subroutine corrmat_gellmann_unopt(, , , ) computes the correlation matrix . Now, let us notice that if calculated directly from the equations above, for , the computational complexity (CC) to obtain the coherence vectors and or the correlation matrix is
The remainder of this paper is structured as follows. In Section 2, we obtain formulas for , , and which are amenable for more efficient numerical computations. In Section 3, we test these formulas by applying them in the calculation of Hilbert-Schmidt quantum discords. In Section 4, we make some final remarks about the usefulness and possible applications of the results reported here.
2. Computing Coherence Vectors and Correlation Matrices
The partial trace function  can be used in order to obtain the reduced states and and to write the components of Bloch’s vectors in the formThus, when computing the coherence vectors of the parties and , we shall solve a similar problem; so let us consider it separately. That is to say, we shall regard a generic density operator written aswhere .
Now, and for the remainder of this paper, we assume that the matrix elements of regarded density operator in the standard computational basis are given. We want to compute Bloch’s vector : . For the sake of doing that, a particular basis must be chosen. Here we pick the generalized Gell-Mann matrices, which are presented as follows in three groups :which are named the diagonal, symmetric, and antisymmetric groups, respectively. The last two groups possess generators each. Any one of these matrices can be obtained by calling the subroutine gellmann(, , , , ). For the first group, , we make and, in this case, one can set to any integer.
It is straightforward seeing that, for the generators above, the corresponding components of Bloch’s vector can be expressed directly in terms of the matrix elements of the density operator as follows:
These expressions were implemented in the Fortran subroutine bloch_vector_gellmann(, , ). With this subroutine and the partial trace function , we can compute the coherence vectors and .
We observe that after these simple algebraic manipulations the computational complexity of Bloch’s vector turns out to be basically the CC for the partial trace function. Hence, from , we have that, for ,
One detail we should keep in mind when making use of the codes linked to this paper is the convention we apply to the indexes of the components of . For the first group of generators, , naturally, . We continue with the second group of generators, , by setting , , . The same convention is used for the third group of generators, , but here we begin with .
Next we address the computation of the correlation matrix , which is a matrix that we write in the following form:with the submatrices given as shown below. For convenience, we define the auxiliary variables:
The matrix elements of , whose dimension is , correspond to the diagonal generators for and diagonal generators for :
The matrix elements of , whose dimension is , correspond to the diagonal generators for and symmetric generators for :
The matrix elements of , whose dimension is , correspond to the diagonal generators for and antisymmetric generators for :
The matrix elements of , whose dimension is , correspond to the symmetric generators for and diagonal generators for :
The matrix elements of , whose dimension is , correspond to the symmetric generators for and symmetric generators for :
The matrix elements of , whose dimension is , correspond to the symmetric generators for and antisymmetric generators for :
The matrix elements of , whose dimension is , correspond to the antisymmetric generators for and diagonal generators for :
The matrix elements of , whose dimension is , correspond to the antisymmetric generators for and symmetric generators for :
The matrix elements of , whose dimension is , correspond to antisymmetric generators for and antisymmetric generators for :
We remark that when implementing these expressions numerically, for the sake of mapping the local to the global computational basis, we utilize, for example,
The subroutine corrmat_gellmann(, , , ) returns the correlation matrix , as written in (11), associated with the bipartite density operator and computed using the Gell-Mann basis, as described in this section. The convention for the indexes of the matrix elements is defined in the same way as for the coherence vectors. The computational complexity for , computed via the optimized expressions obtained in this section, is, for ,By generating some random density matrices , we checked that the expressions and the corresponding code for the unoptimized and optimized versions of , , and agree. Additional tests shall be presented in the next section, where we calculate some quantum discord quantifiers.
3. Computing Hilbert-Schmidt Quantum Discords
The calculation of quantum discord functions (QD) usually involves hard optimization problems [46, 47]. In the last few years, a large amount of effort has been dedicated towards computing QD analytically, with success being obtained mostly for low-dimensional quantum systems [48–65]. Although not meeting all the required properties for a bona fide QD quantifier , the Hilbert-Schmidt discord (HSD) ,is drawing much attention due to its amenability for analytical computations when compared with most other QD measures. In the last equation, the minimization is performed over the classical-quantum states:with being a probability distribution, being an orthonormal basis for , being generic density operators defined in , and being the Hilbert-Schmidt norm of the linear operator , with being the transpose conjugate of .
In this paper, as a basic test for the Fortran code provided to obtain coherence vectors and correlation matrices, we shall compute the following lower bound for the HSD :where are the eigenvalues, sorted in nonincreasing order, of the matrix:In the above equation, stands for the transpose of a vector or matrix. We observe that the other version of the HSD, , can be obtained from the above equations simply by exchanging and and using instead of .
It is interesting that, as was proved in , a bipartite state , with polarization vectors and and correlation matrix , is classical-quantum if and only if there exists a -dimensional projector in the space such thatBased on this fact, an ameliorated version for the Hilbert-Schmidt quantum discord (AHSD) was proposed : with the matrix defined aswhere and are arbitrary functions of . Then, by setting and using the purity,to address the problem of noncontractivity of the Hilbert-Schmidt distance, the following analytical formula was presented :
Thus both discord quantifiers and are, at the end of the day, obtained from the eigenvalues . And the computation of these eigenvalues requires the knowledge of the coherence vector (or ) and of the correlation matrix . These QD measures were implemented in the Fortran functions discord_hs(ssys, , , ) and discord_hsa(ssys, , , ), where ssys = ”s”, with , specifies which version of the quantum discord is to be computed. As an example, let us use the formulas provided in this paper and the associated code to compute the HSD and AHSD of Werner states in (with ): where and The reduced states of are , whose purity is . The results for the HSD and AHSD of are presented in Figure 1.
4. Concluding Remarks
In this paper, we addressed the problem of computing coherence vectors and correlations matrices. We obtained formulas for these functions which make possible considerably more efficient numerical implementation when compared with the direct use of their definitions. We provided Fortran code to calculate all the quantities regarded in this paper. As a test for our formulas and code, we computed Hilbert-Schmidt quantum discords of Werner states. It is important to observe that although our focus here was in quantum information science, the tools provided can find application in other areas, such as the calculation of order parameters and correlations functions for the study of phase transitions in discrete quantum or classical systems.
The author declares that there are no competing interests regarding the publication of this paper.
This work was supported by the following Brazilian funding agencies: Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Processes 441875/2014-9 and 303496/2014-2), Instituto Nacional de Ciência e Tecnologia de Informação Quântica (INCT-IQ) (Process 2008/57856-6), and Coordenação de Desenvolvimento de Pessoal de Nível Superior (CAPES) (Process 6531/2014-08). The author gratefully acknowledges the hospitality of the Physics Institute and Laser Spectroscopy Group at Universidad de la República, Uruguay.
- E. T. Jaynes, Probability Theory, the Logic of Science, Cambridge University Press, New York, NY, USA, 2003.
- M. A. Carlton and J. L. Devore, Probability with Applications in Engineering, Science, and Technology, Springer Texts in Statistics, Springer, New York, NY, USA, 2014.
- C. Gerry and P. Knight, Introductory Quantum Optics, Cambridge University Press, New York, NY, USA, 2005.
- H.-P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press, New York, NY, USA, 2002.
- S. Sachdev, Quantum Phase Transitions, Cambridge University Press, New York, NY, USA, 2011.
- J. Hutchinson, J. P. Keating, and F. Mezzadri, “On relations between one-dimensional quantum and two-dimensional classical spin systems,” Advances in Mathematical Physics, vol. 2015, Article ID 652026, 18 pages, 2015.
- L. Álvarez-Gaumé and M. Vázquez-Mozo, An Invitation to Quantum Field Theory, Springer, Berlin, Germany, 2012.
- M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, Cambridge, UK, 2000.
- J. Preskill, Quantum Information and Computation, http://theory.caltech.edu/people/preskill/ph229/.
- M. M. Wilde, Quantum Information Theory, Cambridge University Press, Cambridge, UK, 2013.
- M. D. Reid, P. D. Drummond, W. P. Bowen et al., “The Einstein-Podolsky-Rosen paradox: from concepts to applications,” Reviews of Modern Physics, vol. 81, no. 4, pp. 1727–1751, 2009.
- R. Horodecki, P. Horodecki, and M. H. Horodecki, “Quantum entanglement,” Reviews of Modern Physics, vol. 81, no. 2, pp. 865–942, 2009.
- M. Li, S.-M. Fei, and X. Li-Jost, “Quantum entanglement: separability, measure, fidelity of teleportation, and distillation,” Advances in Mathematical Physics, vol. 2010, Article ID 301072, 57 pages, 2010.
- L. C. Céleri, J. Maziero, and R. M. Serra, “Theoretical and experimental aspects of quantum discord and related measures,” International Journal of Quantum Information, vol. 9, no. 7-8, pp. 1837–1873, 2011.
- K. Modi, A. Brodutch, H. Cable, T. Paterek, and V. Vedral, “The classical-quantum boundary for correlations: discord and related measures,” Reviews of Modern Physics, vol. 84, no. 4, article 1655, 2012.
- D. O. Soares-Pinto, R. Auccaise, J. Maziero, A. Gavini-Viana, R. M. Serra, and L. C. Céleri, “On the quantumness of correlations in nuclear magnetic resonance,” Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 370, no. 1976, pp. 4821–4836, 2012.
- M. Ávila, G. H. Sun, and A. L. Salas-Brito, “Scales of time where the quantum discord allows an efficient execution of the DQC1 algorithm,” Advances in Mathematical Physics, vol. 2014, Article ID 367905, 4 pages, 2014.
- A. Streltsov, Quantum Correlations Beyond Entanglement, and Their Role in Quantum Information Theory, Springer Briefs in Physics, 2015.
- A. Datta, A. Shaji, and C. M. Caves, “Quantum discord and the power of one qubit,” Physical Review Letters, vol. 100, no. 5, Article ID 050502, 2008.
- D. Cavalcanti, L. Aolita, S. Boixo, K. Modi, M. Piani, and A. Winter, “Operational interpretations of quantum discord,” Physical Review A, vol. 83, Article ID 032324, 2011.
- S. Pirandola, “Quantum discord as a resource for quantum cryptography,” Scientific Reports, vol. 4, article 6956, 2014.
- L. Roa, J. C. Retamal, and M. Alid-Vaccarezza, “Dissonance is required for assisted optimal state discrimination,” Physical Review Letters, vol. 107, no. 8, Article ID 080401, 2011.
- M. Piani, V. Narasimhachar, and J. Calsamiglia, “Quantumness of correlations, quantumness of ensembles and quantum data hiding,” New Journal of Physics, vol. 16, Article ID 113001, 2014.
- B. Dakić, Y. O. Lipp, X. Ma et al., “Quantum discord as resource for remote state preparation,” Nature Physics, vol. 8, no. 9, pp. 666–670, 2012.
- D. Girolami, A. M. Souza, V. Giovannetti et al., “Quantum discord determines the interferometric power of quantum states,” Physical Review Letters, vol. 112, no. 21, Article ID 210401, 2014.
- T. K. Chuan, J. Maillard, K. Modi, T. Paterek, M. Paternostro, and M. Piani, “Quantum discord bounds the amount of distributed entanglement,” Physical Review Letters, vol. 109, no. 7, Article ID 070501, 2012.
- J. Maziero and R. M. Serra, “Classicality witness for two-qubit states,” International Journal of Quantum Information, vol. 10, no. 3, Article ID 1250028, 8 pages, 2012.
- R. Auccaise, J. Maziero, L. C. Céleri et al., “Experimentally witnessing the quantumness of correlations,” Physical Review Letters, vol. 107, no. 7, Article ID 070501, 2011.
- F. M. Paula, I. A. Silva, J. D. Montealegre et al., “Observation of environment-induced double sudden transitions in geometric quantum correlations,” Physical Review Letters, vol. 111, no. 25, Article ID 250401, 2013.
- I. A. Silva, D. Girolami, R. Auccaise et al., “Measuring bipartite quantum correlations of an unknown state,” Physical Review Letters, vol. 110, no. 14, Article ID 140501, 2013.
- G. H. Aguilar, O. Jiménez Farías, J. Maziero, R. M. Serra, P. H. Souto Ribeiro, and S. P. Walborn, “Experimental estimate of a classicality witness via a single measurement,” Physical Review Letters, vol. 108, no. 6, Article ID 063601, 2012.
- A. Bayat and S. Bose, “Entanglement transfer through an antiferromagnetic spin chain,” Advances in Mathematical Physics, vol. 2010, Article ID 127182, 11 pages, 2010.
- M. Cianciaruso, T. R. Bromley, W. Roga, R. Lo Franco, and G. Adesso, “Universal freezing of quantum correlations within the geometric approach,” Scientific Reports, vol. 5, Article ID 10177, 2015.
- T. R. Bromley, M. Cianciaruso, R. Lo Franco, and G. Adesso, “Unifying approach to the quantification of bipartite correlations by Bures distance,” Journal of Physics A: Mathematical and Theoretical, vol. 47, no. 40, Article ID 405302, 20 pages, 2014.
- B. Aaronson, R. L. Franco, G. Compagno, and G. Adesso, “Hierarchy and dynamics of trace distance correlations,” New Journal of Physics, vol. 15, Article ID 093022, 2013.
- B. Aaronson, R. Lo Franco, and G. Adesso, “Comparative investigation of the freezing phenomena for quantum correlations under nondissipative decoherence,” Physical Review A, vol. 88, no. 1, Article ID 012120, 2013.
- B. Bellomo, G. L. Giorgi, F. Galve, R. Lo Franco, G. Compagno, and R. Zambrini, “Unified view of correlations using the square-norm distance,” Physical Review A—Atomic, Molecular, and Optical Physics, vol. 85, no. 3, Article ID 032104, 2012.
- B. Bellomo, R. Lo Franco, and G. Compagno, “Dynamics of geometric and entropic quantifiers of correlations in open quantum systems,” Physical Review A, vol. 86, no. 1, Article ID 012312, 2012.
- R. Lo Franco, B. Bellomo, E. Andersson, and G. Compagno, “Revival of quantum correlations without system-environment back-action,” Physical Review A, vol. 85, no. 3, Article ID 032318, 2012.
- J.-S. Xu, K. Sun, C.-F. Li et al., “Experimental recovery of quantum correlations in absence of system-environment back-action,” Nature Communications, vol. 4, article 2851, 2013.
- R. L. Franco, B. Bellomo, S. Maniscalco, and G. Compagno, “Dynamics of quantum correlations in two-qubit systems within non-Markovian environments,” International Journal of Modern Physics B: Condensed Matter Physics; Statistical Physics; Atomic, Molecular and Optical Physics, vol. 27, no. 1–3, Article ID 1345053, 20 pages, 2013.
- J. Maziero, “Computing partial traces and reduced density matrices,” http://arxiv.org/abs/1601.07458.
- E. Brüning, H. Mäkelä, A. Messina, and F. Petruccione, “Parametrizations of density matrices,” Journal of Modern Optics, vol. 59, no. 1, pp. 1–20, 2012.
- R. A. Bertlmann and P. Krammer, “Bloch vectors for qudits,” Journal of Physics A: Mathematical and Theoretical, vol. 41, no. 23, Article ID 235303, 21 pages, 2008.
- J. Maziero, “Fortran code for generating random probability vectors, unitaries, and quantum states,” Frontiers in ICT, vol. 3, article 4, 2016.
- Y. Huang, “Computing quantum discord is NP-complete,” New Journal of Physics, vol. 16, Article ID 033027, 2014.
- T. Chanda, T. Das, D. Sadhukhan, A. K. Pal, A. Sen De, and U. Sen, “Reducing computational complexity of quantum correlations,” Physical Review A, vol. 92, no. 6, Article ID 062301, 2015.
- S. Luo, “Quantum discord for two-qubit systems,” Physical Review A, vol. 77, Article ID 042303, 2008.
- S. Luo and Q. Zhang, “Observable correlations in two-qubit states,” Journal of Statistical Physics, vol. 136, no. 1, pp. 165–177, 2009.
- M. Ali, A. R. P. Rau, and G. Alber, “Quantum discord for two-qubit X states,” Physical Review A, vol. 81, Article ID 042105, 2010.
- J. Maziero, L. C. Celeri, and R. M. Serra, “Symmetry aspects of quantum discord,” https://arxiv.org/abs/1004.2082.
- D. Girolami and G. Adesso, “Quantum discord for general two-qubit states: analytical progress,” Physical Review A, vol. 83, no. 5, Article ID 052108, 2011.
- Q. Chen, C. Zhang, S. Yu, X. X. Yi, and C. H. Oh, “Quantum discord of two-qubit X state,” Physical Review A, vol. 84, no. 4, Article ID 042313, 2011.
- A. S. M. Hassan, B. Lari, and P. S. Joag, “Tight lower bound to the geometric measure of quantum discord,” Physical Review A, vol. 85, no. 2, Article ID 024302, 2012.
- F. Ciccarello, T. Tufarelli, and V. Giovannetti, “Toward computability of trace distance discord,” New Journal of Physics, vol. 16, Article ID 013038, 2014.
- D. Spehner and M. Orszag, “Geometric quantum discord with Bures distance: the qubit case,” Journal of Physics A: Mathematical and Theoretical, vol. 47, no. 3, Article ID 035302, 20 pages, 2014.
- M. A. Jafarizadeh, N. Karimi, and H. Zahir, “Quantum discord for generalized bloch sphere states,” European Physical Journal D, vol. 68, no. 5, article 136, 2014.
- M. Daoud, R. Ahl Laamara, and S. Seddik, “Hilbert-Schmidt measure of pairwise quantum discord for three-qubit X states,” Reports on Mathematical Physics, vol. 76, no. 2, pp. 207–230, 2015.
- A. Beggi, F. Buscemi, and P. Bordone, “Analytical expression of genuine tripartite quantum discord for symmetrical X-states,” Quantum Information Processing, vol. 14, no. 2, pp. 573–592, 2015.
- Y.-K. Wang, N. Jing, S.-M. Fei, Z.-X. Wang, J.-P. Cao, and H. Fan, “One-way deficit of two-qubit X states,” Quantum Information Processing, vol. 14, no. 7, pp. 2487–2497, 2015.
- Z. Ma, Z. Chen, F. F. Fanchini, and S.-M. Fei, “Quantum discord for d 2 systems,” Scientific Reports, vol. 5, Article ID 10262, 2015.
- N. Jing, X. Zhang, and Y.-K. Wang, “Comment on: oneway deficit of two qubit X states,” Quantum Information Processing, vol. 14, pp. 4511–4521, 2015.
- G. Li, Y. Liu, H. Tang, X. Yin, and Z. Zhang, “Analytic expression of quantum correlations in qutrit Werner states undergoing local and nonlocal unitary operations,” Quantum Information Processing, vol. 14, no. 2, pp. 559–572, 2015.
- P. C. Obando, F. M. Paula, and M. S. Sarandy, “Trace-distance correlations for X states and the emergence of the pointer basis in Markovian and non-Markovian regimes,” Physical Review A, vol. 92, no. 3, Article ID 032307, 2015.
- Y. Huang, “Quantum discord for two-qubit X states: analytical formula with very small worst-case error,” Physical Review A, vol. 88, no. 1, Article ID 014302, 2013.
- A. Brodutch and K. Modi, “Criteria for measures of quantum correlations,” Quantum Information & Computation, vol. 12, no. 9-10, article 0721, 2012.
- B. Dakić, V. Vedral, and C. Brukner, “Necessary and sufficient condition for nonzero quantum discord,” Physical Review Letters, vol. 105, no. 19, Article ID 190502, 2010.
- S. Luo and S. Fu, “Geometric measure of quantum discord,” Physical Review A, vol. 82, no. 3, Article ID 034302, 4 pages, 2010.
- S. J. Akhtarshenas, H. Mohammadi, S. Karimi, and Z. Azmi, “Computable measure of quantum correlation,” Quantum Information Processing, vol. 14, no. 1, pp. 247–267, 2015.
Copyright © 2016 Jonas Maziero. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.