量子计算基础整理（上）

量子计算基础整理

写在前面

年间本来是计划玩电动的，突然觉得有一点点索然无味。遂整理一下先前看过的《量子计算与编程入门》一书的一些知识点，算是读书笔记。

本文的知识点较为简单，主要分为以下两个部分：

1. 量子力学基础理论
2. 量子程序

量子力学部分主要为搭建一个电路所需要的最基础的知识点，量子程序也同样如此。

本人并没有量子力学的基础，但量子计算机关注的重点是计算机，而非量子。若对线性代数或矩阵论的内容较为熟悉，应该没有特别大的困难。

量子力学基础

量子的四个特性

”量子态就是一个微观粒子的状态。“

而由薛定谔方程，可以得到量子的第一个特性——一定会遇到不可分割的最小单位，这种最小单位统称为量子。
$$-\frac{{\hslash }^2}{2\mu }{\nabla }^2\Psi +U\Psi =E\Psi$$

“量子化的属性有很多种，这里优先考虑一种——能量。”

能量的量子化的表现形式为能级。

当我们谈到能级，就可以引出量子的第二个特性——跃迁。

《量子计算与编程入门》一书中对于这两个特性举了一个例子：把一层楼比作量子，对于量子而言，能级就如同拆掉楼梯和电梯的楼层一样，人们只能位于楼层中的某一层，而不能处于两个楼层的中间某个状态（能级）。

但这样并不代表我们没办法上下楼了，相反，如果我们需要上楼的话，我们可以直接超级跳跳上去（跃迁）：）

量子的第三个特性——量子叠加性的例子就非常有名了，那就是我们非常熟悉的薛定谔的猫。即当我们不去观测的时候，猫会处于生与死的叠加态。在我们先前的例子中，如果我们不去观测，每栋楼里的每个人同时存在于所有楼层。

PS：虽然量子有叠加性，但随观测次数的增加，会趋于一个稳定的结果（大数定律）。而宏观物体内的粒子数量非常庞大，所以宏观物体——大量粒子的集合体是绝无可能存在叠加性的。

“薛定谔宣称，不打开盒子，猫就处于生和死的’叠加态‘，又称‘当我们打开盒子，经过了我们的观察，猫就会坍缩到一个确定的生、死状态上’。”

得到这个态的概率是叠加态和测量态的内积的平方，且测量后叠加态就会坍缩到这个确定的态上。

这里便是量子的第四个特性——“测量与坍缩假设”。对于一个叠加态，测量的结果一定是量子化后确定、分立的态中的一个。

量子态的描述

下面就是各种矩阵运算了。

定义

量子态可以用线代中的向量描述，量子态的向量称为态矢，可以分为左矢（ket）与右矢（bra）。
右矢：$$|\psi ⟩=[c_1,c_2,...,c_n{]}^T$$
左矢：$$⟨\psi |=[c_1^{\ast },c_2^{\ast },...,c_n^{\ast }]$$
这样的组合就可以描述一个量子态，其中每一个分量都是复数，T表示转置，如果左右矢在括号内的描述相同，那么称为这两个矢量互为转置共轭。

有了量子态的矩阵表示，我们就可以定义内外积啦。

假设两个量子态矩阵（坐标）表示如下
$$|\alpha ⟩=[a_1,a_2,...,a_n{]}^T$$
$$|\beta ⟩=[b_1,b_2,...,b_n{]}^T$$
那么内积就可以定义为
$$⟨\alpha |\beta ⟩=\sum _{i=1}^n{a}_i\ast b_i$$
外积就可以定义为
$$|\alpha ⟩⟨\beta |=[a_i{b}_i^{\ast }{]}_{n\times n}$$
表示这是一个 nxn 的矩阵

同时，对于二能级量子，能量较高的称为激发态（excited state）记为
$$|e⟩=|1⟩=\left[\begin{array}{c}1\\ 0\end{array}\right]or⟨e|=⟨1|=\left[\begin{array}{cc}1& 0\end{array}\right]$$
能量较低的称为基态（ground state）记为

$$|g⟩=|0⟩=\left[\begin{array}{c}0\\ 1\end{array}\right]or⟨g|=⟨0|=\left[\begin{array}{cc}0& 1\end{array}\right]$$

为了与经典的比特进行类比，通常我们称激发态与基态为量子比特（qubit）。

此时，任意叠加态可以写作激发态与基态量子比特的线性组合
$$|\psi ⟩=\alpha |0⟩+\beta |1⟩$$
并满足归一化条件
$$|\alpha {|}^2+|\beta {|}^2=1$$

状态演化

假设封闭的量子系统演化由与时刻$t_1$与$t_2$酉变换$U$描述。那么在$t_1$与$t_2$时刻的量子态$|\psi ⟩$满足

$$|{\psi }_2⟩=U|{\psi }_1⟩$$

且满足
$$U{U}^{†}=I$$
$U^{†}$表示对$U$取转置共轭，由此可知$U$是可逆矩阵，酉变换也是一个可逆变换。

在量子计算中，酉矩阵也被称为量子门（类比逻辑门）

这里给大家介绍一下Pauli矩阵，Pauli矩阵就是一组量子门

$${\sigma }_0\equiv I\equiv \left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right],{\sigma }_1\equiv {\sigma }_x\equiv X\equiv \left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right]$$
$${\sigma }_2\equiv {\sigma }_y\equiv Y\equiv \left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right],{\sigma }_3\equiv {\sigma }_z\equiv Z\equiv \left[\begin{array}{cc}1& 0\\ 0& -1\end{array}\right]$$

其中X门对于任意量子态
$$X|\psi ⟩=\left[\begin{array}{cc}0& 1\\ 1& 0\end{array}\right]\left[\begin{array}{c}\alpha \\ \beta \end{array}\right]=\left[\begin{array}{c}\beta \\ \alpha \end{array}\right]=\beta |0⟩+\alpha |1⟩$$
所以X门也叫量子非门（quantum NOT gate）

叠加态与测量

上述的两个矢量可以构成一个二维空间的基。任何一个态都可以写为这两个基在复数空间上的线性组合
$$|\psi ⟩=a|0⟩+b{e}^{i\theta }|1⟩$$
$a$.$b$为满足归一化的实数；$e^{i\theta }$表示模为1、幅角为$\theta$的复数（欧拉公式）。

测量的定义为将量子态$|\psi ⟩$投影到另一个态$|\gamma ⟩$上。获得这个态的$|\gamma ⟩$概率是它们内积的平方，即
$$P_{\gamma }=|⟨\psi |\gamma ⟩{|}^2$$
其他概率下会将量子态投影到它的正交态上（由于没能投影到态$|\gamma ⟩$，故只能投影到与$|\gamma ⟩$不构成线性相关的态上）。
$$P_{\gamma \perp }=1-P_{\gamma }$$

相态，纯态和混合态

我们来看看下面的两个量子态：

$$|{\psi }_1⟩=\frac{1}{\sqrt{2}}(|0⟩+|1⟩)$$

$$|{\psi }_2⟩=\frac{1}{\sqrt{2}}(|0⟩-|1⟩)$$
然后算一算在$|0⟩$或者$|1⟩$方向上测量的概率

选取$|1⟩$作为方向，对于$|{\psi }_1⟩$
$$P_1=|⟨\psi |1⟩{|}^2=0.5$$
对于$|{\psi }_2⟩$，同样的
$$P_1=|⟨\psi |1⟩{|}^2=0.5$$

根本分不出来嗷，这是因为缺了个相位信息$\theta$

还有一种情况，假设有两个袋子，一个装满了$|{\psi }_1⟩=\frac{1}{\sqrt{2}}(|0⟩+|1⟩)$放在左侧，另一个是空的。

不断从左侧拿出一个态，测量（量子态塌缩），然后扔到右边的袋子，经过多次实验后，右侧袋子应该具有大致等量的$|0⟩$和$|1⟩$.

此时从右侧取出一个态，在不知道手上是什么态的情况下，能说手上的态是$\frac{1}{\sqrt{2}}(|0⟩+|1⟩)$吗？

答案是不行的，此时右侧的袋子已经变成了

{ ∣ ψ 0 ⟩ = ∣ 0 ⟩ : P 0 = 0.5 , ∣ ψ 1 ⟩ = ∣ 1 ⟩ : P 1 = 0.5 } \{|\psi_0\rangle=|0\rangle: P_0=0.5,|\psi_1\rangle=|1\rangle:P_1=0.5 \} {∣ψ0​⟩=∣0⟩:P0​=0.5,∣ψ1​⟩=∣1⟩:P1​=0.5}

这样的集合形式，同样没有了相位信息。

由此，引出定义，纯态不仅具有概率还有相位，混合态是纯态的概率性叠加，往往失去了全部或部分的相位信息。

混合态的表示

混合态可以描述成态集合和概率的列表，但这样太麻烦了，所以可以改为使用密度矩阵的形式进行描述。

对于纯态，密度矩阵的形式是
$$\rho =|\psi ⟩⟨\psi |$$
对于混合态，密度矩阵的形式是
$$\rho =\sum _i{P}_i|\psi ⟩⟨\psi |$$
其中， { P i , ∣ ψ i ⟩ } \{ P_i,|\psi_i\rangle \} {Pi​,∣ψi​⟩}是系统所处的态及概率

密度矩阵具有以下性质：

1. “对于一个两能级系统表述的态，不论是纯的还是混合的，都可以用密度矩阵$\rho$描述”
2. “$\rho ={\rho }^2$当且仅当量子态为纯态时成立”
3. “$\rho$对角线上的分量表示整个系统如果经历一次测量，可以得到的是这个态的概率。如果只操作和测量一个两能级系统，分辨不出相同的密度矩阵”

密度矩阵可以完备的描述一个两能级系统可能出现的任何状态。为了方便理解，这里引入布洛赫球的概念，方便百事一个量子比特的任意状态。

如果量子态是纯态，那么它是球上的点，z坐标衡量$|0⟩$和$|1⟩$，公式表示
$$P_0(|\psi ⟩)=\frac{1+z}{2}$$
$$P_1(|\psi ⟩)=\frac{1-z}{2}$$
沿平行于平面$xy$的方向，穿过圆面上一点的$z$坐标得到的圆就代表了相位的复平面；这个点于$x$的夹角就是单位复数的幅角。

由此，每个纯态都与球面上的点一一对应了起来。

由于混合态是多个纯态的经典统计概率叠加，每一个纯态分量（连接球面到球心的点）的加权平均即可得到混合态矢量。由此可知，混合态的点在布洛赫球的球内。

最大混合态是球心，不存在任何量子叠加性。

下面我们就拿这个点计算一下，由于布洛赫球上的（1，0，0）和（-1，0，0）是$x$方向的顶点与$-x$方向上的顶点。量子态概率分布就是$z$坐标，即$0$。
$$P_0(|{\psi }_1⟩)=P_0(|{\psi }_2⟩)=0.5$$
幅角${\theta }_1=0,{\theta }_2=\pi$，由此可推断
$$|{\psi }_1⟩=\frac{1}{\sqrt{2}}(|0⟩+|1⟩)$$

$$|{\psi }_2⟩=\frac{1}{\sqrt{2}}(|0⟩-|1⟩)$$

矩阵密度为
$$\rho =\frac{1}{2}|{\psi }_1⟩⟨{\psi }_1|+\frac{1}{2}|{\psi }_2⟩⟨{\psi }_2|=\left[\begin{array}{cc}0.5& 0\\ 0& 0.5\end{array}\right]$$
即最大混合态。

可观测量与量子观测

量子比特的信息无法直接获取，但是我们可以通过测量得到可观测的信息。

这时又要引出量子理论中自拌算子（self-adjoint operators）的概念，自拌有时也称为厄米（Hermitian）。

量子理论的可观测量与经典力学中的动力学量（位置、动量、角动量）对应，其他量如质量、电荷等作为参数引入到系统的哈密顿量（Hamiltonian）。

下面我们就来测量一下量子吧：

假设量子测量是由测量算子的集合 { M i } \{M_i\} {Mi​}来描述，这些算子可以作用在待测量系统的状态空间上，索引$i$表示实验上可能发生的结果。如果测量前量子系统处于最新状态$|\psi ⟩$，则
$$p(i)=⟨\psi |M_i^{†}{M}_i|\psi ⟩$$
测量后系统状态变为
$$\frac{M_i|\psi ⟩}{\sqrt{⟨\psi |M_i^{†}{M}_i|\psi ⟩}}$$

由于

$$\sum _{i}p(i)=\sum _i⟨\psi |M_i^{†}{M}_i|\psi ⟩=1$$

易推得
$$\sum _i{M}_i^{†}{M}_i=I$$
这个方程也叫完备性方程（completeness equation）

单量子比特计算基下有两个测量算子，分别是$M_0=|0⟩⟨0|$，$M_1=|1⟩⟨1|$

这两个测量算子都是自拌的，也就是转置共轭等于其本身，这两个算子同时也是幂等矩阵。

所以 $M_0^{†}{M}_0+M_0^{†}{M}_0=M_0+M_1=I$

设系统被测量时的状态是$|\psi ⟩=\alpha |0⟩+\beta |1⟩$

测量结果为0的概率是
$$p(0)=⟨\psi |M_0^{†}{M}_0|\psi ⟩=⟨\psi |M_0|\psi ⟩=|\alpha {|}^2$$
测量后系统状态变为
$$\frac{M_0|\psi ⟩}{\sqrt{⟨\psi |M_0^{†}{M}_0|\psi ⟩}}=\frac{M_0|\psi ⟩}{|\alpha {|}^2}=\frac{\alpha }{|\alpha |}|0⟩$$
测量结果为1的概率与测量后状态也同理。

下面我们来介绍一下量子观测。

量子观测方式有很多种，比如投影测量（projective measurements）、POVM（positive operator-valued measure）。

本文的最后我们来介绍一下投影变换。当测量算子具有酉变换性质，投影测量与一般测量等价。

投影测量描述一个可观测量$\Lambda$，$\Lambda$是带观测系统的状态空间上的自拌算子。同时$\Lambda$可以分解为

$$\Lambda =\sum _i{\lambda }_i{P}_i$$

这里$P_i$是$\Lambda$的特征值${\lambda }_i$对应特征空间上的投影。测量可能结果对应${\lambda }_i$。在对状态$|\psi ⟩$测量后，得到$i$的结果的概率为
$$p_i=p(\lambda ={\lambda }_i)=⟨\psi |P_i|\psi ⟩$$

测量后，结果$i$发生，量子态坍缩，量子系统的最新状态为
$$\frac{P_i|\psi ⟩}{\sqrt{p_i}}$$

投影测量有一个重要的特征就是很容易计算投影测量的平均值$E(\Lambda )$
$$\begin{gathered} E(\Lambda )=\sum _i{\lambda }_i{p}_i \\ =\sum _i{\lambda }_i⟨\psi |P_i|\psi ⟩ \\ =⟨\psi |\left(\sum _i{\lambda }_i{P}_i\right)|\psi ⟩ \\ =⟨\psi |\Lambda |\psi ⟩ \end{gathered}$$
这样计算就简单多了，观测量$\Lambda$也记作$⟨\Lambda ⟩\equiv ⟨\psi |\Lambda |\psi ⟩$

那么，标准差$\Delta (\Lambda )$也可以计算了
$${\left[\Delta (\Lambda )\right]}^2=⟨(\Lambda -⟨\Lambda ⟩{)}^2⟩=⟨{\Lambda }^2⟩-⟨\Lambda {⟩}^2$$

标准差是一个刻画典型分散程度的度量。

复合系统与联合测量

张量积

张量积是两个向量空间形成的一个更大向量空间的运算。在量子力学中，量子的状态由希尔伯特空间（Hilbert space）中的单位向量来描述。

设 $H_1$ 和 $H_2$ 分别为 $n_1$ 和 $n_2$ 维的希尔伯特空间。$H_1$ 和 $H_2$ 的张量积为一个 $mn$ 维的希尔伯特空间 $H\equiv H_1\otimes H_2$，对于 $H_1$ 中的每一个向量 $|h_1⟩$ 和 $H_2$ 中的每一个向量 $|h_2⟩$ 在 $H$ 中都要有唯一的向量 $|h_1⟩\otimes |h_2⟩$，并在 $H$ 中的向量都可以表示为 $|h_1⟩\otimes |h_2⟩$ 的线性叠加。还要满足以下基本性质：

1. 对于任意 $|h_1⟩\in H_1$，$|h_2⟩\in H_2$，以及任意复数$c\in \mathbb{C}$，都有
   $$c(|h_1⟩\otimes |h_2⟩)=(c|h_1⟩)\otimes |h_2⟩-|h_1⟩\otimes (c|h_2⟩)$$
2. 对于任意 $|h_1^1⟩,|h_1^2⟩\in H_1$，任意$|h_2⟩\in H_2$，都有
   $$(|h_1^1⟩+|h_1^2⟩)\otimes |h_2⟩=h_1^1⟩\otimes |h_2⟩+|h_1^2⟩\otimes |h_2⟩$$
3. 对于任意$|h_1⟩\in H_1$，任意 $|h_2^1⟩,|h_2^2⟩\in H_2$，都有
   $$|h_1⟩\otimes (|h_2^1⟩+|h_2^2⟩)=h_1⟩\otimes |h_2^1⟩+|h_1⟩\otimes |h_2^2⟩$$

复合系统的状态演化

（未完待续）