(57)【要約】  本発明では、パラメトリック発振器１が、パラメトリック発振を用いて、同一の発振周波数を有する複数の擬似的スピンパルスＳＰｉを発振し、相互作用実装部５が、複数の擬似的スピンパルスＳＰｉの発振位相φｉ（暫定）の暫定的な測定結果を用いて、各擬似的スピンパルスＳＰｉが関わる相互作用（σｉに対する比例係数λｉ＋ΣＪｉｊσｊ＋ΣＫｉｊｋσｊσｋ）の大きさ及び符号をフィードバック実装し、擬似的スピン測定部６が、複数の擬似的スピンパルスＳＰｉの発振位相φｉ（定常）の最終的な測定結果に基づいて、複数の擬似的スピンパルスＳＰｉの擬似的なスピンσｉを測定する。
2015156126.tif 000002
【特許請求の範囲】
【請求項１】
イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振器と、
前記複数の擬似的スピンパルスを周回伝搬させるリング共振器と、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定部と、
ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定部が暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算部と、
前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算部が暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装部と、
前記暫定的スピン測定部、前記相互作用計算部及び前記相互作用実装部により構成されるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定部と、
を備えることを特徴とするイジングモデルの量子計算装置。
【請求項２】
前記相互作用計算部は、前記ある擬似的スピンパルスが関わる前記イジングモデルの３体以上の結合係数及び前記暫定的スピン測定部が暫定的に測定した前記他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる３体以上の相互作用を暫定的に計算し、
前記相互作用実装部は、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算部が暫定的に計算した前記ある擬似的スピンパルスが関わる３体以上の相互作用の大きさ及び符号を暫定的に実装する
ことを特徴とする請求項１に記載のイジングモデルの量子計算装置。
【請求項３】
前記パラメトリック発振器は、前記複数の擬似的スピンパルスと同一の発振周波数を有し１対１のペアを組む複数の局部発振パルスをパラメトリック発振させ、
前記リング共振器は、前記複数の局部発振パルスを周回伝搬させ、
前記複数の局部発振パルスの位相について、正相及び逆相の両方を含む状態から、正相及び逆相の一方を含む状態へと、同相化させるパルス同相化部、をさらに備え、
前記暫定的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、１対１のペアを組む前記複数の局部発振パルスの一部を用いてホモダイン検波を行い、
前記相互作用実装部は、前記ある擬似的スピンパルスに対して、１対１のペアを組みその一部の振幅及び位相を制御された局部発振パルスを注入し、
前記擬似的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、１対１のペアを組む前記複数の局部発振パルスの一部を用いてホモダイン検波を行う
ことを特徴とする請求項１又は２に記載のイジングモデルの量子計算装置。
【請求項４】
角周波数ωを有する局部発振パルスを発生させるパルス発生器と、
前記角周波数ωを有する局部発振パルスを用いて、角周波数２ωを有するパルスを発生させる第二高調波発生器と、をさらに備え、
前記パラメトリック発振器は、前記角周波数２ωを有するパルスを用いて、前記複数の擬似的スピンパルスをパラメトリック発振させ、
前記暫定的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、前記角周波数ωを有する局部発振パルスを用いてホモダイン検波を行い、
前記相互作用実装部は、前記ある擬似的スピンパルスに対して、振幅及び位相を制御された前記角周波数ωを有する局部発振パルスを注入し、
前記擬似的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、前記角周波数ωを有する局部発振パルスを用いてホモダイン検波を行う
ことを特徴とする請求項１又は２に記載のイジングモデルの量子計算装置。
【請求項５】
前記相互作用実装部は、前記ある擬似的スピンパルスに対して注入される光の振幅を、計算処理の初期段階ほど大きく制御し、計算処理の終期段階ほど小さく制御する
ことを特徴とする請求項１から４のいずれかに記載のイジングモデルの量子計算装置。
【請求項６】
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの振幅を暫定的に測定する暫定的振幅測定部、をさらに備え、
前記パラメトリック発振器は、前記暫定的振幅測定部が測定した前記複数の擬似的スピンパルスの振幅に基づいて、前記複数の擬似的スピンパルスの振幅が等しくなるように、パラメトリック発振に用いるポンプパルスの振幅をフィードバック制御する
ことを特徴とする請求項１から５のいずれかに記載のイジングモデルの量子計算装置。
【請求項７】
前記イジングモデルの複数のスピンのうち第ｎ群（ｎは１以上の整数）のスピンに擬似的に対応する第ｎ群の擬似的スピンパルスが前記リング共振器を周回伝搬する、請求項１から６のいずれかに記載の第ｎ番のイジングモデルの量子計算装置、を並列に備え、
前記第ｎ番のイジングモデルの量子計算装置が備える前記暫定的スピン測定部が暫定的に測定した前記第ｎ群の擬似的スピンパルスの擬似的なスピンの情報を、並列に備わるイジングモデルの量子計算装置の間で共有させる暫定的スピン共有部、をさらに備える
ことを特徴とするイジングモデルの量子並列計算装置。
【請求項８】
イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振ステップと、
前記複数の擬似的スピンパルスがリング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定ステップが暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算ステップが暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスをパラメトリック増幅させるパラメトリック増幅ステップと、
前記暫定的スピン測定ステップ、前記相互作用実装ステップ及び前記パラメトリック増幅ステップをこの順序で備えるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定ステップと、
を備えることを特徴とするイジングモデルの量子計算方法。
【請求項９】
イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振ステップと、
前記複数の擬似的スピンパルスがリング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定ステップが暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスをパラメトリック増幅させるパラメトリック増幅ステップと、
前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算ステップが暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装ステップと、
前記暫定的スピン測定ステップ、前記パラメトリック増幅ステップ及び前記相互作用実装ステップをこの順序で備えるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定ステップと、
を備えることを特徴とするイジングモデルの量子計算方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、イジングモデルを容易に解くことにより、イジングモデルにマッピングされるＮＰ完全問題などを容易に解くことができる量子計算装置を提供する。
【背景技術】
【０００２】
イジングモデルは、元来は磁性材料のモデルとして研究されてきたが、最近はＮＰ完全問題などからマッピングされるモデルとして注目されている。しかし、イジングモデルは、サイト数が大きいときには、解くことが非常に困難になる。そこで、イジングモデルを実装する量子アニールマシンや量子断熱マシンが提案されている。【０００３】
量子アニールマシンでは、イジング相互作用及びゼーマンエネルギーを物理的に実装してから、系を十分に冷却して基底状態を実現して、基底状態を観測することにより、イジングモデルを解いている。しかし、サイト数が大きいときには、系が冷却の過程で準安定状態にトラップされ、また準安定状態の数はサイト数に対して指数関数的に増大するため、系が準安定状態から基底状態になかなか緩和されないという問題があった。【０００４】
量子断熱マシンでは、横磁場ゼーマンエネルギーを物理的に実装してから、系を十分に冷却して横磁場ゼーマンエネルギーのみの基底状態を実現する。そして、横磁場ゼーマンエネルギーを徐々に下げ、またイジング相互作用を徐々に物理的に実装していき、イジング相互作用及び縦磁場ゼーマンエネルギーを含む系の基底状態を実現して、その基底状態を観測することにより、イジングモデルを解いている。しかし、サイトの数が大きいときには、横磁場ゼーマンエネルギーを徐々に下げ、またイジング相互作用を徐々に物理的に実装する速度はサイト数に対して指数関数的に遅くする必要があるという問題があった。【０００５】
ＮＰ完全問題などをイジングモデルにマッピングし、そのイジングモデルを物理的なスピン系で実装するときには、物理的に近くに位置するサイト間のイジング相互作用は大きく、物理的に遠くに位置するサイト間のイジング相互作用は小さいという自然法則が問題となる。ＮＰ完全問題をマッピングした人工的なイジングモデルでは、物理的に近くに位置するサイト間のイジング相互作用が小さいことがあり、物理的に遠くに位置するサイト間のイジング相互作用が大きいことがありえるからである。この自然なスピン系へのマッピングの難しさも、ＮＰ完全問題などを容易に解くことを困難にしていた。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第５３５４２３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上記の問題を解決するための第１の従来技術（特許文献１を参照）及び第２の従来技術について説明する。ＮＰ完全問題は、磁性体のイジングモデルに置き換え可能であり、磁性体のイジングモデルは、レーザーのネットワークに置き換え可能である。【０００８】
ここで、磁性体のイジングモデルでは、相互作用する原子ペアにおいて、スピン配列のエネルギーが最低となるように、スピンの方向は、逆方向（反強磁性の相互作用の場合）又は同方向（強磁性の相互作用の場合）を指向しようとする。【０００９】
一方で、レーザーのネットワークでは、相互作用するレーザーペアにおいて、発振モードの閥値利得が最低となるように、発振の偏光（第１の従来技術の場合）若しくは位相（第２の従来技術の場合）は、逆回転若しくは逆位相（反強磁性の相互作用の場合）又は同回転若しくは同位相（強磁性の相互作用の場合）を指向しようとする。【００１０】
つまり、１つのレーザーペアからなるシステムでは、発振モードの閥値利得が最低となるように、発振の偏光又は位相を最適化することができる。そして、多くのレーザーペアからなるシステムでは、「ある」レーザーペアで発振の偏光又は位相を最適化しようとすれば、「他の」レーザーペアで発振の偏光又は位相を最適化できないところ、レーザーのネットワークの「全体」として発振の偏光又は位相の「妥協点」を探索することになる。【００１１】
ただし、レーザーのネットワーク全体で発振の偏光又は位相を最適化するときには、各々のレーザーペアで別個の発振モードを立ち上げるのではなく、レーザーのネットワーク全体で１つの発振モードを立ち上げるように、各レーザー間で同期を図る必要がある。【００１２】
このように、第１の従来技術及び第２の従来技術では、各レーザーでポンピング電流を漸増制御し、レーザーのネットワーク全体で閾値利得が最低となる１つの発振モードを立ち上げ、各レーザーの発振の偏光又は位相を測定し、各原子のスピンの方向を測定する。よって、量子アニールマシンにおける準安定状態へのトラップの問題及び量子断熱マシンにおけるイジング相互作用の実装速度の問題を解決することができる。【００１３】
そして、第１の従来技術及び第２の従来技術では、図１及び図２を用いてそれぞれ後述するように、物理的に近くに位置するサイト間のイジング相互作用の大きさのみならず、物理的に遠くに位置するサイト間のイジング相互作用の大きさも自由に制御することができる。よって、サイト間の物理的距離とは無関係に、ＮＰ完全問題などからマッピングされた人工的なイジングモデルを解くことができる。【００１４】
第１の従来技術のイジングモデルの量子計算装置の概要を図１に示す。第２の従来技術のイジングモデルの量子計算装置の概要を図２に示す。【００１５】
イジングハミルトニアンを数式１のようにする。
【数１】
2015156126.tif 000003
【００１６】
イジング相互作用実装部Ｉ１２は、２つの面発光レーザーＶ１、Ｖ２の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ１、Ｖ２の間の擬似的なイジング相互作用Ｊ１２の大きさ及び符号を実装する。【００１７】
イジング相互作用実装部Ｉ１３は、２つの面発光レーザーＶ１、Ｖ３の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ１、Ｖ３の間の擬似的なイジング相互作用Ｊ１３の大きさ及び符号を実装する。【００１８】
イジング相互作用実装部Ｉ１４は、２つの面発光レーザーＶ１、Ｖ４の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ１、Ｖ４の間の擬似的なイジング相互作用Ｊ１４の大きさ及び符号を実装する。【００１９】
イジング相互作用実装部Ｉ２３は、２つの面発光レーザーＶ２、Ｖ３の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ２、Ｖ３の間の擬似的なイジング相互作用Ｊ２３の大きさ及び符号を実装する。【００２０】
イジング相互作用実装部Ｉ２４は、２つの面発光レーザーＶ２、Ｖ４の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ２、Ｖ４の間の擬似的なイジング相互作用Ｊ２４の大きさ及び符号を実装する。【００２１】
イジング相互作用実装部Ｉ３４は、２つの面発光レーザーＶ３、Ｖ４の間で交換される光の振幅及び位相を制御することにより、２つの面発光レーザーＶ３、Ｖ４の間の擬似的なイジング相互作用Ｊ３４の大きさ及び符号を実装する。【００２２】
マスターレーザーＭは、面発光レーザーＶ１〜Ｖ４に対して注入同期を行ない、面発光レーザーＶ１〜Ｖ４の発振周波数を同一周波数に揃える。面発光レーザーＶ１〜Ｖ４の間で同期を図ることにより、面発光レーザーＶ１〜Ｖ４のネットワーク全体で発振の偏光又は位相を最適化するにあたり、面発光レーザーＶ１〜Ｖ４のネットワーク全体で１つの発振モードを立ち上げることができる。【００２３】
第１の従来技術では、不図示のイジングスピン測定部は、面発光レーザーＶ１〜Ｖ４が光を交換する過程で定常状態に到達した後に、面発光レーザーＶ１〜Ｖ４の発振の円偏光の左回り／右回りを測定することにより、面発光レーザーＶ１〜Ｖ４の擬似的なイジングスピンσ１〜σ４の上向き／下向きを測定する。【００２４】
しかし、面発光レーザーＶは、面内異方性を有するため、左回り／右回りの円偏光のいずれについても、同一周波数及び同一閾値利得で発振することは困難である。よって、ある面発光レーザーＶが、単体のレーザーとしては、右回り（又は左回り）の円偏光を有する光を発振するよりも、左回り（又は右回り）の円偏光を有する光を発振しやすいことがあり得る。そして、その面発光レーザーＶは、レーザーのネットワーク全体としては、右回り（又は左回り）の円偏光を有する光を発振することが正答であるところ、左回り（又は右回り）の円偏光を有する光を発振してしまう誤答を生じさせ得る。【００２５】
第２の従来技術では、不図示のイジングスピン測定部は、面発光レーザーＶ１〜Ｖ４が光を交換する過程で定常状態に到達した後に、面発光レーザーＶ１〜Ｖ４の発振の直線偏光の位相の進み／遅れを測定することにより、面発光レーザーＶ１〜Ｖ４の擬似的なイジングスピンσ１〜σ４の上向き／下向きを測定する。【００２６】
ここで、左回り／右回りの円偏光は、水平偏光及び垂直偏光を位相差±π／２で同一の重みで重ね合わせたものである。つまり、イジングスピンの上向き／下向きの情報は、円偏光の左回り／右回りを測定するまでもなく、水平偏光を測定するまでもなく、垂直偏光の位相の進み／遅れを測定することにより、得ることができるのである。よって、第１の従来技術における面発光レーザーＶの面内異方性の問題を解決することができる。【００２７】
第２の従来技術のイジングモデルの量子計算装置の原理を図３に示す。マスターレーザーＭの直線偏光の発振位相０は、初期状態から定常状態まで変化しない。各面発光レーザーＶの直線偏光の発振位相φ（ｔ）は、初期状態においては、マスターレーザーＭの直線偏光の発振位相０と同一の０であることが理想であり、定常状態においては、マスターレーザーＭの直線偏光の発振位相０からずれた±π／２である。定常状態におけるφ（定常）＝±π／２は、σ＝±１に対応する（複号同順）。【００２８】
面発光レーザーＶの各ペアについて、イジング相互作用Ｊｉｊが正であるときには、２つの面発光レーザーＶの擬似的なスピンσが異符号であることが、エネルギー的に有利である。よって、各イジング相互作用実装部は、２つの面発光レーザーＶの発振位相φ（定常）が異符号でありずれをπとするような、発振モードが立ち上がりやすいようにする。【００２９】
面発光レーザーＶの各ペアについて、イジング相互作用Ｊｉｊが負であるときには、２つの面発光レーザーＶの擬似的なスピンσが同符号であることが、エネルギー的に有利である。よって、各イジング相互作用実装部は、２つの面発光レーザーＶの発振位相φ（定常）が同符号でありずれを０とするような、発振モードが立ち上がりやすいようにする。【００３０】
もっとも、イジングモデルの量子計算装置の全体において、一体として１つの発振モードが立ち上がるようにするのであり、面発光レーザーＶの各ペアにおいて、上述の発振モードが実際に立ち上がることもあれば、必ずしも立ち上がらないこともある。【００３１】
ところで、各面発光レーザーＶの直線偏光の発振位相φ（ｔ）は、初期状態においては、理想ではマスターレーザーＭの直線偏光の発振位相０と同一の０であることが望ましいが、実際にはマスターレーザーＭの直線偏光の発振位相０から若干ずれてしまう。【００３２】
各面発光レーザーＶの直線偏光の初期状態の発振位相φ（ｔ＝０）は、各面発光レーザーＶの自走周波数ω、マスターレーザーＭの発振周波数ω０及び注入同期幅ΔωＬ（ωがω０にどの程度近ければ注入同期を図れるか）を用いて、数式２のように表わされる。
【数２】
2015156126.tif 000004
【００３３】
つまり、各面発光レーザーＶの自走周波数ωを、マスターレーザーＭの発振周波数ω０に揃えることができるならば、各面発光レーザーＶの直線偏光の初期状態の発振位相φ（ｔ＝０）は、マスターレーザーＭの直線偏光の発振位相０と同一の０となる。しかし、各面発光レーザーＶの自走周波数ωを、マスターレーザーＭの発振周波数ω０に揃えることが困難であるため、各面発光レーザーＶの直線偏光の初期状態の発振位相φ（ｔ＝０）は、マスターレーザーＭの直線偏光の発振位相０から若干ずれてしまう。【００３４】
よって、ある面発光レーザーＶが、単体のレーザーとしては、マスターレーザーＭの直線偏光の発振位相０より遅れた（又は進んだ）発振位相を有する光を発振するよりも、マスターレーザーＭの直線偏光の発振位相０より進んだ（又は遅れた）発振位相を有する光を発振しやすいことがあり得る。そして、その面発光レーザーＶは、レーザーのネットワーク全体としては、マスターレーザーＭの直線偏光の発振位相０より遅れた（又は進んだ）発振位相を有する光を発振することが正答であるところ、マスターレーザーＭの直線偏光の発振位相０より進んだ（又は遅れた）発振位相を有する光を発振してしまう誤答を生じさせ得る。【００３５】
また、図２において、イジングサイトがＭ個であるとき、面発光レーザーＶはＭ個必要であり、イジング相互作用実装部はＭ（Ｍ−１）／２個必要である。そして、イジングサイトが多数になると、イジングモデルの量子計算装置が大規模かつ複雑になる。【００３６】
そこで、前記課題を解決するために、本発明は、イジングモデルの量子計算装置において、読み出しエラーを防止するとともに、回路構成を簡易にすることを目的とする。
【課題を解決するための手段】
【００３７】
上記目的を達成するために、パラメトリック発振を用いて、同一の発振周波数を有する複数の擬似的スピンパルスを発振し、複数の擬似的スピンパルスの発振位相の暫定的な測定結果を用いて、各擬似的スピンパルスが関わる相互作用の大きさ及び符号をフィードバック実装し、複数の擬似的スピンパルスの発振位相の最終的な測定結果に基づいて、複数の擬似的スピンパルスの擬似的なスピンを測定することとした。【００３８】
具体的には、本発明は、イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振器と、前記複数の擬似的スピンパルスを周回伝搬させるリング共振器と、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定部と、ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定部が暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算部と、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算部が暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装部と、前記暫定的スピン測定部、前記相互作用計算部及び前記相互作用実装部により構成されるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定部と、を備えることを特徴とするイジングモデルの量子計算装置である。【００３９】
この構成によれば、複数の擬似的スピンパルスは、同一の発振周波数を有するため、各擬似的スピンパルスの初期状態の発振位相が、各擬似的スピンパルスの定常状態の２種類の発振位相のうち、一方の位相へは近く他方の位相へは遠い、ということがない。よって、イジングモデルの量子計算装置において、読み出しエラーを防止することができる。【００４０】
そして、イジングサイトがＭ個であるとき、第１及び第２の従来技術では、面発光レーザーをＭ台も必要とするのに対して、本発明では、パラメトリック発振器を１台のみ準備すれば足りる。さらに、イジングサイトがＭ個であるとき、第１及び第２の従来技術では、イジング相互作用実装部をＭ（Ｍ−１）／２個も必要とするのに対して、本発明では、フィードバックループを１系統のみ準備すれば足りる。よって、イジングモデルの量子計算装置において、回路構成を簡易にすることができる。【００４１】
また、本発明は、前記相互作用計算部は、前記ある擬似的スピンパルスが関わる前記イジングモデルの３体以上の結合係数及び前記暫定的スピン測定部が暫定的に測定した前記他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる３体以上の相互作用を暫定的に計算し、前記相互作用実装部は、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算部が暫定的に計算した前記ある擬似的スピンパルスが関わる３体以上の相互作用の大きさ及び符号を暫定的に実装することを特徴とするイジングモデルの量子計算装置である。【００４２】
この構成によれば、各擬似的スピンパルス及び各注入光パルスの線形の重ね合わせの範囲内で、イジングモデルの３体以上の相互作用を実装することができる。【００４３】
また、本発明は、前記パラメトリック発振器は、前記複数の擬似的スピンパルスと同一の発振周波数を有し１対１のペアを組む複数の局部発振パルスをパラメトリック発振させ、前記リング共振器は、前記複数の局部発振パルスを周回伝搬させ、前記複数の局部発振パルスの位相について、正相及び逆相の両方を含む状態から、正相及び逆相の一方を含む状態へと、同相化させるパルス同相化部、をさらに備え、前記暫定的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、１対１のペアを組む前記複数の局部発振パルスの一部を用いてホモダイン検波を行い、前記相互作用実装部は、前記ある擬似的スピンパルスに対して、１対１のペアを組みその一部の振幅及び位相を制御された局部発振パルスを注入し、前記擬似的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、１対１のペアを組む前記複数の局部発振パルスの一部を用いてホモダイン検波を行うことを特徴とするイジングモデルの量子計算装置である。【００４４】
この構成によれば、複数の擬似的スピンパルスの発振位相の暫定的な測定結果を用いた、各擬似的スピンパルスが関わる相互作用の大きさ及び符号のフィードバック実装を、具体的な構成を用いて行うことができる。そして、この構成によれば、次の構成と比べて、擬似的スピンパルス及び局部発振パルスがペアとなってリング共振器を周回伝搬するため、パルス発生器からスピン測定部への光路長の揺らぎの問題をなくすことができる。よって、イジングモデルの量子計算装置において、読み出しエラーを防止することができる。【００４５】
また、本発明は、角周波数ωを有する局部発振パルスを発生させるパルス発生器と、前記角周波数ωを有する局部発振パルスを用いて、角周波数２ωを有するパルスを発生させる第二高調波発生器と、をさらに備え、前記パラメトリック発振器は、前記角周波数２ωを有するパルスを用いて、前記複数の擬似的スピンパルスをパラメトリック発振させ、前記暫定的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、前記角周波数ωを有する局部発振パルスを用いてホモダイン検波を行い、前記相互作用実装部は、前記ある擬似的スピンパルスに対して、振幅及び位相を制御された前記角周波数ωを有する局部発振パルスを注入し、前記擬似的スピン測定部は、前記複数の擬似的スピンパルスの一部に対して、前記角周波数ωを有する局部発振パルスを用いてホモダイン検波を行うことを特徴とするイジングモデルの量子計算装置である。【００４６】
この構成によれば、複数の擬似的スピンパルスの発振位相の暫定的な測定結果を用いた、各擬似的スピンパルスが関わる相互作用の大きさ及び符号のフィードバック実装を、具体的な構成を用いて行うことができる。そして、この構成によれば、先の構成と比べて、擬似的スピンパルス及び局部発振パルスがペアとなってリング共振器を周回伝搬しないため、擬似的スピンパルス及び局部発振パルス間のクロストークをなくすことができ、全局部発振パルスを同相化させるパルス同相化部をなくすことができる。よって、イジングモデルの量子計算装置において、読み出しエラーを防止することができる。【００４７】
また、本発明は、前記相互作用実装部は、前記ある擬似的スピンパルスに対して注入される光の振幅を、計算処理の初期段階ほど大きく制御し、計算処理の終期段階ほど小さく制御することを特徴とするイジングモデルの量子計算装置である。【００４８】
この構成によれば、計算処理の初期段階において、各擬似的スピンパルスの読出結果をできるだけ正しい答えにすることができる。たとえ、計算処理の初期段階において、ある擬似的スピンパルスの読出結果が正しい答えでないとしても、他の擬似的スピンパルスについての正しい答えをフィードバックすれば、計算処理の終期段階において、全擬似的スピンパルスについての正しい答えを得ることができる。よって、イジングモデルの量子計算装置において、読み出しエラーを防止することができる。【００４９】
また、本発明は、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの振幅を暫定的に測定する暫定的振幅測定部、をさらに備え、前記パラメトリック発振器は、前記暫定的振幅測定部が測定した前記複数の擬似的スピンパルスの振幅に基づいて、前記複数の擬似的スピンパルスの振幅が等しくなるように、パラメトリック発振に用いるポンプパルスの振幅をフィードバック制御することを特徴とするイジングモデルの量子計算装置である。【００５０】
この構成によれば、各擬似的スピンパルスの振幅の不均衡による、イジングモデルの結合係数の実質的な書き換えの問題をなくすことができる。よって、イジングモデルの量子計算装置において、読み出しエラーを防止することができる。【００５１】
また、本発明は、前記イジングモデルの複数のスピンのうち第ｎ群（ｎは１以上の整数）のスピンに擬似的に対応する第ｎ群の擬似的スピンパルスが前記リング共振器を周回伝搬する、請求項１から６のいずれかに記載の第ｎ番のイジングモデルの量子計算装置、を並列に備え、前記第ｎ番のイジングモデルの量子計算装置が備える前記暫定的スピン測定部が暫定的に測定した前記第ｎ群の擬似的スピンパルスの擬似的なスピンの情報を、並列に備わるイジングモデルの量子計算装置の間で共有させる暫定的スピン共有部、をさらに備えることを特徴とするイジングモデルの量子並列計算装置である。【００５２】
この構成によれば、イジングモデルのサイト数が多いときでも、複数のイジングモデルの量子計算装置が並列分散処理を実行することにより、各々のイジングモデルの量子計算装置は計算処理負担を軽減することができる。【００５３】
また、本発明は、イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振ステップと、前記複数の擬似的スピンパルスがリング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定ステップが暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算ステップが暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスをパラメトリック増幅させるパラメトリック増幅ステップと、前記暫定的スピン測定ステップ、前記相互作用実装ステップ及び前記パラメトリック増幅ステップをこの順序で備えるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定ステップと、を備えることを特徴とするイジングモデルの量子計算方法である。【００５４】
この構成によれば、暫定的スピン測定ステップと相互作用実装ステップの間に、パラメトリック増幅ステップが入らず、タイムラグがほとんど生じないため、イジングモデルのサイト間のほとんど遅延のない相互作用を実装することができる。【００５５】
また、本発明は、イジングモデルの複数のスピンに擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスをパラメトリック発振させるパラメトリック発振ステップと、前記複数の擬似的スピンパルスがリング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスの位相を暫定的に測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを暫定的に測定する暫定的スピン測定ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、ある擬似的スピンパルスが関わる前記イジングモデルの結合係数及び前記暫定的スピン測定ステップが暫定的に測定した他の擬似的スピンパルスの擬似的なスピンに基づいて、前記ある擬似的スピンパルスが関わる相互作用を暫定的に計算する相互作用計算ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記複数の擬似的スピンパルスをパラメトリック増幅させるパラメトリック増幅ステップと、前記複数の擬似的スピンパルスが前記リング共振器を周回伝搬するたびに、前記ある擬似的スピンパルスに対して注入される光の振幅及び位相を制御することにより、前記相互作用計算ステップが暫定的に計算した前記ある擬似的スピンパルスが関わる相互作用の大きさ及び符号を暫定的に実装する相互作用実装ステップと、前記暫定的スピン測定ステップ、前記パラメトリック増幅ステップ及び前記相互作用実装ステップをこの順序で備えるフィードバックループが繰り返される過程で、前記複数の擬似的スピンパルスが定常状態に到達した後に、前記複数の擬似的スピンパルスの位相を測定することにより、前記複数の擬似的スピンパルスの擬似的なスピンを測定する擬似的スピン測定ステップと、を備えることを特徴とするイジングモデルの量子計算方法である。【００５６】
この構成によれば、暫定的スピン測定ステップと相互作用実装ステップの間に、パラメトリック増幅ステップが入って、タイムラグがある程度生じるものの、イジングモデルのサイト間の実質には遅延のない相互作用を実装することができる。
【発明の効果】
【００５７】
本発明は、イジングモデルの量子計算装置において、読み出しエラーを防止するとともに、回路構成を簡易にすることができる。
【図面の簡単な説明】
【００５８】
【図１】第１の従来技術のイジングモデルの量子計算装置の概要を示す図である。
【図２】第２の従来技術のイジングモデルの量子計算装置の概要を示す図である。
【図３】第２の従来技術のイジングモデルの量子計算装置の原理を示す図である。
【図４】本発明のイジングモデルの量子計算装置の構成を示す図である。
【図５】本発明のイジングモデルの量子計算装置の原理を示す図である。
【図６】第１の実施形態のイジングモデルの量子計算装置の構成を示す図である。
【図７】第２の実施形態のイジングモデルの量子計算装置の構成を示す図である。
【図８】本発明のイジングモデルの量子計算装置の計算結果を示す図である。
【図９】本発明のイジングモデルの量子計算装置の時間発展を示す図である。
【図１０】本発明のイジングモデルの量子計算装置の計算結果を示す図である。
【図１１】本発明のイジングモデルの量子計算装置の時間発展を示す図である。
【図１２】計算処理の全段階で光注入強度を一定にした場合の時間発展を示す図である。
【図１３】計算処理の初期段階に光注入強度を高くした場合の時間発展を示す図である。
【図１４】複数の擬似的スピンパルスの強度を揃えない場合の時間発展を示す図である。
【図１５】複数の擬似的スピンパルスの強度を揃えない場合の計算精度を示す図である。
【図１６】複数の擬似的スピンパルスの強度を揃えた場合の時間発展を示す図である。
【図１７】複数の擬似的スピンパルスの強度を揃えた場合の計算精度を示す図である。
【図１８】本発明のイジングモデルの量子計算方法のループ手順を示す図である。
【図１９】本発明のイジングモデルの量子計算方法のループ手順を示す図である。
【図２０】本発明のイジングモデルの量子並列計算装置の構成を示す図である。
【図２１】本発明のイジングモデルの量子並列計算装置の構成を示す図である。
【発明を実施するための形態】
【００５９】
添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施の例であり、本発明は以下の実施形態に制限されるものではない。なお、本明細書及び図面において符号が同じ構成要素は、相互に同一のものを示すものとする。
【００６０】
（本発明のイジングモデルの量子計算装置の構成及び原理）
本発明のイジングモデルの量子計算装置Ｑの構成を図４に示す。本発明では、イジングハミルトニアンを、１体〜３体の相互作用を含むとして、数式３のようにする。
【数３】
2015156126.tif 000005
【００６１】
パラメトリック発振器１は、イジングモデルの複数のスピンσ１〜σ４に擬似的に対応し同一の発振周波数を有する複数の擬似的スピンパルスＳＰ１〜ＳＰ４をパラメトリック発振させる。リング共振器２は、複数の擬似的スピンパルスＳＰ１〜ＳＰ４を周回伝搬させる。複数の擬似的スピンパルスＳＰ１〜ＳＰ４は、ＳＰ１、ＳＰ２、ＳＰ３、ＳＰ４、ＳＰ１、ＳＰ２、ＳＰ３、ＳＰ４、・・・の順序で、後述のフィードバックループに入る。
【００６２】
暫定的スピン測定部３は、複数の擬似的スピンパルスＳＰ１〜ＳＰ４がリング共振器２を周回伝搬するたびに、複数の擬似的スピンパルスＳＰ１〜ＳＰ４の位相を暫定的に測定することにより、複数の擬似的スピンパルスＳＰ１〜ＳＰ４の擬似的なスピンσ１〜σ４を暫定的に測定する。具体的には、暫定的スピン測定部３は、図６及び図７を用いて後述する局部発振パルスＬＯを用いて、ホモダイン検波を行う。
【００６３】
相互作用計算部４は、ある擬似的スピンパルスＳＰｉが関わるイジングモデルの結合係数λｉ、Ｊｉｊ、Ｋｉｊｋ及び暫定的スピン測定部３が暫定的に測定した他の擬似的スピンパルスＳＰｊ、ＳＰｋの擬似的なスピンσｊ、σｋに基づいて、ある擬似的スピンパルスＳＰｉが関わる相互作用（σｉに対する比例係数λｉ＋ΣＪｉｊσｊ＋ΣＫｉｊｋσｊσｋ）を暫定的に計算する。図４では、ｉ、ｊ、ｋ＝１〜４の場合を示している。
【００６４】
ここで、ＮＰ完全問題などが、イジングモデルにマッピングされた後、相互作用計算部４は、イジングモデルの結合係数λｉ、Ｊｉｊ、Ｋｉｊｋを入力する。
【００６５】
相互作用実装部５は、ある擬似的スピンパルスＳＰｉに対して注入される光の振幅及び位相を制御することにより、相互作用計算部４が暫定的に計算したある擬似的スピンパルスＳＰｉが関わる相互作用（σｉに対する比例係数λｉ＋ΣＪｉｊσｊ＋ΣＫｉｊｋσｊσｋ）の大きさ及び符号を暫定的に実装する。具体的には、相互作用実装部５は、図６及び図７を用いて後述する局部発振パルスＬＯを用いて、注入光パルスを生成する。
【００６６】
擬似的スピン測定部６は、暫定的スピン測定部３、相互作用計算部４及び相互作用実装部５により構成されるフィードバックループが繰り返される過程で、複数の擬似的スピンパルスＳＰ１〜ＳＰ４が定常状態に到達した後に、複数の擬似的スピンパルスＳＰ１〜ＳＰ４の位相を測定することにより、複数の擬似的スピンパルスＳＰ１〜ＳＰ４の擬似的なスピンσ１〜σ４を測定する。具体的には、擬似的スピン測定部６は、図６及び図７を用いて後述する局部発振パルスＬＯを用いて、ホモダイン検波を行う。
【００６７】
ここで、擬似的スピン測定部６が、イジングモデルのスピンσ１〜σ４を出力した後、イジングモデルは、ＮＰ完全問題などにデマッピングされる。
【００６８】
このように、パラメトリック発振器１でポンピング電流を漸増制御し、複数の擬似的スピンパルスＳＰ１〜ＳＰ４のネットワーク全体で閾値利得が最低となる１つの発振モードを立ち上げ、複数の擬似的スピンパルスＳＰ１〜ＳＰ４の発振位相を測定し、複数の擬似的スピンパルスＳＰ１〜ＳＰ４に対応する各原子のスピンの方向を測定する。
【００６９】
本発明のイジングモデルの量子計算装置Ｑの原理を図５に示す。局部発振パルスＬＯの発振位相０は、初期状態から定常状態まで変化しない。各擬似的スピンパルスＳＰの発振位相φ（ｔ）は、初期状態においては、０又は±πであり（各擬似的スピンパルスＳＰは、パラメトリック発振器１により、パラメトリック発振されて、スクイーズド状態にある。）、定常状態においては、初期状態の発振位相０又は±πからずれた±π／２である。定常状態におけるφ（定常）＝±π／２は、σ＝±１に対応する（複号同順）。
【００７０】
各擬似的スピンパルスＳＰについて、１体の相互作用の結合係数λｉが正であるときには、当該擬似的スピンパルスＳＰの擬似的なスピンσが−１であることが、エネルギー的に有利である。よって、相互作用実装部５は、当該擬似的スピンパルスＳＰの発振位相φ（定常）が−π／２であるような、発振モードが立ち上がりやすいようにする。
【００７１】
各擬似的スピンパルスＳＰについて、１体の相互作用の結合係数λｉが負であるときには、当該擬似的スピンパルスＳＰの擬似的なスピンσが＋１であることが、エネルギー的に有利である。よって、相互作用実装部５は、当該擬似的スピンパルスＳＰの発振位相φ（定常）が＋π／２であるような、発振モードが立ち上がりやすいようにする。
【００７２】
２つの擬似的スピンパルスＳＰについて、２体の相互作用の結合係数Ｊｉｊが正であるときには、２つの擬似的スピンパルスＳＰの擬似的なスピンσが異符号であることが、エネルギー的に有利である。よって、相互作用実装部５は、２つの擬似的スピンパルスＳＰの発振位相φ（定常）が異符号であるような、発振モードが立ち上がりやすいようにする。
【００７３】
２つの擬似的スピンパルスＳＰについて、２体の相互作用の結合係数Ｊｉｊが負であるときには、２つの擬似的スピンパルスＳＰの擬似的なスピンσが同符号であることが、エネルギー的に有利である。よって、相互作用実装部５は、２つの擬似的スピンパルスＳＰの発振位相φ（定常）が同符号であるような、発振モードが立ち上がりやすいようにする。
【００７４】
３つの擬似的スピンパルスＳＰについて、３体の相互作用の結合係数Ｋｉｊｋが正であるときには、（１）３つの擬似的スピンパルスＳＰの擬似的なスピンσが−１であること、又は、（２）２つの擬似的スピンパルスＳＰの擬似的なスピンσが＋１であり、１つの擬似的スピンパルスＳＰの擬似的なスピンσが−１であることが、エネルギー的に有利である。よって、相互作用実装部５は、（１）３つの擬似的スピンパルスＳＰの発振位相φ（定常）が−π／２であるような、又は、（２）２つの擬似的スピンパルスＳＰの発振位相φ（定常）が＋π／２であり、１つの擬似的スピンパルスＳＰの発振位相φ（定常）が−π／２であるような、発振モードが立ち上がりやすいようにする。
【００７５】
３つの擬似的スピンパルスＳＰについて、３体の相互作用の結合係数Ｋｉｊｋが負であるときには、（１）３つの擬似的スピンパルスＳＰの擬似的なスピンσが＋１であること、又は、（２）２つの擬似的スピンパルスＳＰの擬似的なスピンσが−１であり、１つの擬似的スピンパルスＳＰの擬似的なスピンσが＋１であることが、エネルギー的に有利である。よって、相互作用実装部５は、（１）３つの擬似的スピンパルスＳＰの発振位相φ（定常）が＋π／２であるような、又は、（２）２つの擬似的スピンパルスＳＰの発振位相φ（定常）が−π／２であり、１つの擬似的スピンパルスＳＰの発振位相φ（定常）が＋π／２であるような、発振モードが立ち上がりやすいようにする。
【００７６】
もっとも、イジングモデルの量子計算装置Ｑの全体において、一体として１つの発振モードが立ち上がるようにするのであり、各擬似的スピンパルスＳＰにおいて、上述の発振モードが実際に立ち上がることもあれば、必ずしも立ち上がらないこともある。
【００７７】
図４及び図５で示した計算原理について詳述する。各擬似的スピンパルスＳＰ１、ＳＰ２、ＳＰ３、ＳＰ４において、発振強度Ａｉ（ｔ）、発振位相φｉ（ｔ）及びキャリアの反転分布数差ＮＣｉ（ｔ）について、レート方程式は、数式４−７のようになる。
【数４】
2015156126.tif 000006
【数５】
2015156126.tif 000007
【数６】
2015156126.tif 000008
【数７】
2015156126.tif 000009
【００７８】
ωは、発振周波数である。Ｑは、各擬似的スピンパルスＳＰの共振器Ｑ値である。Ｐは、反転分布を実現するために各擬似的スピンパルスＳＰについて毎秒注入される電子数、すなわちポンピングレートである。数式４の−（１／２）（ω／Ｑ）Ａｉ（ｔ）は、発振強度Ａｉ（ｔ）が共振器損失に起因して時間経過につれて減少するレートを示している。
【００７９】
τｓｐは、レーザー発振モード以外の他の発振モードへの自然放出に起因する電子寿命である。βは、全自然放出光のうちレーザー発振モードへの結合定数である。数式４の（１／２）ＥＣｉ（ｔ）Ａｉ（ｔ）は、発振強度Ａｉ（ｔ）が誘導放出に起因して時間経過につれて増加するレートを示している。数式４のＥＣｉ（ｔ）は、発振強度Ａｉ（ｔ）が自然放出に起因して時間経過につれて増加するレートを示している。
【００８０】
数式４、５のλｉが関わる項は、１体の相互作用に関わる項である。相互作用実装部５が、擬似的スピンパルスＳＰｉに対して、１体の相互作用（σｉに対する比例係数λｉ）を実装するための注入光パルスを生成する方法を説明する。
【００８１】
相互作用計算部４は、１体の相互作用（σｉに対する比例係数λｉ）を計算する。λｉが正であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をπ／２だけ遅くする位相変調を行い、｜λｉ｜に比例する振幅変調を行い、注入光パルスを生成する。λｉが負であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をπ／２だけ早くする位相変調を行い、｜λｉ｜に比例する振幅変調を行い、注入光パルスを生成する。
【００８２】
数式４の（ω／Ｑ）Ａ｛−λｉｓｉｎφｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、１体の相互作用（σｉに対する比例係数λｉ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振強度Ａｉ（ｔ）が時間経過につれて変化するレートを示す。なお、数式４において、Ａは、比例定数である。
【００８３】
数式５の（１／Ａｉ（ｔ））（ω／Ｑ）Ａ｛−λｉｃｏｓφｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、１体の相互作用（σｉに対する比例係数λｉ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振位相φｉ（ｔ）が時間経過につれて変化するレートを示す。なお、数式５において、Ａは、比例定数である。
【００８４】
数式４、５のＪｉｊが関わる項は、２体の相互作用に関わる項である。相互作用実装部５が、擬似的スピンパルスＳＰｉに対して、２体の相互作用（σｉに対する比例係数ΣＪｉｊσｊ）を実装するための注入光パルスを生成する方法を説明する。
【００８５】
暫定的スピン測定部３は、本周回前に、擬似的スピンパルスＳＰｊの発振位相φｊ（ｔ）及び擬似的なスピンσｊを測定している。相互作用計算部４は、２体の相互作用（σｉに対する比例係数ΣＪｉｊσｊ）を計算する。ｉ、ｊ番目のサイト間について、Ｊｉｊが正であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をφｊ（ｔ）に移して更なる逆相化を施す位相変調を行い、｜Ｊｉｊ｜に比例する振幅変調を行い、注入光パルスを生成する。ｉ、ｊ番目のサイト間について、Ｊｉｊが負であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をφｊ（ｔ）に移すが更なる逆相化を施さない位相変調を行い、｜Ｊｉｊ｜に比例する振幅変調を行い、注入光パルスを生成する。相互作用実装部５は、ｉ、ｊ番目のサイト間の全組み合わせについて、上述のように注入光パルスを生成する。
【００８６】
数式４の−（ω／Ｑ）（１／２）ＪｉｊＡｃｏｓ｛φｊ（ｔ）−φｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、２体の相互作用（σｉに対する比例係数ΣＪｉｊσｊ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振強度Ａｉ（ｔ）が時間経過につれて変化するレートを示す。数式４のΣ（ｊ≠ｉ）は、ｉ番目のサイトにおける、ｉ番目以外の他の全てのサイト（ｊ番目）からの寄与を示す。
【００８７】
数式５の−（１／Ａｉ（ｔ））（ω／Ｑ）（１／２）ＪｉｊＡｓｉｎ｛φｊ（ｔ）−φｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、２体の相互作用（σｉに対する比例係数ΣＪｉｊσｊ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振位相φｉ（ｔ）が時間経過につれて変化するレートを示す。数式５のΣ（ｊ≠ｉ）は、ｉ番目のサイトにおける、ｉ番目以外の他の全てのサイト（ｊ番目）からの寄与を示す。
【００８８】
数式４、５のＫｉｊｋが関わる項は、３体の相互作用に関わる項である。相互作用実装部５が、擬似的スピンパルスＳＰｉに対して、３体の相互作用（σｉに対する比例係数ΣＫｉｊｋσｊσｋ）を実装するための注入光パルスを生成する方法を説明する。
【００８９】
暫定的スピン測定部３は、本周回前に、擬似的スピンパルスＳＰｊ、ＳＰｋの発振位相φｊ（ｔ）、φｋ（ｔ）及び擬似的なスピンσｊ、σｋを測定している。相互作用計算部４は、３体の相互作用（σｉに対する比例係数ΣＫｉｊｋσｊσｋ）を計算する。ｉ、ｊ、ｋ番目のサイト間について、Ｋｉｊｋが正であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をφｊｋ（ｔ）に移して更なる逆相化を施す位相変調を行い、｜Ｋｉｊｋ｜に比例する振幅変調を行い、注入光パルスを生成する。ｉ、ｊ、ｋ番目のサイト間について、Ｋｉｊｋが負であるときには、相互作用実装部５は、局部発振パルスＬＯ（発振位相０）に対して、発振位相をφｊｋ（ｔ）に移すが更なる逆相化を施さない位相変調を行い、｜Ｋｉｊｋ｜に比例する振幅変調を行い、注入光パルスを生成する。相互作用実装部５は、ｉ、ｊ、ｋ番目のサイト間の全組み合わせについて、上述のように注入光パルスを生成する。なお、φｊｋ（ｔ）は、数式１２を用いて後述する。
【００９０】
数式４の−（ω／Ｑ）（１／２）ＫｉｊｋＡｃｏｓ｛φｊｋ（ｔ）−φｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、３体の相互作用（σｉに対する比例係数ΣＫｉｊｋσｊσｋ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振強度Ａｉ（ｔ）が時間経過につれて変化するレートを示す。数式４のΣ（ｊ、ｋ≠ｉ）は、ｉ番目のサイトにおける、ｉ番目以外の他の全てのサイト（ｊ、ｋ番目）からの寄与を示す。
【００９１】
数式５の−（１／Ａｉ（ｔ））（ω／Ｑ）（１／２）ＫｉｊｋＡｓｉｎ｛φｊｋ（ｔ）−φｉ（ｔ）｝は、擬似的スピンパルスＳＰｉに対して、３体の相互作用（σｉに対する比例係数ΣＫｉｊｋσｊσｋ）を実装するための注入光パルスが生成されたときに、ｉ番目のサイトにおける発振位相φｉ（ｔ）が時間経過につれて変化するレートを示す。数式５のΣ（ｊ、ｋ≠ｉ）は、ｉ番目のサイトにおける、ｉ番目以外の他の全てのサイト（ｊ、ｋ番目）からの寄与を示す。
【００９２】
ＦＡ、Ｆφ及びＦＮは、それぞれ、ｉ番目のサイトにおける、発振強度Ａｉ（ｔ）、発振位相φｉ（ｔ）及びキャリアの反転分布数差ＮＣｉ（ｔ）に対する雑音を示す。
【００９３】
定常状態において、数式４は、数式８のようになる。
【数８】
2015156126.tif 000010
ＦＡを無視して、数式８を変形すると、数式９のようになる。
【数９】
2015156126.tif 000011
【００９４】
ここで、数式１０、１１が成立している。そこで、イジングモデル及びレーザーシステムの類似に着目して、数式１２のようにおくと、数式９は、数式１３のようになる。
【数１０】
2015156126.tif 000012
【数１１】
2015156126.tif 000013
【数１２】
2015156126.tif 000014
【数１３】
2015156126.tif 000015
【００９５】
ここで、擬似的スピンパルスＳＰｉに対して、２体の相互作用（σｉに対する比例係数ΣＪｉｊσｊ）を実装するためには、擬似的スピンパルスＳＰｉ、ＳＰｊを線形に重ね合わせるのみで足りる。しかし、擬似的スピンパルスＳＰｉに対して、３体の相互作用（σｉに対する比例係数ΣＫｉｊｋσｊσｋ）を実装するためには、擬似的スピンパルスＳＰｉ、ＳＰｊ、ＳＰｋを線形に重ね合わせるのみでは足らない。
【００９６】
しかし、擬似的スピンパルスＳＰｉ、ＳＰｊ、ＳＰｋの間の非線形効果を利用すれば、イジングモデルの量子計算装置Ｑの回路構成が複雑になる。そこで、数式１２のようにσｊσｋ＝ｓｉｎφｊｋ（ｔ）とおけば、擬似的スピンパルスＳＰｉ及び注入光パルスの間の線形の重ね合わせが利用できて、イジングモデルの量子計算装置Ｑの回路構成が簡易になる。
【００９７】
数式１３をＭ個の全部のサイトについて加算すると、数式１４のようになり、レーザーシステムの全体の閾値利得ΣＥＣｉとして表わされる。
【数１４】
2015156126.tif 000016
【００９８】
ここで、仮定として、数式１５が成立している。そして、定常状態では、数式１６、１７が成立している。このとき、数式１４は、数式１８のようになる。
【数１５】
2015156126.tif 000017
【数１６】
2015156126.tif 000018
【数１７】
2015156126.tif 000019
【数１８】
2015156126.tif 000020
【００９９】
ここで、レーザーの媒質が均一な媒質であるときには、レーザーシステム全体として、最小の閾値利得ΣＥＣｉを実現する発振位相状態｛σｉ｝が選択される。つまり、レーザーシステム全体として、１個の特定の発振モードが選択される。そして、発振モードの間の競合に起因して、１個の特定の発振モードは、他の発振モードを抑制する。つまり、レーザーシステム全体として、数式１８のΣＥＣｉは最小化される。一方で、レーザーシステム全体として、数式１８の（ω／Ｑ）Ｍは一定である。よって、レーザーシステム全体として、数式１８のΣλｉσｉ＋ΣＪｉｊσｉσｊ＋ΣＫｉｊｋσｉσｊσｋは最小化される。つまり、数式３のイジングハミルトニアンを最小化する基底状態が実現されたことになる。
【０１００】
ここで、計算精度を向上させるためには、レーザーシステム全体でのレーザー発振モードについて、最小の閾値利得及びその次に小さい閾値利得の差を、自然放出レートで決定される飽和利得ＥＣ及び光子減衰率ω／Ｑの差であるβ（ω／Ｑ）（１／Ｒ）より十分に大きくする必要がある。ここで、Ｒ＝Ｉ／Ｉｔｈ−１は、規格化ポンプレートであり、Ｉ及びＩｔｈは、それぞれ注入電流及びそのレーザー発振閾値である。よって、βを小さくしＲを高くすることにより、計算精度を向上させることができる。
【０１０１】
図４及び図５を用いて説明したように、複数の擬似的スピンパルスＳＰｉは、同一の発振周波数を有するため、各擬似的スピンパルスＳＰｉの初期状態の発振位相φｉ（ｔ＝０）が、各擬似的スピンパルスＳＰｉの定常状態の２種類の発振位相φｉ（定常）＝±π／２のうち、一方の位相へは近く他方の位相へは遠い、ということがない。よって、イジングモデルの量子計算装置Ｑにおいて、読み出しエラーを防止することができる。
【０１０２】
そして、イジングサイトがＭ個であるとき、第１及び第２の従来技術では、面発光レーザーＶをＭ台も必要とするのに対して、本発明では、パラメトリック発振器１を１台のみ準備すれば足りる。さらに、イジングサイトがＭ個であるとき、第１及び第２の従来技術では、イジング相互作用実装部をＭ（Ｍ−１）／２個も必要とするのに対して、本発明では、フィードバックループを１系統のみ準備すれば足りる。よって、イジングモデルの量子計算装置Ｑにおいて、回路構成を簡易にすることができる。
【０１０３】
さらに、各擬似的スピンパルスＳＰｉ及び各注入光パルスの線形の重ね合わせの範囲内で、イジングモデルの３体以上の相互作用を実装することができる。ここで、イジングモデルの３体の相互作用を実装するときには、数式１２のようにσｊσｋ＝ｓｉｎφｊｋ（ｔ）とおいて、擬似的スピンパルスＳＰｉ及び注入光パルスの線形の重ね合わせを行う。そして、イジングモデルの３体以上の相互作用を実装するときには、数式１２と同様にσｊσｋσｌ・・・＝ｓｉｎφｊｋｌ・・・（ｔ）（Ｎ体の相互作用について、φは（Ｎ−１）個のσの積）とおいて、擬似的スピンパルスＳＰｉ及び注入光パルスの線形の重ね合わせを行う。
【０１０４】
（第１の実施形態のイジングモデルの量子計算装置の構成）
第１の実施形態のイジングモデルの量子計算装置Ｑの構成を図６に示す。
【０１０５】
パラメトリック発振器１は、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）と同一の発振周波数を有し１対１のペアを組む複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）をパラメトリック発振させる。リング共振器２は、複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）を周回伝搬させる。複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）及び複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）は、ＬＯ１、ＳＰ１、・・・、ＬＯ（Ｎ）、ＳＰ（Ｎ）、ＬＯ１、ＳＰ１、・・・、ＬＯ（Ｎ）、ＳＰ（Ｎ）、・・・の順序で、フィードバックループに入る。
【０１０６】
ここで、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）は、発振閾値の近傍でポンプパルスを用いてパラメトリック発振される。そして、複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）は、発振閾値より十分上でポンプパルスを用いてパラメトリック発振される。
【０１０７】
パルス同相化部７は、複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）の位相について、正相及び逆相の両方を含む状態から、正相及び逆相の一方を含む状態へと、同相化させる。よって、パルス同相化部７は、複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）の位相について、図５で示したように０とすることができる。例えば、パルス同相化部７は、２パルス分（隣接の局部発振パルスＬＯは、１つの擬似的スピンパルスＳＰを挟む。）の遅延線及び位相変調器を備える。ただし、パルス同相化部７は、複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）のみ位相変調すべきであり、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）まで位相変調すべきでなく、位相変調器のＯＮ／ＯＦＦの切り替えを行う。
【０１０８】
暫定的スピン測定部３は、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）の一部に対して、１対１のペアを組む複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）の一部を用いてホモダイン検波を行う。ここで、暫定的スピン測定部３は、ｃｏｓ成分及びｓｉｎ成分を測定しており、各成分の測定部に１パルス分（１対１のペアを組む擬似的スピンパルスＳＰ及び局部発振パルスＬＯは、隣接している。）の遅延線を備える。
【０１０９】
相互作用計算部４は、例えば、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）であり、イジングモデルの結合係数λｉ、Ｊｉｊ、Ｋｉｊｋを入力する。相互作用実装部５は、ある擬似的スピンパルスＳＰｉに対して、１対１のペアを組みその一部の振幅及び位相を制御された局部発振パルスＬＯｉを注入する。ここで、相互作用実装部５は、１パルス分（１対１のペアを組む擬似的スピンパルスＳＰ及び局部発振パルスＬＯは、隣接している。）の遅延線上に置かれる。
【０１１０】
擬似的スピン測定部６は、暫定的スピン測定部３と共用であり、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）の一部に対して、１対１のペアを組む複数の局部発振パルスＬＯ１〜ＬＯ（Ｎ）の一部を用いてホモダイン検波を行い、イジングモデルのスピンσを出力する。
【０１１１】
図６を用いて説明したように、複数の擬似的スピンパルスＳＰｉの発振位相φｉ（ｔ）の暫定的な測定結果を用いた、各擬似的スピンパルスＳＰｉが関わる相互作用の大きさ及び符号のフィードバック実装を、具体的な構成を用いて行うことができる。そして、第１の実施形態によれば、第２の実施形態と比べて、擬似的スピンパルスＳＰｉ及び局部発振パルスＬＯｉがペアとなってリング共振器２を周回伝搬するため、パルス発生器８からスピン測定部３、６への光路長の揺らぎの問題をなくすことができる。よって、イジングモデルの量子計算装置Ｑにおいて、読み出しエラーを防止することができる。
【０１１２】
（第２の実施形態のイジングモデルの量子計算装置の構成）
第２の実施形態のイジングモデルの量子計算装置Ｑの構成を図７に示す。
【０１１３】
パルス発生器８は、角周波数ωを有する局部発振パルスＬＯを発生させる。第二高調波発生器９は、角周波数ωを有する局部発振パルスＬＯを用いて、角周波数２ωを有するパルスを発生させる。パラメトリック発振器１は、角周波数２ωを有するパルスを用いて、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）をパラメトリック発振させる。複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）は、ＳＰ１、・・・、ＳＰ（Ｎ）、ＳＰ１、・・・、ＳＰ（Ｎ）、・・・の順序で、フィードバックループに入る。
【０１１４】
ここで、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）は、発振閾値の近傍でポンプパルスを用いてパラメトリック発振される。そして、パルス発生器８は、例えば、モード同期レーザーを備える。よって、パルス発生器８は、角周波数ωを有する局部発振パルスＬＯの位相について、図５で示したように０とすることができる。
【０１１５】
暫定的スピン測定部３は、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）の一部に対して、角周波数ωを有する局部発振パルスＬＯを用いてホモダイン検波を行う。ここで、暫定的スピン測定部３は、ｃｏｓ成分及びｓｉｎ成分を測定している。
【０１１６】
相互作用計算部４は、例えば、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）であり、イジングモデルの結合係数λｉ、Ｊｉｊ、Ｋｉｊｋを入力する。相互作用実装部５は、ある擬似的スピンパルスＳＰｉに対して、振幅及び位相を制御された角周波数ωを有する局部発振パルスＬＯを注入する。
【０１１７】
擬似的スピン測定部６は、暫定的スピン測定部３と共用であり、複数の擬似的スピンパルスＳＰ１〜ＳＰ（Ｎ）の一部に対して、角周波数ωを有する局部発振パルスＬＯを用いてホモダイン検波を行い、イジングモデルのスピンσを出力する。
【０１１８】
図７を用いて説明したように、複数の擬似的スピンパルスＳＰｉの発振位相φｉ（ｔ）の暫定的な測定結果を用いた、各擬似的スピンパルスＳＰｉが関わる相互作用の大きさ及び符号のフィードバック実装を、具体的な構成を用いて行うことができる。そして、第２の実施形態によれば、第１の実施形態と比べて、擬似的スピンパルスＳＰｉ及び局部発振パルスＬＯがペアとなってリング共振器２を周回伝搬しないため、擬似的スピンパルスＳＰｉ及び局部発振パルスＬＯｉ間のクロストークをなくすことができ、全局部発振パルスＬＯｉを同相化させるパルス同相化部７をなくすことができる。よって、イジングモデルの量子計算装置Ｑにおいて、読み出しエラーを防止することができる。
【０１１９】
（本発明のイジングモデルの量子計算装置の計算結果）
数式４−７のレート方程式のシミュレーション結果として、２体の相互作用を含むイジングモデルの量子計算装置Ｑの計算結果及び時間発展を図８及び図９に示す。イジングモデルのハミルトニアンは、図８の上段に示されている。
【０１２０】
このハミルトニアンの最小エネルギー状態では、ある２つのサイトのスピンσが上向きであり、他の２つのサイトのスピンσが下向きであり、これらの最小エネルギー状態は６通り存在することが、数式４−７のレート方程式によらず分かっている。
【０１２１】
図８では、本発明の計算処理を１０００回繰り返した結果、６通りの最小エネルギー状態がほぼ等確率で導き出された一方、その他のエネルギー状態は導き出されなかった。つまり、エラーレートは、１０−３以下であることが分かった。
【０１２２】
図９では、ホモダイン検波結果の時間発展は、σ１＝σ２＝＋１、σ３＝σ４＝−１の最小エネルギー状態を指向している。初期状態から定常状態への所要時間は、〜４０×１００μｓであった。ここで、図９の横軸の時間は、共振器の光寿命（２ｋｍの光ファイバにより構成されるリング共振器２では１００μｓ）により規格化された時間である。
【０１２３】
数式４−７のレート方程式のシミュレーション結果として、４体の相互作用を含むイジングモデルの量子計算装置Ｑの計算結果及び時間発展を図１０及び図１１に示す。イジングモデルのハミルトニアンは、図１０の上段に示されている。
【０１２４】
このハミルトニアンの最小エネルギー状態では、ある１つのサイトのスピンσと他の３つのサイトのスピンσが異符号であり、これらの最小エネルギー状態は８通り存在することが、数式４−７のレート方程式によらず分かっている。
【０１２５】
図１０では、本発明の計算処理を１０００回繰り返した結果、８通りの最小エネルギー状態がほぼ等確率で導き出された一方、その他のエネルギー状態は導き出されなかった。つまり、エラーレートは、１０−３以下であることが分かった。
【０１２６】
図１１では、ホモダイン検波結果の時間発展は、σ１＝＋１、σ２＝σ３＝σ４＝−１の最小エネルギー状態を指向している。初期状態から定常状態への所要時間は、〜８０×１００μｓであった。ここで、図１１の横軸の時間は、共振器の光寿命（２ｋｍの光ファイバにより構成されるリング共振器２では１００μｓ）により規格化された時間である。
【０１２７】
ここで、計算処理の初期段階において、各擬似的スピンパルスＳＰの読出結果は正しい答えであるとは限らず、計算処理の初期段階において、正しくない答えをフィードバックすれば、計算処理の終期段階において、正しい答えを得られない、という疑義が生ずる。
【０１２８】
そこで、相互作用実装部５は、ある擬似的スピンパルスＳＰに対して注入される光の振幅を、計算処理の初期段階ほど大きく制御し、計算処理の終期段階ほど小さく制御する。つまり、相互作用実装部５は、数式４、５における比例定数Ａを、計算処理の初期段階ほど大きく制御し、計算処理の終期段階ほど小さく制御する。
【０１２９】
よって、計算処理の初期段階において、各擬似的スピンパルスＳＰの読出結果をできるだけ正しい答えにすることができる。たとえ、計算処理の初期段階において、ある擬似的スピンパルスＳＰの読出結果が正しい答えでないとしても、他の擬似的スピンパルスＳＰについての正しい答えをフィードバックすれば、計算処理の終期段階において、全擬似的スピンパルスＳＰについての正しい答えを得ることができる。つまり、イジングモデルの量子計算装置Ｑにおいて、読み出しエラーを防止することができる。
【０１３０】
計算処理の全段階で光注入強度を一定にした場合の時間発展を図１２に示す。図１２の場合には、定常状態において、不安定な発振状態が生じている。
【０１３１】
計算処理の初期段階に光注入強度を高くした場合の時間発展を図１３に示す。図１３の場合には、発振閾値直前において、光注入強度を減衰させることにより、定常状態において、不安定な発振状態がなくなっている。そして、初期段階において、光注入強度を高くすることにより、ある擬似的スピンパルスＳＰについて、初期段階においては、誤答が得られたとしても、終期段階においては、正答が得られるようになる。
【０１３２】
そして、数式１５のＡｉ（ｔ）＝Ａｊ（ｔ）を仮定することで、数式１８のイジングモデルの結合係数としてλｉ、Ｊｉｊ、Ｋｉｊｋを得ているところ、数式１５のＡｉ（ｔ）＝Ａｊ（ｔ）が成立しなければ、数式１８のイジングモデルの結合係数としてλｉ、Ｊｉｊ、Ｋｉｊｋと異なるλ’ｉ、Ｊ’ｉｊ、Ｋ’ｉｊｋが得られてしまう、という問題が生ずる。
【０１３３】
そこで、スピン測定部３、６は、複数の擬似的スピンパルスＳＰがリング共振器２を周回伝搬するたびに、複数の擬似的スピンパルスＳＰの振幅を暫定的に測定する。そして、パラメトリック発振器１は、スピン測定部３、６が測定した複数の擬似的スピンパルスＳＰの振幅に基づいて、複数の擬似的スピンパルスＳＰの振幅が等しくなるように、パラメトリック発振に用いるポンプパルスの振幅をフィードバック制御する。
【０１３４】
よって、各擬似的スピンパルスＳＰの振幅の不均衡による、イジングモデルの結合係数λｉ、Ｊｉｊ、Ｋｉｊｋの実質的な書き換えの問題をなくすことができる。つまり、イジングモデルの量子計算装置Ｑにおいて、読み出しエラーを防止することができる。
【０１３５】
複数の擬似的スピンパルスの強度を揃えない場合の時間発展を図１４に示す。図１４の場合には、各擬似的スピンパルスＳＰの強度によらず、各擬似的スピンパルスＳＰについて、ポンプパルスの強度が同様になるように制御されている。よって、各擬似的スピンパルスＳＰの強度が確実に揃えられるとは限らない。
【０１３６】
複数の擬似的スピンパルスの強度を揃えない場合の計算精度を図１５に示す。イジングエネルギー（折れ線グラフで表す。）が低い状態の存在確率（ヒストグラムで示す。）は、本来は高いはずであるところ、図１５の場合には、低くなっている。
【０１３７】
複数の擬似的スピンパルスの強度を揃えた場合の時間発展を図１６に示す。図１６の場合には、各擬似的スピンパルスＳＰの強度に応じて、各擬似的スピンパルスＳＰについて、ポンプパルスの強度がそれぞれに（同様でもよく、異なってもよい。）制御されている。よって、各擬似的スピンパルスＳＰの強度が確実に揃えられることができる。
【０１３８】
複数の擬似的スピンパルスの強度を揃えた場合の計算精度を図１７に示す。イジングエネルギー（折れ線グラフで表す。）が低い状態の存在確率（ヒストグラムで示す。）は、本来は高いはずであるところ、図１７の場合には、高くなっている。
【０１３９】
（本発明のイジングモデルの量子計算方法のループ手順）
本発明のイジングモデルの量子計算方法のループ手順を図１８、１９に示す。
【０１４０】
図１８に示したイジングモデルの量子計算方法のループ手順を実現するために、図４、６、７に示したイジングモデルの量子計算装置Ｑのように、暫定的スピン測定部３（又は擬似的スピン測定部６）、パラメトリック発振器１及び相互作用実装部５が、リング共振器２中の擬似的スピンパルスＳＰの周回伝搬方向にこの順序で配置される。
【０１４１】
パラメトリック発振器１は、複数の擬似的スピンパルスＳＰをパラメトリック発振させる（ステップＳ１）。暫定的スピン測定部３は、擬似的なスピンσを暫定的に測定する（ステップＳ２）。相互作用計算部４は、イジング相互作用を暫定的に計算する（ステップＳ３）。パラメトリック発振器１は、複数の擬似的スピンパルスＳＰをパラメトリック増幅させる（ステップＳ４）。相互作用実装部５は、イジング相互作用を暫定的に実装する（ステップＳ５）。ステップＳ３、Ｓ４は、順序が入れ替わることがある。
【０１４２】
複数の擬似的スピンパルスＳＰが定常状態に到達しなければ（ステップＳ６においてＮＯ）、ステップＳ２〜Ｓ６のループ手順を繰り返す。複数の擬似的スピンパルスＳＰが定常状態に到達したならば（ステップＳ６においてＹＥＳ）、擬似的スピン測定部６は、擬似的なスピンσを最終的に測定する（ステップＳ７）。
【０１４３】
このように、暫定的スピン測定ステップＳ２と相互作用実装ステップＳ５の間に、パラメトリック増幅ステップＳ４が入って、タイムラグがある程度生じるものの、イジングモデルのサイト間の実質には遅延のない相互作用を実装することができる。
【０１４４】
図１９に示したイジングモデルの量子計算方法のループ手順を実現するために、図４、６、７に示したイジングモデルの量子計算装置Ｑと異なり、パラメトリック発振器１、暫定的スピン測定部３（又は擬似的スピン測定部６）及び相互作用実装部５が、リング共振器２中の擬似的スピンパルスＳＰの周回伝搬方向にこの順序で配置される。
【０１４５】
パラメトリック発振器１は、複数の擬似的スピンパルスＳＰをパラメトリック発振させる（ステップＳ１１）。暫定的スピン測定部３は、擬似的なスピンσを暫定的に測定する（ステップＳ１２）。相互作用計算部４は、イジング相互作用を暫定的に計算する（ステップＳ１３）。相互作用実装部５は、イジング相互作用を暫定的に実装する（ステップＳ１４）。パラメトリック発振器１は、複数の擬似的スピンパルスＳＰをパラメトリック増幅させる（ステップＳ１５）。ステップＳ１３、Ｓ１５は、順序が入れ替わることがない。
【０１４６】
複数の擬似的スピンパルスＳＰが定常状態に到達しなければ（ステップＳ１６においてＮＯ）、ステップＳ１２〜Ｓ１６のループ手順を繰り返す。複数の擬似的スピンパルスＳＰが定常状態に到達したならば（ステップＳ１６においてＹＥＳ）、擬似的スピン測定部６は、擬似的なスピンσを最終的に測定する（ステップＳ１７）。
【０１４７】
このように、暫定的スピン測定ステップＳ１２と相互作用実装ステップＳ１４の間に、パラメトリック増幅ステップＳ１５が入らず、タイムラグがほとんど生じないため、イジングモデルのサイト間のほとんど遅延のない相互作用を実装することができる。
【０１４８】
（本発明のイジングモデルの量子並列計算装置の構成）
本発明のイジングモデルの量子並列計算装置Ｐの構成を図２０、２１に示す。
【０１４９】
イジングモデルの量子並列計算装置Ｐは、第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎを並列に備える。第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎとして、図４、６、７に示したイジングモデルの量子計算装置Ｑ及び図１８、１９に示したイジングモデルの量子計算方法を適用することができる。第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎでは、イジングモデルの複数のスピンσのうち、第１群、第２群、・・・、第ｎ群のスピンσに擬似的に対応する、第１群、第２群、・・・、第ｎ群の擬似的スピンパルスＳＰが、各リング共振器２を周回伝搬する。
【０１５０】
暫定的スピン共有部１０は、第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎが備える各暫定的スピン測定部３が暫定的に測定した、第１群、第２群、・・・、第ｎ群の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、並列に備わるイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの間で共有させる。
【０１５１】
図２０に示したイジングモデルの量子並列計算装置Ｐについて説明する。第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの各暫定的スピン測定部３は、第１群、第２群、・・・、第ｎ群の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、暫定的スピン共有部１０及び各相互作用計算部４に出力する。
【０１５２】
暫定的スピン共有部１０は、第１群より他の群（第２群、第ｎ群等）の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、第１番のイジングモデルの量子計算装置Ｑ１の相互作用計算部４に出力する。暫定的スピン共有部１０は、第２群より他の群（第１群、第ｎ群等）の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、第２番のイジングモデルの量子計算装置Ｑ２の相互作用計算部４に出力する。・・・暫定的スピン共有部１０は、第ｎ群より他の群（第１群、第２群等）の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、第ｎ番のイジングモデルの量子計算装置Ｑｎの相互作用計算部４に出力する。
【０１５３】
第１番のイジングモデルの量子計算装置Ｑ１の相互作用計算部４は、第１群の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第１群の擬似的スピンパルスＳＰが関わる相互作用を暫定的に計算する。第２番のイジングモデルの量子計算装置Ｑ２の相互作用計算部４は、第２群の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第２群の擬似的スピンパルスＳＰが関わる相互作用を暫定的に計算する。・・・第ｎ番のイジングモデルの量子計算装置Ｑｎの相互作用計算部４は、第ｎ群の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第ｎ群の擬似的スピンパルスＳＰが関わる相互作用を暫定的に計算する。
【０１５４】
このように、イジングモデルのサイト数が多いときでも、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎが並列分散処理を実行することにより、各々のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎは計算処理負担を軽減することができる。
【０１５５】
図２０に示したイジングモデルの量子並列計算装置Ｐでは、暫定的スピン共有部１０は、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの相互作用計算部４に、擬似的スピンパルスＳＰの擬似的なスピンσの情報を出力する。よって、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの相互作用計算部４は、擬似的スピンパルスＳＰの間の相互作用の計算処理負担を軽減されない。しかし、図４、６、７に示したイジングモデルの量子計算装置Ｑと比べて、図２０に示した各々のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎは、計算処理負担を軽減される。
【０１５６】
図２１に示したイジングモデルの量子並列計算装置Ｐについて説明する。第１番、第２番、・・・、第ｎ番のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの各暫定的スピン測定部３は、第１群、第２群、・・・、第ｎ群の擬似的スピンパルスＳＰの擬似的なスピンσの情報を、暫定的スピン共有部１０及び各相互作用計算部４に出力する。
【０１５７】
暫定的スピン共有部１０は、第１群と他の群（第２群、第ｎ群等）の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第１群と他の群の擬似的スピンパルスＳＰの間の相互作用の情報を、第１番のイジングモデルの量子計算装置Ｑ１の相互作用計算部４に出力する。暫定的スピン共有部１０は、第２群と他の群（第１群、第ｎ群等）の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第２群と他の群の擬似的スピンパルスＳＰの間の相互作用の情報を、第２番のイジングモデルの量子計算装置Ｑ２の相互作用計算部４に出力する。・・・暫定的スピン共有部１０は、第ｎ群と他の群（第１群、第２群等）の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第ｎ群と他の群の擬似的スピンパルスＳＰの間の相互作用の情報を、第ｎ番のイジングモデルの量子計算装置Ｑｎの相互作用計算部４に出力する。
【０１５８】
第１番のイジングモデルの量子計算装置Ｑ１の相互作用計算部４は、第１群に属する複数の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第１群に属する複数の擬似的スピンパルスＳＰの間の相互作用を暫定的に計算する。第２番のイジングモデルの量子計算装置Ｑ２の相互作用計算部４は、第２群に属する複数の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第２群に属する複数の擬似的スピンパルスＳＰの間の相互作用を暫定的に計算する。・・・第ｎ番のイジングモデルの量子計算装置Ｑｎの相互作用計算部４は、第ｎ群に属する複数の擬似的スピンパルスＳＰが関わるイジングモデルの結合係数を保持しており、第ｎ群に属する複数の擬似的スピンパルスＳＰの間の相互作用を暫定的に計算する。
【０１５９】
このように、イジングモデルのサイト数が多いときでも、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎが並列分散処理を実行することにより、各々のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎは計算処理負担を軽減することができる。
【０１６０】
図２１に示したイジングモデルの量子並列計算装置Ｐでは、暫定的スピン共有部１０は、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの相互作用計算部４に、擬似的スピンパルスＳＰの間の相互作用の情報を出力する。よって、複数のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎの相互作用計算部４は、擬似的スピンパルスＳＰの間の相互作用の計算処理負担を軽減される。つまり、図２０に示した各々のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎと比べて、図２１に示した各々のイジングモデルの量子計算装置Ｑ１、Ｑ２、・・・、Ｑｎは、計算処理負担を軽減される。
【産業上の利用可能性】
【０１６１】
本発明のイジングモデルの量子計算装置、イジングモデルの量子並列計算装置及びイジングモデルの量子計算方法は、イジングモデルにマッピングされるＮＰ完全問題などを高速かつ容易に解くのに適している。
【符号の説明】
【０１６２】
Ｖ、Ｖ１、Ｖ２、Ｖ３、Ｖ４：面発光レーザー
Ｍ：マスターレーザー
Ｉ１２、Ｉ１３、Ｉ１４、Ｉ２３、Ｉ２４、Ｉ３４：イジング相互作用実装部
Ｑ、Ｑ１、Ｑ２、Ｑｎ：イジングモデルの量子計算装置
ＳＰ、ＳＰ１、ＳＰ２、ＳＰ３、ＳＰ４、ＳＰ（Ｎ）：擬似的スピンパルス
ＬＯ、ＬＯ１〜ＬＯＮ：局部発振パルス
Ｐ：イジングモデルの量子並列計算装置
１：パラメトリック発振器
２：リング共振器
３：暫定的スピン測定部
４：相互作用計算部
５：相互作用実装部
６：擬似的スピン測定部
７：パルス同相化部
８：パルス発生器
９：第二高調波発生器
１０：暫定的スピン共有部
【図１】
2015156126.tif 000021
【図２】
2015156126.tif 000022
【図３】
2015156126.tif 000023
【図４】
2015156126.tif 000024
【図５】
2015156126.tif 000025
【図６】
2015156126.tif 000026
【図７】
2015156126.tif 000027
【図８】
2015156126.tif 000028
【図９】
2015156126.tif 000029
【図１０】
2015156126.tif 000030
【図１１】
2015156126.tif 000031
【図１２】
2015156126.tif 000032
【図１３】
2015156126.tif 000033
【図１４】
2015156126.tif 000034
【図１５】
2015156126.tif 000035
【図１６】
2015156126.tif 000036
【図１７】
2015156126.tif 000037
【図１８】
2015156126.tif 000038
【図１９】
2015156126.tif 000039
【図２０】
2015156126.tif 000040
【図２１】
2015156126.tif 000041
【国際調査報告】
2015156126.tif 000042
2015156126.tif 000043
2015156126.tif 000044
2015156126.tif 000045