单量子门分解(使用有限的单量子门来近似模拟任意的单量子门)

UTF8gbsn

Introduction

首先我们介绍一下在量子计算集中单量子比特U门。它为一个定义在复数域上的酉矩阵。

$$\left(\begin{array}{cc}{c}_{00}& c_{01}\\ c_{10}& c_{11}\end{array}\right),U{U}^{†}=I$$

通过定义我们可以发现，这样的U门有无数个。但是在实现量子计算机的时候，工程上我
们只能实现有限个单qubit的U门。这样就出现了，使用有限的单量子门来近似模拟任意的
单量子门。其中$U^{†}$是U的共轭转置矩阵。

本文主要根据《Quantum Computation and Quantum Information》
书籍中的内容总结而
来。是一个简化的读书笔记。其中部分细节不明之处还望见谅。可以参考书籍或留言。

Brief Review

Question

本文假设，在工程上实现了两种单qubit门，H和T。这两种门的如下， $$H=\left(\begin{array}{cc}1/\sqrt{2}& 1/\sqrt{2}\\ 1/\sqrt{2}& -1/\sqrt{2}\end{array}\right),T=\left(\begin{array}{cc}{e}^{-i\frac{\pi }{8}}& 0\\ 0& e^{i\frac{\pi }{8}}\end{array}\right)$$

而我们需要使用这两种门来近似任意的U门。当然这种近似模拟，是有一定的误差的。
$$\left(\begin{array}{cc}{c}_{00}& c_{01}\\ c_{10}& c_{11}\end{array}\right),U{U}^{†}=I$$

Steps Overview

1. 首先，我们先来陈述一个定理。我们先不去证明它。只是把它作为已知的内容。这
   样可以把脉络梳理得简单清楚。

   $$U=R_{\hat{n}}(\beta )R_{\hat{m}}(\gamma )R_{\hat{n}}(\delta )$$

   我们来解释一下，这个定理。给定三维空间中的任意两个不平行的轴$\hat{m},\hat{n}$，那么任意的U门可以分解为绕这个两个轴的三个旋转门。其中
   $R_{\hat{n}}(\beta )$表示绕$\hat{n}$,旋转$\beta$。且$R$是一些U门,也就
   是说他们都是一些单量子门。

2. 我们使用H，T来近似构建，任意旋转角度的$R_{\hat{n}}(\theta ),R_{\hat{m}}(\beta )$,其中
   $\hat{n}\simeq (cos(\pi /8),sin(\pi /8),cos(\pi /8)),\hat{m}\simeq (cos(\pi /8),-sin(\pi /8),cos(\pi /8))$，后面我们会讲为什么是近似。如果我们能
   近似模拟任意旋转角度的$R_{\hat{n}}(\theta ),R_{\hat{m}}(\beta )$。我们
   就可一个根据上面的定理，实现任意的U。这就是基本的逻辑，那么下面我们来看如何
   实现任意旋转角度的$R_{\hat{n}}(\theta ),R_{\hat{m}}(\beta )$。

Rotate

这一小节我们来看看如何近似实现任意角度的$R_{\hat{n}}(\theta ),R_{\hat{m}}(\beta )$

Paulli Matrix

首先我们先来介绍一下几个常用的U门

$$X=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right),Y=\left(\begin{array}{cc}0& -i\\ i& 0\end{array}\right),Z=\left(\begin{array}{cc}1& 0\\ 0& -1\end{array}\right)$$

根据这些U门，我们可以来写出绕$x,y,z$轴旋转的U矩阵，如下(还是先把它们作为已知的
公式，有空我们在写证明过程。)

$$R_x(\theta )=e^{-i\theta X/2}=cos(\theta /2)I-i\cdot sin(\theta /2)X=\left(\begin{array}{cc}cos(\frac{\theta }{2})& -i\cdot sin(\frac{\theta }{2})\\ -i\cdot sin(\frac{\theta }{2})& cos(\frac{\theta }{2})\end{array}\right)$$
$$R_y(\theta )=e^{-i\theta Y/2}=cos(\theta /2)I-i\cdot sin(\theta /2)Y=\left(\begin{array}{cc}cos(\frac{\theta }{2})& -sin(\frac{\theta }{2})\\ sin(\frac{\theta }{2})& cos(\frac{\theta }{2})\end{array}\right)$$
$$R_z(\theta )=e^{-i\theta Z/2}=cos(\theta /2)I-i\cdot sin(\theta /2)Z=\left(\begin{array}{cc}{e}^{-i\theta /2}& 0\\ 0& e^{i\theta /2}\end{array}\right)$$

有了这些概念之后，我们还需要写出一个公式它表示绕任意一个轴$\hat{n}=(n_x,n_y,n_z)$旋转$\theta$的U门。我们仍然
把它作为已知的东西。证明放在以后来写。

$$R_{\hat{n}}=e^{-i\theta \hat{n}\cdot \overrightarrow{\sigma }/2}=cos(\frac{\theta }{2})I-i\cdot sin(\frac{\theta }{2})(n_{x}X+n_{y}Y+n_{z}Z)$$
其中，$\overrightarrow{\sigma }=(X,Y,Z)$，它是一个向量。

n Rotate

我们首先来构建$R_{\hat{n}}(\theta )$,让我们来一步一步的算。怎么靠H，T来构
建$R_{\hat{n}}(\theta )$。

1. 首先让我们先来看看两个辅助的U门。

   • $$exp(-i\frac{\pi }{8}Z)=cos(\frac{\pi }{8})I-i\cdot sin(\frac{\pi }{8})Z=\left(\begin{array}{cc}{e}^{-i\frac{\pi }{8}}& 0\\ 0& e^{i\frac{\pi }{8}}\end{array}\right)$$

     这个门就是绕Z轴旋转了$\frac{\pi }{4}$[1],这个门实际上就是T门

   • $$exp(-i\frac{\pi }{8}X)=cos(\frac{\pi }{8})I-i\cdot sin(\frac{\pi }{8})X=HTH$$

2. 有了上面这连个门。我们可以进一步做下面的变换

   $$exp(-i\frac{\pi }{8}Z)exp(-i\frac{\pi }{8}X)=(cos(\frac{\pi }{8})I-i\cdot sin(\frac{\pi }{8})Z)(cos(\frac{\pi }{8})I-i\cdot sin(\frac{\pi }{8})X)$$
   $$=co{s}^2(\frac{\pi }{8})I-i[cos(\frac{\pi }{8})X+sin(\frac{\pi }{8})Y+cos(\frac{\pi }{8})Z]sin(\frac{\pi }{8})$$
   $$exp(-i\frac{\pi }{8}Z)exp(-i\frac{\pi }{8}X)=T(HTH)=THTH$$

   参考3.1中我们发现实际上$exp(-i\frac{\pi }{8}Z)exp(-i\frac{\pi }{8}X)$是一
   个绕某个轴转动一定角度的U门这个轴近似是$(cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$,旋转的角度$\theta$为
   $cos(\theta /2)=co{s}^2(\frac{\pi }{8})$。注意我们准确的表示出了$\theta$但是没有
   准确的表示出$\hat{n}$,是因为系数
   $co{s}^2(\frac{\pi }{8}),sin(\frac{\pi }{8})$，不是一对正弦余弦值$(co{s}^2+si{n}^2=1)$。
   所以，后面$sin(\frac{\pi }{8})$需要同$(cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$,配凑出
   $sin(\frac{\theta }{2})=\sqrt{1-co{s}^4(\frac{\pi }{8})}$所以这条轴不是$(cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$，但是和这条轴接近。

3. 我们总结一下，$U=THTH$门的作用。它可以绕近似$\hat{n}=(cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$,
   一个角度$\theta =2arccos(co{s}^2(\pi /8))$。

4. 我们现在来使用$R_{\hat{n}}(\theta )$近似$R_{\hat{n}}(\alpha )$。其
   中$\theta =2arccos(co{s}^2(\pi /8))$,$\hat{n}=(cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$,$\alpha$是任意角度。注意这个部分是
   关键。可能我这里写得不是很清楚。如果不是很明白，请参考书籍和相关其它资料或留
   言。

   • 对于任意的$\delta >0,N=\frac{2\pi }{\delta }$,我们定义${\theta }_k\in [0,2\pi ),{\theta }_k=(k\theta )mod2\pi$.
     那么由抽屉原理可得，存在不同的$j,k$使得
     $$|{\theta }_k-{\theta }_j|⩽\frac{2\pi }{N}<\delta$$
     不失去一般性令$k>j$，我们可以从这里可以推导出${\theta }_{k-j}<\delta$。这里如
     果看不明白，可以稍微画图思考下或留言。

   • 这样，我们可以看到${\theta }_{l(k-j)}$的序列，其中$l\in (0,1,2,...)$。这
     个序列将会覆盖掉整个$2\pi$,而序列的相邻两项之间的差不大于$\delta$。即
     ${\theta }_{l(k-j)}$可以以精度$\delta$逼近$2\pi$上的所有位置。

   • 也就是说存在某个N,使得$R_{\hat{n}}^N(\theta )$,在误差限
     $ϵ$下,逼近$R_{\hat{n}}(\alpha )$

到这里我们就完成了最重要的工作，使用T,H,来逼近任意的$R_{\hat{n}}(\alpha )$。

m Rotate

我们再来构建$R_{\hat{m}}(\alpha )$,我们稍微的使用H门来变换一下
$R_{\hat{n}}(\alpha )$,就可以得到$R_{\hat{m}}(\alpha )$.

$$H{R}_{\hat{n}}(\alpha )H=R_{\hat{m}}(\alpha )$$

而近似轴为，$\hat{m}\simeq (cos(\frac{\pi }{8}),sin(\frac{\pi }{8}),cos(\frac{\pi }{8}))$

Arbitrary Rotation

我们有了$R_{\hat{n}}(\alpha ),R_{\hat{m}}(\alpha )$,而且
$\hat{n},\hat{m}$不平行。那么我们就可以利用2.2中的定理来构建,任意的U门

$$U=R_{\hat{n}}(\beta )R_{\hat{m}}(\gamma )R_{\hat{n}}(\delta )$$

我们来描述一下具体步骤

1. 把U在数学上分解为
   $$U=R_{\hat{n}}(\beta )R_{\hat{m}}(\gamma )R_{\hat{n}}(\delta )$$

2. 然后根据部门前面提到的分别实现，我们用T，H来近似下面的U门。
   $$R_{\hat{n}}(\beta ),R_{\hat{m}}(\gamma ),R_{\hat{n}}(\delta )$$

3. 当然具体的步骤会很复杂。也可能有的们没有多项式复杂度的替换。所以关于任意
   门分解的问题还有很多其它内容。需要扩展阅读更多的书籍和论文。

The End

我们这里留下了很多坑，需要后面慢慢的填。因为很多公式定理都是直接给出。没有证明，
但是对于理解单U门的近似逼近已经做到了比较清楚的解释。

感谢

[1]weixin_41409816