這篇文章提出了一種可以均勻隨機採樣Clifford算子的方法,相比於此前的採樣算法,該方法實現更為簡單,且可以直接生成Clifford算子對應的量子電路,同時在時間效率上也達到了目前最好的水平。
題目 : A simple method for sampling random Clifford operators
作者:Ewout van den Berg
Scite: 41次
ArXiv: 2008.06011
Clifford群在量子信息與量子計算領域中都具有重要地位,在量子計算的經典模擬[1]、量子信道噪聲刻畫[2]等方面都有著重要的應用,而其中某些應用要求我們能夠高效地從Clifford群中進行均勻隨機採樣。在這篇文章中,作者基於Pauli矩陣的Tableau表示方法,提出了一種能夠對Clifford群進行均勻隨機採樣的算法。下面我們將介紹相關的一些數學概念與文章算法的具體過程。
對於單個qubit而言,Pauli群被定義為所有帶符號( )的Pauli矩陣的集合,
這一集合在矩陣乘法運算下構成一個群;
-qubit Pauli群則由Pauli矩陣的張量積定義,即
-qubit Cliffird群 則被定義為所有能夠正規化(normalize)Pauli群的酉矩陣,即Clifford群中的任何一個元素都可以由一個僅包含Hadamard門、相位門以及控制非門的量子電路 實現,反之亦然[3]。由於任意兩個Pauli算子 , 的乘積可以表示為對應qubit上Pauli算子乘積的張量積,即 。於是可以令 表示張量積 ,其中 而其餘的 ;類似地,我們可以定義 ;那麼就可以用兩個 維向量 來表示Pauli算子,
其中,因子 是來自變換 。使用這種方式,我們就可以用一個 的二元矩陣來表示 個Pauli算子,例如 , 與 的Tableau表示就如圖一所示,另外表中也可以添加一列符號位 來表示Pauli算子的相位 ,文章為了行文簡潔忽略了這一列,但在之後的所有討論中相位列實際都存在。
圖1. , 與 的Tableau表示
前文提到了Clifford群可以通過單比特Hadamard門、相位門以及雙比特控制非門生成,而Tableau表示的好處就在於它能夠很方便地表示Pauli算子在Clifford算子作用下的變化:正如圖2所示,作用在qubit 上的Hadamard門會交換表中對應的兩列(因為Hadamard門會將同一個qubit上的X和Z互換),而相位門則會將 模2加到 上,最後作用在qubit 與 上的控制非門會把 加到 上,把 加到 上(這裡的加都是模2加)。
圖2. Hadamard門、相位門與控制非門作用圖示算法的主要思路是將隨機Clifford算子的採樣問題轉換為採樣一個隨機Pauli算子的Tableau表示。
給定Clifford算子 以及任意兩個Pauli算子 與 ,有
而任意一個Pauli算子都可以寫成 與 乘積的形式,所以只要知道了所有的 與 ,就可以確定 作用在其餘的Pauli算子上的結果;換言之,只需要確定Clifford算子 在這 個基項上的映射,就可以完全描述 。
假設 為初始表格(如圖3所示),那麼任意一張 的表就表示了某個對應的Clifford算子在 這 個基項上作用的結果,於是就把均勻隨機採樣Clifford算子轉換成了均勻隨機採樣Pauli算子的Tableau表示。
圖3. 3-qubit Pauli算子的基項
在以上轉換的基礎上,算法首先採樣一個隨機Pauli算子的Tableau表示,然後通過上文提到的Hadamard門,相位門與控制非門等操作將這張表變換為初始的Tableau表示,那麼由這些量子比特門構成的電路就是我們想要採樣的Clifford算子。
為了降低採樣的複雜度。本文採用了和文獻[5]類似的層次採樣技術 ,即每次只採樣表中兩行,先將這部分還原為初始的表(類似於高斯消元),下一次則在此基礎上再採樣兩行,重複以上操作,經過 次採樣就可以得到結果,且每經過一次採樣,樣本空間都會隨之縮小,操作過程如圖4所示(這裡採樣的兩行只需要滿足反對易即可,而X、Z反對易天然滿足條件)
圖4. 層次採樣過程圖示
算法中最重要的環節就是如何進行高效地採樣,作者提出了如下算法
對於第一行,隨機採樣一個非單位元的Pauli算子。由於每個Pauli算子都可以由一個 位的二元向量表示,所以可以從 中隨機採樣一個整數,並將它的二進位表示賦給表,同時再隨機採樣一個符號。令 為第一行中第一個非零元的下標。對第二行中除下標為 外的其餘所有元素進行隨機採樣,並令第 個元素為0,同時也隨機採樣一個符號。為了保證這兩行是反對易的,需要作如下映射:若這兩行的Pauli算子已經是反對易的,則保持所有元素不變;若這兩行不是反對易的,不妨設第二行的第 個元素為 ,則將 與 取反,即可改變這兩行Pauli算子的對易性(圖5)。
相比於此前的一些工作[4,5],本文工作的主要優勢在於實現簡單,並且可以按照流式方式直接生成採樣得到的Clifford算子對應的量子電路,電路中的不同部分在算法中也可以並行生成;同時算法的時間複雜度為 ,不輸目前最好的算法;此外,可以證明算法生成的電路中最多包含 個量子比特門,且電路的最大深度為 。
[1] Aaronson S, Gottesman D. Improved simulation of stabilizer circuits[J]. Physical Review A, 2004, 70(5).
[2] Emerson J, Alicki R, Życzkowski K, et al. Scalable noise estimation with random unitary operators[J]. Journal of Optics B-quantum and Semiclassical Optics, 2005, 7(10).
[3] Nielsen M A, Chuang I. Quantum computation and quantum information[J]. 2002.
[4] Koenig R, Smolin J A. How to efficiently select an arbitrary Clifford group element[J]. Journal of Mathematical Physics, 2014, 55(12).
[5] Bravyi S, Maslov D. Hadamard-free circuits expose the structure of the clifford group[J]. arXiv preprint arXiv:2003.09412, 2020.
封面圖片來源於: https://www.inverse.com/article/30729-quantum-code-programming-quantum-mechanics-quantum-computers