0
応用数学解説
文献あり

量子コンピュータにおける素因数分解(6):量子ゲートの実装

64
0

はじめに

この記事は量子コンピュータによる素因数分解のアルゴリズムである「Shorのアルゴリズム」の解説記事第6回目です。1-5回目は以下:

  1. 量子コンピュータにおける素因数分解(1): 加算器
  2. 量子コンピュータにおける素因数分解(2): べき剰余
  3. 量子コンピュータにおける素因数分解(3): 量子フーリエ変換
  4. 量子コンピュータにおける素因数分解(4): アルゴリズムの方針 & 関連定理
  5. 量子コンピュータにおける素因数分解(5): 位数推定と量子干渉

本記事では本シリーズの締めくくりとして、今まで紹介した基本的な量子ゲートをいかに物理系で実現するかに関して触れたいと思います。物理学の基礎的な知識が必要な部分がありますし、また議論がいい加減な部分があります(筆者の理解不足によります)。数学徒の方々にはもしかしたら読みにくいかもしれませんがご容赦ください。雰囲気だけ感じてくれれば、と思います。

ここでは2つの実装方法を紹介します。

ひとつは原子のラビ振動を用いた量子ゲートです。これは他の量子ゲートの実装の理解の助けにもなるので紹介します。ただしこのような方法を用いた量子コンピュータが実際に稼働しているかはわかりません。

もうひとつは超伝導ループ・ジョセフソン接合を用いた実現方法ですNielsenHosoyaChoiMukai。これはGoogleのSycamoreプロセッサやIBMの量子コンピュータチップHeronでも採用されている方法であり、現実的な実現方法ですChoiMukai。ただし細かいセットアップは違うと思います。

どちらも発想は同じで、近似的な2準位系を作り、何らかのハミルトニアンを印加することで2準位間の相対位相・重みを変化させ、必要なユニタリー変換を実現します。

シュレーディンガー方程式とBornの規則

本記事ではシュレーディンガー方程式(Schrödinger Equation, SE)を解きます。そのため本章ではその解き方に関して簡単に述べておきます。すでに馴染みがある方は飛ばして「量子ゲートの実装」に進んでください。

シュレーディンガー方程式(SE)

SEは以下の方程式です:

シュレーディンガー方程式(SE)

iψt=H^ψ,   H^:=22m2x2+V(x)

ψは波動関数と呼ばれ、場所xと時刻tに依存します。この波動関数で記述される物体が場所xと時刻tに存在する確率密度は|ψ(x,t)|2で表されます。H^はハミルトニアンと呼ばれ、エネルギーを固有値にもつ演算子です。ハミルトニアンの第1項は運動エネルギー、第2項はポテンシャルエネルギーに対応します。

波動関数は物体に関するすべての情報を担います。時刻t=0における波動関数ψ(x,t=0)が与えられれば、SEを解くことにより任意の時刻での波動関数ψ(x,t)がわかります。すなわちSEは物体の量子力学的状態の時間発展を記述する方程式です。

SEを解くためには、エネルギーの固有値・固有状態を求めればよいです。すなわち以下の方程式
H^ϕn(x)=Enϕn(x)
を解けばよいです。簡単のためここではnは整数とします。これを解くことで{En,ϕn(x)}nZのセットを得ます。Hはエルミート演算子であり、ϕnは完全系をなします。

t=0での波動関数ψ0(x)が与えられたとします。このときψ0を完全系ϕnで展開します:
(1)ψ0(x)=ncnϕn(x)
cnは展開係数であり、一般に複素数です。SEより、個々のエネルギーの固有関数ϕnは以下の時間依存性を持ちます:
exp(iEnt/)ϕn(x)
SEは線形方程式なので、Eq.(1)の時間発展は以下のように与えられます:
ψ(x,t)=ncnexp(iEnt/)ϕn(x)

ということで、量子系を解く際には、方程式H^ϕn(x)=Enϕn(x)を解いてエネルギーの固有値・固有関数を求めることと、波動関数の初期値をこの固有関数で展開することが主な作業となります。

Bornの規則

ある時刻で波動関数がψ=cnϕnだったとします。ただしn|cn|2=1のように規格化されているとします。このとき系のエネルギーを観測すると、Emというエネルギーが確率|cm|2で得られます。Emを観測したのち、波動関数はϕmになります。

これはもっと一般に成立します。ある物理量に対応するエルミート演算子A^の固有状態をϕn(A)、対応する固有値をanとします。系の波動関数がψであるとき、あるamが観測される確率は、波動関数をψ=ncn(A)ϕn(A)のように展開したときの展開係数により|cm(A)|2で与えられます(※n|cn(A)|2=1のように規格化されているとする)。amを観測したのち、波動関数はϕm(A)になります。これをBornの規則と呼びます。

Diracのブラケット表記

以下|0,|1のような記号が現れます。これはDiracのブラケット表記と呼ばれる量子状態の表示方法です。ここでは|0は基底状態(=いちばんエネルギーの低い状態)ϕ0に、|1は第1励起状態(=次にエネルギーの低い状態)ϕ1に対応する記号だと思ってください。すなわち|0,|1はエネルギーの下から2つの状態の基底です。よって
H^|i=Ei|i,   i=0,1
が成立します。

正確に言うと、|0,|1自体はベクトル空間(というかHilbert空間)の元であり、それ以上でも以下でもないです。抽象的な量子状態を表します。これと位置の固有状態|xの双対x|との内積を取った量が波動関数です。すなわち例えばϕ0(x)=x|0です。これは位置表示の波動関数とも呼ばれます。運動量の固有状態の双対p|との内積をとったϕ~(p)=p|0等は運動量表示の波動関数と呼ばれます。でも以下の議論ではそんなに気にしなくてよいです。

行列としてのハミルトニアン

量子論では本来Hilbert空間、すなわち無限次元の空間を扱います。しかしながら量子コンピュータでは基本的には有限次元の空間を扱います。いま相互作用を印加する前のハミルトニアンH^0の固有状態|0,|1があるとします。この系に何らかの相互作用H^Iを印加したとします。全ハミルトニアンをH^=H^0+H^Iとします。このとき
(0|H^|00|H^|11|H^|01|H^|1)
のような行列を作れます。これは|0,|1をそれぞれ(1,0)T,(0,1)Tに対応させてH^を表現したものです。全ハミルトニアンのエネルギーの固有値・固有ベクトルを求めることは、この2行2列の行列の固有値・固有ベクトルを求めることに帰着します。

次章以降に現れる行列としてのハミルトニアンはこのようなものです。

量子ゲートで実現したいこと

本記事の「量子ゲートの実装」とは何を実現することなのかについて述べておきます。

2準位系に関し、以下のような任意の状態の表現法があります:
(2)|ψ=cosθ2|0+eiφsinθ2|1,   0θπ, 0φ<2π
ただしここではoverallの位相ファクターを|0の係数が実になるように固定しています。全体にかかる位相は観測量ではないのでこれで十分です。この状態はθ, φを角度として、半径1の球面上の点で表現できます。これをブロッホ球(Bloch sphere)と呼びます(図1)。

ブロッホ球 ブロッホ球

以下量子ゲートの実装で行うのは、Eq.(2)で表される任意の状態を量子系で実現することです。θ,φを何らかの方法でコントロールできるようにします。これができれば代表的なゲートは以下のように実現できます(量子ゲートに関しては この記事 参照のこと):

  • NOTゲートはθ=0 or π,  φ=0
  • Hadamardゲートはθ=π/2,  φ=0 or π
  • Rkゲートはθ=π/2としてφを変化させる

ということで、以下Eq.(2)の任意の状態を実現する量子系を構成します。

量子ゲートの実装

本章では量子ゲートの実装の具体例を説明します。

実装その1 : 原子によるラビ振動

以下Ref.Hosoyaに沿って説明します。

|0,|1の2準位系を考えます。それは例えば原子の2準位系だとしましょう。本来原子の準位はたくさん存在するのですが、ここでは2準位のみが重要な状況を考えます(それが正当化される状況を作る)。

以下のハミルトニアンを考えます:
H=(0λeiωtλeiωtϵ)
ここで|0,|1を縦ベクトル(1,0)T,(0,1)Tに対応させました。このハミルトニアンは以下のような状況に対応します。基底状態のエネルギーは0、第1励起状態のエネルギーはϵの2つの状態があります。これに角振動数ωの光を時間tだけ照射します。これにより状態は基底状態と励起状態の間を振動します。これはラビ振動と呼ばれます(コヒーレント振動とも呼ばれます)。以下ラビ振動が起こることを具体的に見ていきます。

注入された光によるラビ振動。!FORMULA[80][-221780350][0]のエネルギーをもつ光が系に注入されると、エネルギー差が!FORMULA[81][139816409][0]である2準位系のレベル間で振動が起こる 注入された光によるラビ振動。ωのエネルギーをもつ光が系に注入されると、エネルギー差がϵである2準位系のレベル間で振動が起こる

上記ハミルトニアンによる波動関数の時間発展を求めます。このハミルトニアンは時間依存性を持っているので、まずそれを取り除きます。そこで
ϕ(t):=(a(t)b(t)),   ψ(t):=(a(t)eiωtb(t))
とします。ψ(t)に関するSEは、時間依存性のないハミルトニアンH0を用いて
(3)itϕ(t)=H0ϕ(t),   H0:=(0λλ(ϵω))
のように書き直せます。H0の固有値μ±、および固有ベクトルは

固有値 : μ±=(ϵω)±(ϵω)2+(2λ/)22,固有ベクトル : (λ/μ±1)

です。

証明はここをクリック

H^0の固有値μ|H0μ1|=0
を満たします(||の行列式)。これは
μ=(ϵω)±(ϵω)2+(2λ/)22
で与えられます。ふたつの固有値を符号に対応してμ±とします(μ<μ+, μ<0,μ+>0)。固有関数(=固有ベクトル)は、これを(α,β)Tとすると
(0λλ(ϵω))(αβ)=μ±(αβ)
を満たします。この連立方程式は一次従属なので上成分か下成分かどちらかの方程式のみで固有ベクトルは決定されます。上成分を見れば
α=(λ/μ±)β
なので、固有ベクトルはμ±に対応して
(λ/μ±1)
となります。

ゆえにΦ±:=(λ/μ±,1)Tとすると、SEの一般解Φは、前章の議論より
Φ(t)=A+Φ+eiμ+t/+AΦeiμt/,(4)μ±=(ϵω)±(ϵω)2+(2λ/)22
となります。すなわち角振動数ωの光が系に注入されている間、エネルギー固有値はμ±になり、また波動関数の時間発展はEq.(4)で与えられます。A±は波動関数の初期条件と規格化より定まります。

ここでt=0で波動関数は基底状態|0にあるとします。すなわち
Φ(t=0)=(10)
とします。このときA±はEq.(4)でt=0として
A+Φ++AΦ=(10)
なので
A+=A=μ+μλ(μ+μ)=λ(ϵω)2+(2λ/)2
となります。以上から本セットアップでの波動関数の時間発展は
Φ(t)=λ(ϵω)2+(2λ/)2{(λ/μ+1)eiμ+t/(λ/μ1)eiμt/}
となります。

ここで相互作用がないときのエネルギー差に等しい振動数ω=ϵをもつ光を注入します。このときμ±=±λであり、また
Φ(t)=(cosλtisinλt),ψ(t)=(cosλtieiωtsinλt)
です。いま|0(1,0)T,|1(0,1)Tに対応させているので、これは
(5)ψ(t)cos(λt)|0ieiωtsin(λt)|1
のように書けます。λ,ωをうまく調整できれば、|0|1の相対位相を任意に調整できます。ωは注入する光の角振動数、λは光の振幅で調整します。

Eq.(5)はEq.(2)と本質的に同じです。よってラビ振動を用いれば量子ゲートを作ることができます。ただ筆者はこのような量子ゲートが実際にどこまで作成されたか知りません。量子ゲートではなくて恐縮ですが、現実の原子におけるラビ振動として、例えばセシウム原子のラビ振動を電磁波の強度を測定するために応用する研究が存在しますRabi

その2 : 超伝導量子回路

今度は超伝導体を用いた量子ゲートの実装法を紹介します。最初にも述べたとおり、この方法は実際にGoogleやIBMで使用されている量子コンピュータの実装に用いられていますChoiMukai(※詳細は異なると思います)。この方法が量子ゲートに用いられる理由をRef.Mukaiから引用します:

  1. 制作に半導体加工技術を使うことができる。すなわち現在CPUで使われている技術を応用できる
    ...CPUのような「量子チップ」を作成できる。小さなデバイスで量子コンピュータを実現できる
  2. 回路上の量子系は通信分野で培われてきたマイクロ波技術により精密な制御を実現できる
  3. 超伝導素子のコヒーレンス時間が飛躍的に伸びた
    ...量子コンピュータでは量子系の干渉性が保たれなければならない。これが保たれる時間をコヒーレンス時間と言う。かつてはナノ秒にも満たなかったコヒーレンス時間が現在(※2018年)ではミリ秒に達している

以下超伝導素子を用いた量子ゲートの実装に関して見ていきます。ここではRef.Mukaiを参考にして、実装の概要を記します。 詳細はRefs.KrantzKockumClarke等をご参照ください。

超伝導体、ジョセフソン電流

超伝導体では電流が無抵抗で流れます。超伝導体を流れる電流はクーパー対 (Cooper pair) と呼ばれる2つの電子が相関した状態が担います。電子2つの組なので、クーパー対の電荷は2eです(電子の電荷をeとする)。

無抵抗 & クーパー対に関する注

超伝導体でなぜ無抵抗なのかというと、超伝導体の基底状態と励起状態の間のエネルギー差が大きく(ギャップがある)、その間のエネルギー状態をとれないため、このギャップ内にあるような電子と正イオン(原子核)の散乱状態は許されなくなるからです。
超伝導体でクーパー対が生じる理由は、結晶をなす正イオン(原子核)の格子が電子により歪み、これにより(格子振動により)電子同士に引力が働き、ペアを組んだほうが通常の自由フェルミ粒子模型における基底状態よりエネルギーが低くなることによります。
これらに関しては例えばRef.Niwa等を読んでください。


ここで2つの超伝導体の間に薄い絶縁体を挟みます。これをジョセフソン接合(Josephson junction)と呼びます。

ジョセフソン接合。黄色い部分が超伝導体。黒い部分が絶縁体 ジョセフソン接合。黄色い部分が超伝導体。黒い部分が絶縁体

すると2つの超伝導体の間にジョセフソン電流と呼ばれる電流が流れます。これをジョセフソン効果(Josephson effect)と呼びます。ジョセフソン効果がなぜ起こるかはAppendix 1にゆずります。ここで重要なのは

  1. クーパー対が電荷を担う(電荷2e)
  2. ジョセフソン電流が流れる

の2つです。これは認めてください。

超伝導量子回路の実装

以下Ref.Mukaiに沿って説明します。図3のような超電導量子回路を考えます。

SQUIDを用いた電荷型の超伝導量子回路 SQUIDを用いた電荷型の超伝導量子回路

図の黄色の部分は超伝導体です。左右に1つづつジョセフソン接合があります。ジョセフソン接合2つを超伝導体でつないでループにしたデバイスをSQUID(Superconducting Quantum Interference Device)と呼びます。boxと書いてある部分にはクーパー対が溜まります。この箱に貯まるクーパー対の個数の固有状態が本超伝導量子回路における2準位系をなします。ΦはSQUID回路の中を通る磁場を表します。後で述べますが、Φにより2準位間の遷移の周期を調整します。黒い2本線(||の部分)はキャパシタ(コンデンサー)です。Vg,Vb,Vpはそれぞれゲート、ベース、プローブにかける電位です。

boxの部分のエネルギー、すなわちそれに対するハミルトニアンHは「ジョセフソン電流によるエネルギー+箱に蓄えられたクーパー対の静電エネルギー」で与えられます。電荷による静電エネルギーはコンデンサーと同様で、電荷をQ、静電容量をCとするとQ2/2Cで与えられます。さらにAppendix 1より、ジョセフソン電流によるエネルギーはEJcosθで与えられます。θはジョセフソン接合における超伝導体の巨視的波動関数の位相差です。以上からboxにおけるハミルトニアンは
(6)H=Ec(nng)2EJcosθ
で与えられます。nは箱の電荷、ngはゲートに溜まる電荷です。ここでEcは熱エネルギー換算で数百mKです。Mukai。よって数十mK程度の低温下では、熱浴によるデコヒーレンスを避けることができます。

系が量子力学で記述される場合、クーパー対がnコ存在する量子状態|nで箱の状態を記述できます。ここでは|0がクーパー対0コの状態、|1が1コの状態だと思ってください。|1より上の状態は無視できて、|0,|1のみで記述できる環境が実現できます。このときEq.(3)の第1項は|0,|1の固有状態が独立に存在する場合のエネルギー、第2項はその間の遷移を起こす相互作用とみなせます(電流が流れることでクーパー対の個数を変化させる)。そのようなハミルトニアンは

H=Ec2(ng1)σz+EJ2σx=(Ec2(ng1)EJ/2EJ/2Ec2(ng1))

のように書くことができます。これは|0(0,1)T, |1(1,0)Tと対応させ、また|0,|1のエネルギーの基準を適切にずらせば(Ec(ng1)/2をゼロとする)、Eq.(2)のハミルトニアンH0と同じ形をしています。ゆえにジョセフソン電流が流れると|0,|1の間でラビ振動を起こし、それらの相対位相を任意にコントロールすることができます。

Eq.(2)のλにあたるEJを調整することで|0,|1の間を振動する周期を調整できます。これはSQUIDに通す外部磁束によってコントロールできます。SQUIDの超伝導体中のクーパー対を記述する波動関数の位相は、SQUIDを一周したとき、波動関数が安定的に存在するためにmod 2πで元に戻らなければなりません。またAharonov-Bohm効果より、一周したときの波動関数の位相差は、SQUID内に外場として加えた磁場Φextに依存します。ジョセフソン接合部分での波動関数の位相差θは、波動関数がSQUIDを一周してmod 2πで元に戻る条件より、Φextに依存して変化します。これはEJを変化させるのと同様の効果を生みます。具体的にはEJcos(πΦext/Φ0)のような依存性を持ちますMukaiKrantzKockum。SQUIDおよびAharonov-Bohm効果に関してはAppendix 2をご参照ください。

このようにして超伝導量子回路で量子ゲートを作成することができます。

シリーズ全体のまとめ

今回が「量子コンピュータにおける素因数分解」の最終回なので、全体の総括をしておきます。

量子コンピュータで最も有名であろうアルゴリズム「Shorのアルゴリズム」に関して説明してきました。このアルゴリズムで重要な量子回路は、べき剰余を計算する回路:
ax mod N (a,xはゼロ以上の整数,Nは正の整数. 0a<N)
および量子フーリエ変換:
|j12n/2k=02n1exp(2πijk/2n)|k,   nN, 
を行う回路です。これら回路を用いて
xr1 mod N
を満たすrを見つけます。このrが偶数かつ「適切」であれば、gcd(xr/2±1,N)を古典コンピュータで計算することによりNの因数を見つけることができます。量子回路による計算量、古典コンピュータによる計算量、および試行錯誤が必要な確率的要素を加味しても、log2Nの多項式程度の計算量で因数分解を行うことができます。現在まで、古典コンピュータではこのような高速なアルゴリズムは発見されていません。

上記量子回路で用いられている基本的な量子ゲートの実装方法に関しては、本記事で説明したように、いまのところ超伝導量子回路を用いた方法が主流です。クーパー対が0または1コの量子状態間の遷移をSQUIDで調整することにより、任意の1 qubitの状態を実現することができます。

現在のところ、Shorのアルゴリズムによる素因数分解は2桁程度の整数に留まっているようです。現在最も安全な暗号化技術のひとつである、素因数分解の擬似非可逆性を用いたRSA暗号を破るのに有用な桁数の素因数分解はまだまだ先のようです。しかしながら、Shorのアルゴリズムはその古典計算機に対する本質的な優位性、暗号解読という非常に重要な分野への応用、また様々な量子アルゴリズム開発の指針となるような発想を与えたという意味で、理論的なレベルでも十分価値のあるものだと思います。

おしまい。



Appendix 1 : ジョセフソン効果

本章では超伝導量子回路において基本的な現象であるジョセフソン効果について述べます。特にジョセフソン電流がどうして流れるか、またジョセフソン電流によるエネルギーであるジョセフソンエネルギーについて述べます。参考文献としてRefs.NiwaWikiJosTsunetoを挙げておきます。

ジョセフソン電流

以下Ref.Niwaに沿って説明します。

超伝導体は1つの巨大な(人間のスケールからすると小さいですが)波動関数ψ(r,t)で記述できます。ただし2つの超伝導体ではその位相が異なるとしましょう。いま空間1次元で考えます。2つの超伝導体を1,2とラベルし、それぞれの波動関数をψ1(r,t),ψ2(r,t)として、さらにこれらを極座標で表示します:
ψ1(r,t)=|ψ1(r,t)|exp(iϕ1(r,t)),ψ2(r,t)=|ψ2(r,t)|exp(iϕ2(r,t))
これらの波動関数はクーパー対に対する波動関数です。

この系のシュレーディンガー方程式がどうなるか考えます。いま2つの超伝導体は薄い絶縁体により弱く結合しています。全系のハミルトニアンは
ψ(r,t)=c1ψ1(r,t)+c2ψ2(r,t)
のように書かれます。全形のハミルトニアンは
H^=H^1+H^2+H^T
のように書かれます。ここでH^1,H^2はそれぞれ超伝導体1,2が孤立している場合のハミルトニアン、H^Tはそれらが弱く結合させるハミルトニアンです。いま両方の超伝導体の間に時間に依存しない電圧Vを印加したとします。H^1,H^2のそれぞれのエネルギー固有値をE1,E2として、接合の両端におけるE1,E2の差は
E1E2=2eV
となります。2eはクーパー対の電荷です。いまエネルギーの原点を|E1E2|/2のようにエネルギー差の半分のところにとれば、エネルギーの固有値方程式:H^ψ=Eψはそれぞれの超伝導体で
H^1ψ1(r,t)=eVψ1(r,t),   H^2ψ2(r,t)=eVψ2(r,t)
となります。さらにH^Tに対する仮定として、l=1,2に対するシュレーディンガー方程式に対し、H^Tψl(r,t)Kψl¯(r,t)に置き換えます。ただしl=1,2に対しそれぞれl¯=2,1とします。すなわち1,2の超伝導体に関する波動関数に対するSEを作るときには、物理系のセットアップから定まるあるKという定数を「相手の波動関数」にかけた相互作用ハミルトニアンが加わるとします。すると
iψ˙1(r,t)=eVψ1(r,t)+Kψ2(r,t)iψ˙2(r,t)=eVψ2(r,t)+Kψ1(r,t)
が成立します。ここでドットは時間による偏微分です。

このSEに極座標に分解したψ1,ψ2を入れて解きます。位相差をΔϕ(r,t):=ϕ2(r,t)ϕ1(r,t)とします。上記SEを実部と虚部に分けたうち、虚部のほうをみると
2e|ψ1(r,t)|2t=2eK|ψ1(r,t)||ψ2(r,t)|sin(Δϕ(r,t)),2e|ψ2(r,t)|2t=2eK|ψ1(r,t)||ψ2(r,t)|sin(Δϕ(r,t))
を得ます。ここでも2eはクーパー対の電荷です。クーパー対がこの2つの超伝導体の間を移動すると電流が流れますが、このときの電流は、面積Sの接合面を流れるとすれば
I=2eS|ψ1(r,t)|2t=2eS|ψ2(r,t)|2t(A1)=2eKS|ψ1(r,t)||ψ2(r,t)|sin(Δϕ(r,t))
となります。

ここで注目すべきは、この方程式はVには依存していないことです。すなわち2つの超伝導体に位相差Δϕさえ存在していれば電流が流れます。もし超伝導体1,2の初期の位相差が時間・空間に関し一定ならば、Iは直流電流です。これは直流ジョセフソン効果 (DC Josephson effect) と呼ばれます。ふつう電流と電圧の関係は抵抗をRとしてオームの法則V=IRで記述されますが、ジョセフソン電流に電位差は必要ないし抵抗も関係ありません。

次に実部を解析します。実部の方程式より
Δϕ(r,t)t=2eV+(|ψ2(r,t)|2|ψ1(r,t)|2|ψ1(r,t)||ψ2(r,t)|)Kcos(Δϕ(r,t))
を得ます。1,2の超伝導体は同素材で、位相だけが異なるとすれば|ψ1(r,t)|2=|ψ2(r,t)|2が成立します。このとき上の式の右辺第2項は消えて
(A2)Δϕ(r,t)t=2eV
と非常に簡単になります。これより電圧Vが印加された場合、Δϕ(r,t)は時間tに関して線形の依存性をもつことがわかります:
Δϕ(r,t)=Δϕ(r,t=0)2eVt
これとEq.(1)より、電圧が印加された場合は周期ω=2eV/の交流電流が流れることがわかります。これを交流ジョセフソン効果 (AC Josephson effect) とよび、またこの周期をジョセフソン周波数 (Josephson frequency) と呼びます。

ジョセフソン効果は、超伝導の微視的理論であるBCS理論からも、また対称性の自発的破れが生じた際に現れる秩序変数に関する一般論であるGinzburg-Landau理論からも導出できますNiwa

ジョセフソンエネルギー

ここではジョセフソン電流によるエネルギーを計算します。これはジョセフソン接合における位相差が0からΔϕまで変化した際に流れるジョセフソン電流のエネルギーです。Ref.WikiJosを参考にしています。

位相差をΦとして、Eq.(A1)を
I=IcsinΦ
と書いておきます。またEq.(A2)より
V=2eΦt
です。Φが時刻t=00tΦ=Δϕとします。この間に流れたジョセフソン電流によるエネルギーは
E=0tIVdt=0ΔϕIcsinΦ d(2eΦ)=2eIc[cosΦ]0Δϕ
以上より
E(Δϕ)=EJ(1cosΔϕ),   EJ=2eIc
を得ます。

Appendix 2 : SQUIDとAharonov-Bohm (AB) 効果

SQUIDとは超伝導量子干渉計(Superconducting Quantum Interference Device)の略です。本文図3にあるように、超伝導ループの中に2つのジョセフソン接合を持つデバイスです。通常極めて弱い磁場の検出に用いられる磁気センサーとして用いられますがWiki、本文で述べたように、量子コンピュータでは基底状態と励起状態の振動周期の調整に用いられます。

ここではなぜSQUIDに通す磁場によりその中を流れる電流を制御できるかについて説明します。これはAharonov-Bohm効果と関係します。

Aharonov-Bohm効果

以下Ref.Igiに沿って議論します。Ref.Wachterも参考文献として挙げておきます。

電荷eの粒子の波動関数をψとします。電場・磁場があるとき、SEは以下のようになります:
(A3)iψt={22m(xiecA)2+eϕ}ψ

ここでA,ϕはそれぞれベクトルポテンシャル、スカラーポテンシャルと呼ばれます。これらの量を組にしてAμ=(ϕ,A)のように4次元のベクトルとしたものはゲージ場と呼ばれます。ゲージ場と電場Eと磁場Bとの関係は以下です:

(A4)E=ϕ+1cAt,   B=rotA

Eq.(A3)はχ(x,t)を場所・時刻に依存する任意の関数として

AAχ,   ϕϕ+1ctχ,   ψexp[iecχ]ψ

に対して不変です。この変換をゲージ変換、この変換に対して不変なことをゲージ不変性、χをゲージ関数と呼びます。観測可能な量はすべてゲージ不変です。電場・磁場は観測可能なのでゲージ不変です。「物理量はゲージ不変であれ」という原理はゲージ原理と呼ばれます。ゲージ原理を指針として、SEのゲージ不変性を要請することでEq.(A3)を導くこともできます。A,ϕを結合させ、かつゲージ不変にするには、時間および空間微分を

ixixecA,   ititeϕ
のように置き換えればよいことが知られています。このようにゲージ場と粒子との結合を得る方法はminimal couplingなどと呼ばれます。

ここで重要なのは、量子論においては電場・磁場ではなくゲージ場が本質的に現れることです。局所的な場の理論では電場・磁場のみで運動方程式を書くことは不可能です。

量子論におけるゲージ場の重要性が如実に現れるのがAharanov-Bohm効果です。いま電場は無視してϕ=0とします。Eq.(A3)でϕ=0とした方程式

(A5)iψt=22m(xiecA)2ψ

となりますが、ここでψ

(A6)ψ=exp[iecχ]ψ

とします。ここでχA=χを満たす関数とします。このときB=rotA=×Aはゼロです。これは直接計算でももちろんわかりますが、次のような考察からもわかります:A=χA=0からゲージ関数χによるゲージ変換をしたベクトルポテンシャルですが、A=0で磁場がゼロなのは明らかなので、A=χによる磁場もそのゲージ不変性からゼロになります。

ψでEq.(A5)を書き直すと

itψ=22m2ψ

となります。すなわちAがノンゼロでも、それが与えるBがゼロならば、ψは自由なSEに従います。波動関数のoverallな位相は物理的ではないので、結局この状況では粒子は自由運動をすることがわかります。

さて、A=χを満たすχは線積分を用いて
(A7)χ(x)=C x0xA(x)dx    (Cは粒子の経路)
と書けることは、これをA=χに代入すればすぐにわかります。Eq.(A7)とEq.(A6)と合わせると、ベクトルポテンシャルの存在の元で運動する粒子には、その粒子の軌道上に存在するベクトルポテンシャルA(x)に依存した位相がかかることがわかります。

上記したように位相そのものは観測量ではありません。しかし位相差は観測にかかる場合があります。ゲージ場が波動関数に及ぼす影響を見るために、以下の図のような2重スリットの状況を考えます。

磁場の存在する2重スリット実験。荷電粒子(電荷!FORMULA[264][37763][0])が!FORMULA[265][891513128][0]から!FORMULA[266][1565171705][0]に到達する様子。!FORMULA[267][477447113][0]はそれぞれスリットA,Bを通る場合の到達経路 磁場の存在する2重スリット実験。荷電粒子(電荷e)がx0からxに到達する様子。CA,CBはそれぞれスリットA,Bを通る場合の到達経路

ゲージ場が存在しない場合、経路A,Bを通りx0からxまで移動した粒子のxにおける波動関数をそれぞれψfreeA,B(x)とします。一方ゲージ場が存在する場合の同様の波動関数をそれぞれψgA,B(x)とします。するとEq.(A6)(A7)より

ψgA(x)=exp[iecCA x0xA(x)dx]ψfreeA,   ψgB(x)=exp[iecCB x0xA(x)dx]ψfreeB

となります。よってゲージ場が存在する下で、xにおける波動関数は、経路A,Bを通った波動関数の和になるので

ψtot(x)=ψgA(x)+ψgB(x)=exp[iecCA x0xA(x)dx]{ψfreeA(x)+exp[iec(CB x0xA(x)dxCA x0xA(x)dx)]ψfreeB(x)}=exp[iecCA x0xA(x)dx]{ψfreeA(x)+exp[iecCA(x)dx]ψfreeB(x)}

となります。ここでC=CBCAです。ストークスの定理を用いれば

CA(x)dx=SrotAdS(A8)=SBdS=:Φ

となります。ここでSCを境界とする2次元面です。ΦSを貫く磁束です。よって、粒子がx0からxに到達する確率は、その経路を境界とする2次元面を通る磁場により変化します。Eq.(A8)はゲージ変換AAχで不変であることに注意してください。

粒子がCAまたはCBを通りx0からxに至る確率密度は|ψtot|2です。軌道Cの中に磁場がなければ

|ψtot|2=|ψfreeA+ψfreeB|2

です。一方C内に磁場が存在すれば

|ψtot|2=|ψfreeA(x)+exp[iecCA(x)dx]ψfreeB(x)|2=|ψfreeA|2+|ψfreeB|2+(ψfreeA)exp[ieΦc]ψfreeB(x)+ψfreeA(exp[ieΦc]ψfreeB(x))

となります。最後の2項がSを貫く磁束Φに依存します。よってxに粒子が至る確率もΦに依存します。これがAharanov-Bohm(AB)効果です。上記したようにAB効果はゲージ不変な現象です。

超伝導ループとSQUID

ここまでは電子を2重スリットに通す場合を考えましたが、同様なことを超伝導体のループにあてはめます。超伝導体は1つの巨視的な(クーパー対を記述する)波動関数で表せます。波動関数が安定して存在するには、それが超伝導体を一周すると元に戻らなければなりません。ところが磁場をループ内に通すと、波動関数が一周する間にAB効果により波動関数に2eΦ/cという余計な位相がつきます。Bohr-Sommerfeldの量子化条件と同様、超伝導体内を一周した波動関数が元に戻る条件:

2eΦc=2nπ, nZ

を満たす必要があります。Sは超伝導体ループを境界とする2次元面です。このように超伝導体ループではループを貫く磁場は離散化された値を取ります。いわゆる磁束の量子化が起こります。

SQUIDでは超伝導体ループの途中に存在するジョセフソン接合にさらに位相差θが存在するので
θ+2eΦc=2nπ,   nZ

が成立します。故に超伝導体の位相差θを外部からSQUIDにかける磁束Φにより調節できます。これによりジョセフソン電流はΦの依存性をもち、ひいてはEq.(8)のEJは(有効的に)Φにより変化します。

余談ですが、Aharanov-Bohm効果は、その存在を多くの物理学者が認めるまでに理論的提案から30年ほどかかりました。その間に幾多の反論を受けました。AB効果の存在を示す決定的な実験が日立の基礎研において行われたことは有名です。それに至るまでの経緯およびその実験に関しては外村彰さんの書かれたRef.Tonomuraがわかりやすいです。なぜAB効果に関してそれほど批判があり認め難いのかに関してもこの論文を読まれるとよいと思います。

参考文献

[1]
Nielsen, M.A., Chuang, I.L., Quantum Computation and Quantum Information (10th Anniversary Edition), Cambridge University Press, 2010
[2]
細谷 暁夫, 量子コンピュータの基礎 [第2版], SGCライブラリ69, サイエンス社, 2009, 39-41
[4]
向井 寛人, 友永 顕成, 蔡 兆申, 超伝導量子コンピュータの基礎と最先端, 低温工学 (J. Cryo Super. Soc. Jpn), 2004
[6]
Krantz, P. et al., A quantum engineer's guide to superconducting qubits, Appl. Phys. Rev. 6, 2019, 6
[7]
Kockum, A. F., Nori, F., Quantum bits with Josephson junctions, Fundamentals and Frontiers of the Josephson Effect (Springer International Publishing), 2019, 703-741
[8]
Clarke, John, Wilhelm, Frank K., Superconducting quantum bits, Nature, 2008, 1031-1042
[9]
丹羽,雅昭, 超電導の基礎 [第3版], 東京電機大学出版, 2008, 403-429
[11]
恒藤 敏彦, 超伝導・超流動, 岩波講座 現代の物理学17, 岩波書店, 1993, 55-63
[13]
猪木 慶治, 川合 光, 量子力学 II, 講談社サイエンティフィック, 1994, 305-311
[15]
外村 彰, アハラノフ-ボーム効果とその検証実験 , 日本物理学会誌, 1987, 616–624
投稿日:2024920
更新日:2024921
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

bisaitama
bisaitama
143
67387

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. シュレーディンガー方程式とBornの規則
  3. シュレーディンガー方程式(SE)
  4. Bornの規則
  5. Diracのブラケット表記
  6. 行列としてのハミルトニアン
  7. 量子ゲートで実現したいこと
  8. 量子ゲートの実装
  9. 実装その1 : 原子によるラビ振動
  10. その2 : 超伝導量子回路
  11. シリーズ全体のまとめ
  12. Appendix 1 : ジョセフソン効果
  13. ジョセフソン電流
  14. ジョセフソンエネルギー
  15. Appendix 2 : SQUIDとAharonov-Bohm (AB) 効果
  16. Aharonov-Bohm効果
  17. 超伝導ループとSQUID
  18. 参考文献