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