(57)【要約】【課題】  量子ゲートを高効率で実行できる量子計算機を提供する。【解決手段】  一実施形態に係る量子計算機は、共振器の中に配置された物理系Ｘｉ（ｉ＝１，…,Ｎ）及び物理系Ｙｊ（ｊ＝１，…,Ｎ２）並びに光源部を備える。各物理系Ｘｉでは、状態｜ｅ＞ｉは量子ビットに用いる状態｜０＞ｉ、｜１＞ｉ及びゲート操作の補助に用いる状態｜２＞ｉよりもエネルギーが高く、｜２＞ｉ−｜ｅ＞ｉ遷移は共振器の共振器モードに共鳴している。各物理系Ｙｊでは、状態｜ｅ＞′ｊは状態｜２＞′ｊよりもエネルギーが高く、｜２＞′ｊ−｜ｅ＞′ｊ遷移は前記共振器モードに共鳴している。光源部は、物理系Ｘｓ及び物理系Ｘｔを操作するために、｜１＞ｓ−｜ｅ＞ｓ遷移及び｜１＞ｔ−｜ｅ＞ｔ遷移それぞれに共鳴する第１及び第２のレーザ光と物理系Ｙｊに関して｜２＞′ｊにポピュレーションを集めるための第３のレーザ光とを共振器に照射する。【選択図】図４
2015135377.tif 000002
【特許請求の範囲】
【請求項１】
共振器の中に配置された複数の物理系Ｘｉ（ｉ＝１，２，…,Ｎ、Ｎは２以上の整数）であって、各物理系Ｘｉは、状態｜０＞ｉ、｜１＞ｉ、｜２＞ｉ、｜ｅ＞ｉを含む少なくとも４つの状態を有し、状態｜ｅ＞ｉのエネルギーは量子ビットに用いる状態｜０＞ｉ、｜１＞ｉそれぞれのエネルギー及びゲート操作の補助に用いる状態｜２＞ｉのエネルギーよりも高く、｜２＞ｉ−｜ｅ＞ｉ遷移は前記共振器の共振器モードに共鳴している、複数の物理系Ｘｉと、
前記共振器の中に配置され、前記複数の物理系Ｘｉと異なる物理系Ｙｊ（ｊ＝１，２，…,Ｎ２、Ｎ２は１以上の整数）であって、各物理系Ｙｊは、状態｜２＞′ｊ、｜ｅ＞′ｊを含む少なくとも２つのエネルギー状態を有し、状態｜ｅ＞′ｊのエネルギーは状態｜２＞′ｊのエネルギーよりも高く、｜２＞′ｊ−｜ｅ＞′ｊ遷移は前記共振器モードに共鳴している、物理系Ｙｊと、
前記複数の物理系Ｘｉのうちの２つの物理系Ｘｓ（ｓはＮ以下の自然数）及び物理系Ｘｔ（ｔはｓと異なるＮ以下の自然数）の状態を操作するために、｜１＞ｓ−｜ｅ＞ｓ遷移に共鳴する第１のレーザ光と、｜１＞ｔ−｜ｅ＞ｔ遷移に共鳴する第２のレーザ光と、前記物理系Ｙｊに関して｜２＞′ｊ−｜ｅ＞′ｊ遷移における状態｜２＞′ｊにポピュレーションを集めるための第３のレーザ光と、を前記共振器に照射する光源部と、
を具備することを特徴とする量子計算機。
【請求項２】
前記光源部は、レーザ光を発する光源と、前記レーザ光を第１の分割レーザ光、第２の分割レーザ光、第３の分割レーザ光に分割するための２つのビームスプリッタと、前記第１の分割レーザ光を変調することにより前記第１のレーザ光を生成する第１の光変調器と、前記第２の分割レーザ光を変調することにより前記第２のレーザ光を生成する第２の光変調器と、前記第３の分割レーザ光を変調することにより前記第３のレーザ光を生成する第３の光変調器と、を備えることを特徴とする請求項１に記載の量子計算機。
【請求項３】
前記光源部は、第１の光源と、前記第１の光源から発せられるレーザ光を変調することにより前記第１のレーザ光を生成する第１の光変調器と、第２の光源と、前記第２の光源から発せられるレーザ光を変調することにより前記第２のレーザ光を生成する第２の光変調器と、第３の光源と、前記第３の光源から発せられるレーザ光を変調することにより前記第３のレーザ光を生成する第３の光変調器と、を備えることを特徴とする請求項１に記載の量子計算機。
【請求項４】
前記共振器は、前記複数の物理系Ｘｉを含む第１の空間領域と、前記物理系Ｙｉを含む、前記第１の空間領域と異なる第２の空間領域と、を備え、前記光源部は、前記第２の空間領域に前記第３のレーザ光を照射することを特徴とする請求項１に記載の量子計算機。
【請求項５】
共振器の中に配置された複数の物理系Ｘｉ（ｉ＝１，２，…,Ｎ、Ｎは２以上の整数）であって、各物理系Ｘｉは、状態｜０＞ｉ、｜１＞ｉ、｜２＞ｉ、｜ｅ＞ｉを含む少なくとも４つの状態を有し、状態｜ｅ＞ｉのエネルギーは量子ビットに用いる状態｜０＞ｉ、｜１＞ｉそれぞれのエネルギー及びゲート操作の補助に用いる状態｜２＞ｉのエネルギーよりも高く、｜２＞ｉ−｜ｅ＞ｉ遷移は前記共振器の共振器モードに共鳴している、複数の物理系Ｘｉを用意し、
前記共振器の中に配置され、前記複数の物理系Ｘｉと異なる物理系Ｙｊ（ｊ＝１，２，…,Ｎ２、Ｎ２は１以上の整数）であって、各物理系Ｙｊは、状態｜２＞′ｊ、｜ｅ＞′ｊを含む少なくとも２つのエネルギー状態を有し、状態｜ｅ＞′ｊのエネルギーは状態｜２＞′ｊのエネルギーよりも高く、｜２＞′ｊ−｜ｅ＞′ｊ遷移は前記共振器モードに共鳴している、複数の物理系Ｙｊを用意し、
前記複数の物理系Ｘｉのうちの２つの物理系Ｘｓ（ｓはＮ以下の自然数）及び物理系Ｘｔ（ｔはｓと異なるＮ以下の自然数）の状態を操作するために、｜１＞ｓ−｜ｅ＞ｓ遷移に共鳴する第１のレーザ光と、｜１＞ｔ−｜ｅ＞ｔ遷移に共鳴する第２のレーザ光と、前記複数の物理系Ｙｊに関して｜２＞′ｊ−｜ｅ＞′ｊ遷移における状態｜２＞′ｊにポピュレーションを集めるための第３のレーザ光と、を前記共振器に照射することを特徴とする量子計算方法。
【請求項６】
前記第３の光ビームの照射によって前記複数の物理系Ｙｊに関して状態｜２＞′ｊにポピュレーションを集めた後に、前記第１の光ビーム及び前記第２の光ビームを前記共振器に照射することを備えることを特徴とする請求項５に記載の量子計算方法。
【請求項７】
前記複数の物理系Ｘｉを｜１＞ｉ−｜ｅ＞ｉ遷移の周波数が小さい順にＸ１、Ｘ２、…、ＸＮとし、｜１＞１−｜ｅ＞１遷移と｜１＞２−｜ｅ＞２遷移との間の周波数差をΔ、｜１＞１−｜ｅ＞１遷移と｜１＞ｋ−｜ｅ＞ｋ遷移との間の周波数差をΔｋとした場合に、前記複数の物理系ＸｉはΩ０＜Δ＜Ｎ２ｇ／Ｎ及びΩ０＜Δｊ＜Ｎ２ｇ／Ｎを満たし、ここで、ｇは前記共振器モードと各物理系Ｘｉとの結合定数であり、Ω０は前記第１のレーザ光及び前記第２のレーザ光のラビ周波数の最大値を定めるパラメータであることを特徴とする請求項５に記載の量子計算方法。
【請求項８】
Ｎ２≧Ｎ２を満たすことを特徴とする請求項５に記載の量子計算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、共振器と物理系の結合を利用した量子計算機に関する。
【背景技術】
【０００２】
近年、量子力学的な重ね合わせの状態を用いて演算を行う量子計算機が盛んに研究されている。量子計算機の１つとして、周波数領域で量子ビットを区別する周波数領域量子計算に基づく量子計算機が知られている。周波数領域量子計算においては、位置で量子ビットを区別しないため、操作するつもりのない量子ビットにも離調付きで操作光が作用し、それによって不要な相互作用が生じる。この不要な相互作用はゲートエラーの原因になり得る。量子ビットに用いる遷移の周波数差が非常に大きければこの不要な相互作用の影響は小さくなる。しかしながら、有限の周波数領域に分布する遷移を利用する場合には、より多くの量子ビットを実装するために、なるべく周波数差の小さい遷移を利用することになる。周波数差の小さい遷移を利用する場合においても、不要な相互作用の影響を抑制しながら量子ゲートを実行できる周波数領域量子計算に基づく量子計算機が求められている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００１−２０９０８３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、量子ゲートを高効率で実行することができる量子計算機及び量子計算方法を提供することである。
【課題を解決するための手段】
【０００５】
一実施形態に係る量子計算機は、共振器の中に配置された複数の物理系Ｘｉ（ｉ＝１，２，…,Ｎ、Ｎは２以上の整数）、前記共振器の中に配置され、前記複数の物理系Ｘｉと異なる物理系Ｙｊ（ｊ＝１，２，…,Ｎ２、Ｎ２は１以上の整数）、及び光源部を備える。各物理系Ｘｉは、状態｜０＞ｉ、｜１＞ｉ、｜２＞ｉ、｜ｅ＞ｉを含む少なくとも４つの状態を有し、状態｜ｅ＞ｉのエネルギーは量子ビットに用いる状態｜０＞ｉ、｜１＞ｉそれぞれのエネルギー及びゲート操作の補助に用いる状態｜２＞ｉのエネルギーよりも高く、｜２＞ｉ−｜ｅ＞ｉ遷移は前記共振器の共振器モードに共鳴している。各物理系Ｙｊは、状態｜２＞′ｊ、｜ｅ＞′ｊを含む少なくとも２つのエネルギー状態を有し、状態｜ｅ＞′ｊのエネルギーは状態｜２＞′ｊのエネルギーよりも高く、｜２＞′ｊ−｜ｅ＞′ｊ遷移は前記共振器モードに共鳴している。光源部は、前記複数の物理系Ｘｉのうちの２つの物理系Ｘｓ（ｓはＮ以下の自然数）及び物理系Ｘｔ（ｔはｓと異なるＮ以下の自然数）の状態を操作するために、｜１＞ｓ−｜ｅ＞ｓ遷移に共鳴する第１のレーザ光と、｜１＞ｔ−｜ｅ＞ｔ遷移に共鳴する第２のレーザ光と、前記複数の物理系Ｙｊに関して｜２＞′ｊ−｜ｅ＞′ｊ遷移における状態｜２＞′ｊにポピュレーションを集めるための第３のレーザ光と、を前記共振器に照射する。
【図面の簡単な説明】
【０００６】
【図１】周波数領域量子計算を説明する図。
【図２】周波数領域量子計算に用いる物理系における不要な相互作用を示す図。
【図３】実施形態に係る量子計算機を示す図。
【図４】実施形態に係る周波数領域量子計算に用いる物理系を示す図。
【図５】実施形態において利用するＹ２ＳｉＯ５結晶中のＰｒ３＋イオンが有するエネルギー状態の一部を示す図。
【図６】実施形態に係る２つの物理系を操作する方法を示すフローチャート。
【図７】実施形態に係る、共振器を介したアディアバティックパッセージに用いるパルス波形を示す図。
【図８】他の実施形態に係る量子計算機を示す図。
【図９】さらに他の実施形態に係る量子計算機を示す図。
【発明を実施するための形態】
【０００７】
以下、図面を参照しながら種々の実施形態を説明する。以下の実施形態では、同一の構成要素に同一の参照符号を付して、重複する説明を適宜省略する。
【０００８】
まず、周波数領域量子計算とその量子ゲートに存在する「不要な相互作用」について説明する。また、その不要な相互作用の影響が特に強くなる条件である「共鳴条件」について説明する。さらに、その共鳴条件を制御し不要な相互作用の影響を受けない量子ゲートの操作方法及び構成を説明する。
【０００９】
［周波数領域量子計算の方法］
周波数領域量子計算では、共振器中に配置された複数の物理系であって、各物理系の１つの遷移が共通の共振器モード（固有モード）に共鳴しており且つそれ以外の遷移の周波数が物理系ごとに異なる複数の物理系それぞれを量子ビットとして利用する。物理系としては、例えば、イオンや原子などを利用することができる。周波数領域量子計算は、各物理系の遷移周波数に共鳴するレーザ光を照射することによりその物理系を選択的に操作することによって演算を行うものである。
【００１０】
物理系としてＮ個の四準位系Ｘｉ（ｉ＝１，２，３，…，Ｎ）を用いる場合について説明する。ここで、Ｎは２以上の整数である。各四準位系Ｘｉは４つのエネルギー状態を有する。これらの４つのエネルギー状態をエネルギーの低い順に｜０＞ｉ、｜１＞ｉ、｜２＞ｉ、｜ｅ＞ｉと表示する。各状態（各ケットベクトル）に付した添え字ｉは、その状態を有する四準位系Ｘｉを識別するためのものである。以下では、添え字ｉは省略することもある。状態｜０＞ｉと｜１＞ｉを量子ビットに用い、状態｜２＞ｉをゲート操作の補助に用いる。励起状態である状態｜ｅ＞ｉは、状態｜０＞ｉ、｜１＞ｉ、｜２＞ｉに比べてエネルギーが高い状態である。｜２＞ｉ−｜ｅ＞ｉ遷移（状態｜２＞ｉと状態｜ｅ＞ｉとの間の遷移）は共通の共振器モードに共鳴する遷移である。｜１＞ｉ−｜ｅ＞ｉ遷移の周波数は四準位系Ｘｉごとに異なる。
【００１１】
Ｎ＝３の場合における周波数領域量子計算に使用される物理系を図１に示す。図１に示されるように、物理系Ｘ１、Ｘ２、Ｘ３は共通の共振器モードに結合している。具体的には、｜２＞１−｜ｅ＞１遷移、｜２＞２−｜ｅ＞２遷移、及び｜２＞３−｜ｅ＞３遷移は、共通の共振器モードに共鳴している。すなわち、｜２＞１−｜ｅ＞１遷移、｜２＞２−｜ｅ＞２遷移、及び｜２＞３−｜ｅ＞３遷移の周波数は、共通の共振器モードの共振周波数に等しい。｜１＞１−｜ｅ＞１遷移、｜１＞２−｜ｅ＞２遷移、及び｜１＞３−｜ｅ＞３遷移の周波数は、互いに異なっている。
【００１２】
これらの四準位系Ｘｉのうちの２つの四準位系、例えば、四準位系Ｘ１及びＸ２の状態を操作する場合には、｜１＞１−｜ｅ＞１遷移に共鳴する操作光Ｌ１及び｜１＞２−｜ｅ＞２遷移に共鳴する操作光Ｌ２を四準位系Ｘｉ全体に照射する。操作光Ｌ１及びＬ２の照射によって、理想的には四準位系Ｘ１及びＸ２の状態を選択的に操作することができる。量子ゲートに用いられる状態の操作方法としては、共振器を介したアディアバティックパッセージの方法が知られている。この方法では、例えば、四準位系Ｘ１及びＸ２の状態を初期状態｜１＞１及び｜２＞２から状態｜２＞１及び｜１＞２に操作する場合、操作光Ｌ１及びＬ２のラビ周波数Ω１及びΩ２がτ１＞τ２の条件下で下記数式（１）に従うように操作光Ｌ１及びＬ２の強度を制御する。
【００１３】
【数１】
2015135377.tif 000003
【００１４】
［不要な相互作用］
周波数領域量子計算においては、厳密には下記のような「不要な相互作用」が生じる。不要な相互作用は、位置で四準位系Ｘｉを区別しないために、操作するつもりのない四準位系Ｘｉにも離調付きで操作光が作用することで生じる。具体的には、操作光Ｌ１が四準位系Ｘ２、Ｘｊ（ｊ＝３，４，…，Ｎ）に与える作用と操作光Ｌ２が四準位系Ｘ１、Ｘｊ（ｊ＝３，４，…，Ｎ）に与える作用がある。Ｎ＝３の場合における不要な相互作用を含む物理系を図２に示す。図２において、状態｜１＞ｉと状態｜ｅ＞ｉとを結ぶ実線の両方向矢印が必要な相互作用を表し、破線の両方向矢印が不要な相互作用を表し、状態｜２＞ｉと状態｜ｅ＞ｉとを結ぶ太い実線の両方向矢印が共振器モードを表している。
【００１５】
四準位系Ｘｉの｜１＞ｉ−｜ｅ＞ｉ遷移間の周波数差が非常に大きければ不要な相互作用の影響は非常に小さくなる。しかしながら、一般的には、遷移の周波数は有限の周波数領域の中で分布しているため、多くの量子ビットを利用するためにはなるべく周波数差の小さい遷移を利用することになる。このため、周波数差が小さい場合でも量子ゲートを高効率で実行できることが望ましい。
【００１６】
［共鳴条件］
不要な相互作用の性質について詳しく説明する。図２に示される物理系のような不要な相互作用を含む周波数領域量子計算のための物理系を記述するハミルトニアンは下記数式（２）のように表される。
【００１７】
【数２】
2015135377.tif 000004
【００１８】
ここで、σ（ｉ）ａｂは四準位系Ｘｉの状態｜ｂ＞ｉを状態｜ａ＞ｉに遷移させる演算子であり、ａ及びａ†はそれぞれ共振器モードの消滅演算子及び生成演算子である。ｇは共振器モードと物理系との結合定数であり、γは遷移の緩和速度であり、κは共振器の減衰定数である。Ｈ．ｃ．はエルミート共役である。
【００１９】
ハミルトニアンの各項について説明する。第１項は各イオンの各状態のエネルギー項と各イオンのエネルギー緩和項を含む。第２項は共振器モードのエネルギー項と共振器緩和項を含む。以下では、理想的な強結合系を仮定し、γ＝κ＝０の場合について説明する。第３項は、相互作用項であり、共振器モードと各四準位系Ｘｉの｜２＞ｉ−｜ｅ＞ｉ遷移の相互作用、操作光Ｌ１と各四準位系Ｘｉの｜１＞ｉ−｜ｅ＞ｉ遷移の相互作用、操作光Ｌ２と各四準位系Ｘｉの｜１＞ｉ−｜ｅ＞ｉ遷移の相互作用を含む。
【００２０】
数式（２）のハミルトニアンの相互作用ハミルトニアンは、下記数式（３）に示されるＨ０を用いることで、下記数式（４）のように表される。
【００２１】
【数３】
2015135377.tif 000005
【００２２】
【数４】
2015135377.tif 000006
【００２３】
Δは｜ｅ＞１−｜１＞１遷移と｜ｅ＞２−｜１＞２遷移との間の周波数差を表し、Δｊは｜ｅ＞１−｜１＞１遷移と｜ｅ＞ｊ−｜１＞ｊ遷移との間の周波数差を表す。
【００２４】
数式（４）の相互作用ハミルトニアンＨ′は、下記数式（５）のように必要な相互作用Ｈ１と不要な相互作用Ｖに分けられる。
【００２５】
【数５】
2015135377.tif 000007
【００２６】
数式（５）を用いて、共振器を介したアディアバティックパッセージにおいて、不要な相互作用によって生じるエラー確率を摂動計算によって計算する。初期状態｜ψ（０）＞をＨ１の固有状態の１つであるダークステート｜ψ０（０）＞であるとする。アディアバティックパッセージのエラー確率は、時刻ｔで別の固有状態｜ψｎ（ｔ）＞（ｎ≠０）へ遷移する確率である。この操作を高効率で実行するためには、Ω１，２＜＜ｇとなるような操作光Ｌ１及びＬ２を用いるのが一般的である。そのような場合は、ＶはＨ１より小さいので、下記数式（６）に示すように、Ｖについての摂動計算により時間発展を計算することができる。
【００２７】
【数６】
2015135377.tif 000008
【００２８】
数式（６）中のＥｎはＨ１の固有状態｜ψｎ＞に対応する固有値である。ＥｎはΩ１，２によってのみ時間変化するのでΩ１，２＜＜ｇの場合、その変化量は絶対値に比べて十分小さい。従って、係数Ｃｎ（１）、Ｃｎ（２）の中の各指数関数の指数がゼロのとき、摂動計算が有効でなくなるほどエラー確率は増大する。Ｖ（ｔ′）に振動項があることを考慮すると、１次の係数Ｃｎ（１）から共振器を介したアディアバティックパッセージのエラー確率が増大する条件として下記数式（７）に示す共鳴条件を得ることができる。
【００２９】
【数７】
2015135377.tif 000009
【００３０】
同様に、２次の係数Ｃｎ（２）から下記数式（８）に示す共鳴条件を得る。
【００３１】
【数８】
2015135377.tif 000010
【００３２】
２次の係数Ｃｎ（２）から得られた共鳴条件は１次の係数Ｃｎ（１）から得られた共鳴条件と比べると一般には寄与は小さい。しかしながら、複数の条件が同時に満たされる場合は、その寄与は大きくなる。数式（８）の共鳴条件において複数の条件を同時に満たすための条件として、下記数式（９）に示す条件が得られる。
【００３３】
【数９】
2015135377.tif 000011
【００３４】
共鳴条件の解析解はＨ１の固有値Ｅｎを求めることで得られる。Ω１，２＜＜ｇの場合を考えＨ１のΩ１，２の項をゼロとすると、固有値Ｅｎは良く知られた真空ラビ分裂のアナロジーで求めることができる。状態｜２＞ｉ又は｜ｅ＞ｉにポピュレーションがある四準位系Ｘｉの数をＮ２、状態｜ｅ＞ｉにポピュレーションがある四準位系Ｘｉの数をｎｅ、共振器モードの光子数をｎｃと表記する。Ｈ１の固有値は全励起子数Ｎｅ（Ｎｅ＝ｎｅ＋ｎｃ）とｎｅの最大値ｎｅ｜ｍａｘを指定することによって分類することができる。Ｈ１の固有値の一部を下記数式（１０）に示す。
【００３５】
【数１０】
2015135377.tif 000012
【００３６】
数式（７）、数式（９）及び数式（１０）によって不要な相互作用が増大する条件である共鳴条件を解析的に求めることができる。これは、共振器を介したアディアバティックパッセージを高効率で実行するために、つまり周波数領域量子計算において量子ゲートを高効率で実行するために、避けなければならない条件である。
【００３７】
［共鳴条件の制御］
数式（１０）に示されるように、状態｜２＞ｉ又は｜ｅ＞ｉにポピュレーションがある四準位系Ｘｉの数Ｎ２によってＨ１の固有値が変化する。これは、共振器に結合する遷移に確率振幅が存在する四準位系Ｘｉの数によって共鳴条件を制御できることを意味する。このような制御を共振器に結合する遷移の追加による制御と呼ぶ。
【００３８】
例えば３量子ビットを利用する量子計算機において、共振器を介したアディアバティックパッセージを実行する場合には、四準位系Ｘ１、Ｘ２、Ｘ３の初期状態を｜１＞１、｜２＞２、｜１＞３とすると、状態｜２＞１、｜１＞２、｜１＞３になるように操作する。その操作の中でＮｅ＝１の場合はＮ２＝２であり、Ｎｅ＝２の場合はＮ２＝３である。従って、数式（１０）からＨ１の固有値は下記数式（１１）に示すように得られる。数式（１１）に示される固有値を用いて数式（７）及び数式（９）から共鳴条件が得られる。
【００３９】
【数１１】
2015135377.tif 000013
【００４０】
これに対して、例えば共振器に結合する遷移に遷移確率がある物理系をさらに３個追加した場合には、下記数式（１２）に示すように、Ｈ１の固有値が得られる。
【００４１】
【数１２】
2015135377.tif 000014
【００４２】
この場合、｜Δ｜,｜Δｊ｜＜ｇの領域に共鳴条件を避けられる領域が大幅に増える。従って、確率的に四準位系Ｘｉの周波数分布が与えられたときに、量子ゲートを高効率で実行できる確率が高くなる。さらに、より多くの量子ビットを利用することが可能になる。
【００４３】
このように共振器に結合する遷移の追加により共鳴条件を制御することによって、不要な相互作用を抑制することが可能となり、その結果、高効率な周波数領域量子計算を実行することができる。ここで示した例は一例であり、実際に与えられた物理系の周波数分布に対して、共鳴条件を避けられるように共振器に結合する遷移を追加する数を変えることで、より最適な条件へと制御することができる。
【００４４】
次に、図３から図９を参照しながら実施形態を説明する。
【００４５】
（第１の実施形態）
図３は、第１の実施形態に係る量子計算機３００を概略的に示している。本実施形態では、四準位系Ｘｉとみなせる具体的な物理系として、Ｙ２ＳｉＯ５結晶中にドープされたＰｒ３＋イオンを使用する。具体的には、量子計算機３００は、図３に示されるように、Ｐｒ３＋イオンをドープしたＹ２ＳｉＯ５結晶（Ｐｒ：ＹＳＯ）の側面に誘電体多層膜ミラーを配置し共振器に加工した試料３１５を備える。試料３１５は、クライオスタット３１６内に設置され、低温に（例えば４Ｋに）保たれる。
【００４６】
ここでは、図４に示される例のように、６つのＰｒ３＋イオンを周波数領域量子計算に用い、これらのＰｒ３＋イオンのうちの３つのＰｒ３＋イオンを量子ビットとして利用する場合を想定する。量子ビットとして利用する３つのＰｒ３＋イオンをイオンＸ１、Ｘ２、Ｘ３と表示し、残り３つのＰｒ３＋イオンをイオンＹ１、Ｙ２、Ｙ３と表示する。イオンＹ１、Ｙ２、Ｙ３は共鳴条件の制御のために使用される。
【００４７】
具体的には、イオンＸ１、Ｘ２、Ｘ３の各々は、エネルギーの低い順に状態｜０＞、｜１＞、｜２＞、｜ｅ＞を有する。状態｜０＞、｜１＞は量子ビットに用い、状態｜２＞はゲート操作の補助に用いる。｜２＞−｜ｅ＞遷移は共振器の共通の共振器モードに共鳴している。イオンＹ１、Ｙ２、Ｙ３の各々は、エネルギーの低い順に状態｜０＞′、｜１＞′、｜２＞′、｜ｅ＞′を有する。図４では、イオンＹ１、Ｙ２、Ｙ３の状態｜０＞′、｜１＞′は省略されている。イオンＹ１、Ｙ２、Ｙ３の｜２＞′−｜ｅ＞′遷移は共振器の共通の共振器モードに共鳴している。ここで、イオンＹ１、Ｙ２、Ｙ３に関連する状態にダッシュ符号を付しているが、これは、共鳴状態の制御のために使用される物理系の種類が量子ビットに利用される物理系の種類と異なっていてもよいことを表す。イオンＹ１、Ｙ２、Ｙ３は、イオンＸ１、Ｘ２、Ｘ３と同じ種類のイオン（この例ではＰｒ３＋イオン）、他の種類のイオン、又はその組み合わせであってもよい。イオンＹ１、Ｙ２、Ｙ３がイオンＸ１、Ｘ２、Ｘ３と異なる種類のイオンである場合、イオンＹ１、Ｙ２、Ｙ３は状態｜０＞′、｜１＞′を有していなくてもよい。
【００４８】
図５は、Ｙ２ＳｉＯ５結晶中のＰｒ３＋イオンのエネルギー状態の一部を示している。図５には、基底状態３Ｈ４の超微細構造状態のうちの３つ（｜−１／２＞、｜−３／２＞、｜−５／２＞）と励起状態１Ｄ２の超微細構造状態のうちの３つ（｜＋１／２＞、｜＋３／２＞、｜＋５／２＞）が示されている。励起状態１Ｄ２と基底状態３Ｈ４との間の遷移周波数に対応する波長は約６０６ｎｍである。基底状態３Ｈ４の超微細構造状態｜−１／２＞、｜−３／２＞、｜−５／２＞が図４に示される状態｜０＞、｜１＞、｜２＞にそれぞれ対応し、励起状態１Ｄ２の超微細構造状態｜＋５／２＞が図４に示される状態｜ｅ＞に対応する。この場合、｜０＞−｜ｅ＞遷移、｜１＞−｜ｅ＞遷移、｜２＞−｜ｅ＞遷移は、光学的に遷移可能であり、｜０＞−｜１＞遷移、｜０＞−｜２＞遷移、｜１＞−｜２＞遷移は、光学的に禁制な遷移である。
【００４９】
図３に示される量子計算機３００では、アルゴンイオンレーザ３０１で励起したリング色素レーザ３０２を光源として使用する。リング色素レーザ３０２から発せられたレーザ光は、２つのビームスプリッタ３０３及び３０４によって３つのレーザ光に分割され、これら３つのレーザ光は、音響光学素子３０６、３０７及び３０８にそれぞれ導かれる。具体的には、ビームスプリッタ３０３は、リング色素レーザ３０２からのレーザ光を２つのレーザ光に分割し、これら２つのレーザ光の一方は音響光学素子３０６に入射し、他方はビームスプリッタ３０４に入射する。ビームスプリッタ３０４は、入射レーザ光を２つのレーザ光に分割し、これら２つのレーザ光の一方は音響光学素子３０７に入射し、他方はミラー３０５によって反射され音響光学素子３０８に入射する。
【００５０】
音響光学素子３０６、３０７及び３０８は、制御装置３０９によって生成された信号に従って、入射レーザ光を変調して変調レーザ光３５１、３５２及び３５３をそれぞれ生成する。変調レーザ光３５１は、ミラー３１０及び３１１並びにレンズ３１４によって試料３１５に導かれる。変調レーザ光３５２は、レンズ３１４によって試料３１５に導かれる。変調レーザ光３５３は、ミラー３１２及び３１３並びにレンズ３１４によって試料３１５に導かれる。本実施形態では、アルゴンイオンレーザ３０１、リング色素レーザ３０２、ビームスプリッタ３０３及び３０４、ミラー３０５、音響光学素子３０６〜３０８、ミラー３１０〜３１３、並びに、レンズ３１４によって光源部３２０が形成されている。
【００５１】
図３、図６及び図７を参照して、３つの量子ビットのうちのイオンＸ１及びＸ２の量子ビットを操作する方法を具体的に説明する。まず、光源部３２０は、イオンＹ１、Ｙ２、Ｙ３に関して｜２＞−｜ｅ＞遷移における状態｜２＞にポピュレーションを集めるための変調レーザ光３５３を試料３１５に照射する（図６のステップＳ６０１）。具体的には、音響光学素子３０８は、変調レーザ光３５３が｜０＞−｜ｅ＞遷移又は｜１＞−｜ｅ＞遷移に共鳴するように入射レーザ光を変調する。
【００５２】
続いて、光源部３２０は、変調レーザ光３５３を照射した状態で、イオンＸ１及びＸ２を操作するための変調レーザ光３５１及び３５２を試料３１５に同時に照射する（図６のステップＳ６０２）。ここで、同時に照射するとは、変調レーザ光３５１及び３５２の照射時間が少なくとも一部重なっていることを意味する。具体的には、音響光学素子３０６は、変調レーザ光３５１のラビ周波数Ω１が数式（１）に従って変化するように入射レーザ光を変調し、音響光学素子３０７は、変調レーザ光３５２のラビ周波数Ω２が数式（１）に従って変化するように入射レーザ光を変調する。イオンＸ１及びＸ２の状態を初期状態｜１＞１｜２＞２から状態｜２＞１｜１＞２に操作する場合、τ１＞τ２に設定される。イオンＸ１及びＸ２の状態を初期状態｜１＞１｜２＞２から状態｜２＞１｜１＞２に操作する場合におけるラビ周波数Ω１及びΩ２の時間変化を図７に示す。図７において、横軸は時間を表し、縦軸はラビ周波数を表す。一点鎖線で示される波形がラビ周波数Ω１を示し、実線で示される波形がラビ周波数Ω２を示す。一例では、レーザ光３５１及び３５２のラビ周波数Ω１及びΩ２のパラメータはΩ０＝１ｋＨｚ、τ１＝６４．１ｍｓ、τ２＝５５．９ｍｓ、σ＝２０ｍｓに設定される。
【００５３】
このようにして変調レーザ光３５１、３５２及び３５３を試料３１５に照射することで、共鳴条件を避けながらイオンＸ３の状態（例えば状態｜１＞３）を変化させずに、イオンＸ１及びＸ２の状態を初期状態｜１＞１｜２＞２から状態｜２＞１｜１＞２に操作することができる。
【００５４】
以上のように、第１の実施形態に係る量子計算機では、量子ビットに利用する物理系と異なる、共通の共振器モードに結合する遷移を含む物理系を利用することによって、不要な相互作用を抑制しながら、飛躍的に高効率な量子ゲートを実行することができる。
【００５５】
（第２の実施形態）
第１の実施形態では、状態操作に用いるレーザ光を生成するために、１つの光源を使用している。第２の実施形態では、状態操作に用いるレーザ光ごとに光源を用意する。
図８は、第２の実施形態に係る量子計算機８００を概略的に示している。量子計算機８００は、図８に示されるように、半導体レーザ８０１、８０２及び８０３を備える。半導体レーザ８０１、８０２及び８０３から発生したレーザ光は、音響光学素子３０６、３０７及び３０８にそれぞれ向けられる。音響光学素子３０６、３０７及び３０８などの動作は第１の実施形態と同様であるので、その説明を省略する。本実施形態では、半導体レーザ８０１〜８０３、音響光学素子３０６〜３０８、ミラー３１０〜３１３、及びレンズ３１４によって光源部８２０が形成されている。
第２の実施形態によれば、第１の実施形態と同様に、高効率な量子ゲートを実行することができる。
【００５６】
（第３の実施形態）
図９は、第３の実施形態に係る量子計算機９００を概略的に示している。図８に示される量子計算機９００は、図３に示される第１の実施形態の量子計算機３００と同様の構成を備える。第３の実施形態では、第１の実施形態と異なる部分及び動作について説明する。
本実施形態では、量子ビットとして利用されるイオンＸ１、Ｘ２、Ｘ３と共鳴条件の制御に利用されるイオンＹ１、Ｙ２、Ｙ３とが異なる空間領域に位置する。例えば、イオンＸ１、Ｘ２、Ｘ３が試料３１５の上部に配置されたＰｒ３＋イオンから選択され、イオンＹ１、Ｙ２、Ｙ３が試料３１５の下部に配置されたＰｒ３＋イオンから選択される。この場合、図９に示されるように、光源部９２０は、変調レーザ光３５１及び３５２を試料３１５の上部に照射し、変調レーザ光３５３を試料３１５の下部に照射する。
第３の実施形態によれば、第１の実施形態と同様に、高効率な量子ゲートを実行することができる。
【００５７】
（第４の実施形態）
第４の実施形態では、共通の共振器モードに結合するＰｒ３＋イオンの中から量子ビットに利用する複数のＰｒ３＋イオンを選択する方法例を説明する。なお、第４の実施形態に係る量子計算機は、第１の実施形態に係る量子計算機３００（図３）と同様の構成を有し、同様の方法で量子ゲートを実行するものとし、重複する説明を省略する。
【００５８】
試料３１５の中に含まれ共通の共振器モードに結合する遷移を含むＭ個のＰｒ３＋イオンを、｜１＞−｜ｅ＞遷移の周波数が小さい順に、イオン１、イオン２、…、イオンＭと表示する。ここで、Ｍは２以上の整数である。さらに、｜１＞１−｜ｅ＞１遷移と｜１＞２−｜ｅ＞２遷移との間の周波数差をΔとし、｜１＞１−｜ｅ＞１遷移と｜１＞ｊ−｜ｅ＞ｊ遷移（ｊは３以上Ｍ以下の整数）との間の周波数差をΔｊとする。この場合において、Ω０＜Δ＜Ｎ２ｇ／Ｎ及びイオンΩ０＜Δｊ＜Ｎ２ｇ／Ｎを満たすイオンを量子ビットに利用する。さらに、共鳴状態の制御のために使用するイオンは、周波数差Δｊが大きい順に所定数だけ選択される。
【００５９】
共鳴状態の制御のために使用する物理系（この例ではＰｒ３＋イオン）の数Ｎ２は、１以上の任意の値であってもよいが、量子ビットの数Ｎを二乗した値以上にすることが望ましい。これは、数式（１０）に示されるようにＨ１の固有値がＮ２の平方根に実質的に依存するためである。例えば、３つの量子ビットを利用する場合、共鳴状態の制御のための物理系として、共通の共振器モードに結合する遷移を含む９つのＰｒ３＋イオンを使用する。これら９つのＰｒ３＋イオン（イオンＹ１〜Ｙ９と記載する）は、上述した方法のように周波数差Δｉが大きい順にイオン１〜Ｍから選択される。この場合、第１の実施形態においてイオンＹ１、Ｙ２、Ｙ３に対して行った操作と同様の操作をイオンＹ１〜Ｙ９に対して行う。具体的には、イオンＹ１〜Ｙ９に関して状態｜２＞にポピュレーションが集まるように変調したレーザ光を試料３１５に照射する。
第４の実施形態によれば、第１の実施形態と同様に、高効率な量子ゲートを実行することができる。
【００６０】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００６１】
３００…量子計算機、３０１…アルゴンイオンレーザ、３０２…リング色素レーザ、３０３，３０４…ビームスプリッタ、３０５…ミラー、３０６〜３０８…音響光学素子、３０９…制御装置、３１０〜３１３…ミラー、３１４…レンズ、３１５…試料、３１６…クライオスタット、３２０…光源部、８００…量子計算機、８０１〜８０３…半導体レーザ、８２０…光源部、９００…量子計算機、９２０…光源部。
【図１】
2015135377.tif 000015
【図２】
2015135377.tif 000016
【図３】
2015135377.tif 000017
【図４】
2015135377.tif 000018
【図５】
2015135377.tif 000019
【図６】
2015135377.tif 000020
【図７】
2015135377.tif 000021
【図８】
2015135377.tif 000022
【図９】
2015135377.tif 000023