(57)【特許請求の範囲】
【請求項１】
分岐パラメータを変化させる量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせをそれぞれ実現し、損失を伴わない非散逸的な結合によって相互に結合する複数の量子非線形発振器と、
前記量子非線形発振器の前記分岐パラメータを個別に制御する制御装置と、
前記量子非線形発振器の出力を測定する測定器と、
を具備する量子計算装置。
【請求項２】
前記量子非線形発振器間を接続し、前記の損失を伴わない非散逸的な結合を実現する複数の結合器をさらに具備する請求項１に記載の量子計算装置。
【請求項３】
前記制御装置は、さらに、計算の初期段階において系の時間発展を規定するハミルトニアンの最小または最大の固有値に対応する固有状態が真空状態となるように、前記量子非線形発振器の固有周波数を制御する請求項１または請求項２に記載の量子計算装置。
【請求項４】
前記測定器は、前記量子非線形発振器の出力の正負を測定するパリティ測定器を含む請求項１から３のいずれか１項に記載の量子計算装置。
【請求項５】
前記量子非線形発振器がパラメトリック増幅とカー効果を有するカー・パラメトリック発振器である請求項１から４のいずれか１項に記載の量子計算装置。
【請求項６】
前記カー・パラメトリック発振器が、ジョセフソン接合を有する超伝導回路で実装される請求項５に記載の量子計算装置。
【請求項７】
前記量子非線形発振器が、ジョセフソン接合を有する超伝導回路で実装される請求項１から６のいずれか１項に記載の量子計算装置。
【請求項８】
前記非散逸的な結合が、キャパシタを介した結合であることを特徴とする請求項１から７のいずれか１項に記載の量子計算装置。
【請求項９】
前記非散逸的な結合が、前記量子非線形発振器の固有周波数から離調された共振周波数を持つ結合共振器を用いて実装される請求項１から７のいずれか１項に記載の量子計算装置。
【請求項１０】
分岐パラメータを変化させる量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせをそれぞれ実現する複数の量子非線形発振器を、損失を伴わない非散逸的な結合によって相互に結合し、
前記量子非線形発振器の前記分岐パラメータを制御し、
前記量子非線形発振器の出力を測定すること
を具備する量子計算方法。
【請求項１１】
計算の初期段階において系の時間発展を規定するハミルトニアンの最小または最大の固有値に対応する固有状態が真空状態となるように、前記量子非線形発振器の固有周波数を設定することをさらに具備する請求項１０に記載の量子計算方法。
【請求項１２】
前記量子非線形発振器の固有周波数の初期値を設定し、初期段階以降の固有周波数を制御する離調制御装置と、
計算の初期段階において系の時間発展を規定するハミルトニアンの最小または最大の固有値に対応する固有状態が真空状態となる、という初期条件を満たす範囲で固有周波数を異なる設定に制御して計算することを前記離調制御装置に繰り返えさせ、各計算結果で最良のものを選択して出力する解選択装置と、を具備する請求項１から９のいずれか１項に記載の量子計算装置。
【請求項１３】
前記離調制御装置が前記異なる設定では異なる初期値を選んでその後一定値とする場合を含む請求項１２に記載の量子計算装置。
【請求項１４】
前記離調制御装置は、前記固有周波数の初期値は前記初期条件を満たし、計算終了時にはパラメトリック励起に対する離調がゼロとなるように制御する場合を含む請求項１２または１３に記載の量子計算装置。
【請求項１５】
計算の初期段階において系の時間発展を規定するハミルトニアンの最小または最大の固有値に対応する固有状態が真空状態となる、という初期条件を満たす範囲で、前記量子非線形発振器の固有周波数を異なる設定に制御して計算することを繰り返し、各計算結果で最良のものを選択して出力することを特徴とする請求項１０または１１に記載の量子計算方法。
【請求項１６】
前記異なる設定では異なる初期値を選んでその後一定値とする場合を含むことを特徴とする請求項１５に記載の量子計算方法。
【請求項１７】
前記異なる設定が、初期値は前記初期条件を満たすように設定され、計算終了時にはパラメトリック励起に対する離調がゼロになるように変化させる場合を含むことを特徴とする請求項１５または１６に記載の量子計算方法。
【請求項１８】
量子非線形発振器の分岐後の２つの量子状態を量子ビットとして用い、
量子断熱変化によって量子ビットの初期状態を準備するための分岐パラメータの制御と、
量子断熱変化によって位相シフトゲートを実行するための量子非線形発振器への外部励起の制御と、
量子断熱変化によって回転ゲートを実行するための量子非線形発振器の固有周波数のパラメトリック励起に対する離調の制御と、
量子断熱変化によってエンタングリングゲートを実行するための量子非線形発振器間の結合係数の制御と、
を実行するゲート操作用制御装置をさらに具備する請求項１，２，４，５，６，７，８，９のいずれかに記載の量子計算装置。
【請求項１９】
分岐パラメータを変化させる量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせをそれぞれ実現する複数の量子非線形発振器を用い、
量子非線形発振器の分岐後の２つの量子状態を量子ビットとして用い、
量子非線形発振器への外部励起パルスの照射による量子断熱変化によって位相シフトゲートを実行し、
量子非線形発振器の固有周波数のパラメトリック励起に対する離調をパルス的に変化させる量子断熱変化によって回転ゲートを実行し、
量子非線形発振器間の結合係数をパルス的に変化させる量子断熱変化によってエンタングリングゲートを実行すること、
を具備する量子計算方法。
【請求項２０】
前記量子非線形発振器を初め真空状態に準備し、パラメトリック励起に対する離調を有限値からゼロへと変化させるとともに分岐パラメータの大きさをゼロから分岐点以上へ変化させる量子断熱変化によって、量子ビットの初期状態を準備することを具備する請求項１９に記載の量子計算方法。
【請求項２１】
分岐パラメータを変化させる量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせをそれぞれ実現する複数の量子非線形発振器と、
前記量子非線形発振器が生成する分岐後の２つの量子状態を量子ビットとして用い、前記量子非線形発振器へ外部励起パルスを照射して量子断熱変化によって位相シフトゲートを実行し、前記量子非線形発振器の固有周波数のパラメトリック励起に対する離調をパルス的に変化させる量子断熱変化によって回転ゲートを実行し、前記量子非線形発振器間の結合係数をパルス的に変化させる量子断熱変化によってエンタングリングゲートを実行する制御装置と、
前記量子非線形発振器の出力を測定して前記量子ビットを読み出す測定装置と、
を具備する量子計算装置。
【請求項２２】
前記制御装置は、前記量子非線形発振器を初め真空状態に準備し、パラメトリック励起に対する離調を有限値からゼロへと変化させるとともに分岐パラメータの大きさをゼロから分岐点以上へ変化させる量子断熱変化によって、量子ビットの初期状態を準備する請求項２１に記載の量子計算装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、量子計算装置及び方法に関し、例えば組み合わせ最適化問題を解くことに利用される。
【背景技術】
【０００２】
多くの組み合わせパターンの中から、目的に合った量（目的関数と呼ばれる）が最大、または、最小となるパターンを見つけ出す問題は、社会の様々な場面で現れる。このような問題を組み合わせ最適化問題という。組み合わせ最適化問題は、問題のサイズに対して解の候補の数が指数関数的に増えるため、多くの場合、現在の逐次計算を行うノイマン型コンピュータを用いて解くことが困難であることが知られている。
【０００３】
組み合わせ最適化問題の中で、ビット値を取る変数からなる２次関数を最大、または、最小にする問題をQUBO（quadratic unconstrained binary optimization）問題という。多くの組み合わせ最適化問題がQUBO問題に定式化できる。QUBO問題は、多くの場合解くことが困難であることが知られている。ビット値をスピンの古典的なモデルであるイジングスピン±１に対応させると、QUBO問題は統計力学のイジングモデルのエネルギー最小状態（基底状態）を求める問題（イジング問題と呼ぶ）に帰着される。
【０００４】
近年、このQUBO問題（イジング問題）を高速に解くことができるとされる計算装置が注目されている。例えば、レーザーネットワークやパラメトリック発振器ネットワークの発振現象を利用して、イジングモデルの基底状態を見つけ出す方式が提案されている。この発振現象を利用する方式は、それ以前のQUBO問題の解法であるシミュレーティッドアニーリング（熱揺らぎを利用）や量子アニーリング（量子揺らぎを利用）と異なる原理的に新しい方式であり、途中で局所最適解にトラップされにくいと考えられている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開第２０１２／１１８０６４号パンフレット
【特許文献２】国際公開第２０１５／００６４９４号パンフレット
【非特許文献】
【０００６】
【非特許文献１】arXiv:1503.00194（http://arxiv.org/abs/1503.00194）. "Holonomic quantum computing with cat-codes"; Albert et al.; 2015.
【非特許文献２】New. J. Phys. 16, 045014. "Dynamically protected cat-qubits: a new paradigm for universal quantum computation"; Mirrahimi et al.; 2014.
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところが、この発振現象を利用する方式において、実際に発振しきい値で得られる状態が近似解になる場合がある。その場合、真の最適解へ到達するには、シミュレーティッドアニーリングのように、何らかの雑音によって局所最適解を抜け出しながらより良い解を探す必要がある。この際、上記文献の原理では、原理的に損失が不可欠であり、損失は量子力学的な重ね合わせの状態を破壊するため、量子アニーリングのような量子力学的な効果は期待できない。その結果として、解の精度をあまり上げることができない。
【０００８】
本発明が解決しようとする課題は、上記事情を考慮してなされたものであって、発振現象と量子力学的な効果をともに利用する新しい原理によって、組み合わせ最適化問題を高速に解くことができる量子計算装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
実施形態によれば、量子計算装置は、複数の量子非線形発振器と、制御装置と、測定器と、を含む。複数の量子非線形発振器は、分岐パラメータを変化させる量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせをそれぞれ実現し、損失を伴わない非散逸的な結合によって相互に結合する。制御装置は、量子非線形発振器の分岐パラメータを個別に制御する。測定器は、量子非線形発振器の出力を測定する。
【図面の簡単な説明】
【００１０】
【図１】実施形態の量子計算装置を示す図。
【図２】実施形態の量子計算方法を示すフローチャート。
【図３】実施形態の量子計算装置の概念図。
【図４Ａ】量子モデルでの数値シミュレーションのエラー確率の結果を表す図。
【図４Ｂ】古典モデルでの数値シミュレーションのエラー確率の結果を表す図。
【図５Ａ】量子モデルでの数値シミュレーションの残留エネルギーの結果を表す図。
【図５Ｂ】古典モデルでの数値シミュレーションの残留エネルギーの結果を表す図。
【図６Ａ】図４Ａ及び図４Ｂと図５Ａ及び図５Ｂのシミュレーションにおいて、古典モデルでは１００％誤る問題のエネルギー図。
【図６Ｂ】図６Ａの問題に対する量子モデルでの各スピン状態の確率の時間発展を表す図。
【図６Ｃ】図６Ａの問題に対する古典モデルでの各スピン状態の確率の時間発展を表す図。
【図７】超伝導回路によって実装されたカー・パラメトリック発振器を表す図。
【図８】３つのスピンの場合の量子計算装置の実施例を表す図。
【図９】カー・パラメトリック発振器と読み出しラインの結合をオン・オフできる回路を表す図。
【図１０】結合共振器を用いて非散逸的な結合を実装する量子計算装置の実施例を表す図。
【図１１】結合共振器のパラメータを調整できる回路を表す図。
【図１２】カー・パラメトリック発振器の基底状態・第２励起状態間のエネルギーギャップの離調依存性を表す図。
【図１３】カー・パラメトリック発振器の分岐点におけるウィグナー関数の離調依存性を表す図。
【図１４】離調を３つの一定値に設定した場合のシミュレーション結果を示す図。
【図１５】図１６のシミュレーションにおける離調とポンプ振幅の設定を表す図。
【図１６】離調を図１５のように最終的にゼロにするように制御する場合と一定値とした場合（図１４）のシミュレーション結果を比較する図。
【図１７】離調を制御することで精度を改善する量子計算装置の実施例を表す図。
【図１８】ゲート型万能量子計算機に利用する際の初期化のシミュレーション結果を表す図。
【図１９】ゲート型万能量子計算機に利用する際の位相シフトゲートのシミュレーション結果を表す図。
【図２０】ゲート型万能量子計算機に利用する際の回転ゲートのシミュレーション結果を表す図。
【図２１】ゲート型万能量子計算機に利用する際のエンタングリングゲートのシミュレーション結果を表す図。
【図２２】量子計算装置を用いたゲート型万能量子計算機の概念図。
【図２３】ゲート型万能量子計算機に利用する際の量子計算装置の実施例を表す図。
【発明を実施するための形態】
【００１１】
以下、図面を参照しながら本発明の実施形態について詳細に説明する。
【００１２】
本実施形態の量子計算装置の基本構成要素は、図１に示すように、複数の量子非線形発振器１０１、量子非線形発振器１０１同士を結合する複数の結合器１０２、系のパラメータを制御する制御装置１０３、及び量子非線形発振器１０１を測定する測定器１０４を備えている。まずは、図１の構成図、図２の量子計算方法を示すフローチャートを参照して説明する。
【００１３】
量子非線形発振器１０１は、単独である場合、系のパラメータ（分岐パラメータと呼ぶ）を制御パラメータとする量子断熱変化によって、１つの量子状態を分岐させ区別可能な量子状態の重ね合わせを生成することが可能である。また、量子非線形発振器１０１の損失は、その影響が無視できるほど小さいとする。（理想的な場合、損失はないとする。）
結合器１０２はそれぞれ、複数の量子非線形発振器１０１のうちの２つを、損失なく非散逸的に結合する。
【００１４】
結合器１０２は、装置の実体はなく単に量子非線形発振器１０１同士が直接結合している場合も示すとする。
【００１５】
制御装置１０３は、全ての量子非線形発振器１０１の分岐パラメータを個別に制御する。
【００１６】
測定器１０４は、全ての量子非線形発振器１０１の出力を測定する。量子計算装置は、この測定結果によって計算結果を得ることができる。測定器１０４は、例えば電場振幅の位相のパリティを測定する。
【００１７】
量子計算装置及び方法についてより詳細に以下に記述する。
【００１８】
まず、単独の量子非線形発振器１０１について詳しく説明する。量子非線形発振器１０１は、分岐パラメータを制御パラメータとする量子断熱変化によって、区別可能な量子状態の重ね合わせへと変化する「量子力学的な分岐現象」を示す。イジングスピンの±１は、分岐後の区別可能な２つの量子状態に対応させられる。以下ではこのような発振器として、カー効果を有するパラメトリック発振器（以下、「カー・パラメトリック発振器」）を具体例に適宜用いて説明するが、量子力学的な分岐現象を示す量子非線形発振器であれば同様の効果が期待できる。
【００１９】
カー・パラメトリック発振器の時間発展は次のハミルトニアンＨによって与えられる。
【００２０】
【数１】
0006530326.tif 000002
【００２１】
ここで、
【００２２】
【数２】
0006530326.tif 000003
【００２３】
はそれぞれ発振器の消滅演算子と生成演算子を表す。また、Ｈｐａｒａはパラメトリック増幅、ＨＫｅｒｒはカー効果、Ｈｄｅｔは離調にそれぞれ対応するハミルトニアンであり、ｐはパラメトリック増幅のポンプ振幅、Ｋはカー係数、Δは離調である。（離調は、発振器の固有周波数をω、パラメトリック増幅のためのポンプ周波数をωｐとして、Δ＝ω−ωｐ／２で定義される。）特に断らない限り、以下ではｐ、Ｋ、Δを非負実数とする。
この本実施形態のカー・パラメトリック発振器は、特許文献２で用いられた通常のパラメトリック発振器とは異なり、２光子損失の代わりにカー効果を、１光子損失の代わりに離調を有しており、理想的には損失を有していない点に特長がある。
【００２４】
この量子非線形発振器１０１は、初期状態を真空状態（消滅演算子のゼロ固有値状態）として、パラメトリック増幅のポンプ振幅ｐ（カー・パラメトリック発振器の分岐パラメータ）をゼロから徐々に大きくしていくと、ｐ＝Δで分岐（発振）し始め、ｐ＞Δでは振幅ｘの期待値の大きさが等しく符号が逆符号の２つの発振状態の量子力学的な重ね合わせになる。ここで振幅ｘは以下である。
【００２５】
【数３】
0006530326.tif 000004
【００２６】
これが、前述の量子力学的な分岐現象である。符号が逆符号の２つの発振状態がイジングスピンの±１に対応させられる。ここで、分岐点（発振しきい値）が離調で決まる点に注意する。（これに対し、特許文献２で用いられた通常のパラメトリック発振器の発振しきい値は１光子損失で決まる。）
以下、上で説明した量子非線形発振器１０１を複数用いる本実施形態の量子計算装置について説明する。本実施形態の量子計算装置では、初期状態を真空とした量子断熱変化を利用して最適解を見つけ出す。そのため、初期のハミルトニアンは、真空を最小固有値の固有状態として持つ必要がある。（今はＫを正としているが、Ｋが負の場合はすべての符号を反転し、「正」を「負」に、「最小」を「最大」などと適宜読みかえる。）ここではこの条件を「量子断熱変化の初期条件」と呼ぶこととする。この条件を満たすには、例えば、離調を以下の条件を満たすように設定すればよい（例えば、制御装置１０３が設定する）。（初期設定段階：図２のＳ２０１）
離調の条件（Δｉはｉ番目の量子非線形発振器の離調）：
【００２７】
【数４】
0006530326.tif 000005
【００２８】
つまり、
【００２９】
【数５】
0006530326.tif 000006
【００３０】
で定義される行列Ｍが半正定値（つまり、すべての固有値が非負）となる。
【００３１】
初期値設定段階の次は、解きたい問題に合せて量子非線形発振器の間で、非散逸的な結合を設定する（結合段階：図２のＳ２０２）。本実施形態の量子計算装置は、図３の概念図に示すように、問題に応じて複数の量子非線形発振器３０１を適切に結合３０２したネットワークで構成される。以下では、エネルギー関数が下記のＥで与えられるイジングモデルを例に説明する。以下、Ｅを「イジングエネルギー」と呼ぶ。
【００３２】
【数６】
0006530326.tif 000007
【００３３】
ここで、ｓｉはｉ番目のイジングスピンの値を表し、±１を取る。また、結合定数｛Ｊｉ，ｊ｝はＪｉ，ｉ＝０、及び、Ｊｉ，ｊ＝Ｊｊ，ｉを満たす。
【００３４】
上記イジングモデルに対し、本実施形態の量子計算装置では、Ｎ個の量子非線形発振器を適切に結合し、以下の結合ハミルトニアンＨＩを実装する。
【００３５】
【数７】
0006530326.tif 000008
【００３６】
ξ０は結合エネルギーの大きさを定める正の定数である。
【００３７】
ここで、上記の結合が損失を全く伴っていない点に注意する。そこで、この結合を「非散逸的な結合」と呼ぶ。（これに対し、特許文献２で用いられた相互注入による結合方式はこのようなハミルトニアンでは記述できず、散逸的である。）また、上記の結合は標準的な線形の結合であり、物理的実装が容易である。
【００３８】
上記の初期条件が満たされたとき、ネットワークの分岐点（発振しきい値）は、結合エネルギーが小さい発振パターンほど低くなる。結合エネルギーの最小化は近似的にイジングエネルギーの最小化になるため、初期状態を真空状態として分岐パラメータをゼロから徐々に大きくして発振させると、イジングエネルギーが低い発振パターンの重ね合わせ状態が得られる（分岐段階：図２のＳ２０３）。このようにイジングエネルギーの低い状態が自動的に得られる点が発振現象を利用する方式のメリットであり、これにより、イジングエネルギーが高い局所最適解にトラップされることを避けられる。また分岐点付近では、イジングエネルギーの低い状態が１つだけ選ばれるのではなく、量子揺らぎによってイジングエネルギーが低い複数の状態の量子力学的な重ね合わせとなる。これにより、大域最適解が高い確率で見つかることが期待される。
【００３９】
さらに分岐パラメータを大きくしていくと、最終的にネットワークの状態は大域最適解へと収束していく。これは、断熱条件が満たされているという仮定の下で、量子断熱定理から保証される。初め、分岐パラメータはゼロであり、初期状態である真空状態は全系の基底状態である。分岐パラメータを十分ゆっくり変化させると、量子断熱定理により、全系の状態は最終的なハミルトニアンの基底状態へと変化する。ここで、分岐パラメータが十分大きい場合、非線形項（カー・パラメトリック発振器の場合、パラメトリック増幅とカー効果）が主要となり、発振振幅はどの発振器も同程度となり、位相（上記のｘの符号）だけが異なる状況になる。そして、各発振器からこの位相を測定して測定結果から計算結果を得る（測定段階：図２のＳ２０４）。全系のハミルトニアンの中で位相に依存するのは結合ハミルトニアンだけであり、振幅がすべて等しく位相だけが異なるという条件下では、これはイジングエネルギーに比例する。よって、最終状態である全系の基底状態はイジングエネルギーを最小にする状態、つまり、大域最適解に対応する状態となる。このように、分岐点付近では量子揺らぎによってイジングエネルギーが低い複数の状態の重ね合わせになり、最終的には量子断熱変化によって大域最適解へ収束することは、量子効果が有効に働いていると考えることができる。
【００４０】
以上のように、本実施形態の量子計算装置では、理想的には損失を全く伴わない量子非線形発振器を、やはり損失を伴わない非散逸的な結合でつないだネットワークによって、発振現象（分岐現象）と量子効果をともに有効に利用して組み合わせ最適化問題を解くことができる。
【００４１】
次に、この量子効果を確認するために、カー・パラメトリック発振器を用いた本実施形態の量子計算装置の量子モデルと古典モデルのシミュレーションを行い、比較したことについて図４Ａ及び図４Ｂと図５Ａ及び図５Ｂを参照して説明する。ここで、量子モデルでは上記のハミルトニアンを用いたシュレディンガー方程式を数値的に解き、古典モデルでは次の常微分方程式系を数値的に解いた。
【００４２】
【数８】
0006530326.tif 000009
【００４３】
この方程式は、下記の消滅演算子のハイゼンベルグ方程式から導出される。
【００４４】
【数９】
0006530326.tif 000010
【００４５】
そして、この消滅演算子を下記の式で表される古典近似をして得られる。
【００４６】
【数１０】
0006530326.tif 000011
【００４７】
ここで、このｘｉ，ｙｉは実変数である。
【００４８】
本シミュレーションでは４つのスピンのイジング問題を解いた。結合定数｛Ｊｉ，ｊ｝を−１〜１の一様乱数でランダムに設定し、１０００問の問題を作成し、それぞれについてエラー確率（大域最適解が得られない確率）と残留エネルギー（得られた解のイジングエネルギーと真のイジングエネルギー最小値の差）の期待値を見積もった。その結果をヒストグラムにまとめたものをそれぞれ図４Ａ及び図４Ｂと図５Ａ及び図５Ｂに示す。
【００４９】
図４Ｂを参照すると、古典モデルでも多くの問題で大域最適解が見つかることがわかる。これはイジングエネルギーが低いものほど早く発振する上記の機構が古典モデルでも成り立つからである。一方、図４Ｂからわかるように、古典モデルでは高い確率で誤る問題もかなりある。これらの問題で誤る原因は、局所最適解にトラップされてしまうためと思われる。それに対し、図４Ａを参照すると、量子モデルでは高い確率で誤る問題の数が劇的に減っている。これは、量子揺らぎによって局所最適解にトラップされるのを回避し、量子断熱定理によって大域最適解へ収束したためと考えられる。
【００５０】
以上を確認するため、古典モデルで１００％誤る問題について詳しく見る。図６Ａに、この問題のエネルギー図（energy landscape）を示し、図６Ｂ及び図６Ｃにそれぞれ量子モデル及び古典モデルの各スピン状態の確率の時間変化を示した。ここではスピンの１と−１をビットの０と１に対応させて表示する。イジングエネルギーＥＩｓｉｎｇが最小となる大域最適解は０１００と１０１１であり、エネルギー図の「最適解からの距離」はビット表示における大域最適解からのハミング距離を表す。０００１と１１１０は局所最適解になっている。古典モデルでは、図６Ｃを参照すると、まず局所最適解である０００１と１１１０の確率が１に上がり、その後、その隣の００１１と１１００との間で振動し、大域最適解の確率は最終的にほぼゼロになる。これに対し、量子モデルでは、図６Ｂを参照すると、初めにイジングエネルギーが比較的低いスピン状態の確率が高い重ね合わせの状態となり、最終的には大域最適解の確率が１へと収束している。この結果は、量子揺らぎの効果で局所最適解へのトラップを回避し、量子断熱定理のために大域最適解へ収束した、と解釈できる。
【００５１】
以上から、上述の本計算装置が期待通りに動作することが確認された。
【００５２】
以下、離調の調整による計算の高速化（高精度化）を提案する。（ここで、「高精度」とは、目的関数の最適値からのずれが小さいことを意味する。イジング問題の場合、残留エネルギーが小さいことに相当する。）
量子断熱変化は、基底状態と励起状態のエネルギーギャップの最小値が大きいほど高速化できる。（計算時間が同じであればより高精度になる。）図１２に１つのカー・パラメトリック発振器の基底状態・第２励起状態間のエネルギーギャップの離調依存性を示す。（カー・パラメトリック発振器では、基底状態と第１励起状態の間の遷移はハミルトニアンのパリティ反転対称性から禁止されているため、基底状態・第２励起状態間のエネルギーギャップを考えれば十分。）図１２は、離調が大きいほどエネルギーギャップの最小値が大きくなることを示している。図１２によれば、Δ＝０の場合は約１（ｐ／Ｋ）、Δ＝Ｋの場合は約３（ｐ／Ｋ）、Δ＝２Ｋの場合は約４．５（ｐ／Ｋ）、Δ＝５Ｋの場合は約７．７（ｐ／Ｋ）になる。そこで、ネットワークによって最適化問題を解く場合も、離調を付加することでエネルギーギャップを開き、高速化（高精度化）できることが期待される。
【００５３】
離調の付加による精度の改善は、量子揺らぎの観点からも期待される。図１３に、図１２の各離調に対する分岐点（ｐ＝Δ）におけるウィグナー関数を示す。（ウィグナー関数とは、位相空間（前記の古典モデルの２つの変数ｘとｙで張られる空間）上の擬確率密度関数として知られるもの。）図１３からわかるように、離調が大きいほど分岐点におけるｘ方向の量子揺らぎが大きくなる。よって、離調が大きいほど局所最適解にトラップされにくいと期待される。特に、離調がＫに比べて大きい場合に量子揺らぎの効果が大きくなる。
【００５４】
以上のことを確認するため、上記と同様の４スピンの１０００問を解くシミュレーションを行った。（ここでは量子モデルのみ扱う。）離調は一定ですべての発振器で同じとする。その離調の値として、次の３つの値を選んだ。
1.  Δ＝３ξ０
2.  Δ＝１０ξ０
3.  Δ＝２０ξ０
初めの値は量子断熱変化の初期条件（上式（４）の行列Ｍが半正定値になるという条件）を満たす中で比較的小さい値、他の２つは比較的大きい値である。（なお、本シミュレーションではξ０＝０．２Ｋとした。）また、ｐはゼロから５Ｋ＋Δまで時刻に比例するように増加させた。（最終的な平均光子数をおよそ５にそろえるために５Ｋ＋Δまでとしている。）計算時間（ｐの掃引時間）は１００／Ｋとした。
【００５５】
シミュレーションの結果を図１４に示す。図１４は、残留エネルギーの分布関数Ｐ（Ｅ）（残留エネルギーがＥ以下となる確率）を１０００問の残留エネルギーの結果から見積もり、それを１から引いたもの（つまり、１−Ｐ（Ｅ））を両対数で表している。同じＥに対してグラフが下にあるほど精度が高いことを意味する。図１４から、期待通り、離調を比較的大きく取ったほうが精度が高くなることが確認された。一方、Δ＝１０ξ０の場合のほうがΔ＝２０ξ０の場合よりも精度が高くなったことから、大きいほど良いというものでもないことがわかる。以上から、離調は、単に数式（４）の行列Ｍが正定値になるように取れば何でも良いというものではなく、付加的な離調をあえて付け、その値を変えながら何度か計算し、残留エネルギーが最小となる解を選択する、という方法が望ましい。
【００５６】
離調が大きすぎると、前述の「分岐パラメータが十分大きい場合、非線形項（カー・パラメトリック発振器の場合、パラメトリック増幅とカー効果）が主要となり、発振振幅はどの発振器も同程度となり、位相（上記のｘの符号）だけが異なる状況になる。」という仮定が成り立たなくなることに起因するモデル誤差が大きくなる。このモデル誤差は、離調を最終的にゼロに変化させることで消すことができると考えられる。（完全にゼロでなくてもＫ以下にすれば効果があるため、以下では「離調をゼロにする」というのはＫ以下にすることを意味することとする。）
以上のことを確認するために、初めは離調を大きく取ることでエネルギーギャップを開き、ｐがある程度大きくなった後、離調を最終的にゼロに変化させる場合のシミュレーションを行った。離調の初期値として、図１４のシミュレーションと比較するために、Δ＝１０ξ０とΔ＝２０ξ０の２つの場合を調べた。また、ｐの最終値は、図１３のシミュレーションと比較するために５Ｋとした。（つまり、最終的な平均光子数がおよそ５とした。）以上のように設定したΔとｐの時間変化を図１５に示す。
【００５７】
このシミュレーション結果と図１４のシミュレーション結果との比較を図１６に示す。離調を一定とする場合と異なり、図１５のように離調を制御する場合は、残留エネルギーが０．００１以下と非常に小さくなる確率が増えている。これは、上記のモデル誤差がなくなったためと考えられる。一方、図１５の制御ありでΔ（０）＝２０ξ０の場合、残留エネルギーが大きくなる確率が増えている。これは、図１５の制御の結果、断熱条件が悪くなったためと考えられる。それに比べ、図１５の制御ありでΔ（０）＝１０ξ０の場合は、離調一定の場合よりも全体として高精度になっている。これは、断熱条件を悪化させることなく、モデル誤差を消すことができたため、と考えられる。以上から、離調を一定とする場合と最終的にゼロへ変化させる場合の両方を、値を変えながら何度か計算し、残留エネルギーが最小となる解を選択する、という方法が望ましい。
【００５８】
なお、これまで相互作用項のみのイジングモデルを扱ったが、局所外部磁場項も含んだＥの場合も扱うことができる。ここでＥは以下である。
【００５９】
【数１１】
0006530326.tif 000012
【００６０】
以下、カー・パラメトリック発振器を例にこの場合を説明する。
【００６１】
各カー・パラメトリック発振器に周波数ωｐ／２の電磁波を入力し、次のハミルトニアンを実装する：
【００６２】
【数１２】
0006530326.tif 000013
【００６３】
ここで、α（ｔ）は発振時の振幅の大きさに相当し、カー・パラメトリック発振器の場合は以下で表せる。
【００６４】
【数１３】
0006530326.tif 000014
【００６５】
α（ｔ）は初めゼロであり、ｐ（ｔ）とともにゆっくりと変化する。ｐ（ｔ）が十分大きくなったとき、発振振幅はどの発振器も同程度となり、位相だけが異なる状況になる。全系のハミルトニアンの中で位相に依存するのは結合ハミルトニアンと上記のＨＭａｇだけであり、振幅がすべて等しく位相だけが異なるという条件下では、これはイジングエネルギーに比例する。こうして、前述と同様、量子断熱定理により、最終的にイジングエネルギーを最小にする状態が得られる。
なお、ここでは局所外部磁場項のために入力される電磁波の振幅を実数としたが、あえて上記の振幅に虚部を追加することでエネルギーギャップが開き、精度が改善する場合がある。（虚部は最終的にゼロにする。）
【００６６】
本実施形態の量子計算装置によれば、発振現象と量子力学的な効果をともに利用する新しい原理によって、組み合わせ最適化問題を高速に解くことができる。
本量子計算装置は、各量子非線形発振器の２つの発振状態を量子ビットと見なして、量子ゲート操作に基づく万能量子計算機として用いることもできる。万能量子計算機は、現在の計算機では計算量が多くて解けない問題、例えば、素因数分解を高速に解くことができる。以下、本量子計算装置を万能量子計算機として用いる方法を説明する。
【００６７】
量子ゲート操作に基づく万能量子計算機に必要な処理は、量子ビットの読み出し、量子ビットの初期化、位相シフトゲート、回転ゲート、エンタングリングゲートの５つである。ここで、各ゲートは次のように定義される。（ゲートの定義には任意性があるが、ここでは本装置に合うものを選んだ。）
位相シフトゲート： 
【００６８】
【数１４】
0006530326.tif 000015
【００６９】
回転ゲート： 
【００７０】
【数１５】
0006530326.tif 000016
【００７１】
エンタングリングゲート： 
【００７２】
【数１６】
0006530326.tif 000017
【００７３】
以下、各操作の方法を順に説明するが、エンタングリングゲートの場合を除き、発振器間の結合はすべて切れた状態にあるとする。
【００７４】
量子ビットの読み出しは、前記のイジングスピンの読み出しと同じである。
【００７５】
量子ビットの初期化は、ポンプ振幅をゼロから分岐点（発振しきい値）以上の値ｐ０へ徐々に大きくしていく量子断熱変化によって、各発振器の状態を真空状態から２つの発振状態の重ね合わせ状態
【００７６】
【数１７】
0006530326.tif 000018
【００７７】
（規格化因子は略）へ変化させることで達成される。ここで、
【００７８】
【数１８】
0006530326.tif 000019
【００７９】
は振幅が
【００８０】
【数１９】
0006530326.tif 000020
【００８１】
のコヒーレント状態であり、
【００８２】
【数２０】
0006530326.tif 000021
【００８３】
を量子ビットの｜０＞と見なし、
【００８４】
【数２１】
0006530326.tif 000022
【００８５】
を量子ビットの｜１＞と見なす。（ｐ０は十分大きいとして、これらは直交する。）この初期化で準備される状態を量子ビットで表すと、｜０＞＋｜１＞（規格化因子は略）となる。これは多くの量子アルゴリズム（例えば、ショアの素因数分解アルゴリズム）の初期状態と一致する。なお、この初期化において、各発振器の離調は最終的にゼロであると仮定しているが、前述のエネルギーギャップの考察から、初めは離調を比較的大きな値（Ｋと同程度以上）とし、徐々に減少させ、最終的にゼロにすることで、離調を常にゼロにする場合に比べ、初期化に必要な量子断熱変化を高速化、または、高精度化できる。このことを数値シミュレーションで確認した。初期化に要する時間Ｔを２０／Ｋとし、ポンプ振幅は
【００８６】
【数２２】
0006530326.tif 000023
【００８７】
離調は
【００８８】
【数２３】
0006530326.tif 000024
【００８９】
と設定し、各Δ０に対して最終的に得られた状態と理想的な初期状態
【００９０】
【数２４】
0006530326.tif 000025
【００９１】
（規格化因子は略）の間の忠実度（Fidelity）Ｆを計算した。（２つの状態ベクトルの間のＦはその内積の２乗で定義される。）その結果を図１８に示す。期待した通り、Δ０がＫと同程度以上で高い忠実度が得られることがわかった。
【００９２】
次に、位相シフトゲートの実行方法を説明する。ポンプ振幅はｐ０のまま、離調はゼロのまま、外部励起パルスＥｉｎ（ｔ）によって以下のハミルトニアンを実現する。
【００９３】
【数２５】
0006530326.tif 000026
【００９４】
なお、このパルスは量子断熱変化が起こるように十分ゆっくりと変化するエンベロープを持つものとする。Ｅｉｎ（ｔ）が常にＫに比べて十分小さいとき、各量子ビットのエネルギーは状態に応じて次のようにシフトする。
【００９５】
【数２６】
0006530326.tif 000027
【００９６】
これによって、次の角φの位相シフトゲートが実現される（Ｔはゲート時間）。
【００９７】
【数２７】
0006530326.tif 000028
【００９８】
以上を数値シミュレーションで確認した。前述の初期化のシミュレーションと同じ設定で初期化した後、外部励起パルスＥｉｎ（ｔ）を
【００９９】
【数２８】
0006530326.tif 000029
【０１００】
と変化させた（上の結果から、これで角φの位相シフトゲートが実行される）。ゲート時間Ｔは５／Ｋとした。最終的に得られた状態と理想的な状態
【０１０１】
【数２９】
0006530326.tif 000030
【０１０２】
（規格化因子は略）の間の忠実度（Fidelity）Ｆを計算した。図１９にその結果を示す。−π〜πのφに対して非常に高い忠実度が得られることがわかる。
【０１０３】
次に、回転ゲートの実行方法を説明する。ポンプ振幅はｐ０のまま、離調をパルス的に変化させる。なお、このパルスは量子断熱変化が起こるように十分ゆっくりと変化するエンベロープを持つものとする。離調がある程度大きくなると、離調がゼロのとき縮退していた光子数パリティが偶と奇の状態が分裂し、そのエネルギー差に依存した相対位相θが生じる。これによって、離調のパルス的変化の後、次のように量子ビットの状態が変化する。
【０１０４】
【数３０】
0006530326.tif 000031
【０１０５】
これは、全体にかかる位相因子を除き、角θの回転ゲートである。全体にかかる位相因子は無視できるため、これで回転ゲートが達成される。以上を数値シミュレーションで確認した。本シミュレーションでは、まず、前述のシミュレーションと同じ設定で初期化とφ＝π／２の位相シフトゲートを行い、本シミュレーションの初期状態
【０１０６】
【数３１】
0006530326.tif 000032
【０１０７】
を準備する。その後、離調を
【０１０８】
【数３２】
0006530326.tif 000033
【０１０９】
とパルス的に変化させた（Ｔはゲート時間）。本シミュレーションではＴ=１０／Ｋとした。各Δ０に対して、最終的に得られた状態と角θの回転ゲートの理想的な出力状態の間の忠実度を計算し、忠実度が最大となるθおよび対応する忠実度の最大値を求めた。図２０にその結果を示す。−π〜０のθが得られ、忠実度も十分高い値が得られることがわかった。
【０１１０】
最後に、エンタングリングゲートの実行方法を説明する。ポンプ振幅はｐ０のまま、離調はゼロのまま、エンタングリングゲートを実行したい２つの発振器間の結合係数ｇをパルス的に変化させ、以下のハミルトニアンを実現する。
【０１１１】
【数３３】
0006530326.tif 000034
【０１１２】
ここで、
【０１１３】
【数３４】
0006530326.tif 000035
【０１１４】
は各量子ビットの消滅演算子である。なお、このパルスは量子断熱変化が起こるように十分ゆっくりと変化するエンベロープを持つものとする。ｇ（ｔ）が常にＫに比べて十分小さいとき、２つの量子ビットのエネルギーは状態に応じて次のようにシフトする。
【０１１５】
【数３５】
0006530326.tif 000036
【０１１６】
これによって、次の角Θのエンタングリングゲートが実現される（Ｔはゲート時間）。
【０１１７】
【数３６】
0006530326.tif 000037
【０１１８】
以上を数値シミュレーションで確認した。前述の初期化のシミュレーションと同じ設定で２個の量子ビットを初期化した後、結合係数ｇ（ｔ）を
【０１１９】
【数３７】
0006530326.tif 000038
【０１２０】
と変化させた（上の結果から、これで角Θのエンタングリングゲートが実行される）。ゲート時間Ｔは５／Ｋとした。最終的に得られた状態と理想的な状態
【０１２１】
【数３８】
0006530326.tif 000039
【０１２２】
（規格化因子は略）の間の忠実度（Fidelity）Ｆを計算した。図２１にその結果を示す。０〜πのΘに対して非常に高い忠実度が得られることがわかる。
【０１２３】
ここで、以上の３種類の量子ゲート操作はすべて、量子断熱変化に伴う力学的位相（エネルギー差に依存する位相シフト）を制御することで実現されており、いわゆる幾何学的位相（ベリー位相やホロノミーなど）とは異なる点に注意する。（ここでの量子断熱変化は、各発振器のほぼ縮退した２つの基底状態で張られるヒルベルト空間内に各発振器の状態が閉じ込められている過程を指す。）また、回転ゲートでは、光子数パリティの偶奇の状態の縮退を解く際に、ポンプ振幅を十分小さくするのではなく、離調を十分大きくする点に注意する。（ポンプ振幅を小さくするよりも、離調を大きくする方が、エネルギーギャップを大きく保てるため、より高速な、または、より高精度な量子断熱変化が可能である。）
以上の処理を図２２の概念図にまとめた。量子ゲート操作に必要なパラメータ制御はゲート操作用制御装置２２０１が行う。
【０１２４】
以上のように、量子非線形発振器のネットワークを用いる本量子計算装置は、組み合わせ最適化問題を高速に解くことができるだけでなく、量子ゲート操作に基づく万能量子計算も実行することができる。
【実施例】
【０１２５】
以下、本発明の実施例について説明する。なお、以下では３つのスピンの場合について説明するが、他の場合も同様にして実施できる。また、カー・パラメトリック発振器を具体例に説明するが、量子力学的な分岐現象を示す量子非線形発振器であれば、同様の実施が可能である。
（実施例１）
まず、カー・パラメトリック発振器は、図７に示した超伝導回路によって実装される（Z. R. Lin et al., Nat. Commun. 5, 4480 (2014)）。超伝導マイクロ波共振器７０１、及びdc SQUID７０２からなる。なお、光カー効果と光パラメトリック増幅を利用することも考えられるが、ここでは小さな損失と大きな非線形性が両立できる超伝導回路を用いることとする。
【０１２６】
この回路は、超伝導マイクロ波共振器７０１をdc SQUID７０２で終端する構造になっており、制御装置１０３が外部電流を流しdc SQUID７０２内の磁束を外部電流で変調することで、発振器の固有周波数（離調）、及び、パラメトリック増幅のポンプ振幅を調整できる。これにより、量子断熱変化の初期条件を満たすための離調の調整、および、前述の計算原理を実行するためのパラメトリック増幅率の制御ができ、これらが制御装置１０３によって実行される。また、dc SQUID７０２のジョセフソン接合による非線形性のため、カー効果が得られる。
【０１２７】
上記以外にも、ジョセフソン接合の非線形性を利用した超伝導回路により、量子力学的な分岐現象を示す量子非線形発振器は様々な方法で実装できる。
【０１２８】
カー・パラメトリック発振器間の非散逸的な結合は、キャパシタを介した結合で実装される。図８に３つのカー・パラメトリック発振器８０１のネットワークの概念図を示す。
【０１２９】
制御装置８０３が分岐パラメータをゼロから分岐点を超えるように徐々に変化させた後、図８に示したように、結果の読み出しは、各カー・パラメトリック発振器８０１に読み出しラインをキャパシタ８０２を介して結合し、読み出しラインへの出力をアンプで増幅後、ヘテロダイン検出によって位相を測定することで行われる。得られた位相は、ポンプとの相対位相で決まる１８０°ずつの２つの範囲に対応して±１（パリティ）に変換される。これが計算結果となる。（図８では、以上の処理をパリティ測定器８０４が行う。）
なお、カー・パラメトリック発振器８０１と読み出しラインの結合は、計算中はオフにするのが望ましい。これは、図９に示したように、読み出しラインをdc SQUID９０１で終端することにより実現できる（M. Wallquis et al., Phys. Rev. B 74, 224506 (2006)）。
【０１３０】
局所外部磁場項も含むイジングモデルを扱う場合は、各カー・パラメトリック発振器に電磁波を入力できるポートを作り、そこから電磁波を入力する。
（実施例２）
（実施例１）では、非散逸的な結合を、カー・パラメトリック発振器８０１をキャパシタ８０２を介して直接結合することで実装された。しかし、直接結合する場合、発振器の配置の自由度が小さい上、結合定数を変更することが困難である。そこで、本実施例では、カー・パラメトリック発振器８０１間を結合共振器１００２で結合する（図１０）。これにより、配置の自由度が広がる上、結合共振器１００２のパラメータ調整によって、解きたい問題に応じた結合定数の変更が容易になる。以下ではまず、結合共振器１００２によって所望の結合が実装できることを説明する。
【０１３１】
ｉ番目の発振器とｊ番目の発振器（ｉ＜ｊ）を結合する結合共振器（ｉ，ｊ）の消滅演算子をｂｉ，ｊ、固有周波数をωｉ，ｊ、ｉ番目の発振器、及び、ｊ番目の発振器との結合強度をそれぞれｇｉ，ｊ、及び、ｇｊ，ｉと表す。このとき、結合共振器に対応するハミルトニアンは次のようになる：
【０１３２】
【数３９】
0006530326.tif 000040
【０１３３】
ここで、離調Δｉ，ｊは、Δｉ，ｊ＝ωｉ，ｊ−ωｐ／２で定義される。（便宜上、ｉ＞ｊに対しても離調をΔｊ，ｉ＝ωｉ，ｊ−ωｐ／２と定義する。）
離調Δｉ，ｊの絶対値が他のパラメータ（カー係数Ｋ、パラメトリック増幅のポンプ振幅ｐ、各発振器の離調Δｉ、結合強度ξ０Ｊｉ，ｊ、局所外部磁場項のためのポンプ振幅ξ０α）の絶対値に比べ十分大きいとき、結合共振器１００２の自由度を断熱消去することができ、その結果得られる有効ハミルトニアンは次のようになる。
【０１３４】
【数４０】
0006530326.tif 000041
【０１３５】
よって、下記の式を満たすようにｇｉ，ｊとΔｉ，ｊを選べば、上記の所望のハミルトニアンが実現できる。
【０１３６】
【数４１】
0006530326.tif 000042
【０１３７】
ここで、Δｉ，ｊをすべて負とすると、上記の定義から明らかにＨ’Ｉは非負演算子であり、量子断熱変化の初期条件を満たすための上記の行列Ｍの半正定値条件が自動的に満たされる。なお、行列Ｍの半正定値条件が満たされる範囲で各発振器に離調を追加してもかまわない。
【０１３８】
結合共振器１００２を変形することにより、量子非線形発振器１０１の配置の自由度を広げられる。また、発振器と結合共振器の結合強度ｇｉ，ｊや結合共振器の離調Δｉ，ｊを調整することで、解きたい問題に応じて結合定数｛Ｊｉ，ｊ｝を変更することができる。
【０１３９】
カー・パラメトリック発振器８０１と結合共振器１１０１の結合強度ｇｉ，ｊや結合共振器１１０１の離調Δｉ，ｊの調整は、前記の読み出しラインの場合と同様、結合共振器１１０１をdc SQUID１１０２で終端することにより実現できる。その様子を図１１に示す。これによって、ｇｉ，ｊやΔｉ，ｊの絶対値のみならず符号も調整できる。こうして、解きたい問題に応じて結合定数｛Ｊｉ，ｊ｝を変更することができる。
（実施例３）
離調を制御することによって精度を向上させる実施例について説明する。本実施例では、（実施例１）または（実施例２）の構成に加え、離調を制御する制御装置（１０３、８０３）へ様々な離調設定を送り、測定器（１０４、８０４）から計算結果を受け取り、異なる離調設定で得られた計算結果の中から最良のものを選択して出力する解選択装置を備える。本実施例では、非線形発振器１０１、８０１としてdcSQUIDを有する超伝導マイクロ波共振器を用いるため、離調の制御は外部電流によるdcSQUID内の磁束の制御で実現される。
【０１４０】
本実施例の構成を図１７に示す。解選択装置１７０１から離調設定が制御装置１０３、８０３へ送られ、制御装置１０３、８０３が受け取った離調設定に合わせて各発振器１０１、８０１の離調を制御する。制御として、以下の２つの場合を含むものとする。
１．  量子断熱変化の初期条件を満たす一定値とする。
２．  初期値は量子断熱変化の初期条件を満たす値とし、最終的にはゼロにする。
【０１４１】
上記１の場合なら離調の一定値を、上記２の場合なら離調の初期値をいろいろ変えて計算を繰り返す。解選択装置１７０１はパリティ測定器１０４、８０４から各計算結果を受け取り、それぞれの計算結果に対する残留エネルギーを計算し、残留エネルギーが最小となるものを選び出し、それを出力する。
（実施例４）
量子ゲート操作に基づく万能量子計算機の実施例について説明する。本実施例の構成を図２３に示す。ここでは簡単のため２個の量子ビットだけを示したが、３個以上の場合も同様である。結合係数を制御するために、実施例２の結合共振器による発振器間結合を用いる。結合係数の制御は、結合共振器のdc SQUID内の磁束を外部電流で制御することによって可能である。（数式（４１）からわかるように、結合共振器を用いると余計な離調が生じるが、それがポンプ振幅ｐ０に比べて十分に小さければ、その影響は無視できる。無視できない場合は、各発振器７０１、７０２の離調を調整することで消すことができる。）
本実施例では、量子ゲート操作を実行するためのゲート操作用制御装置２２０１を備える。ゲート操作用制御装置２２０１は、位相シフトゲートを実行する際は外部励起用電磁波を、回転ゲートを実行する際は各発振器のdc SQUID７０２の磁束を制御する外部電流によってポンプ振幅と離調を、エンタングリングゲートを実行する際は結合共振器が有するdc SQUIDの磁束を制御する外部電流によって発振器間の結合係数を、それぞれ制御する。
【０１４２】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【０１４３】
１０１、３０１…量子非線形発振器、１０２…結合器、１０３、８０３…制御装置、１０４…測定器、３０２…結合、７０１…超伝導マイクロ波共振器、７０２、９０１、１１０２…dc SQUID、８０１…カー・パラメトリック発振器、８０２…キャパシタ、８０４…パリティ測定器、１００２、１１０１…結合共振器、１７０１…解選択装置、２２０１…ゲート操作用制御装置。
【図１】
0006530326.tif 000043
【図２】
0006530326.tif 000044
【図３】
0006530326.tif 000045
【図４Ａ】
0006530326.tif 000046
【図４Ｂ】
0006530326.tif 000047
【図５Ａ】
0006530326.tif 000048
【図５Ｂ】
0006530326.tif 000049
【図６Ａ】
0006530326.tif 000050
【図６Ｂ】
0006530326.tif 000051
【図６Ｃ】
0006530326.tif 000052
【図７】
0006530326.tif 000053
【図８】
0006530326.tif 000054
【図９】
0006530326.tif 000055
【図１０】
0006530326.tif 000056
【図１１】
0006530326.tif 000057
【図１２】
0006530326.tif 000058
【図１３】
0006530326.tif 000059
【図１４】
0006530326.tif 000060
【図１５】
0006530326.tif 000061
【図１６】
0006530326.tif 000062
【図１７】
0006530326.tif 000063
【図１８】
0006530326.tif 000064
【図１９】
0006530326.tif 000065
【図２０】
0006530326.tif 000066
【図２１】
0006530326.tif 000067
【図２２】
0006530326.tif 000068
【図２３】
0006530326.tif 000069