はじめに
この記事は量子コンピュータによる素因数分解のアルゴリズムである「Shorのアルゴリズム」の解説記事第6回目です。1-5回目は以下:
-
量子コンピュータにおける素因数分解(1): 加算器
-
量子コンピュータにおける素因数分解(2): べき剰余
-
量子コンピュータにおける素因数分解(3): 量子フーリエ変換
-
量子コンピュータにおける素因数分解(4): アルゴリズムの方針 & 関連定理
-
量子コンピュータにおける素因数分解(5): 位数推定と量子干渉
本記事では本シリーズの締めくくりとして、今まで紹介した基本的な量子ゲートをいかに物理系で実現するかに関して触れたいと思います。物理学の基礎的な知識が必要な部分がありますし、また議論がいい加減な部分があります(筆者の理解不足によります)。数学徒の方々にはもしかしたら読みにくいかもしれませんがご容赦ください。雰囲気だけ感じてくれれば、と思います。
ここでは2つの実装方法を紹介します。
ひとつは原子のラビ振動を用いた量子ゲートです。これは他の量子ゲートの実装の理解の助けにもなるので紹介します。ただしこのような方法を用いた量子コンピュータが実際に稼働しているかはわかりません。
もうひとつは超伝導ループ・ジョセフソン接合を用いた実現方法ですNielsenHosoyaChoiMukai。これはGoogleのSycamoreプロセッサやIBMの量子コンピュータチップHeronでも採用されている方法であり、現実的な実現方法ですChoiMukai。ただし細かいセットアップは違うと思います。
どちらも発想は同じで、近似的な2準位系を作り、何らかのハミルトニアンを印加することで2準位間の相対位相・重みを変化させ、必要なユニタリー変換を実現します。
シュレーディンガー方程式とBornの規則
本記事ではシュレーディンガー方程式(Schrödinger Equation, SE)を解きます。そのため本章ではその解き方に関して簡単に述べておきます。すでに馴染みがある方は飛ばして「量子ゲートの実装」に進んでください。
シュレーディンガー方程式(SE)
SEは以下の方程式です:
は波動関数と呼ばれ、場所と時刻に依存します。この波動関数で記述される物体が場所と時刻に存在する確率密度はで表されます。はハミルトニアンと呼ばれ、エネルギーを固有値にもつ演算子です。ハミルトニアンの第1項は運動エネルギー、第2項はポテンシャルエネルギーに対応します。
波動関数は物体に関するすべての情報を担います。時刻における波動関数が与えられれば、SEを解くことにより任意の時刻での波動関数がわかります。すなわちSEは物体の量子力学的状態の時間発展を記述する方程式です。
SEを解くためには、エネルギーの固有値・固有状態を求めればよいです。すなわち以下の方程式
を解けばよいです。簡単のためここではは整数とします。これを解くことでのセットを得ます。はエルミート演算子であり、は完全系をなします。
での波動関数が与えられたとします。このときを完全系で展開します:
は展開係数であり、一般に複素数です。SEより、個々のエネルギーの固有関数は以下の時間依存性を持ちます:
SEは線形方程式なので、Eq.(1)の時間発展は以下のように与えられます:
ということで、量子系を解く際には、方程式を解いてエネルギーの固有値・固有関数を求めることと、波動関数の初期値をこの固有関数で展開することが主な作業となります。
Bornの規則
ある時刻で波動関数がだったとします。ただしのように規格化されているとします。このとき系のエネルギーを観測すると、というエネルギーが確率で得られます。を観測したのち、波動関数はになります。
これはもっと一般に成立します。ある物理量に対応するエルミート演算子の固有状態を、対応する固有値をとします。系の波動関数がであるとき、あるが観測される確率は、波動関数をのように展開したときの展開係数によりで与えられます(※のように規格化されているとする)。を観測したのち、波動関数はになります。これをBornの規則と呼びます。
Diracのブラケット表記
以下のような記号が現れます。これはDiracのブラケット表記と呼ばれる量子状態の表示方法です。ここではは基底状態(=いちばんエネルギーの低い状態)に、は第1励起状態(=次にエネルギーの低い状態)に対応する記号だと思ってください。すなわちはエネルギーの下から2つの状態の基底です。よって
が成立します。
注
正確に言うと、自体はベクトル空間(というかHilbert空間)の元であり、それ以上でも以下でもないです。抽象的な量子状態を表します。これと位置の固有状態の双対との内積を取った量が波動関数です。すなわち例えばです。これは位置表示の波動関数とも呼ばれます。運動量の固有状態の双対との内積をとった等は運動量表示の波動関数と呼ばれます。でも以下の議論ではそんなに気にしなくてよいです。
行列としてのハミルトニアン
量子論では本来Hilbert空間、すなわち無限次元の空間を扱います。しかしながら量子コンピュータでは基本的には有限次元の空間を扱います。いま相互作用を印加する前のハミルトニアンの固有状態があるとします。この系に何らかの相互作用を印加したとします。全ハミルトニアンをとします。このとき
のような行列を作れます。これはをそれぞれに対応させてを表現したものです。全ハミルトニアンのエネルギーの固有値・固有ベクトルを求めることは、この2行2列の行列の固有値・固有ベクトルを求めることに帰着します。
次章以降に現れる行列としてのハミルトニアンはこのようなものです。
量子ゲートで実現したいこと
本記事の「量子ゲートの実装」とは何を実現することなのかについて述べておきます。
2準位系に関し、以下のような任意の状態の表現法があります:
ただしここではoverallの位相ファクターをの係数が実になるように固定しています。全体にかかる位相は観測量ではないのでこれで十分です。この状態はを角度として、半径1の球面上の点で表現できます。これをブロッホ球(Bloch sphere)と呼びます(図1)。
ブロッホ球
以下量子ゲートの実装で行うのは、Eq.(2)で表される任意の状態を量子系で実現することです。を何らかの方法でコントロールできるようにします。これができれば代表的なゲートは以下のように実現できます(量子ゲートに関しては
この記事
参照のこと):
- NOTゲートは
- Hadamardゲートは
- ゲートはとしてを変化させる
ということで、以下Eq.(2)の任意の状態を実現する量子系を構成します。
量子ゲートの実装
本章では量子ゲートの実装の具体例を説明します。
実装その1 : 原子によるラビ振動
以下Ref.Hosoyaに沿って説明します。
の2準位系を考えます。それは例えば原子の2準位系だとしましょう。本来原子の準位はたくさん存在するのですが、ここでは2準位のみが重要な状況を考えます(それが正当化される状況を作る)。
以下のハミルトニアンを考えます:
ここでを縦ベクトルに対応させました。このハミルトニアンは以下のような状況に対応します。基底状態のエネルギーは、第1励起状態のエネルギーはの2つの状態があります。これに角振動数の光を時間だけ照射します。これにより状態は基底状態と励起状態の間を振動します。これはラビ振動と呼ばれます(コヒーレント振動とも呼ばれます)。以下ラビ振動が起こることを具体的に見ていきます。
注入された光によるラビ振動。のエネルギーをもつ光が系に注入されると、エネルギー差がである2準位系のレベル間で振動が起こる
上記ハミルトニアンによる波動関数の時間発展を求めます。このハミルトニアンは時間依存性を持っているので、まずそれを取り除きます。そこで
とします。に関するSEは、時間依存性のないハミルトニアンを用いて
のように書き直せます。の固有値、および固有ベクトルは
です。
証明はここをクリック
の固有値は
を満たします(はの行列式)。これは
で与えられます。ふたつの固有値を符号に対応してとします(, )。固有関数(=固有ベクトル)は、これをとすると
を満たします。この連立方程式は一次従属なので上成分か下成分かどちらかの方程式のみで固有ベクトルは決定されます。上成分を見れば
なので、固有ベクトルはに対応して
となります。
ゆえにとすると、SEの一般解は、前章の議論より
となります。すなわち角振動数の光が系に注入されている間、エネルギー固有値はになり、また波動関数の時間発展はEq.(4)で与えられます。は波動関数の初期条件と規格化より定まります。
ここでで波動関数は基底状態にあるとします。すなわち
とします。このときはEq.(4)でとして
なので
となります。以上から本セットアップでの波動関数の時間発展は
となります。
ここで相互作用がないときのエネルギー差に等しい振動数をもつ光を注入します。このときであり、また
です。いまに対応させているので、これは
のように書けます。をうまく調整できれば、との相対位相を任意に調整できます。は注入する光の角振動数、は光の振幅で調整します。
Eq.(5)はEq.(2)と本質的に同じです。よってラビ振動を用いれば量子ゲートを作ることができます。ただ筆者はこのような量子ゲートが実際にどこまで作成されたか知りません。量子ゲートではなくて恐縮ですが、現実の原子におけるラビ振動として、例えばセシウム原子のラビ振動を電磁波の強度を測定するために応用する研究が存在しますRabi。
その2 : 超伝導量子回路
今度は超伝導体を用いた量子ゲートの実装法を紹介します。最初にも述べたとおり、この方法は実際にGoogleやIBMで使用されている量子コンピュータの実装に用いられていますChoiMukai(※詳細は異なると思います)。この方法が量子ゲートに用いられる理由をRef.Mukaiから引用します:
- 制作に半導体加工技術を使うことができる。すなわち現在CPUで使われている技術を応用できる
...CPUのような「量子チップ」を作成できる。小さなデバイスで量子コンピュータを実現できる - 回路上の量子系は通信分野で培われてきたマイクロ波技術により精密な制御を実現できる
- 超伝導素子のコヒーレンス時間が飛躍的に伸びた
...量子コンピュータでは量子系の干渉性が保たれなければならない。これが保たれる時間をコヒーレンス時間と言う。かつてはナノ秒にも満たなかったコヒーレンス時間が現在(※2018年)ではミリ秒に達している
以下超伝導素子を用いた量子ゲートの実装に関して見ていきます。ここではRef.Mukaiを参考にして、実装の概要を記します。 詳細はRefs.KrantzKockumClarke等をご参照ください。
超伝導体、ジョセフソン電流
超伝導体では電流が無抵抗で流れます。超伝導体を流れる電流はクーパー対 (Cooper pair) と呼ばれる2つの電子が相関した状態が担います。電子2つの組なので、クーパー対の電荷は2eです(電子の電荷をeとする)。
無抵抗 & クーパー対に関する注
超伝導体でなぜ無抵抗なのかというと、超伝導体の基底状態と励起状態の間のエネルギー差が大きく(ギャップがある)、その間のエネルギー状態をとれないため、このギャップ内にあるような電子と正イオン(原子核)の散乱状態は許されなくなるからです。
超伝導体でクーパー対が生じる理由は、結晶をなす正イオン(原子核)の格子が電子により歪み、これにより(格子振動により)電子同士に引力が働き、ペアを組んだほうが通常の自由フェルミ粒子模型における基底状態よりエネルギーが低くなることによります。
これらに関しては例えばRef.Niwa等を読んでください。
ここで2つの超伝導体の間に薄い絶縁体を挟みます。これをジョセフソン接合(Josephson junction)と呼びます。
ジョセフソン接合。黄色い部分が超伝導体。黒い部分が絶縁体
すると2つの超伝導体の間にジョセフソン電流と呼ばれる電流が流れます。これをジョセフソン効果(Josephson effect)と呼びます。ジョセフソン効果がなぜ起こるかはAppendix 1にゆずります。ここで重要なのは
- クーパー対が電荷を担う(電荷2e)
- ジョセフソン電流が流れる
の2つです。これは認めてください。
超伝導量子回路の実装
以下Ref.Mukaiに沿って説明します。図3のような超電導量子回路を考えます。
SQUIDを用いた電荷型の超伝導量子回路
図の黄色の部分は超伝導体です。左右に1つづつジョセフソン接合があります。ジョセフソン接合2つを超伝導体でつないでループにしたデバイスをSQUID(Superconducting Quantum Interference Device)と呼びます。boxと書いてある部分にはクーパー対が溜まります。この箱に貯まるクーパー対の個数の固有状態が本超伝導量子回路における2準位系をなします。はSQUID回路の中を通る磁場を表します。後で述べますが、により2準位間の遷移の周期を調整します。黒い2本線(||の部分)はキャパシタ(コンデンサー)です。はそれぞれゲート、ベース、プローブにかける電位です。
boxの部分のエネルギー、すなわちそれに対するハミルトニアンは「ジョセフソン電流によるエネルギー+箱に蓄えられたクーパー対の静電エネルギー」で与えられます。電荷による静電エネルギーはコンデンサーと同様で、電荷を、静電容量をとするとで与えられます。さらにAppendix 1より、ジョセフソン電流によるエネルギーはで与えられます。はジョセフソン接合における超伝導体の巨視的波動関数の位相差です。以上からboxにおけるハミルトニアンは
で与えられます。は箱の電荷、はゲートに溜まる電荷です。ここでは熱エネルギー換算で数百mKです。Mukai。よって数十mK程度の低温下では、熱浴によるデコヒーレンスを避けることができます。
系が量子力学で記述される場合、クーパー対がコ存在する量子状態で箱の状態を記述できます。ここではがクーパー対0コの状態、が1コの状態だと思ってください。より上の状態は無視できて、のみで記述できる環境が実現できます。このときEq.(3)の第1項はの固有状態が独立に存在する場合のエネルギー、第2項はその間の遷移を起こす相互作用とみなせます(電流が流れることでクーパー対の個数を変化させる)。そのようなハミルトニアンは
のように書くことができます。これはと対応させ、またのエネルギーの基準を適切にずらせば(をゼロとする)、Eq.(2)のハミルトニアンと同じ形をしています。ゆえにジョセフソン電流が流れるとの間でラビ振動を起こし、それらの相対位相を任意にコントロールすることができます。
Eq.(2)のにあたるを調整することでの間を振動する周期を調整できます。これはSQUIDに通す外部磁束によってコントロールできます。SQUIDの超伝導体中のクーパー対を記述する波動関数の位相は、SQUIDを一周したとき、波動関数が安定的に存在するためにで元に戻らなければなりません。またAharonov-Bohm効果より、一周したときの波動関数の位相差は、SQUID内に外場として加えた磁場に依存します。ジョセフソン接合部分での波動関数の位相差は、波動関数がSQUIDを一周してで元に戻る条件より、に依存して変化します。これはを変化させるのと同様の効果を生みます。具体的にははのような依存性を持ちますMukaiKrantzKockum。SQUIDおよびAharonov-Bohm効果に関してはAppendix 2をご参照ください。
このようにして超伝導量子回路で量子ゲートを作成することができます。
シリーズ全体のまとめ
今回が「量子コンピュータにおける素因数分解」の最終回なので、全体の総括をしておきます。
量子コンピュータで最も有名であろうアルゴリズム「Shorのアルゴリズム」に関して説明してきました。このアルゴリズムで重要な量子回路は、べき剰余を計算する回路:
および量子フーリエ変換:
を行う回路です。これら回路を用いて
を満たすを見つけます。このが偶数かつ「適切」であれば、を古典コンピュータで計算することによりの因数を見つけることができます。量子回路による計算量、古典コンピュータによる計算量、および試行錯誤が必要な確率的要素を加味しても、の多項式程度の計算量で因数分解を行うことができます。現在まで、古典コンピュータではこのような高速なアルゴリズムは発見されていません。
上記量子回路で用いられている基本的な量子ゲートの実装方法に関しては、本記事で説明したように、いまのところ超伝導量子回路を用いた方法が主流です。クーパー対が0または1コの量子状態間の遷移をSQUIDで調整することにより、任意の1 qubitの状態を実現することができます。
現在のところ、Shorのアルゴリズムによる素因数分解は2桁程度の整数に留まっているようです。現在最も安全な暗号化技術のひとつである、素因数分解の擬似非可逆性を用いたRSA暗号を破るのに有用な桁数の素因数分解はまだまだ先のようです。しかしながら、Shorのアルゴリズムはその古典計算機に対する本質的な優位性、暗号解読という非常に重要な分野への応用、また様々な量子アルゴリズム開発の指針となるような発想を与えたという意味で、理論的なレベルでも十分価値のあるものだと思います。
おしまい。
Appendix 1 : ジョセフソン効果
本章では超伝導量子回路において基本的な現象であるジョセフソン効果について述べます。特にジョセフソン電流がどうして流れるか、またジョセフソン電流によるエネルギーであるジョセフソンエネルギーについて述べます。参考文献としてRefs.NiwaWikiJosTsunetoを挙げておきます。
ジョセフソン電流
以下Ref.Niwaに沿って説明します。
超伝導体は1つの巨大な(人間のスケールからすると小さいですが)波動関数で記述できます。ただし2つの超伝導体ではその位相が異なるとしましょう。いま空間1次元で考えます。2つの超伝導体を1,2とラベルし、それぞれの波動関数をとして、さらにこれらを極座標で表示します:
これらの波動関数はクーパー対に対する波動関数です。
この系のシュレーディンガー方程式がどうなるか考えます。いま2つの超伝導体は薄い絶縁体により弱く結合しています。全系のハミルトニアンは
のように書かれます。全形のハミルトニアンは
のように書かれます。ここではそれぞれ超伝導体1,2が孤立している場合のハミルトニアン、はそれらが弱く結合させるハミルトニアンです。いま両方の超伝導体の間に時間に依存しない電圧を印加したとします。のそれぞれのエネルギー固有値をとして、接合の両端におけるの差は
となります。2eはクーパー対の電荷です。いまエネルギーの原点をのようにエネルギー差の半分のところにとれば、エネルギーの固有値方程式:はそれぞれの超伝導体で
となります。さらにに対する仮定として、に対するシュレーディンガー方程式に対し、をに置き換えます。ただしに対しそれぞれとします。すなわちの超伝導体に関する波動関数に対するSEを作るときには、物理系のセットアップから定まるあるという定数を「相手の波動関数」にかけた相互作用ハミルトニアンが加わるとします。すると
が成立します。ここでドットは時間による偏微分です。
このSEに極座標に分解したを入れて解きます。位相差をとします。上記SEを実部と虚部に分けたうち、虚部のほうをみると
を得ます。ここでも2eはクーパー対の電荷です。クーパー対がこの2つの超伝導体の間を移動すると電流が流れますが、このときの電流は、面積の接合面を流れるとすれば
となります。
ここで注目すべきは、この方程式はには依存していないことです。すなわち2つの超伝導体に位相差さえ存在していれば電流が流れます。もし超伝導体1,2の初期の位相差が時間・空間に関し一定ならば、は直流電流です。これは直流ジョセフソン効果 (DC Josephson effect) と呼ばれます。ふつう電流と電圧の関係は抵抗をとしてオームの法則で記述されますが、ジョセフソン電流に電位差は必要ないし抵抗も関係ありません。
次に実部を解析します。実部の方程式より
を得ます。1,2の超伝導体は同素材で、位相だけが異なるとすればが成立します。このとき上の式の右辺第2項は消えて
と非常に簡単になります。これより電圧が印加された場合、は時間に関して線形の依存性をもつことがわかります:
これとEq.(1)より、電圧が印加された場合は周期の交流電流が流れることがわかります。これを交流ジョセフソン効果 (AC Josephson effect) とよび、またこの周期をジョセフソン周波数 (Josephson frequency) と呼びます。
ジョセフソン効果は、超伝導の微視的理論であるBCS理論からも、また対称性の自発的破れが生じた際に現れる秩序変数に関する一般論であるGinzburg-Landau理論からも導出できますNiwa。
ジョセフソンエネルギー
ここではジョセフソン電流によるエネルギーを計算します。これはジョセフソン接合における位相差がからまで変化した際に流れるジョセフソン電流のエネルギーです。Ref.WikiJosを参考にしています。
位相差をとして、Eq.(A1)を
と書いておきます。またEq.(A2)より
です。が時刻で、でとします。この間に流れたジョセフソン電流によるエネルギーは
以上より
を得ます。
Appendix 2 : SQUIDとAharonov-Bohm (AB) 効果
SQUIDとは超伝導量子干渉計(Superconducting Quantum Interference Device)の略です。本文図3にあるように、超伝導ループの中に2つのジョセフソン接合を持つデバイスです。通常極めて弱い磁場の検出に用いられる磁気センサーとして用いられますがWiki、本文で述べたように、量子コンピュータでは基底状態と励起状態の振動周期の調整に用いられます。
ここではなぜSQUIDに通す磁場によりその中を流れる電流を制御できるかについて説明します。これはAharonov-Bohm効果と関係します。
Aharonov-Bohm効果
以下Ref.Igiに沿って議論します。Ref.Wachterも参考文献として挙げておきます。
電荷の粒子の波動関数をとします。電場・磁場があるとき、SEは以下のようになります:
ここではそれぞれベクトルポテンシャル、スカラーポテンシャルと呼ばれます。これらの量を組にしてのように4次元のベクトルとしたものはゲージ場と呼ばれます。ゲージ場と電場と磁場との関係は以下です:
Eq.(A3)はを場所・時刻に依存する任意の関数として
に対して不変です。この変換をゲージ変換、この変換に対して不変なことをゲージ不変性、をゲージ関数と呼びます。観測可能な量はすべてゲージ不変です。電場・磁場は観測可能なのでゲージ不変です。「物理量はゲージ不変であれ」という原理はゲージ原理と呼ばれます。ゲージ原理を指針として、SEのゲージ不変性を要請することでEq.(A3)を導くこともできます。を結合させ、かつゲージ不変にするには、時間および空間微分を
のように置き換えればよいことが知られています。このようにゲージ場と粒子との結合を得る方法はminimal couplingなどと呼ばれます。
ここで重要なのは、量子論においては電場・磁場ではなくゲージ場が本質的に現れることです。局所的な場の理論では電場・磁場のみで運動方程式を書くことは不可能です。
量子論におけるゲージ場の重要性が如実に現れるのがAharanov-Bohm効果です。いま電場は無視してとします。Eq.(A3)でとした方程式
となりますが、ここでを
とします。ここではを満たす関数とします。このときはゼロです。これは直接計算でももちろんわかりますが、次のような考察からもわかります:はからゲージ関数によるゲージ変換をしたベクトルポテンシャルですが、で磁場がゼロなのは明らかなので、による磁場もそのゲージ不変性からゼロになります。
でEq.(A5)を書き直すと
となります。すなわちがノンゼロでも、それが与えるがゼロならば、は自由なSEに従います。波動関数のoverallな位相は物理的ではないので、結局この状況では粒子は自由運動をすることがわかります。
さて、を満たすは線積分を用いて
と書けることは、これをに代入すればすぐにわかります。Eq.(A7)とEq.(A6)と合わせると、ベクトルポテンシャルの存在の元で運動する粒子には、その粒子の軌道上に存在するベクトルポテンシャルに依存した位相がかかることがわかります。
上記したように位相そのものは観測量ではありません。しかし位相差は観測にかかる場合があります。ゲージ場が波動関数に及ぼす影響を見るために、以下の図のような2重スリットの状況を考えます。
磁場の存在する2重スリット実験。荷電粒子(電荷)がからに到達する様子。はそれぞれスリットA,Bを通る場合の到達経路
ゲージ場が存在しない場合、経路を通りからまで移動した粒子のにおける波動関数をそれぞれとします。一方ゲージ場が存在する場合の同様の波動関数をそれぞれとします。するとEq.(A6)(A7)より
となります。よってゲージ場が存在する下で、における波動関数は、経路を通った波動関数の和になるので
となります。ここでです。ストークスの定理を用いれば
となります。ここではを境界とする2次元面です。はを貫く磁束です。よって、粒子がからに到達する確率は、その経路を境界とする2次元面を通る磁場により変化します。Eq.(A8)はゲージ変換で不変であることに注意してください。
粒子がまたはを通りからに至る確率密度はです。軌道の中に磁場がなければ
です。一方内に磁場が存在すれば
となります。最後の2項がを貫く磁束に依存します。よってに粒子が至る確率もに依存します。これがAharanov-Bohm(AB)効果です。上記したようにAB効果はゲージ不変な現象です。
超伝導ループとSQUID
ここまでは電子を2重スリットに通す場合を考えましたが、同様なことを超伝導体のループにあてはめます。超伝導体は1つの巨視的な(クーパー対を記述する)波動関数で表せます。波動関数が安定して存在するには、それが超伝導体を一周すると元に戻らなければなりません。ところが磁場をループ内に通すと、波動関数が一周する間にAB効果により波動関数にという余計な位相がつきます。Bohr-Sommerfeldの量子化条件と同様、超伝導体内を一周した波動関数が元に戻る条件:
を満たす必要があります。は超伝導体ループを境界とする2次元面です。このように超伝導体ループではループを貫く磁場は離散化された値を取ります。いわゆる磁束の量子化が起こります。
SQUIDでは超伝導体ループの途中に存在するジョセフソン接合にさらに位相差が存在するので
が成立します。故に超伝導体の位相差を外部からSQUIDにかける磁束により調節できます。これによりジョセフソン電流はの依存性をもち、ひいてはEq.(8)のは(有効的に)により変化します。
余談ですが、Aharanov-Bohm効果は、その存在を多くの物理学者が認めるまでに理論的提案から30年ほどかかりました。その間に幾多の反論を受けました。AB効果の存在を示す決定的な実験が日立の基礎研において行われたことは有名です。それに至るまでの経緯およびその実験に関しては外村彰さんの書かれたRef.Tonomuraがわかりやすいです。なぜAB効果に関してそれほど批判があり認め難いのかに関してもこの論文を読まれるとよいと思います。