Bures Metric

Helstrom">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.

The Bures metric may be defined as

${\displaystyle [D_{B}(\rho ,\rho +d\rho )]^{2}={\frac {1}{2}}{\mbox{tr}}(d\rho G),}$ 

where ${\displaystyle G}$  is the Hermitian 1-form operator implicitly given by

${\displaystyle \rho G+G\rho =d\rho ,}$ 

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.

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

${\displaystyle D_{B}(\rho _{1},\rho _{2})^{2}=2(1-{\sqrt {F(\rho _{1},\rho _{2})}}),}$ 

where the fidelity function is defined as

${\displaystyle F(\rho _{1},\rho _{2})=\left[{\mbox{tr}}({\sqrt {{\sqrt {\rho _{1}}}\rho _{2}{\sqrt {\rho _{1}}}}})\right]^{2}.}$ 

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

${\displaystyle D_{A}(\rho _{1},\rho _{2})=\arccos {\sqrt {F(\rho _{1},\rho _{2})}},}$ 

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

When both density operators are diagonal (so that they are just classical probability distributions), then let ${\displaystyle \rho _{1}=diag(p_{1},...)}$  and similarly ${\displaystyle \rho _{2}=diag(q_{1},...)}$ , then the fidelity is

$${\displaystyle {\sqrt {F}}=\sum _{i}{\sqrt {p_{i}q_{i}}}}$$

 

${\displaystyle \arccos \left(\sum _{i}{\sqrt {p_{i}q_{i}}}\right)}$ 

${\displaystyle p,q}$ 

${\displaystyle ds^{2}={\frac {1}{2}}\sum _{i}{\frac {dp_{i}^{2}}{p_{i}}}}$ 

Perform a change of variables with ${\displaystyle x_{i}:={\sqrt {p_{i}}}}$ , then the chi-squared metric becomes ${\displaystyle ds^{2}=\sum _{i}dx_{i}^{2}}$ . Since ${\displaystyle \sum _{i}x_{i}^{2}=\sum _{i}p_{i}=1}$ , the points ${\displaystyle x}$  are restricted to move on the positive quadrant of a unit hypersphere. So, the geodesics are just the great circles on the hypersphere, and we also obtain the Wootters distance formula.

If both density operators are pure states, ${\displaystyle \psi ,\phi }$ , then the fidelity is ${\displaystyle {\sqrt {F}}=|\langle \psi |\phi \rangle |}$ , and we obtain the quantum version of Wootters distance

${\displaystyle \arccos(|\langle \psi |\phi \rangle |)}$ .

In particular, the direct Bures distance between any two orthogonal states is ${\displaystyle {\sqrt {2}}}$ , while the Bures distance summed along the geodesic path connecting them is ${\displaystyle \pi /2}$ .

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

${\displaystyle [D_{B}(\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 },}$ 

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

${\displaystyle {\frac {\rho L_{\mu }+L_{\mu }\rho }{2}}={\frac {d\rho ^{\,}}{d\theta ^{\mu }}}.}$ 

In this way, one has

${\displaystyle [D_{B}(\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 },}$ 

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

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

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

${\displaystyle J_{\mu \nu }^{}=4g_{\mu \nu }.}$ 

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.

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

${\displaystyle [D_{B}(\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],}$ 

${\displaystyle [D_{B}(\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].}$ 

For general systems, the Bures metric can be written in terms of the eigenvectors and eigenvalues of the density matrix ${\displaystyle \rho =\sum _{j=1}^{n}\lambda _{j}|j\rangle \langle j|}$  as

${\displaystyle [D_{B}(\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}}},}$ 

as an integral,

${\displaystyle [D_{B}(\rho ,\rho +d\rho )]^{2}={\frac {1}{2}}\int _{0}^{\infty }{\text{tr}}[e^{-\rho t}d\rho e^{-\rho t}d\rho ]\ dt,}$ 

or in terms of Kronecker product and vectorization,

${\displaystyle [D_{B}(\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 ],}$ 

where ${\displaystyle ^{*}}$  denotes complex conjugate, and ${\displaystyle ^{\dagger }}$  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.

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

${\displaystyle \rho ={\frac {1}{2}}(I+{\boldsymbol {r\cdot \sigma }}),}$ 

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

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

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

${\displaystyle dV_{B}={\frac {d^{3}{\boldsymbol {r}}}{8{\sqrt {1-r^{2}}}}},}$ 

which can be used to calculate the Bures volume as

${\displaystyle V_{B}=\iiint _{r^{2}\leq 1}{\frac {d^{3}{\boldsymbol {r}}}{8{\sqrt {1-r^{2}}}}}={\frac {\pi ^{2}}{8}}.}$ 

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

${\displaystyle \rho ={\frac {1}{3}}(I+{\sqrt {3}}\sum _{\nu =1}^{8}\xi _{\nu }\lambda _{\nu }),}$ 

where ${\displaystyle \lambda _{\nu }}$  are the eight Gell-Mann matrices and ${\displaystyle {\boldsymbol {\xi }}\in \mathbb {R} ^{8}}$  the 8-dimensional Bloch vector satisfying certain constraints.

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