Bures metric

In mathematics, in the area of quantum information geometry, the Bures metric (named after Donald Bures) or Helstrom metric (named after Carl W. Helstrom) defines an infinitesimal distance between density matrix operators defining quantum states. It is a quantum generalization of the Fisher information metric, and is identical to the Fubini–Study metric when restricted to the pure states alone.

Definition

The Bures metric [math]\displaystyle{ G }[/math] may be defined as

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2}\mbox{tr}( d \rho G ), where \lt math\gt G }[/math] is Hermitian 1-form operator implicitly given by

[math]\displaystyle{ \rho G + G \rho = d \rho, }[/math]

which is a special case of a continuous Lyapunov equation.

Some of the applications of the Bures metric include that given a target error, it allows the calculation of the minimum number of measurements to distinguish two different states and the use of the volume element as a candidate for the Jeffreys prior probability density for mixed quantum states.

Bures distance

The Bures distance is the finite version of the infinitesimal square distance described above and is given by

[math]\displaystyle{ D_B(\rho_1,\rho_2)^2 = 2(1-\sqrt{F(\rho_1,\rho_2)}), }[/math]

where the fidelity function is defined as

[math]\displaystyle{ F(\rho_1,\rho_2) = \left[ \mbox{tr}( \sqrt{ \sqrt{\rho_1}\rho_2\sqrt{\rho_1}})\right]^2. }[/math]

Another associated function is the Bures arc also known as Bures angle, Bures length or quantum angle, defined as

[math]\displaystyle{ D_A(\rho_1,\rho_2) = \arccos \sqrt{F(\rho_1,\rho_2)}, }[/math]

which is a measure of the statistical distance between quantum states.

Quantum Fisher information

The Bures metric can be seen as the quantum equivalent of the Fisher information metric and can be rewritten in terms of the variation of coordinate parameters as

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2} \mbox{tr}\left( \frac{d \rho}{d \theta^{\mu}} L_{\nu} \right) d \theta^{\mu} d\theta^{\nu}, }[/math]

which holds as long as [math]\displaystyle{ \rho }[/math] and [math]\displaystyle{ \rho+d\rho }[/math] have the same rank. In cases where they do not have the same rank, there is an additional term on the right hand side. [math]\displaystyle{ L_\mu }[/math] is the Symmetric logarithmic derivative operator (SLD) defined from

[math]\displaystyle{ \frac{\rho L_{\mu} + L_{\mu} \rho}{2} = \frac{d \rho^{\,}}{d \theta^{\mu}}. }[/math]

In this way, one has

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2} \mbox{tr}\left[ \rho \frac{L_{\mu} L_{\nu} + L_{\nu} L_{\mu}}{2} \right] d \theta^{\mu} d\theta^{\nu}, }[/math]

where the quantum Fisher metric (tensor components) is identified as

[math]\displaystyle{ J_{\mu \nu} = \mbox{tr}\left[ \rho \frac{L_{\mu} L_{\nu} + L_{\nu} L_{\mu}}{2}\right]. }[/math]

The definition of the SLD implies that the quantum Fisher metric is 4 times the Bures metric. In other words, given that [math]\displaystyle{ g_{\mu\nu} }[/math] are components of the Bures metric tensor, one has

[math]\displaystyle{ J_{\mu\nu}^{ } = 4 g_{\mu \nu}. }[/math]

As it happens with the classical Fisher information metric, the quantum Fisher metric can be used to find the Cramér–Rao bound of the covariance.

Explicit formulas

The actual computation of the Bures metric is not evident from the definition, so, some formulas were developed for that purpose. For 2x2 and 3x3 systems, respectively, the quadratic form of the Bures metric is calculated as

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{4}\mbox{tr}\left[ d \rho d \rho + \frac{1}{\det(\rho)}(\mathbf{1}-\rho)d\rho (\mathbf{1}-\rho)d\rho \right], }[/math]

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{4}\mbox{tr}\left[ d \rho d \rho + \frac{3}{1-\mbox{tr} \rho^3} (\mathbf{1}-\rho)d\rho (\mathbf{1}-\rho)d\rho + \frac{3 \det{\rho} }{1-\mbox{tr} \rho^3} (\mathbf{1}-\rho^{-1})d\rho (\mathbf{1}-\rho^{-1})d\rho \right]. }[/math]

For general systems, the Bures metric can be written in terms of the eigenvectors and eigenvalues of the density matrix [math]\displaystyle{ \rho=\sum_{j=1}^n\lambda_j|j\rangle\langle j| }[/math] as

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2} \sum_{j,k=1}^{n} \frac{|\langle j| d\rho | k\rangle |^2}{\lambda_j+\lambda_k}, }[/math]

as an integral,

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2}\int_0^\infty\text{tr}[e^{-\rho t}d\rho e^{-\rho t}d\rho]\ dt, }[/math]

or in terms of Kronecker product and vectorization,

[math]\displaystyle{ [d(\rho, \rho+d\rho)]^2 = \frac{1}{2}\text{vec}[d\rho]^\dagger\big(\rho^*\otimes \mathbf{1}+\mathbf{1}\otimes\rho\big)^{-1}\text{vec}[d\rho], }[/math]

where [math]\displaystyle{ ^* }[/math] denotes complex conjugate, and [math]\displaystyle{ ^\dagger }[/math] denotes conjugate transpose. This formula holds for invertible density matrices. For non-invertible density matrices, the inverse above is substituted by the Moore-Penrose pseudoinverse. Alternatively, the expression can be also computed by performing a limit on a certain mixed and thus invertible state.

Two-level system

The state of a two-level system can be parametrized with three variables as

[math]\displaystyle{ \rho = \frac{1}{2}( I + \boldsymbol{r\cdot\sigma} ), }[/math]

where [math]\displaystyle{ \boldsymbol{\sigma} }[/math] is the vector of Pauli matrices and [math]\displaystyle{ \boldsymbol{r} }[/math] is the (three-dimensional) Bloch vector satisfying [math]\displaystyle{ r^2\stackrel{\mathrm{def}}{=} \boldsymbol{r\cdot r} \le 1 }[/math]. The components of the Bures metric in this parametrization can be calculated as

[math]\displaystyle{ \mathsf{g} = \frac{\mathsf{I}}{4}+\frac{\boldsymbol{r\otimes r}}{4(1-r^2)} }[/math].

The Bures measure can be calculated by taking the square root of the determinant to find

[math]\displaystyle{ dV_B = \frac{d^3\boldsymbol{r}}{8\sqrt{ 1 - r^2}}, }[/math]

which can be used to calculate the Bures volume as

[math]\displaystyle{ V_B = \iiint_{r^2\leq 1}\frac{d^3\boldsymbol{r}}{8\sqrt{1-r^2}} = \frac{\pi^2}{8} . }[/math]

Three-level system

The state of a three-level system can be parametrized with eight variables as

[math]\displaystyle{ \rho = \frac{1}{3}( I + \sqrt{3} \sum_{\nu=1}^8\xi_\nu\lambda_\nu), }[/math]

where [math]\displaystyle{ \lambda_\nu }[/math] are the eight Gell-Mann matrices and [math]\displaystyle{ \boldsymbol \xi \in\mathbb{R}^8 }[/math] the 8-dimensional Bloch vector satisfying certain constraints.

See also

• Fubini–Study metric
• Fidelity of quantum states
• Fisher information
• Fisher information metric

Further reading

• Uhlmann, A. (1992). "The Metric of Bures and the Geometric Phase". in Gielerak, R.; Lukierski, J.; Popowicz, Z.. Groups and Related Topics. Proceedings of the First Max Born Symposium. pp. 267–274. doi:10.1007/978-94-011-2801-8_23. ISBN 94-010-5244-1. 
• Sommers, H. J.; Zyczkowski, K. (2003). "Bures volume of the set of mixed quantum states". Journal of Physics A 36 (39): 10083–10100. doi:10.1088/0305-4470/36/39/308. Bibcode: 2003JPhA...3610083S. 
• Dittmann, J. (1993). "On the Riemannian Geometry of Finite Dimensional Mixed States". Seminar Sophus Lie 73. http://www.heldermann-verlag.de/jlt/jlt03/DITTLAT.PDF. 
• Slater, Paul B. (1996). "Quantum Fisher-Bures information of two-level systems and a three-level extension". J. Phys. A: Math. Gen. 29 (10): L271–L275. doi:10.1088/0305-4470/29/10/008. 
• Nielsen, M. A.; Chuang, I. L. (2000). Quantum Computation and Quantum Information. Cambridge University Press. ISBN 0-521-63235-8.