量子計算理論(森前 著) の演習問題を解く part1
こんにちは。Kumaです。
最近、量子コンピュータについて勉強しています。
今回は有名な以下の本の演習問題について、解答が載っていないので一部書いてみたいとおもいます。
量子計算理論 量子コンピュータの原理 | 森北出版株式会社
二章はチューリングマシンから始まりますが、あまり詳しくないのでpp.15- (古典的確率状態)から書いていきます。
pp.15.1 Xゲートの確認
としたときに、 および を示せ。
- 解答
行列演算で示す。
および
なので
と示される。同様に
pp.16.1 トフォリゲートの確認
トフォリゲートを「三量子ビットに作用する演算であって、第一および第二量子ビットが1のときに限り第三量子ビットを反転させるもの」とする。即ち
ただしは恒等演算子であって である。
また加算記号はmod 2 の加算である。
このとき、任意のに対して を示せ。
- 解答
確認しながらやっていく。まずは
”は恒等演算子であって である。”
を理解しよう。
であれば納得できそうである。これを示そう。
ここで、 を思い出そう。(直交性)
すると
となり示された。
同様にして も示せるので、がたしかに恒等写像になっていることがわかった。
*1
次にトフォリゲートを理解しよう。
"
というゲートが「三量子ビットに作用する演算であって、第一も第二量子ビットも1のときに限り第三量子ビットを反転させるもの」といえることを確認する。"
具体的にみて納得しよう。例えば, などを確認したい。*2
初めに を確認する。
ここで、出てきた という演算子がちょっと難しいので解説する。
まず、テンソル積 の左右で独立と思って良い。つまりこの場合はテンソル積の左と右 でわけて考える。
テンソル積の左は第一量子ビットと第二量子ビットで作られる空間に対して作用する二量子ビット演算子である。
テンソル積の右は、第三量子ビットだけで作られる空間 に対して作用する一量子ビット演算子である。
つまり は、 *3かつ より
となる。
*4
よって
全く同様にして他の入力に対する演算も確かめられる。
例えば、
確かに入力の第一量子ビットも第二量子ビットも1,つまり のときだけ第三量子ビットが反転しそうである。
これを納得できれば*5、題意は簡単に示せる。
まずのとき、ab = 1 であり、 aかbのいずれかが0のときは ab = 0 である]
そしてあるbit に対して1とmod 2加算()を行うと、bitは反転する()。0であればbitは変化しない。
そのため、明らかに
"任意のに対して " を満たす事(トフォリゲートの別表現)がわかる。
今回はここまで。
光の偏波とパウリ行列の関係
この記事は日曜数学アドベントカレンダーの12/16ぶんです。
adventar.org
こんにちは。Kumaです。
この記事ではいよいよ「光の偏波」と「パウリ行列」の関係をみます。
光の偏波を記述するためのフレームとして
射影演算子とパウリ行列がとても自然であることを示します。
- エルミート行列のパウリ行列展開
- 射影演算子のパウリ行列展開の係数(”成分”)の性質
- 複素ベクトルのユニタリ変換
- 奇妙なことが起こっている
- 光の偏波との関係
- なぜスピノールというのか
- なぜいま偏波か
- まとめ
エルミート行列のパウリ行列展開
2x2エルミート行列は、パウリ行列 と単位行列 の結合で書くことができます。
ここで はすべて実数であって次の式を満たします。
これは以前の記事で証明しています。
electrodynamics.hatenablog.com
他方、二成分の複素ベクトル
に対して
は射影演算子といいます。
をと書くこともあります。
射影演算子は明らかにエルミート行列なので、パウリ行列展開が可能です。
射影演算子のパウリ行列展開の係数(”成分”)の性質
さらに射影演算子のパウリ行列展開における係数(”成分”)である実数を考えます。
ここでのノルムであるとします。するとがによらず成り立つ*1ので
は定数です。自由度ではないので無視しましょう。
結局、射影演算子を考えることでノルム1の2成分複素ベクトル は 実数のベクトル
に対応しました。
さらにこれらは独立ではありません。
それを確認するために、成分表示を定義から求めてみると
となります。上記からも
という関係式を満たすことがわかります。
すなわち、 は三次元空間における半径1の球面と1:1に対応します。
複素ベクトルのユニタリ変換
二成分複素ベクトルに戻って、のユニタリ変換を考えてみましょう。
ユニタリ変換とは、内積を保つ変換です。2つのベクトルの成す角度を変えないような変換であり、直感的には”回転”と思えます。
これはの世界ではどのように表現されるでしょうか。
ユニタリ変換は行列でかけます。一般には複素行列ですが簡単のために実数のものを取ってきましょう。
”回転”というイメージからもわかるように、次のものはユニタリ変換の一種です。
これをに作用させると、
となります。射影演算子の成分表示を求めると
となります。ここで倍角の公式を使っています。
これはよくみると、
とかけます。
は三次元幾何においてz軸まわりの回転操作に相当します。
すなわち
という変換を考えたときに対応するは
と移されます。
奇妙なことが起こっている
実数ベクトルの世界では回転が2倍になってみえているようでした。
すると、実数ベクトルの世界で一周、すなわち回転したときに複素ベクトルの世界ではしか回転していないわけです。
つまりある1つの実数ベクトルに対応する複素ベクトルは2つ存在することになります。
(と が同じ実数ベクトル を与えるという意味です。)
このような2:1の回転の関係にあるようなものをスピノール(Spin-or)といいます。
名前の由来は後で説明します。
こんなへんてこな対応を持っている世界を考えてなにかいいことがあるのでしょうか・・・
光の偏波との関係
その”へんてこな対応の世界”を実際にお見せします。
それこそが光の偏波状態といわれる概念です。
electrodynamics.hatenablog.com
偏光 - Wikipedia
以前の記事では時間依存性をで書いていますが、オイラーの公式からと書いても構いません。
すると光のベクトル表現は
とかけます。光の偏波状態とはtを媒介変数とみたときのベクトルの軌跡(直線、円、楕円)のことでした。
偏波状態に着目するのであれば、振幅の比だけに意味があり、位相の差だけに意味があります。
(軌跡自体の”大きさ”や”始点”には興味がないからです)
そこで、としても一般性を失いません。*2
さらには省略して書くことにしましょう。
すると、偏波状態は以下の手続きで書けます。
まずノルム1の二成分複素ベクトルを考えて
とします。これは光学の分野ではJonesベクトルと呼ばれます。*3
ジョーンズ計算法 - Wikipedia
更に
という行列を作ります。この行列の嬉しい点は、に任意の位相回転を与えても、それが複素共役でキャンセルされて
影響しないことです。すなわち”位相差にだけ着目したい”という偏波状態の要請を自然に満たします。
この行列をウォルフのコヒーレンシー行列といいます。
我々はこれとそっくり同じものを知っています。
そう、コヒーレンシー行列とは光のベクトルからできる射影演算子です!
射影演算子であれば、パウリ行列展開を考えることで実数のベクトルに対応させることができます。
それは、光学の世界ではStokesベクトルと呼ばれています。*4
その定義は
です。
Stokes parameters - Wikipedia
(Representations in fixed bases など)
確かに射影演算子のパウリ行列展開からできる実数ベクトルと全く同じですね!
これが球面と1:1対応することもすでにみました。
この球面は光学の世界ではポアンカレ球と呼ばれています!
ちなみに の作用の物理的な意味としては、例えば
- 光の測定器のほうを 回してから測定する
- xに対してy方向だけに位相差を与える
に相当します。
後者の操作は異方性を有する物体を通すことで可能です。異方性の物体中ではx方向とy方向で光の速度が異なり、位相差が生まれます。
これで、すべてが繋がってきました。光の偏波状態を書くものとしては光のベクトルそのものを考えるよりも
射影演算子を考えるべきで、それはパウリ行列展開という道具を介して実数のベクトルと球面座標でかけるのです。
これをJonesベクトル、コヒーレンシー行列、Stokesベクトル、ポアンカレ球などと呼ぶわけです。
しかし、複素ベクトル⇔実ベクトル対応は1つ難点があるのでした。それは2:1対応だということです。
つまり、とは同じ偏波状態(コヒーレンシー行列)を与え、もちろんStokesベクトルも同じになるのです。
偏波の場合、これは何を意味するのでしょうか?
以下の図をみるとわかります。
Jonesベクトルが回転すると正負がひっくり返ります。しかしこれは振動の初期位相がずれることと同じです。
すなわち、偏波状態としては変わらないのです。(偏波状態は軌跡=曲線のことで、その開始位置はどうでもよいのです)
だからと はその偏波状態を考えるときには同一視するべきで、2:1対応が自然になるのです。
逆に、回転以外は描く直線の角度自体が変わってしまうので偏波状態としても異なります。
なぜスピノールというのか
量子力学をご存知でしょうか。
量子力学では、粒子の状態は複素数で表され、実際に測定できるのはその絶対値のみです。
電子のスピンと呼ばれる状態も考慮したい場合、状態は2成分の複素数に拡張されます。
実際に測定できるのはその絶対値のみであり、かつ確率解釈するためにはそのノルムは1でなければなりません。
そして、絶対値のみを測定できるということは位相回転は物理に影響しない(扱う上で無視したい)ということです。
・・・・なんだか、我々はこれを記述するのに便利そうなものをつい最近みた気がしました。
そう!射影演算子とパウリ行列展開、そしてスピノールです!
詳細はEMANさんのサイトなどを読み進めると参考になるでしょう。
EMANの物理学・物理数学・SU(2)
スピノール(Spin-or)は電子のスピンを記述する言語でもあることから、その名前がついたというわけです。
実数三次元世界で回すと二成分複素ベクトルで表されるスピノールは回ります。(これらは射影演算子を介してつながっています)
実数三次元世界でぐるっと一周してきても、スピノールはまだ半周しかしていない。
つまり最初と異なったスピンの状態(正負反転)になっているのです。
しかし量子力学において、我々は絶対値しか測定できないので正負の符号は物理的には区別できないのです!
(1週回ってきてからスピンを測定しても何事も起きていないように見える)
自然はスピノールの存在を”うまいこと隠してきた”に思えます。
なぜいま偏波か
偏波とパウリ行列が実は思いの外、相性がいいということがわかりました。
偏波というのは実は我々が意識せずに使っています。いま、この瞬間も。
それは、光通信です。
我々はふだん意識せずに通信回線を使っています。
そのような通信回線は、長距離(数十~数百km)を伝送するために必ず光ファイバを経由しています。*5
そして、現在の光ファイバ通信技術は大きく2つの技術で支えられています。
そう!偏波!ここに偏波が使われているのです!
偏波多重信号を生成したり、それを分離したりするためにはJones表現やStokes表現が必須となっています。
現在でもなおこの応用研究は続いており、例えば情報をStokes空間上に埋め込んで送る方法*6なども出てきています。
まとめ
それでは!!
-------------------------------------------------------------------
*1:過去記事参照
*2:これはエネルギーを規格化したことに相当します。更にエネルギー保存系を考える限りは、これはどんな時刻でも満たされていなければなりません
*5:よく「光ファイバってどこで使われているの?」と聞かれます。 光ファイバを経由しない通信は基本的にはないと思って良いです。無線であっても、最終的には集約されて光へ変換され、光ファイバ伝送となります。無線も光も電磁波ですが、この業界では 光 は100THz級の高周波を指す用語です。特に光ファイバでは193THz程度を使います。
*6:実はStokesベクトルは測定が容易という特徴があります。は強度測定器だけで測定可能で、その他の成分については既知かつ実装が容易な座標変換(偏光板)でのように変換したあとでを測れば得られるからです。そこで、安価な通信にとってはJonesベクトルに情報を埋め込むよりも嬉しいのです
*7:このような研究も存在します。その場合、多くの性質はより複雑な等式を満たします。「数学をやっていたら工学もわかっちゃった」、なんて素晴らしいことでしょう!数学やろうぜ!
(番外編)偏波の基礎
こんにちは。Kumaです。
この記事では偏波の基礎を説明します。
この内容はあとあと、これまで書いてきたパウリ行列のお話とリンクしますので
紹介したかったのです。
偏波とはなにか
世の中の波動には縦波と横波が存在します。
- 縦波
媒質の振動方向が波動の進行方向と一致する(平行である)もの。
- 横波、偏波
媒質の振動方向が波動の進行方向と直交するもの。*1
例えば音波は縦波であり、電磁波は横波です。*2
電磁波を考えましょう。横波だということは、波動の進行方向をz軸方向としたときに
媒質の振動はx方向とy方向があるということです。
これをx方向の成分とy方向の成分をもつベクトルで表しましょう。
波動というのは媒質の振動がz方向に”伝わっていく”現象ですから、一般に観測地点と観測時刻の関数です。
いま、観測地点は固定しましょう。そして時間変化は正弦波的であると仮定します。*3
すなわち、
とします。
いま、を変化させたときに
ベクトルの軌跡を考えることができます。これは何らかの曲線だろうと思われます。
この曲線を波動ベクトルの偏波状態といいます。
図でイメージすると、以下のような感じです。
偏光 - Wikipedia
z軸方向から眺めたときの”かたち”と表現する人もいますが、その意味は明らかですね。
偏波の分類
偏波とは、つまるところtで媒介変数表示された
の軌跡でした。軌跡は”振幅”であるの相対関係と”位相”であるの相対関係で決まります。
たとえば簡単な例として位相差がない場合はです。
この場合、は軌跡においてただの初期位相なので無視しましょう。
これは
方向の線分に沿って振動しているだけです。このようなものを直線偏光といいます。線分の方向はの比に依存します。
次にの比が1で逆に位相が90度違うものを考えてみましょう。sinのかわりにcosにすればよいです。
これは円の方程式ですね。
このようなものを円偏光といいます。
時計回りと反時計回りが作れます。
さて、一般の偏光の軌跡はどんなものになりそうでしょうか?
少し考えてみるとわかりますが、一般には楕円を描きます。
desmosでシミュレータを作ってみました。
sとRを変えるといろいろな軌跡ができます。しかしそのいずれも楕円(と円、直線)にしかならないことがわかりますね。
www.desmos.com
話はこれで終わらない。
偏光とはベクトルの軌跡でした。さて、この軌跡はベクトルと【一対一に】対応しているでしょうか。
一見、対応していそうに思えます。
しかし、それは正しくないのです。
直線偏光に対応するようなベクトルは実は2つあるのです。
と
です。これらはベクトルの軌跡の”初期位置”が違うだけで、描く軌跡だけをみると区別できません。
一方で、 の片方だけにマイナスをつけたものは、描く直線が異なるのでこれとは区別できます。
つまり直線偏光は、ベクトルと偏光が【2対1で】対応するのです。
一般の偏光に対してはどうでしょうか?このような”縮退したペア”は存在するでしょうか?
難しくなってきました。
今回はここまでにしましょう!
射影演算子のパウリ行列展開における特殊な性質
こんにちは。Kumaです。
この記事では射影演算子のパウリ行列展開における特殊な性質を紹介します。
ざっくりいうと、パウリ行列で複素数ベクトルを実数のベクトルに埋め込んだとき、
もともとの複素数ベクトルの持っていた性質がどのように”遺伝”するかを与えます。
エルミート行列のパウリ行列展開
2x2エルミート行列は、パウリ行列 と単位行列 の結合で書くことができます。
ここで はすべて実数であって次の式を満たします。
これは以前の記事で証明しています。
electrodynamics.hatenablog.com
射影演算子のパウリ行列展開
天下り的ですが、いま二成分の複素ベクトル
に対して
は射影演算子というのでした。射影演算子はエルミート演算子なので、パウリ行列展開が可能です。
ここで はすべて実数であって次の式を満たします。
複素ベクトルのノルムの遺伝
複素ベクトルのノルム(の二乗)は、パウリ行列展開の世界では次のように”遺伝”しています。
これは の成分表示から明らかです。
すなわち、複素ベクトルが規格化されていることを要請することと を要請することが同じになっています。
さらに、パウリ行列展開の式においてパウリ行列たちのトレースは0だったので、この式は
とも同じ意味です。
つまり、規格化されただけで成り立つ世界を考える場合はがなんであっても常に になるということです。
以降では規格化されただけを考えるとします。*1
すると、自由度はだけということになります。
ここから実数ベクトルを
と定義します。
以降では複素ベクトルと実数ベクトルの関係を調べてみましょう。
複素ベクトルのグローバル位相無視の遺伝
さらに
ですから、を に変えたとしても同じ が得られます。
つまり複素ベクトルから射影演算子 を作ると、のグローバル位相が消えてしまいます。
の自由度はグローバル位相を無視した場合、1つ減っています。射影演算子、ひいては の自由度も、一般の3成分実数ベクトルよりも1つ減っているはずです。
実際、について次の式が成立します。
左辺は複素ベクトルのノルム、右辺はそのパウリ行列展開したときの係数からなる実数ベクトルのノルムであることに注意してください。
いま、は規格化されていて左辺は1ですから、 が成り立ちます。つまり
です。
これが複素ベクトルのグローバル位相無視の遺伝となっています。
射影演算子という特殊なエルミート演算子がもつもう1つの特徴です。
射影演算子のパウリ行列展開
こんにちは。Kumaです。
この記事では射影演算子のパウリ行列展開を紹介します。
ざっくりいうと、複素数ベクトルを実数のベクトルに埋め込む方法になっています。
エルミート行列のパウリ行列展開
2x2エルミート行列は、パウリ行列 と単位行列 の結合で書くことができます。
ここで はすべて実数であって次の式を満たします。
これは以前の記事で証明しています。
electrodynamics.hatenablog.com
複素ベクトル同士の内積と外積
天下り的ですが、いま二成分の複素ベクトル2つを考えます。
さらに、「内積」を定義します。(片方に複素共役を取る以外はふつうのベクトルの内積ですね)
ふたつのベクトルから複素数1つを作り出します。
更に「外積」というものを定義します。
ふたつのベクトルから行列1つを作り出します。
同じベクトル同士の内積をとったとき、それをその複素ベクトルのノルム(の二乗)といいます。
同様に、同じベクトル同士の外積をとったとき、それをその複素ベクトル方向への射影演算子といいます。
(特に、ノルムが1のベクトル同士の場合にいいます)
なぜ同じベクトル同士の外積を射影演算子と呼ぶのでしょうか?確かに行列はベクトルに作用する”演算子”ですが・・・
この理由は読者の練習問題ということで。。
エルミート行列のパウリ行列展開とその座標成分(2)
こんにちは。Kumaです。
この記事はエルミート行列のパウリ行列展開とその座標成分(1)の続きです。
electrodynamics.hatenablog.com
エルミート行列のパウリ行列展開
2x2エルミート行列は、パウリ行列 と単位行列 の結合で書くことができます。
ここで はすべて実数であって次の式を満たします。
本当か?と思うので計算して確かめてみましょう!
証明
------------------
まずtraceの中身 を計算する。
次にtraceを取ります。
よって
となり、確かに行列 に一致します。
-----------------------------------------------
まとめ
という分解公式が成り立ちます。
おまけの性質
この表式
をもう少し掘ってみましょう。
両辺のtraceを考えます。
たちはトレースゼロという性質を持ち、かつtraceは線形演算でtrace(A+B)=trace(A)+trace(B)ですから
が成り立ちます。のトレースはすべて が担っています。つまりはトレースのある部分とない部分に分割されたわけですね。
さらに成分の公式
も少し掘ってみましょう。
ベクトルの場合、ある正規直交な基底ベクトルに関する成分(係数)を計算するには、そのベクトルとの内積を取れば良いのでした。
今回はまるで行列の積をとってtraceを計算することがベクトルの内積に相当するかのようです。実はこれは正しい類推です。
たとえば、たちは次の意味で”正規直交”になっています。
ここで はクロネッカーのデルタと呼ばれ、i=jのときに1, そうでない時には0を返す関数です。
”ベクトルの内積”が"行列積か~ら~の~trace"にすり替わっていると思うと、たしかにこれはベクトルの世界で言う正規直交系に相当します。
それでは!!
エルミート行列のパウリ行列展開とその座標成分(1)
こんにちは。Kumaです。
今回はエルミート行列を展開する方法を紹介します。
これはAdvent Calendarの布石になっています。
エルミート行列の定義
エルミート行列とは、NxNの複素数の正方行列Hであって次の性質を満たすものを指します。
ここで演算子 は複素共役を取ってから転置をすることを意味します。
(簡単のために N = 2 としましょう。)
すなわち、行列の成分を
としたときに、
です。
エルミート行列の実数パラメータ表示
2x2の複素行列は個の複素数(同じことですが 個の実数)を持ちます。
エルミート行列は、その性質 から、以下のように自由度が縛られます。
- 1. から対角要素, は実数でなければならない。
- 2. から、非対角要素, はどちらか一方を決めるともう片方が決まってしまう。
条件1,により, は虚数成分が0と確定するので、を構成する実数の自由度の数は2減ります。
条件2,により,例えばは指定する必要がないので、を構成する実数の自由度の数は2減ります。
以上によりを構成する実数の自由度は となることがわかります。
どうやら、エルミート行列は4個の実数で表示できそうです。
例えば、
特別なエルミート行列としてパウリ行列
もう少しかっこよくを実数で分解する表式があります。そのために
次のような特別なエルミート行列 , , を定義します。
こいつらのことをパウリ行列といいます。
達は何が特別なのかというと、(この記事では使いませんが)変な性質を持っています。
これは添字(1,2,3)をサイクリックに入れ替えても成り立ちます。 は行列なので
かける順番によって結果が異なるわけですが、その”差”は相棒の に定数倍を除いて一致するらしい。
ちょっと便利さがまだよくわからないですね。
次の性質はもっと直感的にも”便利そう”です。
ここで、というのは行列の対角要素の和を取る演算であり、 です。*1
パウリ行列たちは「トレースがゼロな行列」なんですね。
実は行列は行列を特徴づける”固有値”とよばれるN個の大事な数値を持つのですが、固有値の和は
トレースに一致するという定理があります。トレースがわかると、固有値の和はわかるわけです。*2
パウリ行列たちは2個の固有値をもつが、その和は0である ことがわかりました。
まとめ
- となるような2x2複素正方行列をエルミート行列といい、実数4つぶんの自由度がある。
- パウリ行列 , , という特別なものがあり、これらはトレースゼロである。
次に、 をトレースがゼロでない行列とトレースがゼロな行列で分解してみたい・・・のですが
長くなってしまったので、続きは別記事にします。
それでは!!