(57)【要約】【課題】多量子ビットゲートを実行できる、光共振器を用いた量子計算方法を提供する。【解決手段】(N+1)個の物理系（5つのエネルギー準位|0>、|1>、|2>、|3>、|4>を有し、量子ビットが|0>、|1>で表現される）を|2>と|3>の遷移周波数（各物理系で実質的に等しい）に共鳴する共振器モードを有する光共振器内に配置し、N個の制御系と1個の標的系を準備し、標的系に|0>、|1>、|2>と|4>に共鳴する光パルスを照射し、標的系の重ね合わせの状態|c>を|2>へ変化させ、(N+1)個の全ての物理系に|0>と|3>、|1>と|3>に共鳴する光パルスを照射して、制御系に共鳴する光パルスが照射されている間に標的系に共鳴する光パルスの位相をユニタリー変換Uに依存した特定の値だけシフトさせ、標的系に|0>、|1>、|2>と|4>に共鳴し、互いの位相差がユニタリー変換Uに依存した特定の値に設定された光パルスを照射して|2>を|c>へ戻すことにより、多制御ユニタリーゲートを実行する。【選択図】      図２
2005134761.tif 000002
【特許請求の範囲】
【請求項１】
N個（Nは自然数）の制御ビットに依存して他の1個の標的ビットにユニタリー変換Uを施す多制御ユニタリーゲートを実行できる量子計算方法であって、
(N+1)個の物理系（ここで、各々の物理系は、5つのエネルギー準位|0>、|1>、|2>、|3>、|4>を有し、量子ビットが下3準位のうち2つの準位|0>、|1>で表現され、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移が光学的に可能であり、それぞれの物理系の|2>と|3>の2準位間の遷移周波数が実質的に等しい）を、前記|2>と|3>の2準位間の遷移周波数に共鳴する共振器モードを有する光共振器内に配置し、前記(N+1)個の物理系のうちN個の物理系をN個の制御系として用い、他の1個の物理系を標的系として用いるように準備し、
前記標的系に、|0>と|4>、|1>と|4>、および|2>と|4>に共鳴する光パルスを照射し、ユニタリー変換Uに依存する前記標的系の|0>と|1>との重ね合わせの状態|c>を|2>へ変化させ、
(N+1)個の全ての制御系および標的系に|0>と|3>、および|1>と|3>に共鳴する光パルスを照射して、前記制御系に共鳴する光パルスが照射されている間に前記標的系に共鳴する光パルスの位相をユニタリー変換Uに依存した特定の値だけシフトさせ、
前記標的系に、|0>と|4>、|1>と|4>、および|2>と|4>に共鳴し、互いの位相差がユニタリー変換Uに依存した特定の値に設定された光パルスを照射して、|2>を|c>へ戻す
ことを特徴とする量子計算方法。
【請求項２】
前記5つのエネルギー準位は原子またはイオンの核スピンに起因する超微細構造分裂によって生じる準位であることを特徴とする請求項１に記載の量子計算方法。
【請求項３】
原子またはイオンに磁場または電場を印加してエネルギー準位の分裂を生じさせることによって、前記5つのエネルギー準位を設定することを特徴とする請求項１に記載の量子計算方法。
【請求項４】
|0>と|4>、|1>と|4>、および|2>と|4>に共鳴する光を、その強度および位相を制御して照射することによって1量子ビットゲートを実行することを特徴とする請求項１に記載の量子計算方法。
【請求項５】
(N+1)個の物理系（ここで、各々の物理系は、5つのエネルギー準位|0>、|1>、|2>、|3>、|4>を有し、量子ビットが下3準位のうち2つの準位|0>、|1>で表現され、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移が光学的に可能であり、それぞれの物理系の|2>と|3>の2準位間の遷移周波数が実質的に等しい）と、
その内部に前記物理系を配置させ、前記|2>と|3>の2準位間の遷移周波数に共鳴する光共振器と、
各物理系に対し、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移に共鳴する光パルスを生成する光源および周波数・振幅変調器と、
前記光パルスの位相を制御する位相変調器と、
前記周波数・振幅変調器および位相変調器を同期して動作させる制御装置と
を具備したことを特徴とする量子計算機。
【請求項６】
前記物理系のエネルギー準位間の遷移に共鳴する光が照射された物理系からの発光を検出する光検出器を有することを特徴とする請求項５記載の量子計算機。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光を用いた量子計算方法および量子計算機に関する。
【背景技術】
【０００２】
近年、量子力学の原理を使って計算を行う量子計算機の研究が盛んに行われている。量子計算機では、2つの量子状態|0>、|1>が情報の基本となる。これらは現在の計算機において情報の基本となるビットの0と1に相当しているが、量子計算機では計算過程においてα|0>＋β|1>（α、β：複素数）のような重ね合わせの状態が使われるので、現在の計算機のビットと区別してこれを量子ビットと呼ぶ。
【０００３】
量子計算における任意の演算は、1つの量子ビットを変換する1量子ビットゲートと、2つの量子ビットのうち一方の状態に依存して他方の状態を変換する2量子ビットゲートの組み合わせで実行できることが知られている。特に、2量子ビットゲートとしては、一方が|0>のときに他方をそのままとし|1>のときに他方の|0>と|1>を交換するCNOTゲートがあれば十分であることが知られている（非特許文献１参照）。これらの基本的なゲートをユニバーサルゲートと呼ぶ。ユニバーサルゲートを用いれば任意の演算を行うことが可能ではあるが、3量子ビット以上に作用するゲート（2つの量子ビットが|1>のときにのみ他の1つの量子ビットの|0>と|1>を交換する量子トフォリゲートなど）をユニバーサルゲートの組み合わせで行おうとすると多数のゲートが必要になり、操作も複雑になる。よって、Shorの素因数分解のアルゴリズムやGroverのデータベース検索のアルゴリズムなど、3量子ビット以上に作用するゲートを多用する実際の量子アルゴリズムを実現するには、ユニバーサルゲートだけでなく3量子ビット以上に作用するゲートもユニバーサルゲートに分解せずに実行できたほうがよい。
【０００４】
現在までに、いくつかの2量子ビットゲートの実現方法が発表されている。光共振器を用いて、離れた物理系を結合して2量子ビットゲートを行う方法の最初のものとして、Pellizzariの方法が知られている（非特許文献２参照）。この方法で用いられる物理系は1つ1つがそれぞれ2重に縮退した3つのエネルギー準位を持つ3準位系であり、下2準位のうち一方の縮退した2状態で量子ビットを表現し、下2準位の残りの準位と上準位が共振器モードに共鳴している。この方法では縮退した準位を使うので、縮退した2状態に個別の操作をしたい場合には、磁場または電場をかけて縮退を解いてエネルギーで区別するか、物理系として用いる材料固有の光学的選択則を利用するか、いずれかをしなければならないという問題があった。このような問題を持たない共振器を用いた提案もある（非特許文献３参照）。しかし、Pellizzariの方法も含めて従来の方法はすべて2量子ビットゲートの提案であり、3量子ビット以上に作用するゲートを行う共振器を用いた量子計算方法の提案は今のところない。
【非特許文献１】M. A. Nielsen and I. L. Chuang, Quantum Information and Computation, Cambridge Univ. Press （2000）
【非特許文献２】T. Pellizzari, S. A. Gardiner, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 75, 3788 （1995）
【非特許文献３】L. You, X. Y. Yi, and X. H. Su, Phys. Rev. A 67, 032308 （2003）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
本発明の目的は、2量子ビットゲートを含む多量子ビットゲートを実行できる、光共振器を用いた量子計算方法およびそれを実現する量子計算機を提供することにある。
【課題を解決するための手段】
【０００６】
本発明の一態様に係る量子計算方法は、N個（Nは自然数）の制御ビットに依存して他の1個の標的系量子ビットにユニタリー変換Uを施す多制御ユニタリーゲートを実行できる量子計算方法であって、
(N+1)個の物理系（ここで、各々の物理系は、5つのエネルギー準位|0>、|1>、|2>、|3>、|4>を有し、量子ビットが下3準位のうち2つの準位|0>、|1>で表現され、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移が光学的に可能であり、それぞれの物理系の|2>と|3>の2準位間の遷移周波数が実質的に等しい）を、前記|2>と|3>の2準位間の遷移周波数に共鳴する共振器モードを有する光共振器内に配置し、前記(N+1)個の物理系のうちN個の物理系をN個の制御系として用い、他の1個の物理系を標的系として用いるように準備し、
前記標的系に、|0>と|4>、|1>と|4>、および|2>と|4>に共鳴する光パルスを照射し、ユニタリー変換Uに依存する前記標的系の|0>と|1>との重ね合わせの状態|c>を|2>へ変化させ、
(N+1)個の全ての制御系および標的系に|0>と|3>、および|1>と|3>に共鳴する光パルスを照射して、前記制御系に共鳴する光パルスが照射されている間に前記標的系に共鳴する光パルスの位相をユニタリー変換Uに依存した特定の値だけシフトさせ、
前記標的系に、|0>と|4>、|1>と|4>、および|2>と|4>に共鳴し、互いの位相差がユニタリー変換Uに依存した特定の値に設定された光パルスを照射して、|2>を|c>へ戻すことを特徴とする。
【０００７】
本発明の他の態様に係る量子計算機は、
(N+1)個の物理系（ここで、各々の物理系は、5つのエネルギー準位|0>、|1>、|2>、|3>、|4>を有し、量子ビットが下3準位のうち2つの準位|0>、|1>で表現され、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移が光学的に可能であり、それぞれの物理系の|2>と|3>の2準位間の遷移周波数が実質的に等しい）と、
その内部に前記物理系を配置させ、前記|2>と|3>の2準位間の遷移周波数に共鳴する光共振器と、
各物理系に対し、下3準位|0>、|1>、|2>のそれぞれと上2準位|3>、|4>のそれぞれとの間の遷移に共鳴する光パルスを生成する光源および周波数・振幅変調器と、
前記光パルスの位相を制御する位相変調器と、
前記周波数・振幅変調器および位相変調器を同期して動作させる制御装置と
を具備したことを特徴とする。
【発明の効果】
【０００８】
本発明によれば、制御ビットの条件を表す基底を自由に選択して、単一制御ユニタリーゲートだけでなく、多制御ユニタリーゲートも実行することができる。この結果、基本的なユニバーサルゲートのみを使った場合に比べて、実際の量子アルゴリズムの実現で必要なゲートの数が格段に減り、操作も簡略化されるという効果が得られる。
【発明を実施するための最良の形態】
【０００９】
以下、本発明を詳しく説明する。  
本明細書において、制御系とは制御ビットを保持するために用いる物理系のことを呼び、標的系とは標的ビットを保持するために用いる物理系のことを呼ぶ。簡単のため、2つの物理系を制御系として用い、他の1つの物理系を標的系として用い、各制御ビットに依存して標的ビットにユニタリー変換Uを施す多制御ユニタリーゲート（以下、M-C-Uゲートという）をどのように実現するかについて説明する。なお、後述するように、制御ビットがN個（N:2以外の自然数）の場合も同様にして実現できる。
【００１０】
以下、「量子ビットが|d>の場合〜をする」というのは、「量子ビットを表す状態を|d>とそれと直交する状態|c>の重ね合わせで表したとき、その重なり合った状態のうち|d>を含む状態に対して〜を行う」、という意味である。
【００１１】
図１に、ここで説明するM-C-Uゲートの量子回路図を示す。図１に示すように、2つの制御系を物理系1と2、標的系を物理系3とし、2つの制御ビットを量子ビット1と2、標的ビットを量子ビット3とする。このゲートは、各制御ビットがそれぞれ|d>1=U1-1|1>1と|d>2=U2-1|1>2の場合にのみ、標的ビットにユニタリー変換U3を施す（図１では、量子ビット3を標的ビットとしているので、上記のユニタリー変換UをU3と表している）。この場合、このM-C-Uゲートにおける「制御ビットkの条件を表す基底」とは、|d>kとそれと直交する状態|c>kの組{|d>k ,|c>k }のことを意味する。
【００１２】
この特別な場合として、各制御ビットがそれぞれ|1>1と|1>2の場合にのみ標的ビットの|0>3と|1>3を交換する量子トフォリゲートがある。量子トフォリゲートはショアが発見した素因数分解の量子アルゴリズムにおいて多用される非常に重要なゲートである(V. Vedral, A. Barenco, and A. Ekert, Phys. Rev. A 54, 147 (1996))。
【００１３】
3つのユニタリー変換は一般に次のように表せる（任意のユニタリー変換は以下の形で表すことができる）。
【数１】
2005134761.tif 000003
【００１４】
ここで、I(3), σx(3), σy(3), σz(3)は量子ビット3に対する恒等演算子およびパウリ演算子であり、次式で定義される。
【数２】
2005134761.tif 000004
【００１５】
また、Uk(k=1,2)の位相因子は、後に元に戻すユニタリー変換Uk-1を行うことから無視できるので、Θk=φk=0(k=1,2)とする。
【００１６】
図２に、3つの物理系のエネルギー準位とそれらを含む光共振器および共振器モードを示す。物理系k(k=1,2,3)の状態をエネルギーの低いほうから順に|0>k、|1>k、|2>k、|3>k、|4>kとする。|0>kと|1>kによって量子ビットを表す。ゲート操作前は、全ての物理系の状態は|0>kと|1>kで表される状態であるとする。|2>kと|3>kが共振器モードと共鳴している。上記のM-C-Uゲートは以下に説明する3つのステップにより実現される。図３に各ステップにおいて照射する光を示す。
【００１７】
始状態|Ψ0>kを、|d>k、|c>k(k=1,2,3)と共振器モードの光子状態を用いて下記のように表す。
【数３】
2005134761.tif 000005
【００１８】
ここで、|lmn>≡|l>1|m>2|n>3であり、|0>は共振器モードの光子数がゼロの状態（真空）を表す。
【００１９】
（ステップ1）
図３（ステップ1）に示すように、初めに、標的系に|0>3と|4>3、|1>3と|4>3、|2>3と|4>3に共鳴する光パルス列を照射する。量子ビットは位置または周波数で区別する。それぞれの光パルス列の周波数をω04(3)、ω14(3)、ω24(3)とする。物理系3（標的系）に作用する光パルス列の電場は次式で表される。
【数４】
2005134761.tif 000006
【００２０】
エンベロープは次式で与えられる。
【数５】
2005134761.tif 000007
【００２１】
ここで、A(3)は定数であり、fc(t), f2(t)はパルス波形を表し、次式で与えられる。
【数６】
2005134761.tif 000008
【００２２】
上式において2τはパルス幅、t0は時刻0において電場がゼロとなるようにするための正定数である。パルスの間隔を幅の半分に等しくしたのは、以下で説明するアディアバティック・パッセッジという物理過程が効率よく行われることを考慮しているためである。光パルスの具体例については後述する（実施例1、図８参照）。
【００２３】
相互作用表示における相互作用ハミルトニアンは、回転波近似のもとで次式のように表される。
【数７】
2005134761.tif 000009
【００２４】
ハミルトニアンのゼロ固有値に対する固有状態をダークステート（dark state）と呼ぶ。上記のハミルトニアンに対するダークステートで始状態と関係のあるものは、
|l>1|m>2|d>3|0>, |l>1|m>2(Ω2(3)(t)|c>3-Ωc(3)(t)|2>3)|0> (l,m=0,1)
である。ラビ周波数の時間変化が十分遅いとき、始状態のうち|c>3は上記の光パルス列によってダークステートの変化とともに変化する。それ以外は変化しない。|2>3に共鳴する光パルスが照射され始めたときは、ダークステートは|c>3なので|c>3は変化しない。その後、|c>3に共鳴する光パルスが照射され始めるとダークステートは|c>3と|2>3の重ね合わせとなるので、始状態のうち|c>3も同じ重ね合わせの状態へと変化する。光パルス列の照射が終了したとき、ダークステートの変化に伴い|c>3は最終的に-|2>3となる。このようにダークステートを通して状態が変化する現象を「アディアバティック・パッセッジ（adiabatic passage）」と呼ぶ(K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. Phys. 70, 1003 (1998))。上記のアディアバティック・パッセッジの後の状態は、
|Ψ1>=(αddd|ddd>-αddc|dd2>+αdcd|dcd>-αdcc|dc2>
+αcdd|cdd>-αcdc|cd2>+αccd|ccd>-αccc|cc2>)|0>
となる。
【００２５】
（ステップ2）
次に、共振器効果を用いたアディアバティック・パッセッジを行う。図３（ステップ2）に示すように、|0>1と|3>1、|1>1と|3>1、|0>2と|3>2、|1>2と|3>2、|0>3と|3>3、|1>3と|3>3に共鳴する光パルス列を照射する。それぞれの光パルス列の周波数をω03(1)、ω13(1)、ω03(2)、ω13(2)、ω03(3)、ω13(3)とする。物理系k(k=1,2,3)に作用するパルス列の電場E'(k)(t)は次式で与えられる。
【数８】
2005134761.tif 000010
【００２６】
エンベロープは次式で与えられる。
【数９】
2005134761.tif 000011
【００２７】
ここで、A'(k)(1,2,3)は定数であり、パルス波形は次式で与えられる。
【数１０】
2005134761.tif 000012
【００２８】
ステップ1のパルスと重ならないように、ステップ1とステップ2の隣接するパルスはたとえば6τ離すように設定する。また、f1(3)(t)とf2(3)(t)も、重ならないように6τ離すように設定する。以下においても、離すべきパルスは6τ離すように設定する。相互作用ハミルトニアンは次式のように表される。
【数１１】
2005134761.tif 000013
【００２９】
状態の変化に関係のあるダークステートは次のようになる。
【数１２】
2005134761.tif 000014
【００３０】
上記の光パルス列によるアディアバティック・パッセッジにおいて、Ω'c(3)の位相がφ3だけシフトするために、状態は次式のように変わる。
【数１３】
2005134761.tif 000015
【００３１】
（ステップ3）
図３（ステップ3）に示すように、ステップ1のアディアバティック・パッセッジの逆を行うことによって-|2>3→e-iφ3|c>3と戻す。ステップ1との違いは、エンベロープおよび光パルス波形を次のようにすることである。なお、定数A(3)はステップ1と同じでなくてもよい。
【数１４】
2005134761.tif 000016
状態は、
|Ψ3>=(αddd|ddd>+αddce-iφ3|ddc>+αdcd|dcd>+αdcc|dcc>
+αcdd|cdd>+αcdc|cdc>+αccd|ccd>+αccc|ccc>)|0>
となる。
【００３２】
これを演算子Pで始状態を変換する形で表すと、次のようになる。
【数１５】
2005134761.tif 000017
【００３３】
関係式：|d>33<d|+e-iφ3|c>33<c|=e-i(Θ3+φ3/2)U3を用いれば、演算子Pは次のように表される。
【数１６】
2005134761.tif 000018
【００３４】
これは2つの制御ビットがそれぞれ|d>1=U1-1|1>1、|d>2=U2-1|1>2のときにのみ標的ビットにe-i(Θ3+φ3/2)U3をかけることを意味しているので、位相を除いて実現したかったM-C-Uゲートになっている。Θ3+φ3/2=0の場合は位相まで含めて正しく実行できる。最も重要なM-C-Uゲートの1つである量子トフォリゲートは、パラメータの値を
φ3=π, ψ1=ψ2=ψ3=π, θ1=θ2=0, θ3=π/2, Θ3=-π/2
とすることで、位相を含めて正しく実行できる（Θ3+φ3/2=0を満たしている）。
【００３５】
次に、制御ビットが1つであるC-Uゲートを実行する場合について説明する。ここでは、制御ビット1と標的ビットのみを考え、制御ビット2はないものとする。この場合、ステップ3で-|2>3→e-iφ3|c>3を行っている間に、制御ビット1に|d>1→ei[Θ3+φ3/2]|d>1（|c>1はそのまま）という1量子ビットゲートを行えば位相まで含めて常に正しく実行できる。1量子ビットを変換する1量子ビットゲートは、下3準位|0>k、|1>k、|2>kと上準位|4>kに共鳴する光を使ったアディアバティック・パッセッジによって実現できることが知られている(Z. Kis and F. Renzoni, Phys. Rev. A 65, 032318 (2002))。この方法を使って、アディアバティック・パッセッジによって|d>1→ei[Θ3+φ3/2]|d>1を行う。同様にして、任意の1量子ビットゲートは下3準位|0>k、|1>k、|2>kと上準位|4>kに共鳴する光を使ったアディアバティック・パッセッジによって実現できる。このようにした場合、本発明に係る量子計算方法において、演算に用いられる物理的な過程は全てアディアバティック・パッセッジとなる。アディアバティック・パッセッジは、励起状態からの緩和を無視でき、かつ制御すべき物理量の誤差に鈍感である、といったメリットがある。このため、演算が全てアディアバティック・パッセッジで行えることは本発明の利点の1つである。
【００３６】
|d>1→ei[Θ3+φ3/2]|d>1という制御ビット1に対する1量子ビットゲートを実行するには、図４に示すように、制御系1に|0>k、|1>k、|2>kと|4>kに共鳴する光パルス列を照射する。それぞれの光パルス列の周波数をω04(1)、ω14(1)、ω24(1)とする。光パルス列の電場は次式で表される。
【数１７】
2005134761.tif 000019
【００３７】
エンベロープは次式で与えられる。
【数１８】
2005134761.tif 000020
【００３８】
ここで、A"(1)は定数であり、パルス波形は次式で与えられる。
【数１９】
2005134761.tif 000021
【００３９】
次いで、制御ビットがN個（N:自然数）の多制御ユニタリーゲートを実行する場合について説明する。ここでは、1〜N番目の量子ビットを制御ビット、(N+1)番目の量子ビットを標的ビットとする。ステップ1では上で説明したM-C-Uゲートの場合と同様に|c>N+1→-|2>N+1とする。今の場合、ステップ2におけるハミルトニアンは次式で与えられる。
【数２０】
2005134761.tif 000022
【００４０】
対応するダークステートのうち時間に依存するものは、次にように表される。
【数２１】
2005134761.tif 000023
【００４１】
ステップ1の後の状態に含まれる状態のうち、標的ビットが|2>で制御ビットの少なくとも1つが|c>である状態は、上記ダークステートに含まれ、アディアバティック・パッセッジにより変化する。上記の式より、それらの状態のダークステート内での係数はgkΩ'c(N+1)(k=1,…,N)になることがわかる。ステップ2のアディアバティック・パッセッジでは、ラビ周波数Ω'c(N+1)をeiφ(N+1)Ω'c(N+1)へ変換することを行うので、これらの状態の係数にはeiφ(N+1)がかかる。ステップ1の後の状態に含まれる状態のうち、制御ビットが全て|d>である状態は時間に依存しないダークステートとなっており、変化しない。残りの状態（標的ビットが|d>で制御ビットの少なくとも1つが|c>である状態）は上記ダークステートに含まれるが、その係数はgのみのべき乗であり、時間によらないので、ステップ2のアディアバティック・パッセッジにより変化しない。ステップ3において、上で説明したM-C-Uゲートの場合と同様に-|2>N+1→e-iφ(N+1)|c>N+1とすれば、最終的に得られる状態は制御ビットがすべて|d>の場合にのみ標的ビットにe-i(Θ(N+1)+φ(N+1)/2)U3を施したものとなるので、以上の3つのステップによりN制御ビットのM-C-Uゲートも（位相を除いて）実現できることがわかる。
【００４２】
制御ビットがN個（N：2以上の自然数）の場合でΘN+1+φN+1/2≠0であり位相まで正しく実行したい場合は、上で説明したM-C-Uゲートを実行した後、制御ビット1を標的ビット、それ以外の(N-1)個の制御ビットを制御ビットとしたM-C-Uゲートを同じように実行して|d>1|d>2…|d>N→ei[Θ(N+1)+φ(N+1)/2]|d>1|d>2…|d>N（その他はそのまま）とすればよい。
【００４３】
以上のように、N個の制御ビットと他の1つの標的ビットに対し、各制御ビットが|0>と|1>のある重ね合わせの状態|d>（|d>は各制御ビットで異なっていてもよい）の場合にのみ標的ビットにユニタリー変換Uを施すM-C-Uゲートはアディアバティック・パッセッジを3回行うことで（位相を除いて）実現できる。
【００４４】
（実施例）
以下、本発明の実施例を説明する。
【００４５】
（実施例1）
図５に、本実施例における量子計算機の構成を示す。試料20としてPr3+をドープしたY2SiO5結晶（Pr3+:Y2SiO5結晶）を用意し、物理系としてPr3+の超微細準位を用いる。図６に、Pr3+の超微細準位を示す。これはPr3+の核スピンに起因する超微細構造分裂によるもので、低温において非常に長い緩和時間を実現できる。そのため、コヒーレントな重ね合わせの状態を、レーザーを用いて作れることが確認されている(K. Ichimura, K. Yamamoto, and N. Gemma, Phys. Rev. A 58, 4116 (1998))。このPr3+の超微細準位は、下3準位、上2準位を使う本発明の量子計算方法に最適である。
【００４６】
なお、適切なエネルギー準位のイオンがない場合には、各準位が核スピンについて縮退しているので、後述するように（実施例3）磁場を印加してゼーマン分裂を起こし、分裂の大きさを調整したゼーマン準位を用いてもよい。また、使用する物理系によっては、電場を印加することにより所望のエネルギー準位を準備してもよい。
【００４７】
光源10としては周波数ジッターが数キロヘルツにまで周波数安定化したリング色素レーザーを用いる。本発明の量子計算方法に必要な光パルス列は、周波数シフトと強度・位相変調によって実現できるので、これらを行うために音響光学効果素子（AOM）11と電気光学効果素子（EOM）12を使う。光パルスの振幅をAOM11で、その位相をEOM12で適切に設定できるように、AOM11とEOM12を同期させて動作させる制御系15を用いる。
【００４８】
試料20（Pr3+:Y2SiO5結晶）は比較的小さく（1×1×1mm）、ファブリ・ペロ共振器21の焦点に固定され、共振器21ごとクライオスタット13中に設置される。クライオスタット13内は液体ヘリウムを用いて絶対温度1.4Kに保つ。
【００４９】
量子計算終了後の計算結果の読み出しにおいて、読み出したい準位に共鳴する光を照射してその発光を検出するために、クライオスタット13の外部に感度のよい光検出器14を配置する。
【００５０】
図５の装置を用いた量子トフォリゲートの動作を説明する。量子トフォリゲートは、制御ビット1と2がともに|1>の場合のみ標的ビットの|0>と|1>を入れ換えるというゲートである。
【００５１】
まず、ゲート操作を行うのに適切なイオンを選び出す。Pr3+の濃度が十分薄い試料を使い、さらに不均一幅の中心から外れた周波数帯で動作させることで、共振器に共鳴するイオンが数個であるようにする。これらのイオンのうち、図７に示された準位を持つ3つのイオンを使う。共振器モードに共鳴する2準位を|2>k、|3>k(k=1,2,3)とし、それ以外の準位のうち、2つの下準位を|0>k、|1>k、1つの上準位を|4>kと表す。ここでは周波数で分けやすいように共振器に共鳴する下準位を図７のように異なる核スピンの超微細準位になるようにイオンを選んだが、同じ核スピンの超微細準位の不均一幅の中で選んでもよい。
【００５２】
次に、適切な初期状態を用意する作業を行う。|2>k(k=1,2,3)と|3>k、|4>k以外の残りの上準位に共鳴する光を照射し、共振器モードに共鳴する全てのイオンの下準位の状態を共振器モードに共鳴しない下準位に移す。その|2>k(k=1,2,3)に共鳴する光が照射されている状態で、|1>1と|3>1、|1>2と|4>2、|1>3と|4>3に共鳴する光も照射し、3つのイオンの状態を|0>1、|0>2、|0>3に移す。こうして初期化が行われる。同様にして、3つのイオンの状態を|1>1、|0>2、|0>3または|0>1、|1>2、|0>3または|1>1、|1>2、|0>3に移すこともできる。以下においては、初期状態が(|0>1、|0>2、|0>3）、(|1>1、|0>2、|0>3）、(|0>1、|1>2、|0>3)、(|1>1、|1>2、|0>3)の4通りの場合で量子トフォリゲートを行い、その後測定をすることで正しい動作をすることを確認する。初期状態が(|1>1、|1>2、|0>3)の場合だけ量子ビット3の状態が|1>3となり、その他の場合に|0>3のままであれば正しい動作をしたといえる。
【００５３】
量子トフォリゲートを実行する具体的な操作を説明する。量子トフォリゲートと上述した一般的なM-C-Uゲートを比較すると、M-C-Uゲートの場合に使われた状態|d>1、|c>1、|d>2、|c>2とユニタリー変換U3は、量子トフォリゲートではそれぞれ|1>1、|0>1、|1>2、|0>2とσx(3)であるから、照射するパルス列を決めるパラメータの値は以下のようになる。
【００５４】
φ3=π, ψ1=ψ2=ψ3=π, θ1=θ2=0, θ3=π/2, Θ3=-π/2。
【００５５】
これより、|d>3=(|0>+|1>)*2-1/2、|c>3=(|0>-|1>)*2-1/2となる。
【００５６】
以下、始状態が|Ψ0>=|000>|0>=(|00d>+|00c>)*2-1/2|0>の場合を詳しく説明する。トフォリゲートを行う過程は、上述した一般的なM-C-Uゲートの場合と同様に、3つのステップからなる。図８に、3つのステップで照射する光パルスのエンベロープを示す。図８では、以下の式中に現れる定数因子をすべて1とした。
【００５７】
（ステップ1）
初めに、標的ビットに|0>3と|4>3、|1>3と|4>3、|2>3と|4>3に共鳴する光パルス列を照射する。物理系3に作用するパルス列の電場は上で決めたパラメータにより決まり、次式のように表される。
【数２２】
2005134761.tif 000024
【００５８】
エンベロープは次式で与えられる（図８）。
【数２３】
2005134761.tif 000025
【００５９】
ここで、A(3)は定数であり、fc(t), f2(t)はパルス波形を表し、次式で与えられる。
【数２４】
2005134761.tif 000026
【００６０】
本実施例ではτ=10μs、t0=30μsとした。相互作用表示における相互作用ハミルトニアンは回転波近似のもとで次式のように表される。
【数２５】
2005134761.tif 000027
【００６１】
上記のハミルトニアンに対するダークステート（ハミルトニアンのゼロ固有値に対応する固有状態）で始状態と関係のあるものは
|l>1|m>2|d>3|0>, |l>1|m>2(Ω2(3)(t)|c>3-Ωc(3)(t)|2>3)|0> (l,m=0,1)
である。ラビ周波数の時間変化が十分遅いとき、始状態のうち|c>3は上記のパルス列によってダークステートの変化とともに変化する。それ以外は変化しない。|2>3に共鳴する光パルスが照射され始めたときは、ダークステートは|c>3なので|c>3は変化しない。その後、|c>3に共鳴する光パルスが照射され始めるとダークステートは|c>3と|2>3の重ね合わせとなるので、始状態のうち|c>3も同じ重ね合わせの状態へと変化する。光パルス列の照射が終了したとき、ダークステートの変化に伴い、|c>3は最終的に-|2>3となる（アディアバティック・パッセッジ）。上記のアディアバティック・パッセッジの後の状態は、|Ψ1>=(|00d>-|002>)*2-1/2|0>となる。
【００６２】
（ステップ2）
次に、共振器効果を用いたアディアバティック・パッセッジを行う。|0>1と|3>1、|0>2と|3>2、|0>3と|3>3、|1>3と|3>3に共鳴する光パルス列を照射する。それぞれの光パルス列の周波数をω03(1)、ω03(2)、ω03(3)、ω13(3)とする。物理系k(k=1,2,3)に作用するパルス列の電場E'(k)(t)は次式で与えられる。
【数２６】
2005134761.tif 000028
【００６３】
エンベロープは次式で与えられる。
【数２７】
2005134761.tif 000029
【００６４】
ここで、A'(k)(1,2,3)は定数であり、パルス波形f(k)(t)(k=1,2,3)は次式で与えられる。
【数２８】
2005134761.tif 000030
【００６５】
相互作用ハミルトニアンは次式のように表される。
【数２９】
2005134761.tif 000031
【００６６】
状態の変化に関係のあるダークステートは次のようになる。
【数３０】
2005134761.tif 000032
【００６７】
上記のパルス列によるアディアバティック・パッセッジにより、Ω'c(3)の符号が反転するので状態は、|Ψ2>=(|00d>+|002>)*2-1/2|0>に変わる。
【００６８】
（ステップ3）
ステップ1のアディアバティック・パッセッジの逆を行うことによって-|2>3→-|c>3と戻す。ステップ1との違いは、ステップ1との違いは、エンベロープおよび光パルス波形を次のようにすることである。なお、定数A(3)はステップ1と同じでなくてもよい。
【数３１】
2005134761.tif 000033
【００６９】
状態は、|Ψ3>=(|00d>+|00c>)*2-1/2|0>=|000>|0>=|Ψ0>となる。すなわち、始状態が|000>|0>の場合には、以上の3ステップで元に戻り変化しない。
【００７０】
始状態が|Ψ0>=|100>|0>=(|10d>+|10c>)*2-1/2|0>の場合に図８のパルス列を照射すると、各ステップで|Ψ1>=(|10d>-|102>)*2-1/2|0>、|Ψ2>=(|10d>+|102>)*2-1/2|0>、|Ψ3>=(|10d>+|10c>)*2-1/2|0>=|100>|0>=|Ψ0>となり、3ステップで元に戻り変化しない。
【００７１】
始状態が|Ψ0>=|110>|0>=(|11d>+|11c>)*2-1/2|0>の場合に図８のパルス列を照射すると、各ステップで|Ψ1>=(|11d>-|112>)*2-1/2|0>、|Ψ2>=(|11d>+|112>)*2-1/2|0>、|Ψ3>=(|11d>+|11c>)*2-1/2|0>=|111>|0>となり、3ステップで量子ビット3がフリップする。
【００７２】
以上の動作はトフォリゲートの動作と一致しているので、図８のパルス列でトフォリゲートが実現できる。
【００７３】
AOM11とEOM12を用いて図８のエンベロープを持つパルス列を生成し、初期化後の試料に照射した。光パルスの強度はラビ周波数がピーク値で数MHzになるように設定した。その直後に|1>3と|4>3に共鳴する光パルスを照射して発光を光検出器14で検出した。この測定を繰り返し行ったところ、初期状態が(|0>1、|0>2、|0>3）、(|1>1、|0>2、|0>3）、(|0>1、|1>2、|0>3)の場合には発光がほとんど観測されず、初期状態が(|1>1、|1>2、|0>3)の場合には発光が観測された。これはゲートが正しく動作したことを意味する。
【００７４】
（実施例2）
本実施例では、物理系に磁場を印加してエネルギー準位を設定する場合について説明する。図９に、本実施例における量子計算機の構成を示す。図９の装置では、実施例1で用いた図５の装置に加え、試料20に対して磁場を印加するコイル22が設けられている。
【００７５】
試料20に磁場を印加することによって、適切なエネルギー準位が設定できることを示すために、実施例1で使用した周波数帯から共振器の自由スペクトル域(Free Spectral Range, FSR)だけ離れた周波数帯で実験を行った。この場合には、図７に示したような都合のよいエネルギー準位を持った3つのイオンが存在しなかった。そこで、印加する磁場の強度を変化させて調べたところ、磁場が数100Gのところで図７に示したような共振器の共振周波数に等しい遷移周波数のエネルギー準位を持った3つのイオンが見つかった。これらに対し、実施例2と同様にして量子トフォリゲートを実行したところ、正しい動作をすることが確認された。
【図面の簡単な説明】
【００７６】
【図１】本発明の実施形態に係る量子回路図。
【図２】3つの物理系のエネルギー準位とそれらを含む光共振器および共振器モードを示す図。
【図３】M-C-Uゲートを実行する3つのステップで照射する光を示す図。
【図４】C-Uゲートを実行する場合に位相を調整するために制御ビット1に対する1量子ビットゲートを行う際に照射する光を示す図。
【図５】本発明の実施例1において用いられた量子計算機の構成を示す図。
【図６】Pr3+:Y2SiO5結晶中のPr3+イオンのエネルギー準位を示す図。
【図７】本発明の実施例1における量子トフォリゲートに用いられた3つのPr3+イオンのエネルギー準位を示す図。
【図８】本発明の実施例1における量子トフォリゲートを実現するのに用いられた光パルス列のエンベロープを示す図。
【図９】本発明の実施例2において用いられた量子計算機の構成を示す図。
【符号の説明】
【００７７】
10…周波数安定化リング色素レーザー、11…音響光学効果素子（AOM）、12…電気光学効果素子（EOM）、13…クライオスタット、14…光検出器、15…制御装置、20…Pr3+:Y2SiO5結晶、21…光共振器、22…コイル。
【図１】
2005134761.tif 000034
【図２】
2005134761.tif 000035
【図３】
2005134761.tif 000036
【図４】
2005134761.tif 000037
【図５】
2005134761.tif 000038
【図６】
2005134761.tif 000039
【図７】
2005134761.tif 000040
【図８】
2005134761.tif 000041
【図９】
2005134761.tif 000042