贡献者： certain_pineapple
12
在很多时候我们需要去探讨如何计算两个量子态究竟有多么相近，在态矢量的语言下这个答案是很简单的，我们可以使用两个态矢量的内积的模方来描述其二者是否相同或有多么不同。但对于两个密度矩阵来说这个问题并没有那么简单，迹距离和保真度给出了两种度量方法，我们在这里首先介绍迹距离。
定义 1 迹距离
两密度矩阵 $\rho$ 和 $\sigma$ 之间的迹距离 $T\left(\rho,~\sigma\right)$ 定义如下：
\begin{equation}
T\left(\rho,~\sigma\right) \equiv \frac{1}{2} \operatorname {tr} \left\lvert \rho - \sigma \right\rvert ~.
\end{equation}
其中，$ \left\lvert A \right\rvert = \sqrt{A^\dagger A}$
1. 概率分布的全变差距离
为了引入密度矩阵之间的迹距离的概念，我们首先从概率分布的全变差距离开始介绍。
定义 2 概率分布的全变差距离
对于两个样本空间是可数集的概率分布 $\left\{ p_x \right\}$ 和 $\left\{q_x\right\}$，定义其全变差距离为：
\begin{equation}
T\left( p_x,~q_x \right) = \frac{1}{2}\sum_x \left\lvert p_x - q_x \right\rvert ~.
\end{equation}
全变差距离是一种度量
对于定义 1 ，从定义式中可以看到 $T\left(p_x,~q_x\right) = T\left(q_x,~p_x\right)\geqslant 0$，即其是对称且非负的。且对三个概率分布 $\left\{p_x^{(1)}\right\}$，$\left\{p_x^{(2)}\right\}$ 和 $\left\{p_x^{(3)}\right\}$ 有：
\begin{equation}
\begin{aligned}
T\left(p_x^{(1)},~p_x^{(2)}\right) &= \frac{1}{2}\sum_x \operatorname {tr} \left\lvert p_x^{(1)} - p_x^{(2)} \right\rvert \\
&= \frac{1}{2}\sum_x \operatorname {tr} \left\lvert \left(p_x^{(1)} - p_x^{(3)}\right) - \left(p_x^{(2)} - p_x^{(3)}\right) \right\rvert \\
&\leqslant\frac{1}{2}\sum_x \operatorname {tr} \left\lvert p_x^{(1)} - p_x^{(3)} \right\rvert + \frac{1}{2}\sum_x \operatorname {tr} \left\lvert p_x^{(2)} - p_x^{(3)} \right\rvert \\
&=T\left(p_x^{(1)},~p_x^{(3)}\right)+T\left(p_x^{(2)},~p_x^{(3)}\right)~,
\end{aligned}~
\end{equation}
由此可见其的确构成一个度量。
全变差距离的另一种等价形式
接下来，我们将证明对于指标集 $U = \left\{x\right\}$ 和概率分布 $\left\{p_x\right\}$ 和 $\left\{q_x\right\}$，有
\begin{equation}
T\left(p_x,~q_x\right) = \max_{S \subseteq U}\left(\sum_{x\in S}p_x - \sum_{x\in S}q_x\right)~.
\end{equation}
定义 $r_x = p_x - q_x$，由于概率的归一性，有
\begin{equation}
\sum\limits_{x\in U}r_x = 0~.
\end{equation}
记 $S_+ = \left\{x \in U | r_x \geqslant 0\right\}$，$S_- = \left\{x \in U | r_x < 0\right\}$。则有：
\begin{equation}
\sum_{x\in S_+}r_x + \sum_{x\in S_-}r_x = 0~.
\end{equation}
进而：
\begin{equation}
\sum_{x\in S_+}r_x = -\sum_{x\in S_-}r_x~.
\end{equation}
综上所述，有：
\begin{equation}
\begin{aligned}
T\left(p_x,~q_x\right) &= \frac{1}{2}\sum_{x\in U} \left\lvert p_x - q_x \right\rvert \\
&= \frac{1}{2}\sum_{x\in U} \left\lvert r_x \right\rvert \\
&= \frac{1}{2}\sum_{x\in S_+}r_x - \frac{1}{2}\sum_{x\in S_-}r_x \\
&= \sum_{S\in S_+}r_x \\
&= \max_{S \subseteq U}\sum_{x\in S}r_x \\
&= \max_{S\subseteq U}\left(\sum_{x\in S}p_x - \sum_{x\in S}q_x\right)
\end{aligned}~
\end{equation}
2. 密度矩阵的迹距离
从式 2 和式 4 都可以过度得到密度矩阵的迹距离的表达形式：
\begin{equation}
T\left(\rho,~\sigma\right) = \frac{1}{2} \operatorname {tr} \left\lvert \rho - \sigma \right\rvert ~.
\end{equation}
\begin{equation}
T\left(\rho,~\sigma\right) = \max_{ P } \operatorname {tr}\left[P\left(\rho - \sigma\right)\right]~.
\end{equation}
其中 $P$ 代表任意本征值小于 $1$ 的半正定算子，我们首先证明式 10 与式 9 的等价性。
$\left(\rho - \sigma\right)$ 是一个厄米矩阵，那么有：
\begin{equation}
\begin{aligned}
\rho - \sigma &= \sum_i \lambda_i \left| i \right\rangle \left\langle i \right| \\
&= \sum_{\lambda_i \geqslant 0}\lambda_i \left| i \right\rangle \left\langle i \right| + \sum_{\lambda_i<0}\lambda_i \left| i \right\rangle \left\langle i \right| \\
&= M_1 - M_2~.
\end{aligned}~
\end{equation}
其中 $M_1$，$M_2$ 都是半正定算符。
同时：
\begin{equation}
\begin{aligned}
\left\lvert \rho - \sigma \right\rvert &= \sqrt{\left( \rho - \sigma \right)^\dagger \left( \rho - \sigma \right)} \\
&= \sqrt{\sum_i \lambda_i^2 \left| i \right\rangle \left\langle i \right| } \\
&= \sum_i \left\lvert \lambda_i \right\rvert \left| i \right\rangle \left\langle i \right| \\
&= M_1 + M_2~,
\end{aligned}~
\end{equation}
则有 $T\left(\rho,~\sigma\right) = \frac{1}{2}\left( \operatorname {tr}M_1 + \operatorname {tr}M_2\right)$。
考虑 $ \operatorname {tr}\rho = \operatorname {tr}\sigma = 1$，则 $ \operatorname {tr}M_1 - \operatorname {tr}M_2 = \operatorname {tr}\left(\rho - \sigma\right) = 0$，即 $ \operatorname {tr}M_1 = \operatorname {tr}M_2$，$T\left(\rho,~\sigma\right) = \operatorname {tr}M_1$。
设 $P$ 为任意本征值小于 $1$ 的半正定算子，则：
\begin{equation}
\begin{aligned}
\operatorname {tr} \left[ P\left(\rho - \sigma\right) \right] &= \operatorname {tr}\left[ P\left(M_1 - M_2\right) \right] \\
&= \operatorname {tr}\left(PM_1\right) - \operatorname {tr}\left(PM_2\right) \\
&\leqslant \operatorname {tr}\left(PM_1\right) \\
&\leqslant \operatorname {tr}M_1 \\
&= T\left( \rho,~\sigma \right)~.
\end{aligned}~
\end{equation}
所以有 $T\left(\rho,~\sigma\right) = \max_{ P } \operatorname {tr}\left[P\left(\rho - \sigma\right)\right]$。
单量子比特的密度矩阵的迹距离
考虑单 qubit 的密度矩阵的好处是此时密度矩阵可以写成式 8 $\rho = \frac{1}{2}I + \frac{1}{2} \boldsymbol{\mathbf{n}} \cdot \boldsymbol{\mathbf{\sigma}} $。
设 $\rho = \frac{1}{2}I + \frac{1}{2} \boldsymbol{\mathbf{n}} _1\cdot \boldsymbol{\mathbf{\sigma}} $，$\delta = \frac{1}{2}I + \frac{1}{2} \boldsymbol{\mathbf{n}} _2\cdot \boldsymbol{\mathbf{\sigma}} $，那么有：
\begin{equation}
\begin{aligned}
T\left(\rho,~\delta\right) &= \frac{1}{2} \operatorname {tr} \left\lvert \rho - \delta \right\rvert \\
&= \frac{1}{4} \operatorname {tr} \left\lvert \left( \boldsymbol{\mathbf{n}} _1 - \boldsymbol{\mathbf{n}} _2\right) \cdot \boldsymbol{\mathbf{\sigma}} \right\rvert ~,
\end{aligned}~
\end{equation}
考虑 $\left( \boldsymbol{\mathbf{n}} _1 - \boldsymbol{\mathbf{n}} _2\right) \cdot \boldsymbol{\mathbf{\sigma}} $ 的本征值是 $\pm \left\lvert \boldsymbol{\mathbf{n}} _1 - \boldsymbol{\mathbf{n}} _2 \right\rvert $，则：
\begin{equation}
\begin{aligned}
T\left(\rho,~\delta\right)
&= \frac{1}{4} \operatorname {tr} \left\lvert \left( \boldsymbol{\mathbf{n}} _1 - \boldsymbol{\mathbf{n}} _2\right) \cdot \boldsymbol{\mathbf{\sigma}} \right\rvert \\
&= \frac{1}{2} \left\lvert \boldsymbol{\mathbf{n}} _1 - \boldsymbol{\mathbf{n}} _2 \right\rvert ~.
\end{aligned}~
\end{equation}
这也就是说，对于单量子比特的密度矩阵，迹距离就是其 Bloch 矢量的差的模长的一半。
迹距离的性质
- 迹距离是密度矩阵空间的度量，满足非负，对称和三角不等式，且当且仅当 $\rho$ 和 $\sigma$ 相等的时候 $T\left(\rho,~\sigma\right) = 0$。
- $0 \leqslant T\left(\rho,~\sigma\right) \leqslant 1$。
- $T\left(\rho,~\sigma\right)$ 在酉元作用下保持不变，$T\left(U^\dagger\rho U,~U^\dagger\sigma U\right) = T\left(\rho,~\sigma\right)$。
- $T\left(\varepsilon(\rho),~\varepsilon(\sigma)\right) \leqslant T\left(\rho,~\sigma\right)$，其中 $\varepsilon$ 是一个保迹算子，$ \operatorname {tr}\varepsilon(\rho) = \operatorname {tr}\rho$。
- 对输入保持下凸，即 $T\left(\sum_i p_i\rho_i,~\sigma\right) \leqslant \sum_i p_iT\left(\rho_i,~\sigma\right)$。
- 对于纯态满足 $T( \left| \psi \right\rangle \left\langle \psi \right| ,~ \left| \phi \right\rangle \left\langle \phi \right| ) = \sqrt{1 - \left\lvert \left\langle \psi \middle| \phi \right\rangle \right\rvert ^2}$。
1. ^ 本节参考 [1] 和Wikipedia 界面
2. ^ 本节中部分证明参考项目
[1] ^ Michael A. Nielsen，Isaac L. Chuang 著，郑大钟 赵千川 译 量子计算和量子信息（二）——量子信息部分 清华大学出版社