贡献者： certain_pineapple; addis
1在很多时候我们需要去探讨如何计算两个量子态究竟有多么相近，在态矢量的语言下这个答案是很简单的，我们可以使用两个态矢量的内积的模方来描述其二者是否相同或有多么不同。但对于两个密度矩阵来说这个问题并没有那么简单，迹距离和保真度给出了两种度量方法，我们在这里介绍保真度。
1. 定义
定义 1 保真度
对于状态 $\rho$ 和 $\sigma$，保真度的定义为：
\begin{equation}
F\left(\rho, \sigma\right)\equiv \operatorname {tr}\sqrt{\rho^{\frac{1}{2}}\sigma\rho^{\frac{1}{2}}}~.
\end{equation}
在一些文献中的保真度会采用另外一种定义:
定义 2 保真度
对于状态 $\rho$ 和 $\sigma$，保真度的定义为：
\begin{equation}
F_1\left(\rho, \sigma\right)\equiv\left( \operatorname {tr}\sqrt{\rho^{\frac{1}{2}}\sigma\rho^{\frac{1}{2}}}\right)^2~.
\end{equation}
由于定义 1 有更好的经典对应，因此我们此处选取第一种定义，在使用和阅读保真度相关内容时，应注意究竟采用的是哪种定义。
值得注意的是，保真度并不是数学意义上的度量，因为当 $\sigma = \rho$ 时候，$F\left(\rho, \sigma\right) = \operatorname {tr}\sqrt{\rho^{\frac{1}{2}}\rho\rho^{\frac{1}{2}}} = \operatorname {tr}\rho = 1 \neq 0$。
我们从经典信息的角度定义来阐释保真度的含义，考虑经典概率的归一性限制，对于两个概率分布 $\left\{p_i\right\}$ 和 $\left\{q_i\right\}$，有 $\sum\limits_i p_i = \sum\limits_i q_i = 1$。由于概率的性质 $p_i,q_i\geqslant 0$，考虑向量 $ \boldsymbol{\mathbf{v}} _p = \left( \sqrt{p_1},\sqrt{p_2},\cdots\sqrt{p_N}\right)$ 和 $ \boldsymbol{\mathbf{v}} _q = \left( \sqrt{q_1},\sqrt{q_2},\cdots\sqrt{q_N}\right)$，那么他们就是模长为 $1$ 的单位向量，可以视为 $N$ 维单位球面上的单位向量。而他们的内积为 $\left( \boldsymbol{\mathbf{v}} _p, \boldsymbol{\mathbf{v}} _q\right) = \sum\limits_i\sqrt{p_iq_i}$。
接下来我们从保真度的定义出发，如果 $\rho$ 和 $\sigma$ 可以在同一组基 $\left\{ \left\lvert i \right\rangle \right\}$ 下写成对角矩阵的形式，$\rho = \sum\limits_ip_i \left| i \right\rangle \left\langle i \right| $，$\sigma = \sum\limits_iq_i \left| i \right\rangle \left\langle i \right| $，那么可以计算：
\begin{equation}
\begin{aligned}
F\left(\rho,\sigma\right) &= \operatorname {tr}\sqrt{\sum_{i} \sqrt{p_i} \left| i \right\rangle \left\langle i \right| \sum_{j}q_j \left| j \right\rangle \left\langle j \right| \sum_k\sqrt{p_k} \left| k \right\rangle \left\langle k \right| } \\
&= \operatorname {tr}\sqrt{\sum_{i,j,k}\sqrt{p_ip_k}q_j \left\lvert i \right\rangle \left\langle i \middle| j \right\rangle \left\langle j \middle| k \right\rangle \left\langle k \right\rvert } \\
&= \operatorname {tr}\sqrt{\sum_ip_iq_i \left| i \right\rangle \left\langle i \right| } \\
&= \sum_i\sqrt{p_iq_i}~.
\end{aligned}~
\end{equation}
由此可见保真度概念实际上是我们在经典信息中考虑的 $ \boldsymbol{\mathbf{v}} _p$ 和 $ \boldsymbol{\mathbf{v}} _q$ 的量子力学推广。
保真度的另一种表达形式
保真度可以由另一种表达形式给出:
\begin{equation}
F\left(\rho, \sigma\right) = \operatorname {tr} \left\lvert \sqrt{\rho}\sqrt{\sigma} \right\rvert ~,
\end{equation}
其中，$ \left\lvert A \right\rvert = \sqrt{AA^\dagger}$，证明十分简单，直接将绝对值的定义式代入
式 4 中即可得到原本的定义式。
讨论：为什么定义不是 $F\left(\rho,\sigma\right) = \operatorname {tr}\sqrt{\rho\sigma}~$？
2
我们可以很容易的验证如果 $\rho = \sum\limits_i p_i \left| i \right\rangle \left\langle i \right| $，$\sigma = \sum\limits_i q_i \left| i \right\rangle \left\langle i \right| $，有:
\begin{equation}
\begin{aligned}
\operatorname {tr}\sqrt{\rho\sigma} &= \operatorname {tr}\sqrt{\sum_{i,j}p_iq_j \left\lvert i \right\rangle \left\langle i \middle| j \right\rangle \left\langle j \right\rvert } \\
&= \operatorname {tr}\sqrt{\sum_ip_iq_i \left| i \right\rangle \left\langle i \right| } \\&=\sum_i\sqrt{p_iq_i}~.
\end{aligned}~
\end{equation}
貌似通过这种定义也可以实现 $ \boldsymbol{\mathbf{v}} _p$ 和 $ \boldsymbol{\mathbf{v}} _q$ 的量子力学推广，那我们为什么一定要选择看起来较为麻烦的定义 1 呢？
回忆矩阵的函数的定义，我们是通过幂级数的方法定义的，而 $\sqrt{x}$ 在 $0$ 处没有泰勒展开，那我们为什么可以写 $\sqrt{\rho}$ 呢？这是因为 $\rho$ 是密度矩阵，所以是一个厄密矩阵，是必定可对角化且本征值为正数的，我们可以通过 $\rho = \sum\limits_ip_i \left| i \right\rangle \left\langle i \right| ,\sqrt{\rho} = \sum\limits_i\sqrt{p_i} \left| i \right\rangle \left\langle i \right| $ 来定义（注意，这个定义针对厄密矩阵，因为任意矩阵本征值不一定为正实数，这会引入根号的多值问题，而且甚至任意矩阵不一定可以相似对角化）。考虑 $\rho\sigma$，$\left(\rho\sigma\right)^\dagger = \sigma^\dagger \rho^\dagger = \sigma\rho$，除非 $\rho$ 和 $\sigma$ 对易，否则其乘积就不是厄密矩阵，因此我们直接对其乘积开根号是一个很危险的行为。因此我们选择定义 1 的形式，其中所有在根号下的内容的厄密性都很容易验证。
注：事实上也可以拓展开根号的定义来使得 $\sqrt{\rho\sigma}$ 合法化，这样得到的结果是和定义 1 等价的。
2. 性质
保真度在幺正变换下是不变的
保真度在幺正变换下是不变的，也就是对于任意给定幺正矩阵 $U$，有 $F\left(\rho,\sigma\right) = F\left(U \rho U^\dagger, U\sigma U^\dagger\right)$。
证明是较为直接的，考虑密度矩阵都是半正定矩阵，考虑半正定矩阵的性质：$\sqrt{U \rho U^\dagger} = U \sqrt{\rho} U^\dagger$。
那么可以写出：
\begin{equation}
\begin{aligned}
F\left(U \rho U^\dagger, U \sigma U^\dagger\right) &= \operatorname {tr}\sqrt{\sqrt{U \rho U^\dagger}U \sigma U^\dagger \sqrt{U \rho U^\dagger}} \\
&= \operatorname {tr}\sqrt{U\sqrt{\rho}U^\dagger U \sigma U^\dagger U \sqrt{\rho}U^\dagger} \\
&= \operatorname {tr}\left(U\sqrt{\sqrt{\rho}\sigma\sqrt{\rho}}U^\dagger\right) \\
&= \operatorname {tr}\left(U^\dagger U \sqrt{\sqrt{\rho}\sigma\sqrt{\rho}}\right) \\
&= F\left(\rho, \sigma\right)~.
\end{aligned}~
\end{equation}
保真度的对称性
保真度对其输入具有对称性，即 $F\left(\rho,\sigma\right) = F\left(\sigma,\rho\right)$。
一种证明方法是通过 Uhlmann 定理的推论来进行证明的，但我们在此给出一个更加简洁的证明方法，在后续介绍 Uhlmann 定理时，我们将重新提及这个问题。
考虑保真度 $F\left(\rho,\sigma\right)$ 的表达式式 4 ，考虑 $\sqrt{\rho}$ 和 $\sqrt{\sigma}$ 均为半正定矩阵，而两半正定矩阵的乘积可相似对角化，而通过矩阵的绝对值的运算性质可知对于一个可以相似对角化的矩阵 $M$，有 $ \operatorname {tr} \left\lvert M \right\rvert = \operatorname {tr} \left\lvert M^\dagger \right\rvert $。那么有：
\begin{equation}
\begin{aligned}
F\left(\rho,\sigma\right) &= \operatorname {tr} \left\lvert \sqrt{\rho}\sqrt{\sigma} \right\rvert \\
&= \operatorname {tr} \left\lvert \left(\sqrt{\rho}\sqrt{\sigma}\right)^\dagger \right\rvert \\
&= \operatorname {tr} \left\lvert \left(\sqrt{\sigma}\right)^\dagger\left(\sqrt{\rho}\right)^\dagger \right\rvert \\
&= \operatorname {tr} \left\lvert \sqrt{\sigma}\sqrt{\rho} \right\rvert \\
&= F\left(\sigma,\rho\right)~.
\end{aligned}~
\end{equation}
$\rho$ 和 $\sigma$ 的半正定性可以保证最后一步等式的成立。
1. ^ 本文参考 [1] 和 Wikipedia 相关页面。
2. ^ 参考网站
[1] ^ Michael A. Nielsen，Isaac L. Chuang 著，郑大钟 赵千川 译 量子计算和量子信息（二）——量子信息部分 清华大学出版社
致读者： 小时百科一直以来坚持所有内容免费，这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此，我们请求广大读者
热心打赏 ，使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元，我们一个星期内就能脱离亏损， 并保证在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款， 他们的付出帮助了 99% 的读者免费获取知识， 我们在此表示感谢。