量子情報入門 第2章 §2.6
CSS符号よりもさらに広いクラスの符号を統一的に扱うのがスタビライザー形式 である。その中心的な着想は「状態ベクトルそのものではなく、状態を固定する演算子の集合で量子状態を記述する 」というものだ。この視点の転換は単に符号の記述を簡潔にするだけでなく、誤り・ゲート・測定のすべてを統一的な代数の言葉で扱うことを可能にし、さらにGottesman–Knill定理という量子計算の本性に関わる副産物まで生む。本章では、必要な群論の言葉を最小限だけ準備したうえで、形式の全体を一歩ずつ組み上げる。
準備:群論の最小限
以下で使う群論は次の3つの概念だけである。
(i) 群と部分群 。集合 G G G に結合的な積が定義され、単位元と逆元を持つとき G G G を群 という。G G G の部分集合 S S S が同じ積で群をなすとき部分群 という。本章の群はすべて行列の積に関する行列の集合である。
(ii) 生成元 。群 G G G の元の列 g 1 , … , g l g_1,\ldots,g_l g 1 , … , g l の積(繰り返し・順不同を許す)で G G G のすべての元が書けるとき、g 1 , … , g l g_1,\ldots,g_l g 1 , … , g l は G G G を生成する といい、G = ⟨ g 1 , … , g l ⟩ G=\langle g_1,\ldots,g_l\rangle G = ⟨ g 1 , … , g l ⟩ と書く。生成元による記述の利点は圧縮 である:大きさ ∣ G ∣ |G| ∣ G ∣ の群は高々 log 2 ∣ G ∣ \log_2|G| log 2 ∣ G ∣ 個の生成元で生成できる(新しい生成元を1つ加えるたびに生成される群のサイズは少なくとも2倍になるため)。さらに、あるベクトルが S S S のすべての元で固定されることを確かめるには生成元で固定されることだけ 確かめれば十分である——積で書ける元による固定は自動的に従う。この2点が、以下で 2 n 2^n 2 n 次元の状態を n n n 個程度の演算子で記述できるカラクリの正体である。
(iii) 共役 。U g U † U g U^\dagger U g U † の形の変換を(U U U による)共役 という。ユニタリ発展のもとでの演算子の変換がちょうどこの形になることを、次章で徹底的に使う。
Pauli群
主役となる群はPauli群 である。1量子ビットのPauli群は、Pauli行列に ± 1 , ± i \pm1,\pm i ± 1 , ± i の係数を許した16個の行列の集合
G 1 ≡ { ± I , ± i I , ± X , ± i X , ± Y , ± i Y , ± Z , ± i Z } G_1 \equiv \{\pm I,\pm iI,\ \pm X,\pm iX,\ \pm Y,\pm iY,\ \pm Z,\pm iZ\} G 1 ≡ { ± I , ± i I , ± X , ± i X , ± Y , ± iY , ± Z , ± i Z }
である。「なぜ ± 1 , ± i \pm1,\pm i ± 1 , ± i を付けるのか」と思うかもしれないが、これは積で閉じるために必要である:たとえば X Y = i Z XY=iZ X Y = i Z 、Z X ⋅ X Z = − Y ⋅ ( − Y ) ZX \cdot XZ = -Y\cdot(-Y) Z X ⋅ X Z = − Y ⋅ ( − Y ) 等、Pauli行列同士の積には ± i \pm i ± i が自然に発生するので、係数を許さないと群にならない。n n n 量子ビットのPauli群 G n G_n G n は、Pauli行列の n n n 重テンソル積に同じく ± 1 , ± i \pm1,\pm i ± 1 , ± i の係数を許したものとして定義する。G n G_n G n の重要な性質は、任意の2元が可換または反可換のいずれかである こと(各テンソル因子の可換・反可換の符号の積で決まる)、および係数 ± i \pm i ± i の付かない元は g 2 = I g^2=I g 2 = I を満たすことである。
スタビライザーと固定部分空間
定義 (スタビライザー)
S S S を G n G_n G n の部分群とする。S S S のすべての元 g g g に対して g ∣ ψ ⟩ = ∣ ψ ⟩ g\ket\psi=\ket\psi g ∣ ψ ⟩ = ∣ ψ ⟩ を満たす状態 ∣ ψ ⟩ \ket\psi ∣ ψ ⟩ の全体を V S V_S V S と書き、S S S にスタビライズされる部分空間 と呼ぶ。S S S を V S V_S V S のスタビライザー という。
V S V_S V S が部分空間になること(∣ ψ ⟩ , ∣ ϕ ⟩ ∈ V S \ket\psi,\ket\phi\in V_S ∣ ψ ⟩ , ∣ ϕ ⟩ ∈ V S なら任意の線形結合も V S V_S V S )、および V S V_S V S が各元の固有値 + 1 +1 + 1 固有空間の共通部分であることは、定義から直ちに確かめられる。
例1:EPR状態
∣ ψ ⟩ = ( ∣ 00 ⟩ + ∣ 11 ⟩ ) / 2 \ket\psi=(\ket{00}+\ket{11})/\sqrt2 ∣ ψ ⟩ = ( ∣ 00 ⟩ + ∣ 11 ⟩ ) / 2 は X 1 X 2 ∣ ψ ⟩ = ∣ ψ ⟩ X_1X_2\ket\psi=\ket\psi X 1 X 2 ∣ ψ ⟩ = ∣ ψ ⟩ かつ Z 1 Z 2 ∣ ψ ⟩ = ∣ ψ ⟩ Z_1Z_2\ket\psi=\ket\psi Z 1 Z 2 ∣ ψ ⟩ = ∣ ψ ⟩ を満たす(直接計算)。実はこの2つの演算子で固定される状態は(全体位相を除き)∣ ψ ⟩ \ket\psi ∣ ψ ⟩ ただ一つ である。すなわち「2 2 2^2 2 2 個の複素振幅」の代わりに「2個の演算子 X 1 X 2 , Z 1 Z 2 X_1X_2, Z_1Z_2 X 1 X 2 , Z 1 Z 2 」でこの状態は完全に指定できる。
例2:V S V_S V S を生成元から決める
n = 3 n=3 n = 3 、S = ⟨ Z 1 Z 2 , Z 2 Z 3 ⟩ S=\langle Z_1Z_2,\ Z_2Z_3\rangle S = ⟨ Z 1 Z 2 , Z 2 Z 3 ⟩ とする(S S S の元は I , Z 1 Z 2 , Z 2 Z 3 , Z 1 Z 3 I, Z_1Z_2, Z_2Z_3, Z_1Z_3 I , Z 1 Z 2 , Z 2 Z 3 , Z 1 Z 3 の4つ;たとえば Z 1 Z 3 = ( Z 1 Z 2 ) ( Z 2 Z 3 ) Z_1Z_3=(Z_1Z_2)(Z_2Z_3) Z 1 Z 3 = ( Z 1 Z 2 ) ( Z 2 Z 3 ) )。Z 1 Z 2 Z_1Z_2 Z 1 Z 2 の + 1 +1 + 1 固有空間は「第1・2ビットのパリティが偶」の空間で ∣ 000 ⟩ , ∣ 001 ⟩ , ∣ 110 ⟩ , ∣ 111 ⟩ \ket{000},\ket{001},\ket{110},\ket{111} ∣ 000 ⟩ , ∣ 001 ⟩ , ∣ 110 ⟩ , ∣ 111 ⟩ が張る。Z 2 Z 3 Z_2Z_3 Z 2 Z 3 の + 1 +1 + 1 固有空間は ∣ 000 ⟩ , ∣ 100 ⟩ , ∣ 011 ⟩ , ∣ 111 ⟩ \ket{000},\ket{100},\ket{011},\ket{111} ∣ 000 ⟩ , ∣ 100 ⟩ , ∣ 011 ⟩ , ∣ 111 ⟩ が張る。共通部分は
V S = s p a n { ∣ 000 ⟩ , ∣ 111 ⟩ } V_S = \mathrm{span}\{\ket{000},\ \ket{111}\} V S = span { ∣ 000 ⟩ , ∣ 111 ⟩ }
——3量子ビットビット反転符号の符号空間そのもの である。第10章で「Z 1 Z 2 Z_1Z_2 Z 1 Z 2 と Z 2 Z 3 Z_2Z_3 Z 2 Z 3 を測る」と述べた症候測定の正体は、スタビライザー生成元の測定だったのである。
V S V_S V S が非自明であるための条件
どんな部分群 S S S でも役に立つわけではない。たとえば S = { ± I , ± X } S=\{\pm I,\pm X\} S = { ± I , ± X } を考えると、− I ∣ ψ ⟩ = ∣ ψ ⟩ -I\ket\psi=\ket\psi − I ∣ ψ ⟩ = ∣ ψ ⟩ を満たすのは ∣ ψ ⟩ = 0 \ket\psi=0 ∣ ψ ⟩ = 0 だけなので V S = { 0 } V_S=\{0\} V S = { 0 } (自明)。V S ≠ { 0 } V_S\neq\{0\} V S = { 0 } のために必要な条件が次の2つである:
(a) S S S の元はすべて可換 である。
(b) − I ∉ S -I\notin S − I ∈ / S 。
必要性の証明は短い。V S V_S V S に零でない ∣ ψ ⟩ \ket\psi ∣ ψ ⟩ があるとし、M , N ∈ S M,N\in S M , N ∈ S が反可換だと仮定すると
− ∣ ψ ⟩ = − N M ∣ ψ ⟩ = M N ∣ ψ ⟩ = ∣ ψ ⟩ -\ket\psi = -NM\ket\psi = MN\ket\psi = \ket\psi − ∣ ψ ⟩ = − N M ∣ ψ ⟩ = M N ∣ ψ ⟩ = ∣ ψ ⟩
(M , N M,N M , N がともに ∣ ψ ⟩ \ket\psi ∣ ψ ⟩ を固定することを両端で使った)。よって ∣ ψ ⟩ = 0 \ket\psi=0 ∣ ψ ⟩ = 0 となり矛盾——G n G_n G n の2元は可換か反可換しかないので、(a)が従う。(b)は − I ∣ ψ ⟩ = ∣ ψ ⟩ -I\ket\psi=\ket\psi − I ∣ ψ ⟩ = ∣ ψ ⟩ が直ちに矛盾を与える。実はこの2条件は十分 でもあり、次小節で次元まで込めて示す。なお − I ∉ S -I\notin S − I ∈ / S から ± i I ∉ S \pm iI\notin S ± i I ∈ / S 、さらに S S S のすべての元が g 2 = I g^2=I g 2 = I (したがって g † = g g^\dagger=g g † = g 、Hermite)を満たすことも従う。以降、スタビライザーは常に「独立で可換な生成元により生成され − I ∉ S -I\notin S − I ∈ / S 」という規約で扱う。
検査行列表現
生成元の独立性(g i g_i g i をどれか1つ除くと生成される群が真に小さくなること)や可換性を機械的に判定する道具が検査行列 である。g ∈ G n g\in G_n g ∈ G n (係数因子は無視)に対し 2 n 2n 2 n 次元の行ベクトル r ( g ) r(g) r ( g ) を、
第 j 成分 = { 1 g の第 j ビットが X か Y 0 それ以外 , 第 n + j 成分 = { 1 g の第 j ビットが Z か Y 0 それ以外 \text{第 } j \text{ 成分} = \begin{cases}1 & g \text{ の第 } j \text{ ビットが } X \text{ か } Y\\ 0&\text{それ以外}\end{cases},\qquad
\text{第 } n{+}j \text{ 成分} = \begin{cases}1 & g \text{ の第 } j \text{ ビットが } Z \text{ か } Y\\ 0&\text{それ以外}\end{cases} 第 j 成分 = { 1 0 g の第 j ビットが X か Y それ以外 , 第 n + j 成分 = { 1 0 g の第 j ビットが Z か Y それ以外
で定める(左半分が「X X X 成分」、右半分が「Z Z Z 成分」;Y = i X Z Y=iXZ Y = i X Z は両側に 1 1 1 )。生成元 g 1 , … , g l g_1,\ldots,g_l g 1 , … , g l を並べた l × 2 n l\times2n l × 2 n 行列が検査行列である。この表現の要点は2つ:
(i) 積は和になる :r ( g ) + r ( g ′ ) = r ( g g ′ ) r(g)+r(g')=r(gg') r ( g ) + r ( g ′ ) = r ( g g ′ ) ( m o d 2 \bmod\,2 mod 2 ;係数因子を除く)。群の積がベクトルの加法に写る。
(ii) 可換性は「ねじれた内積」で判る :2 n × 2 n 2n\times2n 2 n × 2 n 行列
Λ = ( 0 I I 0 ) ( I は n × n ) \Lambda = \begin{pmatrix}0&I\\ I&0\end{pmatrix}
\qquad(I \text{ は } n\times n) Λ = ( 0 I I 0 ) ( I は n × n )
を使うと、g , g ′ g,g' g , g ′ が可換 ⟺ r ( g ) Λ r ( g ′ ) T = 0 ( m o d 2 ) \iff r(g)\,\Lambda\,r(g')^T = 0 \pmod 2 ⟺ r ( g ) Λ r ( g ′ ) T = 0 ( mod 2 ) 。左辺は「g g g の X X X と g ′ g' g ′ の Z Z Z が同じ場所で衝突する回数+その逆」の偶奇を数えており、X X X と Z Z Z の衝突1回ごとに反可換の符号が1つ出ることの帳簿である(Steane符号の6生成元・全ペアについてこの判定則が実際の行列の可換性と一致することを数値確認済み)。
これらから、独立性の判定が線形代数に落ちる:
命題
− I ∉ S = ⟨ g 1 , … , g l ⟩ -I\notin S=\langle g_1,\ldots,g_l\rangle − I ∈ / S = ⟨ g 1 , … , g l ⟩ のとき、g 1 , … , g l g_1,\ldots,g_l g 1 , … , g l が独立 ⟺ \iff ⟺ 検査行列の行が m o d 2 \bmod\,2 mod 2 で線形独立。
(i)より行の従属 ∑ i a i r ( g i ) = 0 \sum_ia_ir(g_i)=0 ∑ i a i r ( g i ) = 0 は ∏ i g i a i = ± I \prod_ig_i^{a_i}=\pm I ∏ i g i a i = ± I を意味し、− I ∉ S -I\notin S − I ∈ / S より = I =I = I 、すなわちある生成元が他の積で書けることと同値だからである。Steane符号の検査行列(第13章の H H H を左右ブロックに配したもの)は m o d 2 \bmod\,2 mod 2 でランク6であり、6生成元の独立性が機械的に確認できる(数値確認済み)。
固定空間の次元:2 k 2^k 2 k 次元定理
命題 2.4
S = ⟨ g 1 , … , g l ⟩ S=\langle g_1,\ldots,g_l\rangle S = ⟨ g 1 , … , g l ⟩ 、生成元は独立、− I ∉ S -I\notin S − I ∈ / S とする。任意の i i i に対し、g g i g † = − g i gg_ig^\dagger=-g_i g g i g † = − g i かつ g g j g † = g j gg_jg^\dagger=g_j g g j g † = g j (j ≠ i j\neq i j = i )を満たす g ∈ G n g\in G_n g ∈ G n が存在する。
つまり「狙った生成元1つだけの符号を反転させる 」Pauli演算子が必ず見つかる。証明は検査行列の面目躍如で、行が独立なら G Λ x = e i G\Lambda x=e_i G Λ x = e i (e i e_i e i は第 i i i 単位ベクトル)を満たす x x x が解け、r ( g ) = x T r(g)=x^T r ( g ) = x T なる g g g が所望のものになる。この命題は後で、症候の符号合わせ・論理演算子の構成・次元の証明と、繰り返し使われる。
命題
S = ⟨ g 1 , … , g n − k ⟩ S=\langle g_1,\ldots,g_{n-k}\rangle S = ⟨ g 1 , … , g n − k ⟩ (n − k n-k n − k 個の独立・可換な生成元、− I ∉ S -I\notin S − I ∈ / S )のとき、V S V_S V S は 2 k 2^k 2 k 次元である。
証明の鍵は射影子の明示式である。( I + g j ) / 2 (I+g_j)/2 ( I + g j ) /2 は g j g_j g j の + 1 +1 + 1 固有空間への射影子なので、V S V_S V S への射影子は
P S = ∏ j = 1 n − k I + g j 2 P_S = \prod_{j=1}^{n-k}\frac{I+g_j}{2} P S = j = 1 ∏ n − k 2 I + g j
と書ける。より一般に符号ベクトル x ∈ { 0 , 1 } n − k x\in\{0,1\}^{n-k} x ∈ { 0 , 1 } n − k に対し P S x = ∏ j I + ( − 1 ) x j g j 2 P_S^x=\prod_j\frac{I+(-1)^{x_j}g_j}{2} P S x = ∏ j 2 I + ( − 1 ) x j g j とおくと、命題2.4 により各 P S x P_S^x P S x は P S = P S 0 ⋯ 0 P_S=P_S^{0\cdots0} P S = P S 0 ⋯ 0 とユニタリ共役(同じ次元)で、異なる x x x どうしは直交し、しかも ∑ x P S x = I \sum_xP_S^x=I ∑ x P S x = I 。2 n − k 2^{n-k} 2 n − k 個の同次元の直交射影子が 2 n 2^n 2 n 次元を埋め尽くすのだから、各々の次元は 2 n / 2 n − k = 2 k 2^n/2^{n-k}=2^k 2 n / 2 n − k = 2 k 。直観的には「独立な生成元を1つ課すたびに空間が半分になる 」(± 1 \pm1 ± 1 固有空間が半々)ということである。
これで役者が揃った:n − k n-k n − k 個の生成元は 2 k 2^k 2 k 次元の符号空間を定め、そこに k k k 量子ビットを符号化できる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com